北京地区大气温度及重力波活动的季节变化

利用瑞利激光雷达观测数据,分析了北京地区35~70km高度范围内大气温度和重力波活动的季节变化.发现北京地区30~70km高度范围内的大气温度有明显的年周期变化:平流层顶最高温度出现在6,7月份,大约为270K;中间层70km高度最低温度也出现在6,7月份,大约为200K.以2014年10月14日晚数据为例,分析重力波势能密度,发现50km以下重力波势能存在耗散,而在50km以上重力波近乎无耗散地向上传播.通过对比35~50km高度范围内的平均势能密度,对北京地区重力波活动强弱的季节变化进行了研究.研究结果表明,北京上空重力波活动强度具有明显的年周期变化,冬季平均势能密度为18J·kg-1,夏季为8J·kg-1,且冬季重力波活动强度约为夏季的两倍.此外,还分析了春夏秋冬四个季节重力波势能密度随高度的变化.结果表明,不同季节和不同高度的重力波势能密度不同.      

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

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

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

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

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

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

550～600 km高度上热层大气密度的中长期变化

选用了2005年8月20日至2006年7月28日高度550～600 km附近的热层大气密度探测数据,对表征太阳活动的F10.7值和表征地磁活动强度的Ap指数进行了相关特性的统计.分析结果表明,在无明显地磁扰动时热层大气密度日平均值的涨落呈现27日和准半年的周期性变化,但在地磁扰动期间这种变化的周期性会被削弱,且大气密度的周日变化幅度与F10.7值呈正相关关系.      

"神舟3号"运行高度上大气密度的变化

"神舟3号"(SZ-3)大气密度探测器搭载在SZ-3留轨舱上于2002年3月发射入轨,在轨运行期间获得了轨道舱运行高度范围(330-410km)内的大气密度数据.数据分析表明,无明显太阳和地磁扰动时,热层大气密度的主要变化之一是日照和阴影区域之间的涨落变化,最大涨落变化比约为3.0,变化比与太阳和地磁活动程度有关.在2002-04-17和2002-04-19的强地磁扰动时,全球热层大气密度上涨,同时在磁扰峰期探测获得30°N-40°N区域出现密度扰动异常现象.对强地磁扰动在运行轨道高度上大气密度最大涨幅约为60%左右,响应过程在时间上要比地磁扰动过程滞后6-7h,日照和阴影区域中大气密度的响应变化程度明显不同.在太阳活动程度发生变化时,热层大气密度会呈现出明显的正相关变化关系.     

平流层爆发性增温期间大气变化特性的COSMIC掩星观测

利用高精度和高垂直分辨率的COSMIC掩星观测资料, 详细深入分析了2007年冬---2008年春平流层爆发性增温(SSW)期间10~60 km高度范围内大气的变化特性, 尤其是上平流层和低中间层大气的变化特性. 结果表明, 在SSW过程中, 温度场、风场和剩余环流都发生了明显的变化. 根据温度在主增温前和主增温盛期的变化特性, 在水平方向, 大约以55oN为界, 在垂直方向, 大约以42 km为界, 可以将温度场在纬度-高度的分布分为4个区域: 高纬下层增温区, 增温幅度约高达25 K; 高纬上层降温区, 降温幅度约达30\,K; 中纬下层降温区, 降温幅度约为几K; 中纬上层增温区, 增温也约为几K. SSW期间上下层大气纬向风场的变化规律基本相同. 在纬度方向以45oN为界, 45oN以北地区的西风减弱东风增强, 风场变化高达50 m/s; 45oN以南地区西风增强东风减弱, 变化幅度比较小, 约10 m/s. 在2008年1月下旬到2月底, 大气温度和纬向风有明显的振荡现象, 周期约为12天. 剩余环流的环流圈在SSW期间会发生反转, 由此也表明, SSW期间大气中物质的输运方向也会发生改变.      

TIMED/SABER与AURA/MLS临近空间探测温度数据比较

利用AURA/MLS数据（V4.2）和TIMED/SABER数据（V2.0）对20~92km高度的大气温度进行比较分析,计算AURA/MLS数据与TIMED/SABER数据的温度绝对偏差,并分析平均温度偏差在不同季节中随经度、纬度和高度的变化特征.结果表明:20~80km高度的平均温度偏差在±6K以内,相对偏差在3%以内;80~90km高度平均温度偏差减小至-10K以下,相对偏差在9%以内.中低纬度地区平均温度偏差廓线的变化趋势一致,从20km高度的-3K左右的负偏差逐渐增加,在45~50km高度的平流层顶处有较明显的3K左右的正偏差峰值.平均温度偏差随纬度变化明显,随经度变化很小.研究结果可为卫星数据的应用提供参考依据.      

大气密度变化对卫星返回轨道参数的影响

本文在考虑大气密度(90km高度以下) 随地理纬度、高度及时间(月份)随机变化的基础上,沿卫星的返回轨道建立了随机大气密度的统计模型。并应用该模型产生的随机大气密度样本,进行返回轨道的Monte Corlo法模拟计算。通过对模拟计算得到的返回轨道参数样本的统计处理,本文分析了大气密度变化对卫星返回轨道参数的影响。另外,本文还介绍了近似估计大气密度变化对返回轨道参数影响的影响系数法。     

大功率无线电波与低电离层的相互作用

地面入射的大功率无线电波能加热电离层等离子体，引起电离层电子温度和密度的扰动，实现电离层的地面人工变态．本文中，着重考虑电波和电离层相互作用过程中的自吸收，构造一个自治的相互作用模型，在一定功率和频率的加热电波作用下，利用该模型计算了白天低电离层电子温度和由温度的变化而引起的电子密度的变化．计算结果表明，在白天低电离层，电波的自吸收在90km以下比较显著，而最大温度变化在70km高度上，大约增加了2倍．在α复合的假设下，电子密度变化幅度随高度的增加而减少，在70km处，大约增加55％、120km处则为4％左右．     

Periodic variations of ionization rate in the low and middle atmosphere

The cosmic ray ionization source functions which were obtained using a simplified extensive air shower model are used to calculate the eleven year cycle, seasonal and diurnal variations of ionization rate in the low and middle atmosphere. The ionization source function, as a function of the penetrating depth and the energy of cosmic ray particles, is the ionization rate per unit depth for a unit flux of incoming cosmic ray particles with certain energy.The calculation of the eleven year cycle variation of ionization rate in the low and middle atmosphere due to the modulation of galactic cosmic ray intensity by solar activity shows that the amplitude is larger at a higher magnetic latitude and is generally larger at higher altitudes. The relative amplitude of fluctuation of the ionization peak value (at altitudes near 15 km) is up to 45% in the magnetic polar region. The ionization rate, due to the seasonal variation of the atmospheric density, varies from several per cent below the ionization peak to several tens per cent above the peak. This seasonal variation of ionization rate reaches 35% at 70 km. The diurnal variation of atmospheric densities caused by atmospheric tidal oscillation can produce a diurnal variation of the ionization rate to an amplitude of several per cent at altitudes above 40 km. The diurnal oscillation is less than 1% below 35 km.     

Simulating planetary wave propagation to the upper atmosphere during stratospheric warming events at different mountain wave scenarios

Parameterization schemes of atmospheric normal modes (NMs) and orographic gravity waves (OGWs) have been implemented into the mechanistic Middle and Upper Atmosphere Model (MUAM) simulating atmospheric general circulation. Based on the 12-members ensemble of runs with the MUAM, a composite of the stratospheric warming (SW) has been constructed using the UK Met Office data as the lower boundary conditions. The simulation results show that OGW amplitudes increase at altitudes above 30 km in the Northern Hemisphere after the SW event. At altitudes of about 50 km, OGWs have largest amplitudes over North American and European mountain systems before and during the composite SW, and over Himalayas after the SW. Simulations demonstrate substantial (up to 50–70%) variations of amplitudes of stationary planetary waves (PWs) during and after the SW in the mesosphere-lower thermosphere of the Northern Hemisphere. Westward travelling NMs have amplitude maxima not only in the Northern, but also in the Southern Hemisphere, where these modes have waveguides in the middle and upper atmosphere. Simulated variations of PW and NM amplitudes correspond to changes in the mean zonal wind, EP-fluxes and wave refractive index at different phases of the composite SW events. Inclusion of the parameterization of OGW effects leads to decreases in amplitudes (up to 15%) of almost all SPWs before and after the SW event and their increase (up to 40–60%) after the SW in the stratosphere and mesosphere at middle and high northern latitudes. It is suggested that observed changes in NM amplitudes in the Southern Hemisphere during SW could be caused by divergence of increased southward EP-flux. This EP-flux increases due to OGW drag before SW and extends into the Southern Hemisphere.     

Circulation of the atmosphere from the surface to 100 km

A reference model of the atmospheric circulation on Venus based on available observations is presented. The reference atmosphere has the following main features: (i) the entire atmosphere below 85 km moves predominantly from east to west in the planet's reference system (in the same direction as the rotation of the solid planet itself) with the possible exception of the lowest 10 km where velocities are low, (ii) a jet is present near the cloud-top level at 45° latitude in both hemispheres, with a magnitude of approximately 100 ms^?1, (iii) a weak meridional (north-south component) flow directed towards either pole is superimposed on the zonal (east-west) motion at cloud-top level (about 68 km) altitudes, and, (iv) eddies or wave motions are present in the atmosphere, with amplitudes of less than 15 ms^?1 in the upper atmosphere.     

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.     

Ion production and ionization effect in the atmosphere during the Bastille day GLE 59 due to high energy SEPs

The influence of high energy particles, specifically cosmic rays, on atmospheric physics and chemistry is highly discussed. In most of the proposed models the role of ionization in the atmosphere due to cosmic rays is not negligible. Moreover, effect(s) on minor constituents and aerosols are recently observed, specifically over the polar regions during strong solar particle events. According to the recent findings for such effects it is necessary an essential increase of ion production, specifically during the winter period. The galactic cosmic rays are the main source of ionization in the Earth’s stratosphere and troposphere. Occasionally, the atmospheric ionization is significantly enhanced during strong solar energetic particles events, specifically over the polar caps. During the solar cycle 23 several strong ground level enhancements were observed. One of the strongest was the Bastille day event occurred on 14 July 2000. Using a full Monte Carlo 3-D model, we compute the atmospheric ionization, considering explicitly the contribution of cosmic rays with galactic and solar origin, focusing on high energy particles. The model is based on atmospheric cascade simulation with the PLANETOCOSMICS code. The ion production rate is computed as a function of the altitude above the sea level. The ion production rate is computed on a step ranging from 10 to 30?min throughout the event, considering explicitly the spectral and angular characteristics of the high energy part of solar protons as well as their time evolution. The corresponding event averaged ionization effect relative to the average due to galactic cosmic rays is computed in lower stratosphere and upper troposphere at various altitudes, namely 20?km, 15?km, 12?km and 8?km above the sea level in a sub-polar and polar regions. The 24^h and the weekly ionization effects are also computed in the troposphere and low stratosphere. Several applications are discussed.     

Atmospheric dynamics on Venus and Mars

New results from Pioneer Orbiter observations indicate a continued vortex organization of the cloud level atmosphere in either hemisphere, centered over respective poles. Significant changes in the magnitude of the cloud level zonal circulation over a period of several years have been detected. A strong signature of the solar tidal circulation has been detected in the atmospheric circulation with the lowest speeds occurring in equatorial latitudes about 20° upstream of the sub-solar point. Finally, a solar-locked persistent spatial structure has been discovered in the variance of the ultraviolet brightness measured from brightness normalized images of Venus. Vega balloons (drifting at about 53 km altitude near 7°N and 7°S latitudes) have also provided some unique observations of atmospheric circulation, significant among them being the strong vertical motions, the zonality of their drift speeds as well as a significant temperature difference between the two balloons. The temperature difference which amounts to 6.5°K on average is currently being interpreted as a temperature variation with longitude or time.

Diagnostic modelling efforts towards simulating the atmospheric circulation on Venus are continuing and have provided some clues about the processes that maintain them but have not yet been successful in explaining the superrotation of the atmosphere.

Knowledge of the Martian atmospheric dynamics on the other hand is still limited by lack of adequate observations. Numerical modelling of the Martian atmosphere continues to provide most of the information about the atmospheric circulation. The situation regarding the paucity of observations should improve with the completion of the proposed Mars Observer mission. The low circular polar orbit planned provides an excellent opportunity to study the Martian atmosphere.     

日全食期间中间层大气臭氧光化过程试验

利用包括活O_x、HO_x、和NO_x等成份在内的一维时变光化模式,研究了日全食时上中层大气臭氧和O₂(1△g)的变化特性,并通过与观测资料的比较进行了光化模式的试验.结果表明,日全食时太阳辐照变化对1.27μm气晖影响的高度比曙暮期间明显偏低,利用1.27μm气晖光度计测量大气臭氧的高度范围约为52-67km.1997年3月9日漠河日全食时测量的大气臭氧含量变化幅度为50％.这与考虑O_x、HO_x、和NO_x等成份的光化模式理论模拟结果基本一致.它表明在52-67km高度范围不存在臭氧亏损现象.      

太赫兹大气临边探测仪遥感中高层大气风仿真

太赫兹大气临边探测仪（TALIS）是中国正在预研的第一台THz频段的临边探测仪，主要用于高精度、高分辨率的大气遥感测量.TALIS的观测目标主要包括大气温度、大气压强、大气成分（例如H₂O，O₃，HCl，ClO，N₂O，HNO₃等）的垂直分布以及长期变化趋势.由于TALIS的频段覆盖了许多重要的吸收谱线，其观测数据中包含大气风的多普勒信息，因此可以用于反演中高层风的廓线.本文针对TALIS视线多普勒风的观测进行仿真，利用辐射传输模型（ARTS）评估了TALIS测风的潜力和相应的反演精度.结果表明，TALIS的118GHz谱仪具有较好的测量精度，在70km处的精度为12m·s-1.183GHz，633GHz和658GHz谱仪也有一定的测量信息，反演精度分别为19m·s-1（60km），19m·s-1（50km），16m·s-1（50km）.TALIS有一个候选的测风谱仪位于655GHz频段，其在55km处的反演精度为11m·s-1.此外，虽然降低谱分辨率能有效提高系统灵敏度，但并不能提高反演精度，需要通过降低系统噪声来提高测风的精度.