Intensification of El Niño-Southern Oscillation (ENSO)-rainfall variability in response to global warming is a robust feature across Coupled Model Intercomparison Project (CMIP) iterations, regardless of a lack of robust projected changes in ENSO-sea-surface temperature (SST) variability. Previous studies attributed this intensification to an increase in mean SST and moisture convergence over the central-to-eastern Pacific, without explicitly considering underlying nonlinear SST–rainfall relationship changes. Here, by analyzing changes of the tropical SST–rainfall relationship of CMIP6 models, we present a mechanism linking the mean SST rise to amplifying ENSO–rainfall variability. We show that the slope of the SST–rainfall function over Niño3 region becomes steeper in a warmer climate, ~42.1% increase in 2050–2099 relative to 1950–1999, due to the increase in Clausius–Clapeyron-driven moisture sensitivity, ~16.2%, and dynamic contributions, ~25.9%. A theoretical reconstruction of ENSO–rainfall variability further supports this mechanism. Our results imply ENSO’s hydrological impacts increase nonlinearly in response to global warming.

To investigate the role of vegetation and ecosystem diversity on hominin adaptation and migration, we identify past human habitat preferences over time using a transient 3-million-year earth system-biome model simulation and an extensive hominin fossil and archaeological database. Our analysis shows that early African hominins predominantly lived in open environments such as grassland and dry shrubland. Migrating into Eurasia, hominins adapted to a broader range of biomes over time. By linking the location and age of hominin sites with corresponding simulated regional biomes, we also find that our ancestors actively selected for spatially diverse environments. The quantitative results lead to a new diversity hypothesis: Homo species, in particular Homo sapiens , were specially equipped to adapt to landscape mosaics.

The strength of the El Niño‐Southern Oscillation (ENSO)‐Indian summer monsoon rainfall (ISMR) relationship shows considerable decadal fluctuations, which have been previously linked to low‐frequency climatic processes such as shifts in ENSO's center of action or the Atlantic Multidecadal Oscillation. However, random variability can also cause similar variations in the relationship between climate phenomena. Here we propose a statistical test to determine whether the observed time‐evolving correlations between ENSO and ISMR are different from those expected from a simple stochastic null hypothesis model. The analysis focuses on the time evolution of moving correlations, their expected variance, and probabilities for rapid transitions. The results indicate that the time evolution of the observed running correlation between these climate modes is indistinguishable from a system in which the ISMR signal can be expressed as a stochastically perturbed ENSO signal. This challenges previous deterministic interpretations. Our results are further corroborated by the analysis of climate model simulations.

The oldest known hominin remains in Europe [~1.5 to ~1.1 million years ago (Ma)] have been recovered from Iberia, where paleoenvironmental reconstructions have indicated warm and wet interglacials and mild glacials, supporting the view that once established, hominin populations persisted continuously. We report analyses of marine and terrestrial proxies from a deep-sea core on the Portugese margin that show the presence of pronounced millennial-scale climate variability during a glacial period ~1.154 to ~1.123 Ma, culminating in a terminal stadial cooling comparable to the most extreme events of the last 400,000 years. Climate envelope–model simulations reveal a drastic decrease in early hominin habitat suitability around the Mediterranean during the terminal stadial. We suggest that these extreme conditions led to the depopulation of Europe, perhaps lasting for several successive glacial-interglacial cycles.

Driven primarily by variations in the earth's axis wobble, tilt, and orbit eccentricity, our planet experienced massive glacial/interglacial reorganizations of climate and atmospheric CO₂ concentrations during the Pleistocene (2.58 million years ago (Ma)–11.7 thousand years ago (ka)). Even after decades of research, the underlying climate response mechanisms to these astronomical forcings have not been fully understood. To further quantify the sensitivity of the earth system to orbital-scale forcings, we conducted an unprecedented quasi-continuous coupled general climate model simulation with the Community Earth System Model version 1.2 (CESM1.2, ∼3.75 ^∘ horizontal resolution), which covers the climatic history of the past 3 million years (3 Myr). In addition to the astronomical insolation changes, CESM1.2 is forced by estimates of CO₂ and ice-sheet topography which were obtained from a simulation previously conducted with the CLIMBER-2 earth system model of intermediate complexity. Our 3 Ma simulation consists of 42 transient interglacial/glacial simulation chunks, which were partly run in parallel to save computing time. The chunks were subsequently merged, accounting for spin-up and overlap effects to yield a quasi-continuous trajectory. The computer model data were compared against a plethora of paleo-proxy data and large-scale climate reconstructions. For the period from the Mid-Pleistocene Transition (MPT, ∼1 Ma) to the late Pleistocene we find good agreement between simulated and reconstructed temperatures in terms of phase and amplitude (−5.7 ^∘C temperature difference between Last Glacial Maximum and Holocene). For the earlier part (3–1 Ma), differences in orbital-scale variability occur between model simulation and the reconstructions, indicating potential biases in the applied CO₂ forcing. Our model-proxy data comparison also extends to the westerlies, which show unexpectedly large variance on precessional timescales, and hydroclimate variables in major monsoon regions. Eccentricity-modulated precessional variability is also responsible for the simulated changes in the amplitude and flavors of the El Niño–Southern Oscillation. We further identify two major modes of planetary energy transport, which played a crucial role in Pleistocene climate variability: the first obliquity and CO₂-driven mode is linked to changes in the Equator-to-pole temperature gradient; the second mode regulates the interhemispheric heat imbalance in unison with the eccentricity-modulated precession cycle. During the MPT, a pronounced qualitative shift occurs in the second mode of planetary energy transport: the post-MPT eccentricity-paced variability synchronizes with the CO₂ forced signal. This synchronized feature is coherent with changes in global atmospheric and ocean circulations, which might contribute to an intensification of glacial cycle feedbacks and amplitudes. Comparison of this paleo-simulation with greenhouse warming simulations reveals that for an RCP8.5 greenhouse gas emission scenario, the projected global mean surface temperature changes over the next 7 decades would be comparable to the late Pleistocene glacial-interglacial range; but the anthropogenic warming rate will exceed any previous ones by a factor of ∼100.