An objective frequency domain method for quantifying confined aquifer compressible storage using Earth and atmospheric tides

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 abstractio

Vertical groundwater storage properties and changes in confinement determined using hydraulic head response to atmospheric tides

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

Modelling of C/Cl isotopic behaviour during chloroethene biotic reductive dechlorination: Capabilities and limitations of simplified and comprehensive models.

Predicting the fate of chloroethenes in groundwater is essential when evaluating remediation strategies. Such predictions are expected to be more accurate when incorporating isotopic parameters. Although secondary chlorine isotope effects have been observed during reductive dechlorination of chloroethenes, development of modelling frameworks and simulation has thus far been limited. We have developed a novel mathematical framework to simulate the C/Cl isotopic fractionation during reductive dechlorination of chloroethenes. This framework differs from the existing state of the art by incorporating secondary isotopic effects and considering both C and Cl isotopes simultaneously. A comprehensive general model (GM), which is expected to be the closest representation of reality thus far investigated, was implemented. A less computationally intensive simplified model (SM), with the potential for use in modelling of complex reactive transport scenarios, was subsequently validated based on its comparison to GM. The approach of GM considers all isotopocules (i.e. molecules differing in number and position of heavy and light isotopes) of each chloroethene as individual species, of which each is degraded at a different rate. Both models GM and SM simulated plausible C/Cl isotopic compositions of tetrachloroethene (PCE), trichloroethene (TCE) and cis-1,2-dichloroethene (cDCE) during sequential dechlorination when using experimentally relevant kinetic and isotopic parameters. The only major difference occurred in the case where different secondary isotopic effects occur at the different non-reacting positions when PCE is dechlorinated down to cDCE. This observed discrepancy stems from the unequal Cl isotope distribution in TCE that arises due to the occurrence of differential secondary Cl isotopic effects during transformation of PCE to TCE. Additionally, these models are shown to accurately reproduce experimental data obtained during reductive dechlorination by bacterial enrichments harbouring Sulfurospirillum spp. where secondary isotope effects are known to have occurred. These findings underscore a promising future for the development of reactive transport models that incorporate isotopic parameters

Assessing the accuracy of 1-D analytical heat tracing for estimating near-surface sediment thermal diffusivity and water flux under transient conditions

Amplitude decay and phase delay of oscillating temperature records measured at two vertical locations in near-surface sediments can be used to infer water fluxes, thermal diffusivity, and sediment scour/deposition. While methods that rely on the harmonics-based analytical heat transport solution assume a steady state water flux, many applications have reported transient fluxes but ignored the possible violation of this assumption in the method. Here we use natural heat tracing as an example to investigate the extent to which changes in the water flux, and associated temperature signal nonstationarity, can be separated from other influences. We systematically scrutinize the assumption of steady state flow in analytical heat tracing and test the capabilities of the method to detect the timing and magnitude of flux transients. A numerical model was used to synthesize the temperature response to different step and ramp changes in advective thermal velocity magnitude and direction for both a single-frequency and multifrequency temperature boundary. Time-variable temperature amplitude and phase information were extracted from the model output with different signal-processing methods. We show that a worst-case transient flux induces a temperature nonstationarity, the duration of which is less than 1 cycle for realistic sediment thermal diffusivities between 0.02 and 0.13 m2/d. However, common signal-processing methods introduce erroneous temporal spreading of advective thermal velocities and significant anomalies in thermal diffusivities or sensor spacing, which is used as an analogue for streambed scour/deposition. The most time-variant spectral filter can introduce errors of up to 57% in velocity and 33% in thermal diffusivity values with artifacts spanning ±2 days around the occurrence of rapid changes in flux. Further, our results show that analytical heat tracing is unable to accurately resolve highly time-variant fluxes and thermal diffusivities and does not allow for the inference of scour/depositional processes due to the limitations of signal processing in disentangling flux-related signal nonstationarities from those stemming from other sources. To prevent erroneous interpretations, hydrometric data should always be acquired in combination with temperature records