Rainfall recharge thresholds in a subtropical climate determined using a regional cave drip water monitoring network

Quantifying the combination of climatic and hydrological conditions required to generate groundwater recharge is challenging, yet of fundamental importance for groundwater resource management. Here we demonstrate a new unsaturated zone physical method of determining rainfall-recharge thresholds in karst using a regional cave drip water monitoring network. For limestones of the Upper and Lower Macleay Valley, eastern Australia, set in a subtropical climate, we observe thirty-one cave drip water recharge events over a five-year monitoring period. Comparison to antecedent precipitation demonstrates a median observed recharge threshold of 76 mm/week precipitation (Lower Macleay) and 79 mm/week precipitation (Upper Macleay), with lower precipitation thresholds (down to 30 mm/week) possible. We use a simple water budget model to quantify soil and epikarst water storage volumes and to test hypotheses of the hydrological controls. Modelled soil and epikarst water storage capacities of about 65 mm (Lower Macleay) and 80 mm (Upper Macleay) confirm a correspondence between observed weekly precipitation thresholds and soil and epikarst capacities. However, discrepancies between observed and simulated recharge events helps elucidate the likely recharge processes including focussed recharge bypassing the soil and epikarst store, overflow and drainage between multiple karst stores, and tree water use from depth. Our observed recharge thresholds and modelled soil and epikarst storage capacities are comparable to recharge thresholds estimated across a range of water-limited environments globally. The method is readily applicable to any karst region where drip loggers can be installed in a cave system in close proximity to surface climate data

Modern speleothem oxygen isotope hydroclimate records in water-limited SE Australia

Dryland regions are generally projected to become drier under future climate change scenarios. Understanding the long-term natural variability of dryland regions via paleo-reconstructions is therefore highly desirable. The δ18O of two coeval modern speleothems from Cathedral Cave, Wellington, in semi-arid SE Australia are compared to the instrumental record to assess its efficacy as a proxy of past hydrological variability. Stalagmite δ18O was modulated by the frequency of recharge events and epikarst evaporation of storage water. Prolonged intervals between recharge events, such as droughts, resulted in higher stalagmite δ18O. Conversely, periods with more frequent recharge events and a positive water balance, resulted in lower δ18O. Disequilibrium cave processes are likely to be enhanced during dry conditions, although it is argued that these will modulate δ18Ospel in the same direction as epikarst evaporation, effectively amplifying the response of δ18Ospel. Extreme events, such as floods and droughts, were also captured in the stalagmite records, although potentially with a lag of several years. We verify that modern speleothems from semi-arid regions can be used to reconstruct hydroclimate due to variations in δ18Ospel modulated by karst processes. Such records are archives of past changes in recharge rather than precipitation amount or surface temperature, as is commonly applied to speleothem records from non-water-limited regions

Hydroclimate of the Last Glacial Maximum and deglaciation in southern Australia's arid margin interpreted from speleothem records (23-15 ka)

Terrestrial data spanning the Last Glacial Maximum (LGM) and deglaciation from the southern Australian region are sparse and limited to discontinuous sedimentological and geomorphological records with relatively large chronological uncertainties. This dearth of records has hindered a critical assessment of the role of the Southern Hemisphere mid-latitude westerly winds on the region's climate during this time period. In this study, two precisely dated speleothem records for Mairs Cave, Flinders Ranges, are presented, providing for the first time a detailed terrestrial hydroclimatic record for the southern Australian drylands during 23-15 ka. Recharge to Mairs Cave is interpreted from the speleothem record by the activation of growth, physical flood layering, and δ18O and δ13C minima. Periods of lowered recharge are indicated by 18O and 13C enrichment, primarily affecting δ18O, argued to be driven by evaporation of shallow soil/epikarst water in this water-limited environment. A hydrological driver is supported by calcite fabric changes. These include the presence of laminae, visible organic colloids, and occasional dissolution features, related to recharge, as well as the presence of sediment bands representing cave floor flooding. A shift to slower-growing, more compact calcite and an absence of lamination is interpreted to represent reduced recharge. The Mairs Cave record indicates that the Flinders Ranges were relatively wet during the LGM and early deglaciation, particularly over the interval 18.9-15.8 ka. This wetter phase ended abruptly with a shift to drier conditions at 15.8 ka. These findings are in agreement with the geomorphic archives for this region, as well as the timing of events in records from the broader Australasian region. The recharge phases identified in the Mairs Cave record are correlated with, but antiphase to, the position of the westerly winds interpreted from marine core MD03-2611, located 550 km south of Mairs Cave in the Murray Canyons region. The implication is that the mid-latitude westerlies are located further south during the period of enhanced recharge in the Mairs Cave record (18.9-16 ka) and conversely are located further north when greater aridity is interpreted in the speleothem record. A further comparison with speleothem records from the northern Australasian region reveals that the availability of tropical moisture is the most likely explanation driving enhanced recharge, with further amplification of recharge occurring during the early half of Heinrich Stadial 1 (HS1), possibly influenced by a more southerly displaced Intertropical Convergence Zone (ITCZ). A rapid transition to aridity at 15.8 ka is consistent with a retraction of this tropical moisture source.The stable isotope measurements on MC-S1 and MC-S2 were supported by Australian Research Council LIEF grants LE0989624 and LE0668400, and U-Th age measurements by Australian Research Council LIEF grant LE098906

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

Semi-arid zone caves:Evaporation and hydrological controls on δ¹⁸O drip water composition and implications for speleothem paleoclimate reconstructions

Oxygen isotope ratios in speleothems may be affected by external processes that are independent of climate, such as karst hydrology and kinetic fractionation. Consequently, there has been a shift towards characterising and understanding these processes through cave monitoring studies, particularly focussing on temperate zones where precipitation exceeds evapotranspiration. Here, we investigate oxygen isotope systematics at Wellington Caves in semi-arid, SE Australia, where evapotranspiration exceeds precipitation. We use a novel D2O isotopic tracer in a series of artificial irrigations, supplemented by pre-irrigation data comprised four years of drip monitoring and three years of stable isotope analysis of both drip waters and rainfall. This study reveals that: (1) evaporative processes in the unsaturated zone dominate the isotopic composition of drip waters; (2) significant soil zone ‘wetting up’ is required to overcome soil moisture deficits in order to achieve infiltration, which is highly dependent on antecedent hydro-climatic conditions; (3) lateral flow, preferential flow and sorption in the soil zone are important in redistributing subsurface zone water; (4) isotopic breakthrough curves suggest clear evidence of piston-flow at some drip sites where an older front of water discharged prior to artificial irrigation water; and (5) water residence times in a shallow vadose zone (<2 m) are highly variable and can exceed six months. Oxygen isotope speleothem records from semi-arid regions are therefore more likely to contain archives of alternating paleo-aridity and paleo-recharge, rather than paleo-rainfall e.g. the amount effect or mean annual. Speleothem-forming drip waters will be dominated by evaporative enrichment, up to ∼3‰ in the context of this study, relative to precipitation-weighted mean annual rainfall. The oxygen isotope variability of such coeval records may further be influenced by flow path and storage in the unsaturated zone that is not only drip specific but also influenced by internal cave climatic conditions, which may vary spatially in the cave

Calibrating climate-d18O regression models for the interpretation of high-resolution speleothem d18O time series

Providing estimates of past climate changes on interannual-millenial timescales requires suitable regression models between climate and climate proxies. Many proxies appear to show relationships with climate that are timescale dependent. Any proxy-climate model should be able to replicate the major patterns that are observed at multiple timescales. Here we develop a new climate-isotope regression model for speleothems from a middle latitude site. In the low to middle latitudes, daily variation in precipitation isotopes (within individual months) is largely negatively correlated with daily rainfall amount. On interdecadal timescales, though, this relationship appears to be nonstationary. These two points provide a theoretical basis for a new climate-isotope regression model in which the slope and the intercept of a δ18 Oday-Pday line for a given month are modulated by organized patterns of climate variability, such as the extratropical zonal waves (including the annular modes). In constructing this new regression model, we show how daily precipitation-δ18O relationships can be estimated using only monthly δ18O data and daily rainfall amounts. The new regression model provides a consistent picture of 18O variability over a range of timescales, and this has not been the case with any previous climate-isotope regression model

Complex speleothem growth processes revealed by trace element mapping and scanning electron microscopy of annual layers

Closely-spaced transects measured by excimer laser ablation inductively-coupled plasma mass spectrometry (ELA-ICPMS) at 5 and 32 μm spatial resolution are used to generate trace element composition maps (Ba, Sr, Mg, U, Na, P and Al) from MND-S1, a previ