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

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

基于CHAMP卫星观测数据对热层大气密度的经验正交分析

对2003年(太阳活动较高年)至2007年(太阳活动低年) CHAMP卫星的热层大气密度观测数据进行了经验正交函数(EOF)分析, 得到了400 km高度上白天平均大气密度ρ的太阳活动周变化与年度变化等不同变化分量. 研究结果表明, ρ受太阳活动影响较大, 其太阳周变化分量与F10.7指数变化之间的相关系数可高达94.5 %; ρ的太阳周变化分量随纬度增加而减小, 且在中高纬地区, 南半球的值明显大于北半球的值, 在低纬地区则出现基本对称的双峰分布, 即赤道质量密度异常(EMA)结构. 在ρ的年变化中, 呈现出明显的季节变化, 即夏季低冬季高; 同时ρ的年变化幅度随太阳活动增加而增强, 随纬度增加而增强. 将本文结果与经验模式NRLMSISE00在观测条件下的输出数据进行对比, 发现两者的太阳周变化与年变化分量基本一致, 但本文观测数据的太阳周成分随纬度变化略小, 年变化幅度略大, 且NRLMSISE00模式不能再现EMA结构. 研究结果对揭示热层气候学变化特征具有重要意义.      

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

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

低纬100—200km间电子数密度的变化

利用光化平衡模式计算了低纬100—200km间白天电子数密度的变化。求得E-F₁谷区的谷深,谷宽、谷高的变化特征。获得如下结果:a.太阳活动明显影响电子数密度随高度及太阳天顶角的变化,发现太阳活动指数与电子数密度间不仅存在正相关,而且存在负相关;b.太阳活动明显影响E-F₁谷区的形态。在一定太阳活动条件下,对同一太阳赤纬和地理纬度,谷深、谷宽与太阳天顶角的关系难以用一简单函数来表示;c.太阳耀斑、地磁活动对该区电子密度有明显影响;d.在讨论100—200km间电子密度时不能忽略O+（2P）和NO的光电离率。      

基于卫星加速度数据反演的热层大气密度与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%左右.      

北京和东北地区雷暴与太阳活动的关系

太阳活动对于大气电性能可能会产生影响,而大气电性能在雷暴形成过程中可能会起着重要作用.因此,太阳活动对于大气电性能的影响,最终将有可能影响到雷暴的形成.本文用1957年到1978年北京地区和东北地区的10个气象站的13000多个雷暴资料,用时序迭加法分析了雷暴数与太阳磁扇形界面通过、太阳耀斑及太阳黑子11年周期变化等的关系.结果表明,有些年分的雷暴数当行星际磁场方向由指向太阳变为离开太阳的磁扇形界面扫过地球时有相关影响,在界面通过前后3天内或7天后雷暴数明显增大.从季节来看,上半年的界面通过比下半年的界面通过对雷暴的影响大.在各种太阳耀斑分类统计分析中,看到出现在日面上西区特别是西一区(0°-30°)的耀斑对雷暴的影响显著.雷暴数与太阳黑子11年周期变化则没有很显著的相关关系.      

一种电离层垂直剖面物理模型及其与IRI-90模型的比较

在80-500km范围内考虑了3种中性成份的4种离子,从严格的电子和离子密度连续方程出发,对中性风和扩散效应进行了全面、连续的考虑,由此建立了一种电离层的物理模式;在此模式的基础上针对北京地区分别对太阳活动低年(F_10.7=60)、高年(F_10.7=300)的春(DOY=90)、夏(DOY=180)、秋(DOY=270)、冬(DOY=365)进行计算,并将所得结果与IRI-90进行了比较.结果表明: E层为典型的Chapman层: E-F谷区深度一般为0.2-0.5之间,比IRI要深;F₁缘在太阳活动低年的四季都出现,其中夏天最明显,已形成了一个F₁层,冬天最不明显,仅表现为一个轻微的凸缘,在太阳活动高年只有夏天出现了F₁凸缘,这与现有理论相符合,而IRI-90较少出现明显的F₁缘;F₂层的电子密度是活动高年比低年大,平均冬天比夏天大,这与观测结果也基本符合.      

银河宇宙线在电离层D层中电离的全球分布

本文从带电粒子对D层大气电离出发, 给出了宇宙线相对论粒子、非相对论粒子及低能粒子在地球大气中的电离公式, 从而给出了宇宙线在电离层D层中电子产生率q(h)和电子密度N(h)的全球分布.结果表明, 宇宙线产生的q(h)和N(h)具有明显的纬度效应, 在极区产生的q(h)和N(h)要比低纬高得多, 当截止刚度R_c=10—18GV时, q(h)的变化相差很小.太阳活动11年调制对q(h)的影响是明显的, 但远小于R_c对q(h)的影响.大气密度ρ(h)对q(h)的影响主要是随高度的变化.      

极区顶部电离层离子上行的太阳活动依赖性研究

利用第23太阳活动周DMSP F12,F13和F15卫星数据,分别对南北半球极区顶部电离层离子上行的太阳活动依赖性进行了研究.结果表明,南北半球上行事件对太阳活动的响应特征基本一致,即高（低）太阳活动时,离子上行通量以及上行数密度较大（小）,但是上行速度及上行发生率较低（高）.以南半球高纬为例,计算得到离子上行通量、数密度、速度及发生率在高低太阳活动条件下的比值分别约为2.26,3.35,0.71,0.51.对离子上行太阳活动依赖性的可能原因进行了分析.不同太阳活动水平下,光致电离及高能粒子沉降的差异会导致电离层离子密度的不同,而电离层离子密度的变化会改变离子elax-elax中性大气之间的碰撞频率,这是影响离子上行发生率的一个重要原因.      

NRLMSISE-00大气模型与GRACE和CHAMP卫星大气密度数据的对比分析

利用GRACE(Gravity Recovery And Climate Experiment)和CHAMP(Challenging Mini-Satellite Payload)卫星2002-2008年的大气密度数据与NRLMSISE-00大气模型密度结果进行比较,分析了模型密度误差及其特点.结果显示,NRLMSISE-00大气模型计算的密度值普遍偏大,其相对误差随经纬度变化,在高纬度相对较小;相对误差随地方时变化,在02:00LT和15:00LT左右较大,10:00LT和20:00LT左右较小.通过模型密度相对误差与太阳F_10.7指数的对比分析发现,在太阳活动低年模型相对误差最大,而在太阳活动高年相对误差较小;将模型结果分别与GRACEA/B双星和CHAMP卫星的密度数据进行比较,发现对于轨道高度更高的GRACE卫星轨道,模型相对误差更大;在地磁平静期,相对误差与地磁ap指数(当前3h)相关性不强,但是在大磁暴发生时,误差急剧增大.      

Indications of a quasi-20-year cycle of middle atmosphere temperatures

A homogeneous series of 25 years, 1959–1983, of daily measurements of low-frequency radio wave reflection heights in the lower ionosphere (around 80 km), at constant zenith distance of the Sun, has been analysed. After removing the 11-yr solar cycle variation from these data by means of empirical regression coefficients with the solar activity index, F_10.7, a significant residual variation remains with a maximum in 1965 and a minimum in 1975. This residual can be interpreted in terms of a corresponding non-solar variation of neutral air pressure at 80 km, thus indicating that recent climatic temperature changes in the middle atmosphere are of quasi-cyclic character rather than a monotonous trend.     

一氧化氮冷却率的变化及其作用

利用最新的NO经验公式,计算了近两个太阳活动周期100-200km间NO的冷却率,研究了太阳和地磁活动对NO和O冷却率的影响。进一步证明了NO冷却是120km以上热层的主要冷却过程。      

太阳活动突发过程对临近空间大气影响的模拟

空间天气对地球及近地空间具有重要影响,大的空间天气事件对中上层大气动力学和成分具有不同的影响。利用全大气耦合模式WACCM,针对太阳耀斑、太阳质子、地磁暴三类事件,以太阳活动平静期2015年5月10-14日的GEOS-5数据为模式背景场,通过F_10.7、离子产生率、Kp及Ap指数设置,分别模拟三类事件对临近空间大气温度、密度和臭氧的影响。结果表明耀斑事件在三类事件中对临近空间大气温度和密度的影响最为显著。平流层大气温度增加是由耀斑辐射增强引起平流层臭氧吸收紫外辐射发生的光化学反应所致,耀斑事件引起平流层和低热层温度增加约为2～3 K,低热层大气相对密度增加在6%以内;太阳质子事件及磁暴事件主要影响低热层,但太阳质子事件和磁暴事件对低热层温度扰动不大于1 K。     

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.     

On the diurnal variations in the temperature and composition: a three-dimensional model with superrotation

Based on a simplified theoretical interpretation of the composition measurements with the ONMS and OIMS experiments on Pioneer Venus, the conclusion was drawn that the rotation rate of the thermosphere should be close (within a factor of two) to that of the lower atmosphere. A more realistic three-dimensional model of the thermosphere dynamics is now being developed, considering non-linear processes, higher order modes and collisional momentum exchange between the major species CO₂, CO and O, which describes the diurnal variations in temperature and composition (Niemann et al., $\text{JGR}$, 1980). The computed horizontal winds are about 300 m/sec near the terminators and poles. Results are also presented from a two-dimensional (quasi-axisymmetric) spectral model which describes the four day superrotation in the lower atmosphere of Venus.     

The response of the high-latitude thermosphere to geomagnetic substorms

Enhanced ionization and electric fields associated with geomagnetic substorms lead to time-varying momentum (ion drag) and heat (Joule dissipation) sources in the high latitude thermosphere. These momentum and heat sources, in turn, cause changes in the structure of the wind, pressure and temperature fields, and, in addition, provide a mechanism for the generation of atmospheric gravity waves. We have developed a model which has been used to simulate, on a spatial scale of tens of kilometers, the response of the high latitude thermosphere to these momentum and heat sources, which constitute a major factor in the coupling between the magnetosphere, the ionosphere and the upper atmosphere. Simulations reported herein indicate that, because the coupling between ions and neutrals is highly dependent on the neutral-ion collision frequency, the response of the neutral atmosphere varies, not only as a function of total ionization, but also as a function of the details of the electron density profile.     

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.     

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.