skip to main content

Title: LGM Paleoclimate Constraints Inform Cloud Parameterizations and Equilibrium Climate Sensitivity in CESM2

The Community Earth System Model version 2 (CESM2) simulates a high equilibrium climate sensitivity (ECS > 5°C) and a Last Glacial Maximum (LGM) that is substantially colder than proxy temperatures. In this study, we examine the role of cloud parameterizations in simulating the LGM cooling in CESM2. Through substituting different versions of cloud schemes in the atmosphere model, we attribute the excessive LGM cooling to the new CESM2 schemes of cloud microphysics and ice nucleation. Further exploration suggests that removing an inappropriate limiter on cloud ice number (NoNimax) and decreasing the time‐step size (substepping) in cloud microphysics largely eliminate the excessive LGM cooling. NoNimax produces a more physically consistent treatment of mixed‐phase clouds, which leads to an increase in cloud ice content and a weaker shortwave cloud feedback over mid‐to‐high latitudes and the Southern Hemisphere subtropics. Microphysical substepping further weakens the shortwave cloud feedback. Based on NoNimax and microphysical substepping, we have developed a paleoclimate‐calibrated CESM2 (PaleoCalibr), which simulates well the observed twentieth century warming and spatial characteristics of key cloud and climate variables. PaleoCalibr has a lower ECS (∼4°C) and a 20% weaker aerosol‐cloud interaction than CESM2. PaleoCalibr represents a physically more consistent treatment of cloud microphysics than more » CESM2 and is a valuable tool in climate change studies, especially when a large climate forcing is involved. Our study highlights the unique value of paleoclimate constraints in informing the cloud parameterizations and ultimately the future climate projection.

« less
 ;  ;  ;  ;  ;  ;  ;  ;  ;  
Award ID(s):
Publication Date:
Journal Name:
Journal of Advances in Modeling Earth Systems
DOI PREFIX: 10.1029
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The upper end of the equilibrium climate sensitivity (ECS) has increased substantially in the latest Coupled Model Intercomparison Projects phase 6 with eight models (as of this writing) reporting an ECS > 5°C. The Community Earth System Model version 2 (CESM2) is one such high‐ECS model. Here we perform paleoclimate simulations of the Last Glacial Maximum (LGM) using CESM2 to examine whether its high ECS is realistic. We find that the simulated LGM global mean temperature decrease exceeds 11°C, greater than both the cooling estimated from proxies and simulated by an earlier model version (CESM1). The large LGM cooling in CESM2 is attributed to a strong shortwave cloud feedback in the newest atmosphere model. Our results indicate that the high ECS of CESM2 is incompatible with LGM constraints and that the projected future warming in CESM2, and models with a similarly high ECS, is thus likely too large.

  2. Abstract

    Global climate models (GCMs) are challenged by difficulties in simulating cloud phase and cloud radiative effect over the Southern Ocean (SO). Some of the new‐generation GCMs predict too much liquid and too little ice in mixed‐phase clouds. This misrepresentation of cloud phase in GCMs results in weaker negative cloud feedback over the SO and a higher climate sensitivity. Based on a model comparison with observational data obtained during the Southern Ocean Cloud Radiation and Aerosol Transport Experimental Study, this study addresses a key uncertainty in the Community Earth System Model version 2 (CESM2) related to cloud phase, namely ice formation in pristine remote SO clouds. It is found that sea spray organic aerosols (SSOAs) are the most important type of ice nucleating particles (INPs) over the SO with concentrations 1 order of magnitude higher than those of dust INPs based on measurements and CESM2 simulations. Secondary ice production (SIP) which includes riming splintering, rain droplet shattering, and ice‐ice collisional fragmentation as implemented in CESM2 is the dominant ice production process in moderately cold clouds with cloud temperatures greater than −20°C. SIP enhances the in‐cloud ice number concentrations (Ni) by 1–3 orders of magnitude and predicts more mixed‐phase (with percentagemore »occurrence increased from 15% to 21%), in better agreement with the observations. This study highlights the importance of accurately representing the cloud phase over the pristine remote SO by considering the ice nucleation of SSOA and SIP processes, which are currently missing in most GCM cloud microphysics parameterizations.

    « less
  3. Abstract. Equilibrium climate sensitivity (ECS) has been directly estimated using reconstructions of past climates that are different than today's. A challenge to this approach is that temperature proxies integrate over the timescales of the fast feedback processes (e.g., changes in water vapor, snow, and clouds) that are captured in ECS as well as the slower feedback processes (e.g., changes in ice sheets and ocean circulation) that are not. A way around this issue is to treat the slow feedbacks as climate forcings and independently account for their impact on global temperature. Here we conduct a suite of Last Glacial Maximum (LGM) simulations using the Community Earth System Model version 1.2 (CESM1.2) to quantify the forcingand efficacy of land ice sheets (LISs) and greenhouse gases (GHGs) in order to estimate ECS. Our forcing and efficacy quantification adopts the effective radiative forcing (ERF) and adjustment framework and provides a complete accounting for the radiative, topographic, and dynamical impacts of LIS on surface temperatures. ERF and efficacy of LGM LIS are −3.2 W m−2 and 1.1, respectively. The larger-than-unity efficacy is caused by the temperature changes over land and the Northern Hemisphere subtropical oceans which are relatively larger than those in response to a doubling of atmospheric CO2. The subtropical sea-surface temperature (SST) response ismore »linked to LIS-induced wind changes and feedbacks in ocean–atmosphere coupling and clouds. ERF and efficacy of LGM GHG are −2.8 W m−2 and 0.9, respectively. The lower efficacy is primarily attributed to a smaller cloud feedback at colder temperatures. Our simulations further demonstrate that the direct ECS calculation using the forcing, efficacy, and temperature response in CESM1.2 overestimates the true value in the model by approximately 25 % due to the neglect of slow ocean dynamical feedback. This is supported by the greater cooling (6.8 ∘C) in a fully coupled LGM simulation than that (5.3 ∘C) in a slab ocean model simulation with ocean dynamics disabled. The majority (67 %) of the ocean dynamical feedback is attributed to dynamical changes in the Southern Ocean, where interactions between upper-ocean stratification, heat transport, and sea-ice cover are found to amplify the LGM cooling. Our study demonstrates the value of climate models in the quantification of climate forcings and the ocean dynamical feedback, which is necessary for an accurate direct ECS estimation.« less
  4. Abstract

    A framework is introduced to investigate the indirect effect of aerosol loading on tropical deep convection using three-dimensional limited-domain idealized cloud-system-resolving model simulations coupled with large-scale dynamics over fixed sea surface temperature. The large-scale circulation is parameterized using the spectral weak temperature gradient (WTG) approximation that utilizes the dominant balance between adiabatic cooling and diabatic heating in the tropics. The aerosol loading effect is examined by varying the number of cloud condensation nuclei (CCN) available to form cloud droplets in the two-moment bulk microphysics scheme over a wide range of environments from 30 to 5000 cm−3. The radiative heating is held at a constant prescribed rate in order to isolate the microphysical effects. Analyses are performed over the period after equilibrium is achieved between convection and the large-scale environment. Mean precipitation is found to decrease modestly and monotonically when the aerosol number concentration increases as convection gets weaker, despite the increase in cloud liquid water in the warm-rain region and ice crystals aloft. This reduction is traced down to the reduction in surface enthalpy fluxes as an energy source to the atmospheric column induced by the coupling of the large-scale motion, though the gross moist stability remains constant. Increasingmore »CCN concentration leads to 1) a cooler free troposphere because of a reduction in the diabatic heating and 2) a warmer boundary layer because of suppressed evaporative cooling. This dipole temperature structure is associated with anomalously descending large-scale vertical motion above the boundary layer and ascending motion at lower levels. Sensitivity tests suggest that changes in convection and mean precipitation are unlikely to be caused by the impact of aerosols on cloud droplets and microphysical properties but rather by accounting for the feedback from convective adjustment with the large-scale dynamics. Furthermore, a simple scaling argument is derived based on the vertically integrated moist static energy budget, which enables estimation of changes in precipitation given known changes in surfaces enthalpy fluxes and the constant gross moist stability. The impact on cloud hydrometeors and microphysical properties is also examined, and it is consistent with the macrophysical picture.

    « less
  5. Abstract

    Modeling and paleoclimate proxy‐based studies suggest that equilibrium climate sensitivity (ECS) depends on the background climate state, though the reason is not thoroughly understood. Here we study the state dependence of ECS over a large range of global mean surface temperature (GMST) in the Community Atmosphere Model (CAM) Versions 4, 5, and 6 by varying atmospheric CO2concentrations. We find a robust increase of ECS with GMST in all three models, albeit at different rates, which is primarily attributed to strengthening of the shortwave cloud feedback (λcld) at both high and low latitudes. Over high latitudes, increasing GMST leads to a reduction in the cloud ice fraction, weakening the (negative) cloud‐phase feedback due to the phase transition of cloud ice to liquid and thereby strengtheningλcld. Over low‐latitude regions, increasing GMST strengthensλcldlikely through the nonlinear increase in water vapor, which causes low‐cloud thinning through thermodynamic and radiative processes.