Modelling the effects of porous media deformation on the propagation of water-table waves in a sandy unconfined aquifer

This paper examines the influence of porous media deformation on water-table wave dispersion in an unconfined aquifer using a numerical model which couples Richards’ equation to the poro-elastic model. The study was motivated by the findings of Shoushtari et al. (J Hydrol 533:412–440, 2016) who were unable to reproduce the observed wave dispersion in their sand flume data with either numerical Richards’ equation models (assuming rigid porous media) or existing analytic solutions. The water-table wave dispersion is quantified via the complex wave number extracted from the predicted amplitude and phase profiles. A sensitivity analysis was performed to establish the influence of the main parameters in the poro-elastic model, namely Young’s modulus (E) and Poisson’s ratio (ν). For a short oscillation period (T = 16.4 s), the phase lag increase rate (k_i) is sensitive to the chosen values of E and ν, demonstrating an inverse relationship with both parameters. Changes in the amplitude decay rate (k_r), however, were negligible. For a longer oscillation period (T = 908.6 s), variations in the values of E and ν resulted in only small changes in both k_r and k_i. In both the short and long period cases, the poro-elastic model is unable to reproduce the observed wave dispersion in the existing laboratory data. Hence porous media deformation cannot explain the additional energy dissipation in the laboratory data. Shoushtari SMH, Cartwright N, Perrochet P, Nielsen P (2016) The effects of oscillation period on groundwater wave dispersion in a sandy unconfined aquifer: sand flume experiments and modelling. J Hydrol 533:412–440.