Ground Deformation Measurements Research Articles

Oilfield Reservoir Parameter Inversion Based on 2D Ground Deformation Measurements Acquired by a Time-Series MSBAS-InSAR Method

Time-series ground deformation monitoring and reservoir parameter inversion are crucial for the dynamic assessment of oilfield resources and sustainable exploitation in oilfields. As some of the regions with the richest oil reserves in China, the oilfield areas in the western Qaidam Basin were selected as a typical study area. Firstly, we used SAR images collected by the Sentinel-1A satellite from January 2021 to December 2022 and applied the multidimensional small baseline subset (MSBAS) method to obtain vertical and east–west deformation measurements. On this basis, a nonlinear Bayesian inversion method was applied to model the shallow reservoir in a series of complex deformation areas, based on a single-source model and a multi-source model, respectively. As a result, the ground deformation monitoring results obtained by long time-series InSAR clearly reflect the uneven ground deformation caused by the oil extraction and water injection operation processes. There was slight subsidence in the Huatugou oilfield, while significant uplift deformation occurred in the Ganchaigou oilfield and the Youshashan oilfield, with a maximum uplift rate of 48 mm/year. Further analysis indicated that the introduction of the 2D deformation field helps to improve the robustness of oilfield reservoir parameter inversion. Moreover, the dual-source model is more suitable than the single-source model for inverting reservoir parameters of complex deformation. This study not only fills the gap of InSAR deformation monitoring for the oilfields in the western Qaidam Basin but also provides a theoretical reference for the model and method selection of reservoir parameter inversion in other oilfields.

Analysis of Land subsidencein coastal and urban areas by using various techniques– Literature Review

Over the past several decades, subsidence has emerged as a significant issue affecting urban, coastal, and mining areas worldwide. This phenomenon has been exacerbated by a decrease in groundwater resources, which has become a major problem on a global scale. Insufficient availability of surface water has led to excessive exploitation of the earth's groundwater system through unplanned processes. As a result, serious problems have arisen, including economic challenges, hazardous activities, and damage to the Earth's surface. One of the most significant and impactful consequences of excessive groundwater extraction is land subsidence. During a comprehensive literature review focusing on land displacement, 40 research papers were analyzed and discussed. These studies explored various techniques employed to monitor subsidence, such as synthetic aperture radar (SAR), interferometry synthetic aperture radar (InSAR), persistent scatterer interferometry synthetic aperture radar (PS-InSAR), and differential synthetic aperture radar (D-InSAR). Researchers have consistently found a strong correlation between the over extraction of underground fluids, such as water, oil, and gas, and land subsidence. One prominent cause of land subsidence identified by researchers is the over extraction of underground fluids. When excessive amounts of water, oil, or gas are extracted from underground reservoirs, the resulting voids and reduced pressure cause the overlying land to sink and compact. This process can have severe consequences for human settlements, infrastructure, and ecosystems in affected areas. The depletion of groundwater reserves, in particular, has become a pressing concern because of its vital role in supporting agricultural, industrial, and domestic water needs.Furthermore, the underground construction of tunnels has also been identified as a significant contributor to land subsidence. The excavation of tunnels creates voids and alters the natural balance of underground structures, leading to ground settlement and subsidence. This issue is particularly relevant in urban areas where extensive tunneling projects are undertaken to facilitate transportation, utilities, and infrastructure development. To monitor and assess land subsidence, various methods and processes are employed. SAR, InSAR, PS-InSAR, and D-InSAR techniques are valuable tools in this regard. SAR uses satellite-based radar systems to capture high-resolution images of the Earth's surface, enabling the detection of subtle changes over time. In contrast, InSAR employs two or more radar images to measure ground deformation by analyzing interference patterns. PS-InSAR focuses on persistent scatterers, which are specific targets that exhibit stable radar reflections over time, providing precise measurements of subsidence. D-InSAR combines two radar images acquired at different times to accurately calculate differential ground movements. By using advanced monitoring methods, scientists can assess the extent of subsidence and develop strategies to mitigate its negative impacts. It is crucial to adopt sustainable practices and ensure careful management of groundwater resources to effectively address this pressing issue.

Ionospheric compensation in L-band InSAR time-series: Performance evaluation for slow deformation contexts in equatorial regions

Multi-temporal Synthetic Aperture Radar Interferometry (MT-InSAR) is the only geodetic technique allowing to measure ground deformation down to mm/yr over continuous areas. Vegetation cover in equatorial regions favors the use of L-band SAR data to improve interferometric coherence. However, the electron content of ionosphere, affecting the propagation of the SAR signal, shows particularly strong spatio-temporal variations near the equator, while the dispersive nature of the ionosphere makes its effect stronger on low-frequencies, such as L-band signals. To tackle this problem, range split-spectrum method can be implemented to compensate the ionospheric phase contribution. Here, we apply this technique for time-series of ALOS-PALSAR data, and propose optimizations for low-coherence areas. To evaluate the efficiency of this method to retrieve subtle deformation rates in equatorial regions, we compute time-series using four ALOS-PALSAR datasets in contexts of low to medium coherence, showing slow deformation rates (mm/yr to cm/yr). The processed tracks are located in Ecuador, Trinidad and Sumatra, and feature 15 to 19 acquisitions including very high, dominating ionospheric noise, corresponding to equivalent displacements of up to 2 m. The correction method performs well and allows to reduce drastically the noise level due to ionosphere, with significant improvement compared with a simple plane fitting method. This is due to frequent highly non-linear patterns of perturbation, characterizing equatorial TEC distribution. We use semivariograms to quantify the uncertainty of the corrected time-series, highlighting its dependence on spatial distance. Thus, using ALOS-PALSAR-like archive, one can expect a detection threshold on the Line-of-Sight velocity ranging between 3 and 6 mm/yr, depending on the spatial wavelength of the signal to be observed. These values are consistent with the accuracy derived from the comparison of velocities between two tracks in their overlapping area. In the case studies that we processed, the time-series corrected from ionosphere allows to retrieve accurately fault creep and volcanic signal but it is still too noisy for retrieving tiny long-wavelength signals such as slow (mm/yr) interseismic strain accumulation.

Analysis of Aquifer‐System Deformation in the Doñana Natural Space (Spain) Using Unsupervised Cloud‐Computed InSAR Data and Wavelet Analysis

AbstractThis work is focused on the Almonte‐Marismas aquifer located within Doñana Natural Space (SW Spain); this aquifer is threatened by droughts, irrigation‐driven groundwater overexploitation, urban use, and the potential reactivation of gas extraction and storage projects. We present ground deformation measurements from Sentinel‐1 Interferometric Synthetic Aperture Radar (InSAR) data ranging from 1 to 2.5 cm, covering ∼2,500 km2 from 2014 to 2020. Detecting spatially distributed ground deformation over this agricultural area is challenging due to the low‐coherence radar signals; consequently, the ground movement results are on the same order of magnitude as the errors associated with InSAR data. We have approached this issue by considering auxiliary information such as groundwater levels, climatic time series, and pumping rates and analyzing their spatio‐temporal connections to ground displacements. We correlate InSAR and hydrogeological information through wavelet analysis, developing a Python package that allows applying the same methodology in other similar study areas. For the first time in the Doñana area, a significant relationship between distances to pumping wells and the displacement extent is detected. Moreover, other subsidence‐related triggering factors are identified, such as the soil moisture balance, clay shrinking‐swelling processes, and creep of geological formations. These results are highly valuable to support aquifer management decisions in the Doñana Natural Space; in this border region, three groundwater bodies were officially declared overexploited in 2020. Our findings provide a ground motion baseline assessment to help differentiate historical variations from any future anthropogenic effect in this complex marsh land ecosystem.

Strategies for improving the communication of satellite-derived InSAR data for geohazards through the analysis of Twitter and online data portals

Abstract. Satellite-based earth observation sensors are increasingly able to monitor geophysical signals related to natural hazards, and many groups are working on rapid data acquisition, processing, and dissemination to data users with a wide range of expertise and goals. A particular challenge in the meaningful dissemination of Interferometric Synthetic Aperture Radar (InSAR) data to non-expert users is its unique differential data structure and sometimes low signal-to-noise ratio. In this study, we evaluate the online dissemination of ground deformation measurements from InSAR through Twitter, alongside the provision of open-access InSAR data from the Centre for Observation and Modelling of Earthquakes, Volcanoes and Tectonics (COMET) Looking Into Continents from Space with Synthetic Aperture Radar (LiCSAR) processing system. Our aim is to evaluate (1) who interacts with disseminated InSAR data, (2) how the data are used, and (3) to discuss strategies for meaningful communication and dissemination of open InSAR data. We found that the InSAR Twitter community was primarily composed of non-scientists (62 %), although this grouping included earth observation experts in applications such as commercial industries. Twitter activity was primarily associated with natural hazard response, specifically following earthquakes and volcanic activity, where users disseminated InSAR measurements of ground deformation, often using wrapped and unwrapped interferograms. For earthquake events, Sentinel-1 data were acquired, processed, and tweeted within 4.7±2.8 d (the shortest was 1 d). Open-access Sentinel-1 data dominated the InSAR tweets and were applied to volcanic and earthquake events in the most engaged-with (retweeted) content. Open-access InSAR data provided by LiCSAR were widely accessed, including automatically processed and tweeted interferograms and interactive event pages revealing ground deformation following earthquake events. The further work required to integrate dissemination of InSAR data into longer-term disaster risk-reduction strategies is highly specific, to both hazard type and international community of practice, as well as to local political setting and civil protection mandates. Notably, communication of uncertainties and processing methodologies are still lacking. We conclude by outlining the future direction of COMET LiCSAR products to maximize their useability.

The 2013–2020 seismic activity at Sabancaya Volcano (Peru): Long lasting unrest and eruption

Sabancaya volcano is the youngest and second most active volcano in Peru. It is part of the Ampato-Sabancaya volcanic complex which sits to the south of the ancient Hualca Hualca volcano and several frequently active faults, thus resulting in complex volcano-tectonic interactions. After 15 years of repose, in 2013, a series of 4 earthquakes with magnitude >4.5 occurred within 24 h, marking the beginning of a new episode of unrest. Several additional swarms of earthquakes occurred in the following years until magmatic eruptive activity started on 6 November 2016. This activity is ongoing as of this writing, with an average of 50 explosions per day. In this study, we present results of multiparametric monitoring of Sabancaya's activity observed during 2013–2020. Seismic data are used to create a one-dimensional seismic velocity model, to catalog, locate, and characterize earthquakes, to detect repeating earthquake families, and to monitor seismic velocity variations by ambient noise cross-correlation. These analyses are complemented by visual and remote sensing observations and ground deformation measurements. All monitored parameters showed significant changes on 6 November 2016, the day of eruption onset, thus dividing the eruptive activity into pre-eruptive and eruptive stages.The unrest is characterized by high levels of seismic activity with hundreds of events detected per day. Volcano-tectonic (VT) earthquakes were dominant during the pre-eruptive period while long-period (LP) events and explosions have been most numerous since the eruption onset. Earthquake locations highlight long-lasting seismogenic zones along multiple previously active regional faults, as well as along newly identified faults. This VT seismicity is mainly distributed in a sector from the northwest to the east of the volcanic complex at distances of up to 30 km from the crater. We focus our analysis on two eruptive episodes: the eruption onset and subsequent crater migration from south to north, and the increase of lava dome extrusion rate in 2019. Both episodes are accompanied by seismic velocity decreases of up to 0.2% and are preceded by a few weeks by bursts of distal VT activity, including numerous repeating earthquakes. These repeated events were located on several remote tectonic faults (5–25 km from the vent). We suggest that these phenomena could be due to the injection of a batch of magma in the deep reservoir and/or conduit, which would generate 1) a pressure wave propagating in the hydrothermal system, triggering the bursts of seismic activity and 2) slow rising of magma by melting old material filling the conduit that eventually produced the eruptive and dome growth acceleration events.

Monitoring of Water and Tillage Soil Erosion in Agricultural Basins, a Comparison of Measurements Acquired by Differential Interferometric Analysis of Sentinel TopSAR Images and a Terrestrial LIDAR System

Agricultural production, the main pillar of food security, is highly dependent on soil quality, and threatened by erosion processes that degrade soil quality. This article is part of a research to verify the usefulness of differential interferometric analysis on TopSAR (Terrain Observation with Progressive Scans SAR, Synthetic Aperture Radar) images to measure water and tillage erosion in small agricultural basins. For this, images from the Sentinel 1 mission are used, analyzing the deformations on the earth’s surface. The purpose of this research is to verify the accuracy of the proposed method by comparing its measures with the ones taken with the gold standard laser terrestrial LIDAR (Light Detection and Ranging) system, as well as to establish a basic step period framework that guarantees an admissible loss of coherence. The results on a pilot plot in El Molar (north of Madrid, Spain) showed that the differences lay within the range of the error associated with the very LIDAR system and showed that coherence losses correspond with the deformations measured. Given the economic and labor advantages of the differential interferometric analysis, this method could be regarded as an excellent alternative to the use of LIDAR in large-scale studies for measuring ground deformation caused by water and tillage erosion.

Remote monitoring of minewater rebound and environmental risk using satellite radar interferometry

The cessation of dewatering following coalfield abandonment results in the rise of minewater, which can create significant changes in the local and regional hydrogeological regime. Monitoring such change is challenging but essential to avoiding detrimental consequences such as groundwater contamination and surface flooding. Inverse modelling methods using satellite radar interferometry (InSAR) have proven capable for retrospectively mapping minewater level changes, however, there is a need for the capability to remotely monitor changes as they occur. In this study, ground deformation measurements obtained from InSAR are used to develop a method to remotely monitor the spatio-temporal rise of minewater, which could be implemented in near real-time. The approach is demonstrated over the Horlivka mining agglomeration, Ukraine, where there is no other feasible approach possible due to a lack of safe ground access. The results were blindly validated against in-situ measurements before being used to forecast the time until minewater will reach the natural water table and Earth's surface. The findings reveal that, as a result of military conflict in Donbas, an environmental catastrophe could occur where potentially radioactive minewater is forecast to reach the natural water table between May and August of 2024.

Field measurement and numerical investigation of artificial ground freezing for the construction of a subway cross passage under groundwater flow

Artificial ground freezing is an environmentally friendly and reliable method to provide temporary ground improvement and groundwater isolation under adverse geological and hydrological conditions. Groundwater flow supplies a substantial source of heat which has an undesirable impact on the closure of frozen soil body. Furthermore, the phenomenon of frost heave and thaw settlement may threaten the surrounding structures. In terms of a ground freezing case in Beijing Metro Line 19, field measurement of ground temperature and deformation were conducted for the prediction of the formation of frozen curtain and the site stability. A hydro-thermal coupling model considering the effect of groundwater flow and phase transition was proposed to quantify the temperature distribution. The numerical model considers the soil particles, crystal ice and liquid water as separate phases, and thermal parameters as variable values related to temperature. Through two basic physical laws and corresponding governing equations, the model captures the relevant coupled relationship in heat and mass transfer. The proposed model is validated by means of the field measured data. Furthermore, the influence of groundwater velocities on the closure of frozen curtain was investigated. The results indicate groundwater flow has an adverse influence on the formation of freezing curtain, even leading to a failure to the closure. The thickness and temperature of frozen curtain with a groundwater velocity of 2 m/d can attain the design values after freezing for 50 days. Moreover, the ground and tunnel deformation in the process of freezing construction were well controlled.

InSAR and numerical modelling for tailings dam monitoring – the Cadia failure case study

In this paper, ground deformation measurements from synthetic aperture radar interferometry (InSAR) and ground-based prism monitoring are compared to finite-element (FE) simulation results for a recent tailings dam collapse. The InSAR monitoring demonstrated the complex spatial and temporal variability of the tailings dam deformation that is not captured in traditional, point-based monitoring approaches. Potentially anomalous deformation behaviour could have been detected from InSAR 1 year preceding the dam failure. A two-dimensional coupled-consolidation analysis was used to predict the dam behaviour during construction stages. The FE modelling results were in broad agreement with the deformation measurement sources, both in terms of magnitude and trends preceding failure. The results, however, deviate significantly immediately before failure, following the construction of the buttresses. Moreover, the FE modelling is sensitive to parameter uncertainties, such as spatial variability of foundation soil properties. The FE results revealed that deformation in both the upstream and downstream parts of the dam are sensitive to changing foundation properties. However, this impact is potentially much more significant on the downstream compared to the upstream parts of the dam, which are dominated by the behaviour of the tailings. This study highlights the efficacy and complementarity of geotechnical and remote sensing techniques for the monitoring of tailings dams.

Yield estimate for the January 2016 DPRK nuclear test based on InSAR and numerical modelling with mechanical heterogeneity

SUMMARY Interferometric Synthetic Aperture Radar (InSAR) provides a way of measuring ground deformation caused by underground nuclear testing via a satellite remote sensing platform. We present a finite-element model (FEM) of the nuclear test site in the Democratic People's Republic of Korea (DPRK, formerly North Korea) constrained by InSAR data. The model simulates elastic deformation caused by an underground explosion and includes high-resolution topography and mechanical heterogeneity. To the best of our knowledge, we are the first to investigate the effect of mechanical heterogeneity (a layered structure in our case) on deformation modelling at the DPRK test site. Sensitivity testing shows that inclusion of topography, mechanical layering and the combined effect results in 48.2, 22.6 and 91.5 per cent increase in the maximum predicted deformation, respectively. This suggests that previous models that do not account for these complexities may overestimate the nuclear yield and/or underestimate the depth of burial (DoB). We performed a parameter search to solve for the cavity size and DoB that best reproduce ground deformation observed in ALOS-2 InSAR data covering the DPRK's 4th nuclear test (DPRK4). Using these best-fit values with a cavity-yield scaling law, we estimate the yield of DPRK4 to be 10–20 kt at a DoB of 542–826 m. Our method provides an estimate that is independent from and complementary to traditional seismic methods, which is useful for a location that is data-limited.

Polarimetric SAR pixel offset tracking for large-gradient landslide displacement mapping

Interferometry and pixel offset tracking (POT) are two major technologies to measure ground deformation with repeat-pass synthetic aperture radar (SAR) observations. The SAR POT method can tolerate a much larger deformation gradient than interferometric SAR (InSAR), and thus it has been widely applied in monitoring earthquakes, volcanoes, landslides, glaciers, etc. However, there is an intrinsic trade-off between matching accuracy and detail preserving in terms of matching window size for the traditional SAR POT method using only single-polarization intensity. Taking advantage of polarimetric scattering properties, this study develops a novel POT method, named PCM-POT, towards mitigating the dependence of offset estimations on matching window size. It employs a generalized likelihood ratio test (GLRT) upon polarimetric covariance matrices (PCMs) to evaluate the similarity of two candidate pixels and estimate high-resolution surface displacements with high accuracy. The effectiveness and advantage of PCM-POT are demonstrated by an experiment of retrieving the displacements of the Slumgullion landslide from a pair of full-polarimetric UAVSAR images. Comparisons with the traditional POT indicate that using polarimetric scattering properties can reduce the uncertainty in offset estimations over heterogeneous sliding areas. Further experiments suggest that PCM-POT is also suitable for dual-polarimetric satellite SAR data.

High resolution C-band SAR backscatter response to peatland water table depth and soil moisture: a laboratory experiment

ABSTRACT Carbon storage and active carbon sequestration within peatlands strongly depend on water table depth and soil moisture availability. With increasing efforts to protect and restore peatland ecosystems, the assessment of their hydrological condition is highly necessary but remains challenging. Synthetic aperture radar (SAR) satellite observations likely offer an efficient way to obtain regular information with complete spatial coverage over northern peatlands. Studies have indicated that both radar backscatter amplitude and phase are sensitive to peatland condition. Very recently, Differential Interferometric Synthetic Aperture Radar (DInSAR) has been reported as being capable of monitoring ground deformation patterns at the millimetre scale, which are a response to peatland hydrological condition. To further investigate the promise of SAR for peatland monitoring, a laboratory-based polarimetric C-band SAR system was used to acquire the dynamic radar behaviour of a 4 m (l) ×1 m (w) × 0.25 m (d) reconstructed peatland. A forced 4-month drought was introduced with very-high-resolution imagery taken every 2 hours, capturing details of the vertical backscatter patterning through the peat at the centimetric scale. The results showed a clear coherent response both in radar backscatter amplitude and phase to change in water table level and soil moisture. Similar responses were seen across all polarizations. Phase demonstrated a coherent and deterministic change across the experiment; the average differential phase increase across all polarizations was 118° for 17 cm of water table drawdown. Interpreted as the physical movement of the surface, this corresponded to 8.3 mm of surface subsidence. Both phase and amplitude changes were near-linear with changes in the water table depth; amplitude showed a correspondingly strong concomitant mean decrease of 7 dB across all polarizations during the experiment. The results demonstrate the close sensitivity of radar backscatter to hydrological patterns in a peatland ecosystem. The phase result, in particular, strongly supports the notion that differential phase from satellites can be utilized to measure ground deformation as a proxy for the hydrological state.

ASSESSING LEVELLING AND DINSAR FOR DEFORMATION MONITORING IN SEISMIC REGION

Abstract. Carrying out monitoring surveys in seismic regions is good practice both for the assessment of land deformation and the evaluation of building structures standing on it. In this work, topographic levelling and DInSAR techniques have been used for displacement measurement. These geomatic techniques are rarely applied in the same context and attempts are made to combine the results obtained for having a complete analysis of the site. The proposed work analyses, compares and discusses topographic levelling and advanced multi-temporal DInSAR techniques used to detect and measure ground deformation when the occurrence of seismic events might have played a role in the displacement. The area of interest had already been under observation through ground-based monitoring surveys, by means of metal bolts attached to façades of buildings detected by topographic level, from 1998 to 2021. The DInSAR analysis was carried out exploiting Sentinel-1A/B data acquired during the period 2014–2021. The goal of the DInSAR processing stage of the procedure is to derive the deformation map of the area of interest from SAR data. A zero date has been set for both survey methods in order to define similar time series for comparison analysis. The results showed that ground displacements measured by levelling and DInSAR have similar trends. On the geomorphological aspect, the same distribution map of terrain subsidence is found in both techniques.

A new approach for 2-D and 3-D precise measurements of ground deformation from optimized registration and correlation of optical images and ICA-based filtering of image geometry artifacts

High resolution satellite images with improved spatial and temporal resolution provide unprecedented opportunities to monitor Earth Surface changes in 2D and 3D due, for example, to earthquakes, sand dune migration, ice flow, or landslides. The volume of imagery available for such measurements is rapidly growing but the exploitation of these data is challenging due to the various sources of geometric distortions of the satellite imagery. Here we propose a new approach to extract high-quality surface displacement in 3D based on the correlation of multi-date and multi-platform high resolution optical imagery. We additionally show that when a large enough volume of data is available, it is possible to separate the deformation signal from the artifacts due to the satellite jitter and misalignment of the CCDs, which, together with topographic artifacts, are the main source of noise in the measurements. Our method makes use of a reference DEM, but the outcome is independent of the characteristics of the chosen DEM. We use the case-example of the ground deformation caused by the Ridgecrest earthquake sequence to assess the performance of our proposed approach. We show that it outperforms the more standard approach which combines 2-D correlation and DEM differencing. With our technique, we were able to generate high quality measurements of coseismic ground displacement with GSD of 2.4 m, and uncertainties at the 90% confidence level on the NS, EW and vertical displacement measurements of 0.6 m, 0.7 m, and 0.6 m respectively.

Terrain deformation measurements from optical satellite imagery: The MPIC-OPT processing services for geohazards monitoring

Measuring terrain deformation over several spatial and temporal scales is relevant for many applications in Earth Sciences (i.e. active faults, volcanoes, landslides or glaciers understanding). The growing volume of freely available data represents nowadays a challenge in terms of storage capacity and computing resources which, together with the complexity of the processing (code parameterization, combination of the image sequences, co-registration of the images) may prevent the exploitation of long time series. We propose here a new version of the Multiple-Pairwise Image Correlation toolbox for processing OPTical images (MPIC-OPT). The toolbox proposes an end-to-end solution to compute the horizontal sub-pixel ground deformation time series from large Sentinel-2 datasets. In addition to time series inversion, several corrections and filtering options are integrated to reduce the noise and improve the accuracy and precision of the measurements. In particular, an automatic jitter correction based on wavelet filtering is proposed. Moreover, the MPIC-OPT service is deployed on the Tier 1.5 High-Performance Computing cluster (e.g. Datacentre/EOST-A2S) of the University of Strasbourg and is accessible on-line through the ESA Geohazards Exploitation Platform (GEP) and the ForM@Ter Solid Earth computing infrastructure with a user-friendly environment to query the satellite data catalogues, parameterize the processing and visualize the outputs. We test the performances of MPIC-OPT on several use cases: the measurement of the co-seismic ground deformation of the 2019 Ridgecrest earthquake sequence (USA), of the rapid motion of the Slumgullion landslide (USA) and of the glaciers of the Mont-Blanc massif (France/Italy). We show that the results of MPIC-OPT are in agreement with in-situ data. The jitter correction significantly improves the precision (RMSEjitter=0.3m vs. RMSEnojitter=0.5m) and the accuracy (RMSEjitter=0.3m vs. RMSEnojitter=1.3m) of the measurement of the co-seismic displacement of the Rigdecrest seismic deformation. We show that the precision and accuracy of the terrain deformation estimation depend mainly on the correlation threshold and the temporal matching range parameters and we quantify and discuss their impacts. This work opens new perspectives to monitor automatically surface displacements/velocities of natural hazards over large scales and large periods of time.

Multi-Temporal SAR Interferometry for Vertical Displacement Monitoring from Space of Tengiz Oil Reservoir Using SENTINEL-1 and COSMO-SKYMED Satellite Missions

This study focused on the quantitative assessment of the vertical displacement velocities retrieved using Sentinel-1 and Cosmo-SkyMed synthetic aperture radar images for the Tengiz oilfield. Tengiz oilfield was selected as a study area because of its historically reported continuous subsidence and limited up-to-date studies during recent years. The small baseline subset time-series technique was used for the interferometric processing of radar images acquired for the period of 2018–2020. The geospatial and statistical analyses allowed to determine the existing hotspots of the subsidence processes induced by oil extraction in the study area. Ground deformation measurements derived from the Sentinel-1 and COSMO-SkyMed satellite missions showed that the Tengiz oilfield continuously subsided during 2018–2020 with the maximum annual vertical displacement velocity around −77.4 mm/y and −71.5 mm/y, respectively. The vertical displacement velocities derived from the Sentinel-1 and the COSMO-SkyMed images showed a good statistical relationship with R2≥0.73 and RMSE ≤3.68 mm. The cumulative vertical displacement derived from both satellites for the most subsiding location also showed a good statistical relationship with R2 equal to 0.97 and RMSE = ± 4.69. The observed relative differences of measurements by both satellites were acceptable to determine the ongoing vertical surface displacement processes in the study area. These studies demonstrated a practical novelty for the petroleum industry in terms of the comparative assessment of surface displacement measurements using time-series of medium-resolution Sentinel-1 and high-resolution COSMO-SkyMed radar images.

The Submarine Boundaries of Mount Etna’s Unstable Southeastern Flank

Coastal and ocean island volcanoes are renowned for having unstable flanks. This can lead to flank deformation on a variety of temporal and spatial scales ranging from slow creep to catastrophic sector collapse. A large section of these unstable flanks is often below sea level, where information on the volcano-tectonic structure and ground deformation is limited. Consequently, kinematic models that attempt to explain measured ground deformation onshore associated with flank instability are poorly constrained in the offshore area. Here, we attempt to determine the locations and the morpho-tectonic structures of the boundaries of the submerged unstable southeastern flank of Mount Etna (Italy). The integration of new marine data (bathymetry, microbathymetry, offshore seismicity, reflection seismic lines) and published marine data (bathymetry, seafloor geodesy, reflection seismic lines) allows identifying the lineament north of Catania Canyon as the southern lateral boundary with a high level of confidence. The northern and the distal (seaward) boundaries are less clear because no microbathymetric or seafloor geodetic data are available. Hypotheses for their locations are presented. Geophysical imaging suggests that the offshore Timpe Fault System is a shallow second-order structure that likely results from extensional deformation within the moving flank. Evidence for active uplift and compression upslope of the amphitheater-shaped depression from seismic data along with subsidence of the onshore Giarre Wedge block observed in ground deformation data leads us to propose that this block is a rotational slump, which moves on top of the large-scale instability. The new shoreline-crossing structural assessment may now inform and improve kinematic models.

Construction of “Space-Sky-Ground” Integrated Collaborative Monitoring Framework for Surface Deformation in Mining Area

Ground deformation measurements in mining areas play a key role in revealing the surface subsidence law, retrieving the subsidence parameters, warning of geological disasters and restoring the surface ecology. With the development of science and technology, there have emerged a great number of monitoring techniques and buildings of diverse protection levels. The diversity of monitoring techniques and the multiplicity of monitoring objects have brought challenges for surface deformation monitoring in the coal industry. Based on the existing deformation monitoring techniques, this paper established a framework of “space-sky-ground” collaborative monitoring system in mining area. We also constructed an AHP-TOPSIS (Analytic Hierarchy Process method- Technique for Order Preference by Similarity to an Ideal Solution) preference model of “space-sky-ground” collaborative monitoring of surface deformation in mining area, and carried out engineering application. Our study shows that the framework of the “space-sky-ground” collaborative monitoring system for surface subsidence in mining areas established in this paper, combined with the AHP-TOPSIS monitoring preference model, which can fully combine the advantages of each monitoring technique, overcome the limitations of a single monitoring technique, comprehensively obtain the surface subsidence data and work out the surface deformation subsidence pattern. This information provides a data and technical support for surface environment management.

The 2021 Mw 7.4 Madoi Earthquake: An Archetype Bilateral Slip‐Pulse Rupture Arrested at a Splay Fault

AbstractWe combine measurements of ground deformation from Synthetic Aperture Radar images, high‐rate Global Navigation Satellite System and tele‐seismic waveforms to study the rupture kinematics of the Madoi Earthquake, which occurred in eastern Tibet on 21 May 2021 and reached a moment magnitude Mw 7.4. The data show nearly pure left‐lateral motion along a 170 km long rupture and a total duration of 36 s. The earthquake initiated near the middle of the main segment and evolved in a bilateral slip pulse rupture which propagated at a sub‐Rayleigh speed of 2.6–2.8 km/s. In our model, slip is concentrated at depth of less than ∼15 km and reaches a maximum of 4.2 m. The rupture arrested ∼10 s after branching on the extensional splay faults at both extremities. The branching onto the splay faults and the eventual arrest of the rupture is used to provide constraints on the fault frictional properties.