A window into Africa’s past hydroclimates: The sisal_v1 database contribution

Africa spans the hemispheres from temperate region to temperate region and has a long history of hominin evolution. Although the number of Quaternary palaeoclimatic records from the continent is increasing, much of the history of spatial and temporal climatic variability is still debated. Speleothems, as archives of terrestrial hydroclimate variability, can help reveal this history. Here we review the progress made to date, with a focus on the first version of the Speleothem Isotopes Synthesis and AnaLysis (SISAL) database. The geology of Africa has limited development of large karst regions to four areas: along the northern coast bordering the Mediterranean, eastern Africa and the Horn of Africa, southwestern Africa and southern Africa. Exploitation of the speleothem palaeoclimate archives in these regions is uneven, with long histories of research, e.g., in South Africa, but large areas with no investigations such as West Africa. Consequently, the evidence of past climate change reviewed here is irregularly sampled in both time and space. Nevertheless, we show evidence of migration of the monsoon belt, with enhanced rainfall during interglacials observed in northeast Africa, southern Arabia and the northern part of southern Africa. Evidence from eastern Africa indicates significant decadal and centennial scale rainfall variability. In northwestern and southern Africa, precession and eccentricity influence speleothem growth, largely through changing synoptic storm activity

Tropical Indian Ocean basin hydroclimate at the Mid- to Late-Holocene transition and the double drying hypothesis

The spatial pattern of Holocene climate anomalies is crucial to determining the mechanisms of change, distinguishing between unforced and forced climate variability, and understanding potential impacts on past and future human societies. The 4.2 ka event is often regarded as one of the largest and best documented abrupt climate disturbances of the Holocene. Yet outside the data-rich Northern Hemisphere mid-latitudes, the global pattern of climate anomalies is uncertain. In this study we investigate the spatial and temporal variability of the tropical Indian Ocean hydroclimate at the Mid- to Late-Holocene transition. We conducted Monte-Carlo principal component analysis, considering full age uncertainty, on ten high-resolution, precisely dated paleohydroclimate records from around the tropical Indian Ocean basin, all growing continuously or almost continuously between 5 and 3 kyr BP. The results indicate the dominant mode of variability in the region was a drying between 3.97 kyr BP (±0.08 kyr standard error) and 3.76 kyr BP (±0.07 kyr standard error) with dry conditions lasting for an additional 300 years in some records, and a permanent change in others. This drying in PC1, which we interpret as a proxy of summer monsoon variability, fits with a previously recognised tropic wide change in hydroclimate around 4.0 kyr BP. An abrupt event from 4.2 to 3.9 kyr BP is seen locally in individual records but lacks regional coherence. A lack of apparent 4.2 ka event in tropical Indian Ocean hydroclimate has ramifications for climate variability in the Indus valley, and for the Harappan civilization. Through a comparison of existing Indian subcontinent paleoclimate records, upstream climatic variability in the Indian Summer Monsoon and winter Westerly Disturbances source regions, and modern climatology, we present the “Double Drying hypothesis”. A winter rainfall drying between 4.2 and 3.9 kyr BP was followed by a summer rainfall drying between 3.97 kyr BP and at least 3.4 kyr BP. The Double Drying hypothesis provides more detailed climatic context for the Harappan civilization, resolves the cropping paradox, and fits the spatial-temporal pattern of urban abandonment. The consequences for the new Mid- to Late-Holocene Global Boundary Stratotype Section and Point in a stalagmite from Meghalaya are explored

Evaluating model outputs using integrated global speleothem records of climate change since the last glacial

Although quantitative isotope data from speleothems has been used to evaluate isotope-enabled model simulations, currently no consensus exists regarding the most appropriate methodology through which to achieve this. A number of modelling groups will be running isotope-enabled palaeoclimate simulations in the framework of the Coupled Model Intercomparison Project Phase 6, so it is timely to evaluate different approaches to using the speleothem data for data–model comparisons. Here, we illustrate this using 456 globally distributed speleothem δ18O records from an updated version of the Speleothem Isotopes Synthesis and Analysis (SISAL) database and palaeoclimate simulations generated using the ECHAM5-wiso isotope-enabled atmospheric circulation model. We show that the SISAL records reproduce the first-order spatial patterns of isotopic variability in the modern day, strongly supporting the application of this dataset for evaluating model-derived isotope variability into the past. However, the discontinuous nature of many speleothem records complicates the process of procuring large numbers of records if data–model comparisons are made using the traditional approach of comparing anomalies between a control period and a given palaeoclimate experiment. To circumvent this issue, we illustrate techniques through which the absolute isotope values during any time period could be used for model evaluation. Specifically, we show that speleothem isotope records allow an assessment of a model's ability to simulate spatial isotopic trends. Our analyses provide a protocol for using speleothem isotope data for model evaluation, including screening the observations to take into account the impact of speleothem mineralogy on δ18O values, the optimum period for the modern observational baseline and the selection of an appropriate time window for creating means of the isotope data for palaeo-time-slices.Financial support for SISAL activities that have lead to this research has been provided by the Past Global Changes (PAGES) programme; the European Geosciences Union (grant no. W2017/413); the Irish Centre for Research in Applied Geosciences (iCRAG); the European Association of Geochemistry (Early Career Ambassadors program 2017); the Quaternary Research Association UK; the Navarino Environmental Observatory, Stockholm University; University College Dublin (grant no. SF1428), Savillex (UK); John Cantle; Ibn Zohr University, Morocco; the University of Reading; the European Research Council (grant no. 694481); the Natural Environment Research Council (JPI-Belmont project “PAleao-Constraints on Monsoon Evolution and Dynamics (PACMEDY)”); the Geological Survey Ireland (grant no. 2017-SC-056); the Royal Irish Academy (Charlemont Scholar award 2018); the Portuguese Science Foundation (grant no. UID/MAR/00350/2013); and the Deutsche Forschungsgemeinschaft (grant no. RE3994/2-1)

Southern Hemisphere controls on ITCZ variability in southwest Madagascar over the past 117,000 years

Migration of the inter-tropical convergence zone, driven by changes in seasonal insolation and high northern latitude temperatures, is the primary control on tropical rainfall on geologic timescales. We test this paradigm using the timing of growth of stalagmites from southwest Madagascar to infer the timing of expansion of the ITCZ to the south at its southern limit. Over the past 117 ky, speleothems grew in the study area primarily when two conditions are met: summer insolation greater than the mean and relatively high Southern Hemisphere temperatures as indicted by maxima in Antarctic ice core oxygen isotope ratios. We observe little influence of Northern Hemisphere, millennial scale temperature variability on the pluvial periods. Further, we observe periods during which the ITCZ simultaneously expands or contracts in both hemispheres. Because Antarctic isotope maxima are periods of increased atmospheric CO2, our results have implications for how tropical rainfall in the Southern Hemisphere might respond to global warming

Hydroclimate variability in the Madagascar and Southeast African summer monsoons at the Mid- to Late-Holocene transition

The 4.2 ka event at the Mid- to Late-Holocene transition is often regarded as one of the largest and best documented abrupt climate disturbances of the Holocene. The event is most clearly manifested in the Mediterranean and Middle East as a regional dry anomaly beginning abruptly at 4.26 kyr BP and extending until 3.97 kyr BP. Yet the impacts of this regional drought are often extended to other regions and sometimes globally. In particular, the nature and spatial extent of the 4.2 ka event in the tropics have not been established. Here, we present a new stalagmite stable isotope record from Anjohikely, northwest Madagascar. Growing between 5.22 and 2.00 kyr BP, stalagmite AK1 shows a hiatus between 4.31 and 3.93 kyr BP (±40 and ± 35 yrs), replicating a hiatus in another stalagmite from nearby Anjohibe, and therefore indicating a significant drying at the Mid- to Late-Holocene transition. This result is the opposite to wet conditions at the 8.2 ka event, suggesting fundamentally different forcing mechanisms. Dry conditions are also recorded in sediment cores in Lake Malawi, Lake Masoko and the Tatos Basin on Mauritius, also in the southeast African monsoon domain. However, no notable event is recorded at the northern (equatorial East Africa) and eastern (Rodrigues) peripheries of the monsoon domain, while a wet event is recorded in sediment cores at Lake Muzi and Mkhuze Delta to the south. The spatial pattern is largely consistent with the modern rainfall anomaly pattern associated a with weak Mozambique Channel Trough and a northerly austral summer Intertropical Convergence Zone position. Within age error, the observed peak climate anomalies overlap with the 4.2 ka event. However regional hydrological change consistently begins earlier than a 4.26 kyr BP event onset. Gradual hydrological change frequently begins around 4.5 kyr BP, raising doubt as to whether any coherent regional hydrological change is merely coincident with the 4.2 ka event or part of a global climatic anomaly

Comparing the paleoclimates of northwestern and southwestern Madagascar during the late Holocene: Implications for the role of climate in megafaunal extinction.

The relative importance of climate and humans in the disappearance of the Malagasy megafauna remains under debate. Data from southwestern Madagascar imply aridifcation contributed substantially to the late Holocene decline of the megafauna (the Aridifcation Hypothesis). Evidence for aridifcation includes carbon isotopes from tree rings, lacustrine charcoal concentrations and pollen assemblages, and changes in fossil vertebrate assemblages indicative of a local loss of pluvial conditions. In contrast, speleothem records from northwestern Madagascar suggest that megafaunal decline and habitat change resulted primarily from human activity including agropastoralism (the Subsistence Shift Hypothesis). Could there have been contrasting mechanisms of decline in different parts of Madagascar? Or are we lacking the precisely dated, high resolution records needed to fully understand the complex processes behind megafaunal decline? Reconciling these contrasting hypotheses requires additional climate records from southwestern Madagascar. We recovered a stalagmite (AF2) from Asafora Cave in the spiny thicket ecoregion, ~10 km from the southwest coast and just southeast of the Velondriake Marine Reserve. U-series and 14C dating of samples taken from the core of this stalagmite provide a highly precise chronology of the changes in hydroclimate and vegetation in this region over the past 3000 years. Speleothem stable oxygen and carbon isotope analyses provide insight into past rainfall variability and vegetation changes respectively. We compare these records with those for a stalagmite (AB2) from Anjohibe Cave in northwestern Madagascar. Lastly, odds ratio analyses of radiocarbon dates for extinct and extant subfossils allow us to describe and compare the temporal trajectories of megafaunal decline in the southwest and the northwest. Combined, these analyses allow us to test the Aridifcation Hypothesis for megafaunal extinction. The trajectories of megafaunal decline differed in northwestern and southwestern Madagascar. In the southwest, unlike the northwest, there is no evidence of decoupling of speleothem stable carbon and oxygen isotopes. Instead, habitat changes in the southwest were largely related to variation in hydroclimate (including a prolonged drought). The megafaunal collapse here occurred in tandem with the drought, and agropastoralism likely contributed to that demise only after the megafauna had already suffered drought-related population reduction. Our results offer some support for the Aridifcation Hypothesis, but with three caveats: frst, that there was no island-wide aridifcation; second, that aridifcation likely impacted megafaunal decline only in the driest parts of Madagascar; and third, that aridifcation was not the sole factor promoting Comparing the paleoclimates of northwestern and southwestern Madagascar during the late Holocene: Implications for the role of climate in megafaunal extinction Faina et al.: Comparing the paleoclimates of northwestern and southwestern Madagascar 109 megafaunal decline even in the dry southwest. A number of megafaunal species survived the prolonged drought of the first millennium, and then likely succumbed to the activities of agropastoralists

SISALv2: A comprehensive speleothem isotope database with multiple age-depth models

Characterizing the temporal uncertainty in palaeoclimate records is crucial for analysing past climate change, correlating climate events between records, assessing climate periodicities, identifying potential triggers and evaluating climate model simulations. The first global compilation of speleothem isotope records by the SISAL (Speleothem Isotope Synthesis and Analysis) working group showed that age model uncertainties are not systematically reported in the published literature, and these are only available for a limited number of records (ca. 15 %, n = 107=691). To improve the usefulness of the SISAL database, we have (i) improved the database's spatiooral coverage and (ii) created new chronologies using seven different approaches for age depth modelling. We have applied these alternative chronologies to the records from the first version of the SISAL database (SISALv1) and to new records compiled since the release of SISALv1. This paper documents the necessary changes in the structure of the SISAL database to accommodate the inclusion of the new age models and their uncertainties as well as the expansion of the database to include new records and the qualitycontrol measures applied. This paper also documents the age depth model approaches used to calculate the new chronologies. The updated version of the SISAL database (SISALv2) contains isotopic data from 691 speleothem records from 294 cave sites and new age depth models, including age depth temporal uncertainties for 512 speleothems. SISALv2 is available at https://doi.org/10.17864/1947.256 (Comas-Bru et al., 2020a). © 2020 Author(s)

Evaluating model outputs using integrated global speleothem records of climate change since the last glacial

Although quantitative isotopic data from speleothems has been used to evaluate isotope-enabled model simulations, currently no consensus exists regarding the most appropriate methodology through which to achieve this. A number of modelling groups will be running isotope-enabled palaeoclimate simulations in the framework of the Coupled Model Intercomparison Project Phase 6, so it is timely to evaluate different approaches to use the speleothem data for data-model comparisons. Here, we illustrate this using 456 globally-distributed speleothem δ18O records from an updated version of the Speleothem Isotopes Synthesis and Analysis (SISAL) database and palaeoclimate simulations generated using the ECHAM5-wiso isotope-enabled atmospheric circulation model. We show that the SISAL records reproduce the first-order spatial patterns of isotopic variability in the modern day, strongly supporting the application of this dataset for evaluating model-derived isotope variability into the past. However, the discontinuous nature of many speleothem records complicates procuring large numbers of records if data-model comparisons are made using the traditional approach of comparing anomalies between a control period and a given palaeoclimate experiment. To circumvent this issue, we illustrate techniques through which the absolute isotopic values during any time period could be used for model evaluation. Specifically, we show that speleothem isotope records allow an assessment of a model’s ability to simulate spatial isotopic trends. Our analyses provide a protocol for using speleothem isotopic data for model evaluation, including screening the observations to take into account the impact of speleothem mineralogy on 18O values, the optimum period for the modern observational baseline, and the selection of an appropriate time-window for creating means of the isotope data for palaeo time slices

The SISAL database: a global resource to document oxygen and carbon isotope records from speleothems

Stable isotope records from speleothems provide information on past climate changes, most particularly information that can be used to reconstruct past changes in precipitation and atmospheric circulation. These records are increasingly being used to provide “out-of-sample” evaluations of isotope-enabled climate models. SISAL (Speleothem Isotope Synthesis and Analysis) is an international working group of the Past Global Changes (PAGES) project. The working group aims to provide a comprehensive compilation of speleothem isotope records for climate reconstruction and model evaluation. The SISAL database contains data for individual speleothems, grouped by cave system. Stable isotopes of oxygen and carbon (δ 18O, δ 13C) measurements are referenced by distance from the top or bottom of the speleothem. Additional tables provide information on dating, including information on the dates used to construct the original age model and sufficient information to assess the quality of each data set and to erect a standardized chronology across different speleothems. The metadata table provides location information, information on the full range of measurements carried out on each speleothem and information on the cave system that is relevant to the interpretation of the records, as well as citations for both publications and archived data. The compiled data are available at https://doi.org/10.17864/1947.147