Comparative study of high-latitude,mid-latitude and low-latitude ionosphere on basis of local empirical models

The analysis of the regular features of the high-, mid- and low-latitude ionosphere characteristics has been carried out using local empirical models. The local empirical models were derived from the manual scaled ionogram data recorded by DPS-4 Digisondes located at Norilsk (69 N, 88E), Irkutsk (52 N, 104E) and Hainan (19 N, 109E) for a 6-year period from December, 2002 to December, 2008. The technique used to build the local empirical model is described. Primary focuses are diurnal, seasonal and solar cycle variations of the peak electron density and the peak height under low solar activity and their changes with increasing solar activity. The main objective of the paper is to reveal both common and specific features of high-, mid- and low-latitude ionosphere. Based on earlier comparisons with the International Reference Ionosphere model, we analyze how the common and specific features are reproduced by this model.     

Comparison of Diurnal,Seasonal and Solar Cycle Variations of High-latitude,Mid-latitude and Low-latitude Ionosphere

Comparison of regular (diurnal, seasonal and solar cycle) variations of high-latitude,mid-latitude and low-latitude ionospheric characteristics has been provided on basis of local empirical models of the peak electron density and the peak height. The local empirical models were derived from the hand-scaled ionogram data recorded by DPS-4 digisondes located at Norilsk (69°N, 88°E), Irkutsk (52°N, 104°E) and Hainan (19°N, 109°E) for a 6-year period from December, 2002 to December, 2008. The technique used to build the local empirical model is described. The primary focus is diurnal-seasonal behavior under low solar activity and its change with increasing solar activity. Both common and specific features of the high-latitude (Norilsk), mid-latitude (Irkutsk) and low-latitude (Hainan) regular variations were revealed using their local empirical models.      

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.     

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.     

Ionospheric variation at Thailand equatorial latitude station: Comparison between observations and IRI-2001 model predictions

In this paper, the F2-layer critical frequency (foF2) and peak height (hmF2) measured by the FM/CW ionosonde at Thailand equatorial latitude station, namely Chumphon (10.72°N, 99.37°E, dip 3.22) are presented. The measurement data during low solar activity from January 2004 to December 2006 are analyzed based on the diurnal, seasonal variation. The results are then compared with IRI-2001 model predictions. Our study shows that: (1) In general, both the URSI and CCIR options of the IRI model give foF2 close to the measured ones, but the CCIR option produces a smaller range of deviation than the URSI option. The agreement during daytime is generally better than during nighttime. Overestimation mostly occurs in 2004 and 2006, while underestimation is during pre-sunrise hours in June solstice in 2005. The peak foF2 around sunset is higher during March equinox and September equinox than the other seasons, with longer duration of maximum levels in March equinox than September equinox. Large coefficients of variability foF2 occur during pre-sunrise hours. Meanwhile, the best agreement between the observed foF2 and the IRI model is obtained in June solstice. (2) In general, The IRI (CCIR) model predicts the observed hmF2 well during daytime in June solstice from 2004–2006, but it overestimates during March equinox, September equinox and December solstice. For nighttime, the model overestimates hmF2 values for all seasons especially during March equinox and September equinox. However, the model underestimates hmF2 values during September equinox and for some cases during June solstice and December solstice at pre-sunrise. The agreement between the IRI model and the hmF2(M3000_OBS) is worst around noontime, post-sunset and pre-sunrise hours. All comparative studies give feedback for new improvements of CCIR and URSI IRI models.     

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.     

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.     

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 between ionospheric peak parameters retrieved from COSMIC measurement and ionosonde observation over Sanya

In this paper we compared the ionospheric peak parameters (peak electron density of the F region, NmF2, and peak height of the F region, hmF2) retrieved from the FORMASAT-3/COSMIC (COSMIC for short) satellite measurement with those from ionosonde observation at Sanya (18.3°N, 109.6°E) during the period of 2008–2013. Although COSMIC NmF2 (hmF2) tends to be lower (higher) than ionosonde NmF2 (hmF2), the results show that the ionospheric peak parameters retrieved from COSMIC measurement generally agree well with ionosonde observation. For NmF2 the correlation between the COSMIC measurement and the ionosonde observation is higher than 0.89, and for hmF2 the correlation is higher than 0.80. The correlation of the ionospheric peak parameters decreases when solar activity increases. The performance of COSMIC measurement is acceptable under geomagnetic disturbed condition. The correlation of NmF2 between COSMIC and ionosonde measurements is higher (lower) during the nighttime (daytime), while the correlation of hmF2 is lower (higher) during the nighttime (daytime).     

Diurnal and seasonal variation of F2-layer ionospheric parameters at equatorial ionization anomaly crest region and their comparison with IRI-2001

Diurnal and seasonal variations of critical frequency of ionospheric F2-region ‘foF2’ and the height of peak density ‘hmF2’ are studied using modern digital ionosonde observations of equatorial ionization anomaly (EIA) crest region, Bhopal (23.2°N, 77.6°E, dip 18.5°N), during solar minimum period 2007. Median values of these parameters are obtained at each hour using manually scaled data during different seasons and compared with the International Reference Ionosphere-2001 model predictions. The observations suggest that on seasonal basis, the highest values of foF2 are observed during equinox months, whereas highest values of hmF2 are obtained in summer and lowest values of both foF2 and hmF2 are observed during winter. The observed median and IRI predicted values of foF2 and hmF2 are analyzed with upper and lower bound of inter-quartile range (IQR) and it is find out that the observed median values are well inside the inter-quartile range during the period of 2007. Comparison of the recorded foF2 and hmF2 values with the IRI-2001 output reveals that IRI predicted values exhibit better agreement with hmF2 as compared to foF2. In general, the IRI model predictions show some agreement with the observations during the year 2007. Therefore it is still necessary to implement improvements in order to obtain better predictions for EIA regions.     

The variability and IRI2007-predictability of hmF2 over South Africa

This paper presents an investigation into the variability and predictability of the maximum height of the ionospheric F2 layer, hmF2 over the South African region. Data from three South African stations, namely Madimbo (22.4°S, 26.5°E, dip angle: −61.47°), Grahamstown (33.3°S, 26.5°E, dip angle: −64.08°) and Louisvale (28.5°S, 21.2°E, dip angle: −65.44°) were used in this study. The results indicate that hmF2 shows a larger variability around midnight than during the daytime for all seasons. Monthly median hmF2 values were used in all cases and were compared with predictions from the IRI-2007 model, using the URSI (Union Radio-Scientifique Internationale) coefficient option. The analysis covers the diurnal and seasonal hourly hmF2 values for the selected months and time sectors e.g. January, July, April and October for 2003 and 2005. The time ranges between (03h00–23h00 UT; LT = UT + 2h) representing the local sunrise, midday, sunset and midnight hours. The time covers sunrise, midday, sunrise, and midnight hours (03–06h00 UT, 07–11h00 UT, sunrise 16–18h00 UT and 22–23h00 UT; LT = UT + 2h). The dependence of the results on solar activity levels was also investigated. The IRI-2007 predictions follow fairly well the diurnal and seasonal variation patterns of the observed hmF2 values at all the stations. However, the IRI-2007 model overestimates and underestimates the hmF2 value during different months for all the solar activity periods.     

Ionospheric variability at taiwan low latitude station: Comparison between observations and IRI-2001 model

This paper presents the observed ionospheric F-region critical frequency, foF2, and peak height, hmF2, at northern crest of equatorial ionization anomaly (EIA) area station, namely Chung-Li (24.9°N, 121.1°E, dip 35°), and to be compared with International Reference Ionosphere model (IRI-2001) predictions for the period from 1994 to 1999, corresponding to half of the 23rd solar cycle. The diurnal and seasonal variation of foF2 and hmF2 are analyzed for different solar phases, respectively. The result shows the largest discrepancies were observed during nighttime for foF2 and hmF2, respectively. The value of foF2 both CCIR and URSI selected in the IRI model produced a good agreement during the daytime and underestimated during the noon time for high solar activities. The underestimation at noon time is mainly caused by the fountain effect from equator. Further, the peak height hmF2 shows a larger variability around the midnight than daytime in the equinox and winter seasons and reserved in summer, respectively. The study shows that the monthly median values of observed hmF2 is somewhat lower than those predicated by the IRI model, at night time in all the seasons except the period of 04:00–06:00 LT and reverse at daytime in summer. In general the IRI model predictions with respect to the observed in hmF2 is much better than foF2. The percentage deviation of the observed foF2 (hmF2) values with respect to the IRI model varies from 5% to 80% (0–25%) during nighttime and 2–17% (0–20%) at daytime, respectively. In general, the model generates good results, although some improvements are still necessary to implement in order to obtain better simulations for ionospheric low-latitudes region.     

Linkage of the ionospheric peak electron density and height deduced from the topside sounding data

Topside sounding electron density profiles are analyzed to explore interrelations of the F2 layer critical frequency and the peak height for a representative set of conditions provided by ISIS1, ISIS2, IK19 and Cosmos-1809 satellites for the period of 1969–1987. The foF2 and hmF2 are delivered with exponential extrapolation of electron density profile to zero of its 1st derivative. It is shown that the linear regression exists between foF2 and hmF2 under different conditions. The linkage between the two parameters amended to the empirical model of the peak height [Gulyaeva, T.L., Bradley, P.A., Stanislawska, I., Juchnikowski, G. Towards a new reference model of hmF2 for IRI. Adv. Space Res. 42, 666–672, doi:10.1016/j.asr.2008.02.021, 2008] results in an empirical model of the both foF2 and hmF2 expressed by superposition of functions in terms of local-time, season, geodetic longitude, modified dip latitude and solar activity. For the solar activity we use a proxy Fsp index averaged from the mean solar radio flux F10.7s for the past 81 days (3 solar rotations) and F10.7 value for 1 day prior the day of observation. Impact of geomagnetic activity is not discernible with the topside sounding data due to mixed positive and negative storm-time effects. Appreciable differences have been revealed between IRI-CCIR predictions and outcome of the new model which might be attributed to the different techniques of the peak electron density and height derivation, different epochs and different global distribution of the source data as well as the different mathematical functions involved in the maps and the model presentation.     

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.     

Variability of the ionosphere over Irkutsk at low solar activity

Statistical and spectral analyses are performed to investigate variations of two ionosphere F2 layer key parameters, the critical frequency (foF2) and the peak height (hmF2), that were measured over Irkutsk (52.5°N, 104.0°E) from December 2006 to January 2008 under solar minimum. The analyses showed that both parameters contain quasi-harmonic oscillations with periods of T_n = 24/n hours (n = 1–7), among which the diurnal (n = 1) and semidiurnal (n = 2) ones are the strongest. Seasonal variations are explored of mean and median values, spectrum, amplitude, and phase of the diurnal and semidiurnal components of foF2 and hmF2.     

Towards a new reference model of hmF2 for IRI

A numerical model of the peak height of the F2 layer, hmF2_top, is derived from the topside sounding database of 90,000 electron density profiles for a representative set of conditions provided by ISIS1, ISIS2, IK19 and Cosmos-1809 satellites for the period of 1969–1987. The model of regular hmF2 variations is produced in terms of local time, season, geomagnetic latitude, geodetic longitude and solar radio flux. No geomagnetic activity trends were discernible in the topside sounding data. The nighttime peak of hmF2_top evident for mid-latitudes disappears near the geomagnetic equator where a maximum of hmF2_top occurs at sunset hours when it can exceed 500 km at solar maximum. The hmF2 given by the IRI exceeds hmF2_top at the low solar activities. The hmF2_top, obtained by extrapolation of the first derivative of the topside profile to zero shows saturation similar to foF2 the greater the solar activity. The proposed model differs from hmF2 given by IRI based on M(3000)F2 to hmF2 conversion by empirical relationships in terms of foF2, foE and R12 with these quantities mapped globally by the ITU-R (former CCIR) from ground-based ionosonde data. The differences can be attributed to the different techniques of the peak height derivation, different epochs and different global distribution of the source data as well as the different mathematical functions involved in the maps and the model presentation.     

Deformation of the ionosphere structure during the space weather events on October–November 2003

We present the spatial maps of the ionosphere–plasmasphere slab thickness τ (ratio of the vertical total electron content, TEC, to the F-region peak electron density, NmF2) during the intense ionospheric storms of October–November 2003. The model-assisted technology for estimate of the upper boundary of the ionosphere, hup, from the slab thickness components in the bottomside and topside ionosphere – eliminating the plasmasphere contribution of τ – is applied at latitudes 35° to 70°N and longitudes −10° to 40°E, from the data of 20 observatories of GPS-TEC and ionosonde networks, for selected days and hours of October and November 2003. The daily–hourly values of NmF2, hmF2 and TECgps are used as the constrained parameters for the International Reference Ionosphere extended to the plasmasphere, IRI-Plas, during the ionospheric quiet days, positive and negative storm phases for estimate of τ and hup. Good correlation has been found between the slab thickness and the upper boundary of the ionosphere for the intense ionospheric storms at October–November 2003. During the negative phase of the ionospheric storm, when the ionospheric plasma density is exhausted, the nighttime upper boundary of the ionosphere is greatly uplifted towards the magnetosphere tail, while the daytime upper boundary of the ionosphere is reduced below 500 km over the Earth.     

Longitudinal behaviors of the IRI-B parameters of the equatorial electron density profiles retrieved from FORMOSAT-3/COSMIC radio occultation measurements

The electron density profiles in the bottomside F₂-layer ionosphere are described by the thickness parameter B0 and the shape parameter B1 in the International Reference Ionosphere (IRI) model. We collected the ionospheric electron density (N_e) profiles from the FORMOSAT-3/COSMIC (F3/C) radio occultation measurements from DoY (day number of year) 194, 2006 to DoY 293, 2008 to investigate the daytime behaviors of IRI-B parameters (B0 and B1) in the equatorial regions. Our fittings confirm that the IRI bottomside profile function can well describe the averaged profiles in the bottomside ionosphere. Analysis of the equatorial electron density profile datasets provides unprecedented detail of the behaviors of B0 and B1 parameters in equatorial regions at low solar activity. The longitudinal averaged B1 has values comparable with IRI-2007 while it shows little seasonal variation. In contrast, the observed B0 presents semiannual variation with maxima in solstice months and minima in equinox months, which is not reproduced by IRI-2007. Moreover, there are complicated longitudinal variations of B0 with patterns varying with seasons. Peaks are distinct in the wave-like longitudinal structure of B0 in equinox months. An outstanding feature is that a stable peak appears around 100°E in four seasons. The significant longitudinal variation of B0 provides challenges for further improving the presentations of the bottomside ionosphere in IRI.     

An analysis of the topside ionosphere parameters based on the long-duration Irkutsk incoherent scatter radar measurements

The topside ionosphere parameters are studied based on the long-duration Irkutsk incoherent scatter radar (52.9N, 103.3E) measurements conducted in September 2005, June and December 2007. As a topside ionosphere parameter we chose the vertical scale height (VSH) related to the gradient of the electron density logarithm above the peak height. For morphological studies we used median electron density profiles. Besides the median behavior we also studied VSH disturbances (deviations from median values) during the magnetic storm of September 11th 2005. We compared the Irkutsk incoherent scatter radar data with the Millstone Hill and Arecibo incoherent scatter radar observations, the IRI-2007 prediction (using the two topside options) and VSH derived from the Irkutsk DPS-4 Digisonde bottomside measurements.     

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.