Turbulence Spectra and Cospectra Under the Influence of Large Eddies in the Energy Balance EXperiment (EBEX)

During the Energy Balance EXperiment, the patch-by-patch, flood irrigation in a flat cotton field created an underlying surface with heterogeneous soil moisture, leading to a dry (warm)-to-wet (cool) transition within the cotton field under northerly winds. Moreover, the existence of an extremely dry, large bare soil area upstream beyond the cotton field created an even larger step transition from the bare soil region to the cotton field. We investigated the turbulence spectra and cospectra in the atmospheric surface layer (ASL) that was disturbed by large eddies generated over regions upstream and also influenced by horizontal advection. In the morning, the ASL was unstable while in the afternoon a stable internal boundary layer was observed at the site. Therefore, the turbulence data at 2.7 and 8.7 m are interpreted and compared in terms of interactions between large eddies and locally generated turbulence under two atmospheric conditions: the unstable ASL beneath the convective boundary layer (CBL) (hereafter the unstable condition) and the stable ASL beneath the CBL (hereafter the stable condition). We identified the influences of multiple sizes of large eddies on ASL turbulence under both stratifications; these large eddies with multiple sizes were produced over the dry patches and dry, large bare soil areas upstream. As a consequence of the disturbance of large eddies, the broadening, erratic variability, and deviation of spectra and cospectra, relative to those described by Monin–Obukhov similarity theory, are evident in the low- to mid-frequencies. Transfer of momentum, heat, and water vapour by large eddies is distinctly observed from the turbulence cospectra and leads to significant run-to-run variations of residuals of the surface energy balance closure. Our results indicate that these large eddies have greater influences on turbulence at higher levels compared to lower levels, and in the unstable ASL compared to the stable ASL.     

干湿地表的湍流特征及其对深对流影响的大涡模拟

利用大涡模式模拟了对流边界层结构演变以及深对流触发过程.通过改变鲍恩比的敏感性试验研究不同大气初始状况下湿润和干旱下垫面湍流特征及其对深对流触发过程的影响.结果表明:干旱下垫面的混合层干而暖,厚度较大;湿润下垫面相反.由于地表感热通量对热力湍流形成的作用更大,干旱下垫面上湍流混合和夹卷作用更强,使得水汽和相当位温在边界...     

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).     

非均匀网格的过渡湍流理论及其在大气边界层数值模拟中的应用

本文将Stull提出的均匀网格下的过渡湍流理论推广到非均匀网格情形，推广的非均匀网格的过渡湍流理论满足Stull提出的对过渡矩阵系数的要求并具有清晰的物理意义。然后，将非均匀网格的过渡湍流理论应用于一维大气边界层数值模式中，对Wangara资料进行了模拟，并与均匀网格情形进行了对比。计算表明，非均匀网格的过渡湍流模式能很好地模拟Wangara大气边界层平均量与湍流量的变化；本文非均匀网格的过渡湍流理论的推广是可行的，它可能会在大气边界层数值模拟及其他方面（如：中尺度模式）得到应用。     

Variability of precipitation intensity: sensitivity to treatment of moist convection in an RCM and a GCM

The present study investigates the sensitivity of the frequency distribution of precipitation rates to the closure employed in the penetrative mass flux cumulus parameterization of Zhang and McFarlane in the Canadian regional climate model (CRCM) and in the Canadian Centre for Climate Modelling and Analysis third generation global atmospheric general circulation model (AGCM3). The effects of an alternative prognostic closure for mass flux cumulus parameterization in place of the original diagnostic closure are investigated. A set of experiments is performed in which changes in the frequency distribution of precipitation rates and cloud base mass-flux are examined as a function of the parameters that define each closure scheme. The relationship between the frequency distribution of precipitation and cloud base mass flux is examined and a self-consistent relation is found when the depth of convection is taken into account. Experiments performed with the prognostic closure favor relatively strong cloud base mass-flux and deep penetrative convection with relatively more intense convective precipitation. The mean of the frequency distribution of convective precipitation is larger and the heavier events become more intense. Also, experiments performed with the prognostic closure favor less frequent convective activity. However these changes in the distribution of convective component of precipitation are generally offset by opposite changes in the distribution of the resolved large-scale component of precipitation, resulting in relatively smaller changes in total precipitation. The altered partition of precipitation between convective and large-scale components is found to alter the energy balance and the thermodynamic equilibrium structure of the troposphere. The robustness found in the CRCM results regarding the sensitivity of the frequency distribution of precipitation to changes in the closure of the deep convection parameterization is investigated by performing a similar analysis of AGCM3 simulations. A remarkable similarity of AGCM3 and CRCM results is found suggesting that the closure sensitivity identified in this study is robust.     

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.     

Laboratory and Numerical Studies of the Convective Boundary Layer Capped by a Strong Inversion

The convective boundary layer (CBL) with a wide range of stability is simulated experimentally using a thermally stratified wind tunnel, and numerically by direct numerical simulation (DNS). The turbulence structures and flow characteristics of various CBL flows, capped by a strong temperature inversion and affected by surface shear, are investigated. The various vertical profiles of turbulence statistics similar to those from the observed CBL in the field are successfully simulated in both the wind-tunnel experiment and in DNS. The comparison of the wind-tunnel data and DNS results with those of atmospheric observations and water-tank studies shows the crucial dependence of the turbulence statistics in the upper part of the layer on the strength of the inversion layer, as well as the modification of the CBL turbulence regime by the surface shear.     

复杂下垫面模拟域大气边界层非局地闭合模拟研究

在非局地闭合大气边界层模拟研究基础上,选取黑河实验研究资料（HEIFE）,利用二维模式对绿洲和沙漠交界地区的内边界层结构进行了模拟研究,模式考虑了两种下垫面条件下动力、热力特性的不同,并加入了湿度及水汽通量的预报。模拟结果表明:二维非局地闭合模式能得到合理的TIBL结构的基本特征,包括平均场、湍流场等,并能模拟给出一些由下垫面特性形成的大气现象,如沙漠地区逆湿现象的生消过程、绿洲的冷岛效应、绿洲与沙漠因热力差异而形成的局部环流形势等细致结构,符合沙漠站白天观测事实的典型特征。     

Numerical Simulation of Roll Vortices in the Convective Boundary Layer

Roll vortices,which often appear when cold air outbreaks over warm ocean surfaces,are an important system for energy and substance exchange between the land surface and atmosphere.Numerical simulations were carried out in the study to simulate roll vortices in the convective boundary layer(CBL).The results indicate,that with proper atmospheric conditions,such as thermal instability in the CBL,stable stratification in the overlying layer and suitable wind shear,and a temperature jump between the two layers in a two-layer atmosphere,convective bands appear after adding initial pulses in the atmosphere.The simulated flow and temperature fields presented convective bands in the horizontal and roll vortices in the crosswind section. The structure of the roll vortices were similar to those observed in the cloud streets,as well as those from analytical solutions.Simulations also showed the influence of depth and instability strength of the CBL, as well as the stratification above the top of the CBL,on the orientation spacing and strength of the roll vortices.The fluxes caused by the convective rolls were also investigated,and should perhaps be taken into account when explaining the surface energy closure gap in the CBL.     

一个对流边界层中的随机扩散模式

通过对对流边界层(CBL)湍流结构的分析,首次提出用两种不同尺度的湍流模拟CBL中的铅直扩散.在此基础上发展了一个随机扩散模式,并用它模拟了典型对流条件下两种高架连续点源的扩散.与Willis水槽模拟和Lamb等人数值模拟以及CONDORS计划外场试验的结果的比较表明,本模式能成功地模拟CBL中的横向积分浓度.与其它数值模式相比,还具有输入参数少、计算量极小和更加简单实用的优点.     

Evaluation of the Weather Research and Forecast/Urban Model Over Greater Paris

Meteorological modelling in the planetary boundary layer (PBL) over Greater Paris is performed using the Weather Research and Forecast (WRF) numerical model. The simulated meteorological fields are evaluated by comparison with mean diurnal observational data or mean vertical profiles of temperature, wind speed, humidity and boundary-layer height from 6 to 27 May 2005. Different PBL schemes, which parametrize the atmospheric turbulence in the PBL using different turbulence closure schemes, may be used in the WRF model. The sensitivity of the results to four PBL schemes (two non-local closure schemes and two local closure schemes) is estimated. Uncertainties in the PBL schemes are compared to the influence of the urban canopy model (UCM) and the updated Coordination of Information on the Environment (CORINE) land-use data. Using the UCM and the CORINE land-use data produces more realistic modelled meteorological fields. The wind speed, which is overestimated in the simulations without the UCM, is improved below 1,000 m height. Furthermore, the modelled PBL heights during nighttime are strongly modified, with an increase that may be as high as 200 %. At night, the impact of changing the PBL scheme is lower than the impact of using the UCM and the CORINE land-use data.     

The Adaptation of the Atmospheric Boundary Layer to a Change in Surface Roughness

The adaptation of the atmospheric boundary layer to a change in the underlying surface roughness is an interesting problem and hence much research, theoretical, experimental, and numerical, has been undertaken. Within the atmospheric boundary layer an accurate numerical model for the turbulent properties of the atmospheric boundary layer needs to be implemented if physically realistic results are to be obtained. Here, the adaptation of the atmospheric boundary layer to a change in surface roughness is investigated using a first-order turbulence closure model, a one-and-a-half-order turbulence closure model and a second-order turbulence closure model. Perturbations to the geostrophic wind and the pressure gradients are included and it is shown that the second-order turbulence closure model, namely the standard k - model, is inferior to a lower-order closure model if a modification to limit the turbulent eddy size within the atmospheric boundary layer is not included within the model.     

A comparison of local and non-local turbulence closure methods for the case of a cold air outbreak

Numerical experiments have shown that large-eddy-simulation models (LES) are able to reproduce the common features of convective boundary layers (CBL) quite well. Models which cannot resolve the convective motions due to their grid structure (1D-models or models with coarse horizontal and/or vertical resolution) have to take into account the effects of large eddies within their subgrid diffusion terms. Turbulent fluxes are frequently parameterized through first-order-closure methods (K-theory). Recently, non-local closure schemes have also been developed. In this paper we compare 1D-and 2D-models using different local and non-local first-order closure methods. The analysis is carried out for the case of an idealized cold air outbreak (CAO). One of the non-local closures is based on the so-called transilient turbulence theory. The reference states are given by a bulk-model and a 2D-model which resolves the large eddies explicitly. A comparison of the results is presented for characteristic quantities such as evolution of boundary-layer height and surface heat flux as well as mean wind and temperature profiles. It is found that simple local first-order closure does not give good agreement with the reference models. The results of the transilient turbulence model shows that a non-local closure is able to parameterize the effects of the large eddies. Comparable results are produced by a local closure where eddy diffusivities are parameterized by dimensionless gradient-functions.     

Application of the E-ε closure model to simulations of mesoscale topographic effects

A mesoscale Planetary Boundary Layer (PBL) model with a simple turbulence closure scheme based on the turbulence kinetic energy (TKE) equation and the dissipation () equation is used to simulate atmospheric flow over mesoscale topography. Comparative studies with different parameterizations suggest that with a proper closure assumption for turbulence dissipation, the E-model can simulate the circulation induced by the mesoscale topography with results similar to those obtained using the E- model. On the other hand, the first-order closure using O'Brien's cubic interpolation for eddy diffusivities (K) generally produces much larger K profiles in the stable and the unstable regions, which is believed to be due to the overprediction of the height of the PBL. All models with the TKE equation yield quite similar ensemble mean fields, which are found to be little sensitive to the closure assumption for turbulence dissipation, though their predicted magnitudes of TKE and K may differ appreciably. A discussion on the diurnal evolution of the mesoscale topography-induced circulation and the spatial variations of the turbulence fluxes in the surface layer is also given based on the E- model results.     

Toward increasing the accuracy and realism of future optical turbulence calculations

Summary Due to the increased use of laser and ground-to-satellite communications the need for reliable optical turbulence information is growing. Optical turbulence information is important because it describes an atmospheric effect that can degrade the performance of electromagnetic systems and sensors, e.g., free-space optical communications and infrared imaging. However, analysis of selected past research indicates that there are some areas (i.e., data and models) in which optical turbulence information is lacking. For example, line-of-sight optical turbulence data coupled with atmospheric models in hilly terrain, coastal areas, and within cities are few in number or non-existent. In addition, the bulk of existing atmospheric computer models being used to provide estimates of optical turbulence are basically one-dimensional in nature and assume uniform turbulence conditions over large areas. As a result, current optical turbulence theory and models may be deficient and in error for inhomogeneous (nonuniform) turbulence conditions, such as those that occur in urban environments or environments with changing topography and energy budgets. While it is anticipated that theoretical advances in environmental physics (and like disciplines) will be a catalyst for much new work this area, in the interim, we suggest that some very practical computational research can be performed to extend existing low-atmospheric turbulence and micrometeorological calculations beyond current limitations.     

Large-Eddy Simulation Comparison of Neutral Flow Over a Canopy: Sensitivities to Physical and Numerical Conditions,and Similarity to Other Representations

The representation of a neutral atmospheric flow over roughness elements simulating a vegetation canopy is compared between two large-eddy simulation models, wind-tunnel data and recently updated empirical flux-gradient relationships. Special attention is devoted to the dynamics in the roughness sublayer above the canopy layer, where turbulence is most intense. By demonstrating that the flow properties are consistent across these different approaches, confidence in the individual independent representations is bolstered. Systematic sensitivity analyses with the Dutch Atmospheric Large-Eddy Simulation model show that the transition in the one-sided plant-area density from the canopy layer to unobstructed air potentially alters the flow in the canopy and roughness sublayer. Anomalously induced fluctuations can be fully suppressed by spreading the transition over four steps. Finer vertical resolutions only serve to reduce the magnitude of these fluctuations, but do not prevent them. To capture the general dynamics of the flow, a resolution of 10 % of the canopy height is found to suffice, while a finer resolution still improves the representation of the turbulent kinetic energy. Finally, quadrant analyses indicate that momentum transport is dominated by the mean velocity components within each quadrant. Consequently, a mass-flux approach can be applied to represent the momentum flux.     

An Optimal Inverse Method Using Doppler Lidar Measurements to Estimate the Surface Sensible Heat Flux

Inverse methods are widely used in various fields of atmospheric science. However, such methods are not commonly used within the boundary-layer community, where robust observations of surface fluxes are a particular concern. We present a new technique for deriving surface sensible heat fluxes from boundary-layer turbulence observations using an inverse method. Doppler lidar observations of vertical velocity variance are combined with two well-known mixed-layer scaling forward models for a convective boundary layer (CBL). The inverse method is validated using large-eddy simulations of a CBL with increasing wind speed. The majority of the estimated heat fluxes agree within error with the proscribed heat flux, across all wind speeds tested. The method is then applied to Doppler lidar data from the Chilbolton Observatory, UK. Heat fluxes are compared with those from a mast-mounted sonic anemometer. Errors in estimated heat fluxes are on average 18 %, an improvement on previous techniques. However, a significant negative bias is observed (on average $-63\,\%$ ) that is more pronounced in the morning. Results are improved for the fully-developed CBL later in the day, which suggests that the bias is largely related to the choice of forward model, which is kept deliberately simple for this study. Overall, the inverse method provided reasonable flux estimates for the simple case of a CBL. Results shown here demonstrate that this method has promise in utilizing ground-based remote sensing to derive surface fluxes. Extension of the method is relatively straight-forward, and could include more complex forward models, or other measurements.     

Statistical Distribution of Atmospheric Dust Devils on Earth and Mars

Boundary-Layer Meteorology - Numerical weather prediction is often used to supply the mean wind and turbulence fields for atmospheric transport and dispersion plume models as they provide dense...     

Analytical Model of Roll Vortices in the Convective Boundary Layer

Analytical solutions of convective waves in the convective boundary layer (CBL) were obtained with two-layer linearized atmospheric equations including Rayleigh friction, which represents the turbulent viscosity in the lower CBL. The analytical model shows that the interaction between the convection in the lower layer and gravity waves in the upper layer is one of the causes for the formation of convective bands. The flow and temperature fields obtained by the analytical model present the main characteristics of convective bands found in field observations. We have also investigated the influences of atmospheric conditions on the characteristics of the bands. Results accord with previous knowledge about these phenomena.