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

LMDZ5B: the atmospheric component of the IPSL climate model with revisited parameterizations for clouds and convection

Based on a decade of research on cloud processes, a new version of the LMDZ atmospheric general circulation model has been developed that corresponds to a complete recasting of the parameterization of turbulence, convection and clouds. This LMDZ5B version includes a mass-flux representation of the thermal plumes or rolls of the convective boundary layer, coupled to a bi-Gaussian statistical cloud scheme, as well as a parameterization of the cold pools generated below cumulonimbus by re-evaporation of convective precipitation. The triggering and closure of deep convection are now controlled by lifting processes in the sub-cloud layer. An available lifting energy and lifting power are provided both by the thermal plumes and by the spread of cold pools. The individual parameterizations were carefully validated against the results of explicit high resolution simulations. Here we present the work done to go from those new concepts and developments to a full 3D atmospheric model, used in particular for climate change projections with the IPSL-CM5B coupled model. Based on a series of sensitivity experiments, we document the differences with the previous LMDZ5A version distinguishing the role of parameterization changes from that of model tuning. Improvements found previously in single-column simulations of case studies are confirmed in the 3D model: (1) the convective boundary layer and cumulus clouds are better represented and (2) the diurnal cycle of convective rainfall over continents is delayed by several hours, solving a longstanding problem in climate modeling. The variability of tropical rainfall is also larger in LMDZ5B at intraseasonal time-scales. Significant biases of the LMDZ5A model however remain, or are even sometimes amplified. The paper emphasizes the importance of parameterization improvements and model tuning in the frame of climate change studies as well as the new paradigm that represents the improvement of 3D climate models under the control of single-column case studies simulations.     

The role of thermals in the convective boundary layer

Detailed measurements of the structure of thermals throughout the convective boundary layer were obtained from the NCAR Electra aircraft over the ocean during the Air Mass Transformation Experiment (AMTEX). Humidity was used as an indicator of thermals. The variables were first high-pass filtered with a 5 km cutoff digital filter to eliminate mesoscale variations. Segments of the 5 min (30 km length) horizontal flight legs with humidity greater than half the standard deviation of humidity fluctuations for that leg were defined as thermals. This was found to be a better indicator of thermals than temperature in the upper part of the boundary layer since the temperature in a thermal is cooler than its environment in the upper part of the boundary layer. Using mixed-layer scaling, the normalized length and number of thermals were found to scale with the 1/3 and -1/3 powers, respectively, of normalized height, while vertical velocity and temperature scaled according to similarity predictions in the free convection region of the surface layer. The observational results presented here extend throughout the entire mixed layer. Using these results in the equation for mean updraft velocity of a field of thermals, the sum of the vertical pressure gradient and edge-effect terms can be estimated. This residual term is found to be important throughout most of the boundary layer. The magnitude of the divergence of vertical velocity variance within a thermal is found to be larger than the magnitude of the mean updraft velocity term throughout most of the mixed layer.Part of this work was completed while visiting Risø National Laboratory, Denmark.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

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.     

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

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

An ensemble of dynamically identical thermals and vertical profiles of turbulent moments in the convective surface layer of atmosphere

An ensemble of convective thermals is distinguished from the surface layer of penetrative turbulent convection over a heated horizontally uniform surface. For an isolated convective element, an integral model of a quasi-stationary spontaneous jet is proposed which admits an exact analytical solution. A simple statistical model is constructed for an ensemble of dynamically identical thermals. In this model, convective thermals ascend in a static environment, their dynamic parameters are described by the equations of an isolated quasi-stationary jet, and their diameters are stochastic. It is shown that the ensemble of thermals rising in a horizontally homogeneous environment forms surface-layer turbulent moments. The analytical relationships for higher turbulent moments of vertical velocity and temperature are compared with experimental data from the second to the forth order inclusive.     

Resolved Versus Parametrized Boundary-Layer Plumes. Part II: Continuous Formulations of Mixing Rates for Mass-Flux Schemes

The conditional sampling of coherent structures in large-eddy simulations of the convective boundary layer (Couvreux et al. Boundary-layer Meteorol 134:441–458, 2010) is used to propose and evaluate formulations of fractional entrainment and detrainment rates for mass-flux schemes. The proposed formulations are physically-based and continuous from the surface to the top of clouds. Entrainment is related to the updraft vertical velocity divergence, while detrainment depends on the thermal vertical velocity, on buoyancy and on the moisture contrast between the mean plume and its environment. The proposed formulations are first directly evaluated in simulations of shallow clouds. They are then tested in single-column simulations with the thermal plume model, a mass-flux representation of boundary-layer thermals.     

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.     

冬、夏季热带及副热带穿透性对流气候特征分析

文中利用热带测雨卫星 (TRMM) 搭载的测雨雷达 (PR) 1998~2007年的探测结果, 就热带及副热带地区穿透性对流的频次、条件降水强度及垂直廓线等特征进行了分析。研究结果表明: 深对流和穿透性对流都主要发生在热带辐合带 ( ITCZ)、南太平洋辐合带 (SPCZ)、亚洲季风区、20°N以南的非洲以及美洲等地区, 它们的空间分布具有明显的地域性和季节变化特征, 而且陆地深对流更容易发展成为穿透性对流, 但绝大部分地区的穿透性对流频次不超过0.2%。对穿透性对流条件降水强度的分析表明, 热带及副热带大部分地区的穿透性对流条件降水强度在10 mm/h以上, 且洋面的条件降水强度要比陆地大, 但由于其频次较小导致其对总降水的贡献并不大。尽管深对流和穿透性对流降水廓线的外形比较相似, 但相同的高度, 深对流的降水强度要比穿透性对流偏小, 而且这种差异随海陆和纬度的不同而有所区别。此外, 热带地区 (15°S~15°N) 冬、夏季深对流和穿透性对流降水廓线都只存在较小差异, 并没有显示出明显的季节变化。     

The Moist Boundary Layer under a Mid-latitude Weather System

Mid-latitude weather systems are key contributors to the transport of atmospheric water vapour, but less is known about the role of the boundary layer in this transport. We expand a conceptual model of dry boundary-layer structure under synoptic systems to include moist processes, using idealised simulations of cyclone waves to investigate the three-way interaction between the boundary layer, atmospheric moisture and large-scale dynamics. Forced by large-scale thermal advection, boundary-layer structures develop over large areas, analogous to the daytime convective boundary layer, the nocturnal stable boundary layer and transitional regimes between these extremes. A budgeting technique demonstrates the key role of boundary-layer processes in the transport of moisture. Moisture is evaporated from the ocean behind the cold front and in the high-pressure part of the wave, and transported large distances within the boundary layer into the footprint of the warm-conveyor belt. The warm-conveyor belt forms one of the two main processes of boundary-layer ventilation, with shallow cumulus convection being of similar importance.     

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

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

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.     

台风生成阶段周期性对流爆发机制及辐射的影响

最新观测统计分析表明台风生成前对流活动存在周期性特征,周期集中在18~26 h,峰值多出现在夜间。本文通过一系列理想的云分辨数值控制试验和辐射敏感性对比试验讨论了台风生成阶段周期性对流爆发的成因以及辐射日变化在此过程中产生的影响。结果表明,在给定无日变化的恒定太阳辐射情况下,台风生成过程中也会出现与观测一致的周期性对流爆发,说明辐射日变化不是该周期性形成的控制因素。对湿静力能扰动的收支分析发现边界层中湿静力能的周期性增减是导致周期性深对流爆发的关键因素。深对流爆发形成的冷池通过平流作用降低了边界层内的湿静力能,需要通过辐射与海表面通量过程使其恢复,才能再次触发对流。辐射日变化对于一些较弱的热带低压涡旋有调制作用,能使其周期性对流爆发的峰值时间相位改变并趋于在夜晚出现;但对于较强的涡旋,辐射的这种调制作用并不显著,涡旋会保持其最初对流爆发的周期和相位。本研究对观测到的台风生成阶段对流峰值出现在白天或夜间的情况都可以作出合理解释,也深化了对台风生成过程对流和辐射的作用和影响科学问题的理解。     

Bivariate conditional sampling of buoyancy flux during an intense cold-air outbreak

Buoyancy fluxes in the marine atmospheric boundary layer (MABL) for the cloud street regime, observed during the Genesis of Atlantic Lows Experiment (GALE), have been analyzed using the technique of joint frequency distribution. For the lower half of the MABL, the results suggest that the buoyancy flux is mainly generated by the rising thermals and the sinking compensating ambient air, and is mainly consumed by the entrainment and detrainment of thermals, penetrative convection, and the entrainment from the MABL top.The results are compared to those from previous studies of mesoscale cellular convection (Air-Mass Transformation Experiment, AMTEX), the dry convective boundary layer, and the trade-wind MABL. For the lower MABL, the quadrant buoyancy fluxes, fractional coverages, and flux intensities are in good agreement with those of mesoscale cellular convection (AMTEX) and the dry convective boundary layer. The results suggest that, if the buoyancy flux is primarily driven by the temperature flux, the physical processes for generating buoyancy flux mentioned above are about the same for the lower boundary layers over land and ocean, even with different convective regimes. For the trade-wind MABL, the buoyancy flux is mainly driven by the moisture flux; the quadrant flux intensities are stronger than those of the other three studies except for the buoyant updrafts (thermals). These results suggest that the entrainment and detrainment of thermals are more effective in the trade-wind MABL than in the boundary layers driven by the temperature flux.Scale analysis of the buoyancy flux is in good agreement with that of AMTEX. For the lower half of the MABL, the buoyancy flux is mainly generated by the intermediate scale (200 m to 2 km), which includes the dominant convective thermals in the surface layer and the mixed layer. The scale smaller than 200 m is important only in the surface layer. The scale larger than 2 km, which includes the roll vortices, increases its significance upward. While most of the positive and negative fluxes are associated with the updrafts for the intermediate scale, the downdrafts are as important as updrafts for the larger scale.ST Systems Corporation, Lanham, MD, 20706, U.S.A.     

Boundary-layer evolution and nocturnal inversion dispersal—Part II

Observations from an instrumented aircraft are used to study the small-scale structure of turbulence and convection in well-mixed boundary layers and the erosion processes in the nocturnally-formed inversions above them. The ways in which turbulence statistics for temperature, humidity and vertical velocity scale with height in the mixed layer are compared with the results of a three-dimensional model by Deardorff (1974a, b), and agreement is found in many aspects. Conditional sampling enables the statistics of thermals and their environment to be considered separately and, in particular, shows that the mode of the vertical velocity in thermals markedly decreases with height in the upper half of the mixed layer. Thermals may be recognized equally readily by either their excess of temperature or humidity. Transfers of heat and moisture through the nocturnal inversions influence the structure of the upper region of the mixed layer and there is strong evidence that these transfer processes are turbulent and not organized on scales similar to convective thermals.     

A mass flux closure function in a GCM based on the Richardson number

A mass flux closure in a general circulation model (GCM) was developed in terms of the mean gradient Richardson number (GRN), which is defined as the ratio between the buoyancy and the shear-driven kinetic energy in the planetary boundary layer. The cloud resolving model (CRM) simulations using the tropical ocean and global atmosphere-coupled ocean–atmosphere response experiment forcing show that cloud-base mass flux is well correlated with the GRN. Using the CRM simulations, a mass flux closure function is formulated as an exponential function of the GRN and it is implemented in the Arakawa–Schubert convective scheme. The GCM simulations with the new mass flux closure are compared to those of the GCM with the conventional mass flux closure based on convective available potential energy. Because of the exponential function, the new closure permits convective precipitation only when the GRN has a sufficiently large value. When the GRN has a relatively small value, the convection is suppressed while the convective instability is released by large-scale precipitation. As a result, the ratio of convective precipitation to total precipitation is reduced and there is an increase in the frequency of heavy precipitation, more similar to the observations. The new closure also improves the diurnal cycle of precipitation due to a time delay of the large GRN with respect to convective instability.     

一次热带海洋对流云微物理过程的数值模拟

为了理解微物理过程以及次生对流云对热带对流云团发展演变过程的影响,利用二维云分辨模式对2006年1月20日靠近澳大利亚西北部的印度尼西亚暖池附近的一次对流云降水过程进行数值模拟,分析了对流云发展演变不同阶段各种微物理过程的特征,尤其对次生对流的形成以及并合过程进行了深入分析.结果表明:本个例所研究的热带对流云团具有发展...     

20世纪北大西洋温盐环流的年代际变化试评估

根据相对丰富的大气器测资料,综合前人对有限的海洋资料的诊断分析,从北大西洋涛动(NAO)变率、表层海温(SST)变率、格陵兰海和拉布拉多海的深对流活动长期变化等不同角度,对20世纪大洋温盐环流(Thermohaline Circulation,THC)变率进行了试评估.结果表明:(1)19世纪末以来,大西洋温盐环流的变化可分为4个时期:1900年以前的一段时期,THC较强;1904年到1930年,THC较弱;1931年到1972年,THC较强;1973年至1995年,THC较弱,目前则又有所增强.(2)与THC的变化相联系,大西洋主要气候要素的变化,相互间存在着某种协调关系,THC强,NAO弱,北大西洋北部SST升高,格陵兰海的对流活动增强,拉布拉多海的对流活动则减弱.     

Cloud-radiation feedback and atmosphere-ocean coupling in a stochastic multicloud model

Despite recent advances in supercomputing, current general circulation models (GCMs) have significant problems in representing the variability associated with organized tropical convection. Furthermore, due to high sensitivity of the simulations to the cloud radiation feedback, the tropical convection remains a major source of uncertainty in long-term weather and climate forecasts. In a series of recent studies, it has been shown, in paradigm two-baroclinic-mode systems and in aquaplanet GCMs, that a stochastic multicloud convective parameterization based on three cloud types (congestus, deep and stratiform) can be used to improve the variability and the dynamical structure of tropical convection, including intermittent coherent structures such as synoptic and mesoscale convective systems. Here, the stochastic multicloud model is modified with a parameterized cloud radiation feedback mechanism and atmosphere-ocean coupling. The radiative convective feedback mechanism is shown to increase the mean and variability of the Walker circulation. The corresponding intensification of the circulation is associated with propagating synoptic scale systems originating inside of the enhanced sea surface temperature area. In column simulations, the atmosphere ocean coupling introduces pronounced low frequency convective features on the time scale associated with the depth of the mixed ocean layer. However, in the presence of the gravity wave mixing of spatially extended simulations, these features are not as prominent. This highlights the deficiency of the column model approach at predicting the behavior of multiscale spatially extended systems. Overall, the study develops a systematic framework for incorporating parameterized radiative cloud feedback and ocean coupling which may be used to improve representation of intraseasonal and seasonal variability in GCMs.     

The effect of surface friction on the development of tropical cyclones

When tropical cyclones (hereafter referred as TCs) are over the ocean, surface friction plays a dual role in the development of TCs. From the viewpoint of water vapor supply, frictional convergence and Ekman pumping provide a source of moisture for organized cumulus convection and is propitious to the spin-up of TCs. On the other hand, surface friction leads to a dissipation of kinetic energy that impedes the intensification of TCs. Which role is dominant in the developing stage of TCs is a controversial issue. In the present work, the influence of surface friction on the growth of TCs is re-examined in detail by conducting two sets of numerical experiments initialized with different cyclonic disturbances. Results indicate that, because of the inherent complexities of TCs, the impact of surface friction on the evolution of TCs can not be simply boiled down to being positive or negative. In the case that a TC starts from a low-level vortex with a warm core, surface friction and the resultant vertical motion makes an important contribution to the convection in the early developing stage of the TC by accelerating the build-up of convective available potential energy (CAPE) and ensuring moisture supply and the lifting of air parcels. This effect is so prominent that it dominates the friction-induced dissipation and makes surface friction a facilitative factor in the spin-up of the TC. However, for a TC formed from a mesoscale convective vortex (MCV) spawned in a long-lasting mesoscale convective system (MCS), the initial fields, and especially the low-level humidity and cold core, enable the prerequisites of convection (i.e., conditional instability, moisture, and lifting), to be easily achieved even without the help of boundary-layer pumping induced by surface friction. Accordingly, the reliance of the development of TCs on surface friction is not as heavy as that derived from a low-level vortex. The positive effect of surface friction on the development of TCs realized through facilitating favorable conditions for convection is nearly cancelled out by the friction-induced dissipation. However, as SST is enhanced in the latter case, the situation may be changed, and different development speeds may emerge between model TCs with and without surface friction considered. In short, owing to the fact that TC development is a complicated process affected by many factors such as initial perturbations, SST, etc., the importance of surface friction to the intensification of TCs may vary enormously from case to case.