Barometric response functions from borehole water level records and quantification of aquifer vulnerability
Abstract
In confined and semi-confined aquifers, borehole water levels respond to fluctuations in barometric pressure and this response can be used to estimate the properties of aquifer confining layers. We use this response as indicator of groundwater vulnerability for the semi-confined Chalk aquifer in East Yorkshire, UK. Time series data of borehole water levels are corrected for Earth tides and recharge, and barometric response functions (BRFs) estimated using cross-spectral deconvolution-averaging techniques. The resulting BRFs are fitted using a theoretical model of the BRF gain and phase for a semi-confined aquifer (Rojstaczer, 1988) to obtain confining layer properties. For all of the boreholes, non-zero hydraulic diffusivities for the confining layer were found indicating that the aquifer is semi-confined. A ';characteristic time scale' based on the hydraulic and pneumatic diffusivities of the confining layer is introduced as a measure of the degree of aquifer confinement and therefore groundwater vulnerability. The analytical model assumes that the confining layer and aquifer are homogeneous. However, in nature, confining layers are heterogeneous and groundwater vulnerability dominated by the presence of high diffusivity, high flow pathways through the confining layer to the aquifer. A transient numerical model (MODFLOW) was constructed to test the impact of such heterogeneities on the BRF. In the model, an observed barometric pressure time series is used as a boundary condition applied to the upper surface of the top unit of the model (representing the confining layer) and BRFs determined from the time series of model heads in the bottom unit (representing the aquifer). The results from a numerical model with a homogeneous confining layer were found to accurately reproduce the BRFs from a modified version of the analytical model. The introduction of a localized, high diffusive block in the confining layer was found to modify the BRF, reducing the gain amplitude while having limited impact on the phase. It was found that the BRF reflects the presence of a fully penetrating, high diffusivity heterogeneity up to several hundred meters distant from the observation borehole, and shows little sensitivity to the heterogeneity's horizontal dimension. Heterogeneities that are 50% partially penetrating do not significantly impact on the BRF and 90% penetrating heterogeneities can only be detected when large and close to the observation borehole. These results show that BRF gain may be particularly useful in detecting the presence of fully penetrating heterogeneities of high diffusivity within confining layers that potentially enhance groundwater vulnerability. This research has been funded in part through the EU ITN ';IMVUL' (PITN-GA-2008-212298). Reference: Rojstaczer, S. (1988) Determination of fluid-flow properties from the response of water levels in wells to atmospheric loading, Water Resources Research, 24(11), 1927-1938.
- Publication:
-
AGU Fall Meeting Abstracts
- Pub Date:
- December 2013
- Bibcode:
- 2013AGUFM.H21O..06O
- Keywords:
-
- 1829 HYDROLOGY Groundwater hydrology;
- 1872 HYDROLOGY Time series analysis;
- 1847 HYDROLOGY Modeling