Modelling bottom and topside electron density and TEC with profile data from topside ionograms

The paper describes the technique that has been implemented to model the electron density distribution above and below the F2 peak making use of only the profiles obtained from the INTERCOSMOS-19 topside ionograms. Each single profile from the satellite height to the ionosphere peak has been fitted by a semi-Epstein layer function of the type used in the DGR model with shape factor variable with altitude. The topside above the satellite height has been extrapolated to match given values of plasmaspheric electron densities to obtain the full topside profile. The bottomside electron density has been calculated by using the maximum electron density and its altitude estimated from the topside ionogram as input for a modified version of the DGR derived profiler that uses model values for the foF1 and foE layers of the ionosphere. Total electron content has also been calculated. Longitudinal cross sections of vertical profiles from latitudes 50° N to 50° S latitude are shown for low and high geomagnetic activity. These cross sections indicate the equatorial anomaly effect and the changes of the shape of low latitude topside ionosphere during geomagnetic active periods. These results and the potentiality of the technique are discussed.     

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

Ionospheric electron density profiles at sunrise-sunset

In order to improve its representation of the dependence on time and space of the ionospheric parameters, the International Reference Ionosphere ought to take account of realistic sunrise and sunset conditions in the upper atmosphere. Such input is needed for quite a few parameters for which only day and night values were taken as input in the present IRI. Of the 24 hours of a day, true nighttime comprises a fraction of 37% at an altitude of 300 km and only 26% at 1000 km. In order to demarcate the day/night/day transition periods, the present IRI proposes solar zenith angles of 98° to 120°, depending on the altitude.Electron density profiles, obtained during these periods, have been studied with two data sources: 10 vertical-incidence sounding data observed during the meridional voyages of the research vessel “Akademik Korolev” in the Pacific Ocean; 2° data observed at the South Pole. It is shown that the height of the turning point in the sub-peak F2-layer profile and also the corresponding minimum scale height appear to be independent of latitude, season and index of geomagnetic activity. A method is discussed by which the IRI electron density profiles might be improved, in particular during these hours.     

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.     

Observation of the north-south asymmetry in the equatorial spread F intensity at the Brazilian longitude

This study presents an analysis of the observed north-south asymmetry of the range spread F (RSF) intensity at the low latitude region during an equinoctial month of different solar epochs (2002, 2015 and 2017). The ionospheric parameters were obtained during geomagnetic quiet days from four digisonde stations located along the Brazilian longitude, which include a dip equator station (Sao Luiz (SL: 2.33 S, 44.2 W)), conjugate stations (Campo Grande (CG: 20.5°S, 55°W) and Boa Vista (BV: 2.8°N, 60.7°W)) and another low latitude station (Cachoeira Paulista (CP: 22.7°S, 45°W)). The results highlight the competing effect of the post-sunset electric field strength and the trans-equatorial wind on the latitudinal distribution of the irregularity intensity at both hemispheres under varying background ionospheric condition. The RSF intensity was seen to reduce as the solar flux index decreased and the latitudinal peak shifted closer to the dip equator. This was dependent on the variation of the field line mapped irregularity spectrum and the density gradient. Likewise, the north-south asymmetry in the irregularity occurrence was seen to become more significant as a denser ionosphere was observed at the hemisphere with the equatorward meridional wind. This has further proven that the non-linear cascading of the plasma irregularity across the low latitude region is strongly influenced by the local electric field.     

Optical studies of the ionospheric irregularities over the brazilian region by nocturnal images of the OI 630 nm emission

This study is an extension of previous statistical studies (Sobral et al., 1990, 1991, 1999) of both the local time and latitude variations of the zonal drift velocities of ionospheric plasma depletions, over the Brazilian low latitude station Cachoeira Paulista — CP (22.54°S, 45.00°W). The past studies were based on OI 630 nm scanning photometer data and the present one is based on digital OI 630 nm airglow images obtained by an all-sky imager system. These data were gathered between October 1998 and October 1999, at CP. The present results show that, in general, the velocities clearly tended to decrease with local time. Such a decrease should be associated with decreasing intensity of the vertical component of the ambient electric field which, in turn can be accounted for by recombination. All zonal drifts obtained for the 18 nights were eastwards. During equinox, the velocities clearly tended to decrease with local time at lower rates as compared with spring and summer. The highest and lowest zonal drift velocities, from all three seasons considered here, were observed to be in the summer ≈180 ms⁻¹ at 21:45 LT, and in the spring ≈25 ms⁻¹ at 03:15 LT, respectively. Ionospheric plasma bubbles were detected out to the maximum extra-tropical geographical latitude of ≈28° S, which was the highest latitude position analyzed in this study.     

Thermal net flux measurements on the pioneer Venus entry probes

Corrected thermal net radiation measurements from the four Pioneer Venus entry probes at latitudes of 60°N, 31°S, 27°S, and 4°N are presented. Three main conclusions can be drawn from comparisons of the corrected fluxes with radiative transfer calculations: (1) sounder probe net fluxes are consistent with the number density of large cloud particles (mode 3) measured on the same probe, but the IR measurements as a whole are most consistent with a significantly reduced mode 3 contribution to the cloud opacity; (2) at all probe sites, the fluxes imply that the upper cloud contains a yet undetected source of IR opacity; and (3) beneath the clouds the fluxes at a given altitude increase with latitude, suggesting greater IR cooling below the clouds at high latitudes and water vapor mixing ratios of about 2–5×10^?5 near 60°, 2–5×10^?4 near 30°, and >5×10^?4 near the equator.     

Rocket observations of the atomic and molecular oxygen emissions in the equatorial region

A Brazilian sounding rocket, SONDA III, with two airglow photometers and two ionospheric electron density probes, was launched successfully from Natal (5.8°S, 35.2°W), Brazil, on December 11, 1985, at 23:30 GMT. The observed height profiles of the atomic oxygen OI 5577Å and molecular oxygen Atmospheric (0,0) band at 7619Å emissions are discussed. This is the first simultaneous measurement of these emissions in the equatorial region. A preliminary analysis shows that the two emissions have peak emission heights located between 95 and 96 km, and their half widths are about 6 km. The O₂A 7619Å emission peak, however, is located slightly lower, less than 1 km, than that of the OI 5577Å emission.     

Variation of hmF2 and NmF2 deduced from DPS-4 over Multan (Pakistan) and their comparisons with IRI-2012 & IRI-2016 during the deep solar minimum between cycles 23 & 24

We report the results of ionospheric measurements from DPS-4 installed at Multan (Geog coord. 30.18°N, 71.48°E, dip 47.4°). The variations in F2-layer maximum electron density NmF2 and its peak height hmF2 are studied during the deep solar minimum between cycles 23 & 24 i.e 2008–2009 with comparisons conducted with the International Reference Ionosphere (IRI) versions 2012 & 2016. We find that the hmF2 observations peak around the pre-sunrise and sunrise hours depending on the month. Seasonally, the daytime variation of NmF2 is higher in the Equinox and Summer, while daytime hmF2 are slightly higher in the Equinox and Winter. High values of hmF2 around midnight are caused by an increase of upward drifts produced by meridional winds. The ionosphere over Multan, which lies at the verge of low and mid latitude, is affected by both $E\times B$ drifts and thermospheric winds as evident from mid-night peaks and near-sunrise dips in hmF2. The results of the comparison of the observed NmF2 and hmF2 for the year 2008–2009 with the IRI-2012 (both NmF2 and hmF2) and IRI-2016 (only hmF2) estimates indicate that for NmF2, IRI-2012 with Consultative Committee International Radio (CCIR) option produces values in better agreement with observed data. Whereas, for hmF2, IRI-2016 with both International Union of Radio Science (URSI) and CCIR SHU-2015 options, predicts well for nighttime hours throughout the year. However, the IRI-2012 with CCIR option produces better agreement with data during daytime hours. Furthermore, IRI-2012 with CCIR option gives better results during Equinox months, whereas, IRI-2016 with both URSI and CCIR SHU-2015 options predict well for Winter and Summer.     

Mesospheric ionisation over dip equator at sunrise

The electron density profile in the equatorial mesosphere was measured during sunrise time over Thumba(dip lat= 0.6°S). The measurements were carried out in the altitude range 60 to 100 km using rocketborne probes. A sharp layer of ionisation was observed around 80 km with electron density about 10⁸m^?3. It is suggested that hydrated ions are the main constituents of this layer.     

Statistical parameters of nonisothermal lower ionospheric plasma in the electrically active mesosphere

The large V/m electric fields inherent in the lower mesosphere play an essential role in lower ionospheric electrodynamics. They must be the cause of large variations in the electron temperature and the electron collision frequency and consequently of the transition of the ionospheric plasma in the lower part of the D region into a nonisothermal state. This study is based on the datasets on large mesospheric electric fields collected with the 2.2-MHz radar of the Institute of Space and Atmospheric Studies, University of Saskatchewan, Canada (52°N geographic latitude, 60.4°N geomagnetic latitude), and with the 2.3-MHz radar of the Kharkiv V. Karazin National University, Ukraine (49.6°N geographic latitude, 45.6°N geomagnetic latitude). The statistical analysis of these data is presented by [Meek, C.E., Manson, A.H., Martynenko, S.I., Rozumenko, V.T., Tyrnov, O.F. Remote sensing of mesospheric electric fields using MF radars. J. Atmos. Solar-Terr. Phys. 66, 881–890, 2004. 10.1016/j.jastp.2004.02.002]. The large mesospheric electric fields in the 60–67-km altitude range are experimentally established to follow a Rayleigh distribution in the 0 < E < 2.5 V/m interval. These data have permitted the resulting differential distributions of relative disturbances in the electron temperature, θ, and the effective electron collision frequency, η, to be determined. The most probable θ and η values are found to be in the 1.4–2.2 interval, and hence the nonstationary state of the lower part of the D region needs to be accounted for in studying processes coupling the electrically active mesosphere and the lower ionospheric plasma.     

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 comparative study of ionospheric IRIEup and ISP assimilative models during some intense and severe geomagnetic storms

Three-dimensional (3-D) electron density matrices, computed in the Mediterranean area by the IRI climatological model and IRIEup and ISP nowcasting models, during some intense and severe geomagnetic-ionospheric storms, were ingested by the ray tracing software tool IONORT, to synthesize quasi-vertical ionograms. IRIEup model was run in different operational modes: (1) assimilating validated autoscaled electron density profiles only from a limited area which, in our case, is the Mediterranean sector (IRIEup_re(V) mode); (2) assimilating electron density profiles from a larger region including several stations spread across Europe: (a) without taking care of validating the autoscaled data in the assimilation process (IRIEup(NV)); (b) validating carefully the autoscaled electron density profiles before their assimilation (IRIEup(V)).The comparative analysis was carried out comparing IRI, IRIEup_re(V), ISP, IRIEup(NV), and IRIEup(V) foF2 synthesized values, with corresponding foF2 measurements autoscaled by ARTIST, and then validated, at the truth sites of Roquetes (40.80°N, 0.50°E, Spain), San Vito (40.60°N, 17.80°E, Italy), Athens (38.00°N, 23.50°E, Greece), and Nicosia, (35.03°N, 33.16°E, Cyprus). The outcomes demonstrate that: (1) IRIEup_re(V), performs better than ISP in the western Mediterranean (around Roquetes); (2) ISP performs slightly better than IRIEup_re(V) in the central part of Mediterranean (around Athens and San Vito); (3) ISP performance is better than the IRIEup_re(V) one in the eastern Mediterranean (around Nicosia); (4) IRIEup(NV) performance is worse than the IRIEup(V) one; (5) in the central Mediterranean area, IRIEup(V) performance is better than the IRIEup_re(V) one, and it is practically the same for the western and eastern sectors.Concerning the overall performance, nowcasting models proved to be considerably more reliable than the climatological IRI model to represent the ionosphere behaviour during geomagnetic-ionospheric storm conditions; ISP and IRIEup(V) provided the best performance, but neither of them has clearly prevailed over the other one.     

Statistical ionospheric E layer properties measured with the Cyprus Digisonde and comparisons with IRI predictions

This paper reports the diurnal, seasonal, and long term variability of the E layer critical frequency (foE) and peak height (hmE) derived from Digisonde measurements from 2009 to 2016 at the low-middle latitude European station of Nicosia, Cyprus (geographical coordinates: 35°N, 33°E, geomagnetic lat. 29.38°N, I = 51.7°). Manually scaled monthly median values of foE and hmE are compared with IRI-2012 predictions with a view to assess the predictability of IRI. Results show that in general, IRI slightly overestimates foE values both at low and high solar activity. At low solar activity, overestimations are mostly limited to 0.25?MHz (equivalent electron density, 0.775?×?10³?el/m^?3) but can go as high as 0.5?MHz (equivalent electron density, 3.1?×?10³?el/m^?3, during noon) around equinox. In some months, underestimations, though sporadic in nature, up to 0.25?MHz are noted (mostly during sunrise and sunset). At high solar activity, a similar pattern of over-/underestimation is evident. During the entire period of study, over-/under estimations are mostly limited to 0.25?MHz. In very few cases, these exceed 0.25?MHz but are limited to 0.5?MHz. Analysis of hmE reveals that: (1) hmE remains almost constant during ±2 to ±4?h around local noon, (2) hmE values are higher in winter than in spring, summer and autumn, (3) there are two maxima near sunrise and sunset with a noontime minimum in between. During the entire period of study, significant differences between observed hmE and the IRI predictions have been noted. IRI fails to predict hmE and outputs a constant value of 110?km, which is higher than most of the observed values. Over- and under estimations range from 3 to 13?km and from 0 to 3?km respectively.     

Spread-F occurrence during geomagnetic storms near the southern crest of the EIA in Argentina

This work presents, for the first time, the analysis of the occurrence of ionospheric irregularities during geomagnetic storms at Tucumán, Argentina, a low latitude station in the Southern American longitudinal sector (26.9°S, 294.6°E; magnetic latitude 15.5°S) near the southern crest of the equatorial ionization anomaly (EIA). Three geomagnetic storms occurred on May 27, 2017 (a month of low occurrence rates of spread-F), October 12, 2016 (a month of transition from low to high occurrence rates of spread-F) and November 7, 2017 (a month of high occurrence rates of spread-F) are analyzed using Global Positioning System (GPS) receivers and ionosondes. The rate of change of total electron content (TEC) Index (ROTI), GPS Ionospheric L-band scintillation, the virtual height of the F-layer bottom side (h'F) and the critical frequency of the F2 layer (foF2) are considered. Furthermore, each ionogram is manually examined for the presence of spread-F signatures.The results show that, for the three events studied, geomagnetic activity creates favorable conditions for the initiation of ionospheric irregularities, manifested by ionogram spread-F and TEC fluctuation. Post-midnight irregularities may have occurred due to the presence of eastward disturbance dynamo electric fields (DDEF). For the May storm, an eastward over-shielding prompt penetration electric field, (PPEF) is also acting. A possibility is that the PPEF is added to the DDEF and produces the uplifting of the F region that helps trigger the irregularities. Finally, during October and November, strong GPS L band scintillation is observed associated with strong range spread-F (SSF), that is, irregularities extending from the bottom-side to the topside of the F region.     

Characterization of the low latitude plasma density irregularities observed using C/NOFS and SCINDA data

Complex electrodynamic processes over the low latitude region often result in post sunset plasma density irregularities which degrade satellite communication and navigation. In order to forecast the density irregularities, their occurrence time, duration and location need to be quantified. Data from the Communication/Navigation Outage Forecasting System (C/NOFS) satellite was used to characterize the low latitude ion density irregularities from 2011 to 2013. This was supported by ground based data from the SCIntillation Network Decision Aid (SCINDA) receivers at Makerere (Geographic coordinate 32.6°E, 0.3°N, and dip latitude ?9.3°N) and Nairobi (Geographic coordinate 36.8°E, ?1.3°N, and dip latitude ?10.8°N). The results show that irregularities in ion density have a daily pattern with peaks from 20:00 to 24:00 Local Time (LT). Scintillation activity at L band and VHF over East Africa peaked in 2011 and 2012 from 20:00 to 24:00 LT, though in many cases scintillation at VHF persisted longer than that at L band. A longitudinal pattern in ion density irregularity occurrence was observed with peaks over 135–180°E and 270–300°E. The likelihood of ion density irregularity occurrence decreased with increasing altitude. Analysis of C/NOFS zonal ion drift velocities showed that the largest nighttime and daytime drifts were in 270–300°E and 300–330°E longitude regions respectively. Zonal irregularity drift velocities over East Africa were for the first time estimated from L-band scintillation indices. The results show that the velocity of plasma density irregularities in 2011 and 2012 varied daily, and hourly in the range of 50–150 m s^?1. The zonal drift velocity estimates from the L-band scintillation indices had good positive correlation with the zonal drift velocities derived from VHF receivers by the spaced receiver technique.     

NeQuick 2 and IRI Plas VTEC predictions for low latitude and South American sector

Using vertical total electron content (VTEC) measurements obtained from GPS satellite signals the capability of the NeQuick 2 and IRI Plas models to predict VTEC over the low latitude and South American sector is analyzed. In the present work both models were used to calculate VTEC up to the height of GPS satellites. Also, comparisons between the performance of IRI Plas and IRI 2007 have been done. The data correspond to June solstice and September equinox 1999 (high solar activity) and they were obtained at nine stations. The considered latitude range extends from 18.4°N to ?64.7°N and the longitude ranges from 281.3°E to 295.9°E in the South American sector. The greatest discrepancies among model predictions and the measured VTEC are obtained at low latitudes stations placed in the equatorial anomaly region. Underestimations as strong as 40?TECU [1?TECU?=?10¹⁶?m^?2] can be observed at BOGT station for September equinox, when NeQuick2 model is used. The obtained results also show that: (a) for June solstice, in general the performance of IRI Plas for low latitude stations is better than that of NeQuick2 and, vice versa, for highest latitudes the performance of NeQuick2 is better than that of IRI Plas. For the stations TUCU and SANT both models have good performance; (b) for September equinox the performances of the models do not follow a clearly defined pattern as in the other season. However, it can be seen that for the region placed between the Northern peak and the valley of the equatorial anomaly, in general, the performance of IRI Plas is better than that of NeQuick2 for hours of maximum ionization. From TUCU to the South, the best TEC predictions are given by NeQuick2.The source of the observed deviations of the models has been explored in terms of CCIR foF2 determination in the available ionosonde stations in the region. Discrepancies can be also related to an unrealistic shape of the vertical electron density profile and or an erroneous prediction of the plasmaspheric contribution to the vertical total electron content. Moreover, the results of this study could be suggesting that in the case of NeQuick, the underestimation trend could be due to the lack of a proper plasmaspheric model in its topside representation. In contrast, the plasmaspheric model included in IRI, leads to clear overestimations of GPS derived TEC.     

Effect of the seed perturbation amplitude on the equatorial spread F initiation during solar minimum

The contribution of gravity wave (GW) to the initiation/development of spread F during a solar minimum year was investigated through the comparison of the observed precursory parameters and characteristics of the corresponding equatorial spread F (ESF) events. The ionospheric parameters were recorded at the magnetic equatorial station Sao Luis (2.3°S, 44°W, dip latitude 2°S) during March and October 2010. These data were used to estimate the influence of the relative gravity wave amplitude and the ambient ionospheric condition on the diurnal variation of the spread F initiation. The vertical velocity drift indicated a clear control and defines the threshold for the seasonal variability of the ESF occurrence. However, it was insufficient to solely determine or predict the day to day variation of ESF occurrence. Thus, few days with contrasting ambient ionospheric condition and magnitude of GW amplitude were analysed in order to investigate the role of the different precursory factors in the observed diurnal variation of the plasma irregularity development. The density scale length and gravity wave amplitude were shown to immensely contribute to the linear instability growth rate, especially during the days with a low post-sunset rise. Thus, the experimental observations have demonstrated the strong inter-dependence between the precursory factors and they have also highlighted the probable control of the ESF morphology.     

Effects of geomagnetic storms on the bottomside ionospheric F region

The geomagnetic storm is a complex process of solar wind/magnetospheric origin. The variability of the ionospheric parameters increases substantially during geomagnetic storms initiated by solar disturbances. Various features of geomagnetic storm act at various altitudes in the ionosphere and neutral atmosphere. The paper deals with variability of the electron density of the ionospheric bottomside F region at every 10 km of altitude during intense geomagnetic storms with attention paid mainly to the distribution of the F1 region daytime ionisation. We have analysed all available electron density profiles from some European middle latitude stations (Chilton, Pruhonice, Ebro, Arenosillo, Athens) for 36 events that occurred in different seasons and under different levels of solar activity (1995–2003). Selected events consist of both depletion and increase of the F2 region electron density. For European higher middle and middle latitude the F1 region response to geomagnetic storm was found to be negative (decrease of electron density) independent on the storm effect on the F2 region. For lower middle latitude the F1 response is weaker and less regular. Results of the analysis also show that the maximum of the storm effect may sometimes occur below the height of the maximum of electron density (NmF2).     

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.