大气边界层湍流温度序列的信息熵分析

利用大气边界层内近地面的大气湍流温度时间序列，运用功率谱分析、信息熵分析等方法，分析了大气边界层内近地面的大气湍流特点，并对稳定层和不稳定层的大气湍流进行了对比。结果表明，信息熵和功率谱指数是区别稳定层结和不稳定层结大气边界层湍流特征的指标，对造成两者之间的差别做出了对应的解释。     

A Numerically Efficient Parametrization of Turbulent Wind-Turbine Flows for Different Thermal Stratifications

The wake characteristics of a wind turbine in a turbulent atmospheric boundary layer under different thermal stratifications are investigated by means of large-eddy simulation with the geophysical flow solver EULAG. The turbulent inflow is based on a method that imposes the spectral energy distribution of a neutral boundary-layer precursor simulation, the turbulence-preserving method. This method is extended herein to make it applicable for different thermal stratification regimes (convective, stable, neutral) by including suitable turbulence assumptions, which are deduced from velocity fields of a diurnal-cycle precursor simulation. The wind-turbine-wake characteristics derived from simulations that include the parametrization result in good agreement with diurnal-cycle-driven wind-turbine simulations. Furthermore, different levels of accuracy are tested in the parametrization assumptions, representing the thermal stratification. These range from three-dimensional matrices of the precursor-simulation wind field to individual values. The resulting wake characteristics are similar, even for the simplest parametrization set-up, making the diurnal-cycle precursor simulation non-essential for the wind-turbine simulations. Therefore, the proposed parametrization results in a computationally fast, simple, and efficient tool for analyzing the effects of different thermal stratifications on wind-turbine wakes by means of large-eddy simulation.     

Different Turbulent Regimes and Vertical Turbulence Structures of the Urban Nocturnal Stable Boundary Layer

Turbulence in the nocturnal boundary layer(NBL) is still not well characterized, especially over complex underlying surfaces. Herein, gradient tower data and eddy covariance data collected by the Beijing 325-m tower were used to better understand the differentiating characteristics of turbulence regimes and vertical turbulence structure of urban the NBL. As for heights above the urban canopy layer(UCL), the relationship between turbulence velocity scale(V_TKE) and wind speed(V) was con...     

Stability Dependence of Canopy Flows over a Flat Larch Forest

The dependence on atmospheric stability of flow characteristics adjacent to a very rough surface was investigated in a larch forest in Japan. Micrometeorological measurements of three-dimensional wind velocity and air temperature were taken at two heights above the forest, namely 1.7 and 1.2 times the mean canopy height h. Under near-neutral and stable conditions, the observed turbulence statistics suggest that the flow was likely to be that of the atmospheric surface layer (ASL) at 1.7h, and of the roughness sublayer (RSL) at 1.2h. However, in turbulence spectra, canopy-induced large coherent motions appeared clearly at both heights. Even under strongly stable conditions, the large-scale motions were retained at 1.2h, whereas they were overwhelmed by small-scale motions at 1.7h. This phenomenon was probably due to the enhanced contribution of the ASL turbulence associated with nocturnal decay of the RSL depth, because the small-scale motions appeared at frequencies close to the peak frequencies of well-known ASL spectra. This result supports the relatively recent concept that canopy flow is a superimposition of coherent motions and the ASL turbulence. The large-scale motions were retained in temperature spectra over a wider region of stability compared to streamwise wind spectra, suggesting that a canopy effect extended higher up for temperature than wind. The streamwise spacing of dominant eddies according to the plane mixing-layer analogy was only valid in a narrow range at near neutral, and it was stabilised at nearly half its value under stable conditions.     

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

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

Temporal Oscillations in the Convective Boundary Layer Forced by Mesoscale Surface Heat-Flux Variations

A theoretical approach suggests that the surface heterogeneity on a scale of tens of kilometres can generate mesoscale motions that are not in a quasi-stationary state. The starting point of the theoretical approach is the equations of horizontal velocity and potential temperature that are low-pass filtered with a mesoscale cut-off wavelength. The transition of the generated mesoscale motions from a quasi-stationary state to a non-stationary state occurs when horizontal advection is strong enough to level out the potential temperature gradient on the surface heterogeneity scale. Large-eddy simulations (LES) suggest that the convective boundary layer (CBL) changes to a non-stationary state when forced by a surface heat-flux variation of amplitude of 100W m⁻² or higher and a wavelength of the order of 10 km. Spectral analysis of the LES reveals that when the mesoscale motions are in a quasi-stationary state, the energy provided by the surface heat-flux variation remains in organized mesoscale motions on the scale of the surface variation itself. However, in a non-stationary state, the energy cascades to smaller scales, with the cascade extending down into the turbulence scale when the wavelength of the surface heat-flux variation is on a scale smaller than 100 times the CBL height. The energy transfer from the generated mesoscale motions to the CBL turbulence results in the absence of a spectral gap between the two scales. The absence of an obvious spectral gap between the generated mesoscale motions and the turbulence raises questions about the applicability of mesoscale models for studies on the effect of high-amplitude surface heterogeneity on a scale of tens of kilometres. The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

50 Years of the Monin–Obukhov Similarity Theory

In weak wind stable conditions, eddy-correlation fluxes calculated using conventional averaging times of 5 min or longer to define the perturbations are severely contaminated by poorly sampled mesoscale motions. A method is developed to identify the averaging time for each individual data record that captures the turbulence while excluding most of the mesoscale motions. The method is based on multiresolution decomposition of the heat flux, and provides an objective procedure for selecting the averaging time for calculating eddy-correlation fluxes. Eddy-correlation data collected in weak turbulence conditions over grass, snow, a pine forest and the ocean are used to demonstrate the approach.When the small-scale turbulence and mesoscale motions are clearly separated by a gap region in the heat flux cospectra, the variable window width reduces the influence of nonstationarity by more effectively filtering out mesoscale motions compared to traditional methods using constant averaging time. For records where turbulence and mesoscale motions overlap in scale, the method is not well posed, although such records occur infrequently for our datasets. These ambiguous cases correspond to significant nonstationarity at scales that overlap with turbulence scales. The improved turbulence fluxes calculated with the proposed method are the appropriate fluxes for evaluating flux-gradient relationships and Monin–Obukov similarity theory for developing improved model parameterizations of turbulence for weakly turbulent flows     

3mm多普勒云雷达测量反演云内空气垂直速度的研究

垂直指向的W-波段大气云雷达(WACR)不仅测量云粒子的反射率因子Z, 而且测量多普勒速度谱。本文利用寿县气象站WACR在2008年11月1日一个时段的测量, 进行云内空气垂直速度反演试验。首先讨论个例中反射率、 Doppler平均速度及谱宽的分布特征及原因, 然后利用小粒子示踪法和改进的以小粒子示踪法为基础的粒子下降速度w0—反射率因子 (w0－Z) 关系法反演了云内空气垂直速度, 进而结合Doppler谱数据分析反演结果。结果显示: 小粒子示踪法在湍流较弱时能比较精确地反演空气垂直速度, 而湍流较强时, 湍流造成的误差不可忽略; 改进算法在湍流较强时能够减少湍流对反演结果的影响。     

Variability and Maintenance of Turbulence in the Very Stable Boundary Layer

The relationship of turbulence quantities to mean flow quantities, such as the Richardson number, degenerates substantially for strong stability, at least in those studies that do not place restrictions on minimum turbulence or non-stationarity. This study examines the large variability of the turbulence for very stable conditions by analyzing four months of turbulence data from a site with short grass. Brief comparisons are made with three additional sites, one over short grass on flat terrain and two with tall vegetation in complex terrain. For very stable conditions, any dependence of the turbulence quantities on the mean wind speed or bulk Richardson number becomes masked by large scatter, as found in some previous studies. The large variability of the turbulence quantities is due to random variations and other physical influences not represented by the bulk Richardson number. There is no critical Richardson number above which the turbulence vanishes. For very stable conditions, the record-averaged vertical velocity variance and the drag coefficient increase with the strength of the submeso motions (wave motions, solitary waves, horizontal modes and numerous more complex signatures). The submeso motions are on time scales of minutes and not normally considered part of the mean flow. The generation of turbulence by such unpredictable motions appears to preclude universal similarity theory for predicting the surface stress for very stable conditions. Large variation of the stress direction with respect to the wind direction for the very stable regime is also examined. Needed additional work is noted.     

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.     

The Effects of Thermal Stratification on Clustering Properties of Canopy Turbulence

The intermittent structure of turbulence within the canopy sublayer (CSL) is sensitive to the presence of foliage and to the atmospheric stability regime. How much of this intermittency originates from amplitude variability or clustering properties remains a vexing research problem for CSL flows. Using a five-level set of measurements collected within a dense hardwood canopy, the clustering properties of CSL turbulence and their dependence on atmospheric stability are explored using the telegraphic approximation (TA). The binary structure of the TA removes any amplitude variability from turbulent excursions but retains their zero-crossing behaviour, and thereby isolating the role of clustering in intermittency. A relationship between the spectral exponents of the actual and the TA series is derived across a wide range of atmospheric stability regimes and for several flow variables. This relationship is shown to be consistent with a relationship derived for long-memory and monofractal processes such as fractional Brownian motion (fBm). Moreover, it is demonstrated that for the longitudinal and vertical velocity components, the vegetation does not appreciably alter fine-scale clustering but atmospheric stability does. Stable atmospheric stability conditions is characterized by more fine scale clustering when compared to other atmospheric stability regimes. For scalars, fine-scale clustering above the canopy is similar to its velocity counterpart but is significantly increased inside the canopy, especially under stable stratification. Using simplified scaling analysis, it is demonstrated that clustering is much more connected to space than to time within the CSL. When comparing intermittency for flow variables and their TA series, it is shown that for velocity, amplitude variations modulate intermittency for all stability regimes. However, amplitude variations play only a minor role in scalar intermittency. Within the crown region of the canopy, a ‘double regime’ emerges in the inter-pulse duration probability distributions not observed in classical turbulence studies away from boundaries. The double regime is characterized by a power-law distribution for shorter inter-pulse periods and a log-normal distribution for large inter-pulse periods. The co-existence of these two regimes is shown to be consistent with near-field/far-field scaling arguments. In the near-field, short inter-pulse periods are controlled by the source strength, while in the far-field long inter-pulse periods are less affected by the precise source strength details and more affected by the transport properties of the background turbulence.     

Coherent structures and flux contribution over an inhomogeneously irrigated cotton field

The turbulence data measured at two levels (i.e., 8.7 and 2.7?m) in the Energy Balance Experiment (EBEX), which was conducted in San Joaquin Valley in California during the period from July 20 to August 24, 2000, are used to study the characteristics of coherent structures over an irrigated cotton field. Patch-to-patch irrigation in the field generated the dry-to-wet horizontal advection and the oasis effects, leading to the development of a stably internal boundary layer (SIBL) in the late mornings or the early afternoons. The SIBL persisted in the rest of the afternoons. Under this circumstance, a near-neutral atmospheric surface layer (ASL) developed during the period with a stratification transition from the unstable to stable conditions during the daytime. Therefore, EBEX provides us with unique datasets to investigate the features of coherent structures that were generated over the patches upstream and passed by our site in the unstable ASL, the near-neutral ASL, and the SIBL. We use an objective detection technique and the conditional average method that is developed based on the wavelet analysis. Our data reveal some consistencies and inconsistencies in the characteristics of coherent structures as compared with previous studies. Ramp-like structures and sweep?Cejection cycles under the daytime SIBL have similar patterns to those under the nocturnal stable ASL. However, some features (i.e., intermittence) are different from those under the nocturnal stable ASL. Under the three stratifications, thermal and mechanical factors in the ASL perform differently in affecting the ramp intensity for different quantities (i.e., velocity components, temperature, and specific humidity), leading to coherent structures that modulate turbulence flow and alter turbulent transfer differently. It is also found that coherent structures contribute about 10?C20% to the total fluxes in our case with different flux contributions under three stratifications and with higher transporting efficiency in sensible heat flux than latent heat and momentum fluxes.     

On the influence of variations in Monin-Obukhov length on horizontal roll vortices in an inversion-capped planetary boundary layer

Observations of the planetary boundary layer under conditions with strong or moderate winds often show the large-scale boundary-layer motions to be highly elongated in a direction close to that of the geostrophic wind. The properties of such large-scale motions are examined by means of a two-dimensional numerical model. The small-scale turbulence is parameterized using a buoyancy-dependent mixing-length hypothesis. The objective is to understand how the properties and dynamics of such rolls depend upon the relative importance of shear and buoyancy forces.     

Organised Motion and Radiative Perturbations in the Nocturnal Canopy Sublayer above an Even-Aged Pine Forest

Using time series measurements of velocity, carbon dioxideand water vapour concentration, and temperature collected justabove a 15 m tall even-aged pine forest, we quantify the roleof organized motion on scalar and momentum transport withinthe nocturnal canopy sublayer (CSL). We propose a frameworkin which the nocturnal CSL has two end-members, bothdominated by organised motion. These end-members representfully developed turbulent flows at near-neutral or slightly stablestratification and no turbulence for very stable stratification.Our analysis suggests that ramps dominate scalar transport fornear-neutral and slightly stable conditions, while linear canopywaves dominate the flow dynamics for very stable conditions.For intermediate stability, the turbulence is highly damped andoften dominated by fine scale motions. Co-spectral analysissuggests that ramps are the most efficient net scalar mass-transportingagent while linear canopy waves contribute little to net scalartransport between the canopy and atmosphere for averagingintervals that include complete wave cycles. However, canopywaves significantly contribute to the spectral properties of thescalar time series. Ramps are the most frequently occurringorganised motion in the nocturnal CSL for this site.Numerous night-time runs, however, resided between thesetwo end-members. Our analysis suggests that whenradiative perturbations are sufficient large (>20 W m^-2 innet radiation), the flow can switch from being highly dampedfine-scale turbulence to being organized with ramp-like properties. We also found that when ramps are already the dominant eddymotion in the nocturnal CSL, radiative perturbations have aminor impact on scalar transport. Finally, in agreement withprevious studies, we found that ramps and canopy waves havecomparable length scales of about 30–60 metres. Consequencesto night-time flux averaging are also discussed.     

Attached Eddies and Production Spectra in the Atmospheric Logarithmic Layer

We investigate the production components of turbulent spectra within logarithmic layers over flat ground. This assumes that the turbulence giving rise to these spectra consists of active coherent structures (eddies) that are attached to the ground, and whose properties display perfect statistical self-similarity under inner scaling. That is, we take the extreme view that active coherent structures not only contribute to turbulence production spectra but explain the whole of them, so that neither detached eddies nor unstructured motions make any significant contribution. Perfect self-similarityis held to apply only to eddies that are themselves formed totally within the log layer, so the theory applies in the limit of spectra obtained at the hearts of very deep log layers. The model predicts that spectral variance and covariance should become independent of wavenumber at small wavenumbers. This asymptotic behaviour is observed in all neutral spectra from the Kansas experiment. The model also interprets the various positions of the spectral peaks observed at Kansas and in aircraft flights over the sea as consequences of the eddies being aggregated into files aligned with the wind. The observed spectra are therefore consistent with large-scale wedge-like structures being the principal component of active turbulence in the neutral atmospheric surface layer.     

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.     

A Three-Dimensional Backward Lagrangian Footprint Model For A Wide Range Of Boundary-Layer Stratifications

We present a three-dimensional Lagrangian footprint model with the ability to predict the area of influence (footprint) of a measurement within a wide range of boundary-layer stratifications and receptor heights. The model approach uses stochastic backward trajectories of particles and satisfies the well-mixed condition in inhomogeneous turbulence for continuous transitions from stable to convective stratification. We introduce a spin-up procedure of the model and a statistical treatment of particle touchdowns which leads to a significant reduction of CPU time compared to conventional footprint modelling approaches. A comparison with other footprint models (of the analytical and Lagrangian type) suggests that the present backward Lagrangian model provides valid footprint predictions under any stratification and, moreover, for applications that reach across different similarity scaling domains (e.g., surface layer to mixed layer, for use in connection with aircraft measurements or with observations on high towers).     

The Near-Calm Stable Boundary Layer

For the near-calm stable boundary layer, nominally 2-m mean wind speed <0.5 ms⁻¹, the time-average turbulent flux is dominated by infrequent mixing events. These events are related to accelerations associated with wave-like motions and other more complex small-scale motions. In this regime, the relationship between the fluxes and the weak mean flow breaks down. Such near-calm conditions are common at some sites. For very weak winds and strong stratification, the characteristics of the fluctuating quantities change slowly with increasing scale and the separation between the turbulence and non-turbulent motions can become ambiguous. Therefore, a new analysis strategy is developed based on the scale dependence of selected flow characteristics, such as the ratio of the fluctuating potential energy to the kinetic energy. In contrast to more developed turbulence, correlations between fluctuating quantities are small, and a significant heat flux is sometimes carried by very weak vertical motions with large temperature fluctuations. The relation of the flux events to small-scale increases of wind speed is examined. Large remaining uncertainties are noted.     

Global Intermittency and Collapsing Turbulence in the Stratified Planetary Boundary Layer

Direct numerical simulation of the turbulent Ekman layer over a smooth wall is used to investigate bulk properties of a planetary boundary layer under stable stratification. Our simplified configuration depends on two non-dimensional parameters: a Richardson number characterizing the stratification and a Reynolds number characterizing the turbulence scale separation. This simplified configuration is sufficient to reproduce global intermittency, a turbulence collapse, and the decoupling of the surface from the outer region of the boundary layer. Global intermittency appears even in the absence of local perturbations at the surface; the only requirement is that large-scale structures several times wider than the boundary-layer height have enough space to develop. Analysis of the mean velocity, turbulence kinetic energy, and external intermittency is used to investigate the large-scale structures and corresponding differences between stably stratified Ekman flow and channel flow. Both configurations show a similar transition to the turbulence collapse, overshoot of turbulence kinetic energy, and spectral properties. Differences in the outer region resulting from the rotation of the system lead, however, to the generation of enstrophy in the non-turbulent patches of the Ekman flow. The coefficient of the stability correction function from Monin–Obukhov similarity theory is estimated as \(\beta \approx 5.7\) in agreement with atmospheric observations, theoretical considerations, and results from stably stratified channel flows. Our results demonstrate the applicability of this set-up to atmospheric problems despite the intermediate Reynolds number achieved in our simulations.     

Turbulence Spectra and Cospectra Under the Influence of Large Eddies in the Energy Balance EXperiment (EBEX)

During the Energy Balance EXperiment, the patch-by-patch, flood irrigation in a flat cotton field created an underlying surface with heterogeneous soil moisture, leading to a dry (warm)-to-wet (cool) transition within the cotton field under northerly winds. Moreover, the existence of an extremely dry, large bare soil area upstream beyond the cotton field created an even larger step transition from the bare soil region to the cotton field. We investigated the turbulence spectra and cospectra in the atmospheric surface layer (ASL) that was disturbed by large eddies generated over regions upstream and also influenced by horizontal advection. In the morning, the ASL was unstable while in the afternoon a stable internal boundary layer was observed at the site. Therefore, the turbulence data at 2.7 and 8.7 m are interpreted and compared in terms of interactions between large eddies and locally generated turbulence under two atmospheric conditions: the unstable ASL beneath the convective boundary layer (CBL) (hereafter the unstable condition) and the stable ASL beneath the CBL (hereafter the stable condition). We identified the influences of multiple sizes of large eddies on ASL turbulence under both stratifications; these large eddies with multiple sizes were produced over the dry patches and dry, large bare soil areas upstream. As a consequence of the disturbance of large eddies, the broadening, erratic variability, and deviation of spectra and cospectra, relative to those described by Monin–Obukhov similarity theory, are evident in the low- to mid-frequencies. Transfer of momentum, heat, and water vapour by large eddies is distinctly observed from the turbulence cospectra and leads to significant run-to-run variations of residuals of the surface energy balance closure. Our results indicate that these large eddies have greater influences on turbulence at higher levels compared to lower levels, and in the unstable ASL compared to the stable ASL.