基于钠激光雷达观测数据的中间层顶大气温度分布特性研究

作为中间层和热层的边界层,中间层顶存在多种能量交换方式,是大气能量耦合的重要区域。本文利用部署于中国科学院廊坊临近空间大气探测站的钠荧光多普勒激光雷达2013年的观测数据,研究了廊坊上空中间层顶区域大气温度的年度和季节分布特性,并分析了影响温度分布的多种因素。年平均温度廓线图显示,中间层顶位于约97.5 km高度处,温度约191.2 K。受放热化学反应的影响,年平均温度廓线91 km高度处出现了一个198 K的相对温度高点。中间层顶区域大气温度的季节分布受太阳辐射和大气动力学因素综合影响,夏季在大气动力学影响下,中间层顶高度较低,位于88 km高度处,温度也较低,约177 K;冬季太阳辐射起主导作用,中间层顶位于99 km高度处,温度为181 K。通过拟合月平均温度分析了中间层顶区域大气温度年变化和半年变化的振幅和相位特征。结果显示,中间层顶区域上部温度分布主要受太阳辐射的影响;在中间层顶区域下部,大气波动主导了温度分布。      

Semi-empirical model of middle atmosphere wind from the ground to the lower thermosphere

During recent years, special attention has been paid to understanding the background circulation of the middle atmosphere. Particularly in the mesosphere/lower thermosphere (MLT) region, this has involved including data from a range of new radar measurements. It has also involved the comparison of existing empirical middle atmosphere wind models, such as CIRA-86 and HWM-93 to the new data. This has led to the construction of empirical models of MLT winds such as the Global Empirical Wind Model (GEWM). Further investigations are aimed at the construction of new empirical and semi-empirical wind models of the entire middle atmosphere including these new experimental results. The results of a new wind climatology (0–100 km) are presented here, based upon the GEWM, a reanalysis of stratospheric data, and a numerical model which is used to fill the gap between data from the stratospheric and MLT regions.     

Model calculations of middle atmospheric wind and temperature fields

The zonally-averaged circulation of the middle atmosphere between 10 and 110 km is calculated from a two-dimensional grid point model for steady solstice conditions. Using the temperatures of CIRA 72 in a diagnostic study of the vertical-meridional circulation we have investigated the influence of a latitude dependent static stability and the consequences of meridional temperature advection. In the summer mesopause region we found an area of subsidence which is located near the equatorward boundary of the noctilucent cloud region. In a self-consistent dynamical study we have assumed an empirical momentum source in the lower summer thermosphere in addition to the well-known Rayleigh friction to obtain an eastward zonal flow in this region.     

Chapman函数在实际模式大气中的数值积分计算

本文对Chapman阳光掠射函数[Ch(z_p,χ)]进行了数值积分,求得了其在实际模式大气中随观测高度z_p及天顶角χ的变化.计算并讨论了低热层高温度梯度、分子与湍流扩散、重力场及太阳活动对 Ch(z_p,χ)的影响.结果表明,在150 km以下,Ch(z_p,χ)与前人用等标高模式及等标高梯度模式的计算结果差别较大.其中高温度梯度的影响起主导作用.特别是太阳活动对Chapman函数影响较显著,高年与低年之间可变化10—40％(在大天顶角时),这有可能推动热层大气中辐射-光化学-动力学耦合关系的变化.     

Effects of sector structure of the interplanetary magnetic field on the upper mesosphere–lower thermosphere dynamics

In this paper we study the influence of the interplanetary magnetic field (IMF) polarity changes caused by the Earth passing through the IMF sector boundary on the dynamic processes taking place in neutral atmosphere within the altitude interval of the upper mesosphere–lower thermosphere (83–101 km). The analysis has revealed the influence of the IMF sector structure on dynamics of the upper mesosphere–lower thermosphere. There has been a significant seasonal variation of the wind reaction to the IMF polarity changes observed. The influence of the IMF polarity changes on neutral atmosphere dynamics within the altitude range of 83–101 km is most pronounced in the zonal component of neural wind when the IMF polarity changes from negative to positive in all the seasons except for spring and when IMF polarity changes from positive to negative – in spring only.     

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.     

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.     

利用HRDI/UARS资料分析东亚区域中层大气纬向风气候特征

利用美国高层大气研究卫星(UARS)搭载的高分辨率多普勒测风仪(HRDI)获得的中层大气风场观测资料,对东亚区域中层大气纬向风的垂直分布与变化特征进行了分析研究.多年平均结果显示东亚区域中层大气纬向风具有显著的区域特征,与当前普遍使用的参考大气CIRA-86相比存在显著的不同.在冬季,东亚区域中间层西风急流中心位于25°-35°N之间的75 km高度,与CIRA-86相比,该中心纬度偏南5°,高度偏高10 km;在秋季,东亚区域低热层高度存在一个显著的从赤道到高纬度的东风带,而CIRA-86不存在.分析结果还表明,除了夏季中纬度地区,在东亚区域上空中高层大气各高度上均存在相当显著的区域尺度扰动结构.在热带,低热层高度纬向风无论冬夏,沿纬圈方向都表现出相当显著的不均匀性,夏季这种不均匀性进一步向下扩展到55 km高度.与上述热带扰动特征相比,中纬度地区夏季的纬向风在各个高度沿纬圈相当均匀,但是在冬季,中间层和低热层高度都存在沿纬圈方向显著的纬向风扰动结构.      

Neutral upper atmospheres of Venus and Mars

Numerous measurements of the neutral upper atmosphere above 100 km have been made from spacecraft over Venus and over Mars. The Venus exospheric temperatures are unexpectedly low (less than 300°K near noon and less than 130°K near midnight). These very low temperatures may be partially caused by collisional excitation of CO₂ vibrational states by atomic oxygen and partially by eddy cooling. The Venus atmosphere is unexpectedly insensitive to solar EUV variability. On the other hand, the Martian dayside exospheric temperature varies from 150°K to 400°K over the 11-year solar cycle, where CO₂ 15-μm cooling may be less effective because of lower atomic oxygen mixing ratios. On Venus, temperature increases with altitude on the dayside (thermosphere), but decreases with altitude from 100 to 150 km on the nightside (cryosphere). However, dayside Martian temperatures near solar minimum for maximum planet-sun distance and low solar activity are essentially isothermal from 40 km to 200 km. During high solar activity, the thermospheric temperatures of Mars sharply increase. The Venus neutral upper atmosphere contains CO₂, O, CO, C, N₂, N, He, H, D and hot nonthermal H, O, C, and N, while the dayside Mars neutral upper atmosphere contains CO₂, O, O₂, CO, C, N₂, He, H, and Ar. There is evidence on Venus for inhibited day-to-night transport as well as superrotation of the upper atmosphere. Both atmospheres have substantial wave activity. Various theoretical models used to interpret the planetary atmospheric data are discussed.     

Properties of the mesosphere and lower thermosphere

The variability and systematic variations of the properties of the upper mesosphere and lower thermosphere are probably the least well known aspects of the terrestrial atmosphere. Satellite measurements of this region are very limited and rocket and remote sounding techniques do not provide comprehensive coverage. Progress is being made in theoretical studies of this region, primarily with regard to tidal effects, and some progress is being made in analyzing the relatively sparse experimental data that are available. Turbulence dynamics of the region has been studied by analyzing structure measurements at Kwajalein, wind data from Natal and systematic variations of the turbopause altitude determined from measurements of the diffusive separation of argon. One question that is being raised at this time, and it is appropriate at a time near solar maximum, is the extent of solar activity control of the properties of this region of the atmosphere. The occurrence rates and magnitudes of the turbulent diffusivity in the 70 to 90 km altitude region appear to correlate with solar activity with a time lag, as do also the incidence of aurora and the atomic oxygen green line intensity. Solar cycle dependence has been identified in mean zonal wind speeds in the 65 to 110 km altitude region above Saskatoon and in lower thermosphere temperatures measured at Heiss Island and at St. Santin. Millstone Hill data show that the mean meridional wind changes during a solar cycle. Solar cycle variations have also been detected in the stratosphere and troposphere.     

Reevaluation of thermosphere heating by auroral electrons

This paper presents a new calculation of neutral gas heating by precipitating auroral electrons. It is found that the heating rate of the neutral gas is significantly lower than previous determinations below 200 km altitude. The neutral gas heating arises from the many exothermic chemical reactions that take place from the ions and excited species created by the energetic electrons. The calculations show that less than half the energy initially deposited ends up heating the neutral gases. The rest is radiated or lost in the dissociation of O₂ because the O atoms do not recombine in the thermosphere. This paper also presents a new way of calculating the heating rate per ionization that can be used for efficient determination of the overall neutral gas heating for global thermosphere models. The heating rates are relatively insensitive to the neutral atmosphere when plotted against pressure rather than altitude coordinates. At high altitudes, the heating rates are sensitive to the thermal electron density and long-lived species. The calculations were performed with the Field Line Interhemispheric Plasma (FLIP) model using a 2-stream auroral electron precipitation model. The heating rate calculations in this paper differ from previous heating rate calculations in the treatment of backscattered electrons to produce better agreement with observed flux spectra. This paper shows that more realistic model auroral electron spectra can be obtained by reflecting the up going flux back to the ionosphere at the upper boundary of the model. In this case, the neutral gas heating rates are 20%–25% higher than when the backscattered flux escapes from the ionosphere.     

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.     

Numerical modeling of propagation of breaking nonlinear acoustic-gravity waves from the lower to the upper atmosphere

Acoustic-gravity waves (AGWs) observed in the upper atmosphere may be generated near the Earth’s surface due to a variety of meteorological sources. Two-dimensional simulations of vertical propagation and breaking of nonlinear AGWs in the atmosphere are performed. Forcing near the Earth’s surface is used as the AGW source in the model. We use a numerical method based on finite-difference analogues of fundamental conservation laws for solving atmospheric hydrodynamic equations. This approach selects physically correct generalized solutions of the wave hydrodynamic equations. Numerical simulations are performed in a representative region of the Earth’s atmosphere up to altitude 500 km. Vertical profiles of temperature, density, molecular viscosity and heat conductivity were taken from the standard atmosphere model MSIS-90 for January. Calculations were made for different amplitudes and frequencies of lower boundary wave forcing. It is shown that after activating the tropospheric wave forcing, the initial pulse of AGWs may very quickly propagate to altitudes of 100 km and above and relatively slowly dissipate due to molecular viscosity and heat conduction. This may increase the role of transient nonstationary waves in effective energy transport and variations of atmospheric parameters and gas admixtures in a broad altitude range.     

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

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

由电离层电位所建立的大气电状态

本文在已知电离层电位分布下, 解析地计算了大气电位, 电场和电流强度的全球分布.结果表明, 在大气导电率随高度呈指数增加的情况下, 100km高度上的电离层电位, 几乎无衰减地扫到25km以下.大气电场较强的区域主要在20km以下的低层大气区, 其垂直分量比水平分量大4个数量级.而中高层大气电场较弱, 且两分量量级相当.本文还提出了一种考虑地面形状对大气电场影响的解析方法.      

雷暴云准静电场和夜间低电离层的电离

用点电荷模型计算雷暴云突然放电后形成的准静电后形成的准静电场随高度的分布，以E/N（E的电场大小，N为大气密度）为输入参量，在一定条件下，对Boltzmann方程数值求解，计算电离层电子数密度的扰动。计算结果表明，在约70-90km之间，在约放电后的10ms内，准静电场大于中性大气的击穿电场，将引起大气的雪崩电离，从而引起夜间低电离层电子密度的显著增加，但这种电子密度的增加是暂的，在很短的时期内就恢复到平静时的水平，恢复时间随高度的变化而不同。     

Thermospheric storm effects

A brief review is given of our current understanding of the atmospheric perturbations in the thermosphere and exosphere that are related to geomagnetic disturbances and of current efforts to represent these in empirical models of the upper atmosphere. A particular model, based on ESRO4 mass spectrometer observations of neutral composition and density, is presented in detail. This model gives the effects on the principal constituents of the upper atmosphere as a function of the geomagnetic coordinates and the K_p geomagnetic index. It is a modification of an earlier model, the most important difference being the inclusion of the variation with magnetic local time.     

The optical parameters of the atmosphere of Venus

This work is devoted to the derivation of the optical properties of the Venus atmosphere from “Venera-10” optical measurements. Within the framework of a two-layer model of Venus atmosphere it is found that in the spectral interval 0.52 – 0.85 μm the optical thickness of the upper cloud layer is ≈ 50 and the optical parameters of the lower layer are similar to the Rayleigh ones. Comparison is made between the measurements of radiation field within the atmosphere and the results of strict calculations. A preliminary conclusion is suggested that there are considerable numbers of aerosol particles with a radius ? 0.03 μm in the lower layer. The results of the upper boundary of the cloud layer is estimated to be ≈ 70 km.     

Mesospheric nitric oxide and ozone measurements in polar winter at 69°N

Two rocket experiments were carried out just before and after the polar night at Andoya (69°N), Norway to investigate transport of nitric oxide produced by auroral processes into the middle atmosphere and its influence on the ozone chemistry. Nitric oxide densities of (2–5) × 10⁸cm⁻³ found in the 70–90 km region are one to two orders of magnitude larger than those at middle latitudes. The influence on ozone densities in the 70–90 km region due to such enhanced nitric oxide abundance is found to be insignificant as compared to that due to transport in the middle of February. The larger ozone densities found in February (in spite of longer sunlit duration) than in November in the 40–60 km region again support predominance of transport over photochemical loss.     

Numerical modeling of inhomogeneous orographic wave influence on planetary waves in the middle atmosphere

Parameterization of dynamical and thermal effects of stationary orographic gravity waves (OGWs) generated by the Earth’s surface topography is incorporated into a numerical model of general circulation of the middle and upper atmosphere. Responses of atmospheric general circulation and characteristics of planetary waves at altitudes from the troposphere up to the thermosphere to the effects of OGWs propagating from the earth surface are studied. Changes in atmospheric circulation and amplitudes of planetary waves due to variations of OGW generation and propagation in different seasons are considered. It is shown that during solstices the main OGW dynamical and heat effects occur in the middle atmosphere of winter hemispheres, where changes in planetary wave amplitudes due to OGWs may reach up to 50%. During equinoxes OGW effects are distributed more homogeneously between northern and southern hemispheres.