Contributions of incoherent-scatter radar measurements at Arecibo for the improvement of IRI

In the years 1974 to 1977 the incoherent-scatter radar at Arecibo made high resolution measurements of electron density in time and altitude, under the I29 program. This program ran on monthly basis and each observing run was for about 36 hours, consisting of two daytime periods and the intervening night. Electron density was measured from 100 to 500 km with a single 24 μs pulse which provided an altitude resolution of 3.6 km. These measurements provided very accurate values of hmF2, h0.5, E-F valley parameters and topside Ne gradient thereby providing important inputs for the improvement of IRI. This paper gives a summary of these inputs.     

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.     

Comparison of measured and predicted F2 peak altitude

Different models for the F2 peak altitude hmF2 are compared with mean values determined from incoherent scatter measurements of the radar stations Jicamarca/Peru and Arecibo/Puerto Rico at all local times. Our investigation shows that the daytime hmF2 is fairly well represented by the most recent models, whereas at nighttime, higher peak altitudes are measured than predicted by either model. The observed diurnal structure is slightly misinterpreted by the models above Jicamarca and is strongly disagreeing above Arecibo.     

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.     

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.     

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.     

Reconstruction of topside density profile by using the topside sounder model profiler and digisonde data

To improve the accuracy of the real time topside electron density profiles given by the Digisonde software a new model-assisted technique is used. This technique uses the Topside Sounder Model (TSM), which provides the plasma scale height (H_s), O⁺–H⁺ transition height (H_T), and their ratio Rt = H_s/H_T, derived from topside sounder data of Alouette and ISIS satellites. The Topside Sounder Model Profiler (TSMP) incorporates TSM and uses the model quantities as anchor points in construction of topside density (Ne) profiles. For any particular location, TSMP calculates topside Ne profiles by specifying the values of foF2 and hmF2. In the present version, TSMP takes the F2 peak characteristics – foF2, hmF2, and the scale height at hmF2 – from the Digisonde measurements. The paper shows results for the Digisonde stations Athens and Juliusruh. It is found that the topside scale height used in Digisonde reconstruction is less than that extracted from topside sounder profiles. Rough comparison of their bulk distributions showed that they differ by an average factor of 1.25 for locations of Athens and Juliusruh. When the Digisonde scale heights are adjusted by this factor, the reconstructed topside profiles are close to those provided by TSM. Compared with CHAMP reconstruction profiles in two cases, TSMP/Digisonde profiles show lower density between 400 and 2000 km.     

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.     

A global survey of COSMIC ionospheric peak electron density and its height: A comparison with ground-based ionosonde measurements

With a network of ground-based ionosondes distributed around the world, the ionospheric peak electron density and its height measured by FORMOSAT-3/COSMIC satellites in terms of GPS radio occultation technique are extensively examined in this article. It is found that, in spite of the latitude, the mean values of the peak electron density measured by COSMIC satellites are systematically smaller than those observed by ground-based ionosondes. The discrepancy between them is dependent on the latitude, namely, it is small in low and mid-latitudes and large in high-latitude region. Moreover, statistical analysis shows that the slopes of the regression line that is best fitted to the scatter diagram of occultation-retrieved peak electron density (ordinate axis) versus ionosonde-observed peak density (abscissa axis) are universally less than one. This feature is believed to be the result of path average effect of non-uniform distribution of the electron density along the GSP ray during the occultation. A comparison between COSMIC-measured peak height and ionosonde-derived peak height hmF2 indicates that the former is systematically higher than the latter. The difference in the two can be as large as 20% or more in equatorial and low-latitude regions. This result implies that the peak height hmF2 derived from the virtual height through true height analysis based on Titheridge method seems to underestimate the true peak height. The correlation between COSMIC and ionosonde peak electron densities is analyzed and the result reveals that correlation coefficient seems to be dependent on the fluctuation of the occultation-retrieved electron density profile. The correlation will be higher (lower) for the electron density profiles with smaller (larger) fluctuations. This feature suggests that the inhomogeneous distribution of the electron density along the GPS ray path during the occultation plays an important role affecting the correlation between COSMIC and ionosonde measurements.     

磁暴期间中纬度电离层剖面结构变化的数值模拟

利用电离层理论模型模拟了磁暴期间热层大气温度、成分、中性风和电场扰动对电离层电子密度剖面结构，特别是峰值密度和峰值高度变化的影响，结果表明，热层大气温度变化所引起光化反应系数的改变对电离层剖面结构影响不大；热层大气成分特别是Ｎ２／Ｏ的变化能有效地引起密度剖面变化，Ｎ２增加足以使峰值密度产生所观测到的负相暴；由中性风和电场引起等离子体漂移是峰值高度ｈｍＦ２变化的主要原因，但对电子密度的影响不足以抵消     

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

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.     

Initial results of the evaluation of IRI hmF2 performance for minima 22–23 and 23–24

The performance of the International Reference Ionosphere (IRI) in predicting the height of the maximum of electron density (hmF2) has been evaluated for similar geomagnetic latitudes stations in the northern hemisphere (NH) and southern hemisphere (SH), and for the last two minima. As truth-sites, the digisonde stations of Millstone Hill (42.6°N, 288.5°E), USA, and Grahamstown (33.3°S, 26.5°E), South Africa, were considered. A monthly averaged diurnal variation was obtained from all the observations and model output in the months studied, and the corresponding difference was also calculated. For this initial study data from summer and winter in the NH and SH were selected for the solstice comparison, and October data for both stations were used to represent equinox conditions. The choice of these periods depended on data availability and quality. The results show that for the earlier minimum in 1996, in general IRI hmF2 values are in reasonable agreement with the observations. The exceptions are October and December in the SH, where IRI hmF2 tends to high, particularly on the dayside, and also July for which the daytime measured values tend to be larger than the IRI ones. For the recent minimum in 2008, IRI tends to over-estimate the hmF2 in most of the observations. The results support the general assertion that thermospheric temperatures were cooler during the last solar minimum as a consequence of an unusually low, and extended, minimum in solar extreme-ultraviolet flux, and in response to continually increasing long-term trend in anthropogenic carbon dioxide. The cooler temperatures not only decrease density at a fixed height, but also make the corresponding contraction of the atmosphere lower the height of the F-region peak.     

Incoherent-scatter measurements of E-F valley and comparisons with theoretical and empirical models

High resolution electron density measurements from Arecibo incoherent scatter (i.s.) radar are used to make a detailed study of the E-F valley. Features of the important valley parameters like height, width and depth are examined. These features are then compared with the available theoretical and empirical models. The depth of the valley obtained from the empirical models agrees with i.s. measurements for near-noon periods, but disagrees with these measurements for pre-noon and post-noon periods. Further, the i.s. measurements indicate that E-F valley is rather small during daytime as compared to models which give larger width. During the night, the valley is quite wide and deep but the presence of sporadic-E (Es) contaminates the Ne-h profiles observed with the i.s. radar. As a result the valley parameters cannot be determined unambiguously during the night.     

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.     

Comparison of peak characteristics of F2 ionospheric layer over Tehran region at a low solar activity period with IRI-2001 and IRI-2007 models predictions

In the present work values of peak electron density (NmF2) and height of F2 ionospheric layer (hmF2) over Tehran region at a low solar activity period are compared with the predictions of the International Reference Ionosphere models (IRI-2001 and IRI-2007). Data measured by a digital ionosonde at the ionospheric station of the Institute of Geophysics, University of Tehran from July 2006 to June 2007 are used to perform the calculations. Formulations proposed by  and  are utilized to calculate the hmF2. The International Union of Radio Science (URSI) and International Radio Consultative Committee (CCIR) options are employed to run the IRI-2001 and IRI-2007 models. Results show that both IRI-2007 and IRI-2001 can successfully predict the NmF2 and hmF2 over Tehran region. In addition, the study shows that predictions of IRI-2007 model with CCIR coefficient has closer values to the observations. Furthermore, it is found that the monthly average of the percentage deviation between the IRI models predictions and the values of hmF2 and NmF2 parameters are less than 10% and 21%, respectively.     

非相干雷达探测低电离层的结果分析

位于波多黎各的Arecibo非相干雷达可以获得低电离层电子和离子密度, 利用此非相干雷达数据对中纬度低电离层的运动特征进行研究. 得到了电子密度随时间和高度的变化 情况, 结果显著呈现出周日变化特征, 并分析了电子密度随高度的变化规律. 进一步对数据进行频谱分析, 深入研究低电离层电子密度的周日变化效应. 得到电子密度的高度剖面, 发现从F层底部到E层有明显的等离子体沉降. 低电离层的层结构特征及电子密度变化表明, 在该区域还存在不同程度的等离子体扰动, 由此对低电离层的作用因素 进行分析, 认为大气潮汐或声重波可能对低电离层产生扰动, 即低电离层与大气存在一定程度的耦合作用.      

Hemispheric,seasonal and latitudinal dependence of storm-time ionosphere during low solar activity period

Analysis of the seasonal, hemispheric and latitudinal variation of the ionospheric F2 peak during periods of disturbed geomagnetic conditions in 2011, a year of low solar activity, had been studied using hourly data obtained from low- and mid-latitude ionosonde stations. Our results showed an enhancement in F2-layer maximum electron density (NmF2) at daytime over low latitudes. For the mid-latitude stations, NmF2 depletion pre-dominates the daytime and overturned at nighttime. In general, the variation in terms of magnitude is higher in the low-latitude than at mid-latitude. The nighttime decrease in NmF2 is accompanied by a corresponding F2 peak height (hmF2) increase and overturned at daytime. The hmF2 response during the equinoctial months is lower than the solstices. NmF2 shows distinct seasonal, hemispheric and latitudinal dependence in its response. Appearance of a significant ionospheric effect in southern hemisphere is higher than in the northern hemisphere, and is more pronounced in the equinoxes at low latitudes. At mid-latitudes, the ionospheric effect is insignificant at both hemispheres. A negative ionospheric response dominates the whole seasons at the mid-latitude except for March equinox. The reverse is the case for the hmF2 observation. The amplitudes of both the NmF2 and hmF2 increase with increasing latitude and maximize in the southern hemisphere in terms of longitude.     

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.     

Experimental determination of the main F-layer parameters on a global scale using data from a remote optical probe onboard the satellite Interkosmos-Bulgaria-1300

The optical complex of the satellite Interkosmos-Bulgaria-1300 includes the UV-spectrometer “Photon-1”, and the scanning spectrometer EMO-5 operating in the visible spectrum. Using the intensities of the OI 1356 and 6300 nm emission as measured by these two spectrometers, a theoretical evaluation of the two main F-layer parameters has been made: hmF, the height of the F2-maximum and NmF, the electron density at this height, together with the constant Ao characterizing the Ne(h) distribution in the IRI. The MSIS and CIRA-72 models are used for determining the distribution of the neutral component. If we had used CIRA-72, lower NmF values and higher hmF values had been obtained than with MSIS. The results obtained show that the method is very sensitive to the assumed neutral model, particularly when determining hmF. This work demonstrates the possibility of probing the F-region near its peak on a global scale using comperatively cheap and easily available optical methods.