Journal cover Journal topic
Atmospheric Chemistry and Physics An interactive open-access journal of the European Geosciences Union
Journal topic
Atmos. Chem. Phys., 18, 8807-8828, 2018
https://doi.org/10.5194/acp-18-8807-2018
Atmos. Chem. Phys., 18, 8807-8828, 2018
https://doi.org/10.5194/acp-18-8807-2018

Research article 22 Jun 2018

Research article | 22 Jun 2018

# The importance of mixed-phase and ice clouds for climate sensitivity in the global aerosol–climate model ECHAM6-HAM2

Ice containing clouds and climate sensitivity
Ulrike Lohmann and David Neubauer Ulrike Lohmann and David Neubauer
• Institute for Atmospheric and Climate Science, ETH Zurich, Zurich, Switzerland
Abstract

How clouds change in a warmer climate remains one of the largest uncertainties for the equilibrium climate sensitivity (ECS). While a large spread in the cloud feedback arises from low-level clouds, it was recently shown that mixed-phase clouds are also important for ECS. If mixed-phase clouds in the current climate contain too few supercooled cloud droplets, too much ice will change to liquid water in a warmer climate. As shown by , this overestimates the negative cloud-phase feedback and underestimates ECS in the CAM global climate model (GCM). Here we use the newest version of the ECHAM6-HAM2 GCM to investigate the importance of mixed-phase and ice clouds for ECS.

Although we also considerably underestimate the fraction of supercooled liquid water globally in the reference version of the ECHAM6-HAM2 GCM, we do not obtain increases in ECS in simulations with more supercooled liquid water in the present-day climate, different from the findings by . We hypothesize that it is not the global supercooled liquid water fraction that matters, but only how well low- and mid-level mixed-phase clouds with cloud-top temperatures in the mixed-phase temperature range between 0 and 35 C that are not shielded by higher-lying ice clouds are simulated. These occur most frequently in midlatitudes, in particular over the Southern Ocean where they determine the amount of absorbed shortwave radiation. In ECHAM6-HAM2 the amount of absorbed shortwave radiation over the Southern Ocean is only significantly overestimated if all clouds below 0 C consist exclusively of ice. Only in this simulation is ECS significantly smaller than in all other simulations and the cloud optical depth feedback is the dominant cloud feedback. In all other simulations, the cloud optical depth feedback is weak and changes in cloud feedbacks associated with cloud amount and cloud-top pressure dominate the overall cloud feedback. However, apart from the simulation with only ice below 0 C, differences in the overall cloud feedback are not translated into differences in ECS in our model. This insensitivity to the cloud feedback in our model is explained with compensating effects in the clear sky.

1 Introduction

Changes in clouds remain one of the largest uncertainties for the calculation of the response of the climate system to a given radiative forcing ΔF . This response can be described with the following equation:

$\begin{array}{}\text{(1)}& \mathrm{\Delta }F=\mathrm{\Delta }R+\mathrm{\Delta }H+\mathit{\lambda }\mathrm{\Delta }{T}_{\mathrm{s}}.\end{array}$

Here ΔR represents the net radiative imbalance at the top of the atmosphere (TOA), ΔH is the heat taken up by the ocean, ΔTs is the net change in global annual mean surface temperature, and λ is the climate feedback parameter in W m−2 K−1. ΔR and ΔTs change over time. In addition, some studies suggest that the heat uptake by the ocean is also time dependent and even λ is not a constant parameter (e.g., Knutti and Rugenstein2015).

Here we evaluate the increase in ΔTs, which results from a doubling of carbon dioxide (CO2) with respect to preindustrial concentrations, i.e., from a $\mathrm{\Delta }{F}_{\mathrm{2}×{\mathrm{CO}}_{\mathrm{2}}}$ on the order of 3.7 W m−2 . There are two metrics that describe the temperature response to a doubling of CO2, the transient climate response (TCR) and the equilibrium climate sensitivity (ECS). TCR is estimated at the time of CO2 doubling from atmosphere global climate model (GCM) simulations in which the CO2 concentration increases by 1 % yr−1 and that are coupled to a dynamic ocean (e.g., Flato et al.2013). ECS is obtained from coupled atmosphere–ocean simulations that are abruptly exposed to a CO2 doubling relative to preindustrial concentrations and then run until a new equilibrium of the climate system has been established. This requires coupling of the atmosphere GCM to a fully coupled dynamic ocean model (e.g., Gregory et al.2004; Flato et al.2013). As such a coupled system takes a long time to reach a new equilibrium, sometimes ECS is approximated from coupled atmosphere–mixed-layer ocean (MLO) climate models (e.g., Meehl et al.2007; Randall et al.2007).

ECS can also be determined from the so-called Gregory method , in which the TOA radiative flux is regressed against the annual global averaged surface air temperature change. While TCR is more relevant to present-day climate change because it is obtained from transient simulations that at the time of CO2 doubling have not reached a new equilibrium climate, it is more expensive to calculate because it requires fully coupled atmosphere–ocean model simulations over the whole period. Therefore we focus on ECS from coupled atmosphere–MLO simulations in this paper. At the time of the new equilibrium in the, ΔR and ΔH vanish and Eq. (1) reduces to

$\begin{array}{}\text{(2)}& \mathrm{\Delta }{F}_{\mathrm{2}×{\mathrm{CO}}_{\mathrm{2}}}=\mathit{\lambda }\mathrm{\Delta }{T}_{\mathrm{s}},\end{array}$

and ΔTs equals the ECS. λ can be decomposed into different feedbacks (e.g., Soden et al.2008; Shell et al.2008). λ is the sum of the Planck feedback, the water vapor feedback, the lapse rate feedback, the surface albedo feedback, the cloud feedback λc, that will be discussed further below, and a residual term . The residual term considers the interactions among the different components, nonlinear dependencies of ECS from the forcing and errors in the radiative kernel method . This term needs to be small for the decomposition into individual feedbacks to explain the vast majority of the feedback. The average λ value from the GCMs that participated in the Coupled Model Intercomparison Project Phase 5 (CMIP5) is 1.1 W m−2 K−1 with a 90 % uncertainty range of ±0.5 W m−2 K−1 . The uncertainty in λ induces the uncertainty in ECS. Neither the average ECS of approximately 3 K nor its uncertainty has changed much over time.

In a first model intercomparison paper by , the estimates of ECS in GCMs varied between 1.4 and 4.1 K. Estimates from the GCMs used in CMIP5 are somewhat higher with 2.1 and 4.7 K , but the range in estimates remains similar. One of the prime contributors to this range in uncertainty are inter-model differences in the low-level shortwave cloud feedback . While the overall cloud feedback is positive with 0.3 W m−2 K−1, the spread in λc is larger than that of the overall climate feedback parameter and varies by ±0.7 W m−2 K−1 between the CMIP5 models . This means a negative λc cannot be ruled out and the careful conclusion from the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (IPCC AR5) was that λc is likely (with a probability of 66 %) positive. Contributors to the positive net cloud feedback are a poleward shift of the midlatitude storms and decreases in the coverage of low-level clouds, both of which result in less scattering of solar radiation . Another positive cloud feedback that operates mainly in the longwave radiation regime is the increase in the height of deep convective outflow. Because the troposphere warms but the cloud-top temperature remains at roughly the same temperature (fixed anvil temperature; ), the rate of longwave emission to space from these cloud tops does not keep pace with that of the underlying atmosphere, causing more longwave radiation to remain in the Earth–atmosphere system.

One of the more uncertain cloud feedbacks is the negative cloud optical depth (COD) feedback . The COD (τ) is proportional to the cloud water path over the effective cloud particle size. Therefore changes in τ can occur due to changes in the cloud water path, which in turn originate from changes in the hydrological cycle, changes in the size of the cloud droplets and ice crystals and changes in cloud phase from ice to water or vice versa. The negative COD feedback is related to a shift from cloud ice in the present-day climate to cloud liquid water at the same altitude in the warmer climate. Because cloud droplets are generally smaller and more numerous than ice crystals and they have a different refractive index, the optical depth of liquid clouds is larger than that of ice clouds for a constant cloud water path. Precipitation formation is also less efficient for liquid clouds than for ice clouds (e.g., Lohmann1996; Hoose et al.2008), further increasing cloud lifetime and τ of liquid clouds. If, in addition, the cloud water path increases in a warmer climate, τ will be further enlarged. All of these aspects result in a negative COD feedback.

analyzed the ratio of supercooled liquid water to the sum of cloud liquid water and cloud ice (supercooled liquid fraction, SLF) between 0 and 35 C in the CAM5 GCM. They found SLF to be systematically underestimated with respect to CALIOP observations, i.e., too much condensate to be in the form of ice at these temperatures. This led to a too negative cloud-phase feedback and in the CAM5 GCM to a too low ECS. When constraining the present-day SLF by satellite observations, their ECS increased by up to 1.3 C. A similar increase in ECS of 1.5 C was found by when increasing the fraction of supercooled liquid clouds over the Southern Ocean. analyzed the low-COD feedback from eight CMIP5 models with that inferred from ISCCP, MODIS, and PATMOS satellite data. They also concluded that the low-COD feedback is likely too negative at mid and high latitudes in climate models. Motivated by these studies and to understand how universal the findings of the change in ECS using extreme assumptions about the liquid water/ice phase in mixed-phase clouds of are, here we use the ECHAM6-HAM2 GCM to calculate the impact of mixed-phase and ice clouds on ECS.

2 Description of the model and the sensitivity studies

In this paper we use the latest version of the aerosol–climate model ECHAM6.3-HAM2.3, which assembles the most recent versions of the atmospheric general circulation model ECHAM (namely ECHAM6, as described in with the most recent changes described below) and the aerosol module HAM based on the HAM2 version as described in , , and and summarized below. As our simulations are conducted with monthly mean oxidant fields instead of with interactive chemistry as in HAMMOZ , we refer to it as ECHAM6-HAM2.

## 2.1 ECHAM6

ECHAM6 is the latest generation of the atmospheric general circulation model developed by the Max Planck Institute for Meteorology . Like its predecessors, ECHAM6 employs a spectral transform dynamical core and a flux-form semi-Lagrangian tracer transport algorithm from . Vertical mixing occurs through turbulent mixing, moist convection (including shallow, deep, and mid-level convection), and momentum transport by gravity waves arising from boundary effects or atmospheric disturbances. Sub-grid-scale cloudiness (stratiform clouds) is represented using the scheme of , which diagnostically calculates the grid cell cloud fraction as a function of the relative humidity in a given grid cell, once a threshold value is exceeded. Cloud liquid water and cloud ice mixing ratios are treated prognostically following . In ECHAM6-HAM2, additional prognostic equations for the number concentrations of cloud droplets and ice crystals are included . The two-moment cloud microphysics scheme and the aerosol scheme are described in Sect. 2.2. Radiative transfer in ECHAM6 is represented using the radiation transfer broadband model PSrad , which considers 14 and 16 bands for the shortwave (820 to 50 000 cm−1) and longwave (10 to 3000 cm−1) parts of the spectrum, respectively .

Radiative transfer is computed based on the amount of gases, aerosols, and clouds in the atmosphere and their related optical properties. Trace gas concentrations of long-lived greenhouse gases are specified in the model. Optical properties of aerosol particles and clouds are precalculated for each band of the RRTMG scheme using Mie theory offline and they read from a look-up table based on the concentration of cloud droplets and ice crystals as computed by the two-moment scheme.

ECHAM6 includes a new land-surface model JSBACH . JSBACH assumes that each grid box is composed of two fractions, one representing bare soil and the other being covered with vegetation, this one being further subdivided into tiles, one for each of the 11 plant functional types distinguished in JSBACH. Soil hydrology is represented with a single-layer bucket model. A new treatment of the surface albedo is included, which accounts for the different sections of bare and vegetated areas.

In addition to the improved representation of solar radiative transfer by the RRTMG scheme and the improved surface albedo, smaller changes are included. The vertical discretization within the troposphere (in particular in the upper troposphere and lower stratosphere) is slightly different, the representation of convective triggering has been improved, and the tuning of various model parameters was adjusted. In contrast to ECHAM5, ECHAM6 is more commonly used in a middle-atmosphere configuration, i.e., with the two verticals grids L47 and L95 that resolve the atmosphere from the surface up to 0.01 hPa (roughly 80 km).

Changes from ECHAM6.1, as coupled to HAM2 in , to ECHAM6.3 used here include small changes in convection, in the middle atmosphere, improvements in snow and sea ice coverage, an improved aerosol climatology combined with a simple plume implementation of anthropogenic aerosols , and an improved submodel interface. ECHAM6.3 uses the Monte Carlo independent column approximation radiation scheme with the option for spectral sampling in time . The land model JSBACH has been updated with an improved hydrology and soil model. Changes in the cloud cover scheme were made to improve cloud cover in stratocumulus regions. ECHAM6.1 suffered from artificial cloud blinking that occurred due to a bug in the cloud cover scheme, causing the (fractional) cloud cover to be either 0 or 1. This cloud blinking has been removed in ECHAM6.3.

As in previous versions of ECHAM-HAM, ECHAM drives the aerosol and chemistry modules through the generic sub-model interface by providing meteorological conditions such as wind, temperature, pressure, specific humidity, and conditions pertaining to the land surface (taken from JSBACH) such as leaf area index. Aerosol particles and their precursors are transported in the same way as water vapor and cloud species.

## 2.2 Aerosol and cloud microphysics scheme

The aerosol module HAM predicts the evolution of an aerosol ensemble considering five components: sulfate, black carbon, particulate organic matter, sea salt, and mineral dust. The aerosol spectrum is described by the superposition of seven log-normal modes ranging from 0.005 to > 0.5 µm. Aerosol particles within a mode are assumed to be internally mixed, in the sense that each particle can consist of multiple components. Aerosols of different modes are externally mixed, meaning that they coexist in the atmosphere as independent particles. The seven modes of the aerosol number distribution are grouped into four geometrical size classes, including the nucleation, Aitken, accumulation, and coarse modes. Particles in four of the modes contain at least one soluble compound; thus they can take up water and are referred to as soluble. The particles in the other three modes consist of compounds with no or low water solubility and are referred to as insoluble. Through aging processes, insoluble particles can become soluble. Each mode of the aerosol size number distribution is described by the three moments, the aerosol number N, the number median radius r, and the standard deviation σ. The standard deviation is assumed to be constant and is set to 1.59 for the nucleation, Aitken, and accumulation modes and to 2 for the coarse modes. The median radius of each mode is calculated from the aerosol number and aerosol mass, which are transported as tracers. For more details, please refer to the description of HAM in and .

The two-moment cloud microphysics scheme that predicts the number concentrations and mass mixing ratios of cloud droplets and ice crystals in stratiform clouds as implemented in ECHAM5 is described in and . The microphysics scheme includes all phase changes between the water components (condensation, evaporation, freezing, melting, deposition, and sublimation) and precipitation processes (autoconversion of cloud droplets, accretion of raindrops with cloud droplets and snow flakes with cloud droplets and ice crystals, aggregation of ice crystals). Moreover, evaporation of raindrops and melting and sublimation of snowflakes are considered, as well as sedimentation of cloud ice. The cloud microphysics scheme is coupled to the aerosol scheme HAM through the processes of cloud droplet activation and ice crystal nucleation as well as through in-cloud and below-cloud scavenging . Convective clouds are not radiatively active and have much simpler conversion rates (Tiedtke1989). The detrained condensate from convective clouds is a source for the stratiform cloud scheme.

Heterogeneous freezing in mixed-phase clouds occurs by contact and immersion freezing as discussed below. In the standard configuration, cirrus clouds are assumed to form by homogeneous freezing of supercooled solution droplets . Homogeneous and heterogeneous nucleation strongly depends on the vertical velocity that results in supersaturation. The vertical velocity is a superposition of the grid-scale vertical velocity and a sub-grid-scale component, which is linked to the turbulent kinetic energy .

Since the ECHAM-HAM model validation in , the cloud microphysics scheme has undergone the following scientific improvements.

1. While in the previous ECHAM6-HAM version cloud droplet activation followed the empirical scheme by , it now uses a parameterization of cloud droplet activation based on Köhler theory by as discussed in .

2. Previously, the increase in cloud droplet number concentration (CDNC) from cloud droplets that were detrained from convective clouds, was applied everywhere in the grid box. We now weight the detrained CDNC by the detrained mass and add it to the mass-weighted CDNC of the stratiform part. In addition, the split of the detrained cloud water mass into liquid water and ice was made consistent between the number concentrations and mass mixing ratios of cloud droplets and ice crystals.

3. We now assume the shape of ice crystals to be hexagonal plates of crystal type P1a following . Before we used the empirical mass-size relationships from their Table 2.4, which are only valid for crystals between 0.3 and 1.5 mm in size. We replaced that with their empirical relationship in Table 2.2, which is valid from 10 µm to 3 mm and thus covers the whole ice crystal size range. Also, while we assumed that all ice crystals are hexagonal plates in the cloud microphysics scheme, cirrus crystals were treated as spheres for the calculation of the effective ice crystal radius and ice cloud optical properties in . This inconsistency has been eliminated and ice crystals are now hexagonal plates everywhere in the module.

4. The heterogeneous nucleation scheme in the mixed-phase cloud regime among temperatures between 273.15 and 238.15 K considers contact nucleation by mineral dust and immersion nucleation by black carbon and mineral dust following . However, because there is no evidence that Aitken mode particles contribute to freezing , we now limit the immersion freezing of black carbon to particles in the accumulation mode or larger.

5. The temperature dependence of sticking efficiency used for accretion of ice crystals by snow has been changed to the expression used in .

6. In previous versions, the minimum CDNC (CDNCmin) was set to 40 cm−3. The justification for this choice was twofold: first, while observations of clouds with a lower CDNC exist (e.g, Terai et al.2014; Wood et al.2018), these smaller concentrations normally occur in clouds or pockets in clouds that are much smaller than our grid boxes. Second, so far we do not account for nitrate aerosols and our treatment of secondary organic aerosols is rather simplistic and likely underestimates the organic aerosol concentration . Therefore we are likely to underestimate CDNC, which we partly buffer by using CDNCmin = 40 cm−3. However, CDNCmin has a large impact on the aerosol radiative forcing . Therefore we introduced the option to have two ECHAM6.3-HAM2.3 versions, one that keeps the CDNCmin at 40 cm−3 and one in which we lower CDNCmin to 10 cm−3.

In addition to the scientific improvements, we removed smaller inconsistencies, such that cloud droplets and ice crystals could both grow and evaporate/sublimate in one time step and that there is nonzero CDNC below 238.15 K, nonzero ice crystal number concentrations above 273.15 K, and inconsistencies in the calculation of cloud cover, condensation, and the ice crystal number concentration in cirrus clouds. Moreover, the two-moment cloud microphysics scheme is now energy conserving and has been modularized.

3 Model set-up and experiments

In this paper, we compare results from the release version of ECHAM6-HAM2 with CDNCmin = 40 cm−3 (simulation REF) and the sensitivity simulations discussed below with observations. To address the impact of SLF on ECS in ECHAM6-HAM2, we conducted one simulation with no supercooled liquid water below 0 C (simulation ALL_ICE) and one simulation with no ice formation at temperatures >35 C (simulation ALL_LIQ), similar to what has been done in and . The simulation ALL_ICE is set up such that all cloud droplets that are advected to colder temperatures are forced to freeze instantaneously and all the detrained cloud condensate is in the form of ice at temperatures < 0 C. The simulation ALL_LIQ is set up such that heterogeneous freezing is turned off in the temperature range between 0 and 35 C and detrainment of ice crystals from convective clouds is restricted to temperatures below 35 C. Ice crystals sedimenting from colder into warmer cloud layers melt at 35 C. The number of cloud droplets which freeze at temperatures below 35 C had to be reduced by a factor of 100 in order to keep the ice crystal number concentration realistic in simulation ALL_LIQ and to reduce tuning to a minimum.

In this paper we go beyond the study by and also test the impact of other aspects of clouds on ECS, such as the type of nucleation in cirrus. There are discussions if cirrus clouds are mainly formed by homogeneous nucleation as we assume in our reference simulation or are nucleated heterogeneously . We investigate the impact of homogeneous vs. heterogeneous freezing in cirrus clouds by performing one simulation in which all cirrus clouds form heterogeneously (simulation HET) instead of homogeneous nucleation in simulation REF. Heterogeneously nucleated cirrus clouds are optically thinner (Lohmann2008). Thus, with this simulation we aim to address the impact of cirrus on ECS that so far remains uncertain .

Last but not least we investigate the impact of parameterized convection on ECS by performing a simulation in which we completely switched off convection (simulation NOCONV) as in . This is normally performed only in simulations run at horizontal resolutions of less than 10 km, at which the vertical motions associated with deep convection start to be resolved. While our horizontal resolution is much coarser and thus a convective parameterization is needed, there are some inconsistencies in terms of microphysics between convective and stratiform clouds. Therefore we evaluate how the cloud fields, the climate in general, and ECS are simulated if only large-scale clouds are allowed to form.

All simulations were performed in T63 spectral resolution, which corresponds to 1.875× 1.875 and 31 vertical layers with a top at 10 hPa. The present-day atmosphere-only simulations were run over 25 years, after a 3-month spin-up with fixed sea surface temperatures (SSTs) and sea-ice cover (averaged over the atmosphere model intercomparison project (AMIP) climatology for the years 2000–2015). To calculate ECS, ECHAM6-HAM2 has been coupled to a MLO of 50 m depth. The deep ocean heat flux is computed from the atmosphere-only simulations. Because the deep ocean heat flux adjusts to the different set-ups, we computed it individually for each of the five different set-ups used in this study. These simulations were spun up for 25 years after which the simulations were in equilibrium. We then ran them for another 25 years, over which the results were averaged.

The sensitivity studies conducted for our studies are summarized in Table 1. For the calculation of ECS, all simulations need to be in radiative equilibrium at the TOA. This might require retuning. To keep the different simulations as comparable as possible, we only adjusted two parameters inside the two-moment cloud microphysics scheme. γr speeds up the autoconversion rate of cloud droplets to grow to raindrops by collision–coalescence and accounts for the missing sub-grid-scale variability in cloud water and CDNC (Wood2002). γs is the corresponding process in the ice phase, which enhances the aggregation of ice crystals to form snow flakes . The values of these two parameters are also included in Table 1. In addition to the changes in the tuning parameters, we needed to decrease the time step from 7.5 to 5 min in simulation NOCONV because of numerical stability.

Table 1Set-up of the simulations together with the two tuning parameters that differ among the simulations.

All of the simulations listed in Table 1 were run in three different set-ups: (i) atmosphere-only simulations with prescribed SST and sea ice cover using the AMIP 2000–2015 climatology with present-day aerosol emissions and greenhouse gas concentrations for comparison with observations, (ii) atmosphere–MLO simulations at 1×CO2 concentrations with preindustrial aerosol emissions and preindustrial levels of the other greenhouse gas, and (iii) atmosphere–MLO simulations at 2×CO2 concentrations, keeping everything else as in (ii). For simulations (iii) the CO2 concentration was doubled before the 25-year spin-up period.

4 Comparison of ECHAM6-HAM2 with observations

An overview of the annual, global mean state of the climate in the different simulations is given in Table 2 together with observations, where available. The observations of global liquid water path (LWP) range between 30 and 90 g m−2 . restricted the LWP retrievals to ocean regions (LWPoc) and estimated LWPoc from a combination of the SSM/I, TMI, AMSR-E, WindSat, SSMIS, AMSR-2, and GMI satellite data (MAC-LWP) as 81.4 g m−2. Limiting the other retrievals to ocean regions yields a global LWPoc of 42.9 g m−2 from MODIS , 43.9 g m−2 from ATSR-2-AATSR , and 41 g m−2 from AVHRR-PM , illustrating the huge uncertainty of estimating LWPoc. All simulations except for ALL_LIQ fall within this range. In simulation ALL_LIQ, LWPoc amounts to 110 g m−2 because no ice is formed at temperatures >35 C.

The annual zonal mean LWP and ice water path (IWP) from all simulations as well as from satellite observations of LWP from multiple satellite sensors , the MODIS satellite , and the ATSR-2-AATSR satellites and IWP from CALIPSO/CloudSat are shown in Fig. 1. Because the MAC-LWP observations are only available over oceans, we limit the comparison of LWP to ocean regions. Both retrievals from visible–near-infrared sensors and microwave sensors have biases in retrieving LWP . mentioned four biases in LWP retrievals for MODIS (visible–near-infrared sensor) and AMSR-E (microwave sensor): a bias at large solar zenith angles (MODIS), missing pixels of low-lying clouds that are not detected (MODIS), LWP retrievals in cloud-free scenes (AMSR-E), and the partitioning of the microwave signal between cloud and precipitation signals (AMSR-E). The LWP bias of microwave-sensor-based products in clear-sky scenes has been recently corrected by . The LWP estimates of microwave retrieval products are most reliable in regions where LWP is much larger than the rainwater path (RWP). Thus, recommend restricting the LWP satellite data to regions in which LWP(LWP + RWP) > 0.8. This LWPoc,low_pr covers areas dominated by stratocumulus, trade wind cumulus regions in the subtropics, the Southern Ocean, high latitudes in the Northern Hemisphere, and parts of the north Atlantic and north Pacific. LWPoc,low_pr from MAC-LWP should thus have the smallest retrieval biases. Therefore we use it to evaluate the model LWP as shown in Fig. 1.

Table 2Global annual mean of oceanic LWPoc; LWPoc,low_pr; IWP; vertically integrated global cloud droplet (Nl); CDNC at cloud top over oceans (${N}_{\mathrm{l},\mathrm{oc},\mathrm{top}}$); ice crystal number concentration (Ni); total cloud cover (CC); precipitation rate (P); shortwave (SW), longwave (LW), and net cloud radiative effect (CRE) from observations; and 20-year atmosphere-only present-day model simulations with prescribed SST as described in Table 1 for 2003–2012. The satellite observations of LWPoc are taken from , , , and ; of LWPoc,low_pr from ; of IWP from the CloudSat/CALIPSO satellites ; of ${N}_{\mathrm{l},\mathrm{oc},\mathrm{top}}$ from of total cloud cover from and ; of precipitation rate from the GPCP version 2.3 data set . The CRE satellite data are averaged over the period 2001–2011 from the Clouds and the Earth's Radiant Energy System Energy Balanced and Filled Ed2.6r data set and are described in . For the uncertainty range of the SW CRE data from CloudSat and CALIPSO are also considered and for LW CRE data from the TOVS satellite are also considered . See legend of Fig. 1 for details of the observations.

The new multi-satellite observations of LWPoc by , averaged over the years 1988–2016, show a maximum LWP in the tropics and, on average, their estimate is twice as high as the retrievals from MODIS and ATSR-2-AATSR both for the years 2003–2012. LWPoc from MODIS and ATSR-2-AATSR, conversely, increases towards the poles and only a weak secondary maximum is indicated near the Equator. All observations show high values of LWPoc in the storm tracks in both hemispheres, but the exact location varies among the different observations. The most noticeable difference between LWPoc and LWPoc,low_pr is the absence of data in the intertropical convergence zone (ITCZ). In fact, the agreement of simulations REF and HET in terms of the zonal mean structure of LWPoc,low_pr with observations is very good everywhere and better than with LWPoc. LWPoc,low_pr is too low at high latitudes in simulation ALL_ICE and the maxima in the extratropics are overestimated in simulation ALL_LIQ. LWPoc,low_pr is underestimated everywhere in simulation NOCONV because of the large speedup of the autoconversion rate required to reach radiative equilibrium at the TOA in this set-up.

Figure 1Annual zonal mean LWPoc (a), IWP (b), LWPoc,low_pr (c), cloud-top CDNC of clouds between 268 and 300 K (d), SW CRE (e), and LW CRE (f) from the atmosphere-only present-day simulations for the years 2003–2012 as described in Table 1 and observations of LWP from multiple satellite sensors for the years 1988–2016 by (solid line), from MODIS-AQUA collection 6.1 for the years 2003–2012 (dotted line), from ATSR-2-AATSR for the years 2003–2012 (dashed line), CloudSat + CALIPSO satellite observations of IWP for the years 2006–2010 from , low-precipitating oceanic LWP from for the years 2003–2012, oceanic cloud-top CDNC for the years 2003–2015 from , SW and LW CRE. The solid SW and LW CRE lines are from CERES 2003–2012 , the dashed ones from ERBE for the years 1985–1989, and the dotted ones for LW CRE from TOVS satellite data for the years 1985–1993 .

The cloud-top CDNC for oceanic clouds (${N}_{\mathrm{l},\mathrm{oc},\mathrm{top}}$) with temperatures between 268 and 300 K has been obtained from satellite data for the years 2003–2015. From the model simulations, oceanic cloud-top CDNCs were used for temperatures > 273.2 K. ${N}_{\mathrm{l},\mathrm{oc},\mathrm{top}}$ is less zonally symmetric than LWP as a result of the higher emissions in the Northern Hemisphere. ${N}_{\mathrm{l},\mathrm{oc},\mathrm{top}}$ reaches values of > 100 cm−3 between 30 and 60 N in the observations (Fig. 1). All model simulations fail to produce such high values north of 40 N, probably because of an insufficient transport of aerosol particles to the Arctic , missing nitrate aerosols, and underestimation of particulate organic aerosols , all of which contribute to underestimating the concentration of cloud condensation nuclei. ${N}_{\mathrm{l},\mathrm{oc},\mathrm{top}}$ is well simulated south of 30 N, especially in simulation ALL_LIQ. ${N}_{\mathrm{l},\mathrm{oc},\mathrm{top}}$ is highest in simulation ALL_ICE (Table 2) and overestimated with respect to the observations. This is caused by the lower speedup of the autoconversion rate in simulation ALL_ICE (Table 1), which reduces the sink of cloud droplets in this simulation and causes higher CDNC at all altitudes compared to simulation REF (not shown).

The zonal distribution of ${N}_{\mathrm{l},\mathrm{oc},\mathrm{top}}$ in simulation NOCONV differs from the other simulations. It peaks near the Equator because the clouds here cannot form by convection and instead the stratiform cloud scheme needs to take over. As the rain formation depends on CDNC only in the stratiform cloud microphysics scheme but not in the convective one, warm rain is less efficiently formed in the tropics in the stratiform scheme, causing a buildup of cloud droplets in simulation NOCONV as portrayed in Fig. 1. The CDNC peak in the tropics in simulation NOCONV deviates strongly from observations, as do the lower-than-observed values in the extratropics.

ECHAM6-HAM2 in general has problems simulating IWP values in the observed range. This deficiency can largely be attributed to an underestimation in the ice crystal size in simulation REF at least for clouds with a COD < 3 . The only simulation in which IWP in the global mean agrees with the observations is simulation HET (Table 2) because the heterogeneous nucleation is limited by the number concentration of dust aerosols. The dust number concentration is smaller than the number concentration of soluble aerosols that can freeze homogeneously and freezing commences at a lower relative humidity and a higher temperature in simulation HET. Because of the limited number concentration of dust aerosols and the faster depositional growth at higher temperatures, the heterogeneously nucleated ice crystals are larger than the ones nucleated homogeneously. Because of this, they sediment and aggregate more rapidly, causing simulation HET to be initially out of radiative balance. Here retuning was necessary to slow down the rate of aggregation and at the same time speed up the rate of autoconversion. This resulted in one of the lowest LWPs of all simulations and the highest IWPs with the highest ice crystal number concentration.

The IWP peak in the tropics in Fig. 1 is related to liquid-origin cirrus clouds forming the anvils from deep convection in the ITCZ as well as in the warm conveyor belt in the storm tracks. The peak in the ITCZ is not captured in any of the model simulations, suggesting that the model severely underestimates the detrained cloud ice in deep convective clouds. Cloud ice in the extratropics is underestimated in all simulations but simulations HET and ALL_ICE. IWP is largest in simulation HET due to retuning as discussed above. Its zonal distribution shows that the underestimation of IWP in the tropics is compensated for by an overestimation in the extratropics. Simulation ALL_ICE best matches the magnitude of IWP in the extratropics but shifts the peaks slightly polewards of the observations.

The total cloud cover in ECHAM6-HAM2 is that of large-scale clouds because convective clouds are considered to be short lived and to decay within one model time step except for the detrained condensate in the anvils that is taken as a source for the stratiform cloud scheme. The global mean total cloud cover does not vary much among the different simulations and it falls within the observed range of 68 ±5  % in all simulations. It is second-largest in simulation ALL_LIQ because of its large LWP. It is highest in simulation NOCONV because here all clouds are large scale and contribute to the cloud cover.

The observed precipitation rate from GPCP of 2.7 mm d−1 averaged over 1981–2010 is overestimated by 6–13 % in all simulations, a feature that ECHAM6-HAM2 shares with its host model ECHAM6 . As discussed in , GPCP seems to underreport precipitation, and higher precipitation rates are more consistent with the best observational estimates of the surface energy budget .

The cloud radiative effect (CRE) is the difference between the all-sky radiation at the TOA and the clear-sky radiation. The shortwave (SW) CRE amounts to 47.3 W m−2 with a range from 46 to 53.3 W m−2 in the satellite data . In all simulations, except ALL_LIQ, SW CRE amounts to about 50 W m−2, which lies within the observed range. SW CRE is most negative and outside the observational range in simulation ALL_LIQ because of its high LWP and total cloud cover (see also Fig. 1). Conversely, it is least negative in simulation ALL_ICE, in which LWP is rather small and hence τ is smallest. Here SW CRE is severely underestimated south of 40 S due to the lack of supercooled liquid water. While too much absorption of shortwave radiation over the Southern Ocean has been a problem in the many GCMs because of insufficient supercooled liquid water , in our model this problem only arises in the extreme simulation ALL_ICE.

LW CRE amounts to 26.2 W m−2 with a range from 22 to 30.5 W m−2 in the satellite data . As for SW CRE, all simulations but ALL_LIQ calculate LW CRE to be about 24–25 W m−2 in good agreement with observations. Again LW CRE is smallest in simulation ALL_ICE and highest and outside the observational range in simulation ALL_LIQ (see also Fig. 1). Due to the severe underestimation of cloud ice in the tropics in all simulations, LW CRE is systematically underestimated in the tropics. The exception is simulation ALL_LIQ in which the underestimation in tropical cloud ice is compensated for by the thickest high level clouds (Fig. 2).

Figure 2Annual global mean cloud liquid water content (a), ice water content (b) in milligrams per kilogram, and cloud cover (c) in percent as a function of pressure in hectopascals from the various sensitivity atmosphere-only present-day simulations for the years 2003–2012 described in Table 1 and observations of the ice water content from CALIPSO-CloudSat for the years 2006–2010 by .

The global mean observed net CRE ranges between 17.1 and 22.8 W m−2. Here all simulations are more negative than observed, i.e., they overestimate the net negative radiative effect of clouds.

The vertical distribution of the globally and annually averaged cloud liquid water and cloud ice is shown in Fig. 2. Cloud liquid water has its maximum at around 800 hPa associated with low clouds where it varies between 20 and 30 mg kg−1 in the different simulations. It is lowest in simulation NOCONV because of the drastically enhanced autoconversion rate (Table 1). Cloud liquid water decreases to zero at around 400 hPa except in simulation ALL_LIQ in which ice formation is limited to temperatures <35 C. Unfortunately no observational data of the vertical distribution of cloud liquid water are available.

The vertical distribution of cloud ice has been derived from CALIPSO/CloudSat . It peaks at 400 hPa with 6 mg kg−1 and becomes negligible at altitudes above 100 hPa and below 900 hPa (Fig. 2). Here the differences among the sensitivity experiments are more pronounced. In simulations REF and NOCONV the shape of the distribution looks similar to the observed distribution, but the peak value is underestimated by 30–50 % and in simulation REF a secondary maximum is present at around 780 hPa. This peak is related to cloud ice over the Southern Ocean and Antarctica (not shown). In simulation ALL_LIQ, in which the global annual mean IWP is smallest due to suppressed ice formation in mixed-phase clouds, cloud ice peaks at 250 hPa and drops to zero at 600 hPa. This simulation differs most from the observations. It also has the lowest cloud cover in the lower troposphere and the highest cloud fraction between 200 and 400 hPa (Fig. 2). Conversely, simulation NOCONV has the highest coverage of low-level clouds and the smallest coverage of cirrus clouds.

In simulation ALL_ICE the comparison of cloud ice with observations is also less favorable compared to simulation REF because more cloud ice than observed is simulated in the lower atmosphere while the underestimation of cloud ice at higher altitudes has not noticeably improved. The overall best agreement with observations is seen in simulation HET. This is the only simulation in which the global annual mean IWP lies within the observational uncertainty (Table 2). However, the peak in cloud ice at 400 hPa is 25 % too large and too narrow. Although the freezing mechanism was only changed for cirrus clouds, the overall higher amount of cloud ice extends to lower altitudes, causing an overestimation in cloud ice between 600 and 900 hPa. This is mainly a result of the reduced speedup of the aggregation rate that also affects cloud ice in mixed-phase clouds.

SLF has been obtained from the CALIOP satellite and has been compared to an earlier version of ECHAM6-HAM2 . At that time ECHAM6-HAM2 seriously underestimated SLF. This can partly be explained by differences in the definition of SLF in the satellite data and in the model. In addition, some of the model improvements mentioned above were added after the study by . As the lidar signal becomes attenuated at COD > 3, the CALIOP satellite data are only representative for all cloud tops not overlaid by clouds with COD > 3. SLF of the model data is diagnosed similarly to CALIOP starting from the highest cloud layers as long as the cumulative COD of these cloud layers is less than or equal to 3. As shown in Fig. 3, ECHAM6-HAM2 still underestimates SLF but on average simulates SLFs twice as high as presented in . At 10 C SLF amounts to 63 % in CALIOP but only 37 % in simulation REF. The difference between the observed and simulated SLF decreases at colder temperatures because of the general decrease in SLF with decreasing temperature.

SLF from all ECHAM simulations is significantly underestimated except for simulations ALL_LIQ and NOCONV. Simulation NOCONV actually matches the observed SLF rather well, mainly because supercooled liquid water exists at higher altitudes than in simulation REF (Fig. 2). This points to a potential deficiency as to how we handle detrainment from convective clouds or convection itself. In simulation REF we assume that if we are in the Wegener–Bergeron–Findeisen (WBF) regime, following the definition of as described in , and cloud ice is already present, the detrained condensate will be in the form of ice. It is only detrained as supercooled liquid water if the vertical velocity is sufficiently high to exceed saturation with respect to liquid water. This assumption seems to cause a too efficient WBF process, thus depleting the supercooled liquid water too rapidly. We will rethink our approach in the future.

The control simulation in CESM from simulates a SLF of only 20 % at 10 C. Based on their hypothesis that an underestimation in SLF translates into an underestimation in ECS, a smaller underestimation of SLF in simulation REF should lead to a smaller underestimation of ECS in ECHAM6-HAM2. Based on the simulated SLFs, we expect ECS to be lowest in simulation ALL_ICE and to successively increase in simulations HET, REF, NOCONV, and ALL_LIQ.

Figure 3Frequency of occurrence of SLF in percent for clouds with COD < 3 in different temperature bins from CALIOP observations for the years 2006–2010 and the various sensitivity simulations from the atmosphere-only present-day simulations for the years 2003–2012 described in Table 1.

5 Equilibrium climate sensitivity

In a warmer climate, the saturation-specific humidity increases by 7 % C−1 of temperature increase according to the Clausius–Clapeyron equation. Because the relative humidity has been found to remain rather constant , this causes an increase in the specific humidity and speeds up the hydrological cycle leading to higher liquid water contents in clouds and higher precipitation rates.

ECS from all our sensitivity simulations is shown in Table 3. It ranges between 1.8 and 2.6 K and with that lies on the lower side of the range of ECS estimated in the last IPCC report and from CMIP5 models . Similar to , ECS increases by  30 % when increasing SLF from almost zero to one (simulations ALL_ICE and ALL_LIQ in this paper vs. Low-SLF and High-SLF in ). Contrary to their large increase in ECS from their reference simulation to simulation High-SLF, ECS does not change among simulations REF, NOCONV, and ALL_LIQ but remains at  2.5–2.6 K despite the increase in SLF at 10 C from  30 % in HET and  40 % in REF to  70 % in NOCONV and 100 % in ALL_LIQ, i.e., only ECS in simulation ALL_ICE is noticeably different. This suggests that ECHAM6-HAM2 is more sensitive at low SLF. These results seem to support the hypothesis that an increase in ECS is only sensitive to SLF if too little shortwave radiation is reflected back to space from mixed-phase midlatitude clouds, i.e., when they are composed of ice instead of liquid water. This has also been suggested by , who modified CAM5 to detrain more liquid water at colder temperatures in shallow convective clouds, which occur in the cold sector of midlatitude cyclones, for example. It is important that these mixed-phase clouds are not shielded by ice clouds in order for a change in their cloud phase to affect ECS . Since in ECHAM6-HAM2 this shortwave bias is most pronounced in simulation ALL_ICE, this is the only simulation with a distinctively lower ECS. In the reference version of CAM5, too much shortwave radiation is absorbed over the Southern Ocean in their present-day climate , which explains why ECS is also sensitive to an increase in SLF in this model . In addition, did not use a MLO, but a fully coupled dynamic ocean, which could impact the comparison with our results because different methods can lead to differing ECS estimates .

Table 3Global annual mean changes of surface temperature (ΔTs), liquid water path (ΔLWP), ice water path (ΔIWP), vertically integrated cloud droplet (ΔNl) and ice crystal number concentration (ΔNi), total cloud cover (ΔCC), SW CRE (ΔSW CRE), LW CRE (ΔLW CRE), and net CRE from the radiative kernel (RK) method and as normally diagnosed in equilibrium after a CO2 doubling averaged over the years 26 to 50 of the coupled atmosphere – MLO model simulations as well as the net clear-sky feedback parameter λcs diagnosed according to . The individual simulations are described in Table 1. Different from Table 2, here the changes in LWP are global and not limited to ocean regions.

In order to prove our hypothesis for the different relationship between SLF and ECS in ECHAM6-HAM2, we evaluate the changes in cloud fraction as a function of cloud-top pressure and COD between the 2×CO2 and 1×CO2 climates. Because mixed-phase clouds are more prevalent in midlatitudes and high latitudes than in the tropics and subtropics , we expected the largest differences there. The histograms in Figs. 4 and 5 are taken polewards of 40 latitude in both hemispheres. They are taken from a satellite perspective, by using the COSP-ISCCP simulator . This means that the cloud-top pressure refers to the highest cloud in each single sub-column used by the simulator. Sub-columns are used in the simulator to account for sub-grid-scale variability in the cloud (top) distribution. From Fig. 4 we see that only in simulation ALL_ICE is there a large abundance of low- and mid-level optically thin extratropical clouds (with τ<1.3) in the present climate. These clouds show a marked decrease accompanied with a marked increase in medium optical depth clouds ($\mathrm{3.6}<\mathit{\tau }<\mathrm{60}$) between 560 and 800 hPa in the 2×CO2 climate (Fig. 5). These clouds and their changes are absent in all other simulations with a higher SLF (REF, HET, NOCONV, and ALL_LIQ). Only in simulation ALL_ICE are these optically thin ice clouds not shielded by higher clouds prevalent in the present-day climate (not shown) and are converted to low- and mid-level liquid water clouds of medium optical depth in the 2×CO2 climate.

Figure 4Distribution of cloud fraction in the extratropics (> 40 S and N) as a function of cloud optical depth and cloud-top pressure between the 2×CO2 and the 1×CO2 climate in simulations ALL_ICE, HET, REF, NOCONV, and ALL_LIQ from the last 25 years of the MLO simulations described in Table 1.

Figure 5Changes in the distribution of cloud fraction in the extratropics (> 40 S and N) as a function of cloud optical depth and cloud-top pressure between the 2×CO2 and the 1×CO2 climates in simulations ALL_ICE, HET, REF, NOCONV, and ALL_LIQ from the last 25 years of the MLO simulations described in Table 1.

To estimate the radiative effect of this increase in COD, we calculate the different components of the global cloud feedback parameter λc using the radiative kernel decomposition method described in . This method decomposes λc into feedbacks that are associated with changes in cloud amount λamt, cloud-top pressure λctp, and COD λτ as shown in Fig. 6 for all simulations. λctp is positive because of the shift of clouds to higher altitudes in the warmer climate (see also Fig. 7) that enhances their LW CRE. At the same time, the total cloud amount decreases in a warmer climate, most noticeable at lower altitudes (Fig. 7). This decrease in low/mid-level clouds reduces the negative net CRE and also constitutes a positive feedback. λτ is negative in most simulations and smaller than other cloud feedbacks so that λc is positive in most simulations. λc is negative (0.4 W m−2 K−1) only in simulation ALL_ICE because of its large negative λτ. In all other simulations λc ranges between 0.3 W m−2 K−1 and 0.6 W m−2 K−1. The value of 0.6 W m−2 K−1 corresponds to the mean of the analyzed GCMs in AR5 . The estimates from the other simulations (except ALL_ICE) fall within the 90 % range of the cloud feedback of 0.2 to 0.6 W m−2 K−1 as assessed in AR5 .

Figure 6Components of the globally averaged cloud feedback parameter λc in W m−2 K−1 for the simulations described in Table 1. λτ accounts for changes in cloud optical depth, λamt accounts for changes in cloud amount, and λctp accounts for changes in cloud-top pressure.

Figure 7Change in the annual global mean cloud liquid water content (LWC, a) and ice water content (IWC, b) in milligrams per kilogram as a function of pressure in hectopascals between the 2×CO2 and 1×CO2 climates for the simulations described in Table 1.

The largest differences among the simulations are associated with changes in λτ, which varies between 0.6 W m−2 K−1 in simulation ALL_ICE and 0.1 W m−2 K−1 in simulation NOCONV. In simulation ALL_ICE, the large negative value of λτ can be explained with the largest increase in the LWP of 4.8 g m−2 (Table 3) due to the phase change from cloud ice to cloud liquid water mainly in low- and mid-level midlatitude clouds (Fig. 5). This causes a negative optical depth feedback in all regions, but most pronounced polewards of 40 (Fig. 8). Note that simulation ALL_ICE is the only simulation in which the overall cloud feedback parameter is dominated by changes in COD. The results of simulation ALL_ICE agree qualitatively with the results from for regions with large-scale subsidence in which optically thin ice clouds in the present-day climate are not shielded by higher clouds and hence for which the COD feedback is most important. In all other simulations, the cloud feedback is dominated by changes in cloud amount and cloud-top pressure (Fig. 6).

Figure 8Tropical, subtropical, extratropical, and global cloud optical depth feedback parameter λτ in W m−2 K−1 for the sensitivity simulations described in Table 1.

In the other simulations, shown in Fig. 5 for midlatitudes and high latitudes, a rather different picture emerges. Here τ of low clouds hardly changes or decreases as the majority of them are already composed of cloud droplets in the present climate and only a small phase change from ice to liquid occurs due to the doubling of CO2. In all simulations, except for simulation NOCONV, a negative cloud-phase feedback is also visible in the tropics. This negative feedback arises because some tropical clouds that glaciate at the tops, either due to the presence of ice-nucleating particles in the mixed-phase temperature regime or because their tops extend to altitudes with temperatures <35 C in the present climate, will remain supercooled in the warmer climate and hence become optically thicker (Figs. 8 and 9).

Figure 9As Fig. 5 but for the tropics (15 S–15 N).

In the warmer climate, the static stability increases in areas with marine stratus, stratocumulus, and trade wind cumuli in all simulations. This is accompanied by a stronger moisture gradient, which promotes stronger drying by more entrainment , which in turn causes the marine subtropical clouds to become thinner and their optical depth to decrease. In most simulations this leads to a positive COD feedback in the subtropics as shown in Fig. 8.

The parameterization of convective clouds in ECHAM6-HAM2 assumes that they form and dissipate in one time step. During their lifetime, they can detrain cloud water and ice in the environment. Thus the negative COD feedback in the tropics and midlatitudes indicates that more cloud condensate is detrained in the form of liquid or supercooled liquid water rather than as cloud ice in the warmer climate. Thus, next to the rise of the melting level in the warmer climate, the rise of the homogeneous freezing level is also important for the negative cloud-phase feedback. This is best seen in simulation ALL_LIQ in which no ice exists at mixed-phase temperatures, yet the negative cloud-phase feedback still operates in the tropics (Figs. 8 and 9). In the global mean, this negative cloud-phase feedback is visible in all our simulations as a simultaneous increase in cloud liquid water and decrease in cloud ice in the warmer climate at a given pressure level (Fig. 7). In midlatitudes and the tropics it shows up as an increase in optically thicker clouds and decrease in optically thinner clouds for the same cloud-top pressure (Figs. 5 and 9).

The second most important contributor to the differences in the overall cloud feedback between simulations ALL_ICE and ALL_LIQ is the change in cloud-top pressure (λctp). λctp varies between 0.1 W m−2 K−1 in simulation ALL_ICE and 0.4 W m−2 K−1 in simulation ALL_LIQ. The cloud-top pressure feedback is mainly related to cirrus clouds and has been estimated from satellite observations to amount to 0.2 W m−2 K−1 . We obtain a cloud-top pressure feedback of about this value in all other simulations. We also computed the cloud feedbacks separately for low and non-low clouds following the decomposition by , which leads to a weaker cloud-top pressure feedback but in general qualitatively similar results for all simulations (not shown).

λctp is largest in simulation ALL_LIQ in which the global mean changes in cloud liquid water and cloud ice are distinctively different from all other simulations (Fig. 7). It is the only simulation in which cloud liquid water decreases throughout the lower and mid-troposphere because here the poleward shift of the storm tracks and the upward shift of convective clouds are least compensated for by changes from cloud ice to cloud liquid water. The changes in simulation ALL_LIQ are accompanied with the largest decrease in cloud cover between 750 and 370 hPa. Cloud liquid water increases between 250 and 400 hPa at the expense of cloud ice in simulation ALL_LIQ. In all other simulations, the concurrent increase in cloud liquid water and decrease in cloud ice occurs at altitudes below 400 hPa. The different vertical structure of cloud liquid water and cloud ice in simulation ALL_LIQ is a consequence of the design of this simulation. Because cloud ice only forms at temperatures below 35 C, cloud droplets are carried upwards to this temperature, such that the optically thick clouds extend to higher altitudes in the 2×CO2 climate (see Figs. 5 and 9), causing the largest cloud-top pressure feedback as shown in Fig. 6.

The opposite effect can be seen in simulation ALL_ICE, in which all cloud water is converted to ice already at 0 C. The increase in cloud liquid water in low and mid-level clouds in the warmer climate is largest in this simulation mainly because of the negative cloud-phase feedback (Fig. 7). This combined with only a modest increase in cloud ice and cloud cover between 250 and 100 hPa (Fig. 7) causes the cloud-top pressure feedback to be smallest in this simulation (Fig. 6).

Simulation NOCONV is the only simulation in which λτ is slightly positive. λτ is negative in all other simulations in the tropics (Fig. 8) because more cloud condensate is detrained in the form of liquid or supercooled liquid water in the warmer climate. The absence of a convection parameterization and hence detrainment in simulation NOCONV on average leads to a positive tropical λτ, i.e., the tropical and subtropical clouds in simulation NOCONV become optically thinner.

To summarize, the negative COD feedback only dominates the overall cloud feedback if the low- and mid-level midlatitude clouds consist of ice in the present climate and change to liquid in the 2×CO2 climate as in simulation ALL_ICE. However, this does not imply that the overall cloud feedback remains more or less constant in the other simulations. In fact, λc almost doubles from simulations REF and HET to simulations NOCONV and ALL_LIQ, caused by the large increase in cloud-top pressure feedback in simulation ALL_LIQ and a positive λτ in simulation NOCONV as explained above. This marked increase in the overall cloud feedback from simulations HET and REF to simulations NOCONV and ALL_LIQ is, however, not reflected in an increase in ECS because the net clear-sky feedback parameter, i.e., the sum of the clear-sky Planck, water vapor, lapse rate, and surface albedo feedbacks, is largest in simulation ALL_LIQ followed by simulation NOCONV (Table 3).

We hypothesize that different processes contribute to the rather similar ECS in all simulations but ALL_ICE: in simulation ALL_LIQ, this seems to be caused by the optically thick clouds that tend to cluster in the tropics and become optically thicker there. This is apparent from Fig. 9, in which the upward shift in cirrus clouds with an increase in cirrus optical depth between 180 and 310 hPa is much more pronounced in simulation ALL_LIQ than in simulation REF. In addition, the negative cloud-phase feedback is restricted to ice clouds that formed at altitudes with temperatures <35 C in the 1×CO2 climate that are now liquid clouds (Fig. 7).

Figure 10Changes in the annual zonal SW and LW net clear-sky radiation in watts per square meter at the top-of-the-atmosphere and surface temperature (K) (a); changes in SW, LW, and net cloud radiative effects between the 2×CO2 and 1×CO2 climates calculated based on the radiative kernel method (b) and normally diagnosed (c) for the simulations described in Table 1. Negative values for the change in clear-sky outgoing longwave radiation denote higher values in the 2×CO2 climate. Please note the different vertical axes.

This increase in optical depth of high clouds causes their SW CRE to be more negative and completely offsets the reduction in LW CRE from the decrease in cloud cover in the tropics in simulation ALL_LIQ (Fig. 10). Such a clustering of convective clouds was found by using the ECHAM4 model when the emissivity of cirrus clouds was set to 1 (black cirrus), i.e., when the infrared optical depth was increased. No convective clustering was seen in the reference simulation in which the cirrus emissivity was calculated as a function of ice water path and in a simulation with transparent cirrus in the infrared. In a warmer climate, convective clouds have been found to cluster further . The clustering of convective clouds in the warmer climate arises because cloud tops remain at nearly the same temperature in a warmer climate (e.g., Hartmann and Larson2002) but are in a more stable environment, which decreases anvil outflow in the upper troposphere and decreases the anvil cloud fraction . In ECHAM this decrease is most pronounced at altitudes below 250 hPa (Fig. 7) from which more longwave radiation is emitted to space. This leads to an overall negative feedback because of the much larger clear-sky area and offsets the positive cloud feedback. In the present-day climate of simulation ALL_LIQ, the least amount of outgoing longwave radiation is emitted to space because the clouds are optically thickest in this simulation. In a warmer climate, these clouds cluster further in the intertropical convergence zone as can be seen from the largest increase in LW CRE as diagnosed from the radiative kernel (RK) method (Fig. 10). Because these clouds are optically thick, this large increase in LW CRE (RK) is overcompensated for by a decrease in SW CRE (RK), causing the net CRE (RK) in the intertropical convergence zone to be as negative in simulation ALL_LIQ as in simulation REF.

The cloud-top pressure feedback in the extratropics that is significantly positive in all simulations but ALL_ICE where it is close to zero (not shown) is overcompensated for by decreases in cloud cover. As shown in Fig. 10, SW CRE and SW CRE (RK) become less negative in the extratropics and dominate over the decrease in LW CRE and LW CRE (RK), causing a positive Δnet CRE at midlatitudes. The generally smaller positive and larger negative changes in net CRE than in net CRE (RK) are caused by compensating changes in non-cloud climate components as described in detail in . In high latitudes, the decrease in surface albedo due to melting of Arctic sea ice causes especially large negative ΔSW CRE and Δnet CRE values polewards of 70 N. In the clear-sky, water vapor and CO2 absorb more longwave radiation in the warmer climate than in the present day. This decreases the difference between all-sky and clear-sky fluxes and causes the LW CRE to be less positive. Differences in surface temperature operate in opposite ways. While more longwave radiation is emitted to space in the warmer climate (negative Planck feedback), more longwave radiation is also absorbed by clouds and re-emitted to the surface (positive cloud feedback). This clear-sky compensation in the longwave is most pronounced in simulation ALL_LIQ, which is the only simulation in which the global mean ΔLW CRE (RK) is positive (Table 3).

The changes between simulation REF and ALL_LIQ are qualitatively comparable to the ones from , for which the difference in climate sensitivity was also rather small between the reference simulation and the black cirrus simulation, both of which also had a positive cloud feedback and an increase in net CRE. There was a noticeable difference between the transparent cirrus and the reference simulation because only the transparent cirrus simulation produced a negative cloud feedback and a decrease in net CRE. While in our simulation ΔCRE is always negative or close to zero, it is most negative in simulation ALL_ICE, which corresponds to the transparent cirrus simulation.

ECS in simulation HET is also similar to simulation REF, but differences in ΔLW CRE (RK) exist. As in simulation ALL_LIQ, in simulation HET the increase in ΔLW CRE (RK) in the tropics is also larger than in simulation REF (Fig. 10) because of the upward shift of ice clouds (Fig. 7) and their increase in optical depth (Fig. 9). Because these ice clouds are of medium optical depth, the large ΔLW CRE (RK) is overcompensated for by a decrease in ΔSW CRE (RK), causing Δnet CRE (RK) and Δnet CRE to be similar to in simulation REF (Table 3).

In simulation NOCONV the outgoing longwave clear-sky radiation is increased the most. While it is also slightly decreased in the tropics, the decrease is smallest in this simulation because of the absence of convection (Fig. 10). The associated absence of detrainment slows down the Wegener–Bergeron–Findeisen process on the one hand, but required the largest speedup of the autoconversion rate on the other hand (Table 1). The latter causes the present-day cloud liquid water and its increase in the warmer climate to be smallest (Tables 2 and 3). This explains the overall positive λτ and the second largest λc. The most negative clear-sky LW radiation changes in NOCONV explain why ECS does not increase compared to in simulation REF. The absence of convection also seems to limit the increase in cirrus clouds (increase in cloud cover at altitudes above 250 hPa; Fig. 7).

6 Conclusions

In this study we used the newly developed ECHAM6.3-HAM2.3 coupled global aerosol–climate model to assess the influence of different cloud processes for ECS. This work was motivated by the findings of using CAM5, who showed a large influence of correcting the underestimated SLF in present-day mixed-phase clouds on ECS. An underestimate of SLF has also been found in other models and in a different version of CAM5 . The SLF was found to be most sensitive to the glaciation rate, which in turn is influenced by cloud ice microphysics and the model's vertical resolution . In ECHAM6-HAM2, SLF could be improved when switching off the convective parameterization entirely. The absence of convection and hence detrained cloud ice slows down the glaciation rate in ECHAM6-HAM2, which corroborates the findings by on the importance of the glaciation rate for SLF.

In ECHAM6-HAM2, ECS is much smaller than in the CAM5 GCM used by . It varies between 1.8 and 2.5 K in ECHAM6-HAM2 versus between 3.9 and 5.7 K in CAM5 . Thus, while ECS is on the low side of the ECS range between 2.1 and 4.7 K in IPCC AR5 in ECHAM6-HAM2, it is on the high side in CAM5. Note that the percentage increase in ECS of 30 % between the extreme scenarios (ALL_ICE to ALL_LIQ in ECHAM6-HAM2 vs. Low-SLF to High-SLF in CAM5) is similar, indicating that the relative contribution of the negative cloud-phase feedback to the overall ECS is comparable in both GCMs.

However, important differences between the two GCM studies are apparent: increasing SLF in ECHAM6-HAM2 from its calculated values in the reference simulation to one as in simulation ALL_LIQ does not increase ECS in contrast to the findings by . Part of this difference could be caused by the overestimation of shortwave absorption over the Southern Ocean due to too efficient freezing in low- and mid-level shallow convective clouds in CAM5 . We hypothesize that it is the SLF of these optically thin low- and mid-level midlatitude clouds not shielded by overlying clouds that matters for ECS because the cloud-top temperature of these clouds is in the mixed-phase temperature range. Hence, the cloud phase at the cloud top of these clouds plays an important role in the TOA radiation budget. If the absorption of shortwave radiation over the Southern Ocean is correctly simulated, then SLF in other clouds does not matter for ECS. At least this is what the ECHAM6-HAM2 results show because only in simulation ALL_ICE is significantly less shortwave radiation reflected back to space from Southern Ocean clouds and this is the only simulation in which ECS is significantly smaller than in all other simulations. In all other simulations sufficient or even too much shortwave radiation is reflected back to space from Southern Ocean clouds and all of them have similar values of ECS.

The reason why an underestimation of SLF in cloud types other than thin midlatitude clouds with cloud-top temperatures between 0 and 35 C does not seem to matter is because radiative changes in clouds in the warm sector of extratropical cyclones are masked by ice clouds above them . In addition, the tops of tropical deep convective and deep frontal clouds consist of ice and that will not change in a warmer climate. Low-level clouds in the tropics and subtropics already consist of liquid water and therefore their cloud phase will also not change in the future climate.

In our model, it is not only the cloud-phase feedback and the overall cloud feedback that matters for ECS. If this were the case, ECS should be largest in simulation ALL_LIQ, in which the cloud feedback parameter is highest. It seems that in simulation ALL_LIQ tropical deep convective clouds tend to aggregate, which causes a negative feedback by increasing the clear-sky area and with that the longwave emission from clear-sky regions to space . This negative feedback reduces the positive cloud feedback and causes the changes in the net cloud radiative effect and ECS to be the same in simulations REF and ALL_LIQ. Also in simulation NOCONV, in which we switched off the convection parameterization, SLF is higher than in simulation REF because of a reduced Wegener–Bergeron–Findeisen process. Again, ECS remains the same because of the smallest increase in cirrus cloud cover in this simulation that allows more clear-sky longwave radiation emission to space than in the reference simulation.

As discussed by , while cloud-phase improvements in the extratropics affect ECS in their model, they do not seem to matter for the warming during the 21st century in the Community Earth System Model (CESM) because of compensating responses in ocean circulation. Whether this is also the case in other Earth system models, will be subject to future investigations.

Data availability
Data availability.

Model data are available from the authors upon request.

Competing interests
Competing interests.

The authors declare that they have no conflict of interest.

Special issue statement
Special issue statement.

This article is part of the special issue “BACCHUS – Impact of Biogenic versus Anthropogenic emissions on Clouds and Climate: towards a Holistic UnderStanding (ACP/AMT/GMD inter-journal SI)”. It is not associated with a conference.

Acknowledgements
Acknowledgements.

The ECHAM-HAMMOZ model is developed by a consortium composed of ETH Zurich, Max Planck Institut für Meteorologie, Forschungszentrum Jülich, University of Oxford, the Finnish Meteorological Institute, and the Leibniz Institute for Tropospheric Research, and it is managed by the Center for Climate Systems Modeling (C2SM) at ETH Zurich. The research leading to these results has partly received funding from the Center for Climate System Modelling (C2SM) at ETH Zurich, the European Union's Seventh Framework Programme (FP7/2007–2013) project BACCHUS under grant agreement no. 603445, and the Swiss National Science Foundation (project number 200021_160177). We thank Thorsten Mauritsen, Trude Storelvmo, Alejandro Bodas-Salcedo, and the anonymous reviewers for useful comments and suggestions, Sylvaine Ferrachat for the help with earlier simulations on this topic, and Mark Zelinka and Thorsten Mauritsen for providing the radiative kernel used in this study.

The computing time for this work was supported by a grant from the Swiss National Supercomputing Centre (CSCS) under project ID s652 and from ETH Zurich. The data used are listed in the references and/or can be found here. CDNC climatology: ; MAC-LWP: Goddard Earth Sciences Data and Information Services Center (GES DISC, current hosting: http://disc.sci.gsfc.nasa.gov, last access: 20 June 2018); ATSR-2-AATSR v2.0: ; AVHRR-PM v3.0 (prototype): ; MODIS-AQUA collection 6.1: NASA Goddard (https://ladsweb.nascom.nasa.gov, last access: 20 June 2018).

Edited by: Holger Tost
Reviewed by: Trude Storelvmo and two anonymous referees

References

Abdul-Razzak, H. and Ghan, S. J.: A parameterization of aerosol activation: 2. Multiple aerosol types, J. Geophys. Res., 105, 6837–6844, 2000. a

Adler, R. F., Huffman, G. J., Chang, A., Ferraro, R., Xie, P. P., Janowiak, J., Rudolf, B., Schneider, U., Curtis, S., Bolvin, D., Gruber, A., Susskind, J., Arkin, P., and Nelkin, E.: The version-2 global precipitation climatology project (GPCP) monthly precipitation analysis (1979–present), J. Hydrometeorol., 4, 1147–1167, 2003. a, b

Adler, R. F., Gu, G., and Huffman, G. J.: Estimating Climatological Bias Errors for the Global Precipitation Climatology Project (GPCP), J. Appl. Meteor. Clim., 51, 84–99, 2012. a

Barrett, A. I., Hogan, R. J., and Forbes, R. M.: Why are mixed-phase altocumulus clouds poorly predicted by large-scale models? Part 1. Physical processes, J. Geophys. Res., 122, 9903–9926, 2017a. a, b, c

Barrett, A. I., Hogan, R. J., and Forbes, R. M.: Why are mixed-phase altocumulus clouds poorly predicted by large-scale models? Part 2. Vertical resolution sensitivity and parameterization, J. Geophys. Res., 122, 9927–9944, 2017b. a, b

Bennartz, R. and Rausch, J.: Cloud Droplet Number Concentration Climatology, available at: http://hdl.handle.net/1803/8374 (last access: 21 June 2018), 2016. a

Bennartz, R. and Rausch, J.: Global and regional estimates of warm cloud droplet number concentration based on 13 years of AQUA-MODIS observations, Atmos. Chem. Phys., 17, 9815–9836, https://doi.org/10.5194/acp-17-9815-2017, 2017. a, b, c

Block, K. and Mauritsen, T.: Forcing and feedback in the MPI-ESM-LR coupled model under abruptly quadrupled CO2, J. Adv. Model. Earth Syst., 5, 676–691, 2013. a

Bodas-Salcedo, A.: Cloud condensate and radiative feedbacks at midlatitudes in an aquaplanet, Geophys. Res. Lett., 45, 3635–3643, https://doi.org/10.1002/2018GL077217, 2018. a, b

Bodas-Salcedo, A., Webb, M. J., Bony, S., Chepfer, H., Dufresne, J. L., Klein, S. A., Zhang, Y., Marchand, R., Haynes, J. M., Pincus, R., and John, V. O.: COSP Satellite simulation software for model assessment, B. Am. Meteorol. Soc., 92, 1023–1043, 2011. a

Bodas-Salcedo, A., Williams, K. D., Ringer, M. A., Beau, I., Cole, J. N. S., Dufresne, J.-L., Koshiro, T., Stevens, B., Wang, Z., and Yokohata, T.: Origins of the Solar Radiation Biases over the Southern Ocean in CFMIP2 Models, J. Climate, 27, 41–56, 2014. a

Bodas-Salcedo, A., Andrews, T., Karmalkar, A. V., and Ringer, M. A.: Cloud liquid water path and radiative feedbacks over the Southern Ocean, Geophys. Res. Lett., 43, 10938–10946, 2016. a

Bony, S., Stevens, B., Coppin, D., Becker, T., Reed, K. A., Voigt, A., and Medeiros, B.: Thermodynamic control of anvil cloud amount, P. Natl. Acad. Sci. USA, 113, 8927–8932, 2016. a, b

Boucher, O., Randall, D., Artaxo, P., Bretherton, C., Feingold, G., Forster, P., Kerminen, V.-M., Kondo, Y., Liao, H., Lohmann, U., Rasch, P., Satheesh, S. K., Sherwood, S., Stevens, B., and Zhang, X.-Y.: Clouds and Aerosols, in: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Stocker, T., Qin, D., Plattner, G.-K., Tignor, M., Allen, S. K., Boschung, J., Nauels, A., Xia, Y., Bex, V., and Midgley, P. M., 571–657, Cambridge Univ. Press, Cambridge, UK and New York, NY, USA, 2013. a, b, c, d, e

Bourgeois, Q. and Bey, I.: Pollution transport efficiency toward the Arctic: Sensitivity to aerosol scavenging and source regions, J. Geophys. Res., 116, D08213, https://doi.org/10.1029/2010JD015096, 2011. a

Brovkin, V., Boysen, L., Raddatz, T., Gayler, V., Loew, A., and Claussen, M.: Evaluation of vegetation cover and land-surface albedo in MPI-ESM CMIP5 simulations, J. Adv. Model. Earth Syst., 5, 48–57, 2013. a

Ceppi, P., Hartmann, D. L., and Webb, M. J.: Mechanisms of the Negative Shortwave Cloud Feedback in Middle to High Latitudes, J. Climate, 29, 139–157, 2016. a

Cess, R. D., Potter, G. L., Blanchet, J. P., Boer, G. J., Ghan, S. J., Kiehl, J. T., Letreut, H., Li, Z. X., Liang, X. Z., Mitchell, J. F. B., Morcrette, J. J., Randall, D. A., Riches, M. R., Roeckner, E., Schlese, U., Slingo, A., Taylor, K. E., Washington, W. M., Wetherald, R. T., and Yagai, I.: Interpretation of Cloud-Climate Feedback as Produced by 14 Atmospheric General Circulation Models, Science, 245, 513–516, 1989. a

Collins, M., Knutti, R., Arblaster, J., Dufresne, J.-L., Fichefet, T., Friedlingstein, P., Gao, X., Gutowski, W. J., Johns, T., Krinner, G., Shongwe, M., Tebaldi, C., Weaver, A. J., and Wehner, M.: Long-term Climate Change: Projections, Commitments and Irreversibility, in: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Stocker, T., Qin, D., Plattner, G.-K., Tignor, M., Allen, S. K., Boschung, J., Nauels, A., Xia, Y., Bex, V., and Midgley, P. M., 1029–1136, Cambridge Univ. Press, Cambridge, UK and New York, NY, USA, 2013. a, b

Croft, B., Lohmann, U., Martin, R. V., Stier, P., Wurzler, S., Feichter, J., Posselt, R., and Ferrachat, S.: Aerosol size-dependent below-cloud scavenging by rain and snow in the ECHAM5-HAM, Atmos. Chem. Phys., 9, 4653–4675, https://doi.org/10.5194/acp-9-4653-2009, 2009. a

Croft, B., Lohmann, U., Martin, R. V., Stier, P., Wurzler, S., Feichter, J., Hoose, C., Heikkilä, U., van Donkelaar, A., and Ferrachat, S.: Influences of in-cloud aerosol scavenging parameterizations on aerosol concentrations and wet deposition in ECHAM5-HAM, Atmos. Chem. Phys., 10, 1511–1543, https://doi.org/10.5194/acp-10-1511-2010, 2010. a

Cziczo, D. J., Froyd, K. D., Hoose, C., Jensen, E. J., Diao, M. H., Zondlo, M. A., Smith, J. B., Twohy, C. H., and Murphy, D. M.: Clarifying the Dominant Sources and Mechanisms of Cirrus Cloud Formation, Science, 340, 1320–1324, 2013. a

Dufresne, J.-L. and Bony, S.: An assessment of the primary sources of spread of global warming estimates from coupled atmosphere-ocean models, J. Climate, 21, 5135–5144, 2008. a

Elsaesser, G. W., O'Dell, C. W., Lebsock, M. D., Bennartz, R., Greenwald, T. J., and Wentz, F. J.: The Multi-Sensor Advanced Climatology of Liquid Water Path (MAC-LWP), J. Climate, 30, 10193–10210, https://doi.org/10.1175/JCLI-D-16-0902.1, 2018. a, b, c, d, e, f, g, h, i

Flato, G., Marotzke, J., Abiodun, B., Braconnot, P., Chou, S. C., Collins, W., Cox, P., Driouech, F., Emori, S., Eyring, V., Forest, C., Gleckler, P., Guilyardi, E., Jakob, C., Kattsov, V., Reason, C., and Rummukainen, M.: Evaluation of Climate Models, in: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Stocker, T., Qin, D., Plattner, G.-K., Tignor, M., Allen, S. K., Boschung, J., Nauels, A., Xia, Y., Bex, V., and Midgley, P. M.,741–866, Cambridge Univ. Press, Cambridge, UK and New York, NY, USA, 2013. a, b, c, d, e, f, g

Forster, P. M., Andrews, T., Good, P., Gregory, J. M., Jackson, L. S., and Zelinka, M.: Evaluating adjusted forcing and model spread for historical and future scenarios in the CMIP5 generation of climate models, J. Geophys. Res., 118, 1139–1150, 2013. a, b, c

Frey, W. R. and Kay, J. E.: The influence of extratropical cloud phase and amount feedbacks on climate sensitivity, Clim. Dynam., 50, 3097–3116, https://doi.org/10.1007/s00382-017-3796-5, 2017. a, b

Frey, W. R., Maroon, E. A., Pendergrass, A. G., and Kay, J. E.: Do Southern Ocean Cloud Feedbacks Matter for 21st Century Warming?, Geophys. Res. Lett., 44, 12447–12456, 2017. a, b

Gasparini, B., Meyer, A., Neubauer, D., Münch, S., and Lohmann, U.: Cirrus cloud properties as seen by the CALIPSO satellite and ECHAM-HAM global climate model, J. Climate, 31, 1983–2003, https://doi.org/10.1175/JCLI-D-16-0608.1, 2018. a, b

Gettelman, A. and Sherwood, S. C.: Processes Responsible for Cloud Feedback, Current Climate Change Reports, 2, 179–189, 2016. a

Gordon, N. D. and Klein, S. A.: Low-cloud optical depth feedback in climate models, J. Geophys. Res., 119, 6052–6065, 2014. a

Gregory, J. M., Ingram, W. J., Palmer, M. A., Jones, G. S., Stott, P. A., Thorpe, R. B., Lowe, J. A., Johns, T. C., and Williams, K. D.: A new method for diagnosing radiative forcing and climate sensitivity, Geophys. Res. Lett., 31, L03205, https://doi.org/10.1029/2003GL018747, 2004. a, b

Hartmann, D. L. and Larson, K.: An important constaint on tropical cloud – climate feedback, Geophys. Res. Lett., 29, 1951, https://doi.org/10.1029/2002GL015835, 2002. a, b, c, d

Hoose, C., Lohmann, U., Bennartz, R., Croft, B., and Lesins, G.: Global simulations of aerosol processing in clouds, Atmos. Chem. Phys., 8, 6939–6963, https://doi.org/10.5194/acp-8-6939-2008, 2008. a

Hoose, C., Kristjánsson, J. E., Iversen, T., Kirkevåg, A., Seland, Ø., and Gettelman, A.: Constraining cloud droplet number concentration in GCMs suppresses the aerosol indirect effect, Geophys. Res. Lett., 36, L12807, https://doi.org/10.1029/2009GL038568, 2009. a

Iacono, M. J., Delamere, J. S., Mlawer, E. J., Shephard, M. W., Clough, S. A., and Collins, W. D.: Radiative forcing by long-lived greenhouse gases: Calculations with the AER radiative transfer models, J. Geophys. Res., 113, D13103, https://doi.org/10.1029/2008jd009944, 2008. a

Kärcher, B.: Cirrus Clouds and Their Response to Anthropogenic Activities, Current Climate Change Reports, 3, 45–57, 2017. a

Kay, J. E., Wall, C., Yettella, V., Medeiros, B., Hannay, C., Caldwell, P., and Bitz, C.: Global Climate Impacts of Fixing the Southern Ocean Shortwave Radiation Bias in the Community Earth System Model (CESM), J. Climate, 29, 4617–4636, 2016. a, b, c, d

Kinne, S., O'Donnel, D., Stier, P., Kloster, S., Zhang, K., Schmidt, H., Rast, S., Giorgetta, M., Eck, T., and Stevens, B.: MAC-v1: A new global aerosol climatology for climate studies, J. Adv. Model. Earth Sy., 5, 704–740, 2013. a

Knutti, R. and Rugenstein, M. A. A.: Feedbacks, climate sensitivity and the limits of linear models, Philos. T. R. Soc. A, 373, 2054, https://doi.org/10.1098/rsta.2015.0146, 2015. a

Komurcu, M., Storelvmo, T., Tan, I., Lohmann, U., Yun, Y., Penner, J. E., Wang, Y., Liu, X., and Takemura, T.: Intercomparison of the cloud water phase among global climate models, J. Geophys. Res., 119, 3372–3400, 2014. a, b, c, d

Korolev, A.: Limitations of the Wegener-Bergeron-Findeisen mechanism in the evolution of mixed-phase clouds, J. Atmos. Sci., 64, 3372–3375, 2007. a

Krämer, M., Rolf, C., Luebke, A., Afchine, A., Spelten, N., Costa, A., Meyer, J., Zöger, M., Smith, J., Herman, R. L., Buchholz, B., Ebert, V., Baumgardner, D., Borrmann, S., Klingebiel, M., and Avallone, L.: A microphysics guide to cirrus clouds – Part 1: Cirrus types, Atmos. Chem. Phys., 16, 3463–3483, https://doi.org/10.5194/acp-16-3463-2016, 2016. a

Lebsock, M. and Su, H.: Application of active spaceborne remote sensing for understanding biases between passive cloud water path retrievals, J. Geophys. Res., 119, 8962–8979, 2014. a, b

Li, J.-L. F., Waliser, D. E., Chen, W.-T., Guan, B., Kubar, T., Stephens, G., Ma, H.-Y., Deng, M., Donner, L., Seman, C., and Horowitz, L.: An observationally based evaluation of cloud ice water in CMIP3 and CMIP5 GCMs and contemporary reanalyses using contemporary satellite data, J. Geophys. Res., 117, D16105, https://doi.org/10.1029/2012JD017640, 2012. a, b, c, d, e

Lin, H. and Leaitch, W. R.: Development of an in-cloud aerosol activation parameterization for climate modelling, in: Proceedings of the WMO Workshop on Measurement of Cloud Properties for Forecasts of Weather, Air Quality and Climate, 328–335, World Meteorol. Organ., Geneva, 1997. a

Lin, S. J. and Rood, R. B.: Multidimensional flux-form semi-Lagrangian transport schemes, Mon. Weather Rev., 124, 2046–2070, 1996. a

Loeb, N. G., Wielicki, B. A., Doelling, D. R., Smith, G. L., Keyes, D. F., Kato, S., Manalo-Smith, N., and Wong, T.: Toward Optimal Closure of the Earth's Top-of-Atmosphere Radiation Budget, J. Climate, 22, 748–766, 2009. a, b, c, d

Loeb, N. G., Doelling, D. R., Wang, H. L., Su, W. Y., Nguyen, C., Corbett, J. G., Liang, L. S., Mitrescu, C., Rose, F. G., and Kato, S.: Clouds and the Earth's Radiant Energy System (CERES) Energy Balanced and Filled (EBAF) Top-of-Atmosphere (TOA) Edition-4.0 Data Product, J. Climate, 31, 895–918, 2018. a, b, c, d

Lohmann, U.: Sensitivität des Modellklimas eines globalen Zirkulationsmodells der Atmosphäre gegenüber Änderungen der Wolkenmikrophysik, PhD thesis, Max Planck Institute for Meteorology, Hamburg, Germary, 1996. a

Lohmann, U.: Possible aerosol effects on ice clouds via contact nucleation, J. Atmos. Sci., 59, 647–656, 2002. a

Lohmann, U.: Global anthropogenic aerosol effects on convective clouds in ECHAM5-HAM, Atmos. Chem. Phys., 8, 2115–2131, https://doi.org/10.5194/acp-8-2115-2008, 2008. a

Lohmann, U. and Ferrachat, S.: Impact of parametric uncertainties on the present-day climate and on the anthropogenic aerosol effect, Atmos. Chem. Phys., 10, 11373–11383, https://doi.org/10.5194/acp-10-11373-2010, 2010. a

Lohmann, U. and Hoose, C.: Sensitivity studies of different aerosol indirect effects in mixed-phase clouds, Atmos. Chem. Phys., 9, 8917–8934, https://doi.org/10.5194/acp-9-8917-2009, 2009. a, b, c, d

Lohmann, U. and Kärcher, B.: First interactive simulations of cirrus clouds formed by homogeneous freezing in the ECHAM GCM, J. Geophys. Res., 107, AAC 8-1–AAC 8-13, https://doi.org/10.1029/2001JD000767, 2002. a, b

Lohmann, U. and Roeckner, E.: The influence of cirrus cloud-radiative forcing on climate and climate sensitivity in a general circulation model, J. Geophys. Res., 100, 16305–16323, 1995. a, b

Lohmann, U. and Roeckner, E.: Design and performance of a new cloud microphysics scheme developed for the ECHAM general circulation model, Clim. Dynam., 12, 557–572, 1996. a

Lohmann, U., Stier, P., Hoose, C., Ferrachat, S., Kloster, S., Roeckner, E., and Zhang, J.: Cloud microphysics and aerosol indirect effects in the global climate model ECHAM5-HAM, Atmos. Chem. Phys., 7, 3425–3446, https://doi.org/10.5194/acp-7-3425-2007, 2007. a, b, c

Marcolli, C., Gedamke, S., Peter, T., and Zobrist, B.: Efficiency of immersion mode ice nucleation on surrogates of mineral dust, Atmos. Chem. Phys., 7, 5081–5091, https://doi.org/10.5194/acp-7-5081-2007, 2007. a

Matus, A. V. and L'Ecuyer, T. S.: The role of cloud phase in Earth's radiation budget, J. Geophys. Res., 122, 2559–2578, 2017. a, b, c

Mauritsen, T. and Stevens, B.: Missing iris effect as a possible cause of muted hydrological change and high climate sensitivity in models, Nat. Geosci., 8, 346–351, 2015. a, b

McCoy, D. T., Hartmann, D. L., Zelinka, M. D., Ceppi, P., and Grosvenor, D. P.: Mixed-phase cloud physics and Southern Ocean cloud feedback in climate models, J. Geophys. Res., 120, 9539–9554, 2015. a

Meehl, G. A., Stocker, T. F., Collins, W. D., Friedlingstein, P., Gaye, A. T., Gregory, J. M., Kitoh, A., Knutti, R., Murphy, J. M., Noda, A., Raper, S. C. B., Watterson, I. G., Weaver, A. J., and Zhao, Z.-C.: Global climate projections, in: Climate Change 2007: The Physical Science Basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Solomon, S., Qin, D., Manning, M., Chen, Z., Marquis, M., Averyt, K. B., Tignor, M., and Miller, H. L., 747–846, Cambridge Univ. Press, Cambridge, UK and New York, NY, USA, 2007. a

Mitchell, J. F. B., Senior, C. A., and Ingram, W. J.: CO2 and Climate – A Missing Feedback, Nature, 341, 132–134, 1989. a

Neubauer, D., Lohmann, U., Hoose, C., and Frontoso, M. G.: Impact of the representation of marine stratocumulus clouds on the anthropogenic aerosol effect, Atmos. Chem. Phys., 14, 11997–12022, https://doi.org/10.5194/acp-14-11997-2014, 2014. a, b

Pincus, R. and Stevens, B.: Paths to accuracy for radiation parameterizations in atmospheric models, J. Adv. Model. Earth Sy., 5, 225–233, 2013. a, b

Platnick, S., King, M. D., Meyer, K. G., Wind, G., Amarasinghe, N., Marchant, B., Arnold, G. T., Zhang, Z. B., Hubanks, P. A., Ridgway, B., and Riedi, J.: MODIS Cloud Optical Properties: User Guide for the Collection 6 Level-2 MOD06/MYD06 Product and Associated Level-3 Datasets, NASA, available at: https://modis-atmos.gsfc.nasa.gov/_docs/C6MOD06OPUserGuide.pdf (last access: 20 June 2018), 2015. a, b, c, d, e

Platnick, S., Meyer, K. G., King, M. D., Wind, G., Amarasinghe, N., Marchant, B., Arnold, G. T., Zhang, Z., Hubanks, P. A., Holz, R. E., Yang, P., Ridgway, W. L., and Riedi, J.: The MODIS Cloud Optical and Microphysical Products: Collection 6 Updates and Examples From Terra and Aqua, IEEE T. Geosci. Remote Sens., 55, 502–525, 2017. a, b, c, d, e

Poulsen, C., McGarragh, G., Thomas, G., Christensen, M., Povey, A., Grainger, D., Proud, S., and Hollmann, R.: ESA Cloud Climate Change Initiative (ESA Cloud_cci) data: Cloud_cci ATSR2-AATSR L3C/L3U/L2 CLD_PRODUCTS v2.0, Deutscher Wetterdienst (DWD) and Rutherford Appleton Laboratory (Dataset Producer), 2017. a, b, c, d, e

Pruppacher, H. R. and Klett, J. D.: Microphysics of Clouds and Precipitation, Kluwer Acad., Norwell, Mass., 1997. a

Randall, D., Wood, R., Bony, S., Colman, R., Fichefet, T., Fyfe, J., Kattsov, V., Pitman, A., Shukla, J., Srinivasan, J., Stouffer, R., Sumi, A., and Tayler, K.: Climate Models and Their Evaluation, in: Climate Change 2007: The Scientific Basis. Contribution of working group I to the Forth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Solomon, S., Qin, D., Manning, M., Chen, Z., Marquis, M., Averyt, K. B., Tignor, M., and Miller, H. L., 589–662, Cambridge Univ. Press, Cambridge, UK and New York, NY, USA, 2007. a

Reick, C. H., Raddatz, T., Brovkin, V., and Gayler, V.: Representation of natural and anthropogenic land cover change in MPI-ESM, J. Adv. Model. Earth Sy., 5, 459–482, 2013. a

Schultz, M. G., Stadtler, S., Schröder, S., Taraborrelli, D., Franco, B., Krefting, J., Henrot, A., Ferrachat, S., Lohmann, U., Neubauer, D., Siegenthaler-Le Drian, C., Wahl, S., Kokkola, H., Kühn, T., Rast, S., Schmidt, H., Stier, P., Kinnison, D., Tyndall, G. S., Orlando, J. J., and Wespes, C.: The chemistry-climate model ECHAM6.3-HAM2.3-MOZ1.0, Geosci. Model Dev., 11, 1695–1723, https://doi.org/10.5194/gmd-11-1695-2018, 2018. a, b

Seethala, C. and Horvath, A.: Global assessment of AMSR-E and MODIS cloud liquid water path retrievals in warm oceanic clouds, J. Geophys. Res., 115, D13202, https://doi.org/10.1029/2009JD012662, 2010. a

Seifert, A. and Beheng, K. D.: A two-moment cloud microphysics parameterization for mixed-phase clouds. Part 1: Model description, Meteorol. Atmos. Phys., 92, 45–66, 2006. a

Shell, K. M., Kiehl, J. T., and Shields, C. A.: Using the radiative kernel technique to calculate climate feedbacks in NCAR's Community Atmospheric Model, J. Climate, 21, 2269–2282, 2008. a, b, c, d

Soden, B. J., Wetherald, R. T., Stenchikov, G. L., and Robock, A.: Clobal cooling after the eruption of Mount Pinatubo: A test of climate feedback by water vapor, Science, 296, 727–730, 2002. a

Soden, B. J., Held, I. M., Colman, R., Shell, K. M., Kiehl, J. T., and Shields, C. A.: Quantifying Climate Feedbacks Using Radiative Kernels, J. Climate, 21, 3504–3520, 2008. a

Solomon, S., Qin, D., Manning, M., Alley, R. B., Berntsen, T., Bindoff, N. L., Chen, Z., Chidthaisong, A., Gregory, J. M., Hegerl, G. C., Heimann, H., Hewitson, B., Hoskins, B. J., Joos, F., Jouzel, J., Kattsov, V., Lohmann, U., Matsuno, T., Molina, M., Nicholls, N., Overpeck, J., Raga, G., Ramaswamy, V., Ren, J., Rusticucci, M., Somerville, R., Stocker, T. F., Whetton, P., A., W. R., and Wratt, D.: Technical Summary, in: Climate Change 2007: The Physical Science Basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Solomon, S., Qin, D., Manning, M., Chen, Z., Marquis, M., Averyt, K. B., Tignor, M., and Miller, H. L., Cambridge Univ. Press, Cambridge, UK and New York, NY, USA, 2007. a

Spichtinger, P. and Krämer, M.: Tropical tropopause ice clouds: a dynamic approach to the mystery of low crystal numbers, Atmos. Chem. Phys., 13, 9801–9818, https://doi.org/10.5194/acp-13-9801-2013, 2013. a

Stengel, M., Stapelberg, S., Sus, O., Schlundt, C., Poulsen, C., Thomas, G., Christensen, M., Carbajal Henken, C., Preusker, R., Fischer, J., Devasthale, A., Willén, U., Karlsson, K.-G., McGarragh, G. R., Proud, S., Povey, A. C., Grainger, R. G., Meirink, J. F., Feofilov, A., Bennartz, R., Bojanowski, J. S., and Hollmann, R.: Cloud property datasets retrieved from AVHRR, MODIS, AATSR and MERIS in the framework of the Cloud_cci project, Earth Syst. Sci. Data, 9, 881–904, https://doi.org/10.5194/essd-9-881-2017, 2017a. a, b, c, d, e

Stengel, M., Sus, O., Stapelberg, S., Schlundt, C., Poulsen, C., and Hollmann, R.: ESA Cloud Climate Change Initiative (ESA Cloud_cci), Deutscher Wetterdienst (DWD), Offenbach, Germany, 2017b. a, b

Stephens, G. L., Li, J., Wild, M., Clayson, C. A., Loeb, N., Kato, S., L'Ecuyer, T., Stackhouse Jr., P. W., Lebsock, M., and Andrews, T.: An update on Earth's energy balance in light of the latest global observations, Nat. Geosci., 5, 691–696, https://doi.org/10.1038/ngeo1580, 2012. a

Stevens, B., Giorgetta, M., Esch, M., Mauritsen, T., Crueger, T., Rast, S., Salzmann, M., Schmidt, H., Bader, J., Block, K., Brokopf, R., Fast, I., Kinne, S., Kornblueh, L., Lohmann, U., Pincus, R., Reichler, T., and Roeckner, E.: Atmospheric component of the MPI-M Earth System Model: ECHAM6, J. Adv. Model. Earth Sy., 5, 146–172, 2013. a, b, c, d

Stevens, B., Fiedler, S., Kinne, S., Peters, K., Rast, S., Müsse, J., Smith, S. J., and Mauritsen, T.: MACv2-SP: a parameterization of anthropogenic aerosol optical properties and an associated Twomey effect for use in CMIP6, Geosci. Model Dev., 10, 433–452, https://doi.org/10.5194/gmd-10-433-2017, 2017. a

Stier, P.: Limitations of passive remote sensing to constrain global cloud condensation nuclei, Atmos. Chem. Phys., 16, 6595–6607, https://doi.org/10.5194/acp-16-6595-2016, 2016. a

Stier, P., Feichter, J., Kinne, S., Kloster, S., Vignati, E., Wilson, J., Ganzeveld, L., Tegen, I., Werner, M., Balkanski, Y., Schulz, M., Boucher, O., Minikin, A., and Petzold, A.: The aerosol-climate model ECHAM5-HAM, Atmos. Chem. Phys., 5, 1125–1156, https://doi.org/10.5194/acp-5-1125-2005, 2005. a

Stubenrauch, C. J., Rossow, W. B., Kinne, S., Ackerman, S., Cesana, G., Chepfer, H., Girolamo, L. D., Getzewich, B., Guignard, A., Heidinger, A., Maddux, B. C., Menzel, W. P., Minnis, P., Pearl, C., Platnick, S., Poulsen, C., Riedi, J., Sun-Mack, S., Walther, A., Winker, D., Zeng, S., and Zhao, G.: Assessment of Global Cloud Datasets from Satellites: Project and Database Initiated by the GEWEX Radiation Panel, B. Am. Meteorol. Soc., 94, 1031–1049, 2013. a, b, c

Sundqvist, H., Berge, E., and Kristjánsson, J. E.: Condensation and Cloud Parameterization Studies with a Mesoscale Numerical Weather Prediction Model, Mon. Weather Rev., 117, 1641–1657, 1989. a

Susskind, J., Piraino, P., Rokke, L., Iredell, T., and Mehta, A.: Characteristics of the TOVS Pathfinder Path A dataset, B. Am. Meteorol. Soc., 78, 1449–1472, 1997. a, b, c

Tan, I., Storelvmo, T., and Zelinka, M. D.: Observational constraints on mixed-phase clouds imply higher climate sensitivity, Science, 352, 224–227, 2016. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o

Terai, C. R., Bretherton, C. S., Wood, R., and Painter, G.: Aircraft observations of aerosol, cloud, precipitation, and boundary layer properties in pockets of open cells over the southeast Pacific, Atmos. Chem. Phys., 14, 8071–8088, https://doi.org/10.5194/acp-14-8071-2014, 2014. a

Terai, C. R., Klein, S. A., and Zelinka, M. D.: Constraining the low-cloud optical depth feedback at middle and high latitudes using satellite observations, J. Geophys. Res., 121, 9696–9716, 2016. a

Tiedtke, M.: A Comprehensive Mass Flux Scheme for Cumulus Parameterization in Large-Scale Models, Mon. Weather Rev., 117, 3040–3061, 1989. a

Vial, J., Dufresne, J.-L., and Bony, S.: On the interpretation of inter-model spread in CMIP5 climate sensitivity estimates, Clim. Dynam., 41, 3339–3362, 2013. a

Webb, M. J., Lock, A. P., Bretherton, C. S., Bony, S., Cole, J. N. S., Idelkadi, A., Kang, S. M., Koshiro, T., Kawai, H., Ogura, T., Roehrig, R., Shin, Y., Mauritsen, T., Sherwood, S. C., Vial, J., Watanabe, M., Woelfle, M. D., and Zhao, M.: The impact of parametrized convection on cloud feedback, Philos. T. R. Soc. A, 373, 2054, https://doi.org/10.1098/rsta.2014.0414, 2015. a

Wernli, H., Boettcher, M., Joos, H., Miltenberger, A. K., and Spichtinger, P.: A trajectory-based classification of ERA-Interim ice clouds in the region of the North Atlantic storm track, Geophys. Res. Lett., 43, 6657–6664, 2016. a

Williams, K. D., Bodas-Salcedo, A., Déqué, M., Fermepin, S., Medeiros, B., Watanabe, M., Jakob, C., Klein, S. A., Senior, C. A., and Williamson, D. L.: The Transpose-AMIP II Experiment and Its Application to the Understanding of Southern Ocean Cloud Biases in Climate Models, J. Climate, 26, 3258–3274, 2013. a

Wood, R.: Drizzle in stratocumulus: Observations and model results, in: 11th AMS conference on cloud physics, Amer. Meteorol. Soc., Ogden, Utah, 2002. a

Wood, R., O, K.-T., Bretherton, C. S., Mohrmann, J., Albrecht, B. A., Zuidema, P., Ghate, V., Schwartz, C., Eloranta, E., Glienke, S., Shaw, R. A., Fugal, J., and Minnis, P.: Ultraclean Layers and Optically Thin Clouds in the Stratocumulus-to-Cumulus Transition. Part I: Observations, J. Atmos. Sci., 75, 1631–1652, https://doi.org/10.1175/JAS-D-17-0213.1, 2018.  a

Yokohata, T., Webb, M. J., Collins, M., Williams, K. D., Yoshimori, M., Hargreaves, J. C., and Annan, J. D.: Structural Similarities and Differences in Climate Responses to CO2 Increase between Two Perturbed Physics Ensembles, J. Climate, 23, 1392–1410, https://doi.org/10.1175/2009JCLI2917.1, 2010. a

Zelinka, M. D., Klein, S. A., and Hartmann, D. L.: Computing and Partitioning Cloud Feedbacks Using Cloud Property Histograms. Part I: Cloud Radiative Kernels, J. Climate, 25, 3715–3735, https://doi.org/10.1175/JCLI-D-11-00248.1, 2012a. a

Zelinka, M. D., Klein, S. A., and Hartmann, D. L.: Computing and Partitioning Cloud Feedbacks Using Cloud Property Histograms. Part II: Attribution to Changes in Cloud Amount, Altitude, and Optical Depth, J. Climate, 25, 3736–3754, 2012b. a

Zelinka, M. D., Zhou, C., and Klein, S. A.: Insights from a refined decomposition of cloud feedbacks, Geophys. Res. Lett., 43, 9259–9269, 2016. a

Zelinka, M. D., Randall, D. A., Webb, M. J., and Klein, S. A.: Clearing clouds of uncertainty, Nat. Clim. Change, 7, 674–678, https://doi.org/10.1038/nclimate3402, 2017. a, b

Zhang, K., O'Donnell, D., Kazil, J., Stier, P., Kinne, S., Lohmann, U., Ferrachat, S., Croft, B., Quaas, J., Wan, H., Rast, S., and Feichter, J.: The global aerosol-climate model ECHAM-HAM, version 2: sensitivity to improvements in process representations, Atmos. Chem. Phys., 12, 8911–8949, https://doi.org/10.5194/acp-12-8911-2012, 2012. a, b, c, d, e

Zhou, C., Dessler, A. E., Zelinka, M. D., Yang, P., and Wang, T.: Cirrus feedback on inter-annual climate fluctuations, Geophys. Res. Lett., 41, 9166–9173, 2014. a

Short summary
The climate is warming, at the current rate so much so that the 2 ºC target is likely to be exceeded. Uncertainty remains as to when the 2 ºC of warming will be reached. One factor contributing to this uncertainty is how clouds change in the warmer climate. While previously most emphasis was placed on how low clouds change in the warmer climate, here we investigate the importance of mixed-phase and ice clouds.
The climate is warming, at the current rate so much so that the 2 ºC target is likely to be...
Citation