Comparison of the CHAMP radio occultation data with the Canadian advanced digital ionosonde in the Polar Regions

In this paper we present the results of the comparison of the retrieved electron density profiles of the Ionospheric Radio Occultation (IRO) experiment on board CHAMP (CHAllenging Minisatellite Payload), with the ground ionosonde profiles for the Polar Regions. IRO retrieved electron density profiles from CHAMP are compared with Canadian Advanced Digital Ionosonde (CADI) measurements at two vertical sounding stations well within the Polar Cap, Eureka (geog. 80°13′ N; 86°11′ W) and Resolute Bay (geog. 74°41′ N; 94°54′ W). We compared the ionospheric parameters such as the peak electron density of the F-layer (NmF₂) and the peak height of the F-layer (h_mF₂) for a 3-year period, 2004–2006. CHAMP derived NmF₂ shows reasonable agreement with the ionosonde retrieved NmF₂ for both the stations (0.76 and 0.71 correlation coefficient, for Eureka and Resolute Bay, respectively) whereas the h_mF₂ agreement is not that acceptable (0.25 and 0.37 correlation coefficient, respectively). The h_mF₂ from vertical sounding showed less spread than the CHAMP h_mF₂.     

An alternative Shell inversion technique - Analysis and validation based on COSMIC and ionosonde data

Multi-channel Global Positioning System (GPS) carrier phase signals, received by the six low Earth orbiting (LEO) satellites from the Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC) program, were used to undertake active limb sounding of the Earth’s atmosphere and ionosphere via radio occultation. In the ionospheric radio occultation (IRO) data processing, the standard Shell inversion technique (SIT), transformed from the traditional Abel inversion technique (AIT), is widely used, and can retrieve good electron density profiles. In this paper, an alternative SIT method is proposed. The comparison between different inversion techniques will be discussed, taking advantage of the availability of COSMIC datasets. Moreover, the occultation results obtained from the SIT and alternative SIT at 500 km and 800 km, are compared with ionosonde measurements. The electron densities from the alternative SIT show excellent consistency to those from the SIT, with strong correlations over 0.996 and 0.999 at altitudes of 500 km and 800 km, respectively, and the peak electron densities (N_mF₂) from the alternative SIT are equivalent to the SIT, with 0.839 vs. 0.844, and 0.907 vs. 0.909 correlation coefficients when comparing to those by the ionosondes. These results show that: (1) the N_mF₂ and h_mF₂ retrieved from the SIT and alternative SIT are highly consistent, and in a good agreement with those measured by ionosondes, (2) no matter which inversion technique is used, the occultation results at the higher orbits (∼800 km) are better than those at the lower orbits (∼500 km).     

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.     

Ionospheric responses to the 21 August 2017 great American solar eclipse – A multi-instrument study

Herein, we report on the ionospheric responses to a total solar eclipse that occurred on 21 August 2017 over the US region. Ground-based GPS total electron content (TEC) data along with ground-based measurements (Millstone Hill Observatory (MHO) and digital ionosondes) and space-based measurements (COSMIC radio occultation (RO) technique) allowed us to identify eclipse-associated ionospheric responses. TEC data at ~20°, ~30°, and ~40°N latitudes from the west to east longitudes show not only considerable depression but also wave-like characteristics in TEC both in the path of totality and away from it, exclusively on the day of eclipse. Interestingly, the observed depressions are associated with lesser (higher) magnitudes at stations over which the solar obscuration percentage was meager (significant), a clear indication of bow-wave-like features. The MHO observes a 30% reduction in F₂-layer electron densities between 180 and 220 km on eclipse day. Ionosonde-scaled parameters over Boulder (40.4°N, 100°E) and Austin (30.4°N, 94.4°E) show a significant decrease in critical frequencies while an altitude elevation is seen in the virtual heights of the F-layer only during the eclipse day and that decreases are associated with wave-like signatures, which could be attributed to eclipse-generated waves. The estimated vertical electron density profile from the COSMIC RO-based technique shows a maximum depletion of 40%. Relatively intense and moderate depths of TEC depression, considerable reductions in the F₂-layer electron densities measured by the MHO and COSMIC RO-measured densities at the F₂-layer peak, and elevations in virtual heights and reduction in the critical frequencies measured by ionosondes during the eclipse day could be due to the eclipse-induced dynamical effects such as gravity waves (GWs) and their associated electro-dynamical effects (modification of ionospheric electric fields due to GWs).     

Large-scale traveling ionospheric disturbances of auroral origin according to the data of the GPS network and ionosondes

The intensity of large-scale traveling ionospheric disturbances (LS TIDs), registered using measurements of total electron content (TEC) during the magnetic storms on October 29–31, 2003, and on November 7–11, 2004, had been compared with that of local electron density disturbances. The data of TEC measurements at ground-based GPS receivers located near the ionospheric stations and the corresponding values of the critical frequency of the ionospheric F region f_oF2 were used for this purpose. The variations of TEC and f_oF2 were similar for all events mentioned above. The previous assumption that the ionospheric region with vertical extension from 150 to 200 km located near the F-layer maximum mainly contributes to the TEC variations was confirmed for the cases when the electron density disturbance at the F region maximum was not more than 50%. However, this region probably becomes vertically more extended when the electron density disturbance in the ionospheric F region is about 85%.     

Comparison of subpeak electron density profiles deduced from ionograms with the International Reference Ionosphere (IRI)

A comparison is made between the subpeak electron density profiles, obtained at selected local hours by vertical ionospheric sounding at the ionospheric station at Sofia (42.6°N; 23.3°E) and the IRI profiles for spring, summer, winter and two levels of solar activity (R = 10 and 100). It is demonstrated that the ionospheric profiles above Sofia are in rather good agreement with the values computed with IRI.     

Testing the three-dimensional IRI-SIRMUP-P mapping of the ionosphere for disturbed periods

This paper describes the three-dimensional (3-D) electron density mapping of the ionosphere given as output by the assimilative IRI-SIRMUP-P (ISP) model for three different geomagnetic storms. Results of the 3-D model are shown by comparing the electron density profiles given by the model with the ones measured at two testing ionospheric stations: Roquetes (40.8°N, 0.5°E), Spain, and San Vito (40.6°N, 17.8°E), Italy. The reference ionospheric stations from which the autoscaled foF2 and M(3000)F2 data as well as the real-time vertical electron density profiles are assimilated by the ISP model are those of El Arenosillo (37.1°N, 353.3°E), Spain, Rome (41.8°N, 12.5°E), and Gibilmanna (37.9°N, 14.0°E), Italy. Overall, the representation of the ionosphere made by the ISP model is better than the climatological representation made by only the IRI-URSI and the IRI-CCIR models. However, there are few cases for which the assimilation of the autoscaled data from the reference stations causes either a strong underestimation or a strong overestimation of the real conditions of the ionosphere, which is in these cases better represented by only the IRI-URSI model. This ISP misrepresentation is mainly due to the fact that the reference ionospheric stations covering the region mapped by the model turn out to be few, especially for disturbed periods when the ionosphere is very variable both in time and in space and hence a larger number of stations would be required. The inclusion of new additional reference ionospheric stations could surely smooth out this concern.     

Comparing topside and bottomside-measured characteristics of the F2 layer peak

The ionospheric characteristics of the F2 layer peak have been measured with ionosondes from the ground or with satellites from space. The most common characteristics are the F2-peak density NmF2 and peak height hmF2. In addition to these two parameters this paper studies the F2-peak scale height. Comparing the median values of hmF2 and NmF2 obtained from topside and bottomside sounding shows good agreement in general. The Chapman scale height values for the F2 layer peak derived from topside profiles, H_m,top, are generally several times larger than H_m,bot derived from bottomside profiles.     

A numerical study of noontime bite-outs in F2-region electron density

Noontime bite-outs in ionospheric F₂-region electron density in the geomagnetic equatorial, low, and middle latitudes have been reproduced in this study using the Thermosphere Ionosphere Electrodynamic General Circulation Model (TIEGCM). The different driving mechanisms of F₂-region noontime bite-outs at different latitudes were further verified by modeling simulation. In the geomagnetic equator, the upward ExB plasma drifts are the main process to produce the noontime bite-outs in ionospheric electron density. In the geomagnetic low latitudes, both the electrical fields and poleward meridional winds play a crucial role in forming the noontime bite-outs. In contrast to the case at the geomagnetic equator, a weaker fountain effect might be an favorable condition for the noontime bite-outs to occur at low latitudes. For geomagnetic middle latitudes, an upward plasma flux and poleward meridional winds are the dominated drivers in producing the noontime bite-outs. Modeling results show that a large upward plasma flux and poleward meridional winds make the noontime bite-outs to occur and observable over middle latitudes.     

Study of maximum electron density NmF2 at Karachi and Islamabad during solar minimum (1996) and solar maximum (2000) and its comparison with IRI

The monthly hourly medians of maximum electron density, N_mF₂, at two Pakistani ionospheric stations, Karachi and Islamabad, have been determined for solar minimum (1996) and solar maximum (2000) and compared with IRI predictions using the URSI coefficients. At night and pre-noon period the N_mF₂ values at both stations are almost equal during the 2 years. However, at post-noon the values at Karachi are considerably larger than those at Islamabad due to the equatorial or geomagnetic anomaly. Karachi (geomag. coord. 16.44°N, 139.08°E) lies near the region of the equatorial anomaly (+20 and −20 geomagnetic latitude), so most of the N_mF₂ values at Karachi are larger than those at Islamabad (geomag. coord. 24.46°N, 145.67°E). The maximum monthly values of N_mF₂ show a semi-annual variation at Karachi and Islamabad both during 1996 and 2000 as predicted by IRI.     

A realistic simulation framework to evaluate ionospheric tomography

Observation of the 3-dimensional (3-D) electron density of the ionosphere is useful to study large-scale physical processes in space weather events. Ionospheric data assimilation and ionospheric tomography are methods that can create an image of the 3-D electron density distribution. While multiple techniques have been developed over the past 30 years, there are relatively few studies that show the accuracy of the algorithms. This paper outlines a novel simulation approach to test the quality of an ionospheric tomographic inversion. The approach uses observations from incoherent scatter radar (ISR) scans and extrapolates them spatially to create a realistic ionospheric representation. A set of total electron content (TEC) measurements can then be simulated using real geometries from satellites and ground receivers. This data set, for which the ‘truth’ ionosphere is known, is used as input for a tomographic inversion algorithm to estimate the spatial distribution of electron density. The reconstructed ionospheric maps are compared with the truth ionosphere to calculate the difference between the images and the truth.To demonstrate the effectiveness of this simulation framework, an inversion algorithm called MIDAS (Multi-Instrument Data Analysis Software) is evaluated for three geographic regions with differing receiver networks. The results show the importance of the distribution and density of GPS receivers and the use of a realistic prior conditioning of the vertical electron density profile. This paper demonstrates that when these requirements are met, MIDAS can reliably estimate the ionospheric electron density. When the region under study is well covered by GPS receivers, as in mainland Europe or North America, the errors in vertical total electron content (vTEC) are smaller than 1 TECu (2–4%) . In regions with fewer and more sparsely distributed receivers, the errors can be as high as 20–40%. This is caused by poor data coverage and poor spatial resolution of the reconstruction, which has an important effect on the calibration process of the algorithm.     

Ionospheric storm of September 2017 observed at ionospheric station Pruhonice,the Czech Republic

The ionospheric effects induced by the September 2017 storm have been exceptional compared to other events in the solar cycle 24. This paper gathers results of the ionospheric observations at the European middle latitude station Pruhonice. It consists of evaluation of ionospheric vertical and oblique sounding, Digisonde drift measurement, and data obtained from the Continuous Doppler Sounding System. We observed strong ionospheric response with an unusual stratification of ionospheric layers, Large Scale Traveling ionospheric disturbances, changes in electron density, and increase and oscillations in plasma drift velocity.     

Estimation of ground range on the sweep frequency backscatter leading edge

The high frequency management system with backscatter radar supplies the real time ionosphere channel conditions to high frequency users, which leads to the demand for the ground range between the radar location and the scatters on the distant ground. The ionosphere electron density profile is usually inversed to obtain the ground range. An inversion algorithm, with which the ground range on the leading edge of the backscatter ionograms can be obtained without electron density, is presented in this paper. The ray path geometry of the backscatter sounding and the change in the group path on the leading edge with operating frequency are used to derive the ground range. Synthesized backscatter ionogram and experimental backscatter ionograms are processed to validate the algorithm. The results indicate that the algorithm is usable for high frequency management system.     

Planetary wave coupling of the atmosphere–ionosphere system during the Northern winter of 2008/2009

This paper presents the global spatial (latitude and altitude) structure and temporal variability of the ∼23-day ionospheric zonally symmetric (s = 0) planetary wave (PW) seen in the Northern winter of 2008/2009 (October 2008–March 2009). It is shown that these ∼23-day ionospheric oscillations are forced from PWs propagating from below. The COSMIC ionospheric parameters f_oF2 and h_mF2 and electron density at fixed altitudes and the SABER temperatures were utilized in order to define the waves which are present simultaneously in the atmosphere and ionosphere. The long-period PWs from the two data sets have been extracted through the same data analysis method. The similarity between the lower thermospheric ∼23-day (s = 0) temperature PW and its ionospheric electron density response provides valuable and strong experimental evidence for confirming the paradigm of atmosphere–ionosphere coupling.     

An overview of radar soundings of the martian ionosphere from the Mars Express spacecraft

The Mars Express spacecraft carries a low-frequency radar called MARSIS (Mars Advanced Radar for Subsurface and Ionosphere Sounding) that is designed to study the subsurface and ionosphere of Mars. In this paper, we give an overview of the ionospheric sounding results after approximately one year of operation in orbit around Mars. Several types of ionospheric echoes are commonly observed. These include vertical echoes caused by specular reflection from the horizontally stratified ionosphere; echoes from a second layer in the topside ionosphere, possibly associated with O⁺ ions; oblique echoes from upward bulges in the ionosphere; and a variety of other echoes that are poorly understood. The vertical echoes provide electron density profiles that are in reasonable agreement with the Chapman photo-equilibrium model of planetary ionospheres. On the dayside of Mars the maximum electron density is approximately 2 × 10⁵ cm⁻³. On the nightside the echoes are often very diffuse and highly irregular, with maximum electron densities less than 10⁴ cm⁻³. Surface reflections are sometimes observed in the same frequency range as the diffuse echoes, suggesting that small isolated holes exist in the nightside ionosphere, possibly similar to those that occur on the nightside of Venus. The oblique echoes arise from upward bulges in the ionosphere in regions where the crustal magnetic field of Mars is strong and nearly vertical. The bulges tend to be elongated in the horizontal direction and located in regions between oppositely directed arch-like structures in the crustal magnetic field. The nearly vertical magnetic field lines in the region between the arches are thought to connect into the solar wind, thereby allowing solar wind electrons to heat the lower levels of the ionosphere, with an attendant increase in the scale height and electron density.     

An empirical model of electron density in low latitude at 600 km obtained by Hinotori satellite

An empirical model of electron density (Ne) was constructed by using the data obtained with an impedance probe on board Japanese Hinotori satellite. The satellite was in circular orbit of the height of 600 km with the inclination of 31 degrees from February 1981 to June 1982. The constructed model gives Ne at any local time with the time resolution of 90 min and between −25 and 25 degrees in magnetic latitude with its resolution of 5 degrees in the range of F_10.7 from 150 to 250 under the condition of K_p < 4. Spline interpolations are applied to the functions of day of year, geomagnetic latitude and solar local time, and linear interpolation is applied to the function of F_10.7. Longitude dependence of Ne is not taken into account. Our density model can reproduce solar local time variation of electron density at 600 km altitude better than current International Reference Ionosphere (IRI2001) model which overestimates Ne in night time and underestimates Ne in day time. Our density model together with electron temperature model which has been constructed before will enable more understanding of upper ionospheric phenomenon in the equatorial region.     

Longitudinal behaviors of the IRI-B parameters of the equatorial electron density profiles retrieved from FORMOSAT-3/COSMIC radio occultation measurements

The electron density profiles in the bottomside F₂-layer ionosphere are described by the thickness parameter B0 and the shape parameter B1 in the International Reference Ionosphere (IRI) model. We collected the ionospheric electron density (N_e) profiles from the FORMOSAT-3/COSMIC (F3/C) radio occultation measurements from DoY (day number of year) 194, 2006 to DoY 293, 2008 to investigate the daytime behaviors of IRI-B parameters (B0 and B1) in the equatorial regions. Our fittings confirm that the IRI bottomside profile function can well describe the averaged profiles in the bottomside ionosphere. Analysis of the equatorial electron density profile datasets provides unprecedented detail of the behaviors of B0 and B1 parameters in equatorial regions at low solar activity. The longitudinal averaged B1 has values comparable with IRI-2007 while it shows little seasonal variation. In contrast, the observed B0 presents semiannual variation with maxima in solstice months and minima in equinox months, which is not reproduced by IRI-2007. Moreover, there are complicated longitudinal variations of B0 with patterns varying with seasons. Peaks are distinct in the wave-like longitudinal structure of B0 in equinox months. An outstanding feature is that a stable peak appears around 100°E in four seasons. The significant longitudinal variation of B0 provides challenges for further improving the presentations of the bottomside ionosphere in IRI.     

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).     

The global response of terrestrial ionosphere to the December 2015 space weather event

This paper investigates the ionospheric storm of December 19–21, 2015, which was initiated by two successive CME eruptions that caused a G3 space weather event. We used the in situ electron density (N_e) and electron temperature (T_e) and the Total Electron Content (TEC) measurements from SWARM-A satellite, as well as the O/N₂ observations from TIMED/GUVI to study the ionospheric impact. The observations reveal the longitudinal and hemispherical differences in the ionospheric response to the storm event. A positive ionospheric storm was observed over the American, African and Asian regions on 20 December, and the next day showed a negative storm. Both these exhibited hemispheric differences. A positive storm was observed over the East Pacific region on 21 December. It is seen that the net effect of both the disturbance dynamo electric field and composition differences become important in explaining the observed variability in topside ionospheric densities. In addition, we also discuss the T_e variations that occurred as a consequence of the space weather event.