Conditional Wavelet Technique Applied to Aircraft Data Measured in the Thermal Internal Boundary Layer During Sea-Breeze Events

We describe a wavelet-based technique to determine the spectral turbulentcontribution to the vertical flux of sensible heat in a position-wavelength representation. This technique combines awavelet transform (Morlet wavelet) with conditional sampling. We apply this methodto aircraft datacollected during a sea-breeze circulation (BEMA97 experiment) with heterogeneousturbulence conditions horizontally and vertically as well. The turbulent fluxes are analysed with the conditional wavelet techniqueas a function of the wavelength and the horizontal distance.The turbulent processes within the thermal internal boundary layer associatedwith the sea breeze are clearly identified. The results exhibit the wavelength bands through which the upward flux (originating fromthe surface) and the downward flux (originating from the top of theboundary layer) are important.     

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

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

The local atmosphere and the turbulent heat transfer in the Eastern Himalayas

To understand the local atmosphere and heat transfer and to facilitate the boundary-layer parameterization of numerical simulation and prediction,an observational campaign was conducted in the Eastern Himalayas in June 2010.The local atmospheric properties and near-surface turbulent heat transfers were analyzed.The local atmosphere in this region is warmer,more humid and less windy,with weaker solar radiation and surface radiate heating than in the Middle Himalayas.The near-surface turbulent heat transfer in the Eastern Himalayas is weaker than that in the Middle Himalayas.The total heat transfer is mainly contributed by the latent heat transfer with a Bowen ratio of 0.36,which is essentially different from that in the Middle Himalayas and the other Tibetan regions.     

Experimental investigation of atmospheric boundary layer turbulence

This paper describes how to measure turbulence in the atmospheric boundary layer (ABL) in order to address certain problems in modern atmospheric physics. These problems mainly relate to the Earth's energy budget (including the hydrological cycle) and biogeochemical cycles. Starting from the main characteristic numbers and the basic equations of atmospheric turbulent flow, we show what turbulence parameters are important to measure. Special attention is given to the various methods used to compute the turbulent fluxes. We analyse the range of scales which has to be measured to properly capture the eddies contributing to the turbulent transfers. This range of scales determines what sensors can be used in the atmospheric surface layer and in the ABL. We describe the most widely used instruments and their performances. The principal platforms used to deploy these instruments are examined. Aircraft are described in more details, because they allow a thorough exploration of the ABL. In the last section, some examples of ABL turbulence signals measured in various conditions are presented. These examples illustrate horizontally homogeneous turbulence as well as inhomogeneous signals for which standard analysis techniques cannot be used. We show how some recent techniques, like wavelet transforms, can help to investigate this kind of signal. At the end, we present what would be interesting to do in the near future for the study of ABL turbulence.     

Conditional Analysis of an Internal Boundary Layer

Data from a convective internal boundary layer (IBL) are analyzed by focusing on the instantaneousstructure of the top of the IBL instead of the time-average structure.A conditional averaging technique is developed todiscriminate between air from above the IBL and air from below the IBL , which alternately invade some instrument levels due tosubstantial variation of the top of the convective IBL.Sensitivity to the conditional sampling criteria is examined.Inside the IBL , buoyant and mechanicalproduction and dissipation dominate the turbulent kinetic energy budget.The horizontal advection and turbulent transport terms are smaller, but not negligible. The inferred pressure correlation term is negligible.Above the IBL , buoyant production and dissipation, although weak,dominate the turbulent kinetic energy budget. Shear generation andturbulent transport are smaller but significant.     

Surface-layer intermittency investigated with conditional sampling

A conditional sampling technique is used to provide statistics of surface-layer plume properties. A selection criterion based on the high-frequency variance of the horizontal wind component enables an accurate division of plume and nonplume states. The intermittency factor derived with this technique closely matches values obtained using other techniques at various heights in the atmospheric boundary layer. The intermittency factor in addition to other plume statistics are found to be stability dependent. Conditional averages are used to produce scatter diagrams from which the interrelationships between properties of both the plume and nonplume states can be examined. Several provocative relationships discovered in this way are discussed.An extensive investigation into the bimodal nature of the fine structure of turbulence is described. These results provide the most compelling support for the division of surface-layer turbulence into separate states. Length scales derived from the second moments of distributions fitted to conditionally sampled data are found to correlate with external parameters of the flow.Department of Atmospheric Sciences contribution number 514.     

Use of obliquely rotated principal component analysis to identify coherent structures

Principal component analysis (PCA) with oblique rotation is applied to Large Eddy Simulation (LES) results to discern and quantify coherent structures within the convective boundary layer (CBL). Sensitivity tests are first conducted on a moderately convective LES run. Once the ability of PCA to generate robust results is verified, the method is applied to LES runs spanning a range of stability regimes. Interregime similarities and differences in the coherent structures are discussed. For the moderately convective LES run, three-dimensional convective cells are arrayed in two-dimensional bands aligned with the geostrophic wind. The resulting gravity waves in the free atmosphere and convective inflow and outflow in the boundary layer are also captured by the PCA. Convective modes are more sensitive to the ratio of w _* to u _* than are the dynamic modes.PCA has demonstrated advantages over previous analysis methods. PCA score maps provide information on the spatial distribution of phenomena that has not been available from traditional conditional sampling studies. Principal components provide information on the vertical structures of phenomena that would be obscured by life-cycle effects or erratic tilts from the vertical in the conventional approaches to either conditional sampling or composite analysis. Future work includes application of this technique to multi-level observational time series from a surface-layer tower for the Risø Air/Sea Experiment (RASEX).     

Heterogeneity-Induced Heat-Flux Patterns in the Convective Boundary Layer: Can they be Detected from Observations and is There a Blending Height?—A Large-Eddy Simulation Study for the LITFASS-2003 Experiment

An understanding of how the convective boundary layer (CBL) is mixed under heterogeneous surface forcing is crucial for the interpretation of area-averaged turbulence measurements. To determine the height and degree to which a complex heterogeneous surface affects the CBL, large-eddy simulations (LES) for two days of the LITFASS-2003 experiment representing two different wind regimes were undertaken. Spatially-lagged correlation analysis revealed the turbulent heat fluxes to be dependent on the prescribed surface flux pattern throughout the entire CBL including the entrainment layer. These findings prompted the question of whether signals induced by surface heterogeneity can be measured by airborne systems. To examine this question, an ensemble of virtual flights was conducted using LES, according to Helipod flight measurements made during LITFASS-2003. The resulting ensemble-averaged heat fluxes indicated a clear dependence on the underlying surface up to the top of the CBL. However, a large scatter between the flux measurements in different ensemble runs was observed, which was the result of insufficient sampling of the largest turbulent eddies. The random and systematic errors based on the integral length scale did not indicate such a large scatter. For the given flight leg lengths, at least 10–15 statistically independent flight measurements were necessary to give a significant estimate of heterogeneity-induced signals in the CBL. The need for ensemble averaging suggests that the observed blending of heterogeneity-induced signals in the CBL can be partly attributed to insufficient averaging.     

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.     

High-Resolution Fibre-Optic Temperature Sensing: A New Tool to Study the Two-Dimensional Structure of Atmospheric Surface-Layer Flow

We present a novel approach based on fibre-optic distributed temperature sensing (DTS) to measure the two-dimensional thermal structure of the surface layer at high resolution (0.25 m, ≈0.5 Hz). Air temperature observations obtained from a vertically-oriented fibre-optics array of approximate dimensions 8 m × 8 m and sonic anemometer data from two levels were collected over a short grass field located in the flat bottom of a wide valley with moderate surface heterogeneity. The objectives of the study were to evaluate the potential of the DTS technique to study small-scale processes in the surface layer over a wide range of atmospheric stability, and to analyze the space–time dynamics of transient cold-air pools in the calm boundary layer. The time response and precision of the fibre-based temperatures were adequate to resolve individual sub-metre sized turbulent and non-turbulent structures, of time scales of seconds, in the convective, neutral, and stable surface layer. Meaningful sensible heat fluxes were computed using the eddy-covariance technique when combined with vertical wind observations. We present a framework that determines the optimal environmental conditions for applying the fibre-optics technique in the surface layer and identifies areas for potentially significant improvements of the DTS performance. The top of the transient cold-air pool was highly non-stationary indicating a superposition of perturbations of different time and length scales. Vertical eddy scales in the strongly stratified transient cold-air pool derived from the DTS data agreed well with the buoyancy length scale computed using the vertical velocity variance and the Brunt–Vaisala frequency, while scales for weak stratification disagreed. The high-resolution DTS technique opens a new window into spatially sampling geophysical fluid flows including turbulent energy exchange.     

The Influence of Local Stability on Heat and Momentum Transfer within Open Canopies

Eddy-covariance data have been analyzed to investigate the influence of local stability on heat transfer within open canopies. The flux–gradient relationship for heat is derived from the temperature variance equation, and the stability dependence of the flux–gradient relationship is examined and discussed. The results indicate that the strong stability dependence of the nondimensional standard deviation of temperature, and the small contributions of turbulent transport to the temperature variance, lead to a strong stability dependence of the nondimensional temperature gradient within open canopies. Quadrant analysis and hole size analysis were performed for momentum and heat fluxes in the subcanopy, and the results indicate that the contribution of each quadrant to the total flux depends on both the local stability and canopy depth. The intermittency of the turbulent flux does not show a clear dependence on local stability. As the contribution of ejections to the heat flux increases, the vertical flux of the temperature variance changes sign from negative to positive, leading to small temperature variance transport in unstable conditions. Multi-resolution analysis indicates that heat and momentum are transported with different dominant time scales in very unstable conditions, suggesting a different role of local buoyancy in heat and momentum transfer.     

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.     

Estimating the Random Error in Eddy-Covariance Based Fluxes and Other Turbulence Statistics: The Filtering Method

A spatially local decomposition of turbulent fluxes based on properties of spatial filters is used to develop a new method of estimating random error in turbulent moments of any order. The proposed error estimation method does not require an estimate of the integral time scale, which can be highly sensitive to the method used to calculate it. The error estimation method is validated using synthetic flux data with a known ensemble mean and intercompared with existing methods using data from the Advection Horizontal Array Turbulence Study (AHATS). Typical errors for a 27.3-min block of data collected at a height of 8 m are found to be approximately 10% for the heat flux and 7–15% for variances. The error in the momentum flux increases rapidly with increasing atmospheric instability, reaching values of 40% or greater for unstable conditions. A new method based on filtering is also proposed to estimate integral time scales of turbulent quantities.     

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

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

Maximal Overlap Wavelet Statistical Analysis With Application to Atmospheric Turbulence

Statistical tools based on the maximal overlap discrete wavelet transform (MODWT) are reviewed, and then applied to a dataset of aircraft observations of the atmospheric boundary layer from the tropical eastern Pacific, which includes quasi-stationary and non-stationary segments. The wavelet methods provide decompositions of variances and covariances, e.g. fluxes, between time scales that effectively describe a broadband process like atmospheric turbulence. Easily understood statistical confidence bounds are discussed and applied to these scale decompositions, and results are compared to Fourier methods for quasi-stationary turbulence. The least asymmetric LA(8) wavelet filter yields coefficients that exhibit better uncorrelatedness across scales than the Haar filter and is better suited for decomposition of broadband turbulent signals. An application to a non-stationary segment of our dataset, namely vertical profiles of the turbulent dissipation rate, highlights the flexibility of wavelet methods.     

Observations Of The Turbulence Structure Within Two Stratocumulus-Topped, Marine Boundary Layers

Two situations observed during the second Aerosol Characterization Experiment (ACE-2) are analysed from aircraft measurements in the broken stratocumulus (Sc)-topped marine boundary layer. The first one (26 June 1997), characterized by a non-polluted, oceanic air mass, presents a decoupling between the Sc layer (1400–1520 m) and the turbulent mixed layer, this latter extending from the surface up to 580 m. In contrast, the second case (9 July 1997), during which continental air had been advected over the experimental area, presents a well-coupled layer extending from the surface up to the top of the Sc layer(910 m). This coupling, uncommon in this area in the middle of the day, isrelated to the relative shallowness of the boundary layer. For both situations,it is shown that the turbulent fluxes can be computed with reasonably goodaccuracy (better than 10 %), taking into account both the random and thesystematic errors involved in the eddy-correlation technique. Estimationof random error is based on the computation of the integral scale of thecovariance, and systematic error is estimated from the parameterizationof Mann and Lenschow. The fluxes show that the buoyancy, as a sourceof turbulence, is due to latent heat flux rather than sensible heat flux,with values comparable to previous experiments in the Azores-Canariesbasin. In addition, we propose a method to analyse, for coupled situations,the relationship between the fractional cloudiness and the organization ofthe turbulent field below the clouds. This method is based on a conditionalsampling technique. It is shown that this organization cannot be deducedfrom the analysis of the velocity signal, which is dominated by turbulence.However, when the signals are conditionally sampled according to thepresence or absence of clouds, a weak cloud-related organization can beshown, and the cloud-related transports quantified; the values found areof the order of 10 % of the total transfers, i.e. the same order of magnitude asthe errors on the total flux computation. The method developed is thereforepromising, provided that the uncertainties can be reduced by analyzing a highamount of data.     

The measurement of turbulent surface-layer fluxes by use of bichromatic scintillation

Scintillation measurements with a HeNe and a CO₂ laser were used to derive turbulent fluxes of heat and momentum in the surface layer. This was achieved by the structure constant or dissipation technique, i.e., by relating the measured structure constants and inner scales of refractive index fluctuations to structure constants of temperature fluctuations and dissipation rates of turbulent kinetic energy, respectively, and then assuming Monin-Obukhov similarity.The resulting heat fluxes agree well with measurements using the eddy correlation technique but for averaging periods of 10 min, the optical data show a much smoother and physically more plausible behaviour. The optically derived friction velocities are in good agreement with estimates derived from wind velocity and surface roughness. It was also observed that for stationary conditions, 1-min averaged optical measurements already provide good estimates for longer averaged heat and momentum fluxes.Even though some uncertainty remains about the empirical constants and Monin-Obukhov similarity expressions used, the method clearly proves to be of great value for monitoring surface-layer turbulence.     

Coherent Structures Detected in Atmospheric Boundary-layer Turbulence using Wavelet transforms at Huaihe River Basin, China

The presence of coherent structures in turbulent shear flows suggests order in apparently random flows. These coherent structures play an important dynamical role in momentum and scalar transport. To develop dynamical models describing the evolution of such motion, it is necessary to detect and isolate the coherent structures from the background fluctuations. In this paper, we decomposed atmospheric turbulence time series into large-scale eddies, which include coherent structures and small eddies, which are stochastic by using Fourier digital filtering. The wavelet energy computed for the three components of the velocity fluctuations in the large-scale eddies appears to have local maximum values at certain time scales, which correspond to the scales or frequencies of coherent structures. We extract coherent signals from large-scale vortices at this scale by inverse wavelet transform formulae. This method provides an objective technique for examining the turbulence signal associated with coherent structures in the atmospheric boundary layer. The average duration of coherent structures in three directions based on Mexican hat wavelets are 33 s, 34 s and 25 s respectively. Symmetric andanti-symmetric wavelet basis functions give almost the same results. The main features of the structures during the day and night have little difference. The dimensionless durations for u, v and w have linear correlations with each other. These relationships are insensitive to the wavelet basis.     

Turbulent carbon exchange in very stable conditions

Turbulent fluxes obtained using the conventional eddy covariance approach result in erratic results with large time fluctuations in extremely stable conditions. This can limit efforts to estimate components of the nocturnal energy budget and respiratory CO₂ fluxes. Well-organized fluxes that show a clear dependence on turbulent intensity were obtained when multiresolution decomposition was used to estimate turbulent exchanges. CO₂, heat and water vapour fluxes were observed at a site in the eastern Amazon basin that had been cleared for agricultural purposes. Temporal scales of the carbon transfer were determined and shown to be similar to those of latent heat, but as much as three times larger than those of sensible heat. CO₂ eddy diffusivities at the temporal scales on which most of the vertical CO₂ exchange occurs are shown to be 50 times larger than the eddy diffusivity for heat. A process associated with the vertical scale of the scalar accumulation layer is suggested to explain these different scales and turbulent diffusivities of carbon and sensible heat transfer. For an appreciable range of turbulence intensities, the observed vertical turbulent carbon exchange is insufficient to account for the locally respired CO₂ estimated independently. Evidence that shallow drainage currents may account for this is given.