The improvement of IRI profiles of O2+, NO+ and O+ by means of airglow measurements

Methods were suggested in /1,2,3/ to find the maximum electron density N_mF (or N_mF2) and the) corresponding height h_mF (h_mF2) in the ionosphere through a combination of measurements on two airglow lines: the UV oxygen line at λ135.6 nm and the red oxygen line near λ630 nm.     

A different method to update monthly median hmF2 values

The height, h_mF2, and the electron density, N_mF2, of the F2 peak are key model parameters to characterize the actual state of the ionosphere. These parameters, or alternatively the propagation factor, M3000F2, and the critical frequency, f_oF2, of the F2 peak, which are related to h_mF2 and N_mF2, are used to anchor the electron density vertical profile computed with different models such as the International Reference Ionosphere ( Bilitza, 2002), as well as for radio propagation forecast purposes. Long time series of these parameters only exist in an inhomogeneous distribution of points over the surface of Earth, where dedicated instruments (typically ionosondes) have been working for many years. A commonly used procedure for representing median values of the aforementioned parameters all over the globe is the one recommended by the ITU-R ( ITU-R, 1997). This procedure, known as the Jones and Gallet mapping technique, was based on ionosondes measurements gathered from 1954 to 1958 by a global network of around 150 ionospheric stations (  and ). Even though several decades have passed since the development of that innovative work, only few efforts have been dedicated to establish a new mapping technique for computing h_mF2 and N_mF2 median values at global scale or to improve the old method using the increased observational database. Therefore, in this work three different procedures to describe the daily and global behavior of the height of the F2 peak are presented. All of them represent a different and simplified method to estimate h_mF2 and are based on different mathematical expressions. The advantages and disadvantages of these three techniques are analyzed, leading to the conclusion that the recommended procedure to represent h_mF2 is best characterized by a Spherical Harmonics expansion of degree and order equal to 15, since the differences between the h_mF2 values obtained with the Jones and Gallet technique and those obtained using the abovementioned procedure are of only 1%.     

Some effects in the upper atmosphere during geomagnetic storms

This paper examines the response of the high latitude ionosphere–thermosphere system during two intense geomagnetic storms. For that, data taken by instruments on board Dynamic Explorer 2 taken at heights of the F2-layer are used. These results represent a comparison of simultaneous measurements of storm disturbances in gas composition, electron density and temperature in common local time sectors. Documented are an increase in electron temperature and a decrease in electron density; increases both in electron temperature and electron density; and the correlation between electron density decreases and increases in the ratio N₂/O. It is noticed that the decrease in electron density is sometimes due to an increase in the molecular nitrogen density N₂ and not always is attributed both to the increase in N₂ density and the simultaneous decrease in the atomic oxygen density.     

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.     

The variability and predictability of the IRI B0, B1 parameters over Grahamstown,South Africa

The International Reference Ionosphere (IRI) parameters B₀ and B₁ provide a representation of the thickness and shape, respectively, of the F2 layer of the bottomside ionosphere. These parameters can be derived from electron density profiles that are determined from vertical incidence ionograms. This paper aims to illustrate the variability of these parameters for a single mid latitude station and demonstrate the ability of the Neural Network (NN) modeling technique for developing a predictive model for these parameters. Grahamstown, South Africa (33.3°S, 26.5°E) was chosen as the mid latitude station used in this study and the B₀ and B₁ parameters for an 11 year period were determined from electron density profiles recorded at that station with a University of Massachusetts Lowell Center for Atmospheric Research (UMLCAR) Digisonde. A preliminary single station NN model was then developed using the Grahamstown data from 1996 to 2005 as a training database, and input parameters known to affect the behaviour of the F2 layer, such as day number, hour, solar and magnetic indices. An analysis of the diurnal, seasonal and solar variations of these parameters was undertaken for the years 2000, 2005 and 2006 using hourly monthly median values. Comparisons between the values derived from measured data and those predicted using the two available IRI-2001 methods (IRI tables and Gulyaeva, T. Progress in ionospheric informatics based on electron density profile analysis of ionograms. Adv. Space Res. 7(6), 39–48, 1987.) and the newly developed NN model are also shown in this paper. The preliminary NN model showed that it is feasible to use the NN technique to develop a prediction tool for the IRI thickness and shape parameters and first results from this model reveal that for the mid latitude location used in this study the NN model provides a more accurate prediction than the current IRI model options.     

Effect of vertical plasma transport on ionospheric F2-region parameters at equatorial latitude

The variability of the F2-layer even during magnetically quiet times are fairly complex owing to the effects of plasma transport. The vertical E × B drift velocities (estimated from simplified electron density continuity equation) were used to investigate the seasonal effects of the vertical ion drifts on the bottomside daytime ionospheric parameters over an equatorial latitude in West Africa, Ibadan, Nigeria (Geographic: 7.4°N, 3.9°E, dip angle: 6°S) using 1 year of ionsonde data during International Geophysical Year (IGY) of 1958, that correspond to a period of high solar activity for quiet conditions. The variation patterns between the changes of the vertical ion drifts and the ionospheric F2-layer parameters, especially; f_oF₂ and h_mF₂ are seen remarkable. On the other hand, we observed strong anti-correlation between vertical drift velocities and h′F in all the seasons. We found no clear trend between N_mF₂ and h_mF₂ variations. The yearly average value of upward daytime drift at 300 km altitude was a little less than the generally reported magnitude of 20 ms⁻¹ for equatorial F-region in published literature, and the largest upward velocity was roughly 32 ms⁻¹. Our results indicate that vertical plasma drifts; ionospheric F2-layer peak height, and the critical frequency of F2-layer appear to be somewhat interconnected.     

磁暴期间中纬度电离层剖面结构变化的数值模拟

利用电离层理论模型模拟了磁暴期间热层大气温度、成分、中性风和电场扰动对电离层电子密度剖面结构，特别是峰值密度和峰值高度变化的影响，结果表明，热层大气温度变化所引起光化反应系数的改变对电离层剖面结构影响不大；热层大气成分特别是Ｎ２／Ｏ的变化能有效地引起密度剖面变化，Ｎ２增加足以使峰值密度产生所观测到的负相暴；由中性风和电场引起等离子体漂移是峰值高度ｈｍＦ２变化的主要原因，但对电子密度的影响不足以抵消     

A new inversion algorithm for HF sky-wave backscatter ionograms

HF sky-wave backscatter sounding system is capable of measuring the large-scale, two-dimensional (2-D) distributions of ionospheric electron density. The leading edge (LE) of a backscatter ionogram (BSI) is widely used for ionospheric inversion since it is hardly affected by any factors other than ionospheric electron density. Traditional BSI inversion methods have failed to distinguish LEs associated with different ionospheric layers, and simply utilize the minimum group path of each operating frequency, which generally corresponds to the LE associated with the F₂ layer. Consequently, while the inversion results can provide accurate profiles of the F region below the F₂ peak, the diagnostics may not be so effective for other ionospheric layers. In order to resolve this issue, we present a new BSI inversion method using LEs associated with different layers, which can further improve the accuracy of electron density distribution, especially the profile of the ionospheric layers below the F₂ region. The efficiency of the algorithm is evaluated by computing the mean and the standard deviation of the differences between inverted parameter values and true values obtained from both vertical and oblique incidence sounding. Test results clearly manifest that the method we have developed outputs more accurate electron density profiles due to improvements to acquire the profiles of the layers below the F₂ region. Our study can further improve the current BSI inversion methods on the reconstruction of 2-D electron density distribution in a vertical plane aligned with the direction of sounding.     

Validation of B2bot in the NeQuick model during high solar activity at Hainan station

NeQuick ionospheric electron density model, which has been developed to version 2, produces the full electron density profile in the ionosphere. Each part of the profile is modeled using Epstein layer formalism. Simple empirical relations are used to compute the thicknesses of each layer. In order to validate the B2_bot parameter in the NeQuick model during high solar activity, we use the data at Hainan, China (109.1°E, 19.5°N; Geomagnetic coordinates: 178.95°E, 8.1°N), measured with DPS-4, and study the diurnal and seasonal variations of B2_bot, ΔB2 (B2_best − B2_NeQuick 2) and the seasonal median values of B2_best/B2_NeQuick 2 at that region. The results show that, (1) The differences between B2_best and B2_NeQuick 2 have diurnal and seasonal variations. (2) The diurnal variations of B2_NeQuick 2 are smaller than those of B2_best. (3) Generally, except for early morning the experimental values are properly reproduced. (4) Generally, during morning the NeQuick model has an underestimation. The magnitude of underestimation varies with LT and season.     

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.     

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

Improvements in bottomside electron density definition in the Autoscala program

Some improvements introduced in the Autoscala program are presented. They include improvements in E valley modeling of the electron density profile N_e(h), and in the link between the E valley and bottom-side F regions. An abrupt variation in N_e(h) generated by the previous version of Autoscala under night conditions has been eliminated.A series of ionograms recorded by the Millstone Hill digisonde (42.6°, 288.5°) were automatically interpreted by the previous version of Autoscala and by the new one. Data from Incoherent Scatter Radar (ISR) were used to comparatively assess the performance of the two versions. For this purpose, the root mean square errors (RMSEs) of the N_e(h) provided by Autoscala were calculated relative to the corresponding values provided by ISR.A more accurate overall modeling of N_e(h) was achieved by the new Autoscala version (RMSE = 0.51 MHz for the new version against RMSE = 0.67 MHz for the previous one).     

Adaptation of the bottomside electron density profile of the IRI model to data of topside radiosounding

A method is proposed for reconstructing the electron density profiles N(h) of the IRI model from ionograms of topside satellite sounding of the ionosphere. An ionograms feature is the presence of traces of signal reflection from the Earth's surface. The profile reconstruction is carried out in two stages. At the first stage, the N(h) –profile is calculated from the lower boundary of the ionosphere to the satellite height (total profile) by the method presented in this paper using the ionogram. In this case, the monotonic profile of the topside ionosphere is calculated by the classical method. The profile of the inner ionosphere is represented by analytical functions, the parameters of which are calculated by optimization methods using traces of signal reflection, both from the topside ionosphere and from the Earth. At the second stage, the profile calculated from the ionogram is used to obtain the key parameters: the height of the maximum hmF2 of the F2 layer, the critical frequency foF2, the values of B0 and B1, which determine the profile shape in the F region in the IRI model. The input of key parameters, time of observation, and coordinates of sounding into the IRI model allows obtaining the IRI-profile corrected to real experimental conditions. The results of using the data of the ISIS-2 satellite show that the profiles calculated from the ionograms and the IRI profiles corrected from them are close to each other in the inner ionosphere and can differ significantly in the topside ionosphere. This indicates the possibility of obtaining a profile in the inner ionosphere close to the real distribution, which can significantly expand the information database useful for the IRTAM (IRI Realmax Assimilative Modeling) model. The calculated profiles can be used independently for local ionospheric research.     

Thermosphere–ionosphere response to a severe magnetic storm: A case study

This paper reports the response of the ionosphere–thermosphere system to an intense geomagnetic storm. For that, data taken by instruments on board Dynamic Explorer 2 at heights of the F2-layer (molecular nitrogen N₂ and atomic oxygen O compositions, neutral temperature Tg and electron density Ne) were used. The ionospheric response is characterized by a negative storm effect expanding from mid–high to low latitude. It is observed during this severe geomagnetic storm that negative effects were caused mainly by an increase in molecular nitrogen composition N₂ and almost no changes in atomic oxygen composition O.     

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

Bottomside N(h) profiles over Vietnam and their comparison with IRI

Over 70 ionograms were selected from the ground-based ionosonde in Hanoi, Vietnam at times when the AE-E satellite passed above the station. N(h) profiles were calculated from ionograms and thus compared with the IRI model and local AE-E ion densities. For the purpose of the presentation, 15 profiles were selected to cover all seasons and local times. The comparison shows that the observed daytime N(h) profiles have significantly higher gradients below the F2 peak density while at night the observed and the IRI profiles almost coincide. The difference is more pronounced in summer than in equinox. Wintertime comparisons are quite limited and do not lead to reliable conclusions. N_mF2 values taken from the CCIR programme are compared with those observed. The comparison shows that the deviations are not so large.     

Variation of rain drop size distribution with rain rate at a few coastal and high altitude stations in southern peninsular India

Rain drop size distribution (DSD) was measured at four places in Southern India {Thiruvananthapuram, Kochi, Munnar and Sriharikota (SHAR)} using a Joss–Waldvogel (JW) impact type disdrometer. The data for each minute were corrected for dead time errors and rain rate was computed from the corrected data. The data for a whole month were then sorted according to rain rate (R) into several classes ranging from 0.1 to >100 mm/h. The average DSD in each class was computed, and the lognormal distribution function was fitted to the average. In all the cases, the function fitted the data very well. The fit parameters were found to have dependence on rain rate. The total number of drops (N_T), the geometric mean diameter (D_g) and the standard geometric deviation (σ) were also computed from the fit parameters. The standard geometric deviation (σ) was found to be more or less constant with rain rate at all the sites and in all months. The other two parameters (N_T and D_g) were found to vary exponentially with rain rate except in Munnar, a high altitude station. At Thiruvananthapuram, in most of the months, N_T increased exponentially with rain rate up to some value of R, which was different in different months, and then remained more or less constant or decrease slightly. In all cases, the variation of N_T and D_g was such that N_TD_g³ increased linearly with rain rate.     

Comparing the improved Di Giovanni/Radicella model with sounding-based electron density profiles and with the IRI model

The Di Giovanni/Radicella model (DGR) /1/ determines a bottom side electron densty profile alone from the set of routinely scaled ionogram parameters foE, foF1, foF2 and M(3000)F2 and the total electron content; the smoothed sunspot number R12 appears in the calculation. Present designations are DGR2/2/ and DRR3 /3/ [see Appendix]; they are valid in the northern hemisphere. DGR is compared with electron density profiles derived from ionograms obtained at Juliusruh (54.6°N, 13.4°E), and with the (URSI-based) IRI90 at different conditiones. Experimental total electron content (TEC) data are compared to both models. At the considered station, the profiles obtained by both models are reasonably in agreement amongst themselves and with the experimental data.

The TEC derived from the DGR3 model is in good agreement with experimental TEC, whereas, at high solar activity, IRI90 gives too high TEC values, especially during daytime.     

Bottom side profiles for two close stations at the southern crest of the EIA: Differences and comparison with IRI-2012 and NeQuick2 for low and high solar activity

Bottom side electron density profiles for two stations at the southern crest of the Equatorial Ionization Anomaly (EIA), São José dos Campos (23.1°S, 314.5°E, dip latitude 19.8°S; Brazil) and Tucumán (26.9°S, 294.6°E, dip latitude 14.0°S; Argentina), located at similar latitude and separated by only 20° in longitude, have been compared during equinoctial, winter and summer months under low (year 2008, minimum of the solar cycle 23/24) and high solar activity (years 2013–2014, maximum of the solar cycle 24) conditions. An analysis of parameters describing the bottom side part of the electron density profile, namely the peak electron density NmF2, the height hmF2 at which it is reached, the thickness parameter B₀ and the shape parameter B₁, is carried out. Further, a comparison of bottom side profiles and F-layer parameters with the corresponding outputs of IRI-2012 and NeQuick2 models is also reported. The variations of NmF2 at both stations reveal the absence of semi-annual anomaly for low solar activity (LSA), evidencing the anomalous activity of the last solar minimum, while those related to hmF2 show an uplift of the ionosphere for high solar activity (HSA). As expected, the EIA is particularly visible at both stations during equinox for HSA, when its strength is at maximum in the South American sector. Despite the similar latitude of the two stations upon the southern crest of the EIA, the anomaly effect is more pronounced at Tucumán than at São José dos Campos. The differences encountered between these very close stations suggest that in this sector relevant longitudinal-dependent variations could occur, with the longitudinal gradient of the Equatorial Electrojet that plays a key role to explain such differences together with the 5.8° separation in dip latitude between the two ionosondes. Furthermore at Tucumán, the daily peak value of NmF2 around 21:00 LT during equinox for HSA is in temporal coincidence with an impulsive enhancement of hmF2, showing a kind of “elastic rebound” under the action of the EIA. IRI-2012 and NeQuick2 bottom side profiles show significant deviations from ionosonde observations. In particular, both models provide a clear underestimation of the EIA strength at both stations, with more pronounced differences for Tucumán. Large discrepancies are obtained for the parameter hmF2 for HSA during daytime at São José dos Campos, where clear underestimations made by both models are observed. The shape parameter B₀ is quite well described by the IRI-2012 model, with very good agreement in particular during equinox for both stations for both LSA and HSA. On the contrary, the two models show poor agreements with ionosonde data concerning the shape parameter B₁.