Computerized ionospheric tomography with the IRI model

Computerized ionospheric tomography (CIT) is a method to estimate ionospheric electron density distribution by using the global positioning system (GPS) signals recorded by the GPS receivers. Ionospheric electron density is a function of latitude, longitude, height and time. A general approach in CIT is to represent the ionosphere as a linear combination of basis functions. In this study, the model of the ionosphere is obtained from the IRI in latitude and height only. The goal is to determine the best representing basis function from the set of Squeezed Legendre polynomials, truncated Legendre polynomials, Haar Wavelets and singular value decomposition (SVD). The reconstruction algorithms used in this study can be listed as total least squares (TLS), regularized least squares, algebraic reconstruction technique (ART) and a hybrid algorithm where the reconstruction from the TLS algorithm is used as the initial estimate for the ART. The error performance of the reconstruction algorithms are compared with respect to the electron density generated by the IRI-2001 model. In the investigated scenario, the measurements are obtained from the IRI-2001 as the line integral of the electron density profiles, imitating the total electron content estimated from GPS measurements. It has been observed that the minimum error between the reconstructed and model ionospheres depends on both the reconstruction algorithm and the basis functions where the best results have been obtained for the basis functions from the model itself through SVD.     

Empirical model of the TEC response to the geomagnetic activity over the North American region

The paper presents an empirical model of the total electron content (TEC) response to the geomagnetic activity described by the K_p-index. The model is built on the basis of TEC measurements covering the region of North America (50°W–150°W, 10°N–60°N) for the period of time between October 2004 and December 2009. By using a 2D (latitude-time) cross-correlation analysis it is found that the ionospheric response to the geomagnetic activity over the considered geographic region and at low solar activity revealed both positive and negative phases of response. The both phases of the ionospheric response have different duration and time delay with respect to the geomagnetic storm. It was found that these two parameters of the ionospheric response depend on the season and geographical latitude. The presence of two phases, positive and negative, of the ionospheric response imposed the implementation of two different time delay constants in order to properly describe the two different delayed reactions. The seasonal dependence of the TEC response to geomagnetic storms is characterized by predominantly positive response in winter with a short (usually ∼5–6 h) time delay as well as mainly negative response in summer with a long (larger than 15 h) time delay. While the TEC response in March and October is more close to the winter one the response in April and September is similar to the summer one.     

Hierarchical Bayesian modeling of ionospheric TEC disturbances as non-stationary processes

We model regular and irregular variation of ionospheric total electron content as stationary and non-stationary processes, respectively. We apply the method developed to SCINDA GPS data set observed at Bahir Dar, Ethiopia $\left(11.6°\text{N,}37.4°\text{E}\right)$. We use hierarchical Bayesian inversion with Gaussian Markov random process priors, and we model the prior parameters in the hyperprior. We use Matérn priors via stochastic partial differential equations, and use scaled $\text{Inv}-{\chi }^2$ hyperpriors for the hyperparameters. For drawing posterior estimates, we use Markov Chain Monte Carlo methods: Gibbs sampling and Metropolis-within-Gibbs for parameter and hyperparameter estimations, respectively. This allows us to quantify model parameter estimation uncertainties as well. We demonstrate the applicability of the method proposed using a synthetic test case. Finally, we apply the method to real GPS data set, which we decompose to regular and irregular variation components. The result shows that the approach can be used as an accurate ionospheric disturbance characterization technique that quantifies the total electron content variability with corresponding error uncertainties.     

Ionospheric F-region response to geomagnetic disturbances

The F2-region reaction to geomagnetic storms usually called as an ionospheric storm is a rather complicated event. It consists of so called positive and negative phases, which have very complicated spatial and temporal behavior. The main morphological features of ionospheric storms and the main processes governing their behavior were understood at the end of the 1900s and described in a series of review papers. During the recent decade there were many publications dedicated to the problem of ionospheric storms. In this paper a concept of ionospheric storm morphology and physics formulated at the end of the 1990s is briefly summarized and the most interesting results obtained in the 2000s are described. It is shown that the main features of the studies of the previous decade were: the use of GPS TEC data for analyzing the ionospheric storm morphology, attraction of sophisticated theoretical models for studying the processes governing ionospheric behavior in disturbed conditions, and accent to analysis of ionospheric behavior during prominent events (very strong and great geomagnetic storms). Also a special attention was paid to the pre-storm enhancements in foF2 and TEC.     

Time delay of ionospheric TEC storms to geomagnetic storms and pre-storm disturbance events in East Asia

Using the TEC data at Beijing (39.61°N, 115.89°E)/Yakutsk (62.03°N, 129.68°E) stations of East Asia regions and relevant geomagnetic data from 2010 to 2017, we have studied the time delay of ionospheric storms to geomagnetic storms and compare it with our previous results of Taoyuan (25.02°N, 121.21°E) station (Zhang et al., 2020). The data shows a well-known local time dependence of the time delay, and seasonal dependences are different at these stations. In addition, there is no correlation between the time delay and the magnetic storm intensity /solar activity, except the time delay of negative storms has weakly negative dependence on the solar activity. Comparing with the results of Taoyuan station which is located at EIA region in East Asia, we find that the time delay increases nonlinearly as the latitude decreases due to different ionospheric backgrounds at these places. Moreover, the pre-storm disturbance events are found to have similar statistical characteristics as the pre-storm enhancement in Europe middle latitudes (Bure?ová and La?tovi?ka, 2007). By subtracting the common features of the pre-storm disturbance events, we preliminarily infer that auroral activity might be main driver of the pre-storm disturbance events.     

Modeling of ionospheric parameter variations in East Asia during the moderate geomagnetic disturbances

The results of modeling of ionospheric disturbances observed in the East Asian region during moderate storms are presented. The numerical model for ionosphere–plasmasphere coupling developed at the ISTP SB RAS is used to interpret the data of observations at ionospheric stations located in the longitudinal sector of 90–130°E at latitudes from auroral zone to equator. There is obtained a reasonable agreement between measurements and modeling results for winter and equinox. In the summer ionosphere, at the background of high ionization by the solar EUV radiation in the quiet geomagnetic period the meridional thermospheric wind strongly impacts the electron concentration in the middle and auroral ionosphere. The consistent calculations of the thermospheric wind permit to obtain the model results which are closer to summer observations. The actual information about the space-time variations of thermosphere and magnetosphere parameters should be taken into account during storms.     

Midlatitude ionospheric F2-layer response to eruptive solar events-caused geomagnetic disturbances over Hungary during the maximum of the solar cycle 24: A case study

In our study we analyze and compare the response and behavior of the ionospheric F2 and of the sporadic E-layer during three strong (i.e., Dst?<??100nT) individual geomagnetic storms from years 2012, 2013 and 2015, winter time period. The data was provided by the state-of the art digital ionosonde of the Széchenyi István Geophysical Observatory located at midlatitude, Nagycenk, Hungary (IAGA code: NCK, geomagnetic latitude: 46.17° geomagnetic longitude: 98.85°). The local time of the sudden commencement (SC) was used to characterize the type of the ionospheric storm (after Mendillo and Narvaez, 2010). This way two regular positive phase (RPP) ionospheric storms and one no-positive phase (NPP) storm have been analyzed. In all three cases a significant increase in electron density of the foF2 layer can be observed at dawn/early morning (around 6:00 UT, 07:00 LT). Also we can observe the fade-out of the ionospheric layers at night during the geomagnetically disturbed time periods. Our results suggest that the fade-out effect is not connected to the occurrence of the sporadic E-layers.     

Analysis of global TEC prediction performance of IRI-PLAS model

The International Reference Ionosphere (IRI) empirical model provides valuable data for many fields including space and navigation applications. Since the IRI model gives the ionospheric parameters in the altitude range from 50?km to 2000?km, researchers focused on the IRI-PLAS model which is the plasmasphere extension of the IRI model. In this study, Total Electron Content (TEC) prediction performance of the IRI-PLAS model was examined at a global scale using the location of globally distributed 9 IGS stations. Besides the long term (01.01.2015–31.12.2015) behavior of the model, TEC predictions during the equinox and solstice days of 2014–2017 were also tested. IRI-PLAS-TEC values were examined in comparison with GPS-TEC data. Hourly interval of yearly profile exhibits that when the geomagnetic and solar active days are ignored, differences between IRI-PLAS-TEC and GPS-TEC are rather small (～2–3 TECU) at stations in the northern hemisphere, generally ～4–5 TECU level at the southern hemisphere stations and reaching above 10 TECU for few hours. While the IRI-PLAS-TEC generally overestimates the GPS-TEC at southern hemisphere stations during quiet days, the model-derived TEC underestimates GPS-TEC during solar active days. IRI-PLAS-TEC and GPS-TEC values exhibit similar trend for the equinoxes 21 March and 23 September which refer equivalent conditions.     

Study of TEC changes during geomagnetic storms occurred near the crest of the equatorial ionospheric ionization anomaly in the Indian sector

GPS satellites data obtained at Bhopal (23.16° N, 77.36° E, geomagnetic latitude 14.23° N) India were analyzed to study the TEC changes during several geomagnetic storms (−300 nT < Dst < −50 nT) occurred in 2005–2007. We had segregated the storms according to the Dst value, i.e. moderate storms (−100 nT < Dst ? −50 nT), strong storms (−150 nT < Dst < −100 nT), and severe storms (Dst less than −150 nT). Total of 21 geomagnetic storms (10 moderate, 9 strong, 2 severe) are considered for the present study. Deviation in vertical total electron content (VTEC) during the main phase of the storm was found to be associated with the prompt penetration of electric field originated due to the under-shielding and over-shielding conditions for almost all geomagnetic storms discussed in this paper. For most of the storms VTEC shows the positive percentage deviation during the main phase while it shows positive as well as the negative deviation during the recovery phase of the storms. The −80% deviation in VTEC was found for geomagnetic storm occurred on July 17, 2005 and the negative trend continued for recovery phase of the storm. This was mainly due to the thermospheric composition changes by Joule heating effect at auroral latitudes that generate electric field disturbance at low latitudes. Traveling ionospheric disturbances (TIDs) were responsible for the formation of wave like nature in VTEC for the storms occurred on May 15, 2005, whereas it was not observed for storm occurred on August 24, 2005.     

Climatology of ionospheric scintillations and TEC trend over the Ugandan region

This study presents results on the investigation of the diurnal, monthly and seasonal variability of Total Electron Content (TEC), phase (_σΦ${\sigma }_{\mathrm{\Phi }}$) and amplitude (S4) scintillation indices over Ugandan (Low latitude) region. Scintillation Network Decision Aid (SCINDA) data was obtained from Makerere (0.34°N, 32.57°E) station, Uganda for two years (2011 and 2012). Data from two dual frequency GPS receivers at Mbarara (0.60°S, 30.74°E) and Entebbe (0.04°N, 32.44°E) was used to study TEC climatology during the same period of scintillation study. The results show that peak TEC values were recorded during the months of October–November, and the lowest values during the months of July–August. The diurnal peak of TEC occurs between 10:00 and 14:00 UT hours. Seasonally, the ascending and descending phases of TEC were observed during the equinoxes (March and September) and solstice (June and December), respectively. The scintillations observed during the study were classified as weak (0.1≤$\le$S4,_σΦ≤${\sigma }_{\mathrm{\Phi }}\le$0.3) and strong (0.3<$<$S4,_σΦ≤${\sigma }_{\mathrm{\Phi }}\le$1.0). The diurnal scintillation pattern showed peaks between 17:00 and 22:00 UT hour, while the seasonal pattern follows the TEC pattern mentioned above. Amplitude scintillation was more dominant than phase scintillation during the two years of the study. Scintillation peaks occur during the months of March–April and September–October, while the least scintillations occur during the months of June–July. Therefore, the contribution of this study is filling the gap in the current documentation of amplitude scintillation without phase scintillation over the Ugandan region. The scintillations observed have been attributed to wave-like structures which have periods of about 2–3 h, in the range of that of large scale travelling ionospheric disturbances (LSTIDs).     

F2 layer response to geomagnetic disturbances in Eastern Asia under the low solar activity

In this paper, the peculiarities of ionospheric response to geomagnetic disturbances observed at the decay and minimum of solar activity (SA) in the period 2004–2007 are investigated with respect to different geomagnetic conditions. Data from ionospheric stations and results of total electron content (TEC) measurements made at the network of GPS ground-based receivers located within the latitude–longitude sector (20–70°N, 90–160°Е) are used in this study. Three groups of anomalous ionospheric response to geomagnetic disturbances have been observed during low solar activity. At daytime, the large-scale traveling ionospheric disturbances (LSTIDs) could generally be related to the main phase of magnetic storm. Quasi-two-days wavelike disturbances (WLDs) have been also observed in the main phase independent of the geomagnetic storm intensity. Sharp electron density oscillations of short duration (OSD) occurred in the response to the onset of both main and recovery phases of the magnetic storm in the daytime at middle latitudes. A numerical model for ionosphere–plasmasphere coupling was used to interpret the occurrence of LS TIDs. Results showed that the LSTIDs might be associated with the unexpected lifting of F2 layer to the region with the lower recombination rate by reinforced meridional winds that produces the increase of the electron density in the F2 layer maximum.     

Responses of ionospheric foF2 to geomagnetic activities in Hainan

Responses of low-latitude ionospheric critical frequency of F2 layer to geomagnetic activities in different seasons and under different levels of solar activity are investigated by analyzing the ionospheric foF2 data from DPS-4 Digisonde in Hainan Observatory during 2002–2005. The results are as follows: (1) the response of foF2 to geomagnetic activity in Hainan shows obvious diurnal variation except for the summer in low solar activity period. Generally, geomagnetic activity will cause foF2 to increase at daytime and decrease at nighttime. The intensity of response of foF2 is stronger at nighttime than that at daytime; (2) seasonal dependence of the response of foF2 to geomagnetic activity is very obvious. The negative ionospheric storm effect is the strongest in summer and the positive ionospheric storm effect is the strongest in winter; (3) the solar cycle has important effect on the response of foF2 to geomagnetic activity in Hainan. In high solar activity period, the diurnal variation of the response of foF2 is very pronounced in each season, and the strong ionospheric response can last several days. In low solar activity period, ionospheric response has very pronounced diurnal variation in winter only; (4) the local time of geomagnetic activities occurring also has important effect on the responses of foF2 in Hainan. Generally, geomagnetic activities occurred at nighttime can cause stronger and longer responses of foF2 in Hainan.     

Global validation of IRI TEC for high and medium solar activity conditions

The IRI model offers a choice of options for the computation of the electron density profile and electron content (TEC). Recently new options for the topside electron density profile have been developed, which have a strong impact on TEC. Therefore it is important to test massively the IRI and the new options with experimental data. A large number of permanent stations record dual frequency GPS data from which it is possible to obtain TEC values. Thirty-one worldwide distributed stations have been selected to investigate the capabilities of the IRI to reproduce experimental TEC. Data for years 2000 (high solar activity) and 2004 (medium solar activity) have been analyzed computing modeled values with the IRI-2001 and the IRI-2007-NeQuick topside options. It is found that IRI-2007-NeQuick option generally improves the estimate of the slant TEC, especially in the case of high latitudes stations during high solar activity.     

An investigation of ionospheric F region response in the Brazilian sector to the super geomagnetic storm of May 2005

In this paper, we have investigated the responses of the ionospheric F region at equatorial and low latitude regions in the Brazilian sector during the super geomagnetic storm on 15–16 May 2005. The geomagnetic storm reached a minimum Dst of −263 nT at 0900 UT on 15 May. In this paper, we present vertical total electron content (vTEC) and phase fluctuations (in TECU/min) from Global Positioning System (GPS) observations obtained at Belém, Brasília, Presidente Prudente, and Porto Alegre, Brazil, during the period 14–17 May 2005. Also, we present ionospheric parameters h’F, hpF2, and foF2, using the Canadian Advanced Digital Ionosonde (CADI) obtained at Palmas and São José dos Campos, Brazil, for the same period. The super geomagnetic storm has fast decrease in the Dst index soon after SSC at 0239 UT on 15 May. It is a good possibility of prompt penetration of electric field of magnetospheric origin resulting in uplifting of the F region. The vTEC observations show a trough at BELE and a crest above UEPP, soon after SSC, indicating strengthening of nighttime equatorial anomaly. During the daytime on 15 and 16 May, in the recovery phase, the variations in foF2 at SJC and the vTEC observations, particularly at BRAZ, UEPP, and POAL, show large positive ionospheric storm. There is ESF on the all nights at PAL, in the post-midnight (UT) sector, and phase fluctuations only on the night of 14–15 May at BRAZ, after the SSC. No phase fluctuations are observed at the equatorial station BELE and low latitude stations (BRAZ, UEPP, and POAL) at all other times. This indicates that the plasma bubbles are generated and confined on this magnetically disturbed night only up to the low magnetic latitude and drifted possibly to west.     

Low solar activity variability and IRI 2007 predictability of equatorial Africa GPS TEC

Diurnal, seasonal and latitudinal variations of Vertical Total Electron Content (VTEC) over the equatorial region of the African continent and a comparison with IRI-2007 derived TEC (IRI-TEC), using all three options (namely; NeQuick, IRI01-corr and IRI-2001), are presented in this paper. The variability and comparison are presented for 2009, a year of low solar activity, using data from thirteen Global Positioning System (GPS) receivers. VTEC values were grouped into four seasons namely March Equinox (February, March, April), June Solstice (May, June, July), September Equinox (August, September, October), and December Solstice (November, December, January). VTEC generally increases from 06h00 LT and reaches its maximum value at approximately 15h00–17h00 LT during all seasons and at all locations. The NeQuick and IRI01-corr options of the IRI model predict reasonably well the observed diurnal and seasonal variation patterns of VTEC values. However, the IRI-2001 option gave a relatively poor prediction when compared with the other options. The post-midnight and post-sunset deviations between modeled and observed VTEC could arise because NmF2 or the shape of the electron density profile, or both, are not well predicted by the model; hence some improvements are still required in order to obtain improved predictions of TEC over the equatorial region of the Africa sector.     

Spatial distribution of TEC across India in 2005: Seasonal asymmetries and IRI prediction

Total electron content measured simultaneously at 10 locations over India during the low solar activity year 2005 is used to examine the temporal and spatial asymmetries and also to assess the predictability of the International Reference Ionosphere in respect of the observed asymmetrical distribution. The stations are distributed in latitude along 77°E and in longitude along 23°N forming a meridional and a zonal chain respectively. A longitudinal gradient positive towards east was observed in the daytime hours of equinox and summer. Equinoctial asymmetry was prevalent across India during this year. Within the crest and equator, winter anomaly has been observed. It is found that IRI 2012 (with Ne Quick option, URSI coefficients) is unable to fully capture the temporal variation and spatial gradients of the ionization density in the Indian sector during 2005. The amount of offset between the model and measurement varies with local time and location.     

Investigation of ionospheric irregularities and scintillation using TEC at high latitude

The occurrence of ionospheric irregularities at high latitudes, with dimensions of several kms down to decameter scale size shows strong correlation with geomagnetic disturbance, season and solar activity. Transionospheric radio waves propagating through these irregularities experience rapid random fluctuations in phase and/or amplitude of the signal at the receiver, termed scintillation, which can degrade GNSS services. Thus, investigation and prediction of this scintillation effect is very important. To investigate such scintillation effects, a GISTM (GPS Ionospheric Scintillation and TEC Monitoring) NovAtel dual frequency (L1/L2) GPS receiver has been installed at Trondheim, Norway (63.41°$63.41°$ N, 10.4°$10.4°$ E), capable of collecting scintillation indices at a 1 min rate as well as the raw data (phase and intensity) of the satellite signals at a 50 Hz sampling rate and TEC (Total Electron Content) at a 1 Hz rate. Many researchers have reported that both phase and amplitude scintillation is closely associated with TEC fluctuations or associated with a significant developing enhancement or depletion in the TEC. In this study, a novel analogous phase index is developed which provides samples at a 1 min rate. Generally the scintillation indices can help in estimating the irregularity scintillation effect at a one minute rate, but such procedures are time consuming if DFTs of the phase and/or amplitude at a 50 Hz data are required. In this study, instead, this analogous phase index is estimated from 1 Hz rate TEC values obtained from the raw signals and is then compared for weak, moderate and strong scintillation at Trondheim for one year of data collected from the installed GPS receiver. The spectral index of the irregularities (that is the inverse power law of their spatial spectrum) is determined from the resultant phase scintillation psd. The correlations of the scintillation indices and spectral indices with the analogous phase index have been investigated under different geomagnetic conditions (represented by the Kp index) and an approximate linear correlation of phase scintillation with the analogous phase index was found. Then a principal advantage of this index is that it achieves this correlation without requiring a high sampling data rate and the need for DFTs. Thus, the index seems a good candidate for developing a simple means of ionospheric scintillation prediction which could also be utilized in the development of alerts using regional mappings.     

Effects of geomagnetic storms on the bottomside ionospheric F region

The geomagnetic storm is a complex process of solar wind/magnetospheric origin. The variability of the ionospheric parameters increases substantially during geomagnetic storms initiated by solar disturbances. Various features of geomagnetic storm act at various altitudes in the ionosphere and neutral atmosphere. The paper deals with variability of the electron density of the ionospheric bottomside F region at every 10 km of altitude during intense geomagnetic storms with attention paid mainly to the distribution of the F1 region daytime ionisation. We have analysed all available electron density profiles from some European middle latitude stations (Chilton, Pruhonice, Ebro, Arenosillo, Athens) for 36 events that occurred in different seasons and under different levels of solar activity (1995–2003). Selected events consist of both depletion and increase of the F2 region electron density. For European higher middle and middle latitude the F1 region response to geomagnetic storm was found to be negative (decrease of electron density) independent on the storm effect on the F2 region. For lower middle latitude the F1 response is weaker and less regular. Results of the analysis also show that the maximum of the storm effect may sometimes occur below the height of the maximum of electron density (NmF2).     

Regional 4-D modeling of the ionospheric electron density from satellite data and IRI

Accurate knowledge of the electron density is the key point in correcting ionospheric delays of electromagnetic measurements and in studying ionosphere physics. During the last decade Global Navigation Satellite Systems (GNSS) have become a promising tool for monitoring ionospheric parameters such as the total electron content (TEC). In this contribution we present a four-dimensional (4-D) model of the electron density consisting of a given reference part, i.e., the International Reference Ionosphere (IRI), and an unknown correction term expanded in terms of multi-dimensional base functions. The corresponding series coefficients are calculable from the satellite measurements by applying parameter estimation procedures. Since satellite data are usually sampled between GPS satellites and ground stations, finer structures of the electron density are modelable just in regions with a sufficient number of ground stations. The proposed method is applied to simulated geometry-free GPS phase measurements. The procedure can be used, for example, to study the equatorial anomaly.     

TEC derived from some GPS stations in Nigeria and comparison with the IRI and NeQuick models

Total electron content (TEC) measured simultaneously using Global Positioning System (GPS) ionospheric monitors installed at some locations in Nigeria during the year 2011 (Rz = 55.7) was used to study the diurnal, seasonal, and annual TEC variations. The TEC exhibits daytime maximum, seasonal variation and semiannual variations. Measured TEC were compared with those predicted by the improved versions of the International Reference Ionosphere (IRI) and NeQuick models. The models followed the diurnal and seasonal variation patterns of the observed values of TEC. However, IRI model produced better estimates of TEC than NeQuick at all locations.