Simulating Clouds with Global Climate Models - AMS Journals

4 downloads 36 Views 11MB Size Report
Jun 1, 2013 - (201. 0); http://w ww. lasg.ac.cn/F. GOAL. S/CMIP5. CC. C. C. GFDL. High. Resolu tion. A tmosph eric. Model. C180. (GFDL-HiRAM-C180). NOA.
1 JUNE 2013

LAUER AND HAMILTON

3823

Simulating Clouds with Global Climate Models: A Comparison of CMIP5 Results with CMIP3 and Satellite Data AXEL LAUER International Pacific Research Center, University of Hawaii at Manoa, Honolulu, Hawaii

KEVIN HAMILTON International Pacific Research Center, and Department of Meteorology, University of Hawaii at Manoa, Honolulu, Hawaii (Manuscript received 16 July 2012, in final form 1 November 2012) ABSTRACT Clouds are a key component of the climate system affecting radiative balances and the hydrological cycle. Previous studies from the Coupled Model Intercomparison Project phase 3 (CMIP3) showed quite large biases in the simulated cloud climatology affecting all GCMs as well as a remarkable degree of variation among the models that represented the state of the art circa 2005. Here the progress that has been made in recent years is measured by comparing mean cloud properties, interannual variability, and the climatological seasonal cycle from the CMIP5 models with satellite observations and with results from comparable CMIP3 experiments. The focus is on three climate-relevant cloud parameters: cloud amount, liquid water path, and cloud radiative forcing. The comparison shows that intermodel differences are still large in the Coupled Model Intercomparison Project phase 5 (CMIP5) simulations, and reveals some small improvements of particular cloud properties in some regions in the CMIP5 ensemble over CMIP3. In CMIP5 there is an improved agreement of the modeled interannual variability of liquid water path and of the modeled longwave cloud forcing over mid- and high-latitude oceans with observations. However, the differences in the simulated cloud climatology from CMIP3 and CMIP5 are generally small, and there is very little to no improvement apparent in the tropical and subtropical regions in CMIP5. Comparisons of the results from the coupled CMIP5 models with their atmosphere-only versions run with observed SSTs show remarkably similar biases in the simulated cloud climatologies. This suggests the treatments of subgrid-scale cloud and boundary layer processes are directly implicated in the poor performance of current GCMs in simulating realistic cloud fields.

1. Introduction The simulation of clouds with global climate models (GCMs) involves many nonlinear processes spanning a large range of spatial and temporal scales. While cloud systems such as the subtropical stratocumulus decks can extend over thousands of kilometers, cloud droplet formation and droplet growth occur on the micrometer scale. Similarly, time scales relevant to clouds and cloud microphysics range from weeks to fractions of seconds. All of this makes the simulation of clouds with climate models very difficult (Solomon et al. 2007). Because of the large impact of clouds on the radiation budget and

Corresponding author address: Axel Lauer, IPRC, SOEST, University of Hawaii at Manoa, 1680 East–West Rd., POST Bldg. 401, Honolulu, HI 96822. E-mail: [email protected] DOI: 10.1175/JCLI-D-12-00451.1 Ó 2013 American Meteorological Society

their pivotal role in the hydrological cycle, even small changes in cloud properties could have a significant impact on climate (e.g., Hartmann and Doelling 1991). Clouds and their response to climate change therefore remain a major source of uncertainty for projections of the climate response to anticipated anthropogenic forcing (e.g., Cess et al. 1990; Bony and Dufresne 2005; Stowasser et al. 2006; Solomon et al. 2007; Medeiros et al. 2008; Lauer et al. 2010). Even for simulations of the long-term mean cloud fields under present-day conditions, GCMs generally display rather large deviations from observations, and there are quite large disagreements among GCMs in various aspects of their cloud climatology simulations (Weare 2004; Zhang et al. 2005; Waliser et al. 2007, 2009; Lauer et al. 2010; Chen et al. 2011). Deficiencies in GCM representation of cloud fields will affect the application of such models to simulate chemistry–climate interactions (e.g., by causing biases in photolysis rates)

Second-generation Canadian Earth System Model Community Climate System Model, version 4 Me´te´o-France/Centre National de Recherches Me´te´orologiques Coupled Global Climate Model, version 5

CanESM2*

GISS-E2-R*

NASA Goddard Institute for Space Studies (GISS), United States

Me´te´o-France/Centre National de Recherches Me´te´orologiques (CNRM) and Centre Europe´en de Recherches et de Formation Avance´e en Calcul Scientifique (CERFACS), France Australian Commonwealth Scientific and Industrial Research Organization (CSIRO) Marine and Atmospheric Research, Queensland Climate Change Centre of Excellence (QCCCE), Australia National Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics (LASG)/Institute of Atmospheric Physics, China NOAA/Geophysical Fluid Dynamics Laboratory (GFDL), United States

Canadian Center for Atmospheric Research, Canada NCAR, United States

Beijing Climate Center, China Meteorological Administration, China

Host institute

C

C

C C

C C

Schmidt et al. (2006)

2.58 3 28, L40

C

C

2.58 3 28, L24

C

C

s

C

C

C

C

C

2.58 3 28, L24

s

2.88 3 38, L26

Donner et al. (2011); http://nomads.gfdl. noaa.gov/

Li et al. (2010); http://www.lasg.ac.cn/ FGOALS/CMIP5

1.98 3 1.98, L18

C

2.58 3 2.58, L48

Rotstayn et al. (2010)

1.48 3 1.48, L31

C

C

Voldoire et al. (2013)

1.38 3 0.98, L26

C

C

C

C

CLT

LWP

2.88 3 1.78, L26

Gent et al. (2011)

2.88 3 2.88, L35

References Wu et al. (2010); Wu et al. (2012, manuscript submitted to Climate Dyn.) Arora et al. (2011)

2.88 3 2.88, L26

Resolution (atmosphere)

C

C

C

C

C

C

C

C

C

C

s

C

s C

C

C

C

C

C

LCF

C

C

C

C

C

SCF

C

C

C

C

C

s

s

C

C

C

C

s

SST

JOURNAL OF CLIMATE

GISS-E2-H

GFDL-ESM2M

GFDL-ESM2G

GFDL-CM3

Geophysical Fluid Dynamics Laboratory Climate Model, version 3 Geophysical Fluid Dynamics Laboratory Earth System Model 2G Geophysical Fluid Dynamics Laboratory Earth System Model 2M Goddard Institute for Space Studies Model E2-H Goddard Institute for Space Studies Model E2-R

Flexible Global Ocean–Atmosphere– Land System Model, gridpoint version 2 Flexible Global Ocean–Atmosphere– Land System Model, version 2

FGOALS-g2

FGOALS-s2*

Commonwealth Scientific and Industrial Research Organisation, Mark version 3.6.0

CSIRO-Mk3.6.0*

CNRM-CM5*

CCSM4

Beijing Climate Center Climate System Model, version 1.1

Full model name

BCC-CSM1.1*

Model(s)

TABLE 1. CMIP5 coupled models providing data for this study (WCRP CMIP5 multimodel dataset). The black circles specify available data and the open circles specify where no data are available. The models marked with an asterisk also provide data from corresponding runs with prescribed SSTs (AMIP experiments).

3824 VOLUME 26

NorESM1-M*

MRI-CGCM3*

MPI-ESM-P

MPI-ESM-LR*

MIROC5*

MIROC4h

MIROC-ESM-CHEM

MIROC-ESM

IPSL-CM5B-LR

IPSL-CM5A-MR

IPSL-CM5A-LR*

INM-CM4*

HadGEM2-ES*

HadGEM2-CC

HadCM3

Model(s)

Hadley Centre Coupled Model, version 3 Hadley Centre Global Environmental Model, version 2–Carbon Cycle Hadley Centre Global Environmental Model, version 2–Earth System Institute of Numerical Mathematics Coupled Model, version 4 L’Institut Pierre-Simon Laplace Coupled Model version 5A, low resolution L’Institut Pierre-Simon Laplace Coupled Model version 5A, medium resolution L’Institut Pierre-Simon Laplace Coupled Model version 5B, low resolution Model for Interdisciplinary Research on Climate Earth System Model Model for Interdisciplinary Research on Climate Earth System Model, atmospheric chemistry coupled version Model for Interdisciplinary Research on Climate 4, high-resolution version Model for Interdisciplinary Research on Climate 5 Max Planck Institute Earth System Model, low resolution Max Planck Institute Earth System Model, low-resolution grid and paleo mode Meteorological Research Institute Coupled General Circulation Model version 3 Norwegian Earth System Model 1, medium resolution

Full model name

Norwegian Climate Centre, Norway

Meteorological Research Institute (MRI), Japan

Max Planck Institute for Meteorology, Germany

Japan Agency for Marine-Earth Science and Technology (JAMSTEC), Atmosphere and Ocean Research Institute (AORI), University of Tokyo, and National Institute for Environmental Studies (NIES), Japan

Institute for Numerical Mathematics (INM), Russia L’Institut Pierre-Simon Laplace (IPSL), France

Met Office Hadley Centre, United Kingdom

Host institute

TABLE 1. (Continued)

Hourdin et al. (2013)

2.58 3 0.68, L39

Yukimoto et al. (2011)

Kirkeva˚g et al. (2013)

2.58 3 1.98, L26

C

C

C

C 1.18 3 1.18, L48

C

C Roeckner et al. (2006); Wetzel et al. (2006)

1.98 3 1.98, L47

C

C

C

C

1.48 3 1.48, L40

C

C

C

C

C

C

C

C

C

C

C

CLT

0.68 3 0.68, L56

C

C

2.88 3 2.88, L80 2.88 3 2.88, L80

C

C

3.88 3 1.98, L39

Watanabe et al. (2011); Sakamoto et al. (2012) Watanabe et al. (2010)

Dufresne et al. (2013)

3.88 3 1.98, L39

C

C

Volodin et al. (2010)

28 3 1.58, L21

C

s

LWP

s

Collins et al. (2001); Collins et al. (2008) Collins et al. (2011)

References

1.88 3 1.38, L38

1.98 3 1.38, L60

3.88 3 2.58, L19

Resolution (atmosphere)

C

C

C

C

C

C

C

C

C

C

C

C

C

C

C

SCF

C

C

C

C

C

C

C

C

C

C

C

C

C

C

C

LCF

C

C

C

C

C

s

s

s

C

C

C

C

C

C

C

SST

1 JUNE 2013 LAUER AND HAMILTON

3825

3826

JOURNAL OF CLIMATE

and, specifically, to simulate aerosol indirect climate forcing (e.g., underestimation in cloud frequency could underestimate indirect aerosol effect). Significant biases in the simulation of present-day clouds also raise concerns about the accurate representation of cloud feedback processes in climate change projections. In this study we investigate the performance of stateof-the-art global coupled GCMs from the fifth phase of the Coupled Model Intercomparison Project (CMIP5) (Taylor et al. 2012) by comparing simulated cloud properties with observations. To address how much progress has been made over the recent years, we compare the CMIP5 models to results from the previous generation of global coupled models (CMIP3). Four cloud parameters are of particular interest as they largely determine the impact of clouds on the radiation budget and climate in the model simulations: total cloud amount (CA), liquid water path (LWP), ice water path (IWP), and top of the atmosphere (ToA) cloud forcing (CF). Both LWP and IWP contribute to CA and CF. Here we focus on CA, LWP, and CF, which we will partition between shortwave cloud forcing (SCF) and longwave cloud forcing (LCF). All three of these cloud parameters are also available from the CMIP3 models, allowing for a side-by-side comparison of CMIP5 with the previous-generation models. This enables us to assess the progress of simulating clouds with global coupled models. As we will demonstrate, both the CMIP3 and CMIP5 models display a great deal of disagreement among models in their cloud climatology simulations. Models developed at individual centers may have changed greatly between their CMIP3 and CMIP5 versions, so we do not try to track improvements of individual models. Results from the ensemble of all CMIP models are used to document the general progress in simulating clouds in comprehensive GCMs. We will also concentrate on the geographical patterns of vertically integrated cloud parameters. A more detailed evaluation of vertical profiles of water vapor and cloud water content from individual CMIP5 models can be found in Jiang et al. (2012). We will show that even the state-of-the-art CMIP5 models display substantial regional biases in their cloud climatology simulations. These coupled GCMs also all display considerable regional biases in the basic aspects of the climate, as evidenced particularly by their simulation of sea surface temperature (SST). It is natural to ask how much the deficiencies in the cloud simulations are directly related to the biases in SST simulations. We performed a side-by-side comparison of results from the CMIP5 models’ coupled ocean–atmosphere runs with their counterparts using prescribed SSTs based on observations. This comparison allowed us to assess the role of biases in simulated SSTs in mean biases and the

VOLUME 26

intermodel spread of the cloud properties investigated here. Section 2 describes the models, the present-day model simulations, and the satellite data used for this intercomparison. The performance metrics and results from the comparison of the CMIP5 models with CMIP3 results and satellite observations are presented in section 3. Also discussed in section 3 is the comparison between coupled model simulations and those with prescribed SSTs. The main conclusions are given in section 4.

2. Models, model simulations, and satellite data We analyze the results of 27 CMIP5 models (Table 1) from the ‘‘historical’’ model runs (twentieth-century simulations for 1850–2005 conducted with the best record of natural and anthropogenic climate forcing), the results from 18 CMIP5 models run with prescribed observed SSTs [Atmospheric Model Intercomparison Project (AMIP) runs; Table 2] (Taylor et al. 2012), and the results from 23 CMIP3 models (Table 3) for the twentieth-century runs with natural and anthropogenic forcings (20C3M experiments) (Meehl et al. 2007). The model data were obtained from the World Climate Research Programme’s (WCRP) CMIP3 and CMIP5 data archive, which is operated by the Program for Climate Model Diagnosis and Intercomparison (PCMDI). The model data used for this comparison cover 20-yr time periods from 1986 to 2005 for CMIP5 and from 1980 to 1999 for CMIP3. These time periods have been chosen to allow for a maximum overlap with available satellite data, particularly with global observations of liquid water path (O’Dell et al. 2008), which are available for 1988– 2007. The satellite data, instruments, and time periods used for this comparison are summarized in Table 4. We picked the three satellite datasets for LWP, CA, and CF listed in Table 4 particularly for their long time coverage of at least 20 years. A long record is needed for an evaluation of cloud properties on a statistical basis because individual years from the global coupled models do not correspond to any specific observed year. Multiyear climate modes such as ENSO can have a strong impact on cloud properties (e.g., Park and Leovy 2004). For the purpose of this study, the long time series of global radiative fluxes from the International Satellite Cloud Climatology Project flux data (ISCCP-FD), which are available since 1983, seem preferable to shorter time series like those from the Earth Radiation Budget Experiment (ERBE, 1984–1990) (Barkstrom 1984) or the Clouds and Earth’s Radiant Energy System (CERES, 2000–11) (Wielicki et al. 1996). (Note that the CMIP5 historical runs are only required to extend through 2005

Voldoire et al. (2013) Rotstayn et al. (2010)

Li et al. (2010); http://www. lasg.ac.cn/FGOALS/CMIP5 Zhao et al. (2009)

1.48 3 1.48, L31 1.98 3 1.98, L18 2.88 3 1.78, L26

Volodin et al. (2010) Dufresne et al. (2013); Hourdin et al. (2013) Watanabe et al. (2011); Sakamoto et al. (2012); Watanabe et al. (2010) Roeckner et al. (2006); Wetzel et al. (2006)

28 3 1.58, L21 3.88 3 1.98, L39

Norwegian Climate Centre, Norway

C C Kirkeva˚g et al. (2013)

1.18 3 1.18, L48 2.58 3 l.98, L26

s s

0.68 3 0.68, L64

MRI, Japan

C C

C

C C

0.28 3 0.28, L64

1.98 3 1.98, L47; 1.98 3 1.98, L47

Max Planck Institute for Meteorology, Germany

MPI-ESM-LR* Max Planck Institute Earth System Model, medium resolution (MPI-ESM-MR) Meteorological Research Institute Atmospheric General Circulation Model 3.2H (MRI-AGCM3.2H) Meteorological Research Institute Atmospheric General Circulation Model 3.2S (MRI-AGCM3.2S) MRI-CGCM3* NorESMl-M*

Mizuta et al. (2012); Yukimoto et al. (2011)

Collins et al. (2001); Collins et al. (2008)

1.98 3 1.28, L38

1.48 3 1.48, L40

s

Schmidt et al. (2006)

2.58 3 28, L40

C

C

C

C

C

C

C

C

LWP

0.38 3 0.258, L32

0.68 3 0.58, L32

2.88 3 2.88, L35

Wu et al. (2010); Wu et al. (2012, manuscript submitted to Climate Dyn.) Arora et al. (2011)

References

2.88 3 2.88, L26

Resolution (atmosphere)

JAMSTEC, AORI, University of Tokyo, and NIES, Japan

INM, Russia IPSL, France

Goddard Institute for Space Studies (NASA/GISS), United States Met Office Hadley Centre, United Kingdom

Beijing Climate Center, China Meteorological Administration, China Canadian Center for Atmospheric Research, Canada Me´te´o-France/CNRM and CERFACS, France CSIRO Marine and Atmospheric Research, Queensland Climate Change Centre of Excellence LASG/Institute of Atmospheric Physics, China NOAA/GFDL, United States

Host institute

MIROC5*

Hadley Centre Global Environmental Model, version 2–Atmosphere (HadGEM2-A)* INM-CM4* IPSL-CM5A-LR*

GFDL High Resolution Atmospheric Model C180 (GFDL-HiRAM-C180) GFDL High Resolution Atmospheric Model C360 (GFDL-HiRAM-C360) GISS-E2-R*

FGOALS-s2*

CSIRO-Mk3.6.0*

Canadian Atmospheric Model, version 4 (CanAM4)* CNRM-CM5*

BCC-CSM1.1*

Model(s)

C C

C

C

C C

C

C C

C

C

C

C

C

C

C

C

C

CLT

C C

C

C

C C

C

C C

C

C

C

C

C

C

C

C

C

SCT

C C

C

C

C C

C

C C

C

C

C

C

C

C

C

C

C

LCF

TABLE 2. CMIP5 models providing data from simulations with prescribed SSTs (AMIP experiments). The closed circles specify available data, and the open circles specify where no data are available. The models marked with an asterisk also provide data from corresponding coupled runs.

1 JUNE 2013 LAUER AND HAMILTON

3827

Community Climate System Model, version 3 Parallel Climate Model 1 Coupled Global Climate Model, version 2.3.2 Third-generation Coupled Global Climate Model with T-47 spectral resolution Third-generation Coupled Global Climate Model with T-63 spectral resolution Me´te´o-France/Centre National de Recherches Me´te´orologiques Coupled Global Climate Model Commonwealth Scientific and Industrial Research Organisation, Mark version 3.0 Commonwealth Scientific and Industrial Research Organisation, Mark version 3.5 ECHAM5/Max Planck Institute Ocean Model Flexible Global Ocean–Atmosphere– Land System Model gridpoint version 1.0 Geophysical Fluid Dynamics Laboratory Climate Model, version 2.0 Geophysical Fluid Dynamics Laboratory Climate Model, version 2.1 Goddard Institute for Space Studies Atmosphere–Ocean Model Goddard Institute for Space Studies Model E-H Goddard Institute for Space Studies Model E-R

CCSM3

GISS-E-R

NASA GISS, United States

NOAA/GFDL, United States

Max Planck Institute for Meteorology, Germany LASG/Institute of Atmospheric Physics, China

CSIRO Marine and Atmospheric Research, Australia

Me´te´o-France/CNRM, France

Canadian Centre for Climate Modeling and Analysis, Canada

MRI, Japan

Bjerknes Centre for Climate Research and Geophysical Institute, University of Bergen, Nansen Environmental and Remote Sensing Center, Norway NCAR, United States

Host institute

GFDL Global Atmospheric Model Development Team (2004)

2.08 3 2.58, L24

48 3 58, L20

Schmidt et al. (2006); Bleck (2002); Russell et al. (1995)

Yu et al. (2004)

2.88 3 2.88, L26

48 3 38, L12

Roeckner et al. (2003)

1.98 3 1.98, L31

C C

C C

C

C

C s

C

C C

C

C

C

C C

C Gordon et al. (2002)

1.98 3 1.98, L18

C

C

C

C C

C

C

CLT

C

C Salas-Me´lia et al. (2005)

1.98 3 1.98, L45

C

s s

C

McFarlane et al. (1992); Flato (2005)

Yukimoto et al. (2006)

C

C

LWP

1.98 3 1.98, L31

3.88 3 3.88, L31

2.88 3 2.88, L26 2.88 3 2.88, L30

Collins et al. (2006); Washington et al. (2000)

1.48 3 1.48, L26

References http://bjerknes.uib.no

2.88 3 2.88, L31

Resolution (atmosphere)

C

C

s

C

C

C

C

C

C

C

C

C

C C

C

s

SCP

C

C

s

C

C

C

C

C

C

C

C

C

C C

C

s

LCF

JOURNAL OF CLIMATE

GISS-E-H

GISS-AOM

GFDL-CM2.1

GFDL-CM2.0

FGOALS-g1.0

ECHAM5/MPI-OM

CSIRO-Mk3.5

CSIRO-Mk3.0

CNRM-CM3

CGCM3.1(T63)

CGCM3.1(T47)

PCM1 CGCM2.3.2

Bjerknes Centre for Climate Research Bergen Climate Model version 2.0

Full model name

BCCR-BCM2.0

Model(s)

TABLE 3. CMIP3 models providing data for this study (WCRP CMIP3 multimodel dataset; Meehl et al. 2007). The closed circles specify available data, and the open circles specify where no data are available.

3828 VOLUME 26

C C C

Met Office Hadley Centre, United Kingdom

1.38 3 1.98, L38

Gordon et al. (2000); Pope et al. (2000); Martin et al. (2006)

C

C C C 2.58 3 3.88, L19

C

C C C C 1.18 3 1.18, L56

C C C Hasumi and Emori (2004) 2.88 3 2.88, L20

Center for Climate System Research (University of Tokyo), NIES, and Frontier Research Center for Global Change (JAMSTEC), Japan

s

C C C 2.58 3 3.88, L19 IPSL, France

C

C C C

Alekseev et al. (1998); Galin et al. (2003) Hourdin et al. (2006) 58 3 48, L21 INM, Russia

C

s s

UKMO-HadGEM1

UKMO-HadCM3

MIROC3.2(hires)

MIROC3.2(medres)

IPSL-CM4

INM-CM3.0

SCP CLT

C s

LWP References

Roeckner et al. (1996) 1.18 3 1.18, L19 Instituto Nazionale di Geofisica e Vulcanologia, Italy

Instituto Nazionale di Geofisica e Vulcanologia Scale Interaction Experiment Institute of Numerical Mathematics Coupled Model, version 3.0 L’Institut Pierre-Simon Laplace Coupled Model, version 4 Model for Interdisciplinary Research on Climate 3.2, medium-resolution version Model for Interdisciplinary Research on Climate 3.2, high-resolution version Met Office Hadley Centre Coupled Model, version 3 Met Office Hadley Centre Global Environmental Model version 1

Resolution (atmosphere) Host institute Full model name Model(s)

TABLE 3. (Continued)

LAUER AND HAMILTON

INGV-SXG

LCF

1 JUNE 2013

3829

and the CMIP3 20C3M runs only through 1999.) The LWP climatology from O’Dell et al. (2008) has been chosen not only for the number of years with observations available but also because O’Dell et al. adjusted their results to be representative of the diurnal mean. The diurnal cycle of LWP can be important when comparing modeled monthly means with satellite data available only at certain times of the day (e.g., overpass time at 1030 local time), as neglecting it might introduce errors of unknown magnitude. The uncertainties in the observational data given in Table 4 are estimated as the systematic errors, assuming that statistical errors are small when averaging over 20 years of data. O’Dell et al. (2008) give potential systematic errors in the LWP climatology of the order of 15%–30%. The uncertainty estimate for monthly regional average ToA fluxes from ISCCP-FD is 5–10 W m22 (Zhang et al. 2004), translating into 25%–50% for the ToA cloud forcing. Total cloud amount from ISCCP is estimated to be about 10% too low over land and ‘‘about right’’ over the ocean (Rossow and Schiffer 1999). With exception of the data shown in Figs. 1 (and 4) (which show results on the native model grids), all model data have been regridded to a common l8 3 l8 grid for comparison with observations and for calculating multimodel means using bilinear interpolation. If there are multiple ensemble members available for any given model and model experiment, we only consider the first ensemble member in our analysis: CMIP3 is runl and CMIP5 is rlilpl. This facilitates the comparison as the number of available ensemble members varies greatly among the individual models and among CMIP3 and CMIP5. An analysis of the spread in the 20-yr-mean LWP among the ensemble members of individual models shows that its relative standard deviation is typically small and ranges between 5% and 15% in regions with frequent deep convection and between 2% and 5% otherwise. The interensemble spread for LCF and SCF is even smaller. This interensemble spread is much smaller than the intermodel spread. We therefore do not expect our results to change significantly when using ensemble members other than runl and rlilpl or when using ensemble means.

3. Comparison of CMIP5 with CMIP3 and satellite data a. Multiyear annual mean Figure 1 shows the 20-yr annual mean liquid water path averaged over the years 1986–2005 from 24 CMIP5 models that have LWP as an available variable. These maps are compared with the multimodel mean and with the University of Wisconsin LWP climatology (UWisc).

3830

JOURNAL OF CLIMATE

VOLUME 26

TABLE 4. Satellite data used for comparison. Parameter

Satellite instrument(s)

Years

Estimate of systematic errors

Reference

LWP CA SCF, LCF

SSM/I, TMI, AMSR-E ISCCP ISCCP-FD

1988–2007 1988–2007 1988–2007

15%–30% 10% 25%–50% (5–10 W m22)

O’Dell et al. (2008) Rossow et al. (1996) Zhang et al. (2004)

UWisc is based on satellite observations from the Special Sensor Microwave Imager (SSM/I), the Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI), and the Advanced Microwave Scanning Radiometer for Earth Observing System (EOS) (AMSR-E) in the years 1988–2007. The model LWP values are obtained by subtracting the vertically integrated ice water path (clivi) from the vertically integrated cloud water path (ice 1 liquid) clwvi. However, some models provided only

the vertically integrated liquid water path as clwvi. Some of these models are listed on the PCMDI CMIP5 errata page (http://cmip-pcmdi.llnl.gov/cmip5/errata/ cmip5errata.html). In addition to the models listed there, we found that subtracting the ice water path from the cloud water path for EPSL-CM5B-LR results in significant negative LWP values (monthly means of individual grid cells range up to 2300 g m22). We therefore assume that clwvi from IPSL-CM5B-LR is also liquid water only.

FIG. 1. The 20-yr average LWP (1986–2005) from the CMIP5 historical model runs and the multimodel mean in comparison with the UWisc satellite climatology (1988–2007) based on SSM/I, TMI, and AMSR-E (O’Dell et al. 2008).

1 JUNE 2013

LAUER AND HAMILTON

We treat clwvi from the following 9 out of 24 CMIP5 models as liquid water only: CCSM4, CSIRO-Mk3.6.0, IPSL-CM5A-MR, IPSL-CM5A-LR, IPSL-CM5B-LR, MIROC-ESM, MIROC-ESM-CHEM, MPI-ESM-LR, and MPI-ESM-P. Similarly, clwvi from the following CMIP3 models is assumed to be liquid water only: BCCRBCM2.0, CSIRO-Mk3.0, and CSIRO-Mk3.5 (Jiang et al. 2012). The models show a large spread in the overall LWP amounts as well as in the geographical pattern of the simulated LWP. Averaged over all grid cells with available satellite observations, the mean LWP values in the models range from 37 to 167 g m22 with a median of 83 g m22. The mean LWP from satellite observations is 85 g m22 (625 g m22 assuming an uncertainty of 630%) and that from the multimodel mean is 87 g m22. The intermodel spread among the CMIP5 models has narrowed down compared with the CMIP3 models, which give global annual average LWP values ranging between 28 and 209 g m22. The CMIP3 median is 73 g m22, and the multimodel mean is 78 g m22. Interestingly, BCCCSM1.1, CCSM4, and NorESMl-M show similar biases in the simulated annual mean LWP, with a strong overestimation in middle and higher latitudes. All three of these models have in common that the atmosphere component of the model is based on National Center for Atmospheric Research (NCAR)’s Community Atmosphere Model version 5 (CCSM4; Gent et al. 2011) or its earlier version CAM3 [BCC-CSM1.1 (Wu et al. 2010); NorESMl-M (Kirkeva˚g et al. 2013)]. This suggests that the biases in the simulated LWP climatology are determined to a large degree by the atmosphere component of the model. The linear pattern correlation of annual mean LWP between the individual CMIP5 models and the satellite observations ranges from a low of 0.03 to a high of 0.77. It is remarkable that some models simulate annual mean LWP fields with almost no geographical correlation with observations. The correlations of individual CMIP3 model simulations with observations range over 0.22– 0.70. The correlation of the CMIP5 multimodel mean (0.59) did not improve over CMIP3 (0.64). A measure of the performance of the CMIP model ensemble in reproducing observed mean cloud properties is obtained by calculating the differences in modeled (xmod) and observed (xobs) 20-yr means. These differences are then averaged over all N models in the CMIP3 or CMIP5 ensemble to calculate the multimodel ensemble mean bias Dmm, which is defined at each grid point as Dmm 5

1 N

N

2 xobs ) . å (xmod i

i51

(1)

3831

Figure 2 shows 20-yr annual means for liquid water path, total cloud amount, and ToA CF from satellite observations and the ensemble mean bias of the CMIP3 and CMIP5 models. The most striking overall impression from this figure is how similar the geographical distribution of biases are in the CMIP3 and CMIP5 ensemble means for all four cloud properties shown. Both the CMIP3 and the CMIP5 ensemble means overestimate LWP in the midlatitude storm track regions, a bias that is apparent in many of the individual models as well (Fig. 1). This overestimation is not reduced in CMIP5 compared with CMIP3. Expressed as a fraction of the observed values in most parts of the CMIP5 storm track regions, the LWP biases are 30%–50% compared with 20%–40% in CMIP3. The negative bias in the Pacific intertropical convergence zone (ITCZ) in the CMIP3 models (230% to 250% of the observed) is slightly improved in CMIP5 (220% to 250%), particularly in the central Pacific. The systematic underestimation of LWP in the stratocumulus regions, particularly off the coasts of South America (240% to 250%), South Africa (230%), and west Australia (210% to 220%), however, did not improve in CMIP5 with mean biases similar to those in CMIP3. The representation of LWP in the South Pacific convergence zone (SPCZ) did not improve in CMIP5 compared with the previous generation of coupled climate models. In contrast to the overestimation of LWP, the models generally underestimate total CA in low and middle latitudes. While there are small improvements in reproducing the observed total CA in high latitudes in the CMIP5, underestimation of CA, particularly in midlatitudes, is found to be slightly larger in CMIP5 (220%) than in CMIP3 (210% to 220%). The underestimation of the cloud amount in the stratocumulus regions off the subtropical west coasts of the continents is similar in the CMIP5 and CMIP3 ensemble mean results, and the bias ranges between 230% and 250% of the observed CA. It is noteworthy that the ISCCP data used here are known to often show an artifact over the Indian Ocean between 508 and 1008E. This artifact in total cloud amount is related to a dependence of the retrieved data on the satellite zenith angle (Rossow and Garder 1993). The CF is defined as the difference between ToA allsky and clear-sky outgoing radiation in the solar spectral range (SCF) and in the thermal spectral range (LCF). A negative CF corresponds to an energy loss and a cooling effect, and a positive CF corresponds to an energy gain and a warming effect. Biases in annual average SCF are slightly reduced in CMIP5 over the subtropical central North Pacific (30%–50%) and the North Atlantic (20%– 40%), as well as over North Africa (210% to 250%) and South Asia (220% to 260%), compared with CMIP3.

3832

JOURNAL OF CLIMATE

VOLUME 26

FIG. 2. Differences in (top to bottom) 20-yr annual averages of LWP, CA, SCF, and LCF from the (left) CMIP3 and (middle) CMIP5 multimodel means compared with (right) satellite observations. For details, see text.

Very little change between CMIP5 and CMIP3 is found in the stratocumulus regions (230% to 260%), the ITCZ (25% to 225%), the Southern Hemisphere oceans, and the Americas. In contrast to SCF, LCF improved significantly over the oceans, particularly in middle and high

latitudes. The large biases in LCF found in CMIP3 in the Pacific and Atlantic south of the ITCZ and in the SPCZ are still present in CMIP5 even though the magnitude of the biases is reduced (CMIP3: 100%–200%; CMIP5: 50%–150%). The CMIP5 models show a larger

1 JUNE 2013

LAUER AND HAMILTON

underestimation in LCF, particularly over the Americas, Australia, and Asia (220% to 230%), than do the CMIP3 models (210% to 220%). The overall comparisons of the annual mean cloud properties with observations are summarized for individual models and for the ensemble means by the Taylor diagrams for CA, LWP, SCF, and LCF shown in Fig. 3. These give the standard deviation and linear correlation with satellite observations of the total spatial variability calculated from 20-yr annual means. The standard deviations are normalized by the observed values, so the observed climatology is represented in each panel by the black stars on the x axis at x 5 1. In this polar coordinate system, the linear distance between the observations and each model is proportional to the rootmean-square error (rmse) (Taylor 2001) and can be gauged using the green circles centered on the observational dots in Fig. 3 (contour labels in standard deviation). The linear correlation coefficients for total CA among the individual CMIP3 models range from 0.12 to 0.87 (multimodel mean 5 0.76); in CMIP5 this range is from 0.11 to 0.83 (multimodel mean 5 0.80). The rmse of the simulated 20-yr-mean CA ranges from 11% to 20% among the CMIP3 models (multimodel mean 5 11%) and from 10% to 23% (multimodel mean 5 12%) in CMIP5. Total CA is strongly determined by LWP and IWP. Similar to Jiang et al. (2012), we find that the model spread in total CA is smaller than that in LWP (and IWP) (Jiang et al. 2012). A possible reason for that could be that biases in LWP and IWP are partly compensating each other. In addition to the standard model variables, some CMIP5 models provide diagnostics resembling satellite observations of, for instance, total cloud amount, as it would be seen from ISCCP: CanESM2, IPSL-CM5ALR, IPSL-CM5A-MR, MIROC-ESM, and MIROCESM-CHEM. These diagnostics are derived with the Cloud Feedback Model Intercomparison Project (CFMIP) (Bony et al. 2011) integrated satellite simulator, the CFMIP Observation Simulator Package (COSP) (BodasSalcedo et al. 2011). We compared the performance of these five models in reproducing the 20-yr-mean ISCCP observed total cloud amount with and without application of COSP. We found that, while the pattern correlation of CanESM2 improves considerably from 0.67 to 0.81 when using COSP, the performance of the four other models measured by the rmse hardly changes. We therefore do not expect the large intermodel spread in total cloud amount to be reduced considerably when comparing the simulated ISCCP total cloud amount instead of the models’ native total cloud amount, as done here.

3833

Just as for CA, the performance in reproducing the observed multiyear annual mean LWP did not improve considerably in CMIP5 compared with CMIP3. The rmse ranges between 20 and 129 g m22 in CMIP3 (multimodel mean 5 22 g m22) and between 23 and 95 g m22 in CMIP5 (multimodel mean 5 24 g m22). For SCF and LCF, the spread among the models is much smaller compared with CA and LWP. The agreement of modeled SCF and LCF with observations is also better than that of CA and LWP. The linear correlations for SCF range between 0.83 and 0.94 (multimodel mean 5 0.95) in CMIP3 and between 0.80 and 0.94 (multimodel mean 5 0.95) in CMIP5. The rmse of the multimodel mean for SCF is 8 W m22 in both CMIP3 and CMIP5. As in the case for SCF, the CMIP3 and CMIP5 model performance in simulating a realistic LCF climatology is clearly better than for CA and LWP. The correlation of the multimodel mean LCF from CMIP3 is 0.92 (rmse 5 5 W m22), with that of the individual models ranging between 0.66 and 0.90 (rmse 5 6–11 W m22). For CMIP5, the correlation of the multimodel mean LCF is 0.93 (rmse 5 4 W m22) and ranges between 0.70 and 0.92 (rmse 5 4–11 W m22) for the individual models. In both CMIP3 and CMIP5, the large intermodal spread and biases in CA and LWP contrast strikingly with a much smaller spread and better agreement of global average SCF and LCF with observations. The SCF and LCF directly affect the global mean radiative balance of the earth, so it is reasonable to suppose that modelers have focused on ‘‘tuning’’ their results to reproduce aspects of SCF and LCF as the global energy balance is of crucial importance for long climate integrations. The problem of compensating biases in simulated cloud properties is not new and has been reported in previous studies. For example, Zhang et al. (2005) find that many CMIP3 models underestimate cloud amount while overestimating optically thick clouds, explaining why models simulate the ToA cloud forcing reasonably well while showing large biases in simulated LWP (Weare 2004). The problem of too few but too bright clouds has also been identified in some CMIP5 models (e.g., Kay et al. 2012; Nam and Quaas 2012).

b. Coupled models versus prescribed SSTs A first-principles simulation of even the main features of regional climate in coupled atmosphere–ocean models is challenging, as errors can be amplified by positive air– sea coupling feedbacks. Current coupled GCMs typically simulate mean present-day climate with substantial biases from that observed. A well-known example affecting many coupled models is a cold bias in the eastern Pacific equatorial SST accompanied by deficiencies in the overlying atmospheric circulation, along with the appearance

FIG. 3. Taylor diagrams showing the 20-yr annual average performance of the (top) CMIP3 and the (bottom) CMIP5 models for (left to right) CA, LWP, and ToA SCF and LCF as compared to satellite observations. For details, see text.

3834 JOURNAL OF CLIMATE VOLUME 26

1 JUNE 2013

LAUER AND HAMILTON

3835

FIG. 4. Bias in simulated 20-yr-mean SST (1986–2005) from the CMIP5 historical model runs in comparison with the National Oceanic and Atmospheric Administration (NOAA) optimum interpolation daily SST analysis (Reynolds et al. 2007). Also shown are the multimodel means averaged over all 27 models (Table 1) and averaged only over 13 models that also provide data from corresponding AMIP runs (models marked with an asterisk in Table 1).

of an unrealistic ‘‘double ITCZ’’ feature (e.g., Lin 2007). Figure 4 shows the biases in simulated annual mean climatological SST in 26 CMIP5 models. The multimodel mean SST has positive biases of 18–38C in the stratocumulus regions off the west coasts of North and South America, as well as off Africa, and in middle and high latitudes south of 458S. SSTs in the multimodel mean are too low compared with observations in the North Pacific (228C) and the North Atlantic (228 to 248C). The large SST biases in the stratocumulus regions are of particular concern as ocean temperature there may significantly impact the boundary layer cloud fields (e.g., Bony 1997). It is reasonable to ask how much of the large deficiencies we have documented in cloud fields in the CMIP models may be ascribed to the overall biases in the simulated climate in the coupled model runs. Here we compare results from near present day in the CMIP5 coupled model historical runs (Table 1) with results from experiments using the atmosphere components of the CMIP5 models and prescribed observed SSTs (the socalled AMIP runs). As shown in Table 2, 18 of the

CMIP5 models have provided AMIP results we can use. The performance of the models in reproducing observed multiyear mean cloud properties is summarized in the Taylor diagrams for CA, LWP, SCF, and LCF, shown in Fig. 5. As in the case of the coupled CMIP5 models discussed above, we find a large intermodel spread for CA and LWP, while the spread is smaller for SCF and LCF. The corresponding biases and rmse are smaller for the modeled SCF and LCF than for CA and LWP. The intermodel spread in both pattern correlation and rmse does not narrow down when going from the coupled model runs to the AMIP experiments. For LWP, the range of disagreement with observations among individual models is similar in the AMIP runs (linear pattern correlation 5 0.35–0.81, rmse 5 23–96 g m22) than in the coupled simulations (correlation 5 0.03– 0.77, rmse 5 23–95 g m22). The median correlation of the AMIP models with observations is 0.49 (median rmse 5 42 g m22) and that of the coupled models is 0.49 (median rmse 5 43 g m22). Similar results are seen for the other cloud parameters investigated here (CA, SCF,

3836

JOURNAL OF CLIMATE

VOLUME 26

FIG. 5. Taylor diagrams for CA, LWP, and ToA SCF and LCF showing the 20-yr annual average performance of the coupled CMIP5 models (blue and gray dots) in comparison to their counterparts with prescribed SSTs (AMIP runs, red dots) and to satellite observations (black dots). For details, see text.

and LCF). This suggests that the AMIP models do not systematically outperform the coupled models in reproducing observed mean cloud properties. Furthermore, this suggests that the large intermodel spread in CA and LWP is attributable to the representation of

cloud processes rather than to biases in the SST and related aspects of the circulation in the coupled models. Figure 6 shows the biases in simulated cloud properties from the CMIP5 multimodel means computed for 13 coupled models that also provide data from AMIP runs

1 JUNE 2013

LAUER AND HAMILTON

3837

FIG. 6. Differences in (top to bottom) 20-yr annual averages of LWP, CA, SCF, and LCF from (left) the CMIP5 multimodel means averaged over 13 coupled models that also provide data from AMIP runs and (middle) their corresponding AMIP experiments. (right) The differences in bias amplitudes between the coupled and AMIP model runs. Positive values refer to larger biases in the coupled runs, negative values to larger bias amplitudes in the AMIP runs.

(models marked with an asterisk in Table 1) and are averaged over the corresponding AMIP experiments (Table 2). Also shown in the rightmost panels are the differences in the magnitude of the biases between the coupled and the AMIP model runs (jcoupledj 2 jAMIPj). Positive values refer to larger biases in the coupled runs,

and negative values refer to larger bias amplitudes in the AMIP runs. The differences in the bias amplitudes between the coupled and the AMIP model runs are generally small and are typically less than 10% of the modeled values. The geographical patterns of the biases are also very similar in the coupled and AMIP runs.

3838

JOURNAL OF CLIMATE

VOLUME 26

FIG. 7. Biases in simulated 20-yr-mean LWP from (left) the (top to bottom) four individual coupled CMIP5 models and (middle) their AMIP counterparts, with the smallest global average rmse in LWP. (right) The biases in annual mean SST in the coupled runs.

Figure 7 shows the coupled model SST biases and the biases in simulated 20-yr-mean LWP in the coupled and AMIP simulations for four selected individual CMIP5 models. For this purpose, we chose the models that have the smallest rmse in annual mean LWP in both the coupled and the AMIP runs: MIROC5 (rmsecoupled 5 23 g m22, rmseAMIP 5 23 g m22), MRI-CGCM3 (24 g m22, 24 g m22), CNRM-CM5 (30 g m22, 34 g m22), and CSIRO-Mk3.6.0 (35 g m22, 39 g m22). As for the multimodel means, we find the differences in simulated LWP between the coupled and the AMIP simulations to be remarkably small, with differences in the bias magnitudes typically around 10%–20% of the modeled values or less. In the subtropical stratocumulus regions off the coast of North and South America and off the coast of Africa, three of the four coupled models show positive SST biases but similar corresponding negative biases in simulated LWP. Negative SST biases are found over most of the other parts of the subtropical oceans, with smaller bias

magnitudes in LWP in the coupled model runs than in the AMIP experiments. This is particularly the case for the MIROC5 and the MRI-CGCM3 models. Warm biases in SST in southern mid-to-high latitudes also seem to be related to smaller LWP biases in the coupled version of CNRM-CM5 compared with the AMIP run.

c. Interannual variability The interannual variability is estimated as the temporal standard deviation of the deseasonalized monthly means (smod) calculated at each grid cell. The monthly means are normalized by the average over the entire time series (xmod ). The interannual variability of the multimodel ensemble smm is calculated by averaging over all N individual CMIP5 models: smm 5

1 N

N

smod i 3 100%. mod i51 xi

å

(2)

1 JUNE 2013

LAUER AND HAMILTON

3839

FIG. 8. Differences in (top to bottom) the relative temporal standard deviation [Eq. (2)] of 20 years of LWP, CA, SCF, and LCF from the (left) CMIP3 and (middle) CMIP5 models compared with (right) satellite observations. The modeled temporal standard deviations are calculated from monthly anomalies after subtracting the climatological mean seasonal cycle and have been averaged over all models. For details, see text.

Figure 8 shows the relative temporal standard deviation [Eq. (2)] from satellite observations (Table 4) and the deviation of the modeled interannual variability from that of the observations for CA, LWP, SCF, and LCF. The observations show the largest interannual variability in LWP in the ITCZ and SPCZ regions, as well as in the subtropical regions of the central and west Pacific and the Indian Ocean. The modeled interannual variability of LWP in CMIP5 is overestimated throughout most of the globe, ranging between 5% and 15% in the midlatitudes and between 15% and 25% in the subtropics, with the largest overestimation in the subtropical east Pacific. In contrast, the models underestimate the interannual variability in LWP in the equatorial east and central Pacific (210% to 230%). This qualitative pattern is also found in the CMIP3 results, but there is improved agreement with the observed interannual variability of LWP, particularly in midlatitudes and the subtropical

west Pacific and Indian Ocean. Differences between modeled and observed interannual variability in CA, LCF, and SCF are smaller than those in LWP, with the main features of the geographical distribution of the deviations being similar to those for LWP. In contrast to LWP, the differences between the interannual variability in CMIP3 and CMIP5 for CA, SCF, and LCF are only small. The overestimation in interannual variability of CA in CMIP5 is somewhat reduced over the continents compared with CMIP3, but slightly increased over the subtropical oceans. The differences in observed and modeled interannual variability in SCF and LCF in CMIP5 and their changes compared with CMIP3 are qualitatively similar to those in CA.

d. Seasonal cycle We analyze the ability of the models to reproduce the observed geographical patterns of the seasonal cycle of

3840

JOURNAL OF CLIMATE

VOLUME 26

FIG. 9. Differences between the 20-yr seasonal averages for JJA and DJF from (left) CMIP3 and (middle) CMIP5 multimodel means in comparison with (right) satellite observations. Shown are (top to bottom) LWP, CA, SCF, and LCF.

LWP, CA, SCF, and LCF by calculating the differences in seasonal averages over the months June–August (JJA) and December–February (DJF). Figure 9 shows a comparison of the differences in JJA and DJF averages from the CMIP3 and CMIP5 models with satellite observations. Both the CMIP3 and the CMIP5 multimodel means capture the geographical pattern of the JJA–DJF amplitudes reasonably well. The differences in JJA–DJF amplitudes for LWP, CA, and LCF between CMIP3 and CMIP5 are mostly small throughout the globe, with little to no improvement in CMIP5 over CMIP3. In contrast, the amplitudes of the differences in summer and winter averages of SCF from CMIP5

increased in midlatitudes and high latitudes by 10%– 15% compared with CMIP3. The SCF amplitudes from CMIP5 in these regions are in better agreement with satellite observations than those from CMIP3. For each point and each model, a mean seasonal cycle was determined by averaging the data for January over all years, for February over all years, etc. We then normalized this mean seasonal cycle (xj) for each model and point by dividing by the average over all months and years (x). As a performance measure of the shape of the seasonal cycle, we calculate the rmse of the normalized /xmod ) compared modeled mean seasonal cycle (xmod j obs obs with that from observations (xj /x ). The rmse of the

1 JUNE 2013

LAUER AND HAMILTON

3841

FIG. 10. Multimodel mean of the rmse [Eq. (3)] of the mean seasonal cycle of (top to bottom) LWP, CA, SCF, and LCF from comparisons with satellite observations. Results from (left) the CMIP3 models and (right) from CMIP5 models. For details, see text.

multimodel ensemble (rmsemm) is calculated by averaging over all rmse values for the N individual CMIP3 or CMIP5 models:

vffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi !ffi u obs 2 12 xmod u x 1 1 i,j j rmsemm 5 å t å mod 2 obs 3 100%. N i51 12 j51 xi x N

(3)

The largest errors in the multimodel mean seasonal cycle of LWP are found in the tropical Atlantic and east Pacific as well as in the Indian Ocean and the west Pacific (Fig. 10). Here the average rmse ranges between 50% and 80%. The errors in the tropics are largely

unchanged in CMIP5 compared with CMIP3. The errors in the LWP mean seasonal cycle in midlatitudes are much smaller than in the tropics and typically range between 10% and 20%. Again, there is only very little change between CMIP3 and CMIP5 in the models’ ability to reproduce the observed LWP mean seasonal cycle. Consistent with the large errors in reproducing the observed LWP seasonal cycle in the tropics and small errors in midlatitudes, CA, SCF, and LCF show a similar behavior. In general, the error amplitudes are smallest for cloud amount and largest for LCF and LWP. In general, there is little to no improvement in reproducing

3842

JOURNAL OF CLIMATE

the observed seasonal cycles of LWP, CA, SCF, and LCF in CMIP5 compared with CMIP3.

4. Conclusions Comprehensive global models of the kind used for long-term climate projections are known to display significant deficiencies in their cloud simulations. The simulated cloud climate feedbacks activated in global warming projections differ enormously among state-ofthe-art models, and this large degree of disagreement has been a constant feature documented for successive generations of GCMs from the time of the first Intergovernmental Panel on Climate Change (IPCC) assessment (Cess et al. 1990) through the CMIP3 generation models used in the fourth IPCC assessment (Bony and Dufresne 2005; Lauer et al. 2010). Much of the difference among state-of-the-art global models in their simulated climate sensitivity is due to the simulated cloud feedbacks, and the recalcitrance of the cloud simulation issue accounts for the often-lamented fact that the range of GCM-estimated climate sensitivity to increased greenhouse gas concentrations scarcely narrowed from the first IPCC assessment to the fourth IPCC assessment. Even the model-simulated cloud climatologies for present-day conditions are known to depart significantly from observations and, once again, the variation among models is quite remarkable (e.g., Weare 2004; Zhang et al. 2005; Waliser et al. 2007, 2009; Lauer et al. 2010; Chen et al. 2011). In the project reported here, we have evaluated the quality of the cloud climatology simulations in the CMIP5 coupled models, which represent the contemporary state of the art in models to be used for long-term climate projections. The evaluation focused on the geographical patterns of vertically integrated cloud properties, notably, liquid water path, total cloud amount, and shortwave and longwave cloud radiative forcing. We have also tried to assess whether an overall improvement in model performance for these fields can be found relative to that of the CMIP3 models, which represented the state of the art circa 2005. Of course, such a comparison of model performance over different generations is complicated by the fact that many separate changes may be introduced by developers at any particular climate research center over such a long period and improving cloud simulations will be only one of many motivations for such changes. However, given the first-order deficiencies seen in the critically important cloud fields in the global GCMs, it is worthwhile to see if an overall improvement in the stateof-the-art models is apparent. A simple approach to characterizing this overall progress is through use of multimodel mean ensemble results. Possible improvements in performance in the multimodel mean cloud fields

VOLUME 26

are also of interest since the multimodel means are sometimes suggested as, intuitively, the most reliable or probable projections from the CMIP model intercomparison exercises. We found that the long-term mean vertically integrated cloud fields have quite significant deficiencies in all the CMIP5 model simulations. For example, the highest spatial correlation of annual mean LWP with satellite climatology for any of the CMIP5 models is 0.77, and for one of the models the correlation is as small as 0.03. The correlation of the multimodel mean for CMIP5 in this case is 0.59, which is similar to the results for the CMIP3 ensemble of 0.64. Both the CMIP5 and CMIP3 models display a clear bias to simulating too high LWP in midlatitudes. This bias is not reduced in the CMIP5 models. We also investigated the seasonal cycle of cloud properties in the CMIP3 and CMIP5 runs and found significant discrepancies from observations, but only a very small improvement in the regard in CMIP5. Our analysis of the root-mean-square error of simulated LWP, CA, SCF, and LCF supports our findings on little to no changes in the skill of reproducing the observed LWP and CA. In general, the models have higher skills in reproducing observed SCF and LCF than CA and LWP. Of course, since the SCF and LCF directly affect the global mean radiative balance of the earth, efforts by model developers to ‘‘tune’’ their subgridscale parameterizations may focus on these variables. Our study included a detailed comparison of presentday cloud climatologies in the CMIP5 coupled historical runs and in AMIP simulations, in which the atmospheric components are forced with observed SSTs. We found that, even though the historical coupled runs displayed significant mean SST biases, the overall cloud simulations in the coupled and AMIP runs were quite similar. This result seems quite robust, and it is noteworthy the AMIP simulations do not outperform the coupled runs in reproducing observed mean cloud properties. This suggests that the deficiencies in state-of-the-art GCM cloud climatology simulations may be rather directly attributable to the cloud, convection, and boundary layer parameterizations employed. A positive implication of our findings is that some aspects of cloud-related subgridscale parameterizations may be tested adequately through AMIP type simulations and that coupled feedbacks with the ocean may not be first-order significant in this regard. We also conducted an analysis of the interannual variability of the cloud fields in the coupled CMIP3 and CMIP5 simulations. The CMIP5 models showed improvements, particularly in LWP in the midlatitudes and the subtropical west Pacific and Indian Ocean compared with CMIP3. The CMIP5 versus CMIP3 differences in

1 JUNE 2013

LAUER AND HAMILTON

the statistics of interannual variability of SCF and LCF are quite modest, although a systematic overestimation in interannual variability of CA in CMIP3 is slightly improved over the continents in CMIP5. Our analysis showed that there is still a wide intermodel spread, particularly in the modeled LWP and CA. There is generally only very modest improvement in the simulated cloud climatology in CMIP5 compared with CMIP3. The better performance of the models in reproducing observed annual mean SCF and LCF therefore suggests that this good agreement is mainly a result of careful model tuning rather than an accurate fundamental representation of cloud processes in the models. Simulating clouds in comprehensive climate models remains a very challenging problem, and, unfortunately, over recent years the pace of improvement has been modest at best. An important issue not addressed by this study is cloud ice. Cloud ice contributes to the total cloud amount as well as to the cloud forcing. Biases in LWP and IWP can partly compensate each other, potentially masking more problems in the simulated cloud climatology. A closer look at the models’ capability in reproducing observed cloud ice and the contribution of biases in cloud liquid and cloud ice to deficiencies in the simulated climatology of cloud amount and cloud radiative forcing is the logical next step in understanding which cloud processes need particular attention for future model improvements. Acknowledgments. This research was supported by the Japan Agency for Marine-Earth Science and Technology (JAMSTEC), by NASA through Grant NNX07AG53G, and by NOAA through Grant NA09OAR4320075, which sponsored research at the International Pacific Research Center. We acknowledge the World Climate Research Programme’s Working Group on Coupled Modelling, which is responsible for CMIP, and we thank the climate modeling groups (listed in Tables 1–3 of this paper) for producing and making available their model output. For CMIP, the U.S. Department of Energy’s Program for Climate Model Diagnosis and Intercomparison provides coordinating support and led development of software infrastructure in partnership with the Global Organization for Earth System Science Portals. We thank the journal reviewers for valuable comments on the manuscript. REFERENCES Alekseev, V. A., E. M. Volodin, V. Ya. Galin, V. P. Dymnikov, and V. N. Lykossov, 1998: Modelling of the present-day climate by the atmospheric model of INM RAS ‘‘DNM GCM.’’ INM Tech. Rep. N2086-B98, Institute of Numerical Mathematics, Russian Academy of Sciences, 208 pp. Arora, V. K., and Coauthors, 2011: Carbon emission limits required to satisfy future representative concentration pathways of

3843

greenhouse gases. Geophys. Res. Lett., 38, L05805, doi:10.1029/ 2010GL046270. Barkstrom, B. R., 1984: The Earth Radiation Budget Experiment (ERBE). Bull. Amer. Meteor. Soc., 65, 1170–1185. Bleck, R., 2002: An oceanic general circulation model framed in hybrid isopycnic-Cartesian coordinates. Ocean Modell., 4, 55–88. Bodas-Salcedo, A., and Coauthors, 2011: COSP: Satellite simulation software for model assessment. Bull. Amer. Meteor. Soc., 92, 1023–1043. Bony, S., 1997: Sea surface temperature and large-scale circulation influences on tropical greenhouse effect and cloud radiative forcing. J. Climate, 10, 2055–2077. ——, and J.-L. Dufresne, 2005: Marine boundary layer clouds at the heart of tropical cloud feedback uncertainties in climate models. Geophys. Res. Lett., 32, L20806, doi:10.1029/ 2005GL023851. ——, M. Webb, C. Bretherton, S. A. Klein, P. Siebesma, G. Tselioudis, and M. Zhang, 2011: CFMIP: Towards a better evaluation and understanding of clouds and cloud feedbacks in CMIP5 models. CLIVAR Exchanges, Vol. 2, No. 56, International CLIVAR Project Office, Southampton, United Kingdom, 20–24. Cess, R. D., and Coauthors, 1990: Intercomparison and interpretation of climate feedback processes in 19 atmospheric general circulation models. J. Geophys. Res., 95 (D10), 16 601–16 615. Chen, W.-T., C. P. Woods, J.-L. F. Li, D. E. Waliser, J.-D. Chern, W.-K. Tao, J. H. Jiang, and A. M. Tompkins, 2011: Partitioning CloudSat ice water content for comparison with upper tropospheric ice in global atmospheric models. J. Geophys. Res., 116, D19206, doi:10.1029/2010JD015179. Collins, M., S. F. B. Tett, and C. Cooper, 2001: The internal climate variability of HadCM3, a version of the Hadley Centre Coupled Model without flux adjustments. Climate Dyn., 17, 61–81. Collins, W. D., and Coauthors, 2006: The Community Climate System Model: CCSM3. J. Climate, 19, 2122–2143. Collins, W. J., and Coauthors, 2008: Evaluation of the HadGEM2 model. Hadley Center Tech. Note 74, 47 pp. [Available online at http://www.metoffice.gov.uk/media/pdf/8/7/HCTN_74.pdf.] ——, and Coauthors, 2011: Development and evaluation of an earth-system model—HadGEM2. Geosci. Model Dev., 4, 1051–1075, doi:10.5194/gmd-4-1051-2011. Donner, L. J., and Coauthors, 2011: The dynamical core, physical parameterizations, and basic simulation characteristics of the atmospheric component AM3 of the GFDL global coupled model CM3. J. Climate, 24, 3483–3519. Dufresne J.-L., and Coauthors, 2013: Climate change projections using the IPSL-CM5 earth system model: From CMIP3 to CMIP5. Climate Dyn., 40, 2123–2165, doi:10.1007/s00382-012-1636-1. Flato, G. M., 2005: The third generation Coupled Global Climate Model. [Available online at http://www.ec.gc.ca/ccmac-cccma/ default.asp?lang5En&n51299529F-1.] Galin, V. Ya., E. M. Volodin, and S. P. Smyshliaev, 2003: Atmosphere general circulation model of INM RAS with ozone dynamics. Russ. Meteor. Hydrol., 5, 13–22. Gent, P. R., and Coauthors, 2011: The Community Climate System Model version 4. J. Climate, 24, 4973–4991. GFDL Global Atmospheric Model Development Team, 2004: The new GFDL global atmosphere and land model AM2–LM2: Evaluation with prescribed SST simulations. J. Climate, 17, 4641–4673. Gordon, C., C. Cooper, C. A. Senior, H. T. Banks, J. M. Gregory, T. C. Johns, J. F. B. Mitchell, and R. A. Wood, 2000: The simulation of SST, sea ice extents and ocean heat transports in

3844

JOURNAL OF CLIMATE

a version of the Hadley Centre coupled model without flux adjustments. Climate Dyn., 16, 147–168. Gordon, H. B., and Coauthors, 2002: The CSIRO Mk3 Climate System Model. CSIRO Atmospheric Research Tech. Paper 60, 130 pp. [Available online at http://www.cmar.csiro.au/ e-print/open/gordon_2002a.pdf.] Hartmann, D. L., and D. Doelling, 1991: On the net radiative effectiveness of clouds. J. Geophys. Res., 96 (D1), 869–891. Hasumi, H., and S. Emori, Eds., 2004: K-l coupled GCM (MIROC) description. K-l Tech. Rep. 1, Center for Climate System Research, University of Tokyo, Tokyo, Japan, 34 pp. [Available online at http://www.ccsr.u-tokyo.ac.jp/kyosei/hasumi/ MIROC/tech-repo.pdf.] Hourdin, F., and Coauthors, 2006: The LMDZ4 general circulation model: Climate performance and sensitivity to parameterized physics with emphasis on tropical convection. Climate Dyn., 27, 787–813. ——, and Coauthors, 2013: Impact of the LMDZ atmospheric grid configuration on the climate and sensitivity of the IPSL-CM5A coupled model. Climate Dyn., 40, 2167–2192, doi:10.1007/s00382-012-1411-3, in press. Jiang, J. H., and Coauthors, 2012: Evaluation of cloud and water vapor simulations in CMIP5 climate models using NASA ‘‘A-Train’’ satellite observations. J. Geophys. Res., 117, D14105, doi:10.1029/2011JD017237. Kay, J. E., and Coauthors, 2012: Exposing global cloud biases in the Community Atmosphere Model (CAM) using satellite observations and their corresponding instrument simulators. J. Climate, 25, 5190–5207. Kirkeva˚g, A., and Coauthors, 2013: Aerosol–climate interactions in the Norwegian earth system model—NorESM1-M. Geosci. Model Dev., 6, 207–244, doi:10.5194/gmd-6-207-2013. Lauer, A., K. Hamilton, Y. Wang, V. T. J. Phillips, and R. Bennartz, 2010: The impact of global warming on marine boundary layer clouds over the eastern Pacific—A regional model study. J. Climate, 23, 5844–5863. Li, L., and Coauthors, 2010: Development and evaluation of GAMIL2.0 and FGOALS2.0-g. Preprints, Fifth C20C Workshop, Beijing, China, IAP, Chinese Academy of Sciences, 29 pp. [Available online at http://www.lasg.ac.cn/C20C/UserFiles/ File/c20c-ljli.pdf.] Lin, J.-L., 2007: The double-ITCZ problem in IPCC AR4 coupled GCMs: Ocean–atmosphere feedback analysis. J. Climate, 20, 4497–4525. Martin, G. M., M. A. Ringer, V. D. Pope, A. Jones, C. Dearden, and T. J. Hinton, 2006: The physical properties of the atmosphere in the new Hadley Centre Global Environmental Model, HadGEMl. Part I: Model description and global climatology. J. Climate, 19, 1274–1301. McFarlane, N. A., G. J. Boer, J.-P. Blanchet, and M. Lazare, 1992: The Canadian Climate Centre second-generation general circulation model and its equilibrium climate. J. Climate, 5, 1013–1044. Medeiros, B., B. Stevens, I. M. Held, M. Zhao, D. L. Williamson, J. G. Olson, and C. S. Bretherton, 2008: Aquaplanets, climate sensitivity, and low clouds. J. Climate, 21, 4974–4991. Meehl, G. A., C. Covery, T. Delworth, M. Latif, B. McAvaney, J. F. B. Mitchell, R. J. Stouffer, and K. E. Taylor, 2007: The WCRP CMIP3 multimodel dataset: A new era in climate change research. Bull. Amer. Meteor. Soc., 88, 1383–1394. Mizuta, R., and Coauthors, 2012: Climate simulations using MRIAGCM3.2 with 20-km grid. J. Meteor. Soc. Japan, 90, 233–258. Nam, C. C. W., and J. Quaas, 2012: Evaluation of clouds and precipitation in the ECHAM5 general circulation model using

VOLUME 26

CALIPSO and CloudSat satellite data. J. Climate, 25, 4975– 4992. O’Dell, C. W., F. J. Wentz, and R. Bennartz, 2008: Cloud liquid water path from satellite-based passive microwave observations: A new climatology over the global oceans. J. Climate, 21, 1721–1739. Park, S., and C. B. Leovy, 2004: Marine low-cloud anomalies associated with ENSO. J. Climate, 17, 3448–3469. Pope, V., M. L. Gallani, P. R. Rowntree, and R. A. Stratton, 2000: The impact of new physical parameterizations in the Hadley Centre climate model: HadAM3. Climate Dyn., 16, 123–146. Reynolds, R. W., T. M. Smith, C. Liu, D. B. Chelton, K. S. Casey, and M. S. Schlax, 2007: Daily high-resolution-blended analyses for sea surface temperature. J. Climate, 20, 5473–5496. Roeckner, E., and Coauthors, 1996: The atmospheric general circulation model ECHAM4: Model description and simulation of present-day climate. Rep. 218, Max Planck Institute for Meteorology, Hamburg, Germany, 90 pp. [Available online at http://www.mpimet.mpg.de/fileadmin/publikationen/Reports/ MPI-Report_218.pdf.] ——, and Coauthors, 2003: The atmospheric general circulation model ECHAM5. Part I: Model description. Rep. 349, Max Planck Institute, Hamburg, Germany, 127 pp. [Available online at http://www.mpimet.mpg.de/fileadmin/publikationen/ Reports/max_scirep_349.pdf.] ——, and Coauthors, 2006: Sensitivity of simulated climate to horizontal and vertical resolution in the ECHAM5 atmosphere model. J. Climate, 19, 3771–3791. Rossow, W. B., and L. C. Garder, 1993: Validation of ISCCP cloud detections. J. Climate, 6, 2370–2393. ——, and R. A. Schiffer, 1999: Advances in understanding clouds from ISCCP. Bull. Amer. Meteor. Soc., 80, 2261–2287. ——, A. W. Walker, D. E. Beuschel, and M. D. Roiter, 1996: International Cloud Climatology Project (ISCCP) documentation of new cloud datasets. WMO/TD-737, 115 pp. Rotstayn, L. D., M. A. Collier, M. R. Dix, Y. Feng, H. B. Gordon, S. P. O’Farrell, I. N. Smith, and J. Syktus, 2010: Improved simulation of Australian climate and ENSO-related climate variability in a GCM with an interactive aerosol treatment. Int. J. Climatol., 30, 1067–1088, doi:10.1002/joc.1952. Russell, G. L., J. R. Miller, and D. Rind, 1995: A coupled atmosphere–ocean model for transient climate change studies. Atmos.–Ocean, 33, 683–730. Sakamoto, T. T., and Coauthors, 2012: MIROC4h—A new highresolution atmosphere–ocean coupled general circulation model. J. Meteor. Soc. Japan, 90, 325–359. Salas-Me´lia, D., and Coauthors, 2005: Description and validation of the CNRM-CM3 global coupled model. CNRM Working Note 103, 36 pp. [Available online at http://www.cnrm.meteo. fr/scenario2004/paper_cm3.pdf.] Schmidt, G. A., and Coauthors, 2006: Present-day atmospheric simulations using GISS ModelE: Comparison to in situ, satellite, and reanalysis data. J. Climate, 19, 153–192. Solomon, S., D. Qin, M. Manning, M. Marquis, K. Averyt, M. M. B. Tignor, H. L. Miller Jr., and Z. Chen, Eds., 2007: Climate Change 2007: The Physical Science Basis. Cambridge University Press, 996 pp. Stowasser, M., K. Hamilton, and G. J. Boer, 2006: Local and global climate feedbacks in models with differing climate sensitivities. J. Climate, 19, 193–209. Taylor, K. E., 2001: Summarizing multiple aspects of model performance in a single diagram. J. Geophys. Res., 106 (D7), 7183– 7192.

1 JUNE 2013

LAUER AND HAMILTON

——, R. Stouffer, and G. A. Meehl, 2012: An overview of CMIP5 and the experiment design. Bull. Amer. Meteor. Soc., 93, 485–498. Voldoire, A., and Coauthors, 2013: The CNRM-CM5.1 global climate model: Description and basic evaluation. Climate Dyn., 40, 2091–2121, doi:10.1007/s00382-011-1259-y. Volodin, E. M., N. A. Dianskii, and A. V. Gusev, 2010: Simulating present-day climate with the INMCM4.0 coupled model of the atmospheric and oceanic general circulations. Atmos. Oceanic Phys., 46, 414–431, doi:10.1134/S000143381004002X. Waliser, D. E., K.-W. Seo, S. Schubert, and E. Njoku, 2007: Global water cycle agreement in the climate models assessed in the IPCC AR4. Geophys. Res. Lett., 34, L16705, doi:10.1029/ 2007GL030675. ——, and Coauthors, 2009: Cloud ice: A climate model challenge with signs and expectations of progress. J. Geophys. Res., 114, D00A21, doi:10.1029/2008JD010015. Washington, W. M., and Coauthors, 2000: Parallel climate model (PCM) control and transient simulations. Climate Dyn., 16, 755–774. Watanabe, M., and Coauthors, 2010: Improved climate simulation by MIROC5: Mean states, variability, and climate sensitivity. J. Climate, 23, 6312–6335. Watanabe, S., and Coauthors, 2011: MIROC-ESM 2010: Model description and basic results of CMIP5-20c3m experiments. Geosci. Model Dev., 4, 845–872. Weare, B. C., 2004: A comparison of AMIP II model cloud layer properties with ISCCP D2 estimates. Climate Dyn., 22, 281–292. Wetzel, P., E. Maier-Reimer, M. Botzet, J. Jungclaus, N. Keenlyside, and M. Latif, 2006: Effects of ocean biology on the penetrative radiation in a coupled climate model. J. Climate, 19, 3973–3987.

3845

Wielicki, B. A., B. R. Barkstrom, E. F. Harrison, R. B. Lee III, G. L. Smith, and J. E. Cooper, 1996: Clouds and the Earth’s Radiant Energy System (CERES): An Earth Observing System experiment. Bull. Amer. Meteor. Soc., 77, 853–868. Wu, T., and Coauthors, 2010: The Beijing Climate Center for Atmospheric General Circulation Model (BCC-AGCM2.0.1): Description and its performance for the present-day climate. Climate Dyn., 34, 123–147. Yu, Y., X. Zhang, and Y. Guo, 2004: Global coupled ocean– atmosphere general circulation models in LASG/IAP. Adv. Atmos. Sci., 21, 444–455. Yukimoto, S., and Coauthors, 2006: Present-day climate and climate sensitivity in the Meteorological Research Institute Coupled GCM version 2.3 (MRI-CGCM2.3). J. Meteor. Soc. Japan, 84, 333–363. ——, and Coauthors, 2011: Meteorological Research Institute– Earth System Model version 1 (MRI-ESM1): Model description. Tech. Rep. 64, Meteorological Research Institute, Japan Meteorological Agency, Tsukuba, Japan, 83 pp. Zhang, M. H., and Coauthors, 2005: Comparing clouds and their seasonal variations in 10 atmospheric general circulation models with satellite measurements. J. Geophys. Res., 110, D15S02, doi:10.1029/2004JD005021. Zhang, Y., W. B. Rossow, A. A. Lacis, V. Oinas, and M. I. Mishchenko, 2004: Calculation of radiative fluxes from the surface to top of atmosphere based on ISCCP and other global data sets: Refinements of the radiative transfer model and the input data. J. Geophys. Res., 109, D19105, doi:10.1029/2003JD004457. Zhao, M., I. M. Held, S.-J. Lin, and G. A. Vecchi, 2009: Simulations of global hurricane climatology, interannual variability, and response to global warming using a 50-km resolution GCM. J. Climate, 22, 6653–6678.