skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Multiple Equilibria in a Coupled Climate–Carbon Model
Abstract Multiple stable equilibria are intrinsic to many complex dynamical systems, and have been identified in a hierarchy of climate models. Motivated by the idea that the Quaternary glacial–interglacial cycles could have resulted from orbitally forced transitions between multiple stable states mediated by internal feedbacks, this study investigates the existence and mechanisms of multiple equilibria in an idealized, energy-conserving atmosphere–ocean–sea ice general circulation model with a fully coupled carbon cycle. Four stable climates are found for identical insolation and global carbon inventory: an ice-free Warm climate, two intermediate climates (Cold and Waterbelt), and a fully ice-covered Snowball climate. A fifth state, a small ice cap state between Warm and Cold, is found to be barely unstable. Using custom radiative kernels and a thorough sampling of the model’s internal variability, three equilibria are investigated through the state dependence of radiative feedback processes. For fast feedbacks, the systematic decrease in surface albedo feedback from Cold to Warm states is offset by a similar increase in longwave water vapor feedback. At longer time scales, the key role of the carbon cycle is a dramatic lengthening of the adjustment time comparable to orbital forcings near the Warm state. The dynamics of the coupled climate–carbon system are thus not well separated in time from orbital forcings, raising interesting possibilities for nonlinear triggers for large climate changes. Significance Statement How do carbon cycle and other physical processes affect the physical and mathematical properties of the climate system? We use a complex climate model coupled with a carbon cycle to simulate the climate evolution under different initial conditions. Four stable climate states are possible, from the Snowball Earth, in which ice covers the whole planet, to the Warm state, an ice-free world. The carbon cycle drives the global climate change at an extremely slower pace after sea ice retreats. Sea ice and water vapor, on the other hand, constitute the major contributing factors that accelerate faster climate change.  more » « less
Award ID(s):
1952745
PAR ID:
10447866
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Journal of Climate
Volume:
36
Issue:
2
ISSN:
0894-8755
Page Range / eLocation ID:
547 to 564
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Paleoclimate records have been used to estimate the modern equilibrium climate sensitivity. However, this requires understanding how the feedbacks governing the climate response vary with the climate itself. Here we warm and cool a state-of-the-art climate model to simulate a continuum of climates ranging from a nearly ice-covered Snowball Earth to a nearly ice-free hothouse. We find that the pre-industrial (PI) climate is near a stability optimum: warming leads to a less-stable (more-sensitive) climate, as does cooling of more than 2K. Physically interpreting the results, we find that the decrease in stability for climates colder than the PI occurs mainly due to the albedo and lapse-rate feedbacks, and the decrease in stability for warmer climates occurs mainly due to the cloud feedback. These results imply that paleoclimate records provide a stronger constraint than has been calculated in previous studies, suggesting a reduction in the uncertainty range of the climate sensitivity. 
    more » « less
  2. Abstract The ice–albedo feedback associated with sea ice loss contributes to polar amplification, while the water vapor feedback contributes to tropical amplification of surface warming. However, these feedbacks are not independent of atmospheric energy transport, raising the possibility of complex interactions that may obscure the drivers of polar amplification, in particular its manifestation across the seasonal cycle. Here, we apply a radiative transfer hierarchy to an idealized aquaplanet global climate model coupled to a thermodynamic sea ice model. The climate responses and radiative feedbacks are decomposed into the contributions from sea ice loss, including both retreat and thinning, and the radiative effect of water vapor changes. We find that summer sea ice retreat causes winter polar amplification through ocean heat uptake and release, and the resulting decrease in dry energy transport weakens the magnitude of warming. Moreover, sea ice thinning is found to suppress summer warming and enhance winter warming, additionally contributing to winter amplification. The water vapor radiative effect produces seasonally symmetric polar warming via offsetting effects: enhanced moisture in the summer hemisphere induces the summer water vapor feedback and simultaneously strengthens the winter latent energy transport in the winter hemisphere by increasing the meridional moisture gradient. These results reveal the importance of changes in atmospheric energy transport induced by sea ice retreat and increased water vapor to seasonal polar amplification, elucidating the interactions among these physical processes. 
    more » « less
  3. Abstract During the Middle Miocene Climate Transition (MMCT; ∼14.7–13.8 Ma), the global climate experienced rapid cooling, leading to modern‐like temperatures, precipitation patterns, and permanent ice sheets. However, proxy records indicate that atmospheric pCO2and regional climate conditions (SST, ice volume) were highly variable from 17 to 12.5 Ma and these changes were not always synchronous. Here, we report on a series of middle Miocene (∼16–12.5 Ma) simulations using the water isotope enabled earth system model (iCESM1.2) to explore the potential for multiple equilibrium states to explain the observed decoupling between pCO2and regional climates. Our simulations indicate that initial ocean conditions can significantly influence deep water formation in the North Atlantic and lead to multiple ocean equilibria. When the model is initiated from a cold state, residual cool surface water temperatures in the North Atlantic intensify Atlantic Meridional Ocean Circulation (AMOC) and inhibit Arctic sea‐ice formation. When initiated from a warm state, the AMOC remains weak. The different ocean states drive differences in equator‐to‐pole sea surface temperature gradients and sea ice distributions through heat redistribution changes. These equilibria cause variations in temperature gradients and sea ice distribution due to changes in heat redistribution. Additionally, changes in ocean circulation and a reduced temperature gradient in the North Atlantic increase North Atlantic precipitation when the AMOC is strong. These findings underscore the importance of the ocean's initial state in shaping regional climate responses to atmospheric pCO2, potentially explaining regional climate pattern variability observed during the Miocene. 
    more » « less
  4. Abstract Radiative climate feedbacks in the Arctic have been extensively studied, but their spatial and seasonal variations have not been thoroughly examined. Using ERA5 reanalysis data, we examine seasonal variations in Arctic climate feedbacks and their relationship to sea‐ice loss based on changes from 1950–1979 to 1990–2019. The spring and summer seasons experienced large sea‐ice loss, strong surface albedo feedback, and large oceanic heat uptake. Arctic clouds exerted small net cooling in May‐June‐July but moderate warming during the cold season, especially over areas with large sea‐ice loss where cloud liquid and ice water content increased. Arctic water vapor feedback peaked in summer but was weak and uncorrelated with sea‐ice loss. Arctic positive lapse rate feedback (LRF) was strongest in winter over areas with large sea‐ice loss and weak inversion but uncorrelated with atmospheric stability, suggesting that oceanic heating from sea‐ice loss led to enhanced surface warming and the positive LRF. 
    more » « less
  5. null (Ed.)
    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 is 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. 
    more » « less