Evolution and forcing mechanisms of El Niño over the past 21,000 years

Future increase in extreme El Niño supported by past glacial changes

El Niño events, the warm phase of the El Niño-Southern Oscillation (ENSO) phenomenon, amplify climate variability throughout the world1. Uncertain climate model predictions limit our ability to assess whether these climatic events could become more extreme under anthropogenic greenhouse warming2. Palaeoclimate records provide estimates of past changes, but it is unclear if they can constrain mechanisms underlying future predictions3-5. Here we uncover a mechanism using numerical simulations that drives consistent changes in response to past and future forcings, allowing model validation against palaeoclimate data. The simulated mechanism is consistent with the dynamics of observed extreme El Niño events, which develop when western Pacific warm pool waters expand rapidly eastwards because of strongly coupled ocean currents and winds6,7. These coupled interactions weaken under glacial conditions because of a deeper mixed layer driven by a stronger Walker circulation. The resulting decrease in ENSO variability and extreme El Niño occurrence is supported by a series of tropical Pacific palaeoceanographic records showing reduced glacial temperature variability within key ENSO-sensitive oceanic regions, including new data from the central equatorial Pacific. The model-data agreement on past variability, together with the consistent mechanism across climatic states, supports the prediction of a shallower mixed layer and weaker Walker circulation driving more frequent extreme El Niño genesis under greenhouse warming.

Small populations of Palaeolithic humans in Cyprus hunted endemic megafauna to extinction

The hypothesized main drivers of megafauna extinctions in the late Quaternary have wavered between over-exploitation by humans and environmental change, with recent investigations demonstrating more nuanced synergies between these drivers depending on taxon, spatial scale, and region. However, most studies still rely on comparing archaeologically based chronologies of timing of initial human arrival into naïve ecosystems and palaeontologically inferred dates of megafauna extinctions. Conclusions arising from comparing chronologies also depend on the reliability of dated evidence, dating uncertainties, and correcting for the low probability of preservation (Signor-Lipps effect). While some models have been developed to test the susceptibility of megafauna to theoretical offtake rates, none has explicitly linked human energetic needs, prey choice, and hunting efficiency to examine the plausibility of human-driven extinctions. Using the island of Cyprus in the terminal Pleistocene as an ideal test case because of its late human settlement (~14.2-13.2 ka), small area (~11 000 km2), and low megafauna diversity (2 species), we developed stochastic models of megafauna population dynamics, with offtake dictated by human energetic requirements, prey choice, and hunting-efficiency functions to test whether the human population at the end of the Pleistocene could have caused the extinction of dwarf hippopotamus (Phanourios minor) and dwarf elephants (Palaeoloxodon cypriotes). Our models reveal not only that the estimated human population sizes (n = 3000-7000) in Late Pleistocene Cyprus could have easily driven both species to extinction within < 1000 years, the model predictions match the observed, Signor-Lipps-corrected chronological sequence of megafauna extinctions inferred from the palaeontological record (P. minor at ~12-11.1 ka, followed by P. cypriotes at ~10.3-9.1 ka).

A model of long-term population growth with an application to Central West Argentina

We propose an Ideal Specialization Model to help explain the diversity of population growth trajectories exhibited across archaeological regions over thousands of years. The model provides a general set of expectations useful for guiding empirical research, and we provide a concrete example by conducting a preliminary evaluation of three expectations in Central West Argentina. We use kernel density estimates of archaeological radiocarbon, estimates of paleoclimate, and human bone stable isotopes from archaeological remains to evaluate three expectations drawn from the model's dynamics. Based on our results, we suggest that innovations in the production of food and social organization drove demographic transitions and population expansion in the region. The consistency of population expansion in the region positively associates with changes in diet and, potentially, innovations in settlement and social integration.

Demographic models predict end-Pleistocene arrival and rapid expansion of pre-agropastoralist humans in Cyprus

The antiquity of human dispersal into Mediterranean islands and ensuing coastal adaptation have remained largely unexplored due to the prevailing assumption that the sea was a barrier to movement and that islands were hostile environments to early hunter-gatherers [J. F. Cherry, T. P. Leppard, J. Isl. Coast. Archaeol. 13, 191-205 (2018), 10.1080/15564894.2016.1276489]. Using the latest archaeological data, hindcasted climate projections, and age-structured demographic models, we demonstrate evidence for early arrival (14,257 to 13,182 calendar years ago) to Cyprus and predicted that large groups of people (~1,000 to 1,375) arrived in 2 to 3 main events occurring within <100 y to ensure low extinction risk. These results indicate that the postglacial settlement of Cyprus involved only a few large-scale, organized events requiring advanced watercraft technology. Our spatially debiased and Signor-Lipps-corrected estimates indicate rapid settlement of the island within <200 y, and expansion to a median of 4,000 to 5,000 people (0.36 to 0.46 km-2) in <11 human generations (<300 y). Our results do not support the hypothesis of inaccessible and inhospitable islands in the Mediterranean for pre-agropastoralists, agreeing with analogous conclusions for other parts of the world [M. I. Bird et al., Sci. Rep. 9, 8220 (2019), 10.1038/s41598-019-42946-9]. Our results also highlight the need to revisit these questions in the Mediterranean and test their validity with new technologies, field methods, and data. By applying stochastic models to the Mediterranean region, we can place Cyprus and large islands in general as attractive and favorable destinations for paleolithic peoples.

Reconstruction of ENSO variability using the standardized growth index of a Tridacna shell from Yongshu Reef, South China Sea

El Niño-Southern Oscillation (ENSO) is the strongest signal of global interannual climate anomaly and reconstructing past ENSO variations using high-resolution paleoclimate archives can improve our understanding of ENSO variability, as well as improve our ability to predict future climate changes. Here, a daily resolution standardized growth index (SGI) was established using a giant clam (Tridacna spp.) shell specimen MD2 (life span: 1994-2013 CE), collected from the Yongshu Reef, southern South China Sea (SCS). The cross-spectral and correlation analysis indicated that the SGI variation of MD2 was strongly influenced by ENSO variability on an interannual timescale. Tridacna spp. is in symbiosis with zooxanthellae, and its growth index is usually modulated by the photosynthetic efficiency of zooxanthellae. During the El Niño (La Niña) period, the convective anomalies stimulated in western Pacific would increase (decrease) the effective solar radiation on Yongshu Reef, and in turn influence the photosynthesis rate of zooxanthellae and enzyme activity for the calcification site and thus the SGI of giant clam MD2. The SGI can explain 54.7 % of ENSO variance, demonstrating the potential for Tridacna SGI in ENSO reconstruction. Compared with conventional ENSO reconstruction using high-resolution geochemical proxies, the method of giant clam SGI is rapid and economical.

A contracting Intertropical Convergence Zone during the Early Heinrich Stadial 1

Despite the fact that the response of tropical hydroclimate to North Atlantic cooling events during the Heinrich Stadial 1 (HS1) has been extensively studied in African, South American and Indonesia, the nature of such responses remains debated. Here we investigate the tropical hydroclimate pattern over the Indo-Asian-Australian monsoon region during the HS1 by integrating hydroclimatic records, and examining a δ18Oseawater record from Globigerinoides ruber (white) in the tropical Indian Ocean. Our findings indicate that tropical hydrological conditions were synchronously arid in both hemispheres during the early HS1 (~18.3-16.3 ka) in the Indo-Asian-Australian monsoon region, except for a narrow, wet hydrological belt in northern low latitudes, suggesting the existence of a contracted tropical precipitation belt at that time. This study reveals that the meltwater discharge and resulting changes in global temperatures and El Niño exerted a profound influence on the tropical hydroclimate in the Indo-Asian-Australian monsoon region during the early HS1.

Deglacial increase of seasonal temperature variability in the tropical ocean

The relatively stable Holocene climate was preceded by a pronounced event of abrupt warming in the Northern Hemisphere, the termination of the Younger Dryas (YD) cold period1,2. Although this transition has been intensively studied, its imprint on low-latitude ocean temperature is still controversial and its effects on sub-annual to decadal climate variability remain poorly understood1,3,4. Sea surface temperature (SST) variability at these timescales in the tropical Atlantic is expected to intensify under current and future global warming and has considerable consequences for environmental conditions in Africa and South America, and for tropical Pacific climate5–8. Here we present a 100-µm-resolution record obtained by mass spectrometry imaging (MSI) of long-chain alkenones in sediments from the Cariaco Basin9–11 and find that annually averaged SST remained stable during the transition into the Holocene. However, seasonality increased more than twofold and approached modern values of 1.6 °C, probably driven by the position and/or annual range of the Intertropical Convergence Zone (ITCZ). We further observe that interannual variability intensified during the early Holocene. Our results demonstrate that sub-decadal-scale SST variability in the tropical Atlantic is sensitive to abrupt changes in climate background, such as those witnessed during the most recent glacial to interglacial transition.

Thermal coupling of the Indo-Pacific warm pool and Southern Ocean over the past 30,000 years

The role of the tropical Pacific Ocean and its linkages to the southern hemisphere during the last deglacial warming remain highly controversial. Here we explore the evolution of Pacific horizontal and vertical thermal gradients over the past 30 kyr by compiling 340 sea surface and 7 subsurface temperature records, as well as one new ocean heat content record. Our records reveal that La Niña-like conditions dominated during the deglaciation as a result of the more intense warming in the western Pacific warm pool. Both the subsurface temperature and ocean heat content in the warm pool rose earlier than the sea surface temperature, and in phase with South Pacific subsurface temperature and orbital precession, implying that heat exchange between the tropical upper water column and the extratropical Southern Ocean facilitated faster warming in the western Pacific. Our study underscores the key role of the thermal coupling between the warm pool and the Southern Ocean and its relevance for future global warming.

The mechanism of the last deglacial global warming is key for future climate. Here, the authors shed light on the pivotal role of the thermal coupling between the western Pacific warm pool and the Southern Ocean.

El Niño frequency threshold controls coastal biotic communities

El Niño has profound influences on ecosystem dynamics. However, we know little about how it shapes vertebrate faunal community composition over centennial time scales, and this limits our ability to forecast change under projections of future El Niño events. On the basis of correlations between geological records of past El Niño frequency and the species composition of bird and fish remains from a Baja California bone deposit that spans the past 12,000 years, we documented marked faunal restructuring when major El Niño events occurred more than five times per century. This tipping point has implications for the past and future ecology of eastern Pacific coastal environments.

Unraveling forced responses of extreme El Niño variability over the Holocene

Uncertainty surrounding the future response of El Niño-Southern Oscillation (ENSO) variability to anthropogenic warming necessitates the study of past ENSO sensitivity to substantial climate forcings over geological history. Here, we focus on the Holocene epoch and show that ENSO amplitude and frequency intensified over this period, driven by an increase in extreme El Niño events. Our study combines new climate model simulations, advances in coral proxy system modeling, and coral proxy data from the central tropical Pacific. Although the model diverges from the observed coral data regarding the exact magnitude of change, both indicate that modern ENSO variance eclipsed paleo-estimates over the Holocene, albeit against the backdrop of wide-ranging natural variability. Toward further constraining paleo-ENSO, our work underscores the need for multimodel investigations of additional Holocene intervals alongside more coral data from periods with larger climate forcing. Our findings implicate extreme El Niño events as an important rectifier of mean ENSO intensity.

Climate and demography drive 7000 years of dietary change in the Central Andes

Explaining the factors that influence past dietary variation is critically important for understanding changes in subsistence, health, and status in past societies; yet systematic studies comparing possible driving factors remain scarce. Here we compile the largest dataset of past diet derived from stable isotope δ¹³C‰ and δ¹⁵N‰ values in the Americas to quantitatively evaluate the impact of 7000 years of climatic and demographic change on dietary variation in the Central Andes. Specifically, we couple paleoclimatic data from a general circulation model with estimates of relative past population inferred from archaeologically derived radiocarbon dates to assess the influence of climate and population on spatiotemporal dietary variation using an ensemble machine learning model capable of accounting for interactions among predictors. Results reveal that climate and population strongly predict diet (80% of δ15N‰ and 66% of δ13C‰) and that Central Andean diets correlate much more strongly with local climatic conditions than regional population size, indicating that the past 7000 years of dietary change was influenced more by climatic than socio-demographic processes. Visually, the temporal pattern suggests decreasing dietary variation across elevation zones during the Late Horizon, raising the possibility that sociopolitical factors overrode the influence of local climatic conditions on diet during that time. The overall findings and approach establish a general framework for understanding the influence of local climate and demography on dietary change across human history.

Enhanced North Pacific subtropical gyre circulation during the late Holocene

The North Pacific Subtropical Gyre circulation redistributes heat from the Western Pacific Warm Pool towards the mid- to high-latitude North Pacific. However, the driving mechanisms of this circulation and how it changed over the Holocene remain poorly understood. Here, we present alkenone-based sea surface temperature reconstructions along the Kuroshio, California and Alaska currents that cover the past ~7,000 years. These and other paleorecords collectively demonstrate a coherent intensification of the boundary currents, and thereby the basin-scale subtropical gyre circulation, since ~3,000-4,000 years ago. Such enhanced circulation during the late Holocene appears to have resulted from a long-term southward migration of the Intertropical Convergence Zone, associated with Holocene ocean cooling. Our results imply that the North Pacific Subtropical Gyre circulation could be weakened under future global warming.

High- and low-latitude forcings drive Atacama Desert rainfall variations over the past 16,000 years

Late Quaternary precipitation dynamics in the central Andes have been linked to both high- and low-latitude atmospheric teleconnections. We use present-day relationships between fecal pellet diameters from ashy chinchilla rats (Abrocoma cinerea) and mean annual rainfall to reconstruct the timing and magnitude of pluvials (wet episodes) spanning the past 16,000 years in the Atacama Desert based on 81 ¹⁴C-dated A. cinerea paleomiddens. A transient climate simulation shows that pluvials identified at 15.9 to 14.8, 13.0 to 8.6, and 8.1 to 7.6 ka B.P. can be linked to North Atlantic (high-latitude) forcing (e.g., Heinrich Stadial 1, Younger Dryas, and Bond cold events). Holocene pluvials at 5.0 to 4.6, 3.2 to 2.1, and 1.4 to 0.7 ka B.P. are not simulated, implying low-latitude internal variability forcing (i.e., ENSO regime shifts). These results help constrain future central Andean hydroclimatic variability and hold promise for reconstructing past climates from rodent middens in desert ecosystems worldwide.

Hydrological variation recorded in a subalpine peatland of Northeast Asia since the Little Ice Age and its possible driving mechanisms

Hydrological characteristics since the Little Ice Age (LIA) could provide a good reference for current climate analysis and future climate prediction. However, the hydrological variation since the LIA and its driving mechanisms in Northeast Asia remain unclear, which has severely restricted our understanding on the past, present and future hydroclimate changes in these regions. Here we reconstruct the hydrological dynamics over the past 700 years using samples from the Hani peatland a subalpine peatland of Changbai Mountains to reveal these issues. The analytical results from plant macrofossil and grain-size of the HN-1 core and the integrated moisture/precipitation records across the entire Northeast Asia indicate that the hydrological environments in Northeast Asia were wetter conditions during the period of 1300-1700 AD, dry conditions during the period of 1700-1850 AD, and wet conditions during the period of 1850-2018 AD, respectively. The possible driving mechanisms for the hydrological variations in Northeast Asia since the LIA can be divided into three models. La Niña-like conditions induced wetter conditions in Northeast Asia from 1300 to 1700 AD. From 1700 to 1850 AD, strong volcanic aerosol effects superimposed on weaker La Niña-like conditions, resulting in dry conditions in Northeast Asia. However, El Niño-like conditions induced wet conditions in Northeast Asia from 1850 to 2018 AD. These driving models suggest that the teleconnected influence of solar activity/sunspot could control the hydrological dynamics in Northeast Asia on a decadal-centennial scale through the ENSO activities and Walker Circulation changes since the LIA. Based on the periodicity relationship between hydrological conditions and sunspot, it can be predicted that the moisutre conditions in Northeast Asia would gradually decrease from 2030 to 2085 AD, and gradually increase from 2085 to 2140 AD.

Vegetation feedback causes delayed ecosystem response to East Asian Summer Monsoon Rainfall during the Holocene

One long-standing issue in the paleoclimate records is whether East Asian Summer Monsoon peaked in the early Holocene or mid-Holocene. Here, combining a set of transient earth system model simulations with proxy records, we propose that, over northern China, monsoon rainfall peaked in the early Holocene, while soil moisture and tree cover peaked in the mid-Holocene. The delayed ecosystem (soil moisture and tree cover) response to rainfall is caused by the vegetation response to winter warming and the subsequent feedback with soil moisture. Our study provides a mechanism for reconciling different evolution behaviors of monsoon proxy records; it sheds light on the driving mechanism of the monsoon evolution and monsoon-ecosystem feedback over northern China, with implications to climate changes in other high climate sensitivity regions over the globe.

Modulation of late Pleistocene ENSO strength by the tropical Pacific thermocline

The El Niño Southern Oscillation (ENSO) is highly dependent on coupled atmosphere-ocean interactions and feedbacks, suggesting a tight relationship between ENSO strength and background climate conditions. However, the extent to which background climate state determines ENSO behavior remains in question. Here we present reconstructions of total variability and El Niño amplitude from individual foraminifera distributions at discrete time intervals over the past ~285,000 years across varying atmospheric CO₂ levels, global ice volume and sea level, and orbital insolation forcing. Our results show a strong correlation between eastern tropical Pacific Ocean mixed-layer thickness and both El Niño amplitude and central Pacific variability. This ENSO-thermocline relationship implicates upwelling feedbacks as the major factor controlling ENSO strength on millennial time scales. The primacy of the upwelling feedback in shaping ENSO behavior across many different background states suggests accurate quantification and modeling of this feedback is essential for predicting ENSO’s behavior under future climate conditions.

How the El Niño Southern Oscillation depends on the background conditions is not well known. Here, the authors present individual foraminifera distributions which show that central Pacific variability is related to the warmth and depth of the thermocline across varying climate background conditions over the past ~285,000 years.

Different precipitation response over land and ocean to orbital and greenhouse gas forcing

Various proxies suggest a nearly in-phase variation of monsoons with local summer insolation. Oceanic proxies of monsoons document a more complex response. Climate model simulations also indicate that the response is different over land and ocean. Here using a transient simulation by a climate model over the last 22,000 years we have unraveled the factors that lead to these differences within the Indian subcontinent. We show that during the deglacial (22–12 ka) precipitation over India and the Bay of Bengal (BoB) are in phase, whereas they are out of phase across the Holocene (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim$$\end{document}∼ 12 ka to 0 ka). During the deglacial, water vapor amplifies the effect of solar forcing on precipitation over both the regions, whereas contributions from surface latent heat fluxes over the BoB drive an opposite response across the Holocene. We find that greenhouse gas forcing drives similar precipitation response over land and ocean, whereas orbital forcing produces a different response over land and ocean. We have further demonstrated that during periods of abrupt climate change [such as the Bølling–Allerød (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim$$\end{document}∼ 14 ka)], water vapor affects precipitation mainly through its influence on the vertical stability of the atmosphere. These results highlight the complex nature of precipitation over the BoB and thus has implications for the interpretation of monsoon proxies.

The Role of El Niño in Driving Drought Conditions over the Last 2000 Years in Thailand

Irregular climate events frequently occur in Southeast Asia due to the numerous climate patterns combining. Thailand sits at the confluence of these interactions, and consequently experiences major hydrological events, such as droughts. Proxy data, speleothem records, lake sediment sequences and tree ring chronologies were used to reconstruct paleo drought conditions. These trends were compared with modelled and historic El Niño Southern Oscillation (ENSO) data to assess if the ENSO climate phenomena is causing droughts in Thailand. Drought periods were found to occur both during El Niño events and ENSO neutral conditions. This indicates droughts are not a product of one climate pattern, but likely the result of numerous patterns interacting. There is uncertainty regarding how climate patterns will evolve under climate change, but changes in amplitude and variability could potentially lead to more frequent and wider reaching hydrological disasters. It is vital that policies are implemented to cope with the resulting social and economic repercussions, including diversification of crops and reorganisation of water consumption behaviour in Thailand.

The spatiotemporal spread of human migrations during the European Holocene

The European continent was subject to two major migrations of peoples during the Holocene: the northwestward movement of Anatolian farmer populations during the Neolithic and the westward movement of Yamnaya steppe peoples during the Bronze Age. These movements changed the genetic composition of the continent's inhabitants. The Holocene was also characterized by major changes in vegetation composition, which altered the environment occupied by the original hunter-gatherer populations. We aim to test to what extent vegetation change through time is associated with changes in population composition as a consequence of these migrations, or with changes in climate. Using ancient DNA in combination with geostatistical techniques, we produce detailed maps of ancient population movements, which allow us to visualize how these migrations unfolded through time and space. We find that the spread of Neolithic farmer ancestry had a two-pronged wavefront, in agreement with similar findings on the cultural spread of farming from radiocarbon-dated archaeological sites. This movement, however, did not have a strong association with changes in the vegetational landscape. In contrast, the Yamnaya migration speed was at least twice as fast and coincided with a reduction in the amount of broad-leaf forest and an increase in the amount of pasture and natural grasslands in the continent. We demonstrate the utility of integrating ancient genomes with archaeometric datasets in a spatiotemporal statistical framework, which we foresee will enable future studies of ancient populations' movements, and their putative effects on local fauna and flora.

Half-precessional cycle of thermocline temperature in the western equatorial Pacific and its bihemispheric dynamics

The El Niño-Southern Oscillation (ENSO), which is tightly coupled to the equatorial thermocline in the Pacific, is the dominant source of interannual climate variability, but its long-term evolution in response to climate change remains highly uncertain. This study uses Mg/Ca in planktonic foraminiferal shells to reconstruct sea surface and thermocline water temperatures (SST and TWT) for the past 142 ky in a western equatorial Pacific (WEP) core MD01-2386. Unlike the dominant 100-ky glacial-interglacial cycle recorded by SST and δ¹⁸O, which echoes the pattern seen in other WEP sites, the upper ocean thermal gradient shows a clear half-precessional (9.4 ky or 12.7 ky) cycle as indicated by the reconstructed and simulated temperature (ΔT) and δ¹⁸O differences between the surface and thermocline waters. This phenomenon is attributed to the interplay of subtropical-to-tropical thermocline anomalies forced by the antiphased meridional insolation gradients in the two hemispheres at the precessional band. In particular, the TWT shows greater variability than SST, and dominates the ΔT changes which couple with the west-east SST difference in the equatorial Pacific at the half-precessional band, implying a decisive role of the tropical thermocline in orbital-scale climate change.

Modulation of Indian monsoon by water vapor and cloud feedback over the past 22,000 years

To predict how monsoons will evolve in the 21st century, we need to understand how they have changed in the past. In paleoclimate literature, the major focus has been on the role of solar forcing on monsoons but not on the amplification by feedbacks internal to the climate system. Here we have used the results from a transient climate simulation to show that feedbacks amplify the effect of change in insolation on the Indian summer monsoon. We show that during the deglacial (22 ka to 10 ka) monsoons were predominantly influenced by rising water vapor due to increasing sea surface temperature, whereas in the Holocene (10 ka to 0 ka) cloud feedback was more important. These results are consistent with another transient simulation, thus increasing confidence despite potential model biases. We have demonstrated that insolation drives monsoon through different pathways during cold and warm periods, thereby highlighting the changing role of internal factors.

Synchronous 500-year oscillations of monsoon climate and human activity in Northeast Asia

Prehistoric human activities were likely influenced by cyclic monsoon climate changes in East Asia. Here we report a decadal-resolution Holocene pollen record from an annually-laminated Maar Lake in Northeast China, a proxy of monsoon climate, together with a compilation of 627 radiocarbon dates from archeological sites in Northeast China which is a proxy of human activity. The results reveal synchronous ~500-year quasi-periodic changes over the last 8000 years. The warm-humid/cold-dry phases of monsoon cycles correspond closely to the intensification/weakening of human activity and the flourishing/decline of prehistoric cultures. Six prosperous phases of prehistoric cultures, with one exception, correspond approximately to warm-humid phases caused by a strengthened monsoon. This ~500-year cyclicity in the monsoon and thus environmental change triggered the development of prehistoric cultures in Northeast China. The cyclicity is apparently linked to the El Niño-Southern Oscillation, against the background of long-term Holocene climatic evolution. These findings reveal a pronounced relationship between prehistoric human activity and cyclical climate change.

Long-term climate cycles can potentially influence population dynamics, including those of humans. Here, the authors combine climate and archaeological records from Northeast China over the past 8000 years and demonstrate ~500 year cycles in both the monsoon and human activity.

The Asian Summer Monsoon: Teleconnections and Forcing Mechanisms—A Review from Chinese Speleothem δ18O Records

Asian summer monsoon (ASM) variability significantly affects hydro-climate, and thus socio-economics, in the East Asian region, where nearly one-third of the global population resides. Over the last two decades, speleothem δ18O records from China have been utilized to reconstruct ASM variability and its underlying forcing mechanisms on orbital to seasonal timescales. Here, we use the Speleothem Isotopes Synthesis and Analysis database (SISAL_v1) to present an overview of hydro-climate variability related to the ASM during three periods: the late Pleistocene, the Holocene, and the last two millennia. We highlight the possible global teleconnections and forcing mechanisms of the ASM on different timescales. The longest composite stalagmite δ18O record over the past 640 kyr BP from the region demonstrates that ASM variability on orbital timescales is dominated by the 23 kyr precessional cycles, which are in phase with Northern Hemisphere summer insolation (NHSI). During the last glacial, millennial changes in the intensity of the ASM appear to be controlled by North Atlantic climate and oceanic feedbacks. During the Holocene, changes in ASM intensity were primarily controlled by NHSI. However, the spatio-temporal distribution of monsoon rain belts may vary with changes in ASM intensity on decadal to millennial timescales.

Runoff events and related rainfall variability in the Southern Carpathians during the last 2000 years

The occurrence of heavy rainfall events is expected to undergo significant changes under increasing anthropogenic forcing. South-eastern Europe is reacting rapidly to such changes, therefore understanding and forecasting of precipitation variability is vital to better comprehending environmental changes in this area. Here we present a sub-decadal reconstruction of enhanced rainfall events for the past 2000 years from the Southern Carpathians, Romania using peat geochemistry. Five clear periods of enhanced rainfall are identified at 125–250, 600–900, 1050–1300, 1400–1575 and 1725–1980 CE. Significant runoff is observed during the second half of the Medieval Warm Period, whilst the Little Ice Age was characterised by significant variability. The North Atlantic Oscillation appears to be the main control on regional precipitation, but changes in solar irradiance also seem to play a significant role, together with the Siberian High. Comparison of the data presented here with model outputs confirms the ability of models to predict general trends, and major shifts, but highlights the complexity of the region’s hydrological history.

Dynamic Regimes of El Niño Southern Oscillation and Influenza Pandemic Timing

El Niño southern oscillation (ENSO) dynamics has been shown to drive seasonal influenza dynamics. Severe seasonal influenza epidemics and the 2009-2010 pandemic were coincident with chaotic regime of ENSO dynamics. ENSO dynamics from 1876 to 2016 were characterized to determine if influenza pandemics are coupled to chaotic regimes. Time-varying spectra of southern oscillation index (SOI) and sea surface temperature (SST) were compared. SOI and SST were decomposed to components using the algorithm of noise-assisted multivariate empirical mode decomposition. The components were Hilbert transformed to generate instantaneous amplitudes and phases. The trajectories and attractors of components were characterized in polar coordinates and state space. Influenza pandemics were mapped to dynamic regimes of SOI and SST joint recurrence of annual components. State space geometry of El Niños lagged by influenza pandemics were characterized and compared with other El Niños. Timescales of SOI and SST components ranged from sub-annual to multidecadal. The trajectories of SOI and SST components and the joint recurrence of annual components were dissipative toward chaotic attractors. Periodic, quasi-periodic, and chaotic regimes were present in the recurrence of trajectories, but chaos-chaos transitions dominated. Influenza pandemics occurred during chaotic regimes of significantly low transitivity dimension (p < 0.0001). El Niños lagged by influenza pandemics had distinct state space geometry (p < 0.0001). Chaotic dynamics explains the aperiodic timing, and varying duration and strength of El Niños. Coupling of all influenza pandemics of the past 140 years to chaotic regimes of low transitivity indicate that ENSO dynamics drives influenza pandemic dynamics. Forecasts models from ENSO dynamics should compliment surveillance for novel influenza viruses.

Greening of the Sahara suppressed ENSO activity during the mid-Holocene

The evolution of the El Niño-Southern Oscillation (ENSO) during the Holocene remains uncertain. In particular, a host of new paleoclimate records suggest that ENSO internal variability or other external forcings may have dwarfed the fairly modest ENSO response to precessional insolation changes simulated in climate models. Here, using fully coupled ocean-atmosphere model simulations, we show that accounting for a vegetated and less dusty Sahara during the mid-Holocene relative to preindustrial climate can reduce ENSO variability by 25%, more than twice the decrease obtained using orbital forcing alone. We identify changes in tropical Atlantic mean state and variability caused by the momentous strengthening of the West Africa Monsoon (WAM) as critical factors in amplifying ENSO’s response to insolation forcing through changes in the Walker circulation. Our results thus suggest that potential changes in the WAM due to anthropogenic warming may influence ENSO variability in the future as well.

Evolution of the El Niño-Southern Oscillation through the Holocene remains uncertain. Here, via fully coupled model simulations, the authors show that increased Saharan vegetation and reduced dust emissions 6 kyr BP significantly affect ENSO variability through changes in the West African Monsoon strength.

Holocene ENSO-related cyclic storms recorded by magnetic minerals in speleothems of central China

Extreme hydrologic events such as storms and floods have the potential to severely impact modern human society. However, the frequency of storms and their underlying mechanisms are limited by a paucity of suitable proxies, especially in inland areas. Here we present a record of speleothem magnetic minerals to reconstruct paleoprecipitation, including storms, in the eastern Asian monsoon area over the last 8.6 ky. The geophysical parameter IRM_soft-flux represents the flux of soil-derived magnetic minerals preserved in stalagmite HS4, which we correlate with rainfall amount and intensity. IRM_soft-flux exhibits relatively higher values before 6.7 ky and after 3.4 ky and lower values in the intervening period, consistent with regional hydrological changes observed in independent records. Abrupt enhancements in the flux of pedogenic magnetite in the stalagmite agree well with the timing of known regional paleofloods and with equatorial El Niño-Southern Oscillation (ENSO) patterns, documenting the occurrence of ENSO-related storms in the Holocene. Spectral power analyses reveal that the storms occur on a significant 500-y cycle, coincident with periodic solar activity and ENSO variance, showing that reinforced (subdued) storms in central China correspond to reduced (increased) solar activity and amplified (damped) ENSO. Thus, the magnetic minerals in speleothem HS4 preserve a record of the cyclic storms controlled by the coupled atmosphere-oceanic circulation driven by solar activity.

Equatorial Pacific forcing of western Amazonian precipitation during Heinrich Stadial 1

Abundant hydroclimatic evidence from western Amazonia and the adjacent Andes documents wet conditions during Heinrich Stadial 1 (HS1, 18–15 ka), a cold period in the high latitudes of the North Atlantic. This precipitation anomaly was attributed to a strengthening of the South American summer monsoon due to a change in the Atlantic interhemispheric sea surface temperature (SST) gradient. However, the physical viability of this mechanism has never been rigorously tested. We address this issue by combining a thorough compilation of tropical South American paleorecords and a set of atmosphere model sensitivity experiments. Our results show that the Atlantic SST variations alone, although leading to dry conditions in northern South America and wet conditions in northeastern Brazil, cannot produce increased precipitation over western Amazonia and the adjacent Andes during HS1. Instead, an eastern equatorial Pacific SST increase (i.e., 0.5–1.5 °C), in response to the slowdown of the Atlantic Meridional Overturning Circulation during HS1, is crucial to generate the wet conditions in these regions. The mechanism works via anomalous low sea level pressure over the eastern equatorial Pacific, which promotes a regional easterly low-level wind anomaly and moisture recycling from central Amazonia towards the Andes.

Using seafaring simulations and shortest-hop trajectories to model the prehistoric colonization of Remote Oceania

The prehistoric colonization of islands in Remote Oceania that began ∼3400 B.P. represents what was arguably the most expansive and ambitious maritime dispersal of humans across any of the world's seas or oceans. Though archaeological evidence has provided a relatively clear picture of when many of the major island groups were colonized, there is still considerable debate as to where these settlers originated from and their strategies/trajectories used to reach habitable land that other datasets (genetic, linguistic) are also still trying to resolve. To address these issues, we have harnessed the power of high-resolution climatic and oceanographic datasets in multiple seafaring simulation platforms to examine major pulses of colonization in the region. Our analysis, which takes into consideration currents, land distribution, wind periodicity, the influence of El Niño Southern Oscillation (ENSO) events, and "shortest-hop" trajectories, demonstrate that (i) seasonal and semiannual climatic changes were highly influential in structuring ancient Pacific voyaging; (ii) western Micronesia was likely settled from somewhere around the Maluku (Molucca) Islands; (iii) Samoa was the most probable staging area for the colonization of East Polynesia; and (iv) although there are major differences in success rates depending on time of year and the occurrence of ENSO events, settlement of Hawai'i and New Zealand is possible from the Marquesas or Society Islands, the same being the case for settlement of Easter Island from Mangareva or the Marquesas.

Correlation and anti-correlation of the East Asian summer and winter monsoons during the last 21,000 years

Understanding the past significant changes of the East Asia Summer Monsoon (EASM) and Winter Monsoon (EAWM) is critical for improving the projections of future climate over East Asia. One key issue that has remained outstanding from the paleo-climatic records is whether the evolution of the EASM and EAWM are correlated. Here, using a set of long-term transient simulations of the climate evolution of the last 21,000 years, we show that the EASM and EAWM are positively correlated on the orbital timescale in response to the precessional forcing, but are anti-correlated on millennial timescales in response to North Atlantic melt water forcing. The relation between EASM and EAWM can differ dramatically for different timescales because of the different response mechanisms, highlighting the complex dynamics of the East Asian monsoon system and the challenges for future projection.

Future projection of changes in the East Asia Summer and Winter Monsoon are hindered by a lack of understanding of past variability. Here, using longterm transient simulations, the authors show that the monsoons respond in phase to precessional forcing, yet out of phase millennial-scale North Atlantic forcing.

Recent Shift in Climate Relationship Enables Prediction of the Timing of Bird Breeding

Large-scale climate processes influence many aspects of ecology including breeding phenology, reproductive success and survival across a wide range of taxa. Some effects are direct, for example, in temperate-zone birds, ambient temperature is an important cue enabling breeding effort to coincide with maximum food availability, and earlier breeding in response to warmer springs has been documented in many species. In other cases, time-lags of up to several years in ecological responses have been reported, with effects mediated through biotic mechanisms such as growth rates or abundance of food supplies. Here we use 23 years of data for a temperate woodland bird species, the great tit (Parus major), breeding in deciduous woodland in eastern England to demonstrate a time-lagged linear relationship between the on-set of egg laying and the winter index of the North Atlantic Oscillation such that timing can be predicted from the winter index for the previous year. Thus the timing of bird breeding (and, by inference, the timing of spring events in general) can be predicted one year in advance. We also show that the relationship with the winter index appears to arise through an abiotic time-lag with local spring warmth in our study area. Examining this link between local conditions and larger-scale processes in the longer-term showed that, in the past, significant relationships with the immediately preceding winter index were more common than those with the time-lagged index, and especially so from the late 1930s to the early 1970s. However, from the mid 1970s onwards, the time-lagged relationship has become the most significant, suggesting a recent change in climate patterns. The strength of the current time-lagged relationship suggests that it might have relevance for other temperature-dependent ecological relationships.

Heat Transport Compensation in Atmosphere and Ocean over the Past 22,000 Years

The Earth’s climate has experienced dramatic changes over the past 22,000 years; however, the total meridional heat transport (MHT) of the climate system remains stable. A 22,000-year-long simulation using an ocean-atmosphere coupled model shows that the changes in atmosphere and ocean MHT are significant but tend to be out of phase in most regions, mitigating the total MHT change, which helps to maintain the stability of the Earth’s overall climate. A simple conceptual model is used to understand the compensation mechanism. The simple model can reproduce qualitatively the evolution and compensation features of the MHT over the past 22,000 years. We find that the global energy conservation requires the compensation changes in the atmosphere and ocean heat transports. The degree of compensation is mainly determined by the local climate feedback between surface temperature and net radiation flux at the top of the atmosphere. This study suggests that an internal mechanism may exist in the climate system, which might have played a role in constraining the global climate change over the past 22,000 years.

Regional and global forcing of glacier retreat during the last deglaciation

The ongoing retreat of glaciers globally is one of the clearest manifestations of recent global warming associated with rising greenhouse gas concentrations. By comparison, the importance of greenhouse gases in driving glacier retreat during the most recent deglaciation, the last major interval of global warming, is unclear due to uncertainties in the timing of retreat around the world. Here we use recently improved cosmogenic-nuclide production-rate calibrations to recalculate the ages of 1,116 glacial boulders from 195 moraines that provide broad coverage of retreat in mid-to-low-latitude regions. This revised history, in conjunction with transient climate model simulations, suggests that while several regional-scale forcings, including insolation, ice sheets and ocean circulation, modulated glacier responses regionally, they are unable to account for global-scale retreat, which is most likely related to increasing greenhouse gas concentrations.

The extent to which greenhouse gases forced glacier retreat during the last deglaciation remains unclear. Here, the authors recalculate cosmogenic nuclide ages for 195 glacier moraines and show that deglacial glacier retreat was broadly globally synchronous with rising levels of atmospheric CO₂.

Regional and global forcing of glacier retreat during the last deglaciation

The ongoing retreat of glaciers globally is one of the clearest manifestations of recent global warming associated with rising greenhouse gas concentrations. By comparison, the importance of greenhouse gases in driving glacier retreat during the most recent deglaciation, the last major interval of global warming, is unclear due to uncertainties in the timing of retreat around the world. Here we use recently improved cosmogenic-nuclide production-rate calibrations to recalculate the ages of 1,116 glacial boulders from 195 moraines that provide broad coverage of retreat in mid-to-low-latitude regions. This revised history, in conjunction with transient climate model simulations, suggests that while several regional-scale forcings, including insolation, ice sheets and ocean circulation, modulated glacier responses regionally, they are unable to account for global-scale retreat, which is most likely related to increasing greenhouse gas concentrations.