Dr Gabriel Rau

Knowledge of fracture properties and associated flow processes is important for geoscience applications such as nuclear waste disposal, geothermal energy and hydrocarbons. An important tool established in recent years are hydro-mechanical models which provide a useful alternative to experimental methods determining single fracture parameters such as hydraulic aperture. A crucial issue for meaningful numerical modeling is precise imaging of the fracture surfaces to capture geometrical information. Hence, we apply and compare three distinct fracture surface imaging methods: (1) handheld laser scanner (HLS), (2) mounted laser scanner (MLS) and (3) Structure from Motion (SfM) to a bedding plane fracture of sandstone. The imaging reveals that the resolution of the fracture surface obtained from handheld laser scanner (HLS) is insufficient for any numerical simulations, which was therefore rejected. The remaining surfaces are subsequently matched and the resulting fracture dataset is used for detailed fracture flow simulations. The resulting hydraulic aperture is calibrated with laboratory measurements using a handheld air permeameter. The air permeameter data provide a hydraulic aperture of 81 ± 1¿µm. For calibration, mechanical aperture fields are calculated using stepwise increasing contact areas up to 15%. At 5% contact area, the average hydraulic aperture obtained by MLS (85¿µm) is close to the measurement. For SfM, the measurements are fitted at 7% contact area (83¿µm). The flow simulations reveal preferential flow through major channels that are structurally and geometrically predefined. Thus, this study illustrates that resolution and accuracy of the imaging device strongly affect the quality of fluid flow simulations and that SfM provides a promising low-cost method for fracture imaging on cores or even outcrops.

When modelling heat transport in hydrogeological systems, a standard assumption is local thermal equilibrium (LTE), which implies that the porous medium temperature at the interface between the solid and fluid instantaneously reaches equilibrium. Few studies have investigated the validity of the LTE assumption and its violation, also known as local thermal non-equilibrium (LTNE), as well as its impact on the accuracy of heat transport models. While the theoretical conditions under which LTNE occurs in natural groundwater flow have recently been revealed, these have not yet been experimentally verified. We examined the applicability of the LTE assumption by conducting systematic laboratory experiments using eight distinct flow velocities (Reynolds number, Re < 0.37) through sand (0.76 mm grain size), injecting both heat and solute as tracers and observing the response at multiple points downstream. Theoretical heat transport parameters were calculated from solute tracer experiments by applying a retardation factor, and the results were subsequently compared to the actual parameters estimated from heat tracing. Our experimental results confirm that the LTE assumption can be violated under natural groundwater flow conditions, and that LTNE impacts thermal dispersion coefficients more significantly than thermal front velocities. The largest differences between the predicted and estimated thermal transport parameters were 112% and 13% for the dispersion coefficients and velocities, respectively. Our results demonstrate that LTNE effects need to be considered in heat transport modelling especially when analyzing thermal dispersion coefficients.

Urban greening can potentially help mitigate heat-related mortality and flooding facing the >4 billion urban population worldwide. However, the geographical variation of the relative combined hydrological and thermal performance benefits of such interventions are unknown. Here we quantify globally, using a hydrological model, how climate-driven trade-offs exist between hydrological retention and cooling potential of urban greening such as green roofs and parks. Using a Budyko framework, we show that water retention generally increases with aridity in water-limited environments, while cooling potential favors energy-limited climates. Our models suggest that common urban greening strategies cannot yield high performance simultaneously for addressing both urban heat-island and urban flooding problems in most cities globally. Irrigation, if sustainable, may enhance cooling while maintaining retention performance in more arid locations. Increased precipitation variability with climate change may reduce performance of thinner green-infrastructure more quickly compared to greened areas with thicker soils and root systems. Our results provide a conceptual framework and first-order quantitative guide for urban development, renewal and policymaking.

A groundwater recharge investigation in the arid zone of Australia is presented. The investigation used a wide range of hydrogeological techniques including geological mapping, surface and borehole geophysics, groundwater hydraulics, streambed temperature and pressure monitoring, and hydrogeochemical and environmental tracer sampling, and it was complemented by analysis of rainfall intensity from 18 tipping-bucked rain gauges, climate data and stream runoff measurements. Run-off and recharge from a 200-mm rainfall event in January 2015, the largest daily rainfall in the local 50-year record, were investigated in detail. While this major storm provided substantial run-off as a potential source for focused, indirect recharge, it only produced enough actual recharge to the shallow aquifer to temporarily halt a long-term groundwater recession. A series of smaller rainfall-runoff events in 2016 produced a similar recharge response. The results suggest that the total magnitude of a flood event is not the main control on indirect groundwater recharge at this location. A deeper aquifer shows no hydraulic response to surface-water flow events and is isolated from the shallow system, consistent with its Pleistocene groundwater age. This supports a growing body of evidence indicating that attributing or predicting generalised changes in recharge to changes in climate in dryland environments should not be attempted without first unravelling the dynamic processes governing groundwater recharge in the locality of interest. The results should prompt more detailed and long-term field investigation in other arid zone locations to further understand the episodic and nonlinear nature of recharge in such environments.

The groundwater pressure response to the ubiquitous Earth and atmospheric tides provides a largely untapped opportunity to passively characterize and quantify subsurface hydro-geomechanical properties. However, this requires reliable extraction of closely spaced harmonic components with relatively subtle amplitudes but well-known tidal periods from noisy measurements. The minimum requirements for the suitability of existing groundwater records for analysis are unknown. This work systematically tests and compares the ability of two common signal processing methods, the discrete Fourier transform (DFT) and harmonic least squares (HALS), to extract harmonic component properties. First, realistic conditions are simulated by analyzing a large number of synthetic data sets with variable sampling frequencies, record durations, sensor resolutions, noise levels and data gaps. Second, a model of two real-world data sets with different characteristics is validated. The results reveal that HALS outperforms the DFT in all aspects, including the ability to handle data gaps. While there is a clear trade-off between sampling frequency and record duration, sampling rates should not be less than six samples per day and records should not be shorter than¿20 days when simultaneously extracting tidal constituents. The accuracy of detection is degraded by increasing noise levels and decreasing sensor resolution. However, a resolution of the same magnitude as the expected component amplitude is sufficient in the absence of excessive noise. The results provide a practical framework to determine the suitability of existing groundwater level records and can optimize future groundwater monitoring strategies to improve passive characterization using tidal signatures.

Heat transport in natural porous media, such as aquifers or streambeds, is generally modeled assuming local thermal equilibrium (LTE) between the fluid and solid phases. Yet, the mathematical and hydrogeological conditions and implications of this simplification have not been fully established for natural porous media. To quantify the occurrence and effects of local thermal disequilibrium during heat transport, we systematically compared thermal breakthrough curves from a LTE with those calculated using a local thermal nonequilibrium (LTNE) model, explicitly allowing for different temperatures in the fluid and solid phases. For the LTNE model, we developed a new correlation for the heat transfer coefficient representative of the conditions in natural porous aquifers using six published experimental results. By conducting an extensive parameter study (>50,000 simulations), we show that LTNE effects do not occur for grain sizes smaller than 7¿mm or for groundwater flow velocities that are slower than 1.6¿m¿day-1. The limits of LTE are likely exceeded in gravel aquifers or in the vicinity of pumped bores. For such aquifers, the use of a LTE model can lead to an underestimation of the effective thermal dispersion by a factor of up to 30 or higher, while the advective thermal velocity remains unaffected for most conditions. Based on a regression analysis of the simulation results, we provide a criterion which can be used to determine if LTNE effects are expected for particular conditions.

The groundwater response to Earth tides and atmospheric pressure changes can be used to understand subsurface processes and estimate hydraulic and hydro-mechanical properties. We develop a generalised frequency domain approach to disentangle the impacts of Earth and atmospheric tides on groundwater level responses. By considering the complex harmonic properties of the signal, we improve upon a previous method for quantifying barometric efficiency (BE), while simultaneously assessing system confinement and estimating hydraulic conductivity and specific storage. We demonstrate and validate this novel approach using an example barometric and groundwater pressure record with strong Earth tide influences. Our method enables improved and rapid assessment of subsurface processes and properties using standard pressure measurements.

Groundwater is an important global resource and its sustainable use faces major challenges. New methods and advances in computational science could lead to much improved understanding of groundwater processes and subsurface properties. A closer look at current groundwater monitoring practice reveals the need for updates with a special focus on the benefits of high-frequency and high-resolution datasets. To future-proof hydrogeology, this technical note raises awareness about the necessity for improvement, provides initial recommendations and advocates for the development of universal guidelines.

Study Region: Accra Plains. Study Focus: We conducted a field geology mapping, a well inventory exercise, used ERT, drilled boreholes at 8 locations (15¿60 m depth), took drill core samples which we subjected to dilute acid leaching experiments, installed piezometers and equipped them with pressure transducers, analyzed tidal signals in high frequency groundwater hydrographs, carried out pumping tests, and, finally, took 49 groundwater samples. New Insights: Our results indicated a diverse groundwater system. On the one hand, groundwater was found at shallow depths in the saprolite of the Togo Structural Unit (TSU), which, in unweathered state, is composed of phyllites, schists, and quartzites. This system was shallow and predominantly unconfined, as revealed by tidal analysis. In addition, transmissivities of TSU saprolite, all in the order of < 6e-5 m2/s, reduced with depth, which indicated either the lack of a stratiform fractured layer or the presence of such layer beyond drilled depths. On the other hand, groundwater was found in fractures of the gneisses of the Dahomeyan Structural Unit (DSU). This system was potentially deeper, but DSU transmissivities were significantly lower than those of TSU saprolite. Hydrochemically, groundwater was mainly characterized by infiltration of waste-water, evidenced by elevated ion concentrations, including nitrate. Due to the thinly weathered basement, groundwater system development in the area is local and restricted to the Dodowa area.

Groundwater responses to barometric pressure fluctuations are characterized using the concept of barometric efficiency (BE). For semiconfined and confined aquifers, BE values can be used to provide efficient, low-cost estimates of specific storage. This study compares, for the first time, eight existing methods of BE estimation. Comparisons were undertaken using data from the Peel region of Western Australia. Fourier analysis and regression deconvolution methods were used to estimate aquifer confinement status. The former approach was found to be robust and provided a quantitative basis for spatial comparisons of the degree of confinement. The latter approach was confounded by the presence of diurnal and/or semidiurnal signals. For wells at which semiconfined or confined responses were identified, frequency and time domain methods were used to estimate BE values. Most BE estimation methods were similarly confounded by diurnal and/or semidiurnal signals, with the exception of the Acworth et al. (2016) method. Specific storage values calculated from BE values were order-of-magnitude consistent with the results of four historical pumping tests. The methods implemented in this research provide efficient, low-cost alternatives to hydraulic testing for estimating aquifer confinement, as well as the BE and specific storage of semiconfined and confined aquifers. The frequency and duration of observations required by these methods are minimal; for example, typically requiring a minimum of four observations per day over a four month period. In some locations they may allow additional insights to be derived from existing groundwater hydrograph data.

Hydraulic head and gradient measurements underpin practically all investigations in hydrogeology. There is sufficient information in the literature to suggest that head measurement errors can impede the reliable detection of flow directions and significantly increase the uncertainty of groundwater flow rate calculations. Yet educational textbooks contain limited content regarding measurement techniques, and studies rarely report on measurement errors. The objective of our study is to review currently accepted standard operating procedures in hydrological research and to determine the smallest head gradients that can be resolved. To this aim, we first systematically investigate the systematic and random measurement errors involved in collecting time-series information on hydraulic head at a given location: (1) geospatial position, (2) point of head, (3) depth to water, and (4) water level time series. Then, by propagating the random errors, we find that with current standard practice, horizontal head gradients < 10-4 are resolvable at distances =170 m. Further, it takes extraordinary effort to measure hydraulic head gradients< 10-3 over distances < 10 m. In reality, accuracy will be worse than our theoretical estimates because of the many possible systematic errors. Regional flow on a scale of kilometres or more can be inferred with current best-practice methods, but processes such as vertical flow within an aquifer cannot be determined until more accurate and precise measurement methods are developed. Finally, we offer a concise set of recommendations for water level, hydraulic head and gradient time-series measurements. We anticipate that our work contributes to progressing the quality of head time-series data in the hydrogeological sciences and provides a starting point for the development of universal measurement protocols for water level data collection.

Groundwater extraction is increasing rapidly in many areas of the world, causing serious impacts such as falling water tables, ground surface subsidence, water quality degradation, and reduction of stream baseflow on which many ecosystems depend. Methods for understanding and predicting the impacts of groundwater extraction generally lack detailed spatial and temporal knowledge of the subsurface hydrogeomechanical properties. This review provides a comprehensive understanding of Earth and atmospheric tides and their impact on subsurface pore pressure. First, we evaluate the global occurrence of Earth and atmospheric tides. Then, we illustrate their impact on the groundwater response and connect this with the theory of poroelasticity, which underpins quantitative analyses. Finally, we review methods that utilize these impacts to characterize groundwater systems and to quantify their hydrogeomechanical properties. We conclude by highlighting their potential as passive and low-cost investigation techniques and by outlining the research and developments required to progress and make analyses readily available. Thus, hydrogeomechanical properties of subsurface systems could be obtained at unprecedented spatial and temporal resolution, adding additional value to commonly acquired groundwater and atmospheric pressure data.

Clay-rich deposits are commonly assumed to be aquitards which act as natural hydraulic barriers due to their low hydraulic connectivity. Postdepositional weathering processes are known to increase the permeability of aquitards in the near surface but not impact on deeper parts of relatively thick formations. However, syndepositional processes affecting the hydraulic properties of aquitards have previously received little attention in the literature. Here, we analyze a 31 m deep sediment core recovered from an inland clay-rich sedimentary sequence using a combination of techniques including particle size distribution and microscopy, centrifuge dye tracer testing and micro X-ray CT imaging. Subaerial deposition of soils within these fine grained alluvial deposits has led to the preservation of considerable macropores (root channels or animal burrows). Connected pores and macropores thus account for vertical hydraulic conductivity (K) of 4.2×10-1 m/s (geometric mean of 13 samples) throughout the thick aquitard, compared to a matrix K that is likely < 10-10 m/s, the minimum K value that was measured. Our testing demonstrates that such syndepositional features may compromise the hydraulic integrity of what otherwise appears to have the characteristics of a much lower permeability aquitard. Heterogeneity within a clay-rich matrix could also enhance vertical connectivity, as indicated by digital analysis of pore morphology in CT images. We highlight that the paleo-environment under which the sediment was deposited must be considered when aquitards are investigated as potential natural hydraulic barriers and illustrate the value of combining multiple investigation techniques for characterizing clay-rich deposits.

Solute dispersion and mixing in beach aquifers is strongly influenced by highly transient flow induced from wave forcing. While transport at the groundwater-ocean interface has been modelled, little is known about the quantitative effect of wave forcing on solute dispersion and mixing in beach aquifers. We use a prototype-scale laboratory flume experiment to conduct tracer transport experiments at two locations within a wave-forced beach aquifer. For the first time we demonstrate by systematic laboratory experimentation that transient conditions in the re-circulation zone due to run-up and beach face infiltration as well as head oscillations caused by wave forcing strongly disperse and mix subsurface solute plumes. Wave forcing can increase the apparent dispersion by an order of magnitude depending on conditions, compared to solute transport without waves. Our findings illustrate that beach aquifer transport models need to consider the additional dispersion to correctly quantify mixing and biogeochemical processes in this highly dynamic zone.

Groundwater specific storage varies by orders of magnitude, is difficult to quantify, and prone to significant uncertainty. Estimating specific storage using aquifer testing is hampered by the nonuniqueness in the inversion of head data and the assumptions of the underlying conceptual model. We revisit confined poroelastic theory and reveal that the uniaxial specific storage can be calculated mainly from undrained poroelastic properties, namely, uniaxial bulk modulus, loading efficiency, and the Biot-Willis coefficient. In addition, literature estimates of the solid grain compressibility enables quantification of subsurface poroelastic parameters using field techniques such as cross-hole seismic surveys and loading efficiency from the groundwater responses to atmospheric tides. We quantify and compare specific storage depth profiles for two field sites, one with deep aeolian sands and another with smectitic clays. Our new results require bulk density and agree well when compared to previous approaches that rely on porosity estimates. While water in clays responds to stress, detailed sediment characterization from a core illustrates that the majority of water is adsorbed onto minerals leaving only a small fraction free to drain. This, in conjunction with a thorough analysis using our new method, demonstrates that specific storage has a physical upper limit of (Formula presented.) m-1. Consequently, if larger values are derived using aquifer hydraulic testing, then the conceptual model that has been used needs reappraisal. Our method can be used to improve confined groundwater storage estimates and refine the conceptual models used to interpret hydraulic aquifer tests.

Accurate determination of groundwater state of confinement and compressible storage properties at vertical resolution over depth is notoriously difficult. We use the hydraulic head response to atmospheric tides at 2 cpd frequency as a tracer to quantify barometric efficiency (BE) and specific storage (Ss) over depth. Records of synthesized Earth tides, atmospheric pressure, and hydraulic heads measured in nine piezometers completed at depths between 5 and 55 m into unconsolidated smectitic clay and silt, sand and gravel were examined in the frequency domain. The barometric efficiency increased over depth from ~0.05 in silty clay to ~0.15 in sands and gravels. BE for silty clay was confirmed by calculating the loading efficiency as 0.95 using rainfall at the surface. Specific storage was calculated using effective rather than total moisture. The differences in phase between atmospheric pressure and hydraulic heads at 2 cpd were ~180° below 10 m indicating confined conditions despite the low BE. Heads in the sediment above a fine sand and silt layer at 12 m exhibited a time variable phase difference between 0° and 180° indicating varying confinement. Our results illustrate that the atmospheric tide at 2 cpd is a powerful natural tracer for quantifying groundwater state of confinement and compressible storage properties in layered formations from hydraulic heads and atmospheric pressure records without the need for externally induced hydraulic stress. This approach could significantly improve the development of conceptual hydrogeological model used for groundwater resource development and management.

Ephemeral and intermittent flow in dryland stream channels infiltrates into sediments, replenishes groundwater resources and underpins riparian ecosystems. However, the spatiotemporal complexity of the transitory flow processes that occur beneath such stream channels are poorly observed and understood. We develop a new approach to characterise the dynamics of surface water-groundwater interactions in dryland streams using pairs of temperature records measured at different depths within the streambed. The approach exploits the fact that the downward propagation of the diel temperature fluctuation from the surface depends on the sediment thermal diffusivity. This is controlled by time-varying fractions of air and water contained in streambed sediments causing a contrast in thermal properties. We demonstrate the usefulness of this method with multi-level temperature and pressure records of a flow event acquired using 12 streambed arrays deployed along a ~ 12 km dryland channel section. Thermal signatures clearly indicate the presence of water and characterise the vertical flow component as well as the occurrence of horizontal hyporheic flow. We jointly interpret thermal signatures as well as surface and groundwater levels to distinguish four different hydrological regimes: [A] dry channel, [B] surface run-off, [C] pool-riffle sequence, and [D] isolated pools. The occurrence and duration of the regimes depends on the rate at which the infiltrated water redistributes in the subsurface which, in turn, is controlled by the hydraulic properties of the variably saturated sediment. Our results have significant implications for understanding how transitory flows recharge alluvial sediments, influence water quality and underpin dryland ecosystems.

Temperature and moisture content in the variably saturated subsurface are two of the most important physical parameters that govern a wide variety of geochemical and ecological processes. An understanding of thermal and hydraulic processes and properties of transient vadose zones is therefore fundamental in the evaluation of such processes. Here, an investigation of the thermal regime and subsurface properties of a tidally affected, variably saturated streambed is presented. Field and laboratory measurements, as well as a forward numerical model, are jointly employed in the investigation. Temperature, soil moisture, surface level, and water level data were recorded in a transect perpendicular to a tidally driven stream. Frequency-domain analysis of the subsurface temperature measurements revealed the rapid decay of the tidal temperature driver within the top ~30¿cm of sediment. Several techniques were used to evaluate subsurface thermal and hydraulic properties, including thermal conductivity and the soil water retention curve. These properties were used to constrain a forward numerical model that included coupled treatment of relevant variable saturation thermal and hydraulic physics. Even though the investigated vadose zone is intermittent and relatively shallow 20¿cm, the results illustrate how error can be introduced into heat-transport calculations if unsaturated conditions are not taken into account.

Thermal dispersion, caused by fluid velocity and temperature fluctuations in the pore space and the effects of hydrodynamic mixing on the temperature field, controls convective heat transport in saturated porous media. While the thermal dispersion coefficient, a governing parameter in the thermal equilibrium model (TEM), has been investigated for natural systems, the dependence of the thermal dispersion coefficient on particle size remains undetermined. Previous research found that the relationship between the thermal dispersion coefficient and flow velocity follows a power law and that there may be a temperature difference between the solid and fluid phase (thermal non-equilibrium). However, experiments are limited to discrete particle sizes and comparison of the dispersion-velocity relationship is impeded by different experimental approaches. We conducted a series of separate heat and solute transport experiments in a column filled with uniform porous media consisting of different sized glass spheres for a range of flow velocities. The thermal and solute dispersion coefficients obtained from experimental measurements were correlated with flow velocities through the thermal or solute Péclet number. Our results demonstrate that, while solute dispersion is independent of particle size, the dependence of the TEM based thermal dispersion coefficient on flow rates is influenced by the particle size. This is caused by the fact that, unlike solute transport, heat exchanges between fluid and particles and that this induces thermal non-equilibrium between both phases. The results have significant implications for quantifying forced convective heat transport in natural porous media because thermal non-equilibrium between the phases is not considered. The porous media particle size must be considered when selecting appropriate values for the thermal dispersion coefficient.

Soil moisture and temperature are some of the most important controls for a wide variety of geochemical and ecological processes in the vadose zone (VZ). Soil moisture is highly variable both spatially and temporally. The development of methods to measure it on various scales has been the subject of much activity. Recently, geoscientists have been increasingly interested in measuring temperature as a proxy for hydrologic properties and parameters, including soil moisture. Here, we discuss the motivation, primary concepts, equipment, and fundamental thermal and hydraulic models related to heat and water transport in variably saturated porous media. A large variety of methods for heat tracing, including both passive and active-heating methodologies, are detailed. Heat tracing methods offer the capacity to measure soil moisture on a scale from ~¿1¿cm up to several km using temperature, a parameter whose measurement in VZ studies is often required anyway due to its effect on many subsurface processes. Furthermore, heat-tracing methods are not affected by high salinity pore water that can limit electromagnetic soil moisture methods. We also review coupled thermo-hydro VZ modelling software and VZ thermal regime studies and identify several knowledge gaps. With the intention to serve as an introduction to VZ heat-tracing, this review consolidates recent advances and outlines potential themes for future research.

Understanding and managing groundwater resources in drylands is a challenging task, but one that is globally important. The dominant process for dryland groundwater recharge is thought to be as focused, indirect recharge from ephemeral stream losses. However, there is a global paucity of data for understanding and quantifying this process and transferable techniques for quantifying groundwater recharge in such contexts are lacking. Here we develop a generalized conceptual model for understanding water table and groundwater head fluctuations due to recharge from episodic events within ephemeral streams. By accounting for the recession characteristics of a groundwater hydrograph, we present a simple but powerful new water table fluctuation approach to quantify focused, indirect recharge over both long term and event time scales. The technique is demonstrated using a new, and globally unparalleled, set of groundwater observations from an ephemeral stream catchment located in NSW, Australia. We find that, following episodic streamflow events down a predominantly dry channel system, groundwater head fluctuations are controlled by pressure redistribution operating at three time scales from vertical flow (days to weeks), transverse flow perpendicular to the stream (weeks to months), and longitudinal flow parallel to the stream (years to decades). In relative terms, indirect recharge decreases almost linearly away from the mountain front, both in discrete monitored events as well as in the long-term average. In absolute terms, the estimated indirect recharge varies from 80 to 30 mm/a with the main uncertainty in these values stemming from uncertainty in the catchment-scale hydraulic properties.

The groundwater hydraulic head response to the worldwide and ubiquitous atmospheric tide at 2 cycles per day (cpd) is a direct function of confined aquifer compressible storage. The ratio of the responses of hydraulic head to the atmospheric pressure change is a measure of aquifer barometric efficiency, from which formation compressibility and aquifer specific storage can be determined in situ rather than resorting to laboratory or aquifer pumping tests. The Earth tide also impacts the hydraulic head response at the same frequency, and a method is developed here to quantify and remove this interference. As a result, the barometric efficiency can be routinely calculated from 6-hourly hydraulic head, atmospheric pressure, and modeled Earth tide records where available for a minimum of 15¿days duration. This new approach will be of critical importance in assessing worldwide problems of land subsidence or groundwater resource evaluation that both occur due to groundwater abstraction.

Measurement of glycerol dialkyl glycerol tetraethers (GDGTs) preserved in speleothems offers a potential proxy for past temperature but, in general, their origin is unknown. To understand the source of speleothem GDGTs, we undertook an irrigation experiment to activate drip sites within a hydrogeochemically well characterised cave. The cave drip water was analysed for GDGTs, inorganic elements (major ions and trace elements), stable isotopes and dissolved organic matter concentration and character. Published speleothem GDGT records from the site have been observed to be dominated by isoprenoid GDGTs and interpreted as deriving from in situ microbial communities within the cave or vadose zone. The drip water in our irrigation experiment had a GDGT distribution distinct from that of soil and speleothem samples, providing direct evidence that the distinctive GDGT signature in speleothems is derived from a subsurface source. Analysis of GDGTs in this context allowed further elucidation of their source and transport in cave systems, enhancing our understanding of how they might be used as a temperature proxy.

A novel semi-analytical model for the calculation of water saturation levels in the near subsurface using passive temperature measurements is derived. The amplitude and phase of dominant natural diel temperature variations are exploited, although the solution is general so that a cyclical temperature signal of any period could be used. The model is based on the first-principles advection-conduction-dispersion equation, which is fully general for porous media. It requires a single independent soil moisture estimate, but directly considers the spatially variable saturation dependency of thermal conductivity which has been avoided in previous studies. An established empirical model for the thermal conductivity of variably saturated porous media is incorporated and two solutions for saturation are derived. Using data from numerical models, a spatially sequential implementation of one of these solutions is shown to predict the vertical saturation profile to within 2% for a hydraulically stable case and to within the saturation range observed over a single day for percolation rates up to 10 cm/day. The developed model and methodology can aid in the analysis of archived temperature data from the vadose zone and will serve as a powerful tool in future heat-tracing experiments in variably saturated conditions.

A new method was developed for analysing and delineating streambed water fluxes, flow conditions and hydraulic properties using coiled fibre-optic distributed temperature sensing or closely spaced discrete temperature sensors. This method allows for a thorough treatment of the spatial information embedded in temperature data by creating a matrix visualization of all possible sensor pairs. Application of the method to a 5-day field dataset reveals the complexity of shallow streambed thermal regimes. To understand how velocity estimates are affected by violations of assumptions of one-dimensional, saturated, homogeneous flow and to aid in the interpretation of field observations, the method was also applied to temperature data generated by numerical models of common field conditions: horizontal layering, presence of lateral flow and variable streambed saturation. The results show that each condition creates a distinct signature visible in the triangular matrices. The matrices are used to perform a comparison of the behaviour of one-dimensional analytical heat-tracing models. The results show that the amplitude ratio-based method of velocity calculation leads to the most reliable estimates. The minimum sensor spacing required to obtain reliable velocity estimates with discrete sensors is also investigated using field data. The developed method will aid future heat-tracing studies by providing a technique for visualizing and comparing results from fibre-optic distributed temperature sensing installations and testing the robustness of analytical heat-tracing models. Copyright © 2016 John Wiley & Sons, Ltd.

We present results of a detailed study of drip rate variations at 12 drip discharge sites in Glory Hole Cave, New South Wales, Australia. Our novel time series analysis, using the wavelet synchrosqueezed transform, reveals pronounced oscillations at daily and sub-daily frequencies occurring in 8 out of the 12 monitored sites. These oscillations were not spatially or temporally homogenous, with different drip sites exhibiting such behaviour at different times of year in different parts of the cave. We test several hypotheses for the cause of the oscillations, including variations in pressure gradients between karst and cave due to cave breathing effects or atmospheric and earth tides, variations in hydraulic conductivity due to changes in viscosity of water with daily temperature oscillations, and solar-driven daily cycles of vegetative (phreatophytic) transpiration. We conclude that the only hypothesis consistent with the data and hydrologic theory is that daily oscillations are caused by solar-driven pumping by phreatophytic trees which are abundant at the site. The daily oscillations are not continuous and occur sporadically in short bursts (2-14 days) throughout the year due to non-linear modification of the solar signal via complex karst architecture. This is the first indirect observation leading to the hypothesis of tree water use in cave drip water. It has important implications for karst hydrology in regards to developing a new protocol to determine the relative importance of trends in drip rate, such as diurnal oscillations, and how these trends change over timescales of weeks to years. This information can also be used to infer karst architecture. This study demonstrates the importance of vegetation on recharge dynamics, information that will inform both process-based karst models and empirical estimation approaches. Our findings support a growing body of research exploring the impact of trees on speleothem paleoclimate proxies.

Quantifying dryland groundwater recharge as a function of climate variability is becoming increasingly important in the face of a globally depleted resource, yet this remains a major challenge due to lack of adequate monitoring and the complexity of processes involved. A previously unpublished and unique dataset of high density and frequency rainfall measurements is presented, from the Fowlers Gap Arid Zone Research Station in western New South Wales (Australia). The dataset confirms extreme spatial and temporal variability in rainfall distribution which has been observed in other dryland areas and is generally explained by the dominance of individual storm cells. Contrary to previous observations, however, this dataset contains only a few localised storm cells despite the variability. The implications of spatiotemporal rainfall variability on the estimation of groundwater recharge is assessed and show that the most likely recharge mechanism is through indirect and localised, rather than direct, recharge. Examples of rainfall and stream gauge height illustrate runoff generation when a spatially averaged threshold of 15¿25¿mm (depending on the antecedent moisture conditions) is exceeded. Preliminary assessment of groundwater levels illustrates that the regional water table is much deeper than anticipated, especially considering the expected magnitude of indirect and localised recharge. A possible explanation is that pathways for indirect and localised recharge are inhibited by the large quantities of Aeolian dust observed at the site. Runoff readily occurs with water collecting in surface lakes which slowly dry and disappear. Assuming direct groundwater recharge under these conditions will significantly overestimate actual recharge.

There is insufficient information on the movement of air in karst environments to constrain the uncertainty associated with quantifying sources and sinks of methane (CH4) and carbon dioxide (CO2) within caves for global carbon accounting. We analysed cave-air samples for their CO2 and CH4 concentrations ([CO2] and [CH4]) and carbon isotopic compositions from sampling campaigns in winter (August 2014) and summer (February 2015) at numerous heights and locations throughout Gaden and Cathedral caves, in a semi-arid environment, Wellington Caves, NSW, Australia. Ventilation is the dominant control on cave-air CO2 and CH4, with the highest cave-air CO2 concentrations ([CO2]cave) occurring in summer, in association with the lowest cave-air CH4 concentrations ([CH4]cave). Analyses show that the cave-air CO2 has both atmospheric and soil sources. Soil air and cave air in both caves undergo methanogenesis and methanotrophy, but we identify cave-specific differences in cave-air CH4 and CO2. [CH4]cave in Cathedral Cave shows an inverse relationship to [CO2]cave, particularly in areas separated from the main cave passage. In contrast, Gaden Cave has near-atmospheric [CH4]cave and isotopic ratios present at all locations sampled in winter. Where no ventilation is occurring in summer, [CH4]cave in Gaden Cave decreases, but remains reasonably high compared to Cathedral Cave. Our research shows adjacent caves vary in their ability to act as a net sink for CH4, and highlights the need for further studies before global generalisations can be made about the carbon budget of karst environments.