Numerical Study of the Response of an Atmospheric Surface Layer to a Spatially Nonuniform Plant Canopy

High-accuracy large-eddy simulations of neutral atmospheric surface-layer flow over a gapped plant canopy strip have been performed. Subgrid-scale (SGS) motions are parameterized by the Sagaut mixed length SGS model, with a modification to compute the SGS characteristic length self-adaptively. Shaw’s plant canopy model, taking the vertical variation of leaf area density into account, is applied to study the response of the atmospheric surface layer to the gapped dense forest strip. Differences in the region far away from the gap and in the middle of the gap are investigated, according to the instantaneous velocity magnitude, the zero-plane displacement, the potential temperature and the streamlines. The large-scale vortex structure, in the form of a roll vortex, is revealed in the region far away from the gap. The nonuniform spatial distribution of plants appears to cause the formation of the coherent structure. The roll vortex starts in the wake of the canopy, and results in strong fluctuations throughout the entire canopy region. Wind sweeps and ejections in the plant canopy are also attributed to the large vortex structure.     

Time Scales in the Unstable Atmospheric Surface Layer

Calculation of eddy covariances in the atmospheric surface layer (ASL) requires separating the instantaneous signal into mean and fluctuating components. Since the ASL is not statistically stationary, an inherent ambiguity exists in defining the mean quantities. The present study compares four methods of calculating physically relevant time scales in the unstable ASL that may be used to remove the unsteady mean components of instantaneous time signals, in order to yield local turbulent fluxes that appear to be statistically stationary. The four mean-removal time scales are: (t _c ) based on the location of the maximum in the ogive of the heat flux cospectra, () the location of the zero crossing in the multiresolution decomposition of the heat flux, (t ^*) the ratio of the mixed-layer depth over the convective velocity, and () the convergence time of the vertical velocity and temperature variances. The four time scales are evaluated using high quality, three-dimensional sonic anemometry data acquired at the Surface Layer Turbulence and Environmental Science Test (SLTEST) facility located on the salt flats of Utah’s western desert. Results indicate that and , with t _c achieving values about 2–3 times greater than t ^*. The sensitivity of the eddy covariances to the mean-removal time scale (given a fixed 4-h averaging period during midday) is also demonstrated.     

Large-Eddy Simulation of the Stable Atmospheric Boundary Layer using Dynamic Models with Different Averaging Schemes

Large-eddy simulation (LES) of a stable atmospheric boundary layer is performed using recently developed dynamic subgrid-scale (SGS) models. These models not only calculate the Smagorinsky coefficient and SGS Prandtl number dynamically based on the smallest resolved motions in the flow, they also allow for scale dependence of those coefficients. This dynamic calculation requires statistical averaging for numerical stability. Here, we evaluate three commonly used averaging schemes in stable atmospheric boundary-layer simulations: averaging over horizontal planes, over adjacent grid points, and following fluid particle trajectories. Particular attention is focused on assessing the effect of the different averaging methods on resolved flow statistics and SGS model coefficients. Our results indicate that averaging schemes that allow the coefficients to fluctuate locally give results that are in better agreement with boundary-layer similarity theory and previous LES studies. Even among models that are local, the averaging method is found to affect model coefficient probability density function distributions and turbulent spectra of the resolved velocity and temperature fields. Overall, averaging along fluid pathlines is found to produce the best combination of self consistent model coefficients, first- and second-order flow statistics and insensitivity to grid resolution.     

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.     

Geometric Alignments of the Subgrid-Scale Force in the Atmospheric Boundary Layer

In recent years field experiments have been undertaken in the lower atmosphere to perform a priori tests of subgrid-scale (SGS) models for large-eddy simulations (LES). The experimental arrangements and data collected have facilitated studies of variables such as the filtered strain rate, SGS stress and dissipation, and the eddy viscosity coefficient. However, the experimental set-ups did not permit analysis of the divergence of the SGS stress (the SGS force vector), which is the term that enters directly in the LES momentum balance equations. Data from a field experiment (SGS2002) in the west desert of Utah, allows the calculation of the SGS force due to the unique 4 × 4 sonic anemometer array. The vector alignment of the SGS force is investigated under a range of atmospheric stabilities. The eddy viscosity model is likely aligned with the measured SGS force under near-neutral and unstable conditions, while its performance is unsatisfactory under stable conditions.     

Measurements and Modelling of the Wind Speed Profile in the Marine Atmospheric Boundary Layer

We present measurements from 2006 of the marine wind speed profile at a site located 18 km from the west coast of Denmark in the North Sea. Measurements from mast-mounted cup anemometers up to a height of 45 m are extended to 161 m using LiDAR observations. Atmospheric turbulent flux measurements performed in 2004 with a sonic anemometer are compared to a bulk Richardson number formulation of the atmospheric stability. This is used to classify the LiDAR/cup wind speed profiles into atmospheric stability classes. The observations are compared to a simplified model for the wind speed profile that accounts for the effect of the boundary-layer height. For unstable and neutral atmospheric conditions the boundary-layer height could be neglected, whereas for stable conditions it is comparable to the measuring heights and therefore essential to include. It is interesting to note that, although it is derived from a different physical approach, the simplified wind speed profile conforms to the traditional expressions of the surface layer when the effect of the boundary-layer height is neglected.     

Large-eddy Simulation of the Dispersion of Solid Particles in a Turbulent Boundary Layer

A large-eddy simulation (LES) with the dynamic Smagorinsky-Germano subgrid-scale (SGS) model is used to study the dispersion of solid particles in a turbulent boundary layer. Solid particles are tracked in a Lagrangian way. The instantaneous velocity of the surrounding fluid is considered to have a large-scale part (directly computed by the LES) and a small-scale part. The SGS velocity of the surrounding fluid is given by a three-dimensional Langevin model written in terms of SGS statistics at a mesh level. An appropriate Lagrangian correlation time scale is considered in order to include the influences of gravity and inertia of the solid particle. Inter-particle collisions and the influence of particles on the mean flow are also taken into account. The results of the LES are compared with the wind-tunnel experiments of Nalpanis et al. (1993 J Fluid Mech 251: 661–685) and of Tanière et al. (1997 Exp in Fluids 23:463–471) on sand particles in saltation and in modified saltation, respectively.     

An Investigation on Temperature Variance Scaling in the Atmospheric Surface Layer

We present surface-layer measurements of temperature fluctuation variance from a site characterized by small-scale inhomogeneities. Periods of marked radiative forcing are selected. The data characterized by diabatic conditions and vertical heat flux larger than some threshold (here, chosen to be 0.01 K ms⁻¹) agree quite well with the convective scaling in unstable cases, and with the z-less parameterisation (with a large scatter) in stable cases. For near-neutral cases, the similarity function diverges because of the loss of significance of the temperature scale. Departures from similarity are highlighted in cases with smaller thermal fluxes, because horizontal heterogeneity and unsteadiness become important as production terms.     

Scaling Range Exponents From X-Wire Measurements In The Atmospheric Surface Layer

This paper focuses on the behaviour of moments, up to order 6, of longitudinal and vertical velocity increments, measured in the atmospheric surface layer, at a height of 1.7 m. The local derivatives of these moments with respect to the spatial separation between two points indicate that inertial range power-law exponents cannot be determined unambiguously. This is supported by the local slopes of moments of the locally averaged energy dissipation rate but contrasts with the extensive power-law ranges indicated by spectra of longitudinal and vertical velocity fluctuations. The third-order longitudinal velocity structure function provides stronger evidence of anisotropy in the inertial range than either the second-order velocity structure functions or the velocity spectra.     

Statistics of Scalar Fields in the Atmospheric Boundary Layer Based on Large-eddy Simulations. Part 1: Free Convection

Convection in a quasi-steady, cloud-free, shear-free atmospheric boundary layer is investigated based on a large-eddy simulation model. The performed tests indicate that the characteristic (peak) values of statistical moments at the top of the mixed layer are proportional to the interfacial scales (from gradients of scalars in the interfacial layer). Based on this finding a parameterization is proposed for profiles of scalar variances. The parameterization employs two, semi-empirical similarity functions F_m(z/z_i) andF_i(z/z_i), multiplied by a combination of the mixed-layer scales and the interfacial scales.     

Near-Surface Coherent Structures and The Vertical Momentum Flux in a Large-Eddy Simulation of the Neutrally-Stratified Boundary Layer

The near-surface flow of a well-resolved large-eddy simulation of the neutrally-stratified planetary boundary layer is used to explore the relationships between coherent structures and the vertical momentum flux. The near-surface flow is characterized by transient streaks, which are alternating bands of relatively higher and lower speed flow that form parallel to the mean shear direction in the lower part of the boundary layer. Although individual streaks are transient, the overall flow is in a quasi-equilibrium state in which the streaks form, grow, decay and regenerate over lifetimes on the order of tens of minutes. Coupled with the streaky flow is an overturning circulation with alternating bands of updrafts and downdrafts approximately centered on the streaks. The surface stress is dominated by upward ejections of slower moving near-surface air and downward sweeps of higher speed air from higher in the boundary layer. Conditional sampling of the ejection and sweep events shows that they are compact, coherent structures and are intimately related to the streaks: ejections (sweeps) preferentially form in the updrafts (downdrafts) of the three-dimensional streak flow. Hence, consistent with other recent studies, we propose that the streak motion plays an important role in the maintenance of the surface stress by establishing the preferential conditions for the ejections and sweeps that dominate the surface stress. The velocity fluctuation spectra in the model near the surface have a k ⁻¹ spectral slope over an intermediate range of wavenumbers. This behaviour is consistent with recent theoretical predictions that attempt to evaluate the effects of organized flow, such as near-surface streaks, on the variance spectra.     

水、陆不均匀条件下大气边界层结构的模拟研究

应用2004年白洋淀野外观测实验资料,对白洋淀地区水、陆不均匀分布条件下的大气边界层开展了数值模拟试验,结果表明:由于下垫面的热力差异,在360多km2的白洋淀地区,地表特征分布的不均匀可以引发弱的局地环流,影响大气边界层内温度和湿度的空间分布;在500m高度以下,下垫面不均匀分布对湍流动能有明显影响.最后将模拟结果与观测资料进行了对比检验.     

An Improvement of the Louis Scheme for the Surface Layer in an Atmospheric Modelling System

The Louis scheme for calculating the vertical eddy fluxes within the atmospheric surface layer is improved by broadening the original assumptions. In our approach, the momentum and heat transfer roughness lengths (z₀ and z_T respectively) can be different, and z₀ need not be negligibly small compared with the lowest height (z) in modelling. For these conditions, we choose more consistent wind and potential temperature profile forms, then derive new algorithms for calculating fluxes. Improvement is demonstrated for a wide range of z/L (L is the Obukhov length), z/z₀ and z₀ z_T, by comparing these fluxes with those derived from a theoretical surface-layer model. The improved algorithms can be used in atmospheric modelling systems for more varied surfaces and a wide range of atmospheric stability.     

Coherent Structures in the Atmospheric Surface Layer under Stable and Unstable Conditions

Simultaneous measurements of the instantaneous values of absolute temperatureat seven heights within the lower 36 m of the atmospheric boundary layer underdifferent stability conditions were carried out, accompanied by measurements ofthe wind velocity components at two levels and of solar radiation flux at the surface.The data obtained allow one to investigate individual convective cells known ascoherent structures (CS). Outside the CS, i.e., during quiet periods, an instanttemperature profile is in close agreement with the dry-adiabatic lapse rate, butwithin CS the temperature changes much faster with height, and the shape ofthe profile varies significantly.A method was developed to transform temperature records from sensors atseveral heights into an other form, namely, into temporal variations of theheights of isothermal surfaces. Since coherent structures were found to advectwith the mean wind velocity, these temporal height variations may be transformedinto the spatial ones, i.e., into the xoz-plane section of the temperature field.In such a pictorial presentation coherent structures look like asymmetric columnsof heat, penetrating the whole atmospheric surface layer.Coherent structures also exist in the stable stratified surface layer, but they have aninverse asymmetry and occupy only the lower several metres. Wavelike activitydominates in the upper part of the stable surface layer.The characteristic time of surface-layer adjustment to the rapid changes of solarradiation (due to cloud shadows or cloud gaps) was found to be on the order ofone minute. Such a time interval is required for coherent structure to reach the topof surface layer.     

Numerical Simulation of Drag Partition over Rough Surfaces

We present a numerical simulation of drag partition over rough surfaces. A computational fluid dynamics model is applied with high resolution to simulatingturbulent flows over arrays of roughness elements positioned on asmooth surface. The skin drag on the surface and the pressure drag on the roughnesselements are computed. The simulated drag partition compares well with wind-tunnelmeasurements and theoretical estimates for similar rough surfaces. This confirms that the computational approach offers an alternative to wind-tunnel and field experiments in studying drag and drag partition. The model is then applied to studying drag partition over rough surfaces with various roughness configurations. It is shown that drag partition depends not only on the magnitude of the roughness frontal area but also on the sizes and arrangement of roughness elements, because (1) the pressure drag coefficient is sensitive to roughness-element dimensions and (2) the arragement of roughness elements lead to different interferences of turbulent wakes. The impact ofthe latter factor is not insignificant.     

The Coherent Structure of Water Vapour Transfer in the Unstable Atmospheric Surface Layer

Observations of water vapour fluctuations over arice field show vapour ramps. Coherent structuresare first revealed by the frequently occurring ramp pattern in the vapourtrace. Wavelet and pseudo-wavelet analysis techniques were used inconditional sampling, and more than 100 hr of data have been analyzedto determine coherent structure characteristics. The most probablecoherent structure duration was in the range 2–12 sec andthe duration range of the most effective coherent structures shows somedifference between heat and water vapour transfers. Coherent structurescontribute to the major part of the total flux.     

Statistics of Scalar Fields in the Atmospheric Boundary Layer Based on Large-Eddy Simulations. Part II: Forced Convection

Forced convection in a quasi-steady atmospheric boundary layer is investigated based on a large-eddy simulation (LES) model. The performed simulations show that in the upper portion of the mixed layer the dimensionless (in terms of mixed layer scales) vertical gradients of temperature, humidity, and wind velocity depend on the dimensionless height z/z _i and the Reech number Rn. The peak values of variances and covariances at the top of the mixed layer, scaled in terms of the interfacial scales, are functions of the interfacial Richardson number Ri. As a result expressions for the entrainment rates, in the case when the interfacial layer has a finite depth, and a condition for the presence of moistening or drying regimes in the mixed layer, are derived. Profiles of dimensionless scalar moments in the mixed layer are proposed to be expressed in terms of two empirical similarity functions F _m and F _i , dependent on dimensionless height z/z _i , and the interfacial Richardson number Ri. The obtained similarity expressions adequately approximate the LES profiles of scalar statistics, and properly represent the impact of stability, shear, and entrainment. They are also consistent with the parameterization proposed for free convection in the first part of this paper.     

On Modelling the One-Dimensional Atmospheric Boundary Layer

Several commonly used turbulence closure schemes forthe atmospheric boundarylayer (ABL) are applied to simulate neutral, nocturnal and diurnal cyclesituations in a one-dimensional ABL. Results obtained with the differentschemes, E-, E- and its modified versions, and twoversions ofthe q ² Level 2.5, are compared and discussed.     

Characterization of Coherent structures in the Atmospheric Surface Layer

The ramplike coherent structures, observed in the temporal series of temperature and humidity in the atmospheric surface layer, are analyzed using the intermittency function and the wavelet transforms, with Haar, D4 and Mexican Hat functions as mother wavelets, in order to find the most efficient conditional sampling technique. It was found that the intermittency function and the wavelet transform, using Mexican Hat as mother wavelet, are the only ones that sample structures that fulfill the ramplike coherent structures definition of a slow rise followed by a sudden drop in the temporal series. The conditionally averaged structures detected by both techniques were similar for temperature, humidity, and vertical velocity at heights of 3, 5, and 9.4 m. Significant discrepancies were found among the conditional averaged structures detected by both techniques for zonal and meridional components of the wind at 11.5 m. Considering both techniques, it was observed that the averagedcoherent-structure duration ranged from 23.7 ± 0.5 s to 37.8 ± 3.0 s. Furthermore, the averaged number of events per 20-minute period ranged from 20.0 ± 1.0 to 28.5 ± 1.1, and the averaged intermittency factor from 45.0 ± 0.4% to 59.1 ± 1.3%. It was also observed that the averaged duration of the ramplike coherent structures increases with height, while their intensity, number, and intermittency factor decrease. Despite the good matching obtained for temperature and humidity, the coherent-structure properties did not show the expected variation with wind speed, stability parameter, and friction velocity. The Kolmogorov–Smirnov test indicated that the intermittent function and the wavelet transform did not detect coherent structures belonging to the same population.     

The Adaptation of the Atmospheric Boundary Layer to a Change in Surface Roughness

The adaptation of the atmospheric boundary layer to a change in the underlying surface roughness is an interesting problem and hence much research, theoretical, experimental, and numerical, has been undertaken. Within the atmospheric boundary layer an accurate numerical model for the turbulent properties of the atmospheric boundary layer needs to be implemented if physically realistic results are to be obtained. Here, the adaptation of the atmospheric boundary layer to a change in surface roughness is investigated using a first-order turbulence closure model, a one-and-a-half-order turbulence closure model and a second-order turbulence closure model. Perturbations to the geostrophic wind and the pressure gradients are included and it is shown that the second-order turbulence closure model, namely the standard k - model, is inferior to a lower-order closure model if a modification to limit the turbulent eddy size within the atmospheric boundary layer is not included within the model.