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.     

A comprehensive study on middle atmospheric thermal structure over a tropic and sub-tropic stations

The general characteristics of middle atmospheric thermal structure have been studied by making use of the Rayleigh lidar data collected over the period of about four years (1998–2001). Here, the data has been used from two different stations in the Indian sub-continent in tropics (Gadanki; 13.5°N, 79.2°E) and in sub-tropics (Mt. Abu; 24.5°N, 72.7°E). The observed monthly mean temperature profiles are compared with different model atmospheres (CIRA-86 and MSISE-90). We observed, the mean temperature profiles have closer agreement with MSISE-90 than CIRA-86. The temperature profiles measured by lidar and HALOE satellite overpass nearby lidar site are generally in agreement with each other. The systematic and statistical errors in deriving temperature are found to be uniform for both the stations, as 1 K at 50 km, 3 K at 60 km and 10 K at 70 km. The special features of mesospheric temperature inversion (MTI) and double stratopause structure (DBS) are also addressed for both the stations.     

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.     

Total electron content at low latitudes and its comparison with the IRI90

We used total electron content (TEC) data measured by Faraday rotation technique over Cachoeira Paulista (22.5°S, 45°W), in Brazil, to study the TEC variations with the solar flux at 10.7 cm (F10.7) and to compare the results with the IRI90 predictions. The data were divided into summer, equinox and winter. During the analysed period F10.7 varied from 66 up to 330. Our data shows that the observed TEC at 1600 LT (around the diurnal maximum) and at 0500 LT (around the diurnal minimum) increases with F10.7 until saturation is reached which occurs at F10.7≈210 to 220 for equinox and summer, and at F10.7≈180 for winter months. Comparison with the IRI90 predictions shows that IRI overestimates the TEC at 0500 LT for all solar flux values. At 1600 LT, IRI overestimates the observed TEC for low solar flux but underestimates it for high solar flux values.     

On the variabilities of the Total Electron Content (TEC) over the Indian low latitude sector

The Total Electron Content (TEC) from four locations in the Indian sector namely, Trivandrum (8.47°N, 76.91°E, Geomag.0.63°S, 0.3° dip), Waltair (17.7° N, 83.3°E, Geomag. 6.4°N, 20° dip), Bhopal (23.28°N, 77.34°E, Geomag.14.26°N, 33.2° dip), and Delhi (28.58°N, 77.21°E, Geomag.19.2°N, 43.4° dip) during a low sunspot year of 2004 are used to study the variabilities of the TEC. The day time TEC values are higher over Waltair and Bhopal compared to those at Trivandrum and Delhi. Considerable day-to-day variations in the diurnal values of TEC are observed at the anomaly crest locations. The observed GPS-TEC has been compared with the IRI-2007 model derived TEC considering three different options (IRI-2001, IRI-2001 corrected and Ne-Quick) available in the model for the topside electron density. The TEC derived with Ne-Quick and IRI-01 corrected options show better agreement with GPS-TEC while the TEC from IRI-01 method shows larger deviations. From the correlation analysis carried out between TEC value at 1300 h LT and solar indices parameters namely sunspot number (SSN), F10.7 and EUV, it is observed that the correlation is more during equinoctial months and less during summer months. The correlation coefficients observed over the anomaly locations, Bhopal and Delhi are lower compared to those at Trivandrum and Waltair.     

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.     

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.     

Structure of the middle atmosphere of Venus from analysis of Fourier spectrometer measurements aboard Venera 15

Infrared spectra of Venus measured by means of the Fourier spectrometer aboard Venera 15 orbiter were used for retrievals of temperature profiles of the atmosphere in the altitude range from 60 to 95 km. Monotonous profiles are typical for latitudes lower than 60° latitude, but on the dayside near the equator some traces of the upper atmospheric inversion were found in the profiles. At latitudes greater than 60°N the profiles contain an isothermal or inversional part between 10 and 100 hPa pressure levels. Temperature profiles inside the “warm dipole” maxima are the same as outside of this region and consequently these optical properties should be explained only by peculiarities of cloud structure there.     

Behavior of the equivalent slab thickness over three European stations

The total electron content (TEC) derived from the global positioning system (GPS) and the F2-layer peak electron density obtained from Digisonde data have been used to study the diurnal, seasonal and solar activity variations of the ionospheric equivalent slab thickness (τ) over three European stations located at Pruhonice (50.0°N, 15.0°E), Ebro (40.8°N, 0.5°E) and El Arenosillo (37.1°N, 353.3°E). The diurnal variation of the τ is characterized by daytime values lower than nighttime ones for all seasons at low solar activity while daytime values larger than nighttime characterizes the diurnal variation for summer at high solar activity. A double peak is noticeable at dusk and at dawn, better expressed for winter at low solar activity. The seasonal variations of τ depend on local time and solar activity, the daytime values of τ increases from winter to summer whereas nighttime values of τ show the opposite. The effect of the solar activity on τ depends on local time and season, there being very sensitive for winter nighttime values of τ. The results of this study are compared with those presented by other authors.     

Towards validation of SCIAMACHY lunar occultation NO2 vertical profiles

Vertical profiles of stratospheric nitrogen dioxide (NO₂) have been retrieved from moderate resolution lunar occultation transmission spectra measured by Scanning Imaging Absorption spectroMeter for Atmospheric CHartographY (SCIAMACHY) on board the European Environmental Satellite (ENVISAT). These measurements were taken over the high southern latitude of 50°–90° during the period of 2003–2005. To assess the accuracy of the retrieved NO₂ profiles, the SCIAMACHY nighttime NO₂ profiles were compared with NO₂ profiles retrieved from sunrise solar occultation spectra measured by the Halogen Occultation Experiment (HALOE) and the Stratospheric Aerosol and Gas Experiments II (SAGE II) using a photochemical correction model. The validation results show good agreement of SCIAMACHY lunar occultation NO₂ with scaled HALOE and SAGE II profiles. The relative mean differences (rmd) with scaled HALOE profiles are within −13% to +5% and standard deviations (rms) of the relative differences are within 3–19% between 25 and 38 km. The rmd and rms with scaled SAGE II NO₂ profiles are in the range of −9 to +7 and 10–17% respectively between 22 and 39 km.     

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.     

Bayesian approach to retrieval of vertical ozone profile from radiometry data

A technique for retrieving vertical distributions (profiles) of atmospheric gas constituents from data of passive remote sensing of the atmosphere is proposed. The goal of the technique based on the statistical (Bayesian) approach to solution of inverse problems is construction of probability distribution for a sought quantity throughout the interval of the studied heights. It is assumed that initial data contain measurement noise, and a priori information about properties of the profile is used. It is proposed to approximate the sought profile by a function in the form of an artificial neural network. This approximation allows optimal inclusion of a priori information into retrieval procedure, thus ensuring the most effective regularization of the problem. Efficiency of the proposed technique is demonstrated on an example of retrieval of vertical ozone profile from data of ground-based sounding of the atmosphere in the millimeter wavelength range. Results of profile retrieval from model data and from spectra of radiation temperature of the atmosphere measured in the Apatity (67° N, 33° E) in the winter of 2002–2003 are presented.     

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.     

Diurnal variations of ion flux intensity on a synchronous orbit

The data from the synchronous-orbit satellites of the Gorizont series are used to study the dependences of the ion flux variation amplitudes in the synchronous altitude region (the diurnal behaviour) on particle energies and on the form and rigidity of the particle energy spectrum. The proton fluxes were measured in the energy range E 60–120 keV, and the [N,0]²⁺ and [C,N,0]⁴⁺ ion fluxes in the energy range E 60–70 keV/e.

The ratio of the diurnal variation amplitudes of the studied ions is shown to correspond to the similarity of their energy spectra in the E/Q representation. The magnetically-quiet time gradient of the distribution function F(μ,J,L) in the synchronous-orbit region is shown to be (∂F/∂L)=0 for the H⁺ and [N,0]²⁺ ions and (∂F/∂L) > 0 for the [C,N,0]⁴⁺ ions (at the values of μ corresponding to the examined energy ranges). During magnetically-disturbed periods the inner boundary of the (∂F/∂L)=0 region shifts to lower L and (∂ F/∂L) = O in the synchronous altitude region must be also for the [C,N,O]⁴⁺ ions.     

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.     

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.     

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.     

Diurnal and seasonal variations of F2 layer characteristics over Irkutsk during the decrease in solar activity in 2003–2006: Observations and IRI-2001 model predictions

The diurnal and seasonal variations of F2 layer characteristics (critical frequency, peak height and bottomside thickness) over Irkutsk, Russia (52.3 N and 104.3 E) are studied by the method of running medians. The comparison with the IRI-2001 model during the decrease in solar activity in 2003–2006 revealed cases of both close agreement and systematic differences between predictions and observations. The systematic difference is not the only reason for disagreement between IRI and observations; there are also intrayear variations which are not associated with seasonal behavior. The period of observation was too short to make conclusions about solar activity dependence of the noon bottomside thickness and the modification of its diurnal behavior with decreasing solar activity.     

Regular features of the polar ionosphere characteristics from Digisonde measurements over Norilsk

Regular features of the polar ionosphere have been studied using its local empirical model of the electron density distribution in the bottomside ionosphere. The local empirical model was derived from the hand-scaled ionogram data recorded by DPS-4 Digisonde at Norilsk, Russia (69.4N, 88.1E; 60N GLAT, 166E GLON) for a 6-year period from December, 2002 to December, 2008. The paper describes the technique used to build the local empirical model and discusses its diurnal, seasonal, and solar activity specifications in comparison with the standard IRI-2007 climatological model for the same period of time, long-term observations from the European Incoherent Scatter UHF radar (1988–1999), and the high-latitude ionosondes data. Primary focus of the paper is behavior of the three F2 layer parameters: the F2 peak density (NmF2), the peak height (hmF2) and the bottomside thickness (B0). Special emphasis of the paper is the analysis of the winter anomaly manifestation at Norilsk and the peculiar diurnal–seasonal behavior of hmF2 under low solar activity, named as a “polar day effect”.     

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.