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.     

Isotopic composition and origin of polar precipitation in present and glacial climate simulations

The Hamburg atmospheric general circulation model (AGCM) ECHAM‐4 is used to identify the main source regions of precipitation falling on Greenland and Antarctica. Both water isotopes H₂¹⁸O and HDO are explicitly built into the water cycle of the AGCM, and in addition the capability to trace water from different source regions was added to the model. Present and LGM climate simulations show that water from the most important source regions has an isotopic signature similar to the mean isotope values of the total precipitation amount. But water from other source regions (with very different isotopic signatures) contributes an additional, non‐negligible part of the total precipitation amount on both Greenland and Antarctica. Analyses of the temperature‐isotope‐relations for both polar regions reveal a solely bias of the glacial isotope signal on Greenland, which is caused by a strong change in the seasonal deposition of precipitation originating from nearby polar seas and the northern Atlantic. Although the performed simulations under LGM boundary conditions show a decrease of the δ ¹⁸O values in precipitation in agreement with ice core measurements, the AGCM fails to reproduce the observed simultaneous decrease of the deuterium excess signal.     

Coupled climate modelling of ocean circulation changes during ice age inception

Freshening of high latitude surface waters can change the large-scale oceanic transport of heat and salt. Consequently, atmospheric and sea ice perturbations over the deep water production sites excite a large-scale response establishing an oceanic "teleconnection" with time scales of years to centuries. To study these feedbacks, a coupled atmosphere-ocean-sea ice model consisting of a two dimensional atmospheric energy and moisture balance model (EMBM) coupled to a thermodynamic sea ice model and an ocean general circulation model is utilised. The coupled model reproduces many aspects of the present oceanic circulation. We also investigate the climate impact of changes in fresh water balance during an ice age initiation. In this experiment part of the precipitation over continents is stored within continental ice sheets. During the buildup of ice sheets the oceanic stratification in the North Atlantic is weakened by a reduced continental run-off leading to an enhanced thermohaline circulation. Under these conditions salinity is redistributed such that deep water is more saline than under present conditions. Once the ice sheets built up, we simulate an ice age climate without net fresh water storage on the continents. In this case the coupled model reproduces the shallow and weak overturning cell, an ice edge advance insulating the upper ocean, and many other aspects of the glacial circulation.     

Frequency distributions of transient regional climate change from perturbed physics ensembles of general circulation model simulations

Large ensembles of coupled atmosphere–ocean general circulation model (AOGCM) simulations are required to explore modelling uncertainty and make probabilistic predictions of future transient climate change at regional scales. These are not yet computationally feasible so we have developed a technique to emulate the response of such an ensemble by scaling equilibrium patterns of climate change derived from much cheaper “slab” model ensembles in which the atmospheric component of an AOGCM is coupled to a mixed-layer ocean. Climate feedback parameters are diagnosed for each member of a slab model ensemble and used to drive an energy balance model (EBM) to predict the time-dependent response of global surface temperature expected for different combinations of uncertain AOGCM parameters affecting atmospheric, land and sea-ice processes. The EBM projections are then used to scale normalised patterns of change derived for each slab member, and hence emulate the response of the relevant atmospheric model version when coupled to a dynamic ocean, in response to a 1% per annum increase in CO₂. The emulated responses are validated by comparison with predictions from a 17 member ensemble of AOGCM simulations, constructed from variants of HadCM3 using the same parameter combinations as 17 members of the slab model ensemble. Cross-validation permits estimation of the spatial and temporal dependence of emulation error, and also allows estimation of a correction field to correct discrepancies between the scaled equilibrium patterns and the transient response, reducing the emulation error. Emulated transient responses and their associated errors are obtained from the slab ensemble for 129 pseudo-HadCM3 versions containing multiple atmospheric parameter perturbations. These are combined to produce regional frequency distributions for the transient response of annual surface temperature change and boreal winter precipitation change. The technique can be extended to any surface climate variable demonstrating a scaleable, approximately linear response to forcing.     

A global ocean biogeochemistry general circulation model and its simulations

An ocean biogeochemistry model was developed and incorporated into a global ocean general circulation model (LICOM) to form an ocean biogeochemistry general circulation model (OBGCM). The model was used to study the natural carbon cycle and the uptake and storage of anthropogenic CO2 in the ocean. A global export production of 12.5 Pg C yr-1 was obtained. The model estimated that in the pre-industrial era the global equatorial region within 15o of the equator released 0.97 Pg C yr-1 to the atmosphere, which was balanced by the gain of CO2 in other regions. The post-industrial air-sea CO2 flux indicated the oceanic uptake of CO2 emitted by human activities. An increase of 20－50 mol kg-1 for surface dissolved inorganic carbon (DIC) concentrations in the 1990s relative to pre-industrial times was obtained in the simulation, which was consistent with data-based estimates. The model generated a total anthropogenic carbon inventory of 105 Pg C as of 1994, which was within the range of estimates by other researchers. Various transports of both natural and anthropogenic DIC as well as labile dissolved organic carbon (LDOC) were estimated from the simulation. It was realized that the Southern Ocean and the high-latitude region of the North Pacific are important export regions where accumulative air-sea CO2 fluxes are larger than the DIC inventory, whereas the subtropical regions are acceptance regions. The interhemispheric transport of total natural carbon (DIC+LDOC) was found to be northward (0.11 Pg C yr-1), which was just balanced by the gain of carbon from the atmosphere in the Southern Hemisphere.     

Characteristics of sub-daily precipitation extremes in observed data and regional climate model simulations

Theoretical and Applied Climatology - The study compares characteristics of observed sub-daily precipitation extremes in the Czech Republic with those simulated by Hadley Centre Regional Model...     

Ocean-atmosphere interactions and climate drift in a coupled general circulation model

 We have analysed numerical simulations performed with a global 3D coupled atmosphere-ocean model to focus on the role of atmospheric processes leading to sea surface temperature (SST) drift in the tropics. Negative SST errors occur coherently in space and time with large positive errors in latent heat and momentum fluxes at the tropical air-sea interface, as diagnosed from forced SST simulations. The warm pool in the western Pacific disappears after a few years of simulation. Strong SST gradients enforce regions of high precipitation that are thin and stationary north of the equator. We detail the implications for the ocean-atmosphere system of such upheaval in the deep convection location. A sensitivity experiment to empirically formulate air-sea drag coefficient shows that the rapid warm pool erosion is not sensitive to changes in the formulation of the surface drag coefficient over the oceans because the corresponding changes in turbulent heat fluxes and LW cooling approximately cancel one another. In the eastern Pacific, the improvement in SST is striking and caused by feedbacks between SST, surface turbulent fluxes and boundary layer cloud fraction, which decreases as SST warms. Received: 8 December 1998 / Accepted: 6 January 2000     

Convective and stratiform precipitation characteristics in an ensemble of regional climate model simulations

We apply a recently proposed algorithm for disaggregating observed precipitation data into predominantly convective and stratiform, and evaluate biases in characteristics of parameterized convective (subgrid) and stratiform (large-scale) precipitation in an ensemble of 11 RCM simulations for recent climate in Central Europe. All RCMs have a resolution of 25 km and are driven by the ERA-40 reanalysis. We focus on mean annual cycle, proportion of convective precipitation, dependence on altitude, and extremes. The results show that characteristics of total precipitation are often better simulated than are those of convective and stratiform precipitation evaluated separately. While annual cycles of convective and stratiform precipitation are reproduced reasonably well in most RCMs, some of them consistently and substantially overestimate or underestimate the proportion of convective precipitation throughout the year. Intensity of convective precipitation is underestimated in all RCMs. Dependence on altitude is also simulated better for stratiform and total precipitation than for convective precipitation, for which several RCMs produce unrealistic slopes. Extremes are underestimated for convective precipitation while they tend to be slightly overestimated for stratiform precipitation, thus resulting in a relatively good reproduction of extremes in total precipitation amounts. The results suggest that the examined ensemble of RCMs suffers from substantial deficiencies in reproducing precipitation processes and support previous findings that climate models’ errors in precipitation characteristics are mainly related to deficiencies in the representation of convection.     

Tropical precipitation regimes and mechanisms of regime transitions: contrasting two aquaplanet general circulation models

The atmospheric general circulation models ARPEGE-climate and LMDz are used in an aquaplanet configuration to study the response of a zonally symmetric atmosphere to a range of sea surface temperature (SST) forcing. We impose zonally-symmetric SST distributions that are also symmetric about the equator, with varying off-equatorial SST gradients. In both models, we obtain the characteristic inter-tropical convergence zone (ITCZ) splitting that separates two regimes of equilibrium (in terms of precipitations): one with one ITCZ over the equator for large SST gradients in the tropics, and one with a double ITCZ for small tropical SST gradients. Transition between these regimes is mainly driven by changes in the low-level convergence that are forced by the SST gradients. Model-dependent, dry and moist feedbacks intervene to reinforce or weaken the effect of the SST forcing. In ARPEGE, dry advective processes reinforce the SST forcing, while a competition between sensible heat flux and convective cooling provides a complex feedback on the SST forcing in the LMDz. It is suggested that these feedbacks influence the location of the transition in the parameter range.     

Decadal climate variability simulated in a coupled general circulation model

A 1000 year integration of the CSIRO coupled ocean-atmosphere general circulation model is used to study low frequency (decadal to centennial) climate variability in precipitation and temperature. The model is shown to exhibit sizeable decadal variability for these fields, generally accounting for approximately 20 to 40% of the variability (greater than one year) in precipitation and up to 80% for temperature. An empirical orthogonal function (EOF) analysis is applied to the model output to show some of the major statistical modes of low frequency variability. The first EOF spatial pattern looks very much like that of the interannual ENSO pattern. It bears considerable resemblance to observational estimates and is centred in the Pacific extending into both hemispheres. It modulates both precipitation and temperature globally. The EOF has a time evolution that appears to be more than just red noise. Finally, the link between SST in the Pacific with Australian rainfall variability seen in observations is also evident in the model. Received: 29 August 1998 / Accepted: 31 July 1999     

Downscaling and projection of precipitation from general circulation model predictors in an equatorial climate region by the automated regression-based statistical method

The authors have applied an automated regression-based statistical method, namely, the automated statistical downscaling (ASD) model, to downscale and project the precipitation climatology in an equatorial climate region (Peninsular Malaysia). Five precipitation indices are, principally, downscaled and projected: mean monthly values of precipitation (Mean), standard deviation (STD), 90th percentile of rain day amount, percentage of wet days (Wet-day), and maximum number of consecutive dry days (CDD). The predictors, National Centers for Environmental Prediction (NCEP) products, are taken from the daily series reanalysis data, while the global climate model (GCM) outputs are from the Hadley Centre Coupled Model, version 3 (HadCM3) in A2/B2 emission scenarios and Third-Generation Coupled Global Climate Model (CGCM3) in A2 emission scenario. Meanwhile, the predictand data are taken from the arithmetically averaged rain gauge information and used as a baseline data for the evaluation. The results reveal, from the calibration and validation periods spanning a period of 40 years (1961–2000), the ASD model is capable to downscale the precipitation with reasonable accuracy. Overall, during the validation period, the model simulations with the NCEP predictors produce mean monthly precipitation of 6.18–6.20 mm/day (root mean squared error 0.78 and 0.82 mm/day), interpolated, respectively, on HadCM3 and CGCM3 grids, in contrast to 6.00 mm/day as observation. Nevertheless, the model suffers to perform reasonably well at the time of extreme precipitation and summer time, more specifically to generate the CDD and STD indices. The future projections of precipitation (2011–2099) exhibit that there would be an increase in the precipitation amount and frequency in most of the months. Taking the 1961–2000 timeline as the base period, overall, the annual mean precipitation would indicate a surplus projection by nearly 14~18 % under both GCM output cases (HadCM3 A2/B2 scenarios and CGCM3 A2 scenario). According to the model simulation, the September–November periods might be the more significant months projecting the increment of the precipitation amount around over 50 %, while the precipitation deficit would be seen in March–May periods.     

Non-linear climate feedback analysis in an atmospheric general circulation model

 A method is described for evaluating the ‘partial derivatives’ of globally averaged top-of-atmosphere (TOA) radiation changes with respect to basic climate model physical parameters. This method is used to analyse feedbacks in the Australian Bureau of Meteorology Research Centre general circulation model. The parameters considered are surface temperature, water vapour, lapse rate and cloud cover. The climate forcing which produces the changes is a globally uniform sea surface temperature (SST) perturbation. The first and second order differentials of model parameters with respect to the forcing (i.e. SST changes) are estimated from quadratic least square fitting. Except for total cloud cover, variables are found to be strong functions of global SST. Strongly non-linear variations of lapse rate and high cloud amount and height appear to relate to the non-linear response in penetrative convection. Globally averaged TOA radiation differentials with respect to model parameters are also evaluated. With the exception of total cloud contributions, a high correlation is generally found to exist, on the global mean level, between TOA radiation and the respective parameter perturbations. The largest non-linear terms contributing to radiative changes are those due to lapse rate and high cloud. The contributions of linear and non-linear terms to the overall radiative response from a 4 K SST perturbation are assessed. Significant non-linear responses are found to be associated with lapse rate, water vapour and cloud changes. Although the exact magnitude of these responses is likely to be a function of the particular model as well as the imposed SST perturbation pattern, the present experiments flag these as processes which cannot properly be understood from linear theory in the evaluation of climate change sensitivity. Received: 16 January 1997/Accepted: 9 May 1997     

Walker circulation controls ENSO atmospheric feedbacks in uncoupled and coupled climate model simulations

Many climate models strongly underestimate the two most important atmospheric feedbacks operating in El Niño/Southern Oscillation (ENSO), the positive (amplifying) zonal surface wind feedback and negative (damping) surface-heat flux feedback (hereafter ENSO atmospheric feedbacks, EAF). This hampers a realistic representation of ENSO dynamics in these models. Here we show that the atmospheric components of climate models participating in the 5th phase of the Coupled Model Intercomparison Project (CMIP5) when forced by observed sea surface temperatures (SST), already underestimate EAF on average by 23%, but less than their coupled counterparts (on average by 54%). There is a pronounced tendency of atmosphere models to simulate stronger EAF, when they exhibit a stronger mean deep convection and enhanced cloud cover over the western equatorial Pacific (WEP), indicative of a stronger rising branch of the Pacific Walker Circulation (PWC). Further, differences in the mean deep convection over the WEP between the coupled and uncoupled models explain a large part of the differences in EAF, with the deep convection in the coupled models strongly depending on the equatorial Pacific SST bias. Experiments with a single atmosphere model support the relation between the equatorial Pacific atmospheric mean state, the SST bias and the EAF. An implemented cold SST bias in the observed SST forcing weakens deep convection and reduces cloud cover in the rising branch of the PWC, causing weaker EAF. A warm SST bias has the opposite effect. Our results elucidate how biases in the mean state of the PWC and equatorial SST hamper a realistic simulation of the EAF.     

Stratospheric climate and variability from a general circulation model and observations

The climate and natural variability of the large-scale stratospheric circulation simulated by a newly developed general circulation model are evaluated against available global observations. The simulation consisted of a 30-year annual cycle integration performed with a comprehensive model of the troposphere and stratosphere. The observations consisted of a 15-year dataset from global operational analyses of the troposphere and stratosphere. The model evaluation concentrates on the simulation of the evolution of the extratropical stratospheric circulation in both hemispheres. The December–February climatology of the observed zonal mean winter circulation is found to be reasonably well captured by the model, although in the Northern Hemisphere upper stratosphere the simulated westerly winds are systematically stronger and a cold bias is apparent in the polar stratosphere. This Northern Hemisphere stratospheric cold bias virtually disappears during spring (March–May), consistent with a realistic simulation of the spring weakening of the mean westerly winds in the model. A considerable amount of monthly interannual variability is also found in the simulation in the Northern Hemisphere in late winter and early spring. The simulated interannual variability is predominantly caused by polar warmings of the stratosphere, in agreement with observations. The breakdown of the Northern Hemisphere stratospheric polar vortex appears therefore to occur in a realistic way in the model. However, in early winter the model severely underestimates the interannual variability, especially in the upper troposphere. The Southern Hemisphere winter (June–August) zonal mean temperature is systematically colder in the model, and the simulated winds are somewhat too strong in the upper stratosphere. Contrary to the results for the Northern Hemisphere spring, this model cold bias worsens during the Southern Hemisphere spring (September–November). Significant discrepancies between the model results and the observations are therefore found during the breakdown of the Southern Hemisphere polar vortex. For instance, the simulated Southern Hemisphere stratosphere westerly jet continuously decreases in intensity more or less in situ from June to November, while the observed stratospheric jet moves downward and poleward.This paper was presented at the Third International Conference on Modelling of Global Climate Change and Variability, held in Hamburg 4–8 Sept. 1995 under the auspice of the Max Planck Institute for Meteorology, Hamburg. Editor for these papers is L. Dümenil.     

The sensitivity of the tropical circulation and Maritime Continent precipitation to climate model resolution

The dependence of the annual mean tropical precipitation on horizontal resolution is investigated in the atmospheric version of the Hadley Centre General Environment Model. Reducing the grid spacing from about 350 km to about 110 km improves the precipitation distribution in most of the tropics. In particular, characteristic dry biases over South and Southeast Asia including the Maritime Continent as well as wet biases over the western tropical oceans are reduced. The annual-mean precipitation bias is reduced by about one third over the Maritime Continent and the neighbouring ocean basins associated with it via the Walker circulation. Sensitivity experiments show that much of the improvement with resolution in the Maritime Continent region is due to the specification of better resolved surface boundary conditions (land fraction, soil and vegetation parameters) at the higher resolution. It is shown that in particular the formulation of the coastal tiling scheme may cause resolution sensitivity of the mean simulated climate. The improvement in the tropical mean precipitation in this region is not primarily associated with the better representation of orography at the higher resolution, nor with changes in the eddy transport of moisture. Sizeable sensitivity to changes in the surface fields may be one of the reasons for the large variation of the mean tropical precipitation distribution seen across climate models.     

Entropy production and coarse graining of the climate fields in a general circulation model

We extend the analysis of the thermodynamics of the climate system by investigating the role played by processes taking place at various spatial and temporal scales through a procedure of coarse graining. We show that the coarser is the graining of the climatic fields, the lower is the resulting estimate of the material entropy production. In other terms, all the spatial and temporal scales of variability of the thermodynamic fields provide a positive contribution to the material entropy production. This may be interpreted also as that, at all scales, the temperature fields and the heating fields resulting from the convergence of turbulent fluxes have a negative correlation, while the opposite holds between the temperature fields and the radiative heating fields. Moreover, we obtain that the latter correlations are stronger, which confirms that radiation acts as primary driver for the climatic processes, while the material fluxes dampen the resulting fluctuations through dissipative processes. We also show, using specific coarse-graining procedures, how one can separate the various contributions to the material entropy production coming from the dissipation of kinetic energy, the vertical sensible and latent heat fluxes, and the large scale horizontal fluxes, without resorting to the full three-dimensional time dependent fields. We find that most of the entropy production is associated to irreversible exchanges occurring along the vertical direction, and that neglecting the horizontal and time variability of the fields has a relatively small impact on the estimate of the material entropy production. The approach presented here seems promising for testing climate models, for assessing the impact of changing their parametrizations and their resolution, as well as for investigating the atmosphere of exoplanets, because it allows for evaluating the error in the estimate of their thermodynamical properties due to the lack of high-resolution data. The findings on the impact of coarse graining on the thermodynamic fields on the estimate of the material entropy production deserve to be explored in a more general context, because they provide a way for understanding the relationship between forced fluctuations and dissipative processes in continuum systems.     

Resolution effects on regional climate model simulations of seasonal precipitation over Europe

We analyze a set of nine regional climate model simulations for the period 1961–2000 performed at 25 and 50 km horizontal grid spacing over a European domain in order to determine the effects of horizontal resolution on the simulation of precipitation. All of the models represent the seasonal mean spatial patterns and amount of precipitation fairly well. Most models exhibit a tendency to over-predict precipitation, resulting in a domain-average total bias for the ensemble mean of about 20% in winter (DJF) and less than 10% in summer (JJA) at both resolutions, although this bias could be artificially enhanced by the lack of a gauge correction in the observations. A majority of the models show increased precipitation at 25 km relative to 50 km over the oceans and inland seas in DJF, JJA, and ANN (annual average), although the response is strongest during JJA. The ratio of convective precipitation to total precipitation decreases over land for most models at 25 km. In addition, there is an increase in interannual variability in many of the models at 25 km grid spacing. Comparison with gridded observations indicates that a majority of models show improved skill in simulating both the spatial pattern and temporal evolution of precipitation at 25 km compared to 50 km during the summer months, but not in winter or on an annual mean basis. Model skill at higher resolution in simulating the spatial and temporal character of seasonal precipitation is found especially for Great Britain. This geographic dependence of the increased skill suggests that observed data of sufficient density are necessary to capture fine-scale climate signals. As climate models increase their horizontal resolution, it is thus a key priority to produce high quality fine scale observations for model evaluation.