Impact of the LMDZ atmospheric grid configuration on the climate and sensitivity of the IPSL-CM5A coupled model

The IPSL-CM5A climate model was used to perform a large number of control, historical and climate change simulations in the frame of CMIP5. The refined horizontal and vertical grid of the atmospheric component, LMDZ, constitutes a major difference compared to the previous IPSL-CM4 version used for CMIP3. From imposed-SST (Sea Surface Temperature) and coupled numerical experiments, we systematically analyze the impact of the horizontal and vertical grid resolution on the simulated climate. The refinement of the horizontal grid results in a systematic reduction of major biases in the mean tropospheric structures and SST. The mid-latitude jets, located too close to the equator with the coarsest grids, move poleward. This robust feature, is accompanied by a drying at mid-latitudes and a reduction of cold biases in mid-latitudes relative to the equator. The model was also extended to the stratosphere by increasing the number of layers on the vertical from 19 to 39 (15 in the stratosphere) and adding relevant parameterizations. The 39-layer version captures the dominant modes of the stratospheric variability and exhibits stratospheric sudden warmings. Changing either the vertical or horizontal resolution modifies the global energy balance in imposed-SST simulations by typically several W/m² which translates in the coupled atmosphere-ocean simulations into a different global-mean SST. The sensitivity is of about 1.2 K per 1 W/m² when varying the horizontal grid. A re-tuning of model parameters was thus required to restore this energy balance in the imposed-SST simulations and reduce the biases in the simulated mean surface temperature and, to some extent, latitudinal SST variations in the coupled experiments for the modern climate. The tuning hardly compensates, however, for robust biases of the coupled model. Despite the wide range of grid configurations explored and their significant impact on the present-day climate, the climate sensitivity remains essentially unchanged.     

The LMDZ4 general circulation model: climate performance and sensitivity to parametrized physics with emphasis on tropical convection

The LMDZ4 general circulation model is the atmospheric component of the IPSL–CM4 coupled model which has been used to perform climate change simulations for the 4th IPCC assessment report. The main aspects of the model climatology (forced by observed sea surface temperature) are documented here, as well as the major improvements with respect to the previous versions, which mainly come form the parametrization of tropical convection. A methodology is proposed to help analyse the sensitivity of the tropical Hadley–Walker circulation to the parametrization of cumulus convection and clouds. The tropical circulation is characterized using scalar potentials associated with the horizontal wind and horizontal transport of geopotential (the Laplacian of which is proportional to the total vertical momentum in the atmospheric column). The effect of parametrized physics is analysed in a regime sorted framework using the vertical velocity at 500 hPa as a proxy for large scale vertical motion. Compared to Tiedtke’s convection scheme, used in previous versions, the Emanuel’s scheme improves the representation of the Hadley–Walker circulation, with a relatively stronger and deeper large scale vertical ascent over tropical continents, and suppresses the marked patterns of concentrated rainfall over oceans. Thanks to the regime sorted analyses, these differences are attributed to intrinsic differences in the vertical distribution of convective heating, and to the lack of self-inhibition by precipitating downdraughts in Tiedtke’s parametrization. Both the convection and cloud schemes are shown to control the relative importance of large scale convection over land and ocean, an important point for the behaviour of the coupled model.     

Control of deep convection by sub-cloud lifting processes: the ALP closure in the LMDZ5B general circulation model

Recently, a new conceptual framework for deep convection scheme triggering and closure has been developed and implemented in the LMDZ5B general circulation model, based on the idea that deep convection is controlled by sub-cloud lifting processes. Such processes include boundary-layer thermals and evaporatively-driven cold pools (wakes), which provide an available lifting energy that is compared to the convective inhibition to trigger deep convection, and an available lifting power (ALP) at cloud base, which is used to compute the convective mass flux assuming the updraft vertical velocity at the level of free convection. While the ALP closure was shown to delay the local hour of maximum precipitation over land in better agreement with observations, it results in an underestimation of the convection intensity over the tropical ocean both in the 1D and 3D configurations of the model. The specification of the updraft vertical velocity at the level of free convection appears to be a key aspect of the closure formulation, as it is weaker over tropical ocean than over land and weaker in moist mid-latitudes than semi-arid regions. We propose a formulation making this velocity increase with the level of free convection, so that the ALP closure is adapted to various environments. Cloud-resolving model simulations of observed oceanic and continental case studies are used to evaluate the representation of lifting processes and test the assumptions at the basis of the ALP closure formulation. Results favor closures based on the lifting power of sub-grid sub-cloud processes rather than those involving quasi-equilibrium with the large-scale environment. The new version of the model including boundary-layer thermals and cold pools coupled together with the deep convection scheme via the ALP closure significantly improves the representation of various observed case studies in 1D mode. It also substantially modifies precipitation patterns in the full 3D version of the model, including seasonal means, diurnal cycle and intraseasonal variability.     

Mid-Holocene and last glacial maximum climate simulations with the IPSL model: part II: model-data comparisons

The climates of the mid-Holocene (MH, 6,000 years ago) and the Last Glacial Maximum (LGM, 21,000 years ago) have been extensively documented and as such, have become targets for the evaluation of climate models for climate contexts very different from the present. In Part 1 of the present work, we have studied the MH and LGM simulations performed with the last two versions of the IPSL model: IPSL_CM4, run for the PMIP2/CMIP3 (Coupled Model Intercomparion Project) projects and IPSL_CM5A, run for the most recent PMIP3/CMIP5 projets. We have shown that not only are these models different in their simulations of the PI climate, but also in their simulations of the climatic anomalies for the MH and LGM. In the Part 2 of this paper, we first examine whether palaeo-data can help discriminate between the model performances. This is indeed the case for the African monsoon for the MH or for North America south of the Laurentide ice sheet, the South Atlantic or the southern Indian ocean for the LGM. For the LGM, off-line vegetation modelling appears to offer good opportunities to distinguish climate model results because glacial vegetation proves to be very sensitive to even small differences in LGM climate. For other cases such as the LGM North Atlantic or the LGM equatorial Pacific, the large uncertainty on the SST reconstructions, prevents model discrimination. We have examined the use of other proxy-data for model evaluation, which has become possible with the inclusion of the biogeochemistry morel PISCES in the IPSL_CM5A model. We show a broad agreement of the LGM–PI export production changes with reconstructions. These changes are related to the mixed layer depth in most regions and to sea-ice variations in the high latitudes. We have also modelled foraminifer abundances with the FORAMCLIM model and shown that the changes in foraminifer abundance in the equatorial Pacific are mainly forced by changes in SSTs, hence confirming the SST-foraminifer abundance relationship. Yet, this is not the case in all regions in the North Atlantic, where food availability can have a strong impact of foraminifer abundances. Further work will be needed to exhaustively examine the role of factors other than climate in piloting changes in palaeo-indicators.     

Key features of the IPSL ocean atmosphere model and its sensitivity to atmospheric resolution

This paper presents the major characteristics of the Institut Pierre Simon Laplace (IPSL) coupled ocean–atmosphere general circulation model. The model components and the coupling methodology are described, as well as the main characteristics of the climatology and interannual variability. The model results of the standard version used for IPCC climate projections, and for intercomparison projects like the Paleoclimate Modeling Intercomparison Project (PMIP 2) are compared to those with a higher resolution in the atmosphere. A focus on the North Atlantic and on the tropics is used to address the impact of the atmosphere resolution on processes and feedbacks. In the North Atlantic, the resolution change leads to an improved representation of the storm-tracks and the North Atlantic oscillation. The better representation of the wind structure increases the northward salt transports, the deep-water formation and the Atlantic meridional overturning circulation. In the tropics, the ocean–atmosphere dynamical coupling, or Bjerknes feedback, improves with the resolution. The amplitude of ENSO (El Niño-Southern oscillation) consequently increases, as the damping processes are left unchanged.     

Development of an ensemble data assimilation system with LMDZ5 AGCM for regional reanalysis

Climate Dynamics - The present study describes a newly developed LMDZ–DART data assimilation system based on the Ensemble Kalman Filter (EnKF) with a stretched grid atmospheric model and...     

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.     

Resolution of the atmospheric model matters for the Northern Hemisphere Mid-Holocene climate

This study evaluates the dependence of simulated surface air temperatures on model resolution and orography for the mid-Holocene. Sensitivity experiments with the atmospheric general circulation model ECHAM5 are performed with low (∼3.75°, 19 vertical levels) and high (∼1.1°, 31 vertical levels) resolution. Results are compared to the respective preindustrial runs. It is found that the large-scale temperature anomalies for the mid-Holocene (compared to preindustrial) are significantly different in the low- and high-resolution versions. For boreal winter, differences are mainly related to circulation changes caused by the response to thermal forcing in conjunction with orographic resolution. For summer, shortwave cloud radiative forcing emerges as an important factor. The anomaly differences (low minus high resolution version) in the Northern Hemisphere are regionally as large as the anomalous mid-Holocene temperature signals. Furthermore, they depend on the applied surface boundary conditions. We conclude that the resolution matters for the Northern Hemisphere response in mid-Holocene simulations, which should be taken into account in model-model and data-model comparisons.     

Performance of two cumulus convection parameterizations for two heavy precipitation events in the Western Mediterranean

Summary A set of mesoscale numerical simulations using the Emanuel and Kain-Fritsch deep convection schemes has been performed in order to determine the sensitivity of the forecast-especially, the rainfall-to the scheme used. The study is carried out for two cases of heavy precipitation in the coastal zone of the Western Mediterranean, where the topographic forcing is of primary influence. The first one, characterized by an almost stationary synoptic situation, is dominated by warm, moist advection at low levels; the second one, of frontal type, presents a much stronger dynamic forcing at upper levels. Although the comparison attempt is conditioned by the limited number of considered cases, the numerical results provide at least some preliminary conclusions. The inclusion of a convective scheme improves the forecast precipitation, through two actions: directly, producing more realistic rainfall patterns in areas of convection; indirectly, avoiding excessive precipitation in areas with orographic or dynamical upward forcing by drying and stabilizing the atmosphere upstream. In particular, the Kain-Fritsch scheme seems to be more sensitive to the orographic forcing, in agreement with observations.With 21 Figures     

Validation for a tropical belt version of WRF: sensitivity tests on radiation and cumulus convection parameterizations

对基本气候态和降水日变化的分析是检验模式模拟性能、理解模式误差来源的重要手段。为了评估出对热带气候模拟效果较好的物理参数化方案组合,本文应用WRF带状区域模式,主要比较了四种积云对流参数化方案:NewTiedtke、Kain-Fritsch、newSAS、Tiedtke,和两种辐射参数化方案:RRTMG和CAM,对热带带状区域的气候模拟结果。研究表明:使用NewTiedtke积云对流参数化方案和RRTMG辐射方案的试验,表现出对气温、降水及降水日变化等综合性最好的模拟性能;NewTiedtke积云对流参数化方案能模拟出较好的降水空间分布和降水日变化位相分布特征;与RRTMG辐射方案相比,CAM辐射方案会使温度模拟偏低,特别是陆地上更明显,这种陆地上的冷偏差可能主要来源于Tmin的模拟偏冷。     

The Beijing Climate Center atmospheric general circulation model: description and its performance for the present-day climate

The Beijing Climate Center atmospheric general circulation model version 2.0.1 (BCC_AGCM2.0.1) is described and its performance in simulating the present-day climate is assessed. BCC_AGCM2.0.1 originates from the community atmospheric model version 3 (CAM3) developed by the National Center for Atmospheric Research (NCAR). The dynamics in BCC_AGCM2.0.1 is, however, substantially different from the Eulerian spectral formulation of the dynamical equations in CAM3, and several new physical parameterizations have replaced the corresponding original ones. The major modification of the model physics in BCC_AGCM2.0.1 includes a new convection scheme, a dry adiabatic adjustment scheme in which potential temperature is conserved, a modified scheme to calculate the sensible heat and moisture fluxes over the open ocean which takes into account the effect of ocean waves on the latent and sensible heat fluxes, and an empirical equation to compute the snow cover fraction. Specially, the new convection scheme in BCC_AGCM2.0.1, which is generated from the Zhang and McFarlane’s scheme but modified, is tested to have significant improvement in tropical maximum but also the subtropical minimum precipitation, and the modified scheme for turbulent fluxes are validated using EPIC2001 in situ observations and show a large improvement than its original scheme in CAM3. BCC_AGCM2.0.1 is forced by observed monthly varying sea surface temperatures and sea ice concentrations during 1949–2000. The model climatology is compiled for the period 1971–2000 and compared with the ERA-40 reanalysis products. The model performance is evaluated in terms of energy budgets, precipitation, sea level pressure, air temperature, geopotential height, and atmospheric circulation, as well as their seasonal variations. Results show that BCC_AGCM2.0.1 reproduces fairly well the present-day climate. The combined effect of the new dynamical core and the updated physical parameterizations in BCC_AGCM2.0.1 leads to an overall improvement, compared to the original CAM3.     

Decadal-timescale changes of the Atlantic overturning circulation and climate in a coupled climate model with a hybrid-coordinate ocean component

A wide range of statistical tools is used to investigate the decadal variability of the Atlantic Meridional Overturning Circulation (AMOC) and associated key variables in a climate model (CHIME, Coupled Hadley-Isopycnic Model Experiment), which features a novel ocean component. CHIME is as similar as possible to the 3rd Hadley Centre Coupled Model (HadCM3) with the important exception that its ocean component is based on a hybrid vertical coordinate. Power spectral analysis reveals enhanced AMOC variability for periods in the range 15–30 years. Strong AMOC conditions are associated with: (1) a Sea Surface Temperature (SST) anomaly pattern reminiscent of the Atlantic Multi-decadal Oscillation (AMO) response, but associated with variations in a northern tropical-subtropical gradient; (2) a Surface Air Temperature anomaly pattern closely linked to SST; (3) a positive North Atlantic Oscillation (NAO)-like pattern; (4) a northward shift of the Intertropical Convergence Zone. The primary mode of AMOC variability is associated with decadal changes in the Labrador Sea and the Greenland Iceland Norwegian (GIN) Seas, in both cases linked to the tropical activity about 15 years earlier. These decadal changes are controlled by the low-frequency NAO that may be associated with a rapid atmospheric teleconnection from the tropics to the extratropics. Poleward advection of salinity anomalies in the mixed layer also leads to AMOC changes that are linked to processes in the Labrador Sea. A secondary mode of AMOC variability is associated with interannual changes in the Labrador and GIN Seas, through the impact of the NAO on local surface density.     

Towards a better understanding of changes in wintertime cold extremes over Europe: a pilot study with CNRM and IPSL atmospheric models

Recent winter seasons have evidenced that global warming does not exclude the occurrence of exceptionally cold and/or snowy episodes in the Northern mid-latitudes. The expected rarefaction of such events is likely to exacerbate both their societal and environmental impacts. This paper therefore aims to evaluate model uncertainties underlying the fate of wintertime cold extremes over Europe. Understanding why climate models (1) still show deficiencies in simulating present-day features and (2) differ in their responses under future scenarios for the twentyfirst century indeed constitutes a crucial challenge. Here we propose a weather-regime approach in order to separate the contributions of large-scale circulation and non-dynamical processes to biases or changes in the simulated mean and extreme temperatures. We illustrate our methodology from the wintertime occurrence of extremely cold days in idealized atmosphere-only experiments performed with two of the CMIP5 climate models (CNRM-CM5 and IPSL-CM5A-LR). First we find that most of the present-day temperature biases are due to systematic errors in non-dynamical processes, while the main features of the large-scale dynamics are well captured in such experiments driven by observed sea-surface temperatures, with the exception of a generalized underestimation of blocking episodes. Then we show that uncertainties associated with changes in large-scale circulation modulate the depletion in cold extremes under an idealized scenario for the late twentyfirst century. These preliminary results suggest that the original methodology proposed in this paper can be helpful for understanding spreads of larger model-ensembles when simulating the response of temperature extremes to climate change.     

Combining ERBE and ISCCP data to assess clouds in the Hadley Centre, ECMWF and LMD atmospheric climate models

 This study compares radiative fluxes and cloudiness fields from three general circulation models (the HadAM4 version of the Hadley Centre Unified model, cycle 16r2 of the ECMWF model and version LMDZ 2.0 of the LMD GCM), using a combination of satellite observations from the Earth Radiation Budget Experiment (ERBE) and the International Satellite Cloud Climatology Project (ISCCP). To facilitate a meaningful comparison with the ISCCP C1 data, values of column cloud optical thickness and cloud top pressure are diagnosed from the models in a manner consistent with the satellite view from space. Decomposing the cloud radiative effect into contributions from low-medium- and high-level clouds reveals a tendency for the models' low-level clouds to compensate for underestimates in the shortwave cloud radiative effect caused by a lack of high-level or mid-level clouds. The low clouds fail to compensate for the associated errors in the longwave. Consequently, disproportionate errors in the longwave and shortwave cloud radiative effect in models may be taken as an indication that compensating errors are likely to be present. Mid-level cloud errors in the mid-latitudes appear to depend as much on the choice of the convection scheme as on the cloud scheme. Convective and boundary layer mixing schemes require as much consideration as cloud and precipitation schemes when it comes to assessing the simulation of clouds by models. Two distinct types of cloud feedback are discussed. While there is reason to doubt that current models are able to simulate potential `cloud regime' type feedbacks with skill, there is hope that a model capable of simulating potential `cloud amount' type feedbacks will be achievable once the reasons for the remaining differences between the models are understood. Received: 23 January 2000 / Accepted: 24 January 2001     

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.