Sensitivity of climate simulations to the parameterization of cumulus convection in the Canadian climate centre general circulation model

Abstract

A simplified cumulus parameterization scheme, suitable for use in GCMs, is presented. This parameterization is based on a plume ensemble concept similar to that originally proposed by Arakawa and Schubert (1974). However, it employs three assumptions which significantly simplify the formulation and implementation of the scheme. It is assumed that an ensemble of convective‐scale updrafts with associated saturated downdrafts may exist when the atmosphere is locally conditionally unstable in the lower troposphere. However, the updraft ensemble is comprised only of those plumes which are sufficiently buoyant to penetrate through this unstable layer. It is assumed that all such plumes have the same upward mass flux at the base of the convective layer. The third assumption is that moist convection, which occurs only when there is convective available potential energy (CAPE) for reversible ascent of an undiluted parcel from the sub‐cloud layer, acts to remove CAPE at an exponential rate with a specified adjustment time scale.

The performance of the scheme and its sensitivity to choices of disposable parameters is illustrated by presenting results from a series of idealized single‐column model tests. These tests demonstrate that the scheme permits establishment of a quasi‐equilibrium between large‐scale forcing and convective response. However, it is also shown that the strength of convective downdrafts is an important factor in determining the nature of the equilibrium state. Relatively strong down‐drafts give rise to an unsteady irregularly fluctuating state characterized by alternate periods of deep and shallow convection.

The effect of using the scheme for GCM climate simulations is illustrated by presenting selected results of a multi‐year simulation carried out with the Canadian Climate Centre GCM using the new parameterization (the CONV simulation). Comparison of these results with those for a climate simulation made with the standard model (the CONTROL simulation, as documented by McFarlane et al., 1992) reveals the importance of other parameterized processes in determining the ultimate effect of introducing the new convective scheme. The radiative response to changes in the cloudiness regime is particularly important in this regard.     

Design criteria for water tank models of dispersion in the planetary convective boundary layer

Design criteria for laboratory water-analogs of clear-air penetrative convection in the atmosphere are described. Consideration is given to the range of factors relevant to modelling both turbulent penetrative convection and the dispersion of buoyant point-source plumes within the convective boundary layer. Scaling arguments based on mixed-layer and plume scaling show that at typical laboratory scales, saline convection can satisfy the requirements for modelling buoyant plume dispersion under strongly convective (light wind) conditions better than heated water tanks or wind tunnels.     

A saline laboratory model of the planetary convective boundary layer

A laboratory water-analog of clear-air penetrative convection in the atmosphere has been constructed to continue studies of the turbulent dispersion of buoyant plumes in the convective boundary layer (CBL). A unique feature is the utilization of saline rather than thermal convection, which has been made possible by the development of a reliable method for delivering a controllable buoyancy flux through a porous membrane. It has been shown in an earlier paper that at typical laboratory scales, a saline convection tank is well suited to modelling buoyant plume dipersion under strongly convective (light wind) conditions.A range of experiments has clearly demonstrated the validity of the model. Results for density and velocity variances show much less scatter than most comparable measurements because of the greatly improved sampling that is possible in the tank. The results are generally in good agreement with field data and other laboratory simulations but the improved accuracy of the data has highlighted the anomalously low values for the horizontal velocity variances produced by large-eddy simulations of the CBL. The cause of this apparent underprediction remains unresolved.     

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

考虑动量和浮力通量耗散的烟气抬升解析模式

环流湍流对烟流抬升有重要的作用。本文考虑了环境湍流引起的烟流的动量和浮力通量耗散，导出了烟流抬升的控制方程，给出了烟流的轨迹方程和中性条件下热浮力烟流的终极抬升高度。在烟囱附近，轨迹方程接近2/3次律；当下风距离增大时，与2/3次律有明显偏离并逐渐变平。与外场试验资料的比较表明，本模式能较好地模拟浮力烟流的轨迹和终极高度。     

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.     

An improvement of the mass flux convection parameterization scheme and its sensitivity tests for seasonal prediction over China

A modified cumulus parameterization scheme, suitable for use in a seasonal forecast model, is presented. This parameterization scheme is an improvement of the mass flux convection scheme developed by Gregory and Rowntree (1989; 1990). This convection scheme uses a “bulk” cloud model to present an ensemble of convective clouds, and aims to represent shallow, deep, and mid-level convection. At present,this convection scheme is employed in the NCC T63L20 model (National Climate Center, China Meteorological Administration). Simulation results with this scheme have revealed some deficiencies in the scheme,although to some extent, it improves the accuracy of the simulation. In order to alleviate the deficiencies and reflect the effect of cumulus convection in the actual atmosphere, the scheme is modified and improved.The improvements include (i) the full estimation of the effects of the large-scale convergence in the lower layer upon cumulus convection, (ii) the revision of the initial convective mass flux, and (iii) the regulation of convective-scale downdrafts. A comparison of the results obtained by using the original model and the modified one shows that the improvement and modification of the original convection scheme is successful in simulating the precipitation and general circulation field, because the modified scheme provides a good simulation of the main features of seasonal precipitation in China, and an analysis of the anomaly correlation eoetfieient between the simulation and the observations confirms the improved results.     

THE LIMITATION OF CLOUD-BASE MASS FLUX IN CUMULUS PARAMETERIZATION AND ITS APPLICATION IN A HIGH-RESOLUTION MODEL

A large area of unrealized precipitation is produced with the standard convective parameterization scheme in a high-resolution model, while subgrid-scale convection that cannot be explicitly resolved is omitted without convective parameterization. A modified version of the convection scheme with limited mass flux at cloud base is introduced into a south-China regional high-resolution model to alleviate these problems. A strong convection case and a weak convection case are selected to analyze the influence of limited cloud-base mass flux on precipitation forecast. The sensitivity of different limitation on mass flux at cloud base is also discussed. It is found that using instability energy closure for Simplified Arakawa- Schubert Scheme will produce better precipitation forecast than the primary closure based on quasi-equilibrium assumption. The influence of the convection scheme is dependent on the upper limit of mass flux at cloud base. The total rain amount is not so sensitive to the limitation of mass flux in the strong convection case as in the weak one. From the comparison of two different methods for limiting the cloud-base mass flux, it is found that shutting down the cumulus parameterization scheme completely when the cloud-base mass flux exceeds a given limitation is more suitable for the forecast of precipitation.     

Betts对流调整方案在移动套网格台风数值模式中的初步试验

Betts对流调整方案采用的是垂直温湿场向观测的准平衡热力结构松驰的方法。本文在一个考虑了台风初值化的移动套网格台风数值模式中,对Betts方案进行了初步试验。结果表明：Betts方案能较好的预报出台风移动的路径和降水过程,并在台风路径、降水落区预报以及对流增温的垂直分布等方面都略优于Kuo方案的预报结果。     

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

Penetrative convection in rapidly rotating flows: preliminary results from numerical simulation

Turbulent convection forced by a surface heat flux into a stably stratified region is a feature of both the atmospheric and oceanic planetary boundary layers. Of particular interest is the interface between the convective layer and the stable stratification, where the entrainment of fluid into the convective layer by penetrating plumes may lead to a reverse buoyancy flux, and an enhancement of the stable stratification. Whereas in the atmosphere the influence of rotation on this penetrative convection is negligible, oceanic convection may be subjected to lower Rossby numbers and hence greater rotational influence. To isolate the effects of rotation, we present three numerical solutions for turbulent penetrative convection, characterised by different rotation rates, with all other parameters being held constant. Our results indicate that at lower Rossby numbers the lateral scale of the plumes is reduced, whereas the vertical vorticity of the plumes is much enhanced. Vertical transports of buoyancy and kinetic energy across the convective layer are reduced, leading to less efficient penetration at the interface with the stratified layer, and hence less reverse buoyancy flux in this region.     

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.     

Convection model for stratus cloud over a warm water surface

A numerical model is developed to simulate convective stratus cloud formation over the sea. The model is based on quasi-steady state moist plumes advecting over an area of increasing sea-surface temperature with the sea warmer than the air, but it is also a good approximation for non-steady states whenever the model is matched to air-sea temperature differences. Combining the effects of upward transfer of heat and moisture fluxes as well as adiabatic cooling, stratus cloud forms and spreads downward in this field. The depth of the convective field, the sea surface temperature gradient, the liquid water content at cloud top, and the horizontal pressure gradient are the four controlling parameters for the convective field. Alternatively, the wind speed, air-sea temperature difference, and the mixing ratio of the air, derived variables in this treatment, can be taken, with the depth, as the basic parameters.The entity type of convective model used here has the advantage that it models the transport and modification of air parcels and hence provides a method for studying drop size development in stratus clouds. It also uses relationships derived from water tank experiments with plumes and tested in dry convection and so needs no parameters specific to each situation. Its most important feature, however, is that the mean motion of plumes, rather than turbulent diffusion, transports the moisture.The upward growth and erosion of a temperature step increase to produce an inversion can be attributed to liquid water present in haze drops or cloud drops, rather than to turbulent diffusion. Radiative transfer is not a necessary requirement, but may either enhance or slow down the process.     

A bulk mass flux convection scheme for climate model: description and moisture sensitivity

A convection scheme for climate model is developed based on Tiedtke’s (Mon Weather Rev 117:1779–1800, 1989) bulk mass flux framework and is evaluated with observational data and cloud resolving model simulation data. The main differences between the present parameterization and Tiedtke’s parameterization are the convection trigger, fractional entrainment and detrainment rate formulations, and closure method. Convection is triggered if the vertical velocity of a rising parcel is positive at the level at which the parcel is saturated. The fractional entrainment rate depends on the vertical velocity and buoyancy of the parcel as well as the environmental relative humidity. For the fractional detrainment rate, a linear decrease in the updraft mass flux above maximum buoyancy level is assumed. In the closure method, the cloud base mass flux is determined by considering both cloud layer instability and subcloud layer turbulent kinetic energy as controlling factors in the strength of the convection. The convection scheme is examined in a single column framework as well as using a general circulation model. The present bulk mass flux (BMF) scheme is compared with a simplified Relaxed Arakawa-Schubert (RAS) scheme. In contrast to the RAS, which specifies the cloud top, cloud top height in BMF depends on environmental properties, by considering the conditions of both the parcel and its environment in a fractional entrainment and detrainment rate formulations. As a result, BMF shows improved sensitivity in depth and strength of convection on environmental humidity compared to RAS, by strengthening coupling between cloud and environment. When the mid to lower troposphere is dry, the cloud resolving model and BMF produce cloud top around the dry layer and moisten the layer. In the framework of general circulation model, enhanced coupling between convection and environmental humidity in BMF results in improved representation of eastward propagating intraseasonal variability in the tropics—the Madden-Julian oscillation.     

The rise of forced plumes in a stably stratified crossflow

Solutions are found describing the rise of turbulent buoyant plumes in a neutral or stably stratified crossflow from a source having finite fluxes of buoyancy, momentum and volume. Plumes from such sources are known as forced plumes. The solutions arise from a unified theory having one disposable parameter, the entrainment constant, and they describe a wide range of behaviour of plumes and jets in a crossflow. Solutions for buoyant plumes indicate that, for any given source, heights of rise are reduced with increased crossflow velocity and increased stratification while for constant environmental parameters, heights of rise are increased with increased fluxes of source buoyancy and momentum and are decreased with source radius. Maximum dilutions occur in plumes from small sources with relatively large buoyancy fluxes in light crossflows with small stratification.     

半拉格朗日、半隐式欧拉方程组大气数值模式研究

介绍了一个具有较高时间积分效率的三维弹性大气数值模式。其中,控制方程为能够描述大气非静力及可压缩性的欧拉方程组,时间积分采用了高效率的“半拉格朗日、半隐式”方案。通过引用一个简化的“云物理过程”参数化方案,针对大气中深厚湿对流过程的若干问题进行了数值试验,结果表明,模式能够对该过程中的一些基本和复杂现象进行有效地模拟。试验结果还表明,动力学框架具有很好的稳定性,能够实现高效率的时间积分,它还具有较好的频散特性,保证了空间计算的精度,从而完好地描述了风暴的形态。     

A Simulation Study on the Characteristics of Cloud Microphysics of Heavy Rainfall in the Meiyu Front

A heavy rainfall in the Meiyu front during 4--5 July 2003 is simulated by use of the non-hydrostatic mesoscale model MM5 (V3--6) with different explicit cloud microphysical parameterization schemes. The characteristics of microphysical process of convective cloud are studied by the model outputs. The simulation study reveals that: (1) The mesoscale model MM5 with explicit cloud microphysical process is capable of simulating the instant heavy rainfall in the Meiyu front, the rainfall simulation could be improved significantly as the model resolution is increased, and the Goddard scheme is better than the Reisner or Schultz scheme. (2) The convective cloud in the Meiyu front has a comprehensive structure composed of solid, liquid and vapor phases of water, the mass density of water vapor is the largest one in the cloud; the next one is graupel, while those of ice, snow, rain water and the cloud water are almost same. The height at which mass density peaks for different hydrometeors is almost unchangeable during the heavy rainfall period. The mass density variation of rain water, ice, and graupel are consistent with that of ground precipitation, while that of water vapor in the low levels is 1--2 h earlier than the precipitation. (3) The main contribution to the water vapor budget in the atmosphere is the convergence of vapor flux through advection and convection, which provides the main vapor source of the rainfall. Besides the basic process of the auto-conversion of cloud water to rain water, there is an additional cloud microphysical process that is essential to the formation of instant heavy rainfall, the ice-phase crystals are transformed into graupels first and then the increased graupels mix with cloud water and accelerates the conversion of cloud water to rain water. The positive feedback mechanism between latent heat release and convection is the main cause to maintain and develop the heavy precipitation.     

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

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