Turbulence Characteristics of the Shear-Free Convective Boundary Layer Driven by Heterogeneous Surface Heating

Large-eddy simulations (LESs) are employed to investigate the turbulence characteristics in the shear-free convective boundary layer (CBL) driven by heterogeneous surface heating. The patterns of surface heating are arranged as a chessboard with two different surface heat fluxes in the neighbouring patches, and the heterogeneity scale Λ in four different cases is taken as 1.2, 2.5, 5.0 and 10.0 km, respectively. The results are compared with those for the homogeneous case. The impact of the heterogeneity scale on the domain-averaged CBL characteristics, such as the profiles of the potential temperature and the heat flux, is not significant. However, different turbulence characteristics are induced by different heterogeneous surface heating. The greatest turbulent kinetic energy (TKE) is produced in the case with the largest heterogeneity scale, whilst the TKE in the other heterogeneous cases is close to that for the homogeneous case. This result indicates that the TKE is not enhanced unless the scale of the heterogeneous surface heating is large enough. The potential temperature variance is enhanced more significantly by a larger surface heterogeneity scale. But this effect diminishes with increasing CBL height, which implies that the turbulent eddy structures are changed during the CBL development. Analyses show that there are two types of organized turbulent eddies: one relates to the thermal circulations induced by the heterogeneous surface heating, whilst the other identifies with the inherent turbulent eddies (large eddies) induced by the free convection. At the early stage of the CBL development, the dominant scale of the organized turbulent eddies is controlled by the scale of the surface heterogeneity. With time increasing, the original pattern breaks up, and the vertical velocity eventually displays horizontal structures similar to those for the homogeneous heating case. It is found that after this transition, the values of λ/z _i (λ is the dominant horizontal scale of the turbulent eddies, z _i is the boundary-layer height) ≈1.6, which is just the aspect ratio of large eddies in the CBL.     

A laboratory study of the turbulent velocity characteristics in the convective boundary layer

Based on the measurement of the velocity field in the convective boundary layer （CBL） in a convection water tank with the particle image velocimetry （PIV） technique, this paper studies the characteristics of the CBL turbulent velocity in a modified convection tank. The experiment results show that the velocity distribution in the mixed layer clearly possesses the characteristics of the CBL thermals, and the turbulent eddies can be seen obviously. The comparison of the vertical distribution of the turbulent velocity variables indicates that the modeling in the new tank is better than in the old one. The experiment data show that the thermal＇s motion in the entrainment zone sometimes fluctuates obviously due to the intermittence of turbulence. Analyses show that this fluctuation can influence the agreement of the measurement data with the parameterization scheme, in which the convective Richardson number is used to characterize the entrainment zone depth. The normalized square velocity wi^2/w＊^2. at the top of the mixed layer seems to be time-dependent, and has a decreasing trend during the experiments. This implies that the vertical turbulent velocity at the top of the mixed layer may not be proportional to the convective velocity （w＊）.     

涡旋相关法测定湍流通量偏低的研究

针对野外实验所发现的不同观测法测定地表能通量不平衡问题,进行了均匀加热大气边界层的大涡模拟实验.用模拟的湍流风、温度和湿度涨落的时间序列证实,对流边界层低频涡普遍存在,并经常以一簇一簇热泡的形式出现.风速较小时,有限时长的取样不足以捕捉低频涡的贡献,可造成涡旋相关法测量的统计量异常偏低.仿照涡旋相关法的步骤进行数据处理发现,经去除平均或趋势计算的温度和湿度通量偏低程度在边界层下部随观测高度的增高而显著,其中尤以湿度通量为甚.其结果在一定程度上可以解释低风速条件下地表能通量测量的不闭合问题,但是尚不能完全解释诸如青藏高原实验出现的严重不闭合.文中对此作了探讨性的讨论.     

Characteristics of vertical turbulent velocities in the urban convective boundary layer

Turbulence measurements of the vertical velocity component were obtained by an instrumented aircraft under fair weather conditions in the St. Louis, Missouri, metropolitan area. Time series of vertical velocity fluctuations from horizontal flight segments made in the lower part of and near the middle of the convective boundary layer (CBL) over the urban area and surrounding region were subjected to various statistical and objective analyses. Higher order vertical velocity moments, and positive and negative velocity statistics, were computed. The horizontal dimensions of updrafts and downdrafts, and related properties of these turbulent eddies were derived by conditional sampling analysis. Emphasis is on a comparison of the results from urban and selected rural measurements from the lower part of the CBL.The vertical velocity probability density distribution for each segment was positively skewed and the mode was negative. The means and standard deviations of positive and negative velocity fluctuations were greater over the urban area. The urban vertical velocity variance was 50% greater than rural values, and power spectra revealed greater production of vertical turbulent energy in the urban area over a wide frequency range.The mean and maximum widths of downdrafts were generally larger than the corresponding values for updrafts. Differences between urban and rural eddy sizes were not statistically significant. The widths of the largest updraft and downdraft are comparable to the boundary-layer depth, Z _i, and the mean value of the ratio of spectral peak wavelength to Z _iwas about 1.3 and 1.1 for urban and rural areas, respectively. Convective similarity scaling parameters appeared to order both the urban and rural measurements.On assignment from the National Oceanic and Atmospheric Administration, U.S. Dept. of Commerce.     

Aircraft Observations of the Development of Thermal Internal Boundary Layers and Scaling of the Convective Boundary Layer Over Non-Homogeneous Land Surfaces

This study investigates the convective boundary layer (CBL) that develops over anon-homogeneous surface under different thermal and dynamic conditions. Analysesare based on data obtained from a Russian research aircraft equipped with turbulentsensors during the GAME-Siberia experiment over Yakutsk in Siberia, from April to June 2000.Mesoscale thermal internal boundary layers (MTIBLs) that radically modified CBLdevelopment were observed under unstable atmospheric conditions. It was found thatMTIBLs strongly influenced the vertical and horizontal structures of virtual potentialtemperature, specific humidity and, most notably, the vertical sensible and latent heatfluxes. MTIBLs in the vicinity of the Lena River lowlands were confirmed by clouddistributions in satellite pictures.MTIBLs spread through the entire CBL and radically modify its structure if the CBL isunstable, and strong thermal features on the underlying surface have horizontal scalesexceeding 10 km. MTIBL detection is facilitated through the use of special parameterslinking shear stress and convective motion.The turbulent structure of the CBL with and without MTIBLs was scaled usingthe mosaic or flux aggregate approach. A non-dimensional parameterL_Rau/L_hetero (where L_Rau is Raupach's length and L_hetero is the horizontal scale of the surface heterogeneity)estimates the application limit of similarity and local similarity scaling models forthe mosaic parts over the surface. Normalized vertical profiles of wind speed, airtemperature, turbulent sensible and latent heat fluxes for the mosaic parts withL_RauL_hetero < 1 could be estimated by typical scalingcurves for the homogeneous CBL. Traditional similarity scaling models for the CBLcould not be applied for the mosaic parts with L_Rau/L_hetero > 1.For some horizontally non-homogeneous CBLs, horizontal sensible heat fluxes werecomparable with the vertical fluxes. The largest horizontal sensible heat fluxes occurred at the top of the surface layer and below the top of the CBL.Formerly affiliated to the Frontier Observational Research System for Global ChangeFormerly affiliated to the Frontier Observational Research System for Global Change     

日间对流边界层中的非局地动量混合

日间对流边界层最显著的结构特征是在热力作用下所形成的组织化对流。与小尺度湍涡不同的是，组织化对流具有边界层尺度的垂直相干性，可实现垂直贯穿边界层的非局地物质和能量传输。本文针对对流边界层中的动量混合，探究组织化对流对动量输送的贡献。以高精度大涡模拟数据为研究资料，通过傅里叶变换、本征正交分解和经验模态分解3种滤波方法，分离组织化对流和背景湍涡，计算与两者相关的非局地和局地动量通量，发现与组织化对流相关的非局地动量通量是总通量的重要组成部分，并主导混合层中的垂直动量输送。而后，基于协谱和相位谱分析，探究组织化对流的空间结构对动量传输的影响，发现在热力主导的不稳定环境中，单体型环流结构对动量的传输效率较低。而在风切较强的近中性环境中，滚涡型组织化结构可使垂直和水平流向扰动速度的相位差减小，从而提升动量传输效率。研究结果表明，边界层方案需要包含非局地动量通量项，其参数化应考虑整体稳定度对传输效率的影响。     

实验速度场测量技术及对流边界层特征研究

在对流槽中对对流边界层(CBL)温度场实验研究的基础上,进一步尝试通过实验技术测量速度场并分析研究CBL中的速度场特征。在应用PIV测量技术时选用铝粉作示踪粒子。实验证明了在混合层中速度分布明显具有对流边界层热泡特性;混合层顶部的速度分布很好地反应出夹卷层的结构特征;湍流速度特征量的垂直分布合理,与野外实测结果和类似的对流槽实验结果接近;误差分析表明示踪粒子的跟随性良好,粒子速度的测量结果能真实地反应流体的运动特征,从而得证了分析结果的可靠性。     

Resolved Versus Parametrized Boundary-Layer Plumes. Part I: A Parametrization-Oriented Conditional Sampling in Large-Eddy Simulations

A conditional sampling based on the combination of a passive tracer emitted at the surface and thermodynamic variables is proposed to characterise organized structures in large-eddy simulations of cloud-free and cloudy boundary layers. The sampling is evaluated against more traditional sampling of dry thermals or clouds. It enables the characterization of convective updrafts from the surface to the top of the boundary layer (or the top of cumulus clouds), describing in particular the transition from the sub-cloud to the cloud layer, and retrieves plume characteristics, entrainment and detrainment rates, variances and fluxes. This sampling is used to analyze the contribution of boundary-layer thermals to vertical fluxes and variances.     

Measurements Of Thermal Updraft Intensity Over Complex Terrain Using American White Pelicans And A Simple Boundary-Layer Forecast Model

An examination of boundary-layer meteorological and avian aerodynamic theories suggests that soaring birds can be used to measure the magnitude of vertical air motions within the boundary layer. These theories are applied to obtain mixed-layer normalized thermal updraft intensity over both flat and complex terrain from the climb rates of soaring American white pelicans and from diagnostic boundary-layer model-produced estimates of the boundary-layer depth z_i and the convective velocity scale w_*. Comparison of the flatland data with the profiles of normalized updraft velocity obtained from previous studies reveals that the pelican-derived measurements of thermal updraft intensity are in close agreement with those obtained using traditional research aircraft and large eddy simulation (LES) in the height range of 0.2 to 0.8 z_i. Given the success of this method, the profiles of thermal vertical velocity over the flatland and the nearby mountains are compared. This comparison shows that these profiles are statistically indistinguishable over this height range, indicating that the profile for thermal updraft intensity varies little over this sample of complex terrain. These observations support the findings of a recent LES study that explored the turbulent structure of the boundary layer using a range of terrain specifications. For terrain similar in scale to that encountered in this study, results of the LES suggest that the terrain caused less than an 11% variation in the standard deviation of vertical velocity.     

On the Nature of the Transition Between Roll and Cellular Organization in the Convective Boundary Layer

Both observational and numerical studies of the convective boundary layer (CBL) have demonstrated that when surface heat fluxes are small and mean wind shear is strong, convective updrafts tend to organize into horizontal rolls aligned within 10–20\(^\circ \) of the geostrophic wind direction. However, under large surface heat fluxes and weak to negligible shear, convection tends to organize into open cells, similar to turbulent Rayleigh-Bénard convection. Using a suite of 14 large-eddy simulations (LES) spanning a range of \(-z_i/L\) between zero (neutral) and 1041 (highly convective), where \(z_i\) is the CBL depth and L is the Obukhov length, the transition between roll- and cellular-type convection is investigated systematically for the first time using LES. Mean vertical profiles including velocity variances and turbulent transport efficiencies, as well the “roll factor,” which characterizes the rotational symmetry of the vertical velocity field, indicate the transition occurs gradually over a range of \(-z_i/L\); however, the most significant changes in vertical profiles and CBL organization occur from near-neutral conditions up to about \(-z_i/L \approx \) 15–20. Turbulent transport efficiencies and quadrant analysis are used to characterize the turbulent transport of momentum and heat with increasing \(-z_i/L\). It is found that turbulence transports heat efficiently from weakly to highly convective conditions; however, turbulent momentum transport becomes increasingly inefficient as \(-z_i/L\) increases.     

A Simple Model for the Vertical Transport of Reactive Species in the Convective Atmospheric Boundary Layer

We have developed a simple, steady-state, one-dimensional second-order closure model to obtain continuous profiles of turbulent fluxes and mean concentrations of non-conserved scalars in a convective boundary layer without shear. As a basic tool we first set up a model for conserved species with standard parameterizations. This leads to formulations for profiles of the turbulent diffusivity and the ratio of temperature-scalar covariance to the flux of the passive scalar. The model is then extended to solving, in terms of profiles of mean concentrations and fluxes, the NO _x –O₃ triad problem. The chemical reactions involve one first-order reaction, the destruction of NO₂ with decay time τ, and one second-order reaction, the destruction of NO and O₃ with the reaction constant k. Since the fluxes of the sum concentrations of NO _x = NO + NO₂ and O₃ + NO₂ turn out to be constant throughout the boundary layer, the problem reduces to solving two differential equations for the concentration and the flux of NO₂. The boundary conditions are the three surface fluxes and the fluxes at the top of the boundary layer, the last obtained from the entrainment velocity, and the concentration differences between the free troposphere and the top of the boundary layer. The equations are solved in a dimensionless form by using 1/(kτ) as the concentration unit, the depth h of the boundary layer as the length unit, the convective velocity scale w _* as the velocity unit, and the surface temperature flux divided by w _* as the temperature unit. Special care has been devoted to the inclusion of the scalar–scalar covariance between the concentrations of O₃ and NO. Sample calculations show that the fluxes of the reactive species deviate significantly from those of non-reactive species. Further, the diffusivities, defined by minus the flux divided by the concentration gradient may become negative for reactive species in contrast to those of non-reactive species, which in the present model are never negative.     

黄河源区鄂陵湖湖面和湖边草地对流边界层湍流结构特征的大涡模拟研究

为研究黄河源区边界层湍流特征及其对物质和能量输送的影响,本文首次采用大涡模拟的方法,对比分析了黄河源区两种不同下垫面上（鄂陵湖和湖边草地）对流边界层（CBL）中精细的湍流结构特征。使用资料为2012年夏季黄河源区鄂陵湖流域野外观测实验的GPS探空资料、涡动相关观测资料。分析表明,模拟的黄河源区草地和湖上CBL的平均结构与实测结果吻合较好,但草地和湖上CBL的湍流结构特征差异较明显。模拟结果显示,草地CBL内湍能收支、湍流特征量的时空分布和湍涡结构特征均与陆地上热力驱动CBL的研究结果一致;湖上CBL顶部存在明显的对流卷特征,且夹卷层的湍流强度比草地的强,而草地近地面湍强则更大。通过改变水平分辨率的模拟试验,发现两个不同下垫面上模拟结果对模式分辨率的敏感性不同,湖面CBL的模拟要选择较高的水平分辨率（50～100 m）,以提高近湖面和夹卷层对湍流动能和湍流通量模拟的精度,也充分模拟出各种尺度的波对湍流通量的累积贡献。考虑到计算时间等影响,模拟草地边界层精细的湍流结构时建议选择网格距为100～200 m。     

Large-eddy simulations of the effects of hilly terrain on the convective boundary layer

Large-eddy simulations of the convective boundary layer are compared over hilly versus flat surfaces. Moderate values for the height and horizontal spacing of the hills were selected. Thermally-direct hill-valley circulations are induced by the uneven terrain, accounting for a significant fraction of the resolved energy in the boundary-layer eddies. The probability of upward eddy motion reaches up to 70% over the hilltops and down to 15% over the valleys. Above-average values of both subgrid scale turbulent kinetic energy and upward eddy heat transport are found above the higher terrain. Horizontal spectra of vertical motion are strongly biased toward the horizontal scales of the terrain. Vertical profiles of atmospheric variables obtained by horizontal averaging, however, exhibit no significant differences between hilly and flat terrain simulations.     

A Parameterization of Dry Thermals and Shallow Cumuli for Mesoscale Numerical Weather Prediction

For numerical weather prediction models and models resolving deep convection, shallow convective ascents are subgrid processes that are not parameterized by classical local turbulent schemes. The mass flux formulation of convective mixing is now largely accepted as an efficient approach for parameterizing the contribution of larger plumes in convective dry and cloudy boundary layers. We propose a new formulation of the EDMF scheme (for Eddy Diffusivity\Mass Flux) based on a single updraft that improves the representation of dry thermals and shallow convective clouds and conserves a correct representation of stratocumulus in mesoscale models. The definition of entrainment and detrainment in the dry part of the updraft is original, and is specified as proportional to the ratio of buoyancy to vertical velocity. In the cloudy part of the updraft, the classical buoyancy sorting approach is chosen. The main closure of the scheme is based on the mass flux near the surface, which is proportional to the sub-cloud layer convective velocity scale w _*. The link with the prognostic grid-scale cloud content and cloud cover and the projection on the non- conservative variables is processed by the cloud scheme. The validation of this new formulation using large-eddy simulations focused on showing the robustness of the scheme to represent three different boundary layer regimes. For dry convective cases, this parameterization enables a correct representation of the countergradient zone where the mass flux part represents the top entrainment (IHOP case). It can also handle the diurnal cycle of boundary-layer cumulus clouds (EUROCS\ARM) and conserve a realistic evolution of stratocumulus (EUROCS\FIRE).     

Large-Eddy Simulations of Surface Heterogeneity Effects on the Convective Boundary Layer During the LITFASS-2003 Experiment

We investigate the impact of observed surface heterogeneities during the LITFASS-2003 experiment on the convective boundary layer (CBL). Large-eddy simulations (LES), driven by observed near-surface sensible and latent heat fluxes, were performed for the diurnal cycle and compare well with observations. As in former studies of idealized one- and two-dimensional heterogeneities, secondary circulations developed that are superimposed on the turbulent field and that partly take over the vertical transport of heat and moisture. The secondary circulation patterns vary between local and roll-like structures, depending on the background wind conditions. For higher background wind speeds, the flow feels an effective surface heat-flux pattern that derives from the original pattern by streamwise averaging. This effective pattern generates a roll-like secondary circulation with roll axes along the mean boundary-layer wind direction. Mainly the upstream surface conditions control the secondary circulation pattern, where the fetch increases with increasing background wind speed. Unlike the entrainment flux that appears to be slightly decreased compared to the homogeneously-heated CBL, the vertical flux of sensible heat appears not to be modified in the mixed layer, while the vertical flux of latent heat shows different responses to secondary circulations. The study illustrates that sufficient time averaging and ensemble averaging is required to separate the heterogeneity-induced signals from the raw LES turbulence data. This might be an important reason why experiments over heterogeneous terrain in the past did not give any clear evidence of heterogeneity-induced effects.     

Airborne measurements of turbulent fluxes during LITFASS-98: Comparison with ground measurements and remote sensing in a case study

Summary ?Simultaneous flight measurements with the research aircraft Do 128 and the helicopter-borne turbulence probe Helipod were performed on 18 June 1998 during the LITFASS-98 field experiment. The area-averaged turbulent vertical fluxes of momentum, sensible, and latent heat were determined on a 15 km × 15 km and a 10 km × 10 km flight pattern, respectively. The flights were carried out over heterogeneous terrain at different altitudes within a moderately convective boundary layer with Cumulus clouds. Co-spectra-analysis demonstrated that the small scale turbulent transport was completely sampled, while the comparatively small flight patterns were possibly of critical size regarding the large-scale turbulence. The phygoide of the airplane was identified as a significant peak in some co-spectra. The turbulent fluxes of momentum and sensible heat at 80 m above the ground showed systematic dependence on the location of the flight legs above the heterogeneous terrain. This was not observed for the latent heat flux, probably due to the vertical distribution of humidity in the boundary layer. Statistical error analysis of the fluxes F showed that the systematic statistical error ΔF was one order of magnitude smaller than the standard deviation σ _F . The difference between area-averaged fluxes derived from simultaneous Helipod and Do 128 measurements was much smaller than σ _F , indicating that the systematic statistical error was possibly over-estimated by the usual method. In the upper half of the boundary layer the airborne-measured sensible heat flux agreed well with windprofiler/RASS data. A linear fit was the best approximation for the height dependence of all three fluxes. The linear extrapolations of the latent and sensible heat fluxes to the ground were in good agreement with tower, scintillometer, and averaged ground-station measurements on various surface types. Systematic discrepancies between airborne and ground-based measurements were not found. Received June 18, 2001; revised December 21, 2001; accepted June 3, 2002     

Vertical ozone fluxes and related deposition parameters over agricultural and forested landscapes

The spatial variability and temporal behavior of the vertical flux of ozone have been investigated from turbulence measurements collected on aircraft flight legs in the daytime period during two consecutive summer experimental field programs. The data were obtained during horizontal flight legs conducted over agricultural crops and forested land in three different regions of the eastern United States.Results from individual experimental cases and statistics derived from all cases in each region are presented. Ozone flux generally exhibited a significant height dependency. The strongest negative (downward) fluxes in the lowest-level flight legs were primarily attributed to the uptake of ozone by the surface and vegetative cover. Fluxes were near-zero in the middle of the convective boundary layer (CBL) in the afternoon period. As ozone flux was proportional to concentration, slightly stronger fluxes were found in low-level urban plume segments where ozone concentrations were 10–20 ppb higher than in the surrounding area. The derived deposition velocity showed no such bias as a function of position across the urban plume. Ozone flux differences were not apparent between the more heavily forested sections and the primarily agricultural cropland areas in these regions. During the afternoon period when no clear temporal trend was evident, means from values obtained below 0.15Z _i (Z _i being the CBL height) were -0.421 and -0.431 ppb m^-2 s^-1 for ozone flux and 0.81 and 0.82 cm s^-1 for the derived mean deposition velocity in the southeastern Pennsylvania and central Ohio areas, respectively. These experimental results for ozone provide support to a dry deposition parameterization module which computes grid-area averaged deposition velocities for use in regional-scale models.On assignment from the National Oceanic and Atmospheric Administration, U.S. Department of Commerce.     

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.     

Vertical velocity field in the convective boundary layer as observed with an acoustic Doppler sodar

The morning development of the daytime convective boundary layer (CBL) during fine weather has been observed with an acoustic Doppler sodar of the C.R.P.E. In particular, the vertical profile of the vertical velocity third-order statistic W* ³ has been obtained. This quantity is a maximum near 0.3z _I where z _I, is the height of the CBL. The histogram of vertical velocity in the CBL shows a relationship between W ³ and the convective velocity W _* and is useful for convective plume determination.     

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.