Overexploitation of groundwater resources in the faulted basin of Querétaro, Mexico: A 3D deformation and stress analysis

Abstract The City of Queretaro is located on a continental basin filled since the Oligocene with lacustrine and alluvial sediments, pyroclastic deposits, and interbedded fractured basalts. The graben structure of the basin was formed by two major North-South trending normal faults, among which the thicknesses of the filling materials vary many tens of meters in close distances. Hence, important differences of hydraulic and mechanical properties characterize the various geologic units. Groundwater has been strongly withdrawn over the last three decades in the study area, with a decline of the piezometric level exceeding 100 m. Because of the high variability of the geologic deposits, the piezometric decrease and consequently the effective stress increase are characterized by a large spatial variability. Piezometric variations are also due to faults that strongly impact on groundwater flow dynamics. The deformation and effective stress variability has caused large differential subsidence causing ground fracturing that has damaged the urban infrastructure of the City of Queretaro. This complex geological setting has been properly accounted into a three-dimensional (3D) flow and geomechanical modeling approach to quantify the displacement, deformation, and stress fields caused by water withdrawal. The static geologic model was accurately defined using geological logs from extraction wells, field mapping of faults, fractures, and the integration of major structures reported in previous geophysical works. The model has been calibrated using observed groundwater levels and land settlement records. The simulations have spanned the period from 1970 to 2011. The modeling results highlight that the areas where large differential subsidence and horizontal displacements developed correspond to the portions of the city where ground fractures are observed. Normal and shear components of the stress field changes accumulate along the discontinuity surfaces at depth, providing evidence that large piezometric declines can be a key factor triggering the fault reactivation. The spatial relationship between major withdrawals, discontinuities of the geologic structure, and accumulation of large stress and strain fields clearly emerged from the outcomes of the 3D geomechanical model.

[1]  Andrea Franceschini,et al.  A novel Lagrangian approach for the stable numerical simulation of fault and fracture mechanics , 2016, J. Comput. Phys..

[2]  D. Carreon-Freyre,et al.  Groundwater flow in a transboundary fault-dominated aquifer and the importance of regional modeling: the case of the city of Querétaro, Mexico , 2016, Hydrogeology Journal.

[3]  G. Gambolati,et al.  Geomechanics of subsurface water withdrawal and injection , 2015 .

[4]  A. Mayer,et al.  Estimation of fault-zone conductance by calibration of a regional groundwater flow model: Desert Hot Springs, California , 2007 .

[5]  J. Pacheco,et al.  Delimitation of ground failure zones due to land subsidence using gravity data and finite element modeling in the Querétaro valley, México , 2006 .

[6]  M. Biot.  General Theory of Three‐Dimensional Consolidation , 1941 .

[7]  Margarita López Martínez,et al.  Quaternary intra-arc extension in the central Trans-Mexican volcanic belt , 2001 .

[8]  A. F. Nieto-Samaniego,et al.  Structural analysis of a relay ramp in the Querétaro graben, central Mexico: Implications for relay ramp development , 2011 .

[9]  J. Pacheco,et al.  The Sanfandila earthquake sequence of 1998, Queretaro, Mexico: activation of an undocumented fault in the northern edge of central Trans-Mexican Volcanic Belt , 2003 .

[10]  Laura Luna-González,et al.  Influencia de la estratigrafía y estructura geológica en el flujo de agua subterránea del Valle de Querétaro , 2005 .

[11]  Giuseppe Gambolati,et al.  Groundwater pumping and land subsidence in the Emilia‐Romagna coastland, Italy: Modeling the past occurrence and the future trend , 2006 .

[12]  T. Burbey.  Three-dimensional deformation and strain induced by municipal pumping, Part 2: Numerical analysis , 2006 .

[13]  Xun Jiao,et al.  Three-dimensional numerical modeling of land subsidence in Shanghai, China , 2016, Hydrogeology Journal.

[14]  M. Hernández-Marín,et al.  Controls on initiation and propagation of pumping-induced earth fissures: insights from numerical simulations , 2010 .

[15]  Lianyang Zhang.  Evaluation of rock mass deformability using empirical methods – A review , 2017 .

[16]  Á. F. Nieto-Samaniego,et al.  Estratigrafía y deformación extensional en la región San Miguel de Allende-Querétaro, México , 2001 .

[17]  F. Tupin,et al.  Time series analysis of Mexico City subsidence constrained by radar interferometry , 2009 .

[18]  E. Chaussard,et al.  Land subsidence in central Mexico detected by ALOS InSAR time-series , 2014 .

[19]  C. Janna,et al.  A Novel Approach to Model Earth Fissure Caused by Extensive Aquifer Exploitation and its Application to the Wuxi Case, China , 2018 .

[20]  Falk Amelung,et al.  Land Subsidence in Las Vegas, Nevada, 1935–2000: New Geodetic Data Show Evolution, Revised Spatial Patterns, and Reduced Rates , 2002 .

[21]  Howard A. Zebker,et al.  Inverse modeling of interbed storage parameters using land subsidence observations, Antelope Valley, California , 2003 .

[22]  M. López‐Martínez,et al.  The Amazcala caldera, Queretaro, Mexico. Geology and geochronology , 2001 .

[23]  Xianfang Song,et al.  Simulation of fully coupled finite element analysis of nonlinear hydraulic properties in land subsidence due to groundwater pumping , 2015, Environmental Earth Sciences.

[24]  M. Cerca,et al.  Assessment of groundwater flow in volcanic faulted areas. A study case in Queretaro, Mexico , 2015 .

[25]  Timothy H. Dixon,et al.  Space geodetic imaging of rapid ground subsidence in Mexico City , 2008 .

[26]  Mohammad Reza Nikudel,et al.  Numerical simulation and prediction of regional land subsidence caused by groundwater exploitation in the southwest plain of Tehran, Iran , 2016 .

[27]  W. Hung,et al.  Modeling aquifer-system compaction and predicting land subsidence in central Taiwan , 2012 .

[28]  T. Burbey,et al.  Review: Regional land subsidence accompanying groundwater extraction , 2011 .

[29]  D M Allen,et al.  Characterizing a Faulted Aquifer by Field Testing and Numerical Simulation , 1999, Ground water.

[30]  P. Teatini,et al.  Shearing along faults and stratigraphic joints controlled by land subsidence in the Valley of Queretaro, Mexico , 2016, Hydrogeology Journal.

[31]  J. Arzate,et al.  A method to predict the group fissuring and faulting caused by regional groundwater decline , 2002 .

[32]  D. Carreon Freyre.  Land subsidence processes and associated ground fracturing in central Mexico , 2010 .