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.     

The effect of a regional increase in ocean surface roughness on the tropospheric circulation: a GCM experiment

The sensitivity of the atmospheric circulation to an increase in ocean surface roughness in the Southern Hemisphere storm track is investigated in a paired general circulation model experiment. Such a change in sea roughness could be induced by ocean waves generated by storms. Two extended permanent-July runs are made. One with standard sea surface roughness, the other with ten times as a large surface roughness over open sea poleward of 40° S. The regional increase in ocean surface roughness significantly modifies the tropospheric circulation in the Southern Hemisphere. The strongest effect is the reduction of tropospheric winds (by 2 m/s or 10%) above the area with increased roughness. The poleward eddy momentum flux is reduced in the upper troposphere and the meridional eddy sensible heat flux is reduced in the lower troposphere. Zonal mean and eddy kinetic energy are consistently reduced.     

Scaling Variances of Scalars in a Convective Boundary Layer Under Different Entrainment Regimes

For the presentation and analysis of atmospheric boundary-layer (ABL) data, scales are used to non-dimensionalise the observed quantities and independent variables. Usually, the ABL height, surface sensible heat flux and surface scalar flux are used. This works well, so long as the absolute values of the entrainment ratio for both the scalar and temperature are similar. The entrainment ratio for temperature naturally ranges from −0.4 to −0.1. However, the entrainment ratio for passive scalars can vary widely in magnitude and sign. Then the entrainment flux becomes relevant as well. The only customary scalar scale that takes into account both the surface flux and the entrainment flux is the bulk scalar scale, but this scale is not well-behaved for large negative entrainment ratios and for an entrainment ratio equal to −1. We derive a new scalar scale, using previously published large-eddy simulation results for the convective ABL. The scale is derived under the constraint that scaled scalar variance profiles are similar at those heights where the variance producing mechanisms are identical (i.e., either near the entrainment layer or near the surface). The new scale takes into account that scalar variance in the ABL is not only related to the surface flux of that scalar, but to the scalar entrainment flux as well. Furthermore, it takes into account that the production of variance by the entrainment flux is an order of magnitude larger than the production of variance by the surface flux (per unit flux). Other desirable features of the new scale are that it is always positive (which is relevant when scaling standard deviations) and that the scaled variances are always of order 1–10.     

Perturbations to the Spatial and Temporal Characteristics of the Diurnally-Varying Atmospheric Boundary Layer Due to an Extensive Wind Farm

The effect of extensive terrestrial wind farms on the spatio-temporal structure of the diurnally-evolving atmospheric boundary layer is explored. High-resolution large-eddy simulations of a realistic diurnal cycle with an embedded wind farm are performed. Simulations are forced by a constant geostrophic velocity with time-varying surface boundary conditions derived from a selected period of the CASES-99 field campaign. Through analysis of the bulk statistics of the flow as a function of height and time, it is shown that extensive wind farms shift the inertial oscillations and the associated nocturnal low-level jet vertically upwards by approximately 200 m; cause a three times stronger stratification between the surface and the rotor-disk region, and as a consequence, delay the formation and growth of the convective boundary layer (CBL) by approximately 2 h. These perturbations are shown to have a direct impact on the potential power output of an extensive wind farm with the displacement of the low-level jet causing lower power output during the night as compared to the day. The low-power regime at night is shown to persist for almost 2 h beyond the morning transition due to the reduced growth of the CBL. It is shown that the wind farm induces a deeper entrainment region with greater entrainment fluxes. Finally, it is found that the diurnally-averaged effective roughness length for wind farms is much lower than the reference value computed theoretically for neutral conditions.     

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

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

峨眉山地区近地层微气象特征研究

基于2019年12月至2020年11月峨眉山站梯度塔资料、辐射观测资料和地表通量资料,采用涡动相关法对峨眉山地区近地层的地表通量和蒸散发量的变化进行分析,并估算了零平面位移、空气动力粗糙度、空气热力粗糙度、动量通量输送系数和感热通量输送系数等重要的空气动力学和热力学参数.研究表明:近地面风速呈现高层高、低层低的特征,且...     

干旱区绿洲诱发的中尺度运动的模拟及其关键因子的敏感性实验

利用一个非静力平衡的、高分辨的、二维中尺度大气数值模式，并在仅考虑简单过程的情况下，模拟了干旱区中绿洲所诱发的中尺度运动，并进行了这种中尺度大气运动的强度对绿洲水平尺度、绿洲与周围环境的水平热力差异、大尺度背景场水平风速和大尺度地表加热率等一些重要物理参数关系的敏感性实验研究。研究发现：中尺度大气运动强度随水平热力差异的增大而加强，随背景场水平风速和大尺度地表加热率增强而分别减弱。但随绿洲水平尺度的变化并不像前三个因子一样为单调函数，而是在绿洲水平尺度为20km时中尺度大气运动最强，绿洲水平尺度更大或更小时中尺度大气运动强度均会减弱。通过统计甘肃省河西地区的绿洲水平尺度分布规律，发现绿洲分布最集中的尺度在15－25km，与模拟所得到的能源发最强中尺度运动的绿洲水平尺度基本一致。     

Wind profiles,momentum fluxes and roughness lengths at Cabauw revisited

We describe the results of an experiment focusing on wind speed and momentum fluxes in the atmospheric boundary layer up to 200 m. The measurements were conducted in 1996 at the Cabauw site in the Netherlands. Momentum fluxes are measured using the K-Gill Propeller Vane. Estimates of the roughness length are derived using various techniques from the wind speed and flux measurements, and the observed differences are explained by considering the source area of the meteorological parameters. A clear rough-to-smooth transition is found in the wind speed profiles at Cabauw. The internal boundary layer reaches the lowest k-vane (20 m) only in the south-west direction where the obstacle-free fetch is about 2 km. The internal boundary layer is also reflected in the roughness lengths derived from the wind speed profiles. The lower part of the profile (< 40 m) is not in equilibrium and no reliable roughness analysis can be given. The upper part of the profile can be linked to a large-scale roughness length. Roughness lengths derived from the horizontal wind speed variance and gustiness have large footprints and therefore represent a large-scale average roughness. The drag coefficient is more locally determined but still represents a large-scale roughness length when it is measured above the local internal boundary layer. The roughness length at inhomogeneous sites can therefore be determined best from drag coefficient measurements just above the local internal boundary layers directly, or indirectly from horizontal wind speed variance or gustiness. In addition, the momentum and heat fluxes along the tower are analysed and these show significant variation with height related to stability and possibly surface heterogeneity. It appears that the dimensionless wind speed gradients scale well with local fluxes for the variety of conditions considered, including the unstable cases.     

Scalar Fluxes from Urban Street Canyons Part II: Model

A practical model is developed for the vertical flux of a scalar, such as heat, from an urban street canyon that accounts for variations of the flow and turbulence with canyon geometry. The model gives the magnitude and geometric dependence of the flux from each facet of the urban street canyon, and is shown to agree well with wind-tunnel measurements described in Part I. The geometric dependence of the flux from an urban street canyon is shown to be determined by two physical processes. Firstly, as the height-to-width ratio of the street canyon increases, so does the roughness length and displacement height of the surface. This increase leads to a reduction in the wind speed in the inertial sublayer above the street canyons. Since the speed of the circulations in the street are proportional to this inertial sublayer wind speed, the flux then reduces with the inertial sublayer wind speed. This process is dominant at low height-to-width ratios. Secondly, the character of the circulations within the street canyon also varies as the height-to-width ratio increases. The flow in the street is partitioned into a recirculation region and a ventilated region. When the street canyon has high height-to-width ratios the recirculation region occupies the whole street canyon and the wind speeds within the street are low. This tendency decreases the flux at high height-to-width ratios. These processes tend to reduce the flux density from the individual facets of the street canyon, when compared to the flux density from a horizontal surface of the same material. But the street canyon has an increased total surface area, which means that the total flux from the street canyon is larger than from a horizontal surface. The variations in scalar flux from an urban street canyon with geometry is over a factor of two, which means that the physical mechanisms responsible should be incorporated into energy balance models for urban areas.     

半干旱区榆中地表粗糙度年变化及影响机理

利用2006年6月—2010年12月兰州大学半干旱气候与环境观测站 (SACOL) 观测资料,分析了黄土高原自然植被下垫面榆中地表粗糙度时间变化特征,考虑到地形、植被物理特征以及降水和热力条件的影响,分析了东南风向和西北风向粗糙度年变化规律及其影响机理,并分别给出这两个风向归一化粗糙度与时间的拟合关系式。研究发现:对于非均一下垫面,由于地形起伏和下垫面植被差别造成的不同风向粗糙度差异显著。选取东南风向和西北风向,这两个风向的地表粗糙度无论是量级还是年变化特征都有很大差别,且由于地形和植被的差别,东南风向粗糙度年变化趋势与稳定度年变化趋势一致,粗糙度与稳定度存在一定相关关系,而西北风向粗糙度年变化趋势与降水量年变化趋势一致,粗糙度与降水量相关性较好。     

Mosaic方法在非均匀下垫面上的适用性研究

在非均匀下垫面情况下, Mosaic方法是目前国际上广泛运用于模式中计算地表通量的方法.大量的研究表明, 下垫面的非均匀分布会引发局地环流, 非均匀分布的空间尺度较大时, 所引起的环流甚至可以达到海陆风的强度.这种环流的存在直接影响到次网格地表通量的计算.次网格地表非均匀分布, 尤其是大尺度模式中的次网格非均匀分布, 必将影响地表通量的计算.本文针对次网格地表非均匀问题, 设计了高分辨率的Mosaic试验和非均匀试验, 开展了不同背景风情况下的一系列数值试验, 以探讨这种影响的程度.结果表明, 在土壤湿度空间分布不均匀的情况下, 运用Mosaic方法计算得到的地表潜热通量偏小, 背景风较小的时候偏差较大, 背景风增强时偏差减小.     

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.     

Comparative Studies In The Local Circulations Induced By Land-Use And By Topography

A two-dimensional numerical model witha simplified land-use paramaterizationis used to investigate the effects of land usetopography on local circulationsystems. A criterion is presented indicating therelative importance of land use onlocal circulations.Land-use contrast in the numerical model is parameterized by moisture availability and roughness length. Numerical experiments were carried out under various atmospheric stabilities with various dimensions for the mountains.Numerical results show the following: (1) Anabatic winds prevail in high mountains,while the land–land breeze is stronger when the horizontal contrast of the heat fluxfrom each land surface is large. (2) In the early morning, anabatic winds tend to prevailover land–land breezes, but the land–land breeze tends to prevail in the later afternoon.(3) While atmospheric stability has a large influence on the mesoscale circulation, thehorizontal scale of the mountain is not as important for the intensity of the mesoscalecirculation. (4) In strong stable conditions or weak insolation, the anabatic wind tendsto be more active than the land–land breeze; namely, the heat flux from inhomogeneousland use becomes a less important factor for the generation of a mesoscale disturbancethan the orographic forcing in the case of strong atmospheric stability. (5) The predominantmesoscale circulation is predicted by a criterion based upon the horizontal adjustmenttheory of the mixed layer. The criterion proposed in this study is based on the mountainheight, the ratio of the heat flux from the different land use patterns, the atmospheric staticstability, and the time-integrated heat flux. The criterion gives results that generally agreewith the numerical results.     

DERIVATION OF FLUXES FROM PROFILES OVER A MODERATELY HOMOGENEOUS FOREST

Measurements of fluxes and profiles of wind andtemperature are performed in the roughness layer ofa moderately homogeneous forest location. Weinvestigate to what extent vertical scalar fluxescan be derived from profile measurements. Theinfluence of inhomogeneities in the upwind terrainis investigated with footprint analysis and with aninhomogeneous surface-layer model. Four methods toestimate displacement height are suggested, amongthem is a method involving the structure parameterof the vertical wind. All methods give a decrease ofdisplacement height with increasing wind speed,while roughness length is found to increase withincreasing wind speed. For near-neutral conditionsdimensionless temperature gradients are found to besubstantially lower than the surface-layer valuesfound in the literature for homogeneous terrain with lowvegetation. Dimensionless shear however iscomparable with the surface-layer value. The heightof the roughness layer is 20 times the roughnesslength. Two schemes with locally derived surfaceparameters are tested to derive friction velocityand sensible heat flux from the profilemeasurements. These site specific schemes performsatisfactorily. A third scheme based on surface parameters chosen a priorifrom the literatureperforms significantly worse especially for low windspeed and unstable cases.     

A Scale-Dependent Dynamic Model for Scalar Transport in Large-Eddy Simulations of the Atmospheric Boundary Layer

An important challenge in large-eddy simulationsof the atmospheric boundarylayer is the specification of the subgrid-scale(SGS) model coefficient(s)and, in particular, how to account for factorssuch as position in the flow,grid/filter scale and atmospheric stability.A dynamic SGS model (thatassumes scale invariance of the coefficients)is implemented in simulationsof a neutral boundary layer with a constantand uniform surface flux of apassive scalar. Results from our simulationsshow evidence that the lumpedcoefficient in the eddy-diffusion modelcomputed with the dynamic proceduredepends on scale. This scale dependence isstronger near the surface, and itis more important for the scalar than for thevelocity field (Smagorinskycoefficient) due to the stronger anisotropicbehaviour of scalars. Based onthese results, a new scale-dependent dynamicmodel is developed for theeddy-diffusion lumped coefficient. The newmodel, which is similar to theone proposed earlierfor the Smagorinsky coefficient,is fully dynamic, thus not requiring anyparameter specification or tuning.Simulations with the scale-dependent dynamicmodel yield the expected trendsof the coefficients as functions of positionand filter/grid scale.Furthermore, in the surface layer the newmodel gives improved predictionsof mean profiles and turbulence spectra ascompared with the traditionalscale-invariant dynamic model.     

Flux Footprint Simulation Downwind of a Forest Edge

Surface fluxes, originating from forest patches, are commonly calculated from atmospheric flux measurements at some height above that patch using a correction for flux arising from upwind surfaces. Footprint models have been developed to calculate such a correction. These models commonly assume homogeneous turbulence, resulting in a simulated atmospheric flux equal to the average surface flux in the footprint area. However, atmospheric scalar fluxes downwind of a forest edge have been observed to exceed surface fluxes in the footprint area. Variations in atmospheric turbulence downwind of the forest edge, as simulated with an E – model, can explain enhanced atmospheric scalar fluxes. This E – model is used to calculate the footprint of atmospheric measurements downwind of a forest edge. Atmospheric fluxes appear mainly enhanced as a result of a stronger sensitivity to fluxes from the upwind surface. A sensitivity analysis shows that the fetch over forest, necessary to reach equilibrium between atmospheric fluxes and surface fluxes, tends to be longer for scalar fluxes as compared to momentum fluxes. With increasing forest density, atmospheric fluxes deviate even more strongly from surface fluxes, but over shorter fetches. It is concluded that scalar fluxes over forests are commonly affected by inhomogeneous turbulence over large fetches downwind of an edge. It is recommended to take horizontal variations in turbulence into account when the footprint is calculated for atmospheric flux measurements downwind of a forest edge. The spatially integrated footprint is recommended to describe the ratio between the atmospheric flux and the average surface flux in the footprint.     

Turbulent Flow Properties Around a Staggered Wind Farm

The fundamental properties of turbulent flow around a perfectly staggered wind farm are investigated in a wind tunnel. The wind farm consisted of a series of 10 rows by 2–3 columns of miniature wind turbines spaced 5 and 4 rotor diameters in the streamwise and spanwise directions respectively. It was placed in a boundary-layer flow developed over a smooth surface under thermally neutral conditions. Cross-wire anemometry was used to obtain high resolution measurements of streamwise and vertical velocity components at various locations within and above the wind farm. The results show that the staggered configuration is more efficient in terms of momentum transfer from the background flow to the turbines compared to the case of an aligned wind turbine array under similar turbine separations in the streamwise and spanwise directions. This leads to improved power output of the overall wind farm. A simplified analysis suggests that the difference in power output between the two configurations is on the order of 10%. The maximum levels of turbulence intensity in the staggered wind farm were found to be very similar to that observed in the wake of a single wind turbine, differing substantially with that observed in an aligned configuration with similar spacing. The dramatic changes in momentum and turbulence characteristics in the two configurations show the importance of turbine layout in engineering design. Lateral homogenization of the turbulence statistics above the wind farm allows for the development of simple parametrizations for the adjustment of flow properties, similar to the case of a surface roughness transition. The development of an internal boundary layer was observed at the upper edge of the wind farm within which the flow statistics are affected by the superposition of the ambient flow and the flow disturbance induced by the wind turbines. The adjustment of the flow in this layer is much slower in the staggered situation (with respect to its aligned counterpart), implying a change in the momentum/power available at turbine locations. Additionally, power spectra of the streamwise and vertical velocity components indicate that the signature of each turbine-tip vortex structure persists to locations deep within the wind farm.     

Modification of nocturnal drainage flow due to urban surface heat flux

Dense observations and numerical experiments were carried out to estimate the modification of mesoscale circulation, particularly cold drainage wind. It was confirmed that nocturnal drainage flow can develop on clear calm summer day and change due to orographical forcing and the heterogeneity of heat flux induced by the discontinuity of land-use. The temperature of nocturnal drainage flow at Sungji Valley, Busan Korea, tended to increase as it passed over the urban surface due to anthropogenic heat. The increase in temperature reached 2.9 K at night. The roughness associated with the exchange of momentum flux alone and the pass of nocturnal drainage flow is important for modifying the characteristics of flow Numerical simulations carried out under various surface conditions showed good agreement with the observations. Urban heat fluxes from the surface during the day are fundamental causes of the changes in the urban mesoscale circulation. In addition, the impact of a discontinuity of surface heat flux on mesoscale flow modification tends to be greater at night than during the day because the direction of urban surface heat fluxes at night is different from that in rural areas. In addition, the criterion to estimate the increase in temperature nocturnal drainage flow was also proposed, and provided results that generally agreed with the numerical results.     

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.     

Airflow from mustard to fallow

Real terrain is complex and its effects on the atmosphere are not well understood. The two‐dimensional change‐of‐roughness problem is investigated in an experimental study of airflow from mustard to fallow. The surface friction velocity is found to overshoot and return slowly to the equilibrium value. The wind profile in the modified region is more nearly logarithmic than the Peterson model predicts. The simple Elliot theory is in fair agreement with the wind profile data.

Values for the ratio of the standard deviation of the vertical wind to the friction velocity reveal no variation with stability and are in close agreement with the value for laboratory flows. Values for the ratio of the standard deviation of the longitudinal wind component to the friction velocity are comparable to other recently determined atmospheric values and show a variation attributable to large‐scale terrain properties.