Ionospheric Electron Content Research Articles

Improvement of the Ionospheric Radio Occultation Retrievals by Means of Accurate Global Ionospheric Maps

AbstractThe Abel inversion method is the classic approach to retrieve electron density profiles from radio occultation measurements, assuming that the electron density is only dependent upon height. In the early 2000s, the total electron content (TEC)‐aided inversion method was introduced, where horizontal gradients in the electron density distribution are taken into account by using information from an external model of total electron content (TEC). We show how the quality of the chosen external ionospheric model influences the radio occultation retrievals. Using a two‐layer TEC model, the 68% percentiles of the global error distributions of retrieved critical frequency and bottomside ionospheric electron content for 2014 improve more than 30%, with regard to the classic Abel inversion results. In a well‐sounded region like Europe, this improvement raises to about 40–45%, while 10% improvement is achieved using two‐layer maps with regard to the retrievals using single‐layer TEC maps. We point out that using the TEC to describe the horizontal gradient incurs an implicit mismodeling that, until now, has not been taken into account. A new technique is presented that, thanks to TEC modeling by means of two layers, can assess the impacts of such mismodeling. From this technique, a method to select the most accurate radio occultation retrievals is applied to demonstrate that the resulting errors in the retrieved critical frequencies for the European region are smaller than 7% in nearly 82% of the cases, very similar to the relative difference between measurements of two nearby ionosondes.

The 6 September 2017 X-Class Solar Flares and Their Impacts on the Ionosphere, GNSS, and HF Radio Wave Propagation.

On 6 September 2017, the Sun emitted two significant solar flares (SFs). The first SF, classified X2.2, peaked at 09:10 UT. The second one, X9.3, which is the most intensive SF in the current solar cycle, peaked at 12:02 UT and was accompanied by solar radio emission. In this work, we study ionospheric response to the two X‐class SFs and their impact on the Global Navigation Satellite Systems and high‐frequency (HF) propagation. In the ionospheric absolute vertical total electron content (TEC), the X2.2 SF caused an overall increase of 2–4 TECU on the dayside. The X9.3 SF produced a sudden increase of ~8–10 TECU at midlatitudes and of ~15–16 TECU enhancement at low latitudes. These vertical TEC enhancements lasted longer than the duration of the EUV emission. In TEC variations within 2–20 min range, the two SFs provoked sudden increases of ~0.2 TECU and 1.3 TECU. Variations in TEC from geostationary and GPS/GLONASS satellites show similar results with TEC derivative of ~1.3–1.7 TECU/min for X9.3 and 0.18–0.24 TECU/min for X2.2 in the subsolar region. Further, analysis of the impact of the two SFs on the Global Navigation Satellite Systems‐based navigation showed that the SF did not cause losses‐of‐lock in the GPS, GLONASS, or Galileo systems, while the positioning error increased by ~3 times in GPS precise point positioning solution. The two X‐class SFs had an impact on HF radio wave propagation causing blackouts at <30 MHz in the subsolar region and <15 MHz in the postmidday sector.

Precise ionospheric electron content monitoring from single-frequency GPS receivers

The number of existing global positioning system (GPS) single-frequency receivers continues growing. More than 90% of GPS receivers are implemented as low-cost single-frequency chipsets embedded in smartphones. This provides new opportunities, in particular for ionospheric sounding. In this context, we present the new sidereal days ionospheric graphic (SIg) combination of single-frequency GNSS measurements. SIg is able to monitor, for each given GNSS transmitter–receiver pair, the vertical total electron content (VTEC) relative to the previous observation with the same or almost the same line-of-sight (LOS) vector. In such arrangements the SIg multipath error mostly cancels, thus increasing the accuracy of the ΔVTEC significantly. This happens for the GPS constellation after one sidereal day (about 23 h 56 m) and for Galileo after 10 sidereal days approximately. Moreover, we show that the required calibration of the corresponding carrier phase ambiguity can be accurately performed by means of VTEC global ionospheric maps (GIMs). The results appear almost as accurate as those based on the dual-frequency technique, i.e., about 1 TECU or better, and with much more precision and resolution than the GIM values in the ionospheric region sounded by each given single-frequency receiver. The performance is demonstrated using actual data from 9 permanent GPS receivers during a total solar eclipse on August 21, 2017 over North America, where the corresponding ionospheric footprint is clearly detected in agreement with the total solar eclipse predictions. The advantages of extending SIg to lower carrier frequencies and the feasibility of applying it to other global navigation satellite system (GNSS) systems are also studied. This is shown in terms of a fully consistent VTEC depletion signature of the same eclipse phenomena, obtained with Galileo-only data in North America at mid and low latitude. Finally the SIg feasibility, including the cycle slip detection, is shown as well with actual mass-market single frequency GPS receivers at mid and high latitude.

Ionospheric Electron Content During Solar Cycle 23

AbstractThe solar minimum period between solar cycles 23 and 24 has generated extensive study, in part because the changes in the upper atmosphere and ionosphere were so dramatic. Thermospheric density during 2008–2009, at the benchmark 400‐km altitude, was shown to be approximately 30% lower than the previous solar minimum during 1996, in observations and in model simulations. Ionospheric densities were also estimated to be lower, by 10% to 20%. However, earlier analyses using data from the global network of Global Navigation Satellite System receivers did not find any significant reduction in global mean total electron content. Recent work reexamining those data using a limited set of receivers, but with consistent processing, identified a 19% reduction (Emmert et al., 2017, https://doi.org/10.1002/2016JA023680). In this study, we compare model simulations of the thermosphere‐ionosphere system during the 2008–2009 solar minimum to its predecessor and estimate a 16% reduction in global mean total electron content. This is consistent with the neutral thermosphere density change seen in model simulations and in satellite drag observations. In the model simulations, the primary driver is a 10% reduction in solar extreme ultraviolet irradiance, with a smaller contribution from lower geomagnetic activity, and a very small decrement due to increase in atmospheric carbon dioxide.

Using the Subsurface Soil Sounding Radar for Investigating the Structure and Total Electron Content of the Martian Ionosphere

The possibilities of using the Martian soil subsurface sounding radar for investigating the structure of the plasma shell surrounding the planet have been considered. Based on the numerical modeling results and actual soil sounding data, it has been shown that the soil sounding mode of the radio-locating MARSIS radar can be used to assess the structure of the Martian ionosphere. As the emitted signals pass to the planet’s surface, it is possible to use the reflected signals to estimate the total electron content of the Martian ionosphere along the flight track of the spacecraft.

Plasmaspheric Electron Content Inferred from Residuals between GNSS-Derived and TOPEX/JASON Vertical TEC Data

The plasmasphere, which is located above the ionosphere, is a significant component of Earth’s atmosphere, and the plasmasphere electron content (PEC) distribution is determined by different physical mechanisms to those of the ionosphere electron content (IEC). However, the observation for the PEC is very limited. In this study, we introduced a methodology (called zero assumption method, which is based on the assumption that PEC can reach zero) to extract the PEC over TOPEX/JASON (T/J) and global navigation satellite system (GNSS) overlapping areas. Results show that the daily systematic bias (T/J vertical TEC &gt; GNSS-derived vertical TEC) for both low (2009) and high (2011) solar activity condition is consistent, and the systematic bias for JASON2 and JASON1 is different. We suggest that systematic biases predominantly arise from the sea state bias (SSB), especially the tracker bias. After removing the systematic bias, we extracted reliable PEC inferred from differences between GNSS-derived vertical TEC and T/J vertical TEC data. Finally, the characteristics of the plasmaspheric component distribution for different local times, latitudes, and seasons were investigated.

Analysis of temporal\u2013longitudinal\u2013latitudinal characteristics in the global ionosphere based on tensor rank-1 decomposition

Combining analyses of spatial and temporal characteristics of the ionosphere is of great significance for scientific research and engineering applications. Tensor decomposition is performed to explore the temporal–longitudinal–latitudinal characteristics in the ionosphere. Three-dimensional tensors are established based on the time series of ionospheric vertical total electron content maps obtained from the Centre for Orbit Determination in Europe. To obtain large-scale characteristics of the ionosphere, rank-1 decomposition is used to obtain U^{(1)}, U^{(2)}, and U^{(3)}, which are the resulting vectors for the time, longitude, and latitude modes, respectively. Our initial finding is that the correspondence between the frequency spectrum of U^{(1)} and solar variation indicates that rank-1 decomposition primarily describes large-scale temporal variations in the global ionosphere caused by the Sun. Furthermore, the time lags between the maxima of the ionospheric U^{(2)} and solar irradiation range from 1 to 3.7 h without seasonal dependence. The differences in time lags may indicate different interactions between processes in the magnetosphere–ionosphere–thermosphere system. Based on the dataset displayed in the geomagnetic coordinates, the position of the barycenter of U^{(3)} provides evidence for north–south asymmetry (NSA) in the large-scale ionospheric variations. The daily variation in such asymmetry indicates the influences of solar ionization. The diurnal geomagnetic coordinate variations in U^{(3)} show that the large-scale EIA (equatorial ionization anomaly) variations during the day and night have similar characteristics. Considering the influences of geomagnetic disturbance on ionospheric behavior, we select the geomagnetic quiet GIMs to construct the ionospheric tensor. The results indicate that the geomagnetic disturbances have little effect on large-scale ionospheric characteristics.

27-day variation in solar-terrestrial parameters: Global characteristics and an origin based approach of the signals

The Earth and the near interplanetary medium are affected by the Sun in different ways. Those processes generated in the Sun that induce perturbations into the Magnetosphere-Ionosphere system are called geoeffective processes and show a wide range of temporal variations, like the 11-year solar cycle (long term variations), the variation of ∼27 days (recurrent variations), solar storms enduring for some days, particle acceleration events lasting for some hours, etc.In this article, the periodicity of ∼27 days associated with the solar synodic rotation period is investigated. The work is mainly focused on studying the resulting 27-day periodic signal in the magnetic activity, by the analysis of the horizontal component of the magnetic field registered on a set of 103 magnetic observatories distributed around the world. For this a new method to isolate the periodicity of interest has been developed consisting of two main steps: the first one consists of removing the linear trend corresponding to every calendar year from the data series, and the second one of removing from the resulting series a smoothed version of it obtained by applying a 30-day moving average. The result at the end of this process is a data series in which all the signal with periods larger than 30 days are canceled.The most important characteristics observed in the resulting signals are two main amplitude modulations: the first and most prominent related to the 11-year solar cycle and the second one with a semiannual pattern. In addition, the amplitude of the signal shows a dependence on the geomagnetic latitude of the observatory with a significant discontinuity at approx. ±60°.The processing scheme was also applied to other parameters that are widely used to characterize the energy transfer from the Sun to the Earth: F10.7 and Mg II indices and the ionospheric vertical total electron content (vTEC) were considered for radiative interactions; and the solar wind velocity for the non-radiative interactions between the solar wind and the magnetosphere. The 27-day signal obtained in the magnetic activity was compared with the signals found in the other parameters resulting in a series of cross-correlations curves with maximum correlation between 3 and 5 days of delays for the radiative and between 0 and 1 days of delay for the non-radiative parameters. This result supports the idea that the physical process responsible for the 27-day signal in the magnetic activity is related to the solar wind and not to the solar electromagnetic radiation.

The Effect of Subauroral Polarization Streams on Ionosphere and Thermosphere During the 2015 St. Patrick's Day Storm: Global Ionosphere‐Thermosphere Model Simulations

AbstractUsing the Millstone Hill incoherent scatter radar observations during 2015 St. Patrick's Day storm, subauroral polarization streams (SAPSs) have been specified in the nonhydrostatic Global Ionosphere‐Thermosphere Model simulations. The results reveal that the effect of SAPS on the coupled thermosphere‐ionosphere system includes the following: (1) Sudden frictional heating of SAPS results in acoustic‐gravity waves in the thermosphere. The vertical oscillation is localized, while the meridional disturbance propagates poleward and equatorward. (2) The SAPS‐associated horizontal wind field includes an enhanced westward wind within SAPS channel and a twin of vortex‐like winds north (clockwise) and south (anticlockwise) of subauroral latitudes. (3) Due to the neutral‐ion drag, ions in the vicinity of SAPS channel oscillate vertically with neutrals, resulting in a perturbation of ~0.3 TECu in ionospheric total electron content. The SAPS‐induced traveling atmospheric disturbances can elevate the plasma and increase the total electron content in midlatitude ionosphere. (4) It is confirmed that the Coriolis force can contribute to the poleward turning of the neutral wind during the post‐SAPS interval. In addition, the traveling atmospheric disturbance induced by the variation of auroral input and high‐latitude convection is possibly the primary cause of the poleward neutral wind surge during the magnetic storm on 17–18 March 2015. The combination of the two factors can make the northward meridional wind surge reach a magnitude of 100 m/s. This study improves our understanding of the SAPS's effect on neutral dynamics and ion‐neutral coupling processes during geomagnetically disturbed intervals.

Prediction of global ionospheric VTEC maps using an adaptive autoregressive model

In this contribution, an adaptive autoregressive model is proposed and developed to predict global ionospheric vertical total electron content maps (VTEC). Specifically, the spherical harmonic (SH) coefficients are predicted based on the autoregressive model, and the order of the autoregressive model is determined adaptively using the F-test method. To test our method, final CODE and IGS global ionospheric map (GIM) products, as well as altimeter TEC data during low and mid-to-high solar activity period collected by JASON, are used to evaluate the precision of our forecasting products. Results indicate that the predicted products derived from the model proposed in this paper have good consistency with the final GIMs in low solar activity, where the annual mean of the root-mean-square value is approximately 1.5 TECU. However, the performance of predicted vertical TEC in periods of mid-to-high solar activity has less accuracy than that during low solar activity periods, especially in the equatorial ionization anomaly region and the Southern Hemisphere. Additionally, in comparison with forecasting products, the final IGS GIMs have the best consistency with altimeter TEC data. Future work is needed to investigate the performance of forecasting products using the proposed method in an operational environment, rather than using the SH coefficients from the final CODE products, to understand the real-time applicability of the method.Graphical The ionosphere affects national defense, aerospace, economic development, and human life. Short-term predictions of global ionospheric vertical TEC (VTEC) maps are important for the scientific analysis of the ionosphere and practical applications, such as satellite navigation. In this contribution, an adaptive autoregressive model is proposed and developed for the prediction of global ionospheric vertical TEC maps. The forecasting products are compared with final CODE and IGS GIMs, as well as altimeter TEC from JASON data under both low and mid-to-high solar activity. The results indicate that forecasting products exhibit good consistency with final GIMs.

Real-Time Precise Point Positioning (RTPPP) with raw observations and its application in real-time regional ionospheric VTEC modeling

Precise Point Positioning (PPP) is an absolute positioning technology mainly used in post data processing. With the continuously increasing demand for real-time high-precision applications in positioning, timing, retrieval of atmospheric parameters, etc., Real-Time PPP (RTPPP) and its applications have drawn more and more research attention in recent years. This study focuses on the models, algorithms and ionospheric applications of RTPPP on the basis of raw observations, in which high-precision slant ionospheric delays are estimated among others in real time. For this purpose, a robust processing strategy for multi-station RTPPP with raw observations has been proposed and realized, in which real-time data streams and State-Space-Representative (SSR) satellite orbit and clock corrections are used. With the RTPPP-derived slant ionospheric delays from a regional network, a real-time regional ionospheric Vertical Total Electron Content (VTEC) modeling method is proposed based on Adjusted Spherical Harmonic Functions and a Moving-Window Filter. SSR satellite orbit and clock corrections from different IGS analysis centers are evaluated. Ten globally distributed real-time stations are used to evaluate the positioning performances of the proposed RTPPP algorithms in both static and kinematic modes. RMS values of positioning errors in static/kinematic mode are 5.2/15.5, 4.7/17.4 and 12.8/46.6 mm, for north, east and up components, respectively. Real-time slant ionospheric delays from RTPPP are compared with those from the traditional Carrier-to-Code Leveling (CCL) method, in terms of function model, formal precision and between-receiver differences of short baseline. Results show that slant ionospheric delays from RTPPP are more precise and have a much better convergence performance than those from the CCL method in real-time processing. 30 real-time stations from the Asia-Pacific Reference Frame network are used to model the ionospheric VTECs over Australia in real time, with slant ionospheric delays from both RTPPP and CCL methods for comparison. RMS of the VTEC differences between RTPPP/CCL method and CODE final products is 0.91/1.09 TECU, and RMS of the VTEC differences between RTPPP and CCL methods is 0.67 TECU. Slant Total Electron Contents retrieved from different VTEC models are also validated with epoch-differenced Geometry-Free combinations of dual-frequency phase observations, and mean RMS values are 2.14, 2.33 and 2.07 TECU for RTPPP method, CCL method and CODE final products, respectively. This shows the superiority of RTPPP-derived slant ionospheric delays in real-time ionospheric VTEC modeling.

Tsallis non-extensive statistical mechanics in the ionospheric detrended total electron content during quiet and storm periods

In this paper, the Tsallis non-extensive q-statistics in ionospheric dynamics was investigated using the total electron content (TEC) obtained from two Global Positioning System (GPS) receiver stations. This investigation was carried out considering the geomagnetically quiet and storm periods. The micro density variation of the ionospheric total electron content was extracted from the TEC data by method of detrending. The detrended total electron content, which represent the variation in the internal dynamics of the system was further analyzed using for non-extensive statistical mechanics using the q-Gaussian methods. Our results reveals that for all the analyzed data sets the Tsallis Gaussian probability distribution (q-Gaussian) with value q>1 were obtained. It was observed that there is no distinct difference in pattern between the values of qquiet and qstorm. However the values of q varies with geophysical conditions and possibly with local dynamics for the two stations. Also observed are the asymmetric pattern of the q-Gaussian and a highly significant level of correlation for the q-index values obtained for the storm periods compared to the quiet periods between the two GPS receiver stations where the TEC was measured. The factors responsible for this variation can be mostly attributed to the varying mechanisms resulting in the self-reorganization of the system dynamics during the storm periods. The result shows the existence of long range correlation for both quiet and storm periods for the two stations.

Collocated ionosonde and dense GPS/GLONASS network measurements of midlatitude MSTIDs

To analyze midlatitude medium-scale travelling ionospheric disturbances (MSTIDs) over Kazan (55.5°N, 49°E), Russia, the sufficiently dense network of GNSS receivers (more than 150 ground-based stations) were used. For the first time, daytime MSTIDs in the form of their main signature (band structure) on high-resolution two-dimensional maps of the total electron content perturbation (TEC maps) are compared with ionosonde data with a high temporal resolution. For a pair of events, a relationship between southwestward TEC perturbations and evolution of F2 layer traces was established. So F2 peak frequency varied in antiphase to TEC perturbations. The ionograms show that during the movement of plasma depletion band (overhead ionosonde) the F2 peak frequency is the highest, and vice versa, for the plasma enhancement band, the F2 peak frequency is the lowest. One possible explanation may be a greater inclination of the radio beam from the vertical during the placement of a plasma enhancement band above the ionosonde, as evidenced by the absence of multiple reflections and the increased occurrence rate of additional cusp trace. Another possible explanation may be the redistribution of the electron content in the topside ionosphere with a small decrease in the F peak concentration of the layer with a small increase in TEC along the line-of-sight. Analysis of F2 peak frequency variation shows that observed peak-to-peak values of TEC perturbation equal to 0.4 and 1 TECU correspond to the values of ΔN/N equal to 13% and 28%. The need for further research is evident.

Global Ionospheric and Thermospheric Effects of the June 2015 Geomagnetic Disturbances: Multi-Instrumental Observations and Modeling.

By using data from multiple instruments, we investigate ionospheric/thermospheric behavior during the period from 21 to 23 June 2015, when three interplanetary shocks (IS) of different intensities arrived at Earth. The first IS was registered at 16:45 UT on 21 June and caused ~50 nT increase in the SYM‐H index. The second IS arrived at 5:45 UT on 22 June and induced an enhancement of the auroral/substorm activity that led to rapid increase of thermospheric neutral mass density and ionospheric vertical total electron content at high latitudes. Several hours later, topside electron content and electron density increased at low latitudes on the nightside. The third and much larger IS arrived at 18:30 UT on 22 June and initiated a major geomagnetic storm that lasted for many hours. The storm provoked significant effects in the thermosphere and ionosphere on both dayside and nightside. In the thermosphere, the dayside neutral mass density exceeded the quiet time levels by 300–500%, with stronger effects in the summer hemisphere. In the ionosphere, both positive and negative storm effects were observed on both dayside and nightside. We compared the ionospheric observations with simulations by the coupled Sami3 is Also a Model of the Ionosphere/Rice Convection Model (SAMI3/RCM) model. We find rather good agreement between the data and the model for the first phase of the storm, when the prompt penetration electric field (PPEF) was the principal driver. At the end of the storm main phase, when the ionospheric effects were, most likely, driven by a combination of PPEF and thermospheric winds, the modeling results agree less with the observations.

Analysis of ionospheric vertical total electron content before the 1 April 2014 Mw 8.2 Chile earthquake

This paper studies ionospheric vertical total electron content (VTEC) variations before the 1 April 2014 Mw 8.2 Chile earthquake. VTEC derived from 14 global positioning system (GPS) stations and global ionospheric map (GIM) were used to analyze ionospheric variations before the earthquake using the sliding interquartile range method, and the results showed that significant positive VTEC anomalies occurred on 28 March. To explore possible causes of these anomalies, effects of solar and geomagnetic activities were examined, and VTEC variations during 17 March to 31 March in 2009–2013 were cross-compared. Also, VTEC for a full year before the earthquake was investigated. The results indicated that the anomalies were weakly associated with high solar activities and geomagnetic storms and that the anomalies were not normal seasonal and diurnal variations. An analysis of the spatial distribution of the observed anomalies was also presented, and it demonstrated that the anomalies specifically appeared around the epicenter on 28 March. We suggest that the observed anomalies may be associated with the subsequent Chile earthquake. Equatorial anomaly variations were analyzed to discuss the possible physical mechanism, and results showed that the equatorial anomaly unusually increased on 28 March, which indicates that anomalous electric fields generated in the earthquake preparation area and the meridional wind are possible causes of the observed ionospheric anomalies.

Solar activity indices as a proxy for the variation of ionospheric Total Electron Content (TEC) over Bahir Dar, Ethiopia during the year 2010–2014

This paper was aimed at investigating the solar variations of vTEC as a function of solar activity parameters, EUV and F10.7 radio flux. The daily values of ionospheric vertical Total Electron Content (vTEC) were observed using a dual frequency GPS receiver deployed at Bahir Dar (11.6°N and 37.36°E), Ethiopia. Measurements were taken during the period of 2010–2014 for successive five years and analysis was done on only quiet day observations. A quadratic fit was used as a model to describe the daily variation of vTEC in relation to solar parameters. Linear and non-linear coefficients of the vTEC variations were calculated in order to capture the trend of the variation. The variation of vTEC have showed good agreement with the trend of solar parameters in almost all of the days we consider during the period of our observations. We have explicitly observed days with insignificant TECU deviation (eg. modeling with respect to EUV, DOY=49 in 2010 and modeling with respect to F10.7, DOY=125 in 2012 and the like) and days with maximum deviation (about 50 TECU). A maximum deviation were observed, on average, during months of equinox whereas minimum during solstice months. This implies that there is a need to consider more parameters, including EUV and F10.7, that can affect the variation of vTEC during equinox seasons. Relatively, small deviations was observed in modeling vTEC as a function of EUV compared to that of the variation due to F10.7cm flux. This may also tell us that EUV can be more suitable in modeling the solar variation of vTEC especially for longterm trends. Even though, the linear trend of solar variations of vTEC was frequently observed, significant saturation and amplification trends of the solar variations of vTEC were also observed to some extent across the months of the years we have analyzed. This mixed trend of the solar variation of vTEC implies the need for thorough investigation on the effect of solar parameters on TEC. However, based on long-term dataset, we came to conclude that the solar variations of vTEC is dominated by its linear pattern.

Ionospheric TEC modelling for earthquakes precursors from GNSS data

In present study, we measured the Total Electron Content (TEC) variation in the ionosphere from Global Navigation Satellite System (GNSS) data which might have been induced by earthquakes in the Himalayan region. The results were analysed with other inducing factors (geomagnetic storm and solar flare) affecting TEC in order to constraint the causative factor. The study has been performed to understand a relationship between ionospheric electron content and earthquakes occurrences with special emphasis on Himalayan region and provides information on spatio-temporal variation of TEC from GNSS observation stations vis-à-vis prominent earthquakes of the region. The results indicate that the ground based GNSS (GPS) observations show the deviation in vertical total electron content (vTEC) in ionosphere few days prior to the seismic event as evident from our two continuously operating reference stations (CORS) as well as from CORS of UNAVCO data centre. Ionospheric perturbation has also been observed in case of low magnitude earthquakes (Mw 4.9 in present study) whenever recording station lies very close to the epicentre. TEC variation is found to increase as the epicentre distances decreases. In case of Mw 7.8, 2015 Nepal earthquake the TEC variation is found to increase by 15–20 TEC unit recorded at station separated by 60 km apart. This may provide us the avenue for epicentre detection as TEC concentration was found to increase as we move closer to the epicentre. TEC variations (mostly high TEC) have been observed during a period of 0–8 days prior to 4 earthquakes: 1st April 2015 Pipalkoti earthquake (4.9Mw), 25th April 2015 Nepal earthquake (7.8Mw), 26 April 2015 (6.7Mw) and 12th May 2015 (7.3Mw) Nepal earthquakes. Significant low TEC values were also observed before 13–14 days prior to first two earthquake events. Overall the study has revealed that low TEC followed by couple of high TEC values correlate well with the seismic events in Himalayan region.

Methodology and consistency of slant and vertical assessments for ionospheric electron content models

A summary of the main concepts on global ionospheric map(s) [hereinafter GIM(s)] of vertical total electron content (VTEC), with special emphasis on their assessment, is presented in this paper. It is based on the experience accumulated during almost two decades of collaborative work in the context of the international global navigation satellite systems (GNSS) service (IGS) ionosphere working group. A representative comparison of the two main assessments of ionospheric electron content models (VTEC-altimeter and difference of Slant TEC, based on independent global positioning system data GPS, dSTEC-GPS) is performed. It is based on 26 GPS receivers worldwide distributed and mostly placed on islands, from the last quarter of 2010 to the end of 2016. The consistency between dSTEC-GPS and VTEC-altimeter assessments for one of the most accurate IGS GIMs (the tomographic-kriging GIM ‘UQRG’ computed by UPC) is shown. Typical error RMS values of 2 TECU for VTEC-altimeter and 0.5 TECU for dSTEC-GPS assessments are found. And, as expected by following a simple random model, there is a significant correlation between both RMS and specially relative errors, mainly evident when large enough number of observations per pass is considered. The authors expect that this manuscript will be useful for new analysis contributor centres and in general for the scientific and technical community interested in simple and truly external ways of validating electron content models of the ionosphere.

Determinations of ionosphere and plasmasphere electron content for an African chain of GPS stations

Abstract. The confluence of recent instrumentation deployments in Africa with developments for the determination of plasmasphere electron content using Global Positioning System (GPS) receivers has provided new opportunities for investigations in that region. This investigation, using a selected chain of GPS stations, extends the method (SCORPION) previously applied to a chain of GPS stations in North America in order to separate the ionosphere and plasmasphere contributions to the total electron content (TEC) during a day (24 July) in 2011. The results span latitudes from the southern tip of Africa, across the Equator, to the southern Arabian Peninsula, providing a continuous latitudinal profile for both the ionosphere and plasmasphere during this day.The peak diurnal vertical ionosphere electron content (IEC) increases from about 14 TEC units (1 TEC unit = 1016 electrons m−2) at the southernmost station to about 32 TEC units near the geographic equator, then decreases to about 28 TEC units at the Arabian Peninsula. The peak diurnal slant plasmasphere electron content (PEC) varies between about 4 and 7 TEC units among the stations, with a local latitudinal profile that is significantly influenced by the viewing geometry at the station location, relative to the magnetic field configuration. In contrast, the peak vertical PEC varies between about 1 and 6 TEC units among the stations, with a more uniform latitudinal variation.Comparisons to other GPS data analyses are also presented for TEC, indicating the influence of the PEC on the determination of latitudinal TEC variations and also on the absolute TEC levels, by inducing an overestimate of the receiver bias. The derived TEC latitudinal profiles, in comparison to global map profiles, tend to differ from the map results only about as much as the map results differ among themselves. A combination of ionosonde IEC and alternative GPS TEC measurements, which in principle permits a PEC determination through their difference, was compared to the composite and separate ionosphere and plasmasphere contributions derived solely by the SCORPION method for one station. Although there is considerably more scatter in the PEC values derived from the difference of the GPS TEC and ionosonde IEC measurements compared to the PEC values derived by the SCORPION method, the average overhead values for this day are comparable for the two methods, near 2 TEC units, at the South African site examined.This initial investigation provides a basis for day-to-day TEC monitoring for Africa, with separate ionosphere and plasmasphere electron content determinations.

A statistical study of global ionospheric map total electron content changes prior to occurrences of M ≥ 6.0 earthquakes during 2000–2014

AbstractThere are many reports on the occurrence of anomalous changes in the ionosphere prior to large earthquakes. However, whether or not these changes are reliable precursors that could be useful for earthquake prediction is controversial within the scientific community. To test a possible statistical relationship between ionospheric disturbances and earthquakes, we compare changes in the total electron content (TEC) of the ionosphere with occurrences of 1279 M ≥ 6.0 earthquakes globally for 2000–2014. We use TEC data from the global ionosphere map (GIM) and an earthquake list declustered for aftershocks. For each earthquake, we look for anomalous changes in GIM‐TEC within 2.5° latitude and 5.0° longitude of the earthquake location (the spatial resolution of GIM‐TEC). Although case studies of individual earthquakes that used short periods of data sometimes yield GIM‐TEC changes that were considered possible earthquake‐related phenomena, our analysis has not found any statistically significant changes prior to earthquakes when considering all 1279 earthquakes together. Thus, we have found no evidence that would suggest that monitoring changes in GIM‐TEC might be useful for predicting earthquakes.