No Arabic abstract
Cenozoic temperature, sea level and CO2 co-variations provide insights into climate sensitivity to external forcings and sea level sensitivity to climate change. Climate sensitivity depends on the initial climate state, but potentially can be accurately inferred from precise paleoclimate data. Pleistocene climate oscillations yield a fast-feedback climate sensitivity 3 +/- 1{deg}C for 4 W/m2 CO2 forcing if Holocene warming relative to the Last Glacial Maximum (LGM) is used as calibration, but the error (uncertainty) is substantial and partly subjective because of poorly defined LGM global temperature and possible human influences in the Holocene. Glacial-to-interglacial climate change leading to the prior (Eemian) interglacial is less ambiguous and implies a sensitivity in the upper part of the above range, i.e., 3-4{deg}C for 4 W/m2 CO2 forcing. Slow feedbacks, especially change of ice sheet size and atmospheric CO2, amplify total Earth system sensitivity by an amount that depends on the time scale considered. Ice sheet response time is poorly defined, but we show that the slow response and hysteresis in prevailing ice sheet models are exaggerated. We use a global model, simplified to essential processes, to investigate state-dependence of climate sensitivity, finding an increased sensitivity towards warmer climates, as low cloud cover is diminished and increased water vapor elevates the tropopause. Burning all fossil fuels, we conclude, would make much of the planet uninhabitable by humans, thus calling into question strategies that emphasize adaptation to climate change.
Additional material supporting the article Target atmospheric CO2: Where should humanity aim?
A climate state close to a tipping point will have a degenerate linear response to perturbations, which can be associated with extreme values of the equilibrium climate sensitivity (ECS). In this paper we contrast linearized (`instantaneous) with fully nonlinear geometric (`two-point) notions of ECS, in both presence and absence of tipping points. For a stochastic energy balance model of the global mean surface temperature with two stable regimes, we confirm that tipping events cause the appearance of extremes in both notions of ECS. Moreover, multiple regimes with different mean sensitivities are visible in the two-point ECS. We confirm some of our findings in a physics-based multi-box model of the climate system.
Paleoclimate data show that climate sensitivity is ~3 deg-C for doubled CO2, including only fast feedback processes. Equilibrium sensitivity, including slower surface albedo feedbacks, is ~6 deg-C for doubled CO2 for the range of climate states between glacial conditions and ice-free Antarctica. Decreasing CO2 was the main cause of a cooling trend that began 50 million years ago, large scale glaciation occurring when CO2 fell to 450 +/- 100 ppm, a level that will be exceeded within decades, barring prompt policy changes. If humanity wishes to preserve a planet similar to that on which civilization developed and to which life on Earth is adapted, paleoclimate evidence and ongoing climate change suggest that CO2 will need to be reduced from its current 385 ppm to at most 350 ppm. The largest uncertainty in the target arises from possible changes of non-CO2 forcings. An initial 350 ppm CO2 target may be achievable by phasing out coal use except where CO2 is captured and adopting agricultural and forestry practices that sequester carbon. If the present overshoot of this target CO2 is not brief, there is a possibility of seeding irreversible catastrophic effects.
We use numerical climate simulations, paleoclimate data, and modern observations to study the effect of growing ice melt from Antarctica and Greenland. Meltwater tends to stabilize the ocean column, inducing amplifying feedbacks that increase subsurface ocean warming and ice shelf melting. Cold meltwater and induced dynamical effects cause ocean surface cooling in the Southern Ocean and North Atlantic, thus increasing Earths energy imbalance and heat flux into most of the global oceans surface. Southern Ocean surface cooling, while lower latitudes are warming, increases precipitation on the Southern Ocean, increasing ocean stratification, slowing deepwater formation, and increasing ice sheet mass loss. These feedbacks make ice sheets in contact with the ocean vulnerable to accelerating disintegration. We hypothesize that ice mass loss from the most vulnerable ice, sufficient to raise sea level several meters, is better approximated as exponential than by a more linear response. Doubling times of 10, 20 or 40 years yield multi-meter sea level rise in about 50, 100 or 200 years. Recent ice melt doubling times are near the lower end of the 10-40 year range, but the record is too short to confirm the nature of the response. The feedbacks, including subsurface ocean warming, help explain paleoclimate data and point to a dominant Southern Ocean role in controlling atmospheric CO2, which in turn exercised tight control on global temperature and sea level. The millennial (500-2000 year) time scale of deep ocean ventilation affects the time scale for natural CO2 change and thus the time scale for paleo global climate, ice sheet, and sea level changes, but this paleo millennial time scale should not be misinterpreted as the time scale for ice sheet response to a rapid large human-made climate forcing.
The stability properties of intermediate-order climate models are investigated by computing their Lyapunov exponents (LEs). The two models considered are PUMA (Portable University Model of the Atmosphere), a primitive-equation simple general circulation model, and MAOOAM (Modular Arbitrary-Order Ocean-Atmosphere Model), a quasi-geostrophic coupled ocean-atmosphere model on a beta-plane. We wish to investigate the effect of the different levels of filtering on the instabilities and dynamics of the atmospheric flows. Moreover, we assess the impact of the oceanic coupling, the dissipation scheme and the resolution on the spectra of LEs. The PUMA Lyapunov spectrum is computed for two different values of the meridional temperature gradient defining the Newtonian forcing. The increase of the gradient gives rise to a higher baroclinicity and stronger instabilities, corresponding to a larger dimension of the unstable manifold and a larger first LE. The convergence rate of the rate functional for the large deviation law of the finite-time Lyapunov exponents (FTLEs) is fast for all exponents, which can be interpreted as resulting from the absence of a clear-cut atmospheric time-scale separation in such a model. The MAOOAM spectra show that the dominant atmospheric instability is correctly represented even at low resolutions. However, the dynamics of the central manifold, which is mostly associated to the ocean dynamics, is not fully resolved because of its associated long time scales, even at intermediate orders. This paper highlights the need to investigate the natural variability of the atmosphere-ocean coupled dynamics by associating rate of growth and decay of perturbations to the physical modes described using the formalism of the covariant Lyapunov vectors and to consider long integrations in order to disentangle the dynamical processes occurring at all time scales.