Variations of Thermosphere Density Based on Joint Analysis of In-Situ Measurement Data From Shenzhou Spacecrafts

Based on the measurements made by Atmospheric Density Detectors (ADDs) onboard Chinese spacecraft Shenzhou 2-4, the variations of thermosphere density are revealed. During the quiet period, the density at spacecraft altitude of 330～410km exhibited a dominant diurnal variation, with high value on dayside and low value on nightside. The ratio of the diurnal maximum density to the minimum ranged from 1.7 to 2.0. The ratio shows a positive correlation with the level of solar activity and a negative correlation with the level of geomagnetic activity. When a geomagnetic disturbance comes, the atmospheric density at the altitude of 330～410km displayed a global enhancement. For a strong geomagnetic disturbance, the atmospheric density increased by about 56%, and reached its maximum about 6～7 hours after the geomagnetic disturbance peak. The density asymmetry was also observed both in the southern and northern hemisphere during the geomagnetic disturbance peak.      

利用COSMIC数据构建的全球垂直TEC图的验证研究

在日固坐标系(地磁纬度和地方时)下, 累积地方时过去24h的COSMIC(Constellation Observing System for Meteorology Ionosphere andClimate)观测资料, 通过对110$sim$750km高度范围内的电子密度进行数值积分得到各掩星点的垂直TEC值, 进而利用Kriging方法插值产生近实时的全球地方时MAGLat2.5°×2h的COSMIC TEC图. 利用2008年1月1日至2010年6月30日共30个月的COSMIC数据, 逐日构建COSMICTEC图, 将其与全球导航卫星系统服务组织(International GNSS Service,IGS)发布的全球电离层TEC图(Global Ionospheric Maps, GIMs)以及OSTM/JASON-2卫星高度计观测值分别进行比对,证明利用COSMIC掩星资料构建全球电离层垂直TEC图是可行的.      

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.     

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

Effects of geomagnetic storms on the bottomside ionospheric F region

The geomagnetic storm is a complex process of solar wind/magnetospheric origin. The variability of the ionospheric parameters increases substantially during geomagnetic storms initiated by solar disturbances. Various features of geomagnetic storm act at various altitudes in the ionosphere and neutral atmosphere. The paper deals with variability of the electron density of the ionospheric bottomside F region at every 10 km of altitude during intense geomagnetic storms with attention paid mainly to the distribution of the F1 region daytime ionisation. We have analysed all available electron density profiles from some European middle latitude stations (Chilton, Pruhonice, Ebro, Arenosillo, Athens) for 36 events that occurred in different seasons and under different levels of solar activity (1995–2003). Selected events consist of both depletion and increase of the F2 region electron density. For European higher middle and middle latitude the F1 region response to geomagnetic storm was found to be negative (decrease of electron density) independent on the storm effect on the F2 region. For lower middle latitude the F1 response is weaker and less regular. Results of the analysis also show that the maximum of the storm effect may sometimes occur below the height of the maximum of electron density (NmF2).     

Properties of the mesosphere and thermosphere and comparison with CIRA 72

Exospheric temperatures of several reference atmosphere are reviewed and a recommendation is made for the exospheric temperature of a proposed mean CIRA. One of the deficiencies of CIRA 72 and other present thermospheric models is the representation of density changes with geomagnetic activity. This deficiency is illustrated with samples of data. The data show the effects of geomagnetic activity, particle precipitation, a solar proton event, and gravity waves. An empirical model developed from the unique AFGL satellite density data bank using multiple linear regression is reviewed. The present model is for low to moderate solar flux and quiet geomagnetic conditions, but it is planned to extend the model to active conditions. Good progress has been made since CIRA 72 was specified in our knowledge and understanding of the properties of the lower thermosphere, although there are still some unresolved problems. The biggest progress has been made in the theory of tidal effects and of particulate energy deposition and of electrojet heating. On the other hand, it is still not possible to define adequately the systematic variations of the lower boundary conditions of thermospheric models. This is due to lack of knowledge of the systematic variations of the structure properties in the 100 to 120 km altitude region and inadequate information on the mesospheric turbulence profile and variations in the turbopause altitude.     

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.     

F2 layer response to geomagnetic disturbances in Eastern Asia under the low solar activity

In this paper, the peculiarities of ionospheric response to geomagnetic disturbances observed at the decay and minimum of solar activity (SA) in the period 2004–2007 are investigated with respect to different geomagnetic conditions. Data from ionospheric stations and results of total electron content (TEC) measurements made at the network of GPS ground-based receivers located within the latitude–longitude sector (20–70°N, 90–160°Е) are used in this study. Three groups of anomalous ionospheric response to geomagnetic disturbances have been observed during low solar activity. At daytime, the large-scale traveling ionospheric disturbances (LSTIDs) could generally be related to the main phase of magnetic storm. Quasi-two-days wavelike disturbances (WLDs) have been also observed in the main phase independent of the geomagnetic storm intensity. Sharp electron density oscillations of short duration (OSD) occurred in the response to the onset of both main and recovery phases of the magnetic storm in the daytime at middle latitudes. A numerical model for ionosphere–plasmasphere coupling was used to interpret the occurrence of LS TIDs. Results showed that the LSTIDs might be associated with the unexpected lifting of F2 layer to the region with the lower recombination rate by reinforced meridional winds that produces the increase of the electron density in the F2 layer maximum.     

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

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

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.     

Intercomparisons of temperature,density and wind measurements from in situ and satellite techniques

Intercomparisons between satellite retrieved temperatures (TIROS N series) and those derived from radiosonde and rocketsonde profiles have been made covering the years 1980–1984. Differences in the measurement parameters between 100 and 0.4 mbar (～16–55 km) are described; generally radiosonde/satellite differences are less than 1°K, while rocketsonde/satellite differences reach 7–8°K in the upper stratosphere. Comparisons between the various $\text{in situ}$ devices indicate that radiosonde/rocketsonde differrences are less than 1°K while precision studies of the rocketsonde instrument find that the rocketsonde measurements are internally consistent to less than 1°K up to 50 km and to less than 3°K to 60 km. Density data obtained with the small rocketsondes ($\text{in situ}$ thermistors and inflatable spheres) and with the large sounding rocket systems show that density measurements usually agree to within 15 percent up to 85 km. Comparisons of the various atmospheric parameters obtained from different instruments are important, however the usefulness of intermixing the measurements is obvious and increased emphasis should be placed on procedures for intermingling such data. Suggestions are made on how this might be accomplished.     

Middle atmosphere models and comparison with shuttle reentry density data

Several middle atmosphere models will be reviewed, including a new set of models produced by Groves in 1985. The latter models are based on rocket and rawinsonde in situ measurements and satellite remote sounding temperature data. The models are compared with measurements made with instruments on board U. S. Shuttles during their reentry. Very useful atmospheric density data have been obtained in the altitude region from 50 to 80 km. The measurements are unique in that they are made by a vehicle travelling almost horizontally through the atmosphere at a velocity of 6 to 7 km/sec. This results in measurements along a path of approximately 8,000 km in a time interval of about 20 minutes. The results show some unique features.     

Properties of the neutral density and composition in the thermosphere

Measurements of the density and composition of the thermosphere between 150 and 500 km, which were obtained by the S3-1 satellite, have been compared with the Jacchia and MSIS models. The measurements of the densities of O, N₂, N and Ar show some differences from the current models which should be considered during the preparation of the next CIRA model. The Ar measurements are particularly useful in examining the response of the neutral atmosphere to geomagnetic heating. These results are useful in establishing the appropriate lower boundary conditions for modeling of the thermosphere.     

Temperature and altitude of the polar mesopause in summer

A review and summary of 60 in situ experiments is provided which determined the temperature and altitude of the mesopause north of 58°N latitude during the summer months of May through August. These experiments employed 4 experimental techniques; acoustic grenades, rigid and inflatable falling spheres, and Pitot-static tubes. Excellent agreement is found among the results obtained from different techniques. During June and July the average mesopause temperature drops below 130 K, the average mesopause altitude is 88.5 km. The climatological tables of CIRA 1986 indicate, however, a mean mesopause temperature of approximately 140 K at 91 km for corresponding geophysical conditions.     

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.     

Simultaneous investigation of magnetospheric plasma and spacecraft charging

The results of magnetospheric plasma fluxes measurement in geosynchronous orbit are presented. The measurements were done onboard the geosynchronous Russian “Electro” spacecraft (SC) in 1995–1997 years. The diurnal variations of the plasma density and temperature were observed at various levels of geomagnetic activity.     

强磁暴期间TIEGCM模式与CHAMP卫星热层大气密度的比较分析

利用NCAR-TIEGCM模式计算了2003年11月20—21日强磁暴期间410km高度上的大气密度,并与CHAMP/STAR加速度计反演数据进行对比和分析. 结果表明,模式结果能够准确反映磁暴期间大气密度的分布和变化情况,与实测结果在变化趋势和量级上具有较好的一致性,但在精细结构和数值大小上仍存在一定差异. 模式低估了磁暴期间大气密度的增幅,实测大气密度增幅高达250%～400%,而模式结果为100%～125%. 模式结果与实测数据的偏差在高纬地区高于低纬地区,日侧高于夜侧. 通过模式和实测数据的分析发现,磁暴期间大气密度扰动具有日夜侧和南北半球不对称性. 此外,模式能够准确反映磁暴期间大气密度扰动从高纬向低纬的传播以及大气密度对SYM-H指数响应的延迟特性.      

Recent advances in upper atmospheric structure

A possible quantitative explanation of the semi-annual variation in thermospheric density has been obtained in terms of a semi-annual variation in the computed globally averaged vertical energy carried by propagating tides from the lower and middle atmosphere into the thermosphere. The effect is primarily due to seasonal changes in the distribution of water vapor and in the solar declination angle and Sun-Earth distance. An MSIS-83 empirical model of the thermosphere, representing a revision of the earlier MSIS models, has been prepared. The database used covers a wider range of solar activity than previous models and an improved magnetic storm representation is included. Atomic oxygen profiles in the 100 to 160 km altitude region of the auroral thermosphere have been recalculated from measured quenching of N₂(A³∑_u⁺) using the latest laboratory rates and the results are in good agreement with the mean CIRA 1972 profile. A new empirical model of thermospheric variations with geomagnetic activity has been developed incorporating variations with local magnetic time, latitude dependent terms which can vary with the magnitude of the geomagnetic disturbance, and an altitude dependent expression for the equatorial wave. A new index ML, derived from the AL index, has been developed that appears to have promise to represent the variations of thermospheric species with geomagnetic activity. Satellite measured values of solar UV flux, ground-based observations of CaK plages, sunspot numbers and 10.7 cm solar radio flux have been analyzed for temporal variations. Some differences have been identified and the significance to empirical and theoretical upper atmosphere models is discussed.     

太阳活动与热层大气密度的相关性研究

为分析太阳活动对热层大气的影响,使用250km,400km,550km高度处热层大气密度与太阳F_10.7指数数据,研究了二者的周期变化及相关关系. 结果表明,热层大气密度的变化与太阳活动呈现相似的变化趋势;两者均具有显著的27天及11年周期变化特征,热层大气密度还存在7～11天及0.5年和1年的变化特征,且高度越高越明显;热层大气密度对太阳活动的最佳响应滞后为3天,无论何种地磁活动水平下,400km高度处相关性高于250km,550km处相关性最小,且太阳活动下降相期间高于上升相;250km,400km和550km高度处热层大气密度和太阳活动的统计结果分别为饱和、线性和放大关系;高度越高的热层大气密度对太阳活动响应越敏感.