Non-local mixing of momentum in the convective boundary layer

Recently Frech and Mahrt proposed a closure scheme which includes alarge-scale stress term to represent the effects of non-local momentummixing in the convective boundary layer. Here large-eddy simulation (LES)datasets are used to evaluate the performance of this scheme across a rangeof stabilities between neutral and highly convective conditions, and as afunction of baroclinity. Generally the inclusion of the non-local term inthe closure model leads to results in better agreement with LES, althoughsome modifications to the model formulation are suggested.     

Impact of terrain heterogeneity on near-surface turbulence structure

This study investigates the impact of terrain heterogeneity on local turbulence measurements using 18 months of turbulence data taken on a 30 m tower at the SIRTA mixed land-use observatory under varying stability conditions and fetch configurations. These measurements show that turbulence variables such as the turbulent kinetic energy or momentum fluxes are strongly dependent on the upstream complexity of the terrain (presence of trees or buildings, open field). However, using a detection technique based on wavelet transforms which permits the isolation of the large-scale coherent structures from small-scale background fluctuations, the study shows that, for all stability conditions, whatever the upstream complexity of the terrain, the coherent structures display universal properties which are independent of the terrain nature: the frequency of occurrence, time duration of the coherent structures, the time separation between coherent structures and the relative contribution of the coherent structures to the total fluxes (momentum and heat) appear to be independent of the upstream roughness. This is an important result since coherent structures are known to transport a large portion of the total energy. This study extends to all stability conditions a numerical study by Fesquet et al. [Fesquet, C., Dupont, S., Drobinski, P., Barthlott, C., Dubos, T., 2008. Impact of terrain heterogeneities on coherent structures properties: experimental and numerical approaches. In: 18th Symposium on Boundary Layers and Turbulence. No. 11B.1. Stockholm, Sweden., Fesquet, C., Dupont, S., Drobinski, P., Dubos, T., Barthlott, C., in press. Impact of terrain heterogeneity on coherent structure properties: numerical approach. Bound.-Layer Meteorol.] conducted in neutral conditions which shows that a reason for such behavior is that the production of local active turbulence in an internal boundary layer associated with coherent structure originating from the outer layer and impinging onto the ground is not sensitive to the nature of the terrain.     

Long-term study of coherent structures in the atmospheric surface layer

A long-term study of coherent turbulence structures in the atmospheric surface layer has been carried out using 10 months of turbulence data taken on a 30-m tower under varying meteorological conditions. We use an objective detection technique based on wavelet transforms. The applied technique permits the isolation of the coherent structures from small-scale background fluctuations which is necessary for the development of dynamical models describing the evolution and properties of these phenomena. It was observed that coherent structures occupied 36% of the total time with mean turbulent flux contributions of 44% for momentum and 48% for heat. The calculation of a transport efficiency parameter indicates that coherent structures transport heat more efficiently than momentum. Furthermore, the transport efficiency increases with increasing contribution of the structures to the overall transport.     

城市建筑动力学效应对对流边界层影响的敏感性试验

本文将大涡模拟应用于城市对流边界层(CBL)湍流结构和流场特征的研究,在大涡模式中,拖曳系数取与建筑物高度及建筑物高度标准差有关的表达式以考虑次网格建筑物对风速和湍流动能(TKE)的面积平均影响.模拟结果表明,由于城市建筑物对气流的拖曳作用,使建筑物冠层及整个CBL内风速大幅度减小,城市冠层内部风速减小尤为明显,在夹卷层内,风速有一明显的跃变.在边界层中部对流运动已经发展成为较强的热泡,城市建筑物的动力学效应使热泡的水平尺度增大,CBL内平均上升气流速度和下沉气流速度减小,同时使CBL中上升气流所占比例比平坦地面增大.城市建筑物使CBL低层热通量、动量通量、速度方差和位温方差明显增大,但对近地层高度以上的湍流量影响不大.     

Turbulent mixing of reactive gases in the convective boundary layer

We study the interactions of chemistry and turbulent mixing of tracersin the convective boundary layer with a second-order closure model,including higher order chemistry terms. In order to limit the number of predictive equations we prescribe the profiles for ¯w¯, ¯w¯ ¯ and the lengthscale l. However, for model validation we treat temperature and humidity asinert tracers, and compare the results with profiles observed during theAir Mass Transformation Experiment, and with similarity expressions for thesurface layer. We find good agreement of the mean profiles, but the (co-)variances are slightly underpredicted. Furthermore, the model usesdiagnostic equations expressing third moments of concentration in terms ofsecond moments and their vertical derivatives. They are compared withlarge-eddy model results, showing good agreement and, therefore, thesimplifications are justified. The model is applied to the transport of two gases subject to one bimolecular reaction. The importance of concentration correlations on themean transformation rate is studied. For two gases diffusing in oppositedirections we find for moderate and fast chemistry a 50% and90% decreased transformation rate due to the negatively correlatedconcentrations. These values are similar to large-eddy results of Schumannand Sykes et al. For two bottom-up tracers we find that the covariance ofboth reactive species is either positive or negative, increasing or reducingthe effective transformation rate depending on the Damköhler number (the ratio of the turbulent and the chemistry timescale). A significantdirect influence of chemistry on the flux divergence is found in bothcases. According to the model the effective transport to mid-levels of theboundary layer is increased when two reactive tracers diffuse in oppositedirections, and decreased in the case of two bottom-up tracers.     

A numerical study of daily transitions in the convective boundary layer

We examine daily (morning–afternoon) transitions in the atmospheric boundary layer based on large-eddy simulations. Under consideration are the effects of the stratification at the top of the mixed layer and of the wind shear. The results describe the transitory behaviour of temperature and wind velocity, their second moments, the boundary-layer height Z _m (defined by the maximum of the potential temperature gradient) and its standard deviation σ _m , the mixed-layer height z _i (defined by the minimum of the potential temperature flux), entrainment velocity W _e, and the entrainment flux H _i . The entrainment flux and the entrainment velocity are found to lag slightly in time with respect to the surface temperature flux. The simulations imply that the atmospheric values of velocity variances, measured at various instants during the daytime, and normalized in terms of the actual convective scale w_*, are not expected to collapse to a single curve, but to produce a significant scatter of observational points. The measured values of the temperature variance, normalized in terms of the actual convective scale Θ_*, are expected to form a single curve in the mixed layer, and to exhibit a considerable scatter in the interfacial layer.     

Negative water vapour skewness and dry tongues in the convective boundary layer: observations and large-eddy simulation budget analysis

This study focuses on the intrusion of dry air into the convective boundary layer (CBL) originating from the top of the CBL. Aircraft in-situ measurements from the IHOP_2002 field campaign indicate a prevalence of negative skewness of the water vapour distribution within the growing daytime CBL over land. This negative skewness is interpreted according to large-eddy simulations (LES) as the result of descending dry downdrafts originating from above the mixed layer. LES are used to determine the statistical properties of these intrusions: their size and thermodynamical characteristics. A conditional sampling analysis demonstrates their significance in the retrieval of moisture variances and fluxes. The rapid CBL growth explains why greater negative skewness is observed during the growing phase: the large amounts of dry air that are quickly incorporated into the CBL prevent a full homogenisation by turbulent mixing. The boundary-layer warming in this phase also plays a role in the acquisition of negative buoyancy for these dry tongues, and thus possibly explains their kinematics in the lower CBL. Budget analysis helps to identify the processes responsible for the negative skewness. This budget study underlines the main role of turbulent transport, which distributes the skewness produced at the top or the bottom of the CBL into the interior of the CBL. The dry tongues contribute significantly to this turbulent transport.     

Three-Dimensional Scalar Microfront Systems in a Large-Eddy Simulation of Vegetation Canopy Flow

Large-eddy simulation is used to reproduce neutrallystratified airflow inside and immediately above a vegetation canopy. A passive scalaris released from the canopy and the evolution of scalar concentration above the canopyis studied. The most significant characteristic of the scalar concentration is the repeatedformation and dissipation of scalar microfronts, a phenomenon that has been observedin nature. These scalar microfronts consist of downstream-tilted regions of highscalar concentration gradients. Computer visualization tools and a conditional samplingand compositing technique are utilized to analyze these microfronts. Peaks in positivepressure perturbation exceeding an experimental threshold are found to be effectiveindicators of scalar microfronts. Convergence of the streamwise velocity componentand divergence of the cross-stream velocity component are observed in the immediatevicinity of scalar microfronts, which helps explain their relatively longlifetimes. Many of these three-dimensional features have been observedin previous field studies of canopy flow.     

Effects of mesoscale sea-surface temperature fronts on the marine atmospheric boundary layer

A numerical modelling study is presented focusing on the effects of mesoscale sea-surface temperature (SST) variability on surface fluxes and the marine atmospheric boundary-layer structure. A basic scenario is examined having two regions of SST anomaly with alternating warm/cold or cold/warm water regions. Conditions upstream from the anomaly region have SST values equal to the ambient atmosphere temperature, creating an upstream neutrally stratified boundary layer. Downstream from the anomaly region the SST is also set to the ambient atmosphere value. When the warm anomaly is upstream from the cold anomaly, the downstream boundary layer exhibits a more complex structure because of convective forcing and mixed layer deepening upstream from the cold anomaly. An internal boundary layer forms over the cold anomaly in this case, generating two distinct layers over the downstream region. When the cold anomaly is upstream from the warm anomaly, mixing over the warm anomaly quickly destroys the shallow cold layer, yielding a more uniform downstream boundary-layer vertical structure compared with the warm-to- cold case. Analysis of the momentum budget indicates that turbulent momentum flux divergence dominates the velocity field tendency, with pressure forcing accounting for only about 20% of the changes in momentum. Parameterization of surface fluxes and boundary-layer structure at these scales would be very difficult because of their dependence on subgrid-scale SST spatial order. Simulations of similar flow over smaller scale fronts (<5 km) suggest that small-scale SST variability might be parameterized in mesoscale models by relating the effective heat flux to the strength of the SST variance.     

Minimum friction velocity and heat transfer in the rough surface layer of a convective boundary layer

A simple model is deduced for the surface layer of a convective boundary layer for zero mean wind velocity over homogeneous rough ground. The model assumes large-scale convective circulation driven by surface heat flux with a flow pattern as it would be obtained by conditional ensemble averages. The surface layer is defined here such that in this layer horizontal motions dominate relative to vertical components. The model is derived from momentum and heat balances for the surface layer together with closures based on the Monin-Obukhov theory. The motion in the surface layer is driven by horizontal gradients of hydrostatic pressure. The balances account for turbulent fluxes at the surface and fluxes by convective motions to the mixed layer. The latter are the dominant ones. The model contains effectively two empirical coefficients which are determined such that the model's predictions agree with previous experimental results for the horizontal turbulent velocity fluctuations and the temperature fluctuations. The model quantitatively predicts the decrease of the minimum friction velocity and the increase of the temperature difference between the mixed layer and the ground with increasing values of the boundary layer/roughness height ratio. The heat transfer relationship can be expressed in terms of the common Nusselt and Rayleigh numbers, Nu and Ra, as Nu ~ Ra^{% MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaSGbaeaaca% aIXaaabaGaaGOmaaaaaaa!3779!\[{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}\]}. Previous results of the form Nu ~ Ra^{% MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaSGbaeaaca% aIXaaabaGaaG4maaaaaaa!377A!\[{1 \mathord{\left/ {\vphantom {1 3}} \right. \kern-\nulldelimiterspace} 3}\]} are shown to be restricted to Rayleigh-numbers less than a certain value which depends on the boundary layer/roughness height ratio.     

Analysis of turbulent exchange and coherent structures in the stable atmospheric boundary layer based on tower observations

Using data collected at the Spanish low troposphere research centre CIBA (Centro de Investigación de la Baja Atmósfera) and at the Cabauw Experimental Site for Atmospheric Research (CESAR) in the Netherlands, we analyzed the most significant features of different coherent structures occurring in the stable atmospheric boundary layer. In particular, we used both the Reynolds and wavelet methods to analyze a solitary wave, a gravity wave, a density current and a low-level jet. For each of these structures, we found that wavelet analysis had the capacity to distinguish the different scales involved in these events due to the different timing and heights of the thermal instabilities and downdrafts associated with the disturbances. In addition, the wavelet method highlights the different roles of turbulence and coherent structures in the transfer of heat, moisture and CO₂ in the nocturnal boundary layer.     

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

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

The Effects of Vegetation Density on Coherent Turbulent Structures within the Canopy Sublayer: A Large-Eddy Simulation Study

Large-eddy simulation has become an important tool for the study of the atmospheric boundary layer. However, since large-eddy simulation does not simulate small scales, which do interact to some degree with large scales, and does not explicitly resolve the viscous sublayer, it is reasonable to ask if these limitations affect significantly the ability of large-eddy simulation to simulate large-scale coherent structures. This issue is investigated here through the analysis of simulated coherent structures with the proper orthogonal decomposition technique. We compare large-eddy simulation of the atmospheric boundary layer with direct numerical simulation of channel flow. Despite the differences of the two flow types it is expected that the atmospheric boundary layer should exhibit similar structures as those in the channel flow, since these large-scale coherent structures arise from the same primary instability generated by the interaction of the mean flow with the wall surface in both flows. It is shown here that several important similarities are present in the two simulations: (i) coherent structures in the spanwise-vertical plane consist of a strong ejection between a pair of counter-rotating vortices; (ii) each vortex in the pair is inclined from the wall in the spanwise direction with a tilt angle of approximately 45°; (iii) the vortex pair curves up in the streamwise direction. Overall, this comparison adds further confidence in the ability of large-eddy simulation to produce large-scale structures even when wall models are used. Truncated reconstruction of instantaneous turbulent fields is carried out, testing the ability of the proper orthogonal decomposition technique to approximate the original turbulent field with only a few of the most important eigenmodes. It is observed that the proper orthogonal decomposition reconstructs the turbulent kinetic energy more efficiently than the vorticity.     

The role of surface heterogeneity in modelling the stable boundary layer

High-resolution numerical simulations are performed for three nights for each of two areas of the United Kingdom. Area-averaging techniques are then used to calculate effective stability functions for turbulence parametrizations in models using typical mesoscale and global spatial resolutions. Comparisons are made with parametrizations commonly used in numerical weather prediction models. The present results do not suggest that significant enhancement of the stability functions above 50 m is justified. Closer to the surface, more significant enhancement is observed in some regions. It is shown that the amount of enhancement is related to the variability of the orography.     

热带气旋边界层精细结构研究综述

热带气旋边界层中包含许多小尺度环流系统,其精细结构与台风灾害紧密联系。随着观测资料的丰富、理论研究的深入、计算能力的提升及国家防台减灾的迫切需要,近20 a热带气旋边界层精细结构的研究日益受到国内外学者的重视,并取得了显著进展,目前该方向成为当前热带气旋研究的热点领域之一。本文对近20 a热带气旋边界层精细结构的研究从观测和模拟两方面进行了回顾,主要集中在近年来发现的热带气旋边界层滚涡和龙卷尺度涡旋,并对热带气旋边界层精细结构的研究进行展望。     

Exploring the Effects of Microscale Structural Heterogeneity of Forest Canopies Using Large-Eddy Simulations

The Regional Atmospheric Modeling System (RAMS)-based Forest Large-Eddy Simulation (RAFLES), developed and evaluated here, is used to explore the effects of three-dimensional canopy heterogeneity, at the individual tree scale, on the statistical properties of turbulence most pertinent to mass and momentum transfer. In RAFLES, the canopy interacts with air by exerting a drag force, by restricting the open volume and apertures available for flow (i.e. finite porosity), and by acting as a heterogeneous source of heat and moisture. The first and second statistical moments of the velocity and flux profiles computed by RAFLES are compared with turbulent velocity and scalar flux measurements collected during spring and winter days. The observations were made at a meteorological tower situated within a southern hardwood canopy at the Duke Forest site, near Durham, North Carolina, U.S.A. Each of the days analyzed is characterized by distinct regimes of atmospheric stability and canopy foliage distribution conditions. RAFLES results agreed with the 30-min averaged flow statistics profiles measured at this single tower. Following this intercomparison, two case studies are numerically considered representing end-members of foliage and midday atmospheric stability conditions: one representing the winter season with strong winds above a sparse canopy and a slightly unstable boundary layer; the other representing the spring season with a dense canopy, calm conditions, and a strongly convective boundary layer. In each case, results from the control canopy, simulating the observed heterogeneous canopy structure at the Duke Forest hardwood stand, are compared with a test case that also includes heterogeneity commensurate in scale to tree-fall gaps. The effects of such tree-scale canopy heterogeneity on the flow are explored at three levels pertinent to biosphere-atmosphere exchange. The first level (zero-dimensional) considers the effects of such heterogeneity on the common representation of the canopy via length scales such as the zero-plane displacement, the aerodynamic roughness length, the surface-layer depth, and the eddy-penetration depth. The second level (one-dimensional) considers the normalized horizontally-averaged profiles of the first and second moments of the flow to assess how tree-scale heterogeneities disturb the entire planar-averaged profiles from their canonical (and well-studied planar-homogeneous) values inside the canopy and in the surface layer. The third level (three-dimensional) considers the effects of such tree-scale heterogeneities on the spatial variability of the ejection-sweep cycle and its propagation to momentum and mass fluxes. From these comparisons, it is shown that such microscale heterogeneity leads to increased spatial correlations between attributes of the ejection-sweep cycle and measures of canopy heterogeneity, resulting in correlated spatial heterogeneity in fluxes. This heterogeneity persisted up to four times the mean height of the canopy (h _c ) for some variables. Interestingly, this estimate is in agreement with the working definition of the thickness of the canopy roughness sublayer (2h _c –5h _c ).     

Influence of surface heterogeneity on scalar dissimilarity in the roughness sublayer

While it is generally known that surface heterogeneity weakens the application of Monin–Obukhov similarity (MOS), few studies have investigated how seasonal changes in the degree of surface heterogeneity at a particular site may influence the validity of the similarity application. Exploiting seasonal changes in forest function associated with senescence, we conduct a unique evaluation of the effects of surface heterogeneity on the validity of similarity theory at two sites through time. Using high frequency (10 Hz) velocity and scalar time series collected within the roughness sublayer over mixed hardwood deciduous and coniferous forests during both periods of peak leaf area and senescence of deciduous foliage, we examined conformity with proposed universal flux-variance predictions and agreement amongst normalized standard deviations of different scalars (temperature, water vapour and carbon dioxide concentrations). Normalized scalar standard deviations were elevated above MOS flux-variance predictions, with more pronounced deviations observed during and following senescence, particularly in the case of CO₂. Power-law scaling of normalized standard deviations as a function of stability was upheld and robust to seasonal changes in surface heterogeneity. However, dissimilarity of normalized standard deviations for the scalars increased during senescence, as heterogeneity in the source/sink field increased. Scalewise decomposition of scalar time series using wavelet analysis indicated that correlations between scalars were conservative through much of the inertial cascade but decayed for eddies < 10 m. Senescence lowered correlations between scalars over a wide range of eddy sizes. These results demonstrate how seasonal changes in surface physiology can cause a temporal production of heterogeneity in the source/sink field, thus weakening similarity applications in the roughness sublayer.     

Airborne measurements of turbulent trace gas fluxes and analysis of eddy structure in the convective boundary layer over complex terrain

Using the new high-frequency measurement equipment of the research aircraft DO 128, which is described in detail, turbulent vertical fluxes of ozone and nitric oxide have been calculated from data sampled during the ESCOMPTE program in the south of France. Based on airborne turbulence measurements, radiosonde data and surface energy balance measurements, the convective boundary layer (CBL) is examined under two different aspects. The analysis covers boundary-layer convection with respect to (i) the control of CBL depth by surface heating and synoptic scale influences, and (ii) the structure of convective plumes and their vertical transport of ozone and nitric oxides. The orographic structure of the terrain causes significant differences between planetary boundary layer (PBL) heights, which are found to exceed those of terrain height variations on average. A comparison of boundary-layer flux profiles as well as mean quantities over flat and complex terrain and also under different pollution situations and weather conditions shows relationships between vertical gradients and corresponding turbulent fluxes. Generally, NO_x transports are directed upward independent of the terrain, since primary emission sources are located near the ground. For ozone, negative fluxes are common in the lower CBL in accordance with the deposition of O₃ at the surface.The detailed structure of thermals, which largely carry out vertical transports in the boundary layer, are examined with a conditional sampling technique. Updrafts mostly contain warm, moist and NO_x loaded air, while the ozone transport by thermals alternates with the background ozone gradient. Evidence for handover processes of trace gases to the free atmosphere can be found in the case of existing gradients across the boundary-layer top. An analysis of the size of eddies suggests the possibility of some influence of the heterogeneous terrain in mountainous area on the length scales of eddies.