1. Introduction
The Global Electron Content (GEC) can be defined as the total number of free electrons in the ionosphere and plasmasphere up to the Global Positioning System (GPS) satellite altitude of 20,200 km [
1,
2]. The GEC can be computed by integrating the global ionospheric maps (GIMs) of VTEC on the overall ionosphere [
3]. The GIMs are routinely computed from worldwide Global Navigation Satellite System (GNSS) multifrequency measurements by different analysis centers in the context of the International GNSS Service [
4,
5,
6].
In order to be compared with the GEC we can consider the K-index [
7], which quantifies disturbances in the horizontal component of Earth’s magnetic field with an integer in the range 0–9 with 1 being calm and 5 or more indicating a geomagnetic storm. The Kp-index is calculated by combining the data from multiple magnetic observatories around the world to determine the global planetary index.
Motivated by previous studies of specific geomagnetic storms where the GEC seems to evolve correlated with geomagnetic indices [
8], GEC and Kp are compared in this work in both the temporal and frequency domains. The study is focused on the GEC and Kp spectrum evolution on the 10 days time intervals centered at the 34 geomagnetic events with Kp ≥ 7.5 from 2000 to 2020.
Our first motivation for the present work was the apparent correlation found between some space weather indices and GEC time evolutions around the St. Patrick 2015 geomagnetic storm [
8], as can be seen in
Figure 1.
This apparent relationship of the GEC (hereinafter
G) vs. SW indices has been studied in previous works under different approaches. Indeed, Ref. [
9] summarized the GEC storm time modeling, from the JPL GIMs from 1999 to 2014, associated with variability of smoothed and normalized Auroral Electrojet index (positive correlation). The study considers the DGEC values taking the hourly ratio of instant GEC to 7 preceding days median.
Moreover, the authors of [
10] studied 90 storm events from 1999 to 2011 from IGS GIMs in a different way. They conclude that the ratio
is closely correlated with geomagnetic Kp index (positively) and time weighted Dst index (negatively) and also affected by F10.7, where
is the 6 quiet days’ smoothed average GEC (Kp < 3.0).
More recently the authors of [
11], with an approach similar to that used by [
12,
13], studied the geomagnetic storm of 7–9 September 2017 with the UPC-IonSAT UQRG GIMs:
shows a clear correlation with SW indices like Kp and SYM-H (where the quiet time GEC,
is obtained by using the three quiet days before the storm on which the Kp index is below 4).
In order to dispose of a most general view of GEC and Kp variability, not restricted to the given reference time windows for detrending (like the three ones of the mentioned recent works), we have performed a systematic comparison of the overall spectrum time evolution of GEC and Kp.
2. Methodology
The GEC time series for this study has been extracted from the GIMs of VTEC estimated by UPC-IonSAT with TOMION-v1 software from global GPS carrier phase measurements (i.e., not affected by pseudorange noise, multipath or differential code biases, see for instance [
14]) since the end of 1996. An spectral analysis of the GEC time series is presented and, afterwards, the spectrograms have been computed and studied around major geomagnetic events.
2.1. How Is the GEC Computed
In order to compute the so-called UQRG GIM, first, a dual-layer voxel-based tomographic model is solved from dual-frequency carrier GPS data only by means of a forward Kalman filter; and, secondly, the kriging interpolation is applied. It has a resolution of 15 min in time, 5° in longitude and 2.5° in latitude, and is one of the best behaving GIMs in the International GNSS Service, IGS ([
4,
5]).
We have considered the GEC, computed as the surface integral of the GIM VTEC, from 2000 to 2020 (almost two solar cycles) but downdated to the Kp cadence of 3 h. This provided us a total of 34 periods of 10-days to study, centered at the geomagnetic events with Kp ≥ 7.5 in such a time interval, only taking one single event when they are consecutive (separated by less than 3 days). Then the GEC and Kp, but specially their spectra time evolution, by applying a continuous wavelet transformation, are compared.
2.2. How Are the Spectrograms Computed
Spectrograms are obtained through a continuous wavelet transformation, by using the corresponding package under the R statistical programming language. The basic wavelet function is the Morlet form:
To transform a time series
, indexed on time
t, we compute the scalar product of
and the shifted (by a time
) and scaled version of
, with scaling factor
s. The output of the transformation is:
And the final results and corresponding figures shown below correspond to the wavelet power spectrum
.
Since the time series are finite samples, the projection over the wavelet function is less significant for close to the boundary when the periods are large. This is represented in the figures below by a lighter color region.
3. Results
After performing the GEC spectral analysis, the comparisons between Kp and GEC are presented in spectral and time domains, and in both ways: qualitatively, by visual inspection, and quantitatively, by means of the two corresponding Pearson correlation coefficients.
3.1. Spectrum of the GEC Time Evolution
The spectral analysis of the GEC time evolution (
Figure 2), computed every 15 min since the end of 1996 to September 2023, is summarized and compared with the previous results [
15] in this section, following the decreasing order of spectral terms significance. We can see during more than two solar cycles, the annual and semiannual periods, with GEC maxima at the equinoxes, and minima at the solstices (which appear systematically low at the June one). And we can appreciate as well (in the resolution limit of the plot) that the higher frequency oscillations are associated with the solar rotation.
These four periods appear as the first, fourth, third and sixteenth most intense periods in the Fast Fourier transform (hereinafter FFT) of the GEC time series (
Table 1), as it can be seen as well in the semilog and log-log representation of their module vs. period (
Figure 3 and
Figure 4, respectively), confirming the appearance under visual inspection. These GEC periods given by the FFT, after being applied to the almost one million GEC values derived from almost 27 years of 15 min GIMs, were also detected and reported in Figure 22 of [
15], computed since June 1998 to end of 2007 from our UPCG GIMs at 2 h, i.e., a time series of around 40 thousands GIMs. In other words, in this study we analyze a set of GEC values more than 20 times larger than in [
15] and with and extension almost three times larger (around 2.5 solar cycles vs. almost one), reaching to higher periods in the spectrum that were not available in the former study (see NA in
Table 1), but confirming the most part of periods in the lower period interval.
Table 1 shows the (ordered) local maxima for the absolute value of the FFT of GEC, the periods where they occur and a estimation for the error on such periods. Since frequency is an integer index, the error on the period increases with the period value.Indeed, in
Table 1 the complete set of the twenty most prominent periods in the temporal GEC spectrum can be found, and it shows that nine of them are in agreement with the above commented study (see the fifth column in the above-mentioned table). Moreover, it can be seen that other peak periods which are clearly defined in the GEC FFT: the one of
days is very well in agreement with the solar synodic rotation period (27.2753 days [
16]), which affects the sunspot groups to the solar sunspots rotation period, the daily, half-daily and 6 h, the four of them are also observed in [
15].
3.2. Qualitative Comparison of GEC and Kp
Examples of different qualitative levels of GEC and Kp spectra similarity can be seen in
Figure 5, in terms of extremely high (day 2017-09-08, top plot), high (day 2013-10-02, second plot), mid (day 2012-03-09, third plot), low (day 2005-05-15, fourth plot) and extremely low similarity (2002-05-24, fifth plot).
In summary, more than 50% of the GEC and Kp spectrograms during the geomagnetic events appear high or extremely high correlated by visual inspection, and only 21% of them appear low or extremely low correlated (see
Table 2). In order to understand the context of low or extremely low spectral similarity between Kp and GEC, we represent the time evolution of Kp and GEC qualitative spectral similarity, besides F10.7 index, during the solar cycle in
Figure 6. It can be seen that geomagnetic storms with Kp
presenting extremely low similarity between GEC and Kp spectra (four cases of 34) happen around the two extremely high F10.7 peaks (>200) of previous solar cycle.
This might be related with the correlation between the GEC and the F10.7 index, following [
1,
12]), which might affect the comparison with space weather indices, as it was pointed out by [
11]. But before concluding about this potential explanations of this result, we have also performed, in the next subsection, a quantitative assessment of the degree of similarity between both, GEC and Kp distributions.
3.3. Quantitative Comparison of GEC and Kp
It can be seen in
Figure 7 that the quantitative similarity, correlation of both GEC and Kp spectra, confirms the qualitative one previously shown in
Figure 6: indeed, a correlation decrease from predominant positive values up to +0.6 to values within the interval
, can be seen around the 2000–2002 solar cycle peak (coinciding with extremely low similarity in the previous plot). And higher correlation values, from 0.5 to almost 0.8, happen otherwise.
It is interesting to see in the same
Figure 7 the corresponding time evolution of the correlation of the GEC and Kp time series, which is systematically lower (as it is confirmed in the one-to-one spectrum vs. time GEC-Kp correlations represented in
Figure 8), and in some periods beyond the 2000–2002 solar cycle peak, much lower. This result might be related with the time window used in the spectral analysis, less sensitive to latencies between GEC and Kp reactions to space weather conditions.
4. Conclusions
The GEC time series presents clear spectral signatures in periods which are confirmed, with regard to previous studies with a much reduced dataset, and extended in this research, based on +900,000 GIMs with +5000 VTEC pixel global values of each GIM (resolution of 15 min × 5° × 2.5° in time, longitude and latitude since the end of 1996 to September 2023). Motivated by previous works, the GEC and Kp geomagnetic index spectrograms have been systematically compared during time windows of 10 days around the 34 geomagnetic events, defined with Kp , from 2000 to 2020. After visual inspection, the majority of them (18 geomagnetic events) present a high or extremely high similarity, confirming the potentiality of using directly the GEC as SW index. However, there is a minority of GEC and Kp comparisons showing low or extremely low spectra similarity (7 of 34), which happened mostly around the strong solar flux epochs corresponding to the bimodal solar cycle peaks in 2000 and 2002. This result, quantitatively confirmed, is consistent with the potential contribution of high ionization epochs to the GEC spectra differently than the geomagnetic activity, pointed out by other authors in previous works.
In summary, we show in this work the higher correlation of GEC and Kp spectra, vs. the correlation of the GEC and Kp time series, during high geomagnetic activity periods. It can indicate the corresponding correlation of the extraordinary electron currents at high latitudes with the variability of electron content at large scale. Regardless of the potential common physical origin, the results strongly suggest the use of the GEC spectra to characterize space weather conditions.
Author Contributions
Conceptualization, M.H.-P. and J.M.A.-F.; methodology, J.M.A.-F. and M.H.-P.; software, J.M.A.-F.; validation, J.M.A.-F. and M.H.-P.; formal analysis, M.H.-P. and J.M.A.-F.; investigation, J.M.A.-F. and M.H.-P.; resources, M.H.-P.; data curation, J.M.A.-F. and M.H.-P.; writing—original draft preparation, M.H.-P.; writing—review and editing, J.M.A.-F., M.H.-P., H.L., D.R.-D., R.O.-P., A.G.-R., V.G., G.O.-P., E.M.-M., H.Y. and Q.L.; visualization, J.M.A.-F. and M.H.-P.; supervision, M.H.-P.; project administration, M.H.-P.; funding acquisition, M.H.-P. All authors have read and agreed to the published version of the manuscript.
Funding
This research was partially funded by the European Union (PITHIA-NRF project, grant number H2020-INFRAIA-2018-2020101007599).
Data Availability Statement
Acknowledgments
The last part of the work has been performed coinciding with the starting of the DISPEC project (HORIZON-CL4-2023-SPACE-01-71). The authors would like to extend their sincere appreciation to the European Union.
Conflicts of Interest
The authors declare no conflict of interest.
References
- Afraimovich, E.; Astafyeva, E.; Oinats, A.; Yasukevich, Y.V.; Zhivetiev, I. Global electron content: A new conception to track solar activity. Ann. Geophys. 2008, 26, 335–344. [Google Scholar] [CrossRef]
- Gulyaeva, T.; Veselovsky, I. Two-phase storm profile of global electron content in the ionosphere and plasmasphere of the Earth. J. Geophys. Res. Space Phys. 2012, 117, A09324. [Google Scholar] [CrossRef]
- Hernández-Pajares, M. GNSS Ionosphere. Encyclopedia of Geodesy; Springer: Berlin/Heidelberg, Germany, 2022; pp. 1–7. [Google Scholar]
- Hernández-Pajares, M.; Roma-Dollase, D.; Krankowski, A.; García-Rigo, A.; Orús-Pérez, R. Methodology and consistency of slant and vertical assessments for ionospheric electron content models. J. Geod. 2017, 91, 1405–1414. [Google Scholar] [CrossRef]
- Roma-Dollase, D.; Hernández-Pajares, M.; Krankowski, A.; Kotulak, K.; Ghoddousi-Fard, R.; Yuan, Y.; Li, Z.; Zhang, H.; Shi, C.; Wang, C.; et al. Consistency of seven different GNSS global ionospheric map** techniques during one solar cycle. J. Geod. 2018, 92, 691–706. [Google Scholar] [CrossRef]
- Hernández-Pajares, M.; Juan, J.M.; Sanz, J.; Aragón-Àngel, À.; García-Rigo, A.; Salazar, D.; Escudero, M. The ionosphere: Effects, GPS modeling and the benefits for space geodetic techniques. J. Geod. 2011, 85, 887–907. [Google Scholar] [CrossRef]
- Bartels, J.; Heck, N.; Johnston, H. The three-hour-range index measuring geomagnetic activity. Terr. Magn. Atmos. Electr. 1939, 44, 411–454. [Google Scholar] [CrossRef]
- Béniguel, Y.; Cherniak, I.; Garcia-Rigo, A.; Hamel, P.; Hernández-Pajares, M.; Kameni, R.; Kashcheyev, A.; Krankowski, A.; Monnerat, M.; Nava, B.; et al. MONITOR Ionospheric Network: Two case studies on scintillation and electron content variability. Ann. Geophys. 2017, 35, 377–391. [Google Scholar] [CrossRef]
- Gulyaeva, T.; Yenen, S.D. Statistical Investigation of Ionospheric Storm and Substorm Events Using GEC and AE-Index. In Proceedings of the 14th Ionospheric Effects Symposium (IES2015), Old Town Alexandria, VA, USA, 12–14 May 2015. [Google Scholar]
- Li, S.; Galas, R.; Ewert, D.; Peng, J. An empirical model for the ionospheric global electron content storm-time response. Acta Geophys. 2016, 64, 253–269. [Google Scholar] [CrossRef]
- Imtiaz, N.; Younas, W.; Khan, M. Response of the low-to mid-latitude ionosphere to the geomagnetic storm of September 2017. Ann. Geophys. 2020, 38, 359–372. [Google Scholar] [CrossRef]
- Nava, B.; Rodríguez-Zuluaga, J.; Alazo-Cuartas, K.; Kashcheyev, A.; Migoya-Orué, Y.; Radicella, S.; Amory-Mazaudier, C.; Fleury, R. Middle-and low-latitude ionosphere response to 2015 St. Patrick’s Day geomagnetic storm. J. Geophys. Res. Space Phys. 2016, 121, 3421–3438. [Google Scholar] [CrossRef]
- Kashcheyev, A.; Nava, B.; Alazo-Cuartas, K.; Radicella, S. Multivariable comprehensive analysis of two great geomagnetic storms of 2015. J. Geophys. Res. Space Phys. 2018, 123, 5000–5018. [Google Scholar] [CrossRef]
- Hernández-Pajares, M.; Lyu, H.; Garcia-Fernandez, M.; Orus-Perez, R. A new way of improving global ionospheric maps by ionospheric tomography: Consistent combination of multi-GNSS and multi-space geodetic dual-frequency measurements gathered from vessel-, LEO-and ground-based receivers. J. Geod. 2020, 94, 73. [Google Scholar] [CrossRef]
- Hernández-Pajares, M.; Juan, J.; Sanz, J.; Orus, R.; García-Rigo, A.; Feltens, J.; Komjathy, A.; Schaer, S.; Krankowski, A. The IGS VTEC maps: A reliable source of ionospheric information since 1998. J. Geod. 2009, 83, 263–275. [Google Scholar] [CrossRef]
- Feng, X.; Yang, L.; **ang, C.; Jiang, C.; Ma, X.; Wu, S.; Zhong, D.; Zhou, Y. Validation of the 3D AMR SIP–CESE solar wind model for four Carrington rotations. Sol. Phys. 2012, 279, 207–229. [Google Scholar] [CrossRef]
Figure 1.
(
Top plots): Space weather indices,
,
,
,
and
vs. time, in days of March 2015, around the St. Patrick 2015 geomagnetic storm (see Sudden Storm Commencement indicated by the vertical red line). (
Bottom plot): GEC vs. time in days of year 2015, and with the same time scale than top plots (reproduced from [
8]).
Figure 1.
(
Top plots): Space weather indices,
,
,
,
and
vs. time, in days of March 2015, around the St. Patrick 2015 geomagnetic storm (see Sudden Storm Commencement indicated by the vertical red line). (
Bottom plot): GEC vs. time in days of year 2015, and with the same time scale than top plots (reproduced from [
8]).
Figure 2.
GEC, expressed in GECU, vs. time, in years, since the end of 1996 to September 2023. It has been computed from UQRG GIMs, with a time resolution of 15 min.
Figure 2.
GEC, expressed in GECU, vs. time, in years, since the end of 1996 to September 2023. It has been computed from UQRG GIMs, with a time resolution of 15 min.
Figure 3.
Absolute value of the GEC FFT, expressed in GECU, vs. logarithm of period in days, since the end of 1996 to September 2023. It has been computed from UQRG GIMs, with a time resolution of 15 min.
Figure 3.
Absolute value of the GEC FFT, expressed in GECU, vs. logarithm of period in days, since the end of 1996 to September 2023. It has been computed from UQRG GIMs, with a time resolution of 15 min.
Figure 4.
Logarithm of the absolute value of the GEC FFT, expressed in GECU, vs. logarithm of period in days, since the end of 1996 to September 2023. It has been computed from UQRG GIMs, with a time resolution of 15 min.
Figure 4.
Logarithm of the absolute value of the GEC FFT, expressed in GECU, vs. logarithm of period in days, since the end of 1996 to September 2023. It has been computed from UQRG GIMs, with a time resolution of 15 min.
Figure 5.
Representative examples of GEC and Kp time evolution (top and bottom, respectively of left-hand plots), GEC spectrograms (central plots) and Kp spectrograms (right plots), vs. time (horizontal axis, in days) and period (vertical axis, in multiples of three hours), with extremely high similarity (day 8 September 2017, row a), high similarity (day 2 October 2013, row b), mid similarity (day 9 March 2012, row c), low similarity (day 15 May 2005, row d) and extremely low similarity (24 May 2002, row e).
Figure 5.
Representative examples of GEC and Kp time evolution (top and bottom, respectively of left-hand plots), GEC spectrograms (central plots) and Kp spectrograms (right plots), vs. time (horizontal axis, in days) and period (vertical axis, in multiples of three hours), with extremely high similarity (day 8 September 2017, row a), high similarity (day 2 October 2013, row b), mid similarity (day 9 March 2012, row c), low similarity (day 15 May 2005, row d) and extremely low similarity (24 May 2002, row e).
Figure 6.
Evolution during more than one solar cycle of the level of qualitative similarity between GEC and Kp spectrograms (violet points), compared with the evolution of the F10.7 solar radio flux index (green line).
Figure 6.
Evolution during more than one solar cycle of the level of qualitative similarity between GEC and Kp spectrograms (violet points), compared with the evolution of the F10.7 solar radio flux index (green line).
Figure 7.
Evolution during more than one solar cycle of the correlation between GEC and Kp spectrograms (violet points) and the correlation between the GEC and Kp time series (light blue points). The evolution of the F10.7 solar radio flux index is also represented (green line).
Figure 7.
Evolution during more than one solar cycle of the correlation between GEC and Kp spectrograms (violet points) and the correlation between the GEC and Kp time series (light blue points). The evolution of the F10.7 solar radio flux index is also represented (green line).
Figure 8.
Points representing the correlation between GEC and Kp spectra vs. the corresponding correlation of GEC and Kp time evolutions during the geomagnetic events with Kp ≤ 7.5 (the straight line represent equal values for both).
Figure 8.
Points representing the correlation between GEC and Kp spectra vs. the corresponding correlation of GEC and Kp time evolutions during the geomagnetic events with Kp ≤ 7.5 (the straight line represent equal values for both).
Table 1.
The first twenty period values and error estimation in days, corresponding to the twenty highest absolute values of the GEC FFT modulus, are included in decreasing order. They have been computed from the complete GEC time series obtained from UQRG GIMs computed from UPC by the second author since the end of 1996 to September 2023. The best defined maxima from visual inspection of the corresponding plots,
Figure 3 and
Figure 4, are emphasized in bold. The periods also given in years and approximate solar rotation period (27.3 days) time units. The fifth column indicates whether the estimated period was already detected in [
15] or not, from just one solar cycle of GIMs, including the periods too large to be detected in such a work (identified as Non Available, NA).
Table 1.
The first twenty period values and error estimation in days, corresponding to the twenty highest absolute values of the GEC FFT modulus, are included in decreasing order. They have been computed from the complete GEC time series obtained from UQRG GIMs computed from UPC by the second author since the end of 1996 to September 2023. The best defined maxima from visual inspection of the corresponding plots,
Figure 3 and
Figure 4, are emphasized in bold. The periods also given in years and approximate solar rotation period (27.3 days) time units. The fifth column indicates whether the estimated period was already detected in [
15] or not, from just one solar cycle of GIMs, including the periods too large to be detected in such a work (identified as Non Available, NA).
Period for 20 Highest GEC Local Maxima | | Is This Period Detected in [15]? |
---|
/Days | /Years | /27.3 Days |
---|
| | | | NA |
| | | | NA |
| | | | Yes ) days |
| | | | Yes ) days |
| | | | NA |
| | | | NA |
| | | | Yes () days |
| | | | No |
| | | | No |
| | | | No |
| | | | Yes () days |
| | | | Yes () days |
| | | | No |
| | | | No |
| | | | No |
| | | | Yes () days |
| | | | Yes () days |
| | | | Yes () days |
| | | | Yes () days |
| | | | No |
Table 2.
Summary of GEC and Kp spectrogram similarity determined by visual inspection.
Table 2.
Summary of GEC and Kp spectrogram similarity determined by visual inspection.
GEC and Kp Spectograms Similarity | Number of Cases | Percentage of Cases |
---|
Extremely high | 10 | 29% |
High | 8 | 24% |
Mid | 9 | 26% |
Low | 3 | 9% |
Extremely low | 4 | 12% |
| Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).