Regional 4-D modeling of the ionospheric electron density

The knowledge of the electron density is the key point in correcting ionospheric delays of electromagnetic measurements and in studying the ionosphere. During the last decade GNSS, in particular GPS, has become a promising tool for monitoring the total electron content (TEC), i.e., the integral of the electron density along the ray-path between the transmitting satellite and the receiver. Hence, geometry-free GNSS measurements provide informations on the electron density, which is basically a four-dimensional function depending on spatial position and time. In addition, these GNSS measurements can be combined with other available data including nadir, over-ocean TEC observations from dual-frequency radar altimetry (T/P, JASON, ENVISAT), and TECs from GPS-LEO occultation systems (e.g., FORMOSAT-3/COSMIC, CHAMP) with heterogeneous sampling and accuracy.     

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

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

Local ionospheric electron density profile reconstruction in real time from simultaneous ground-based GNSS and ionosonde measurements

The purpose of the LIEDR (local ionospheric electron density profile reconstruction) system is to acquire and process data from simultaneous ground-based total electron content (TEC) and digital ionosonde measurements, and subsequently to deduce the vertical electron density distribution above the ionosonde’s location. LIEDR is primarily designed to operate in real time for service applications and, for research applications and further development of the system, in a post-processing mode. The system is suitable for use at sites where collocated TEC and digital ionosonde measurements are available. Developments, implementations, and some preliminary results are presented and discussed in view of possible applications.     

Variability of the ionospheric electron density at fixed heights and validation of IRI-2007 profile’s prediction at Ilorin

A study on the variability of the equatorial ionospheric electron density was carried out at fixed heights below the F2 peak using one month data for each of high and low solar activity periods. The data used for this study were obtained from ionograms recorded at Ilorin, Nigeria, and the study covers height range from 100 km to the peak of the F2 layer for the daytime hours and height range from 200 km to the peak of the F2 layer for the nighttime hours. The results showed that the deviation of the electron density variation from simple Chapman variation begins from an altitude of about 200 km for the two months investigated. Daytime minimum variability of between 2.7% and 9.0% was observed at the height range of about 160 and 200 km during low solar activity (January 2006) and between 3.7% and 7.8% at the height range of 210 and 260 km during high solar activity (January 2002). The nighttime maximum variability was observed at the height range of 210 and 240 km at low solar activity and at the height range of 200 and 240 km at high solar activity. A validation of IRI-2007 model electron density profile’s prediction was also carried out. The results showed that B0 option gives a better prediction around the noontime.     

The subauroral electron density trough: Comparison between satellite observations and IRI-2007 model estimates

We compare electron density predictions of the International Reference Ionosphere (IRI-2007) model with in situ measurements of the satellites CHAMP and GRACE for the years from 2005 to 2010 over the subauroral regions. The electron density between 58° and 68° Mlat are considered. The trough region Ne peaks during local summers and attain the valley during local winter. Around −100°E and 60°E, two larger electron density sectors features can be seen in both hemispheres during all three seasons, which attributed to the electron extending from middle latitude to trough region. From 2005 to the beginning of 2010, the model overestimates the trough region Ne by 20% on average and the decrease of Ne in this region can also be seen during the last solar minimum. In the southern hemisphere, the model prediction shows quite well consistence with the observation during all three seasons while the huge difference between observations and model estimation implies that the IRI-2007 model needs significant improvement to predict better the trough region in northern hemisphere.     

Electron density inversed by plasma lines induced by suprathermal electron in the ionospheric modification experiment

Incoherent scatter radar (ISR) is the most powerful ground-based measurement facility to study the ionosphere. The plasma lines are not routinely detected by the incoherent scatter radar due to the low intensity, which falls below the measured spectral noise level of the incoherent scatter radar. The plasma lines are occasionally enhanced by suprathermal electrons through the Landau damping process and detectable to the incoherent scatter radar. In this study, by using the European Incoherent Scatter Association (EISCAT) UHF incoherent scatter radar, the experiment observation presents that the enhanced plasma lines were observed. These plasma lines were considered as manifest of the suprathermal electrons generated by the high-frequency heating wave during the ionospheric modification. The electron density profile is also obtained from the enhanced plasma lines. This study can be a promising technique for obtaining the accurate electron density during ionospheric modification experiment.     

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.     

An empirical model of electron density in low latitude at 600 km obtained by Hinotori satellite

An empirical model of electron density (Ne) was constructed by using the data obtained with an impedance probe on board Japanese Hinotori satellite. The satellite was in circular orbit of the height of 600 km with the inclination of 31 degrees from February 1981 to June 1982. The constructed model gives Ne at any local time with the time resolution of 90 min and between −25 and 25 degrees in magnetic latitude with its resolution of 5 degrees in the range of F_10.7 from 150 to 250 under the condition of K_p < 4. Spline interpolations are applied to the functions of day of year, geomagnetic latitude and solar local time, and linear interpolation is applied to the function of F_10.7. Longitude dependence of Ne is not taken into account. Our density model can reproduce solar local time variation of electron density at 600 km altitude better than current International Reference Ionosphere (IRI2001) model which overestimates Ne in night time and underestimates Ne in day time. Our density model together with electron temperature model which has been constructed before will enable more understanding of upper ionospheric phenomenon in the equatorial region.     

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.     

Obtaining more accurate electron density profiles from bending angle with GPS occultation data: FORMOSAT-3/COSMIC constellation

Since 1995, with the first GPS occultation mission on board Low Earth Orbiter (LEO) GPS/MET, inversion techniques were being applied to GPS occultation data to retrieve accurate worldwide distributed refractivity profiles, i.e. electron density profiles in the case of Ionosphere. Important points to guarantee the accuracy is to take into account horizontal gradients and topside electron content above the LEO orbit. This allows improving the accuracy from 20% to 50%, depending on the conditions, latitude and epoch regarding to Solar cycle as reported in previous works.     

Evaluation of the IRI-2007 model options for the topside electron density

The International Reference Ionosphere (IRI) 2007 provides two new options for the topside electron density profile: (a) a correction of the IRI-2001 model, and (b) the NeQuick topside formula. We use the large volume of Alouette 1, 2 and ISIS 1, 2 topside sounder data to evaluate these two new options with special emphasis on the uppermost topside where IRI-2001 showed the largest discrepancies. We will also study the accurate representation of profiles in the equatorial anomaly region where the profile function has to accommodate two latitudinal maxima (crests) at lower altitudes but only a single maximum (at the equator) higher up. In addition to IRI-2001 and the two new IRI-2007 options we also include the Intercosmos-based topside model of Triskova, Truhlik, and Smilauer [Triskova, L., Truhlik, V., Smilauer, J. An empirical topside electron density model for calculation of absolute ion densities in IRI. Adv. Space Res. 37 (5), 928–934, 2006] (TTS model) in our analysis. We find that overall IRI-2007-NeQ gives the best results but IRI-2007-corrected provides a more realistic representation of the altitudinal–latitudinal structure in the equatorial anomaly region. The applicability of the TTS model is limited by the fact that it is not normalized to the F2 peak density and height.     

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.     

Evaluating the IRI topside model for the South African region: An overview of the modelling techniques

The representation of the topside ionosphere (the region above the F2 peak) is critical because of the limited experimental data available. Over the years, a wide range of models have been developed in an effort to represent the behaviour and the shape of the electron density (Ne) profile of the topside ionosphere. Various studies have been centred around calculating the vertical scale height (VSH) and have included (a) obtaining VSH from Global Positioning System (GPS) derived total electron content (TEC), (b) calculating the VSH from ground-based ionosonde measurements, (c) using topside sounder vertical Ne profiles to obtain the VSH. One or a combination of the topside profilers (Chapman function, exponential function, sech-squared (Epstein) function, and/or parabolic function) is then used to reconstruct the topside Ne profile. The different approaches and the modelling techniques are discussed with a view to identifying the most adequate approach to apply to the South African region’s topside modelling efforts. The IRI-2001 topside model is evaluated based on how well it reproduces measured topside profiles over the South African region. This study is a first step in the process of developing a South African topside ionosphere model.     

Evolution of artificial ion clouds in the ionospheric F-region in CRRES satellite experiments from radiophysical sounding data

In the paper the results of the G-09 and G-12 Caribbean CRRES experiments radiophysical probing of artificial barium ion clouds in the ionosphere carried by ionosonde and vertical Doppler sounder are considered.     

Variation of ionospheric electron and ion temperatures during periods of minimum to maximum solar activity by the SROSS-C2 satellite over Indian low and equatorial latitudes

To study the variation of ionospheric electron and ion temperatures with solar activity the data of electron and ion temperatures were recorded with the help of Retarding Potential Analyzer payload aboard Indian SROSS-C2 satellite at an average altitude of ∼500 km. The main focuses of the paper is to see the diurnal, seasonal and latitudinal variations of electron and ion temperatures during periods of minimum to maximum solar activity. The ionospheric temperatures in the topside show strong variations with altitude, latitude, season and solar activity. In present study, the temperature variations with latitude, season and solar activity have been studied at an average altitude ∼500 km. The peak at sunrise has been observed during all seasons, in both electron and ion temperatures. Further, the ionospheric temperatures vary with latitude in day time. The latitudinal variation is more pronounced for low solar activity than for high solar activity.     

Dependence of the F-region peak electron density (foF2) on solar activity observed in the equatorial ionospheric anomaly region in the Brazilian sector

The long-term (solar cycle) changes in the Sun and how it affects the ionospheric F-region observed at São José dos Campos (23.2° S, 45.9° W), Brazil, a location under the southern crest of the equatorial ionospheric anomaly, have been investigated in this paper. The dependence of the F-region peak electron density (foF2) on solar activity during the descending phase of the 23rd solar cycle for the periods of high, medium, and low solar activity has been studied. The ionospheric F-region peak electron densities observed during high and medium solar activity show seasonal variations with maxima close to the equinox periods, whereas during the low solar activity the maxima during the equinox periods is absent. However, during the low solar activity only change observed is a large decrease from summer to winter months. We have further investigated changes in the different ionospheric F-region parameters (minimum virtual height of the F-region (h′F), virtual height at 0.834foF2 (hpF2), and foF2) during summer to winter months in low solar activity periods, 2006–2007 and 2007–2008. Large changes in the two ionospheric parameters (hpF2 and foF2) are observed during summer to winter months in the two low solar activity periods investigated.     

Time series modeling and analysis of trends of daily averaged ionospheric total electron content

A 10.7 cm solar radio flux F10.7, geomagnetic planetary equivalent amplitude (Ap index), and period variations were considered in this paper to construct a linear model for daily averaged ionospheric total electron content (TEC). The correlation coefficient of the modeled results and International GNSS Service (IGS) observables was approximately 0.97, which implied that the model could accurately reflect the realistic variation characteristics of the daily averaged TEC. The influences of the different factors on TEC and its characteristics at different latitudes were examined with this model. Results show that solar activity, annual and semiannual cycles are the three most important factors that affect daily averaged TEC. Solar activity is the primary determinant of TEC during periods with high solar activity, whereas periodic factors primarily contribute to TEC during periods with minimum solar activity. The extent of the influences of the different factors on TEC exhibits obvious differences at varying latitudes. The magnitude of the semiannual variation becomes less significant with the increase in latitude. Furthermore, a geomagnetic storm causes an increase in TEC at low latitudes and a decrease at high latitudes.     

Ionospheric electron density over Resolute Bay according to E-CHAIM model and RISR radar measurements

In this study, predictions of the E-CHAIM ionospheric model are compared with measurements by the incoherent scatter radars RISR at Resolute Bay, Canada, in the northern polar cap. Reasonable coverage was available for all seasons except winter for which no conclusions were drawn. It is shown that ratios of the model-to measured electron densities are close to unity in the central part of the F layer, around its peak. This is particularly evident for summer daytime. Distributions of the ratios are wider for other seasons indicating larger number of cases when the model underestimates or overestimates. E-CHAIM underestimates the electron density at ionospheric topside and bottomside by ～ 10–20 %. At the bottomside, the underestimations are strongest in summer and equinoctial nighttime. At the topside, the underestimations are strongest in autumn nighttime. Model overestimations are noticeable in the middle part of the F layer during dawn hours in autumn. Overall, the model tends to not predict highest-observed peak electron densities and the largest-observed heights of the peak.     

Solar wind electron density and temperature over solar cycle 23: Thermal noise measurements on Wind

We present the solar wind plasma parameters obtained from the Wind spacecraft during more than nine years, encompassing almost the whole solar cycle 23. Since its launch in November 1994 Wind has frequently observed the in-ecliptic solar wind upstream of the Earth’s bow shock. The WIND/WAVES thermal noise receiver was specially designed to measure the in situ plasma thermal noise spectra, from which the electron density and temperature can be accurately determined. We present and discuss histograms of such measurements performed from 1994 to 2003. Using these large data sets, we study the density and core temperature variations with solar activity cycle and with different regimes of the solar wind. We confirm the anticorrelation of the electron density with the sunspot number, and obtain a positive correlation of the core temperature, with the sunspot number.     

A 2-D comparison of ionospheric convection derived from SuperDARN and DMSP measurements

Ion drift vectors measured by the DMSP satellites are compared with plasma convection vectors obtained by the SuperDARN HF radars through the standard Map Potential algorithm of Ruohoniemi and Baker [Ruohoniemi, M., Baker, K.B. Large-scale imaging of high-latitude convection with super dual auroral radar network HF radar observations. J. Geophys. Res. 103, 20797–20811, 1998]. Despite significant data spread, the agreement can be qualified as reasonable for a data set comprising of 149 satellite passes over the Northern Hemisphere at high latitudes. The slope of the best-fit line relating SuperDARN and DMSP velocity magnitudes is of the order of 0.3 with a tendency for the SuperDARN velocities to be smaller. The agreement between the azimuths of the ion drift and convection is better with the slope of the best-fit line being close to 1. It is shown that consistency between the radar and satellite measurements is much better if the SuperDARN line-of-sight velocities are compared with the DMSP cross-track ion drifts for events showing slow spatial and temporal variations of the convection. If areas of strong convection changes are included into comparison, the degree of agreement deteriorates drastically. This result implies that differences in the spatial and temporal resolutions of DMSP and SuperDARN measurements are crucial factors contributing to the observed discrepancies. In addition, some differences are introduced when the SuperDARN line-of-sight velocities are filtered and reprocessed into vectors with the application of a background convection model.