Land Subsidence Due to Groundwater Exploitation in Unconfined Aquifers: Experimental and Numerical Assessment with Computational Fluid Dynamics

Abstract

Land subsidence is a global challenge that enhances the vulnerability of aquifers where climate change and driving forces are occurring simultaneously. To comprehensively analyze this issue, integrated modeling tools are essential. This study advances the simulation of subsidence using Computational Fluid Dynamics (CFD); it assessed the effects of exploitation and recharge of groundwater on the vertical displacement of coarse and fine sands in a laboratory-scale aquifer. A model was developed by integrating the Navier–Stokes equations to study the groundwater flow and Terzaghi’s law for the vertical displacement of sands. The boundary conditions used were Dirichlet based on the changes in the hydraulic head over time. The specific storage coefficient was used to calibrate the model. The findings confirmed that subsidence occurs at slower rates in soil with fine sands with average particle diameters of 0.39 mm than in coarse sands with average particle diameters of 0.67 mm. The maximum discrepancy between the experimental and the numerical reaffirms that CFD platforms can be used to simulate subsidence dynamics and potentially allow the simultaneous simulation of other dynamics. Concluding remarks and recommendations are highlighted considering the up-to-date advances and future work to improve the research on subsidence in unconfined aquifers.

Flow Pattern and Escape Hazards of People from Flood Intrusion into the Staircase of Underground Spaces with Multiple Rest Platforms

“…The VOF method is suitable for solving the free surface tracking description problem for two or more immiscible fluids, wherein only one fluid (e.g., air) is compressible [38][39][40]. The method has been validated and is widely used in wave simulations [41][42][43][44][45], dam-break models [38,46], and multiform step-fall simulations [47][48][49]. The core idea of the VOF method is to set the volume fraction as an α function, to characterize the proportion of fluid occupancy in a cell.…”