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.     

Comparison of different N(h) profiles with mid-latitude ionospheric observations

After inversion, N(h) profiles obtained from ionograms that had been recorded during high solar activity at two mid-latitutde stations have been compared with those derived from IRI90 and DGR ionospheric models. A small data set has been selected such that both geomagnetically quiet and disturbed conditions are represented.     

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.     

New input to IRI from the worldwide digisonde network

Electron density profiles derived from Digisonde ionograms at Argentia, Millstone Hill, Wallops Island, Bermuda, Dourbes and Karachi are compared with IRI model prediction. Four months of data for 1989/90 were analyzed. For a number of station/months N(h) profiles were available every 15 or 30 minutes providing a good statistical database for the evaluation of the IRI model in terms of diurnal and seasonal variations. The data presented here are part of the VIM study (Validation of Ionospheric Models) initiated by the URSI Working Group G3 on Ionospheric Informatics.     

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.     

Validating ionospheric models with measured electron density profiles

On behalf of an URSI Working Group 3 initiated study (VIM), three ionospheric models, IRI, PL/PRISM and FLIP, are compared with electron density profiles derived from ionograms Millstone Hill. Four months of data in 1989/90 were analyzed. For most of the time, N(h) profiles were available every 15 minutes providing a good statistical database for the evaluation of the ionospheric models in terms of diurnal and seasonal variations.     

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.     

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.     

A statistical study of the F2 layer stratification at the northern equatorial ionization anomaly

Ionograms recorded at Puer station (PUR, 22.7°N, 101.05°E, Dip Latitude 12.9°N) in the Southwest of China from January 2015 to December 2016 were used to study characteristics of the F2 layer stratification at the northern equatorial ionization anomaly. Ionosonde observations show that the development of the F2 layer stratification is different under different conditions. Both the upward and downward movement of the F2 layer stratification could be observed. The F2 layer stratification could occur both at daytime and nighttime. The new cusp could originate from different positions on ionograms. Moreover, statistical results indicate that the F2 layer stratification occurred later in the winter than in other seasons at daytime, it occurred frequently in the local spring, and most of ionograms with the F2 layer stratification at post-midnight occurred in March and April. Our results also show that the F2 layer stratification has a correlation with solar activity.     

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.     

Variability of the behavior of the bottomside (B0, B1) parameters obtained from the ground-based ionograms at China’s low latitude station

In this paper, data of (B0, B1) parameters deduced from the electron density profiles that are inverted from the ionograms recorded at Hainan (19.4°N, 109.0°E), China during a three year period from March 2002 to February 2005 are used to study the diurnal and seasonal variation of (B0, B1) parameters at low latitude. The observational results are compared with the IRI2001 model predictions. Variability study of (B0, B1) in terms of percentage ratio of the inter-quartiles to the median values and correlative analysis between (B0, B1) parameters and other ionospheric characteristics such as hmF2 and M(3000)F2 are also made. Our present study showed that: (1) for daytime hours, the IRI2001 model results with new table option (B0_Tab) is in a better agreement with the observational results (B0_Obs) than the IRI2001 model results with Gulyaeva option (B0_Gul) for summer season, whereas B0_Gul is in a better agreement with B0_Obs than B0_Tab for winter season. For nighttime, in general, B0_Gul is in a better agreement with B0_Obs than B0_Tab. For other occasions, both B0_Tab and B0_Gul showed some systematic deviations from the observational ones. Moreover, the deviations of B0_Tab and B0_Gul from B0_Obs showed opposite trends; (2) the monthly upper (lower) quartiles of (B0, B1) parameter showed a good linear relationship with the monthly median values, this makes it possible to do the regression analysis between the monthly upper (lower) quartiles and the monthly median values, which can give a measure of the variability of these parameters. In terms of the percentage ratio of the inter-quartiles to the median values, the variability of B0 showed a diurnal variation ranging between 22% and 36% with maximum value occurring at pre-sunrise hours, whereas the variability of B1 showed a diurnal variation ranging between 15% and 30% with higher value by daytime than at night; (3) B0 shows high linear correlative relationships with hmF2 and M(3000)F2 for most of the local time period of a day except for a few hours around midnight, whereas B1 showed high linear correlations with B0, hmF2 for daytime hours, but not for nighttime hours. This suggests that it maybe is possible to obtain the synthetic database of (B0, B1) parameter or to construct the model of (B0, B1) using database of hmF2 or M(3000)F2 which is much easier to obtain from experimental measurements.     

Comparison of ion composition from the goddard comprehensive ionosphere database with the international reference ionosphere model

We compared the IRI values of T_e, N_e, T_i, O⁺, H⁺, He⁺, O₂₊, and NO⁺ with AE-C values, obtained from the Goddard Comprehensive Ionosphere Database (GCID), composed of data from the satellites, AE-B, OGO-6, ISIS-2, AE-C, AE-D, and AE-E. O⁺ - H⁺ transition levels were derived from the IRI and AE-C altitude profiles. Some discrepancies were found between IRI and the AE-C data. The IRI electron density was found to be about a factor of 2 higher than the data. The H⁺ composition agreed best among the IRI ions, with an average AE-C/IRI ratio of 1.05. The temperatures of both electrons and ions agreed quite well: the average ratios of AE-C/IRI was found to be .99 for electrons and 1.17 for ions.     

Implementation of a new characteristic parameter into the IRI sub-peak electron density profile

The shape of electron density profile in the International Reference Ionosphere could be improved significantly if the height h_g and electron density N_g of the F region sub-peak inflexion point were entered in the set of the profile standard parameters. To study variations of these important parameters, the N(h) analysis of the statistically-summarized ionograms at the latitudes of 40–80°N of the Eastern hemisphere has been carried out for the two-hours intervals of local time, three seasons (winter, equinox and summer) and two levels of solar activities characterized by Covington indicesF_10.7 = 100 and 200. It is shown that the parameters of the inflexion point can be expressed in most cases via the peak parameters of the F2 layer ash_g= 0.8 hmF2 and N_g= 0.5 NmF2.     

N(h) profiles from ground and “Interkosmos-19” satellite data

During the years of high solar activity (1979–1980) simultaneous soundings of the ionosphere were made by the Bulgarian ground stations at Sofia and Michurin and by the satellite Interkosmos 19. The N(h) profiles for the ground stations were computed in association with IZMIRAN, Moscow and those for the satellite by SIB-IZMIRAN.During some of the satellite transit orbits the electron densities at the maximum of the F2 layer coincide precisely with those obtained from ground ionograms. In other cases differences were observed, perhaps because of the differences in time; the ionospheric stations at Sofia and Michurin worked on the ordinary 15-minute regime.     

Physical and model interpretation of HF radio propagation on the St. Petersburg–Longyearbyen (Svalbard) path

HF radio wave observations have been carried out with an oblique ionospheric sounding (OIS) method on the radio path from St. Petersburg to Longyearbyen (Svalbard), and experimental ionograms were obtained for December 2001. These ionograms have been analysed to investigate the impact of the main ionospheric trough (MIT) and magnetic disturbances on the signals on this path. The observations during weakly disturbed (K_р = 2) magnetic conditions on 14–15 December 2001 were compared with predictions from ray-tracing through a numerical model of the ionosphere. The ray-tracing computer program synthesizes the OIS ionograms by means of the “shooting method”. This method calculates trajectories of HF radio waves for different values of elevation angle and transmission frequency. There was a variety of calculated trajectories, from which we choose those which reach the receiver, and the selected paths provide a synthesis of the oblique ionograms. To simulate HF radio wave propagation, we apply a three-dimensional distribution of the electron density calculated with the mathematical model of the high-latitude ionosphere developed in the Polar Geophysical Institute (PGI). These numerical simulations permit us to interpret specific peculiarities of the OIS data such as abnormal propagation modes, increased delays of signals, enhanced MOF (maximum observed frequency) values etc. New results of the study are summarised as follows. (1) An unusual feature of the propagation along the path is the change of propagation mechanism during substorms on entering a path midpoint (or 1-hop reflection point) to the MIT. (2) Even weak substorms, having the distinguished intensities, lead to the appearance of different types of irregularities observed by the CUTLASS radar and therefore to the different propagation modes and F2MOF values. (3) The PGI model of the ionosphere was first used for ray-tracing at high latitudes. The model results are basically in a good qualitative agreement with experimental observations. This model provides the satisfactory agreement between the calculated and experimental F2MOF values while not correctly representing the fine structure of the experimental OIS ionograms at night. An agreement between the calculated and experimental data is better for day and evening hours than at night.     

Cross-hemisphere comparison of mid-latitude ionospheric variability during 1996–2009: Juliusruh vs. Hobart

Analysis of a long-time series of hourly median characteristics of the ionospheric plasma at two mid-latitude locations in the Northern and Southern hemisphere, Juliusruh (54.6N; 13.4E) and Hobart (42.9S; 147.3E), reveals patterns of their synchronous and independent variability. We studied timelines of GPS vTEC, ionogram-derived F2-layer peak electron density NmF2, ionospheric equivalent slab thickness τ, and their ratios at two locations during the complete 23rd solar cycle and its following period of the extremely low solar activity in 2008–2009. This study has also involved the comparative analysis of the observed data versus the model predictions by IRI-2012. During the high solar activity in 2000–2002, seasonal variations show a complicated cross-hemisphere behavior influenced by the winter and semi-annual anomalies, with the largest noon-time values of TEC and NmF2 observed around equinoxes. Strength of the winter anomaly in NmF2 was significantly greater at Juliusruh in comparison with Hobart. The winter anomaly in GPS vTEC values was much weaker than in NmF2 for the Northern hemisphere mid-latitudes and was entirely absent at the Southern hemisphere. Cross-hemisphere analysis of the equivalent slab thickness shows its clear seasonal dependence for all levels of solar activity: the day-time maximum τ_max is observed during local summer, whereas the day-time minimum τ_min is observed during local winter. The night-time values of τ were higher compared to the day-time values during the winter and equinox seasons. Comparative model-data study shows rather good IRI performance of the day-time NmF2 for mid-latitudes of both hemispheres and rather noticeable overestimations for the mid-night NmF2 values during high solar activity. Analysis of IRI vTEC demonstrates the model limitations, related with the absence of the plasmaspheric part, and actual demand in a reliable and standard ionosphere–plasmasphere model for analysis of GPS vTEC.     

Comparison of ionospheric data from Soviet and Chinese stations with the IRI-86 and the faim models

Median values of ionospheric data from Moscow(55.5° N, 37.3° E), Tashkent (41.3° N, 69.6° E), Beijing (40.0° N, 116.3° E) and Guangzhou (23.1° N, 113.3° E) for different years at LT: 00.00, 06.00, 12.00 and 18.00 were used to produce vertical profiles of electron density with the analytic procedure introduced by Giovanni and Radicella. These profiles were compared with the IRI-86 model profiles adjusted to the F2 experimental peak values and with those obtained with the Fully Analytic Ionospheric Model (FAIM). The main differences with the IRI model are related to the shape of the topside profile under certain conditions and to the presence of the F1 layer in the median values used. The differences with the FAIM model appears to be more marked. The results are discussed.     

Comparison of ionospheric F2 peak parameters foF2 and hmF2 with IRI2001 at Hainan

Monthly median values of foF2, hmF2 and M(3000)F2 parameters, with quarter-hourly time interval resolution for the diurnal variation, obtained with DPS4 digisonde at Hainan (19.5°N, 109.1°E; Geomagnetic coordinates: 178.95°E, 8.1°N) are used to investigate the low-latitude ionospheric variations and comparisons with the International Reference Ionosphere (IRI) model predictions. The data used for the present study covers the period from February 2002 to April 2007, which is characterized by a wide range of solar activity, ranging from high solar activity (2002) to low solar activity (2007). The results show that (1) Generally, IRI predictions follow well the diurnal and seasonal variation patterns of the experimental values of foF2, especially in the summer of 2002. However, there are systematic deviation between experimental values and IRI predictions with either CCIR or URSI coefficients. Generally IRI model greatly underestimate the values of foF2 from about noon to sunrise of next day, especially in the afternoon, and slightly overestimate them from sunrise to about noon. It seems that there are bigger deviations between IRI Model predictions and the experimental observations for the moderate solar activity. (2) Generally the IRI-predicted hmF2 values using CCIR M(3000)F2 option shows a poor agreement with the experimental results, but there is a relatively good agreement in summer at low solar activity. The deviation between the IRI-predicted hmF2 using CCIR M(3000)F2 and observed hmF2 is bigger from noon to sunset and around sunrise especially at high solar activity. The occurrence time of hmF2 peak (about 1200 LT) of the IRI model predictions is earlier than that of observations (around 1500 LT). The agreement between the IRI hmF2 obtained with the measured M(3000)F2 and the observed hmF2 is very good except that IRI overestimates slightly hmF2 in the daytime in summer at high solar activity and underestimates it in the nighttime with lower values near sunrise at low solar activity.     

Comparison of ionosphere characteristic parameters obtained by ionosonde with IRI-2007 model over Southeast Asia

In this work, the foF2 and hmF2 parameters at the conjugate points near the magnetic equator of Southeast Asia are studied and compared with the International Reference Ionosphere (IRI) model. Three ionosondes are installed nearly along the magnetic meridian of 100°E; one at the magnetic equator, namely Chumphon (10.72°N, 99.37°E, dip angle 3.0°N), and the other two at the magnetic conjugate points, namely Chiang Mai (18.76°N, 98.93°E, dip angle 12.7°N) and Kototabang (0.2°S, 100.30°E, dip angle 10.1°S). The monthly hourly medians of the foF2 and hmF2 parameters are calculated and compared with the predictions obtained from the IRI-2007 model from January 2004 to February 2007. Our results show that: the variations of foF2 and hmF2 predicted by the IRI-2007 model generally show the similar feature to the observed data. Both parameters generally show better agreement with the IRI predictions during daytime than during nighttime. For foF2, most of the results show that the IRI model overestimates the observed foF2 at the magnetic equator (Chumphon), underestimates at the northern crest (Chiang Mai) and is close to the measured ones at the southern crest of the EIA (Kototabang). For hmF2, the predicted hmF2 values are close to the hmF2(M3000F2_OBS) during daytime. During nighttime, the IRI model gives the underestimation at the magnetic equator and the overestimation at both EIA crests. The results are important for the future improvements of the IRI model for foF2 and hmF2 over Southeast Asia region.     

Parameterized Regional Ionospheric Model and a comparison of its results with experimental data and IRI representations

We describe a Parameterized Regional Ionospheric Model (PARIM) to calculate the spatial and temporal variations of the ionospheric electron density/plasma frequency over the Brazilian sector. The ionospheric plasma frequency values as calculated from an enhanced Sheffield University Plasmasphere–Ionosphere Model (SUPIM) were used to construct the model. PARIM is a time-independent 3D regional model (altitude, longitude/local time, latitude) used to reproduce SUPIM plasma frequencies for geomagnetic quiet condition, for any day of the year and for low to moderately high solar activity. The procedure to obtain the modeled representation uses finite Fourier series so that all plasma frequency dependencies can be represented by Fourier coefficients. PARIM presents very good results, except for the F region peak height (hmF2) near the geomagnetic equator during times of occurrence of the F₃ layer. The plasma frequency calculated by IRI from E region to bottomside of the F region present latitudinal discontinuities during morning and evening times for both solar minimum and solar maximum conditions. Both the results of PARIM and the IRI for the E region peak density show excellent agreement with the observational values obtained during the conjugate point equatorial experiment (COPEX) campaign. The IRI representations significantly underestimate the foF2 and hmF2 compared to the observational results over the COPEX sites, mainly during the evening–nighttime period.