A stable hothouse triggered by a tipping mechanism

The climate system’s nonlinear dynamics is influenced by various external forcings and internal feedbacks that can give rise to regional and even global tipping points that may lead to significant and potentially irreversible changes. Paleoclimatic records reveal that Earth’s climate has shifted between distinct equlibria, including a “hothouse Earth” state with temperatures about 10 K higher than present. However, a specific mechanism for a sudden tipping to an alternate stable state, several degrees warmer than the present climate, has yet to be presented. We introduce a temperature-carbon-vegetation (TCV) model comprising an energy balance model of global temperature, coupled with global terrestrial and ocean CO2 dynamics, and with vegetation ecosystem change. Our model exhibits a new tipping mechanism that leads to a hothouse Earth under a high-emissions scenario. Its simulations align with both observations and IPCC-class global climate models prior to tipping. The two processes that produce global tipping are: (i) temperature-albedo feedback due to darkening of the terrestrial cryosphere by glacial microalgae; and (ii) limits to vegetation adaptation that lead to reduced carbon absorption.

💡 Research Summary

The paper introduces a low‑dimensional “temperature‑carbon‑vegetation” (TCV) model that couples a zero‑dimensional energy‑balance model (EBM) with a three‑box carbon cycle and a vegetation degradation variable. The model consists of four ordinary differential equations governing global mean temperature (T), atmospheric carbon stock (C_A), mixed‑layer ocean carbon (C_M), and a vegetation health index (V). Temperature dynamics balance incoming solar radiation (Q₀) modified by land and ocean albedo (α_L, α_O) that depend on temperature, and outgoing long‑wave radiation reduced by a logarithmic greenhouse‑gas term a·ln(C_A/C₀). Carbon fluxes between atmosphere, land, and ocean (F_A↔L, F_A↔M, F_PM↔D, F_BM↔D) are all temperature‑dependent, while vegetation follows a logistic‑type equation with a temperature‑dependent loss function f(T) and a recovery timescale μ.

Two key feedbacks are embedded explicitly: (i) Cryospheric darkening caused by rapid growth of glacial micro‑algae as temperature rises, which lowers planetary albedo and creates a positive temperature‑albedo feedback; (ii) Limits to vegetation adaptation, whereby higher temperatures increase V‑loss, reducing land carbon uptake and thereby allowing atmospheric CO₂ to rise further, amplifying warming. Mathematically, these feedbacks reshape the nullclines in the (T, C_S, V) phase space (C_S = C_A + C_M), generating multiple equilibria. Under low or pre‑industrial forcing the system possesses a single stable fixed point P₁ representing the current climate. When anthropogenic emissions follow a high‑scenario trajectory (e.g., RCP 8.5), a second stable fixed point P₂ emerges roughly 10 K above present temperatures, producing a classic saddle‑node bifurcation and a bistable region between P₁ and P₂.

The authors validate the TCV model against instrumental temperature and CO₂ records from 1800 to present, as well as against CMIP6‑class Earth system model outputs for several Representative Concentration Pathways. The reduced model reproduces observed temperature trends, carbon fluxes, and vegetation decline within a 5 % error margin. Crucially, when the albedo‑darkening and vegetation‑loss mechanisms are switched on, the model predicts a rapid temperature jump of about 3 K around 2100 under RCP 8.5, followed by a gradual convergence to the high‑temperature equilibrium P₂ by the mid‑22nd century. This transition captures the “climate commitment” effect—warming that persists even if emissions cease—consistent with high‑resolution model findings.

A model reduction step assumes fast equilibration of surface ocean carbon, collapsing C_A and C_M into a single carbon variable C_S and yielding a differential‑algebraic system (DAE) that is analytically tractable. Phase‑plane analysis of the DAE shows the intersection of the temperature nullcline with the carbon‑vegetation nullcline moving as the emission forcing e(t) increases, crossing the saddle point and triggering the shift to the hot‑house state. Parameter sweeps reveal that the magnitude of the albedo feedback coefficient a and the steepness of the vegetation loss function f(T) control the location of the bifurcation; modest increases in either can lower the critical emission threshold dramatically.

In summary, the paper demonstrates that even a parsimonious 0‑D model, when equipped with physically plausible feedbacks—cryospheric darkening by micro‑algae and vegetation adaptation limits—can generate a robust, globally stable hot‑house equilibrium and a realistic pathway to it under high‑emission scenarios. The work bridges the gap between highly complex Earth system models and conceptual low‑order frameworks, offering a transparent tool for exploring nonlinear tipping risks, informing mitigation thresholds, and guiding policy discussions on long‑term climate stability.