Unsaturated zone hydrology and cave drip discharge water response: Implications for speleothem paleoclimate record variability

High-frequency, spatially-dense discharge monitoring was conducted over fifteen months to characterise unsaturated zone flow at Harrie Wood Cave (HWC), in the Snowy Mountains, Yarrangobilly (SE Australia). The cave was formed in the Late Silurian Yarrangobilly Limestone, a fractured rock associated with very low primary porosity due to past diagenesis. Over our monitoring period we simultaneously measured rainfall, soil moisture saturation and drip discharge rate at fourteen sites to characterise infiltration-discharge relationships. All drip discharge sites exhibited non-Gaussian distributions, indicating long periods where low discharge predominates, punctuated by short infrequent periods of high discharge. However, there was significant variability in discharge between sites and consequently no spatial correlation in the cave. We investigated the depth-discharge relationship at HWC and found a moderate relationship between depth and drip discharge lag (response) times to soil moisture content, but only weak relationships between depth and mean and maximum discharge. This highlights that the karst architecture plays an important role in controlling drip discharge dynamics. Principal Component Analysis (PCA) and Agglomerative Hierarchal Clustering (AHC) were used to classify similar drip types, revealing five unique drip regimes. Two-phase flow and non-linear response to recharge behaviour were observed, suggesting secondary porosity is controlling unsaturated zone flow in mature limestone environments with low primary porosity. Using the data presented here, the first coupled conceptual and box hydrological flow model was developed. This study highlights the heterogeneous nature of hydrological flow in karst and the need to understand unsaturated zone hydrology at the individual drip discharge level, to inform speleothem studies for high-resolution paleoclimate reconstruction.