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.     

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     

黄土高原复杂地形受中尺度运动影响的稳定边界层湍流特征

利用兰州大学半干旱气候与环境观测站(Semi-Arid Climate and Environment Observatory of Lanzhou University,简称SACOL)2008年12月观测资料,研究了稳定边界层湍流特征.使用涡动相关资料研究湍流通量时,定义湍流的平均时间τ内的中尺度运动是造成湍流统计量变化范围大的主要原因,稳定情形? τ取几十秒至几分钟.对梯度理查森数大于0.3的强稳定情形的湍流尺度分解(MRD)谱分析表明,感热通量在112.4～449.9 s存在谱隙,尺度大于谱隙的中尺度运动造成了通量观测资料离散性大,甚至有支配性影响.动量通量的谱隙在112.4～224.9 s之间.弱风时,中尺度运动的影响更大,垂直风速标准差以0.1的比率随中尺度风速变化;垂直风速标准差同广义风速表现出很好的相关性,并随着广义风速消失而消失.三维风速标准差与摩擦速度呈很好的线性关系,垂直、水平、横风风速的无量纲标准差分别为1.35、2.54、2.21.对湍流动能的研究发现,在梯度理查森数大于0.3的条件下,仍然存在连续的湍流.以湍动能为依据,分析了湍流的平稳时间长度,其长度随稳定度变化而变化,2008年12月7～11日从133.5 s变化到856.2 s,湍流平稳时间长度反映了中尺度运动的发生频率.     

Large‐eddy simulation of small‐scale surface effects on the convective boundary‐layer structure

Abstract

The effects of small‐scale surface inhomogeneities on the turbulence structure in the convective boundary layer are investigated using a high‐resolution large‐eddy simulation model. Surface heat flux variations are sinusoidal and two‐dimensional, dividing the total domain into a checkerboard‐like pattern of surface hot spots with a 500‐m wavelength in the x and y directions, or 1/4 of the domain size. The selected wind speeds were 1 and 4 m s^‐l, respectively. As a comparison, a simulation of the turbulence structure was performed over a homogeneous surface.

When the wind speed is light, surface heat flux variations influence the horizontally averaged turbulence statistics, including the higher moments despite the small characteristic length of the surface perturbation. Stronger mean wind speeds weaken the effects of inhomogeneous surface conditions on the turbulence structure in the convective boundary layer.

Results from conditional sampling show that when the mean wind speed is small, weak mean circulations occur, with updraft branches above the high heat flux regions and down‐draft branches above the low heat flux regions. The inhomogeneous surface induces significant differences in the turbulence statistics between the high and low heat flux regions. However, the effect of the surface perturbations weaken rapidly when the mean wind speed increases. This research has implications in the explanation of the large‐scale variability commonly encountered in aircraft observations of atmospheric turbulence.     

Aspects Of The Atmospheric Surface Layers On Mars And Earth

The structures of mean flow and turbulence in the atmospheric surface boundary layer have been extensively studied on Earth, and to a far less extent on Mars, where only the Viking missions and the Pathfinder mission have delivered in-situ data. Largely the behaviour of surface-layer turbulence and mean flow on Mars is found to obey the same scaling laws as on Earth. The largest micrometeorological differences between the two atmospheres are associated with the low air density of the Martian atmosphere. Together with the virtual absence of water vapour, it reduces the importance of the atmospheric heat flux in the surface energy budget. This increases the temperature variation of the surface forcing the near-surface temperature gradient and thereby the diabatic heat flux to higher values than are typical on the Earth, resulting in turn in a deeper daytime boundary layer. As wind speed is much like that of the Earth, this larger diabatic heat flux is carried mostly by larger maximal values of T_*, the surface scale temperature. The higher kinematic viscosity yields a Kolmogorov scale of the order of ten times larger than on Earth, influencing the transition between rough and smooth flow for the same surface features.The scaling laws have been validated analysing the Martian surface-layer data for the relations between the power spectra of wind and temperature turbulence and the corresponding mean values of wind speed and temperature. Usual spectral formulations were used based on the scaling laws ruling the Earth atmospheric surface layer, whereby the Earth's atmosphere is used as a standard for the Martian atmosphere.     

Cross‐shore variations of near‐surface wind velocity and atmospheric turbulence at the land‐sea boundary during CASP

Abstract

Airborne measurements of mean wind velocity and turbulence in the atmospheric boundary layer under wintertime conditions of cold offshore advection suggest that at a height of 50 m the mean wind speed increases with offshore distance by roughly 20% over a horizontal scale of order 10 km. Similarly, the vertical gust velocity and turbulent kinetic energy decay on scales of order 3.5 km by factors of 1.5 and 3.2, respectively. The scale of cross‐shore variations in the vertical fluxes of heat and downwind momentum is also 10 km, and the momentum flux is found to be roughly constant to 300 m, whereas the heat flux decreases with height. The stability parameter, z/L (where z = 50 m and L is the local Monin‐Obukhov length), is generally small over land but may reach order one over the warm ocean. The magnitude and horizontal length scales associated with the offshore variations in wind speed and turbulence are reasonably consistent with model results for a simple roughness change, but a more sophisticated model is required to interpret the combined effects of surface roughness and heat flux contrasts between land and sea.

Comparisons between aircraft and profile‐adjusted surface measurements of wind speed indicate that Doppler biases of 1–2 m s^?1 in the aircraft data caused by surface motions must be accounted for. In addition, the wind direction measurements of the Minimet anemometer buoy deployed in CASP are found to be in error by 25 ± 5°, possibly due to a misalignment of the anemometer vane. The vertical fluxes of heat and momentum show reasonably good agreement with surface estimates based on the Minimet data.     

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.     

Determination Of The Surface Drag Coefficient

This study examines the dependence of the surface drag coefficienton stability, wind speed, mesoscale modulation of the turbulent flux and method of calculation of the drag coefficient. Data sets over grassland, sparse grass, heather and two forest sites are analyzed. For significantly unstable conditions, the drag coefficient does not depend systematically on z/L but decreases with wind speed for fixed intervals of z/L, where L is the Obukhov length. Even though the drag coefficient for weak wind conditions is sensitive to the exact method of calculation and choice of averaging time, the decrease of the drag coefficient with wind speed occurs for all of the calculation methods. A classification of flux calculation methods is constructed, which unifies the most common previous approaches.The roughness length corresponding to the usual Monin–Obukhovstability functions decreases with increasing wind speed. This dependence on wind speed cannot be eliminated by adjusting the stability functions. If physical, the decrease of the roughness length with increasing wind speed might be due to the decreasing role of viscous effectsand streamlining of the vegetation, although these effects cannot be isolated from existing atmospheric data.For weak winds, both the mean flow and the stress vector often meander significantly in response to mesoscale motions. The relationship between meandering of the stress and wind vectors is examined. For weak winds, the drag coefficient can be sensitive to the method of calculation, partly due to meandering of the stress vector.     

An Analysis Of Secondary Circulations And Their Effects Caused By Small-Scale Surface Inhomogeneities Using Large-Eddy Simulation

A parallelized large-eddy simulation model has been used to investigate the effects of two-dimensional, discontinuous, small-scale surface heterogeneities on the turbulence structure of the convective boundary layer.Heterogeneities had a typical size of about the boundary-layer heightz_i. They were produced by a surface sensible heat flux pattern ofchessboard-type and of strong amplitude as typical, e.g., for the marginalice zone. The major objectives of this study were to determinethe effects of such strong amplitude heat flux variations and to specify theinfluence of different speeds and directions of the background wind.Special emphasis has been given to investigate the secondary circulations induced by the heterogeneities by means of three-dimensional phase averages.Compared with earlier studies of continuous inhomogeneities, the same sizeddiscontinuous inhomogeneities in this study show similar but stronger effects.Significant changes compared with uniform surface heating are only observedwhen the scale of the inhomogeneities is increased to z_i. Especially the vertical energy transport is much more vigorous and even the mean emperature profile shows a positive lapse rate within the whole mixed layer. However, the effects are not directly caused by the different shape of the inhomogeneities but can mainly be attributed to the large amplitude of the imposed heat flux,as it is typical for the partially ice covered sea during cold air outbreaks.The structure of the secondary flow is found to be very sensitive to the wavelength and shape of the inhomogeneities as well as to the heatflux amplitude, wind speed and wind direction. The main controlling parameter is the near-surface temperature distribution and the related horizontal pressure gradient perpendicular to the main flow direction. The secondary flow varies from a direct circulation with updraughts mainly above the centre of the heated regions to a more indirect circulation with updraughts beneath the centre and downdraughts above it. For background winds larger than 2.5 m s^–1 a roll-like circulation pattern is observed.From previous findings it has often been stated that moderate backgroundwinds of 5 m s^–1 eliminate all impacts of surface inhomogeneitiesthat could potentially be produced in realistic landscapes. However, this studyshows that the effects caused by increasing the wind speed stronglydepend on the wind direction relative to the orientation of theinhomogeneities. Secondary circulations remain strong, even for abackground wind of 7.5 m s^–1, when the wind direction is orientatedalong one of the two diagonals of the chessboard pattern. On the otherhand, the effects of inhomogeneities are considerably reduced, even undera modest background wind of 2.5 m s^–1, if the wind direction isturned by 45°. Mechanisms for the different flow regimesare discussed.     

2009年“莫拉克”台风登陆过程阵风特征分析

利用上海台风研究所移动观测车获取的“莫拉克”台风登陆过程中超声风、温等观测资料对地面阵风特性进行了诊断分析.结果表明,在风速时间序列中叠加有周期为3-7 min的阵风扰动,显现出明显的相干结构,即沿顺风方向阵风风速峰期有下沉运动,谷期有上升运动;阵风扰动的各向异性特征明显,沿顺风方向的阵风扰动能量最大,其次是沿侧向和垂直方向的扰动能量;沿顺风方向的阵风垂直动量通量向下传播,而沿侧风方向阵风扰动动量垂直通量总体贡献接近于0.阵风扰动沿顺风方向的积分空间尺度和时间尺度最大,沿侧风方向和垂直方向其次,均明显大于湍流的积分空间和时间尺度.此外,阵风扰动的其他特征还包括:感热垂直通量极小;当平均风速较大时阵风风向变化幅度较小,而风速较小时阵风风向变化幅度则较大;动力学分析表明,阵风扰动主要表现出重力内波的一些特性.     

A New Scheme for the Simulation of Microscale Flow and Dispersion in Urban Areas by Coupling Large-Eddy Simulation with Mesoscale Models

A coupling scheme is proposed for the simulation of microscale flow and dispersion in which both the mesoscale field and small-scale turbulence are specified at the boundary of a microscale model. The small-scale turbulence is obtained individually in the inner and outer layers by the transformation of pre-computed databases, and then combined in a weighted sum. Validation of the results of a flow over a cluster of model buildings shows that the inner- and outer-layer transition height should be located in the roughness sublayer. Both the new scheme and the previous scheme are applied in the simulation of the flow over the central business district of Oklahoma City (a point source during intensive observation period 3 of the Joint Urban 2003 experimental campaign), with results showing that the wind speed is well predicted in the canopy layer. Compared with the previous scheme, the new scheme improves the prediction of the wind direction and turbulent kinetic energy (TKE) in the canopy layer. The flow field influences the scalar plume in two ways, i.e. the averaged flow field determines the advective flux and the TKE field determines the turbulent flux. Thus, the mean, root-mean-square and maximum of the concentration agree better with the observations with the new scheme. These results indicate that the new scheme is an effective means of simulating the complex flow and dispersion in urban canopies.     

On the generation of sea waves by a turbulent wind

Miles' (1957) theory of wave generation is modified to include the effects of air turbulence. It is shown that the fluctuating wave-induced Reynolds stresses are important only within a critical layer surrounding the critical height where the mean wind speed matches the wave phase speed. When the critical height is such that the critical layer extends to the air-sea interface, the energy transfer from the wind to the wave is affected greatly by the turbulence.     

Local advection of momentum,heat, and moisture in micrometeorology

The local advection of momentum, heat and moisture in micrometeorology due to a horizontal inhomogeneity in surface conditions is numerically investigated by a higher-order turbulence closure model which includes equations for the mean quantities, turbulent fluxes, and the viscous dissipation rate. The application of the two-dimensional model in this paper deals with the simulation of the flow from an extensive smooth dry area to a grassy wet terrain. The mean wind speed, temperature, and humidity distributions in the resulting internal boundary layer downstream of the surface discontinuity are determined such that the energy and moisture balances at the Earth's surface are satisfied.Numerical calculations of the mean temperature and humidity profiles are compared with available observed ones. The results include the advective effects on turbulent flux distributions, surface energy balance, evaporation rate, and Bowen ratio. The sensitivity of the predicted mean profiles and turbulent flux distributions to the surface relative humidity, thermal stratification, and the roughness change is discussed.NRC-NAS Resident Research Associate at AFCRL.     

A wind-tunnel investigation of the wind speed and turbulence characteristics close to the ground over various escarpment shapes

A wind-tunnel investigation of the wind flow over two-dimensional forward-facing escarpments is reported as part of a continuing research programme into the effects of local topography on the wind flow close to the ground. Four sharp-edged escarpments with their slopes varying between a cliff and a 4:1 gradient, were placed normal to a simulated neutrally-stable rural boundary layer which was modelled to a scale of 1:300. The resulting flows close to the surfaces of the escarpments were measured with a hot-wire anemometer. The modifications to the mean wind speed, turbulence intensity and energy spectra over the escarpments are described. The results indicate the extent and magnitude of the modification to the flow and suggest that significant changes in turbulence characteristics only occur in the wake region close behind the crest, where a shift of energy to higher frequencies is evident.     

Impact of waves on air-sea exchange of sensible heat and momentum

The impact of sea waves on sensible heat and momentum fluxes is described. The approach is based on the conservation of heat and momentum in the marine atmospheric surface layer. The experimental fact that the drag coefficient above the sea increases considerably with increasing wind speed, while the exchange coefficient for sensible heat (Stanton number) remains virtually independent of wind speed, is explained by a different balance of the turbulent and the wave-induced parts in the total fluxes of momentum and sensible heat.Organised motions induced by waves support the wave-induced stress which dominates the surface momentum flux. These organised motions do not contribute to the vertical flux of heat. The heat flux above waves is determined, in part, by the influence of waves upon the turbulence diffusivity.The turbulence diffusivity is altered by waves in an indirect way. The wave-induced stress dominates the surface flux and decays rapidly with height. Therefore the turbulent stress above waves is no longer constant with height. That changes the balance of the turbulent kinetic energy and of the dissipation rate and, hence the diffusivity.The dependence of the exchange coefficient for heat on wind speed is usually parameterized in terms of a constant Stanton number. However, an increase of the exchange coefficient with wind speed is not ruled out by field measurements and could be parametrized in terms of a constant temperature roughness length. Because of the large scatter, field data do not allow us to establish the actual dependence. The exchange coefficient for sensible heat, calculated from the model, is virtually independent of wind speed in the range of 3–10 ms^-1. For wind speeds above 10 ms^-1 an increase of 10% is obtained, which is smaller than that following from the constant roughness length parameterization.The investigation was in part supported by the Netherlands Geosciences Foundation (GOA) with financial aid from the Netherlands Organization for Scientific Research (NWO).     

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...     

Turbulence Structure in the Wake Region of a Meteorological Tower

A meteorological tower significantly modifies the air flow, the mean windspeed and wind direction as well as the turbulencestructure of the air. Suchchanges can be noticed in particular in the wake region of the tower.Measurementson the 200 m tower ofForschungszentrum Karlsruhewere carried outusing Solent sonic anemometers in the lee of the towerand cup anemometers on both sides.In the wake region, spectral energydensity is increased in the high-frequency range. Superposition of this disturbance spectrum on the undisturbedspectrum yields a `knee' in the resulting spectrum. In the case of low turbulence intensity with stable stratification,a plateau with a constant energy content is observed in front of the knee.This effect is caused by the new production of turbulence energy from the mean flow as well as by an energy transfer fromlarger to smaller vortices. Power spectra in strongly stable conditionsshow a more rapid decrease of intensity in the region where the inertialsubrange is expected.The relevant scales of wake turbulence are derived from the maximum of the disturbance spectrum.Locations of the high-frequency peak do not depend on atmospheric stability,but are controlled mainly by mean wind speed.Apart from the reduction of the mean wind speed, the spectra and cospectra exhibit a strong anisotropy for such cases.The results demonstrate the significant influence of a tower on turbulence spectra in the wake region.     

Analysis on the Interaction between Turbulence and Secondary Circulation of the Surface Layer at Jinta Oasis in Summer

The kinetic energy variations of mean flow and turbulence at three levels in the surface layer were calculated by using eddy covariance data from observations at Jinta oasis in 2005 summer.It is found that when the mean horizontal flow was stronger,the turbulent kinetic energy was increased at all levels,as well as the downward mean wind at the middle level.Since the mean vertical flow on the top and bottom were both negligible at that time,there was a secondary circulation with convergence in the upper half and divergence in the lower half of the column.After consideration of energy conversion,it was found that the interaction between turbulence and the secondary circulation caused the intensification of each other.The interaction reflected positive feedback between turbulence and the vertical shear of the mean flow.Turbulent sensible and latent heat flux anomaly were also analyzed.The results show that in both daytime and at night,when the surface layer turbulence was intensified as a result of strengthened mean flow,the sensible heat flux was decreased while the latent heat flux was increased.Both anomalous fluxes contributed to the cold island effect and the moisture island effect of the oasis.     

A Simple Method of Estimating Scalar Fluxes Over Forests

A simple aerodynamic-variance method is proposed to fill gaps in continuous CO₂ flux measurements in rainy conditions, when open-path analysers do not function. The method requires turbulent conditions (friction velocity greater than 0.1 ms^–1), and uses measurements of mean wind speed, and standard deviations of temperature and CO₂ concentration fluctuations to complement, and at times replace, eddy-covariance measurements of friction velocity, sensible heat flux and CO₂ flux. Friction velocity is estimated from the mean wind speed with a flux-gradient relationship modified for the roughness sublayer. Since normalised standard deviations do not follow Monin-Obukhov similarity theory in the roughness sublayer, a simple classification scheme according to the scalar turbulence scale was used. This scheme is shown to produce sensible heat and CO₂ flux estimates that are well correlated with the measured values.     

The influence of nonstationarity on the turbulent flux–gradient relationship for stable stratification

Extensive eddy-correlation datasets are analyzed to examine the influence of nonstationarity of the mean flow on the flux–gradient relationship near the surface. This nonstationarity is due to wavelike motions, meandering of the wind vector, and numerous unidentified small-scale mesoscale motions. While the data do not reveal an obvious critical gradient Richardson number, the maximum downward heat flux increases approximately linearly with increasing friction velocity for significant stability. The largest of our datasets is chosen to more closely examine the influence of stability, nonstationarity, distortion of the mean wind profile and self-correlation on the flux-gradient relationship. Stability is expressed in terms of z/L, the gradient Richardson number or the bulk Richardson number over the tower layer. The efficiency of the momentum transport systematically increases with increasing nonstationarity and attendant distortion of the mean wind profile. Enhancement of the turbulent momentum flux associated with nonstationarity is examined in terms of the nondimensional shear, Prandtl number and the eddy diffusivity.