大气边界层阵风相干结构的产生条件

壁湍流相干结构的发现是近代湍流研究的重大进展之一,从20世纪50年代开始,在大气边界层湍流中也发现了相干结构——对流云街,并进行了系统的研究。近些年来,人们发现在近地层湍流中也存在相干结构。利用北京325 m气象塔对城市下垫面中大风和小风天气的风速分析,发现较有规律的周期3~6 min的阵风,且有明显的相干结构,而对不同下垫面的阵风研究,均发现存在这种相干结构,这种阵风相干结构对通量输送有不可忽视的作用。本文利用2012年4月甘肃省民勤县巴丹吉林沙漠观测塔的超声风速和平均场风速、温度观测资料,对阵风相干结构的产生条件进行了分析。采用傅立叶变换,将三维超声风速按频率分成基流(周期10分钟以上)、阵风扰动(周期1到10分钟)、湍流脉动(周期小于1分钟)三部分,结合平均场的资料分析发现:阵风相干结构出现在静力中性、不稳定甚至略微稳定的条件下,或者说机械作用主导的大气边界层,阵风区就会出现相干结构,热力作用对其有抑制和干扰的作用。从而,阵风的相干结构和壁面相干结构都出现在中性条件下,是机械湍流的现象,都主导着动量能量的输运。阵风区的相干结构并不等同于对流云街,他们出现在不同的大气稳定度条件下且尺度不同。     

用小波系数谱方法分析湍流湿度脉动的相干结构

小波系数谱分析方法是结合小波分析和高分辨率谱分析的一种统计方法,可以用来同时识别时间序列中相干结构的生命尺度和出现周期,可以很好地描述相干结构的演变过程。基于此方法,作者分析了2004年11月在河北省白洋淀地区的陆地和岛上两个观测点（分别代表陆地和水面两种不同下垫面）湍流湿度脉动的相干结构特征,结果表明陆地和水上湿度序列的相干结构尺度分布相似,并且尺度与周期之间的关系一致:小于5 s的相干结构不连续出现,而且通常伴有更大尺度的相干结构,而5~30s的相干结构有与其尺度差不多的周期。在寻找更大尺度的相干结构时发现存在一个尺度,当大于某个周期时,在各个周期上这个尺度的相干结构都显著;与正交小波变换识别相干结构主尺度的方法识别的相干结构主尺度一致。另外,小尺度结构不连续出现也可以解释小尺度湍流能量变化比较大。     

大气边界层湍流涡旋结构的小波分解

利用离散正交小波在若干物理判别准则下对边界层湍流脉动信号进行去噪和多尺度分解,从而有效地区分出均匀各向同性小涡成份和大尺度含能涡旋成份。能谱分析发现,小涡的谱动力学行为具有非常好的标度不变性,标度关系满足Kolmogorov的“-5／3”律。       

用连续子波变换提取城市冠层大气湍流的相干结构

切变湍流的相干结构是湍流研究中的重大发现,它表明湍流在表面上看来不规则运动中具有可检测的有序运动,这种相干结构在切变湍流的脉动生成和发展中起着主宰作用.因此识别和提取相干结构对于认识和研究湍流是非常重要的.用数字滤波法将包含相干结构的大尺度信号提取出来以后,再用子波分析,根据子波能量极大值的判别方法,分别确定出大气湍流三个方向上的速度脉动信号相干结构的频率或时间尺度,然后由确定尺度上的连续子波反演公式,提取出大气湍流三个方向上的速度脉动信号相干结构所对应的波形.     

大气边界层湍流多尺度分形特征的研究

运用离散正交小波变换将湍流信号分解为不同尺度，计算其分数维。考察其分数维的变化得出：随着分解层次得增加，提取湍流信号得低频部分趋于简单光滑，分数维不断减小，高频部分呈现复杂，分数维趋于定值，平均为1．70左右。说明大气边界层湍流信号在某些尺度上，存在明显的自相似性特征。     

试探大气边界层湍流运动的若干问题

对湍流运动的几个基本理论问题作了一定程度的探讨，并根据湍流动能方程分析了大气边界层湍流产生和发展的能量机制，以此为判据推导了滴运动的凡主范围和物理含意较为清楚的混合长表达式。     

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

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

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

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

近地层大气湍流微结构对比分析

利用黑河资料分析了戈壁下垫面上的近地层湍流微结构,并与其他下垫面上的端流微结构进行比较。结果表明：黑河戈壁滩上的湍流微结构满足局地相似理论,且有不同于其他下垫面的湍流结构特征;横向湍强较强,大小与纵向湍强相当。对湍流进行多尺度分析发现,黑河地区湍流结构呈现明显的多尺度分布,湍涡尺度量级从１０＾－２到１０＾２不等; 相比于ｖ方向,ｕ方向上占主要地位的湍涡尺度较大,但是所占比例罗小;随着不稳定度增加,     

Long-term study of coherent structures in the atmospheric surface layer

A long-term study of coherent turbulence structures in the atmospheric surface layer has been carried out using 10 months of turbulence data taken on a 30-m tower under varying meteorological conditions. We use an objective detection technique based on wavelet transforms. The applied technique permits the isolation of the coherent structures from small-scale background fluctuations which is necessary for the development of dynamical models describing the evolution and properties of these phenomena. It was observed that coherent structures occupied 36% of the total time with mean turbulent flux contributions of 44% for momentum and 48% for heat. The calculation of a transport efficiency parameter indicates that coherent structures transport heat more efficiently than momentum. Furthermore, the transport efficiency increases with increasing contribution of the structures to the overall transport.     

Analysis of turbulent exchange and coherent structures in the stable atmospheric boundary layer based on tower observations

Using data collected at the Spanish low troposphere research centre CIBA (Centro de Investigación de la Baja Atmósfera) and at the Cabauw Experimental Site for Atmospheric Research (CESAR) in the Netherlands, we analyzed the most significant features of different coherent structures occurring in the stable atmospheric boundary layer. In particular, we used both the Reynolds and wavelet methods to analyze a solitary wave, a gravity wave, a density current and a low-level jet. For each of these structures, we found that wavelet analysis had the capacity to distinguish the different scales involved in these events due to the different timing and heights of the thermal instabilities and downdrafts associated with the disturbances. In addition, the wavelet method highlights the different roles of turbulence and coherent structures in the transfer of heat, moisture and CO₂ in the nocturnal boundary layer.     

The Theoretical Model of Atmospheric Turbulence Spectrum in Surface Layer

It is shown that the slope of energy spectrum obtained from the velocity solution of Kdv-Burgers equation lies between -5/3 and -2 in the dilogarithmic coordinates paper. The spectrum is very close to one of Kolmogorov’s isotropic turbulence and Frisch’s intermittent turbulence in inertial region. In this paper, the Kdv-Burgers equation to describe atmospheric boundary layer turbulence is obtained. In the equation, the 1 / Re, corresponds to dissipative coefficient v, to dispersive coefficient β, then (v/ 2β)2 corresponds to .We prove that the wave number corresponding to maximum energy spectrum decreases with the decrease of stability (i.e., the increase of in eddy-containing region. And the spectrim amplitude decreases with the increase of (i.e., the decrease of stability). These results are consistent with actual turbulence spectrum of atmospheric surface layer from turbulence data.     

Characteristics of atmospheric turbulence in the surface layer over Antarctica

This paper presents meteorological measurements made during the antarctic summer period, on two 9 m and 3 m towers, on the rocky and ice shelf terrains of the Indian antarctic stations Maitri and Dakshin Gangotri, respectively. The measurements of fluctuations in temperature and wind speed made with relatively lesser precision instrumentation pertain to smaller wave numbers ~10^-2 m^-1 appropriate to outer scale L ₀ of the atmospheric turbulence spectrum. Autocorrelation analysis of the fluctuations in temperature and wind speed has been performed. A new autoregressive scheme has been developed to represent the computed autocorrelation functions by a Yule statistical model, and to estimate the correlation period T ₀ of the turbulent medium. Height profiles of outer scale L ₀ of turbulence may be given in terms of T ₀ and mean wind speed u. Further, the similarity theory of Monin-Obukhov has been used to compute height profiles of temperature structure parameter C _T ². At Maitri, values of L ₀ and C _T ² are higher between 03–22 h local time than between 22–03 h. Values of L ₀ and C _T ² are smaller over the ice shelf terrain of the Dakshin Gangotri station, compared to those over the rocky terrain of the Maitri station.     

A model of turbulence spectra in the atmospheric surface layer

The spectral equations of turbulent kinetic energy and temperature variance have been solved by using Onsager's energy cascade model and by extending Onsager's model to closure of terms that embody the interaction of turbulent and mean flow.The spectral model yields the following results: In a stably stratified shear flow, the peak wave numbers of the spectra of energy and temperature variance shift toward larger wave numbers as stability increases. In an unstably stratified flow, the peak wave numbers of energy spectra move toward smaller wave numbers as instability increases, whereas the opposite trend is observed for the peak wave numbers of temperature variance spectra. Hence, the peak wave numbers of temperature spectra show a discontinuity at the transition from stable to unstable stratification. At near neutral stratification, both spectra reveal a bimodal structure.The universal functions of the Monin-Obukhov similarity theory are predicted to behave as _m ~ _H ~ (- Z/L)^-1/3 in an extremely unstable stratification and as _m ~ _H ~ z/L in an extremely stable stratification. For a stably stratified flow, a constant turbulent Prandtl number is expected.     

Characteristics of turbulent structures in the unstable atmospheric surface layer

An atmospheric surface-layer (ASL) experiment conducted at a meteorological site in the Oostelijk-Flevoland polder of the Netherlands is described. Turbulent fluctuations of wind velocity, air temperature and static pressure were measured, using three 10 m towers.Simultaneous turbulent signals at several heights on the towers were used to investigate the properties of the turbulent structures which contribute most significantly to the turbulent vertical transports in the unstable ASL. These turbulent structures produce between 30 and 50% of the mean turbulent vertical transport of horizontal alongwind momentum and they contribute to between 40 and 50% of the mean turbulent vertical heat transport; in both cases this occurs during 15 to 20% of the total observation time.The translation speed of the turbulent structures equals the wind speed averaged over the depth of the ASL, which scales on the surface friction velocity. The inclination angle of the temperature interface at the upstream edge of the turbulent structures to the surface is significantly smaller than that of the internal shear layer, which is associated with the temperature interface. The turbulent structures in the unstable ASL are determined by a large-scale temperature field: Convective motions, which encompass the whole depth of the planetary boundary layer (PBL), penetrate into the ASL. The curvature of the vertical profile of mean horizontal alongwind velocity forces the alignment of the convective cells in the flow direction (Kuettner, 1971), which have an average length of several hundreds of metres and an average width of a few tens of metres. This mechanism leads to the formation of turbulent structures, which extend throughout the depth of the ASL.     

Resolved-scale turbulence in the atmospheric surface layer from a large eddy simulation

Large eddy simulation has encountered difficulties in handling turbulence in the atmospheric surface layer due to deficiencies in sub-grid scale models. This paper addresses the possibility of resolving the turbulence in the upper part of the surface layer by a low-aspect ratio of grid spacing. Results show that resolved-scale shear stresses dominate over the sub-grid scale components so that effects due to the sub-grid scale model can be ignored in this region. The effects of the lower boundary condition on the resolved-scale turbulence in the upper part of the surface layer are discussed. It is concluded that the normalized mean velocity shear and resolved turbulence in the upper part of the surface layer are not affected by the specification of the lower boundary condition. In addition, the present work proposes a new independent model parameter, the Smagorinsky Reynolds Number (Re_SM), and demonstrates that this number determines the resolved turbulence in the upper part of the surface layer.     

0814号强台风“黑格比”风场相干结构统计和演化特征研究

为了探究台风风场相干结构统计特性以及在台风不同空间结构部位的演化规律,采用小波系数谱分析(Spavelet Analysis)方法对0814号强台风“黑格比”风速场的湍流相干结构统计特征及演化特性进行了系统分析。研究结果表明：对于0814号强台风“黑格比”过程,顺风向脉动风速与横风向脉动风速的相干结构呈现出多尺度分布特征;顺风向、横风向和竖向风速分量的相干结构主周期均值分别为70.2 s、50.3 s、22.8 s,主尺度均值分别为54.4 s、38.9 s、18.2 s,三个方向相干结构的周期尺度比均值分别为1.28、1.27、1.29;台风“黑格比”风场水平向空间部位中,前外环流区(FOV)的顺风向、横风向和竖向的相干结构主尺度均值分别为68.8 s、41.7 s、14.4 s;前眼壁区(FEW)的三个方向脉动风速的相干结构主尺度均值分别为69.9 s、32.5 s、19.5 s;风眼区三个方向脉动风速的相干结构主尺度分别为83.5 s、50.3 s、15.3 s;后眼壁区(BEW)其值为55.1 s、42.3 s、11.2 s;后外环流区(BOV)三个方向脉动风速相干结构的主尺度分别为53.6 s、39.1 s、10.9 s;整体上表现为前外环流区(FOV)和前眼壁区(FEW)相干结构主尺度大于后眼壁区(BEW)和后外环流区(BOV),风眼区相干结构主尺度则为最大。     

TIPEX高原东南部近地层湍流参数的分析

The characteristics of the turbulence spectrum, turbulence variance, and turbulence flux at Qamdo overthe Southeastern Tibetan Plateau measured during the TIPEX experiment from 18 May to 30 June 1998are analyzed by the eddy correlation method. During periods of intense convection, most of the spectraof the 3D winds, temperature, and humidity follow a power law with a slope of -2/3. The normalizedvariances of the 3D winds in relation to z/L satisfy the similarity law but the normalized variances oftemperature and humidity are related to z/L by a -1/3 power law only in unstable conditions. In near-neutral stratification, σu/u* and σy/u* at Qamdo and Gerze are nearly constant in rugged terrain andσw/u* at Qamdo and Gerze is similar to the value found in plains, which indicates that the, effects oflandform on vertical turbulence is not significant over the southeast and western Tibetan Plateau. Duringthe dry period, the sensible heat dominates, comprising 81% of the heat intensity, with the other 19%being latent heat; during the wet period on the other hand, the latent heat accounts for about 64% ofthe available energy and the sensible heat only around 36%. The maximum intensity of the heating of airoccurs in the middle of the Plateau during the summertime.