ALBERT

Description: The importance of Antarctic sea ice and Southern Ocean warming has come into the focus of polar research during the last couple of decades. Especially around West Antarctica, where warm water masses approach the continent and where sea ice has declined, the distribution and evolution of sea ice play a critical role in the stability of nearby ice shelves. Organic geochemical analyses of marine seafloor surface sediments from the Antarctic continental margin allow an evaluation of the applicability of biomarker-based sea-ice and ocean temperature reconstructions in these climate-sensitive areas. We analysed highly branched isoprenoids (HBIs), such as the sea-ice proxy IPSO25 and phytoplankton-derived HBI-trienes, as well as phytosterols and isoprenoidal glycerol dialkyl glycerol tetraethers (GDGTs), which are established tools for the assessment of primary productivity and ocean temperatures respectively. The combination of IPSO25 with a phytoplankton marker (i.e. the PIPSO25 index) permits semi-quantitative sea-ice reconstructions and avoids misleading over- or underestimations of sea-ice cover. Comparisons of the PIPSO25-based sea-ice distribution patterns and TEX86L- and RI-OH′-derived ocean temperatures with (1) sea-ice concentrations obtained from satellite observations and (2) instrument measurements of sea surface and subsurface temperatures corroborate the general capability of these proxies to determine oceanic key variables properly. This is further supported by model data. We also highlight specific aspects and limitations that need to be taken into account for the interpretation of such biomarker data and discuss the potential of IPSO25 as an indicator for the former occurrence of platelet ice and/or the export of ice-shelf water.

Description: Frequently occurring mega-droughts under current global climate change have attracted broad social attention. A paleoclimatic perspective is needed to increase our understanding of the causes and effects of droughts. South-western (SW) China has been threatened by severe seasonal droughts. Our current knowledge of millennial-scale dry and wet phases in this region is primarily based on the variability of the Indian summer monsoon. However, water availability over land does not always follow patterns of monsoonal precipitation but also depends on water loss from evaporation and transpiration. Here, we reconstructed precipitation intensity, lake hydrological balance and the soil water stress index (SWSI) for the last 27 000 years. Grain size, geochemical and pollen records from Yilong Lake reveal the long-term relationships and inconsistencies of dry–wet patterns in meteorological, hydrological and soil systems in the central Yunnan region, SW China. Our results show that the long-term trends among precipitation, hydrological balance and soil moisture varied through time. The hydrological balance and soil moisture were primarily controlled by temperature-induced evaporation change during periods of low precipitation such as the Last Glacial Maximum and Younger Dryas. During periods of high precipitation (the early to late Holocene), intensified evaporation from the lake surface offset the effects of increased precipitation on the hydrological balance. However, abundant rainfall and the dense vegetation canopy circumvented a soil moisture deficit that might have resulted from rising temperature. In conclusion, the hydrological balance in the central Yunnan region was more sensitive to temperature change while soil moisture could be further regulated by vegetation changes over millennial timescales. Therefore, under future climate warming, the surface water shortage in the central Yunnan region may become even more serious. Our study suggests that reforestation efforts may provide some relief to soil moisture deficits in this region.

Description: Deglaciations are characterized by relatively fast and near-synchronous changes in ice sheet volume, ocean temperature, and atmospheric greenhouse gas concentrations, but glacial inception occurs more gradually. Understanding the evolution of ice sheet, ocean, and atmosphere conditions from interglacial to glacial maximum provides insight into the interplay of these components of the climate system. Using noble gas measurements in ancient ice samples, we reconstruct mean ocean temperature (MOT) from 74 to 59.7 ka, covering the Marine Isotope Stage (MIS) 5a–4 boundary, MIS 4, and part of the MIS 4–3 transition. Comparing this MOT reconstruction to previously published MOT reconstructions from the last and penultimate deglaciation, we find that the majority of the last interglacial–glacial ocean cooling must have occurred within MIS 5. MOT reached equally cold conditions in MIS 4 as in MIS 2 (−2.7 ± 0.3 ∘C relative to the Holocene, −0.1 ± 0.3 ∘C relative to MIS 2). Using a carbon cycle model to quantify the CO2 solubility pump, we show that ocean cooling can explain most of the CO2 drawdown (32 ± 4 of 40 ppm) across MIS 5. Comparing MOT to contemporaneous records of benthic δ18O, we find that ocean cooling can also explain the majority of the δ18O increase across MIS 5 (0.7 ‰ of 1.3 ‰). The timing of ocean warming and cooling in the record and the comparison to coeval Antarctic isotope data suggest an intimate link between ocean heat content, Southern Hemisphere high-latitude climate, and ocean circulation on orbital and millennial timescales.

Description: The middle Miocene climate transition (MMCT), around 14 Ma, was associated with a significant climatic shift, but the mechanisms triggering the event remain enigmatic. We present a clumped isotope (Δ47) bottom-water temperature (BWT) record from 16.0 to 12.2 Ma from Ocean Drilling Program (ODP) Site 747 in the Southern Ocean and compare it to existing BWT records from different latitudes. We show that BWTs in the Southern Ocean reached 8–10 ∘C during the Miocene climatic optimum. These high BWT values indicate considerably warmer bottom-water conditions than today. Nonetheless, bottom-water δ18O (calculated from foraminiferal δ18O and Δ47) suggests substantial amounts of land ice throughout the interval of the study. Our dataset further demonstrates that BWTs at Site 747 were variable with an overall cooling trend across the MMCT. Notably, a cooling of around 3–5 ∘C preceded the stepped main increase in benthic δ18O, interpreted as global ice volume expansion, and appears to have been followed by a transient bottom-water warming starting during or slightly after the main ice volume increase. We speculate that a regional freshening of the upper water column at this time may have increased stratification and reduced bottom-water heat loss to the atmosphere, counteracting global cooling in the bottom waters of the Southern Ocean and possibly even at larger scales. Feedbacks required for substantial ice growth and/or tectonic processes may have contributed to the observed decoupling of global ice volume and Southern Ocean BWT.

Description: The last deglaciation offers an unique opportunity to understand the climate–ice-sheet interactions in a global warming context. In this paper, to tackle this question, we use an Earth system model of intermediate complexity coupled to an ice sheet model covering the Northern Hemisphere to simulate the last deglaciation and the Holocene (26–0 ka). We use a synchronous coupling every year between the ice sheet and the rest of the climate system and we ensure a closed water cycle considering the release of freshwater flux to the ocean due to ice sheet melting. Our reference experiment displays a gradual warming in response to the forcings, with no abrupt changes. In this case, while the amplitude of the freshwater flux to the ocean induced by ice sheet retreat is realistic, it is sufficient to shut down the Atlantic meridional overturning circulation from which the model does not recover within the time period simulated. However, with reduced freshwater flux we are nonetheless able to obtain different oceanic circulation evolutions, including some abrupt transitions between shut-down and active circulation states in the course of the deglaciation. The inclusion of a parameterisation for the sinking of brines around Antarctica also produces an abrupt recovery of the Atlantic meridional overturning circulation, absent in the reference experiment. The fast oceanic circulation recoveries lead to abrupt warming phases in Greenland. Our simulated ice sheet geometry evolution is in overall good agreement with available global reconstructions, even though the abrupt sea level rise at 14.6 ka is underestimated, possibly because the climate model underestimates the millennial-scale temperature variability. In the course of the deglaciation, large-scale grounding line instabilities are simulated both for the Eurasian and North American ice sheets. The first instability occurs in the Barents–Kara seas for the Eurasian ice sheet at 14.5 ka. A second grounding line instability occurs ca. 12 ka in the proglacial lake that formed at the southern margin of the North American ice sheet. With additional asynchronously coupled experiments, we assess the sensitivity of our results to different ice sheet model choices related to surface and sub-shelf mass balance, ice deformation and grounding line representation. While the ice sheet evolutions differ within this ensemble, the global climate trajectory is only weakly affected by these choices. In our experiments, only the abrupt shifts in the oceanic circulation due to freshwater fluxes are able to produce some millennial-scale variability since no self-generating abrupt transitions are simulated without these fluxes.

Description: Rapid and continuous analysis of radiocarbon (14C) concentration in carbonate samples at spatial resolution down to 100 µm has been made possible with the new LA-AMS (laser ablation accelerator mass spectrometry) technique. This novel approach can provide radiocarbon data at a spatial resolution similar to that of stable carbon (C) isotope measurements by isotope ratio mass spectrometry of micromilled samples and, thus, can help to interpret δ13C signatures, which otherwise are difficult to understand due to numerous processes contributing to changes in the C-isotope ratio. In this work, we analyzed δ13C and 14C on the Holocene stalagmite SPA 127 from the high-alpine Spannagel Cave (Austria). Both proxies respond in a complex manner to climate variability. Combined stable carbon and radiocarbon profiles allow three growth periods characterized by different δ13C signatures to be identified: (i) the period 8.5 to 8.0 ka is characterized by relatively low δ13C values with small variability combined with a comparably high radiocarbon reservoir effect (expressed as dead carbon fraction, dcf) of around 60 %. This points towards C contributions of host rock dissolution and/or from an “old” organic matter (OM) reservoir in the karst potentially mobilized due to the warm climatic conditions of the early Holocene. (ii) Between 8 and 3.8 ka there was a strong variability in δ13C with values ranging from −8 ‰ to +1 ‰ and a generally lower dcf. The δ13C variability is most likely caused by changes in C exchange between cave air CO2 and dissolved inorganic carbon in drip water in the cave, which are induced by reduced drip rates as derived from reduced stalagmite growth rates. Additionally, the lower dcf indicates that the OM reservoir contributed less to stalagmite growth in this period possibly as a result of reduced meteoric precipitation or because it was exhausted. (iii) In the youngest section between 3.8 and 2.4 ka, comparably stable and low δ13C values, combined with an increasing dcf reaching up to 50 % again, hint towards a contribution of an aged OM reservoir in the karst. This study reveals the potential of combining high-resolution 14C profiles in speleothems with δ13C records in order to disentangle climate-related C dynamics in karst systems.

Description: The ∼74 ka Toba eruption was one of the largest volcanic events of the Quaternary. There is much interest in determining the impact of such a large event, particularly on the climate and hominid populations at the time. Although the Toba eruption has been identified in both land and marine archives as the Youngest Toba Tuff, its precise place in the ice core record is ambiguous. Several volcanic sulfate signals have been identified in both Antarctic and Greenland ice cores and span the Toba eruption 40Ar/39Ar age uncertainty. Here, we measure sulfur isotope compositions in Antarctic ice samples from the Dome C (EDC) and Dronning Maud Land (EDML) ice cores at high temporal resolution across 11 of these potential Toba sulfate peaks to identify candidates with sulfur mass-independent fractionation (S-MIF), indicative of an eruption whose plume reached altitudes at or above the stratospheric ozone layer. Using this method, we identify several candidate sulfate peaks that contain stratospheric sulfur. We further narrow down potential candidates based on the isotope signatures by identifying sulfate peaks that are due to a volcanic event at tropical latitudes. In one of these sulfate peaks at 73.67 ka, we find the largest ever reported magnitude of S-MIF in volcanic sulfate in polar ice, with a Δ33S value of −4.75 ‰. As there is a positive correlation between the magnitude of the S-MIF signal recorded in ice cores and eruptive plume height, this could be a likely candidate for the Toba super-eruption, with a plume top height in excess of 45 km. These results support the 73.7±0.3 ka (1σ) 40Ar/39Ar age estimate for the eruption, with ice core ages of our candidates with the largest magnitude S-MIF at 73.67 and 73.74 ka. Finally, since these candidate eruptions occurred on the transition into Greenland Stadial 20, the relative timing suggests that Toba was not the trigger for the large Northern Hemisphere cooling at this time although we cannot rule out an amplifying effect.

Description: Here we present the experimental design and results from a new mid-Pliocene simulation using the latest version of the UK's physical climate model, HadGEM3-GC31-LL, conducted under the auspices of CMIP6/PMIP4/PlioMIP2. Although two other palaeoclimate simulations have been recently run using this model, they both focused on more recent periods within the Quaternary, and therefore this is the first time this version of the UK model has been run this far back in time. The mid-Pliocene Warm Period, ∼3 Ma, is of particular interest because it represents a time period when the Earth was in equilibrium with CO2 concentrations roughly equivalent to those of today, providing a possible analogue for current and future climate change. The implementation of the Pliocene boundary conditions is firstly described in detail, based on the PRISM4 dataset, including CO2, ozone, orography, ice mask, lakes, vegetation fractions and vegetation functional types. These were incrementally added into the model, to change from a pre-industrial setup to a Pliocene setup. The results of the simulation are then presented, which are firstly compared with the model's pre-industrial simulation, secondly with previous versions of the same model and with available proxy data, and thirdly with all other models included in PlioMIP2. Firstly, the comparison with the pre-industrial simulation suggests that the Pliocene simulation is consistent with current understanding and existing work, showing warmer and wetter conditions, and with the greatest warming occurring over high-latitude and polar regions. The global mean surface air temperature anomaly at the end of the Pliocene simulation is 5.1 ∘C, which is the second highest of all models included in PlioMIP2 and is consistent with the fact that HadGEM3-GC31-LL has one of the highest Effective Climate Sensitivities of all CMIP6 models. Secondly, the comparison with previous generation models and with proxy data suggests a clear increase in global sea surface temperatures as the model has undergone development. Up to a certain level of warming, this results in a better agreement with available proxy data, and the “sweet spot” appears to be the previous CMIP5 generation of the model, HadGEM2-AO. The most recent simulation presented here, however, appears to show poorer agreement with the proxy data compared with HadGEM2 and may be overly sensitive to the Pliocene boundary conditions, resulting in a climate that is too warm. Thirdly, the comparison with other models from PlioMIP2 further supports this conclusion, with HadGEM3-GC31-LL being one of the warmest and wettest models in all of PlioMIP2, and if all the models are ordered according to agreement with proxy data, HadGEM3-GC31-LL ranks approximately halfway among them. A caveat to these results is the relatively short run length of the simulation, meaning the model is not in full equilibrium. Given the computational cost of the model it was not possible to run it for a longer period; a Gregory plot analysis indicates that had it been allowed to come to full equilibrium, the final global mean surface temperature could have been approximately 1.5 ∘C higher.

Description: The evolution of the Cenozoic cryosphere from unipolar to bipolar over the past 30 million years (Myr) is broadly known. Highly resolved records of carbonate (CaCO3) content provide insight into the evolution of regional and global climate, cryosphere, and carbon cycle dynamics. Here, we generate the first Southeast Atlantic CaCO3 content record spanning the last 30 Myr, derived from X-ray fluorescence (XRF) ln(Ca / Fe) data collected at Ocean Drilling Program Site 1264 (Walvis Ridge, SE Atlantic Ocean). We present a comprehensive and continuous depth and age model for the entirety of Site 1264 (∼ 316 m; 30 Myr). This constitutes a key reference framework for future palaeoclimatic and palaeoceanographic studies at this location. We identify three phases with distinctly different orbital controls on Southeast Atlantic CaCO3 deposition, corresponding to major developments in climate, the cryosphere and the carbon cycle: (1) strong ∼ 110 kyr eccentricity pacing prevails during Oligocene–Miocene global warmth (∼ 30–13 Ma), (2) increased eccentricity-modulated precession pacing appears after the middle Miocene Climate Transition (mMCT) (∼ 14–8 Ma), and (3) pervasive obliquity pacing appears in the late Miocene (∼ 7.7–3.3 Ma) following greater importance of high-latitude processes, such as increased glacial activity and high-latitude cooling. The lowest CaCO3 content (92 %–94 %) occurs between 18.5 and 14.5 Ma, potentially reflecting dissolution caused by widespread early Miocene warmth and preceding Antarctic deglaciation across the Miocene Climatic Optimum (∼ 17–14.5 Ma) by 1.5 Myr. The emergence of precession pacing of CaCO3 deposition at Site 1264 after ∼ 14 Ma could signal a reorganisation of surface and/or deep-water circulation in this region following Antarctic reglaciation at the mMCT. The increased sensitivity to precession at Site 1264 between 14 and 13 Ma is associated with an increase in mass accumulation rates (MARs) and reflects increased regional CaCO3 productivity and/or recurrent influxes of cooler, less corrosive deep waters. The highest carbonate content (%CaCO3) and MARs indicate that the late Miocene–early Pliocene Biogenic Bloom (LMBB) occurs between ∼ 7.8 and 3.3 Ma at Site 1264; broadly contemporaneous with the LMBB in the equatorial Pacific Ocean. At Site 1264, the onset of the LMBB roughly coincides with appearance of strong obliquity pacing of %CaCO3, reflecting increased high-latitude forcing. The global expression of the LMBB may reflect increased nutrient input into the global ocean resulting from enhanced aeolian dust and/or glacial/chemical weathering fluxes, due to enhanced glacial activity and increased meridional temperature gradients. Regional variability in the timing and amplitude of the LMBB may be driven by regional differences in cooling, continental aridification and/or changes in ocean circulation in the late Miocene.

Description: In this study, we analyze extreme daily precipitation during the pre-industrial period from 1501 BCE to 1849 CE in simulations from the Community Earth System Model version 1.2.2. A peak-over-threshold (POT) extreme value analysis is employed to examine characteristics of extreme precipitation and to identify connections of extreme precipitation with the external forcing and with modes of internal variability. The POT analysis shows that extreme precipitation with similar statistical characteristics, i.e., the probability density distributions, tends to cluster spatially. There are differences in the distribution of extreme precipitation between the Pacific and Atlantic sectors and between the northern high and southern low latitudes. Extreme precipitation during the pre-industrial period is largely influenced by modes of internal variability, such as El Niño–Southern Oscillation (ENSO), the Pacific North American, and Pacific South American patterns, among others, and regional surface temperatures. In general, the modes of variability exhibit a statistically significant connection to extreme precipitation in the vicinity to their regions of action. The exception is ENSO, which shows more widespread influence on extreme precipitation across the Earth. In addition, the regions with which extreme precipitation is more associated, either by a mode of variability or by the regional surface temperature, are distinguished. Regional surface temperatures are associated with extreme precipitation over lands at the extratropical latitudes and over the tropical oceans. In other regions, the influence of modes of variability is still dominant. Effects of the changes in the orbital parameters on extreme precipitation are rather weak compared to those of the modes of internal variability and of the regional surface temperatures. Still, some regions in central Africa, southern Asia, and the tropical Atlantic ocean show statistically significant connections between extreme precipitation and orbital forcing, implying that in these regions, extreme precipitation has increased linearly during the 3351-year pre-industrial period. Tropical volcanic eruptions affect extreme precipitation more clearly in the short term up to a few years, altering both the intensity and frequency of extreme precipitation. However, more apparent changes are found in the frequency than the intensity of extreme precipitation. After eruptions, the return periods of extreme precipitation increase over the extratropical regions and the tropical Pacific, while a decrease is found in other regions. The post-eruption changes in the frequency of extreme precipitation are associated with ENSO, which itself is influenced by tropical eruptions. Overall, the results show that climate simulations are useful to complement the information on pre-industrial extreme precipitation, as they elucidate statistical characteristics and long-term connections of extreme events with natural variability.