MSIS-86和IRI-86模式的扩展应用

本文简要地介绍了MSIS-86和IRI-86两个模式,并用此两模式所给的参数推算出热层中性大气风系、碰撞频率和电离层电导率的时空分布。      

用经验模式参数研究电离层发电机效应

本文利用MSIS-86和IRI-86模式的基本参数,求得热层风系和电导率的三维分布;继而从发电机理论得到了电离层电位及层电流密度在北半球的分布和变化。本研究将大量数据统计平均的大气模式参数与热层、电离层理论研究联系了起来。     

基于电子密度同化的热层大气密度预报方法

采用热层电离层耦合模式TIEGCM和集合卡尔曼滤波同化方法,利用同化COSMIC电离层掩星电子密度数据优化热层电离层参量,并将模式预报的大气密度与CHAMP卫星大气密度数据进行对比,分别开展模拟和实测数据的同化预报实验.在模拟数据同化实验中,状态向量包含温度、风场和离子成分的实验结果表明,仅优化温度即可达到最优的热层大气密度预报效果.在实测数据同化实验中,将温度作为状态向量参数,优化结果表明,循环同化过程中模式预报的大气密度相对偏差的均方根误差在48h内从38%减小到27%,同化稳定时间至少需要30h.预报过程中大气密度预报效果的改善持续时间为34h.这表明电子密度同化能够改善热层大气密度的预报精度,设计的实验方案合理可行,可获得较长的预报时效.      

电离层水平电场与风场的耦合模拟研究

设计了一个将电离层水平电场与风场耦合的模拟方案,研究了电流函数和风场在耦合前后的变化与差异. 研究发现,水平电场与风场相互反馈后,风场的变化比电流函数小. 经向风在白天有较明显的差异,夜晚的差异比白天小,主要出现在中高纬地区,并随高度的增加而增大,300km左右达到最大值,其后几乎保持不变. 纬向风有与经向风相似的变化,但纬向风耦合前后的差异比经向风小. 电流函数在耦合后有较大改变,两个涡旋强度都有较强增加,并且北半球的增强大于南半球,而夜晚差异较小. 结果表明,在研究的高度范围内,风场对电场的控制作用大于电场对风场的影响.      

Mathematical modeling of plasma drifts over equatorial low latitude regions

This paper presents a mathematical model to simulate ionospheric plasma drifts at equatorial low latitude regions by coupling of E- and F-regions. The governing non-linear differential equations (of elliptic and parabolic nature) are solved numerically through finite-difference schemes and obtained neutral winds and electric fields. The temperature and electron density profiles are generated utilizing MSIS-86 atmospheric model. The continuity equation is employed to obtain night-time E-region density profile using measured ionograms at Trivandrum (India). The computed vertical and zonal plasma drifts are comparable with measured Jacamarca plasma drifts with little variations during noon and evening times. The plasma drifts at Trivandrum (8.5° N, 76.5° E, dip 0.5° N) are compared with those of Jicamarca (12° S, 76.9° W, dip 2° N). Neutral wind simulations of present model agree well with those of horizontal wind model (HWM-93). The post-sunset enhancement and its reversal are also discussed.     

太阳活动对EUV辐射在大气中吸收过程的影响

本文利用A-E卫星在太阳活动21周峰年间观测到的EUV辐射资料,高层大气成分的吸收截面,以及MSIS-86热层大气模式,研究了EUV辐射在大气中的吸收过程;在透射比为1/e和0.1/100时分别计算了透射高度随波长及太阳活动的变化。在波长范围50—1050内对37个波段分别求出了透射高度随太阳天顶角的变化。结果表明,当太阳活动增强时各波段的透射高度均升高,而且透射比越大则透射高度随太阳活动的变化也越剧烈。当透射比为一定时,太阳天顶角越大则透射高度随太阳活动的变化也越大。除此之外还存在一个相反效应,即太阳活动会使Chapman函数变小,这反过来又促使透射高度降低。这两种效应的综合作用结果可较好地解释某些电离层观测中的日没效应。     

Thermospheric storms

Severe perturbations of the neutral upper atmosphere caused by the dissipation of solar wind energy have become known as thermospheric storms. This review summarizes the temperature and density changes observed during such events. The following topics are discussed: (1) the morphology and origin of composition changes at high and middle latitudes; (2) low-latitude effects; (3) time delays between magnetic and thermospheric perturbations; (4) systematic variations; and (5) accuracy of model predictions.     

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.     

On the structure and dynamics of the thermosphere

Thermospheric temperature, composition and wind measurements from the Dynamics Explorer satellite (DE-2) are interpreted using a three dimensional, multiconstituent spectral model. The analysis accounts for tides driven by the absorbed solar radiation as well as energy and momentum coupling involving the magnetosphere and lower atmosphere. We discuss phenomena associated with the annual tide, polar circulation, magnetic storms and substorms.     

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

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

Vertical plasma transport in the ionosphere over Irkutsk during St. Patrick’s Day geomagnetic storm: Observations and modeling

The St. Patrick’s Day storm being the strongest geomagnetic storm of Solar Cycle 24 caused strong changes in ionospheric and thermospheric dynamics. The paper presents a study of vertical plasma transport in the ionosphere during the St. Patrick’s Day storm with using both observations and modeling. The observations give the ionospheric peak height obtained with the chirp vertical sounding ionosonde and the neutral wind velocities obtained with the Fabry-Perot interferometer. The ionospheric peak height is an indicator of the total vertical plasma transport, while meridional wind and electromagnetic drift are the two main drivers of the vertical plasma transport. The Global Self-consistent Model of the Thermosphere, Ionosphere, and Protonosphere used in this study gives the total set of ionospheric and thermospheric parameters including F2-layer peak height, neutral wind velocities, electric field, and neutral composition. The model/data comparison allows us to obtain two main results. The first one is an estimation of the model prediction possibilities under storm conditions. The second result is an indirect assessment of the neutral wind and electric field contribution into the changes in the ionospheric peak height in the case of the St. Patrick’s Day geomagnetic storm.     

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.     

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

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

Neutral winds in the polar thermosphere as measured from dynamics explorer

Remote sensing measurements of the meridional thermospheric neutral wind using the Fabry-Perot Interferometer on Dynamics Explorer have been combined within-situ measurements of the zonal component using the Wind and Temperature Spectrometer on the same spacecraft. The two data sets with appropriate spatial phasing and averaging determine the vector wind along the track of the polar orbiting spacecraft. A study of fifty-eight passes over the Southern (sunlit) pole has enabled the average Universal Time dependence of the wind field to be determined for essentially a single solar local time cut. The results show the presence of a “back-ground” wind field driven by solar EUV heating upon which is superposed a circulating wind field driven by high latitude momentum and energy sources.     

基于卫星加速度数据反演的热层大气密度与NRLMSISE-00模式结果的比较研究

利用CHAMP/STAR加速度数据反演的热层大气密度与NRLMSISE-00模式反演的热层大气密度进行比较, 结果表明, 热层大气密度在春秋季期间高于冬夏季, 并且太阳活动高年比低年更加显著; 日照面和阴影区大气密度的比值在低纬地区由太阳活动高年的4下降到太阳活动低年的2左右, 中纬地区大约由3变化到1.5, 高纬地区变化较小; NRLMSISE-00模式能够较好地模拟热层大气密度的变化趋势, 但是磁暴期间模式精度较差. 统计结果表明, 模式整体比反演结果偏高, 2002-2008年相对偏差分别为16.512%, 20.004%, 18.915%, 18.245%, 25.161%, 33.261%和41.980%; NRLMSISE-00模式在高纬地区的相对偏差为27.337%, 高于中低纬地区的24.047%; 模式在中等太阳活动水平相对偏差较为稳定, 基本在15%左右.      

High accuracy satellite drag model (HASDM)

The dominant error source in force models used to predict low-perigee satellite trajectories is atmospheric drag. Errors in operational thermospheric density models cause significant errors in predicted satellite positions, since these models do not account for dynamic changes in atmospheric drag for orbit predictions. The Air Force Space Battlelab’s High Accuracy Satellite Drag Model (HASDM) estimates and predicts (out three days) a dynamically varying global density field. HASDM includes the Dynamic Calibration Atmosphere (DCA) algorithm that solves for the phases and amplitudes of the diurnal and semidiurnal variations of thermospheric density near real-time from the observed drag effects on a set of Low Earth Orbit (LEO) calibration satellites. The density correction is expressed as a function of latitude, local solar time and altitude. In HASDM, a time series prediction filter relates the extreme ultraviolet (EUV) energy index E_10.7 and the geomagnetic storm index a_p, to the DCA density correction parameters. The E_10.7 index is generated by the SOLAR2000 model, the first full spectrum model of solar irradiance. The estimated and predicted density fields will be used operationally to significantly improve the accuracy of predicted trajectories for all low-perigee satellites.     

Horizontal and vertical thermospheric cross-wind from GOCE linear and angular accelerations

Thermospheric wind measurements obtained from linear non-gravitational accelerations of the Gravity field and steady-state Ocean Circulation Explorer (GOCE) satellite show discrepancies when compared to ground-based measurements. In this paper the cross-wind is derived from both the linear and the angular accelerations using a newly developed iterative algorithm. The two resulting data sets are compared to test the validity of wind derived from angular accelerations and quantify the uncertainty in accelerometer-derived wind data. In general the difference is found to be less than 50?m/s vertically after high-pass filtering, and 100?m/s horizontally. A sensitivity analysis reveals that continuous thrusting is a major source of uncertainty in the torque-derived wind, as are the magnetic properties of the satellite. The energy accommodation coefficient is identified as a particularly promising parameter for improving the consistency of thermospheric cross-wind data sets in the future. The algorithm may be applied to obtain density and cross-wind from other satellite missions that lack accelerometer data, provided the attitude and orbit are known with sufficient accuracy.     

A global study of meridional winds in the thermosphere

We present neutral meridional winds derived from existing ground-based measurements of the height of the F2-layer maximum electron density (h_max). The method of calculation uses an ionospheric model to determine the relationship between h_max and the neutral wind along a magnetic meridian. The meridional wind is derived from a comparison of the modeled layer height with the measured height. This method is used in the global mapping of thermospheric winds using radar and ionosonde measurements from the Global Thermospheric Mapping Study at the summer and winter solstices. The diurnal behavior of the meridional wind is calculated for the two seasons from Millstone Hill radar measurements and from several ionosonde stations. The major features of the winds calculated from Millstone Hill radar data are a 200 m/s southward wind between 0200 and 0400 hours local time at the summer solstice and a northward daytime wind in winter. A sampling of winds derived from ionosonde data shows the diurnal pattern to vary with geographic latitude and longitude. Nighttime equatorward winds are found to be larger in Europe than at other locations of similar latitude.     

黄河站和Tromso站亚暴期间热层风场观测结果

中高层大气风场探测对研究大气物理过程具有极为重要的意义，尤其是在极地地区，风场对大气结构的影响更为剧烈.针对亚暴期间中国北极黄河站和日本Tromso站上空OI557.7nm气辉层（低热层）中性风场，利用全天空法布里-珀罗干涉仪（all-sky Fabry-Perot Interferometer，all-sky FPI）探测气辉谱线的多普勒频移，反演气辉层的大气风场信息.结果表明，低热层风场平均水平在100m·s-1左右，热层风场在极地地区更为剧烈，纬度相对较低的Tromso站探测到的风速整体小于同期黄河站上空的风速.结合离子风数据，分析离子拖拽和焦耳加热对中性风的影响过程，发现极光亚暴不仅对低热层风场有增强作用，也有明显的抑制效果，但整体风向都垂直于极光弧变化.