The relations between the linear and bilinear atmospheric tidal fields in geometric height and in log-pressure coordinates

The relationships between the linearized meteorological variables as expressed in geometric height and in log-pressure coordinates are derived from the assumptions of classical atmospheric tidal theory. While the horizontal velocity components are the same to first-order in the two coordinate systems, a linearized vertical velocity differencew-H ₀ occurs because of the periodic vertical displacement of the constant pressure surfaces due to time-dependent, hydrostatic density perturbations; a linearized temperature differenceT- also results when these displacements occur in the presence of a zero-order vertical gradient of temperature. Both of these differences can be expressed in terms of the tidal geopotential field. For a given tidal mode, both differences are generally proportional to the square root of the ratio of the tidal mode's equivalent depth and the atmospheric scale height; the temperature difference is also proportional to the background temperature lapse rate. It is further shown that the two classical tidal vertical structure equations commonly derived in their respective geometric height and log-pressure coordinate systems are in fact identical to first-orderas long as the same thermotidal forcing function is used. Expressions for the zonal-mean components of the tidal bilinear fluxes, formed by zonally averaging the product of two longitudinally varying, linearized tidal fields, are also derived for the two coordinate systems. For the bilinear fields the largest relative differences (a few tens of percent) are for the tidal zonal-mean forcing per unit mass of the zonal wind. For Earth and Mars, differences between the tidal vertical velocity fields are generally less than 25% but may be significantly larger in the Martian atmosphere during one of its episodic planetary-scale dust storms. Tidal temperature differences are generally smaller.     

The form of Newtonian cooling in atmospheric tidal theory

The different forms of the classical atmospheric tidal vertical structure equation derived by Lindzen and McKenzie (1967) and by Dickinson and Geller (1968) are shown to be due solely to their different treatment of the infrared radiative damping of tidal temperature perturbations. While both parameterize the infrared cooling as a Newtonian cooling, the former use the tidal temperature perturbation along a constant geometric height surface while the latter use the perturbation along a constant pressure surface. These two tidal temperature perturbations differ when the background atmospheric temperature varies with height. Scaling arguments show that for most applications the differences between the use of these two forms of Newtonian cooling are negligible. However, relative differences of 30% in amplitude and 15° in phase can occur in the Martian atmosphere, with its relatively short radiative time constants, for a possibly resonant planetary-scale Rossby tidal mode.     

非绝热加热对大气局地扰动位能的影响和机理

本文运用统计相关和奇异值分解方法,系统地考察了非绝热加热对大气局地扰动位能的影响特征和机理问题.分析结果表明,热带地区海表温度异常和扰动位能的耦合相关特征与厄尔尼诺和南方涛动变率的关系密切,扰动位能在热带外地区的耦合模态空间型呈现出与北太平洋—北美大气遥相关型极为相似的分布特征.通过进一步考察大气中视热源与扰动位能的耦合特征,发现在低纬地区仍主要反映了厄尔尼诺和南方涛动的影响机制,在热带外地区,视热源异常的耦合模态则表现出与北半球环状模类似的特点,这表明热带外地区大气的局地能量有效性与热带非绝热加热的遥强迫以及局地性热源加热强迫有着极为密切的联系.     

A coupled oscillator model of shelf and ocean tides

The resonances of tides in the coupled open ocean and shelf are modeled by a mechanical analogue consisting of a damped driven larger mass and spring (the open-ocean) connected to a damped smaller mass and spring (the shelf). When both masses are near resonance, the addition of even a very small mass can significantly affect the oscillations of the larger mass. The influence of the shelf is largest if the shelf is resonant with weak friction. In particular, an increase of friction on a near-resonant shelf can, perhaps surprisingly, lead to an increase in ocean tides. On the other hand, a shelf with large friction has little effect on ocean tides. Comparison of the model predictions with results from numerical models of tides during the ice ages, when lower sea levels led to a much reduced areal extent of shelves, suggests that the predicted larger tidal dissipation then is related to the ocean basins being close to resonance. New numerical simulations with a forward global tide model are used to test expectations from the mechanical analogue. Setting friction to unrealistically large values in Hudson Strait yields larger North Atlantic _M2$M_2$ amplitudes, very similar to those seen in a simulation with the Hudson Strait blocked off. Thus, as anticipated, a shelf with very large friction is nearly equivalent in its effect on the open ocean to the removal of the shelf altogether. Setting friction in shallow waters throughout the globe to unrealistically large values yields even larger open ocean tidal amplitudes, similar to those found in simulations of ice-age tides. It thus appears that larger modeled tides during the ice ages can be a consequence of enhanced friction in shallower water on the shelf in glacial times as well as a reduced shelf area then. Single oscillator and coupled oscillator models for global tides show that the maximum extractable power for human use is a fraction of the present dissipation rate, which is itself a fraction of global human power consumption.     

东经120°E中间层和低热层大气潮汐及其季节变化特征

利用为期一年的卫星遥感温度（SABER/TIMED）资料重建了120°E子午圈内中间层和低热层大气潮汐各主要频率分量（周日、半日和8小时潮汐）.这些主要频率分量随高度振幅增大，在97 km高度达到显著的振幅；其中迁移性周日潮汐在97 km高度出现极大振幅，然后随高度衰减.本文从考察迁移性成分和非迁移性成分各自在总潮汐中贡献角度出发，着重讨论了那些对形成该子午圈中97 km高度上整体潮汐扰动起控制作用的潮汐成分.结果显示，对周日和半日频率这两种潮汐而言，迁移性成分控制了它们的总体时空分布.在春分季节，迁移性周日潮的控制作用最显著，决定了赤道和两半球热带的活动中心；其中北半球副热带地区的季节变化形势与以往利用武汉（30°N，114°E）流星雷达风测量资料开展分析得到的结果是一致的；其他季节受非迁移性成分明显影响，例如，在本文关注的2005年中，夏至季节受（1，0）模、（1，-3）模和（1，-2）模的共同影响形成了从赤道向南延伸的活动中心，极值中心位于赤道附近，振幅达到了20 K以上，是全年的最大值.受迁移性成分控制，半日潮活动主要出现在两半球热带地区，北半球活动中心位于秋分季节（振幅达到13 K），南半球活动中心位于春分和夏至之间.其他季节受非迁移成分的影响，形成若干分布在两半球的活动中心.在本文关注的40°S～40°N范围内，与周日潮和半日潮相比，8小时潮汐具有显著较低的振幅；另外，虽然迁移性成分在一年中的大部分时间系统地分布在两半球热带地区，但是非迁移成分具有与迁移性成分相当或更大的振幅，在整体上控制了这种潮汐的时空分布.     

电离层人工调制激发的VLF波在磁层的传播特性及应用研究

电离层人工调制可以激发甚低频（VLF）波,其中向上传播进入磁层的VLF波,不但能够用来研究磁层中的各种物理现象,且具有人工沉降高能粒子,消除辐射带等实际用途.本文使用射线追踪方法,模拟电离层调制激发的VLF波在磁层的传播路径,分析激发纬度和调制频率对传播路径和传播特性的影响;并基于低频波的色散方程和波粒共振条件,分析VLF波传播路径上与磁层高能粒子的最低共振能及其分布.研究表明,VLF波通过在磁层来回反射向更高的L-shell传播,最终稳定在某一L-shell附近.以较低的调制频率或者从较高的纬度激发的VLF波能够传播到更高的L-shell,但是,当激发纬度过高时,低频波也可能不发生磁层反射而直接进入电离层和大气层.低频波在磁层的传播过程中,在较高的纬度或者较低的L-shell能够与较高能量的电子发生共振相互作用,在较高的L-shell并且低纬地区,能够与较低能量的电子发生共振相互作用.共振谐数越高,能发生波粒共振的电子能量越高.     

机载激光雷达对青岛及周边海域的气溶胶探测

工作波长为532 nm的机载大气环境探测激光雷达AEDAL（Atmospheric Environment Detecting Airborne Lidar）装载在CMS_3807飞机上,于2005年11月7～11日期间在青岛地区及周边海域上空进行了飞行探测.此次实验的目的有两个：验证我国用于大气环境探测的激光雷达技术已经具备从地基向空基乃至天基发展的条件;获得青岛地区及周边海域边界层结构及大气气溶胶时空分布变化的特点.激光雷达的高时空分辨率为获取飞行路径上的边界层结构及气溶胶时空分布提供了可能.为了研究下垫面对边界层及气溶胶时空分布的影响,预定的飞行路径上包含了丰富的地形变化,有城市、丘陵、海区等.通过给出11月8日及11日的探测结果,不仅得到了不同地区边界层结构及气溶胶的时空分布特点,还可以看到冷锋、地形、地面气象场等因素对它们的影响.     

强潮汐激发地震火山活动的新证据

强潮汐与火山地震活动密切相关. 天文资料表明, 2000年6月到8月日本Izu半岛最活跃的火山地震活动正好处于天文大潮时期和日长变化最小值时期. 在厄尔尼诺事件和拉尼娜事件发生前后,东西太平洋海面高度分别升降40cm, 水均衡作用使洋壳反向升降13cm. 由此形成的东西太平洋地壳跷跷板运动加强了强潮汐对地震火山活动的激发作用.     

三维分层流动过双山脉地形激发的大气船舶重力波动力学理论和数值试验

利用我们建立的三维分层线性理论计算模式和中尺度数值模式ARPS, 分别研究了三维分层流动过双山脉地形产生的三维线性和非线性山脉重力波和大气船舶的结构特征及其形成机制.线性理论计算结果表明三维三层流动过双山脉地形时,两个山脉各自强迫出一个发散模态的山脉背风波,在第二个山脉背风面,三维三层流动过双山脉地形可以强迫出两个发散模态的拦截背风波,大大加强了对大气环流的拖曳作用.非线性数值模拟结果表明,流动过山所产生的非线性山脉重力波和大气船舶完全不同于三维分层线性理论计算模式所产生的山脉重力波和大气船舶的结构和特征,由于分层流体之间的非线性相互作用,三维三层流动过双山脉地形时,可在第二个山脉背风面激发4个发散模态的拦截背风波. 三维三层流动过双山脉地形所强迫的山脉重力波和大气船舶,具有同三维三层流动过孤立山脉所产生的山脉重力波和大气船舶完全不同的结构和特征,三维流动过双山脉地形对两个山脉之间的距离表现出极大的敏感性.对于相距较远的两个山脉,流动过双山脉所强迫的山脉重力波表现为4个发散模态的拦截背风波,波动的能量相对于相距较近的两个山脉能传播到更高的高度.     

A study on three-dimensional structures of the ionospheric dynamo currents induced by the neutral winds simulated by the Kyushu-GCM

Three-dimensional structures of the ionospheric dynamo currents are examined using the neutral winds in a general circulation model of the middle atmosphere at Kyushu University. A quasi-three-dimensional ionospheric dynamo model is constructed assuming an infinite parallel conductivity in the ionosphere. This model is able to simulate both the equatorial electrojet and the global Sq current system successfully. The simulated results reveal that the equatorial electrojet is confined in quite narrow latitudes around the equator accompanied with meridional current circulations and satisfies a non-divergent structure mainly within the E region. A vertically stratified double layered structure is seen in the east–west current density near the focus latitude of the global Sq current system. It is shown that the stratified structure mainly consists of the east–west Hall current associated with the eastward wind of zonal wavenumbers 1 and 2 in the lower altitudes and the westward wind of zonal wavenumber 2 in the upper altitudes. The day-to-day variation of the neutral winds can significantly vary the induced ionospheric dynamo current system, which is recognized as changes of the focus latitude and/or the maximum value of the equatorial electrojet.     

Investigation of the time variability of diurnal tides and resonant FCN period

The time variability of diurnal tides was investigated by analyzing gravity observations from global superconducting gravimeter (SG) stations with running time intervals. Through least-square and Bayesian approaches, FCN resonance parameters were estimated for each data section after obtaining the tidal parameters of mainly diurnal tidal waves. The correlation of the time variation in diurnal tidal waves and FCN period was discussed. For comparison, a similar method was used to analyze VLBI observations to study the time variability of nutation terms and FCN period. The variation trend of the FCN period totally depends on the Ψ1 wave in tidal gravity and on the retrograde annual term in nutation. We observed a similar variation trend in the FCN periods obtained from different SG stations worldwide and VLBI observations. The relation between diurnal tides and LOD variations is discussed and the possible mechanisms of the decadal variation in FCN periods were discussed.     

A note on the maintenance of the atmospheric kinetic energy

The winter simulations of the GLAS climate model and the NCAR community climate model are used to examine the maintenance of the atmospheric kinetic energy. It is found that the kinetic energy is generated in the lower latitudes south of the maximum westerlies, transported northward and then, destroyed in the midlatitudes north of the maximum westerlies. Therefore, the atmospheric kinetic energy is maintained by the counterbalance between the divergence (convergence) of kinetic energy flux and generation (destruction) of kinetic energy in lower (middle) latitudes.     

非迁移潮作用下的电离层赤道异常区单峰现象及其经度分布

当喷泉效应较弱而双峰结构发展不充分的时候,可能在赤道异常区仅能够观测到一个电子密度的峰值,称之为单峰现象.本文利用CHAllenging Minisatellite Payload卫星在2001-2010年的电子密度数据给出了单峰的发生规律,单峰在地方时早上08：00-10：00和下午16：00-19：00发生率高,发生位置在经度上呈现多波数分布,尤其在10：00-18：00明显：在分季时多呈现四波,而在冬至季时以三波为主.单峰发生多的经度,正好对应着双峰的结构特征较弱之处.究其原因,是非迁移潮的DE2和DE3分量调制了背景风场和大气发电机电场,在电场和喷泉效应减弱的经度,双峰结构难以形成时,就会表现为单峰结构.本文扩展了对单峰现象的地方时、季节和经度分布等规律的了解,明确了非迁移潮在其中施加的影响,由此,单峰同双峰现象一样可以用于研究非迁移潮对热层-电离层的作用.

     

大气变率对北极地区近期海冰变化趋势影响数值模拟研究

为研究近期21年(1989—2009年)北极地区海冰变化原因,本文利用欧洲中期天气预报中心ERA-Interim数据集资料和美国麻省理工学院MITgcm全球海冰-海洋耦合模式开展了不同大气强迫条件下海冰变化的数值模拟研究.研究工作中共设计了6个数值试验,除1个试验全部采用1989—2009年每日4个时次的大气强迫场外,其余5个试验各有一种大气强迫(地表气温、地表大气比湿、向下短波辐射通量、向下长波辐射通量和地表风)采用1989年月平均结果.分析了各模拟试验结果中3月和9月北极地区海冰面积的年际变化特征及最小二乘拟合意义下的线性变化趋势,并以ERA-Interim结果为参照标准对各模拟试验结果进行了对比和检验,以说明不同大气强迫量变率对海冰变化的作用.结果表明:地表气温变率和向下长波辐射通量变率是造成海冰面积减少的主要原因;向下短波辐射通量变率对海冰面积变化影响几乎可以忽略;地表大气比湿变率对海冰面积线性变化趋势影响较小,但对海冰面积年际变化特征有调制作用;地表风变率对海冰季节变化、海冰面积线性变化趋势及年际变化特征均有明显影响,说明提高大气风应力精度是改善海冰数值模拟结果的重要手段.     

大气变率对北极地区近期海冰变化趋势影响数值模拟研究

为研究近期21年(1989—2009年)北极地区海冰变化原因,本文利用欧洲中期天气预报中心ERA-Interim数据集资料和美国麻省理工学院MITgcm全球海冰-海洋耦合模式开展了不同大气强迫条件下海冰变化的数值模拟研究.研究工作中共设计了6个数值试验,除1个试验全部采用1989—2009年每日4个时次的大气强迫场外,其余5个试验各有一种大气强迫(地表气温、地表大气比湿、向下短波辐射通量、向下长波辐射通量和地表风)采用1989年月平均结果.分析了各模拟试验结果中3月和9月北极地区海冰面积的年际变化特征及最小二乘拟合意义下的线性变化趋势,并以ERA-Interim结果为参照标准对各模拟试验结果进行了对比和检验,以说明不同大气强迫量变率对海冰变化的作用.结果表明:地表气温变率和向下长波辐射通量变率是造成海冰面积减少的主要原因;向下短波辐射通量变率对海冰面积变化影响几乎可以忽略;地表大气比湿变率对海冰面积线性变化趋势影响较小,但对海冰面积年际变化特征有调制作用;地表风变率对海冰季节变化、海冰面积线性变化趋势及年际变化特征均有明显影响,说明提高大气风应力精度是改善海冰数值模拟结果的重要手段.     

On the contributions of atmospheric pressure and wind to daily sea level in the northern Adriatic Sea

Daily sea level variability in the Adriatic Sea is studied from different data sets using Empirical Orthogonal Functions, in connection with atmospheric pressure and wind stress. The first mode explains 56–69% of total variance and consists of uniform sea level variability all over the basin, correlated with atmospheric pressure through the inverse barometer effect. The second mode explains 13–16% of variance and accounts for an along-basin sea level gradient, which is correlated with the meridional wind stress component. The first two Principal Components are used as proxies to pressure- and wind-induced components of storm surges in the northern Adriatic. The analysis of the frequency of the most intense events in the 1957–2005 period shows that the wind contribution to storm surges has decreased, while no significant trends are found in the contribution of atmospheric pressure.     

Estimation of model errors generated by atmospheric forcings for ocean data assimilation: experiments in a regional model of the Bay of Biscay

The characterization of model errors is an essential step for effective data assimilation into open-ocean and shelf-seas models. In this paper, we propose an experimental protocol to properly estimate the error statistics generated by imperfect atmospheric forcings in a regional model of the Bay of Biscay, nested in a basin-scale North Atlantic configuration. The model used is the Hybrid Coordinate Ocean Model (HYCOM), and the experimental protocol involves Monte Carlo (or ensemble) simulations. The spatial structure of the model error is analyzed using the representer technique, which allows us to anticipate the subsequent impact in data assimilation systems. The results show that the error is essentially anisotropic and inhomogeneous, affecting mainly the model layers close to the surface. Even when the forcings errors are centered around zero, a divergence is observed between the central forecast and the mean forecast of the Monte Carlo simulations as a result of nonlinearities. The 3D structure of the representers characterizes the capacity of different types of measurement (sea level, sea surface temperature, surface velocities, subsurface temperature, and salinity) to control the circulation. Finally, data assimilation experiments demonstrate the superiority of the proposed methodology for the implementation of reduced-order Kalman filters.     

Wave-mean flow interaction and its relationship with the atmospheric energy cycle with diabatic heating

In a linear system, the wave characteristics depend strongly on the distributions with height of wind ve-locity and static stability. The simplest case is for con-stant wind and static stability (e.g., isothermal atmos-phere with pressure and density exponentially de-creasing with height). In such circumstances there is no convergence or divergence in wave energy flux, therefore, no energy exchange between the wave and mean flow. In the atmosphere wind speed varies with increasing height, inte…     

海潮误差对GRACE时变重力场解算的影响研究

海潮误差是 GRACE 时变重力场反演中重要的误差源,目前发布的海潮模型中主要包含振幅较大的主潮波分量模型,在时变重力场反演中次潮波的影响也是不可忽略的,因此,GRACE 时变重力场反演中的海潮误差主要包括受限于海潮模型误差和次潮波影响.本文利用轨道模拟方法检测了短周期潮波的混频周期以及次潮波对ΔC20, ΔC30的时序特征,并进一步通过轨道模拟结果分析了海潮误差对时变重力场反演的影响,然后通过实测数据解算分析了海潮误差对当前 GRACE 时变重力场解算的影响,研究发现:(1) 利用轨道模拟能够有效地检测短周期潮波的混频周期;(2)时变重力场解算过程中,次潮波的影响大于海潮模型误差的影响;(3)海潮模型误差以及次潮波影响是当前 GRACE 没有达到基准精度的重要因素之一.     

体应变潮汐观测资料信息提取与映震效果

根据体应变观测资料的特点, 研究了体应变观测资料的处理方法, 提出了在时间域分析的数学模型。在数学模型中, 将潮汐因子、气压系数、非潮汐应变速率和加速度一并考虑, 在资料处理结果中, 这些参数可同时取得。分别用时间域分析方法和频率域分析方法对观测数据进行了处理, 发现潮汐因子、气压系数、应变速率和加速度等参数在大同地震前均有大小不同、持续时间不等的异常, 这些参数可作为物理意义明确的地震中短期预测的指标。