太阳和地磁活动对磁赤道地区低热层NO密度的影响研究

利用SNOE卫星1998年3月11日至2000年9月30日共计935天观测的NO密度和太阳软X射线数据,分析了低热层NO的时空分布特征及其对太阳和地磁活动的响应,得出了以下结论:NO密度从96.67 km开始增加,大约在105～110 km高度达到最大,随后开始减小;同一高度处一般夏季期间最大,冬季次之,春秋分季最小;密度峰值大小变化范围约为(0.5～1.5)×108 mol/cm3,峰值高度基本分布在107 km和113 km高度处,且不随太阳活动变化,平均值约为107 km;NO密度与太阳软X射线及地磁Ap指数的相关系数在不同高度存在0、1和2天的最佳延迟时间,而同太阳软X射线的统计关系在不同高度和季节存在"线性"、"放大"和"饱和"现象;从统计和事件分析结果来看,太阳活动对磁赤道地区低热层NO密度的气候尺度变化的影响远大于地磁活动,但地磁活动对NO短期变化贡献非常明显.     

极点热层密度变化特征及其受太阳风扇形结构调制

本文基于2002年至2010年的GRACE卫星的观测密度统计分析南北极点的热层大气密度的世界时（即磁地方时）变化.研究发现：在9—11月份地球处于行星际磁场为背向太阳的扇区内（背向扇区）时,南极点热层密度在约17:00 UT（13:30 MLT）达到最大值,比日平均值高约22%;而在6—8月份,当地球处于行星际磁场为面向太阳的扇区内（面向扇区）时,北极点热层密度在06:00 UT（12:30 MLT）达到最大值,比日平均值高约13%.南极点的磁纬是-74°,其在15:30 UT处于磁地方时正午,恰与极尖区位置重合.北极点在5:30 UT处于磁地方时正午,此时北极点与极尖区位置最靠近.因此,极点热层大气密度的磁地方时变化可能是其周期性靠近极尖区的结果.南北极点热层密度的磁地方时变化分别在背向和面向扇区内更明显,这可能与行星际磁场B_y分量对南北半球密度的不同影响有关.统计结果还表明,极点热层大气密度的磁地方时变化在冬季半球内不明显.这可能是由于在冬季半球,沉降于极尖区的粒子相比夏季半球少、沉降高度低,因而能量沉降所引起的热层上部的密度增强较小.

     

利用IGS数据分析全球TEC的周年和半年变化特性

利用太阳活动高年（2000年）IGS提供的全球TEC数据,采用傅里叶展开的方法,分析了白天电离层TEC周年和半年变化的全球特征.结果显示：电离层TEC周年变化幅度在南北半球中高纬度地区较大、赤道和低纬地区很小.半年变化的幅度在“远极地区”（远离地球南北地磁极点的东北亚和南美地区） 比“近极地区”（靠近地球南北地磁极点的北美和澳大利亚）大得多.进一步的统计显示,全球大部分地区TEC在春秋月份出现最大值,北半球近极地区最大值在冬季出现.南半球的南美和澳大利亚部分地区,最大值出现在夏季.同样,采用傅里叶方法分析了中性大气模式MSIS90计算的全球大气原子分子浓度比值（［O/N2］）的数据,发现在南北半球中高纬度地区,中性成分［O/N2］周年变化幅度较大且有明显的冬季异常现象,依据Rishbeth等提出的理论,我们认为大气成分［O/N2］可能对TEC周年变化的产生有重要作用,并且也是TEC在近极地区出现冬季异常现象的主要原因.TEC半年变化的全球分布特征形成的原因较复杂,我们初步分析可能是由于中性成分［O/N2］、太阳天顶角控制的电离层光化学产生率变化共同作用而产生的.     

中纬热层大气质量密度经向结构差异研究

本文利用CHAMP卫星以及全球电离层-热层模型（GITM）来研究太阳活动低年（2007-2009年）中纬热层大气质量密度（ρ）的经度结构变化.结果如下：（1）ρ存在明显的经度单波结构（单峰和单谷）,且南北半球反相,波峰和波谷随着地方时增加而向东移动;（2）模拟表明离子拖曳效应在ρ结构差异的形成中起到了重要的作用,欧亚地区电子密度经度差异性较弱,不足以影响ρ经度分布,导致该地区ρ经度差异不明显;（3）在磁中纬地区,太阳天顶角的经度差异可达20°~30°,太阳光加热的经度不均匀性是导致ρ经度差异的另一个主要原因.     

磁暴期间热层大气密度变化

基于CHAMP卫星资料,分析了2002—2008年267个磁暴期间400 km高度大气密度变化对季节、地方时与区域的依赖以及时延的统计学特征,得到暴时大气密度变化的一些新特点,主要结论如下: 1)两半球大气密度绝对变化(δρa)结果在不同强度磁暴、不同地方时不同.受较强的焦耳加热和背景中性风共同作用,在北半球夏季,中等磁暴过程中夜侧和大磁暴中,夏半球的δρa强于冬半球;由于夏季半球盛行风环流造成的扰动传播速度快,北半球夏季日侧30°附近大气,北(夏)半球到达峰值的时间早于南(冬)半球.而可能受半球不对称背景磁场强度所导致的热层能量输送率影响,北半球夏季强磁暴和中磁暴个例的日侧,南半球δρa强于北半球;春秋季个例中日侧30°附近大气,北半球先于南半球1~2 h达到峰值. 2)受叠加在背景环流上的暴时经向环流影响,春秋季暴时赤道大气密度达到峰值的时间最短,日/夜侧大气分别在Dstmin后1 h和2 h达到峰值.至点附近夜侧赤道大气达到峰值时间一致,为Dstmin后3 h;不同季节日侧结果不同,在北半球冬季时赤道地区经过更长的时间达到峰值. 3)日侧赤道峰值时间距离高纬度峰值时间不受季节影响,为3 h左右.在春秋季和北半球冬季夜侧,赤道大气密度先于高纬度达到峰值,且不同纬度大气密度的峰值几乎无差别,表明此时低纬度存在其他加热源起着重要作用.     

Stratospheric ozone response to a solar proton event: Hemispheric asymmetries

Ozone depression in the polar stratosphere during the energetic solar proton event on 4 August 1972 was observed by the backscattered ultraviolet (BUV) experiment on the Nimbus 4 satellite. Distinct asymmetries in the columnar ozone content, the amount of ozone depressions and their temporal variations above 4 mb level (38 km) were observed between the two hemispheres. The ozone destroying solar particles precipitate rather symmetrically into the two polar atmospheres due to the geomagnetic dipole field These asymmetries can be therefore ascribed to the differences mainly in dynamics and partly in the solar illumination and the vertical temperature structure between the summer and the winter polar atmospheres. The polar stratosphere is less disturbed and warmer in the summer hemisphere than the winter hemisphere since the propagation of planetary wave from the troposphere is inhibited by the wind system in the upper troposphere, and the air is heated by the prolonged solar insolation. Correspondingly, the temporal variations of stratospheric ozone depletion and its recovery appear to be smooth functions of time in the (northern) summer hemisphere and the undisturbed ozone amount is slighily, less than that of its counterpart. On the other hand, the tempotal variation of the upper stratospheric ozone in the winter polar atmosphere (southern hemisphere) indicates large amplitudes and irregularities due to the disturbances produced by upward propagating waves which prevail in the polar winter atmosphere. These characteristic differences between the two polar atmospheres are also evident in the vertical distributions of temperature and wind observed by balloons and rocker soundings.     

火星电离层电子密度对太阳辐射变化的响应

太阳辐射是火星电离层变化的重要控制因素.利用火星全球勘探者号(Mars Global Surveyor,MGS)电离层掩星探测数据,并结合一个火星电离层总电子含量(Total Electron Content,TEC)经验模型,研究了火星北半球高纬地区电离层电子密度对太阳辐射变化的响应特性.在考虑了火星掩星数据中电离层电子密度对太阳天顶角的依赖关系后,发现随着太阳辐射增强,火星电离层M₂层峰值密度增大,但增长偏离线性趋势,而M₂层峰值高度和大气标高没有很明显的变化趋势.从100~200 km高度区域掩星电子密度剖面积分得到的TEC及底部和顶部TEC也随太阳辐射增大而增大,但增长率有所减小,表明火星电离层可能存在类似地球电离层的饱和特征.MGS掩星TEC及其底部和顶部剖面的TEC与经验模型TEC的比值均与太阳辐射强度变化呈反相关特性,表明在强太阳辐射情形下200 km以上电子含量在TEC中占比增大.这一特性意味着太阳活动增强,在火星顶部电离层区域,动力学过程对电离层的控制逐渐超过光化学过程.

     

Simulation of the stratospheric ozone and temperature response to the solar irradiance variability during sun rotation cycle

Despite substantial progress in atmospheric modeling, the agreement of the simulated atmospheric response to decadal scale solar variability with the solar signal in different atmospheric quantities obtained from the statistical analysis of the observations cannot be qualified as successful. An alternative way to validate the simulated solar signal is to compare the sensitivity of the model to the solar irradiance variability on shorter time scales. To study atmospheric response to the 28-day solar rotation cycle, we used the chemistry–climate model SOCOL that represents the main physical–chemical processes in the atmosphere from the ground up to the mesopause. An ensemble simulation has been carried out, which is comprised of nine 1-year long runs, driven by the spectral solar irradiance prescribed on a daily basis using UARS SUSIM measurements for the year 1992. The correlation of zonal mean hydroxyl, ozone and temperature averaged over the tropics with solar irradiance time series have been analyzed. The hydroxyl has robust correlations with solar irradiance in the upper stratosphere and mesosphere, because the hydroxyl concentration is defined mostly by the photolysis. The simulated sensitivity of the hydroxyl to the solar irradiance changes is in good agreement with previous estimations. The ozone and temperature correlations are more complicated because their behavior depends on non-linear dynamics and transport in the atmosphere. The model simulates marginally significant ozone response to the solar irradiance variability during the Sun rotation cycle, but the simulated temperature response is not robust. The physical nature of this is not clear yet. It seems likely that the temperature (and partly the ozone) daily fields possess their own internal variability, which is not stable and can differ from year to year reflecting different dynamical states of the system.     

顶部电离层离子密度经度结构的特征及其随季节、太阳活动和倾角的变化

本文利用DMSP卫星测量数据和傅里叶分解和重构方法, 研究了地磁平静期顶部电离层总离子密度(Ni)经度结构的多重波数特征及波数4的年变化、逐年变化、地方时差异和随倾角的变化.傅里叶分解和重构的结果表明, 顶部电离层平均Ni的经度结构中同时含有以波数1至波数4为主的多重分量, 不同波数分量的幅度和相位各不相同.对波数4分量的分析表明, 波数4的幅度在春秋季最强, 北半球夏季高于冬季;随太阳活动水平增强, 波数4分量的幅度增高, 至太阳活动高年幅度达到最高, 此后随太阳活动水平降低而减小, 与F10.7呈正相关;春秋季和北半球夏季波数4分量在傍晚最强, 晚上和上午次之, 黎明最弱, 从09 LT到21 LT, 波数4的相位依次滞后, 暗示向东移动.分析还发现, 日落期间波数4幅度依赖倾角, 春秋季随倾角的变化呈双峰结构, 两个极大出现在倾角±18°附近, 暗示赤道等离子体喷泉效应对顶部电离层经度结构的控制作用.     

Observations of Antarctic precipitable water vapor and its response to the solar activity based on GPS sensing

Predicting global climate change is a great challenge and must be based on a thorough understanding of how the climate system components behave. Precipitable water vapor (PWV) is one of the key components in determining and predicting the global climate system. It is well known that the local surface temperature and pressure have a direct influence on the production of PWV. However, the influence of solar activity on atmospheric dynamics and their physical mechanisms is still an open debate, where past studies are focused at mid-latitude regions. A new method of determining and quantifying the solar influence on PWV based on GPS observations to correlate the GPS PWV and total electron content (TEC) variations is proposed. Observed data from Scott Base (SBA) and McMurdo (MCM) stations from 2003 to 2005 have been used to study the response of PWV to solar activity. In the analysis, the effects of local conditions (wind speed and relative humidity) on the distribution of PWV are investigated. Results show significant correlation between PWV and solar activity for four geomagnetic storms, with correlation coefficients of 0.74, 0.77, 0.64 and 0.69, which are all significant at the 95% confidence level. There was no significant correlation between TEC and PWV changes during the absence of storms. On a monthly analysis, a strong relationship exists between PWV and TEC during storm-affected days, with correlation coefficients of 0.83 and 0.89 (99% confidence level) for SBA and MCM respectively. These indicate a statistically significant seasonal signal in the Antarctic region, which is very active (higher) during the summer and inactive (lower) for the winter periods.     

对应地震反应谱的当量功率谱密度函数的数值算法

提出根据地震反应谱曲线求解当量功率谱密度函数曲线的一种数值算法,为结构随机地震反应分析中确定和计算地面运动功率谱密度提供依据。在相对位移谱与结构相对位移响应峰值的均值之间建立联系;在现有迭代方法的基础上,对算法公式作进一步的推导和整理,将迭代计算转变为线性方程组的求解计算,进而避免迭代方法的不收敛。通过对数值计算结果的分析发现,由于规范反应谱曲线存在斜率突变,采用本文算法计算出的功率谱密度曲线在规范反应谱曲线的两端和分段转折点处会出现倍周期分岔甚至混沌现象。为此建议适当增加低频范围内离散频率点的数目,且采取整体光滑的反应谱曲线,从而提高计算结果的稳定性。     

地球磁场对太阳风的加卸载响应与地震

将磁暴过程作为地球磁场对太阳风的加卸响应，计算分析了北京地磁台１９６５．１－１９７９．１２和１９８９．１－１９９１．１２共１８年垂直分量Ｚ的暴时场Ｄｓｔ加卸载响应比值ｆＤ（Ｚ）的变化。     

The Polar Cap (PC) indices: Relations to solar wind parameters and global magnetic activity

The polar geomagnetic activity resulting from solar wind–magnetosphere interactions can be characterized the Polar Cap (PC) indices, PCN and PCS. PC index values are derived from polar magnetic variations calibrated on a statistical basis such that the index approximate values in units of mV/m of the interplanetary “geo-effective” (or “merging”) electric field (E_M) conveyed by the solar wind. The timing and amplitude relations of the PC index to solar wind plasma and magnetic field parameters are reported. The solar wind effects are parameterized in terms of the geo-effective electric field (E_M) and the dynamical pressure (P_DYN). The PC index has a delayed and damped response to E_M variations and display saturation-like effects for E_M values exceeding 10 mV/m. Steady or slowly varying levels of solar wind dynamical pressure have little or no impact on the PC index above the effects related to E_M for which the solar wind velocity is also a factor. Sharp increases in the dynamical pressure generate impulsive variations in the PC index comprising a initial negative impulse of 5–10 min duration followed by a positive impulse lasting 10–20 min. Typical amplitudes of both the negative and the positive impulses are 0.2–0.5 units. A sharp decrease in the pressure produces the inverse sequence of pulses in the PC index. Auroral substorm activity represented by the AL index level has a marked influence on the average PC/E_M level at the transition from very quiet (AL0 nT) to disturbed conditions while more or less disturbed conditions (AL<−100 nT) have no systematic effect on the average PC/E_M values. At distinct substorm events the PC/E_M ratio has a minimum (0.8) in the pre-onset phase at around 20 min before substorm onset. The average ratio gradually increases in the expansion phase to reach a maximum value (1.1) at around 40 min after substorm onset (or 20 min after the largest (negative) peak in AL). At substorm recovery during the next 2 h the PC/E_M ratio decreases. Finally, we report on the application of polar magnetic variations to model the disturbance storm time (Dst) index development during magnetic storms by using the PC index as a source function to quantify the energy input to the ring current representing accumulated storm energy and characterized by the Dst index.     

台风暴雨过程中位势散度波作用密度分析和预报应用研究

在位涡的基础上,考虑散度效应,把水平风矢量旋转90°后的旋度在广义位温梯度方向上的投影定义为位势散度,并把其二阶扰动量定义为位势散度波作用密度,该波作用密度代表扰动热量的扰动输送,与位势稳定度的发展演变有关.登陆台风Morakot(2009)的诊断分析表明,在台风登陆台湾岛,进入台湾海峡和登陆大陆三个阶段,波作用密度高值区与大陆和台湾降水区相重叠,扰动热量的扰动输送显著;在前两个阶段,大陆和台湾降水区波作用密度逐渐增大,扰动热量的扰动输送逐步增强;在台风登陆福建霞浦市后,大陆地区波作用密度减小,扰动热量的扰动输送减弱.在影响波作用密度局地变化的诸多物理因素中,扰动非地转风位涡项是主要强迫项.统计分析表明,波作用密度与观测降水有一定的相关性,据此发展了位势散度波作用密度降水预报方程,其本质是数值模式预报资料的动力延伸释用,预报应用表明,利用24 h模式预报场计算的波作用密度降水能够较好地再现观测降水的落区,对台风Morakot引发的暴雨有一定的预报技巧.     

The response of CO2 flux to rain pulses at a saline desert

As the substantial component of the ecosystem respiration, soil CO₂ flux is strongly influenced by infrequent and unpredictable precipitation in arid region. In the current study, we investigated the response of soil CO₂ flux to rain pulses at a saline desert in western China. Soil CO₂ flux was measured continuously during the whole growing season of 2009 at six sites. We found that there were remarkable changes in amplitude or diurnal patterns of soil CO₂ flux induced by rainfall events: from bimodal before rain to a single peak after that. Further analysis indicated that there is a significant linear relationship (P < 0.001) between soil CO₂ flux and soil temperature (T_soil). However, a hysteresis between the waveform of diurnal course of CO₂ flux and T_soil was observed: with soil CO₂ flux always peaked earlier than T_soil. Furthermore, a double exponential decay function was fitted to the soil CO₂ flux after rainfall, and total carbon (C) releases were estimated by numerical integration for rainfall events. The relative enhancement and total C release, in association with the rain pulses, was linearly related to the amount of precipitation. According to the size and frequency of rainfall events, the total amount of C release induced by rain pulses was computed as much as 7.88 g C·m^–2 in 2009, equivalent to 10.25% of gross primary production. These results indicated that rain pulses played a significant role in the carbon budget of this saline desert ecosystem, and the size of them was a good indicator of rain‐induced flux enhancement. Copyright © 2012 John Wiley & Sons, Ltd.     

地壳应变场对气压短周期变化的响应

大气压力变化可分为由日、月引力和太阳热辐射产生的潮汐变化和非潮汐频段的短周期气压变化,其引起的地壳应变可达10-9量级.我们选用分钟采样的观测精度达10-9量级的钻孔应变仪,对山东省泰安台2O02～2005年气压及体应变分钟值资料进行深入分析,研究了地壳应变场对短周期气压变化的响应规律.频谱分析结果表明,短周期气压变化对体应变的影响主要集中在周期大于2000 s的低频段.利用小波分解方法,得出了短周期气压系数随不同频段和不同时间的变化规律:高频段(周期小于960 s)气压系数数值较小,随频率变化幅度较大;低频段(周期大于960 s)气压系数数值较大,但变化趋于平稳;在频段一定的情况下,短周期气压系数随时间变化较平稳.从钻孔应变仪观测原理和按余弦函数分布的载荷作用于地壳表面的应力解出发,得出了地壳不同深度对不同频率短周期气压场的响应,理论计算与实际观测基本相符.研究结果表明,短周期气压变化对低频信号(周期大于2000 s)的影响较大,提取低频地球自由振荡信号时要消除气压影响,对苏门答腊大地震激发的地球球型自由振荡的功率谱分析验证了这一结论.     

沉水植物苦草根系功能性状对基质类型和种植密度的响应及其与个体表现关系

根系功能性状体现了植物细根的生长状况及其对外部环境的适应性,然而根系功能性状响应环境变化的时间稳定性如何仍不清楚。苦草(Vallisneria natans)是水体沉水植被修复的先锋种类,细根在其种群重建初期起着重要作用。该研究设计了草垫、草垫+沙子、草垫+黄泥及草垫+底泥4种生长基质和100株/m²以及200株/m²两种种植密度并将其两两组合,以模拟自然生境不同的基质和密度情况,并在不同时期对苦草功能性状指标进行取样测定,通过重复测量二元方差分析研究基质和密度对根系功能指标的影响,并通过线性拟合模型探讨个体水平的生长表现(生物量分配、营养吸收、个体定植)与具体的根系功能性状之间的关系对环境变化的响应。结果表明:所有测量指标均受到基质条件的显著影响,部分根系功能性状指标如根冠比(RSR)、根干重、比根长(SRL)、根比表面积(SRA)、根组织密度(RTD)、根系锚定力和根体积受密度的影响显著。苦草生物量的分配主要受基质条件的影响,对营养的吸收、运输及根系固着能力受基质和密度的共同影响。线性拟合模型表明SRA、RTD、根直径、根表面积以及根体积在个体水平上能一定程度地预测RSR对环境变化的响应,且RTD的预测效果最好;根冠比、SRA、生根数、总根长、根表面积以及根体积在个体水平上能一定程度地预测SRL对环境变化的响应,其中生根数的预测效果最好;SRA、根直径、根表面积以及根体积在个体水平上能一定程度地预测根系锚定力对环境变化的响应,并且SRA的预测效果最好。但是基质和密度在不同时期对苦草根系功能性状的影响方向会发生改变,当用根系功能性状预测植物个体表现时,需要考虑种群重建所处的时期。