Statistical framework for estimating GNSS bias

We present a statistical framework for estimating global navigation satellite system (GNSS) non-ionospheric differential time delay bias. The biases are estimated by examining differences of measured line integrated electron densities (TEC) that are scaled to equivalent vertical integrated densities. The spatio-temporal variability, instrumentation dependent errors, and errors due to inaccurate ionospheric altitude profile assumptions are modeled as structure functions. These structure functions determine how the TEC differences are weighted in the linear least-squares minimization procedure, which is used to produce the bias estimates. A method for automatic detection and removal of outlier measurements that do not fit into a model of receiver bias is also described. The same statistical framework can be used for a single receiver station, but it also scales to a large global network of receivers. In addition to the Global Positioning System (GPS), the method is also applicable to other dual frequency GNSS systems, such as GLONASS (Globalnaya Navigazionnaya Sputnikovaya Sistema). The use of the framework is demonstrated in practice through several examples. A specific implementation of the methods presented here are used to compute GPS receiver biases for measurements in the MIT Haystack Madrigal distributed database system. Results of the new algorithm are compared with the current MIT Haystack Observatory MAPGPS bias determination algorithm. The new method is found to produce estimates of receiver bias that have reduced day-to-day variability and more consistent coincident vertical TEC values.Comment: 18 pages, 5 figures, submitted to AM

Diurnal Variation of TEC and S 4 Index During the Period of Low Geomagnetic Activity at Ile-Ife, Nigeria

Ile-Ife lies on the equatorial anomaly region where the ionospheric current is greatly influenced by the existence of the equatorial electrojet. The dual frequency SCINDA NovAtel GSV 4004B GPS receiver recently installed at Ile-Ife [on geographical latitude 7°33′N and longitude 4°33′E and geomagnetic dipole (coordinate) of latitude 9.84°N and longitude 77.25°E] is currently operational and recording data from the available global positioning system satellites. The receiver provides the data on total electron content (TEC) and the scintillation index (S[subscript 4]). This paper presents the first sets of results from this station. Data records for the month of February 2010 were analyzed using the WinTec-P software program and these were interpreted to discuss the diurnal variation of the TEC and S[subscript 4] index during the period considered, as having low geomagnetic activity. The vertical TEC in this study showed that the values vary widely from as low as 0 TECu about sunrise to about 35 TECu during the day. Depletion in TEC was also noticed about sunset and marked by the occurrence of scintillations with a maximum index value of 0.3. Results of the IRI models and the observed TEC differ considerably; hence, there is the need to improve IRI models for its adaptability to the Africa ionospheric conditions

Hemispherical Shifted Symmetry in Polar Cap Patch Occurrence: A Survey of GPS TEC Maps From 2015–2018

Much theoretical and observational work has been devoted to studying the occurrence of F region polar cap patches in the Northern Hemisphere; considerably less work has been applied to the Southern Hemisphere. In recent years, the Madrigal database of mappings of total electron content (TEC) has improved in Southern Hemisphere coverage, to the point that we can now carry out a study of patch frequency and occurrence. We find that Southern Hemisphere patch occurrence is very similar to that of the Northern Hemisphere with a half‐year offset, plus an offset in universal time of approximately 12 hr. This is further supported by running an ionospheric model for both hemispheres and applying the same patch‐to‐background technique. Further, we present a simple physical mechanism involving a sunlit dayside plasma source concurrent with a dark polar cap, which yields a patch‐to‐background pattern very much like that seen in the TEC mappings for both hemispheres

Ionospheric symmetry caused by geomagnetic declination over North America

We describe variations in total electron content (TEC) in the North American sector exhibiting pronounced longitudinal progression and symmetry with respect to zero magnetic declination. Patterns were uncovered by applying an empirical orthogonal function (EOF) decomposition procedure to a 12 year ground-based American longitude sector GPS TEC data set. The first EOF mode describes overall average TEC, while the strong influence of geomagnetic declination on the midlatitude ionosphere is found in the second EOF mode (or the second most significant component). We find a high degree of correlation between spatial variations in the second EOF mode and vertical drifts driven by thermospheric zonal winds, along with well-organized temporal variation. Results strongly suggest a causative mechanism involving varying declination with longitude along with varying zonal wind climatology with local time, season, and solar cycle. This study highlights the efficiency and key role played by the geomagnetic field effect in influencing mesoscale ionospheric structures over a broad midlatitude range.National Science Foundation (U.S.) (Grant ATM-0733510)National Science Foundation (U.S.) (Grant ATM-0856093)National Science Foundation (U.S.) (Grant AGS-1242204)China Scholarship CouncilHaystack Observator

GPS TEC observations of dynamics of the mid‐latitude trough during substorms

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/95166/1/grl28288.pd

Ionospheric longitudinal variations at midlatitudes: Incoherent scatter radar observation at Millstone Hill

Incoherent scatter radar (ISR) extra-wide coverage experiments during the period of 1978–2011 at Millstone Hill are used to investigate longitudinal differences in electron density. This work is motivated by a recent finding of the US east-west coast difference in TEC suggesting a combined effect of changing geomagnetic declination and zonal winds. The current study provides strong supporting evidence of the longitudinal change and the plausible mechanism by examining the climatology of electron density Ne on both east and west sides of the radar with a longitude separation of up to 40o for different heights within 300–450 km. Main findings include: 1) The east-west difference can be up to 60% and varies over the course of the day, being positive (East side Ne > West side Ne) in the late evening, and negative (West side Ne > East side Ne) in the pre-noon. 2) The east-west difference exists throughout the year. The positive (relative) difference is most pronounced in winter; the negative (relative) difference is most pronounced in early spring and later summer. 3) The east-west difference tends to enhance toward decreasing solar activity, however, with some seasonal dependence; the enhancements in the positive and negative differences do not take place simultaneously. 4) Both times of largest positive and largest negative east-west differences in Ne are earlier in summer and later in winter. The two times differ by 12–13 h, which remains constant throughout the year. 5) Variations at different heights from 300–450 km are similar. Zonal wind climatology above Millstone Hill is found to be perfectly consistent with what is expected based on the electron density difference between the east and west sides of the site. The magnetic declination-zonal wind mechanism is true for other longitude sectors as well, and may be used to understand longitudinal variations elsewhere. It may also be used to derive thermospheric zonal winds.National Natural Science Foundation (China) (Grant 40890164)National Science Foundation (U.S.) (Grants ATM-0733510 and ATM- 6920184

Conjugate ionospheric perturbation during the 2017 solar eclipse

An edited version of this paper was published by AGU. Copyright 2021 American Geophysical Union.We report new findings of total electron content (TEC) perturbations in the southern hemisphere at conjugate locations to the northern eclipse on August 21, 2017. We identified a persistent conjugate TEC depletion by 10%–15% during the eclipse time, elongating along magnetic latitudes with at least ∼5° latitudinal width. As the Moon's shadow swept southward, this conjugate depletion moved northward and became most pronounced at lower magnetic latitudes (>−20°N). This depletion was coincident with a weakening of the southern crest of the equatorial ionization anomaly (EIA), while the northern EIA crest stayed almost undisturbed or was slightly enhanced. We suggest these conjugate perturbations were associated with dramatic eclipse initiated plasma pressure reductions in the flux tubes, with a large portion of shorter tubes located at low latitudes underneath the Moon's shadow. These short L-shell tubes intersect with the F region ionosphere at low and equatorial latitudes. The plasma pressure gradient was markedly skewed northward in the flux tubes at low and equatorial latitudes, as was the neutral pressure. These effects caused a general northward motion tendency for plasma within the flux tubes, and inhibited normal southward diffusion of equatorial fountain plasma into the southern EIA region. We also identified posteclipse ionospheric disturbances likely associated with the global propagation of eclipse-induced traveling atmospheric disturbances in alignment with the Moon's shadow moving direction

An Ionosphere Specification Technique Based on Data Ingestion Algorithm and Empirical Orthogonal Function Analysis Method

A data ingestion method in reproducing ionospheric electron density and total electron content (TEC) was developed to incorporate TEC products from the Madrigal Database into the NeQuick 2 model. The method is based on retrieving an appropriate global distribution of effective ionization parameter (Az) to drive the NeQuick 2 model, which can be implemented through minimizing the difference between the measured and modeled TEC at each grid in the local time‐modified dip latitude coordinates. The performance of this Madrigal TEC‐driven‐NeQuick 2 result is validated through the comparison with various International Global Navigation Satellite Systems Services global ionospheric maps and ionosonde data. The validation results show that a general accuracy improvement of 30–50% can be achieved after data ingestion. In addition, the empirical orthogonal function (EOF) analysis technique is used to construct a parameterized time‐varying global Az model. The quick convergence of EOF decomposition makes it possible to use the first six EOF series to represent over 90% of the total variances. The intrinsic diurnal variation and spatial distribution in the original data set can be well reflected by the constructed EOF base functions. The associated EOF coefficients can be expressed as a set of linear functions of F10.7 and Ap indices, combined with a series of trigonometric functions with annual/seasonal variation components. The NeQuick TEC driven by EOF‐modeled Az shows 10–15% improvement in accuracy over the standard ionosphere correction algorithm in the Galileo navigation system. These preliminary results demonstrate the effectiveness of the combined data ingestion and EOF modeling technique in improving the specifications of ionospheric density variations.Key PointsThe Madrigal TEC data are ingested into the NeQuick 2 model through deriving an effective ionization parameter (Az)The Empirical Orthogonal Function (EOF) analysis technique is used to construct a parameterized time‐varying Az model to make a predictionThe TEC data ingestion and EOF modeling are effective in bringing certain systematic improvement of ionosphere now‐cast/forecastPeer Reviewedhttps://deepblue.lib.umich.edu/bitstream/2027.42/146373/1/swe20760_am.pdfhttps://deepblue.lib.umich.edu/bitstream/2027.42/146373/2/swe20760.pd

Merging of Storm Time Midlatitude Traveling Ionospheric Disturbances and Equatorial Plasma Bubbles

Postsunset midlatitude traveling ionospheric disturbances (TIDs) and equatorial plasma bubbles (EPBs) were simultaneously observed over American sector during the geomagnetic storm on 8 September 2017. The characteristics of TIDs are analyzed by using a combination of the Millstone Hill incoherent scatter radar data and 2‐D detrended total electron content (TEC) from ground‐based Global Navigation Satellite System receivers. The main results associated with EPBs are as follows: (1) stream‐like structures of TEC depletion occurred simultaneously at geomagnetically conjugate points, (2) poleward extension of the TEC irregularities/depletions along the magnetic field lines, (3) severe equatorial and midlatitude electron density (Ne) bite outs observed by Defense Meteorological Satellite Program and Swarm satellites, and (4) enhancements of ionosphere F layer virtual height and vertical drifts observed by equatorial ionosondes near the EPBs initiation region. The stream‐like TEC depletions reached 46° magnetic latitudes that map to an apex altitude of 6,800 km over the magnetic equator using International Geomagnetic Reference Field. The formation of this extended density depletion structure is suggested to be due to the merging between the altitudinal/latitudinal extension of EPBs driven by strong prompt penetration electric field and midlatitude TIDs. Moreover, the poleward portion of the depletion/irregularity drifted westward and reached the equatorward boundary of the ionospheric main trough. This westward drift occurred at the same time as the sudden expansion of the convection pattern and could be attributed to the strong returning westward flow near the subauroral polarization stream region. Other possible mechanisms for the westward tilt are also discussed.Key PointsPostsunset EPBs driven by PPEF were observed to merge with midlatitude TIDs forming stream‐like depletion structures over American sectorDepletions reached 46 MLAT that map to 6,800 km over the equator and drifted westward reaching the equatorward boundary of the main troughStrong convection flow near SAPS region and disturbance thermospheric wind contributed to the westward drift of the midlatitude depletionsPeer Reviewedhttps://deepblue.lib.umich.edu/bitstream/2027.42/148412/1/swe20807.pdfhttps://deepblue.lib.umich.edu/bitstream/2027.42/148412/2/swe20807_am.pd

Coordinated Groundâ Based and Spaceâ Based Observations of Equatorial Plasma Bubbles

This paper presents coordinated and fortuitous groundâ based and spaceborne observations of equatorial plasma bubbles (EPBs) over the South American area on 24 October 2018, combining the following measurements: Globalâ scale Observations of Limb and Disk far ultraviolet emission images, Global Navigation Satellite System total electron content data, Swarm in situ plasma density observations, ionosonde virtual height and drift data, and cloud brightness temperature data. The new observations from the Globalâ scale Observations of Limb and Disk/ultraviolet imaging spectrograph taken at geostationary orbit provide a unique opportunity to image the evolution of plasma bubbles near the F peak height over a large geographic area from a fixed longitude location. The combined multiâ instrument measurements provide a more integrated and comprehensive way to study the morphological structure, development, and seeding mechanism of EPBs. The main results of this study are as follows: (1) The bubbles developed a westward tilted structure with 10â 15° inclination relative to the local geomagnetic field lines, with eastward drift velocity of 80â 120 m/s near the magnetic equator that gradually decreased with increasing altitude/latitude. (2) Waveâ like oscillations in the bottomside F layer and detrended total electron content were observed, which are probably due to upward propagating atmospheric gravity waves. The wavelength based on the mediumâ scale traveling ionospheric disturbance signature was consistent with the interbubble distance of â ¼500â 800 km. (3) The atmospheric gravity waves that originated from tropospheric convective zone are likely to play an important role in seeding the development of this equatorial EPBs event.Plain Language SummaryThis study presents multiâ instrument observations of equatorial plasma density depletions occurred on 24 October 2018 by using Globalâ scale Observations of Limb and Disk far ultraviolet images, Global Navigation Satellite System total electron content data, electron density measurements from Swarm satellite, ionosonde measurements, and cloud temperature data. This multiâ instrument study generated an integrated and detailed image revealing both largeâ scale and mesoscale structures of the equatorial plasma depletion. Our results also suggest that atmospheric gravity waves originating from tropospheric convection activity could play a significant seeding role in the development of equatorial plasma bubbles.Key PointsCombined GOLD/UV spectrograph images and groundâ based TEC data revealed EPB features and development over a large geographic areaBottomside F layer oscillations and traveling ionospheric disturbance were observed by ionosonde and detrended TEC resultsAtmospheric gravity waves likely play an important role in seeding the Râ T instability and the development of this EPB eventPeer Reviewedhttps://deepblue.lib.umich.edu/bitstream/2027.42/153570/1/jgra55456_am.pdfhttps://deepblue.lib.umich.edu/bitstream/2027.42/153570/2/jgra55456.pd