Laboratory Experiments on Convective Entrainment Using a Saline Water Tank

Entrainment fluxes in a shear-free convective boundary layer have been measured with a saline water tank set-up. The experiments were targeted towards measuring the entrainment behaviour for medium to high Richardson numbers and use a two-layer design, i.e. two stacked non-stratified (neutral) layers with different densities. With laser induced fluorescence (LIF), the entrainment flux of a fluorescent dye is measured for bulk Richardson numbers in the range 30–260. It is proposed that a carefully chosen combination of top-down and bottom-up processes improves the accuracy of LIF-based entrainment observations. The observed entrainment fluxes are about an order of magnitude lower than reported for thermal water tanks: the derived buoyancy entrainment ratio, $A$ , is found to be $A \approx 0.02$ , which is to be compared with $A\approx 0.25$ for a thermal convection tank (Deardorff et al., J Fluid Mech 100:41–64, 1980). An extensive discussion is devoted to the influence of the Reynolds and Prandtl numbers in laboratory experiments on entrainment.     

Some characteristics of cumulus convection over the Tibetan Plateau

The diagnostic model of the cumulus convection proposed by Yanai et al. (1973) was applied to the atmosphere over the Tibetan Plateau, and used to estimate the vertical mass flux, entrainment and detrainment, excess temperature and moisture, liquid water content, and condensation and precipitation rates of highland cloud clusters. The re-sults illustrated that in clouds over the Tibetan Plateau, the water vapor condensation rate, liquid water content, and efficiency of the rain generation process are less than those in the tropics (represented by the Marshall Islands region). Therefore, the condensational latent heat released over the Tibetan Plateau, overall, is much smaller than that in the tropics. The water vapor and liquid water detrainment from shallow nonprecipitating cumulus clouds, and their entrainment into deep cumulus clouds, serve as a growing mechanism for the deep precipitating cumulus towers over the Tibetan Plateau. It should be noted that there is a stronger detrainment of liquid water from cumulus clouds and a stronger re-evaporation rate in environment. The process of the condensation-detrainment-re-evaporation-entrainment is repeatedly in progress. It would play an important role in maintaining of cumulus convection on the condition that the supply of moisture is not plentiful over the Tibetan Plateau.The analyses also showed that the cloud mass flux Mc over the Tibetan Plateau is less, and the large-scale av-erage upward motion is much less than those over the Marshall Islands. Stronger compensating downward motion in the cloud environment over the Tibetan Plateau, responsible for the area’s strong environmental heating rate was re-vealed, and would link to the stability of the South Asian High in summer.     

A note on hydraulic theory of internal bores

Mass and momentum conservation across an internal bore, together with an assumption that energy is dissipated in both fluid layers, yields a range of possible bore speeds. The upper speed limit is that given by Wood and Simpson [Wood, I.R., Simpson, J.E., 1984. Jumps in layered miscible fluids. J. Fluid Mech., 140: 215–231.] who assume no energy dissipation in the contracting layer, while the lower limit is that of Klemp et al. [Klemp, J.B., Rotunno, R., Skamarock, W.C., 1997. On the propagation of internal bores. J. Fluid Mech., 331: 81–106] who assume no energy dissipation in the expanding layer. The two bore speeds agree to within a few percent, except when the expanding layer is shallow upstream and the internal bore propagates as a gravity current.     

A mass-flux cumulus parameterization scheme for large-scale models: description and test with observations

A simple mass-flux cumulus parameterization scheme suitable for large-scale atmospheric models is presented. The scheme is based on a bulk-cloud approach and has the following properties: (1) Deep convection is launched at the level of maximum moist static energy above the top of the boundary layer. It is triggered if there is positive convective available potential energy (CAPE) and relative humidity of the air at the lifting level of convection cloud is greater than 75%; (2) Convective updrafts for mass, dry static energy, moisture, cloud liquid water and momentum are parameterized by a one-dimensional entrainment/detrainment bulk-cloud model. The lateral entrainment of the environmental air into the unstable ascending parcel before it rises to the lifting condensation level is considered. The entrainment/detrainment amount for the updraft cloud parcel is separately determined according to the increase/decrease of updraft parcel mass with altitude, and the mass change for the adiabatic ascent cloud parcel with altitude is derived from a total energy conservation equation of the whole adiabatic system in which involves the updraft cloud parcel and the environment; (3) The convective downdraft is assumed saturated and originated from the level of minimum environmental saturated equivalent potential temperature within the updraft cloud; (4) The mass flux at the base of convective cloud is determined by a closure scheme suggested by Zhang (J Geophys Res 107(D14), doi:10.1029/2001JD001005, 2002) in which the increase/decrease of CAPE due to changes of the thermodynamic states in the free troposphere resulting from convection approximately balances the decrease/increase resulting from large-scale processes. Evaluation of the proposed convection scheme is performed by using a single column model (SCM) forced by the Atmospheric Radiation Measurement Program’s (ARM) summer 1995 and 1997 Intensive Observing Period (IOP) observations, and field observations from the Global Atmospheric Research Program’s Atlantic Tropical Experiment (GATE) and the Tropical Ocean and Global Atmosphere Coupled Ocean–Atmosphere Response Experiment (TOGA COARE). The SCM can generally capture the convective events and produce a realistic timing of most events of intense precipitation although there are some biases in the strength of simulated precipitation.     

Note on the generalized thermal theory for gravity currents in the deceleration phase

The generalized thermal theory for gravitational convection, produced from instantaneous buoyancy sources on sloping boundaries, developed in Dai and Garcia (2010) is examined in this note. An assumption implicitly made therein, that detrained fluid carries no momentum, was inappropriate and the solution was not physical in special cases. The generalized thermal theory is now improved by considering the momentum carried away by detrained mixed fluid. An asymptotic velocity–distance relation for gravity currents further downslope in the deceleration phase is provided and agreement with reported experimental data is found.     

Improvements in Climate Simulation with Modifications to the Tiedtke Convective Parameterization in the Grid-Point Atmospheric Model of IAP LASG （GAMIL）

The grid-point atmospheric model of IAP LASG （GAMIL） was developed in and has been evaluated since early 2004. Although the model shows its ability in simulating the global climate, it suffers from some problems in simulating precipitation in the tropics. These biases seem to result mainly from the treatment of the subgrid scale convection, which is parameterized with Tiedtke＇s massflux scheme （or the Zhang-McFarlane scheme, as an option） in the model. In order to reduce the systematic biases, several modifications were made to the Tiedtke scheme used in GAMIL, including （1） an increase in lateral convective entrainment/detrainment rate for shallow convection, （2） inclusion of a relative humidity threshold for the triggering of deep convection, and （3） a reduced efficiency for the conversion of cloud water to rainwater in the convection scheme.
Two experiments, one with the original Tiedtke scheme used in GAMIL and the other with the modified scheme, were conducted to evaluate the performance of the modified scheme in this study. The results show that both the climatological mean state, such as precipitation, temperature and specific humidity, and interannual variability in the model simulation are improved with the use of this modified scheme. Results from several additional experiments show that the improvements in the model performance in different regions mainly result from either the introduction of the relative humidity threshold for triggering of the deep convection or the suppressed shallow convection due to enhanced lateral convective entrainment/detrainment rates.     

Particle dispersion and mixing induced by breaking internal gravity waves

The purpose of this paper is to analyze diapycnal mixing induced by the breaking of an internal gravity wave — the primary wave — either standing or propagating. To achieve this aim we apply two different methods. The first method consists of a direct estimate of vertical eddy diffusion from particle dispersion while the second method relies upon potential energy budgets [Winters, K.B., Lombard, P.N., Riley, J.J., D’Asaro, E.A., 1995. J. Fluid Mech. 289, 115–128; Winters, K.B., D’Asaro, E.A., 1996. J. Fluid Mech. 317, 179–193]. The primary wave we consider is of small amplitude and is statically stable, a case for which the breaking process involves two-dimensional instabilities. The dynamics of the waves have been previously analyzed by means of two-dimensional direct numerical simulations [Bouruet-Aubertot, P., Sommeria, J., Staquet, C., 1995. J. Fluid Mech. 285, 265–301; Bouruet-Aubertot, P., Sommeria, J., Staquet, C., 1996. Dyn. Atmos. Oceans 29, 41–63; Koudella, C., Staquet, C., 1998. In: Davis, P. (Ed.), Proceedings of the IMA Conference on Mixing and Dispersion on Stably-stratified Flows, Dundee, September 1996. IMA Publication]. High resolution three-dimensional calculations of the same wave are also reported here [Koudella, C., 1999].A local estimate of mixing is first inferred from the time evolution of sets of particles released in the flow during the breaking regime. We show that, after an early evolution dominated by shear effects, a diffusion law is reached and the dispersion coefficient is fairly independent of the initial seeding location of the particles in the flow.The eddy diffusion coefficient, K, is then estimated from the diapycnal diffusive flux. A good agreement with the value inferred from particle dispersion is obtained. This finding is of particular interest regarding the interpretation of in situ estimates of K inferred either from tracer dispersion or from microstructure measurements. Computation of the Cox number, equal to the ratio of eddy diffusivity to molecular diffusivity, shows that the Cox number varies within the interval [9, 262], which corresponds to the range of vertical eddy diffusivity measured in the interior of the ocean. The Cox number is found to depend on the turbulent Froude number squared.We show eventually that mixing results in a weak distortion of the initial density profile and we relate this result to observations made at small scale in the ocean.Comparisons between the analysis of the two-dimensional and high resolution (256³) three-dimensional direct numerical simulations of the primary wave were also conducted. We show that the energetics and the amount of mixing are very close when the primary wave is of small amplitude. This results from the fact that, for a statically stable wave, the dynamics of the initially two-dimensional primary wave remains mostly two-dimensional even after the onset of wavebreaking.     

修正的质量通量积云对流方案及其模拟试验研究 Ⅱ:三种积云方案的积云对流活动及MFS方案相关参数的敏感性试验

作者在“修正的质量通量积云对流方案及其模拟试验研究I方案介绍及对1991年洪涝过程的模拟”论文中提出的质量通量方案成功地植入区域气候模式RegCM2的基础上，对比分析质量通量方案MFS，Kuo方案和Grell方案对积云对流活动的模拟，结果表明质量通量方案较好地模拟了积云对流活动过程。针对质量通量方案中的一些参数，如云水向雨水的转换率、混合卷入率、混合卷出率和下沉气流的强度等进行了一系列的敏感性试验，试验结果表明积云对流活动对这些参数比较敏感，因此提高这些参数的准确性是改进积云对流参数化方案的重要内容之一。     

MODIFIED MASS FLUX CUMULUS CONVECTIVE PARAMETERIZATION SCHEME AND ITS SIMULATION EXPERIMENT—PARTⅡ:CUMULUS CONVECTION ACTIVITIES SIMULATED FROM THREE SCHEMES AND SENSITIVITY EXPERIMENTS BY USE OF MASS FLUX SCHEME*

By comparison of simulated cumulus convection processes in RegCM2,using the Kuo scheme,the Grell scheme and the mass flux scheme (MFS),it is found that the MFS can simulate the cumulus heating and moistening very well.A series of sensitivity tests show that the parameters for specifying the conversion coefficient from cloud droplets to raindrops,the turbulent entrainment and detrainment rates in updrafts anddowndrafts,and the intensity of the downdrafts have different degrees of influence upon the cumulus convection.Therefore.it is quite important for cumulus parameterization scheme to define these parameters as accurately as possible.     

MODIFIED MASS FLUX CUMULUS CONVECTIVE PARAMETERIZATION SCHEME AND ITS SIMULATION EXPERIMENT—PART II:CUMULUS CONVECTION ACTIVITIES SIMULATED FROM THREE SCHEMES AND SENSITIVITY EXPERIMENTS BY USE OF MASS FLUX SCHEME

By comparison of simulated cumulus convection processes in RegCM2,using the Kuo scheme,the Grell scheme and the mass flux scheme (MFS),it is found that the MFS can simulate thecumulus heating and moistening very well.A series of sensitivity tests show that the parametersfor specifying the conversion coefficient from cloud droplets to raindrops,the turbulententrainment and detrainment rates in updrafts anddowndrafts,and the intensity of thedowndraftshave different degrees of influence upon the cumulus convection.Therefore.it is quite importantfor cumulus parameterization scheme to define these parameters as accurately as possible.     

The Volumetric Particle Approach for Concentration Fluctuations and Chemical Reactions in Lagrangian Particle and Particle-grid Models

A new approach is proposed to predict concentration fluctuations in the framework of one-particle Lagrangian stochastic models. The approach is innovative since it allows the computation of concentration fluctuations in dispersing plumes using a Lagrangian one-particle model with micromixing but with no need for the simulating of background particles. The extension of the model for the treatment of chemically reactive plumes is also accomplished and allows the computation of plume-related chemical reactions in a Lagrangian one-particle framework separately from the background chemical reactions, accounting for the effect of concentration fluctuations on chemical reactions in a general, albeit approximate, manner. These characteristics should make the proposed approach an ideal tool for plume-in-grid calculations in chemistry transport models. The results are compared to the wind-tunnel experiments of Fackrell and Robins (J Fluid Mech, 117:1–26, 1982) for plume dispersion in a neutral boundary layer and to the measurements of Legg et al. (Boundary-Layer Meteorol, 35:277–302, 1986) for line source dispersion in and above a model canopy. Preliminary reacting plume simulations are also shown comparing the model with the experimental results of Brown and Bilger (J Fluid Mech, 312:373–407, 1996; Atmos Environ, 32:611–628, 1998) to demonstrate the feasibility of computing chemical reactions in the proposed framework.     

The effect of ageostrophy on the stability of thin oceanic vortices

This paper examines the stability of vortices in a two-layer ocean on the f-plane. The mean depth of the upper layer is assumed to be much smaller than the depth of the lower layer. Using the primitive equations, we derive an asymptotic criterion for baroclinic instability of compensated (i.e. confined to the upper layer) vortices. Surprisingly, it coincides exactly with a similar criterion derived from the quasigeostrophic equations [Benilov, E.S., 2003. Instability of quasigeostrophic vortices in a two-layer ocean with thin upper layer. J. Fluid Mech. 475, 303–331]. Thus, to leading order in , ageostrophy does not affect the stability properties of thin compensated vortices. As a result, whether a vortex is stable or not, depends on its shape, not amplitude (although the growth rate of an unstable vortex does depend on its amplitude).     

Application of a Bivariate Gamma Distribution for a Chemically Reacting Plume in the Atmosphere

The joint concentration probability density function of two reactive chemical species is modelled using a bivariate Gamma distribution coupled with a three-dimensional fluctuating plume model able to simulate the diffusion and mixing of turbulent plumes. A wind-tunnel experiment (Brown and Bilger, J Fluid Mech 312:373–407, 1996), carried out in homogeneous unbounded turbulence, in which nitrogen oxide is released from a point source in an ozone doped background and the chemical reactions take place in non-equilibrium conditions, is considered as a test case. The model is based on a stochastic Langevin equation reproducing the barycentre position distribution through a proper low-pass filter for the turbulence length scales. While the meandering large-scale motion of the plume is directly simulated, the internal mixing relative to the centroid is reproduced using a bivariate Gamma density function. The effect of turbulence on the chemical reaction (segregation), which in this case has not yet attained equilibrium, is directly evaluated through the covariance of the tracer concentration fields. The computed mean concentrations and the O₃–NO concentration covariance are also compared with those obtained by the Alessandrini and Ferrero Lagrangian single particle model (Alessandrini and Ferrero, Physica A 388:1375–1387, 2009) that entails an ad hoc parametrization for the segregation coefficient.     

Experiments with density currents on a sloping bottom in a rotating fluid

The behaviour of a density current on a sloping bottom in a rotating system is investigated by laboratory experiments. The main result is that the dense bottom outflow induces cyclonic vortices in the upper fluid layer, which are formed periodically and move to the west parallel to coast. Two regimes of vortex formation have been identified. For strong density currents and weak rotation, vortices are formed by stretching of the upper layer near the source as found also in the experiments by Lane-Serff and Baines (1998) [Lane-Serff, G.F., Baines, P.G., 1998. Eddy formation by dense flows on slopes in a rotating fluid. J. Fluid Mech. 363, 229–253]. For weak density currents and strong rotation vortices are due to instability of the bottom plume itself as found in the numerical simulations of Jiang and Garwood (1996) [Jiang, L., Garwood, W. Jr., 1996. Three-dimensional simulations of overflows on continental slopes. J. Phys. Oceanogr. 26, 1224–1233].     

Large-eddy Simulation of the Dispersion of Solid Particles in a Turbulent Boundary Layer

A large-eddy simulation (LES) with the dynamic Smagorinsky-Germano subgrid-scale (SGS) model is used to study the dispersion of solid particles in a turbulent boundary layer. Solid particles are tracked in a Lagrangian way. The instantaneous velocity of the surrounding fluid is considered to have a large-scale part (directly computed by the LES) and a small-scale part. The SGS velocity of the surrounding fluid is given by a three-dimensional Langevin model written in terms of SGS statistics at a mesh level. An appropriate Lagrangian correlation time scale is considered in order to include the influences of gravity and inertia of the solid particle. Inter-particle collisions and the influence of particles on the mean flow are also taken into account. The results of the LES are compared with the wind-tunnel experiments of Nalpanis et al. (1993 J Fluid Mech 251: 661–685) and of Tanière et al. (1997 Exp in Fluids 23:463–471) on sand particles in saltation and in modified saltation, respectively.     

A Lagrangian stochastic model for the trajectories of particle pairs and its application to the prediction of concentration variance within plant canopies

A simple Lagrangian stochastic model for the trajectories of particle pairs in high Reynolds-number turbulent flows is presented. In this model, the velocities of particle pairs are initially correlated but subsequently each particle moves independently. The independent single-particle trajectories are simulated using Thomson's model [J. Fluid Mech. 180, 529–556, 1987]. This two-particle model exactly satisfies the well-mixed condition for Gaussian turbulence when length scales, characterizing the two-point Eulerian velocity correlation function, vanish. Temperature variances, due to heat released as a passive scalar from an elevated plane source, within a model plant canopy (Coppin et al. Boundary Layer Meteorol. 35, 167–191, 1986) are shown to be well predicted by the model. It is suggested that for strongly inhomogeneous flows, the two-point Eulerian velocity function is of secondary importance in determining the simulated trajectories of particle pairs compared to the importance of ensuring satisfaction of the two-to-one constraint (Borgas and Sawford. J. Fluid Mech. 279, 69–99, 1994); i.e ensuring that one-particle statistics obtained from the two-particle model are the same as those obtained from the corresponding one-particle model. Limitations of this modelling approach are discussed.     

Atmospheric boundary-layer structure observed during a haze event due to forest-fire smoke

During a haze event in Baltimore, U.S.A. from July 6 to 8, 2002, smoke from forest fires in the Québec region (Canada), degraded air quality and impacted upon local climate, decreasing solar radiation and air temperature. The smoke particles in and above the atmospheric boundary layer (ABL) served as a tracer and provided a unique opportunity to investigate the ABL structure, especially entrainment. Elastic backscatter lidar measurements taken during the haze event distinctly reveal the downward sweeps (or wisps) of smoke-laden air from the free atmosphere into the ABL. Visualisations of mechanisms such as dry convection, the entrainment process, detrainment, coherent entrainment structures, and mixing inside the ABL, are presented. Thermals overshooting at the ABL top are shown to create disturbances in the form of gravity waves in the free atmosphere aloft, as evidenced by a corresponding ripple structure at the bottom of the smoke layer. Lidar data, aerosol ground-based measurements and supporting meteorological data are used to link free atmosphere, mixed-layer and ground-level aerosols. During the peak period of the haze event (July 7, 2002), the correlation between time series of elastic backscatter lidar data within the mixed layer and the scattering coefficient from a nephelometer at ground level was found to be high (R=0.96 for z =324 m, and R=0.89 for z=504 m). Ground-level aerosol concentration was at a maximum about 2 h after the smoke layer intersected with the growing ABL, confirming that the wisps do not initially reach the ground.     

修正的质量通量积云对流方案及其模拟试验研究Ⅰ:方案介绍及对1991年洪涝过程的模拟

文中在综合比较各类积云对流参数化方案优缺点的基础上 ,主要参考陈伯民等修正的ECMWF质量通量积云参数化方案 ,对其进行简化和修改 ,发展了一个质量通量积云对流参数化方案 ,文中表示为 MFS(Mass Flux Scheme)。MFS是一种综合型的方案 ,既考虑了大尺度水汽辐合的重要性 ,又考虑了积云中的上升运动、下沉运动、环境中的补偿下沉运动 ,以及卷入、卷出和蒸发等 ,用总体云模式来描述积云与环境的相互作用 ,同时考虑了深对流和浅对流。将 MFS植入 NCAR区域气候模式 Reg CM2中 ,对 1 991年 5～ 7月江淮地区特大降水过程的夏季风气候特征和变化进行了模拟 ,并与 NCAR Reg CM2选用 Kuo方案 ,在同样初、边值条件和其它物理过程选择下的模拟结果进行了对比分析。分析结果表明 ,植入 MFS后的模式能够模拟这次极端的降水气候事件。在某些方面 ,如地表气温 ,降水的模拟上 ,植入 MFS后的模式的模拟结果要比原模式的结果更合理     

Diurnal cycle of precipitation over the Maritime Continent simulated by a spectral cumulus parameterization

The validity of a spectral cumulus parameterization (spectral scheme) for simulating a diurnal cycle of precipitation over the Maritime Continent (MC) was examined using a regional atmospheric model. The impacts of entrainment parameterization and each type of convective closure, i.e., non-equilibrium (or equilibrium) closure for deep convection, mid-level, and shallow convective closures, were also examined. When vertically variable entrainment and appropriate convective closures were employed, the model adequately simulated a diurnal cycle of precipitation over both land and ocean as compared to the observation. Analysis regarding the entrainment parameterization revealed that variable entrainment parameterization was needed not only for simulating better mean patterns of precipitation, but also for more realistic phases of diurnal cycles. The impacts of convective closures appeared in the differences in the precipitation amplitude. Analysis on diurnal cycles of convective properties and tendencies revealed that the cycles between boundary layer forcing and convective heating determined convection strength and were affected by each type of convective closure. It can be concluded that the spectral scheme with appropriate convective closures is able to simulate a realistic diurnal cycle over the MC.     

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