Last Glacial Maximum pattern effects reduce climate sensitivity estimates.

Here, we show that the Last Glacial Maximum (LGM) provides a stronger constraint on equilibrium climate sensitivity (ECS), the global warming from increasing greenhouse gases, after accounting for temperature patterns. Feedbacks governing ECS depend on spatial patterns of surface temperature ("pattern effects"); hence, using the LGM to constrain future warming requires quantifying how temperature patterns produce different feedbacks during LGM cooling versus modern-day warming. Combining data assimilation reconstructions with atmospheric models, we show that the climate is more sensitive to LGM forcing because ice sheets amplify extratropical cooling where feedbacks are destabilizing. Accounting for LGM pattern effects yields a median modern-day ECS of 2.4°C, 66% range 1.7° to 3.5°C (1.4° to 5.0°C, 5 to 95%), from LGM evidence alone. Combining the LGM with other lines of evidence, the best estimate becomes 2.9°C, 66% range 2.4° to 3.5°C (2.1° to 4.1°C, 5 to 95%), substantially narrowing uncertainty compared to recent assessments.

Climate system asymmetries drive eccentricity pacing of hydroclimate during the early Eocene greenhouse.

The early Eocene Climatic Optimum (EECO) represents the peak of Earth's last sustained greenhouse climate interval. To investigate hydroclimate variability in western North America during the EECO, we developed an orbitally resolved leaf wax δ2H record from one of the most well-dated terrestrial paleoclimate archives, the Green River Formation. Our δ2Hwax results show ∼60‰ variation and evidence for eccentricity and precession forcing. iCESM simulations indicate that changes in the Earth's orbit drive large seasonal variations in precipitation and δ2H of precipitation at our study site, primarily during the summer season. Our findings suggest that the astronomical response in δ2Hwax is attributable to an asymmetrical climate response to the seasonal cycle, a "clipping" of precession forcing, and asymmetric carbon cycle dynamics, which further enhance the influence of eccentricity modulation on the hydrological cycle during the EECO. More broadly, our study provides an explanation for how and why eccentricity emerges as a dominant frequency in climate records from ice-free greenhouse worlds.

The PhanSST global database of Phanerozoic sea surface temperature proxy data.

Paleotemperature proxy data form the cornerstone of paleoclimate research and are integral to understanding the evolution of the Earth system across the Phanerozoic Eon. Here, we present PhanSST, a database containing over 150,000 data points from five proxy systems that can be used to estimate past sea surface temperature. The geochemical data have a near-global spatial distribution and temporally span most of the Phanerozoic. Each proxy value is associated with consistent and queryable metadata fields, including information about the location, age, and taxonomy of the organism from which the data derive. To promote transparency and reproducibility, we include all available published data, regardless of interpreted preservation state or vital effects. However, we also provide expert-assigned diagenetic assessments, ecological and environmental flags, and other proxy-specific fields, which facilitate informed and responsible reuse of the database. The data are quality control checked and the foraminiferal taxonomy has been updated. PhanSST will serve as a valuable resource to the paleoclimate community and has myriad applications, including evolutionary, geochemical, diagenetic, and proxy calibration studies.

Spatial patterns of climate change across the Paleocene-Eocene Thermal Maximum.

The Paleocene-Eocene Thermal Maximum (PETM; 56 Ma) is one of our best geological analogs for understanding climate dynamics in a "greenhouse" world. However, proxy data representing the event are only available from select marine and terrestrial sedimentary sequences that are unevenly distributed across Earth's surface, limiting our view of the spatial patterns of climate change. Here, we use paleoclimate data assimilation (DA) to combine climate model and proxy information and create a spatially complete reconstruction of the PETM and the climate state that precedes it ("PETM-DA"). Our data-constrained results support strong polar amplification, which in the absence of an extensive cryosphere, is related to temperature feedbacks and loss of seasonal snow on land. The response of the hydrological cycle to PETM warming consists of a narrowing of the Intertropical Convergence Zone, off-equatorial drying, and an intensification of seasonal monsoons and winter storm tracks. Many of these features are also seen in simulations of future climate change under increasing anthropogenic emissions. Since the PETM-DA yields a spatially complete estimate of surface air temperature, it yields a rigorous estimate of global mean temperature change (5.6 ∘C; 5.4 ∘C to 5.9 ∘C, 95% CI) that can be used to calculate equilibrium climate sensitivity (ECS). We find that PETM ECS was 6.5 ∘C (5.7 ∘C to 7.4 ∘C, 95% CI), which is much higher than the present-day range. This supports the view that climate sensitivity increases substantially when greenhouse gas concentrations are high.

Northern Hemisphere vegetation change drives a Holocene thermal maximum.

The Holocene thermal maximum, a period of global warmth evident in early to mid-Holocene proxy reconstructions, is controversial. Most model simulations of the Holocene have not reproduced this warming, leading to a disagreement known as the Holocene Temperature Conundrum. Pollen records document the expansion of vegetation in the early and mid-Holocene African Sahara and Northern Hemisphere mid- and high latitudes, which has been overlooked in previous modeling studies. Here, we use time slice simulations of the Community Earth System Model to assess the impact of Northern Hemisphere vegetation change on Holocene annual mean temperatures. Our simulations indicate that expansion of Northern Hemisphere vegetation 9000 and 6000 years ago warms Earth's surface by ~0.8° and 0.7°C, respectively, producing a better match with proxy-based reconstructions. Our results suggest that vegetation change is critical for modeling Holocene temperature evolution and highlight its role in driving a mid-Holocene temperature maximum.

Globally resolved surface temperatures since the Last Glacial Maximum.

Climate changes across the past 24,000 years provide key insights into Earth system responses to external forcing. Climate model simulations1,2 and proxy data3-8 have independently allowed for study of this crucial interval; however, they have at times yielded disparate conclusions. Here, we leverage both types of information using paleoclimate data assimilation9,10 to produce the first proxy-constrained, full-field reanalysis of surface temperature change spanning the Last Glacial Maximum to present at 200-year resolution. We demonstrate that temperature variability across the past 24 thousand years was linked to two primary climatic mechanisms: radiative forcing from ice sheets and greenhouse gases; and a superposition of changes in the ocean overturning circulation and seasonal insolation. In contrast with previous proxy-based reconstructions6,7 our results show that global mean temperature has slightly but steadily warmed, by ~0.5 °C, since the early Holocene (around 9 thousand years ago). When compared with recent temperature changes11, our reanalysis indicates that both the rate and magnitude of modern warming are unusual relative to the changes of the past 24 thousand years.

Subseafloor Archaea reflect 139 kyrs of paleodepositional changes in the northern Red Sea.

The vertical distribution of subseafloor archaeal communities is thought to be primarily controlled by in situ conditions in sediments such as the availability of electron acceptors and donors, although sharp community shifts have also been observed at lithological boundaries suggesting that at least a subset of vertically stratified Archaea form a long-term genetic record of coinciding environmental conditions that occurred at the time of sediment deposition. To substantiate this possibility, we performed a highly resolved 16S rRNA gene survey of vertically stratified archaeal communities paired with paleo-oceanographic proxies in a sedimentary record from the northern Red Sea spanning the last glacial-interglacial cycle (i.e., marine isotope stages 1-6; MIS1-6). Our results show a strong significant correlation between subseafloor archaeal communities and drastic paleodepositional changes associated with glacial low vs. interglacial high stands (ANOSIM; R = .73; p = .001) and only a moderately strong correlation with lithological changes. Bathyarchaeota, Lokiarchaeota, MBGA, and DHVEG-1 were the most abundant identified archaeal groups. Whether they represented ancient cell lines from the time of deposition or migrated to the specific sedimentary horizons after deposition remains speculative. However, we show that the majority of sedimentary archaeal tetraether membrane lipids were of allochthonous origin and not produced in situ. Slow post-burial growth under energy-limited conditions would explain why the downcore distribution of these dominant archaeal groups still indirectly reflect changes in the paleodepositional environment that prevailed during the analyzed marine isotope stages. In addition, archaea seeded from the overlying water column such as Thaumarchaeota and group II and III Euryarchaeota, which were likely not have been able to subsist after burial, were identified from a lower abundance of preserved sedimentary DNA signatures, and represented direct markers of paleoenvironmental changes in the Red Sea spanning the last six marine isotope stages.

Past climates inform our future.

As the world warms, there is a profound need to improve projections of climate change. Although the latest Earth system models offer an unprecedented number of features, fundamental uncertainties continue to cloud our view of the future. Past climates provide the only opportunity to observe how the Earth system responds to high carbon dioxide, underlining a fundamental role for paleoclimatology in constraining future climate change. Here, we review the relevancy of paleoclimate information for climate prediction and discuss the prospects for emerging methodologies to further insights gained from past climates. Advances in proxy methods and interpretations pave the way for the use of past climates for model evaluation-a practice that we argue should be widely adopted.

Glacial cooling and climate sensitivity revisited.

The Last Glacial Maximum (LGM), one of the best studied palaeoclimatic intervals, offers an excellent opportunity to investigate how the climate system responds to changes in greenhouse gases and the cryosphere. Previous work has sought to constrain the magnitude and pattern of glacial cooling from palaeothermometers1,2, but the uneven distribution of the proxies, as well as their uncertainties, has challenged the construction of a full-field view of the LGM climate state. Here we combine a large collection of geochemical proxies for sea surface temperature with an isotope-enabled climate model ensemble to produce a field reconstruction of LGM temperatures using data assimilation. The reconstruction is validated with withheld proxies as well as independent ice core and speleothem δ18O measurements. Our assimilated product provides a constraint on global mean LGM cooling of -6.1 degrees Celsius (95 per cent confidence interval: -6.5 to -5.7 degrees Celsius). Given assumptions concerning the radiative forcing of greenhouse gases, ice sheets and mineral dust aerosols, this cooling translates to an equilibrium climate sensitivity of 3.4 degrees Celsius (2.4-4.5 degrees Celsius), a value that is higher than previous LGM-based estimates but consistent with the traditional consensus range of 2-4.5 degrees Celsius3,4.

Emergence of an equatorial mode of climate variability in the Indian Ocean.

Presently, the Indian Ocean (IO) resides in a climate state that prevents strong year-to-year climate variations. This may change under greenhouse warming, but the mechanisms remain uncertain, thus limiting our ability to predict future changes in climate extremes. Using climate model simulations, we uncover the emergence of a mode of climate variability capable of generating unprecedented sea surface temperature and rainfall fluctuations across the IO. This mode, which is inhibited under present-day conditions, becomes active in climate states with a shallow thermocline and vigorous upwelling, consistent with the predictions of continued greenhouse warming. These predictions are supported by modeling and proxy evidence of an active mode during glacial intervals that favored such a state. Because of its impact on hydrological variability, the emergence of such a mode would become a first-order source of climate-related risks for the densely populated IO rim.

Simulation of Eocene extreme warmth and high climate sensitivity through cloud feedbacks.

The Early Eocene, a period of elevated atmospheric CO2 (>1000 ppmv), is considered an analog for future climate. Previous modeling attempts have been unable to reproduce major features of Eocene climate indicated by proxy data without substantial modification to the model physics. Here, we present simulations using a state-of-the-art climate model forced by proxy-estimated CO2 levels that capture the extreme surface warmth and reduced latitudinal temperature gradient of the Early Eocene and the warming of the Paleocene-Eocene Thermal Maximum. Our simulations exhibit increasing equilibrium climate sensitivity with warming and suggest an Eocene sensitivity of more than 6.6°C, much greater than the present-day value (4.2°C). This higher climate sensitivity is mainly attributable to the shortwave cloud feedback, which is linked primarily to cloud microphysical processes. Our findings highlight the role of small-scale cloud processes in determining large-scale climate changes and suggest a potential increase in climate sensitivity with future warming.

Glacial changes in tropical climate amplified by the Indian Ocean.

The mechanisms driving glacial-interglacial changes in the climate of the Indo-Pacific warm pool are poorly understood. Here, we address this question by combining paleoclimate proxies with model simulations of the Last Glacial Maximum climate. We find evidence of two mechanisms explaining key patterns of ocean cooling and rainfall change interpreted from proxy data. Exposure of the Sahul shelf excites a positive ocean-atmosphere feedback involving a stronger surface temperature gradient along the equatorial Indian Ocean and a weaker Walker circulation-a response explaining the drier/wetter dipole across the basin. Northern Hemisphere cooling by ice sheet albedo drives a monsoonal retreat across Africa and the Arabian Peninsula-a response that triggers a weakening of the Indian monsoon via cooling of the Arabian Sea and associated reductions in moisture supply. These results demonstrate the importance of air-sea interactions in the Indian Ocean, amplifying externally forced climate changes over a large part of the tropics.

Comparison of three methods for the methylation of aliphatic and aromatic compounds.

RATIONALE: Methylation protocols commonly call for acidic, hot conditions that are known to promote organic 1 H/2 H exchange in aromatic and aliphatic C-H bonds. Here we tested two such commonly used methods and compared a third that avoids these acidic conditions, to quantify isotope effects with each method and to directly determine acidic-exchange rates relevant to experimental conditions. METHODS: We compared acidic and non-acidic methylation approaches catalyzed by hydrochloric acid, acetyl chloride and EDCI (1-ethyl-3-(3-dimethylaminopropyl)carbodiimide)/DMAP (4-dimethylaminopyridine), respectively. These were applied to two analytes: phthalic acid (an aromatic) and octacosanoic acid (an aliphatic). We analyzed yield by gas chromatography/flame ionization (GC/FID) and hydrogen and carbon isotopic compositions by isotope ratio mass spectrometry (GC/IRMS). We quantified the 1 H/2 H exchange rate on dimethyl phthalate under acidic conditions with proton nuclear magnetic resonance (1 H-NMR) measurements. RESULTS: The δ2 H and δ13 C values and yield were equivalent among the three methods for methyl octacosanoate. The two acidic methods resulted in comparable yield and isotopic composition of dimethyl phthalate; however, the non-acidic method resulted in lower δ2 H and δ13 C values perhaps due to low yields. Concerns over acid-catalyzed 1 H/2 H exchange are unwarranted as the effect was trivial over a 12-h reaction time. CONCLUSIONS: We find product isolation yield and evaporation to be the main concerns in the accurate determination of isotopic composition. 1 H/2 H exchange reactions are too slow to cause measurable isotope fractionation over the typical duration and reaction conditions used in methylation. Thus, we are able to recommend continued use of acidic catalysts in such methylation reactions for both aliphatic and aromatic compounds.

Rainfall regimes of the Green Sahara.

During the "Green Sahara" period (11,000 to 5000 years before the present), the Sahara desert received high amounts of rainfall, supporting diverse vegetation, permanent lakes, and human populations. Our knowledge of rainfall rates and the spatiotemporal extent of wet conditions has suffered from a lack of continuous sedimentary records. We present a quantitative reconstruction of western Saharan precipitation derived from leaf wax isotopes in marine sediments. Our data indicate that the Green Sahara extended to 31°N and likely ended abruptly. We find evidence for a prolonged "pause" in Green Sahara conditions 8000 years ago, coincident with a temporary abandonment of occupational sites by Neolithic humans. The rainfall rates inferred from our data are best explained by strong vegetation and dust feedbacks; without these mechanisms, climate models systematically fail to reproduce the Green Sahara. This study suggests that accurate simulations of future climate change in the Sahara and Sahel will require improvements in our ability to simulate vegetation and dust feedbacks.

Early onset of industrial-era warming across the oceans and continents.

The evolution of industrial-era warming across the continents and oceans provides a context for future climate change and is important for determining climate sensitivity and the processes that control regional warming. Here we use post-ad 1500 palaeoclimate records to show that sustained industrial-era warming of the tropical oceans first developed during the mid-nineteenth century and was nearly synchronous with Northern Hemisphere continental warming. The early onset of sustained, significant warming in palaeoclimate records and model simulations suggests that greenhouse forcing of industrial-era warming commenced as early as the mid-nineteenth century and included an enhanced equatorial ocean response mechanism. The development of Southern Hemisphere warming is delayed in reconstructions, but this apparent delay is not reproduced in climate simulations. Our findings imply that instrumental records are too short to comprehensively assess anthropogenic climate change and that, in some regions, about 180 years of industrial-era warming has already caused surface temperatures to emerge above pre-industrial values, even when taking natural variability into account.

Changes in northeast African hydrology and vegetation associated with Pliocene-Pleistocene sapropel cycles.

East African climate change since the Late Miocene consisted of persistent shorter-term, orbital-scale wet-dry cycles superimposed upon a long-term trend towards more open, grassy landscapes. Either or both of these modes of palaeoclimate variability may have influenced East African mammalian evolution, yet the interrelationship between these secular and orbital palaeoclimate signals remains poorly understood. Here, we explore whether the long-term secular climate change was also accompanied by significant changes at the orbital-scale. We develop northeast African hydroclimate and vegetation proxy data for two 100 kyr-duration windows near 3.05 and 1.75 Ma at ODP Site 967 in the eastern Mediterranean basin, where sedimentation is dominated by eastern Sahara dust input and Nile River run-off. These two windows were selected because they have comparable orbital configurations and bracket an important increase in East African C4 grasslands. We conducted high-resolution (2.5 kyr sampling) multiproxy biomarker, H- and C-isotopic analyses of plant waxes and lignin phenols to document orbital-scale changes in hydrology, vegetation and woody cover for these two intervals. Both intervals are dominated by large-amplitude, precession-scale (approx. 20 kyr) changes in northeast African vegetation and rainfall/run-off. The δ(13)Cwax values and lignin phenol composition record a variable but consistently C4 grass-dominated ecosystem for both intervals (50-80% C4). Precessional δDwax cycles were approximately 20-30‰ in peak-to-peak amplitude, comparable with other δDwax records of the Early Holocene African Humid Period. There were no significant differences in the means or variances of the δDwax or δ(13)Cwax data for the 3.05 and 1.75 Ma intervals studied, suggesting that the palaeohydrology and palaeovegetation responses to precessional forcing were similar for these two periods. Data for these two windows suggest that the eastern Sahara did not experience the significant increase in C4 vegetation that has been observed in East Africa over this time period. This observation would be consistent with a proposed mechanism whereby East African precipitation is reduced, and drier conditions established, in response to the emergence of modern zonal sea surface temperature gradients in the tropical oceans between 3 and 2 Ma.This article is part of the themed issue 'Major transitions in human evolution'.

Past and future rainfall in the Horn of Africa.

The recent decline in Horn of Africa rainfall during the March-May "long rains" season has fomented drought and famine, threatening food security in an already vulnerable region. Some attribute this decline to anthropogenic forcing, whereas others maintain that it is a feature of internal climate variability. We show that the rate of drying in the Horn of Africa during the 20th century is unusual in the context of the last 2000 years, is synchronous with recent global and regional warming, and therefore may have an anthropogenic component. In contrast to 20th century drying, climate models predict that the Horn of Africa will become wetter as global temperatures rise. The projected increase in rainfall mainly occurs during the September-November "short rains" season, in response to large-scale weakening of the Walker circulation. Most of the models overestimate short rains precipitation while underestimating long rains precipitation, causing the Walker circulation response to unrealistically dominate the annual mean. Our results highlight the need for accurate simulation of the seasonal cycle and an improved understanding of the dynamics of the long rains season to predict future rainfall in the Horn of Africa.

A TEX86 surface sediment database and extended Bayesian calibration.

Quantitative estimates of past temperature changes are a cornerstone of paleoclimatology. For a number of marine sediment-based proxies, the accuracy and precision of past temperature reconstructions depends on a spatial calibration of modern surface sediment measurements to overlying water temperatures. Here, we present a database of 1095 surface sediment measurements of TEX86, a temperature proxy based on the relative cyclization of marine archaeal glycerol dialkyl glycerol tetraether (GDGT) lipids. The dataset is archived in a machine-readable format with geospatial information, fractional abundances of lipids (if available), and metadata. We use this new database to update surface and subsurface temperature calibration models for TEX86 and demonstrate the applicability of the TEX86 proxy to past temperature prediction. The TEX86 database confirms that surface sediment GDGT distribution has a strong relationship to temperature, which accounts for over 70% of the variance in the data. Future efforts, made possible by the data presented here, will seek to identify variables with secondary relationships to GDGT distributions, such as archaeal community composition.

Abrupt shifts in Horn of Africa hydroclimate since the Last Glacial Maximum.

The timing and abruptness of the initiation and termination of the Early Holocene African Humid Period are subjects of ongoing debate, with direct consequences for our understanding of abrupt climate change, paleoenvironments, and early human cultural development. Here, we provide proxy evidence from the Horn of Africa region that documents abrupt transitions into and out of the African Humid Period in northeast Africa. Similar and generally synchronous abrupt transitions at other East African sites suggest that rapid shifts in hydroclimate are a regionally coherent feature. Our analysis suggests that the termination of the African Humid Period in the Horn of Africa occurred within centuries, underscoring the nonlinearity of the region's hydroclimate.