Abstract. Improving the understanding of changes in the mean and variability of climate variables as well as their interrelation is crucial for reliable climate change projections. Comparisons between general circulation models and paleoclimate archives using indirect proxies for temperature or precipitation have been used to test and validate the capability of climate models to represent climate changes. The oxygen isotopic ratio δ18O, a proxy for many different climate variables, is routinely measured in speleothem samples at decadal or higher resolution, and single specimens can cover full glacial–interglacial cycles. The calcium carbonate cave deposits are precisely dateable and provide well preserved (semi-)continuous albeit multivariate climate signals in the lower and mid-latitudes, where the measured δ18O in the mineral does not directly represent temperature or precipitation. Therefore, speleothems represent suitable archives to assess climate model abilities to simulate climate variability beyond the timescales covered by meteorological observations (101–102 years). Here, we present three transient isotope-enabled simulations from the Hadley Center Climate Model version 3 (iHadCM3) covering the last millennium (850–1850 CE) and compare them to a large global dataset of speleothem δ18O records from the Speleothem Isotopes Synthesis and AnaLysis (SISAL) database version 2 (Comas-Bru et al., 2020b). We systematically evaluate offsets in mean and variance of simulated δ18O and test for the main climate drivers recorded in δ18O for individual records or regions. The time-mean spatial offsets between the simulated δ18O and the speleothem data are fairly small. However, using robust filters and spectral analysis, we show that the observed archive-based variability of δ18O is lower than simulated by iHadCM3 on decadal and higher on centennial timescales. Most of this difference can likely be attributed to the records' lower temporal resolution and averaging or smoothing processes affecting the δ18O signal, e.g., through soil water residence times. Using cross-correlation analyses at site level and modeled grid-box level, we find evidence for highly variable but generally low signal-to-noise ratios in the proxy data. This points to a high influence of cave-internal processes and regional climate particularities and could suggest low regional representativity of individual sites. Long-range strong positive correlations dominate the speleothem correlation network but are much weaker in the simulation. One reason for this could lie in a lack of long-term internal climate variability in these model simulations, which could be tested by repeating similar comparisons with other isotope-enabled climate models and paleoclimate databases.

Abstract. Global changes in the climate, especially the warming trend in mean temperature, have received increasing public and scientific attention. Improving the understanding of changes in the mean and variability of climate variables as well as their interrelation is crucial for reliable climate change projections. Comparisons between general circulation models and paleoclimate archives using indirect proxies for temperature and/or precipitation have been used to test and validate the capability of climate models to represent climate changes. The oxygen isotopic ratio δ18O is routinely measured in speleothem samples at decadal or higher resolution and single specimens can cover full Glacial-Interglacial cycles. The calcium carbonate cave deposits are precisely dateable and provide well preserved (semi-) continuous, albeit multivariate climate signals in the lower and mid-latitudes, where the measured δ18O in the mineral does not directly represent temperature or precipitation. Therefore, speleothems represent suitable archives to assess simulated climate model abilities for the simulation of climate variability beyond the timescales covered by meteorological observations (10–100 yr). Here, we present three transient isotope enabled simulations from the Hadley Center Climate Model version 3 (iHadCM3) covering the last millennium (850–1850 CE) and compare these to a large global dataset of speleothem δ18O records from the Speleothem Isotopes Synthesis and AnaLysis (SISAL) database version 2 (Comas-Bru et al., 2020). We evaluate systematically offsets in mean and variance of simulated δ18O and test for the main climate drivers for individual records or regions. The time-mean spatial offsets between the simulated δ18O and the speleothem data are fairly small. However, using robust filters and spectral analysis, we show that the observed proxy-based variability of δ18O is lower (higher) than simulated by iHadCM3 on decadal (centennial) timescales. Most of this difference can likely be attributed to the records' lower temporal resolution and averaging processes affecting the δ18O signal. Using cross-correlation analyses at site-level and modeled gridbox level, we find evidence for highly variable but generally low signal-to-noise ratios in the proxy data. This points at a high influence of cave-internal processes and regional climate particularities and could suggest low regional representativity of individual sites. Long-range strong positive correlations dominate the speleothem correlation network but are much weaker in the simulation. One reason for this could lie in a lack of longterm internal climate variability in these model simulations, which could be tested by repeating similar comparisons with other isotope-enabled climate models and paleoclimate databases.

Abstract. The incorporation of water isotopologues into the hydrology of general circulation models (GCMs) facilitates the comparison between modeled and measured proxy data in paleoclimate archives. However, the variability and drivers of measured and modeled water isotopologues, as well as the diversity of their representation in different models, are not well constrained. Improving our understanding of this variability in past and present climates will help to better constrain future climate change projections and decrease their range of uncertainty. Speleothems are a precisely datable terrestrial paleoclimate archives and provide well-preserved (semi-)continuous multivariate isotope time series in the lower latitudes and mid-latitudes and are therefore well suited to assess climate and isotope variability on decadal and longer timescales. However, the relationships of speleothem oxygen and carbon isotopes to climate variables are influenced by site-specific parameters, and their comparison to GCMs is not always straightforward. Here we compare speleothem oxygen and carbon isotopic signatures from the Speleothem Isotopes Synthesis and Analysis database version 2 (SISALv2) to the output of five different water-isotope-enabled GCMs (ECHAM5-wiso, GISS-E2-R, iCESM, iHadCM3, and isoGSM) over the last millennium (850–1850 CE). We systematically evaluate differences and commonalities between the standardized model simulation outputs. The goal is to distinguish climatic drivers of variability for modeled isotopes and compare them to those of measured isotopes. We find strong regional differences in the oxygen isotope signatures between models that can partly be attributed to differences in modeled surface temperature. At low latitudes, precipitation amount is the dominant driver for stable water isotope variability; however, at cave locations the agreement between modeled temperature variability is higher than for precipitation variability. While modeled isotopic signatures at cave locations exhibited extreme events coinciding with changes in volcanic and solar forcing, such fingerprints are not apparent in the speleothem isotopes. This may be attributed to the lower temporal resolution of speleothem records compared to the events that are to be detected. Using spectral analysis, we can show that all models underestimate decadal and longer variability compared to speleothems (albeit to varying extents). We found that no model excels in all analyzed comparisons, although some perform better than the others in either mean or variability. Therefore, we advise a multi-model approach whenever comparing proxy data to modeled data. Considering karst and cave internal processes, e.g., through isotope-enabled karst models, may alter the variability in speleothem isotopes and play an important role in determining the most appropriate model. By exploring new ways of analyzing the relationship between the oxygen and carbon isotopes, their variability, and co-variability across timescales, we provide methods that may serve as a baseline for future studies with different models using, e.g., different isotopes, different climate archives, or different time periods.

<p>Speleothems are terrestrial paleoclimate archives that occur abundantly in the low and mid latitudes. They archive changes in the past hydroclimate in many ways, including the rate of calcium carbonate accumulation – their growth rate. However, determining speleothem growth rates, and in particular growth rate changes, is challenging due to speleothem inherent features such as growth hiatuses, and large and abrupt growth rate changes. Low dating resolution poses an additional problem, as the U/Th measurements that allow for precise dating are time-consuming and expensive.</p> <p>Here, we analyze speleothem growth rates during the Holocene – an ideal period for method testing due to the high abundance of speleothem records in the SISALv2 database. In particular, we compare speleothem growth rates in the early (12-8 kyr BP), mid (8-4 kyr BP) and the late Holocene (4-0 kyr BP). Using synthetically-modelled stalagmites, we test the strengths and weaknesses of state-of-the-art age-depth modelling methods to determine a set of necessary requirements to quantify speleothem growth rates and growth rate changes. Using these, we find slightly higher growth rates in the early Holocene within speleothems that cover at all periods. Comparing growth rates of speleothems that cover only one of the respective periods in the Holocene did not distinguish any period of highest or lowest growth rate. Detailed regional studies and comparison to model data are used to further interpret these results. Reliably determining growth rate changes in the Holocene may help in further understanding and characterizing hydroclimate changes as archived in speleothems also beyond the Holocene.</p>

<div> <p>Data Assimilation in paleoclimatology (PaleoDA) is a method that has been used in several climate reconstructions for the last millennium. By fusing information from both climate proxies and general circulation models (GCMs), PaleoDA provides statistical estimates of climate fields that are dynamically consistent. However, existing reconstructions mostly rely on calibrated tree ring data and assimilate proxy records on a single, annual time scale. Ice cores and speleothems, which record past variations in the oxygen isotope ratio of precipitation, often have a lower and irregular time resolution, but reliably record climate variations on decadal to centennial time scales.<span data-ccp-props="{}"> </span></p> </div> <div> <p>Here, we implemented a computationally efficient DA algorithm that enables the assimilation of proxy records on multiple timescales. The algorithm has been applied to speleothem and ice core records from the SISALv2 and Iso2k database and five isotope-enabled GCMs. Reconstructions of global mean temperature changes during the last millennium compare well in both amplitude and uncertainty to recent studies. The potential of incorporating speleothems is shown with a reconstruction of hydroclimatic changes in tropical South America, where speleothems represent the most abundant type of hydroclimate archive. The experiments performed suggest an increased reconstructed decadal to centennial variability by using proxy records on multiple timescales. Making use of different climate models shows the influence of model biases on the reconstructions. Future PaleoDA reconstructions could be improved from more proxy records and the multiple time scale approach to provide a globally complete picture of past climate changes.<span data-ccp-props="{}"> </span></p> </div>

<p align="justify">The Pantanal is a large region located in the central parts of South America (140,000 km2) with a unique climate and vegetation setting. This region is subjected to seasonal floods, which makes the Pantanal one of the most important wetlands on the planet. In this region occur transitions between different biomes, such as the Amazon Forest, Cerrado (Brazilian Savanna), and Atlantic Forest, located to the North, East, and South of this region, respectively. The area also serves as a moisture pathway for the South American Summer Monsoon (SASM), which connects the Amazon Basin with the La Plata Basin. The two major drainage basins of South America. Due to the complex hydrology of the rivers and lakes of this region, it is necessary to combine multiple proxy archives from different parts of the Pantanal basin to understand its climate and vegetation evolution during the Holocene.</p> <p>Here we present isotope records from stalagmites collected at sites located at the northern and southern borders of the Pantanal. Hiatuses in speleothem deposition during the mid-Holocene identified in several stalagmites indicate overall dry conditions in the region at this period. However, the drier conditions recorded in the northern portion of the basin precede similar conditions in the South by approximately two thousand years. Furthermore, summer insolation seems to drive the intensity of the SASM at the North Pantanal, while its influence is weaker in the southern part. During the late Holocene, this establishes a moisture gradient between a wetter North and drier South. Our record also shows a strong multidecadal to centennial variability, which was probably responsible for the high hydrology complexity of the rivers of the Pantanal, which are subject to seasonal floods and migration of its channels and tributaries.</p>

Abstract. The incorporation of water isotopologues into the hydrology of general circulation models (GCMs) facilitates the comparison between modelled and measured proxy data in paleoclimate archives. However, the variability and drivers of measured and modelled water isotopologues, and indeed the diversity of their representation in different models are not well constrained. Improving our understanding of this variability in past and present climates will help to better constrain future climate change projections and decrease their range of uncertainty. Speleothems are a precisely datable paleoclimate archive and provide well preserved (semi-)continuous multivariate isotope time series in the lower and mid-latitudes, and are, therefore, well suited to assess climate and isotope variability on decadal and longer timescales. However, the relationship between speleothem oxygen and carbon isotopes to climate variables also depends on site-specific parameters, and their comparison to GCMs is not always straightforward. Here we compare speleothem oxygen and carbon isotopic signatures from the Speleothem Isotopes Synthesis and AnaLysis database version 2 (SISALv2) to the output of five different water-isotope-enabled GCMs (ECHAM5-wiso, GISS-E2-R, iCESM, iHadCM3, and isoGSM) over the last millennium (850–1850 common era, CE). We systematically evaluate differences and commonalities between the standardized model simulation outputs. The goal is to distinguish climatic drivers of variability for both modelled and measured isotopes. We find strong regional differences in the oxygen isotope signatures between models that can partly be attributed to differences in modelled temperatures. At low latitudes, precipitation amount is the dominant driver for water isotope variability, however, at cave locations the agreement between modelled temperature variability is higher than for precipitation variability. While modelled isotopic signatures at cave locations exhibited extreme events coinciding with changes in volcanic and solar forcing, such fingerprints are not apparent in the speleothem isotopes, and may be attributed to the lower temporal resolution of speleothem records compared to the events that are to be detected. Using spectral analysis, we can show that all models underestimate decadal and longer variability compared to speleothems, although to varying extent. We found that no model excels in all analyzed comparisons, although some perform better than the others in either mean or variability. Therefore, we advise a multi-model approach, whenever comparing proxy data to modelled data. Considering karst and cave internal processes through e.g. isotope-enabled karst models may alter the variability in speleothem isotopes and play an important role in determining the most appropriate model. By exploring new ways of analyzing the relationship between the oxygen and carbon isotopes, their variability, and co-variability across timescales, we provide methods that may serve as a baseline for future studies with different models using e.g. different isotopes, different climate archives, or time periods.

No abstract