Subsidence Bowl Research Articles

Characteristics of the Residual Surface Deformation of Multiple Abandoned Mined-Out Areas Based on a Field Investigation and SBAS-InSAR: A Case Study in Jilin, China

Residual surface deformation resulting from abandoned mined-out areas can lead to severe damage to ground structures (e.g., buildings and infrastructure in mining areas) and the local ecological environment. Long-term monitoring and analyses of surface deformation characteristics of abandoned mined-out areas are significant for preventing potential disasters. In this study, a detailed field investigation first was conducted in Ying’an and Baoshan coal mines located in Jilin Province, China, to survey mining-induced disasters in the mining areas. Based on the 40 Sentinel-1A images acquired from 14 February 2017 to 17 May 2020, small baseline subset interferometry synthetic aperture radar (SBAS-InSAR) technology was employed to obtain the time-series residual surface deformation. Validation of the SBAS-derived results is performed by comparing with the results obtained via leveling measurements. The root mean square error (RMSE) between SBAS-derived and leveling measurements results was found to be 1.144 mm, reflecting a fairly good agreement. Furthermore, the ordinary Kriging interpolation approach was adopted to obtain information on the deformation across the entire area. The spatial–temporal evolution characteristics of the derived subsidence bowls in multiple mined-out areas were revealed. The deformation characteristics for the abandoned mined-out areas in different periods were not completely consistent. Finally, the potential mechanism underlying the inconsistency in the subsidence associated with underground coal exploitation is analyzed. The findings of this study can provide insights into local construction and ecological improvement as well as guidance for the prediction of deformation in abandoned mined-out areas.

Present-day land subsidence rates, surface faulting hazard and risk in Mexico City with 2014–2020 Sentinel-1 IW InSAR

Among the fastest sinking cities globally, the metropolitan area of Mexico City is the target of an unprecedented satellite investigation based on over 300 Sentinel-1 Synthetic Aperture Radar (SAR) Interferometric Wide swath mode scenes acquired in 2014–2020. Two-pass differential Interferometric SAR (InSAR) and the parallelized Small BAseline Subset (SBAS) repeat-pass InSAR approach provide a complete account of spatial patterns, long-term trend and present-day settlement rates affecting the city. The 3D deformation field reveals that foremost is the role of the vertical velocity VU, with peaks of −38.7 cm/year in Nezahualcóyotl, −32.0 cm/year in Gustavo A. Madero and Venustiano Carranza, and −39.1 cm/year in Iztapalapa. Settlement at the metropolitan Cathedral in Cuauhtémoc is ongoing at up to −8.8 cm/year, consistently with the last six decades. Volcanic edifices mark ground stability “islets” inside the main subsidence bowls. East-west rates are limited, except for some horizontal strain (up to ±5 cm/year) within the subsidence bowls. Comparison with surface geology and geotechnical zoning confirms that aquitard compaction is the predominant process. The power relationship between VU [cm/year] and the thickness of lacustrine clay deposits HC [m] is: VU=176∗HC1.8. The 2019–2020 deformation scenario shows that subsidence still involves most of Nezahualcóyotl, with −3.0 cm/month VU. A well-defined long-term deformation behaviour comes out from 2014–2019 InSAR time series and comparison with 2008–2020 GPS data. RMSE of 0.9 cm is found at MMX1 station deployed within the lacustrine unit, and 0.6 cm at TNGF station onto hard volcanic rocks. The sharpest subsidence gradients and angular distortions β in 2017–2019 (up to over 1/400, i.e. 0.14°) – thus the greatest vulnerability to surface faulting and cracking as a consequence of large tensile stress in the soil caused by differential settlement – are found at the foothills of Sierra de Santa Catarina, Peñón del Marqués, Cerro Chimalhuachi, Peñón de los Baños and Sierra de Chichinautzin, where the transition unit is narrower (or absent). Faults and cracks develop where β > 1/3000, i.e. 0.03%, in 2017–2019. The observed width of the influence zone (i.e. damage band), where β is still significant to induce damage, is 250 m. Differential settlement and surface faulting could compromise the serviceability of housing and utility infrastructure. It is estimated that over 457,000 properties and 1.5 million inhabitants of the Valley of Mexico Metropolitan Area (ZMVM) are in zones at high to very high surface faulting risk, mainly in Iztapalapa, Tláhuac, Chimalhuacán and Chalco. Increased flood exposure due to formation of topographic depressions involves over 751,000 properties and ~2.7 million inhabitants of the ZMVM, mainly in Nezahualcóyotl, Tláhuac, Venustiano Carranza, Iztapalapa, Gustavo A. Madero and Ecatepec de Morelos. These municipalities are often hit by catastrophic floods.

Ground Deformation and Its Causes in Abbottabad City, Pakistan from Sentinel-1A Data and MT-InSAR

Land subsidence, as one of the engineering geological problems in the world, is generally caused by compression of unconsolidated strata due to natural or anthropogenic activities. We employed interferometric point target analysis (IPTA) as a multi-temporal interferometric synthetic aperture radar (MT-InSAR) technique on ascending and descending Sentinel-1A the terrain observation with progressive scans SAR (TOPSAR) images acquired between January 2015 and December 2018 to analyze the spatio-temporal distribution and cause of subsidence in Abbottabad City of Pakistan. The line of sight (LOS) average deformation velocities along ascending and descending orbits were decomposed into vertical velocity fields and compared with geological data, ground water pumping schemes, and precipitation data. The decomposed and averaged vertical velocity results showed significant subsidence in most of the urban areas in the city. The most severe subsidence was observed close to old Karakorum highway, where the subsidence rate varied up to −6.5 cm/year. The subsidence bowl profiles along W–E and S–N transects showed a relationship with the locations of some water pumping stations. The monitored LOS time series histories along an ascending orbit showed a close correlation with the rainfall during the investigation period. Comparative analysis of this uneven prominent subsidence with geological and precipitation data reflected that the subsidence in the Abbottabad city was mainly related to anthropogenic activities, overexploitation of water, and consolidation of soil layer. The study represents the first ever evidence of land subsidence and its causes in the region that will support the local government as well as decision and policy makers for better planning to overcome problems of overflowing drains, sewage system, littered roads/streets, and sinking land in the city.

LAND SUBSIDENCE MONITORING WITH MULTI-TRACK SAR DATA

Abstract. In this study, we investigated wide-area land subsidence in Hebei Plain using 78 C-band Sentinel-1A SAR images acquired from May 2017 to May 2018 covering two tracks. High-precision time series retrieval was performed by NSBAS technology based on a single track. The offset deformation values of two tracks calculated in this paper were −0.09 mm/year, and the multi-track deformation rate was successfully merged. Using the cross-validation of redundant observations in the multi-track overlap area to evaluate the accuracy, and it was found that 90% pixel difference between the two track overlap areas was within 9 mm. The standard deviation was 5.38 mm, and the subsidence trend of the overlap area was consistent. Twenty-four subsidence bowls were extracted, and the maximum subsidence rate reached −62 mm/year.

Quantifying Ground Subsidence Associated with Aquifer Overexploitation Using Space-Borne Radar Interferometry in Kabul, Afghanistan

Rapid population growth combined with recent drought events and decades of political instability have left the residents of Kabul facing water scarcity, significantly relying on groundwater. Groundwater overexploitation might have induced various magnitudes of ground subsidence, however, to date, no comprehensive study of ground subsidence in Kabul has been conducted. In this study, we investigated the spatio-temporal evolution of ground deformation phenomena and its main governing processes in Kabul from 2014 to 2019 using C-Band Sentinel-1 derived Interferometric Synthetic Aperture Radar (InSAR) time-series from both ascending and descending orbits to extract the two-dimensional (2D) surface displacement field. Four subsidence bowls were distinguished with highly variable spatial extents and deformation magnitudes over four separate aquifer basins, with the maximum value of −5.3 cm/year observed in the Upper Kabul aquifer basin. A wavelet analysis suggests that there is a strong correlation between the groundwater level variations and subsidence. Investigation of hydrogeological data further reveals that the observed subsidence could be attributed to the presence of highly compressible clayey soils. This detailed space-borne regional survey provides new insights into the main governing mechanism of land subsidence in Kabul and may direct better mitigation plans of potential hazards.

Identification of the ground movements caused by mining-induced seismicity with the satellite interferometry

Abstract. The assessment of the impact of mining-induced seismicity on the natural environment and infrastructure is often limited to the analysis of terrain surface vibrations. However, similar seismic phenomena, like earthquakes, may also imply dislocations and deformations of the rock mass. Such ground movements may occur in areas which are not directly under the influence of the mining. The study of the displacement field caused by mining-induced seismicity is usually carried out with the use of geodetic methods. Classical geodetic measurements provide discrete information about observed ground movements. As a result, they generally do not provide spatially and temporally relevant estimates of the total range and values of ground movements for specific periods of interest. Moreover, mining-induced seismicity causes a severe threat to buildings. That is why, regarding the complexity of the mechanism of occurrence of mining-induced seismicity and their impact on ground movements, this problem remains a substantial research issue. The presented research aimed to analyse the ground movements caused by mining-induced seismicity. The ground displacements were established based on data from Sentinel-1 satellites applying differential interferometric synthetic aperture radar (DInSAR). The results of the investigation in the copper mining area of the Lower Silesia region of Poland revealed that the observed subsidence caused by mining-induced seismicity usually has a shape of a regular ellipse. The radius of these ground movements does not exceed approximately 2–3 km from the mining-induced tremor's epicenter, and the total subsidence reaches ca. 10–20 cm. More than 50 % of the total subsidence is observed on the surface within a few days after the mining tremor occurrence. Furthermore, the deformations of the surface occur when the energy of mining-induced tremor reaches values of the order of 105 J or higher. The presented research can contribute to better identification and evaluation of the mechanism of the rock mass deformation process caused by mining-induced seismicity. In addition, the use of satellite radar interferometry improves the quality of monitoring of these dynamic phenomena significantly. The data retrieved using this method allow for quasi-continuous monitoring of the local subsidence bowls caused by mining-induced seismicity.

Reduced rate of land subsidence since 2016 in Beijing, China: evidence from Tomo-PSInSAR using RadarSAT-2 and Sentinel-1 datasets

ABSTRACTLand subsidence associated with groundwater extraction in the city of Beijing, China, has been a problem for decades. Remote sensing has been used extensively in prior studies to monitor subsidence in Beijing. However, given recent changes in precipitation and groundwater management, there is an urgent need to update the subsidence record and to evaluate whether the long-term spatiotemporal patterns of subsidence have changed. This study therefore investigates the recent spatiotemporal patterns of land subsidence in Beijing by tomography-based persistent scatterer interferometry SAR (Tomo-PSInSAR) technology, using 39 RadarSAT-2 images from 2012 to 2015 and 33 Sentinel-1 images from 2016 to May-2018, and drawing upon Geographic Information System (GIS) spatial analysis methods. Vertical ground deformation rates in Beijing were found to range from −176.2 to +12.3 mm year−1 from 2012 to 2015, but subsequently decreased to −119 to +8 mm year−1 from 2016 to May 2018. Three spatial scale of subsidences are evident: At the metropolitan regional-scale, the total area of subsidence area is about 1235.2 km2, and comprises four main subsiding regions, located in the northern and eastern parts of the city. More than 85% of the subsiding area is located between the Fifth and Sixth Ring Roads. At a more local scale, eight main subsidence bowls are characterized by different patterns of subsidence. Some of the subsidence bowls are separated by active faults. Time-series data of the displacement show that the decreasing subsidence rate after 2016 could be due to the 1 m rise in mean groundwater level from the end of 2014 to mid-2018. This change in groundwater level is likely due to an increase in precipitation since 2016, and water transfers, which reached 2.3 × 109 m3 by 2017 from the South-North Water Transfer Project. At the scale of individual infrastructure projects, the Beijing subway, main roads and the Capital Airport all show severe uneven subsidence, which is a cause for concern. To our knowledge, this research is the first study using satellite SAR remote sensing methods to document the change in the land subsidence rate of Beijing. Starting in 2016, the rate notably declined, suggesting that subsidence mitigation strategies are beginning to have an effect.

Remotely sensing large- and small-scale ground subsidence: A case study of the Guangdong–Hong Kong–Macao Greater Bay Area of China

As the world's largest city cluster, the Guangdong–Hong Kong–Macao Greater Bay Area (GBA) is vulnerable to significant subsidence. The widely distributed sediments and rapid urbanization in the GBA result in the concurrence of large- and small-scale subsidence. Mono-sensor synthetic aperture radar (SAR) images usually suffer from the limitation of either low resolution or small coverage, and thus, are not applicable to completely monitoring the GBA. In this study, we used Sentinel-1 (S1), COSMO-SkyMed (CSK) and TerraSAR-X (TSX) images jointly to reveal multi-scale subsidence of the GBA. The overall subsidence ranged from 0 to 112.3 mm/yr derived from the 2015–2017 S1 images. Three regional subsidence bowls (Zhuhai–Zhongshan, Jiangmen, and Guangzhou–Zhongshan) formed in the western alluvial plain. The high correlation between regional subsidence and Quaternary sediments confirms that sediment consolidation is the main cause of subsidence. The land use and numerical modeling results suggest that groundwater extraction and artificial loading are the triggering factors. Two representative local subsidence cases were analyzed using high-resolution images: reclamation settlement at the Hong Kong International Airport (HKIA) and sinkhole subsidence due to the excavation of the Shenzhen Wenbo skyscraper foundation pit. The validation at the HKIA showed that the measurements from the CSK and S1 data both agreed well with the leveling data. However, CSK outperformed S1 in the sense that it increased the point density by a factor of 5, improved the height precision by a factor of 4, and showed fewer false alarms. CSK is therefore more applicable to monitoring the local subsidence of key infrastructures. The cross-heading tracks of TSX and CSK images detected precursory subsidence before the sinkhole collapse from two sides, indicating that the cross-heading tracks benefit the comprehensive monitoring of local subsidence in high-rise and high-density built environments. In summary, the synergistic use of multi-sensor SAR images demonstrates the practicability of the operational surveillance of multi-scale subsidence from regional surveying to the fine monitoring of local areas in the GBA.

Current Activity of the Long Point Fault in Houston, Texas Constrained by Continuous GPS Measurements (2013–2018)

The Long Point Fault is one of the most active urban faults in Houston, Texas, which belong to a complex system of normal growth faults along the Texas Gulf Coast. To assess the activity of the Long Point Fault, a GPS array with 12 permanent stations was installed along the two sides of the 16-km-long fault scarp in 2013. GPS datasets were processed with the Precise Point Positioning (PPP) and Double-Difference (DD) methods. The daily PPP solutions with respect to the International Global Navigation Satellite System (GNSS) Reference Frame 2014 (IGS14) were converted to the Stable Houston Reference Frame (Houston16). The six-year continuous GPS observations indicate that the Long Point Fault is currently inactive, with the rates of down-dip-slip and along-strike-slip being below 1 mm/year. The Long Point Fault area is experiencing moderate subsidence varying from 5 to 11 mm/year and a coherent horizontal movement towards the northwest at a rate of approximately 2 to 4 mm/year. The horizontal movement is induced by the subsidence bowl that has been developing since the 1980s in the Jersey Village area. Current surficial damages in the Long Point Fault area are more likely caused by ongoing uneven subsidence and its induced horizontal strains, as well as the significant seasonal ground deformation, rather than deep-seated or tectonic-controlled fault movements. The results from this study suggest a cause-and-effect relationship between groundwater withdrawals and local faulting, which is pertinent to plans for future urban development, use of groundwater resources, and minimization of urban geological hazards.

Research on development characteristics and failure mechanism of land subsidence and ground fissure in Xi'an, monitored by using time-series SAR interferometry

The development characteristics of land subsidence of Xi'an by an advanced InSAR technique are investigated based on Envisat, Radarsat-2 and Sentinel-1 satellites, and it can be found that there still have been five typical land subsidence bowls with maximum subsidence rates beyond 40 mm/year in Xi'an. In 2012, the land subsidence took on slightly slowing down with a ground rebound area of up to 6.3 km2. But in 2003–2010, 2016 and 2017, there was more serious land subsidence and especially in 2017 the maximum velocity reached up to −137 mm/year. Yuhuazhai-Nanyaotou subsidence bowl as the greatest subsidence bowl, its area that rate was below 50 mm/year in 2017 was up to 10.64 km2. Moreover, these land subsidence bowls mainly concentrate upon tectonic depression area due to huge thickness clay soil and are spatially and temporally coincident with confined aquifer depression cones, suggesting that over-exploitation groundwater plays a decisive role for land subsidence. Ground fissure belts are generally controlled by normal faults in Xi'an, but asymmetry settlement of land surface can also give rise to the abnormal activities of ground fissure belts, which have caused serious threat or damage to building structures. In addition, it could also be found that surface uplift area in 2012 involves in both the tectonic depression area and loess ridge in terms of their same thickness sand layer.

Ground surface response to continuous compaction of aquifer system in Tehran, Iran: Results from a long-term multi-sensor InSAR analysis

Most of the developed groundwater basins in Iran are subject to land subsidence hazards resulting from the over-extraction of groundwater. Several areas in Tehran, the capital city and a provincial center in north-central Iran, have been reported to be subsiding at different rates. In this study, we present the results of an Interferometric Synthetic Aperture Radar (InSAR) time series analysis of Tehran using different SAR data between 2003 and 2017. By constructing more than 400 interferograms derived from 39 Envisat ASAR (C-band), 10 ALOS PALSAR (L-band), 48 TerraSAR-X (X-band), and 64 Sentinel-1 (C-band) SAR datasets, we compile displacement time series from interferometric observations using the Small Baseline (SB) technique. Our analysis identifies 3 distinct subsidence features in Tehran with rates exceeding 25 cm/yr in the western Tehran Plain, approximately 5 cm/yr in the immediate vicinity of Tehran international airport, and 22 cm/yr in the Varamin Plain to the southeast of Tehran city. The temporal pattern of land subsidence, which is dominated by a decreasing trend, generally follows the regional decline in groundwater level, which suggests that anthropogenic processes caused by excessive groundwater extraction are the primary cause of land subsidence. Integrating a decadal time series of subsidence constructed from multi-sensor InSAR with in-situ observations suggests that inelastic and permanent compaction dominates the main deformation regime of the Tehran aquifer, and the ratio between elastic and inelastic deformation is approximately 0.4. A geological analysis indicates that the shape of the subsidence bowl in the western Tehran Plain does not follow the trend of major mapped faults in the region. In contrast, the subsidence bowl in Varamin is controlled by the Pishva Fault, which suggests that either this fault acts as a hydrologic barrier to the groundwater flow in this region or that the differences in sediment thickness causes the discontinuity in land subsidence.

Multi-Scale Analysis of the Relationship between Land Subsidence and Buildings: A Case Study in an Eastern Beijing Urban Area Using the PS-InSAR Technique

Beijing is severely affected by land subsidence, and rapid urbanisation and building construction might accelerate the land subsidence process. Based on 39 Envisat Advanced Synthetic Aperture Radar (ASAR) images acquired between 2003–2010, 55 TerraSAR-X images acquired between 2010–2016, and urban building information, we analysed the relationship between land subsidence and buildings at the regional, block, and building scales. The results show that the surface displacement rate in the Beijing urban area ranged from −109 mm/year to +13 mm/year between 2003–2010, and from −151 mm/year to +19 mm/year between 2010–2016; two subsidence bowls were mainly distributed in the eastern part of the Chaoyang District. The displacement rate agreed well with the levelling measurements, with an average bias of less than six mm/year. At the regional scale, the spatial pattern of land subsidence was mainly controlled by groundwater extraction, compressible layer thickness, and geological faults. Subsidence centres were located in the area around ground water funnels with a compressible layer depth of 50–70 m. The block-scale analysis demonstrated a clear correlation between the block construction age and the spatial unevenness of subsidence. The blocks constructed between 1998–2005 and after 2005 showed considerably more subsidence unevenness and temporal instability than the blocks constructed before 1998 during both time periods. The examination of the new blocks showed that the spatial unevenness increased with building volume variability. For the 16 blocks with a high building volume, variability, and subsidence unevenness, the building-scale analysis showed a positive relationship between building volume and settlement in most blocks, although the R2 was lower than 0.5. The results indicate that intense building construction in urban areas could cause differential settlement at the block scale in Beijing, while the settlement of single buildings could be influenced by the integrated effects of building volume, foundation structures, and the hydrogeological background.

Spatiotemporal Evolution of Land Subsidence in the Beijing Plain 2003–2015 Using Persistent Scatterer Interferometry (PSI) with Multi-Source SAR Data

Land subsidence is one of the most important geological hazards in Beijing, China, and its scope and magnitude have been growing rapidly over the past few decades, mainly due to long-term groundwater withdrawal. Interferometric Synthetic Aperture Radar (InSAR) has been used to monitor the deformation in Beijing, but there is a lack of analysis of the long-term spatiotemporal evolution of land subsidence. This study focused on detecting and characterizing spatiotemporal changes in subsidence in the Beijing Plain by using Persistent Scatterer Interferometry (PSI) and geographic spatial analysis. Land subsidence during 2003–2015 was monitored by using ENVISAT ASAR (2003–2010), RADARSAT-2 (2011–2015) and TerraSAR-X (2010–2015) images, with results that are consistent with independent leveling measurements. The radar-based deformation velocity ranged from −136.9 to +15.2 mm/year during 2003–2010, and −149.4 to +8.9 mm/year during 2011–2015 relative to the reference point. The main subsidence areas include Chaoyang, Tongzhou, Shunyi and Changping districts, where seven subsidence bowls were observed between 2003 and 2015. Equal Fan Analysis Method (EFAM) shows that the maximum extensive direction was eastward, with a growing speed of 11.30 km2/year. Areas of differential subsidence were mostly located at the boundaries of the seven subsidence bowls, as indicated by the subsidence rate slope. Notably, the area of greatest subsidence was generally consistent with the patterns of groundwater decline in the Beijing Plain.

Pumping-induced stress and strain in aquifer systems in Wuxi, China

Excessive groundwater withdrawal from an aquifer system leads to three-dimensional displacement, causing changes in the states of stress and strain. Often, land subsidence and sometimes earth fissures ensue. Field investigation indicates that land subsidence and earth fissures in Wuxi, a city in eastern China, are mainly due to excessive groundwater withdrawal, and that they are temporally and spatially related to groundwater pumping. Groundwater withdrawal may cause tensile strain to develop in aquifer systems, but tensile strain does not definitely mean tensile stress. Where earth fissures are concerned, the stress state should be adopted in numerical simulations instead of the strain state and displacement. The numerical simulation undertaken for the Wuxi area shows that the zone of tensile strain occupies a large area on the ground surface; nevertheless, the zone of tensile stress is very limited. The zone of tensile stress often occurs near the ground surface, beneath which the depth to the bedrock surface is relatively small and has considerable variability. Earth fissures often initiate near the ground surface where tensile stress occurs. Tensile stress and earth fissures rarely develop at the centers of land subsidence bowls, where compressive stress is dominant.

Wuhan Surface Subsidence Analysis in 2015–2016 Based on Sentinel-1A Data by SBAS-InSAR

The Terrain Observation with Progressive Scans (TOPS) acquisition mode of Sentinel-1A provides a wide coverage per acquisition and features a repeat cycle of 12 days, making this acquisition mode attractive for surface subsidence monitoring. A few studies have analyzed wide-coverage surface subsidence of Wuhan based on Sentinel-1A data. In this study, we investigated wide-area surface subsidence characteristics in Wuhan using 15 Sentinel-1A TOPS Synthetic Aperture Radar (SAR) images acquired from 11 April 2015 to 29 April 2016 with the Small Baseline Subset Interferometric SAR (SBAS InSAR) technique. The Sentinel-1A SBAS InSAR results were validated by 110 leveling points at an accuracy of 6 mm/year. Based on the verified SBAS InSAR results, prominent uneven subsidence patterns were identified in Wuhan. Specifically, annual average subsidence rates ranged from −82 mm/year to 18 mm/year in Wuhan, and maximum subsidence rate was detected in Houhu areas. Surface subsidence time series presented nonlinear subsidence with pronounced seasonal variations. Comparative analysis of surface subsidence and influencing factors (i.e., urban construction, precipitation, industrial development, carbonate karstification and water level changes in Yangtze River) indicated a relatively high spatial correlation between locations of subsidence bowl and those of engineering construction and industrial areas. Seasonal variations in subsidence were correlated with water level changes and precipitation. Surface subsidence in Wuhan was mainly attributed to anthropogenic activities, compressibility of soil layer, carbonate karstification, and groundwater overexploitation. Finally, the spatial-temporal characteristics of wide-area surface subsidence and the relationship between surface subsidence and influencing factors in Wuhan were determined.

Integrated space geodesy for mapping land deformation over Choushui River Fluvial Plain, Taiwan

ABSTRACTGroundwater is an important part of the precious water resources. As the fresh surface water resources become scarcer because of climate change, population growth, and industrial activities, more and more groundwater has been extracted to meet the demands of various water uses (e.g. municipal, industrial, and agricultural). Excessive groundwater extraction leads to severe ground subsidence which compromises the safety of surface and underground infrastructures. Modelling the effects of groundwater extraction is vital to the management and sustainable use of groundwater. However, results of such modelling have to be validated with inputs such as the field survey of ground subsidence. Levelling and continuous global positioning system (GPS) receiver networks are routinely used to collect these field measurements. Unfortunately, these techniques have limitations in terms of areal coverage and density of survey marks and, as a result, subsidence hot spots can be easily missed out. In order to provide a comprehensive picture of subsidence to aid geotechnical modelling and to assess the effectiveness of measures used to mitigate ground subsidence, satellite imaging radar interferometry techniques (interferometric synthetic aperture radar (InSAR) can be used to complement other deformation monitoring techniques. In this study, 20 Advanced Land Observing Satellite (ALOS) Phased Array L-band Synthetic Aperture Radar (PALSAR) images acquired from 31 December 2006 to 26 February 2011 were used to map the land displacement over the Choushui River Fluvial Plain (CRFP), Taiwan. The GPS measurements acquired at 10 continuously operating reference stations (CORS) were used to refine the orbit error in the each differential interferogram obtained from each radar image pair. The displacement time series over the distributed scatterers and the persistent scatterers were analysed. Several subsidence bowls were identified in CRFP. A quantitative comparison was conducted to compare the radar measurements to the GPS measurements over 36 GPS CORS stations. Good agreement between both measurements was observed with coefficient of determination (R2) of 0.97, absolute mean difference of 3.2 mm year−1, and standard deviation of 4 mm year−1. The InSAR-measured Line-of-Sight displacement and GPS-measured horizontal displacement were integrated to derive the vertical displacement map. Two displacement maps were generated using two ALOS-2 PALSAR-2 pairs acquired between 2015 and 2016. Similar subsidence patterns were found in the two maps compared to the 2006–2011 displacement rate map, suggesting the land over the same region might have continued to fall.

Unwrapping-free interpolation of sparse DInSAR phase data: experimental validation

ABSTRACTIn multi-temporal applications of synthetic aperture radar (SAR) interferometry, differential phase contributions due to atmospheric inhomogeneities, estimated over sparse points, have to be interpolated and removed from the regular-grid interferograms in order to highlight the phase stability of more image pixels, which then add to the available data to infer useful information about terrain displacements or other phenomena of interest. Interpolation is usually done on the phase data after a phase unwrapping (PU) operation. In a previous work, we considered the alternative interpolation step applied directly to the complex phasor derived from the wrapped phase, thus bypassing the error-prone sparse PU operation. In this article, the performances of the proposed methodology are evaluated over atmospheric phase screen (APS) data estimated from a previous processing through persistent scatterers interferometry (PSI) methods. The original persistent scatterer (PS) population is reduced by thresholding their inter-image coherence values, and then further subsampled randomly in a rectangle inside a detected subsidence bowl. Both the classical and the proposed interpolation procedures are applied to the subsampled APS phase values. The interpolated fields are then removed from the rest of the PS, and the residual phase values are compared in terms of inter-image coherence. Results confirm that interpolating complex phasors, thus avoiding PU, gives results equivalent to the standard procedure in good sampling conditions. Moreover, when point sparsity induces phase aliasing, thus hindering the PU operation, the proposed method allows to better recover phase information over unsampled pixels, improving the final results of the PSI processing.

Characteristics of land subsidence, earth fissures and related disaster chain effects with respect to urban hazards in Xi’an, China

The ancient capital of Xi’an is a typical city in China with a shortage of water resources. With an urban population of 7 million, its domestic and industrial use of water depends almost entirely on groundwater. Large-scale long-term extraction of groundwater has caused severe land subsidence and triggered the reactivation of active normal faults in the Quaternary under-consolidated sediments underlying the city. This has further led to ground ruptures that have gradually evolved into earth fissure zones. During this process, there has been a close spatiotemporal relationship between land subsidence and the formation of earth fissures. From field monitoring data, the subsidence bowls are typically located in loess depressions between earth fissure zones within structural blocks of sediment cut by active Quaternary faults and preexisting fault planes. The subsidence bowls are generally elliptical in shape, with the long axes consistent with the preferred orientation of the earth fissures. The underlying Quaternary active faults and preexisting fault planes predispose the intervening sediment blocks to formation of earth fissures, and intense groundwater extraction has led to accelerated subsidence and earth fissuring during the past 50 years. Mining groundwater, land subsidence and earth fissures occur in a sequence and constitute a disaster chain with respect to urban hazards in Xi’an. Thus, it is important to understand the relationships between these events. This paper summarizes the characteristics of land subsidence and earth fissures in Xi’an and discusses the chain of relationships that connects them.

Impact of rock salt creep law choice on subsidence calculations for hydrocarbon reservoirs overlain by evaporite caprocks

AbstractAccurate forward modeling of surface subsidence above producing hydrocarbons reservoirs requires an understanding of the mechanisms determining how ground deformation and subsidence evolve. Here we focus entirely on rock salt, which overlies a large number of reservoirs worldwide, and specifically on the role of creep of rock salt caprocks in response to production‐induced differential stresses. We start by discussing available rock salt creep flow laws. We then present the subsidence evolution above an axisymmetric finite element representation of a generic reservoir that extends over a few kilometers and explore the effects of rock salt flow law choice on the subsidence response. We find that if rock salt creep is linear, as appropriate for steady state flow by pressure solution, the subsidence response to any pressure reduction history contains two distinct components, one that leads to the subsidence bowl becoming narrower and deeper and one that leads to subsidence rebound and becomes dominant at later stages. This subsidence rebound becomes inhibited if rock salt deforms purely through steady state power law creep at low stresses. We also show that an approximate representation of transient creep leads to relatively small differences in subsidence predictions. Most importantly, the results confirm that rock salt flow must be modeled accurately if good subsidence predictions are required. However, in practice, large uncertainties exist in the creep behavior of rock salt, especially at low stresses. These are a consequence of the spatial variability of rock salt physical properties, which is practically impossible to constrain. A conclusion therefore is that modelers can only resort to calculating bounds for the subsidence evolution above producing rock salt‐capped reservoirs.

Monitoring and Analysis of Land Subsidence Along Beijing-Tianjin Inter-City Railway

Land subsidence, as a slow and persistent geological disaster, threatens the safety of urban constructions and linear structures. Beijing-Tianjin intercity railway was built across some of the severe subsidence funnels. By combining with the hydrogeology condition and the InSAR (Interferometric Synthetic Aperture Radar) monitoring result, the corresponding analysis was conducted to reveal the situation and variational rule of the railway subsidence. Several indicators such as gradient, profiles and time-series analysis were used for describing the the land settlement along and across the railway. The result shows that areas with thick compressible layers, frequent over-exploitation and great decline of groundwater has larger subsidence rate. The middle confined groundwater level has much closer trend with subsidence than shallow groundwater. The railway is threatened by the subsidence along it, despite massive efforts are done to avoid crossing the subsidence cones in the railway route selection. The gradient, which is used to describe the uneven subsidence, shows some extrema at the edge area of the subsidence bowl. The maximal gradient value which means the most serious non-uniform settlement of the railway happens at a distance of 14200 m from the beginning. The profiles of the three lines have relatively larger vertical velocity in the middle where the railway goes through. With time goes by, the cumulative non-uniform deformation along the railway should be paid enough attention.