Poroelastic Response to Karst Conduit Pressurization: A Finite Element Modeling Exercise Toward the Use of Tiltmeters in Karst Aquifer Monitoring Applications

Abstract

Tiltmeters have the potential to resolve ground deformation due to changes in hydraulic head induced by conduit pressurization. Conduit pressure variations cause groundwater to be stored or released from storage within the surrounding rock matrix. We modeled this process and infer whether the resulting deformation is measurable with tiltmeters and what behavior to expect by fully coupling porous media flow and solid mechanics in a poroelastic, 2D finite element model. Parameter sets globally representative of Paleozoic, Mesozoic, and Cenozoic confined and unconfined aquifers are considered. Our analysis focuses on the impact of the parameterization on pore pressure, vertical displacement, and tilt. We find that the spatial distribution of the poroelastic signal depends on the hydraulic diffusivity, and its magnitude depends on the mechanical and coupling parameters. Additional analysis of the impacts of conduit radius and depth suggests that tilt polarity could be an indicator of conduit location and relative conduit size. We calibrated the model to data observations acquired at the Santa Fe River Sink-Rise system in north-central Florida, US. We find that an overlying clay-rich layer may act to partially confine the aquifer. Although the observed tilt signal is present in radial and transverse components and polarity reversals occur, we were able to recover the magnitude and general trend of the tilt response.