Surface Temperature and Surface-Layer Turbulence in a Convective Boundary Layer

Previous laboratory and atmospheric experiments have shown that turbulence influences the surface temperature in a convective boundary layer. The main objective of this study is to examine land-atmosphere coupled heat transport mechanism for different stability conditions. High frequency infrared imagery and sonic anemometer measurements were obtained during the boundary layer late afternoon and sunset turbulence (BLLAST) experimental campaign. Temporal turbulence data in the surface-layer are then analyzed jointly with spatial surface-temperature imagery. The surface-temperature structures (identified using surface-temperature fluctuations) are strongly linked to atmospheric turbulence as manifested in several findings. The surface-temperature coherent structures move at an advection speed similar to the upper surface-layer or mixed-layer wind speed, with a decreasing trend with increase in stability. Also, with increasing instability the streamwise surface-temperature structure size decreases and the structures become more circular. The sequencing of surface- and air-temperature patterns is further examined through conditional averaging. Surface heating causes the initiation of warm ejection events followed by cold sweep events that result in surface cooling. The ejection events occur about 25 % of the time, but account for 60–70 % of the total sensible heat flux and cause fluctuations of up to 30 % in the ground heat flux. Cross-correlation analysis between air and surface temperature confirms the validity of a scalar footprint model.     

Surface-Layer Characteristics in the Stable Boundary Layer with Strong and Weak Winds

An extensive meteorological dataset obtained from the plumevalidation experiment conducted by the Electric Power Research Institute (EPRI) atKincaid during 1980–1981 is analysed for studying the characteristic differences in thesurface-layer parameters in strong and weak wind stable conditions. The surface-layerparameters are computed using the similarity functions _m and _h proposed byBeljaars and Holtslag. The weak winds are characterized using the geostrophic wind speedas well as the wind speed at the 10-m level. The surface fluxes are found to be finitein weak wind conditions.Empirical formulations for the eddy diffusivities of momentum(K_M) and heat (K_H), and drag (C_D) and heat exchange (C_H) coefficients, as powerlaw functions of the bulk Richardson number (Ri_B), are proposed under both strong andweak wind conditions. Results are close to those based on observations taken from the IndianInstitute of Technology low wind diffusion experiment, the Land surface processes experiment,the Hanford diffusion experiment, the Cabauw field experiment and the Cooperative Atmospheric SurfaceExchange Study 1999 (CASES-99) experiment. In addition, the fluxes obtained fromthe proposed empirical relations are in good agreement with those based on similarity theory as wellas the turbulence measurements taken from the CASES-99 experiment.     

Intermittency of Turbulence in the Atmospheric Boundary Layer: Scaling Exponents and Stratification Influence

We show the relationship between the intermittency of turbulence and the type of stratification for different atmospheric situations during the SABLES98 field campaign. With this objective, we first demonstrate the scaling behaviour of the velocity structure functions corresponding to these situations; next, we analyze the curvature of the scaling exponents of the velocity structure functions versus the order of these functions (ζ _p vs. p), where ζ _p are the exponents of the power relation for the velocity structure function with respect to the scale. It can be proved that this curve must be concave, under the assumption that the incompressible approximation does not break down at high Reynolds numbers. The physical significance of this kind of curvature is that the energy dissipation rate increases as the scale of the turbulent eddies diminishes (intermittency in the usual sense). However, the constraints imposed by stability, preventing full development of the turbulence, allow the function ζ _p versus p to show any type of curvature. In this case, waves of high frequency trapped by the stability, or bursts of turbulence caused by the breaking up of internal waves, may produce a redistribution of energy throughout the scaling range. Due to this redistribution, the variation with the scale of the energy dissipation rate may be smaller (decreasing the intermittency) and, even in more stable situations, this rate may diminish (instead of increasing) as the scale diminishes (convex form of the curve ζ _p vs. p).     

Scaling Properties of Temperature Spectra and Heat-Flux Cospectra in the Surface Friction Layer Beneath an Unstable Outer Layer

Temperature variance and temperature power spectra in the unstable surface layer have always presented a problem to the standard Monin-Obukhov similarity model. Recently that problem has intensified with the demonstration by Smedman et al. (2007, Q J Roy Meteorol Soc 133: 37–51) that temperature spectra and heat-flux cospectra can have two distinct peaks in slightly unstable conditions, and by McNaughton et al. (2007, Nonlinear Process Geophys 14: 257–271) who showed that the wavenumber of the peak of temperature spectra in a convective boundary layer (CBL), closely above the surface friction layer (SFL), can be sensitive to the CBL depth, z _i. Neither the two-peak form at slight instability nor the dependence of peak position on z _i at large instability is compatible with the Monin-Obukhov model. Here we examine the properties of temperature spectra and heat-flux cospectra from between these extremes, i.e. from within the unstable SFL, in two experiments. The analysis is based on McNaughton’s model of the turbulence structure in the SFL. According to this model, heat is transported through most of the SFL by sheet plumes, created by the action of impinging outer eddies. The smallest and most effective of these outer eddies have sizes that scale on SFL depth, z _s. The z _s-scale eddies and plumes are organised within the overall convection pattern in the CBL, and in turn they organise the motion of smaller eddies within the SFL, whose sizes scale on height, z. The main experimental results are: (1) the peak amplitudes of the temperature spectra in the SFL are collapsed with a scaling factor (z_sz)^1/3e_o^2/3{(z_{\rm s}z)^{1/3}\varepsilon_{\rm o}^{2/3}} divided by the square of the surface temperature flux, where e_o{\varepsilon_{\rm o}} is the dissipation rate of turbulent energy in the outer CBL (above the SFL); (2) the peak wavenumbers of the temperature spectra are collapsed with the mixed length scale (z _i z _s)^1/2; (3) the peak wavenumbers of the heat-flux cospectra are collapsed with the doubly-mixed length scale (z _i z _s)^1/4 z ^1/2; (4) for z/z _s < 0.03, the peak in the cospectrum is replaced by another peak at a wavenumber about a magnitude larger. This peak’s position scales on z; (5) all these findings are consistent with the observations of Smedman et al.     

Resolution Sensitivity and Scaling of Large-Eddy Simulations of the Stable Boundary Layer

Large-eddy simulations (LES) of the continuously turbulent quasi-equilibrium stable boundary layer (SBL) are conducted with grid lengths in the range of 12.5 m to 2 m, in order to explore resolution sensitivity, and determine at what point grid convergence occurs. The structure of the mean potential temperature, winds, and turbulent fluxes varies significantly over this resolution range. The highest resolution simulations show a significant degree of convergence. The dimensionless momentum diffusivity asymptotes to a value of 0.06, corresponding to a limiting flux Richardson number of 0.15.Using the converged simulations, some scaling hypotheses underpinning first-order and second-order closure models are revisited. The effective Richardson number stability functions of the LES are compared with the forms often used in numerical weather prediction (NWP). The mixing implied by the LES is less than that used in NWP. The commonly used similarity profiles for heat and momentum fluxes, and the scalings for dissipation and pressure covariances are compared with the LES. This information could provide guidance for the next generation of SBL parametrization schemes.     

Vertical Profiles of Turbulence Parameters in the Thermal Internal Boundary Layer

Boundary-Layer Meteorology - The correct simulation of pollutant dispersion in coastal regions demands understanding the turbulence structure of the thermal internal boundary layer (TIBL), which...     

Research Progress on Estimation of the Atmospheric Boundary Layer Height

Atmospheric boundary layer height(ABLH) is an important parameter used to depict characteristics of the planetary boundary layer(PBL) in the lower troposphere. The ABLH is strongly associated with the vertical distributions of heat, mass, and energy in the PBL, and it is a key quantity in numerical simulation of the PBL and plays an essential role in atmospheric environmental assessment. In this paper, various definitions and methods for deriving and estimating the ABLH are summarized, from the perspectives of turbulent motion, PBL dynamics and thermodynamics, and distributions of various substances in the PBL. Different methods for determining the ABLH by means of direct observation and remote sensing retrieval are reviewed, and comparisons of the advantages and disadvantages of these methods are presented. The paper also summarizes the ABLH parameterization schemes, discusses current problems in the estimation of ABLH, and finally points out the directions for possible future breakthroughs in the ABLHrelated research and application.     

Impact of the Nocturnal Low-Level Jet and Orographic Waves on Turbulent Motions and Energy Fluxes in the Lower Atmospheric Boundary Layer

The nocturnal low-level jet (LLJ) and orographic (gravity) waves play an important role in the generation of turbulence and pollutant dispersion and can affect the energy production by wind turbines. Additionally, gravity waves have an influence on the local mixing and turbulence within the surface layer and the vertical flux of mass into the lower atmosphere. On 25 September 2017, during a field campaign, a persistent easterly LLJ and gravity waves were observed simultaneously in a coastal area in the north of France. We explore the variability of the wind speed, turbulent eddies, and turbulence kinetic energy in the time–frequency and space domain using an ultrasonic anemometer and a scanning wind lidar. The results reveal a significant enhancement of the turbulence-kinetic-energy dissipation (by?50%) due to gravity waves in the LLJ shear layer (below the jet core) during the period of wave propagation. Large magnitudes of zonal and vertical components of the shear stress (approximately 0.4 and 1.5 m² s^?2, respectively) are found during that period. Large eddies (scales of 110 to 280 m) matching the high-wind-speed regime are found to propagate the momentum downwards, which enhances the mass transport from the LLJ shear layer to the roughness layer. Furthermore, these large-scale eddies are associated with the crests while comparatively small-scale eddies are associated with the troughs of the gravity wave.

     

Derivation of Structure Parameters of Temperature and Humidity in the Convective Boundary Layer from Large-Eddy Simulations and Implications for the Interpretation of Scintillometer Observations

We derive the turbulent structure parameters of temperature $C_{T}^2$ and humidity $C_q^2$ from high-resolution large-eddy simulations (LES) of a homogeneously-heated convective boundary layer. Boundary conditions and model forcing were derived from measurements at Cabauw in The Netherlands. Three different methods to obtain the structure-parameters from LES are investigated. The shape of the vertical structure-parameter profiles from all three methods compare well with former experimental and LES results. Depending on the method, deviations in the magnitude up to a factor of two are found and traced back to the effects of discretization and numerical dissipation of the advection scheme. Furthermore, we validate the LES data with airborne and large-aperture scintillometer (LAS) measurements at Cabauw. Virtual path measurements are used to study the variability of $C_{T}^2$ in the mixed layer and surface layer and its implications for airborne and LAS measurements. A high variability of $C_{T}^2$ along a given horizontal path in the LES data is associated with plumes (high values) and downdrafts (low values). The path average of $C_{T}^2$ varies rapidly in time due to the limited path length. The LES results suggest that measured path averages require sufficient temporal averaging and an adequate ratio of path length to height above the ground for the LAS in order to approach the domain average of $C_{T}^2$ .     

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.     

Techniques for the Eduction of Coherent Structures from Flow Measurements in the Atmospheric Boundary Layer

Two empirical methods to detect coherent motions embedded in the flow field have been compared, namely the variable interval time average (VITA) method and a wavelet-based technique, both with artificial signals as well as velocity measurements from the atmospheric boundary layer over a forest canopy. It has been found that the wavelet method is slightly better than the VITA approach in coherent structure eduction, even if the results of both techniques are comparable. Also the application of the present approach to simultaneous conditionally sampled wind data has highlighted some important features of coherent structures and gust generation in canopy flows.     

Estimation Of Atmospheric Water Vapour Flux Profiles In The Nocturnal Unstable Urban Boundary Layer With Doppler Sodar And Raman Lidar

The vertical wind profiles determined by Doppler sodar and the water vapourmixing ratio profiles obtained by Raman lidar are used to estimate the atmosphericwater vapour flux profiles in the nocturnal urban boundary layer under unstableconditions. The experiment was conducted for several nights in the central areaof Rome under a variety of moisture conditions and different urban boundary-layerflow regimes. Despite some scatter in the profiles, the latent heat flux is found tobe positive throughout the depth of the nocturnal urban boundary-layer. Thelayer-averaged flux shows a variation between -4 to +40 W m^-2, whileindividual values of flux in excess of +150 W m^-2 pertain to a case offree convection during cold air advection caused by the sea breeze. The qualityof flux estimates is found to be highly limited by the low sampling rates employedin the experiment resulting in errors to the order of 60%. Therefore, the results mustbe viewed as estimates rather than precise measurements. The skewness profiles ofthe turbulent fluctuations of vertical velocity and water vapour mixing ratio are alsopositive.     

Radiative Processes in the Stable Boundary Layer: Part II. The Development of the Nocturnal Boundary Layer

The interaction between radiation and turbulence in the stable boundary layer over land is explored using an idealized model, with a focus on the surface layer after the evening transition. It is shown that finer vertical resolution is required in transitional boundary layers than in developed ones. In very light winds radiative cooling determines the temperature profile, even if similarity functions without a critical Richardson number are used; standard surface similarity theory applied over thick layers then yields poor forecasts of near-surface air temperatures. These points are illustrated with field data. Simulations of the developing nocturnal boundary layer are used to explore the wider role of radiation. Comparatively, radiation is less significant within the developed stable boundary layer than during the transition; although, as previous studies have found, it remains important towards the top of the stable layer and in the residual layer. Near the ground, reducing the surface emissivity below one is found to yield modest relative radiative warming rather than intense cooling, which reduces the potential importance of radiation in the developed surface layer. The profile of the radiative heating rate may be strongly dependent on other processes, leading to quite varied behaviour.     

Scaling of the Vertical Spreading of a Plume of a Passive Tracer in a Rough-Wall Neutral Boundary Layer

The influence of surface roughness on the dispersion of a passive scalar in a rough wall turbulent boundary layer has been studied using wind-tunnel experiments. The surface roughness was varied using different sizes of roughness elements, and different spacings between the elements. Vertical profiles of average concentration were measured at different distances downwind of the source, and the vertical spread of the plume was computed by fitting a double Gaussian profile to the data. An estimate of the integral length scale is derived from the turbulence characteristics of the boundary layer and is then used to scale the measured values of plume spread. This scaling reduces the variability in the data, confirming the validity of the model for the Lagrangian integral time scale, but does not remove it entirely. The scaled plume spreading shows significant differences from predictions of theoretical models both in the near and in the far field. In the region immediately downwind of the source this is due to the influence of the wake of the injector for which we have developed a simple model. In the far field we explain that the differences are mainly due to the absence of large-scale motions. Finally, further downwind of the source the scaled values of plume spread fall into two distinct groups. It is suggested that the difference between the two groups may be related to the lack of dynamical similarity between the boundary-layer flows for varying surface roughness or to biased estimates of the plume spread.     

Variability of the Structure Parameters of Temperature and Humidity Observed in the Atmospheric Surface Layer Under Unstable Conditions

The structure parameters of temperature and humidity are important in scintillometry as they determine the structure parameter of the refractive index of air, the primary atmospheric variable obtained with scintillometers. In this study, we investigate the variability of the logarithm of the Monin-Obukhov-scaled structure parameters (denoted as $\log ({\widetilde{C_{s}^2}_{\mathrm {}}})$ ) of temperature and humidity. We use observations from eddy-covariance systems operated at three heights (2.5, 50, and 90 m) within the atmospheric surface layer under unstable conditions. The variability of $\log ({\widetilde{C_{s}^2}_{\mathrm {}}})$ depends on instability and on the size of the averaging window over which $\log ({\widetilde{C_{s}^2}_{\mathrm {}}})$ is calculated. If instability increases, differences in $\log ({\widetilde{C_{s}^2}_{\mathrm {}}})$ between upward motions (large $C_{s}^2$ ) and downward motions (small $C_{s}^2$ ) increase. The differences are, however, not sufficiently large to result in a bimodal probability density function. If the averaging window size increases, the variances of $\log ({\widetilde{C_{s}^2}_{\mathrm {}}})$ decrease. A linear regression of the variances of $\log ({\widetilde{C_{s}^2}_{\mathrm {}}})$ versus the averaging window size for various stability classes shows an increase of both the offset and slope (in absolute sense) with increasing instability. For temperature, data from the three heights show comparable results. For humidity, in contrast, the offset and slope are larger at 50 and 90 m than at 2.5 m. In the end we discuss how these findings could be used to assess whether observed differences in $C_{s}^2$ along a scintillometer path or aircraft flight leg are just within the range of local variability in $C_{s}^2$ or could be attributed to surface heterogeneity. This is important for the interpretation of data measured above a heterogeneous surface.     

一次重污染过程及其边界层气象特征量分析

针对2013年12月14~25日出现的区域性重污染过程,采用Hysplit后推气团轨迹模式分析了污染形成源,利用污染中心邢台的探空、地面数据计算了大气稳定度、混合层高度、逆温等气象特征量,并对混合层高度、相对湿度、能见度与PM2.5浓度进行了相关分析。结果表明:此次重污染以局地排放为主要形成源,河北省中南部地区大气层结偏稳定,逆温层厚(平均230 m)、强度大(平均2.34℃/100 m)、混合层高度低(平均618 m)是影响污染程度的重要因素;PM2.5浓度与混合层高度呈现负相关(R=-0.80),与能见度呈指数相关(R=-0.77),与相对湿度呈弱的正相关(R=0.62)。     

Analytical Solution for the Convectively-Mixed Atmospheric Boundary Layer

Based on the prognostic equations of mixed-layer theory assuming a zeroth order jump at the entrainment zone, analytical solutions for the boundary-layer height evolution are derived with different degrees of accuracy. First, an exact implicit expression for the boundary-layer height for a situation without moisture is analytically derived without assuming any additional relationships or specific initial conditions. It is shown that to expand the solution to include moisture, only minor approximations have to be made. Second, for relatively large boundary-layer heights, the implicit representation is simplified to an explicit function. Third, a hybrid expression is proposed as a reasonable representation for the boundary-layer height evolution during the entire day. Subsequently, the analysis is extended to present the evolution of any boundary-layer averaged scalar, either inert or under idealized chemistry, as an analytical function of time and boundary-layer height. Finally, the analytical solutions are evaluated. This evaluation includes a sensitivity analysis of the boundary-layer height for the entrainment ratio, the free tropospheric lapse rate of the potential temperature, the time-integrated surface flux and the initial boundary-layer height and potential temperature jump.     

Simulating Dispersion in the Evening-Transition Boundary Layer

We investigate dispersion in the evening-transition boundary layer using large-eddy simulation (LES). In the LES, a particle model traces pollutant paths using a combination of the resolved flow velocities and a random displacement model to represent subgrid-scale motions. The LES is forced with both a sudden switch-off of the surface heat flux and also a more gradual observed evolution. The LES shows ‘lofting’ of plumes from near-surface releases in the pre-transition convective boundary layer; it also shows the subsequent ‘trapping’ of releases in the post-transition near-surface stable boundary layer and residual layer above. Given the paucity of observations for pollution dispersion in evening transitions, the LES proves a useful reference. We then use the LES to test and improve a one-dimensional Lagrangian Stochastic Model (LSM) such as is often used in practical dispersion studies. The LSM used here includes both time-varying and skewed turbulence statistics. It is forced with the vertical velocity variance, skewness and dissipation from the LES for particle releases at various heights and times in the evening transition. The LSM plume spreads are significantly larger than those from the LES in the post-transition stable boundary-layer trapping regime. The forcing from the LES was thus insufficient to constrain the plume evolution, and inclusion of the significant stratification effects was required. In the so-called modified LSM, a correction to the vertical velocity variance was included to represent the effect of stable stratification and the consequent presence of wave-like motions. The modified LSM shows improved trapping of particles in the post-transition stable boundary layer.     

An Analytical Slab Model for the Growth of the Coastal Thermal Internal Boundary Layer Under Near-Neutral Onshore Flow Conditions

Analytical parameterisations of the thermal internal boundary-layer (TIBL) height based on the slab approach are widely used in coastal dispersion models. However, they tend to a singular behaviour when the stability of the onshore flow is close to neutral. Assuming that convective turbulence dominates mixing, we derive a more general analytical model that is valid for both stable and neutral onshore flows. The model is based on the existing framework for the slab approach but involving the Zilitinkevich correction (or the spin-up term). The height variation of the onshore flow lapse rate is accounted for in the model by including an initial TIBL height. An algebraic form of the model also includes the mechanical mixing contribution to the TIBL growth and is, therefore, suitable for use when the overland surface heat flux is small and friction velocity large. The new analytical model is tested with field measurements taken under near-neutral onshore flow conditions. The performance of the model is shown to be better than a commonly used TIBL parameterisation scheme.