The state of the art in D-region IRI modelling

While describing the formulation of IRI-79, in so far as it is applicable to the sub E-peak electron density profile, an attempt is made to compare its predictions with the available experimental evidence. The improvements needed for D-region electron density modelling are summarized. The current D-region IRI modelling effort is illustrated with a typical example.     

Improvement of the solar-cycle variation of IRI lower ionosphere models by means of radio wave propagation data

Data on day-time and night-time radio wave absorption in the frequency range 50 to 2614 kHz, obtained in long-term observational programmes in Central Europe, are compared with corresponding absorption values calculated from electron density profiles of the International Reference Ionosphere (IRI-1979) using the full-wave method.Discrepancies between calculated and observed absorption values were found for the diurnal and the solar-cycle variation, the amplitudes of the solar-cycle variation of absorption being considerably larger than the observed variation.A modification of the solar-activity dependence of the D-region electron density parameters is derived, which provides an improvement of the solar-cycle variation as well as the diurnal variation of the IRI electron density profiles.     

Global empirical models of the density peak height and of the equivalent scale height for quiet conditions

Monthly average electron density profiles have been calculated from hourly electron density N(h) recorded in 26 digisonde stations distributed worldwide encompassing the time interval 1998–2006. The ionospheric electron density peak height of the F2 region, hmF2, and the effective scale height at the hmF2, Hm, deduced from average profiles have been analyzed to obtain the quiet-time behavior and have been analytically modeled by the spherical harmonic analysis (SH) technique using the modip latitude as the coordinate of the reference system. The coefficients of the SH models of hmF2 and Hm are bounded to the solar activity, and the temporal and seasonal variations are considered by Fourier expansion of the coefficients. The SH models provide a tool to predict hmF2 and Hm located anywhere in the range of latitudes between of 70°N and 70°S and at any time. The SH analytical model for hmF2 improves the fit to the observations by 10% in average compared to the IRI prediction, and it might improve the IRI prediction of hmF2 by more than 30% at high and low latitudes. The analytical model for Hm predicts the quiet behavior of the effective scale height with accuracy better than 15% in average which enables to obtain a good estimation of vertical profiles. These results could be useful to estimate information for the topside profile formulation.     

Response of the mid-latitude D-region ionosphere to the total solar eclipse of 22 July 2009 studied using VLF signals in South Korean peninsula

In this paper, we analyze VLF signals received at Busan to study the the D-region changes linked with the solar eclipse event of 22 July 2009 for very short (∼390 km) transmitter–receiver great circle path (TRGCP) during local noon time 00:36–03:13 UT (09:36–12:13 KST). The eclipse crossed south of Busan with a maximum obscuration of ∼84%. Observations clearly show a reduction of ∼6.2 dB in the VLF signal strength at the time of maximum solar obscuration (84% at 01:53 UT) as compared to those observed on the control days. Estimated values of change in Wait ionospheric parameters: reflection height (h′) in km and inverse scale height parameter (β) in km⁻¹ from Long Wave Propagation Capability (LWPC) model during the maximum eclipse phase as compared to unperturbed ionosphere are 7 km and 0.055 km⁻¹, respectively. Moreover, the D-region electron density estimated from model computation shows 95% depletion in electron density at the height of ∼71 km. The reflection height is found to increase by ∼7 km in the D-region during the eclipse as compared to those on the control days, implying a depletion in the Lyman-α flux by a factor of ∼7. The present observations are discussed in the light of current understanding on the solar eclipse induced D-region dynamics.     

IRI 2001/90 TEC predictions over a low latitude station

Total electron content (TEC) over Tucumán (26.9°S, 294.6°W) measured with Faraday technique during the high solar activity year 1982, is used to check IRI 2001 TEC predictions at the southern crest of the equatorial anomaly region. Comparisons with IRI 90 are also made. The results show that in general IRI overestimates TEC values around the daily minimum and underestimates it the remaining hours. Better predictions are obtained using ground ionosonde measurements as input coefficients in the IRI model. The results suggest that for hours of maximum TEC values the electron density profile is broader than that assumed by the model. The main reason for the disagreement would be the IRI shape of the electron density profile.     

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.     

Progress in ionospheric informatics based on electron-density profile analysis of ionograms

Evolution of techniques for electron-density profile analysis of ionograms attests to a close coupling with ionospheric modelling. Advanced achievements in both fields are used in the International Reference Ionosphere, IRI, the selection being made by criteria of internal consistency. This paper presents the results of comparisons of modern methods of calculating electron-density profiles from ground-based vertical-incidence ionograms. It is shown that the shapes of the resulting electron-density profiles depend upon the peak height and electron density of the E and F-layers. Formulae for the F2-layer sub-peak electron density, the E-F valley and their relation to the peak height of the ionosphere are presented as functions of solar-zenith angle. Implementation of such relationships allows refinement and considerable simplification of representations in the IRI, depending on the solar and geophysical controlling parameters.     

Calibration of IRI–ITU-R peak density and height over the oceans with topside sounding data

Accuracy of IRI electron density profile depends on the F2 layer peak density and height converted by empirical formulae from the critical frequency and M3000F2 factor provided by the ITU-R (former CCIR). The CCIR/ITU-R maps generated from ground-based ionosonde measurements suffer from model assumptions, in particular, over the oceans where relatively few measurements are available due to a scarcity of ground-based ionosondes. In the present study a grid-point calibration of IRI/ITU-R maps for the foF2 and hmF2 over the oceans is proposed using modeling results based on the topside true-height profiles provided by ISIS1, ISIS2, IK-19 and Cosmos-1809 satellites for the period of 1969–1987. Topside soundings results are compared with IRI and the Russian standard model of ionosphere, SMI, and grouped to provide an empirical calibration coefficient to the peak density and height generated from ITU-R maps. The grid-point calibration coefficients maps are produced in terms of the solar activity, geodetic latitude and longitude, universal time and season allowing update of IRI–ITU-R predictions of the F2 layer peak parameters.     

A topside investigation over a mid-latitude digisonde station in Cyprus

We examine the systematic differences between topside electron density measurements recorded by different techniques over the low-middle latitude operating European station in Nicosia, Cyprus (geographical coordinates: 35.14^oN, 33.2^oE), (magnetic coordinates 31.86^oN, 111.83 ^oE). These techniques include space-based in-situ data by Langmuir probes on board.European Space Agency (ESA) Swarm satellites, radio occultation measurements on board low Earth orbit (LEO) satellites from the COSMIC/FORMOSAT-3 mission and ground-based extrapolated topside electron density profiles from manually scaled ionograms. The measurements are also compared with International Reference Ionosphere Model (IRI-2016) topside estimations and IRI-corrected NeQuick topside formulation (method proposed by Pezzopane and Pignalberi (2019)). The comparison of Swarm and COSMIC observations with digisonde and IRI estimations verifies that in the majority of cases digisonde underestimates while IRI overestimates Swarm observations but in general, IRI provides a better topside representation than the digisonde. For COSMIC and digisonde profiles matched at the F layer peak the digisonde systematically underestimates topside COSMIC electron density values and the relative difference between COSMIC and digisonde increases with altitude (above hmF2), while IRI overestimates the topside COSMIC electron density but after a certain altitude (~150 km above hmF2) this overestimation starts to decrease with altitude. The IRI-corrected NeQuick underestimates the majority of topside COSMIC electron density profiles and relative difference is lower up to approximately 100 km (above the hmF2) and then it increases. The overall performance of IRI-corrected NeQuick improves with respect to IRI and digisonde.     

The interest of information obtained by incoherent scatter technique

Incoherent scatter measurements are of great value for establishing and improving the IRI. This holds, in particular, for the topside electron density profile, for the valley occurring under certain conditions between regions E and F, and for the electron and ion temperature profiles. Extension in time of the observations could be very helpful for future work on IRI.     

Response of low latitude D-region ionosphere to the total solar eclipse of 22 July 2009 deduced from ELF/VLF analysis

Response of the D-region of the ionosphere to the total solar eclipse of 22 July 2009 at low latitude, Varanasi (Geog. lat., 25.27° N; Geog. long., 82.98° E; Geomag. lat. = 14° 55’ N) was investigated using ELF/VLF radio signal. Tweeks, a naturally occurring VLF signal and radio signals from various VLF navigational transmitters are first time used simultaneously to study the effect of total solar eclipse (TSE). Tweeks occurrence is a nighttime phenomena but the obscuration of solar disc during TSE in early morning leads to tweek occurrence. The changes in D-region ionospheric VLF reflection heights (h) and electron density (n_e: 22.6–24.6 cm⁻³) during eclipse have been estimated from tweek analysis. The reflection height increased from ∼89 km from the first occurrence of tweek to about ∼93 km at the totality and then decreased to ∼88 km at the end of the eclipse, suggesting significant increase in tweek reflection height of about 5.5 km during the eclipse. The reflection heights at the time of totality during TSE are found to be less by 2–3 km as compared to the usual nighttime tweek reflection heights. This is due to partial nighttime condition created by TSE. A significant increase of 3 dB in the strength of the amplitude of VLF signal of 22.2 kHz transmitted from JJI-Japan is observed around the time of the total solar eclipse (TSE) as compared to a normal day. The modeled electron density height profile of the lower ionosphere depicts linear variation in the electron density with respect to solar radiation as observed by tweek analysis also. These low latitude ionospheric perturbations on the eclipse day are discussed and compared with other normal days.     

Latitudinal influences on the quiet daytime D-region

Four versions of a steady-state quiet D-region model are presented. They differ from each other as a result of latitudinal differences in total neutral particle concentrations, nitric oxide concentrations and cosmic ray ionization rates. The total ion concentration profiles of all four versions have minima near 70 km which range from about 10⁸ m^?3 at high latitudes to 3.5 × 10⁷ m^?3 at equatorial latitudes for a solar zenith angle of 60°. Neutral density differences among the four cases result in important vertical shifts for the respective D-region profiles relative to one another. A “C-layer” is evident for the high and mild-latitude models at large solar zenith angles. The altitude where the negative ion/electron concentrations ratio is unity varies from about 63 to 67 km. The computed results are compared briefly with the extensive data base in the literature.     

Prediction of total electron content using the International Reference Ionosphere

The International Reference Ionosphere (IRI) is a model of the ionosphere, based on experimental data, which has been proposed as a standard ionospheric model. As such, it should be tested extensively to determine its range of validity. One of the ways in which the electron denisty profile given by the IRI, especially above the peak of the F layer, can be tested is to compare calculated and observed values of total electron content (TEC). We have therefore studied the discrepancies between calculated and observed values of TEC recorded at 15 stations covering a wide range of longitudes and latitudes, mainly in the northern hemisphere, and mainly for high levels of solar activity. W have found that the IRI produces reasonably accurate values of TEC at mid and high latitudes, but that it greatly underestimates the daytime values of TEC at low latitudes. We conclude therefore that the daytime electron density profile given by the IRI is reasonably accurate at mid and high latitudes, at least above the peak of the F2 layer. The situation at low latitudes clearly requires more work, and we have suggested two possible lines of study. The generally low discrepancies at night indicate that the night-time electron density profiles given by the IRI correspond fairly closely to the actual profiles.     

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.     

Accuracy of IRI profiles of ionospheric density and temperatures derived from comparisons to Kharkov incoherent scatter radar measurements

The incoherent scatter radar (ISR) facility in Kharkov, Ukraine (49.6°N, 36.3°E) measures vertical profiles of electron density, electron and ion temperature, and ion composition of the ionospheric plasma up to 1100 km altitude. Acquired measurements constitute an accurate ionospheric reference dataset for validation of the variety of models and alternative measurement techniques. We describe preliminary results of comparing the Kharkov ISR profiles to the international reference ionosphere (IRI), an empirical model recognized for its reliable representation of the monthly-median climatology of the density and temperature profiles during quiet-time conditions, with certain extensions to the storm times. We limited our comparison to only quiet geomagnetic conditions during the autumnal equinoxes of 2007 and 2008. Overall, we observe good qualitative agreement between model and data both in time and with altitude. Magnitude-wise, the measured and modeled electron density and plasma temperatures profiles appear different. We discovered that representation accuracy improves significantly when IRI is driven by observed-averaged values of the solar activity index rather than their predictions. This result motivated us to study IRI performance throughout protracted solar minimum of the 24th cycle. The paper summarizes our observations and recommendations for optimal use of the IRI.     

利用太阳耀斑的电离层效应确定低纬度D层顶一氧化氮浓度

通过氧鎓离子化学计算了一次强太阳耀斑扰动期间低纬度(28°38′N)D层顶NO浓度。该方法适用于O+₂离子化学生存时间较短且负离子浓度与正离子相比可以忽略的80至87.5公里高度。计算结果表明,它相当好地符合Meira的实验观测剖面,只是85公里处的最小值特别弱。      

Midday bottomside electron density profiles during moderate solar activity and comparison with IRI-2001

Bottomside electron density (Ne-h) profiles during midday (10–14 h) are analyzed using modern digital ionosonde observations at a low-middle latitude station, New Delhi (28.6N, 77.2E, dip 42.4N), for the period from January 2003 to December 2003, pertaining to moderate solar activity (MSA). Each individual profile is normalized with respect to the peak height and density (hmF2, NmF2) of the F2-region. These profiles are compared with those obtained from the International Reference Ionosphere (IRI-2001) model. Bilitza [Bilitza, D. International Reference Ionosphere 2000. Radio Sci. 36 (2), 261–275, 2001] using both the options namely: Gulyaeva’s model [Gulyaeva, T.L. Progress in ionospheric informatics based on electron density profile analysis of ionograms. Adv. Space Res. 7 (6) 39–48, 1987] and B0 Tab. option [Bilitza, D., Radicella, S.M., Reinisch, B.W., Adeniyi, J.O., Mosert Gonzalez, M.E., Zhang, S.R., Obrou, O. New B0 and B1 models for IRI. Adv. Space Res. 25 (1), 89–95, 2000]. The study reveals that during summer and equinox, the IRI model with B0 Tab. option in general, produces better agreement with the observed median profiles, while the IRI predictions using Gulyaeva’s option, overestimate the electron density distribution at all the heights below the F2-peak. However, during winter, in general, the IRI model, using both the options, reveals shows fairly good agreement with the observations.     

The quiet nighttime low-latitude ionosphere as observed by TIMED/GUVI

In this paper, we examine the nighttime ionosphere climatology structure in the low latitude region and discrepancies between Global Ultraviolet Imager (GUVI) observations and the IRI model predictions using (1) the magnetic zonal mean of electron number density as a function of altitude and magnetic latitude, (2) vertical electron density profiles at various levels of F10.7 index, (3) nighttime descent and magnitude decrease of the ionosphere, (4) point-to-point comparisons of F-peak height (hmF2) and density (NmF2), and (5) the magnetic longitudinal variations of hmF2 and NmF2. The data collected from the Thermosphere, Ionosphere, Mesosphere, Energetics, and Dynamics (TIMED) mission since its launch in December 2001 have provided great opportunities for many scientific investigations of the ionosphere. In this analysis, we investigate the climatology of the nighttime low-latitude ionosphere under low geomagnetic activity (kp ? 4) using the electron density profiles inferred from the airglow measurements obtained by the GUVI aboard the TIMED spacecraft and compared with the results obtained from IRI (International Reference Ionosphere) model-2001. The observed climatology is an essential tool for further understanding the electrodynamics in the low-latitude region and improving the model’s prediction capability. The time range of the GUVI data used in this study is from 2002 (day 053) to 2006 (day 304), and the IRI model predictions were produced at every GUVI location. The ionosphere observed is generally of greater density than what IRI predicts throughout the night for all four seasons for low and moderate solar activity while the model over-predicts the electron density near the F-region peak at high solar activity before midnight. Observations show that the height of the F-region peak has a steep descent from dusk to midnight and near midnight the height of layer is insensitive to solar conditions, significantly different than what is predicted by IRI. Longitudinal features shown in GUVI data are present in the low-latitude ionosphere after sunset and continue through to midnight after which the low-latitude ionosphere is largely zonally symmetric.