Sinkhole formation and subsidence along the Dead Sea coast,Israel

More than 4,000 sinkholes have formed since the 1980s within a 60-km-long and 1-km-wide strip along the western coast of the Dead Sea (DS) in Israel. Their formation rate accelerated in recent years to >400 sinkholes per year. They cluster mostly in specific sites up to 1,000 m long and 200 m wide, which align parallel to the general direction of the fault systems associated with the DS Rift. The abrupt appearance of the sinkholes reflects changes to the groundwater regime around the shrinking DS. The eastward retreat of the shoreline and the lake-level drop (1 m/year in recent years) cause an eastward and downward migration of the fresh/saline groundwater interface. Consequently, a subsurface salt layer, which was previously enveloped by saline groundwater, is gradually being invaded and submerged by relatively fresh groundwater, and cavities form due to the rapid dissolution of the salt. Collapse of the overlying sediments into these cavities results in sinkholes at the surface. An association between sinkhole sites and land subsidence is revealed by interferometric synthetic aperture radar (InSAR) measurements. On a broad scale (hundreds of meters), subsidence occurs due to compaction of fine-grained sediments as groundwater levels decline along the retreating DS shoreline. At smaller scales (tens of meters), subsidence appears above subsurface cavities in association with the sinkholes, serving in many cases as sinkhole precursors, a few weeks to more than a year before their actual appearance at the surface. This paper overviews the processes of sinkhole formation and their relation to land subsidence.     

Three-dimensional numerical modeling of land subsidence in Shanghai,China

Shanghai, in China, has experienced two periods of rapid land subsidence mainly caused by groundwater exploitation related to economic and population growth. The first period occurred during 1956–1965 and was characterized by an average land subsidence rate of 83 mm/yr, and the second period occurred during 1990–1998 with an average subsidence rate of 16 mm/yr. Owing to the establishment of monitoring networks for groundwater levels and land subsidence, a valuable dataset has been collected since the 1960s and used to develop regional land subsidence models applied to manage groundwater resources and mitigate land subsidence. The previous geomechanical modeling approaches to simulate land subsidence were based on one-dimensional (1D) vertical stress and deformation. In this study, a numerical model of land subsidence is developed to simulate explicitly coupled three-dimensional (3D) groundwater flow and 3D aquifer-system displacements in downtown Shanghai from 30 December 1979 to 30 December 1995. The model is calibrated using piezometric, geodetic-leveling, and borehole extensometer measurements made during the 16-year simulation period. The 3D model satisfactorily reproduces the measured piezometric and deformation observations. For the first time, the capability exists to provide some preliminary estimations on the horizontal displacement field associated with the well-known land subsidence in Shanghai and for which no measurements are available. The simulated horizontal displacements peak at 11 mm, i.e. less than 10 % of the simulated maximum land subsidence, and seems too small to seriously damage infrastructure such as the subways (metro lines) in the center area of Shanghai.     

The effect of stratigraphic heterogeneity on areal distribution of land subsidence at Taiyuan, northern China

Taiyuan city has experienced serious land subsidence since the 1950s, with the maximum accumulated subsidence of 2,960 mm and a total affected area of 585 km² by 2003. Tectonic settlement was found to contribute about 1% and the major cause is over-exploitation of groundwater. The spatial–temporal relationship between the areal distribution of land subsidence and that of groundwater drawdown from 1956 to 2000 indicates that although land subsidence centers overall match groundwater depression cones, there are local deviations, and that although the time series curves are similar, land subsidence always lags behind groundwater level decline, with different lag time at different sites. The major findings of this work on the control of stratigraphic heterogeneity on the spatial pattern of land subsidence at Taiyuan include: (1) land subsidence centers shift from the corresponding groundwater depression cones to the sides with thicker accumulated clay layers; (2) under the same pumping rate, land subsidence at places with more clay interlayers and thinner individual interlayers is greater and the lag time shorter; and (3) land subsidence is closely related to the physical properties of clay soils. The Interbed Storage Package-1, a modular subroutine of MODFLOW was employed to simulate the areal distribution of individual layer compression. The modeling results show that compression of different clay layers has different contribution to land subsidence. Pumping groundwater from water-bearing zones close to the most compressible clay layers should therefore be carefully controlled.     

Numerical simulation of land subsidence induced by groundwater overexploitation in Su-Xi-Chang area,China

Su-Xi-Chang area is one of the typical regions in China which suffers from severe land subsidence. Various field monitoring records were integrated to study the characteristics and mechanisms of land subsidence in this region. The development of the land subsidence in this region shows a tight spatial and temporal correlation with the groundwater pumping. Based on the analysis of the field data, it is found that the deformation patterns of the hydrogeologic units are greatly related to the hydrogeologic properties and groundwater level variations. Some have an elastic behavior, others may have an elastic–plastic rheology. Hence, a 3D finite element numerical model considering the rheological properties of the soil was developed to simulate the groundwater level and land subsidence. Both hydraulic conductivity and specific storage were expected to vary with the porosity during the process of consolidation. Multiscale finite element method (MsFEM) was applied to solve the model during the period from 1996 to 2004. After calibrating the model with the observed groundwater level and subsidence data, the parameters of the multi-layers system were estimated. The calibrated model outputs fit reasonably well with the observed data. Consequently the model can be applied to predict groundwater level and land subsidence in future pumping scenarios. The model predictive results show that land subsidence rate can be controlled and even rebound may occur after the implementation of the groundwater exploitation prohibition. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Land subsidence and earth fissures in south-central and southern Arizona,USA

Land subsidence due to groundwater overdraft has been an ongoing problem in south-central and southern Arizona (USA) since the 1940s. The first earth fissure attributed to excessive groundwater withdrawal was discovered in the early 1950s near Picacho. In some areas of the state, groundwater-level declines of more than 150 m have resulted in extensive land subsidence and earth fissuring. Land subsidence in excess of 5.7 m has been documented in both western metropolitan Phoenix and Eloy. The Arizona Department of Water Resources (ADWR) has been monitoring land subsidence since 2002 using interferometric synthetic aperture radar (InSAR) and since 1998 using a global navigation satellite system (GNSS). The ADWR InSAR program has identified more than 25 individual land subsidence features that cover an area of more than 7,300 km². Using InSAR data in conjunction with groundwater-level datasets, ADWR is able to monitor land subsidence areas as well as identify areas that may require additional monitoring. One area of particular concern is the Willcox groundwater basin in southeastern Arizona, which is the focus of this paper. The area is experiencing rapid groundwater declines, as much as 32.1 m during 2005–2014 (the largest land subsidence rate in Arizona State—up to 12 cm/year), and a large number of earth fissures. The declining groundwater levels in Arizona are a challenge for both future groundwater availability and mitigating land subsidence associated with these declines. ADWR’s InSAR program will continue to be a critical tool for monitoring land subsidence due to excessive groundwater withdrawal.     

Hydrogeochemical and isotopic constraints on the origins of dryland salinity,Murray Basin,Victoria, Australia

Combined hydrogeological and hydrogeochemical data allow flow systems and the origins of solutes in the Honeysuckle Creek area of the southeastern Murray Basin, which is an area affected by dryland salinity, to be constrained. Recharge occurs both on the uplands that are composed of fractured Violet Town Volcanic rocks and the Riverine Plain that comprises sediments of the Shepparton and Coonambidgal Formations. Groundwater from the Violet Town Volcanics has low salinity (<20 mmol/L Cl) and major ion geochemistry that is controlled largely by dissolution of silicate minerals. Low Cl/Br ratios (as low as 281 molar) suggest that this groundwater has not dissolved halite. Groundwater that recharged through the Riverine Plain sediments has higher Cl/Br ratios (up to 1146) and Cl concentrations of <20 mmol/L, consistent with it dissolving minor halite. Higher salinity (>20 mmol/L) groundwater has intermediate Cl/Br ratios (600–1000), which indicate that the high salinities do not simply result from halite dissolution. Rather, mixing of groundwater homogenises Cl/Br ratios, and evaporation as a consequence of a shallow water table is the dominant process that increases salinity. Oxygen and H isotopes also indicate that mixing and evaporation have occurred. These results indicate that land use over the whole region, not just the uplands, needs to be considered in any salinity management plans. Additionally future development of salinity is controlled by depth to the water table on the plains and the efficiency of recharge rather than by salt stores (halite or brines) in the unsaturated zone.     

Groundwater flow and its effect on salt dissolution in Gypsum Canyon watershed,Paradox Basin,southeast Utah,USA

Groundwater flow is an important control on subsurface evaporite (salt) dissolution. Salt dissolution can drive faulting and associated subsidence on the land surface and increase salinity in groundwater. This study aims to understand the groundwater flow system of Gypsum Canyon watershed in the Paradox Basin, Utah, USA, and whether or not groundwater-driven dissolution affects surface deformation. The work characterizes the groundwater flow and solute transport systems of the watershed using a three-dimensional (3D) finite element flow and transport model, SUTRA. Spring samples were analyzed for stable isotopes of water and total dissolved solids. Spring water and hydraulic conductivity data provide constraints for model parameters. Model results indicate that regional groundwater flow is to the northwest towards the Colorado River, and shallow flow systems are influenced by topography. The low permeability obtained from laboratory tests is inconsistent with field observed discharges, supporting the notion that fracture permeability plays a significant role in controlling groundwater flow. Model output implies that groundwater-driven dissolution is small on average, and cannot account for volume changes in the evaporite deposits that could cause surface deformation, but it is speculated that dissolution may be highly localized and/or weaken evaporite deposits, and could lead to surface deformation over time.     

基于地面沉降控制的区域性松散沉积层地下水可采资源规划评价

为了满足区域性松散沉积层地区地面沉降的防控要求,规划评价地下水的可采资源量,根据渗流理论和土力学理论,建立了地下水三维非稳定渗流与地面沉降耦合数学模型,并考虑了含水层孔隙度、渗透系数、储水率随含水层发生固结沉降的变化特征。采用三维有限元数值分析方法,以江苏省南通市地下水开采为例,基于地面沉降的控制要求,规划评价出了各乡镇各含水层的地下水可采资源量。结果表明：地下水开采布局科学规划后,总的可采资源量为17 870.56×10⁴ m³/a,较现状开采量10 902.32×10⁴ m³/a有较大幅度的增加。地下水三维变参数非稳定流与地面沉降耦合模型可以更加精确地刻画三维水文地质体的特征,更加符合实际情况。     

Numerical Simulation of Groundwater Flow and Vulnerability in Wadi El-Natrun Depression and Vicinities,West Nile Delta,Egypt

During the last 25 years, rapid and unplanned land reclamation activity has been carried out in the areas located in both south and east of Wadi El - Natrun Depression of Egypt. Accordingly, negative effects on groundwater levels and vulnerability are frequently caused by localized high levels of abstraction and the return-flow of polluted irrigation water respectively. A groundwater model is a computational method that presents an approximation of an underground water system. In this study the groundwater system is simulated both in quantity and quality by using Mass Balance Transfer Model (NETPATH), Groundwater Modeling System (GMS) and DRASTIC Model to investigate the water - rock interactions, groundwater levels drawdown and vulnerability respectively. Three main geochemical processes namely dedolomitisation, dissolution of halite and silicate weathering were estimated during the flow path. The present over-abstraction of groundwater (105.84 million m³/year) has induced a general head drawdown from 3 to 40 m in years 2015 and 2050 respectively. Best estimate using a 3D GMS hydraulic model was (157000 m³/day) a strategy proposed for the management of groundwater without critical depletion (second scenario). The results document the extent to which a high drawdown can greatly reach 4 m by the end of simulation year 2050. The vulnerability maps of groundwater were constructed using the DRASTIC index method. The results indicated that, the southeastern and central portions of the study area are having high vulnerability rate (> 110). Modified DRASTIC map showed many more dominant high risk areas in the eastern parts of the study area that were low risk, which may be attributed to return flow of polluted irrigation water.     

On causes and impacts of land subsidence in Bandung Basin, Indonesia

The Bandung Basin is a large intra-montane basin surrounded by volcanic highlands, in western Java, Indonesia, inhabited by more than seven million people. The basin, an area of about 2,300 km², is a highland plateau at approximately 650–700 m above sea level and is surrounded by up to 2,400 m high Late Tertiary and Quaternary volcanic terrain. Based on the results of nine GPS surveys conducted since 2000 up to 2011, it was shown that several locations in the Bandung Basin have experienced land subsidence, with an average rate of about ?8 cm/year and can go up to about ?23 cm/year in certain locations. A hypothesis has been proposed by several studies that land subsidence observed in several locations in the Bandung Basin has been caused mainly by excessive groundwater extraction. It is found that there is a strong correlation between the rates of groundwater level lowering with the GPS-derived rates of land subsidence in several locations in Bandung Basin. The GPS results in this study detected significant subsidence in the textile industry area, where very large volumes of groundwater are usually extracted. The impact of land subsidence in Bandung can be seen in several forms, mainly in the cracking and damage of houses, buildings and infrastructure. Land subsidence also aggravates the flooding in Bandung Basin, which has brought huge economic losses and deteriorated the quality of life and environment in the affected areas.     

Threat of land subsidence in and around Kolkata City and East Kolkata Wetlands,West Bengal,India

This paper attempts to estimate the possible rate of land subsidence of Kolkata City including Salt Lake City and the adjoining East Kolkata Wetlands located at the lower part of the deltaic alluvial plain of South Bengal basin. Demand of groundwater for drinking, agricultural and industrial purposes has increased due to rapid urbanization. The subsurface geology consists of Quaternary sediments comprising a succession of clay, silty clay and sand of various grades. Groundwater occurs mostly under confined condition except in those places where the top aquitard has been obliterated due to the scouring action of past channels. Currently, the piezometric head shows a falling trend and it may be accelerated due to further over-withdrawal of groundwater resulting in land subsidence. The estimated mean land subsidence rate is 13.53 mm/year and for 1 m drop in the piezometric head, the mean subsidence is 3.28 cm. The surface expression of the estimated land subsidence is however, cryptic because of a time lag between the settlement of the thick low-permeable aquitard at the top and its surface expression. Therefore, groundwater of the cities and wetland areas should be developed cautiously based on the groundwater potential to minimize the threat of land subsidence.     

弹塑性变形条件下抽水引起的地面沉降三维数值模拟

基于部分耦合原理,采用TOUGH2和FLAC3D建立抽水引起的三维地面沉降弹塑性模型,模型中综合考虑土体的弹塑性变形特征、渗流-应力的双向耦合作用以及参数的非线性,探讨了持续抽水和脉冲抽水两种抽水过程中地面沉降发展演化过程。研究结果表明:（1）集中抽水停止后地面沉降会发生回弹,抽水中心沉降量不断减小。由于水平方向存在水力梯度,地下水继续向地下水位漏斗中心渗流从而导致沉降漏斗的范围仍继续扩大;（2）脉冲抽水导致土体的孔隙水压力、渗透系数以及沉降量均呈周期性波动变化,地面沉降会局部回弹,但总体仍随着抽水的持续,沉降量不断增加;（3）在抽水量相同前提下,对比持续抽水与脉冲抽水两种方式引发的塑形沉降量可知,抽水速率小、脉冲式多次开采导致的塑性沉降量较小,持续抽水的抽水速率越小、脉冲抽水间隔越短越有利于控制地面沉降。研究成果为地面沉降数值模拟提供了一种新方法,其中算例研究能为抽水条件下地面沉降的控制提供参考。     

Groundwater resources management under environmental constraints in Shiroishi of Saga plain, Japan

A sinking of the land surface due to the pumping of groundwater has long been recognized as an environmental issue in the Shiroishi plain of Saga, Japan. Land subsidence can have several negative economic and social implications such as changes in groundwater and surface water flow patterns, restrictions on pumping in land subsidence prone areas, localized flooding, failure of well casings as well as shearing of structures. To minimize such an environmental effect, groundwater management should be considered in this area. In this study, a new integrated numerical model that integrates a three-dimensional numerical groundwater flow model coupled with a one-dimensional soil consolidation model and a groundwater optimization model was developed to simulate groundwater movement, to predict ground settlement and to search for optimal safe yield of groundwater without violating physical, environmental and socio-economic constraints. It is found that groundwater levels in the aquifers greatly vary from season to season in response to the varying climatic and pumping conditions. Consequently, land subsidence has occurred rapidly throughout the area with the Shiroishi plain being the most prone. The predicted optimal safe yield of the pumping amount is about 5 million m³. The study also suggests that pumping with this optimal amount will minimize the rate of land subsidence over the entire area. An erratum to this article can be found at     

华北平原典型地段地面沉降演化特征与机理研究

华北平原是世界上最大的地下水系统之一,地面沉降问题突出。由于沉积环境变化多样、地质条件差异性显著和人类开采活动强烈,使得该地区地面沉降成因机理复杂。本文采用卫星对地观测技术与传统手段相结合,监测地面沉降灾变过程,分析华北平原地面沉降发展历史和现状。结合应力-应变图解法及土工实验研究地面沉降差异性特征及滞后变形成因机理。取得了以下关键认识:(1)华北平原地面沉降空间分布差异性明显,沉降主要分布在平原区第四纪沉积凹陷,呈现东西分带、南北分段特点。地面沉降空间发展部分受到北东向和北西向构造控制。在沧县隆起区,地面沉降也比较发育,主要原因是沧县隆起在第四纪时期构造运动相对不活跃,沉积了较厚的第四系;存在与构造走向一致的3期古河道,该地区赋存丰富的地下水资源并被大量开采。(2)地面沉降发生发展与地下水开采历史密切相关,沉降主要压缩贡献层随地下水开采层位变化而变化。北京平原100 m以深地层对地面沉降贡献呈增加趋势。天津平原目前地面沉降的主要贡献层来自300 m以下地层。(3)气候干旱导致地下水补给量减少,同时增加了地下水的开采,因而是引起地面沉降的重要间接驱动因素。高层建筑荷载、基坑降排水、地热开采对地面沉降的影响应引起足够重视。(4)地面沉降具有很强的滞后性,最大滞后时间可达25年。除了渗透固结成因以外,土体蠕变是另外一个重要原因。更新世地层在不同荷载下,蠕变特征明显。沧县隆起晚更新世地层次固结可达到总变形28.3%。(5)土的物理性质、地下水位变化模式对土层变形特征具有重要影响。不同埋深地层在地下水位变化条件下的变形特征存在较大的差异(弹性、黏弹性、黏弹塑性)。浅部含水组呈现以弹性为主的变形特征。     

The development and control of the land subsidence in the Yangtze Delta,China

Land subsidence in China occurs predominantly in 17 provinces (cities) situated in the eastern and middle regions of the country, including Shanghai, Tianjin and Jiangsu, and Hebei provinces. It is primarily caused by groundwater overpumping. One of the areas most severely affected by land subsidence is the Yangtze Delta, most of which consists of Shanghai City, the Su-Xi-Chang area (Suzhou, Wuxi and Changzhou cities) of Jiangsu Province, and the Hang-Jia-Hu area (Hangzhou, Jiaxing and Huzhou cities) of Zhejiang Province. The excessive exploitation of groundwater forms in a large regional cone of depression and, consequently, land subsidence is also regional, currently centered in the Shanghai and Su-Xi-Chang areas. In 2002, the maximum cumulative subsidence of Shanghai, Su-Xi-Chang and Hang-Jia-Hu were 2.63 m, 2.00 and 1.06 m, respectively. The land subsidence area is continuing to expand throughout the Yangtze Delta. To study the characteristics and the pattern of this land subsidence, the government has implemented a monitoring system involving the placement of 37 groups of extensometers (layers marks) and drilling of more than 1000 observation wells. These provide an invaluable historical record of deformation and pore water pressure and facilitate studies on the special features of soil deformation when the groundwater level changes due to pumping. Several measures have been taken in recent years to control the development of the land subsidence in the different areas; these include groundwater injection, prohibition of pumping deep confined groundwater, and an adjustment of the pumping depth and magnitude of the groundwater withdrawn. At present, although the subsidence area is still increasing slowly, the subsidence rate is controlled.     

黄河三角洲地下水动态变化及其与地面沉降的关系

为了分析黄河三角洲地下水动态及其与地面沉降的关系, 利用多年地下水和地面沉降监测数据, 发现黄河三角洲广饶县和东营区的地下水动态变化剧烈且地面沉降严重, 含水层多处于超采状态, 浅、深层地下水降落漏斗先后出现.深层地下水降落漏斗中心水位下降速度达2~3m/a.近年来, 东营和广饶地面沉降漏斗中心沉降量和速率分别为155.1mm、28.2mm/a和356.0mm、64.7mm/a.借助GIS技术及数理统计法, 发现深层地下水降落漏斗与沉降漏斗空间耦合良好, 深层地下水位与地面高程呈线性正相关, 相关系数为0.92, 深层地下水过度开采已成为影响沉降的最根本因素.井灌区第三粘性压缩层成为地面沉降主要贡献层, 且深层地下水降落漏斗中心的地下水位已低于第三承压含水层临界水位, 沉降趋于严重.      

冲洪积平原地面沉降特征及主控因素——以北京平原为例

北京由于长期过量开采地下水,相继引发了一系列地质环境问题,其中地面沉降问题尤为突出。回顾了北京地面沉降发展历史,从平面和垂向上分析了地面沉降特征,在此基础上对北京冲洪积平原区沉降的主控因素进行了研究。结果表明:(1)平面上,沉降分为南、北2个大区,7个沉降中心。北区已由多个单独沉降中心区扩展成一个大区域,南区北扩明显;(2)垂向上,南区第一压缩层为沉降主贡献层,沉降占比42%,浅部地层沉降速率减小,深部地层沉降速率增加。土体变形特征为塑性变形,包含蠕变变形;北区第二压缩层为沉降主贡献层,沉降占比65%,浅部沉降量值很小且波动平缓,深部沉降量相对较大。土体变形特征为浅部以弹性变形为主,深部以塑性变形为主,包含蠕变变形;(3)沉降受构造作用及基底格架控制,北东方向受冲洪积扇上部单一砂卵砾石的地层条件控制扩展范围有限,沉降整体向北西、南东方向扩张;(4)地层结构决定沉降平面和垂向分布特征,尤其北部冲洪积与南部湖相沉积的差异,是产生深浅部地层沉降贡献率不同的重要因素;(5)地下水开采仍是沉降产生的主因,地下水漏斗的扩展和沉降中心的分布高度吻合,主要沉降层地下水位下降速率与沉降速率成正比。     

苏锡常地区地面沉降地质结构三维可视化模型虚拟现实系统研究

从苏锡常地区地面沉降研究实际需要出发,介绍了利用虚拟现实技术构建地面沉降地质结构三维可视化模型虚拟现实系统的设计思路、开发模式、实现方法和功能特点。该系统真实表达了研究区三维地层结构特征,逼真再现了地面沉降模型所预测的地下水流场、地面沉降发展变化过程,并虚拟表现出地面沉降可能造成的后果,为苏锡常地区地面沉降机理及预测研究提供了一个全新的三维可视化平台。     

三维区域地面沉降数值模拟

由于地裂缝研究及地表变形监测技术（例如GPS,InSAR等）的快速发展,抽取地下水引起的地层水平位移受到关注。传统区域地面沉降模型虽然求解快速但不能模拟水平位移;比奥模型虽然能够模拟土体的三维变形,但模型求解的计算量较大,较少应用于大尺度的区域地面沉降数值模拟。为解决以上问题,推导了解耦三维地面沉降数学模型,模型推导过程显示：比奥模型假设法向总应力和不变,则可简化为解耦三维地面沉降模型;解耦三维地面沉降模型假设土体仅有垂向一维变形,则可简化为传统区域地面沉降模型。同时通过数值试验验证了解耦三维地面沉降模型可以作为比奥模型的替代模型和传统区域地面沉降模型的改进模型,用来模拟抽取地下水引起的三维区域地面沉降。     

唐山沿海地区地面沉降渗流固结耦合模拟研究

唐山沿海地区经济在迅速发展,沿海地区城市化规模在扩大,地下水开采量增大,地面沉降加剧.文中分析了唐山沿海地区的水文地质条件,概化为3个含水层、3个弱透水层,共6个压缩层.建立了三维地下水流和垂向一维压缩完全耦合模型.采用25a的观测资料校正模型,计算值与实测值拟合较好,模型具有较高的仿真性和适用性.预测了10a末的地面沉降;当地下水以现有开采量开采时,沉降中心累计达1192.3mm, 10a沉降352.3mm,沉降速率为35.23mma-1;当地下水的开采量在现有开采量的基础上增加10%时,沉降中心累计达1260.8mm, 10a沉降420.8mm,沉降速率为42.08mma-1; 当地下水的开采量在现有开采量的基础上减小10%时,沉降中心累计达1088.7mm, 10a沉降247.9mm,沉降速率为24.79mma-1.增大10%的地下水开采量, 10a地面沉降量增加68.5mm;减少10%的地下水开采量, 10a地面沉降量减少104.4mm.因此,控制地下水开采量是控制地面沉降的有效方法.