Validation and comparison of different techniques for the derivation of digital elevation models and volcanic monitoring (Vulcano Island,Italy)

High accurate digital elevation models (DEM) acquired periodically over a volcanic area can be used for monitoring crustal deformations. Airborne stereoscopic photography is a powerful tool for the derivation of high resolution DEM, especially when combined with Global Positioning System (GPS). We analyse data acquired on Vulcano Island (Italy) to assess the performance of two photogrammetry methods for DEM generation. The first method is based on automatic digital processing of scanned airborne stereo images from a film camera (Wild RC20). In the second method digital stereo data from the multi-spectral High Resolution Stereo Camera-Airborne (HRSC-A) are used. Accuracy assessment through comparison with kinematic GPS height profiles shows that both DEMs have accuracy on the order of few decimetres. Direct comparison of the two DEMs on the La Fossa volcanic cone provides a standard deviation of the residuals of 78 cm. Residuals greater than two metres between the two DEMs acquired at one year interval are locally evidenced in unstable areas with uneven morphology. The application of photogrammetric DEMs is also discussed within a SAR interferometry study carried out on Vulcano Island to evaluate the potentialities of such techniques for ground deformation monitoring. Although accuracy better than 1 m or 2 m is not required for satellite SAR interferometry, we show how the precise photogrammetric DEMs could still significantly improve SAR interferograms of Vulcano Island.     

Application of the SBAS-DInSAR technique to fault creep: A case study of the Hayward fault, California

We present a quantitative assessment of the capability of the differential SAR interferometry (DInSAR) technique referred to as Small BAseline Subset (SBAS) approach to investigate fault creep phenomena. In particular we have computed, via the SBAS-DInSAR algorithm, time series of the surface displacements relevant to the Hayward fault zone, within the San Francisco Bay Area (California), from the European Space Agency's ERS-1/2 satellite radar data for the 1992 to 2000 time period. Starting from the DInSAR time series we measured the relative displacements across the fault with no need for any atmospheric filtering step. These results have been systematically compared to the measurements available from the alignment arrays that are located along the fault. Our analysis shows that the standard deviation of the differences between the DInSAR and the in situ measurements is on the order of 2 mm. Moreover, the estimated mean deformation rates have an accuracy that is better than 1 mm/year.     

Two‐scale surface deformation analysis using the SBAS‐DInSAR technique: a case study of the city of Rome,Italy

We have exploited the capability of the differential synthetic aperture radar (SAR) interferometry (DInSAR) technique, referred to as Small BAseline Subset (SBAS) approach, to analyse surface deformation at two distinct spatial scales: a low resolution, large scale, and a fine resolution, local scale. At the large scale, the technique investigates DInSAR data with a ground resolution of the order of 100 m×100 m and leads to generate mean deformation velocity maps and associated time series for areas extending to some thousands of square kilometres. At the local scale, the technique exploits the SAR images at full spatial resolution (typically of the order of 5 m×20 m), detecting and analysing localized deformation phenomena. The study is focused on the city of Rome, Italy, and we used the ERS‐1/2 satellite radar data relevant to the 1995–2000 time period. The presented results demonstrate the capability of the SBAS approach to retrieve, from the low‐resolution DInSAR data, large‐scale deformation information leading to identify several sites affected by significant displacements. Our analysis permitted us to conclude that a major contribution to the detected displacements is due to the consolidation of the alluvial soils present in the area, mostly enforced by the buildings' overload. Furthermore, in a selected area, a detailed analysis was carried out by exploiting the full resolution DInSAR data. In this case we investigated deformation phenomena at the scale of single buildings. As key result we showed that differential displacements of few mm a^?1, affecting single man‐made structures or building complexes, could be detected, thus allowing to identify sites that may potentially be involved in critical situations.     

Urban subsidence in the city of Prato (Italy) monitored by satellite radar interferometry

Differential Synthetic Aperture Radar (SAR) interferometry has been widely used to monitor centimetric surface deformations in geophysical applications. In this letter, this technique is applied to study the displacement field in an urban area. A set of six European Remote Sensing satellite (ERS)-1/2 SAR images has been used to detect, map and quantify the subsidence occurring in the city of Prato near Florence (Italy). Four areas which have been affected by strong subsidence during the period 1993-2000 have been spotted within the city. The analysis of three interferograms processed from images taken two years apart shows that the deformation rate appears to be relatively constant, with a maximum value of about 8.3 cm y m 1 .     

Monitoring dewatering induced subsidence and fault reactivation using interferometric synthetic aperture radar

Using multitemporal differential interferometric synthetic aperture radar analysis integrated with pumping and site geologic data we present evidence for hydrologically induced large subsidence in and around an ongoing open‐pit mine with intensive dewatering operations. Analysis of numerous differential synthetic aperture radar interferometry (DInSAR) pairs spanning the period 1993 to 2001 reveals the abrupt appearance of these features to intervals of a few to several months. Along a section through the anomaly, we plotted dewatering associated changes in the groundwater levels at monitoring wells. We also used DInSAR to extract several individual kilometre‐lengths, centimetre amplitude normal fault reactivation events in the alluvial sediments adjacent to the mine dewatering operation. High‐resolution remote sensing analyses provide strong evidence that these features align with faults active in the last several thousand years. We interpret these reactivations as mechanically involving only the upper few hundred metres of the existing fault plane above the alluvial aquifer affected by the mine dewatering.     

The present status of subsiding land vulnerable to roof collapse in the Jharia Coalfield,India, as obtained from shorter temporal baseline C-band DInSAR by smaller spatial subset unwrapped phase profiling

In this paper, we identified recently subsiding areas in Jharia Coalfield, Jharkhand, India from the shorter temporal baseline Radarsat-2 C-band interferometric synthetic aperture radar (InSAR) data pairs of 2012. Although shorter wavelength C-band differential InSAR (DInSAR) is more sensitive to slow deformation and better suited for higher precision land subsidence measurement, the dynamic and adverse land cover in mining areas and resulting temporal decorrelation problem poses a serious problem for DInSAR observation in mining areas. We used smaller temporal baseline data pairs and adopted InSAR coherence-guided incremental filtering with smaller moving windows to highlight the deformation fringes over temporal decorrelation noise. We identified the deformation fringes and validated them based on ground information to prepare the land subsidence map of the coalfield in 2012. Several new, previously unreported subsidence areas were detected in the present study with a total subsiding area of 6.9 km². The recent incidence of roof collapse on 15 November 2014 at Angar Patra village in Katras region of the coalfield where 45 houses collapsed and 10 people were injured is situated in a highly subsiding vulnerable area as obtained from the present study. Due to spatial discontinuities of InSAR coherence, DInSAR phase unwrapping for the entire study area in one go did not appear feasible. To avoid this problem, we performed DInSAR processing in smaller spatial subsets and unwrapping of the subset interferograms by a ‘minimum cost flow’ algorithm. Subsequently, we plotted unwrapped phase profiles across the deformation fringes and retrieved the maximum deformation phase with respect to background phase and translated them into radar line of sight (LOS) displacement rates. For obtaining the average subsidence rates, we adopted InSAR coherence-weighted LOS displacement rates taking into account the contribution of each data pair as a function of DInSAR phase quality of the fringe areas. Ground-based subsidence measurements by precision levelling were conducted in four test sites that had been undergoing active underground mining during the observation period. We compared space-borne DInSAR-based subsidence rates obtained by the adopted technique with precision levelling measurements. Overall, the results are found to agree well. In the four test sites with gentle to flat topography, land subsidence occurs at slow to moderate rates due to compression of in-filled material (resulting from sand stowing in underground mining), without any evidence of roof collapse. In such cases, the horizontal displacement component is less significant, and overall surface displacement occurs essentially in the vertical direction. However, we assessed the nature of subtle horizontal strain to infer relative shrinkage or dilation of the land surface which could be additive or subtractive to vertical displacement in DInSAR-based LOS displacement.     

Insight into ground deformations at Lascar volcano (Chile) from SAR interferometry, photogrammetry and GPS data: Implications on volcano dynamics and future space monitoring

We present a detailed study of Lascar volcano (Chile) based on the combination of satellite, aerial and ground-based data, in order (i) to better characterize the deformation style of Andean explosive volcanoes, and (ii) to provide new insights on the potential of space techniques to monitor active volcanic deformations on such edifices. Lascar is one of the most active volcanoes in Central Andes characterized by a recent cyclic activity. Additionally, it is located in favourable conditions for radar imaging. Lascar thus offers very good conditions for studying large to small scale ground deformations associated with volcano dynamics. The analysis of InSAR (Synthetic Aperture Radar interferometry) time series data from the European and Japanese satellites (ERS, JERS) acquired between 1993 and 2000, encompassing three eruptive events, confirmed the absence of broad far-field deformation signal. Thus during the recent activity of Lascar we discard significant magmatic input at depth. The following approaches were used to improve the InSAR signal / noise ratio in order to detect possible local deformation. We carried out a quantitative evaluation of the potential tropospheric contribution in INSAR interferograms for the Salar de Atacama-Lascar area using radar (ASAR-ENVISAT) and spectrometer (MODIS) data. We also used an accurate aerial photogrammetric and GPS constrained DEM in our InSAR data reprocessing. We find a co-eruptive ground-deformation confined into the summit crater for the 1995 eruption. This deformation has spatial dimension of 500 by 400 m and relates to a subsidence of crater floor up to 17 mm. We interpret it as pressure or volume decrease at subsurface levels below the active crater. Our study made it possible to image a new near-field volcanic deformation confined within the summit crater of the Lascar volcano. It also demonstrates that the combination of precise photogrammetry DEM and INSAR data can significantly improve our ability to remotely sense subtle surface deformation on these explosive volcanoes. This methodology might contribute to better understand volcano dynamics and to complement their monitoring in remote areas.     

Monitoring of slow ground deformation by ERS radar interferometry on the Vauvert salt mine (France): Comparison with ground-based measurement

The differential SAR Interferometry (DInSAR) technique has been applied to a test site near Vauvert (France) to detect and monitor ground deformation. This site corresponds to the location of an industrial exploitation of underground salt using the solution mining technique. An area of subsidence has been observed using in situ measurements. Despite conditions unfavorable for InSAR because of the vegetal cover, we show that radar remote sensing observations provide valuable information which substantially improves our knowledge of the phenomenon. An adaptive phase filtering process has been used to improve the coherence level. In particular, our study shows that the geometry of the subsidence bowl is different to that previously assumed using ground-based techniques only. The size of the subsidence bowl (8 km) is larger than expected. This information will be useful for further modeling of the deformation and to improve the coverage of the in situ measurement networks. It also shows that radar interferometry can be used for the long-term monitoring of such sites and to predict potential environmental issues.     

Surface movements in Bologna (Po Plain — Italy) detected by multitemporal DInSAR

We studied the surface deformations affecting the southeastern sector of the Po Plain sedimentary basin, in particular the area of Bologna. To this aim an advanced DInSAR technique, referred to as DInSAR-SBAS (Small BAseline Subset), has been applied. This technique allows monitoring the temporal evolution of a deformation phenomenon, via the generation of mean deformation velocity maps and displacement time series from a data set of acquired SAR images. In particular, we have processed a set of SAR data acquired by the European Remote Sensing Satellite (ERS) sensors and compared the achieved results with optical levelling measurements, assumed as reference. The surface displacements detected by DInSAR SBAS from 1992 to 2000 are between 10 mm/year in the historical part of Bologna town, and up to 59 mm/year in the NE industrial and agricultural areas. Former measurements from optical levelling referred to 1897 show 2-3 mm/year vertical movements. This trend of displacement increased in the second half of the 20th century and the subsidence rate reached 60 mm/year. We compared the more recent levelling campaigns (in 1992 and late 1999) and DInSAR results from 1992 to 1999. The standard deviation of the difference between levelling data, projected onto the satellite Line Of Sight, and DInSAR results is 2 mm/year. This highlights a good agreement between the measurements provided by two different techniques. The explanation of soil movements based on interferometric results, ground data and geological observations, allowed confirming the anthropogenic cause (surface effect due to the overexploitation of the aquifers) and highlights a natural, tectonic, subsidence.     

Surface deformation from persistent scatterers SAR interferometry and fusion with leveling data: A case study over the Choushui River Alluvial Fan, Taiwan

The Choushui River Alluvial Fan (CRAF), one of the most important agricultural areas in Taiwan, suffers severe subsidence caused by groundwater withdrawal. We use 20 ENVISAT images from 2006 to 2008 to derive vertical deformation over CRAF by persistent scatterer interferometry (PSI), which effectively reduces errors affecting conventional DInSAR techniques. The pixel density over CRAF is 107.6 pixels/km², compared to 0.19 points/km² of the leveling benchmarks. PSI yields vertical displacements matching the leveling result to 0.6 cm/year (RMS), and provides a higher spatial resolution of subsidence than the latter. We also develop a data fusion method that considers the high-precision, low point-density leveling result as a smoothed correction to the PSI result. The combined field is more representative of overall deformation characteristics than the PSI-only or leveling-only field, and it is better suited for the assessment of the impact of subsidence over CRAF.     

Detecting changes in surface soil moisture content using differential SAR interferometry

The feasibility of measuring changes in surface soil moisture content with differential interferometric synthetic aperture radar (DInSAR) has received little attention in comparison with other active microwave techniques. In this study, multi-polarization C- and L-band DInSAR is explored as a potential tool for the measurement of changes in surface soil moisture in agricultural areas. Using 10 ascending phased array L-band SAR (PALSAR) scenes acquired by the Japanese Advanced Land Observing Satellite (ALOS) and 12 descending advanced SAR (ASAR) scenes acquired by the European ENVISAT satellite between July 2007 and November 2009, a series of 27 differential interferograms covering a common study area over southern Ireland were generated to investigate whether small-scale changes in phase are linked to measured soil moisture changes. Comparisons of observed mean surface displacement and in situ mean soil moisture change show that C-band cross-polarization pairs displayed the highest correlation coefficients over both the barley (correlation coefficient, r = 0.51, p = 0.04)- and potato crop (r = 0.81, p = 0.003)-covered fields. Current results support the hypothesis that a soil moisture phase contribution exists within differential interferograms covering agricultural areas.     

Ground deformation due to tectonic, hydrothermal, gravity, hydrogeological, and anthropic processes in the Campania Region (Southern Italy) from Permanent Scatterers Synthetic Aperture Radar Interferometry

We apply the Permanent Scatterers Synthetic Aperture Radar Interferometry (PS-InSAR) technique to the Campania Region (Southern Italy), which includes the Southern Apennines chain and Plio-Quaternary structural depressions, with the aim to detect ground displacements at a regional scale. The study area, which extends for 13,600 km², is characterized by intense urbanization, active volcanoes (Phlegraean Fields, Vesuvius and Ischia), seismogenic structures, landslides, hydrogeological instability. PS-InSAR technique allows us to identify a set of radar benchmarks (PS) where accurate displacement measurements can be carried out. About 1.7·10⁶ PS are identified by processing Synthetic Aperture Radar (SAR) images acquired in ascending and descending orbits from 1992 to 2001 by the European Remote Sensing satellites (ERS). The PS-InSAR application at regional scale detected ground deformations ranging from + 28 to − 39 mm/yr. The calculated velocity values are consistent with the available GPS and levelling data from selected areas. We identify volcanic areas in which the deformation is mainly related to the depressurization of the local hydrothermal systems, and recognize deformations along seismogenic and aseismic NNW-SSE and NW-SE faults. The deformations localized along the Southern Apennines chain are mainly related to landslides while those occurring in the plains are due to subsidence processes induced by intensive drainage from wells, i.e. anthropic activity. The review of 9 years of SAR data shows that tectonic, volcanic/hydrothermal, gravity, and anthropic processes are responsible for the ground deformation of Campania. The proposed joint interpretation of deformation fields related to natural and anthopogenic factors provides a comprehensive view of the dynamics of the Earth’s surface.     

Detecting and quantifying mountain permafrost creep from in situ inventory,space-borne radar interferometry and airborne digital photogrammetry

In this paper three different techniques for detecting and quantifying mountain permafrost creep are compared: (1) rock glacier inventory and characterization from in situ indicators, (2) space-borne radar interferometry, and (3) digital processing of repeated airborne imagery. The specific characteristics of the three methods and their complementarity are investigated for the Fletschhorn mountain range in the Simplon/Saas valley region, Swiss Alps. We found that radar interferometry is suitable to quantify the degree of activity and the order of surface velocity of rock glaciers over large areas in one process, with the possibility to also detect very small movements of inactive and relict rock glaciers. On the other hand, aero-photogrammetry represents a valuable base for additional interpretation of the three-dimensional surface flow field (including speed, direction and change in thickness) of the most active rock glaciers. Results from radar interferometry can also form the basis for further detailed in situ investigations.     

Process models in the practice of distributed software development: A systematic review of the literature

ContextDistributed Software Development (DSD) has recently become an active research area. Although considerable research effort has been made in this area, as yet, no agreement has been reached as to an appropriate process model for DSD.PurposeThis paper is intended to identify and synthesize papers that describe process models for distributed software development in the context of overseas outsourcing, i.e. “offshoring”.MethodWe used a systematic review methodology to search seven digital libraries and one topic-specific conference.ResultsWe found 27 primary studies describing stage-related DSD process models. Only five of such studies looked into outsourcing to a subsidiary company (i.e. “internal offshoring”). Nineteen primary studies addressed the need for DSD process models. Eight primary studies and three literature surveys described stage-based DSD process models, but only three of such models were empirically evaluated.ConclusionWe need more research aimed at internal offshoring. Furthermore, proposed models need to be empirically validated.     

Load distribution to minimise pressure-related pain on foot: a model

The optimal force distribution to minimise pain or discomfort at the foot–shoe interface is still not known. Most shoe-related products attempt to distribute the load uniformly without much consideration to the bony and soft tissue regions. An experiment was conducted to first determine the pressure pain threshold (PPT) and tissue deformation on the plantar surface of the foot. Circular probes of areas 0.5, 1.0 and 2.0 cm² at indentation speeds of 0.5, 1 and 2 mm/s showed that PPT depends on the location stimulated, area of stimulation and the indentation speed. The results also showed that tissue stiffness is quite low for small deformations ( < 4 mm), but significantly higher at large deformations (>4 mm). The stiffness at the larger deformation region was positively correlated with PPT (r = 0.63, p < 0.001). The data were further used to develop a model with PPT, deformation and stimulated area.

Practitioner Summary: Pressure at which there is an onset of pain is higher when a larger area of soft tissue is stimulated. Bony areas may be better suited to bear load on smaller areas to minimise pressure-related pain. Thus, manipulating supporting surface stiffness and surface contours can help minimise pain.     

Coseismic deformation and source model of the 12 November 2017 MW 7.3 Kermanshah Earthquake (Iran–Iraq border) investigated through DInSAR measurements

A large earthquake with a magnitude of M_W 7.3 struck the border of Iran and Iraq at the province of Kermanshah, Iran. In our study, coseismic deformation and source model of the 12 November 2017 Kermanshah Earthquake are investigated using ALOS-2 ScanSAR and Sentinel-1A/B TOPSAR Differential Interferometric Synthetic Aperture Radar (DInSAR) techniques. Geodetic inversion has been performed to constrain source parameters and invert slip distribution on the fault plane. The optimised source model from joint inversion shows a blind reverse fault with a relatively large right-lateral component, striking 353.5° NNW-SSE and dipping 16.3° NE. The maximum slip is up to 3.8 m at 12–14 km depth and the inferred seismic moment is 1.01 × 10²⁰ Nm, corresponding to M_W 7.3, consistent with seismological solutions. The high-resolution optical images from SuperView-1 satellite suggest that most of the linear surface features mapped by DInSAR measurements are landslides or surface cracks triggered by the earthquake. Coulomb stress changes on the source fault indicating consistency between aftershock distribution and high loaded stress zones. Based on the stress change on neighbouring active faults around this area, the Kermanshah Earthquake has brought two segments of the Zagros Mountain Front Fault (MFF), MFF-1 and MFF-2, 0.5–3.1 MPa and 0.5–1.96 MPa closer to failure, respectively, suggesting the risk of future earthquakes. Recent major aftershocks (M_W≥ 5.0) could probably ease the seismic hazard on MFF-2, but the risk of earthquakes on MFF-2 is still increasing.     

广义质心及其在仿射变换参数恢复中的应用

目的 全局的仿射变换配准需估计出仿射变换的参数,现有算法要么效果不佳,要么对二值图像无能为力。本文改造传统质心的定义,提出广义质心的概念。方法 传统的质心以二重积分定义,所提广义质心利用变形累次积分定义,传统质心只是这种广义质心的特例。本文给出了广义质心保持仿射变换前后对应关系的条件,并提出了一种利用这种广义质心进行仿射变换参数恢复的算法。结果 该算法对灰度和二值图像的仿射变换参数恢复都适用,实验结果也表明现有的交叉权重矩方法耗时是本文算法耗时的25倍,但它们的恢复效果相差不大,并且本文算法要比现有的图像矩构造非线性方程组方法恢复效果好。结论 本文提出了广义质心,利用这种广义质心进行仿射变换参数恢复算法,对二值图像和灰度图像均适用,恢复效果较好,并且计算量较小。     

Landslide monitoring by corner reflectors differential interferometry SAR

Although the differential interferometric synthetic aperture radar (DInSAR) technique has the potential for monitoring ground deformation with millimetric accuracy, it still suffers from temporal and geometric decorrelation. In this paper, the corner reflectors differential interferometry synthetic aperture radar (CRDInSAR) technique was used to overcome the limitations of conventional DInSAR. We studied the basic principles of CRDInSAR, discussed the calculation of the flat earth and topographic phases based on the geometry of satellite and corner reflectors, presented the phase unwrapping approach for the sparse grid of corner reflectors, then investigated the construction of, and solution to, the unwrapped phase system equation. Subsequently, we applied CRDInSAR to monitor the displacement of the Shuping landslide in Hubei Province, China. In this study, we processed five SAR images on a descending pass acquired by the Environmental Satellite (ENVISAT) Advanced Synthetic Aperture Radar (ASAR) sensor from September 2005 to March 2006, and compared the achieved results of CRDInSAR with Global Positioning System (GPS) measurements at the same time of the SAR data acquisition assumed as reference. The result indicates a good agreement between the measurements provided by the two different techniques, which shows that CRDInSAR allows monitoring of slow landslide deformation in low coherence areas and provides accurate results.     

Choice of the objective function for pseudogradient measurement of image parameters

Approaches to the choice of objective functions and the form of the pseudogradient are considered for pseudogradient evaluation of interframe geometrical image deformations. The possibilities of reducing the computational cost of calculation of the pseudogradient are analyzed. Minkina Galina Leonidovna. Born in 1983. Fifth-year student at Ulyanovsk State Technical University, specialty applied mathematics. Scientific interests: optimization of parameters of algorithms for estimating interframe geometrical image deformations. Samoilov Mikhail Yurievich. Born in 1981. Graduated from Moscow State University in 2003. PhD student at Ulyanovsk State Technical University. Scientific interests: minimization of the computational cost in estimation of interframe geometrical image deformations. Tashlinskii Alexander Grigorievich. Born in 1954. Graduated from the Ulyanovsk Polytechnical Institute in 1977. Received doctoral degree in 2000. Professor in the Department of Systems of Computer Aided Design. Scientific interests: statistical image processing, in particular, estimation of spatiotemporal deformations of sequences of dynamical images. Author of 200 publications, including more than 60 papers and a monograph on estimation of parameters of spatial deformations of sequences of images. Full member of the International Academy of Authors of Scientific Discoveries and Investigations and of the Russian Academy of Natural Sciences. Awarded medals from these academies.     

Deformation retrieval in large areas based on multibaseline DInSAR algorithm: a case study in Cangzhou,northern China

The DInSAR technique with a multibaseline is becoming popular nowadays to investigate slow urban deformation. In this paper, we focus on deformation retrieval in large areas, including urban and suburban areas. Based on the multibaseline DInSAR algorithm proposed by Mora, three extensions are derived. First, least‐squares adjustment and error‐controlling methods are used to obtain stable deformation velocity and height error estimations. The least‐squares QR factorizaiton algorithm is emphasized to solve large, linear, and sparse functions. Second, a new complex network is presented to limit noise effects on the Delaunay triangular network. Third, by combining complex and Delaunay networks, large‐area deformation is investigated, from centre urban areas to suburban areas. The enhanced algorithm is performed to investigate the subsidence of Cangzhou, Hebei province (northern China) during 1993–1997 by using 9 ERS SLC data. The experimental results show serious subsidence in the region and are validated by levelling data and groundwater wells data. Compared with levelling data, the estimation errors of linear deformation velocity in urban areas are in the range of (?2, 2) mm year^?1, and in suburban areas, the errors are in the range of (?26, 15) mm year^?1, which is sufficiently feasible to determine the status of subsidence relative to the maximum deformation velocity of about ?100 mm year^?1. The subsidence centres in urban areas are consistent with the spatial distribution of groundwater wells, which provides evidence that groundwater overexploitation is the main cause of subsidence in Cangzhou. The closure of wells will be a good way to control subsidence in the future.