On the planetary ionospheric vortex electric fields

The paper presents a physical mechanism of large-scale vortex electric field generation in the ionospheric E- and F-layers. It shows that the planetary-scale, synoptic short-period (from several second to several hours) and fast processes (with propagation velocity higher than 1 km/s) produce a planetary-scale internal vortex electric field. Its value may far exceed that of the dynamo-field generated in the same ionospheric layer by local wind motion. We found, that an ionospheric source of the vortex electric field is spatial inhomogeneity of the geomagnetic field.     

Application of HF Doppler measurements for the investigation of internal atmospheric waves in the ionosphere

We present results of the spectral analysis of data series of Doppler frequency shifted signals reflected from the ionosphere, using experimental data received at Kazan University, Russia. Spectra of variations with periods from 1 min to 60 days have been calculated and analyzed for different scales of periods. The power spectral density for spring and winter differs by a factor of 3–4. Local maxima of variation amplitude are detected, which are statistically significant. The periods of these amplitude increases range from 6 to 12 min for winter, and from 24 to 48 min for autumn. Properties of spectra for variations with the periods of 1–72 h have been analyzed. The maximum of variation intensity for all seasons and frequencies corresponds to the period of 24 h. Spectra of variations with periods from 3 to 60 days have been calculated. The maxima periods of power spectral density have been detected by the MUSIC method for the high spectral resolution. The detected periods correspond to planetary wave periods. Analysis of spectra for days with different level of geomagnetic activity shows that the intensity of variations for days with a high level of geomagnetic activity is higher.     

Effects of X2-class solar flare events on ionospheric GPS-TEC and radio waves over Brazilian sector

In this investigation, we present and discuss the effects of 6 X2-class solar flare events in the ionospheric F region over Brazilian sector that occurred during 2013 to 2015. For this investigation, we present the vertical total electron content (VTEC) observations from nearly 120 Global Positioning System (GPS) receivers all over the Brazilian sector for each event. Also, ionospheric sounding observations obtained in São José dos Campos (23.2°S, 45.9°W, dip latitude 17.6°S; hereafter referred to as SJC), under the southern crest of the equatorial ionospheric anomaly (EIA), Brazil, are presented. The observations show that the greatest TEC impact occurs with the EUV fluxes increases lasting for more than one hour and when the solar active region is located close to the solar disc center. We present a detailed study of the efficiency of the EUV flux with wavelengths ranging from 0.1 to 190?nm for the F region ionization. The largest increase of ΔTEC occurs below the magnetic equator line, covering mainly the central, northeast, southeast and south regions, which includes the equatorial ionospheric anomaly (EIA) region. The ionograms show partial or total fade out in the echoes traces observed causing blackouts of radio signals of up to 60?min, which can have serious consequences to technological systems of public and private agencies around Brazilian sector. This study can help to better understand the effects of solar flares in the ionospheric F region.     

Doppler utilised Kalman estimation (DUKE) of ionospheric delay for satellite navigation

The ionospheric delay experienced by the satellite navigation signals depends upon the Total Electron Content (TEC) and needs to be corrected. While the single frequency receivers always use parametric models to correct this delay, dual frequency receivers, when suffers a loss of lock of one of its signal, also has to resort to these models. Here, an alternative method, based on Doppler, surrogated by range rate variation, has been attempted to estimate the ionospheric delay using a Kalman filter. GPS data have been used for all visible satellites over four days selected around the equinox and solstice with nominal geomagnetic conditions and estimations done in continuous and calibrated modes. Results of continuous estimation, obtained for a mid latitude station, showed moderate accuracy while it was significantly better for the calibrated mode with no seasonal dependence. Estimations done for station within the extent of equatorial anomaly, has not only resulted in relative deterioration in performance, but also shown seasonal dependence. Compared with estimates of Klobuchar model, the Calibrated estimation showed superior performance, conspicuously in the mid latitude station. However, for the continuous mode, performance was at par with the model at higher latitudes but inferior to it in regions within the extent of the equatorial anomaly.     

Comparing ionospheric models with mid-latitude ionosonde observations

The purpose of this research work is to validate the ionospheric models (IRI and CHIU) to assess its suitability and usefulness as an operational tool. The ionospheric model is a computer model designed to predict the state of the global ionosphere for 24 h. The scope was limited to conduct comparisons between the predicted F₂ layer critical frequencies (f₀F₂) against observed ionosonde data. The ionospheric prediction model (IPM) was designed to predict by using monthly median sunspot number, while the observation data are taken from two digital ionospheric sounding stations (Okinawa, 26.28N, 127.8E and Wakkanai, 45.38N, 141.66E) which lies within the mid-latitude region of the globe. Analysis of the f₀F₂ data from stations for year (2001) with high solar activity and year (2004) with low solar activity, four months (March, June, September and December) chosen based primarily on data availability. From results it seen that the ratio between monthly median predicted and observed f₀F₂ values for each model used in this research work and for the chosen months was nonlinear with local time, so the empirical formula for applying correction factors were determined, these formula can be used to correct the error occurred in predicted f₀F₂ value.     

Comparison of model and experimental ionospheric parameters at high latitudes

Modern use and study of the auroral region needs to attract a wider class of models for describing conditions of radio wave propagation in the ionosphere. In this paper the possibilities of the International Reference Ionosphere model, well-proven and widespread in the mid-latitudes, are investigated in the high latitude zone. Model and measured values of the critical frequency foF2 for two mid-latitude stations (Juliusruh and Goosebay) and four high-latitude ones (Loparskaya, Sodankyla, Sondrestrom, Thule) are compared. Deviations of medians, variations from day to day and solar activity trends seemed to be similar for both areas. This similarity is irrespective of the RZ12 index. Special attention is paid to the TEC parameter and its determination using 6 versions of models, a new version of the model IRI2010 (IRI-Plas) among them. It is shown that the IRI-Plas model significantly improves the definition of TEC in contrast to the versions of IRI2007 and the new model NeQuick. The use of the median of the experimental equivalent slab thickness, together with the current values of the TEC, increases by a factor of two the agreement between calculated and measured foF2 values as compared with the variations from day to day. This allows foF2 to be defined in near-real time.     

Ionospheric oscillations caused by geomagnetic Pi2 pulsations and their observations by multipoint continuous Doppler sounding; first results

In 2004, we started operating a continuous Doppler sounding system to investigate ionospheric signatures of infrasonic, short period acoustic gravity waves and geomagnetic pulsations. Since January 2007, four stable 3.59 MHz transmitters have been in operation in the western part of the Czech Republic. Multipoint measurements enable us to investigate horizontal propagation of waves and disturbances in the ionosphere and to estimate horizontal distances over which these waves (disturbances) are correlated.     

Variations in the ionospheric scale height parameter at the F2 peak over Grahamstown,South Africa

The Grahamstown, South Africa (33.3°S, 26.5°E) ionospheric field station operates a Lowell digital pulse ionospheric sounder (Digisonde) whose output includes scaled parameters derived from the measured ionogram. One of these output parameters is the ionospheric scale height parameter (H), and this paper presents an analysis of the seasonal, diurnal, and solar activity variations of this parameter over the Grahamstown station. Ionosonde data from three years 2002, 2003, and 2004 were used in this study. The data was subjected to a general trend analysis to remove any outliers and then the monthly median data were used to explore the different variations. The results of this analysis were found to be similar to what has already been presented in the literature for low latitude stations, and are presented as well as the correlation at this mid-latitude station between the H parameter, the IRI shape parameter (B0), and the peak electron density (NmF2).     

Atmosphere-ionosphere coupling from convectively generated gravity waves

Ionospheric variability impacts operational performances of a variety of technological systems, such as HF communication, Global Positioning System (GPS) navigation, and radar surveillance. The ionosphere is not only perturbed by geomagnetic inputs but is also influenced by atmospheric tides and other wave disturbances propagating from the troposphere to high altitudes. Atmospheric Gravity Waves (AGWs) excited by meteorological sources are one of the largest sources of mesoscale variability in the ionosphere. In this paper, Total Electron Content (TEC) data from networks of GPS receivers in the United States are analyzed to investigate AGWs in the ionosphere generated by convective thunderstorms. Two case studies of convectively generated gravity waves are presented. On April 4, 2014 two distinct large convective systems in Texas and Arkansas generated two sets of concentric AGWs that were observed in the ionosphere as Traveling Ionospheric Disturbances (TIDs). The period of the observed TIDs was 20.8 min, the horizontal wavelength was 182.4 km, and the horizontal phase speed was 146.4 m/s. The second case study shows TIDs generated from an extended squall line on December 23, 2015 stretching from the Gulf of Mexico to the Great Lakes in North America. Unlike the concentric wave features seen in the first case study, the extended squall line generated TIDs, which exhibited almost plane-parallel phase fronts. The TID period was 20.1 min, its horizontal wavelength was 209.6 km, and the horizontal phase speed was 180.1 m/s. The AGWs generated by both of these meteorological events have large vertical wavelength (>100 km), which are larger than the F2 layer thickness, thus allowing them to be discernible in the TEC dataset.     

A prediction model of short-term ionospheric foF2 based on AdaBoost

In this paper, the AdaBoost-BP algorithm is used to construct a new model to predict the critical frequency of the ionospheric F2-layer (foF2) one hour ahead. Different indices were used to characterize ionospheric diurnal and seasonal variations and their dependence on solar and geomagnetic activity. These indices, together with the current observed foF2 value, were input into the prediction model and the foF2 value at one hour ahead was output. We analyzed twenty-two years’ foF2 data from nine ionosonde stations in the East-Asian sector in this work. The first eleven years’ data were used as a training dataset and the second eleven years’ data were used as a testing dataset. The results show that the performance of AdaBoost-BP is better than those of BP Neural Network (BPNN), Support Vector Regression (SVR) and the IRI model. For example, the AdaBoost-BP prediction absolute error of foF2 at Irkutsk station (a middle latitude station) is 0.32 MHz, which is better than 0.34 MHz from BPNN, 0.35 MHz from SVR and also significantly outperforms the IRI model whose absolute error is 0.64 MHz. Meanwhile, AdaBoost-BP prediction absolute error at Taipei station from the low latitude is 0.78 MHz, which is better than 0.81 MHz from BPNN, 0.81 MHz from SVR and 1.37 MHz from the IRI model. Finally, the variety characteristics of the AdaBoost-BP prediction error along with seasonal variation, solar activity and latitude variation were also discussed in the paper.     

Epochwise prediction of GPS single differenced ionospheric delays of formation flying spacecraft

Both single and dual frequency GPS relative navigation filters may benefit from proper predictions of single differenced ionospheric delays. In this article, the single differenced ionospheric delays of GPS observations are predicted for the GRACE formation during the switch manoeuvre.Two prediction methods are considered. The first is based on a Taylor expansion to first order of a mapping function that maps slant total electron content measurements to vertical total electron content estimates. The second method fits a shape profile through undifferenced ionospheric data available. It then raytraces through this profile to estimate the difference in total electron content along the path of the GPS signals.Continuously changing ionospheric conditions hamper the assessment of the quality of the predictions. Comparison of both methods shows that the raytracing method performs better. The difference of predictions and measurements generally shows a smaller RMS than the measurements alone. However, both methods suffer from a number of systematically unpredicted observations, which arise from small unaccounted differential variations in electron densities along the path of the GPS signals. These prediction methods perform better when spacecraft separation is small. Baselines considered here range from tens of kilometres down to several hundred metres. When smallest spacecraft separation occurs (0.4 km), the single differenced ionospheric delays exhibit RMS values of 0.0089 m. The first method shows a difference between measurements and predictions with an RMS of 0.0081 m. For the second method the difference RMS is found to be 0.0067 m.     

Using neural networks to determine the optimum solar input for the prediction of ionospheric parameters

Neural networks (NNs) are proving to be ideal tools for modeling the behaviour of the ionosphere. The NNs are trained using a database of archived data describing the relationship between the output parameter and an input space. The input space is designed from knowledge of those variables that affect the behaviour of the output parameter. For ionospheric parameters this input space would always include a solar variable due to the strong influence that the sun has on ionospheric behaviour.     

Variability of the ionospheric electron density at fixed heights and validation of IRI-2007 profile’s prediction at Ilorin

A study on the variability of the equatorial ionospheric electron density was carried out at fixed heights below the F2 peak using one month data for each of high and low solar activity periods. The data used for this study were obtained from ionograms recorded at Ilorin, Nigeria, and the study covers height range from 100 km to the peak of the F2 layer for the daytime hours and height range from 200 km to the peak of the F2 layer for the nighttime hours. The results showed that the deviation of the electron density variation from simple Chapman variation begins from an altitude of about 200 km for the two months investigated. Daytime minimum variability of between 2.7% and 9.0% was observed at the height range of about 160 and 200 km during low solar activity (January 2006) and between 3.7% and 7.8% at the height range of 210 and 260 km during high solar activity (January 2002). The nighttime maximum variability was observed at the height range of 210 and 240 km at low solar activity and at the height range of 200 and 240 km at high solar activity. A validation of IRI-2007 model electron density profile’s prediction was also carried out. The results showed that B0 option gives a better prediction around the noontime.     

Effects observed in the equatorial and low latitude ionospheric F-region in the Brazilian sector during low solar activity geomagnetic storms and comparison with the COSMIC measurements

The main objective of the present investigation has been to compare the ionospheric parameters (NmF2 and hmF2) observed by two ground-based ionospheric sounders (one at PALMAS- located near the magnetic equator and the other at Sao Jose dos Campos-located in the low-latitude region) in the Brazilian sector with that by the satellite FORMOSAT-3/COSMIC radio occultation (RO) measurements during two geomagnetic storms which occurred in December 2006 and July 2009. It should be pointed out that in spite of increasing the latitude (to 10°) and longitude (to 20°) around the stations; we had very few common observations. It has been observed that both the peak electron density (NmF2) and peak height (hmF2) observed by two different techniques (space-borne COSMIC and ground-based ionosondes) during both the geomagnetic storm events compares fairly well (with high correlation coefficients) at the two stations in the Brazilian sector. It should be pointed out that due to equatorial spread F (ESF) in the first storm (December 2006) and no-reflections from the ionosphere during nighttime in the second storm (July 2009), we had virtually daytime data from the two ionosondes.     

Observations of ionospheric irregularities and its correspondence with sporadic e occurrence over South Korea and Japan

This study analyzed the occurrence of ionospheric irregularities over South Korea and Japan (mid-latitudes) during the years 2010–2015. The irregularities were quantified using the rate of change of total electron content (TEC) index (ROTI), which detects irregularities with scale sizes in the range of 400 m–2.5 km. The ROTI threshold for an ionospheric irregularity to have occurred was set as 0.1 TECU/min. Results showed that ionospheric irregularities mostly occur during night-time and around local noon-time in the seasons of spring and summer. In addition, the percentage of ionospheric irregularities had a high positive correlation with solar flux (F10.7) (r > 0.72). For the first time, we showed good correspondence between ionospheric irregularities measured by the ROTI index and sporadic E (Es). The median ROTI associated with ionospheric irregularities over a South Korea station (DAEJ) and a Japan station (KGNI) were 0.131 and 0.125 TECU/min, respectively. However, in severe cases of ionospheric irregularities, the ROTI values over DAEJ (KGNI) can reach 0.246 (0.217) and 0.314 (0.339) TECU/min during day and night, respectively. These critical ROTI values can be important in interpreting and monitoring ionospheric irregularity occurrence over South Korea and Japan.     

Diurnal and seasonal variation of F2-layer ionospheric parameters at equatorial ionization anomaly crest region and their comparison with IRI-2001

Diurnal and seasonal variations of critical frequency of ionospheric F2-region ‘foF2’ and the height of peak density ‘hmF2’ are studied using modern digital ionosonde observations of equatorial ionization anomaly (EIA) crest region, Bhopal (23.2°N, 77.6°E, dip 18.5°N), during solar minimum period 2007. Median values of these parameters are obtained at each hour using manually scaled data during different seasons and compared with the International Reference Ionosphere-2001 model predictions. The observations suggest that on seasonal basis, the highest values of foF2 are observed during equinox months, whereas highest values of hmF2 are obtained in summer and lowest values of both foF2 and hmF2 are observed during winter. The observed median and IRI predicted values of foF2 and hmF2 are analyzed with upper and lower bound of inter-quartile range (IQR) and it is find out that the observed median values are well inside the inter-quartile range during the period of 2007. Comparison of the recorded foF2 and hmF2 values with the IRI-2001 output reveals that IRI predicted values exhibit better agreement with hmF2 as compared to foF2. In general, the IRI model predictions show some agreement with the observations during the year 2007. Therefore it is still necessary to implement improvements in order to obtain better predictions for EIA regions.     

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.     

Coordinated in situ measurements of plasma irregularities and ground based scintillation observations at the crest of equatorial anomaly

First comparison of in situ density fluctuations measured by the DEMETER satellite with ground based GPS receiver measurements at the equatorial anomaly station Bhopal (geographic coordinates (23.2°N, 77.6°E); geomagnetic coordinates (14.29°N, 151.12°E)) for the low solar activity year 2005, are presented in this paper. Calculation of the diurnal maximum of the strength of the equatorial electrojet, which can serve as precursor to ionospheric scintillations in the anomaly region is also done. The Langmuir Probe experiment and Plasma Analyzer onboard DEMETER measure the electron and ion densities respectively. Irregularities in electron density distribution cause scintillations on transionospheric links and there exists a close relationship between an irregularity and scintillation. In 40% of the cases, DEMETER detects the irregularity structures (dNe/Ne ? 5% and dNi/Ni (O⁺) ? 5%) and GPS L band scintillations (S4 ? 0.2) are also observed around the same time, for the low solar activity period. It is found that maximum irregularity intensity is obtained in the geomagnetic latitude range of 10–20° for both electron density and ion density. As the GPS signals pass through this irregularity structure, scintillations are recorded by the GPS receiver installed at the equatorial anomaly station, Bhopal it is interesting to note that in situ density fluctuations observed on magnetic flux tubes that pass over Bhopal can be used as indicator of ionospheric scintillations at that site. Many cases of density fluctuations and associated scintillations have been observed during the descending low solar activity period. The percentage occurrence of density irregularities and scintillations shows good correspondence with diurnal maximum of the strength of electrojet, however this varies with different seasons with maximum correspondence in summer (up to 66%) followed by equinox (up to 50%) and winter (up to 46%). Also, there is a threshold value of EEJ strength to produce density irregularities ((dNe/Ne)_max ? 5%) and for moderate to strong scintillations (S4 ? 0.3) to occur. For winter this value is found to be ∼40 nT whereas for equinox and summer it is around 50 nT.     

Magnetospheric and ionospheric manifestations of interplanetary foreshocks

The geometry of a typical interplanetary shock front in the vicinity of the Earth’s orbit predicts that the leading edge of the foreshock region comes into contact with the magnetosphere a few hours ahead of geomagnetic sudden impulses (SI). There is reason to believe that the interaction of the magnetosphere with the foreshock leads to magnetic and ionospheric disturbances, which can be detected by ground-based instruments. We searched for specific precursors of SIs in data from the Scandinavian riometer network and in the short period geomagnetic pulsation data from mid-latitude magnetometers. We found that SIs were preceded by the following three features: (1) an increase in riometric absorption, (2) excitation of Pcl magnetic pulsations and (3) a spectral broadening of the Pc3 magnetic pulsations. Our observations may be useful for the study of acceleration processes in the solar wind. These observations are also of potential forecasting interest.     

Comparison of peak characteristics of F2 ionospheric layer over Tehran region at a low solar activity period with IRI-2001 and IRI-2007 models predictions

In the present work values of peak electron density (NmF2) and height of F2 ionospheric layer (hmF2) over Tehran region at a low solar activity period are compared with the predictions of the International Reference Ionosphere models (IRI-2001 and IRI-2007). Data measured by a digital ionosonde at the ionospheric station of the Institute of Geophysics, University of Tehran from July 2006 to June 2007 are used to perform the calculations. Formulations proposed by  and  are utilized to calculate the hmF2. The International Union of Radio Science (URSI) and International Radio Consultative Committee (CCIR) options are employed to run the IRI-2001 and IRI-2007 models. Results show that both IRI-2007 and IRI-2001 can successfully predict the NmF2 and hmF2 over Tehran region. In addition, the study shows that predictions of IRI-2007 model with CCIR coefficient has closer values to the observations. Furthermore, it is found that the monthly average of the percentage deviation between the IRI models predictions and the values of hmF2 and NmF2 parameters are less than 10% and 21%, respectively.