Total Electron Content Research Articles

Latitude Variation of the Post‐Sunset Plasma Density Enhancement During the Minor Geomagnetic Storm on 27 May 2021

AbstractIn this study, multiple instrumental observations including Global Navigation Satellite System total electron content (TEC), plasma drift velocity measured by Sanya (18.3°N, 109.6°E, dip latitude 12.6°N) Incoherent Scatter Radar (SYISR) and F2‐layer peak electron density (NmF2) and peak height (hmF2) from ionosonde and SYISR have been used to investigate ionospheric responses during a minor yet highly geo‐effective geomagnetic storm on 26–27 May 2021. Our findings revealed a significant time delay in the post‐sunset plasma density enhancement peak across different latitudes over East Asia, that is, the lower the geographic latitude, the earlier the peak appeared. The plasma density enhancement was accompanied by a decrease in hmF2 prior to NmF2 peak around sunset. The newly built SYISR measurements around sunset verified that the field‐aligned drift decreased the ionosphere with a notable time delay at latitude, beneficial to electron density enhancements at lower altitudes within the 16–30°N latitudinal band but a small TEC change. While at 30–50°N, it is possible that the competition between storm‐induced equatorward winds and downward field‐aligned drift depressed hmF2 decline and the buildup increased both NmF2 and TEC. The ICON observations suggested that the meridional wind during this minor storm event modulated the direction of plasma transport near sunset, playing a dominant role in post‐sunset plasma density enhancement from low to middle latitudes. These results provide fresh insight into the electrodynamic mechanisms of post‐sunset enhancements at middle and low latitudes over East Asia, and also enhance our understanding of the intricate behaviors within the ionosphere‐thermosphere system in response to a minor storm.

Thermospheric Wind Response to March 2023 Storm: Largest Wind Ever Observed With a Fabry‐Perot Interferometer in Tromsø, Norway Since 2009

AbstractSolar cycles 24–25 were quiet until a geomagnetic storm with a Sym‐H index of −170 nT occurred in late March 2023. On March 23–24, a Fabry‐Perot interferometer (FPI; 630 nm) in Tromsø, Norway, recorded the highest thermospheric wind speed of over 500 m/s since 2009. Comparisons with magnetometer readings in Scandinavia showed that a large amount of electromagnetic energy was transferred to the ionosphere‐thermosphere system. Total electron content maps suggested an enlarged auroral oval and revealed that the FPI observed winds near the polar cap instead of inside the oval for a long period during the storm main phase. The FPI wind had a strong equatorward component during the storm, likely because of the powerful anti‐sunward ionospheric plasma flow in the polar cap. The positive Y‐component of the IMF for 6 days before the storm caused a successive westward component of the FPI‐measured wind during the storm main phase. On March 24, the first day of the storm recovery phase, thermospheric wind disturbed and the ionospheric density decreased significantly at high latitudes. This density depression lasted for several days, and a large amount of electromagnetic energy during the storm modified the thermospheric dynamics and ionospheric plasma density.

Variation in the Reflection Height of VLF/LF Transmitter Signals in the D‐Region Ionosphere and the Possible Source: A 2018 Meteoroid in Hokkaido, Japan

AbstractSeveral studies have examined ionospheric variation associated with meteorites, meteoroids, or meteors based on Global Satellite Navigation System total electron content observations. However, there have been few quantitative studies of the D‐region of the ionosphere (60–90 km), which is associated with meteoroids. We investigated variation in the D‐region during the passage of a meteoroid over northeastern Hokkaido, Japan, at 11:55:55 UT on 18 October 2018, using very low‐frequency (VLF, 3–30 kHz) and low‐frequency (LF, 30–300 kHz) signals observed by three transmitters [JJY (40 kHz), JJY (60 kHz), and JJI (22.2 kHz)], at Rikubetsu, Japan. Periodic variation of 100–200 s was observed in the VLF and LF amplitudes upon arrival of the acoustic wave. The vertical seismic velocity of Hi‐net and F‐net data also showed acoustic waves. Although the main period of the acoustic wave was 0.1–0.5 s in the seismic data, a longer period component (100–200 s) remained during propagation up to the D‐region ionosphere. The estimated velocity of the acoustic waves was ∼340 m/s on the ground according to the Hi‐net seismic data. The acoustic wave originated near the endpoint (25 km altitude) of the meteoroid trajectory. Based on the observed propagation time of the acoustic waves and ray tracing results, the acoustic waves propagated obliquely from near the endpoint of the meteoroid trajectory up to a D‐region height (about ∼90 km altitude), south of the Rikubetsu receiver.

The Hunt for Perpendicular Magnetic Field Measurements in Plasma

The one consistent technique for remotely estimating the magnetic field in plasma has been Faraday rotation. It is only sensitive to the portion of the vector parallel to the propagation path. We show how to remotely detect the portion of the vector that is perpendicular using a modified measurement. Isolating this electromagnetic propagation wave mode to measure the magnetic field enables us to (i) study more about the magnetic field vector in plasma, (ii) reduce error in total electron content measurements, and (iii) discover new magnetic field information from archived data sets. The Appleton–Hartree equation is used to verify a new approach to calculating the phase change to an electromagnetic wave propagating through a plasma at frequencies larger than the gyrofrequency, the cyclotron frequency, and the upper hybrid frequency. Focusing on the perpendicular propagation modes, the simplified equation for the integrated path effect from a perpendicular magnetic field is calculated. The direction of the perpendicular component is unknown, because the magnetic field is squared. Isolating the magnetic term in the equation with dual frequency waves is shown. We also show how to eliminate the magnetic field contribution to total electron content measurements with a similar approach. In combination with Faraday rotation, the degeneracy of the magnetic field vector direction is reduced to a cone configuration.

Seasonal Features of the Ionospheric Total Electron Content Response at Low Latitudes during Three Selected Geomagnetic Storms

In the present paper, the response of the ionospheric Total Electron Content (TEC) at low latitudes during several geomagnetic storms occurring in different seasons of the year is investigated. In the analysis of the ionospheric response, the following three geomagnetic events were selected: (i) 23–24 April 2023; (ii) 22–24 June 2015 and (iii) 16 December 2006. Global TEC data were used, with geographic coordinates recalculated with Rawer’s modified dip (modip) latitude, which improved the accuracy of the representation of the ionospheric response at low and mid-latitudes. By decomposition of the zonal distribution of the relative deviation of the TEC values from the hourly medians, the spatial distribution of the anomalies, the dependence of the response on the local time and their evolution during the selected events were analyzed. As a result of the study, it was found that the positive response (i.e., an increase in electron density relative to quiet conditions) in low latitudes occurs at the modip latitudes 30° N and 30° S. An innovative result related to the observed responses during the considered events is that they turn out to be relatively stationary. The longitude variation in the observed maxima changes insignificantly during the storms. Depending on the season, there is an asymmetry between the two hemispheres, which can be explained by the differences in the meridional neutral circulation in different seasons.

SHARAD Mapping of Mars Dayside Ionosphere Patterns: Relationship to Regional Geology and the Magnetic Field

AbstractThe electron density of the Martian ionosphere is modulated by solar wind forcing and crustal magnetic fields. Sounding observations from the orbital Shallow Radar (SHARAD) map the ionospheric total electron content (TEC) at a spatial resolution of ∼35 km. Averaging over a 250‐km diameter window from data collected weeks to years apart yields the first map of long‐term stable dayside martian TEC features. An extensive region of suppressed TEC in the southern highlands correlates with strong radial magnetic fields, but in other areas no simple correlation is observed. The TEC maps do follow the outlines of exposed Noachian crust and patches of magnetization in Tharsis not reset by volcanic activity. SHARAD TEC mapping may capture magnetic field strength at an intermediate height between the surface and the altitudes of orbital measurements underlying spherical harmonic models. Existing and future data will allow SHARAD TEC mapping to ∼100 km spatial resolution.

Investigating the Global Performance of the BDS-2 and BDS-3 Joint Real-Time Undifferenced and Uncombined Precise Point Positioning Using RTS Products from Different Analysis Centers

Compared to the traditional ionospheric-free (IF) precise point positioning (PPP) model, the undifferenced and uncombined (UU) PPP has the advantages of lower observation noise and the ability to obtain ionospheric information. Thanks to the IGS (International GNSS Service), real-time service (RTS) can provide RT vertical total electron content (VTEC) products, and an enhanced RT UU-PPP based on the RT-VTEC constraints can be achieved. The global performance of the BeiDou Navigation Satellite System-2 (BDS-2) and BDS-3 joint RT UU-PPP using different RTS products was investigated. There is not much difference in the RTS orbit accuracy of medium earth orbit (MEO) satellites among all analysis centers (ACs), and the optimal orbit accuracy is better than 5, 9, and 7 cm in the radial, along-track, and cross-track directions, respectively. The orbit accuracy of inclined geosynchronous orbit (IGSO) satellites is worse than that of MEO satellites. Except for CAS of 0.46 ns, the RTS clock accuracy of MEO satellites for other ACs achieves 0.2–0.27 ns, and the corresponding accuracy is about 0.4 ns for IGSO satellites. In static positioning, due to the limited accuracy of RT-VTEC, the convergence time of the enhanced RT UU-PPP is longer than that of RT IF-PPP for most ACs and can be better than 25 and 20 min in the horizontal and vertical components, respectively. After convergence, the 3D positioning accuracy of the static RT UU-PPP is improved by no more than 8.7%, and the optimal horizontal and vertical positioning accuracy reaches 3.5 and 7.0 cm, respectively. As for the kinematic mode with poor convergence performance, with the introduction of RT-VTEC constraints, the convergence time of RT UU-PPP can be slightly shorter and reaches about 55 and 60 min in the horizontal and vertical components, respectively. Both the horizontal and vertical positioning accuracies of the kinematic RT UU-PPP can be improved and achieve around 7.5 and 10 cm, respectively.

Characteristic analysis of the differences between total electron content (TEC) values in global ionosphere map (GIM) grids

Abstract. Using total electron content (TEC) from a global ionosphere map (GIM) for ionospheric delay correction is a common method of eliminating ionospheric errors in satellite navigation and positioning. On this basis, the TEC of a puncture point can be obtained by GIM grid TEC interpolation. However, in terms of grid, only few studies have analyzed the TEC value size characteristics of its four grid points, that is, the TEC difference characteristics among them. In view of this, by utilizing the GIM data from high solar-activity years (2014) and low solar-activity years (2021) provided by CODE (Center for Orbit Determination in Europe), this paper proposes the grid TEC difference as a way of analyzing TEC variation characteristics within the grid, which is conducive to exploring and analyzing the variation characteristics of the ionosphere TEC in the single-station area. The value is larger in high solar-activity years and generally small in low solar-activity years, and the value of high-latitude areas is always smaller than that of low-latitude areas. Specifically, in high solar-activity years, most of the GIM grid TEC internal differences are within 4 TECu (1 TECu = 1016 electrons m−2) in high-latitude and midlatitude regions, while only 78.17 % are in low-latitude regions. In low solar-activity years, the TEC difference values within a GIM grid are mostly less than 2 TECu, and most of them in the high and middle latitudes are within 1 TECu. The main finding of this analysis is that the grid TEC differences are small for most GIM grids, especially in the midlatitudes to high latitudes of low solar years. This means that relevant extraction methods and processes can be simplified when TEC within these GIM grids is needed.

Geomagnetic storm effect on equatorial ionosphere over Sri Lanka through total electron content observations from continuously operating reference stations network during Mar–Apr 2022

Abstract The technological advancements in the current era have highlighted the increasing significance of satellite-based positioning, navigation, and timing services in a wide range of dynamic and critical applications. This has led to significant efforts towards enhancing the performance of global navigation satellite systems (GNSS) operating under challenging ionospheric conditions. The Sri Lankan ionosphere region is a focal point of equatorial aeronomy scientists, being situated in the southernmost landmass of the Indian longitude sector within the vicinity of the magnetic equator where a combination of electric, wind, and temperature dynamics exerts a substantial influence on the ionosphere but was relatively unexplored in the past. In the present work, we employed a Kriging interpolation technique on the total electron content (TEC) variables from ten GNSS receivers operating under the Continuously Operating Reference Stations (CORS) network in Sri Lanka first ever of its kind to deliver two-dimensional regional ionospheric TEC maps at hourly intervals, both during quiet and disturbed ionospheric conditions in the equinoctial March and April months of 2022. The latitudinal variation patterns are discernable from the hourly TEC maps. Furthermore, a comparative analysis of the performance of GNSS-derived TEC with that of the routinely published Global Ionospheric Maps (GIMs) confirms overestimation characteristics of the latter irrespective of the local time of observation. The generated regional ionospheric maps are fairly responsive to the onset of the storm and the recovery phase thereafter. The extent of nighttime ionospheric irregularity is also probed through the rate of TEC index (ROTI) variations, demonstrating that the irregularities were insignificant during the selected storm event.

Response of the lower and upper ionosphere after the eruption of Shiveluch volcano on april 10, 2023

The disturbances in the lower ionosphere and in the region of the maximum of the ionospheric F2 layer during the Shiveluch volcanic eruption in April 2023 are analyzed based on data from ground-based magnetometers and GPS radio sounding of the ionosphere. The magnetic stations were located at distances of 455 km (Paratunka) and 752 km (Magadan) from the volcano. The variations in the magnetic field and total electron content of the ionosphere were studied as characteristics of the ionospheric response to this event. An analysis of the measurements showed that the impact on the ionosphere is carried out by seismic Rayleigh waves and atmospheric acoustic-gravity waves generated by volcanic explosions. The energy of several explosions was estimated from the amplitude of the ionospheric signal in the total electron content.

The Analysis of Ionospheric TEC Anomalies Prior to the Jiuzhaigou Ms7.0 Earthquake Based on BeiDou GEO Satellite Data

The position between BeiDou geostationary Earth orbit (GEO) satellites and ground-based receiving stations can roughly be considered to be constant with negligible fluctuations; thus, the total electron content (TEC) data over a fixed ionospheric piercing point (IPP) can be continuously acquired, which is advantageous for monitoring ionospheric disturbances. Focused on the Jiuzhaigou Ms7.0 earthquake that occurred on 8 August 2017, the TEC data inverted by the BeiDou GEO satellite were analyzed to extract ionospheric disturbances potentially associated with the earthquake. It was found that significant anomalies in ionospheric TEC occurred 10–11 days, 6–7 days, and 1–9 h prior to the earthquake, which was mainly located in the southeast and southwest directions within about 2500 km distance from the epicenter. Comparing the spatial and temporal characteristics between the ionospheric disturbance and the radon gas near the surface, the atmospheric electric field, and the spectrum of TEC data, it was considered that the chemical and acoustic–gravity wave pathway may play an important role in the lithosphere–atmosphere–ionosphere coupling (LAIC) mechanism.

Relationship between GIX, SIDX, and ROTI ionospheric indices and GNSS precise positioning results under geomagnetic storms

Ionospheric indices give information about ionospheric perturbations, which may cause absorption, diffraction, refraction, and scattering of radio signals, including those from global navigation satellite systems (GNSS). Therefore, there may be a relationship between index values and GNSS positioning results. A thorough understanding of ionospheric indices and their relationship to positioning results can help monitor and forecast the reliability and accuracy of GNSS positioning and support the precision and safety of life applications. In this study, we present the relationship between three indices: Gradient Ionosphere indeX (GIX), Sudden Ionospheric Disturbance indeX (SIDX), and Rate of Total electron content Index (ROTI) in relation to precise positioning results. We used two approaches: precise point positioning (PPP) and relative positioning for long baselines. We focus on GNSS stations located in Europe for two selected geomagnetic storms: March 17, 2015, and May 22, 2015. Our results show that in the case of PPP, positioning degradation occurred mainly at high latitudes and was mostly caused by rapid small-scale changes in ionospheric electron content represented by SIDX and ROTI. We also showed a significant correlation between cycle slips of GNSS signals and ROTI (0.88). The most significant degradations for relative positioning for low and medium latitudes were associated with large spatial gradients reflected by the GIX.

Geomagnetic Storm Variation of Vertical Total Electron Content (VTEC) Over Some Euro-African Stations

Geomagnetic storms are events which have physical effects on some ionospheric parameters that, to some extent,affects the state and dynamics of the ionosphere with important implications on GNSS applications. Here, the total electron content (TEC) of Brussels (50.80oN, 04.37oE), Madrid (40.43oN, 04.25oW) and Irkutsk (52.22oN, 104.32oE), which are all mid-latitude European stations are compared with Libreville (00.35oN, 09.67oE) and Lusaka (15.43oS, 28.32oE) which are equatorial and low-latitude stations respectively. This study is done over two geomagnetic storms that took place in the solstice period of 2004. Deviations of storm time VTEC from solar quiet (Sq) averages are calculated, analysed and presented. Similarities and differences of storm effects are observed in the European stations with enhancements and depressions. Diurnal solar quiet day variations showed high VTEC during the post-noon hours for all the stations. The VTEC deviations during storm time at Libreville lie within, for Lusaka it is. For the mid-latitude European stations, the deviations are lower such that is recorded at Brussels while is recorded for both Irkutsk and Madrid. Enhancement of VTEC during the daytime storm period is attributable to the super-fountain effect caused by the prompt penetration electric fields (PPEFs) into the ionosphere and magnetosphere while low VTEC at night-time is attributed to the process of recombination. Understanding the behaviour of the ionosphere during geomagnetic storms is important and necessary for a better understanding of the applications of GNSS.

Prediction of ionospheric total electron content data using NARX neural network model

Successful prediction of ionospheric total electron content (TEC) data will help in correction of positioning errors in global navigation satellite systems (GNSS) caused by the ionosphere. This research paper proposes a prediction model for ionospheric TEC using a nonlinear autoregressive with exogenous inputs (NARX) neural network that utilizes past TEC data alongwith solar and geomagnetic indices namely F10.7, disturbed storm (Dst), Kp, Ap, and time of the day. We assess the prediction capability of our model at different latitudes during different solar activity years. We compare our results with another NARX model which uses previous TEC data along with time of the day, day of the year and season as exogenous parameters. The results show that for the solar minimum year the TEC prediction accuracy improves by 35.71% and for the solar maximum year it improves by 31.20%. The results using root mean square error (RMSE), mean absolute error (MAE), correlation coefficient, and symmetric mean absolute percentage error (sMAPE) clearly indicate that solar and geomagnetic indices along with time of the day help in enhancing prediction accuracy of TEC across different latitudinal regions during both solar minimum and maximum years.

Localized Plasma Density Peak at Middle Latitudes During the April 2023 Geomagnetic Storm

AbstractThis paper conducts a multi‐instrument analysis of a latitudinal plasma density peak at the middle latitudes during the early recovery phase of the April 2023 geomagnetic storm. The total electron content (TEC), peak density of the F layer, and the in situ plasma density from Swarm and Defense Meteorological Satellite Program (DMSP) satellites all capture this peak feature. This narrow latitudinal peak structure appeared around 50°N and extended from 40°E to 150°E in longitude with a prolonged duration of about 8 hr from sunset to midnight. This mid‐latitude peak reveals a noticeable equatorward motion and a slight westward shift. According to the plasma composition observations from DMSP satellites, this peak structure shows an O+ ions dominance, which means that this peak is more likely to be formed by an internal rather than an external source from the plasmasphere. Meanwhile, the middle latitude Fabry–Perot interferometer (FPI) observed strong equatorward thermospheric winds, and the peak height of the F layer presented a visible elevation, which suggests that the equatorward wind lifting caused the plasma density enhancement. Besides, the O/N2 ratio significantly decreased at lower and middle latitudes, and ion drift observations showed a distinct subauroral westward channel. Based on these simultaneous measurements, this structure's sharp equatorward and poleward boundaries might be related to the O/N2 ratio change and the subauroral polarization stream (SAPS) flow separately.

Ionospheric Disturbances Generated by the 2015 Calbuco Eruption: Comparison of GITM‐R Simulations and GNSS Observations

AbstractVolcanic eruptions provide broad spectral forcing to the atmosphere and understanding the primary mechanisms that are relevant to explain the variety in waveform characteristics in the Ionosphere‐Thermosphere (IT) is still an important open question for the community. In this study, Global Navigation Satellite System (GNSS) Total Electron Content (TEC) data are analyzed and compared to simulations performed by the Global Ionosphere‐Thermosphere Model with Local Mesh Refinement (GITM‐R) for the first phase of the 2015 Calbuco eruption that occurred on 22 April. A simplified source representation and spectral acoustic‐gravity wave (AGW) propagation model are used to specify the perturbation at the lower boundary of GITM‐R at 100 km altitude. Two assumptions on the propagation structure, Direct Spherical (DS) and Ground Coupled (GC), are compared to the GNSS data and these modeling specifications show good agreement with different aspects of the observations for some waveform characteristics. Most notably, GITM‐R is able to reproduce the relative wave amplitude of AGWs as a function of radial distance from the vent, showing acoustic dominant forcing in the near field (<500 km) and gravity dominant forcing in the far‐field (>500 km). The estimated apparent phase speeds from GITM‐R simulations are consistent with observations with ∼10% difference from observation for both acoustic wave packets and a trailing gravity mode. The relevance of the simplifications made in the lower atmosphere to the simulated IT response is then discussed.

Controlling Factors of the Seasonal Variation of the Latitudinal Location of the Equatorial Ionization Anomaly Crest

AbstractThe latitudinal location of the Equatorial Ionization Anomaly (EIA) crest has seasonal variation, and there are disagreements on the interpretation of such seasonal characteristic in previous studies. Some studies suggested that this seasonal characteristic is determined by the seasonal characteristic of the equatorial electric field. Others suggested that this seasonal characteristic is determined by the seasonal changes of the thermospheric wind. The current paper uses Total Electron Content (TEC) data and the Thermosphere Ionosphere Electrodynamics General Circulation Model (TIEGCM) to analyze the seasonal variation of the northern EIA crest in the eastern Asian sector under low solar activity. Our results show that the monthly averaged latitudinal location of the northern EIA crest has a good linear relationship (r = 0.74) with the monthly averaged Equatorial Electrojet (EEJ) intensity, which is a commonly used proxy of the eastward electric field. However, TIEGCM simulations with and without F‐region wind indicate that such a relationship might be attributed to wind effects. Additionally, the linear relationship between the EEJ intensity and the northern‐southern EIA crest distance is not significant (r = 0.47) in the eastern Asian sector. Our results suggest that a good correspondence between the eastward electric field and the latitudinal location of the EIA crest is not assured annually, as the seasonally varying F‐region wind significantly influences EIA evolution.

A novel ionospheric TEC mapping function with azimuth parameters and its application to the Chinese region

The ionospheric mapping function (MF) for Global Navigation Satellite System (GNSS), a mutual projection method for the slant total electron content (STEC) and vertical total electron content, is one of the significant factors affecting the performance of ionospheric models. The commonly used MF assumes isotropic TEC variations and takes into account only the satellite elevation angle, which may result in significant ionospheric projection errors, especially at low elevation angles. Based on the single-layer model, we propose an additional azimuth parameter mapping function (APMF). The APMF was estimated and evaluated by the NeQuick model during the periods of January 2014 and January 2022 from the aspect of simulation and measured STEC during the periods of 2014 and 2022 from the aspect of actual measurements over China, respectively. Compared to the modified single-layer model mapping function (MSLM-MF), the experimental results indicate that (1) The APMF can significantly reduce the ionospheric projection error, and the fluctuation in errors with different azimuth angles is small. (2) According to the evaluation based on the NeQuick simulation during the TEC peak time, when the ionosphere is quite active, the upper and lower quartiles of the absolute projection error boxplot of the APMF relative to the MSLM-MF in January 2014 are reduced by 56.1% and 60.0%, respectively, and in January 2022, they are reduced by 67.7% and 65.2%, respectively. Similarly, the upper whiskers in the boxplot are reduced by 54.7% and 67.5% in January 2014 and January 2022, respectively; the APMF performance in terms of the root mean square error (RMSE) is improved by 47.0% in January 2014 and 58.3% in January 2022. (3) According to the evaluation based on the measured STEC from GNSS raw data during the TEC peak time, the upper and lower quartiles of the absolute mapping error boxplot of the APMF relative to the MSLM-MF in 2014 are reduced by 48.9% and 46.9%, respectively, while in 2022, they are reduced by 48.3% and 41.2%, respectively. The upper whiskers in the boxplot are reduced by 41.8% and 35.2% in 2014 and 2022, respectively; the APMF performance in terms of RMSE is improved by 44.6% in 2014 and 39.2% in 2022.

3‐D Ionospheric Imaging Over the South American Region With a New TEC‐Based Ionospheric Data Assimilation System (TIDAS‐SA)

AbstractThis study has developed a new TEC‐based ionospheric data assimilation system for 3‐D regional ionospheric imaging over the South American sector (TIDAS‐SA) (45°S–15°N, 35°–85°W, and 100–800 km). The TIDAS‐SA data assimilation system utilizes a hybrid Ensemble‐Variational approach to incorporate a diverse set of ionospheric data sources, including dense ground‐based Global Navigation Satellite System (GNSS) line‐of‐sight Total Electron Content (TEC) data, radio occultation data from the Constellation Observing System for Meteorology, Ionosphere, and Climate‐2 (COSMIC‐2), and altimeter TEC data from the JASON‐3 satellite. TIDAS‐SA can produce a reanalyzed three‐dimensional (3‐D) electron density spatial variation with a high time cadence, yielding spatial‐temporal resolution of 1° (latitude) × 1° (longitude) × 20 km (altitude) × 5 min. This allows us to reconstruct and study the 3‐D ionospheric morphology with multi‐scale structures. The performance of the data assimilation system is validated against independent ionosonde and in situ measurements through an experiment for a strong geomagnetic storm event on 03–04 November 2021. The results demonstrate that TIDAS‐SA can provide detailed and altitude‐resolved information that accurately characterizes the storm‐time ionospheric disturbances in vertical and horizontal domains over the equatorial and low‐latitude regions of South America.

Relativistic Electron Precipitation Driven by Mesoscale Transients, Inferred From Ground and Multi‐Spacecraft Platforms

AbstractPrecipitation of relativistic electrons into the Earth's atmosphere regulates the outer radiation belt fluxes and contributes to magnetosphere‐atmosphere coupling. One of the main drivers of such precipitation is electron scattering by whistler‐mode waves. Such waves typically originate at the equator, where they can resonate with and scatter sub‐relativistic (tens to a few hundred keV) electrons. However, they can occasionally propagate far away from the equator along field lines, reaching middle latitudes, where they can resonate with and scatter relativistic (>500 keV) electrons. Such a propagation is typical for the dayside, but statistically has not been found on the nightside where the waves are quickly damped along their propagation due to Landau damping. Here we explore two events of relativistic electron precipitation from low‐altitude observations on the nightside. Combining measurements of whistler‐mode waves from ground observatories, relativistic electron precipitation from low‐altitude satellites, total electron content maps from GPS receivers, and magnetic field and electron flux from equatorial satellites, we show wave ducting by plasma density gradients is the possible channel that allows the waves to reach middle latitudes and scatter relativistic electrons. We suggest that both whistler‐mode wave generation and ducting can be driven by equatorial mesoscale (with spatial scales of about one Earth radius) transient structures during nightside injections. We also compare these nightside events with observations of ducted waves and relativistic electron precipitation at the dayside, where wave generation and ducting are driven by ultra‐low‐frequency waves. This study demonstrates the potential importance of mesoscale transients in relativistic electron precipitation, but does not however unequivocally establish that ducted whistler‐mode waves are the primary cause of the observed electron precipitation.