与东北冷涡相伴的高空急流诱发平流层重力波的数值模拟研究

使用中尺度数值模式WRF-ARW,针对2010年6月发生在中国东北地区一例伴随对流层高空西风急流(位于~9km高度)演变过程出现的平流层重力波活动特征开展了数值模拟.事件发生期间,对流层区域环流处在一个东北冷涡系统的控制之下.模拟结果再现了该东北冷涡的发展和维持过程,以及与之相伴的高空急流的特征.模拟结果揭示出在急流区域上空的平流层中存在显著重力波活动现象.分析结果显示,重力波活动与急流存在紧密联系,在水平方向上,重力波呈显著的二维结构,出现在急流出口区上部并逆背景流向西传播.功率谱分析结果表明盛行波动具有~700km水平尺度、9~12h时间尺度以及4~5km垂直波长.由于急流的存在,造成其与平流层中下部之间存在显著的水平风速垂直切变,与切变相伴的耗散使得上传的重力波动量通量数值随着高度升高而递减.同时,在18~20km高度间出现的西风-东风转换带极大地抑制了波动在垂直方向的传播,形成显著动量通量沉积效应.估算结果表明,在11~20km高度之间,这种效应的整体作用相当于对该层背景流施加强度为0.86m·s-1·day-1的动力阻曳.     

中层大气中重力波饱和机制的数值分析

通过对向上传播的重力波波包在中层大气中的非线性传播过程进行数值模拟,讨论中层大气中重力波的饱和机制.数值模拟结果表明,向上传播的重力波波包的扰动振幅在接近波包的本征水平相速度之前随高度单调增长,而当波振幅接近本征水平相速度时,在对流不稳定区域出现等位温面的翻转,同时波振幅的增长达到饱和（波振幅随高度不发生变化）.小尺度对流在等位温面的翻转和波饱和发生后产生,随后波包开始破碎,这些非线性过程的最终结果产生湍流.表明导致重力波饱和的关键因素是等位温面的翻转而不是诸如波破碎、湍流、波-流相互作用等其他的一些物理过程.     

风垂直切变对中尺度地形对流降水影响的研究

针对长江中下游中尺度地形特点以及暴雨过程发生发展期间风垂直切变的主要观测特征,设计了一系列中尺度地形的三维理想数值试验,分析了干大气地形流和重力波特征,探讨了条件不稳定湿大气地形对流降水的模态分布,在此基础上研究了圆形、直线风垂直切变和切变厚度对中尺度地形对流降水强度和模态分布的影响.结果发现:在 F_r≈1的干大气条件下,气流遇到地形后分支、绕流和爬升现象同时存在,地形激发的重力波在水平和垂直方向上传播,其在迎风坡、背风坡、地形上游和下游的振幅不同,并组织出不同强度的垂直上升运动.在F_r > 1的条件不稳定湿大气下,地形对流降水主要存在三种模态,即迎风坡和背风坡准静止对流降水以及地形下游移动性对流降水,地形对流降水的形成与重力波在低层组织的上升运动密切相关.风垂直切变对地形对流降水的强度和模态分布有重要作用,其中圆形风垂直切变(风随高度旋转)不仅影响地形下游对流降水系统的移动方向,而且影响迎风坡和背风坡山脚处对流降水中心的分布和强度;直线风垂直切变(风随高度无旋转)主要影响地形对流降水的移动速度和强度.风随高度自下而上顺(逆)时针旋转,地形对流系统向下游传播时向右(左)偏移.风垂直切变主要通过影响地形重力波的结构和传播以及对流系统的形成、移动方向和速度,来影响地形对流降水的模态分布,其中对流层中低层的风垂直切变对地形对流降水强度和模态分布有重要影响.     

与东北冷涡相伴的高空急流诱发平流层重力波的数值模拟研究

使用中尺度数值模式WRF-ARW,针对2010年6月发生在中国东北地区一例伴随对流层高空西风急流（位于~9 km高度）演变过程出现的平流层重力波活动特征开展了数值模拟. 事件发生期间,对流层区域环流处在一个东北冷涡系统的控制之下. 模拟结果再现了该东北冷涡的发展和维持过程,以及与之相伴的高空急流的特征. 模拟结果揭示出在急流区域上空的平流层中存在显著重力波活动现象. 分析结果显示,重力波活动与急流存在紧密联系,在水平方向上,重力波呈显著的二维结构,出现在急流出口区上部并逆背景流向西传播. 功率谱分析结果表明盛行波动具有~700 km水平尺度、9~12 h时间尺度以及4~5 km垂直波长. 由于急流的存在,造成其与平流层中下部之间存在显著的水平风速垂直切变,与切变相伴的耗散使得上传的重力波动量通量数值随着高度升高而递减. 同时,在18~20 km高度间出现的西风-东风转换带极大地抑制了波动在垂直方向的传播,形成显著动量通量沉积效应. 估算结果表明,在11~20 km高度之间,这种效应的整体作用相当于对该层背景流施加强度为0.86 m·s^-1·day^-1的动力阻曳.     

北极地区低平流层惯性重力波的观测研究

南极地区重力波活动有大量报道，相对而言，北极地区重力波的研究还很少.本文利用极区Ny-Alesund站点（78.9°N，11.9°E）无线电探空仪从2012年4月1日到2017年3月31日共5年的观测数据，统计分析了北极地区低平流层惯性重力波的特征.观测显示，月平均纬向风在20 km以下盛行东向风，再随着高度增加，逐渐呈现出半年振荡现象.对流层顶高度在5~13 km范围内变化，其月平均高度显示出年循环，最高出现在夏季，约为10 km，最低出现在冬季，约为8.5 km.对流层和低平流层月平均温度都显示出明显的年周期变化，这与中低纬度观测结果有所不同.结合Lomb-Scargle谱分析和矢端曲线方法，估算了准单色惯性重力波参数.个例研究表明，低平流层惯性重力波呈现出远离源区的自由传播性质.统计结果显示，惯性重力波的水平和垂直波长分别集中在50~450 km和1~4 km范围内，本征频率集中在1~2.5倍惯性频率间，这些值都比中低纬度观测值稍小.垂直方向本征相速度主要集中在-0.3~0 m·s^-1，而纬向和经向本征相速度集中在-40~40 m·s^-1之间.在5年的观测中，大约91.5%的惯性重力波向上传播.在冬季和早春，由于极地平流层极涡活动，激发出向下传播的惯性重力波，因此，向下传播的比例上升到相应月份的20%左右.由于低层大气盛行的东向风的滤波效应，低平流层大部分惯性重力波向西传播.波能量呈现出明显的年周期变化，最大值在冬季、最小值在夏季，与北半球中低纬度观测结果一致，表明北半球重力波活动普遍冬季强、夏季弱.     

基于多台站OH全天空气辉成像仪观测的中国中纬地区重力波传播特性

利用位于中国中纬地区6个OH气辉成像仪2012年1月至2013年12月两年的观测数据,我们研究分析了重力波传播特征.结果表明重力波的水平波长、观测周期和水平相速度分别主要分布于10~35km,4~12min和30~100m·s-1范围.夏季,重力波主要沿极向方向传播.然而,冬季,他们有向赤道方向和平行于赤道方向的传播趋势.同时,我们结合TRMM卫星和ECMWF数据,发现在夏季重力波的北向传播趋势可能主要由观测台站南方的对流活动导致.然而,对流层顶附近的急流可能在冬季重力波的主要传播方向方面做出较大贡献.分析结果也表明,低层-中层大气背景风的滤波效应仅在夏季与中国中纬地区重力波纬向传播方向各向异性吻合较好.     

利用瑞利激光雷达观测北京地区上平流层地形重力波活动

本文利用中国科学院空间科学与应用研究中心的瑞利激光雷达首次观测到了平流层地形重力波活动的现象,并结合美国国家环境预报中心(NCEP)的全球预报系统(GFS)的风场数据分析了该地形重力波的基本参数.与惯性重力波相比较,地形重力波的密度扰动没有下传的相位,在同一高度上,其扰动相位保持不变.北京空间科学与应用研究中心瑞利激光雷达自2012年开始观测实验以来,已经观测到多起地形重力波活动事件.本文以2013年11月11日的观测数据为例,研究北京上空的地形重力波活动,并结合GFS风场数据分析了北京上平流层地形重力波的波长、传播方向、传播速度等参量.通过分析得到在2013年11月11日北京上空存在一列传播方向为北偏西52.4°,水平波长为5.5km,平均垂直波长约为6.0km的地形重力波.     

青藏高原上空UTLS区域一次地形重力波过程中的物质上传

利用美国航空航天局MERRA(Modern-Era Retrospective Analysis for Research and Applications)再分析资料和MODIS(Moderate-Resolution Imaging Spectroradiometer)卫星资料以及欧洲气象中心ECMWF-Interim(European Centre for Medium-Range Weather Forecasts)再分析资料,分析了发生于青藏高原北侧上空的一次地形重力波事件,并使用中尺度预报模式WRF-ARW.V3.0(Weather Research and Forecasting model,V3.0)对其进行了数值模拟.在此基础上,诊断分析了此次地形重力波在UTLS(Upper Troposphere and Lower Stratosphere)区域造成的物质和能量垂直传输特征.分析结果表明这一中尺度地形重力波信号的水平波长约为600km,与地形扰动水平尺度接近,重力波在对流层中传播的垂直波长约为3km,在垂直方向上随着高度的增加呈现出由东向西倾斜的结构特征.此次地形重力波上传进入平流层并在150hPa附近破碎,波破碎后动量通量在短时间内发生了强烈的衰减,重力波携带的能量在破碎高度附近释放.重力波破碎的同时垂直方向湍流混合变得异常强烈,湍流交换系数可在短时间内增加到背景值的8倍以上,剧烈湍流混合过程导致了对流层上层的空气进入平流层,使下平流层空气出现了位势涡度和臭氧的低值区,在浮力频率的垂直剖面中也可以看到由于地形重力波过程造成的平流层下层浮力频率异常低值区.     

热源激发重力波特征以及波流作用的数值模拟研究

本文在二维等温可压大气中引入了一个随时间和空间变化的热源扰动,分别以静止风和中纬1月份月平均向东的纬向风急流为背景,对不同背景下热源激发的重力波的传播详细过程及其特性进行了数值模拟研究.热源激发出来的重力波在初始阶段有很宽的频谱范围,随后由于重力波的传播效应,水平波长和垂直波长分布范围随时间都有所减小.顺风传播的重力波...     

对流层上传重力波的非线性演化

利用二维全隐欧拉格式对重力波在可压、非等温大气中的非线性传播过程进行了数值模拟和分析.分析结果表明，从对流层顶激发的重力波能稳定地经平流层传到中层顶，从而将能量和动量从一个区域带到另一个区域；在向上传播过程中，重力波经历了发展、位温翻转、对流直至最终破碎的演变；重力波的破碎是对流和小尺度波动的重要的源，对流不稳定和翻转是非线性现象的一个基本特征.计算还显示，扰动源的大小直接影响着重力波的非线性传播过程，当扰动源足够小时，重力波能稳定传播，而大振幅扰动可以加速重力波的破碎.     

Parametrization of momentum and energy depositions from gravity waves generated by tropospheric hydrodynamic sources

The mechanism of generation of internal gravity waves (IGW) by mesoscale turbulence in the troposphere is considered. The equations that describe the generation of waves by hydrodynamic sources of momentum, heat and mass are derived. Calculations of amplitudes, wave energy fluxes, turbulent viscosities, and accelerations of the mean flow caused by IGWs generated in the troposphere are made. A comparison of different mechanisms of turbulence production in the atmosphere by IGWs shows that the nonlinear destruction of a primary IGW into a spectrum of secondary waves may provide additional dissipation of nonsatu-rated stable waves. The mean wind increases both the effectiveness of generation and dissipation of IGWs propagating in the direction of the wind. Competition of both effects may lead to the dominance of IGWs propagating upstream at long distances from tropospheric wave sources, and to the formation of eastward wave accelerations in summer and westward accelerations in winter near the mesopause.     

Linear mechanism by which internal gravity waves are generated and intensified in the ionosphere when interacting with an inhomogeneous zonal wind: 2. Internal gravity wave generation and intensification during the linear evolution stage

The linear mechanism by which internal gravity waves (IGWs) are generated and subsequently intensified in a stably stratified dissipative ionosphere in the presence of an inhomogeneous zonal wind (shear flow) has been studied. In the case of shear flows, the operators of linear problems are nonself-adjoint and the corresponding eigenfunctions are nonorthogonal; a canonical approach can hardly be used to study such motions. It is more adequate to apply the so-called nonmodal calculation. Dynamic equations and equations of energy transfer of IGW disturbances in the ionosphere with a shear flow have been obtained based on a nonmodal approach. Exact analytical solutions for the constructed dynamic equations have been found. The growth rate of the IGW shear instability has been determined. It has been established that IGW disturbances are intensified in an algebraically power manner rather than exponentially in the course of time. The effectiveness of the linear mechanism by which IGWs are intensified when interacting with an inhomogeneous zonal wind is analyzed. It has been indicated that IGWs effectively obtain the shear flow energy during the linear evolution stage and substantially increase (by an order of magnitude) their amplitude and energy. The frequency and the wave vector of generated IGW modes depend on time; therefore, a wide spectrum of wavelike disturbances, depending on the linear, rather than nonlinear, turbulent effects, is formed in the ionosphere with a shear flow. Thereby, a new degree of freedom appears, and the turbulent state of atmospheric—ionospheric layers can be formed on IGW disturbances.     

高纬电离层特性的实例研究

本文利用EISCAT雷达资料讨论磁层-电离层耦合的高纬电离层效应。研究表明,即使在夏季极昼情况下,磁扰期间的磁层过程对高纬电离层形态的影响也远大于太阳紫外辐射的作用。高能粒子沉降使电离层E层的电子密度大大增加;而磁层对流速度变大会使F层内电离复合加强。因而,磁扰时经常出现n_e（E层）>_e（F层）的情况。此现象不仅与宁静时完全相反,而且与中低纬电离层形态变化也有很大差别。     

Relative role of subinertial and superinertial modes in the coastal long wave response forced by the landfall of a tropical cyclone

A set of numerical experiments has been performed in order to analyze the long-wave response of the coastal ocean to a translating mesoscale atmospheric cyclone approaching the coastline at a normal angle. An idealized two-slope shelf topography is chosen. The model is forced by a radially symmetric atmospheric pressure perturbation with a corresponding gradient wind field. The cyclone's translation speed, radius, and the continental shelf width are considered as parameters whose impact on the long wave period, modal structure, and amplitude is studied. Subinertial continental shelf waves (CSW) dominate the response under typical forcing conditions and on the narrower shelves. They propagate in the downstream (in the sense of Kelvin wave propagation) direction. Superinertial edge wave modes have higher free surface amplitudes and faster phase speeds than the CSW modes. While potentially more dangerous, edge waves are not as common as subinertial shelf waves because their generation requires a wide, gently sloping shelf and a storm system translating at a relatively high (∼10 m s⁻¹ or faster) speed. A relatively smaller size of an atmospheric cyclone also favors edge wave generation. Edge waves with the highest amplitude (up to 60% of the forced storm surge) propagate upstream. They are produced by a storm system with an Eulerian time scale equal to the period of a zero-mode edge wave with the wavelength of the storm spatial scale. Large amplitude edge waves were generated during Hurricane Wilma's landfall (2005) on the West Florida shelf with particularly severe flooding occurring upstream of the landfall site.     

旋转流体中的二维惯性重力波-涡相互作用

本文利用f-平面浅水模型讨论了二维情况下快速惯性重力波与平衡涡旋之间的相互作用问题.当快速的惯性重力波传入涡旋区域后，波－涡之间发生了相互作用，涡旋可被加深，涡旋区域的风场等可以出现非对称分布，当惯性重力波全部传出涡旋区域并远离之后，涡旋恢复到其初始状态，但此时惯性重力波可以出现非对称结构，这与一维情况不同.在波－涡发生相互作用过程中，其物理量变化的最大区域分布基本与波动传播方向一致，在波动进入或移出涡旋时，相对应存在两个最大物理量变化区域，此外，相互作用导致的物理量的最大变化与波动结构、强度、传播方向以及涡旋特性有关.非线性条件下的波－涡相互作用要强于线性过程.     

Seasonal variation of space–time parameters of internal gravity waves at Kharkiv (49°30′N, 36°51′E)

Internal gravity waves (IGWs) over Kharkiv (49°30′N, 36°51′E) have been studied by using an automatic goniometer of a meteor radar (AG MR), with its antenna directed to the East. The AG MR carries out Doppler measurements of the radial drift velocity and the position of the reflecting area of meteor trails. In order to obtain information about IGW parameters, an algorithm was used that processes the AG MR wind data by dividing the measured volume into several sub-volumes, and by carrying out wavelet analysis of wind variations in these areas. Results of 1 year (1987) of AG MR data show a good qualitative correspondence with literature results. There is an increase of IGW activity during the change of the zonal background wind in spring and autumn. It was found that the mean IGW horizontal phase velocity and predominating horizontal propagation direction change with season. However, the mean period (1.5 h), the dominant amplitude (30 m/s) and the vertical phase velocity do not strongly vary in the course of the year. The mean vertical and horizontal wavelengths were found to be 40 and 250 km, respectively.     

Comparison of Early Evolutions of Mimas and Enceladus

Thermal history of Mimas and Enceladus is investigated from the beginning of accretion to 400 Myr. The numerical model of convection combined with the parameterized theory is used. The following heat sources are included: short lived and long lived radioactive isotopes, accretion, serpentinization, and phase changes. The heat transfer processes are: conduction, solid state convection, and liquid state convection. We find that temperature of Mimas’ interior was significantly lower than that of Enceladus. If Mimas accreted 1.8 Myr after CAI then the internal melting and differentiation did not occur at all. Comparison of thermal models of Mimas and Enceladus indicates that conditions favorable for the start of tidal heating lasted for a short time (~10⁷ yr) in Mimas and for ~10⁸ yr in Enceladus. This could explain the Mimas—Enceladus paradox. In fact, in view of the chronology based on cometary impact rate, one cannot discard a possibility that also Mimas was for some time active and it has the interior differentiated on porous core and icy mantle.     

Field observations of infragravity waves and their behaviour on rock shore platforms

Infragravity wave (IGW) transformation was quantified from field measurements on two shore platforms on New Zealand's east coast, making this the first study to describe the presence, characteristics and behaviour of IGWs on rock platform coasts. Data was collected using a cross‐shore array of pressure transducers during a 22 hour experiment on Oraka shore platform and a 36 hour experiment at Rothesay Bay shore platform. A low pass Fourier filter was used to remove gravity wave frequency oscillations, allowing separate analysis of IGWs and the full wave spectrum. Offshore IGW heights were measured to be 7 cm (Oraka) and 9 cm (Rothesay Bay), which were 21% (Oraka) and 7.5% (Rothesay Bay) the height of incident wave height. At the cliff toe, significant IGW height averaged 15 cm at Oraka and 13 cm at Rothesay Bay. This increase in IGW height over the platform during both experiments is attributed to shoaling of 40 to 55% over the last 50–60 m before the cliff toe, respectively. Shoaling across the platform was quantified as the change in IGW height from the platform edge to cliff toe, resulting in a maximum increase of 1·88 and 2·63 on Rothesay Bay and Oraka platforms. IGW height at the cliff toe showed a strong correlation with incident wave height. The proportional increase in IGW height shows a strong correlation to water level on each platform. The rate of shoaling of long period waves on the shallow, horizontal platforms increased at higher water levels resulting in a super elevation in water level at the cliff toe during high tide. Greater IGW shoaling was also observed on the wider (Oraka) shore platform. Results from this study show the first measurements of IGWs on shore platforms and identify long wave motion a significant process in a morphodynamic understanding of rock coast. Copyright © 2011 John Wiley & Sons, Ltd.