大气中的内船舶波

通过波动方程的分析解求得稳定层结大气中孤立地形激发的大气内船舶波,证明这类波是由非拦截波和拦截波两部分构成,非拦截波中的辐散波和横波的形态,王要由地形形状和走向决定,它们都可向大气上层传输,而拦截背风波的波形主要由层结、风速及风速切变等大气条件决定.     

一个大气边界层非线性三维数值模式及其在太湖流域的应用

本模式用以K半经验理论为基础的闭合方程组求解,引进了地形作用、动量、热量和水汽的侧向扩散和垂直湍流交换以及凝结加热等物理过程,其中湍流交换系数K是地表粗糙度、几何高度、大气稳定度及风速切变等因子的函数。本模式还考虑了山地、平原和水面在地形高度、地表粗糙度、辐射、蒸发及热交换等方面的差异,在下垫面上建立热平衡方程与边界层内控制方程耦合。将模式应用于太湖流域,计算得到的该流域边界层内温、压、湿、风的分布特征与实际情况相似,模式具有一定的实用性。     

一个大气边界层非线性三维数值模式及其在太湖流域的应用

本模式用以K半经验理论为基础的闭合方程组求解,引进了地形作用、动量、热量和水汽的侧向扩散和垂直湍流交换以及凝结加热等物理过程,其中湍流交换系数K是地表粗糙度、几何高度、大气稳定度及风速切变等因子的函数。本模式还考虑了山地、平原和水面在地形高度、地表粗糙度、辐射、蒸发及热交换等方面的差异,在下垫面上建立热平衡方程与边界层内控制方程耦合。将模式应用于太湖流域,计算得到的该流域边界层内温、压、湿、风的分布特征与实际情况相似,模式具有一定的实用性。     

《地球物理学报》总目次

0001磁暴期间几种主要磁扰成分的演化特征………吴迎燕,徐文耀,陈耿雄,陈博,…刘晓灿0010由非相干散射雷达数据重建极光沉降粒子能谱…………………………………蔡红涛,…马淑英0018电离层暴时的foF2预报技术研究…………………………陈丹俊,吴健,…王先义0024重力波非共振激发的数值研究…………………………黄开明,张绍东,…易帆0034三维分层流动过双山脉地形激发的大气船舶重力波动力学理论和数值试验………………………李子良0043夜间稳定边界层中小尺度地形激发的形式阻力和波动阻力……………姜瑜君,桑建国,刘辉志,…刘树华0051平流…     

MYJ和YSU方案对WRF边界层气象要素模拟的影响

研究新一代中尺度气象模式WRF中两种大气边界层方案(MYJ,YSU)对沈阳冬季大气边界层结构模拟的影响,重点分析温度层结、低层风场、边界层高度等对污染物扩散有重要影响的气象要素.和观测数据的比对表明WRF基本能够模拟出温度风速的日变化特征,但模拟风速偏大.YSU方案由于模拟的边界层顶卷挟和边界层内混合作用较强,夜间接地逆温强度低于MYJ方案,逆温维持时间比MYJ方案短4小时,同时模拟边界层高度也高于MYJ方案,有利于污染物垂直扩散.边界层高度的3种计算方法中,湍流动能方法计算的边界层高度最高,Richardson数方法次之,位温方法得到的高度最低.Richardson数方法对临界值的选取较敏感.     

大气边界层湍流的动力非平稳性的验证

首次用验证时间序列中是否存在动力非平稳性的一种简单图示方法——space time index法来分析大气边界层湍流的动力平稳性特征.本文以取自淮河流域和威斯康星森林下垫面条件下的三维高精度风速和温度、湿度湍流脉动资料对大气边界层湍流的平稳性特征进行了分析.结果表明space time index方法能有效地检验大气边界层湍流信号中是否存在动力平稳性.另外,均匀下垫面条件（水稻田）及复杂下垫面条件（森林）下的大气边界层湍流信号中几乎都存在动力非平稳性,大气湍流动力学非平稳性可能是边界层湍流信号相当普遍具有的一种特性.大气边界层湍流中的间歇性和相干结构使得其非平稳性图形的特征不同于一般时间序列非平稳性图形的“V”型特征;森林下垫面条件下的湍流信号比相对均匀下垫面（水稻田）下的湍流信号更有组织性,相干结构更强.     

我国干旱区深厚大气边界层与陆面热力过程的关系研究

大气边界层是地球大气重要物理特征之一.干旱气候背景下的大气边界层特征与人们以往对典型大气边界层的认识有很大不同,其形成机制也比较特殊.本文利用冬、夏两季典型时段在极端干旱区敦煌开展的大气边界层和陆面过程综合观测试验,对大气边界层厚度与净辐射、地-气温差和感热通量等陆面热力因素的日变化特征进行了对比分析,研究了大气边界层的发展和维持衰退过程与陆面热力和动力过程的关系,发现发展过程消耗的能量要比维持衰退过程高得多,而且进入残余层后对流边界层发展对陆面热力作用更敏感,发展也更迅速.风速的动力作用对大气边界层发展也有一定的影响,尤其对冬季稳定大气边界层影响较大.极强的陆面热力作用是我国干旱区形成深厚大气边界层的主导因素.     

大气边界层研究进展

大气边界层对云和对流的发展、演变有重要作用.本文回顾了在大气边界层高度计算方法,边界层的时空分布特征、结构和发展机理,以及边界层参数化方案等方面的主要研究进展.大气边界层高度计算方法主要分为基于大气廓线观测数据计算和基于模式参数化方案计算两大类;大气边界层高度频率分布形态具有明显的日变化特征,并且稳定、中性和对流边界层高度的频率分布呈现出不同的Gamma分布特征;地面湿度状况对边界层发展影响明显,对于不同的下垫面热力性质和地形状况,大气边界层高度呈现出明显的空间差异,青藏高原边界层高度明显高于一般平原地区;在强烈的地面加热驱动下,对流边界层与残余层通过正反馈机制循环增长可以形成4000 m以上的超高大气边界层;研制大气边界层、浅对流以及云物理方案的统一参数化框架是未来数值预报模式的发展趋势.     

北冰洋及其邻近海域极昼期间大气边界层结构特征试验研究

1年7-月中国"雪龙"号破冰船首航北极科学考察期间,利用TMT(系留式气象塔)大气边界层探测系统,对北冰洋及其邻近海域不同下垫面进行较长时间的连续多次观测. 观测到的大气边界层廓线,其中包括比湿廓线的多层、逆湿结构(逆湿强度最强为0.7g/(10m·kg));风向、风速廓线的切变结构(风向切变为1.1°/m,风速切变为0.11(m·s^-1)/m. 特别是持续数日、厚达几百米、其平均逆温强度有时高达1.4℃/10m的逆温结构,这种北冰洋海域特有的大气边界层结构,与地球气候系统中其他圈层的大气边界层结构有明显差别. 据此,提出了对该海域大气-海冰-海洋间动量、感热和潜热等湍流通量垂直交换以及热量平衡等有重要影响的大气逆温的屏障过程. 这为研究北极地区对全球气候影响、模拟北极地区现代气候和未来气候等提供一条新的思路. 此外,文中还给出该海域不同下垫面稳定大气边界层的高度.     

广州地区一次严重灰霾过程的垂直探测

利用激光雷达和微波辐射计在广州地区一次严重灰霾过程进行探测,采用小波变换协方差方法反演边界层高度,从热力和动力的角度分析边界层气象因子和能见度之间的关系.结果表明,边界层高度有明显的日变化过程,和地面能见度变化有很好的一致性.在清洁过程,边界层高度超过1 km,发生严重灰霾时,最高仅为500 m.50-100 m的温度梯度滞后30 h和能见度相关性更加显著,相关系数达到0.77.高层能见度（255 m）和下层稳定度有显著的负相关,清洁日负相关系数最大达-0.76,灰霾日负相关系数最大达-0.49.在相关的边界层气象因子中,地表通风系数和地面能见度成线性关系,相关系数最高,达0.88.边界层高度、地面风速、相对湿度和地面能见度的相关系数分别为0.76,0.67和-0.77.各气象因子之间具有较高的相关性.此次灰霾过程的主控气象因子为地表通风系数,在没有边界层高度探测的地区可利用地面能见度和风速估算边界层高度.     

Estimation of orographically induced wave drag in the stable boundary layer during the CASES-99 experimental campaign

This paper addresses the quantification of gravity wave drag due to small hills in the stable boundary layer. A single column atmospheric model is used to forecast wind and temperature profiles in the boundary layer. Next, these profiles are used to calculate vertical profiles of gravity wave drag. Climatology of wave drag magnitude and “wave drag events” is presented for the CASES-99 experimental campaign. It is found that gravity wave drag events occur for several relatively calm nights, and that the wave drag is then of equivalent magnitude as the turbulent drag. We also illustrate that wave drag events modify the wind speed sufficiently to substantially change the surface sensible heat flux.     

Effects of surface drag on low-level frontogenesis within baroclinic waves

Using a three-dimensional nonhydrostatic mesoscale numerical model (MM5), the evolution and structures of baroclinic waves with and without surface drag in case of dry and moist atmosphere are simulated, with special emphases on the effects of surface drag on the low-level frontal structure and frontogenesis. There are two different effects of surface drag on the low-level frontogenesis in the dry case. On one hand, the surface drag weakens the low-level frontogenesis and less inclined to develop the baroclinic wave due to the dissipation. But on the other hand, the surface drag induces a strong ageostrophic flow, which prolongs the low-level frontogenesis and finally leads to the enhancement of cold front. Compared with the no surface drag case, the surface drag increases the frontal slope espe- cially in the boundary layer, where the front is almost vertical to the surface, and then enhances the prefrontal vertical motion. All these conclusions expanded the analytical theory of Tan and Wu (1990). In the moist atmosphere, the influence of surface drag on frontal rainbands is also obvious. The surface drag weakens the convection, and reduces the energy dissipation near the surface when the initial relative humidity is relatively weak. At this time, the confluence induced post-frontal updrafts moves across the cold front and reinforces the prefrontal convection, which is beneficial to the maintenance of the rainband in cold sector. Given the enhancement of relative humidity, the moist convection domi- nates the low-level frontogenesis while the retardation of surface drag on energy dissipation is not obvious, therefore the effects of surface drag on the low-level frontogenesis and precipitation are re- duced.     

On the variability of the Charnock constant and the functional dependence of the drag coefficient on wind speed

A model for the air–sea interface, based on the coupled pair of similarity relations for “aerodynamically” rough flow in both fluids, is presented, which is applied to fetch-limited and high wind speed conditions which occur, for example, in hurricanes. It is shown that the specification of the maximum 10-m drag coefficient and the 10-m wind speed and the peak wave speed at which it occurs are sufficient to uniquely determine the drag law, which asymptotes at low wind speeds to a Charnock constant similar to that for the fully developed wind wave sea and is almost independent of the peak wave speed at the maximum in drag coefficient. A feature of the drag law is that it is of Charnock form, almost independent of the wave age, consistent with the transfer of momentum to the wave spectrum being due to the smaller rather than the dominant wavelengths. The analysis is also applied to a variable sea state in which either the surface wind or the surface Stokes drift vary, but the peak wave speed is kept constant. The corresponding variability in the Charnock constant is in general accord with observations.     

On the atmospheric internal ship waves

The analytical solutions of the atmospheric internal ship waves induced by three-dimensional terrain are obtained by solving the atmospheric wave equation. The solutions show that the waves consist of the untrapped and trapped parts. The patterns of the diverging wave and transverse wave in the untrapped parts are mainly determined by the shape and orientation of the terrain. This kind of wave may transport the wave energy to the upper atmosphere. The patterns of trapped lee waves are decided by the atmospheric conditions such as stratification, mean wind speeds and wind shear.     

A theory of enhanced saturation of the gravity wave spectrum due to increases in atmospheric stability

In this paper we consider a vertical wavenumber spectrum of vertically propagating gravity waves impinging on a rapid increase in atmospheric stability. If the high-wavenumber range is saturated below the increase, as is usually observed, then the compression of vertical scales as the waves enter a region of higher stability results in that range becoming supersaturated, that is, the spectral amplitude becomes larger than the saturation limit. The supersaturated wave energy must then dissipate in a vertical distance of the order of a wavelength, resulting in an enhanced turbulent energy dissipation rate. If the wave spectrum is azimuthally anisotropic, the dissipation also results in an enhanced vertical divergence of the vertical flux of horizontal momentum and enhanced wave drag in the same region. Estimates of the enhanced dissipation rates and radar reflectivities appear to be consistent with the enhancements observed near the high-latitude summer mesopause. Estimates of the enhanced mean flow acceleration appear to be consistent with the wave drag that is needed near the tropopause and the high-latitude summer mesopause in large-scale models of the atmosphere. Thus, this process may play a significant role in determining the global effects of gravity waves on the large-scale circulation.     

黄土高原复杂地形上边界层低空急流对近地层湍流的影响

利用中尺度气象数值模式（Weather Research and Forecasting Model,WRF）模拟风场,结合兰州大学半干旱气候与环境观测站（Semi-Arid Climate and Environment Observatory of Lanzhou University,SACOL）湍流观测资料,分析了黄土高原复杂地形上稳定边界层低空急流对近地层湍流活动的影响.黄土高原复杂地形上稳定边界层低空急流的形成与地形作用引发的局地环流有关.低空急流对近地层湍流活动有强烈影响,剪切作用使小尺度湍涡活动加剧,湍动能增大,同时非平稳运动被压制.低空急流发生时,观测数据有87.3%是弱稳定情形（梯度理查森数小于0.25）;而无低空急流时,对应时段的观测表明65.4%属于强稳定层结（梯度理查森数大于0.3）,非平稳运动造成湍流功率谱在低频端迅速增大.与无低空急流和弱低空急流情形相比,强低空急流发生时,近地层湍动能增大1倍,湍动能在垂直方向上的传递增大1个量级,且方向向下,约为-3 × 10^-3 m³·s^-3,湍流在上层产生并向下传递.     

Experimental investigation of drag coefficients of gobi surfaces

The response of gobi surfaces to the near-surface air flow can be characterized quantitatively by drag coefficients. By using wind tunnel tests, an attempt is made to define the relationship between the drag coefficients of gobi surfaces and gravel size and coverage. It is concluded that the drag coefficients of gobi surfaces tend to be constants when gravel coverage is over 40%-50%. Consequently, we think that the gobi deflation planes expanding vastly in the arid Northwestern China are aerodynamically stable, at least not the supplying sources of current dust storms, and therefore the emphasis on dust storm control should be paid on the so-called "earth gobi" that has low gravel coverage. The prediction model for drag coefficients of gobi surfaces has been developed by regressing drag coefficients on gravel size and coverage, the predicted results are in reasonably good agreement with wind tunnel results (R 2 = 0.94). The change of drag coefficients with gravel friction Reynolds number implies that the development extent of drag effect increases with gravel size and coverage.     

Correlation equation for the marine drag coefficient and wave steepness

This work questions, starting from dimensional considerations, the generality of the belief that the marine drag coefficient levels off with increasing wind speed. Dimensional analysis shows that the drag coefficient scales with the wave steepness as opposed to a wave-age scaling. A correlation equation is employed here that uses wave steepness scaling at low aspect ratios (inverse wave steepnesses) and a constant drag coefficient at high aspect ratios. Invoked in support of the correlation are measurements sourced from the literature and at the FINO1 platform in the North Sea. The correlation equation is then applied to measurements recorded from buoys during the passage of hurricanes Rita, Katrina (2005) and Ike (2008). Results show that the correlation equation anticipates the expected levelling off in deeper water, but a drag coefficient more consistent with a Charnock type relation is also possible in more shallower water. Some suggestions are made for proceeding with a higher-order analysis than that conducted here.     

On the variability of the Charnock constant and the functional dependence of the drag coefficient on wind speed: Part II-Observations

An analytical expression for the 10 m drag law in terms of the 10 m wind speed at the maximum in the 10 m drag coefficient, and the Charnock constant is presented, which is based on the results obtained from a model of the air-sea interface derived in Bye et al. (2010). This drag law is almost independent of wave age and over the mid-range of wind speeds (5?17 ms^?1) is very similar to the drag law based on observed data presented in Foreman and Emeis (2010). The linear fit of the observed data which incorporates a constant into the traditional definition of the drag coefficient is shown to arise to first-order as a consequence of the momentum exchange across the air-sea boundary layer brought about by wave generation and spray production which are explicitly represented in the theoretical model.     

藏南羊卓雍错湖面大气湍流特征观测分析

湍流运动是大气边界层的本质特征,是地表与大气之间能量和物质交换的主要方式.本文利用2016和2017年4-10月藏南羊卓雍错湖泊涡动观测资料,分析了湖面大气湍流方差和湍流特征量的统计和变化特征.结果表明：（1）不稳定层结下,三维风速分量和超声虚温、水汽密度、CO₂密度的无量纲标准差随稳定度变化符合Monin-Obukhov相似理论的"1/3"或"-1/3"次幂律,垂直风速的拟合效果最好;稳定层结下,除CO₂密度无量纲标准差与稳定度无明显关系外,其他量基本上满足相似性规律;中性条件下,以上物理量的无量纲标准差分别趋近常数：3.57、3.93、0.77、20.91、6.35和11.96.（2）水平方向平均湍流强度（0.60和0.58）大于垂直方向（0.13）,三维方向湍流强度与平均风速的变化呈显著负相关,相关系数分别为-0.39、-0.42和-0.34.（3）湖面湍流动能随风速呈线性增长,增长率达0.45 m/s;近中性层结时湍流动能最大,层结越稳定或不稳定湍流动能均减小.（4）湖泊下午到傍晚动量输送较强,13：00-22：30时间段平均动量通量达0.091 kg/（m·s²）;热量输送以潜热为主,潜热通量日平均值（77.3 W/m²）是感热通量（14.6 W/m²）的5.3倍,感热和潜热通量日变化峰值分别出现在5：30（22.4 W/m²）和16：00（106.6 W/m²）.