Total Electron Content Response Research Articles

Response of total electron content to the October 25, 2022 partial solar eclipse from high to low latitudes in the Euro-Asian region

Solar eclipses (SEs) pertain to high-energy sources. They are capable of causing significant disturbances in all geoshells and geophysical fields. Despite the almost century-long history of studying the influence of eclipses on geoshells, the task of observing and analyzing disturbances in them remains relevant. This is explained by the fact that the response of the media significantly depends on the state of atmospheric and space weather, the location and time of observation, the SE magnitude, etc. The aim of this paper is to present the results of a statistical analysis of the features of the total electron content (TEC) variations during the October 25, 2022 SE from high to low latitudes in the Euro-Asian region. Data from Global Navigation Satellite System were used for the analysis. The error in estimating TEC in the vertical column did not exceed 0.1 TECU. For the first time, the response of the ionosphere to an SE has been studied on a global scale (from the Reykjavik station to the Novosibirsk station) using methods of statistical analysis, including the time moments after sunset. The decrease in TEC reached 1.0–19.0 TECU, and the relative decrease was approximately –0.32. Changes in all TEC disturbance parameters depended in a complex manner on the relative obscuration area of the solar disk. This can be explained by non-monotonic changes in electron density during the daytime, by the intensification of the ionosphere–plasmasphere coupling, the disturbance of the ionospheric current system, vertical plasma drift in crossed electric and magnetic fields, as well as the generation of wave disturbances during the SE period. The delay time of the maximum ionospheric response (maximum |δV| value) with respect to the maximum magnitude of the eclipse was 15–25 min. The duration of the ionospheric response to the SE slightly exceeded the duration of the perturbing source influence.

Different Response of the Ionospheric TEC and EEJ to Ultra‐Fast Kelvin Waves in the Mesosphere and Lower Thermosphere

AbstractWe studied the response of ionospheric total electron content (TEC) and equatorial electrojet (EEJ) to the ultra‐fast Kelvin wave (UFKW) at the equator in the mesosphere using zonal wind data obtained from TIMED Doppler Interferometer (TIDI), EEJ data over the monitoring station Jicamarca (12°S, 77°W) and global TEC maps. The least squares fitting method is utilized to perform a spectral analysis of zonal wind, EEJ and TEC. Our analysis results demonstrate that UFKW events can be divided into four categories: (a) UFKW events with both TEC and EEJ response; (b) UFKW events with TEC response but without EEJ response; (c) UFKW events with EEJ response but without TEC response; (d) UFKW events without neither TEC response nor EEJ response. The first type of UFKW events occur the most often and is generally thought to generate a response in EEJ at approximately 105–110 km through the dynamo effect. The polarization electric field associated with EEJ then produces a response in the ionospheric TEC through the fountain effect. The lack of EEJ response in the second type of UFKWs may be due to the influence of eastward background winds. We found that all UFKW events with EEJ response have a response in TEC. The fourth type of UFKWs have smaller amplitudes, shorter vertical wavelengths and longer periods, which make them more likely to dissipate and cannot propagate to higher altitudes. These UFKWs cannot propagate to the altitude of EEJ and produce a response in EEJ, much less in TEC.

2‐D Total Electron Content and 3‐D Ionospheric Electron Density Variations During the 14 October 2023 Annular Solar Eclipse

AbstractThis study investigates the ionospheric total electron content (TEC) responses in the 2‐D spatial domain and electron density variations in the 3‐D spatial domain during the annular solar eclipse on 14 October 2023, using ground‐based Global Navigation Satellite System (GNSS) observations, a novel TEC‐based ionospheric data assimilation system (TIDAS), ionosonde measurements, and satellite in situ data. The main results are summarized as follows: (a) The 2‐D TEC responses exhibited distinct latitudinal differences. The mid‐latitude ionosphere exhibited a more substantial TEC decrease of 25%–40% along with an extended recovery time of 3–4 hr. In contrast, the equatorial and low‐latitude ionosphere experienced a smaller TEC reduction of 10%–25% and a faster recovery time of 20–50 min. The minimal eclipse effect was observed near the northern equatorial ionization anomaly crest region. (b) The ionospheric electron density variations during the eclipse were effectively reconstructed by TIDAS data assimilation in the 3‐D domain, providing important altitude information with validity. (c) The ionospheric electron density variations showed a notable altitude‐dependent feature. The eclipse led to a substantial electron density reduction of 30%–50%, with the maximum depletion occurring around the ionospheric F2‐layer peak height (hmF2) of 250–350 km. The post‐eclipse recovery of electron density exhibited a relatively slower pace near the F2‐layer peak height than that at lower and higher altitudes.

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.

Polar Vortex Strength Impacts on the Longitudinal Structure of Thermospheric Composition and Ionospheric Electron Density

AbstractThis work examines the weak stratospheric polar vortex event in late February 2008 and the strong polar vortex event in late February 2021 to understand the potential influence of polar vortex strength on midlatitude thermospheric composition and ionospheric plasma density. The weak polar vortex event shows a reduction in thermospheric O/N2 at all longitudes, changes in the ionospheric total electron content (TEC) that vary with longitude, and an intensification of the semidiurnal tides. The strong polar vortex event shows O/N2 elevated by up to 25% in the Asian sector and a 25% decrease in semidiurnal tidal amplitude. The TEC response during the strong polar vortex event is complex and shows enhanced TEC in the Asian sector where O/N2 is elevated. These are the first observations of anomalies in the thermospheric composition and ionospheric plasma associated with a strong polar vortex event.

The 6 February 2023 Türkiye Earthquake Sequence as Detected in the Ionosphere

AbstractOn 6 February 2023, a series of large earthquakes struck Turkey and Northern Syria. The main earthquake of Mw 7.8 occurred at 01:17:34 UTC and was followed by the three notable (Mw > 5.5) aftershocks within the next 18 min. Then, ∼9 hr later, the biggest aftershock with magnitude Mw 7.5 and a Mw 6.0 earthquake occurred to the north‐east from the first main earthquake. In this work, we use data of ground‐based Global Navigation Satellite Systems (GNSS) receivers in Turkey, Israel and Cyprus to analyze the ionospheric response to this series of earthquakes. We separate these events in two groups: the first sequence of earthquakes (at 01–02 UTC) and the second sequence (at 10–11 UTC). For the first sequence, we observe a clear N‐shaped total electron content (TEC) response after the Mw 7.8 mainshock earthquake and Mw 6.7 aftershock, and a smaller TEC disturbance that is, most likely, caused by the Mw 5.6 earthquake. The latter is now the smallest earthquake detected by using ionospheric GNSS data. The co‐seismic ionospheric disturbances (CSID) propagated from the epicentral area in the south‐west direction with velocities of about 750–830 m/s. For the second sequence, we observed the response to the Mw 7.5 aftershock earthquake and the Mw 6.0 aftershock. The CSID propagated both to the south‐west and the north‐west to the epicentral area, with velocities of about 950–1,100 m/s.

Analysis of the Ionospheric Response to Sudden Stratospheric Warming and Geomagnetic Forcing over Europe during February and March 2023

A study of the behavior of the main characteristics of the ionosphere over Europe during the 26–28 February 2023 ionospheric storm was carried out in this present work. The additional influence of sudden stratospheric warming on the ionosphere was considered. The behavior of the critical frequency of the ionosphere foF2 (characterizing the maximum electron density), the peak height of the F2-layer (hmF2), and Total Electron Content (TEC) were investigated through their relative deviations from the quiet conditions. The behavior of the TEC over Europe showed the geographic latitudinal dependence of the response. The variability in the ionospheric critical frequency was represented by the data of 10 ionospheric stations for vertical sounding located in two groups: (i) near the prime meridian and (ii) near the 25° E meridian. Some differences were found in the response compared to the TEC response, which was explained by the different responses of the top maximum region and bottom maximum region. The peak height of the F2 layer varied strongly during the storm, which was due to the forced drift of ionospheric plasma induced by additional electric fields. The present detailed analysis of the ionospheric response shows that the considered storm exhibited characteristic features inherent in the winter season but with some manifestations of reactions in equinox conditions.

Low latitude ionospheric TEC response to the solar flares during the peak of solar cycle 24

We describe the newly discovered response of differential vertical total electron content (DVTEC) in the low-equatorial latitude ionosphere to the solar flares during the peak of solar cycle 24 i.e., the year 2014. GPS TEC data for the existing investigation has been acquired from the international GNSS services network station in Bangalore, India (geomagnetic latitude 4.58οN). For the year 2014, we inspected five X-class solar flares, forty-nine M-class solar flares, and three hundred ninety-six C-class solar flares in the current study. To ascertain the relationship between ΔDVTEC and the ΔX-ray and ΔEUV fluxes, the two main ionizing radiation fluxes during flare events are used. Additionally, to identify electron density change in the ionosphere due to solar flares, we used two well-known methodologies, the baseline method and the mean method. The baseline method provides a better correlation between ΔDVTEC and ΔX-ray for X-, M-, and C-class flares. For M- and C- class flares, both methods do not show any meaningful correlation. The Solar disc location effect shows an enhanced (0.573) correlation between ΔDVTEC and ΔX-ray*cos (CMD) for the X-class flares, whereas for M- and C-class flares, no significant change was observed, which indicates feeble or no CMD effect for lower class solar flares while it is opposite for X-class flares. Here we suggest that the high solar activity, and location of the observation station to the EIA trough region might be the possible reason for the observed results.

Response of the Ionospheric Total Electron Content During St. Patrick's Day Geomagnetic Storm in 2015 over the Philippine-Taiwan Region

The ionospheric total electron content (TEC) response during the St. Patrick's Day geomagnetic storm of 2015 (max Kp = 8) is an interest of study in the field of space weather since it is the strongest geomagnetic storm of Solar Cycle 24. The cause of this storm is a coronal mass ejection (CME) from the Sun on 15 Mar 2015 recorded by SOHO/LASCO and arrived at ~ 05:00 UT 17 Mar 2015. The CME speed is estimated at ~ 668 km/s. At the initial phase of the storm, the disturbance storm time (Dst) index rose from 15 nanoteslas nT at 03:21–56 nT at 05:31 UT. The main phase of the storm caused the Dst index to drop to the minimum value of –223 nT from 05:31–22:33 UT on 17 Mar 2015. In the Philippine-Taiwan sector, TEC during the main phase was enhanced by ~ 25 TECU on 10 GNSS receiver stations and was depleted afterward with a maximum depletion of ~ 33 TECU on PBAS. The recovery phase showed a TEC depletion on all stations for the whole day. The minimum dTEC (differential TEC) is observed at PBAS (Basco, Batanes, Philippines), and the maximum percent depletion is observed at TWTF (Taoyuan, Taiwan) during the early hours of photoionization. Results from this study verified that GNSS TEC measurements along Asian sectors dropped significantly during St. Patrick's geomagnetic storm.

A semi-supervised total electron content anomaly detection method using LSTM-auto-encoder

Total electron content (TEC) is one of the important features used in studying of ionospheric properties. In seismo-ionospheric studies, variations/anomalies of the vertical total electron content over a seismic event are used to study the manifestation of lithospheric in the ionosphere. It is significant to design a robust TEC anomaly detection algorithm, which is capable of detecting non-spurious Seismo-induced TEC variations. In this paper, the long short term memory (LSTM) based auto-encoder network is presented for the detection of TEC anomalies recorded by GNSS receivers. The LSTM-auto-encoder is applied as a semi-supervised scheme to learn TEC responses in quiet solar and geomagnetic conditions. It then uses the learned TEC features to identify anomalies in a given TEC input data. The method is implemented to detect TEC anomalies recorded by three different GNSS-TEC receivers (ankr, ista, and tubi) in Türkiye. The model is also used to verify results from a recently published study on Mexico earthquake (magnitude 7.4). In each case, plausible results are obtained, and the relationship between detected anomalies with some lithosphere-atmosphere processes are discussed. The method highlights the significance/applications of AI in studying ionospheric variations.

Analysis of ionospheric TEC response to solar and geomagnetic activities at different solar activity stages

Studying the relationship of total electron content (TEC) to solar or geomagnetic activities at different solar activity stages can provide a reference for ionospheric modeling and prediction. On the basis of solar activity indices, geomagnetic activity parameters, and ionospheric TEC data at different solar activity stages, this study analyzes the overall variation relationships of solar and geomagnetic activities with ionospheric TEC, the characteristics of the quasi-27-day periodic oscillations of the three variables at different stages, and the delayed TEC response of solar activity by conducting correlation analysis, Butterworth band-pass filtering, Fourier transform, and time lag analysis. The following results are obtained. (1) TEC exhibits a significant linear relationship with solar activity at different solar activity stages. The correlation coefficients |R| are arranged as follows: |R|EUV > |R|F10.7 > |R|sunspot number. No significant linear relationship exists between TEC and geomagnetic activity parameters (|R| < 0.35). (2) TEC, solar activity indices, and geomagnetic activity parameters have a period of 10.5 years. The maximum amplitudes of the Fourier spectrum for TEC and solar activity indices are nearly 27 days and those of geomagnetic activity parameters are nearly 27 and 13.5 days. (3) The deviations of the quasi-27-day significant periodic oscillation of TEC and solar activity indices are consistent. (4) No evident relationship exists between the quasi-27-day periodic oscillation of TEC and geomagnetic activity parameters. (5) The delay time of TEC for the 10.7 cm solar radio flux and extreme ultraviolet is always consistent, whereas that for sunspot number varies at each stage.

Near‐Real‐Time Analysis of the Ionospheric Response to the 15 January 2022 Hunga Tonga‐Hunga Ha'apai Volcanic Eruption

AbstractWe present a near‐real‐time (NRT) scenario of analysis of ionospheric response to the 15 January 2022 Hunga Tonga‐Hunga Ha'apai eruption by using Global Navigation Satellite Systems data in the near‐field (in the vicinity of the volcano), and in the far‐field (Japan, North America, and South America). We introduce a new method to determine instantaneous velocities using an interferometric approach and using the time derivative of the total electron content (TEC). Moreover, for the first time, we propose a novel method that automatically estimates the apparent propagation velocity of ionospheric disturbances from NRT travel‐time diagrams. By using our new methods, we analyzed the dynamics of co‐volcanic ionospheric disturbances generated by the Hunga‐Tonga eruption, and we estimated the first propagation velocity in the near‐field to be ∼800–950 m/s, subsequently decreasing to ∼600 m/s. Based on these values, we conclude that in the near‐field, we detect ionospheric signatures of acoustic waves. In the far field, the apparent velocity of ionospheric disturbances was estimated to be between 277 and 365 m/s, which corresponds to the propagation of the Lamb wave. It is important to note that our new methods can successfully perform at low spatial resolution networks and with 30‐s cadence data. Also, they enable NRT spatio‐temporal analysis of ionospheric TEC response to smaller‐amplitude events.

Study of the equatorial and low-latitude total electron content response to plasma bubbles during solar cycle 24–25 over the Brazilian region using a Disturbance Ionosphere indeX

Abstract. This work uses the Disturbance Ionosphere indeX (DIX) to evaluate the ionospheric responses to equatorial plasma bubble (EPB) events from 2013 to 2020 over the Brazilian equatorial and low latitudes. We have compared the DIX variations during EPBs to ionosonde and All-Sky Imager data, aiming to evaluate the physical characteristics of these events. Our results show that the DIX was able to detect EPB-related TEC disturbances in terms of their intensity and occurrence times. Thus, the EPB-related DIX responses agreed with the ionosphere behavior before, during, and after the studied cases. Finally, we found that the magnitude of those disturbances followed most of the trends of solar activity, meaning that the EPB-related total electron content variations tend to be higher (lower) in high (low) solar activity.

The 15 January 2022 Hunga Tonga Eruption History as Inferred From Ionospheric Observations

AbstractOn 15 January 2022, the Hunga Tonga‐Hunga Ha’apai submarine volcano erupted violently and triggered a giant atmospheric shock wave and tsunami. The exact mechanism of this extraordinary eruptive event, its size and magnitude are not well understood yet. In this work, we analyze data from the nearest ground‐based receivers of Global Navigation Satellite System to explore the ionospheric total electron content (TEC) response to this event. We show that the ionospheric response consists of a giant TEC increase followed by a strong long‐lasting depletion. We observe that the explosive event of 15 January 2022 began at 04:05:54UT and consisted of at least five explosions. Based on the ionospheric TEC data, we estimate the energy released during the main major explosion to be between 9 and 37 Megatons in trinitrotoluene equivalent. This is the first detailed analysis of the eruption sequence scenario and the timeline from ionospheric TEC observations.

Response of the Ionospheric TEC to SSW and Associated Geomagnetic Storm Over the American Low Latitudinal Sector

AbstractDuring the sudden stratospheric warming (SSW) event in 2013, we investigated the American low latitude around 75°W. We used 12 Global Positioning System (GPS) receivers, a pair of magnetometers, and the NASA Thermosphere Ionosphere Mesosphere Energetics and Dynamics (TIMED) satellite airglow instrument to unveil the total electron content (TEC), inferred vertical drift, and the changes in the neutral composition, respectively. A major SSW characterized the 2013 SSW event with the main phase (7–27 January 2013) overlapped by a minor geomagnetic storm (17 January 2013). The late morning inferred downward‐directed E X B drift did not support the varying equatorial ionization anomaly (EIA) signature during the SSW onset (7 January 2013). The mid‐January (15–16 January 2013) witnessed enhancement in the varying inferred upward‐directed E X B drift at both hemispheres. On 17 January 2013, there were reductions in the varying inferred upward‐directed E X B drift at both hemispheres. Generally, the SSW effect on TEC around 15–16 January 2013 is more pronounced than the SSW onset. During the mid‐January (15–16 January 2013), the higher northern EIA crests are facilitated majorly by the SSW compared to the photo‐ionization that primarily enabled the southern crests. On 17 January 2013, the combined effect of photo‐ionization and SSW contribution was majorly responsible for the slight reduction in the northern crest. In the southern hemisphere, photo‐ionization played the lead role as the SSW, and the minor geomagnetic storm roles are secondary in enhancing the southern crest.

Interhemispheric Asymmetries in Ionospheric Electron Density Responses During Geomagnetic Storms: A Study Using Space‐Based and Ground‐Based GNSS and AMPERE Observations

AbstractWe utilize Total Electron Content (TEC) measurements and electron density (Ne) retrieval profiles from Global Navigation Satellite System (GNSS) receivers onboard multiple Low Earth Orbit (LEO) satellites to characterize large‐scale ionosphere‐thermosphere system responses during geomagnetic storms. We also analyze TEC measurements from GNSS receivers in a worldwide ground‐based network. Measurements from four storms during June and July 2012 (boreal summer months), December 2015 (austral summer month), and March 2015 (equinoctial month) are analyzed to study global ionospheric responses and the interhemispheric asymmetry of these responses. We find that the space‐based and ground‐based TECs and their responses are consistent in all four geomagnetic storms. The global 3D view from GNSS‐Radio Occultation (RO) Ne observations captures enhancements and the uplifting of Ne structures at high latitudes during the initial and main phases. Subsequently, Ne depletion occurs at high latitudes and starts progressing into midlatitude and low latitude as the storm reaches its recovery phase. A clear time lag is evident in the storm‐induced Ne perturbations at high latitudes between the summer and winter hemispheres. The interhemispheric asymmetry in TEC and Ne appears to be consistent with the magnitudes of the Active Magnetosphere and Planetary Electrodynamics Response Experiment (AMPERE) high latitude integrated field‐aligned currents (FACs), which are 3–4 MA higher in the summer hemisphere than in the winter hemisphere during these storms. The ionospheric TEC and Ne responses combined with the AMPERE‐observed FACs indicate that summer preconditioning in the ionosphere‐thermosphere system plays a key role in the interhemispheric asymmetric storm responses.

Persistent Eastward EEJ Enhancement During the Geomagnetic Storm Recovery Phases

AbstractRecent case studies disclosed that sometimes the dayside ionospheric equatorial electric field underwent persistent long‐duration (&gt;4 hr) eastward enhancements during geomagnetic storm recovery phases, rather than anticipated decrease or reversal. However, it is still unknown the occurrence and morphology of the persistent enhancement as well as its impact on the low‐latitude ionospheric plasma distribution. For the first time, we examined the issue with superposed epoch analyses of equatorial electrojet (EEJ) and GNSS total electron content (TEC) responses in the American (49°W) and Indian (78°E) sectors. The occurrence rates of the persistent enhancement during storm recovery phases are more than 30%, with higher occurrence rates in second to fourth days of the recovery phase, equinoxes, and American sector. The corresponding averaged EEJ and TEC responses were disclosed and compared with those storm samples without the persistent enhancement. On the first day of recovery phase, the dayside equatorial ionosphere was dominated by the prevalent westward electric field and causative inhibited equatorial fountain. In contrast, after the first day, the dayside morphologies showed different patterns: (a) for the persistent cases, the large eastward EEJ and TEC enhancement indicated the enhanced eastward electric field at forenoon and causative intensified equatorial fountain. (b) For the storms without the persistent cases, the general westward EEJ and equatorial TEC enhancements were dominated. Our study indicates that dayside eastward enhancement is an important disturbed pattern of equatorial ionospheric electric field during the late recovery phases, and the enhancement generally causes the intensification of equatorial ionization anomaly.

African and American Equatorial Ionization Anomaly (EIA) Responses to 2013 SSW Event

AbstractThis study investigates the responses of the African and American Equatorial Ionization Anomaly (EIA) regions to 2013 Sudden Stratospheric Warming (SSW) event. The Total Electron Content (TEC) data obtained from chains of Global Positioning System receivers within ±40° geomagnetic latitudes in the African and American sectors were used to construct the EIA structures for both longitudinal sectors. The responses of the EIA structures, constructed from the TEC, ΔTEC, and ionospheric irregularities data to the 2013 SSW event were investigated. During the SSW peak phase, EIA structures in both longitudinal sectors responded significantly, with the pole‐ward flow of plasma from the equator to higher crests' locations (strengthening of the EIA). Furthermore, a clear asymmetry in plasma distribution in the northern and southern crests of the EIA was observed. Generally, for the entire data span, TEC enhancements and ionospheric irregularities occurrences during SSW were more in the American sector than the African sector. The geomagnetic activity of 17 January 2013 caused negative TEC response in the African sector and positive TEC response in the American sector. Moderate storm‐induced TEC enhancements were generally lower than SSW‐induced TEC enhancements. Furthermore, solar flux‐induced TEC of 10 January 2013 was lesser than the SSW‐induced TEC of 15–16 January 2013.

English

The present paper reviews CODG TEC response to the moderate storm of July 26, 2003 and that of July 29 at Niamey station in West African equatorial station (Geo Lat 13&deg; 28&#39;45.3 &quot;N; Geo long: 02&deg; 10&#39;59.5&quot; E) in Niger. These two moderate geomagnetic storms are of solar wind origin. The study showed an increase in CODG TEC values during the initial phase of the storm and a decrease in CODG TEC values during the main and recovery phases of the storm. In general, in the equatorial region, positive and negative storms occurred (increase and decrease in TEC values). Here we highlight the prereversal enhancement (PRE), a particularity of the equatorial ionosphere. Our investigations show M, B and R profiles in addition to the well-known dome profile in the variation of the CODG TEC in Niamey. This work suggests a variation of the equatorial electrojet intensity during the storm recovery phase, the appearance of the counter-electrojet as well as a disturbance of the ExB drift. This study is local, in the equatorial region in the West African sector. Data are needed to corroborate some hypotheses such as changes in profiles associated with the presence or absence of electrojet and counter-electrojet currents. We advocate for a new IHY for Africa and for data exchange between researchers. Key words: Ionosphere, solar wind, disturbance, total electron content (TEC).

Quantifying Contributions of External Drivers to the Global Ionospheric State

AbstractTo advance our understanding of the global response of the ionospheric total electron content (TEC) to various external drivers, we apply entropy‐based information theory to quantify contributions of solar, interplanetary, and lower atmospheric drivers to the global ionospheric state described by the Global Ionospheric Map (GIM). By computing normalized transfer entropy on 18 years of GIM TEC, F10.7, solar wind, and lower atmospheric data, we obtain the predictive information transfer from the present drivers to the global ionospheric state at various future times. We find that the solar extreme ultraviolet (EUV) irradiance dominates the information transfer within 3 days into the future, while the lower atmospheric migrating tidal sources dominate beyond 3 days into the future. The maximum information transfer from individual drivers reveals the maximum contributions from the drivers to the global ionospheric state. Among all drivers considered, the lower atmospheric migrating tidal sources contribute most. The solar EUV irradiance and tropospheric deep convection are secondary drivers with comparable contributions at about half of the contribution from migrating tidal sources. The interplanetary driving contributes the least to the global ionospheric state during both geomagnetic storm and quiet time, even though the contribution from the interplanetary magnetic field can enhance four to six times during geomagnetic storm time than quiet time.