Wind-Tunnel and Numerical Simulations of the Coastal Thermal Internal Boundary Layer

Wind-tunnel experiments in a thermally stratified wind tunnel and direct numerical simulations were performed to simulate the thermal internal boundary layer (TIBL) that developed over a coastal area in a sea-breeze flow. The results of the simulations were analyzed to investigate turbulence structure in the TIBL. To study the effects of the atmospheric stability over the sea on the TIBL, two vertical profiles of temperature were created in the upstream portion of the wind-tunnel experiment and the direct numerical simulation. Turbulence statistics of the TIBL changed significantly according to the temperature profile over the sea, indicating that the stability of the flow over the sea has a significant effect on the structure and turbulence characteristics of the TIBL. Furthermore, the TIBL heights were estimated from the vertical profiles of the local Richardson number. The estimated TIBL heights agreed with those predicted by a pre-existing relation, suggesting that both the wind-tunnel experiment and the direct numerical simulation accurately reproduced the growth of the TIBL.     

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

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

Laboratory Simulation Of Vertical Plume Dispersion Within A Convective Boundary Layer

Mesoscale measurements of the vertical dispersion coefficient ₂ by using a composite turbulence water tank were validated through a comparison with CONDORS (Convective Diffusion Observed with Remote Sensors) field data, and were analysed with respect to the intensity of the thermal flux, mechanical turbulence, and plume release height.It seems possible to correct the plume _z values for different release heights below 0.5z_i (z_i is the mixing height) by applying an equation expressing the height dependency of turbulence intensity. The downwind distance where the plume's mass centre height approaches its final level was also analysed with respect to the above three parameters, and an empirical equation to estimate the downwinddistance derived.     

Improving Non-local Parameterization of the Convective Boundary Layer

We propose improvements in the “non-local” parameterization scheme of the convective boundary layer. The countergradient terms for components of the momentum fluxes are introduced in a form analogous to those for other scalars. The scheme also includes explicit expressions for entrainment fluxes of momentum, temperature, and humidity. A simplified procedure for calculating the boundary-layer height is proposed, consisting of two steps: the evaluation of the convection level, followed by the assessment of the depth of the interfacial layer.     

对流边界层中过山气流的数值模拟

采用ARPS4.0非静力中尺度气象模式模拟了对流边界层中气流过山引起的地形波,讨论了地形及大气条件改变对其的影响.模拟表明,当大气边界层是对流边界层时,气流过山引起的地形强迫,仍能在上部稳定层结中造成足够的垂直扰动,产生向上传播的重力内波,重力内波引起的波动阻力仍不可忽略.     

大气对流边界层发展的模拟研究

室内水槽模拟是大气边界层研究的一种重要手段。利用室内模拟水槽对大气边界层的发展进行了模拟,通过处理平均温度廓线和光斑图像得到了对流边界层顶部位置h2和边界层高度zi。结果表明,不同测量方法得到的结果一致性很好,与实际大气的边界层发展情况也较为接近。同时,根据试验情况确定初始条件和边界条件,使用边界层参数化模型进行了数值模拟,其结果与室内模拟的结果也较吻合。     

The Neutral, Barotropic Planetary Boundary Layer, Capped by a Low-Level Inversion

The presence of a low-level, capping inversion layer will affect the height and structure of the planetary boundary layer (PBL). Results from models of varying levels of sophistication, including analytical, turbulent kinetic energy (TKE), second-order closure (SOC), large-eddy simulation (LES) and direct numerical simulation (DNS) models, are used to investigate this influence for the neutral, barotropic PBL. Predicted and observed profiles of stress and geostrophic departure components, and integral measures, such as the parameters of Rossby-number similarity theory, are compared for the KONTUR, Marine Stratocumulus, JASIN, Leipzig, Pre-Wangara and Upavon field experiments.Analytical models of the equilibrium value of inversion height z_i, which depend on the surface friction velocity u_*, and both the Coriolis parameter f and the free-flow Brunt-Väisälä frequency N, are found to give reasonable estimates of the PBL height. They also indicate that only the KONTUR and Marine Stratocumulus experiments were strongly influenced by N. More quantitative comparisons would require larger, more comprehensive datasets. The effects of the presence of a capping inversion on the profile structure were found to be insignificant for h_* = |f|z_i/u_* > 0.15.The simple analytical model performed quite well over all values of h_*; it predicted the profiles of the longitudinal stress component (in the direction of the surface stress) better than the lateral component. The more advanced models performed well for small values of h_* (for flow over the sea), but systematically underestimated the cross-isobaric angle for flow over land. These models predicted the profiles of the lateral stress component better than the longitudinal component. The profiles of the analytical model agreed with those of the advanced models when the constant eddy viscosity of the outer layer was increased.Agreement with DNS was achieved by increasing the eddyviscosity of the analytical model by a factor of 5.Zilitinkevich and Esau(2002, Boundary-Layer Meteorology, 104, 371–379)suggest that the neutral, barotropic values of A and B of Rossby-numbersimilarity theory are not universal constants, but depend on the ratio N/|f|. The dependence for A and B is calculated using the analytical model and TKE models. Over the sea (h_* 0.1; N/|f| 100, where we have used the Zilitinkevich-Esau relation to convert between h_* and N/|f|) there is agreement between the model predictions and observations; however over land where the equilibrium boundary-layer height is greater (h_* 0.35; N/|f| 10) the inconsistency between the advanced model predictions (TKE, SOC, LES, and DNS) and observations, as noted previously by Hess and Garratt, still exists. We attribute this disagreement to violations of the strict assumptions of steady, horizontally homogeneous, neutral, barotropic conditions implicit in the observations. At small values of z_i and a strongly stable background stratification (h_* 0.04; N/|f| 1000) both the TKE and analytical models predict that A and B depend significantly on h_*, however observations are unavailable to confirm these predictions. Zilitinkevich and Esau call this case the `long-lived near-neutral PBL', and state that it is found in cold weather at high latitudes.     

Discussion and Applications of Slab Models of the Convective Boundary Layer based on Turbulent Kinetic Energy Budget Parameterisations

Some of the most widely used slab model formulations for applications in the convective boundary layer are analysed and discussed. Three main classes are identified based on different approximations of the turbulent kinetic energy equation. The models appear to be quite insensitive to the initial values for boundary-layer height, and temperature discontinuity at the boundary-layer top. The slab models are applied to a case of sea-land transition from the literature, and a case of convective boundary layer time evolution over a homogeneous terrain at San Pietro Capofiume (Bologna, Italy). The different parameterisations turn out to be almost equivalent for the cases studied. The models generally underpredict the value for the height, while all give very good estimates for the mean mixed-layer temperature.     

A Variable Mesh Spacing for Large-Eddy Simulation Models in the Convective Boundary Layer

A variable vertical mesh spacing for large-eddy simulation (LES) models in a convective boundary layer (CBL) is proposed. The argument is based on the fact that in the vertical direction the turbulence near the surface in a CBL is inhomogeneous and therefore the subfilter-scale effects depend on the relative location between the spectral peak of the vertical velocity and the filter cut-off wavelength. From the physical point of view, this lack of homogeneity makes the vertical mesh spacing the principal length scale and, as a consequence, the LES filter cut-off wavenumber is expressed in terms of this characteristic length scale. Assuming that the inertial subrange initial frequency is equal to the LES filter cut-off frequency and employing fitting expressions that describe the observed convective turbulent energy one-dimensional spectra, it is feasible to derive a relation to calculate the variable vertical mesh spacing. The incorporation of this variable vertical grid within a LES model shows that both the mean quantities (and their gradients) and the turbulent statistics quantities are well described near to the ground level, where the LES predictions are known to be a challenging task.     

An Observational Study of Wind Profiles in the Baroclinic Convective Mixed Layer

A comprehensive planetary boundary-layer (PBL) and synoptic data set is used to isolate the mechanisms that determine the vertical shear of the horizontal wind in the convective mixed layer. To do this, we compare a fair-weather convective PBL with no vertical shear through the mixed layer (10 March 1992), with a day with substantial vertical shear in the north-south wind component (27 February). The approach involves evaluating the terms of the budget equations for the two components of the vertical shear of the horizontal wind; namely: the time-rate-of-change or time-tendency term, differential advection, the Coriolis terms (a thermal wind term and a shear term), and the second derivative of the vertical transport of horizontal momentum with respect to height (turbulent-transport term). The data, gathered during the 1992 STorm-scale Operational and Research Meteorology (STORM) Fronts Experiments Systems Test (FEST) field experiment, are from gust-probe aircraft horizontal legs and soundings, 915-MHz wind profilers, a 5-cm Doppler radar, radiosondes, and surface Portable Automated Mesonet (PAM) stations in a roughly 50 × 50 km boundary-layer array in north-eastern Kansas, nested in a mesoscale-to-synoptic array of radiosondes and surface data.We present evidence that the shear on 27 February is related to the rapid growth of the convective boundary layer. Computing the shear budget over a fixed depth (the final depth of the mixed layer), we find that the time-tendency term dominates, reflecting entrainment of high-shear air from above the boundary layer. We suggest that shear within the mixed layer occurs when the time-tendency term is sufficiently large that the shear-reduction terms – namely the turbulent-transport term and differential advection terms – cannot compensate. In contrast, the tendency term is small for the slowly-growing PBL of 10 March, resulting in a balance between the Coriolis terms and the turbulent-transport term. Thus, the thermal wind appears to influence mixed-layer shear only indirectly, through its role in determining the entrained shear.     

Simulated Airborne Flux Measurements in a LES generated Convective Boundary Layer

One aim of past boundary-layer experiments with aircraft was the determination of areally averaged heat fluxes. In spite ofsophisticated instrumentation the measured fluxes extrapolated to the ground differed significantly from fluxes measured directly at ground stations. This studypresents simulated sensible heat flux measurements with aircraft flightsthrough a synthetic convective boundary layer created by a401 × 401 × 42 cubic-grid large eddy simulation (LES) with agrid spacing of 50 m. After some considerations with respect to necessary measurement lengths using results ofLenschow and Stankov (1986 – J. Atmos. Sci. 43, 1198–1209), simulated measurementcampaigns were carried out in three modelruns. During each model run five sets ofmeasurement runs were carried out successively.During each set of runs 10 aircraftflew at 10 altitudes with a ground speedof 100 m s^-1 simultaneously throughtime and space. In total, 150 legs were carried out, 15 at each flight level. The resulting`measured' heat fluxes were compared withthose of the `true' flux profiles obtaineddirectly from the ensemble-averagedLES-generated data. No significant systematic error between `measured' and `true' profiles was observed. Furthermore, the comparison of the resulting relative error with the theory ofLenschow and Stankov showed a good agreement at allmeasurement levels.     

Parameterization of Entrainment in a Sheared Convective Boundary Layer Using a First-order Jump Model

Basic entrainment equations applicable to the sheared convective boundary layer (CBL) are derived by assuming an inversion layer with a finite depth, i.e., the first-order jump model. Large-eddy simulation data are used to determine the constants involved in the parameterizations of the entrainment equations. Based on the integrated turbulent kinetic energy budget from surface to the top of the CBL, the resulting entrainment heat flux normalized by surface heat flux is a function of the inversion layer depth, the velocity jumps across the inversion layer, the friction velocity, and the convection velocity. The developed first-order jump model is tested against large-eddy simulation data of two independent cases with different inversion strengths. In both cases, the model reproduces quite reasonably the evolution of the CBL height, virtual potential temperature, and velocity components in the mixed layer and in the inversion layer.The part of this work was done when the first author visited at NCAR.     

On the Determination of the Neutral Drag Coefficient in the Convective Boundary Layer

Based on the idea that free convection can be considered as a particular case of forced convection, where the gusts driven by the large-scale eddies are scaled with the Deardorff convective velocity scale, a new formulation for the neutral drag coefficient, C_Dn, in the convective boundary layer (CBL) is derived. It is shown that (i) a concept of C_Dn can still be used under strongly unstable conditions including a pure free-convection regime even when no logarithmic portion in the velocity profile exists; (ii) gustiness corrections must be applied for rational calculations of C_Dn; and (iii) the stratification function used in the derivation of C_Dn should satisfy the theoretical free-convection limit. The new formulation is compared with the traditional relationship for C_Dn, and data collected over the sea (during the Tropical Ocean-Global Atmosphere Coupled Ocean-Atmosphere Response Experiment (TOGA COARE) and the San Clemente Ocean Probing Experiment (SCOPE)) and over land (during the BOREX-95 experiment) are used to illustrate the difference between the new and traditional formulations. Compared to the new approach, the traditional formulation strongly overestimates C_Dn and z_o in the CBL for mean wind speed less than about 2 m s^-1. The new approach also clarifies several contradictory results from earlier works. Some aspects related to an alternate definition of the neutral drag coefficient and the wind speed and the stress averaging procedure are considered.     

TRAC98: Detection of Coherent Structures in a Convective Boundary Layer using Airborne Measurements

The TRAC98 experimental campaign (Turbulence Radar Aircraft Cells) devoted to coherent structures analysis took place over the Beauce plain (France) during summer 1998. It allowed us to collect a large dataset of airborne measurements in addition to various ground measurements. This study aims at diagnosing the occurrence of coherent structures within the atmospheric boundary layer (ABL) through airborne measurements. The statistical analysis performed as a first step from turbulent parameters underlined the homogeneity of the ABL over the Beauce plain. However mixed-layer scaling failed at the top of the ABL, even when taking into account the entrainment rate. Coherent structures were detected through the analysis of ABL isotropy, using the opportunity of sampling with two perpendicular crossing planes, one of them being aligned with the wind. This approach allowed us to determine an organization scheme of the ABL for three of the five flights (ARAT30, MIV30 and MIV27). For the ARAT30 flight, the analysis was pursued by focusing on measurements of fluctuations in the inner flight legs. In this way, the low-level cloud cover has been investigated from the downward visible radiation (VISD). The results indicated an anisotropy of the horizontal cloud size. Secondly, the variations of some parameters were analysed through lagged correlation functions. This allowed us to infer relationships between the vertical velocity, VISD, mixing ratio and lifting condensation level. Length scales have also been extracted, and confirmed the ABL organization during the ARAT30 flight. Finally, the anisotropy observed in various flights has been investigated with respect to the underestimation of the latent heat fluxes revealed by the imbalance of the surface energy budget.     

Coherence and Scale of Vertical Velocity in the Convective Boundary Layer from a Doppler Lidar

We utilized a Doppler lidar to measure integral scale and coherence of vertical velocity w in the daytime convective boundary layer (CBL). The high resolution 2 μm wavelength Doppler lidar developed by the NOAA Environmental Technology Laboratory was used to detect the mean radial velocity of aerosol particles. It operated continuously in the zenith-pointing mode for several days in the summer 1996 during the “Lidars in Flat Terrain” experiment over level farmland in central Illinois. We calculated profiles of w integral scales in both the alongwind and vertical directions from about 390 m height to the CBL top. In the middle of the mixed layer we found, from the ratio of the w integral scale in the vertical to that in the horizontal direction, that the w eddies are squashed by a factor of about 0.65 as compared to what would be the case for isotropic turbulence. Furthermore, there is a significant decrease of the vertical integral scale with height. The integral scale profiles and vertical coherence show that vertical velocity fluctuations in the CBL have a predictable anisotropic structure. We found no significant tilt of the thermal structures with height in the middle part of the CBL.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

边界层参数化方案对暴雨数值模拟的影响

选取2003年7月4-5日南京暴雨个例，采用非静力中尺度模式MM5进行模拟，着重研究了不同边界层参数化方案对雨量中心强度、雨区分布的影响。结果表明：对于不同的边界层参数化方案，垂直速度场、水汽通量散度场、涡度场、水平风场的散度以及θse场都表现出不同的特征；合理边界层方案的引入对预报效果有明显的改进；结合边界层和自由大气的动力、热力结构进行了综合分析，给出了边界层作用与自由大气动力、热力结构的配置情况。说明这种配置对暴雨的形成是至关重要的。     

Equilibrium Evaporation and the Convective Boundary Layer

A theory is developed for surface energy exchanges in well-mixed, partlyopen systems, embracing fully open and fully closed systems as limits.Conservation equations for entropy and water vapour are converted intoan exact rate equation for the potential saturation deficit D in a well-mixed, partly open region. The main contributions to changes in D arise from (1) the flux of D at the surface, dependent on a conductance g_q that is a weighted sum of the bulk aerodynamic and surface conductances; and (2) the exchange flux of D with the external environment by entrainment or advection, dependent on a conductance g_e that is identifiable with the entrainment velocity when the partly open region is a growing convective boundary layer (CBL). The system is fully open when g_e/g_q , and fully closed when g_e/g_q 0. The equations determine the steady state surface energy balance (SEB) in a partly open system, the associated steady-state deficit, and the settling time scale needed to reach the steady state. The general result for the steady-state SEB corresponds to the equations of conventional combination theory for the SEB of a vegetated surface, with the surface-layer deficit replaced by the external deficit and with g_q replaced by the series sum (g_q ^-1 + g_e ^-1)^-1. In the fully open limit D is entirely externally prescribed, while in the fully closed limit, D is internally determined and the SEB approaches thermodynamic equilibrium energy partition. In the case of the CBL, the conductances g_q and g_e are themselves functions of D through short-term feedbacks, induced by entrainment in the case of g_e and by both physiological and aerodynamic (thermal stability) processes in the case of g_q. The effects of these feedbacks are evaluated. It is found that a steady-state CBL is physically achievable only over surfaces with at least moderate moisture availability; that entrainment has a significant accelerating effect on equilibration; that the settling time scale is well approximated by h/(g_q + g_e), where h is the CBL depth; and that this scale is short enough to allow a steady state to evolve within a semi-diurnal time scale only when h is around 500 m or less.     

Forward-in-time and Backward-in-time Dispersion in the Convective Boundary Layer: the Concentration Footprint

The turbulence field obtained using a large-eddy simulation model is used to simulate particle dispersion in the convective boundary layer with both forward-in-time and backward-in-time modes. A Lagrangian stochastic model is used to treat subgrid-scale turbulence. Results of forward dispersion match both laboratory experiments and previous numerical studies for different release heights in the convective boundary layer. Results obtained from backward dispersion show obvious asymmetry when directly compared to results from forward dispersion. However, a direct comparison of forward and backward dispersion has no apparent physical meaning and might be misleading. Results of backward dispersion can be interpreted as three-dimensional or generalized concentration footprints, which indicate that sources in the entire boundary layer, not only sources at the surface, may influence a concentration measurement at a point. Footprints at four source heights in the convective boundary layer corresponding to four receptors are derived using forward and backward dispersion methods. The agreement among footprints derived with forward and backward methods illustrates the equivalence between both approaches. The paper shows explicitly that Lagrangian simulations can yield identical footprints using forward and backward methods in horizontally homogeneous turbulence.     

Evidence Of Dynamical Coupling Between The Residual Layer And The Developing Convective Boundary Layer

The diurnal cycle of the atmospheric boundary layer (ABL) hasbeen documented on 8 August 1998 in the framework of the Étude et Simulation de la QUalité de l'air en Ile-de-France (ESQUIF) experiment that took place in the Paris area. The ABL structure was documented by means of a ground-based lidar, surface meteorological stations and soundings. The interaction between the residual layer and the convective boundary layer is investigated using the collected data as well as mesoscale modelling. As opposed to the generally accepted concept, we find evidence of entrainment at the top ofthe residual layer. High temporal simulations of the 8 August 1998 casemade with the mesoscale atmospheric model Meso-NH also show evidenceof mixing at the top of the residual layer (RL). This mixing is believed to be related to the presence of convective (gravity) waves in the RL.     

Application of Dynamic Subgrid-scale Models for Large-eddy Simulation of the Daytime Convective Boundary Layer over Heterogeneous Surfaces

The sensitivity of large-eddy simulation (LES) to the representation of subgrid-scale (SGS) processes is explored for the case of the convective boundary layer (CBL) developing over surfaces with varying degrees of spatial heterogeneity. Three representations of SGS processes are explored: the traditional constant Smagorinsky–Lilly model and two other dynamic models with Lagrangian averaging approaches to calculate the Smagorinsky coefficient (C _S ) and SGS Prandtl number (Pr). With initial data based roughly on the observed meteorology, simulations of daytime CBL growth are performed over surfaces with characteristics (i.e. fluxes and roughness) ranging from homogeneous, to striped heterogeneity, to a realistic representation of heterogeneity as derived from a recent field study. In both idealized tests and the realistic case, SGS sensitivities are mostly manifest near the surface and entrainment zone. However, unlike simulations over complex domains or under neutral or stable conditions, these differences for the CBL simulation, where large eddies dominate, are not significant enough to distinguish the performance of the different SGS models, irrespective of surface heterogeneity.