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.     

Heterogeneity-Induced Heat-Flux Patterns in the Convective Boundary Layer: Can they be Detected from Observations and is There a Blending Height?—A Large-Eddy Simulation Study for the LITFASS-2003 Experiment

An understanding of how the convective boundary layer (CBL) is mixed under heterogeneous surface forcing is crucial for the interpretation of area-averaged turbulence measurements. To determine the height and degree to which a complex heterogeneous surface affects the CBL, large-eddy simulations (LES) for two days of the LITFASS-2003 experiment representing two different wind regimes were undertaken. Spatially-lagged correlation analysis revealed the turbulent heat fluxes to be dependent on the prescribed surface flux pattern throughout the entire CBL including the entrainment layer. These findings prompted the question of whether signals induced by surface heterogeneity can be measured by airborne systems. To examine this question, an ensemble of virtual flights was conducted using LES, according to Helipod flight measurements made during LITFASS-2003. The resulting ensemble-averaged heat fluxes indicated a clear dependence on the underlying surface up to the top of the CBL. However, a large scatter between the flux measurements in different ensemble runs was observed, which was the result of insufficient sampling of the largest turbulent eddies. The random and systematic errors based on the integral length scale did not indicate such a large scatter. For the given flight leg lengths, at least 10–15 statistically independent flight measurements were necessary to give a significant estimate of heterogeneity-induced signals in the CBL. The need for ensemble averaging suggests that the observed blending of heterogeneity-induced signals in the CBL can be partly attributed to insufficient averaging.     

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

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

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

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

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.     

Spectral Structure of Small-Scale Turbulent and Mesoscale Fluxes in the Atmospheric Boundary Layer over a Thermally Inhomogeneous Land Surface

Spectral analysis was performed on aircraft observations of a convective boundary layer (CBL) that developed over a thermally inhomogeneous, well-marked mesoscale land surface. The observations, part of the GAME-Siberia experiment, were recorded between April and June 2000 over the Lena River near Yakutsk City. A special integral parameter termed the ‘reduced depth of the CBL’ was used to scale the height of the mixed layer with variable depth. Analysis of wavelet cospectra and spectra facilitated the separation of fluxes and other variables into small-scale turbulent fluctuations (with scales less than the reduced depth of the CBL, approximately 2 km) and mesoscale fluctuations (up to 20 km). This separation approach allows for independent exploration of the scales. Analyses showed that vertical distributions obeyed different laws for small-scale fluxes and mesoscale fluxes (of sensible heat, water vapour, momentum and carbon dioxide) and for other variables (wind speed and air temperature fluctuations, coherence and degree of anisotropy). Vertical profiles of small-scale turbulent fluxes showed a strong decay that differed from generally accepted similarity models for the CBL. Vertical profiles of mesoscale fluxes and other variables clearly showed sharp inflections at the same relative (with respect to the reduced depth of the CBL) height of approximately 0.55 in the CBL. Conventional similarity models for sensible heat fluxes describe both small-scale turbulent and mesoscale flows. The present results suggest that mesoscale motions that reach up to the relative level of 0.55 could be initiated by thermal surface heterogeneity. Entrainment between the upper part of the CBL and the free atmosphere may cause mesoscale motions in that region of the CBL.     

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.     

The Influence of Thermally-Induced Mesoscale Circulations on Turbulence Statistics Over an Idealized Urban Area Under a Zero Background Wind

The influence of mesoscale circulations induced by urban-rural differential surface sensible heat flux and roughness on convective boundary-layer (CBL) flow statistics over an isolated urban area has been examined using large-eddy simulation (LES). Results are analyzed when the circulations influence the entire urban area under a zero background wind. For comparison, the CBL flow over an infinite urban area with identical urban surface characteristics under the same background meteorological conditions is generated as a control case (without circulations). The turbulent flow over the isolated urban area exhibits a mix of streaky structure and cellular pattern, while the cellular pattern dominates in the control case. The mixed-layer height varies significantly over the isolated urban area, and can be lower near the edge of the urban area than over the rural area. The vertical profiles of turbulence statistics over the isolated urban area vary horizontally and are dramatically different from the control case. The turbulent kinetic energy (TKE) sources include wind shear, convergence, and buoyancy productions, compared to only buoyancy production in the control case. The normalized vertical velocity variance is reduced compared to the control case except in the central urban area where it is little affected. The low-level flow convergence is mainly responsible for the enhanced horizontal velocity variance in the central urban area, while wind shear is responsible for the additional local maximum of the horizontal velocity variance near the middle of the CBL outside the central area. Parameterizations in the prognostic equation for TKE used in mesoscale models are evaluated against the LES results over the isolated urban area. We also discuss conditions under which the urban-induced circulations occur and when they may affect the entire urban area. Given that urban-induced circulations can influence the entire urban area within hours for an urban area of a realistic size, it is inappropriate to directly apply empirical relations of turbulence statistics derived under horizontally-homogenous flow conditions to an urban area.     

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.     

A Comparison of Higher-Order Vertical Velocity Moments in the Convective Boundary Layer from Lidar with In Situ Measurements and Large-Eddy Simulation

We have analyzed measurements of vertical velocity w statistics with the NOAA high resolution Doppler lidar (HRDL) from about 390 m above the surface to the top of the convective boundary layer (CBL) over a relatively flat and uniform agricultural surface during the Lidars-in-Flat-Terrain (LIFT) experiment in 1996. The temporal resolution of the zenith-pointing lidar was about 1 s, and the range-gate resolution about 30 m. Vertical cross-sections of w were used to calculate second- to fourth-moment statistics of w as a function of height throughout most of the CBL. We compare the results with large-eddy simulations (LES) of the CBL and with in situ aircraft measurements. A major cause of the observed case-to-case variability in the vertical profiles of the higher moments is differences in stability. For example, for the most convective cases, the skewness from both LES and observations changes more with height than for cases with more shear, with the observations changing more with stability than the LES. We also found a decrease in skewness, particularly in the upper part of the CBL, with an increase in LES grid resolution.     

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＊）.     

ON THE PDF MODELS FOR ATMOSPHERIC DIFFUSION IN BOUNDARY LAYER

In this paper, taking its turbulent exchange coefficient as a function of the Lagrangian timescale and standard variance of the turbulence in atmosphere, the atmospheric dispersion PDFmodels are obtained on the basis of atmospheric diffusion K-theory. In the model the statistics ofwind speed are directly used as its parameters instead of classic dispersion parameters. The bi-Gaussian PDF is derived in convective boundary layer (CBL), from the statistics of verticalvelocity in both of the downdraft and updraft regions that are investigated theoretically in the otherpart of this paper. Giving the driven parameters of the CBL (including the convective velocity scalew* and the mixing depth h_i) and the time-averaged wind speed at release level, the PDF model isable to simulate the distribution of concentration released at any levels in the CBL. The PDF'ssimulations are fairly consistent with the measurements in CONDORS experiment or the resultsbrought out by some numerical simulations.     

A Three-Dimensional Correlation/Spectral Model for Turbulent Velocities in a Convective Boundary Layer

A three-dimensional model for correlation functions and spectra in theatmospheric, convective boundary layer (CBL) is presented. The modelincludes vertical inhomogeneities introduced by eddy-blocking at the ground.By assuming the disturbance to the turbulent flow resulting from the groundblocking is irrotational, an equation is developed which allows one to writethe inhomogeneous, two-dimensional (2D) cross spectra for the blocked flowin terms of the 2D cross spectra for a homogeneous flow. VonKármán's energy spectrum then is used to determine thehomogeneous, 2D cross spectra. Although there are only two adjustableparameters in the model, the variance and a length scale, the model is shownto agree quite well with a diversity of previous results for the CBL.     

Temporal Oscillations in the Convective Boundary Layer Forced by Mesoscale Surface Heat-Flux Variations

A theoretical approach suggests that the surface heterogeneity on a scale of tens of kilometres can generate mesoscale motions that are not in a quasi-stationary state. The starting point of the theoretical approach is the equations of horizontal velocity and potential temperature that are low-pass filtered with a mesoscale cut-off wavelength. The transition of the generated mesoscale motions from a quasi-stationary state to a non-stationary state occurs when horizontal advection is strong enough to level out the potential temperature gradient on the surface heterogeneity scale. Large-eddy simulations (LES) suggest that the convective boundary layer (CBL) changes to a non-stationary state when forced by a surface heat-flux variation of amplitude of 100W m⁻² or higher and a wavelength of the order of 10 km. Spectral analysis of the LES reveals that when the mesoscale motions are in a quasi-stationary state, the energy provided by the surface heat-flux variation remains in organized mesoscale motions on the scale of the surface variation itself. However, in a non-stationary state, the energy cascades to smaller scales, with the cascade extending down into the turbulence scale when the wavelength of the surface heat-flux variation is on a scale smaller than 100 times the CBL height. The energy transfer from the generated mesoscale motions to the CBL turbulence results in the absence of a spectral gap between the two scales. The absence of an obvious spectral gap between the generated mesoscale motions and the turbulence raises questions about the applicability of mesoscale models for studies on the effect of high-amplitude surface heterogeneity on a scale of tens of kilometres. The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

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     

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.     

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.     

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

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

论边界层中的大气扩散PDF模式

基于大气扩散Ｋ理论，用作为风速脉动均方差和拉氏时间尺度函数的湍流交换系数，得到了直接利用风速脉动几率密度而不用扩散参数的大气扩散ＰＤＦ模式。分别研究了对流边界层上升气流区与下降区垂直速度的统计特征，求得双正态ＰＤＦ模式。在给定ＣＢＬ自身参数如对流特征速度ｗ＊，顶高ｈｉ和源高度上的平均风速时，该模式计算出的无量纲浓度分布与室内外测试结果一致。     

Application of a Large-Eddy Simulation Database to Optimisation of First-Order Closures for Neutral and Stably Stratified Boundary Layers

Large-eddy simulation (LES) is a well-established numerical technique, resolving the most energetic turbulent fluctuations in the planetary boundary layer. By averaging these fluctuations, high-quality profiles of mean quantities and turbulence statistics can be obtained in experiments with well-defined initial and boundary conditions. Hence, LES data can be beneficial for assessment and optimisation of turbulence closure schemes. A database of 80 LES runs (DATABASE64) for neutral and stably stratified planetary boundary layers (PBLs) is applied in this study to optimize first-order turbulence closure (FOC). Approximations for the mixing length scale and stability correction functions have been made to minimise a relative root-mean-square error over the entire database. New stability functions have correct asymptotes describing regimes of strong and weak mixing found in theoretical approaches, atmospheric observations and LES. The correct asymptotes exclude the need for a critical Richardson number in the FOC formulation. Further, we analysed the FOC quality as functions of the integral PBL stability and the vertical model resolution. We show that the FOC is never perfect because the turbulence in the upper half of the PBL is not generated by the local vertical gradients. Accordingly, the parameterised and LES-based fluxes decorrelate in the upper PBL. With this imperfection in mind, we show that there is no systematic quality deterioration of the FOC in the strongly stable PBL provided that the vertical model resolution is better than 10 levels within the PBL. In agreement with previous studies, we found that the quality improves slowly with the vertical resolution refinement, though it is generally wise not to overstretch the mesh in the lowest 500 m of the atmosphere where the observed, simulated and theoretically predicted stably stratified PBL is mostly located. The submission to a special issue of the “Boundary-Layer Meteorology” devoted to the NATO advanced research workshop “Atmospheric Boundary Layers: Modelling and Applications for Environmental Security”.