Global comparison of the model results of GSM TIP with IRI for summer conditions

This paper presents the results of the numerical calculations thermosphere/ionosphere parameters which were executed with using of the Global Self-consistent Model of the Thermosphere, Ionosphere and Protonosphere (GSM TIP)and comparison of these results with empirically-based model IRI-2001. Model GSM TIP was developed in West Department of IZMIRAN and solves self-consistently the time-dependent, 3-D coupled equations of the momentum, energy and continuity for neutral particles (O₂, N₂, O), ions (O⁺, H⁺), molecular ions (M⁺) and electrons and largescale eletric field of the dynamo and magnetospheric origin in the range of height from 80 km to 15 Earth’s radii. The empirically derived IRI model describes the E and F regions of the ionosphere in terms of location, time, solar activity and season. Its output provides a global specification not only of N_e but also on the ion and electron temperatures and the ion composition. These two models represent a unique set of capabilities that reflect major differences in along with a substantial approaches of the first-principles model and global database model for the mapping ionosphere parameters. We focus on global distribution of the N_e, T_i, T_e and TEC for the one moment UT and fixed altitudes: 110 km, h_mF2, 300 km and 1000 km. The calculations were executed with using of GSM TIP and IRI models for August 1999, moderate solar activity and quiet geomagnetic conditions. Results present as the global differences between the IRI and GSM TIP models predictions. The discrepancies between model results are discussed.     

Latitudinal variation of the topside electron temperature at different levels of solar activity

A database of electron temperature (T_e) measurements comprising of most of the available satellite measurements in the topside ionosphere is used for studying the solar activity variations of the electron temperature T_e at different latitudes, altitudes, local times and seasons. The T_e data are grouped into three levels of solar activity (low, medium, high) at four altitude ranges, for day and night, and for equinox and solstices. We find that in general T_e changes with solar activity are small and comparable in magnitude with seasonal changes but much smaller than the changes with altitude, latitude, and from day to night. In all cases, except at low altitude during daytime, T_e increases with increasing solar activity. But this increase is not linear as assumed in most empirical T_e models but requires at least a parabolic approximation. At 550 km during daytime negative as well as positive correlation is found with solar activity. Our global data base allows to quantify the latitude range and seasonal conditions for which these correlations occur. A negative correlation with solar activity is found in the invdip latitude range from 20 to 55 degrees during equinox and from 20 degrees onward during winter. In the low latitude (20 to −20 degrees invdip) F-region there is almost no change with solar activity during solstice and a positive correlation during equinox. A positive correlation is also observed during summer from 30 degrees onward.     

Modelling bottom and topside electron density and TEC with profile data from topside ionograms

The paper describes the technique that has been implemented to model the electron density distribution above and below the F2 peak making use of only the profiles obtained from the INTERCOSMOS-19 topside ionograms. Each single profile from the satellite height to the ionosphere peak has been fitted by a semi-Epstein layer function of the type used in the DGR model with shape factor variable with altitude. The topside above the satellite height has been extrapolated to match given values of plasmaspheric electron densities to obtain the full topside profile. The bottomside electron density has been calculated by using the maximum electron density and its altitude estimated from the topside ionogram as input for a modified version of the DGR derived profiler that uses model values for the foF1 and foE layers of the ionosphere. Total electron content has also been calculated. Longitudinal cross sections of vertical profiles from latitudes 50° N to 50° S latitude are shown for low and high geomagnetic activity. These cross sections indicate the equatorial anomaly effect and the changes of the shape of low latitude topside ionosphere during geomagnetic active periods. These results and the potentiality of the technique are discussed.     

Comparison of topside ionosphere scale height modeled by the Topside Sounder Model and incoherent scatter radar ionospheric model

The topside ionosphere scale height extracted from two empirical models are compared in the paper. The Topside Sounder Model (TSM) provides directly the scale height (H_T), while the incoherent scatter radar ionospheric model (ISRIM) provides electron density profiles and its scale height (H_R) is determined by the lowest gradient in the topside part of the profile. H_T and H_R are presented for 7 ISR locations along with their dependences on season, local time, solar flux F10.7, and geomagnetic index ap. Comparison reveals that H_T values are systematically lower than respective H_R values as the average offset for all 7 stations is 55 km. For the midlatitude stations Arecibo, Shigaraki, and Millstone Hill this difference is reduced to 43 km. The range of variations of H_R is much larger than that of H_T, as the H_T range overlaps the lower part of the H_R range. Dependences on ap, DoY and LT are much stronger in the ISRIM than in TSM. This results in much larger values of H_R at higher ap. Diurnal amplitude of H_R is much larger than that of H_T, with large maximum of H_R at night. The present comparison yields the conclusion that the ISR measurements provide steeper topside Ne profiles than that provided by the topside sounders.     

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

Ionospheric electron temperature at solar maximum

Langmuir probe measurements made at solar maximum from the Dynamics Explorer-2 satellite in 1981 and 1982 are employed to examine the latitudinal variation of electron temperature, T_e, at altitudes between 300 and 400 km and its response to 27 day variations of solar EUV. Comparison of these data with T_e models based on the solar minimum measurements from Atmosphere Explorer-C suggest that the daytime T_e does not change very much during the solar cycle, except at low latitudes where an especially large 27 day variation occurs. The 27 day component decreases from about 7°/F10.7 unit at the equator to 3°/F10.7 unit at 851V 3 middle and higher latitudes. From these DE-2 measurements, and those from AE-C, we conclude that the daytime T_e near the F₂ peak is more responsive to short-term (daily) variations in F10.7 than to any longer term changes that may occur between solar minimum and solar maximum. To investigate this sensitivity of the dayside ionosphere to solar activity we employ the inverse relationship of T_e and N_e, that was found at solar minimum, to see if it can be used to order the T_e behaviour at solar maximum. We introduce a simple quadratic correction for the F10.7 influence on T_e based on the entire daytime AE-C and DE-2 data base between 300 and 400 km. Although this equation may be found useful, the systematic deviations of the DE-2 data suggest that the solar minimum model does not accurately describe the T_e-N_e relationships at solar maximum, at least above 300 km where the DE-2 measurements were made. Future work with this data base should attempt to see if such a relationship exists.     

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.     

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 response to solar and magnetospheric protons during January 15–22, 2005: EAGLE whole atmosphere model results

We present an analysis of the ionosphere and thermosphere response to Solar Proton Events (SPE) and magnetospheric proton precipitation in January 2005, which was carried out using the model of the entire atmosphere EAGLE. The ionization rates for the considered period were acquired from the AIMOS (Atmospheric Ionization Module Osnabrück) dataset. For numerical experiments, we applied only the proton-induced ionization rates of that period, while all the other model input parameters, including the electron precipitations, corresponded to the quiet conditions. In January 2005, two major solar proton events with different energy spectra and proton fluxes occurred on January 17 and January 20. Since two geomagnetic storms and several sub-storms took place during the considered period, not only solar protons but also less energetic magnetospheric protons contributed to the calculated ionization rates. Despite the relative transparency of the thermosphere for high-energy protons, an ionospheric response to the SPE and proton precipitation from the magnetotail was obtained in numerical experiments. In the ionospheric E layer, the maximum increase in the electron concentration is localized at high latitudes, and at heights of the ionospheric F2 layer, the positive perturbations were formed in the near-equatorial region. An analysis of the model-derived results showed that changes in the ionospheric F2 layer were caused by a change in the neutral composition of the thermosphere. We found that in the recovery phase after both solar proton events and the enhancement of magnetospheric proton precipitations associated with geomagnetic disturbances, the TEC and electron density in the F region and in topside ionosphere/plasmasphere increase at low- and mid-latitudes due to an enhancement of atomic oxygen concentration. Our results demonstrate an important role of magnetospheric protons in the formation of negative F-region ionospheric storms. According to our results, the topside ionosphere/plasmasphere and bottom-side ionosphere can react to solar and magnetospheric protons both with the same sign of disturbances or in different way. The same statement is true for TEC and foF2 disturbances. Different disturbances of foF2 and TEC at high and low latitudes can be explained by topside electron temperature disturbances.     

SAMI2 model results for the quiet time low latitude ionosphere over India

Efficacy of SAMI2 model for the Indian low latitude region around 75°E longitudes has been tested for different levels of solar flux. With a slight modification of the plasma drift velocity the SAMI2 model has been successful in reproducing quiet time ionospheric low latitude features like Equatorial Ionization Anomaly. We have also showed the formation of electron hole in the topside equatorial ionosphere in the Indian sector. Simulation results show the formation of electron hole in the altitude range 800–2500?km over the magnetic equator. Indian zone results reveal marked differences with regard to the time of occurrence, seasonal appearances and strength of the electron hole vis-a-vis those reported for the American equatorial region.     

Ionospheric heating due to solar flares as measured by SROSS-C2 satellite

The data on thermal fluctuations of the topside ionosphere have been measured by Retarding Potential Analyser (RPA) payload aboard the SROSS-C2 satellite over the Indian region for half of the solar cycle (1995–2000). The data on solar flare has been obtained from National Geophysical Data Center (NGDC) Boulder, Colorado (USA) and other solar indices (solar radio flux and sunspot number) were download from NGDC website. The ionospheric electron and ion temperatures show a consistent enhancement during the solar flares. The enhancement in the electron temperature is 28–92% and for ion temperature it is 18–39% compared to the normal day’s average temperature. The enhancement of ionospheric temperatures due to solar flares is correlated with the variation of sunspot and solar radio flux (F10.7cm). All the events studied in the present paper fall in the category of subflare with almost same intensity. The ionospheric electron and ion temperatures enhancement have been compared with the IRI model values.     

An empirical model of electron density for low and middle latitudes below 200 km

Our empirical model of electron density (n_e) for quiet and weakly disturbed geomagnetic conditions (Kp not greater 4) takes account of comparative analysis of existing models and of experimental data obtained by rockets and incoherent scatter radar. The model describes the n_e distribution in the 80 to 200 km height range at low and middle latitudes, and to some extent, in the subauroral region. It is presented in analytical form thus allowing one to calculate electron density profiles for any time. The electron density distribution at 140 km depends on the season (day of the year) and on the solar zenith angle. Profile variations during the day are for one season shown. Different from other models, ours specifies the variations during sunrise and sunset and reflects the particular profile shape at night admitting the occurrence of an intermediate layer.     

An experimental and empirical model of electron temperature for altitudes of 500 to 1000 km and for a high solar activity period

On the basis of systematic electron temperature measurements onboard the Interkosmos-19 satellite, an experimental global model of electron temperature T_e has been constructed; namely, a set of samples representing 10 intervals of measured T_e, accompanied by values of the geographic longitude, solar zenith angle, season of the year, Covington index, Dst and Kp, grouped according to the invariant latitude, geomagnetic time and altitude. On the basis of the experimental model, the coefficients of the empirical models for the summer and winter seasons, for geophysically quiet conditions, and for heights of 520, 600, 920 and 1000 km are calculated. For heights of 680, 760 and 840 km with fewer data available, the coefficients are provisional.     

Studies of electron and ion temperatures at 500 km altitude during sunrise using Indian SROSS C2 satellite

Solar dependence of electron and ion temperatures (T_e and T_i) in the ionosphere is studied using RPA data onboard SROSS C2 at an altitude of ∼500 km and 77°E longitude during early morning hours (04:00–07:00 LT) for three solar activities: solar minimum, moderate and maximum during winter, summer and equinox months in 10°S–20°N geomagnetic latitude. In winter the morning overshoot phenomenon is observed around 06:00 LT (T_e enhances to ∼4000 K) during low-solar activity and to T_e ∼ 3800 K, during higher solar activity. In summer, it is observed around 05:30 LT, but the rate of T_e enhancement is higher during moderate solar activity (∼2700 K/hr) than the low-solar activity (∼1700 K/hr). During equinox, this phenomenon is delayed and is observed around 06:00 LT (∼4200 K) during all three activities.     

The global response of terrestrial ionosphere to the December 2015 space weather event

This paper investigates the ionospheric storm of December 19–21, 2015, which was initiated by two successive CME eruptions that caused a G3 space weather event. We used the in situ electron density (N_e) and electron temperature (T_e) and the Total Electron Content (TEC) measurements from SWARM-A satellite, as well as the O/N₂ observations from TIMED/GUVI to study the ionospheric impact. The observations reveal the longitudinal and hemispherical differences in the ionospheric response to the storm event. A positive ionospheric storm was observed over the American, African and Asian regions on 20 December, and the next day showed a negative storm. Both these exhibited hemispheric differences. A positive storm was observed over the East Pacific region on 21 December. It is seen that the net effect of both the disturbance dynamo electric field and composition differences become important in explaining the observed variability in topside ionospheric densities. In addition, we also discuss the T_e variations that occurred as a consequence of the space weather event.     

Ionosphere-plasmasphere electron fluxes at middle latitudes obtained from whistlers

985 whistlers observed between 1970 and 1975 in Hungary have been processed for equatorial plasmaspheric electron density and tube electron content above 1000 km (N_T). The hourly median value of N_T exhibits a diurnal variation with an amplitude of 1×10¹³ electrons/cm²-tube. 75 per cent of the electron flux values obtained from the time variation of N_T are lower than 6×10⁸ el cm^?2s^?1, while in some cases the fluxes reach a value as high as 3×10⁹ el cm^?2s^?1. Between 17 and 04 LT the dominant flux direction is toward the ionosphere. The data also indicate that the day to day filling of the plasmasphere after magnetic disturbances continues through several days without exhibiting saturation, with higher filling rates for lower values of average K_p.     

不同上边界条件下的极区电离层数值模拟

利用一维自洽的极区电离层模型,研究了沿磁力线方向不同电离层-磁层耦合条件下极区电离层的响应.此模型在110-610km的电离层空间区域内,综合求解描述极区电离层的连续性方程、动量方程和能量方程,以得到电离层数值解.研究发现,上边界条件在200 km以上的高度能显著地影响电离层参量的形态.较高的O+上行速度对应较低的F层峰值和较高的电子温度.不同边界O+上行速度对应的温度高度剖面完全不同.200km以上电子温度高度剖面不但由来自磁层的热流通量所控制,同时还受到场向O+速度的影响.对利用电离层模型研究电离层内部物理过程提出了建议.      

On the relative abundance of helium ions in the topside ionosphere

Based on data from satellite INTERCOSMOS-BULGARIA-1300, the latitudinal distribution of oxygen and helium ions in the topside ionosphere is discussed for night-time equinox at high solar activity. A comparison with the corresponding IRI-79 distribution is made. The vertical IRI ion composition profile is checked with measurements made with VERTICAL-10 rocket. Some recommendations are made in order to improve the IRI-modelling of the ion composition in the topside ionosphere.