Initially observed in the semiarid basins of southwestern USA, earth fissures due to aquifer over‐exploitation are presently threatening a large number of subsiding basins in various countries worldwide. Different mechanics have been proposed to explain this process, such as differential compaction, horizontal movements, and fault reactivation. Numerical modeling and prediction of this major geohazard caused by overuse of groundwater resources are challenging because of two main requirements: shifting from the classical continuous to discontinuous geomechanics and incorporating two‐dimensional features (the earth fissures) into large three‐dimensional (3‐D) modeling domain (the subsiding basin). In this work, we proposed a novel modeling approach to simulate earth fissure generation and propagation in 3‐D complex geological settings. A nested two‐scale approach associated with an original nonlinear elastoplastic finite element/interface element simulator allows modeling the mechanics of earth discontinuities, in terms of both sliding and opening. The model is applied on a case study in Wuxi, China, where groundwater pumping between 1985 and 2004 has caused land subsidence larger than 2 m. The model outcomes highlight that the presence of a shallow (∼80 m deep) bedrock ridge crossing the Yangtze River delta is the key factor triggering the earth fissure development in this area. Bending of the alluvial deposits around the ridge tip and shear stress due to the uneven piezometric change and asymmetrical shape of the bedrock have caused the earth fissure to onset at the land surface and propagate downward to a maximum depth of about 20–30 m. Maximum sliding and opening are computed in the range of 10–40 cm, in agreement with the order of magnitude estimated in the field.

A modeling analysis is used to investigate the relative susceptibility of various hydrogeologic configurations to aseismic rupture generation due to deformation of aquifer systems accompanying groundwater pumping. An advanced numerical model (GEPS3D) is used to simulate rupture generation and propagation for three typical processes: (i) reactivation of a preexisting fault, (ii) differential compaction due to variations in thickness of aquifer/aquitard layers constituting the aquifer system, and (iii) tensile fracturing above a bedrock ridge that forms the base of the aquifer system. A sensitivity analysis is developed to address the relative importance of various factors, including aquifer depletion, aquifer thickness, the possible uneven distribution and depth below land surface of the aquifer/aquitard layers susceptible to aquifer‐system compaction, and the height of bedrock ridges beneath the aquifer system which contributes to thinning of the aquifer system. The rupture evolution is classified in two occurrences. In one, the rupture develops at either the top of the aquifer or at land surface and does not propagate. In the other, the developed rupture propagates from the aquifer top toward the land surface and/or from the land surface downward. The aquifer depth is the most important factor controlling rupture evolution. Specifically, the probability of a significant rupture propagation is higher when the aquifer top is near land surface. The numerical results are processed by a statistical regression analysis to provide a general methodology for a preliminary evaluation of possible ruptures development in exploited aquifer systems susceptible to compaction and accompanying land subsidence. A comparison with a few representative case studies in Arizona, USA, China, and Mexico supports the study outcomes.

Abstract. The Suzhou-Wuxi-Changzhou (known as "Su-Xi-Chang") area, located in the southern part of Jiangsu Province, China, experienced serious land subsidence caused by overly exploitation of groundwater. The largest cumulative land subsidence has reached 3 m. With the rapid progress of land subsidence since the late 1980s, more than 20 earth fissures developed in Su-Xi-Chang area, although no pre-existing faults have been detected in the surroundings. The mechanisms of earth fissure generation associated with excessive groundwater pumping are: (i) differential land subsidence, (ii) differences in the thickness of the aquifer system, and (iii) bedrock ridges and cliffs at relatively shallow depths. In this study, the Guangming Village Earth Fissures in Wuxi area are selected as a case study to discuss in details the mechanisms of fissure generation. Aquifer exploitation resulted in a drop of groundwater head at a rate of 5-6 m yr −1 in the 1990s, with a cumulative drawdown of 40 m. The first earth fissure at Guangming Village was observed in 1998. The earth fissures, which developed in a zone characterized by a cumulative land subsidence of approximately 800 mm, are located at the flank of a main subsidence bowl with differential subsidence ranging from 0 to 1600 mm in 2001. The maximum differential subsidence rate amounts to 5 mm yr −1 between the two sides of the fissures. The fissure openings range from 30 to 80 mm, with a cumulative length of 1000 m. Depth of bed rock changes from 60 to 140 m across the earth fissure. The causes of earth fissure generation at Guangming Village includes a decrease in groundwater levels, differences in the thickness of aquifer system, shallow depths of bedrock ridges and cliffs, and subsequent differential land subsidence.

Abstract. Shanghai city has been suffering land subsidence caused by overly exploitation of ground water since 1921, which is a serious problem for this coastal city with altitude of 2.2-4.8 m above mean sea level. The largest cumulative land subsidence amounted to 2.6 m in the downtown area. Measures to decrease the ground water exploitation, change the pumping aquifers, and increase aquifer artificial recharge have been used to mitigate land subsidence since 1961. It is necessary to develop a proper numerical model to simulate and predict land subsidence. In this study, a decoupled three-dimensional (3-D) finite element land subsidence model including a 3-D ground water flow model and a 3-D geo-mechanical model was developed to simulate the 3-D deformation of the aquifer systems in the center area of Shanghai. The area of downtown Shanghai is 660 km 2 , with 10 million inhabitants, dense high buildings, and 11 metro lines. The simulation spans the period from 1979 to 1995. Two different assumptions have been tested on the side boundary, i.e., precluding the three components of the displacement, or assuming a free-displacement condition. The distribution of calculated land subsidence and horizontal displacements in different aquifers was analyzed. The computed vertical displacement fitted well with the available observations. It has been verified that the two different assumptions on the lateral boundaries in the geo-mechanical model caused different results just limited on nodes close to boundary. The developed 3-D land subsidence model is reasonable and can be used to simulate and predict 3-D movement of aquifer systems in the center area of Shanghai, which could provide scientific support to local government in controlling land subsidence and differential movements of the land surface.

Abstract. The stress variation induced by aquifer overdraft in sedimentary basins with shallow bedrock may cause rupture in the form of pre-existing fault activation or earth fissure generation. The process is causing major detrimental effects on a many areas in China and Mexico. Ruptures yield discontinuity in both displacement and stress field that classic continuous finite element (FE) models cannot address. Interface finite elements (IE), typically used in contact mechanics, may be of great help and are implemented herein to simulate the fault geomechanical behaviour. Two main approaches, i.e. Penalty and Lagrangian, are developed to enforce the contact condition on the element interface. The incorporation of IE incorporation into a three-dimensional (3-D) FE geomechanical simulator shows that the Lagrangian approach is numerically more robust and stable than the Penalty, thus providing more reliable solutions. Furthermore, the use of a Newton-Raphson scheme to deal with the non-linear elasto-plastic fault behaviour allows for quadratic convergence. The FE -IE model is applied to investigate the likely ground rupture in realistic 3-D geologic settings. The case studies are representative of the City of Wuxi in the Jiangsu Province (China), and of the City of Queretaro, Mexico, where significant land subsidence has been accompanied by the generation of several earth fissures jeopardizing the stability and integrity of the overland structures and infrastructure.