Verification of IRI plasma temperature at great altitude by satellite data

Electron temperature data obtained aboard two satellites are compared with the IRI model. Our comparison may serve as an input to a future version of IRI. HINOTORI data show IRI90 to predict night time T_e quite well with the exception of a few peculiar phenomena, in particular a morning overshoot at height around 600 km and anomalous T_e observations of AKEBONO for the first time allowed to extend reliably the temperature profile up to − 8000km. This profile clearly depends on local time at all invariant latitudes. The data suggest strong influence of solar activity on electron temperature.     

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.     

Ionospheric electron density profiles at sunrise-sunset

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

An empirical model of electron temperature for low and middle latitudes in the 100–200 km height region

An empirical model of electron temperature (T_e) for low and middle latitudes is proposed in view of IRI. It is constructed on the basis of experimental data obtained at 100 to 200 km by probe and incoherent scatter methods. Below 150 km the model gives two T_e values: one from incoherent scatter data and another from probe measurements. The model can be used for all seasons for quiet geomagnetic conditions (Kp not greater 3) and at almost all levels of solar activity (F_10.7 between 70 and 200). It is presented in an analytical form that allows one to calculate T_e profiles for different latitudes, longitudes and at any season (day). Depending on geomagnetic latitude and solar zenith angle, electron temperature distributions are presented for two heights along with T_e profile variations during the day (at middle latitudes).     

Analysis of global TEC prediction performance of IRI-PLAS model

The International Reference Ionosphere (IRI) empirical model provides valuable data for many fields including space and navigation applications. Since the IRI model gives the ionospheric parameters in the altitude range from 50?km to 2000?km, researchers focused on the IRI-PLAS model which is the plasmasphere extension of the IRI model. In this study, Total Electron Content (TEC) prediction performance of the IRI-PLAS model was examined at a global scale using the location of globally distributed 9 IGS stations. Besides the long term (01.01.2015–31.12.2015) behavior of the model, TEC predictions during the equinox and solstice days of 2014–2017 were also tested. IRI-PLAS-TEC values were examined in comparison with GPS-TEC data. Hourly interval of yearly profile exhibits that when the geomagnetic and solar active days are ignored, differences between IRI-PLAS-TEC and GPS-TEC are rather small (～2–3 TECU) at stations in the northern hemisphere, generally ～4–5 TECU level at the southern hemisphere stations and reaching above 10 TECU for few hours. While the IRI-PLAS-TEC generally overestimates the GPS-TEC at southern hemisphere stations during quiet days, the model-derived TEC underestimates GPS-TEC during solar active days. IRI-PLAS-TEC and GPS-TEC values exhibit similar trend for the equinoxes 21 March and 23 September which refer equivalent conditions.     

Electron density reference profile in the lower ionosphere

The new IRI formula, as accepted at the 1983 Stara Zagora Workshop, prescribes the use of Epstein functions for reproducing logarithmic electron density profiles. In this paper we discuss solutions which might be applicable to the lower ionosphere. The experimental data base is briefly reviewed. It appears that the stratification near 80 km must be accepted as a regular feature of the daytime lower ionosphere. The C-layer problem is left open. In order to reproduce such profiles, one needs three LAY-functions. Examples show that the weighted sum of these does very well represent experimental profiles, the amplitudes being determined by a least square fit. For profile synthesis (as in IRI) a least square determination of the three amplitudes, admitting four linear conditions, is proposed.     

Modelling of ionospheric temperature profiles

The amount of measured temperature data accumulated in recent years allows and asks for improvement and refinement of the rather crude temperature models employed in the International Reference Ionosphere 1979. By combining the mission-oriented models by BRACE, THEIS /2/ for the AE-C satellites and by SPENNER, PLUGGE /1/ for the AEROS-A satellite, a much better diurnal and latitudinal reliability can be obtained. It is also suggested that IRI should have the option to make use of the strong anti-correlation between electron temperature and density in cases where actual measured densities are available. For daytime condition, incorporation of a density dependent model into IRI can significantly enhance the prediction quality of IRI in the altitude range 300 to 600 km.Furtheron, the solar activity dependence of the electron temperature and of the density-temperature-relation are investigated by comparing the low solar activity models with more recent AE-C and -DE data.     

Electron density models for the lower ionosphere

Probably the only reliable method of checking an electron density model below 70 km is to calculate from it what would be obtained by VLF or LF propagation over certain paths, and to compare the results with actual observations. This has been done for the IRI at various frequencies from 16 to 70 kHz; the results agree in places but differ substantially elsewhere. Previous models described by the author give satisfactory results and it is suggested that certain features of them might be incorporated with advantage in the IRI. In particular, it is impossible to get agreement with VLF propagation in all seasons by means of a model varying only with solar zenith angle, such as the IRI from 50–90 km.     

Prediction of total electron content using the International Reference Ionosphere

The International Reference Ionosphere (IRI) is a model of the ionosphere, based on experimental data, which has been proposed as a standard ionospheric model. As such, it should be tested extensively to determine its range of validity. One of the ways in which the electron denisty profile given by the IRI, especially above the peak of the F layer, can be tested is to compare calculated and observed values of total electron content (TEC). We have therefore studied the discrepancies between calculated and observed values of TEC recorded at 15 stations covering a wide range of longitudes and latitudes, mainly in the northern hemisphere, and mainly for high levels of solar activity. W have found that the IRI produces reasonably accurate values of TEC at mid and high latitudes, but that it greatly underestimates the daytime values of TEC at low latitudes. We conclude therefore that the daytime electron density profile given by the IRI is reasonably accurate at mid and high latitudes, at least above the peak of the F2 layer. The situation at low latitudes clearly requires more work, and we have suggested two possible lines of study. The generally low discrepancies at night indicate that the night-time electron density profiles given by the IRI correspond fairly closely to the actual profiles.     

Global validation of IRI TEC for high and medium solar activity conditions

The IRI model offers a choice of options for the computation of the electron density profile and electron content (TEC). Recently new options for the topside electron density profile have been developed, which have a strong impact on TEC. Therefore it is important to test massively the IRI and the new options with experimental data. A large number of permanent stations record dual frequency GPS data from which it is possible to obtain TEC values. Thirty-one worldwide distributed stations have been selected to investigate the capabilities of the IRI to reproduce experimental TEC. Data for years 2000 (high solar activity) and 2004 (medium solar activity) have been analyzed computing modeled values with the IRI-2001 and the IRI-2007-NeQuick topside options. It is found that IRI-2007-NeQuick option generally improves the estimate of the slant TEC, especially in the case of high latitudes stations during high solar activity.     

IRI 2001/90 TEC predictions over a low latitude station

Total electron content (TEC) over Tucumán (26.9°S, 294.6°W) measured with Faraday technique during the high solar activity year 1982, is used to check IRI 2001 TEC predictions at the southern crest of the equatorial anomaly region. Comparisons with IRI 90 are also made. The results show that in general IRI overestimates TEC values around the daily minimum and underestimates it the remaining hours. Better predictions are obtained using ground ionosonde measurements as input coefficients in the IRI model. The results suggest that for hours of maximum TEC values the electron density profile is broader than that assumed by the model. The main reason for the disagreement would be the IRI shape of the electron density profile.     

On the use of NeQuick topside option in IRI-2007

The latest version of IRI includes various options for the computation of the topside electron density profile. One of the possible choices is based on NeQuick model. Its inclusion in IRI has been made transferring all the formulations used in NeQuick model. In details, an Epstein layer function is used to describe the electron density profile and the topside shape is controlled by an empirical parameter, connected to the NeQuick F2 bottomside thickness parameter, B2bot. It is computed also in this IRI topside option in order to maintain self-consistency with its original formulation. This paper analyses the possibility of using the IRI bottomside parameters for this option and its impact on the profile and TEC. The case of experimental peak values given as input is also analysed.     

A 12 year comparison of MIDAS and IRI 2007 ionospheric Total Electron Content

Data assimilation in conventional meteorological applications uses measurements in conjunction with a physical model. In the case of the ionised region of the upper atmosphere, the ionosphere, assimilation techniques are much less mature. The empirical model known as the International Reference Ionosphere (IRI) could be used to augment data-sparse regions in an ionospheric now-cast and forecast system. In doing so, it is important that it does not introduce systematic biases to the result. Here, the IRI model is compared to ionospheric observations from the Global Positioning System satellites over Europe and North America. Global Positioning System data are processed into hour-to-hour monthly averages of vertical Total Electron Content using a tomographic technique. A period of twelve years, from January 1998 to December 2009, is analysed in order to capture variations over the whole solar cycle. The study shows that the IRI model underestimates Total Electron Content in the daytime at solar maximum by up to 37% compared to the monthly average of GPS tomographic images, with the greatest differences occurring at the equinox. IRI shows good agreement at other times. Errors in TEC are likely due to peak height and density inaccuracies. IRI is therefore a suitable model for specification of monthly averages of Total Electron Content and can be used to initialise a data assimilation process at times away from solar maximum. It may be necessary to correct for systematic deviations from IRI at solar maximum, and to incorporate error estimation into a data assimilation scheme.     

A simplified D-region ion chemistry scheme and its possible use for IRI lower ionosphere modelling

Ion composition of the D region is principally characterized by the existence of two distinct regions of predominant molecular ions and predominant cluster ions, separated from each other by a rather sharp ‘transition height’, which is proposed to be included in the IRI as an additional parameter, supplementing the electron density models. It is possible to predict the position of this ‘transition height’ at a given place and time with the aid of a simplified ion chemistry scheme which is shown to be satisfactorily compatible with experimental ion composition data available in the literature. Our suggested method of this prediction makes use of the (IRI or experimental) electron density profile at the location and season in question, together with an effective clustering rate coeeficient calculated from corresponding temperature and density profiles taken from a suitable reference model of the neutral atmosphere.     

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

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

瑞利激光雷达反演中层大气温度算法

瑞利激光雷达是探测中层大气温度的重要设备,具有高时空分辨率且不存在探测盲区的优点.为充分利用激光雷达探测到的原始数据,改进了传统的Chanin-Haunchecorne方法,采用由均匀搜索生成温度初值的方式反演大气温度.位于北京延庆（40.3°N,116.2°E）的瑞利激光雷达具有589nm和532nm两个通道,将589nm通道用于计算,532nm通道作为参考.在150~250K区间等间隔选取多个温度初值,利用589nm通道反演60~70km高度范围内的大气温度廓线和大气密度廓线,利用参考密度廓线选取准确的温度初值,反演得到准确的大气温度廓线.将589nm通道和532nm通道反演的温度廓线相比较,发现二者具有较高的一致性.改进后的方法有效利用了信噪比较差的光信号,使589nm通道的温度探测上限从60km提高到70km.