Anisotropy of Unstably Stratified Near-Surface Turbulence

Boundary-Layer Meteorology - Classic Monin–Obukov similarity scaling states that in a stationary, horizontally homogeneous flow, in the absence of subsidence, turbulence is dictated by the...     

在EBEX-2000实验资料中湍流耗散率、长度尺度和结构参数特征

采用2000年8月在美国加州棉花地两个高度上应用超声三分量仪、快速响应温度和湿度仪进行的EBEX-2000 (International Energy Balance Experiment, 2000, 简称EBEX-2000) 风速三分量、温度和湿度湍流实验观测数据, 计算分析了在不同稳定度下的湍流能量和热量耗散率和湍流结构参数特征.并与Kansas和长白山原始森林湍流实验得到的结果进行了比较, 得到了一些湍流特征量在不同下垫面情况下的一些有意义的特征.     

Role of Shear and the Inversion Strength During Sunset Turbulence Over Land: Characteristic Length Scales

The role of shear and inversion strength on the decay of convective turbulence during sunset over land is systematically studied by means of large-eddy simulations. Different decay rates have been found for the vertical and horizontal velocity fluctuations, resulting in an increase of the anisotropy for all the studied cases. Entrainment, which persists during the decay process, favours the appearance of vertical upward movements associated with a conversion from kinetic to potential energy. Particular attention is paid to the evolution of the characteristic length scale of the various turbulent variables during this process. The length scale evolution is found to depend on the wind shear characteristics, but not on the strength of the inversion. In general the length scales of the variables grow during decay because small-scale fluctuations dissipate faster than large-scale fluctuations. Only the length scale of the vertical velocity component remains nearly constant during decay. Spectral analysis of the variance budgets shows that pressure correlations are responsible for fixing this length scale, effectively compensating the strong but oscillating influence of buoyancy. In the shear cases, after an initial period of growth, the length scales start to decrease once the buoyancy-generated variance has sufficiently subsided. Also here the effect of pressure redistribution is crucial, as it transfers the spectral influence of shear to the other velocity components.     

The Near-Surface Stress Profile and Conditional Averaging of Turbulence

Using the conditional average formulation, we suggest a new explanation for why the stress in the atmospheric surface layer is often observed to vary with height. In essence, because turbulence series are always correlated for small lags, the steady-state equations of motion with negligible viscous terms that traditionally require vertical fluxes to be constant with height accordingly now require the vertical fluxes to vary with height. This result has implications for interpreting and validating Monin–Obukhov similarity theory.     

Local Scales of Turbulence in the Stable Boundary Layer

Local, gradient-based scales, which contain the vertical velocity and temperature variances, as well as the potential temperature gradient, but do not include fluxes, are tested using data collected during the CASES-99 experiment. The observations show that the scaling based on the temperature variance produces relatively smaller scatter of empirical points. The resulting dimensionless statistical moments approach constant values for sufficiently large values of the Richardson number Ri. This allows one to derive predictions for the Monin–Obukhov similarity functions φ _m and φ _h , the Prandtl number Pr and the flux Richardson number Rf in weak turbulence regime.     

The Use of Weather Forecasts to Characterise Near-Surface Optical Turbulence

The propagation of optical and electromagnetic waves is affected by small-scale atmospheric turbulence, quantified by the structure parameter of the refractive index. In the atmospheric surface layer, the mean structure parameter C_n²{C_{n}^{2}} , as averaged over the large-scale turbulence, relates to meteorological forcings through well-documented relationships. Present-day numerical weather forecast models routinely produce these forcings at the global scale. This study introduces a method where the products of such a model are used to calculate the mean optical turbulence near the surface. The method is evaluated against scintillometry measurements over climatologically distinct sites in Western Europe. The diurnal cycle modulation, and regional and seasonal contrasts, are all reproduced by our predictions. Hence, the present method explains and predicts some essential aspects of the meteorological variability of C_n²{C_{n}^{2}} near the surface. The noted discrepancies combine instrumental limitations, site peculiarities, differences related to distinct averaging procedures, and model errors, notably from weather forecasts. The minute-scale fluctuations of the measured scintillation rate are also analysed in the light of the forecast weather conditions. Fair-weather daytime periods consistently show a small short-term variability compared to the nighttime and perturbed weather periods. Thus, this short-term variability appears to have a predictable component.     

Impact of Boundary-Layer Processes on Near-Surface Turbulence Within the West African Monsoon

High frequency measurements of near-surface meteorological data acquired in north Benin during the 2006 West African monsoon seasonal cycle, in the context of the African Monsoon Multidisciplinary Analysis (AMMA) experiment, offer insight into the characteristics of surface turbulence in relation to planetary boundary-layer (PBL) processes. A wide range of conditions is encountered at the lower and upper limits of the PBL: (i) from water-stressed to well-fed vegetation, and (ii) from small to large humidity and temperature jumps at the PBL top inversion, due to the Saharan air layer overlying the monsoonal flow. As a result, buoyant convection at the surface and entrainment at the PBL top play very different roles according to the considered scalar. We show that, when the boundary-layer height reaches the shear level between the monsoonal and Harmattan flows, the temperature source and humidity sink at the boundary-layer top are sufficient to allow the entrainment to affect the entire boundary layer down to the surface. This situation occurs mainly during the drying and moistening periods of the monsoon cycle and affects the humidity statistics in particular. In this case, the humidity turbulent characteristics at the surface are no longer driven solely by buoyant convection, but also by entrainment at the boundary-layer top. Consequently, the Monin–Obukhov similarity theory appears to fail for the parameterisation of humidity-related moments.     

Comparing Estimates of Turbulence Based on Near-Surface Measurements in the Nocturnal Stable Boundary Layer

Tethered Lifting System (TLS) estimates of the dissipation rate of turbulent kinetic energy (e){(\varepsilon)} are reasonably well correlated with concurrent measurements of vertical velocity variance (s_w²){(\sigma_{w}^{2})} obtained from sonic anemometers located on a nearby 60-m tower during the CASES-99 field experiment. Additional results in the first 100 m of the nocturnal stable boundary layer confirm our earlier claim that the presence of weak but persistent background turbulence exists even during the most stable atmospheric conditions, where e{\varepsilon} can exhibit values as low as 10⁻⁷ m² s⁻³. We also present a set of empirical equations that incorporates TLS measurements of temperature, horizontal wind speed, and e{\varepsilon} to provide a proxy measurement for s_w²{\sigma_{w}^{2}} at altitudes higher than tower heights.     

广州市近地层大气的湍流微结构和谱特征

本文利用UVW脉动风速仪资料分析了广州市区近地层大气的湍流强度、相关系数、尺度和速度谱,并获得了不同稳定度条件下的速度谱模式.结果表明,城市近地层大气湍流在惯性副区接近局地各向同性、速度谱符合Kolmogorov相似理论;气流方向上下垫面粗糙度的增加,使沿海城市近地层大气湍流能量(特别是铅直方向)比平坦、均匀下垫面上的增加.     

云南大理干湿季近地层湍流特征对比分析

利用大理国家气候观象台2007年3月至2008年1月观测资料,采用涡动相关法等计算方法分析了该地区湍流强度、湍流方差、湍流通量等特征量的日变化规律和干湿季变化特征。结果表明：湍流强度干季大于湿季;湍流方差与稳定度满足1／3次方定律,风速方差在稳定条件下比不稳定条件下离散,水平方向比垂直方向离散;湍流通量有明显日变化特征,感热、动量通量干季大于湿季,潜热通量湿季大于干季,干湿季热量交换以潜热为主。     

Turbulence Structure of the Unstable Atmospheric Surface Layer and Transition to the Outer Layer

We present a new model of the structure of turbulence in the unstable atmospheric surface layer, and of the structural transition between this and the outer layer. The archetypal element of wall-bounded shear turbulence is the Theodorsen ejection amplifier (TEA) structure, in which an initial ejection of air from near the ground into an ideal laminar and logarithmic flow induces vortical motion about a hairpin-shaped core, which then creates a second ejection that is similar to, but larger than, the first. A series of TEA structures form a TEA cascade. In real turbulent flows TEA structures occur in distorted forms as TEA-like (TEAL) structures. Distortion terminates many TEAL cascades and only the best-formed TEAL structures initiate new cycles. In an extended log layer the resulting shear turbulence is a complex, self-organizing, dissipative system exhibiting self-similar behaviour under inner scaling. Spectral results show that this structure is insensitive to instability. This is contrary to the fundamental hypothesis of Monin--Obukhov similarity theory. All TEAL cascades terminate at the top of the surface layer where they encounter, and are severely distorted by, powerful eddies of similar size from the outer layer. These eddies are products of the breakdown of the large eddies produced by buoyancy in the outer layer. When the outer layer is much deeper than the surface layer the interacting eddies are from the inertial subrange of the outer Richardson cascade. The scale height of the surface layer, z _s, is then found by matching the powers delivered to the creation of emerging TEAL structures to the power passing down the Richardson cascade in the outer layer. It is z _s = u _* ³ /k_s, where u _* is friction velocity, k is the von Kármán constant and _s is the rate of dissipation of turbulence kinetic energy in the outer layer immediately above the surface layer. This height is comparable to the Obukhov length in the fully convective boundary layer. Aircraft and tower observations confirm a strong qualitative change in the structure of the turbulence at about that height. The tallest eddies within the surface layer have height z _s, so z _s is a new basis parameter for similarity models of the surface layer.     

Water-Channel Estimation of Eulerian and Lagrangian Time Scales of the Turbulence in Idealized Two-Dimensional Urban Canopies

Lagrangian and Eulerian statistics are obtained from a water-channel experiment of an idealized two-dimensional urban canopy flow in neutral conditions. The objective is to quantify the Eulerian \((T^{\mathrm{E}})\) and Lagrangian \((T^{\mathrm{L}})\) time scales of the turbulence above the canopy layer as well as to investigate their dependence on the aspect ratio of the canopy, AR, as the latter is the ratio of the width (W) to the height (H) of the canyon. Experiments are also conducted for the case of flat terrain, which can be thought of as equivalent to a classical one-directional shear flow. The values found for the Eulerian time scales on flat terrain are in agreement with previous numerical results found in the literature. It is found that both the streamwise and vertical components of the Lagrangian time scale, \(T_\mathrm{u}^\mathrm{L} \) and \(T_\mathrm{w}^\mathrm{L} \), follow Raupach’s linear law within the constant-flux layer. The same holds true for \(T_\mathrm{w}^\mathrm{L} \) in both the canopies analyzed \((AR= 1\) and \(AR= 2\)) and also for \(T_\mathrm{u}^\mathrm{L} \) when \(AR = 1\). In contrast, for \(AR = 2\), \(T_\mathrm{u}^\mathrm{L} \) follows Raupach’s law only above \(z=2H\). Below that level, \(T_\mathrm{u}^\mathrm{L} \) is nearly constant with height, showing at \(z=H\) a value approximately one order of magnitude greater than that found for \(AR = 1\). It is shown that the assumption usually adopted for flat terrain, that \(\beta =T^{\mathrm{L}}/T^{\mathrm{E}}\) is proportional to the inverse of the turbulence intensity, also holds true even for the canopy flow in the constant-flux layer. In particular, \(\gamma /i_\mathrm{u} \) fits well \(\beta _\mathrm{u} =T_\mathrm{u}^\mathrm{L} /T_\mathrm{u}^\mathrm{E} \) in both the configurations by choosing \(\gamma \) to be 0.35 (here, \(i_\mathrm{u} =\sigma _\mathrm{u} / \bar{u} \), where \(\bar{u} \) and \(\sigma _\mathrm{u} \) are the mean and the root-mean-square of the streamwise velocity component, respectively). On the other hand, \(\beta _\mathrm{w} =T_\mathrm{w}^\mathrm{L} /T_\mathrm{w}^\mathrm{E} \) follows approximately \(\gamma /i_\mathrm{w} =0.65/\left( {\sigma _\mathrm{w} /\bar{u} } \right) \) for \(z > 2H\), irrespective of the AR value. The second main objective is to estimate other parameters of interest in dispersion studies, such as the eddy diffusivity of momentum \((K_\mathrm{{T}})\) and the Kolmogorov constant \((C_0)\). It is found that \(C_0\) depends appreciably on the velocity component both for the flat terrain and canopy flow, even though for the latter case it is insensitive to AR values. In all the three experimental configurations analyzed here, \(K_\mathrm{{T}}\) shows an overall linear growth with height in agreement with the linear trend predicted by Prandtl’s theory.     

On the Lagrangian and Eulerian Time Scales of Turbulence Within a Two-Dimensional Array of Obstacles

Boundary-Layer Meteorology - Fields of Lagrangian ( $$T^{L}$$ ) and Eulerian ( $$T^{E}$$ ) time scales of the turbulence within a regular array of two-dimensional obstacles of unit aspect ratio...     

Length Scales of Scalar Diffusion in the Convective Boundary Layer: Laboratory Observations

A laboratory study of scalar diffusion in the convective boundary layer has found results that are consistent with a 1999 large-eddy simulation (LES) study by Jonker, Duynkerke and Cuijpers. For bottom-up and top-down scalars (introduced as ‘infinite’ area sources of passive tracer at the surface and inversion, respectively) the dominant length scale was found to be much larger than the length scale for density fluctuations, the latter being equal to the boundary-layer depth h. The variance of the normalized passive scalar grew continuously with time and its magnitude was about 3–5 times larger for the top-down case than for the bottom-up case. The vertical profiles of the normalized passive scalar variance were found to be approximately constant through the convective boundary layer (CBL) with a value of about 3–8c_*² for bottom-up and 10–50c_*² for top-down diffusion. Finally, there was some evidence of a minimum in the variance and dominant length scale for scalar flux ratios (top-down to bottom-up flux) close to −0.5. All these convection tank results confirm the LES results and support the hypothesis that there is a distinct difference in behaviour between the dynamic and passive variables in the CBL.     

Hurricane Wind Power Spectra, Cospectra, and Integral Length Scales

Atmospheric turbulence is an important factor in the modelling of wind forces on structures and the losses they produce in extreme wind events. However, while turbulence in non-hurricane winds has been thoroughly researched, turbulence in tropical cyclones and hurricanes that affect the Gulf and Atlantic coasts has only recently been the object of systematic study. In this paper, Florida Coastal Monitoring Program surface wind measurements over the sea surface and open flat terrain are used to estimate tropical cyclone and hurricane wind spectra and cospectra as well as integral length scales. From the analyses of wind speeds obtained from five towers in four hurricanes it can be concluded with high confidence that the turbulent energy at lower frequencies is considerably higher in hurricane than in non-hurricane winds. Estimates of turbulence spectra, cospectra, and integral turbulence scales presented can be used for the development in experimental facilities of hurricane wind flows and the forces they induce on structures.     

Eddy Covariance Tilt Corrections over a Coastal Mountain Area in South-east China: Significance for Near-Surface Turbulence Characteristics

Turbulence characteristics of an atmospheric surface layer over a coastal mountain area were investigated under different coordinate frames.Performances of three methods of coordinate rotation:double rotation(DR),triple rotation(TR),and classic planar-fit rotation(PF) were examined in terms of correction of eddy covariance flux.Using the commonly used DR and TR methods,unreasonable rotation angles are encountered at low wind speeds and cause significant run-to-run errors of some turbulence characteristics.The PF method rotates the coordinate system to an ensemble-averaged plane,and shows large tilt error due to an inaccurate fit plane over variable terrain slopes.In this paper,we propose another coordinate rotation scheme.The observational data were separated into two groups according to wind direction.The PF method was adapted to find an ensemble-averaged streamline plane for each group of hourly runs with wind speed exceeding 1.0 m s 1.Then,the coordinate systems were rotated to their respective bestfit planes for all available hourly observations.We call this the PF10 method.The implications of tilt corrections for the turbulence characteristics are discussed with a focus on integral turbulence characteristics,the spectra of wind-velocity components,and sensible heat and momentum fluxes under various atmospheric stabilities.Our results show that the adapted application of PF provides greatly improved estimates of integral turbulence characteristics in complex terrain and maintains data quality.The comparisons of the sensible heat fluxes for four coordinate rotation methods to fluxes before correction indicate that the PF10 scheme is the best to preserve consistency between fluxes.