阵风因子与大气边界层要素的关系及预报试验

利用2008—2016年舟山4个海岛气象站大风资料和欧洲中心的ERA-interim再分析资料,分析了阵风因子随平均风速、风向、小时、月份等分布的气候特征,统计阵风因子与边界层的大气稳定度、250～1000m风速与10m风速的比值、6h变温等要素的相关后,选取最佳预报因子,利用BP人工神经网络方法,根据不同因子组合对阵风进行循环试报。结果表明:①平均风速较小时阵风因子波动范围大;靠近大陆站点的阵风因子及来自陆地方向气流的阵风因子偏大。②白天11:00—16:00受太阳辐射影响大气湍流相对较强,阵风因子偏大。7—9月沿海受台风影响频繁,其阵风因子要大些,而11—12月阵风因子偏大则与来流方向的地表粗糙度较大有关。③阵风因子和边界层不同高度的风速与10m风速比值,及气温具有明显正相关,与边界层大气稳定度参数具有负相关,相关结果印证了阵风主要物理成因与动量的垂直湍流输送有关。④阵风循环试报表明最佳组模型试报的绝对误差及方差均比对比组模型减少约11%～25%,具有较好的预报效果。     

大气边界层强风的阵性和相干结构

我国北方春季冷锋过境后,常骤发强风,甚至起沙扬尘,持续数小时甚至一二天,通过对边界层超声风温仪的资料分析,可知大风常叠加有周期为3～6 min的阵风,较有规律,且有明显的相干结构:阵风风速峰期有下沉运动,谷期有上升运动;阵风扰动以沿平均流的顺风方向分量为主,横向和垂直方向的分量都较小,其本质是低频次声波和重力波的混合;阵风沿顺风向且向下传播.周期小于1 min的脉动在水平面上基本是各向同性的不规则的湍涡.大风期间,无论是平均流、阵风和湍流脉动,至少在120 m高度以下,主要都有西风和北风动量下传,感热上传.平均流的动量下传强于由脉动下传的量,与一般天气情况不同,而且阵风与湍流的动量下传的量值差不多.平均流和阵风在动量传送上起相当大的作用.     

大气边界层阵风扬尘机理

东亚春季冷锋后的强风及其伴随的系统下沉气流和叠加其上的阵风,都对当地起沙十分有利.但此系统的下沉气流却压制着沙尘,使之不能上扬到大气边界层内.经过仔细计算和分析发现,正是阵风的三维相干结构,使得积聚在风沙边界层的沙尘可以持续地克服系统的下沉气流,而不断上扬到大气边界层的中上层.再藉中上层和其上系统的上升气流,将沙尘带到对流层,扩散到更大范围和传送到远方.这是东亚春季强风起沙扬尘最常见和最有效的机理.     

大气边界层湍流相干结构的识别

首先利用数字滤波方法对淮河流域试验的大气边界层湍流观测资料进行三项分解，将大气边界层湍流的风速信号分解为近似各项同性的小尺度涡和各向异性的大尺度涡。然后再将大尺度涡信号进行离散正交小波分解，寻求相干结构的主要特征尺度。对于大气边界层湍流垂直脉动风速来说，其相干结构的主要特征尺度为16s；对径向与纬向脉动来说，其相干结构的主要特征尺度为32－64s。在此基础上，利用小波的反变换提取出相干结构的信号与非相干结构的信号，并计算两间的相关系数，最大仅有0．02。此外，对原始大气湍流观测信号不进行数字滤波，直接利用本中子波分析法提取湍流相干结构所得结果作比较研究；并探讨了采用对称或似对称离散正交小波对此研究的影响。     

强风天气下边界层结构特征

近地层观测的强风运动表明,叠加在平均流动之上的脉动通常有两种,一种是随机的湍流脉动,还有一种具有相干结构的阵风扰动。分析表明,上层强风的剪切运动产生阵风,并向下传递能量,对近地层的通量传输起到重要作用。本文利用北京325 m气象塔、位于海拔1257 m的妙峰山测风塔和位于海拔1688 m的灵山测风塔的资料,分析了强风天气下,边界层上层出现阵风并向下传递的过程,进一步证实无论在近地层还是边界层上层,强风期间,叠加在平均流动上除了高频湍流脉动之外,还有周期为1~10分钟的阵风,即相干结构。阵风峰期有下沉运动,阵风谷期有上升运动。这些相干结构在边界层上层产生,向下运动和传播过程中受到平均气流梯度的切变作用和地面摩擦,破碎为湍流结构。边界层上层的阵风和湍流产生的动量通量向下传递,使得强风期间,边界层中阵风和湍流对通量具有同样的输送能力,对边界层中沙尘、污染物等气溶胶的传输具有重要作用。本研究为模式中进行通量输送参数化方案的修正提供了观测和理论依据。     

沿海地区大气边界层扩散特征探讨

利用辽东湾近地层三轴风速仪资料分别计算不同稳定度的大气扩散参数及湍流强度。并将不同情况下的大气扩散参数及湍流强度进行了比较。结果表明：在地形相对平坦的辽东湾,在不同的稳定度下,σy和σx值均略大于同级布里格斯公式的计算值。在不稳定条件(C、B类)下,湍强值重大。     

重庆冬季大气边界层湍流扩散能力的初步研究

本文利用1989年12月16日至1990年1月14日在重庆市用三轴风速仪所测得到的三轴风速资料，分别计算了有风、静风、有雾时的大气扩散参数， 以及无雾有风条件下的湍流强度。文中对有风时、静风时、有雾时的扩散参数作了比较和研究，对无雾有风时的湍流强度作了分析。从而有助于了解重庆在有风、 静风或有雾时边界层内湍流扩散能力的大小。     

辽东湾地区大气边界层扩散特征探讨

利用辽东湾地区近地层三轴风速仪资料,分别计算了不同稳定度下的大气扩散参数及湍流强度。并对不同情况下的大气扩散参数及湍流强度进行了比较。结果表明,在地形相对平坦的辽东湾地区,在不同的稳定度下,σy和σx值均略大于同级Briggs公式的计算值,但小于地形相对复杂的本溪地区。在不稳定条件(C、B类)下,湍强值最大。     

近壁湍流边界层相干结构形成的动力学分析

以壁面脉冲扰动来构建湍流边界层近壁区对称与非对称相干结构理论模式,采用直接数值模拟的方法,研究了湍流边界层近壁区对称与非对称相干结构形成的理论机制和演化规律.计算结果显示非对称相干结构比对称相干结构更容易增长,这是由于非对称型结构能在湍流边界层近壁区域诱导产生较大的展向扰动速度,且不同初始结构类型分布对湍流边界层近壁区相干结构形成的动力学机制不同所致,从而加深了对近壁湍流边界层单个相干结构形成的某种动力学机制的认识理解.     

The Coherent Structure of Water Vapour Transfer in the Unstable Atmospheric Surface Layer

Observations of water vapour fluctuations over arice field show vapour ramps. Coherent structuresare first revealed by the frequently occurring ramp pattern in the vapourtrace. Wavelet and pseudo-wavelet analysis techniques were used inconditional sampling, and more than 100 hr of data have been analyzedto determine coherent structure characteristics. The most probablecoherent structure duration was in the range 2–12 sec andthe duration range of the most effective coherent structures shows somedifference between heat and water vapour transfers. Coherent structurescontribute to the major part of the total flux.     

Simulation of a Stably Stratified Atmospheric Boundary Layer Using One-Dimensional Turbulence

One-dimensional turbulence (ODT) is a single-column simulation in which vertical motions are represented by an unsteady advective process, rather than their customary representation by a diffusive process. No space or time averaging of mesh-resolved motions is invoked. Molecular-transport scales can be resolved in ODT simulations of laboratory-scale flows, but this resolution of these scales is prohibitively expensive in ODT simulations of the atmospheric boundary layer (ABL), except possibly in small subregions of a non-uniform mesh.Here, two methods for ODT simulation of the ABL on uniform meshes are described and applied to the GABLS (GEWEX Atmospheric Boundary Layer Study; GEWEX is the Global Energy and Water Cycle Experiment) stable boundary-layer intercomparison case. One method involves resolution of the roughness scale using a fixed eddy viscosity to represent subgrid motions. The other method, which is implemented at lower spatial resolution, involves a variable eddy viscosity determined by the local mesh-resolved flow, as in multi-dimensional large-eddy simulation (LES). When run at typical LES resolution, it reproduces some of the key high-resolution results, but its fidelity is lower in some important respects. It is concluded that a more elaborate empirically based representation of the subgrid physics, closely analogous to closures currently employed in LES of the ABL, might improve its performance substantially, yielding a cost-effective ABL simulation tool. Prospects for further application of ODT to the ABL, including possible use of ODT as a near-surface subgrid closure framework for general circulation modeling, are assessed.     

On the Possible Impact of a Following-Swell on the Atmospheric Boundary Layer

A simple model of the atmospheric boundary layer over the ocean where the swell impact on the atmosphere is explicitly accounted for is suggested. The model is based on Ekman’s equations, where the stress in the wave boundary layer is split into two parts: the turbulent and wave-induced stress. The turbulent stress is parameterized traditionally via the eddy viscosity proportional to the generalized mixing length. The wave-induced stress directed upward (from swell to the atmosphere) is parameterized using the formalism of the wind-over-waves coupling theory. The model can be seen as an extension of the model by Kudryavtsev and Makin (J Phys Oceanogr 34:934–949, 2004) to the scale of the entire atmospheric boundary layer by including the Coriolis force into the momentum conservation equation and generalizing the definition of the mixing length. The regime of low winds for swell propagating along the wind direction is studied. It is shown that the impact of swell on the atmosphere is governed mainly by the swell parameter—the coupling parameter that is the product of the swell steepness and the growth rate coefficient. When the coupling parameter drops below − 1 the impact of swell becomes significant and affects the entire atmospheric boundary layer. The turbulent stress is enhanced near the surface as compared to the no-swell case, and becomes negative above the height of the inner region. The wind profile is characterized by a positive gradient near the surface and a negative gradient above the height of the inner region forming a characteristic bump at the height of the inner region. Results of the model agree at least qualitatively with observations performed in the atmosphere in presence of swell.     

On Modelling the One-Dimensional Atmospheric Boundary Layer

Several commonly used turbulence closure schemes forthe atmospheric boundarylayer (ABL) are applied to simulate neutral, nocturnal and diurnal cyclesituations in a one-dimensional ABL. Results obtained with the differentschemes, E-, E- and its modified versions, and twoversions ofthe q ² Level 2.5, are compared and discussed.     

TRAC98: Detection of Coherent Structures in a Convective Boundary Layer using Airborne Measurements

The TRAC98 experimental campaign (Turbulence Radar Aircraft Cells) devoted to coherent structures analysis took place over the Beauce plain (France) during summer 1998. It allowed us to collect a large dataset of airborne measurements in addition to various ground measurements. This study aims at diagnosing the occurrence of coherent structures within the atmospheric boundary layer (ABL) through airborne measurements. The statistical analysis performed as a first step from turbulent parameters underlined the homogeneity of the ABL over the Beauce plain. However mixed-layer scaling failed at the top of the ABL, even when taking into account the entrainment rate. Coherent structures were detected through the analysis of ABL isotropy, using the opportunity of sampling with two perpendicular crossing planes, one of them being aligned with the wind. This approach allowed us to determine an organization scheme of the ABL for three of the five flights (ARAT30, MIV30 and MIV27). For the ARAT30 flight, the analysis was pursued by focusing on measurements of fluctuations in the inner flight legs. In this way, the low-level cloud cover has been investigated from the downward visible radiation (VISD). The results indicated an anisotropy of the horizontal cloud size. Secondly, the variations of some parameters were analysed through lagged correlation functions. This allowed us to infer relationships between the vertical velocity, VISD, mixing ratio and lifting condensation level. Length scales have also been extracted, and confirmed the ABL organization during the ARAT30 flight. Finally, the anisotropy observed in various flights has been investigated with respect to the underestimation of the latent heat fluxes revealed by the imbalance of the surface energy budget.     

A Scale-Dependent Dynamic Model for Scalar Transport in Large-Eddy Simulations of the Atmospheric Boundary Layer

An important challenge in large-eddy simulationsof the atmospheric boundarylayer is the specification of the subgrid-scale(SGS) model coefficient(s)and, in particular, how to account for factorssuch as position in the flow,grid/filter scale and atmospheric stability.A dynamic SGS model (thatassumes scale invariance of the coefficients)is implemented in simulationsof a neutral boundary layer with a constantand uniform surface flux of apassive scalar. Results from our simulationsshow evidence that the lumpedcoefficient in the eddy-diffusion modelcomputed with the dynamic proceduredepends on scale. This scale dependence isstronger near the surface, and itis more important for the scalar than for thevelocity field (Smagorinskycoefficient) due to the stronger anisotropicbehaviour of scalars. Based onthese results, a new scale-dependent dynamicmodel is developed for theeddy-diffusion lumped coefficient. The newmodel, which is similar to theone proposed earlierfor the Smagorinsky coefficient,is fully dynamic, thus not requiring anyparameter specification or tuning.Simulations with the scale-dependent dynamicmodel yield the expected trendsof the coefficients as functions of positionand filter/grid scale.Furthermore, in the surface layer the newmodel gives improved predictionsof mean profiles and turbulence spectra ascompared with the traditionalscale-invariant dynamic model.     

北京冬季大风过程大气边界层特征分析

利用北京中国科学院大气物理研究所325 m气象观测塔的气象梯度资料和湍流资料,分析了2014年11月29日至12月5日北京两次大风过程中气象要素和湍流输送特征的变化。第一次大风过程的强度和持续时间均高于第二次大风过程。强烈的风速垂直切变主要集中在距地面100 m高度范围内,最强风速垂直切变达到0.31 s~(-1)。大风过程中,阵风系数呈现随高度减小的趋势,越接近地面,阵风系数愈大。阵风强度的变化与阵风系数相似,100 m以下高度时,阵风强度随高度增大而减小。大风过程自上而下改变边界层结构,平均动能、湍流动能和摩擦速度最先从上层(280 m)发生变化且迅速增加。近地层由于风速垂直梯度的显著差异,近地层垂直方向的湍流强度最大。大风时各功率谱在低频区(0.01 s~(-1))达到峰值,大风过后各高度的能量都有所下降。