The partitioning of attached and detached eddy motion in the atmospheric surface layer using Lorentz wavelet filtering

Townsend's attached eddy hypothesis states that the turbulent structure in the constant stress layer can be decomposed into attached and detached eddy motion. This paper proposes and tests a methodology for separating the attached and detached eddy motion from time series measurements of velocity and temperature. The proposed methodology is based on the time-frequency localization and filtering capabilities of the orthonormal wavelet transforms. Using a relative entropy statistical measure, the optimal wavelet basis is identified first. The turbulence time series measurements are then transformed into the wavelet domain where the contribution of specific events in the time-frequency domain is identified. The filtering scheme utilizes a recently constructed Lorentz thresholding methodology that successfully eliminates all wavelet coefficients associated with the detached eddy motion. While this filtering scheme lacks the compression efficiency of the classical Donoho and Johnstone's universal thresholding model, it conserves the higher-order statistics and important turbulence interactions related to the Reynolds stresses. Following the filtering scheme, the attached eddy motion time series is re-constructed by an inverse wavelet transform of the non-zero wavelet coefficients. The proposed partitioning methodology for attached and detached eddy motion is tested using 56 Hz triaxial sonic anemometer velocity and temperature measurements above a uniform dry lake bed in Owens valley, California, for a wide range of atmospheric stability conditions. Validation that the wavelet filtered time series represents the attached eddy motion is also discussed in the context of conservation of turbulence energy and surface fluxes.     

On the downward transport of turbulent kinetic energy in the surface layer over plant canopies

Downward fluxes of turbulent kinetic energy have been frequently observed in the air layer just above plant canopies. In order to investigate the mechanism for such downward transport, analysis of observational data is attempted. Height-dependency of turbulent kinetic energy flux and turbulence statistics including higher order moments is represented as a function of a non-dimensional height z/H, where z is an observational height and H an average height of plant canopies. Downward fluxes and non-Gaussianity of wind velocity fluctuations are predominant just above plant canopies and decrease with increasing height. The downward flux is closely related to the high intensity of turbulence and the non-Gaussianity of wind velocity fluctuations, especially with a positive skewness in the longitudinal wind and a negative skewness in the vertical wind. The analysis method of conditional sampling and averaging is applied to the present observations. The results show that the predominance of the intermittent inrush phase over the intermittent ejection phase leads to the above-mentioned non-Gaussianity. Finally, a simple explanation is given in order to interpret the turbulent flow structure in the air layer near the plant canopies, which is associated with the downward energy transport process.     

郑州近54年降水变化的多时间尺度分析

利用Morlet小波分析方法,对1951~2004年郑州市年降水距平时间序列的小波变化特征进行分析,揭示降水变化的多时间尺度结构,分析其中存在的主要周期振荡和突变点,并与功率谱分析结果进行比较。结果表明:郑州市年降水存在多时间尺度的周期变化特征,郑州市年平均降水存在3 a、准5 a、准9 a和准21 a的主周期。郑州市降水变化存在着明显的多时间尺度,即年代际尺度和年际尺度的周期性变化。郑州市降水变化还显示出自己的突变点分布及其位相结构。在降水量总趋势下降的前提下,2004年后的一段时间还将处在一个偏多期。     

Turbulent transport processes of momentum and sensible heat in the surface layer over a paddy field

A sonic anemometer-thermometer was used to measure turbulent fluxes of momentum and sensible heat and related turbulence statistics just above plant canopies in unstable conditions. The stability dependence of transport processes is presented, using observational data in a wide range of instability. The analysis of joint probability distributions of w – u, w – T, w – wu, and w – wT confirmed that just above plant canopies, downdrafts were remarkably efficient for vertical transport of momentum and scalar quantities in near neutral conditions. Furthermore, it was shown that updrafts became more efficient than downdrafts for vertical transport, especially of scalar quantities, in very unstable conditions.     

Detection and analysis of coherent structures in urban turbulence

Summary The continuous wavelet transform provides a suitable tool to visualize the vertical structure of turbulence and to detect coherent structures in turbulent time series. This is demonstrated with a simple example of an artificially ramp structured time series. In this study turbulence data, i.e. the fluctuations of the horizontal wind components u′ and v′, the vertical component w′ and temperature T′, sampled with 20.83 Hz and measured simultaneously at three levels (z/h=1.5, 2.1 and 3.2, with z as the sensor height and h the height of the roughness elements) over an urban canopy in the inner city of Basel, Switzerland, are analyzed. The detection of the coherent structures was performed using the Mexican hat wavelet and the zero-crossing method. The analysis for unstable conditions shows that organized structures (ejection-sweep cycles) cover about 45% of the total run time. A conditional average from a total of 116 detected ejection-sweep sequences during 7 hours was calculated over a time window of 100 seconds. This dominating time scale was derived from peak frequencies of the wavelet spectra as well as from the Fourier spectra. It is shown that the normalized amplitudes of fluctuations of temperature and longitudinal wind speed during the events are largest at the lowest measurement level just above the canopy and decrease with increasing distance from the roughness elements. A comparison of related studies over different non-urban surfaces (mainly forests) shows that the shape of conditionally averaged ejection-sweep sequences is very similar for all canopies, however, the dominating time scale in general increases the rougher the surface is and the higher the roughness elements are.     

NUMERICAL SIMULATIONS OF TURBULENT FLOW IN AND ABOVE PLANT CANOPIES

The turbulent flow in and above plant canopies is of fundamental importance to the understanding oftransport processes of momentum,heat and mass between plant canopies and atmosphere,and to microme-teorology.The Reynolds stress equation model(RSM)has been applied to calculate the turbulence in cano-pies in this paper.The calculated mean wind velocity profiles,Reynolds stress,turbulent kinetic energy andviscous dissipation rate in a corn canopy and a spruce forest are compared with field observed data and withWilson's and Shaw's model.The velocity profiles and Rynolds stress calculated by both models are in goodagreement,and the length scale of turbulence appears to be similar.     

植被内部及其上方湍流场的数值模拟

植被内部及其上方的湍流流场对于了解植被与大气之间的动量、热量和质量交换过程极其重要。本文把高阶湍流封闭模型的Reynolds应力方程模型(RSM)应用于植被湍流的计算,得到了风速、湍流动能、Reynolds应力及能量耗散率的垂直分布,与现场观测数据比较,甚为满意。     

我国华东地区气候变化与副高指数的小波分析

根据1951～2007年华东地区逐月平均气温和降水资料,在利用连续小波变换对华东地区的降水、气温与副高指数的时频变化特征进行分析的基础上,采用交叉小波变换方法分析了华东地区近57年的气候变化与副高之间的相关关系.结果表明,华东地区气候要素与副高指数距平序列都存在着不同时间尺度的周期分量,并且各周期分量的强度是不同的.其中以2～8年尺度的年际振荡较为明显,10年以上的年代际变化以1980年为突变点.     

Large-Scale Turbulence Structures and Their Contributions to the Momentum Flux and Turbulence in the Near-Neutral Atmospheric Boundary Layer Observed from a 213-m Tall Meteorological Tower

Large-scale turbulence structures in the near-neutral atmospheric boundary layer (ABL) are investigated on the basis of observations made from the 213-m tall meteorological tower at Tsukuba, Japan. Vertical profiles of wind speed and turbulent fluxes in the ABL were obtained with sonic anemometer-thermometers at six levels of the tower. From the archived data, 31 near-neutral cases are selected for the analysis of turbulence structures. For the typical case, event detection by the integral wavelet transform with a large time scale (180 s) from the streamwise velocity component (u) at the highest level (200 m) reveals a descending high-speed structure with a time scale of approximately 100 s (a spatial scale of 1 km at the 200-m height). By applying the wavelet transform to the u velocity component at each level, the intermittent appearance of large-scale high-speed structures extending also in the vertical is detected. These structures usually make a large contribution to the downward momentum transfer and induce the enhancement of turbulent kinetic energy. This behaviour is like that of “active” turbulent motions. From the analysis of the two-point space–time correlation of wavelet coefficients for the u velocity component, the vertical extent and the downward influence of large-scale structures are examined. Large fluctuations in the large-scale range (wavelet variance at the selected time scale) at the 200-m level tend to induce the large correlation between the higher and lower levels.     

Observations of Coherent Turbulence Structures in the Near-Neutral Atmospheric Boundary Layer

Turbulence structures of high Reynolds number flow in the near-neutral atmospheric boundary layer (ABL) are investigated based on observations at Shionomisaki and Shigaraki, Japan. A Doppler sodar measured the vertical profiles of winds in the ABL. Using the integral wavelet transform for the time series of surface wind data, the pattern of a descending high-speed structure with large vertical extent (from the surface to more than 200-m level) is depicted from the Doppler sodar data. Essentially this structure is a specific type of coherent structure that has been previously shown in experiments on turbulent boundary-layer flows. Large-scale high-speed structures in the ABL are extracted using a long time scale (240 s) for the wavelet transform. The non-dimensional interval of time between structures is evaluated as 3.0–6.2 in most cases. These structures make a large contribution to downward momentum transfer in the surface layer. Quadrant analyses of the turbulent motion measured by the sonic anemometer (20-m height) suggest that the sweep motion (high-speed downward motion) plays a substantial role in the downward momentum transfer. In general, the contribution of sweep motions to the momentum flux is nearly equal to that of ejection motions (low-speed upward motions). This contribution of sweep motions is related to the large-scale high-speed structures.     

Turbulent transport of scalar quantities within and above a paddy field

An experiment was conducted to study turbulent transport processes of scalar quantities within and above a rice plant canopy. A sonic anemometer-thermometer and a Lyman- humidiometer were used to measure the turbulent fluxes of sensible and latent heat and related turbulence statistics within a paddy field. The sensible and latent heat fluxes measured at two heights within and above the plant canopy showed that the upper layer of this plant canopy was an active source region and that the source strength of sensible and latent heat depended on the solar radiation and physiology of rice plants. Analysis of joint probability distributions of w and T and of w and q within this plant canopy showed that downdrafts were remarkably efficient for upward transport of sensible and latent heat in the daytime. The vertical fluxes of temperature and humidity variance were also divergent from the upper layer of plant canopies. The power spectra of temperature and humidity within the plant canopy decreased rapidly in the high frequency range, compared with the - 2/3 law relationship of nS(n) vs log n observed above plant canopies.     

北极涛动对河南省降水的可能影响

采用连续小波变换和交叉小波变换方法,分析了北极涛动指数和河南省月降水距平序列的时频变化特征,并讨论了近50多年来北极涛动对河南省降水的可能影响。结果表明,北极涛动指数和河南省降水距平序列中都存在着不同时间尺度的周期振荡,且时域中的分布特征和振荡强度也不尽相同;两者之间在6~8 a、12~16 a和25 a以上尺度的周期振荡呈正相关,而在9~11 a和16~24 a尺度上表现为负相关关系。     

正交小波变换研究复杂下垫面边界层的湍流特征

小波变换方法具有较好的时频局部特性,非常适合于分析非平稳的湍流信号。本文对35 m铁塔的超声风速测量数据进行了离散正交小波变换,计算了各向同性系数(ISO isotropy coefficient)及小波功率谱,以此对实验场所代表的水陆交际复杂下垫面的近地面层湍流特征进行了研究。ISO系数可以很好地描述实际大气在不同尺度的各向同性特征。根据ISO系数的分布,我们可以通过设定一定的阈值(ISO=0.7)来得到各向同性小涡的分离尺度(记为ISO0.7),即湍流的各向同性尺度范围。研究表明,下垫面对于边界层湍流各向同性特性具有相当大的影响,当风从不同下垫面吹过时ISO0.7尺度的均值具有较明显的差异。同时湍流度和稳定度对湍流各向同性也有一定的影响。对湍流小波功率谱研究表明,功率谱谱幂率在分离出的小涡所处的频段接近于-5/3。而当风从陆面吹过且湍流度较弱时,湍流小波功率谱小涡以分离尺度所对应的频率(fms)为界具有明显的两段趋势:高于fms的频段,谱幂率一般接近-5/3;而低于fms的频段的谱幂率接近-3/3。这反映了下垫面特征对湍流功率谱的影响。同时在不同的风向转变过程中下垫面对功率谱的影响具有差异。     

Structure Function Analysis of Two-Scale Scalar Ramps. Part II: Ramp Characteristics and Surface Renewal Flux Estimation

Ramp features in the turbulent scalar field are associated with turbulent coherent structures, which dominate energy and mass fluxes in the atmospheric surface layer. Although finer scale ramp-like shapes embedded within larger scale ramp-like shapes can readily be perceived in turbulent scalar traces, their presence has largely been overlooked in the literature. We demonstrate the signature of more than one ramp scale in structure functions of the turbulent scalar field measured from above bare ground and two types of short plant canopies, using structure-function time lags ranging in scale from isotropic to larger than the characteristic coherent structures. Spectral analysis of structure functions was used to characterize different scales of turbulent structures. By expanding structure function analysis to include two ramp scales, we characterized the intermittency, duration, and surface renewal flux contribution of the smallest (i.e., Scale One) and the dominant (i.e., Scale Two) coherent structure scales. The frequencies of the coherent structure scales increase with mean wind shear, implying that both Scale One and Scale Two are shear-driven. The embedded Scale One turbulent structure scale is ineffectual in the surface-layer energy and mass transport process. The new method reported here for obtaining surface renewal-based scalar exchange works well over bare ground and short canopies under unstable conditions, effectively eliminating the α calibration for these conditions and forming the foundation for analysis over taller and more complex surfaces.     

Development and Validation of a Lagrangian Probability Density Function Model of Horizontally-Homogeneous Turbulence Within and Above Plant Canopies

Second- and third-order turbulence closure models have met with mixed success when applied to the prediction of turbulent flows within and above plant canopies and model predictions are typically no better than those obtained using simpler two-equation (k-e){(k-\varepsilon)} models. This is because the local gradient diffusion approximation (a crucial model requirement) cannot represent accurately turbulent transport that is dominated by the presence of ejections and sweeps whose length scales are comparable with the canopy height. To make progress, turbulent transport must be treated without approximation, as in Lagrangian probability density function (PDF) models. This study is the first to develop and to validate a PDF model of horizontally-homogeneous canopy flow. The model relies upon a prescribed length-scale that has been used elsewhere in the modelling of turbulent flows. Model predictions compare favourably with measurements of neutrally stratified turbulent flows within and above canopies of mature corn and forested eucalypt trees.     

Analysis of oscillations in the stable atmospheric boundary layer using wavelet methods

A wavelet method is used to estimate kinetic energy and fluxes from data collected under stable conditions during the CASES-99 field campaign. Results in the high frequency range are compared with those obtained by the traditional method used to estimate turbulent moments, which is based on the Reynolds decomposition of variables into a mean and a turbulent part. The fact that the wavelet transform performs much better as a filter than the averaging process accounts for most of the disagreements between results. Since the wavelet method can be applied at very different spectral ranges, it is also used to analyse two different coherent structures: a density current and a train of internal gravity waves. The strong burst of turbulence related to the density current reflects the complexity of the first event. The wavelet method discriminates the different scales of motion, which are present in the perturbation, and is therefore an ideal tool for assessing the interactions between them. A method based on the phase difference between wavelet-transformed time series is then applied to the analysis of the horizontal and vertical structure of the gravity waves, and a three-dimensional image of the oscillations is provided.     

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

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

Energy balance and turbulent flux partitioning in a corn–soybean rotation in the Midwestern US

Quantifying the energy balance above plant canopies is critical for better understanding of water balance and changes in regional weather patterns. This study examined temporal variations of energy balance terms for contrasting canopies [corn (Zea mays L.) and soybean (Glycine max L. Merr.)]. We monitored energy balance for 4 years using eddy-covariance systems, net radiometers, and soil heat flux plates in adjacent production fields near Ames, Iowa. On an annual basis, soybean exhibited 20% and 30% lower sensible heat flux (H) and Bowen ratio than corn, respectively. As canopies developed, a gradual shift in turbulent fluxes occurred with decreasing H and increasing latent heat flux (LE), but with a more pronounced effect for corn. Conversely, during mid-growing season and as both canopies progressively senesced, H in general increased and LE decreased; however, soybean exhibited slightly greater LE and much lower H than corn. These temporal variations in magnitude and partitioning of turbulent fluxes translated into a pronounced energy imbalance for soybean (0.80) and an enhanced closure for corn (0.98) in August and September. These discrepancies could be directly associated with differences in momentum transport as shown by friction velocities of 0.34 and 0.28 m s^?1 for corn and soybean, respectively. These results support influential roles of plant canopy on intensity and mode of surface energy exchange processes.     

Hilbert方法在大气边界层湍流特征分析中的有效性

介绍一种新的建立在经验模态分解(EMD)方法基础上的非线性、非平稳数据分析技术一Hilbert分析技术,并首次将其应用于大气边界层(PBL)湍流数据的分析,初步探讨了其在PBL湍流研究中的有效性.通过对城市与森林冠层上湍流资料的能量分布特征和统计平稳度进行分析、比较,结果表明:Hilbert谱分析能有效地对PBL湍流信号进行分析.它的边缘谱分析能够有效地探测PBL湍流信号的能量分布特征,统计平稳度分析也能有效地给出PBL湍流信号平稳性的定量化测量,这些将有助于建立合适的数据质量控制方法,以及对现有空气质量与扩散模式中扩散参数的计算加以改进.文中个例分析中,城市和森林冠层上空的湍流有一定相似性,湍流混合都比较充分,但森林冠层上湍流信号的能量更多地集中在大尺度湍涡,且扰动风速的高频部分具有更强的间歇性.对于相近高度的湍流信号来说,多数情况下,森林冠层上相同尺度的湍涡表现得比城市冠层上更不稳定,但湍涡的含能量要更低.