Zhang-McFarlane Scheme Research Articles

Wave‐Convection Interactions Amplify Convective Parameterization Biases in the South Pacific Convergence Zone

AbstractClimate models have long‐standing difficulties simulating the South Pacific Convergence Zone (SPCZ) and its variability. For example, the default Zhang‐McFarlane (ZM) convection scheme in the Community Atmosphere Model version 5 (CAM5) produces too much light precipitation and too little heavy precipitation in the SPCZ, with this bias toward light precipitation even more pronounced in the SPCZ than in the tropics as a whole. Here, we show that implementing a recently developed convection scheme in the CAM5 yields significant improvements in the simulated SPCZ during austral summer and discuss the reasons behind these improvements. In addition to intensifying both mean rainfall and its variability in the SPCZ, the new scheme produces a larger heavy rainfall fraction that is more consistent with observations and state‐of‐the‐art reanalyses. This shift toward heavier, more variable rainfall increases both the magnitude and altitude of diabatic heating associated with convective precipitation, intensifying lower tropospheric convergence and increasing the influence of convection on the upper‐level circulation. Increased diabatic production of potential vorticity in the upper troposphere intensifies the distortion effect exerted by convection on transient Rossby waves that pass through the SPCZ. Weaker distortion effects in simulations using the ZM scheme allow waves to propagate continuously through the region rather than dissipating locally, further reducing updrafts and weakening convection in the SPCZ. Our results outline a dynamical framework for evaluating model representations of tropical–extratropical interactions within the SPCZ and clarify why convective parameterizations that produce “top‐heavy” profiles of deep convective heating better represent the SPCZ and its variability.

Dynamically computed characteristic adjustment time scale for Zhang–McFarlane convective parameterization scheme

AbstractCharacteristic adjustment time scale τ is an important parameter in convective parameterization scheme. In this paper, we propose a dynamic calculation formula of τ in the Zhang–McFarlane scheme based on convective available potential energy (CAPE) in order to improve the simulation of precipitation. The performance of using the new formula is tested and evaluated using the Community Atmosphere Model (CAM5) for numerical experiments. The results indicate that the root mean square errors of mean precipitation amount over the Northern Hemisphere eastern equatorial Pacific Ocean, the central Indo-Pacific Warm Pool, the Tibetan Plateau, and northern South America are reduced by 20.15, 16.40, 12.56, and 6.45% in summer, respectively, after a dynamic τ is adopted. The dynamically computed τ improves the Pacific Walker circulation by strengthening the western ascending and eastern descending, and increases the intensity of deep convective precipitation in tropical areas by at least 1 mm/day. Dynamic τ also increases the frequency of large deep convective precipitation in the equatorial Pacific by almost 50%, and decreases the frequency of light deep convective precipitation in many other regions. The dynamic τ has positive effects in the areas with large CAPE and spatiotemporal variability, and also leads to an increase in deep convective precipitation and a significant decrease in large-scale precipitation over high-altitude regions of windward slopes. Therefore, the generation of exaggerated precipitation in the highlands is reduced in the CAM5. Overall, the dynamic calculation of characteristic adjustment time scale can improve precipitation simulation.

A Shallow‐Deep Unified Stochastic Mass Flux Cumulus Parameterization in the Single Column Community Climate Model

AbstractCumulus parameterization (CP) in state‐of‐the‐art global climate models is based on the quasi‐equilibrium assumption (QEA), which views convection as the action of an ensemble of cumulus clouds, in a state of equilibrium with respect to a slowly varying atmospheric state. This view is not compatible with the organization and dynamical interactions across multiple scales of cloud systems in the tropics and progress in this research area was slow over decades despite the widely recognized major shortcomings. Novel ideas on how to represent key physical processes of moist convection‐large‐scale interaction to overcome the QEA have surged recently. The stochastic multicloud model (SMCM) CP in particular mimics the dynamical interactions of multiple cloud types that characterize organized tropical convection. Here, the SMCM is used to modify the Zhang‐McFarlane (ZM) CP by changing the way in which the bulk mass flux and bulk entrainment and detrainment rates are calculated. This is done by introducing a stochastic ensemble of plumes characterized by randomly varying detrainment level distributions based on the cloud area fraction of the SMCM. The SMCM is here extended to include shallow cumulus clouds resulting in a unified shallow‐deep CP. The new stochastic multicloud plume CP is validated against the control ZM scheme in the context of the single column Community Climate Model of the National Center for Atmospheric Research using data from both tropical ocean and midlatitude land convection. Some key features of the SMCM CP such as it capability to represent the tri‐modal nature of organized convection are emphasized.

Implementation and Evaluation of a Double-Plume Convective Parameterization in NCAR CAM5

AbstractPerformance of global climate models (GCMs) is strongly affected by the cumulus parameterization (CP) used. Similar to the approach in GFDL AM4, a double-plume CP, which unifies the deep and shallow convection in one framework, is implemented and tested in the NCAR Community Atmospheric Model version 5 (CAM5). Based on the University of Washington (UW) shallow convection scheme, an additional plume was added to represent the deep convection. The shallow and deep plumes share the same cloud model, but use different triggers, fractional mixing rates, and closures. The scheme was tested in single-column, short-term hindcast, and AMIP simulations. Compared with the default combination of the Zhang–McFarlane scheme and UW scheme in CAM5, the new scheme tends to produce a top-heavy mass flux profile during the active monsoon period in the single-column simulations. The scheme increases the intensity of tropical precipitation, closer to TRMM observations. The new scheme increased subtropical marine boundary layer clouds and high clouds over the deep tropics, both in better agreement with observations. Sensitivity tests indicate that regime-dependent fractional entrainment rates of the deep plume are desired to improve tropical precipitation distribution and upper troposphere temperature. This study suggests that a double-plume approach is a promising way to combine shallow and deep convections in a unified framework.

Effects of coupling a stochastic convective parameterization with the Zhang–McFarlane scheme on precipitation simulation in the DOE E3SMv1.0 atmosphere model

Abstract. A stochastic deep convection parameterization is implemented into the US Department of Energy (DOE) Energy Exascale Earth System Model (E3SM) Atmosphere Model version 1.0 (EAMv1). This study evaluates its performance in simulating precipitation. Compared to the default model, the probability distribution function (PDF) of rainfall intensity in the new simulation is greatly improved. The well-known problem of “too much light rain and too little heavy rain” is alleviated, especially over the tropics. As a result, the contribution from different rain rates to the total precipitation amount is shifted toward heavier rain. The less frequent occurrence of convection contributes to suppressed light rain, while more intense large-scale and convective precipitation contributes to enhanced heavy total rain. The synoptic and intraseasonal variabilities of precipitation are enhanced as well to be closer to observations. The sensitivity of the rainfall intensity PDF to the model vertical resolution is examined. The relationship between precipitation and dilute convective available potential energy in the stochastic simulation agrees better with that in the Atmospheric Radiation Measurement (ARM) observations compared with the standard model simulation. The annual mean precipitation is largely unchanged with the use of the stochastic scheme except over the tropical western Pacific, where a moderate increase in precipitation represents a slight improvement. The responses of precipitation and its extremes to climate warming are similar with or without the stochastic deep convection scheme.

Improving diurnal rainfall phase over the Southern Great Plains in warm seasons by using a convective triggering design

AbstractThis study reports the improvements in the diurnal rainfall phase over the Southern Great Plains (SGP) with the Community Atmosphere Model by incorporating two convective initiation designs into the Zhang–Mcfarlane scheme. To capture the convection initiation of the nocturnal rainfall regime, the triggers emphasized the importance of suppression of convective inhibition and identification of elevated buoyancy for initiating convection. Consistent with the single‐column study of Wang et al. (2015), the model rainfall bias reduced with an increase in the heavy nocturnal rainfall and with the suppression of the spurious daytime convection in global domain experiments. We found that, with the climate model resolution of 1°, the nocturnal convection over the SGP can be well represented when the destabilization mechanisms of potentially buoyant layer at the low‐level troposphere is appropriately captured. Such destabilization is caused by the moisture convergence between the Rocky Mountains–SGP circulations and the uplifted low‐level jets from the Gulf of Mexico. Similar improvements are also found in other regions on the lee side of high terrains, such as the Tibetan Plateau.

Assessing the Resolution Adaptability of the Zhang‐McFarlane Cumulus Parameterization With Spatial and Temporal Averaging

AbstractWith increasing computational capabilities, cumulus parameterizations that are adaptable to the smaller grid spacing and temporal interval for high‐resolution climate model simulations are needed. In this study, we propose a method to improve the resolution adaptability of the Zhang‐McFarlane (ZM) scheme, by implementing spatial and temporal averaging to the CAPE tendency. This method allows for a more consistent application of the quasi‐equilibrium (QE) hypothesis at high spatial and temporal resolutions. The resolution adaptability of the original ZM scheme, the scheme with spatial averaging, and the scheme with spatiotemporal averaging at 4–32 km grid spacings are assessed using the Weather Research and Forecasting (WRF) model by comparing to cloud resolving model (CRM) simulation results coarse‐grained to these same grid spacings. We show the original ZM scheme has poor resolution adaptability, with spatiotemporally averaged subgrid convective transport and convective precipitation increasing significantly as the resolution increases. The spatial averaging method improves the resolution adaptability of the ZM scheme and better conserves total transport and total precipitation. Temporal averaging further improves the resolution adaptability of the scheme. With better constrained (although smoothed) convective transport and precipitation, both the spatial distribution and time series of total precipitation at 4 and 8 km grid spacings are improved with the averaging methods. The results could help develop resolution adaptability for other cumulus parameterizations that are based on the QE assumption.

Impacts of the triggering function of cumulus parameterization on warm-season diurnal rainfall cycles at the Atmospheric Radiation Measurement Southern Great Plains site

In this study, we investigated the impacts of the triggering function of the deep convection scheme on diurnal rainfall variation in the middle latitudes by using the single-column version of the Community Atmospheric Model (SCAM). Using the climate statistics of a long-term ensemble analysis of SCAM simulations, we quantified and validated the diurnal rainfall climatological regimes at the Atmospheric Radiation Measurement Southern Great Plains (SGP) site. The results showed that the averaged diurnal rainfall cycle simulated using the default Zhang-Mcfarlane (ZM) scheme of the SCAM peaks near noon, which is far earlier than the observed nighttime peak phase. This bias was due to the ZM scheme, which produced spurious daytime rainfall, even during days in which only light rainfall was observed. By contrast, using a weather-focused scheme, the Simplified Arakawa-Schubert (SAS) scheme, we successfully simulated the nocturnal peak of the diurnal cycle. Experiments conducted on the ZM and SAS schemes featuring different triggering functions revealed that the relaxation of launching parcels above the planetary boundary layer (PBL) and the inclusion of convective inhibition (CIN) were crucial designs for the model to capture the nocturnal rainfall events of the SGP. The inclusion of CIN reduces spurious weak convective events, and the allowance of launching parcels being above the PBL better captures convective cloud base. The results of this study highlight the modulatory effect of low-level inhomogeneity on the diurnal variation of convection over midlatitudes and the importance of the triggering function of the deep convection scheme in capturing those variations.

Development of a toy column model and its application in testing cumulus convection parameterizations

A single-column model is constructed based on parameterizations inherited from the Finite-volume/Spectral Atmospheric Model F/SAMIL and tested in simulations of tropical convective systems. Two representative convection schemes are compared in terms of their performances on precipitation types, individual physical tendencies, and temperature and moisture fields. The main difference between the two selected schemes is in their representation of entraining/detraining process. The Tiedtke scheme assumes bulk entrainment, while the Zhang-McFarlane scheme parameterizes entrainment/detrainment rates under the spectrum concept. Large-scale forcing and verification data are taken from the GATE phase III field campaign, during which abundant convective events were observed. Given the same triggering function and closure assumption, results show that entrainment/detrainment representation remains the dominant factor on the simulation of cumulus mass flux and of temperature and moisture fields. By analyzing sources and sinks of heat and moisture, this study reveals how parameterization components compensate for each other and make model results insensitive to parameterization changes in certain fields, thus suggesting the need to treat parameterizations as systems rather than individual components.

Improvements in WRF simulation skills of southeastern United States summer rainfall: physical parameterization and horizontal resolution

Realistic regional climate simulations are important in understanding the mechanisms of summer rainfall in the southeastern United States (SE US) and in making seasonal predictions. In this study, skills of SE US summer rainfall simulation at a 15-km resolution are evaluated using the weather research and forecasting (WRF) model driven by climate forecast system reanalysis data. Influences of parameterization schemes and model resolution on the rainfall are investigated. It is shown that the WRF simulations for SE US summer rainfall are most sensitive to cumulus schemes, moderately sensitive to planetary boundary layer schemes, and less sensitive to microphysics schemes. Among five WRF cumulus schemes analyzed in this study, the Zhang–McFarlane scheme outperforms the other four. Further analysis suggests that the superior performance of the Zhang–McFarlane scheme is attributable primarily to its capability of representing rainfall-triggering processes over the SE US, especially the positive relationship between convective available potential energy and rainfall. In addition, simulated rainfall using the Zhang–McFarlane scheme at the 15-km resolution is compared with that at a 3-km convection-permitting resolution without cumulus scheme to test whether the increased horizontal resolution can further improve the SE US rainfall simulation. Results indicate that the simulations at the 3-km resolution do not show obvious advantages over those at the 15-km resolution with the Zhang–McFarlane scheme. In conclusion, our study suggests that in order to obtain a satisfactory simulation of SE US summer rainfall, choosing a cumulus scheme that can realistically represent the convective rainfall triggering mechanism may be more effective than solely increasing model resolution.

Sensitivity of simulated tropical intraseasonal oscillations to cumulus schemes

The sensitivity of simulated tropical intraseasonal oscillations (ISO) to different cumulus parameterization schemes was analyzed using an atmospheric general circulation model (latest version-SAMIL2.2.3) developed at the Laboratory for Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics (LASG) at the Institute of Atmospheric Physics (IAP) of the Chinese Academy of Sciences. Results show that the basic features of tropical climatological intraseasonal oscillations (CISO) can be captured using all three cumulus schemes. The CISO simulated by the Tiedtke scheme was found to be more realistic than that of the Manabe and Zhang-McFarlane schemes. The results of simulated transient intraseasonal oscillations (TISO) indicate that although the Tiedtke and the Zhang-McFarlane schemes in the new version SAMIL2.2.3 have been adjusted according to different problems, only the latter can simulate the eastward propagation of the 27–50-day TISO mode. It may be associated with the more realistic diabatic heating profile simulated by the Zhang-McFarlane scheme. In addition, the Manabe scheme in SAMIL2.2.3 is the same as that in the prior version SAMIL2.08. However, some aspects of the physical process, such as the radiation scheme and aerosol condition, have been changed. Conversely the eastward propagation from 100°E to the west of the tropical 27–50-day TISO mode only can be simulated using the Manabe scheme of SAMIL 2.08. Consequently, not all the improvements of physical parameterization schemes work well in every respect. The coordinated developments between dynamic frame and physical processes, and among different physical processes, are important methods that may be used to improve the model.

Convection Parameterization, Tropical Pacific Double ITCZ, and Upper-Ocean Biases in the NCAR CCSM3. Part II: Coupled Feedback and the Role of Ocean Heat Transport

Abstract This study investigates the coupled atmosphere–ocean feedback and the role of ocean dynamic heat transport in the formation of double ITCZ over the tropical Pacific in the NCAR Community Climate System Model, version 3 (CCSM3) and its alleviation when a revised Zhang–McFarlane (ZM) convection scheme is used. A hierarchy of coupling strategy is employed for this purpose. A slab ocean model is coupled with the atmospheric component of the Community Atmosphere Model, version 3 (CAM3) to investigate the local feedback between the atmosphere and the ocean. It is shown that the net surface energy flux differences in the southern ITCZ region between the revised and original ZM scheme seen in the stand-alone CAM3 simulations can cool the SST by up to 1.5°C. However, the simulated SST distribution is very sensitive to the prescribed ocean heat transport required in the slab ocean model. To understand the role of ocean heat transport, the fully coupled CCSM3 model is used. The analysis of CCSM3 simulations shows that the altered ocean dynamic heat transport when the revised ZM scheme is used is largely responsible for the reduction of SST bias in the southern ITCZ region, although surface energy flux also helps to cool the SST in the first few months of the year in seasonal variation. The results, together with those from Part I, suggest that the unrealistic simulation of convection over the southern ITCZ region in the standard CCSM3 leads to the double-ITCZ bias through complex coupled interactions between atmospheric convection, surface winds, latent heat flux, cloud radiative forcing, SST, and upper-ocean circulations. The mitigation of the double-ITCZ bias using the revised ZM scheme is achieved by altering this chain of interactions.

The Role of Convective Moistening in the Madden–Julian Oscillation

Abstract This study compares two models that differ primarily in their cloud parameterizations and produce extremely different simulations of the Madden–Julian oscillation (MJO). The Community Atmosphere Model (CAM) version 3.0 from NCAR uses the Zhang–McFarlane scheme for deep convection and does not produce an MJO. The “superparameterized” version of the CAM (SP-CAM) replaces the cloud parameterizations with a two-dimensional cloud-resolving model (CRM) in each grid column and produces an extremely vigorous MJO. This analysis shows that the CAM is unable to produce high-humidity regions in the mid- to lower troposphere because of a lack of coupling between parameterized convection and environmental relative humidity. The SP-CAM produces an overly moist column due in part to excessive near-surface winds and evaporation during strong convective events. In the real tropics and the SP-CAM, convection within a high-humidity environment produces intense latent heating, which excites the large-scale circulation that is the signature of the MJO. The authors suggest that a model must realistically represent convective processes that moisten the entire tropical troposphere in order to produce a simulation of the MJO.

Convection Parameterization, Tropical Pacific Double ITCZ, and Upper-Ocean Biases in the NCAR CCSM3. Part I: Climatology and Atmospheric Feedback

Abstract The role of convection parameterization in the formation of double ITCZ and associated upper-ocean biases in the NCAR Community Climate System Model, version 3 (CCSM3) is investigated by comparing the simulations using the original and revised Zhang–McFarlane (ZM) convection schemes. Ten-year model climatologies show that the simulation with the original ZM scheme produces a typical double ITCZ bias, whereas all biases related to the spurious double ITCZ and overly strong cold tongue in precipitation, sea surface temperature (SST), wind stress, ocean thermocline, upper-ocean currents, temperature, and salinity are dramatically reduced when the revised ZM scheme is used. These results demonstrate that convection parameterization plays a critical role in the formation of double ITCZ bias in the CCSM3. To understand the physical mechanisms through which the modifications of the convection scheme in the atmospheric model alleviate the double ITCZ bias in the CCSM3, the authors investigate the impacts of convection schemes on the atmospheric forcing and feedback in the uncoupled Community Atmospheric Model, version 3 (CAM3). It is shown that the CAM3 simulation with the original ZM scheme also produces a signature of double ITCZ bias in precipitation, whereas the simulation with the revised ZM scheme does not. Diagnostic analyses have identified three factors on the atmospheric side (i.e., the sensitivity of convection to SST, the convection–shortwave flux–SST feedback, and the convection–wind–evaporation–SST feedback) that may contribute to the differences in the coupled simulations.

The MJO in an AGCM with three different cumulus parameterization schemes

The Madden–Julian oscillation (MJO) is simulated using an AGCM with three different cumulus parameterization schemes: a moist convective adjustment (MCA) scheme, the Zhang–McFarlane (ZM) mass-flux scheme, and the Tiedtke scheme. Results show that the simulated MJO is highly dependent on the cumulus parameterization used. Among the three cumulus parameterizations, only the MCA scheme produces MJO features similar to observations, including the reasonable spatial distribution, intraseasonal time scales and eastward propagation. Meanwhile, the amplitude is too large and the eastward propagation speed too fast than observations and the relationship between precipitation and low-level wind anomaly is unrealistic with enhanced convection occurring within easterly anomalies instead of westerly anomalies as in observations. The over-dependence of precipitation on boundary convergence produced by the MCA scheme is presumably responsible for this unrealistic phase relation in the simulation. The other two schemes produce very poor simulations of the MJO: spectral power of westward propagation is larger than that of eastward propagation in zonal wind and precipitation, indicating a westward propagation of the intraseasonal variability. The mean state and vertical profile of diabatic heating are perhaps responsible for the differences in these simulations. The MCA scheme produces relatively realistic climate background. When either ZM or Tiedtke scheme is used, the observed extension of westerly winds from the western Pacific to the dateline is missing and precipitation over the equatorial region and SPCZ is dramatically underestimated. In addition, diabatic heating produced by both ZM and Tiedtke schemes are very weak and nearly uniform with height. The heating profile produced by the MCA scheme has a middle-heavy structure with much larger magnitude than those produced by the other two schemes. In addition, a very unrealistic boundary layer heating maximum produced by the MCA scheme induces too strong surface convergence, which perhaps contributes to the too strong intraseasonal variability in the simulation.

Improved Forecasts of the Diurnal Cycle in the Tropics Using Multiple Global Models. Part I: Precipitation

Abstract The Tropical Rainfall Measuring Mission (TRMM) satellite supplemented with the Defense Meteorological Satellites Program (DMSP) microwave dataset provides accurate rain-rate estimates. Furthermore, infrared radiances from the geostationary satellites provide the possibility for mapping the diurnal change of tropical rainfall. Modeling of the phase and amplitude of the tropical rainfall is the theme of this paper. The present study utilizes a suite of global multimodels that are identical in all respects except for their cumulus parameterization algorithms. Six different cumulus parameterizations are tested in this study. These include the Florida State University (FSU) Modified Kuo Scheme (KUO), Goddard Space Flight Center (GSFC) Relaxed Arakawa–Schubert Scheme (RAS1), Naval Research Laboratory–Navy Operational Global Atmospheric Prediction System (NRL–NOGAPS) Relaxed Arakawa–Schubert Scheme (RAS2), NCEP Simplified Arakawa–Schubert Scheme (SAS), NCAR Zhang–McFarlane Scheme (ZM), and NRL–NOGAPS Emanuel Scheme (ECS). The authors carried out nearly 600 experiments with these six versions of the T170 Florida State University global spectral model. These are 5-day NWP experiments where the diurnal change datasets were archived at 3-hourly intervals. This study includes the estimation of skills of the phase and amplitudes of the diurnal rain using these member models, their ensemble mean, a multimodel superensemble, and those from a single unified model. Test results are presented for the global tropics and for some specific regions where the member models show difficulty in predicting the diurnal change of rainfall. The main contribution is the considerable improvement of the modeling of diurnal rain by deploying a multimodel superensemble and by constructing a single unified model. The authors also present a comparison of these findings on the modeling of diurnal rain from another suite of multimodels that utilized different versions of cloud radiation algorithms (instead of different cumulus parameterization schemes) toward defining the suite of multimodels. The principal result is that the superensemble does provide a future forecast for the total daily rain and for the diurnal change of rain through day 5 that is superior to forecasts provided by the best model. The training of the superensemble with good observed estimates of rain, such as those from TRMM, is necessary for such forecasts.

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.