Harmonic analysis of the ionospheric electron densities retrieved from FORMOSAT-3/COSMIC radio occultation measurements

In order to investigate the regular variations of the ionosphere, the least-squares harmonic estimation is applied to the time series of ionospheric electron densities derived from about five years of Global Positioning System radio occultation observations by FORMOSAT-3/COSMIC satellites. The analysis is done for different latitudes and altitudes in the region of Iran. The least-squares harmonic estimation is found to be a powerful tool for the frequency analysis of the completely unevenly spaced time series of radio occultation measurements. Although the obtained results are slightly different from the exact expected cycles (i.e. annual and diurnal components with their Fourier decompositions, and the 27-day period) due to the low horizontal resolution of radio occultation measurements, high vertical resolution of the observations enables us to detect not only the total electron content variations but also periodic patterns of electron densities at different altitudes of the ionosphere. The dominant diurnal and annual signals together with their Fourier series decompositions are obtained, which are consistent with the previous analyses on the total electron content. In the equatorial anomaly band, the annual component is weaker than its Fourier decomposition periods. In particular, the semiannual period dominates the annual component, indicating the relationship between the semiannual variation of the electron densities and the ionospheric equatorial anomaly. From detection of the phases of the components, it is revealed that the annual signal generally has its maximum value in summers at high altitudes, and in the winters at low altitudes. This is probably due to the higher [O/N₂] ratios in winter than in the summer in the lower ionosphere. Furthermore, the semiannual component mostly peaks around solstices or about a month before/after them.     

Validation of FORMOSAT-3/COSMIC radio occultation electron density profiles by incoherent scatter radar data

Low Earth Orbiting satellites carrying a dual frequency GPS receiver onboard offer a unique opportunity to remote sensing of the global ionosphere on a continuous basis. No other profiling technique unifies profiling through the entire F2-layer with global coverage. The FORMOSAT-3/COSMIC data can make a positive impact on the global ionosphere study providing essential information about the height electron density distribution and particularly over regions that are not accessible with ground-based measuring instruments such as ionosondes and GPS dual frequency receivers. Therefore, it is important to verify occultation profiles with other techniques and to obtain experience in the reliability of their derivation. In the given study we present results of comparison of the electron density profiles derived from radio occultation measurements on-board FS-3/COSMIC and from the Kharkov incoherent scatter radar sounding.     

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.     

Theoretical simulation of electron density and temperature distribution at Indian equatorial and low latitude ionosphere

The electron density and temperature distribution of the equatorial and low latitude ionosphere in the Indian sector has been investigated by simultaneously solving the continuity, momentum and energy balance equations of ion and electron flux along geomagnetic field lines from the Northern to the Southern hemisphere. Model algorithm is presented and results are compared with the electron density and electron temperature measured in situ by Indian SROSS C2 satellite at an altitude of ∼500 km within 31°S–34°N and 75 ± 10°E that covers the Indian sector during a period of low solar activity. Equatorial Ionization Anomaly (EIA) observed in electron density, morning and afternoon enhancements, equatorial trough in electron temperature have been simulated by the model within reasonable limits of accuracy besides reproducing other normal diurnal features of density and temperature.     

Radio occultation electron density profiles from the FORMOSAT-3/COSMIC satellites over the Brazilian region: A comparison with Digisonde data

This study aims to validate the electron density profiles from the FORMOSAT-3/COSMIC satellites with data from Digisondes in Brazil during the low solar activity period of the years 2006, 2007 and 2008. Data from three Brazilian Digisondes located in Cachoeira Paulista (22.7°S, 45°W), São Luís (2.5°S, 44.2°W) and Fortaleza (3.8°S, 38°W) were used in the comparisons. Only the profiles whose density peak have been obtained near the stations coordinates were chosen for the comparison. Although there is generally good agreement, some cases of discrepancies are observed. Some of these discrepancies cannot be explained simply by the differences in the position and local time of the measurements made by the satellite and the ground-based station. In such cases it is possible that local conditions, such as the presence of a trans-equatorial wind or electron density gradients, could contribute to the observed differences. Comparison of the F2 layer peak parameters, the NmF2 and hmF2, obtained from the two techniques showed that, in general, the agreement for NmF2 is pretty good and the NmF2 has a better correlation than hmF2. Cachoeira Paulista had the worst correlation for hmF2 possibly because this station is situated in the region under the influence of the equatorial ionization anomaly, a region where it is more difficult to apply the RO technique without violating the spherical symmetry condition.     

Plasmaspheric electron content derived from GPS TEC and FORMOSAT-3/COSMIC measurements: Solar minimum condition

The plasmaspheric electron content (PEC) was estimated by comparison of GPS TEC observations and FORMOSAT-3/COSMIC radio occultation measurements at the extended solar minimum of cycle 23/24. Results are retrieved for different seasons (equinoxes and solstices) of the year 2009. COSMIC-derived electron density profiles were integrated up to the height of 700 km in order to retrieve estimates of ionospheric electron content (IEC). Global maps of monthly median values of COSMIC IEC were constructed by use of spherical harmonics expansion. The comparison between two independent measurements was performed by analysis of the global distribution of electron content estimates, as well as by selection specific points corresponded to mid-latitudes of Northern America, Europe, Asia and the Southern Hemisphere. The analysis found that both kinds of observations show rather similar diurnal behavior during all seasons, certainly with GPS TEC estimates larger than corresponded COSMIC IEC values. It was shown that during daytime both GPS TEC and COSMIC IEC values were generally lower at winter than in summer solstice practically over all specific points. The estimates of PEC (h > 700 km) were obtained as a difference between GPS TEC and COSMIC IEC values. Results of comparative study revealed that for mid-latitudinal points PEC estimates varied weakly with the time of a day and reached the value of several TECU for the condition of solar minimum. Percentage contribution of PEC to GPS TEC indicated the clear dependence from the time with maximal values (more than 50–60%) during night-time and lesser values (25–45%) during day-time.     

An update global model of hmF2 from values estimated from ionosonde and COSMIC/FORMOSAT-3 radio occultation

In this paper, we present our recent work on developing an updated global model of the ionospheric F2 peak height hmF2 parameter by combining data from the Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC/FORMOSAT-3) radio occultation (RO) measurements and from the extended global ionosonde stations. In particular, 10 Chinese ionosonde stations’ data are newly introduced into this study. The modeling technique used is based on a two-layer empirical orthogonal function (EOF) expansion. Global distributions of hmF2 maps calculated using the newly constructed global model and the one provided by the International Reference Ionosphere model (IRI-ITU-R) are compared with the global distributions of hmF2 obtained by the COSMIC RO measurements and quantitative statistical analysis of the differences between the model results and those of the COSMIC RO measurements is made for the low (2008) and high (2012) solar activity years. The obtained average root-mean-square differences (RMSEs) for our model are 27.7 km (11.1%) and 31.0 km (9.8%), respectively for the years 2008 and 2012, whereas those for the IRI-ITU-R model are 39.9 km (16.9%) and 35.0 km (11.6%), respectively. Comparison of the results calculated both by our model and the IRI-ITU-R model with the digisonde observation is also made. The comparisons show that the newly constructed global hmF2 model can reproduce reasonably well the observations and perform better than IRI-ITU-R model.     

Assessment of precision in ionospheric electron density profiles retrieved by GPS radio occultations

The Constellation Observing System for Meteorology Ionosphere and Climate (COSMIC) is a six satellite radio occultation mission that was launched in April 2006. The close proximity of these satellites during some months after launch provides a unique opportunity to evaluate the precision of Global Positioning System (GPS) radio occultation (RO) retrievals of ionospheric electron density from nearly collocated and simultaneous observations. RO data from 30 consecutive days during July and August 2006 are divided into ten groups in terms of daytime or nighttime and latitude. In all cases, the best precision values (about 1%) are found at the F peak height and they slightly degrade upwards. For all daytime groups, it is seen that electron density profiles above about 120 km height exhibit a substantial improvement in precision. Nighttime groups are rather diverse: in particular, the precision becomes better than 10% above different levels between 120 and 200 km height. Our overall results show that up to 100–200 km (depending on each group), the uncertainty associated with the precision is in the order of the measured electron density values. Even worse, the retrieved values tend sometimes to be negative. Although we cannot rely directly on electron density values at these altitudes, the shape of the profiles could be indicative of some ionospheric features (e.g. waves and sporadic E layers). Above 200 km, the profiles of precision are qualitatively quite independent from daytime or latitude. From all the nearly collocated pairs studied, only 49 exhibited a difference between line of sight angles of both RO at the F peak height larger than 10°. After analyzing them we find no clear indications of a significant representativeness error in electron density profiles due to the spherical assumption above 120 km height. Differences in precision between setting and rising GPS RO may be attributed to the modification of the processing algorithms applied to rising cases during the initial period of the COSMIC mission.     

Equatorial vertical E × B drift velocities inferred from ionosonde measurements over Ouagadougou and the IRI-2007 vertical ion drift model

F-region vertical plasma drift velocities were deduced from the hourly hmF2 values acquired from ionogram data over a near dip equatorial station Ouagadougou (12.4°N, 358.5°E, dip angle 5.9°N) in Africa. Our results are compared against the global empirical model of Scherliess and Fejer (1999) incorporated in the IRI model (IRI-2007) for 1600 to 0800 LT from 1 year of data during sunspot maximum year of 1989 (yearly average solar flux intensity, F10.7 = 192) corresponding to the peak phase of solar cycle 22, under magnetically quiet conditions. The drifts are entirely downward between 2000 and 0500 LT bin for both techniques and the root mean square error (RMSE) between the modeled and the ionosonde vertical plasma drifts during these periods is 3.80, 4.37, and 4.74 m/s for June solstice, December solstice and equinox, respectively. Ouagadougou average vertical drifts show evening prereversal enhancement (PRE) velocity peaks (V_ZP) of about 16, 14, and 17 m/s in June solstice, December solstice, and equinox, respectively, at 1900–2000 LT; whereas global empirical model average drifts indicate V_ZP of approximately 33 m/s (June solstice), 29 m/s (December solstice), and 50 m/s (equinox) at 1800 LT. We find very weak and positive correlation (+0.10376) between modeled V_ZP versus F10.7, while ionosonde V_ZP against F10.7 gives worst and opposite correlation (−0.05799). The results also show that modeled V_ZP–Ap indicates good and positive correlation (+0.64289), but ionosonde V_ZP–Ap exhibits poor and negative correlation (−0.22477).     

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.     

Variations in statistical parameters of the NmF2 equinoctial asymmetry with latitude and solar activity near noon

Hourly values of the F2-layer peak density, NmF2, measured by 95 ionosondes near noon from 1957 to 2011 at low and middle geomagnetic latitudes of the northern and southern geographic hemispheres are used in a statistical study of the NmF2 equinoctial asymmetry. The ratios, R, of NmF2 measured during 61 days around the March equinox to NmF2 measured during 61 days around the September equinox at the same UT near noon during geomagnetically quiet daytime conditions for approximately the same solar activity conditions over the same ionosonde are analyzed. The conditional probability of the occurrence of R in an interval of R, the most probable value of R, and the mean expected value of R are calculated for the first time for the low, moderate, and high solar activity levels to study variations in these statistical parameters with latitude and solar activity. These statistical parameters are averaged over 5° geomagnetic latitude interval in the northern and southern geographic hemispheres to calculate and to study for the first time trends in latitude and solar activity of these averaged NmF2 equinoctial asymmetry statistical characteristics.     

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.     

A statistical study of the mid-latitude NmF2 winter anomaly

Hourly values of the F2-layer peak density, NmF2, measured by 62 ionosonde stations from 1957 to 2008 at middle geomagnetic latitudes of the northern and southern hemispheres are used in a statistical study of the F2-region winter anomaly. We analyze a maximum value, R, of the winter/summer geomagnetically quiet daytime NmF2 ratio over each ionosonde for approximately the same winter and summer solar activity conditions. The conditional probability of an occurrence of R in an interval of R, the most frequent value of R, the mean expected value of R, and the conditional probability to observe the F2-region winter anomaly are calculated and studied for low, moderate, and high solar activity conditions. It is found that northern hemisphere stations show significant cross-correlations of winter anomaly statistics with geomagnetic latitude for moderate and high solar activities.     

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.     

Evaluation of the IRI-2016 and NeQuick electron content specification by COSMIC GPS radio occultation,ground-based GPS and Jason-2 joint altimeter/GPS observations

We examined performance of two empirical profile-based ionospheric models, namely IRI-2016 and NeQuick-2, in electron content (EC) and total electron content (TEC) representation for different seasons and levels of solar activity. We derived and analyzed EC estimates in several representative altitudinal intervals for the ionosphere and the plasmasphere from the COSMIC GPS radio occultation, ground-based GPS and Jason-2 joint altimeter/GPS observations. It allows us to estimate a quantitative impact of the ionospheric electron density profiles formulation in several altitudinal intervals and to examine the source of the model-data discrepancies of the EC specification from the bottom-side ionosphere towards the GPS orbit altitudes. The most pronounced model-data differences were found at the low latitude region as related to the equatorial ionization anomaly appearance. Both the IRI-2016 and NeQuick-2 models tend to overestimate the daytime ionospheric EC and TEC at low latitudes during all seasons of low solar activity. On the contrary, during high solar activity the model results underestimated the EC/TEC observations at low latitudes. We found that both models underestimated the EC for the topside ionosphere and plasmasphere regions for all levels of solar activity. For low solar activity, the underestimated EC from the topside ionosphere and plasmasphere can compensate the overestimation of the ionospheric EC and, consequently, can slightly decrease the resulted model overestimation of the ground-based TEC. For high solar activity, the underestimated EC from the topside ionosphere and plasmasphere leads to a strengthening of the model underestimation of the ground-based TEC values. We demonstrated that the major source of the model-data discrepancies in the EC/TEC domain comes from the topside ionosphere/plasmasphere system.     

Influence of the equinoctial asymmetry on the semi-annual anomaly in NmF2 near noon in the northern geographic hemisphere: A statistical study

Hourly values of NmF2 measured by 72 ionosondes near noon from 1957 to 2012 at low and middle geomagnetic latitudes of the northern geographic hemisphere are used to study the spring and autumn semi-annual anomalies in NmF2. The spring/summer, autumn/summer, spring/winter, and autumn/winter geomagnetically quiet NmF2 ratios, X_M, X_S, Y_M, and Y_S, respectively, taken near noon over each ionosonde for approximately the same winter, spring, summer, and autumn solar activity conditions are analyzed. The probabilities of occurrences of X_M, X_S, Y_M, and Y_S in intervals of change of these parameters, together with its most frequent and the mean expected values, and the probabilities to observe the F2-region spring and autumn semi-annual anomalies are calculated and studied for the first time for low, moderate, and high solar activity conditions. The geomagnetic latitude range between 10 and 55 degrees, where the ionosondes are located, is divided into 9 intervals of the same length of 5 degrees, the statistical parameters are averaged over each 5 degree interval, and the trends in these averaged statistical characteristics of the NmF2 spring and autumn semi-annual anomalies are calculated and studied for the first time. It is shown that the NmF2 equinoctial asymmetry can significantly affect the statistical parameters of the semi-annual anomaly in NmF2.     

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.     

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.     

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.     

The influence of the spatial and temporal collocation windows on the comparisons of the ionospheric characteristic parameters derived from COSMIC radio occultation and digisondes

GPS radio occultation (RO) ionospheric products obtained by Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC) mission during the year of 2014 and the observations from 3 digisonde stations which are located at different latitudes are used to study the influence of different time and space collocation windows on the comparisons of the ionospheric characteristic parameters (ICPs), including the peak density and peak height, derived from the two techniques. The results show that the correlation coefficients (CC) and the standard deviation of the absolute biases (SDAB) between the ICPs derived from the two techniques vary distinctly under different spatial and time collocation windows. Generally, the CC (SDAB) of the ICPs decrease (increase) as the size of the collocation window increases in time dimension or in space dimension. The rate of change of the statistic parameters with the increase in the size of the collocation window in time dimension and space dimension is analyzed for each digisonde station. It is found that within the collocation window of $\left(60\mathrm{m}\mathrm{i}\mathrm{n},{20}^{°},{20}^{°}\right)$, the influence of the increase of $1^{°}$ in the space window on the statistical comparison is much more significant than that of the increase of 1?min in the time window, and it is suggested that there can be appropriate relaxation on the time window within the threshold of 60?min to get a balance between the quality of the comparison results and the number of the matched pairs. In addition, it is found that the same variations in the longitude window and in the latitude window may have different influences on the comparison results when the horizontal gradients in electron density are distinctly different along different directions at the digisonde station, and strict space collocation window is preferred when comparing the observations from COSMIC RO with those from the digisonde station in such cases.