Nothing Special   »   [go: up one dir, main page]

You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Sign in to use this feature.

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

Article Types

Countries / Regions

Search Results (414)

Search Parameters:
Keywords = time series InSAR

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
18 pages, 12063 KiB  
Article
Deformation Monitoring and Analysis of Beichuan National Earthquake Ruins Museum Based on Time Series InSAR Processing
by Jing Fan, Weihong Wang, Jialun Cai, Zhouhang Wu, Xiaomeng Wang, Hui Feng, Yitong Yao, Hongyao Xiang and Xinlong Luo
Remote Sens. 2024, 16(22), 4249; https://doi.org/10.3390/rs16224249 - 14 Nov 2024
Viewed by 288
Abstract
Since the Wenchuan earthquake in 2008, Old Beichuan County-town has experienced significant subsidence due to the disruption of the geological environment and the concurrent increase in precipitation. The ongoing land surface deformation poses a threat to the preservation and utilization of the Beichuan [...] Read more.
Since the Wenchuan earthquake in 2008, Old Beichuan County-town has experienced significant subsidence due to the disruption of the geological environment and the concurrent increase in precipitation. The ongoing land surface deformation poses a threat to the preservation and utilization of the Beichuan National Earthquake Ruins Museum (BNERM), as well as to the safety of urban residents’ lives. However, the evolutionary characteristics of surface deformation in these areas remain largely unexplored. Here, we focused on the BNERM control zone and employed the small-baseline subset interferometric synthetic aperture radar (SBAS-InSAR) technique to accurately measure land surface deformation and its spatiotemporal changes. Subsequently, we integrated this data with land cover types and precipitation to investigate the driving factors of deformation. The results indicate a slight overall elevation increase in the study area from June 2015 to May 2023, with deformation rates varying between −35.2 mm/year and 22.9 mm/year. Additionally, four unstable slopes were identified within the BNERM control zone. Our analysis indicates that surface deformation in the study area is closely linked to changes in land cover types and precipitation, exhibiting a seasonal cumulative pattern, and active geological activity may also be a cause of deformation. This study provides invaluable insights into the surface deformation characteristics of the BNERM and can serve as a scientific foundation for the protection of earthquake ruins, risk assessment, early warning, and disaster prevention measures. Full article
(This article belongs to the Section Remote Sensing in Geology, Geomorphology and Hydrology)
Show Figures

Figure 1

Figure 1
<p>Overview of the study area. (<b>a</b>) BNERM Control Zone and SAR satellite imagery coverage. (<b>b</b>) The extent of the study area (from Google Earth on 6 September 2020), where RJP for Renjiaping Earthquake Memorial Museum and a comprehensive Service Area, OBC for Earthquake Ruins Protection Area of Old Beichuan, TJS for Tangjiashan Secondary Hazard Demonstration and Natural Recovery Area. (<b>c</b>) The current status of the BNERM.</p>
Full article ">Figure 2
<p>The main framework and flow chart of the methodology.</p>
Full article ">Figure 3
<p>(<b>a</b>) The location of the reference point. (<b>b</b>) The time series of the reference point. (<b>c</b>) Interference baseline map.</p>
Full article ">Figure 4
<p>(<b>a</b>) The deformation rate in the LOS direction of the study area and (<b>b</b>–<b>d</b>) correspond to the deformation rate maps of RJP, OBC, and TJS, respectively.</p>
Full article ">Figure 5
<p>(<b>a</b>–<b>c</b>) are the selected feature points of RJP, OBC, and TJS. (<b>d</b>,<b>e</b>) are the unstable slope boundaries of the feature points B3 and C3. (<b>f</b>,<b>g</b>) are the current status of the collapse in Jingjiashan and landslides in Wangjiyan.</p>
Full article ">Figure 6
<p>Time series cumulative variogram of feature points. (<b>a</b>–<b>c</b>) correspond to A1, B1, C1, which are the feature points in P1 region; (<b>d</b>–<b>f</b>) correspond to A2, B2, C2, which are the feature points in OBC region; (<b>g</b>–<b>i</b>) correspond to A3, B3, C3, which are the feature points in TJS region.</p>
Full article ">Figure 7
<p>(<b>a</b>) Cumulative displacement map. (<b>b</b>) RMSE distribution of cumulative deformation in the study area.</p>
Full article ">Figure 8
<p>RMSE distribution of cumulative deformation.</p>
Full article ">Figure 9
<p>(<b>a</b>) 2017 Land Use Types in the Study Area. (<b>b</b>) 2023 Land Use Types in the Study Area.</p>
Full article ">Figure 10
<p>(<b>a</b>) Map of land cover type changes in the deformed area between 2017 and 2023. (<b>b</b>) Map of land cover types changes in the deformed area that changed between 2017 and 2023.</p>
Full article ">Figure 11
<p>(<b>a</b>) Surface deformation rate of Renjiaping Earthquake Memorial Museum and comprehensive service area, and RJP stands for Renjiaping Earthquake Memorial Museum and a comprehensive Service Area. (<b>b</b>) The optical images of the area in 2015. (<b>c</b>) The optical images of the area in 2020.</p>
Full article ">Figure 12
<p>Precipitation and time series cumulative deformation trends. (<b>a</b>) is the Renjiaping Earthquake Memorial Museum and comprehensive Service Area. (<b>b</b>) is the Earthquake Ruins Protection Area of Old Beichuan, and (<b>c</b>) is the Tangjiashan Secondary Hazard Demonstration and Natural Recovery Area. (<b>d</b>) is the average of the study area.</p>
Full article ">Figure 13
<p>Distribution of faults in the vicinity of the study area and the seismotectonic context (the yellow ball indicates large and small earthquakes that occurred after the 12 May 2008 earthquake and between 17 June 2015, and the red ball indicates earthquakes that occurred between 18 June 2015 and 31 May 2023).</p>
Full article ">
20 pages, 22822 KiB  
Article
Monitoring Aeolian Erosion from Surface Coal Mines in the Mongolian Gobi Using InSAR Time Series Analysis
by Jungrack Kim, Bayasgalan Amgalan and Amanjol Bulkhbai
Remote Sens. 2024, 16(21), 4111; https://doi.org/10.3390/rs16214111 - 3 Nov 2024
Viewed by 932
Abstract
Surface mining in the southeastern Gobi Desert has significant environmental impacts, primarily due to the creation of large coal piles that are highly susceptible to aeolian processes. Using spaceborne remote sensing and numerical simulations, we investigated erosional processes and their environmental impacts. Our [...] Read more.
Surface mining in the southeastern Gobi Desert has significant environmental impacts, primarily due to the creation of large coal piles that are highly susceptible to aeolian processes. Using spaceborne remote sensing and numerical simulations, we investigated erosional processes and their environmental impacts. Our primary tool was Interferometric Synthetic Aperture Radar (InSAR) data from Sentinel-1 imagery collected between 2017 and 2022. We analyzed these data using phase angle information from the Small Baseline InSAR time series framework. The time series analyses revealed intensive aeolian erosion in the coal piles, represented as thin deformation patterns along the potential pathways of aerodynamic transportation. Further analysis of multispectral data, combined with correlations between wind patterns and trajectory simulations, highlighted the detrimental impact of coal dust on the surrounding environment and the mechanism of aeolian erosion. The lack of mitigation measures, such as water spray, appeared to exacerbate erosion and dust generation. This study demonstrates the feasibility of using publicly available remote sensing data to monitor coal mining activities and their environmental hazards. Our findings contribute to a better understanding of coal dust generation processes in surface mining operations as well as the aeolian erosion mechanism in desert environments. Full article
(This article belongs to the Special Issue Remote Sensing and Geophysics Methods for Geomorphology Research)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>(<b>a</b>) Location of study areas, (<b>b</b>) geographical and topographic context of Ail Bayan and Tavan Tolgoi coal mines, (<b>c</b>) surrounding hydrological contexts (source: <a href="https://eic.mn/" target="_blank">https://eic.mn/</a>, accessed on 1 July 2024), (<b>d</b>) coal production and transportation at Tavan Tolgoi (43.625N, 105.474E) and (<b>e</b>) coal dust generation at Ail Bayan (43.717N, 108.946E) (images were taken in March 2019). Note that the transportation vehicles in Tavan Tolgoi are well-confined so as not to produce coal dust.</p>
Full article ">Figure 2
<p>Acquisition times and connection graphs of employed ascending/descending Sentinel-1 InSAR pairs over Ail Bayan (<b>a</b>,<b>b</b>) and Tavan Tolgoi (<b>c</b>–<b>e</b>). Note that the phase coherences of InSAR pairs are always higher than 0.7. In Ail Bayan, the ascending mode InSAR observations were interpolated into the descending mode time domain for decomposition. Similarly, over Tavan Tolgoi, the ascending and descending mode observations in path 62 were interpolated into the descending mode time domain of path 135. While the thresholds for perpendicular and temporal baselines were set to 150 m and 25 days, respectively, some InSAR pairs exceeding these thresholds were included to enhance interferometric coverage.</p>
Full article ">Figure 3
<p>Processing workflow for InSAR time series data, including integration with other satellite and spatial datasets.</p>
Full article ">Figure 4
<p>(<b>a</b>) Google map, which reveals the details of Ail Bayan, (<b>b</b>) topography presented in Copernicus 30 m DEM, (<b>c</b>) ascending LOS velocity and (<b>d</b>) descending LOS deformation velocity.</p>
Full article ">Figure 5
<p>(<b>a</b>) Decomposed horizontal velocity and (<b>b</b>) vertical velocity in Ail Bayan. The overlaid average wind velocities were extracted using GEE and interpolated to a 1 km resolution from the original 11.132 km ERA5-Land data using kriging.</p>
Full article ">Figure 6
<p>(<b>a</b>) Google map which reveals the details of Tavan Tolgoi, (<b>b</b>) topography presented in Copernicus 30 m DEM, (<b>c</b>) ascending LOS velocity, (<b>d</b>) descending LOS deformation velocity of path 135 coverage, and (<b>e</b>) descending LOS deformation velocity of path 62 coverage.</p>
Full article ">Figure 7
<p>Decomposed velocities in Tavan Tolgoi: (<b>a</b>) horizontal velocity and (<b>b</b>) vertical velocity. Note that the wind directions are similar to those in the Ail Bayan area, blowing from west to east.</p>
Full article ">Figure 8
<p>The behavior of seven RoIs along with mean wind velocities in different modes: (<b>a</b>) ascending mode, (<b>b</b>) descending mode, (<b>c</b>) decomposed horizontal deformation velocities, and (<b>d</b>) decomposed vertical deformation velocities.</p>
Full article ">Figure 9
<p>Correlation maps between InSAR deformation velocities and average wind velocities for corresponding periods in (<b>a</b>) ascending mode, (<b>b</b>) descending mode, (<b>c</b>) horizontal component of decomposed InSAR velocities, and (<b>d</b>) vertical component of decomposed InSAR velocities.</p>
Full article ">Figure 10
<p>Spectral signature analyses using Sentinel-2 time series images on (<b>a</b>) 23 March 2018, (<b>b</b>) 17 April 2018, (<b>c</b>) 2 May 2018, and (<b>d</b>) 22 May 2018. Note that a lower SID value indicates greater spectral similarity. (<b>e</b>) Visual band view of Sentinel-2 image (2 May 2018), (<b>f</b>) spectral signatures in Down1 and Down3/Up3 areas representing coal mine and major FMP aeolian sites.</p>
Full article ">Figure 11
<p>Wind factors influencing coal mine dust generation: (<b>a</b>) monthly wind velocity at an altitude of 10 m, (<b>b</b>) friction velocities for different 10 m wind velocities and roughness lengths, (<b>c</b>) trajectory simulations originating from the coal mine during the sand dust season from March to May 2018, (<b>d</b>) trajectory simulations during the summer season of 2018.</p>
Full article ">Figure 12
<p>Average NMDI maps for (<b>a</b>) Ail Bayan from 24 September 2017 to 14 August 2018, (<b>b</b>) the same region from 23 January 2017 to 9 May 2022, (<b>c</b>) the Tavan Tolgoi region from 11 September 2017 to 10 March 2018, and (<b>d</b>) the same region from 4 January 2017 to 26 May 2022.</p>
Full article ">Figure 13
<p>Environmental consequences of coal mine dust generation: (<b>a</b>) HYSPLIT trajectory simulations originating from the coal mine in 8 March 2018, (<b>b</b>) HYSPLIT trajectory simulations originating from the coal mine in 18 April 2018, using an ensemble HYSPLIT model with 150-h forward trajectory options. (<b>c</b>) Ground photos from Down3 area in Ail Bayan showing contaminated soil and vegetation by blown FMP (images were taken in March 2019).</p>
Full article ">
23 pages, 7056 KiB  
Article
Land Subsidence Predictions Based on a Multi-Component Temporal Convolutional Gated Recurrent Unit Model in Kunming City
by Tao Chen, Di Ning and Yuhang Liu
Appl. Sci. 2024, 14(21), 10021; https://doi.org/10.3390/app142110021 - 2 Nov 2024
Viewed by 529
Abstract
Land subsidence (LS) is a geological hazard driven by both natural conditions and human activities. Traditional LS time-series prediction models often struggle to accurately capture nonlinear data characteristics, leading to suboptimal predictions. To address this issue, this paper introduces a multi-component temporal convolutional [...] Read more.
Land subsidence (LS) is a geological hazard driven by both natural conditions and human activities. Traditional LS time-series prediction models often struggle to accurately capture nonlinear data characteristics, leading to suboptimal predictions. To address this issue, this paper introduces a multi-component temporal convolutional gate recurrent unit (MC-TCGRU) model, which integrates a fully adaptive noise-ensemble empirical-mode decomposition algorithm with a deep neural network to account for the complexity of time-series data. The model was validated using typical InSAR subsidence data from Kunming, analyzing the impact of each component on the prediction performance. A comparative analysis with the TCGRU model and models based on seasonal-trend decomposition using LOESS (STL) and empirical-mode decomposition (EMD) revealed that the MC-TCGRU model significantly enhanced the prediction accuracy by reducing the complexity of the original data. The model achieved R² values of 0.90, 0.93, 0.51, 0.93, and 0.96 across five points, outperforming the compared models. Full article
(This article belongs to the Special Issue Advanced Remote Sensing Technologies and Their Applications)
Show Figures

Figure 1

Figure 1
<p>Location of the study area, with an indication of the administrative boundaries of Kunming City.</p>
Full article ">Figure 2
<p>Structure of the GRU network unit.</p>
Full article ">Figure 3
<p>Structure of the TCN residual block.</p>
Full article ">Figure 4
<p>Structure of the TCGRU model.</p>
Full article ">Figure 5
<p>Structure of the proposed MC-TCGRU model.</p>
Full article ">Figure 6
<p>Deformation velocity map of the study area, with indication of five distinct LS regions and points.</p>
Full article ">Figure 7
<p>CEEMDAN decomposition results for 5 LS points. (<b>a</b>) P1; (<b>b</b>) P2; (<b>c</b>) P3; (<b>d</b>) P4; and (<b>e</b>) P5.</p>
Full article ">Figure 8
<p>Prediction results of each component based on CEEMDAN for point P1; (<b>a</b>) IMF<sub>1</sub> component; (<b>b</b>) IMF<sub>2</sub> component; (<b>c</b>) IMF<sub>3</sub> component; and (<b>d</b>) residual component.</p>
Full article ">Figure 8 Cont.
<p>Prediction results of each component based on CEEMDAN for point P1; (<b>a</b>) IMF<sub>1</sub> component; (<b>b</b>) IMF<sub>2</sub> component; (<b>c</b>) IMF<sub>3</sub> component; and (<b>d</b>) residual component.</p>
Full article ">Figure 9
<p>Point P1 LS time-series prediction results (<b>left</b>) and enlarged view of test set (<b>right</b>).</p>
Full article ">Figure 10
<p>Point P2 LS time-series prediction results (<b>left</b>) and enlarged view of test set (<b>right</b>).</p>
Full article ">Figure 11
<p>Point P3 LS time-series prediction results (<b>left</b>) and enlarged view of test set (<b>right</b>).</p>
Full article ">Figure 12
<p>Point P4 LS time-series prediction results (<b>left</b>) and enlarged view of test set (<b>right</b>).</p>
Full article ">Figure 13
<p>Point P5 LS time-series prediction results (<b>left</b>) and enlarged view of test set (<b>right</b>).</p>
Full article ">Figure 14
<p>Comparison of five AMEs on five LS points. (<b>a</b>) P1. (<b>b</b>) P2. (<b>c</b>) P3. (<b>d</b>) P4. (<b>e</b>) P5.</p>
Full article ">Figure 14 Cont.
<p>Comparison of five AMEs on five LS points. (<b>a</b>) P1. (<b>b</b>) P2. (<b>c</b>) P3. (<b>d</b>) P4. (<b>e</b>) P5.</p>
Full article ">
28 pages, 6037 KiB  
Article
Statistical and Independent Component Analysis of Sentinel-1 InSAR Time Series to Assess Land Subsidence Trends
by Celina Anael Farías, Michelle Lenardón Sánchez, Roberta Bonì and Francesca Cigna
Remote Sens. 2024, 16(21), 4066; https://doi.org/10.3390/rs16214066 - 31 Oct 2024
Viewed by 665
Abstract
Advanced statistics can enable the detailed characterization of ground deformation time series, which is a fundamental step for thoroughly understanding the phenomena of land subsidence and their main drivers. This study presents a novel methodological approach based on pre-existing open-access statistical tools to [...] Read more.
Advanced statistics can enable the detailed characterization of ground deformation time series, which is a fundamental step for thoroughly understanding the phenomena of land subsidence and their main drivers. This study presents a novel methodological approach based on pre-existing open-access statistical tools to exploit satellite differential interferometric synthetic aperture radar (DInSAR) data to investigate land subsidence processes, using European Ground Motion Service (EGMS) Sentinel-1 DInSAR 2018−2022 datasets. The workflow involves the implementation of Persistent Scatterers (PS) time series classification through the PS-Time tool, deformation signal decomposition via independent component analysis (ICA), and drivers’ investigation through spatio-temporal correlation with geospatial and monitoring data. Subsidence time series at the three demonstration sites of Bologna, Ravenna and Carpi (Po Plain, Italy) were classified into linear and nonlinear (quadratic, discontinuous, uncorrelated) categories, and the mixed deformation signal of each PS was decomposed into independent components, allowing the identification of new spatial clusters with linear, accelerating/decelerating, and seasonal trends. The relationship between the different independent components and DInSAR-derived displacement velocity, acceleration, and seasonality was also analyzed via regression analysis. Correlation with geological and groundwater monitoring data supported the investigation of the relationship between the observed deformation and subsidence drivers, such as aquifer resource exploitation, local geological setting, and gas extraction/reinjection. Full article
(This article belongs to the Special Issue Monitoring Geohazard from Synthetic Aperture Radar Interferometry)
Show Figures

Figure 1

Figure 1
<p>The study areas of (1) Ravenna, (2) Bologna, and (3) Carpi–Correggio–Soliera: (<b>a</b>) geographical location in Italy; (<b>b</b>) extent of the European Ground Motion Service (EGMS) Level-3 (L3) and Level-2b (L2b) dataset footprints used for the statistical analysis, overlapped onto the Copernicus Global Digital Elevation Model [<a href="#B44-remotesensing-16-04066" class="html-bibr">44</a>]; and (<b>c</b>) detail of the mean vertical deformation velocity from EGMS L3 datasets, overlapped onto a Google satellite imagery basemap.</p>
Full article ">Figure 2
<p>(<b>a</b>) Mean LOS deformation velocity; (<b>b</b>) acceleration; (<b>c</b>) annual seasonality amplitude; and (<b>d</b>) PS-Time classification maps for Ravenna, overlapped onto Google satellite imagery. The area selected for the following ICA analysis is highlighted on (<b>a</b>). DCV = discontinuous with constant velocity; DVV = discontinuous with variable velocity.</p>
Full article ">Figure 3
<p>(<b>a</b>) Mean LOS deformation velocity; (<b>b</b>) acceleration; (<b>c</b>) annual seasonality amplitude; and (<b>d</b>) PS-Time classification maps for Bologna, overlapped onto Google satellite imagery, with indication of the administrative boundary of the city of Bologna (black polygon). The rectangles (i.e., 1 in (<b>c</b>), and 2 in (<b>b</b>)) indicate the testing areas utilized in the following ICA analysis. DCV = discontinuous with constant velocity; DVV = discontinuous with variable velocity.</p>
Full article ">Figure 4
<p>(<b>a</b>) Mean LOS deformation velocity; (<b>b</b>) acceleration; (<b>c</b>) annual seasonality amplitude; and (<b>d</b>) PS-Time classification maps in the Carpi–Correggio–Soliera area, overlapped onto Google satellite imagery. The area selected for the following ICA analysis is highlighted on (<b>a</b>). DCV = discontinuous with constant velocity; DVV = discontinuous with variable velocity.</p>
Full article ">Figure 5
<p>Independent components identified in Ravenna (Ra) testing area, overlapped onto Google satellite imagery.</p>
Full article ">Figure 6
<p>Independent components identified in Bologna (Bo), covering (<b>a</b>) Area 1, and (<b>b</b>) Area 2, overlapped onto Google satellite imagery.</p>
Full article ">Figure 7
<p>Independent components identified in Soliera (So), overlapped onto Google satellite imagery.</p>
Full article ">Figure 8
<p>Correlation between mean deformation velocity, acceleration and amplitude of the APC, from EGMS products and PS-Time analysis, and the ICs retained for the area of Ravenna. Linear or bilinear fitting (red lines) and R<sup>2</sup> values are shown in the graphs that show the best correlation.</p>
Full article ">Figure 9
<p>Correlation between mean deformation velocity, acceleration, and amplitude of the APC, from EGMS products and PS-Time analysis, and the ICs retained for Area 1 in Bologna. Linear or quadratic fitting (red lines) and R2 values are shown in the graphs that show the best correlation.</p>
Full article ">Figure 10
<p>Correlation between mean deformation velocity, acceleration, and amplitude of the APC, from EGMS products and PS-Time analysis, and the ICs retained for Area 2 in Bologna. Linear or quadratic fitting (red lines) and R<sup>2</sup> values are shown in the graphs that show the best correlation.</p>
Full article ">Figure 11
<p>Correlation between mean deformation velocity, acceleration, and amplitude of the APC, from EGMS products and PS-Time analysis, and the ICs retained for the area of Soliera. Linear or quadratic fitting (red lines) and R<sup>2</sup> values are shown in the graphs that show the best correlation.</p>
Full article ">Figure 12
<p>(<b>a</b>) Geological map with the location of the principal gas fields operating near the coast of Ravenna and two ARPAE groundwater monitoring wells, overlapped onto Google satellite imagery; (<b>b</b>) Comparison between piezometric level variations in ARPAE’s monitoring wells RA49-00 and RA29-00 and the deformation time series of contiguous points; (<b>c</b>) Deformation velocities observed within each lithological unit, expressed in [mm/year]. Gas exploitation data in (<b>a</b>) is made available by the Italian Ministry of Environment and Energy Security [<a href="#B58-remotesensing-16-04066" class="html-bibr">58</a>], while the location of the monitoring wells and the geological layers were downloaded from the MinERva Portal, managed by the Emilia-Romagna Region service [<a href="#B46-remotesensing-16-04066" class="html-bibr">46</a>].</p>
Full article ">Figure 13
<p>(<b>a</b>) Position of ARPAE’S groundwater monitoring wells and the recorded change in piezometric levels (Δ<span class="html-italic">h<sub>i</sub></span>) for the area of Bologna during the studied time period (2018−2022), overlapped onto Google Satellite imagery; (<b>b</b>) Comparison between piezometric level variations in three of the wells and the deformation time series of contiguous PS–DS points; (<b>c</b>) Geological map of Bologna; (<b>d</b>) Deformation velocities observed within each lithological unit, ex-pressed in [mm/year]. Geological layers used in (<b>c</b>) were downloaded from the MinERva Portal, managed by the Emilia-Romagna Region service [<a href="#B46-remotesensing-16-04066" class="html-bibr">46</a>].</p>
Full article ">Figure 14
<p>(<b>a</b>) Geological map of Carpi–Correggio–Soliera subsidence hotspot, overlapped onto Google satellite imagery; (<b>b</b>) Comparison between piezometric level variations in MO10-01 ARPAE’s monitoring well and a deformation time series of a contiguous PS–DS point; (<b>c</b>) Deformation velocities observed within each lithology, expressed in [mm/year]. Geological layers used were downloaded from the MinERva Portal, managed by the Emilia-Romagna Region service [<a href="#B46-remotesensing-16-04066" class="html-bibr">46</a>].</p>
Full article ">Figure A1
<p>Example of a time series classified as “Bilinear” by PS-Time automatic classification algorithm, in the southern area of Soliera.</p>
Full article ">Figure A2
<p>Time series of one of the PS–DS points scored positively for Bo2–IC2 seasonal component.</p>
Full article ">Figure A3
<p>Acceleration variations vs. buffer distances from Angela Angelina reinjection well.</p>
Full article ">
20 pages, 6644 KiB  
Article
Refined Coseismic Slip and Afterslip Distributions of the 2021 Mw 6.1 Yangbi Earthquake Based on GNSS and InSAR Observations
by Zheng Liu, Keliang Zhang, Weijun Gan and Shiming Liang
Remote Sens. 2024, 16(21), 3996; https://doi.org/10.3390/rs16213996 - 28 Oct 2024
Viewed by 608
Abstract
On 21 May 2021, an Mw 6.1 earthquake occurred in Yangbi County, Dali Bai Autonomous Prefecture, Yunnan Province, with the epicenter located in an unmapped blind fault approximately 7 km west of the Weixi-Qiaohou fault (WQF) on the southeastern margin of the Qinghai–Tibetan [...] Read more.
On 21 May 2021, an Mw 6.1 earthquake occurred in Yangbi County, Dali Bai Autonomous Prefecture, Yunnan Province, with the epicenter located in an unmapped blind fault approximately 7 km west of the Weixi-Qiaohou fault (WQF) on the southeastern margin of the Qinghai–Tibetan Plateau. While numerous studies have been conducted to map the coseismic slip distribution by using the Global Navigation Satellite System (GNSS), Interferometric Synthetic Aperture Radar (InSAR) and seismic data as well as their combinations, the understanding of deformation characteristics during the postseismic stage remains limited, mostly due to the long revisiting time interval and large uncertainty of most SAR satellites. In this study, we refined coseismic slip and afterslip distributions with nonlinear inversions for both fault geometry and relaxation time. First, we determined the fault geometry and coseismic slip distribution of this earthquake by joint inversion for coseismic offsets in the line-of-sight (LOS) direction of both Sentinel-1A/B ascending and descending track images and GNSS data. Then, the descending track time series of Sentinel-1 were further fitted using nonlinear least squares to extract the coseismic and postseismic deformations. Finally, we obtained the refined coseismic slip and afterslip distributions and investigated the spatiotemporal evolution of fault slip by comparing the afterslip with aftershocks. The refined coseismic moment magnitude, which was of Mw 6.05, was smaller than Mw 6.1 or larger, which was inferred from our joint inversion and previous studies, indicating a significant reduction in early postseismic deformation. In contrast, the afterslip following the mainshock lasted for about six months and was equivalent to a moment release of an Mw 5.8 earthquake. These findings not only offer a novel approach to extracting postseismic deformation from noisy InSAR time series but also provide valuable insights into fault slip mechanisms associated with the Yangbi earthquake, enhancing our understanding of seismic processes. Full article
(This article belongs to the Special Issue Monitoring Geohazard from Synthetic Aperture Radar Interferometry)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Tectonics and SAR tracks of the Yangbi Region. The black lines indicate active faults near the epicenter. The image frames for Sentinel-1A ascending (blue box), Sentinel-1A descending (red box) and Sentinel-1B ascending (yellow box) are shown. The red star represents the epicenter of the 2021 Yangbi Mw6.1 earthquake. The white box indicates the position of (<b>b</b>). Aftershock data were provided by Tian et al. (2023) [<a href="#B21-remotesensing-16-03996" class="html-bibr">21</a>]. Abbreviations: WQF, Weixi-Qiaohou Fault; LCF, Lancangjiang Fault; RRF, Red River Fault. (<b>b</b>) Aftershocks in 60 days following the mainshock. The color bar depicts the days after the mainshock. Solid blue triangles denote GNSS stations.</p>
Full article ">Figure 2
<p>Interferograms of the 2021 Mw 6.1 Yangbi earthquake. The variation in fringe colors represents changes in LOS direction deformations: (<b>a</b>) on the ascending track between 20 May 2021 and 26 May 2021 and (<b>b</b>) on the descending track between 10 and 22 May 2021.</p>
Full article ">Figure 3
<p>Line-of-sight (LOS) deformation maps of the 2021 Mw 6.1 Yangbi earthquake: (<b>a</b>) on the ascending track between 20 and 26 May 2021 and (<b>b</b>) on the descending track between 10 and 22 May 2021, red represents regions moving toward the satellite, while blue represents regions moving away from the satellite(with gray lines representing coseismic contours at 2 cm intervals).</p>
Full article ">Figure 4
<p>The posterior probability density distribution of the margins of the 2021 Yangbi earthquake fault. The bottom row shows the histograms of the marginal probability density distributions for each parameter. The black waved lines represent probability distribution curves and the red line represents the maximum a posteriori probability solution.</p>
Full article ">Figure 5
<p>(<b>a</b>) Selected iterative search results of the fault geometry. The rectangles in beige denote faults during the search process and the rectangle in green represents the optimal fault geometry. (<b>b</b>) Spatial distribution characteristics of the optimal fault.</p>
Full article ">Figure 6
<p>Coseismic deformation field model and inversion residuals: (<b>a</b>) coseismic deformation; (<b>b</b>) modeled result; and (<b>c</b>) residuals. The arrows represent the horizontal displacements of the GNSS stations, with red representing observed values, blue representing model values and black representing residuals (descending images from 10 to 22 May 2021).</p>
Full article ">Figure 7
<p>Coseismic deformation field model and inversion residuals: (<b>a</b>) coseismic deformation; (<b>b</b>) modeled result; and (<b>c</b>) residuals. The arrows represent the horizontal displacements of the GNSS stations, with red representing observed values, blue representing model values and black representing residuals (ascending images from 20 to 26 May 2021).</p>
Full article ">Figure 8
<p>(<b>a</b>) Spatial distribution of coseismic fault slip for the optimal solution. (<b>b</b>) Coseismic slip distribution on the fault plane. The black arrows represent the direction of slip of the subfaults. The aftershocks were provided by Tian et al. [<a href="#B21-remotesensing-16-03996" class="html-bibr">21</a>].</p>
Full article ">Figure 9
<p>Temporal characteristics of the surface deformation in 6 months after the mainshock of the Yangbi earthquake. The red star denotes the epicenter of the mainshock. Black dots (Selected_Points) denote 2191 points used to refine the coseismic offset and extract postseismic deformation. The black triangles represent eight points (A–H) from the “black dots” and their time series are displayed in <a href="#remotesensing-16-03996-f009" class="html-fig">Figure 9</a>.</p>
Full article ">Figure 10
<p>The postseismic fitting curves of the selected points (points <b>A</b>–<b>H</b>) in <a href="#remotesensing-16-03996-f009" class="html-fig">Figure 9</a>. The blue arrows represent the coseismic observation values of the LOS directions at selected points, and the red arrows represent the coseismic values predicted by the postseismic fitting model at the selected points.</p>
Full article ">Figure 10 Cont.
<p>The postseismic fitting curves of the selected points (points <b>A</b>–<b>H</b>) in <a href="#remotesensing-16-03996-f009" class="html-fig">Figure 9</a>. The blue arrows represent the coseismic observation values of the LOS directions at selected points, and the red arrows represent the coseismic values predicted by the postseismic fitting model at the selected points.</p>
Full article ">Figure 11
<p>Comparison between slip inversions from (<b>top</b>) raw descending interferometric pairs and (<b>bottom</b>) refined coseismic offsets. Subfault dimension: 2 km × 2 km.</p>
Full article ">Figure 12
<p>Characteristics of spatial distribution of fault slip due to postseismic afterslip.</p>
Full article ">
21 pages, 10071 KiB  
Article
Deformation Monitoring and Analysis of Baige Landslide (China) Based on the Fusion Monitoring of Multi-Orbit Time-Series InSAR Technology
by Kai Ye, Zhe Wang, Ting Wang, Ying Luo, Yiming Chen, Jiaqian Zhang and Jialun Cai
Sensors 2024, 24(20), 6760; https://doi.org/10.3390/s24206760 - 21 Oct 2024
Viewed by 857
Abstract
Due to the limitations inherent in SAR satellite imaging modes, utilizing time-series InSAR technology to process single-orbit satellite image data typically only yields one-dimensional deformation information along the LOS direction. This constraint impedes a comprehensive representation of the true surface deformation of landslides. [...] Read more.
Due to the limitations inherent in SAR satellite imaging modes, utilizing time-series InSAR technology to process single-orbit satellite image data typically only yields one-dimensional deformation information along the LOS direction. This constraint impedes a comprehensive representation of the true surface deformation of landslides. Consequently, in this paper, after the SBAS-InSAR and PS-InSAR processing of the 30-view ascending and 30-view descending orbit images of the Sentinel-1A satellite, based on the imaging geometric relationship of the SAR satellite, we propose a novel computational method of fusing ascending and descending orbital LOS-direction time-series deformation to extract the landslide’s downslope direction deformation of landslides. By applying this method to Baige landslide monitoring and integrating it with an improved tangential angle warning criterion, we classified the landslide’s trailing edge into a high-speed, a uniform-speed, and a low-speed deformation region, with deformation magnitudes of 7~8 cm, 5~7 cm, and 3~4 cm, respectively. A comparative analysis with measured data for landslide deformation monitoring revealed that the average root mean square error between the fused landslide’s downslope direction deformation and the measured data was a mere 3.62 mm. This represents a reduction of 56.9% and 57.5% in the average root mean square error compared to the single ascending and descending orbit LOS-direction time-series deformations, respectively, indicating higher monitoring accuracy. Finally, based on the analysis of landslide deformation and its inducing factors derived from the calculated time-series deformation results, it was determined that the precipitation, lithology of the strata, and ongoing geological activity are significant contributors to the sliding of the Baige land-slide. This method offers more comprehensive and accurate surface deformation information for dynamic landslide monitoring, aiding relevant departments in landslide surveillance and management, and providing technical recommendations for the fusion of multi-orbital satellite LOS-direction deformations to accurately reconstruct the true surface deformation of landslides. Full article
Show Figures

Figure 1

Figure 1
<p>InSAR technology imaging mode.</p>
Full article ">Figure 2
<p>Main body of the landslide and research area, where (<b>a</b>) is the main body of the landslide and the extent of the study area and (<b>b</b>) is the location of the Baige landslide.</p>
Full article ">Figure 3
<p>The specific imaging modes of satellites in ascending and descending orbits. Where (<b>a</b>) is the ascending orbit satellite imaging mode and (<b>b</b>) is the descending orbit satellite imaging mode.</p>
Full article ">Figure 4
<p>Fusion extraction flow chart.</p>
Full article ">Figure 5
<p>LOS direction deformation transformation model.</p>
Full article ">Figure 6
<p>Spatial and temporal baseline maps, where (<b>a</b>,<b>b</b>) are time–position plots of the descending and ascending orbit images; (<b>c</b>,<b>d</b>) are time–baseline plots of the descending and ascending orbit images. (These diagrams were drawn using ENVI-SARscape5.6.2 software).</p>
Full article ">Figure 7
<p>Time-series deformation results obtained from the SBAS-InSAR. Where (<b>a</b>) is for the descending orbit dataset and (<b>b</b>) is for the ascending orbit dataset.</p>
Full article ">Figure 8
<p>Location of feature points on the trailing edge of the Baige landslide: (<b>a</b>) is the location of the points in the satellite map and (<b>b</b>) is the location of the points in the descending orbit time-series deformation result.</p>
Full article ">Figure 9
<p>Time-series deformation results for each monitoring point. Where (<b>a</b>–<b>i</b>) are the time-series deformations of the corresponding monitoring points (A–I). The blue line is the descending orbit LOS direction time-series deformation, the orange line is the ascending orbit LOS direction time-series deformation, and the green line is the landslide’s downslope direction time-series deformation.</p>
Full article ">Figure 10
<p>Delineation of landslide areas.</p>
Full article ">Figure 11
<p>Distribution of errors between each time-series deformation result dataset and measured data, where (B), (F) and (H) represent the deformation information of points B, F and H respectively, M-D<sub>S</sub> in green is the error between the measured data and the extracted landslide’s downslope direction deformation, M-D<sub>A</sub> in red is the error between the measured data and the ascending time-series deformation, M-D<sub>D</sub> in blue is the error between the measured data and the descending time-series deformation, and the red dashed line is the time when the first landslide occurred.</p>
Full article ">Figure 12
<p>Relationship between tangential angle and monthly mean deformation and rainfall in different regions around the landslide trailing edge. Where (<b>a</b>–<b>c</b>) are the specific information of Area D1, D2, and D3, respectively, and (<b>d</b>) is the location map of the areas.</p>
Full article ">Figure A1
<p>Results from PS-InSAR processing of descending orbit images: (<b>a</b>,<b>b</b>) are the time-position plot and time-baseline plot of the descending orbit images; (<b>c</b>) is a permanent scatterer point position map of the descending orbit images.</p>
Full article ">Figure A2
<p>Results from PS-InSAR processing of ascending orbit images: (<b>a</b>,<b>b</b>) are the time-position plot and time-baseline plot of the ascending orbit images; (<b>c</b>) is a permanent scatterer point position map of the ascending orbit images.</p>
Full article ">Figure A3
<p>Stratigraphic lithology and topographic map of the Baige landslide.</p>
Full article ">
23 pages, 48646 KiB  
Article
Land Subsidence Detection Using SBAS- and Stacking-InSAR with Zonal Statistics and Topographic Correlations in Lakhra Coal Mines, Pakistan
by Tariq Ashraf, Fang Yin, Lei Liu and Qunjia Zhang
Remote Sens. 2024, 16(20), 3815; https://doi.org/10.3390/rs16203815 - 14 Oct 2024
Viewed by 776
Abstract
The adverse combination of excessive mining practices and the resulting land subsidence is a significant obstacle to the sustainable growth and stability of regions associated with mining activities. The Lakhra coal mines, which contain some of Pakistan’s largest coal deposits, have been overlooked [...] Read more.
The adverse combination of excessive mining practices and the resulting land subsidence is a significant obstacle to the sustainable growth and stability of regions associated with mining activities. The Lakhra coal mines, which contain some of Pakistan’s largest coal deposits, have been overlooked in land subsidence monitoring, indicating a considerable oversight in the region. Subsidence in mining areas can be spotted early when using Interferometric Synthetic Aperture Radar (InSAR), which can precisely monitor ground changes over time. This study is the first to employ the Small Baseline Subset (SBAS)-InSAR and stacking-InSAR techniques to identify land subsidence at the Lakhra coal mines. This research offers critical insights into subsidence mechanisms in the study area, which has never been previously investigated for ground deformation monitoring, by utilizing 150 Sentinel-1A (ascending) images obtained between January 2018 and September 2023. A total of 102 deformation spots were identified using SBAS-InSAR, while stacking-InSAR detected 73 deformation locations. The most extensive cumulative subsidence in the Lakhra coal mine was −114 mm, according to SBAS-InSAR, with a standard deviation of 6.63 mm. In comparison, a subsidence rate of −19 mm/year was reported using stacking-InSAR with a standard deviation of 1.17 mm/year. The rangeland covered 88.8% of the total area and exhibited the most significant deformation values, as determined by stacking and SBAS-InSAR techniques. Linear regression showed that there was not a strong correlation between subsidence and topographic factors. As detected by optical remote sensing data, the subsidence locations were near or above the mines in the research area, indicating that widespread mining in Lakhra coal mines was the cause of subsidence. Our findings suggest that SAR interferometric time series analysis is helpful for proactively identifying and controlling subsidence difficulties in mining regions by closely monitoring activities, hence reducing negative consequences on operations and the environment. Full article
Show Figures

Figure 1

Figure 1
<p>Location of the Study Area Lakhra Mines. (<b>A</b>) Top left corner shapefile of Pakistan and area showing Sindh province where Lakhra mines are located. (<b>B</b>) Area of Interest zoomed view.</p>
Full article ">Figure 2
<p>(<b>A</b>) Digital Elevation Model of the Study Area. (<b>B</b>) Land Cover Map.</p>
Full article ">Figure 3
<p>SBAS- and Stacking-InSAR Workflow.</p>
Full article ">Figure 4
<p>Interferogram Network and Average Spatial Coherence.</p>
Full article ">Figure 5
<p>(<b>A</b>) Stacking-InSAR Results. (<b>B</b>,<b>C</b>) Magnified View and Profile Plot, Respectively.</p>
Full article ">Figure 6
<p>(<b>A</b>,<b>B</b>) Magnified view of (Stacking) deformation in Upper and Lower Lakhra overlaid on satellite imagery, respectively.</p>
Full article ">Figure 7
<p>Enlarged Google Earth Image of the Upper (<b>left side</b>) and Lower Lakhra (<b>right side</b>).</p>
Full article ">Figure 8
<p>SBAS-InSAR Results. (<b>A</b>) Displacement Rates in Lakhra Coal Mines. (<b>B</b>,<b>C</b>) Magnified View and Profile Plot, Respectively.</p>
Full article ">Figure 9
<p>(<b>A</b>,<b>B</b>) Magnified View of SBAS Accumulative Deformation in the Upper and Lower Lakhra Overlaid on Satellite Imagery, Respectively.</p>
Full article ">Figure 10
<p>SBAS-InSAR Time Series (2018–2023).</p>
Full article ">Figure 11
<p>Terrain Factors and SBAS-InSAR Deformation Distribution. (<b>A</b>–<b>C</b>) Connection of aspect, slope, and elevation with land deformation, respectively.</p>
Full article ">Figure 12
<p>SBAS- and Stacking-InSAR Standard Deviation Plot.</p>
Full article ">Figure 13
<p>SBAS-InSAR Time-Series Plot of 10 Locations from Upper (<b>A</b>) and Lower (<b>B</b>) Lakhra Mines.</p>
Full article ">Figure 14
<p>Expansion of Mines and Subsidence in the Lower Lakhra in 2018–2023. (<b>A</b>,<b>B</b>) Google Earth images of mines from 2018–2023 (<b>C</b>,<b>D</b>) subsidence observed during 2018–2023.</p>
Full article ">Figure 15
<p>Shaft Mining in Lakhra Coal Mines. Source: Mineral Transformation Plan Vision 2025, Government of Pakistan.</p>
Full article ">Figure 16
<p>Coal Production in Pakistan. Source: BP Statistical Review 2022 and Pakistan Energy Yearbook 2022.</p>
Full article ">
17 pages, 11779 KiB  
Article
InSAR Analysis of Partially Coherent Targets in a Subsidence Deformation: A Case Study of Maceió
by Ana Cláudia Teixeira, Matus Bakon, Daniele Perissin and Joaquim J. Sousa
Remote Sens. 2024, 16(20), 3806; https://doi.org/10.3390/rs16203806 - 13 Oct 2024
Viewed by 728
Abstract
Since the 1970s, extensive halite extraction in Maceió, Brazil, has resulted in significant geological risks, including ground collapses, sinkholes, and infrastructure damage. These risks became particularly evident in 2018, following an earthquake, which prompted the cessation of mining activities in 2019. This study [...] Read more.
Since the 1970s, extensive halite extraction in Maceió, Brazil, has resulted in significant geological risks, including ground collapses, sinkholes, and infrastructure damage. These risks became particularly evident in 2018, following an earthquake, which prompted the cessation of mining activities in 2019. This study investigates subsidence deformation resulting from these mining operations, focusing on the collapse of Mine 18 on 10 December 2023. We utilized the Quasi-Persistent Scatterer Interferometric Synthetic Aperture Radar (QPS-InSAR) technique to analyze a dataset of 145 Sentinel-1A images acquired between June 2019 and April 2024. Our approach enabled the analysis of cumulative displacement, the loss of amplitude stability, the evolution of amplitude time series, and the amplitude change matrix of targets near Mine 18. The study introduces an innovative QPS-InSAR approach that integrates phase and amplitude information using amplitude time series to assess the lifecycle of radar scattering targets throughout the monitoring period. This method allows for effective change detection following sudden events, enabling the identification of affected areas. Our findings indicate a maximum cumulative displacement of −1750 mm, with significant amplitude changes detected between late November and early December 2023, coinciding with the mine collapse. This research provides a comprehensive assessment of deformation trends and ground stability in the affected mining areas, providing valuable insights for future monitoring and risk mitigation efforts. Full article
(This article belongs to the Section Engineering Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>Overview of the study area. (<b>a</b>) Location of the state of Alagoas; (<b>b</b>) area of interest demarcated by the red polygon; (<b>c</b>) border of the neighborhoods, with the marking of the neighborhoods closest to the mining complex, and location of the mines (yellow triangle). Image background: Bing Aerial Map.</p>
Full article ">Figure 2
<p>Representation of the location and identification of halite wells (indicated by colored triangles) and geological faults (indicated by red lines) as identified in [<a href="#B8-remotesensing-16-03806" class="html-bibr">8</a>,<a href="#B29-remotesensing-16-03806" class="html-bibr">29</a>]. Each triangle is numbered by the identification number associated with the wells. The color and size of the triangles indicate the depth of each well: smaller, darker triangles represent shallower wells, while larger, yellow triangles indicate deeper wells. Mines 26 and 34 are excluded due to unknown depths. Image background: Bing Aerial Map.</p>
Full article ">Figure 3
<p>The spatial and temporal baseline distribution and interferogram pairs used. In the first image (image on the left), all SAR images were spatially co-registered in the reference image acquired on 13 July 2021. In the second image (image on the right), each image is linked to other images throughout the time series having multi-references.</p>
Full article ">Figure 4
<p>LOS direction velocity (mm/year) map for each PS point obtained using the traditional PSI approach. Velocity values range from −50 to 10 mm/year, represented by a color scale from red to blue, where red indicates −50 mm/year and blue indicates 10 mm/year. Image background: Bing Aerial Map.</p>
Full article ">Figure 5
<p>Cumulative displacement (mm) of each PS point obtained using the QPS-InSAR approach from June 2019 to April 2024. Cumulative displacement values range from less than −400 to 10 mm, represented by a color scale from red to blue, where red indicates values less than −400 mm, which can reach up to −1750 mm, and blue indicating 10 mm. Image background: Bing Aerial Map.</p>
Full article ">Figure 6
<p>InSAR cumulative displacement profile from the points within the black box. Segment A is located in the lagoon, segment B marks the boundary between the Muntage and Pinheiro neighborhoods, and segment C is located further from the mining area. Image background: Bing Aerial Map.</p>
Full article ">Figure 7
<p>Time series of the closest point to each mine identified by the number and categorized by deformation class: (<b>a</b>) low, (<b>b</b>,<b>c</b>) medium, (<b>d</b>) high, and (<b>e</b>) very high. Image background: Bing Aerial Map.</p>
Full article ">Figure 8
<p>(<b>a</b>) LOS direction velocity (mm/year) map for each PS point obtained using the traditional PSI approach. (<b>b</b>) Cumulative displacement (mm) of each PS point obtained with the QPS-InSAR approach from June 2019 to April 2024. Image background: Bing Aerial Map.</p>
Full article ">Figure 9
<p>Time off (years) of points in the mine area (<b>a</b>) and with emphasis on points close to Mine 18 (<b>b</b>). Image background: Bing Aerial Map.</p>
Full article ">Figure 10
<p>Amplitude time series (<b>a</b>,<b>c</b>) and amplitude change matrix (<b>b</b>,<b>d</b>) from two points near Mine 18.</p>
Full article ">
18 pages, 4210 KiB  
Article
Quantifying Creep on the Laohushan Fault Using Dense Continuous GNSS
by Wenquan Zhuang, Yuhang Li, Ming Hao, Shangwu Song, Baiyun Liu and Lihong Fan
Remote Sens. 2024, 16(19), 3746; https://doi.org/10.3390/rs16193746 - 9 Oct 2024
Viewed by 524
Abstract
The interseismic behavior of faults (whether they are locked or creeping) and their quantitative kinematic constraints are critical for assessing the seismic hazards of faults and their surrounding areas. Currently, the creep of the eastern segment of the Laohushan Fault in the Haiyuan [...] Read more.
The interseismic behavior of faults (whether they are locked or creeping) and their quantitative kinematic constraints are critical for assessing the seismic hazards of faults and their surrounding areas. Currently, the creep of the eastern segment of the Laohushan Fault in the Haiyuan Fault Zone at the northeastern margin of the Tibetan Plateau, as revealed by InSAR observations, lacks confirmation from other observational methods, particularly high-precision GNSS studies. In this study, we utilized nearly seven years of observation data from a dense GNSS continuous monitoring profile (with a minimum station spacing of 2 km) that crosses the eastern segment of the Laohushan Fault. This dataset was integrated with GNSS data from regional continuous stations, such as those from the Crustal Movement Observation Network of China, and multiple campaign measurements to calculate GNSS baseline change time series across the Laohushan Fault and to obtain a high spatial resolution horizontal crustal velocity field for the region. A comprehensive analysis of this primary dataset indicates that the Laohushan Fault is currently experiencing left-lateral creep, characterized by a partially locked shallow segment and a deeper locked segment. The fault creep is predominantly concentrated in the shallow crustal region, within a depth range of 0–5.7 ± 3.4 km, exhibiting a creep rate of 1.5 ± 0.7 mm/yr. Conversely, at depths of 5.7 ± 3.4 km to 16.8 ± 4.2 km, the fault remains locked, with a loading rate of 3.9 ± 1.1 mm/yr. The shallow creep is primarily confined within 3 km on either side of the fault. Over the nearly seven-year observation period, the creep movement within approximately 5 km of the fault’s near field has shown no significant time-dependent variation, instead demonstrating a steady-state behavior. This steady-state creep appears unaffected by postseismic effects from historical large earthquakes in the adjacent region, although the deeper (far-field) tectonic deformation of the Laohushan Fault may have been influenced by the postseismic effects of the 1920 Haiyuan M8.5 earthquake. Full article
(This article belongs to the Special Issue Advances in Multi-GNSS Technology and Applications)
Show Figures

Figure 1

Figure 1
<p>Regional active tectonics map. (<b>a</b>) Gray lines represent major active faults [<a href="#B30-remotesensing-16-03746" class="html-bibr">30</a>]. The red dashed line area indicates the Pull-Apart Basin [<a href="#B31-remotesensing-16-03746" class="html-bibr">31</a>]. Hollow circles represent historical earthquakes [<a href="#B32-remotesensing-16-03746" class="html-bibr">32</a>]. Blue pentagons mark the distribution of GNSS continuous observation profile sites established by the Second Monitoring Center of the China Earthquake Administration across the Laohushan Fault. Red triangles represent GNSS survey stations from the Fifteen Digital Seismic Network Construction Project. The pink triangles denote observation stations from the Crustal Movement Observation Network of China. The blue triangles indicate stations from the National GNSS Geodetic Control Network of China. Green triangles and green squares represent published GNSS observation results. Major fault segments of the Haiyuan Fault are shown as thick solid lines, with the following abbreviations: JQH F., the Jinqianghe Fault; MMS F., the Maomaoshan Fault; LHS F., the Laohushan Fault; HYW F., the western segment of the Haiyuan Fault; HYM F., the central segment of the Haiyuan Fault; HYE F., the eastern segment of the Haiyuan Fault. (<b>b</b>) Details of the solid black line in (<b>a</b>). (<b>c</b>) The orange rectangle in (<b>c</b>) shows our study area.</p>
Full article ">Figure 2
<p>Time series of GNSS short-baseline horizontal vectors. (<b>a</b>) Eastward horizontal baseline vector; (<b>b</b>) northward horizontal baseline vector.</p>
Full article ">Figure 3
<p>Baseline component variation time series for the Laohushan Fault; (<b>a</b>) time series of baseline changes parallel to the Laohushan Fault, where Vpa represents the variation rate of the baseline component parallel to the fault; (<b>b</b>) distribution of GNSS stations; (<b>c</b>) time series of baseline changes perpendicular to the Laohushan Fault, where Vpe represents the variation rate of the baseline component perpendicular to the fault. Green points represent the baseline time series across the fault. Blue points represent the baseline time series on the same side of the fault. Red solid lines represent the linear rate of best fit from the least-squares method.</p>
Full article ">Figure 4
<p>GNSS horizontal velocity field of the Laohushan Fault and surrounding area (relative to the Eurasian reference frame); (<b>a</b>) the pink box indicates the location of the GNSS profile shown in <a href="#remotesensing-16-03746-f005" class="html-fig">Figure 5</a>; (<b>b</b>) an enlarged view showing the denser GNSS velocity field around the Maomaoshan Fault and Laohushan Fault; (<b>c</b>) the errors associated with common points used in the 200 integration of the published velocity field.</p>
Full article ">Figure 5
<p>(<b>a</b>,<b>b</b>) display the GNSS velocity components parallel to the fault, with pink squares indicating the data and gray squares representing outliers (which were not used in the inversion). (<b>c</b>,<b>d</b>) show the GNSS velocity components perpendicular to the fault, with similar color coding: pink squares for data and gray squares for outliers. The black dashed lines in (<b>a</b>–<b>d</b>) indicate the position of the fault as determined by inversion, with yellow shading representing the error margin. The 0 km mark denotes the actual fault position. The blue solid lines in (<b>a</b>,<b>b</b>) represent the inverse tangent curves calculated from the fault slip model, where SS denotes the slip rate, D1 indicates the locking depth, CC represents the shallow creep rate, and D2 denotes the creep depth. The gray shading in (<b>c</b>,<b>d</b>) represents the differential motion perpendicular to the fault on both sides. The black solid line shows the actual fault position.</p>
Full article ">Figure 6
<p>Leveling vertical land motion rates in the Haiyuan Fault and its surrounding regions. Scale of the arrows indicates velocity magnitude, underlined by color. Positive (negative) values denote crustal uplift (subsidence).</p>
Full article ">Figure 7
<p>(<b>a</b>) Rate fitting of the baseline component parallel to the fault in the near-field across the fault; (<b>b</b>) rate fitting of the baseline component parallel to the fault in the far-field across the fault. The red dashed line represents the least-squares linear fit, and the blue dashed line represents the result of the fitting using Equation (2).</p>
Full article ">Figure 8
<p>GNSS velocity field parallel to the Haiyuan Fault. The colors in this Figure represent the magnitude of the velocity, while the direction of the arrows indicates the velocity direction. The red beach ball represents the focal mechanism solution [<a href="#B57-remotesensing-16-03746" class="html-bibr">57</a>].</p>
Full article ">Figure 9
<p>The distribution of small earthquake activity on the Laohushan Fault: (<b>a</b>) The spatial distribution of small earthquake activity on the Laohushan Fault. (<b>b</b>) Red circles represent the projections of small earthquake activity along the depth of the Laohushan Fault, totaling 36 events. (<b>c</b>) The statistical distribution of earthquakes of magnitude 3 or greater with respect to the fault depth.</p>
Full article ">
20 pages, 7449 KiB  
Article
Study on the Relationship between Groundwater and Land Subsidence in Bangladesh Combining GRACE and InSAR
by Liu Ouyang, Zhifang Zhao, Dingyi Zhou, Jingyao Cao, Jingyi Qin, Yifan Cao and Yang He
Remote Sens. 2024, 16(19), 3715; https://doi.org/10.3390/rs16193715 - 6 Oct 2024
Viewed by 1074
Abstract
Due to a heavy reliance on groundwater, Bangladesh is experiencing a severe decline in groundwater storage, with some areas even facing land subsidence. This study aims to investigate the relationship between groundwater storage changes and land subsidence in Bangladesh, utilizing a combination of [...] Read more.
Due to a heavy reliance on groundwater, Bangladesh is experiencing a severe decline in groundwater storage, with some areas even facing land subsidence. This study aims to investigate the relationship between groundwater storage changes and land subsidence in Bangladesh, utilizing a combination of GRACE and InSAR technologies. To clarify this relationship from a macro perspective, the study employs GRACE data merged with GLDAS to analyze changes in groundwater storage and SBAS-InSAR technology to assess land subsidence. The Dynamic Time Warping (DTW) method calculates the similarity between groundwater storage and land subsidence time series, incorporating precipitation and land cover types into the data analysis. The findings reveal the following: (1) Groundwater storage in Bangladesh is declining at an average rate of −5.55 mm/year, with the most significant declines occurring in Rangpur, Mymensingh, and Rajshahi. Notably, subsidence areas closely match regions with deeper groundwater levels; (2) The similarity coefficient between the time series of groundwater storage and land subsidence changes exceeds 0.85. Additionally, land subsidence in different regions shows an average lagged response of 2 to 6 months to changes in groundwater storage. This study confirms a connection between groundwater dynamics and land subsidence in Bangladesh, providing essential knowledge and theoretical support for further research. Full article
Show Figures

Figure 1

Figure 1
<p>Geographic location of the study area.</p>
Full article ">Figure 2
<p>Workflow of this study.</p>
Full article ">Figure 3
<p>Spatial distribution of various types of data: (<b>a</b>) groundwater trend calculated based on GRACE and GLDAS; (<b>b</b>) changes in total surface water storage from GLDAS (including soil moisture, surface runoff, canopy water, and snow water); (<b>c</b>) average precipitation; (<b>d</b>) annual average groundwater storage based on GRACE and GLDAS.</p>
Full article ">Figure 4
<p>Comparison of relevant data: (<b>a</b>) original GRACE TWS; (<b>b</b>) GLDAS data; (<b>c</b>) changes in groundwater processing from three institutions; (<b>d</b>) average groundwater storage.</p>
Full article ">Figure 5
<p>Spatial distribution of land subsidence in Bangladesh. (<b>I</b>) Rangpur; (<b>II</b>) Mymensingh; (<b>III</b>) Rajshahi; (<b>IV</b>) Dhaka; (<b>V</b>) Chittagong; (<b>VI</b>) Khulna.</p>
Full article ">Figure 6
<p>Comparison of InSAR results with GNSS station data in the LOS direction: (<b>a</b>) the VLKA station; (<b>b</b>) the COML station; (<b>c</b>) the BNGM station; (<b>d</b>) the BNTL station; (<b>e</b>) the DHA2 station.</p>
Full article ">Figure 7
<p>Comparison of groundwater storage and precipitation: (<b>a</b>) precipitation from 2014 to 2023; (<b>b</b>) comparison of GRACE-extracted groundwater storage changes and detrended precipitation.</p>
Full article ">Figure 8
<p>The land subsidence and land cover changes in Mymensingh.</p>
Full article ">Figure 9
<p>Comparison of time series curves of groundwater storage and land subsidence in different regions: (<b>I</b>) Rangpur; (<b>II</b>) Mymensingh; (<b>III</b>) Rajshahi; (<b>IV</b>) Dhaka; (<b>V</b>) Chittagong; (<b>VI</b>) Khulna.</p>
Full article ">
25 pages, 34633 KiB  
Article
Identification of Potential Landslides in the Gaizi Valley Section of the Karakorum Highway Coupled with TS-InSAR and Landslide Susceptibility Analysis
by Kaixiong Lin, Guli Jiapaer, Tao Yu, Liancheng Zhang, Hongwu Liang, Bojian Chen and Tongwei Ju
Remote Sens. 2024, 16(19), 3653; https://doi.org/10.3390/rs16193653 - 30 Sep 2024
Viewed by 993
Abstract
Landslides have become a common global concern because of their widespread nature and destructive power. The Gaizi Valley section of the Karakorum Highway is located in an alpine mountainous area with a high degree of geological structure development, steep terrain, and severe regional [...] Read more.
Landslides have become a common global concern because of their widespread nature and destructive power. The Gaizi Valley section of the Karakorum Highway is located in an alpine mountainous area with a high degree of geological structure development, steep terrain, and severe regional soil erosion, and landslide disasters occur frequently along this section, which severely affects the smooth flow of traffic through the China-Pakistan Economic Corridor (CPEC). In this study, 118 views of Sentinel-1 ascending- and descending-orbit data of this highway section are collected, and two time-series interferometric synthetic aperture radar (TS-InSAR) methods, distributed scatter InSAR (DS-InSAR) and small baseline subset InSAR (SBAS-InSAR), are used to jointly determine the surface deformation in this section and identify unstable slopes from 2021 to 2023. Combining these data with data on sites of historical landslide hazards in this section from 1970 to 2020, we constructed 13 disaster-inducing factors affecting the occurrence of landslides as evaluation indices of susceptibility, carried out an evaluation of regional landslide susceptibility, and identified high-susceptibility unstable slopes (i.e., potential landslides). The results show that DS-InSAR and SBAS-InSAR have good agreement in terms of deformation distribution and deformation magnitude and that compared with single-orbit data, double-track SAR data can better identify unstable slopes in steep mountainous areas, providing a spatial advantage. The landslide susceptibility results show that the area under the curve (AUC) value of the artificial neural network (ANN) model (0.987) is larger than that of the logistic regression (LR) model (0.883) and that the ANN model has a higher classification accuracy than the LR model. A total of 116 unstable slopes were identified in the study, 14 of which were determined to be potential landslides after the landslide susceptibility results were combined with optical images and field surveys. These 14 potential landslides were mapped in detail, and the effects of regional natural disturbances (e.g., snowmelt) and anthropogenic disturbances (e.g., mining projects) on the identification of potential landslides using only SAR data were assessed. The results of this research can be directly applied to landslide hazard mitigation and prevention in the Gaizi Valley section of the Karakorum Highway. In addition, our proposed method can also be used to map potential landslides in other areas with the same complex topography and harsh environment. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Location of the study area; (<b>a</b>,<b>b</b>) show where the study area is located; (<b>c</b>) shows the spatial distribution of historical landslide hazards in the study area; and (<b>d</b>,<b>e</b>) are field photographs of the highway surroundings.</p>
Full article ">Figure 2
<p>Pearson correlation test for landslide environmental factors; where * ** *** represents the level of significance.</p>
Full article ">Figure 3
<p>Landslide susceptibility evaluation indicators: (<b>a</b>) Elevation, (<b>b</b>) Curvature, (<b>c</b>) Slope, (<b>d</b>) TWI, (<b>e</b>) Aspect, (<b>f</b>) Soil type, (<b>g</b>) Land cover, (<b>h</b>) Snowmelt, (<b>i</b>) River Distance, (<b>j</b>) MNDWI, (<b>k</b>) Fault Distance, (<b>l</b>) NDBI, (<b>m</b>) Earthquake Distance.</p>
Full article ">Figure 4
<p>Spatial connectivity map of SBAS-InSAR and DS-InSAR interferometric image pairs, where (<b>a</b>) is the ascending orbit SBAS-InSAR connectivity map; (<b>b</b>) is the descending orbit SBAS-InSAR connectivity map; (<b>c</b>) is the ascending orbit DS-InSAR connectivity map; and (<b>d</b>) is the descending orbit DS-InSAR connectivity map. The dots in the figure indicate the numbers of the interferometric images, and the lines indicate the number of image pairs.</p>
Full article ">Figure 5
<p>Technical process.</p>
Full article ">Figure 6
<p>Monitoring deformation in the Gaizi Valley section of the Karakorum Highway using SBAS-InSAR and DS-InSAR; (<b>A</b>,<b>C</b>) are the results of ascending-orbit SBAS-InSAR and DS-InSAR data; (<b>B</b>,<b>D</b>) are the results of descending-orbit SBAS-InSAR and DS-InSAR data; The black rectangular areas (a)–(i) in each subfigure represent aggregated areas of unstable slopes; and the background of the figure is a hillshaded view showing the mountains.</p>
Full article ">Figure 7
<p>Histograms of DS-InSAR and SBAS-InSAR deformation points; (<b>a</b>) ascending SBAS-InSAR; (<b>b</b>) descending SBAS-InSAR; (<b>c</b>) ascending DS-InSAR; (<b>d</b>) descending DS-InSAR.</p>
Full article ">Figure 8
<p>LR and ANN model ROC curves.</p>
Full article ">Figure 9
<p>Landslide susceptibility map of the ANN model (<b>a</b>) vs. that of the LR model (<b>b</b>).</p>
Full article ">Figure 10
<p>Distribution of potential landslides in the study area.</p>
Full article ">Figure 11
<p>Detailed characterization of potential landslides: (1)–(12) potential landslides (white circled lines); (13)–(14) unstable slopes that are not potential landslides (red circled lines); where (<b>a</b>) demonstrates the spatial relationship between potential landslides and TS-InSAR; (<b>b</b>) Demonstrate the spatial relationship between potential landslides and landslide susceptibility; and (<b>c</b>) is a high-resolution image of potential landslides.</p>
Full article ">Figure 12
<p>Time deformation characteristics of potential landslides at GZHG4: (<b>a</b>) SBAS-InSAR deformation rate for the ascending track; (<b>b</b>) SBAS-InSAR deformation rate for the descending track; (<b>c</b>) DS-InSAR deformation rate for the ascending track; (<b>d</b>) DS-InSAR deformation rate for the descending track; (<b>e</b>) cumulative deformation totals at points P1, P2, and P3 in the ascending orbit direction; (<b>f</b>) cumulative deformation totals at points P1, P2, and P3 in the descending orbit direction.</p>
Full article ">Figure 13
<p>Field validation of potential landslide sites (GZHG5): (<b>a</b>) ascending DS-InSAR deformation distribution; (<b>b</b>) cumulative deformation at points P1–P3; (<b>c</b>) map of the overall profile of the slope; (<b>d</b>) manual excavation at the foot of the slope; (<b>e</b>) accumulation of material washed down the slope.</p>
Full article ">Figure 14
<p>Nonpotential landslides in a melting snow gully: (<b>a</b>) location of the study area where the snowmelt zone is located; (<b>b</b>) high-resolution image of the unstable area; (<b>c</b>) distribution of landslide hazard susceptibility; (<b>d</b>) ascending SBAS-InSAR deformation rate map; (<b>e</b>) descending SBAS-InSAR deformation rate map; (<b>f</b>,<b>g</b>) field site photographs.</p>
Full article ">Figure 15
<p>Landslides induced by engineering facilities at a mine site: (<b>a</b>) location of the engineering works at the mine site; (<b>b</b>) high-resolution imagery; (<b>c</b>) landslide susceptibility distribution; (<b>d</b>) ascending SBAS-InSAR deformation rate maps; (<b>e</b>) field photographs.</p>
Full article ">
20 pages, 22765 KiB  
Article
Landslide Susceptibility Assessment Based on Multisource Remote Sensing Considering Inventory Quality and Modeling
by Zhuoyu Lv, Shanshan Wang, Shuhao Yan, Jianyun Han and Gaoqiang Zhang
Sustainability 2024, 16(19), 8466; https://doi.org/10.3390/su16198466 - 29 Sep 2024
Viewed by 664
Abstract
The completeness of landslide inventories and the selection of evaluation models significantly impact the accuracy of landslide susceptibility assessments. Conventional field geological survey methods and single remote-sensing technology struggle to reliably identify landslides under complex environmental conditions. Moreover, prevalent landslide susceptibility evaluation models [...] Read more.
The completeness of landslide inventories and the selection of evaluation models significantly impact the accuracy of landslide susceptibility assessments. Conventional field geological survey methods and single remote-sensing technology struggle to reliably identify landslides under complex environmental conditions. Moreover, prevalent landslide susceptibility evaluation models are often plagued by issues such as subjectivity and overfitting. Therefore, we investigated the uncertainty in susceptibility modeling from the aspects of landslide inventory quality and model selection. The study focused on Luquan County in Yunnan Province, China. Leveraging multisource remote-sensing technologies, particularly emphasizing optical remote sensing and InSAR time-series deformation detection, the existing historical landslide inventory was refined and updated. This updated inventory was subsequently used to serve as samples. Nine evaluation indicators, encompassing factors such as distance to faults and tributaries, lithology, distance to roads, elevation, slope, terrain undulation, distance to the main streams, and average annual precipitation, were selected on the basis of the collation and organization of regional geological data. The information value and two coupled machine-learning models were formulated to evaluate landslide susceptibility. The evaluation results indicate that the two coupled models are more appropriate for susceptibility modeling than the single information value (IV) model, with the random forest model optimized by genetic algorithm in Group I2 exhibiting higher predictive accuracy (AUC = 0.796). Furthermore, comparative evaluation results reveal that, under equivalent model conditions, the incorporation of a remote-sensing landslide inventory significantly enhances the accuracy of landslide susceptibility assessment results. This study not only investigates the impact of landslide inventories and models on susceptibility outcomes but also validates the feasibility and scientific validity of employing multisource remote-sensing technologies in landslide susceptibility assessment. Full article
Show Figures

Figure 1

Figure 1
<p>Some visual information about the study area. (<b>A</b>) The study area is located in the north-central part of Yunnan Province. (<b>B</b>) A sketch of the regional seismotectonic setting and historical seismicity. (<b>C</b>) The black triangle shows the distribution of landslide samples in the historical inventory, and the lower-right corner shows the distribution of engineered lithologies in the study area.</p>
Full article ">Figure 2
<p>Stacking-InSAR technological process.</p>
Full article ">Figure 3
<p>Acquisition process of landslide inventories via multisource remote-sensing technologies.</p>
Full article ">Figure 4
<p>Updated landslide distribution pattern within the study region (the red triangles represent newly identified landslides via multisource remote-sensing technologies, and black triangles represent previous landslides). (A) New landslide-clustering zone along the Pudu River and (B, C) new landslide-clustering zones along the Jinsha River.</p>
Full article ">Figure 5
<p>Typical signs of surface deformation. (<b>A</b>,<b>B</b>) InSAR deformation rate map (ascending) of the new landslide-clustering zone; (<b>a</b>,<b>b</b>) interferometric loop characterization of typical landslide InSAR deformation-monitoring results; and (<b>c</b>,<b>d</b>) signs of deformation damage evident in optical images of landslides.</p>
Full article ">Figure 6
<p>Landslide susceptibility zoning map of Group I1 (<b>a(I)</b>: IV; <b>b(I)</b>: GLR; and <b>c(I)</b>: GWIV).</p>
Full article ">Figure 7
<p>Landslide susceptibility zoning map of Group I2 (<b>a(II)</b>: IV; <b>b(II)</b>: GLR; and <b>c(II)</b>: GWIV).</p>
Full article ">Figure 8
<p>ROC curve.</p>
Full article ">Figure 9
<p>Difference in landslide probability after updating the landslide catalog in the study area. (A) Changes in susceptibility along the Pudu River; (B, C) Changes in susceptibility along the Jinsha River.</p>
Full article ">Figure 10
<p>Feature importance visualization.</p>
Full article ">Figure 11
<p>The spatial distribution of population density in 2020 (Data from WorldPop [<a href="#B25-sustainability-16-08466" class="html-bibr">25</a>]).</p>
Full article ">
26 pages, 29457 KiB  
Article
SSBAS-InSAR: A Spatially Constrained Small Baseline Subset InSAR Technique for Refined Time-Series Deformation Monitoring
by Zhigang Yu, Guanghui Zhang, Guoman Huang, Chunquan Cheng, Zhuopu Zhang and Chenxi Zhang
Remote Sens. 2024, 16(18), 3515; https://doi.org/10.3390/rs16183515 - 22 Sep 2024
Viewed by 1262
Abstract
SBAS-InSAR technology is effective in obtaining surface deformation information and is widely used in monitoring landslides and mining subsidence. However, SBAS-InSAR technology is susceptible to various errors, including atmospheric, orbital, and phase unwrapping errors. These multiple errors pose significant challenges to precise deformation [...] Read more.
SBAS-InSAR technology is effective in obtaining surface deformation information and is widely used in monitoring landslides and mining subsidence. However, SBAS-InSAR technology is susceptible to various errors, including atmospheric, orbital, and phase unwrapping errors. These multiple errors pose significant challenges to precise deformation monitoring over large areas. This paper examines the spatial characteristics of these errors and introduces a spatially constrained SBAS-InSAR method, termed SSBAS-InSAR, which enhances the accuracy of wide-area surface deformation monitoring. The method employs multiple stable ground points to create a control network that limits the propagation of multiple types of errors in the interferometric unwrapped data, thereby reducing the impact of long-wavelength signals on local deformation measurements. The proposed method was applied to Sentinel-1 data from parts of Jining, China. The results indicate that, compared to the traditional SBAS-InSAR method, the SSBAS-InSAR method significantly reduced phase closure errors, deformation rate standard deviations, and phase residues, improved temporal coherence, and provided a clearer representation of deformation in time-series curves. This is crucial for studying surface deformation trends and patterns and for preventing related disasters. Full article
Show Figures

Figure 1

Figure 1
<p>Simulation of phase unwrapping errors. (<b>a</b>) Shows a sample of the simulated interference pattern. (<b>b</b>) Displays a sample of the phase unwrapping result. (<b>c</b>) Illustrates the average absolute phase difference for each block relative to the P1 block, based on 1200 simulations.</p>
Full article ">Figure 2
<p>Flowchart of SSBAS-InSAR method.</p>
Full article ">Figure 3
<p>Schematic diagram of phase correction in a local area.</p>
Full article ">Figure 4
<p>Overview of the study area. (<b>a</b>) Map of the study area. (<b>b</b>) Land cover status, derived from Esri Land Cover 2021 data [<a href="#B36-remotesensing-16-03515" class="html-bibr">36</a>]. (<b>c</b>) Surface elevation of the study area.</p>
Full article ">Figure 5
<p>Spatiotemporal baseline diagram. The red dots represent SAR imagery. Thanks to the excellent orbital control capability of the Sentinel-1 satellite [<a href="#B38-remotesensing-16-03515" class="html-bibr">38</a>], the 200 m spatial baseline threshold did not exclude any image pairs, ensuring redundancy in the baseline connections and avoiding the introduction of extremely low-coherence data. A total of 369 interferometric pairs were formed.</p>
Full article ">Figure 6
<p>Data preparation for point selection. (<b>a</b>) Stacking-InSAR deformation rate results. (<b>b</b>) Average coherence graph.</p>
Full article ">Figure 7
<p>Control point distribution map. The base map shows the Stacking-InSAR deformation rate, with regions of average coherence greater than 0.8. The red rectangle highlights the area within the study region where the absolute value of the Stacking-InSAR deformation rate exceeds 12 cm/year, and the average coherence remains above 0.8.</p>
Full article ">Figure 8
<p>Examples of multi-control point phase constraints. (<b>a</b>,<b>d</b>,<b>g</b>,<b>j</b>,<b>m</b>,<b>p</b>) The original unwrapped images; (<b>b</b>,<b>e</b>,<b>h</b>,<b>k</b>,<b>n</b>,<b>q</b>) the unwrapped images after applying multi-control point constraints; (<b>c</b>,<b>f</b>,<b>i</b>,<b>l</b>,<b>o</b>,<b>r</b>) the differences between the original unwrapped images and those corrected with multi-control point constraints. The area marked by the red rectangle in (<b>l</b>) indicates where the absolute value of the phase exceeds 20 radians.</p>
Full article ">Figure 9
<p>Deformation rates obtained from time-series processing. (<b>a</b>) Deformation rate derived using the SBAS-InSAR method. (<b>b</b>) Deformation rate derived using the SSBAS-InSAR method.</p>
Full article ">Figure 10
<p>Comparison of triangular phase closure. (<b>a</b>,<b>c</b>) The non-zero triangular phase closure numbers of each point of the conventional method and the proposed method, respectively. (<b>b</b>,<b>d</b>) The histograms of the non-zero triangular phase closure numbers of the conventional method and the proposed method, respectively. The unit interval width is 1.</p>
Full article ">Figure 11
<p>Comparison of temporal coherence. (<b>a</b>,<b>c</b>) The temporal coherence for each point in the conventional SBAS-InSAR method and the proposed SSBAS-InSAR method. (<b>b</b>,<b>d</b>) The corresponding histograms, with a unit interval width of 0.01.</p>
Full article ">Figure 12
<p>Comparison of the standard deviation of velocity. (<b>a</b>,<b>c</b>) The standard deviation of velocity for each point using the conventional method and the proposed SSBAS-InSAR method, respectively. (<b>b</b>,<b>d</b>) The histograms of the standard deviation of velocity for the conventional method and the proposed method, respectively. The unit interval width is 0.05 mm/year.</p>
Full article ">Figure 13
<p>Comparison of phase residue. (<b>a</b>,<b>c</b>) The phase residue for each point using the conventional method and the proposed SSBAS-InSAR method, respectively. The red dotted line divides the image into eight regions. (<b>b</b>,<b>d</b>) The histograms of phase residue for the conventional method and the proposed method, respectively. The unit interval width is 0.01 m.</p>
Full article ">Figure 14
<p>Comparison of phase residue for each patch. (<b>a</b>–<b>h</b>) The paches numbered 1 to 8 in <a href="#remotesensing-16-03515-f013" class="html-fig">Figure 13</a> (<b>a</b>,<b>c</b>), respectively. The red dotted line in each figure represents y = x, while the green dotted line indicates the linear regression line.</p>
Full article ">Figure 15
<p>Deformation results for selected areas in Jining. (<b>a</b>) Deformation rate map generated by the SBAS-InSAR method. (<b>b</b>) Deformation rate map generated by the SSBAS-InSAR method. The DEM base map used is GLO-30.</p>
Full article ">Figure 16
<p>Optical images of the locations of selected points. (<b>a</b>) Optical image of the area where point P1 is located. (<b>b</b>) Optical image of the area where point P2 is located. (<b>c</b>) Optical image of the area where point P3 is located. (<b>d</b>) Optical image of the area where the SDJX point is located.</p>
Full article ">Figure 17
<p>Comparison of deformation and temperature at point P1. (<b>a</b>,<b>c</b>,<b>e</b>,<b>g</b>) Comparison of time-series deformation and temperature curves obtained through the SBAS-InSAR, SBAS-InSAR+VEC, SBAS-InSAR+DTLF, and SSBAS-InSAR methods, respectively. (<b>b</b>,<b>d</b>,<b>f</b>,<b>h</b>) The correlation coefficient diagrams between deformation and temperature for the four methods, respectively. The red line indicates the linear regression line.</p>
Full article ">Figure 18
<p>Deformation curve at point P2. (<b>a</b>–<b>d</b>) The deformation curves at point P3 obtained by the SBAS-InSAR, SBAS-InSAR+VEC, SBAS-InSAR+DTLF, and SSBAS-InSAR methods, respectively.</p>
Full article ">Figure 19
<p>Deformation status at point P3. (<b>a</b>) The deformation rate result obtained by the SSBAS-InSAR method. (<b>b</b>) The time-series deformation curves at point P3 obtained by the SBAS-InSAR, SBAS-InSAR+VEC, SBAS-InSAR+DTLF, and SSBAS-InSAR methods.</p>
Full article ">Figure 20
<p>Location and deformation of GNSS stations. (<b>a</b>) Distribution of GNSS stations in the study area and surrounding regions. (<b>b</b>) Three-dimensional deformation of the SDJX point as observed by GNSS.</p>
Full article ">Figure 21
<p>Comparison of GNSS deformation results and time-series InSAR results.</p>
Full article ">Figure A1
<p>Schematic diagram of triangle phase closure.</p>
Full article ">
17 pages, 12325 KiB  
Article
Development and Comparison of InSAR-Based Land Subsidence Prediction Models
by Lianjing Zheng, Qing Wang, Chen Cao, Bo Shan, Tie Jin, Kuanxing Zhu and Zongzheng Li
Remote Sens. 2024, 16(17), 3345; https://doi.org/10.3390/rs16173345 - 9 Sep 2024
Viewed by 732
Abstract
Land subsidence caused by human engineering activities is a serious problem worldwide. We selected Qian’an County as the study area to explore the evolution of land subsidence and predict its deformation trend. This study utilized synthetic aperture radar interferometry (InSAR) technology to process [...] Read more.
Land subsidence caused by human engineering activities is a serious problem worldwide. We selected Qian’an County as the study area to explore the evolution of land subsidence and predict its deformation trend. This study utilized synthetic aperture radar interferometry (InSAR) technology to process 64 Sentinel-1 data covering the area, and high-precision and high-resolution surface deformation data from January 2017 to December 2021 were obtained to analyze the deformation characteristics and evolution of land subsidence. Then, land subsidence was predicted using the intelligence neural network theory, machine learning methods, time-series prediction models, dynamic data processing techniques, and engineering geology of ground subsidence. This study developed three time-series prediction models: a support vector regression (SVR), a Holt Exponential Smoothing (Holt) model, and multi-layer perceptron (MLP) models. A time-series prediction analysis was conducted using the surface deformation data of the subsidence funnel area of Zhouzi Village, Qian’an County. In addition, the advantages and disadvantages of the three models were compared and analyzed. The results show that the three developed time-series data prediction models can effectively capture the time-series-related characteristics of surface deformation in the study area. The SVR and Holt models are suitable for analyzing fewer external interference factors and shorter periods, while the MLP model has high accuracy and universality, making it suitable for predicting both short-term and long-term surface deformation. Ultimately, our results are valuable for further research on land subsidence prediction. Full article
(This article belongs to the Topic Environmental Geology and Engineering)
Show Figures

Figure 1

Figure 1
<p>Location and data: (<b>a</b>,<b>b</b>) location of the study area; (<b>c</b>) coverage of Sentinel-1 data.</p>
Full article ">Figure 2
<p>Structure of multilayer perceptron.</p>
Full article ">Figure 3
<p>The deformation rate and profile position of the land subsidence in Zhouzi Village.</p>
Full article ">Figure 4
<p>Time series cumulative deformation of Zhouzi Village land subsidence.</p>
Full article ">Figure 5
<p>Temporal deformation characteristics of typical profiles (<b>a</b>) A-A′ coherent points deformation velocity; (<b>b</b>) B-B′ coherent points deformation velocity; (<b>c</b>) accumulated deformation of selected coherent points.</p>
Full article ">Figure 6
<p>Time series deformation prediction results based on three models.</p>
Full article ">Figure 7
<p>Distribution of land subsidence prediction errors. The blue points represent the minimum error, and the red points represent the maximum error.</p>
Full article ">Figure 8
<p>Deformation error of typical coherent point prediction.</p>
Full article ">
25 pages, 8750 KiB  
Article
Liaohe Oilfield Reservoir Parameters Inversion Based on Composite Dislocation Model Utilizing Two-Dimensional Time-Series InSAR Observations
by Hang Jiang, Rui Zhang, Bo Zhang, Kangyi Chen, Anmengyun Liu, Ting Wang, Bing Yu and Lin Deng
Remote Sens. 2024, 16(17), 3314; https://doi.org/10.3390/rs16173314 - 6 Sep 2024
Viewed by 657
Abstract
To address the industry’s demand for sustainable oilfield development and safe production, it is crucial to enhance the scientific rigor and accuracy of monitoring ground stability and reservoir parameter inversion. For the above purposes, this paper proposes a technical solution that employs two-dimensional [...] Read more.
To address the industry’s demand for sustainable oilfield development and safe production, it is crucial to enhance the scientific rigor and accuracy of monitoring ground stability and reservoir parameter inversion. For the above purposes, this paper proposes a technical solution that employs two-dimensional time-series ground deformation monitoring based on ascending and descending Interferometric Synthetic Aperture Radar (InSAR) technique first, and the composite dislocation model (CDM) is utilized to achieve high-precision reservoir parameter inversion. To validate the feasibility of this method, the Liaohe Oilfield is selected as a typical study area, and the Sentinel-1 ascending and descending Synthetic Aperture Radar (SAR) images obtained from January 2020 to December 2023 are utilized to acquire the ground deformation in various line of sight (LOS) directions based on Multitemporal Interferometric Synthetic Aperture Radar (MT-InSAR). Subsequently, by integrating the ascending and descending MT-InSAR observations, we solved for two-dimensional ground deformation, deriving a time series of vertical and east-west deformations. Furthermore, reservoir parameter inversion and modeling in the subsidence trough area were conducted using the CDM and nonlinear Bayesian inversion method. The experimental results indicate the presence of uneven subsidence troughs in the Shuguang and Huanxiling oilfields within the study area, with a continuous subsidence trend observed in recent years. Among them, the subsidence of the Shuguang oilfield is more significant and shows prominent characteristics of single-source center subsidence accompanied by centripetal horizontal displacement, the maximum vertical subsidence rate reaches 221 mm/yr, and the maximum eastward and westward deformation is more than 90 mm/yr. Supported by the two-dimensional deformation field, we conducted a comparative analysis between the Mogi, Ellipsoidal, and Okada models in terms of reservoir parameter inversion, model fitting efficacy, and residual distribution. The results confirmed that the CDM offers the best adaptability and highest accuracy in reservoir parameter inversion. The proposed technical methods and experimental results can provide valuable references for scientific planning and production safety assurance in related oilfields. Full article
(This article belongs to the Special Issue Synthetic Aperture Radar Interferometry Symposium 2024)
Show Figures

Figure 1

Figure 1
<p>The overall technical flow char.</p>
Full article ">Figure 2
<p>Diagram of vertical and horizontal deformation and slope of the subsidence trough.</p>
Full article ">Figure 3
<p>Simplified schematics of two<span class="html-italic">-</span>dimensional time series. The blue circles represent ascending and descending SAR data at time <span class="html-italic">t<sub>i</sub></span>. The horizontal solid line <span class="html-italic">I<sub>i</sub></span> between two points indicates the interferogram, while Δ<span class="html-italic">t<sub>i</sub></span> denotes the time interval between adjacent images.</p>
Full article ">Figure 4
<p>Diagram of the composite dislocation model (CDM).</p>
Full article ">Figure 5
<p>Study area and image coverage. The red and blue boxes indicate the coverage areas of the ascending and descending SAR data, respectively.</p>
Full article ">Figure 6
<p>The spatio-temporal baseline of the (<b>a</b>) ascending and (<b>b</b>) descending interferometric pairs.</p>
Full article ">Figure 7
<p>LOS deformation velocities for ascending (<b>a</b>) and descending (<b>b</b>) datasets. The red box highlights areas of significant subsidence, and the red star denotes the location of the chosen reference point.</p>
Full article ">Figure 8
<p>Profile deformation distribution. (<b>a</b>,<b>c</b>) represent the deformation velocity results for the Shuguang and Huanxiling oilfields from the ascending and descending datasets, including the selected profile lines and feature points. (<b>b</b>,<b>d</b>) represent the same for the deformation velocity results and profile lines for the ascending and descending datasets. (<b>e</b>–<b>h</b>) depict the deformation distribution along the profile lines AA’, BB’, CC’, and DD’ for the ascending and descending datasets.</p>
Full article ">Figure 9
<p>Time-series deformation for feature Points P1 (<b>a</b>) and P2 (<b>b</b>) in the ascending and descending datasets.</p>
Full article ">Figure 10
<p>2D deformation velocity maps. (<b>a</b>,<b>b</b>) represent the vertical and horizontal deformation velocities for the Shuguang oilfield area. (<b>c</b>,<b>d</b>) show the vertical and horizontal deformation velocities for the Huanxiling oilfield area. Positive and negative values for vertical deformation velocities indicate uplift and subsidence, respectively, while positive and negative values for horizontal deformation velocities represent eastward and westward deformation, respectively.</p>
Full article ">Figure 11
<p>Vertical and horizontal deformation distributions along the L1 profile line for the Shuguang oilfield (<b>a</b>) and the L2 profile line for the Huanxiling oilfield (<b>b</b>). (<b>c</b>,<b>d</b>) show the 3D effects of vertical deformation velocities for the Shuguang and Huanxiling oilfield areas, respectively.</p>
Full article ">Figure 12
<p>Vertical (<b>a</b>) and horizontal (<b>b</b>) time-series deformation characteristics from January 2021 to December 2021.</p>
Full article ">Figure 13
<p>3D results of reservoir parameter inversion for the Shuguang oilfield using the CDM.</p>
Full article ">Figure 14
<p>The observed deformation field in the (<b>a</b>) vertical and (<b>d</b>) horizontal directions. The modeled deformation field from the CDM parameter inversion in the (<b>b</b>) vertical and (<b>e</b>) horizontal directions. The residuals in the (<b>c</b>) vertical and (<b>f</b>) horizontal directions.</p>
Full article ">Figure 15
<p>The observed deformation field in the (<b>a</b>) vertical and (<b>d</b>) horizontal directions. The modeled deformation field from the Mogi parameter inversion in the (<b>b</b>) vertical and (<b>e</b>) horizontal directions. The residuals in the (<b>c</b>) vertical and (<b>f</b>) horizontal directions.</p>
Full article ">Figure 16
<p>The observed deformation field in the (<b>a</b>) vertical and (<b>d</b>) horizontal directions. The modeled deformation field from the Ellipsoidal parameter inversion in the (<b>b</b>) vertical and (<b>e</b>) horizontal directions; the residuals in the (<b>c</b>) vertical and (<b>f</b>) horizontal directions.</p>
Full article ">Figure 17
<p>The observed deformation field in the (<b>a</b>) vertical and (<b>d</b>) horizontal directions. The modeled deformation field from the Okada parameter inversion in the (<b>b</b>) vertical and (<b>e</b>) horizontal directions. The residuals in the (<b>c</b>) vertical and (<b>f</b>) horizontal directions.</p>
Full article ">Figure 18
<p>Vertical (<b>a</b>) and horizontal (<b>b</b>) observed deformation and the L1 profile line distribution. Fitting results of the four inversion models compared to the observed deformation along the vertical (<b>c</b>) and horizontal (<b>d</b>) profile lines.</p>
Full article ">Figure 19
<p>Histograms and statistics of residuals for the four models in the vertical (<b>a</b>) and horizontal (<b>b</b>) directions. The horizontal axis represents residual values in mm/yr, while the vertical axis indicates the number of residuals.</p>
Full article ">
Back to TopTop