Adaptive foraging behaviours in the Horn of Africa during Toba supereruption.

Although modern humans left Africa multiple times over 100,000 years ago, those broadly ancestral to non-Africans dispersed less than 100,000 years ago1. Most models hold that these events occurred through green corridors created during humid periods because arid intervals constrained population movements2. Here we report an archaeological site-Shinfa-Metema 1, in the lowlands of northwest Ethiopia, with Youngest Toba Tuff cryptotephra dated to around 74,000 years ago-that provides early and rare evidence of intensive riverine-based foraging aided by the likely adoption of the bow and arrow. The diet included a wide range of terrestrial and aquatic animals. Stable oxygen isotopes from fossil mammal teeth and ostrich eggshell show that the site was occupied during a period of high seasonal aridity. The unusual abundance of fish suggests that capture occurred in the ever smaller and shallower waterholes of a seasonal river during a long dry season, revealing flexible adaptations to challenging climatic conditions during the Middle Stone Age. Adaptive foraging along dry-season waterholes would have transformed seasonal rivers into 'blue highway' corridors, potentially facilitating an out-of-Africa dispersal and suggesting that the event was not restricted to times of humid climates. The behavioural flexibility required to survive seasonally arid conditions in general, and the apparent short-term effects of the Toba supereruption in particular were probably key to the most recent dispersal and subsequent worldwide expansion of modern humans.

A unified explanation for the morphology of raised peatlands.

Raised peatlands, or bogs, are gently mounded landforms that are composed entirely of organic matter1-4 and store the most carbon per area of any terrestrial ecosystem5. The shapes of bogs are critically important because their domed morphology4,6,7 accounts for much of the carbon that bogs store and determines how they will respond to interventions8,9 to stop greenhouse gas emissions and fires after anthropogenic drainage10-13. However, a general theory to infer the morphology of bogs is still lacking4,6,7. Here we show that an equation based on the processes universal to bogs explains their morphology across biomes, from Alaska, through the tropics, to New Zealand. In contrast to earlier models of bog morphology that attempted to describe only long-term equilibrium shapes4,6,7 and were, therefore, inapplicable to most bogs14-16, our approach makes no such assumption and makes it possible to infer full shapes of bogs from a sample of elevations, such as a single elevation transect. Our findings provide a foundation for quantitative inference about the morphology, hydrology and carbon storage of bogs through Earth's history, as well as a basis for planning natural climate solutions by rewetting damaged bogs around the world.

One-third of Southern Ocean productivity is supported by dust deposition.

Natural iron fertilization of the Southern Ocean by windblown dust has been suggested to enhance biological productivity and modulate the climate1-3. Yet, this process has never been quantified across the Southern Ocean and at annual timescales4,5. Here we combined 11 years of nitrate observations from autonomous biogeochemical ocean profiling floats with a Southern Hemisphere dust simulation to empirically derive the relationship between dust-iron deposition and annual net community production (ANCP) in the iron-limited Southern Ocean. Using this relationship, we determined the biological response to dust-iron in the pelagic perennially ice-free Southern Ocean at present and during the last glacial maximum (LGM). We estimate that dust-iron now supports 33% ± 15% of Southern Ocean ANCP. During the LGM, when dust deposition was 5-40-fold higher than today, the contribution of dust to Southern Ocean ANCP was much greater, estimated at 64% ± 13%. We provide quantitative evidence of basin-wide dust-iron fertilization of the Southern Ocean and the potential magnitude of its impact on glacial-interglacial timescales, supporting the idea of the important role of dust in the global carbon cycle and climate6-8.

Revisiting the Holocene global temperature conundrum.

Recent global temperature reconstructions for the current interglacial period (the Holocene, beginning 11,700 years ago) have generated contrasting trends. This Review examines evidence from indicators and drivers of global change, as inferred from proxy records and simulated by climate models, to evaluate whether anthropogenic global warming was preceded by a long-term warming trend or by global cooling. Multimillennial-scale cooling before industrialization requires extra climate forcing and major climate feedbacks that are not well represented in most climate models at present. Conversely, global warming before industrialization challenges proxy-based reconstructions of past climate. The resolution of this conundrum has implications for contextualizing post-industrial warming and for understanding climate sensitivity to several forcings and their attendant feedbacks, including greenhouse gases. From a large variety of available evidence, we find support for a relatively mild millennial-scale global thermal maximum during the mid-Holocene, but more research is needed to firmly resolve the conundrum and to advance our understanding of slow-moving climate variability.

Modern temperatures in central-north Greenland warmest in past millennium.

The Greenland Ice Sheet has a central role in the global climate system owing to its size, radiative effects and freshwater storage, and as a potential tipping point1. Weather stations show that the coastal regions are warming2, but the imprint of global warming in the central part of the ice sheet is unclear, owing to missing long-term observations. Current ice-core-based temperature reconstructions3-5 are ambiguous with respect to isolating global warming signatures from natural variability, because they are too noisy and do not include the most recent decades. By systematically redrilling ice cores, we created a high-quality reconstruction of central and north Greenland temperatures from AD 1000 until 2011. Here we show that the warming in the recent reconstructed decade exceeds the range of the pre-industrial temperature variability in the past millennium with virtual certainty (P < 0.001) and is on average 1.5 ± 0.4 degrees Celsius (1 standard error) warmer than the twentieth century. Our findings suggest that these exceptional temperatures arise from the superposition of natural variability with a long-term warming trend, apparent since AD 1800. The disproportionate warming is accompanied by enhanced Greenland meltwater run-off, implying that anthropogenic influence has also arrived in central and north Greenland, which might further accelerate the overall Greenland mass loss.

The geography of climate and the global patterns of species diversity.

Climate's effect on global biodiversity is typically viewed through the lens of temperature, humidity and resulting ecosystem productivity1-6. However, it is not known whether biodiversity depends solely on these climate conditions, or whether the size and fragmentation of these climates are also crucial. Here we shift the common perspective in global biodiversity studies, transitioning from geographic space to a climate-defined multidimensional space. Our findings suggest that larger and more isolated climate conditions tend to harbour higher diversity and species turnover among terrestrial tetrapods, encompassing more than 30,000 species. By considering both the characteristics of climate itself and its geographic attributes, we can explain almost 90% of the variation in global species richness. Half of the explanatory power (45%) may be attributed either to climate itself or to the geography of climate, suggesting a nuanced interplay between them. Our work evolves the conventional idea that larger climate regions, such as the tropics, host more species primarily because of their size7,8. Instead, we underscore the integral roles of both the geographic extent and degree of isolation of climates. This refined understanding presents a more intricate picture of biodiversity distribution, which can guide our approach to biodiversity conservation in an ever-changing world.

Precipitation regime changes in High Mountain Asia driven by cleaner air.

High Mountain Asia (HMA) has experienced a spatial imbalance in water resources in recent decades, partly because of a dipolar pattern of precipitation changes known as South Drying-North Wetting1. These changes can be influenced by both human activities and internal climate variability2,3. Although climate projections indicate a future widespread wetting trend over HMA1,4, the timing and mechanism of the transition from a dipolar to a monopolar pattern remain unknown. Here we demonstrate that the observed dipolar precipitation change in HMA during summer is primarily driven by westerly- and monsoon-associated precipitation patterns. The weakening of the Asian westerly jet, caused by the uneven emission of anthropogenic aerosols, favoured a dipolar precipitation trend from 1951 to 2020. Moreover, the phase transition of the Interdecadal Pacific Oscillation induces an out-of-phase precipitation change between the core region of the South Asian monsoon and southeastern HMA. Under medium- or high-emission scenarios, corresponding to a global warming of 0.6-1.1 °C compared with the present, the dipolar pattern is projected to shift to a monopolar wetting trend in the 2040s. This shift in precipitation patterns is mainly attributed to the intensified jet stream resulting from reduced emissions of anthropogenic aerosols. These findings underscore the importance of considering the impact of aerosol emission reduction in future social planning by policymakers.

Forced changes in the Pacific Walker circulation over the past millennium.

The Pacific Walker circulation (PWC) has an outsized influence on weather and climate worldwide. Yet the PWC response to external forcings is unclear1,2, with empirical data and model simulations often disagreeing on the magnitude and sign of these responses3. Most climate models predict that the PWC will ultimately weaken in response to global warming4. However, the PWC strengthened from 1992 to 2011, suggesting a significant role for anthropogenic and/or volcanic aerosol forcing5, or internal variability. Here we use a new annually resolved, multi-method, palaeoproxy-derived PWC reconstruction ensemble (1200-2000) to show that the 1992-2011 PWC strengthening is anomalous but not unprecedented in the context of the past 800 years. The 1992-2011 PWC strengthening was unlikely to have been a consequence of volcanic forcing and may therefore have resulted from anthropogenic aerosol forcing or natural variability. We find no significant industrial-era (1850-2000) PWC trend, contrasting the PWC weakening simulated by most climate models3. However, an industrial-era shift to lower-frequency variability suggests a subtle anthropogenic influence. The reconstruction also suggests that volcanic eruptions trigger El Niño-like PWC weakening, similar to the response simulated by climate models.

Pathogen spillover driven by rapid changes in bat ecology.

During recent decades, pathogens that originated in bats have become an increasing public health concern. A major challenge is to identify how those pathogens spill over into human populations to generate a pandemic threat1. Many correlational studies associate spillover with changes in land use or other anthropogenic stressors2,3, although the mechanisms underlying the observed correlations have not been identified4. One limitation is the lack of spatially and temporally explicit data on multiple spillovers, and on the connections among spillovers, reservoir host ecology and behaviour and viral dynamics. We present 25 years of data on land-use change, bat behaviour and spillover of Hendra virus from Pteropodid bats to horses in subtropical Australia. These data show that bats are responding to environmental change by persistently adopting behaviours that were previously transient responses to nutritional stress. Interactions between land-use change and climate now lead to persistent bat residency in agricultural areas, where periodic food shortages drive clusters of spillovers. Pulses of winter flowering of trees in remnant forests appeared to prevent spillover. We developed integrative Bayesian network models based on these phenomena that accurately predicted the presence or absence of clusters of spillovers in each of the 25 years. Our long-term study identifies the mechanistic connections between habitat loss, climate and increased spillover risk. It provides a framework for examining causes of bat virus spillover and for developing ecological countermeasures to prevent pandemics.

Natural short-lived halogens exert an indirect cooling effect on climate.

Observational evidence shows the ubiquitous presence of ocean-emitted short-lived halogens in the global atmosphere1-3. Natural emissions of these chemical compounds have been anthropogenically amplified since pre-industrial times4-6, while, in addition, anthropogenic short-lived halocarbons are currently being emitted to the atmosphere7,8. Despite their widespread distribution in the atmosphere, the combined impact of these species on Earth's radiative balance remains unknown. Here we show that short-lived halogens exert a substantial indirect cooling effect at present (-0.13 ± 0.03 watts per square metre) that arises from halogen-mediated radiative perturbations of ozone (-0.24 ± 0.02 watts per square metre), compensated by those from methane (+0.09 ± 0.01 watts per square metre), aerosols (+0.03 ± 0.01 watts per square metre) and stratospheric water vapour (+0.011 ± 0.001 watts per square metre). Importantly, this substantial cooling effect has increased since 1750 by -0.05 ± 0.03 watts per square metre (61 per cent), driven by the anthropogenic amplification of natural halogen emissions, and is projected to change further (18-31 per cent by 2100) depending on climate warming projections and socioeconomic development. We conclude that the indirect radiative effect due to short-lived halogens should now be incorporated into climate models to provide a more realistic natural baseline of Earth's climate system.

A warming-induced reduction in snow fraction amplifies rainfall extremes.

The intensity of extreme precipitation events is projected to increase in a warmer climate1-5, posing a great challenge to water sustainability in natural and built environments. Of particular importance are rainfall (liquid precipitation) extremes owing to their instantaneous triggering of runoff and association with floods6, landslides7-9 and soil erosion10,11. However, so far, the body of literature on intensification of precipitation extremes has not examined the extremes of precipitation phase separately, namely liquid versus solid precipitation. Here we show that the increase in rainfall extremes in high-elevation regions of the Northern Hemisphere is amplified, averaging 15 per cent per degree Celsius of warming-double the rate expected from increases in atmospheric water vapour. We utilize both a climate reanalysis dataset and future model projections to show that the amplified increase is due to a warming-induced shift from snow to rain. Furthermore, we demonstrate that intermodel uncertainty in projections of rainfall extremes can be appreciably explained by changes in snow-rain partitioning (coefficient of determination 0.47). Our findings pinpoint high-altitude regions as 'hotspots' that are vulnerable to future risk of extreme-rainfall-related hazards, thereby requiring robust climate adaptation plans to alleviate potential risk. Moreover, our results offer a pathway towards reducing model uncertainty in projections of rainfall extremes.

Safe and just Earth system boundaries.

The stability and resilience of the Earth system and human well-being are inseparably linked1-3, yet their interdependencies are generally under-recognized; consequently, they are often treated independently4,5. Here, we use modelling and literature assessment to quantify safe and just Earth system boundaries (ESBs) for climate, the biosphere, water and nutrient cycles, and aerosols at global and subglobal scales. We propose ESBs for maintaining the resilience and stability of the Earth system (safe ESBs) and minimizing exposure to significant harm to humans from Earth system change (a necessary but not sufficient condition for justice)4. The stricter of the safe or just boundaries sets the integrated safe and just ESB. Our findings show that justice considerations constrain the integrated ESBs more than safety considerations for climate and atmospheric aerosol loading. Seven of eight globally quantified safe and just ESBs and at least two regional safe and just ESBs in over half of global land area are already exceeded. We propose that our assessment provides a quantitative foundation for safeguarding the global commons for all people now and into the future.

Landscape dynamics and the Phanerozoic diversification of the biosphere.

The long-term diversification of the biosphere responds to changes in the physical environment. Yet, over the continents, the nearly monotonic expansion of life started later in the early part of the Phanerozoic eon1 than the expansion in the marine realm, where instead the number of genera waxed and waned over time2. A comprehensive evaluation of the changes in the geodynamic and climatic forcing fails to provide a unified theory for the long-term pattern of evolution of life on Earth. Here we couple climate and plate tectonics models to numerically reconstruct the evolution of the Earth's landscape over the entire Phanerozoic eon, which we then compare to palaeo-diversity datasets from marine animal and land plant genera. Our results indicate that biodiversity is strongly reliant on landscape dynamics, which at all times determine the carrying capacity of both the continental domain and the oceanic domain. In the oceans, diversity closely adjusted to the riverine sedimentary flux that provides nutrients for primary production. On land, plant expansion was hampered by poor edaphic conditions until widespread endorheic basins resurfaced continents with a sedimentary cover that facilitated the development of soil-dependent rooted flora, and the increasing variety of the landscape additionally promoted their development.

Increased heat risk in wet climate induced by urban humid heat.

Cities are generally warmer than their adjacent rural land, a phenomenon known as the urban heat island (UHI). Often accompanying the UHI effect is another phenomenon called the urban dry island (UDI), whereby the humidity of urban land is lower than that of the surrounding rural land1-3. The UHI exacerbates heat stress on urban residents4,5, whereas the UDI may instead provide relief because the human body can cope with hot conditions better at lower humidity through perspiration6,7. The relative balance between the UHI and the UDI-as measured by changes in the wet-bulb temperature (Tw)-is a key yet largely unknown determinant of human heat stress in urban climates. Here we show that Tw is reduced in cities in dry and moderately wet climates, where the UDI more than offsets the UHI, but increased in wet climates (summer precipitation of more than 570 millimetres). Our results arise from analysis of urban and rural weather station data across the world and calculations with an urban climate model. In wet climates, the urban daytime Tw is 0.17 ± 0.14 degrees Celsius (mean ± 1 standard deviation) higher than rural Tw in the summer, primarily because of a weaker dynamic mixing in urban air. This Tw increment is small, but because of the high background Tw in wet climates, it is enough to cause two to six extra dangerous heat-stress days per summer for urban residents under current climate conditions. The risk of extreme humid heat is projected to increase in the future, and these urban effects may further amplify the risk.

Topography of mountain belts controlled by rheology and surface processes.

It is widely recognized that collisional mountain belt topography is generated by crustal thickening and lowered by river bedrock erosion, linking climate and tectonics1-4. However, whether surface processes or lithospheric strength control mountain belt height, shape and longevity remains uncertain. Additionally, how to reconcile high erosion rates in some active orogens with long-term survival of mountain belts for hundreds of millions of years remains enigmatic. Here we investigate mountain belt growth and decay using a new coupled surface process5,6 and mantle-scale tectonic model7. End-member models and the new non-dimensional Beaumont number, Bm, quantify how surface processes and tectonics control the topographic evolution of mountain belts, and enable the definition of three end-member types of growing orogens: type 1, non-steady state, strength controlled (Bm > 0.5); type 2, flux steady state8, strength controlled (Bm ≈ 0.4-0.5); and type 3, flux steady state, erosion controlled (Bm < 0.4). Our results indicate that tectonics dominate in Himalaya-Tibet and the Central Andes (both type 1), efficient surface processes balance high convergence rates in Taiwan (probably type 2) and surface processes dominate in the Southern Alps of New Zealand (type 3). Orogenic decay is determined by erosional efficiency and can be subdivided into two phases with variable isostatic rebound characteristics and associated timescales. The results presented here provide a unified framework explaining how surface processes and lithospheric strength control the height, shape, and longevity of mountain belts.

Life rather than climate influences diversity at scales greater than 40 million years.

The diversity of life on Earth is controlled by hierarchical processes that interact over wide ranges of timescales1. Here, we consider the megaclimate regime2 at scales ≥1 million years (Myr). We focus on determining the domains of 'wandering' stochastic Earth system processes ('Court Jester'3) and stabilizing biotic interactions that induce diversity dependence of fluctuations in macroevolutionary rates ('Red Queen'4). Using state-of-the-art multiscale Haar and cross-Haar fluctuation analyses, we analysed the global genus-level Phanerozoic marine animal Paleobiology Database record of extinction rates (E), origination rates (O) and diversity (D) as well as sea water palaeotemperatures (T). Over the entire observed range from several million years to several hundred million years, we found that the fluctuations of T, E and O showed time-scaling behaviour. The megaclimate was characterized by positive scaling exponents-it is therefore apparently unstable. E and O are also scaling but with negative exponents-stable behaviour that is biotically mediated. For D, there were two regimes with a crossover at critical timescale [Formula: see text] ≈ 40 Myr. For shorter timescales, D exhibited nearly the same positive scaling as the megaclimate palaeotemperatures, whereas for longer timescales it tracks the scaling of macroevolutionary rates. At scales of at least [Formula: see text] there is onset of diversity dependence of E and O, probably enabled by mixing and synchronization (globalization) of the biota by geodispersal ('Geo-Red Queen').

Small rainfall changes drive substantial changes in plant coexistence.

Although precipitation patterns have long been known to shape plant distributions1, the effect of changing climate on the interactions of species and therefore community composition is far less understood2,3. Here, we explored how changes in precipitation alter competitive dynamics via direct effects on individual species, as well as by the changing strength of competitive interactions between species, using an annual grassland community in California. We grew plants under ambient and reduced precipitation in the field to parameterize a competition model4 with which we quantified the stabilizing niche and fitness differences that determine species coexistence in each rainfall regime. We show that reduced precipitation had little direct effect on species grown alone, but it qualitatively shifted predicted competitive outcomes for 10 of 15 species pairs. In addition, species pairs that were functionally more similar were less likely to experience altered outcomes, indicating that functionally diverse communities may be most threatened by changing interactions. Our results highlight how important it is to account for changes to species interactions when predicting species and community response to global change.

Africa's oldest dinosaurs reveal early suppression of dinosaur distribution.

The vertebrate lineages that would shape Mesozoic and Cenozoic terrestrial ecosystems originated across Triassic Pangaea1-11. By the Late Triassic (Carnian stage, ~235 million years ago), cosmopolitan 'disaster faunas' (refs. 12-14) had given way to highly endemic assemblages12,13 on the supercontinent. Testing the tempo and mode of the establishment of this endemism is challenging-there were few geographic barriers to dispersal across Pangaea during the Late Triassic. Instead, palaeolatitudinal climate belts, and not continental boundaries, are proposed to have controlled distribution15-18. During this time of high endemism, dinosaurs began to disperse and thus offer an opportunity to test the timing and drivers of this biogeographic pattern. Increased sampling can test this prediction: if dinosaurs initially dispersed under palaeolatitudinal-driven endemism, then an assemblage similar to those of South America4,19-21 and India19,22-including the earliest dinosaurs-should be present in Carnian deposits in south-central Africa. Here we report a new Carnian assemblage from Zimbabwe that includes Africa's oldest definitive dinosaurs, including a nearly complete skeleton of the sauropodomorph Mbiresaurus raathi gen. et sp. nov. This assemblage resembles other dinosaur-bearing Carnian assemblages, suggesting that a similar vertebrate fauna ranged high-latitude austral Pangaea. The distribution of the first dinosaurs is correlated with palaeolatitude-linked climatic barriers, and dinosaurian dispersal to the rest of the supercontinent was delayed until these barriers relaxed, suggesting that climatic controls influenced the initial composition of the terrestrial faunas that persist to this day.

Comprehensive evidence implies a higher social cost of CO2.

The social cost of carbon dioxide (SC-CO2) measures the monetized value of the damages to society caused by an incremental metric tonne of CO2 emissions and is a key metric informing climate policy. Used by governments and other decision-makers in benefit-cost analysis for over a decade, SC-CO2 estimates draw on climate science, economics, demography and other disciplines. However, a 2017 report by the US National Academies of Sciences, Engineering, and Medicine1 (NASEM) highlighted that current SC-CO2 estimates no longer reflect the latest research. The report provided a series of recommendations for improving the scientific basis, transparency and uncertainty characterization of SC-CO2 estimates. Here we show that improved probabilistic socioeconomic projections, climate models, damage functions, and discounting methods that collectively reflect theoretically consistent valuation of risk, substantially increase estimates of the SC-CO2. Our preferred mean SC-CO2 estimate is $185 per tonne of CO2 ($44-$413 per tCO2: 5%-95% range, 2020 US dollars) at a near-term risk-free discount rate of 2%, a value 3.6 times higher than the US government's current value of $51 per tCO2. Our estimates incorporate updated scientific understanding throughout all components of SC-CO2 estimation in the new open-source Greenhouse Gas Impact Value Estimator (GIVE) model, in a manner fully responsive to the near-term NASEM recommendations. Our higher SC-CO2 values, compared with estimates currently used in policy evaluation, substantially increase the estimated benefits of greenhouse gas mitigation and thereby increase the expected net benefits of more stringent climate policies.

Deciphering and prediction of transcriptome dynamics under fluctuating field conditions.

Determining the drivers of gene expression patterns is more straightforward in laboratory conditions than in the complex fluctuating environments where organisms typically live. We gathered transcriptome data from the leaves of rice plants in a paddy field along with the corresponding meteorological data and used them to develop statistical models for the endogenous and external influences on gene expression. Our results indicate that the transcriptome dynamics are predominantly governed by endogenous diurnal rhythms, ambient temperature, plant age, and solar radiation. The data revealed diurnal gates for environmental stimuli to influence transcription and pointed to relative influences exerted by circadian and environmental factors on different metabolic genes. The model also generated predictions for the influence of changing temperatures on transcriptome dynamics. We anticipate that our models will help translate the knowledge amassed in laboratories to problems in agriculture and that our approach to deciphering the transcriptome fluctuations in complex environments will be applicable to other organisms.