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 (431)

Search Parameters:
Keywords = time series InSAR

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
19 pages, 2621 KiB  
Article
Multi-Scale Debris Flow Warning Technology Combining GNSS and InSAR Technology
by Xiang Zhao, Linju He, Hai Li, Ling He and Shuaihong Liu
Water 2025, 17(4), 577; https://doi.org/10.3390/w17040577 - 17 Feb 2025
Viewed by 122
Abstract
The dynamic loads of fluid impact and static loads, such as the gravity of a rock mass during the formation of debris flows, exhibit a coupled effect of mutual influence. Under this coupling effect, surface monitoring points in disaster areas experience displacement. However, [...] Read more.
The dynamic loads of fluid impact and static loads, such as the gravity of a rock mass during the formation of debris flows, exhibit a coupled effect of mutual influence. Under this coupling effect, surface monitoring points in disaster areas experience displacement. However, existing methods do not consider the dynamic–static coupling effects of debris flows on the surface. Instead, they rely on GNSS or InSAR technology for dynamic or static single-scale monitoring, leading to high Mean Absolute Percentage Error (MAPE) values and low warning accuracy. To address these limitations and improve debris flow warning accuracy, a multi-scale warning method was proposed based on Global Navigation Satellite System (GNSS) and Synthetic Aperture Radar Interferometry (InSAR) technology. GNSS technology was utilized to correct coordinate errors at monitoring points, thereby enhancing the accuracy of monitoring data. Surface deformation images were generated using InSAR and Small Baseline Subset (SBAS) technology, with time series calculations applied to obtain multi-scale deformation data of the surface in debris flow disaster areas. A debris flow disaster morphology classification model was developed using a support vector mechanism. The actual types of debris flow disasters were employed as training labels. Digital Elevation Model (DEM) files were utilized to extract datasets, including plane curvature, profile curvature, slope, and elevation of the monitoring area, which were then input into the training model for classification training. The model outputted the classification results of the hidden danger areas of debris flow disasters. Finally, the dynamic and static coupling variables of surface deformation were decomposed into valley-type internal factors (rock mass static load) and slope-type triggering factors (fluid impact dynamic load) using the moving average method. Time series prediction models for the variable of the dynamic–static coupling effects on surface deformation were constructed using polynomial regression and particle swarm optimization (PSO)–support vector regression (SVR) algorithms, achieving multi-scale early warning of debris flows. The experimental results showed that the error between the predicted surface deformation results using this method and the actual values is less than 5 mm. The predicted MAPE value reached 6.622%, the RMSE value reached 8.462 mm, the overall warning accuracy reached 85.9%, and the warning time was under 30 ms, indicating that the proposed method delivered high warning accuracy and real-time warning. Full article
(This article belongs to the Special Issue Flowing Mechanism of Debris Flow and Engineering Mitigation)
Show Figures

Figure 1

Figure 1
<p>Annual surface deformation rate chart.</p>
Full article ">Figure 2
<p>Classification model of hidden danger area form.</p>
Full article ">Figure 3
<p>Prediction process of PSO-SVR combination model.</p>
Full article ">Figure 4
<p>Gaussian distribution for solving the optimal estimate.</p>
Full article ">Figure 5
<p>Flowchart of vegetation coverage data collection.</p>
Full article ">Figure 6
<p>Time series curve of surface deformation under external factors [<a href="#B5-water-17-00577" class="html-bibr">5</a>,<a href="#B6-water-17-00577" class="html-bibr">6</a>].</p>
Full article ">Figure 7
<p>Absolute error curve of time series prediction using different prediction methods.</p>
Full article ">Figure 8
<p>Effect of debris flow disaster risk warning [<a href="#B5-water-17-00577" class="html-bibr">5</a>,<a href="#B6-water-17-00577" class="html-bibr">6</a>].</p>
Full article ">Figure 9
<p>Real-time warning test.</p>
Full article ">
22 pages, 24221 KiB  
Article
SBAS-InSAR Monitoring of Landslides and Glaciers Along the Karakoram Highway Between China and Pakistan
by Basit Ali Khan, Chaoying Zhao, Najeebullah Kakar and Xuerong Chen
Remote Sens. 2025, 17(4), 605; https://doi.org/10.3390/rs17040605 - 10 Feb 2025
Viewed by 446
Abstract
Global assessments of landslide impact on critical communication infrastructure have become urgent because of rising occurrences related to human activities and climate change. The landslide and glacial slide susceptibility along the Karakoram Highway poses a significant threat to the infrastructure ecosystem, local communities, [...] Read more.
Global assessments of landslide impact on critical communication infrastructure have become urgent because of rising occurrences related to human activities and climate change. The landslide and glacial slide susceptibility along the Karakoram Highway poses a significant threat to the infrastructure ecosystem, local communities, and the critical China–Pakistan Economic Corridor. This research paper utilized the Small Baseline Subset InSAR technique to monitor the deformation patterns over the past 5 years, yielding high-resolution insights into the terrain instability in this geologically active region. The SBAS time series results reveal that the substantial cumulative deformation in our study area ranges from 203 mm to −486 mm, with annual deformation rates spanning from 62 mm/year to −104 mm/year. Notably, the deformation that occurred is mainly concentrated in the northern section of our study area. The slope’s aspect is responsible for the maximum deformed material flow towards the Karakoram Highway via steep slopes, lost glacial formations, and the climate variations that cause the instability of the terrain. The given pattern suggests that the northern area of the Karakoram Highway is exposed to a greater risk from the combined influence of glacial slides, landslides, and climatic shifts, which call for the increased monitoring of the Karakoram Highway. The SBAS-InSAR method is first-rate in deformation monitoring, and it provides a scientific basis for developing real-time landslide monitoring systems. The line of sight limitations and the complexity and imprecision of weather-induced signal degradation should be balanced through additional data sources, such as field surveys to conduct large slide and glacial slide susceptibility evaluations. These research results support proactive hazard mitigation and infrastructure planning along the China–Pakistan Economic Corridor by incorporating SBAS-InSAR monitoring into the original planning. The country’s trade policymakers and national level engineers can enhance transport resilience, efficiently manage the landslide and glacial slide risks, and guarantee safer infrastructure along this strategic trade route. Full article
(This article belongs to the Section Engineering Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>The study area location map (<b>a</b>) geographic extent of study area between boundaries of China and Pakistan (<b>b</b>) google earth location map.</p>
Full article ">Figure 2
<p>Flowchart of the methodology used in this research.</p>
Full article ">Figure 3
<p>Comparison of the interferograms (<b>a</b>) before and (<b>b</b>) after the atmospheric correction.</p>
Full article ">Figure 4
<p>Deformation time series maps from SBAS technique (2019–2023).</p>
Full article ">Figure 5
<p>Annual deformation rate map along the KKH.</p>
Full article ">Figure 6
<p>Deformation causing landsliding (<b>a</b>) was detected by the SBAS, and (<b>b</b>) Google Earth Maps was used to demonstrate landsliding and geomorphological characteristics along sections of the KKH.</p>
Full article ">Figure 7
<p>Deformation and geomorphological characteristics (<b>a</b>) SBAS deformation map (<b>b</b>) Google Earth map along sections of the KKH.</p>
Full article ">Figure 8
<p>Deformation and geomorphological characteristics (<b>a</b>) landsliding occurred due to the deformation and (<b>b</b>) Google Earth map of the deformed area along sections of the KKH.</p>
Full article ">Figure 9
<p>(<b>a</b>) Glacial deformation and (<b>b</b>) Google Maps image of the glacial sliding area along sections of the KKH.</p>
Full article ">Figure 10
<p>Five-year SBAS deformation map along the KKH.</p>
Full article ">Figure 11
<p>Three example points in red color and KKH in black colored line for the demonstration of SBAS-InSAR time series for the past 5 years.</p>
Full article ">Figure 12
<p>SBAS-InSAR time series deformation for three example points (2019–2023).</p>
Full article ">
22 pages, 15578 KiB  
Article
Analysis of Ground Subsidence Evolution Characteristics and Attribution Along the Beijing–Xiong’an Intercity Railway with Time-Series InSAR and Explainable Machine-Learning Technique
by Xin Liu, Huili Gong, Chaofan Zhou, Beibei Chen, Yanmin Su, Jiajun Zhu and Wei Lu
Land 2025, 14(2), 364; https://doi.org/10.3390/land14020364 - 10 Feb 2025
Viewed by 278
Abstract
The long-term overextraction of groundwater in the Beijing–Tianjin–Hebei region has led to the formation of the world’s largest groundwater depression cone and the most extensive land subsidence zone, posing a potential threat to the operational safety of high-speed railways in the region. As [...] Read more.
The long-term overextraction of groundwater in the Beijing–Tianjin–Hebei region has led to the formation of the world’s largest groundwater depression cone and the most extensive land subsidence zone, posing a potential threat to the operational safety of high-speed railways in the region. As a critical transportation hub connecting Beijing and the Xiong’an New Area, the Beijing–Xiong’an Intercity Railway traverses geologically complex areas with significant ground subsidence issues. Monitoring and analyzing the causes of land subsidence along the railway are essential for ensuring its safe operation. Using Sentinel-1A radar imagery, this study applies PS-InSAR technology to extract the spatiotemporal evolution characteristics of ground subsidence along the railway from 2016 to 2022. By employing a buffer zone analysis and profile analysis, the subsidence patterns at different stages (pre-construction, construction, and operation) are revealed, identifying the major subsidence cones along the Yongding River, Yongqing, Daying, and Shengfang regions, and their impacts on the railway. Furthermore, the XGBoost model and SHAP method are used to quantify the primary influencing factors of land subsidence. The results show that changes in confined water levels are the most significant factor, contributing 34.5%, with strong interactions observed between the compressible layer thickness and confined water levels. The subsidence gradient analysis indicates that the overall subsidence gradient along the Beijing–Xiong’an Intercity Railway currently meets safety standards. This study provides scientific evidence for risk prevention and the control of land subsidence along the railway and holds significant implications for ensuring the safety of high-speed rail operations. Full article
(This article belongs to the Special Issue Assessing Land Subsidence Using Remote Sensing Data)
Show Figures

Figure 1

Figure 1
<p>Study area and the extent of SAR imagery.</p>
Full article ">Figure 2
<p>Ps-InSAR technology roadmap.</p>
Full article ">Figure 3
<p>XGBoost technology roadmap.</p>
Full article ">Figure 4
<p>Time-series evolution trend of subsidence in the background area along the Beijing–Xiong’an Intercity Railway.</p>
Full article ">Figure 5
<p>Accuracy verification of InSAR monitoring subsidence results with leveling precision.</p>
Full article ">Figure 6
<p>Annual average subsidence rates in the buffer zone along the Beijing–Xiong’an Intercity Railway across different periods.</p>
Full article ">Figure 7
<p>Relationship between annual average subsidence rates along the profile of the Beijing–Xiong’an Intercity Railway and the distribution of subsidence funnels.</p>
Full article ">Figure 8
<p>Annual average subsidence rates of the Beijing–Xiong’an Intercity Railway cross-section during different periods (The dashed lines of different colors represent areas where the sedimentation rate varies greatly).</p>
Full article ">Figure 9
<p>Cumulative subsidence at stations along the Beijing–Xiong’an Intercity Railway.</p>
Full article ">Figure 10
<p>Variations in annual average subsidence rates of stations along the Beijing–Xiong’an Intercity Railway across different periods.</p>
Full article ">Figure 11
<p>Distribution of monitoring points on both sides of stations along the Beijing–Xiong’an Intercity Railway.</p>
Full article ">Figure 12
<p>Subsidence rate differences between the eastern and western sides of Xiong’an Station and Bazhou North Station.</p>
Full article ">Figure 13
<p>Subsidence rate differences between the eastern and western sides of Gu’an East Station, Daxing Airport Station, and Beijing Daxing Station.</p>
Full article ">Figure 14
<p>Accuracy validation based on the Extreme Gradient Boosting (XGBoost) model.</p>
Full article ">Figure 15
<p>Distribution of the importance of subsidence-influencing factors for the Beijing–Xiong’an Intercity Railway.</p>
Full article ">Figure 16
<p>SHAP interpretability analysis of overall subsidence characteristics.</p>
Full article ">Figure 17
<p>SHAP interpretability analysis of factor interactions.</p>
Full article ">Figure 18
<p>Variations in subsidence gradient along the Beijing–Xiong’an Intercity Railway.</p>
Full article ">
16 pages, 6946 KiB  
Article
Earthquake Damage Susceptibility Analysis in Barapani Shear Zone Using InSAR, Geological, and Geophysical Data
by Gopal Sharma, M. Somorjit Singh, Karan Nayak, Pritom Pran Dutta, K. K. Sarma and S. P. Aggarwal
Geosciences 2025, 15(2), 45; https://doi.org/10.3390/geosciences15020045 - 1 Feb 2025
Viewed by 652
Abstract
The identification of areas that are susceptible to damage due to earthquakes is of utmost importance in tectonically active regions like Northeast India. This may provide valuable inputs for seismic hazard analysis; however, it poses significant challenges. The present study emphasized the integration [...] Read more.
The identification of areas that are susceptible to damage due to earthquakes is of utmost importance in tectonically active regions like Northeast India. This may provide valuable inputs for seismic hazard analysis; however, it poses significant challenges. The present study emphasized the integration of Interferometric Synthetic Aperture Radar (InSAR) deformation rates with conventional geological and geophysical data to investigate earthquake damage susceptibility in the Barapani Shear Zone (BSZ) region of Northeast India. We used MintPy v1.5.1 (Miami INsar Timeseries software in PYthon) on the OpenSARLab platform to derive time series deformation using the Small Baseline Subset (SBAS) technique. We integrated geology, geomorphology, gravity, magnetic field, lineament density, slope, and historical earthquake records with InSAR deformation rates to derive earthquake damage susceptibility using the weighted overlay analysis technique. InSAR time series analysis revealed distinct patterns of ground deformation across the Barapani Shear Zone, with higher rates in the northern part and lower rates in the southern part. The deformation values ranged from 6 mm/yr to about 18 mm/yr in BSZ. Earthquake damage susceptibility mapping identified areas that are prone to damage in the event of earthquakes. The analysis indicated that about 46.4%, 51.2%, and 2.4% of the area were low, medium, and high-susceptibility zones for earthquake damage zone. The InSAR velocity rates were validated with Global Positioning System (GPS) velocity in the region, which indicated a good correlation (R2 = 0.921; ANOVA p-value = 0.515). Additionally, a field survey in the region suggested evidence of intense deformation in the highly susceptible earthquake damage zone. This integrated approach enhances our scientific understanding of regional tectonic dynamics, mitigating earthquake risks and enhancing community resilience. Full article
(This article belongs to the Special Issue Earthquake Hazard Modelling)
Show Figures

Figure 1

Figure 1
<p>(<b>A</b>) Northeastern India showing study area and GPS station locations (green triangles) used for InSAR data validation. ML: Meghalaya, AS: Assam, NL: Nagaland, MN: Manipur, TR: Tripura, MZ: Mizoram, AR: Arunachal Pradesh, SK: Sikkim. (<b>B</b>) Lineaments extracted from LISS-IV Indian satellite image (background image) in study region within 10 km buffer from Barapani Shear Zone (BSZ). The pie chart to the right represents the orientation of lineaments, whereas the yellow dots are some of the important settlements in the study region.</p>
Full article ">Figure 2
<p>Workflow adopted in present study to derive earthquake damage susceptibility map (MintPy workflow adopted from [<a href="#B28-geosciences-15-00045" class="html-bibr">28</a>]). The color legend indicates tools utilized for executing selected task.</p>
Full article ">Figure 3
<p>InSAR (SBAS)-based deformation time series of BSZ for 2017–2024 derived using MintPy approach. (<b>left</b>) Cumulative displacement in BSZ (in millimeters) during 2017–2024. (<b>right</b>) (<b>A</b>–<b>D</b>) Indicates a few selected locations (in right image) and their corresponding deformation profiles over time within BSZ (average annual velocity). Each dot represents the deformation value over time in centimeters, represented by the slope of a line.</p>
Full article ">Figure 4
<p>Deformation rates (2017–2024) across BSZ at locations X–X′, Y–Y′, Z–Z′, and P–P′ marked in <a href="#geosciences-15-00045-f003" class="html-fig">Figure 3</a> (left).</p>
Full article ">Figure 5
<p>Parameters derived in the present study: (<b>A</b>) InSAR velocity, (<b>B</b>) slope, (<b>C</b>) past earthquake (green dots are past earthquakes for the duration 2012–2023), (<b>D</b>) lineament density. The values corresponding to low, medium, and high classes are provided in <a href="#geosciences-15-00045-t001" class="html-table">Table 1</a>.</p>
Full article ">Figure 6
<p>Parameters from geological and geophysical data products obtained from the Geological Survey of India (GSI): (<b>A</b>) geomorphology, (<b>B</b>) geology, (<b>C</b>) Bouger anomaly, (<b>D</b>) magnetic field. The values corresponding to low, medium, and high classes are provided in <a href="#geosciences-15-00045-t001" class="html-table">Table 1</a>.</p>
Full article ">Figure 7
<p>Earthquake damage susceptibility map derived from a combination of InSAR and geological and geophysical parameters. (<b>a</b>,<b>b</b>) Two locations (among many) of field surveys showing corresponding field photographs on the right of the map.</p>
Full article ">Figure 8
<p>Linear regression and square of the correlation coefficient between GPS Line-of-sight (GPS LOS) and InSAR velocities near a few selected GPS station locations in the northeastern region of India. The <span class="html-italic">x</span>-axis represents the observed velocity (mm/yr) from InSAR measurements, whereas the <span class="html-italic">y</span>-axis represents the GPS Line-of-Sight velocity (mm/yr) computed using GPS velocities.</p>
Full article ">Figure 9
<p>Analysis of Variance (ANOVA) test comparing GPS and InSAR velocities data. The analysis demonstrates no statistically significant differences (<span class="html-italic">p</span> &gt; 0.05) between the GPS and InSAR measurement data.</p>
Full article ">
20 pages, 60234 KiB  
Article
Combining InSAR and Time-Series Clustering to Reveal Deformation Patterns of the Heifangtai Loess Terrace
by Hao Xu, Bao Shu, Qin Zhang, Guohua Xiong and Li Wang
Remote Sens. 2025, 17(3), 429; https://doi.org/10.3390/rs17030429 - 27 Jan 2025
Viewed by 501
Abstract
The Heifangtai Loess terrace in northwest China is frequently affected by landslides due to hydrological factors, establishing it as a significant research area for loess landslides. Advanced time-series InSAR technology facilitates the retrieval of surface deformation information, thereby aiding in the monitoring of [...] Read more.
The Heifangtai Loess terrace in northwest China is frequently affected by landslides due to hydrological factors, establishing it as a significant research area for loess landslides. Advanced time-series InSAR technology facilitates the retrieval of surface deformation information, thereby aiding in the monitoring of landslide deformation status. However, existing methods that analyze deformation patterns do not fully exploit the displacement time series derived from InSAR, which hampers the exploration of potentially coexisting deformation patterns within the area. This study integrates InSAR with time-series clustering methods to reveal the surface deformation patterns and their spatial distribution characteristics in Heifangtai. Initially, utilizing the Sentinel-1 ascending and descending SAR data stack from January 2020 to June 2023, we optimize the interferometric phase based on distributed scatterer characteristics to reduce noise levels and obtain higher spatial density of measurement points. Subsequently, by combining the differential interferometric datasets from both ascending and descending orbits, the multidimensional small baseline subsets technique is employed to calculate the two-dimensional deformation time series. Finally, time-series clustering methods are utilized to extract the deformation patterns present and their spatial distribution from all measurement point time series. The results indicate that the deformation of the Heifangtai is primarily distributed around the surrounding area of the platform, with subsidence deformation being more intense than horizontal deformation. The entire terrace exhibits five deformation patterns: eastward subsidence, westward subsidence, vertical subsidence, westward movement, and eastward movement. The spatial distribution of these patterns suggests that the areas beneath the platform, namely Yanguoxia Town and Dangchuan Village, may be more susceptible to landslide threats in the future. Furthermore, wavelet analysis reveals the response relationship between rainfall and various deformation patterns, further enhancing the interpretability of these patterns. These findings hold significant implications for subsequent landslide monitoring, early warning, and risk prevention. Full article
Show Figures

Figure 1

Figure 1
<p>Flowchart for revealing surface deformation patterns based on InSAR and time series clustering.</p>
Full article ">Figure 2
<p>Schematic diagram of MSBAS. The ascending and descending SAR acquisitions are represented by orange upward triangles and bright blue downward triangles, respectively. The horizontal solid line between the two triangles represents the differential interferogram, labeled as <math display="inline"><semantics> <mrow> <mi>D</mi> <mi>i</mi> <mi>n</mi> <msub> <mi>t</mi> <mi>i</mi> </msub> </mrow> </semantics></math>. The vertical dashed line divides the consecutive acquisitions into time intervals <math display="inline"><semantics> <mrow> <mo>Δ</mo> <msub> <mi>t</mi> <mi>i</mi> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 3
<p>Study area and SAR datasets coverage. (<b>a</b>) Location of the study area indicated by the red polygon. The blue and pink rectangles show the coverage of the ascending and descending Sentinel-1 data. (<b>b</b>) Heifangtai Google satellite image and the distribution of historical landslides (revised after [<a href="#B17-remotesensing-17-00429" class="html-bibr">17</a>]).</p>
Full article ">Figure 4
<p>Average coherence of ascending and descending SAR datasets. (<b>a</b>) Coherent regular grid estimation of ascending SAR datasets. (<b>b</b>) Coherent regular grid estimation of descending SAR datasets. (<b>c</b>) Adaptive coherence estimation of ascending SAR datasets based on distributed scatterers. (<b>d</b>) Adaptive coherence estimation of descending SAR datasets based on distributed scatterers.</p>
Full article ">Figure 5
<p>Spatio-temporal baseline distribution of interferogram. (<b>a</b>) Ascending. (<b>b</b>) Descending.</p>
Full article ">Figure 6
<p>Mean velocity maps of Heifangtai terrace along LOS directions. (<b>a</b>) Ascending. (<b>b</b>) Descending. (<b>c</b>) Ascending based on DS. (<b>d</b>) Descending based on DS.</p>
Full article ">Figure 7
<p>Standard deviation maps of LOS deformation velocity. (<b>a</b>) Ascending. (<b>b</b>) Descending. (<b>c</b>) Ascending based on DS. (<b>d</b>) Descending based on DS.</p>
Full article ">Figure 8
<p>Two-dimensional deformation velocity maps of Heifangtai terrace. (<b>a</b>) Horizontal east–west deformation velocity map. (<b>b</b>) Vertical deformation velocity map.</p>
Full article ">Figure 9
<p>The PCs obtained from PCA and the results of the knee point detection. (<b>a</b>) First four PCs of horizontal displacement. (<b>b</b>) First four PCs of vertical displacement. (<b>c</b>) Knee point detection results of horizontal displacement. (<b>d</b>) Knee point detection results of vertical displacement.</p>
Full article ">Figure 10
<p>Clustering time series of horizontal and vertical displacement datasets. (<b>a</b>–<b>c</b>) Three clustering time series of horizontal east–west displacement. (<b>d</b>,<b>e</b>) Two clustering time series of vertical displacement.</p>
Full article ">Figure 11
<p>Geographic distribution map of deformation patterns of Heifangtai terrace combined with horizontal and vertical displacement clustering results.</p>
Full article ">Figure 12
<p>Cross wavelet transform and wavelet coherence of the fluctuation component of vertical deformation displacement and rainfall.</p>
Full article ">Figure 13
<p>Cross wavelet transform and wavelet coherence of the fluctuation component of horizontal deformation displacement and rainfall.</p>
Full article ">
12 pages, 20046 KiB  
Communication
Time-Series Change Detection Using KOMPSAT-5 Data with Statistical Homogeneous Pixel Selection Algorithm
by Mirza Muhammad Waqar, Heein Yang, Rahmi Sukmawati, Sung-Ho Chae and Kwan-Young Oh
Sensors 2025, 25(2), 583; https://doi.org/10.3390/s25020583 - 20 Jan 2025
Viewed by 563
Abstract
For change detection in synthetic aperture radar (SAR) imagery, amplitude change detection (ACD) and coherent change detection (CCD) are widely employed. However, time-series SAR data often contain noise and variability introduced by system and environmental factors, requiring mitigation. Additionally, the stability of SAR [...] Read more.
For change detection in synthetic aperture radar (SAR) imagery, amplitude change detection (ACD) and coherent change detection (CCD) are widely employed. However, time-series SAR data often contain noise and variability introduced by system and environmental factors, requiring mitigation. Additionally, the stability of SAR signals is preserved when calibration accounts for temporal and environmental variations. Although ACD and CCD techniques can detect changes, spatial variability outside the primary target area introduces complexity into the analysis. This study presents a robust change detection methodology designed to identify urban changes using KOMPSAT-5 time-series data. A comprehensive preprocessing framework—including coregistration, radiometric terrain correction, normalization, and speckle filtering—was implemented to ensure data consistency and accuracy. Statistical homogeneous pixels (SHPs) were extracted to identify stable targets, and coherence-based analysis was employed to quantify temporal decorrelation and detect changes. Adaptive thresholding and morphological operations refined the detected changes, while small-segment removal mitigated noise effects. Experimental results demonstrated high reliability, with an overall accuracy of 92%, validated using confusion matrix analysis. The methodology effectively identified urban changes, highlighting the potential of KOMPSAT-5 data for post-disaster monitoring and urban change detection. Future improvements are suggested, focusing on the stability of InSAR orbits to further enhance detection precision. The findings underscore the potential for broader applications of the developed SAR time-series change detection technology, promoting increased utilization of KOMPSAT SAR data for both domestic and international research and monitoring initiatives. Full article
(This article belongs to the Section Remote Sensors)
Show Figures

Figure 1

Figure 1
<p>Location of study site along with KOMPSAT-5 time-series image footprints.</p>
Full article ">Figure 2
<p>Dataset for change detection analysis: (<b>a</b>) optical footprint of the study site (source: Google Earth; image acquisition date: 4 October 2024), (<b>b</b>) KOMPSAT-5 SAR imagery of the study site, and (<b>c</b>) generated ground truth for accuracy assessment.</p>
Full article ">Figure 3
<p>Preprocessing of KOMPSAT-5 time-series images: (<b>a</b>) KOMPSAT-5 time-series stack, (<b>b</b>) KOMPSAT-5 radiometric terrain-corrected time-series stack, (<b>c</b>) KOMPSAT-5 radiometric terrain-corrected normalized time-series stack.</p>
Full article ">Figure 4
<p>Detailed methodological framework adopted for change detection using KOMPSAT-5 images.</p>
Full article ">Figure 5
<p>Experimental results to obtain appropriate statistical homogeneous pixels (SHPs).</p>
Full article ">Figure 6
<p>Statistical homogeneous pixels (SHPs) selection: (<b>a</b>) KOMPSAT-5 image, (<b>b</b>) resultant SHPs over urban segments. The San Francisco port area, highlighted within the red box, was selected for time-series change detection using the proposed technique.</p>
Full article ">Figure 7
<p>Change detection results by utilizing KOMPSAT-5 time-series images: (<b>a</b>) pre-image, (<b>b</b>) post-image, (<b>c</b>) de-correlation between pre- and post-image, (<b>d</b>) adaptive thresholding results, (<b>e</b>) detected changed area, (<b>f</b>) ground truth data.</p>
Full article ">
22 pages, 7199 KiB  
Article
Three-Dimensional Deformation Prediction Based on the Improved Segmented Knothe–Dynamic Probabilistic Integral–Interferometric Synthetic Aperture Radar Model
by Shuang Wang, Genyuan Liu, Zhihong Song, Keming Yang, Ming Li, Yansi Chen and Minhua Wang
Remote Sens. 2025, 17(2), 261; https://doi.org/10.3390/rs17020261 - 13 Jan 2025
Viewed by 459
Abstract
Coal is the main mineral resource, but over-exploitation will cause a series of geological disasters. Interferometric synthetic aperture radar (InSAR) technology provides a superior monitoring method to compensate for the inadequacy of traditional measurements for mine surface deformation monitoring. In this study, the [...] Read more.
Coal is the main mineral resource, but over-exploitation will cause a series of geological disasters. Interferometric synthetic aperture radar (InSAR) technology provides a superior monitoring method to compensate for the inadequacy of traditional measurements for mine surface deformation monitoring. In this study, the whole process of mining a working face in Huaibei Mining District, Anhui Province, is taken as the object of study. The ALOS PALSAR satellite radar image data and ground measurements were acquired, and the ISK-DPIM-InSAR deformation monitoring model with the dynamic probabilistic integral model (DPIM) was proposed by combining the probabilistic integral method (PIM) and the improved segmented Knothe time function (ISK). The ISK-DPIM-InSAR model constructs the inversion equations of InSAR line-of-sight deformation, north–south and east–west horizontal movement deformation, vertical deformation, inverts the optimal values of the predicted parameters of the workforce through the particle swarm algorithm, and substitutes it into the ISK-DPIM-InSAR model for predicting the three-dimensional dynamic deformation of a mining face. Simulated workface experiments determined the feasibility of the model, and by comparing the level observation results of the working face, it is confirmed that the ISK-DPIM-InSAR model can accurately monitor the three-dimensional deformation of the surface in the mining area. Full article
Show Figures

Figure 1

Figure 1
<p>Dynamic projection model.</p>
Full article ">Figure 2
<p>InSAR side-view imaging geometric relationships and angular parameters.</p>
Full article ">Figure 3
<p>Particle swarm algorithm flow.</p>
Full article ">Figure 4
<p>Simulated working surface plan.</p>
Full article ">Figure 5
<p>Simulated work face of surface deformation between 300 and 350 days. (<b>a</b>) Simulated of work face of surface subsidence between 300 days and 350 days; (<b>b</b>) simulated work face of horizontal ground movement deformation along the east–west direction between 300 days and 350 days; (<b>c</b>) simulated work face of horizontal ground movement deformation along the north–south direction between 300 days and 350 days; and (<b>d</b>) simulated work face surface deformation values along the LOS between 300 days and 350 days.</p>
Full article ">Figure 6
<p>Simulation of LOS deformation results of extracted feature points.</p>
Full article ">Figure 7
<p>DPIM-InSAR model predicts 3D surface deformation values for simulated working surfaces. (<b>a</b>) Simulated working face of vertical surface deformation on day 350. (<b>b</b>) Simulated working face of horizontal ground movement deformation along the north–south direction on day 350. (<b>c</b>) Simulated working face of horizontal ground movement deformation along the east–west direction on day 350.</p>
Full article ">Figure 8
<p>Comparison of predicted and true values for simulated workface. (<b>a</b>) Comparison of subsidence along the AB observation line. (<b>b</b>) Comparison of subsidence along the CD observation line. (<b>c</b>) Comparison of horizontal movement deformation along east–west direction. (<b>d</b>) Comparison of horizontal movement deformation along the north–south direction.</p>
Full article ">Figure 9
<p>Location of the study area.</p>
Full article ">Figure 10
<p>LOS-oriented deformation results and characteristic points monitored by D-InSAR.</p>
Full article ">Figure 11
<p>Predicted vertical deformation results for DPIM-InSAR model inversion parameters.</p>
Full article ">Figure 12
<p>DPIM-InSAR model inversion parameters predicted north–south horizontal movement deformation results.</p>
Full article ">Figure 13
<p>DPIM-InSAR model inversion parameters predicted east–west horizontal movement deformation results.</p>
Full article ">Figure 14
<p>Comparison of DPIM-InSAR model projections and level data.</p>
Full article ">
23 pages, 10005 KiB  
Article
Time-Series InSAR Technology for Monitoring and Analyzing Surface Deformations in Mining Areas Affected by Fault Disturbances
by Kuan He, Youfeng Zou, Zhigang Han and Jilei Huang
Remote Sens. 2024, 16(24), 4811; https://doi.org/10.3390/rs16244811 - 23 Dec 2024
Viewed by 805
Abstract
Faults, as unique geological structures, disrupt the mechanical connections between rock masses. During coal mining, faults in the overlying strata can disturb the original stress balance, leading to fault activation and altering the typical subsidence patterns. This can result in abnormal ground deformation [...] Read more.
Faults, as unique geological structures, disrupt the mechanical connections between rock masses. During coal mining, faults in the overlying strata can disturb the original stress balance, leading to fault activation and altering the typical subsidence patterns. This can result in abnormal ground deformation and significant damage to surface structures, posing a serious geological hazard in mining areas. This study examines the influence of a known fault (F13 fault) on ground subsidence in the Wannian Mine of the Fengfeng Mining Area. We utilized 12 Sentinel-1A images and applied SBAS-InSAR, StaMPS-InSAR, and DS-InSAR time-series InSAR methods, alongside the D-InSAR method, to investigate surface deformations caused by the F13 fault. The monitoring accuracy of these methods was evaluated using leveling measurements from 28 surface movement observation stations. In addition, the density of effective monitoring points and the relative strengths and limitations of the three time-series methods were compared. The findings indicate that, in low deformation areas, DS-InSAR has a monitoring accuracy of 7.7 mm, StaMPS-InSAR has a monitoring accuracy of 16.4 mm, and SBAS-InSAR has an accuracy of 19.3 mm. Full article
Show Figures

Figure 1

Figure 1
<p>Technical roadmap overview.</p>
Full article ">Figure 2
<p>Overview map of the research area.</p>
Full article ">Figure 3
<p>Spatiotemporal baseline combination diagram. (The third image, acquired on 22 January 2020, is the super master image).</p>
Full article ">Figure 4
<p>SBAS−InSAR time−series deformation map.</p>
Full article ">Figure 5
<p>SBAS−InSAR deformation rate map.</p>
Full article ">Figure 6
<p>Deformation rate map of the study area was obtained using the StaMPS−InSAR method.</p>
Full article ">Figure 7
<p>The deformation rate map of the study area was obtained using the DS−InSAR method.</p>
Full article ">Figure 8
<p>Comparison of monitoring data from 28 surface movement observation stations.</p>
Full article ">Figure 9
<p>Mechanism diagram of surface deformation induced by fault disturbance. (<b>a</b>) Mechanical model of key layer before initial fracture without fault; (<b>b</b>) Mechanical model of key layers before initial fracture with faults.</p>
Full article ">Figure 10
<p>Time−series differential interferometry phase map.</p>
Full article ">Figure 11
<p>Three−dimensional deformation map generated by SBAS−InSAR.</p>
Full article ">Figure 12
<p>Comparison of average correlation coefficients before and after phase optimization ((<b>a</b>) average correlation coefficient from time−series D−InSAR; (<b>b</b>) average correlation coefficient after phase optimization; (<b>c</b>) comparison of correlation coefficients before and after phase optimization for the surface movement observation stations).</p>
Full article ">Figure 13
<p>Pearson correlation coefficient analysis of monitoring results in the deformation zone. ((<b>a</b>) Pearson correlation coefficient between DS−InSAR and StaMPS−InSAR; (<b>b</b>) Pearson correlation coefficient between SBAS−InSAR and DS−InSAR; (<b>c</b>) Pearson correlation coefficient between SBAS−InSAR and StaMPS−InSAR).</p>
Full article ">
18 pages, 31612 KiB  
Article
Land Subsidence Velocity and High-Speed Railway Risks in the Coastal Cities of Beijing–Tianjin–Hebei, China, with 2015–2021 ALOS PALSAR-2 Multi-Temporal InSAR Analysis
by Qingli Luo, Mengli Li, Zhiyuan Yin, Peifeng Ma, Daniele Perissin and Yuanzhi Zhang
Remote Sens. 2024, 16(24), 4774; https://doi.org/10.3390/rs16244774 - 21 Dec 2024
Viewed by 576
Abstract
Sea-level rise has important implications for the economic and infrastructure security of coastal cities. Land subsidence further exacerbates relative sea-level rise. The Beijing–Tianjin–Hebei region (BTHR) along the Bohai Bay is one of the areas most severely affected by ground subsidence in the world. [...] Read more.
Sea-level rise has important implications for the economic and infrastructure security of coastal cities. Land subsidence further exacerbates relative sea-level rise. The Beijing–Tianjin–Hebei region (BTHR) along the Bohai Bay is one of the areas most severely affected by ground subsidence in the world. This study applies the Small Baseline Subset Interferometric Synthetic Aperture Radar (SBAS InSAR) method to analyze 47 ALOS PALSAR-2 images with five frames, mapping subsidence across 21,677.7 km2 and revealing spatial patterns and trends over time from 2015 to 2021. This is one of the few published research studies for large-scale and long-term analysis of its kind using ALOS-2 data in this region. The results reveal the existence of six major areas affected by severe subsidence in the study area, with the most pronounced in Jinzhan Town, Beijing, with the maximum subsiding velocity of −94.42 mm/y. Except for the two subsidence areas located in Chaoyang District of Beijing and Guangyang District of Langfang City, the other areas with serious subsidence detected are all located in suburban areas; this means that the strict regulations of controlling urban subsidence for downtown areas in the BTHR have worked. The accumulated subsidence is highly correlated with the time in the time series. Moreover, the subsidence of 161.4 km of the Beijing–Tianjin Inter-City High-Speed Railway (HSR) and 194.5 km of the Beijing–Shanghai HSR (out of a total length of 1318 km) were analyzed. It is the first time that PALSAR-2 data have been used to simultaneously investigate the subsidence along two important HSR lines in China and to analyze relatively long sections of the routes. The above two railways intersect five and seven subsiding areas, respectively. Within the range of the monitored railway line, the percentage of the section with subsidence velocity below −10 mm/y in the monitoring length range is 11.2% and 27.9%; this indicates that the Beijing–Shanghai HSR has suffered more serious subsidence than the Beijing–Tianjin Inter-City HSR within the monitoring period. This research is also beneficial for assessing the subsidence risk associated with different railways. In addition, this study further analyzed the potential reasons for the serious land subsidence of the identified areas. The results of the geological interpretation still indicate that the main cause of subsidence in the area is due to hydrogeological characteristics and underground water withdrawal. Full article
(This article belongs to the Special Issue Synthetic Aperture Radar Interferometry Symposium 2024)
Show Figures

Figure 1

Figure 1
<p>Study area and coverage of the datasets: (<b>a</b>) the location of the study area in China; (<b>b</b>) the spatial coverage of the ALOS PALSAR-2 data in the study area; (<b>c</b>) the coverage of ALOS PALSAR-2 data overlapped with Tianditu maps; the blue rectangles in (<b>b</b>,<b>c</b>) are the spatial coverage of the ALOS PALSAR-2 data; and the orange and red lines in (<b>c</b>) highlight the Beijing–Tianjin Inter-City HSR and the Beijing–Shanghai HSR, respectively.</p>
Full article ">Figure 2
<p>Differential interference results of Frame 770A between 9 May 2019 and 6 May 2021: (<b>a</b>) topographic information; (<b>b</b>) the ALOS PALSAR-2 image from 9 May 2019; (<b>c</b>) the differential interferogram; and (<b>d</b>) the coherence.</p>
Full article ">Figure 3
<p>Processing for PALSAR-2 imagery by SBAS.</p>
Full article ">Figure 4
<p>The overall average subsidence velocity in the BTHR. Stable areas are indicated in green, while orange and red represent subsiding areas with the velocity of −30 and -90 mm/y, respectively. The main subsidence centers are identified as follows: (<b>A</b>) Jinzhan Town, (<b>B</b>) Wenjing Sub-District, (<b>C</b>) Guangyang District, (<b>D</b>) Wangqingtuo Town, (<b>E</b>) Shengfang Town, and (<b>F</b>) Tuanbo Town.</p>
Full article ">Figure 5
<p>Subsidence history of ground points in subsidence funnel.</p>
Full article ">Figure 6
<p>The partially enlarged average deformation map of Beijing–Tianjin Inter-City HSR. A1–A5 are the main five serious subsiding centers along the Beijing–Tianjin Inter-City HSR.</p>
Full article ">Figure 7
<p>The partially enlarged average deformation map of Beijing–Shanghai HSR. B1–B7 are the main seven subsiding centers along the Beijing–Shanghai HSR.</p>
Full article ">Figure 8
<p>Comparisons of SBAS InSAR results and leveling data: (<b>a</b>) the Beijing–Tianjin Inter-City HSR; (<b>b</b>) the Beijing–Shanghai HSR; (<b>c</b>) scatterplot for correlation analysis along the Beijing–Tianjin Inter-City HSR; and (<b>d</b>) scatterplot along the Beijing–Shanghai HSR.</p>
Full article ">Figure 9
<p>Land resources distribution and groundwater extraction in the BTHR: (<b>a</b>) groundwater exploitation situation; and (<b>b</b>) geological map. These two maps are from the 2015 China Geological Survey conducted by the Ministry of Land and Resources. A, B, C, D, E, and F are the main six subsiding centers marked in <a href="#remotesensing-16-04774-f004" class="html-fig">Figure 4</a>.</p>
Full article ">Figure 10
<p>Profiles of average subsidence velocity along railways from SBAS InSAR results: (<b>a</b>) the profile of Beijing–Tianjin Inter-City HSR; (<b>b</b>) the profile of Beijing–Shanghai HSR. A1–A5 and B1–B7 are the sections affected by relatively serious subsidence along the Beijing–Tianjin Inter-City HSR and the Beijing–Shanghai HSR, respectively.</p>
Full article ">Figure 11
<p>The histograms of the distances between the point pairs from the HSRs and the SBAS InSAR results: (<b>a</b>) the histogram of the distances along Beijing–Tianjin Inter-City HSR; and (<b>b</b>) the histogram of the distances along Beijing–Shanghai HSR.</p>
Full article ">
19 pages, 32077 KiB  
Article
Present-Day Tectonic Deformation Characteristics of the Northeastern Pamir Margin Constrained by InSAR and GPS Observations
by Junjie Zhang, Xiaogang Song, Donglin Wu and Xinjian Shan
Remote Sens. 2024, 16(24), 4771; https://doi.org/10.3390/rs16244771 - 21 Dec 2024
Viewed by 599
Abstract
The Pamir is located on the northwestern margin of the Tibetan Plateau, which is an area of intense continental deformation and part of the famous India–Himalaya collision zone. The dominant structural deformation in the eastern Pamir is characterized by a 250 km long [...] Read more.
The Pamir is located on the northwestern margin of the Tibetan Plateau, which is an area of intense continental deformation and part of the famous India–Himalaya collision zone. The dominant structural deformation in the eastern Pamir is characterized by a 250 km long east–west extensional fault system, known as the Kongur Shan extensional system (KSES), which has developed a series of faults with different orientations and characteristics, resulting in highly complex structural deformation and lacking sufficient geodetic constraints. We collected Sentinel-1 SAR data from December 2016 to March 2023, obtained high-resolution ascending and descending LOS velocities and 3D deformation fields, and combined them with GPS data to constrain the current motion characteristics of the northeastern Pamirs for the first time. Based on the two-dimensional screw dislocation model and using the Bayesian Markov chain Monte Carlo (MCMC) inversion method, the kinematic parameters of the fault were calculated, revealing the fault kinematic characteristics in this region. Our results demonstrate that the present-day deformation of the KSES is dominated by nearly E–W extension, with maximum extensional motion concentrated in its central segment, reaching peak extension rates of ~7.59 mm/yr corresponding to the Kongur Shan. The right-lateral Muji fault at the northern end exhibits equivalent rates of extensional motion with a relatively shallow locking depth. The strike-slip rate along the Muji fault gradually increases from west to east, ranging approximately between 4 and 6 mm/yr, significantly influenced by the eastern normal fault. The Tahman fault (TKF) at the southernmost end of the KSES shows an extension rate of ~1.5 mm/yr accompanied by minor strike-slip motion. The Kashi anticline is approaching stability, while the Mushi anticline along the eastern Pamir frontal thrust (PFT) remains active with continuous uplift at ~2 mm/yr, indicating that deformation along the Tarim Basin–Tian Shan boundary has propagated southward from the South Tian Shan thrust (STST). Overall, this study demonstrates the effectiveness of integrated InSAR and GPS data in constraining contemporary deformation patterns along the northeastern Pamir margin, contributing to our understanding of the region’s tectonic characteristics. Full article
Show Figures

Figure 1

Figure 1
<p>Tectonics and seismicity of the study area. (<b>a</b>) The yellow and red rectangle shows the spatial footprint of the Sentinel-1 InSAR coverage. Blue arrows show the GPS horizontal velocity field with respect to the stable Eurasian plate [<a href="#B28-remotesensing-16-04771" class="html-bibr">28</a>]. Circles of different colors represent earthquake events of varying magnitudes. (<b>b</b>) Fault structures in the eastern part of PFT. Red and pink focal mechanisms represent the mainshock and aftershocks of the 1985 Wuqia earthquake. Brown focal mechanisms represent the mainshock of the 2016 Aketao earthquake. (<b>c</b>) Fault segmentation in KESE. S1–S5 correspond to different segments, respectively. STST = southern Tian Shan thrust, PFT = Pamir frontal thrust, KSES = Kongur Shan extensional system, MPT = main Pamir thrust, KKF = Karakax fault, KYTS = Kashgar–Yecheng transfer system, TFF = Talas–Fergana fault, TT = Tuomuluoan thrust, MF = Muji fault, KATF = King Ata Tagh normal fault, KSF = Kongur Shan normal fault, MAF = Muztagh Ata normal fault, TF = Tahman normal fault, TKF = Tashkorgan normal fault.</p>
Full article ">Figure 2
<p>The Sentinel-1 A/B data processing workflow. It consists of three steps, including interferograms generation, SBAS time series analysis, and three-dimensional deformation field solution.</p>
Full article ">Figure 3
<p>Perpendicular and temporal baseline plot showing the network of interferograms on one ascending track (<b>a</b>) and one descending track (<b>b</b>) used in this study. The number of total interferograms are labelled for each track.</p>
Full article ">Figure 4
<p>Interpolated GNSS velocities map. The interpolated GNSS velocities using the method outlined by Shen et al. [<a href="#B47-remotesensing-16-04771" class="html-bibr">47</a>], (<b>a</b>) corresponds to EW and (<b>b</b>) corresponds to NS. GNSS velocities are resampled to a resolution of 0.01 degrees. Different colored circles represent different GPS data, and the color bars for GPS various data and the interpolated velocity field are identical.</p>
Full article ">Figure 5
<p>The satellite line-of-sight (LOS) velocity fields of the northeastern Pamir margin. Red lines represent the fault crossing profiles, each profile for 60 km long and 10 km wide, distributed along six sub-faults of the KSES. (<b>a</b>) corresponds to ascending track and (<b>b</b>) corresponds to descending track.</p>
Full article ">Figure 6
<p>Joint InSAR-GPS three-dimensional deformation field. (<b>a</b>–<b>c</b>) are east–west, north–south, and vertical components, respectively.</p>
Full article ">Figure 7
<p>GPS profiles and results. (<b>a</b>) Four profiles in KSES. Each profile is 300 km long and 50 km wide. (<b>b</b>) The GPS data was projected parallel to and perpendicular to the local fault, respectively, and combined with fault strike.</p>
Full article ">Figure 8
<p>The cross-fault profiles of ascending and descending LOS deformation velocities. (<b>a</b>) represents ascending LOS velocity profiles, and (<b>b</b>) represents descending LOS velocity profiles. (<b>a</b>–<b>f</b>) correspond successively to the six profiles in <a href="#remotesensing-16-04771-f005" class="html-fig">Figure 5</a>. Black dots are binned average values every 1 km along the profile. Gray vertical stripes indicate the mountains on profiles. Red and purple lines are the best-fitting models. The black dotted line indicates the fault location.</p>
Full article ">Figure 9
<p>An example of Bayesian MCMC inversion results for profile aa’. Posterior marginal probability density functions illustrating parameter estimation and uncertainty quantification. (<b>Top</b>): profile aa’ topography from the Copernicus DEM data with 30 m spatial resolution (average elevation: white line; min/max: gray lines). (<b>Middle</b>): InSAR LOS velocities with the best-fitting predicted velocities. (<b>Bottom</b>): model-predicted fault-parallel, fault-normal, and vertical velocities.</p>
Full article ">Figure 10
<p>The LOS velocities and vertical component of the eastern PFT. (<b>a</b>–<b>c</b>) are ascending, descending, and vertical components, respectively. The abnormal deformation area corresponding to the black circle and black rectangle are caused by industrial activity. The satellite images corresponding to these two regions are shown in <a href="#app1-remotesensing-16-04771" class="html-app">Figure S5</a>. (<b>d</b>–<b>f</b>) correspond to the results of profiles aa’, bb’, and cc’ in (<b>a</b>–<b>c</b>), respectively. The pink, blue, and green points correspond to the results of profiles aa’, bb’, and cc’, respectively.</p>
Full article ">
21 pages, 9480 KiB  
Article
Collapse Hotspot Detection in Urban Area Using Sentinel-1 and TerraSAR-X Dataset with SBAS and PSI Techniques
by Niloofar Alizadeh, Yasser Maghsoudi, Tayebe Managhebi and Saeed Azadnejad
Land 2024, 13(12), 2237; https://doi.org/10.3390/land13122237 - 20 Dec 2024
Viewed by 720
Abstract
Urban areas face an imminent risk of collapse due to structural deficiencies and gradual ground subsidence. Therefore, monitoring surface movements is crucial for detecting abnormal behavior, implementing timely preventive measures, and minimizing the detrimental effects of this phenomenon in residential regions. In this [...] Read more.
Urban areas face an imminent risk of collapse due to structural deficiencies and gradual ground subsidence. Therefore, monitoring surface movements is crucial for detecting abnormal behavior, implementing timely preventive measures, and minimizing the detrimental effects of this phenomenon in residential regions. In this context, interferometric synthetic aperture radar (InSAR) has emerged as a highly effective technique for monitoring slow and long-term ground hazards and surface motions. The first goal of this study is to explore the potential applications of persistent scatterer interferometry (PSI) and small baseline subset (SBAS) algorithms in collapse hotspot detection, utilizing a dataset consisting of 144 Sentinel-1 images. The experimental results from three areas with a history of collapses demonstrate that the SBAS algorithm outperforms PSI in uncovering behavior patterns indicative of collapse and accurately pinpointing collapse points near real collapse sites. In the second phase, this research incorporated an additional dataset of 36 TerraSAR-X images alongside the Sentinel-1 data to compare results based on radar images with different spatial resolutions in the C and X bands. The findings reveal a strong correlation between the TerraSAR-X and Sentinel-1 time series. Notably, the analysis of the TerraSAR-X time series for one study area identified additional collapse-prone points near the accident site, attributed to the higher spatial resolution of these data. By leveraging the capabilities of InSAR and advanced algorithms, like SBAS, this study highlights the potential to identify areas at risk of collapse, enabling the implementation of preventive measures and reducing potential harm to residential communities. Full article
(This article belongs to the Special Issue Assessing Land Subsidence Using Remote Sensing Data)
Show Figures

Figure 1

Figure 1
<p>Study area location; (<b>a</b>) Map of Tehran with InSAR-derived subsidence data overlay; (<b>b</b>) Abshenasan Blvd.; (<b>c</b>) Moghadam Sq.; (<b>d</b>) Shariati St.; (<b>e</b>) Karimkhan St.</p>
Full article ">Figure 2
<p>Workflow of the proposed algorithm, comprising three main stages: preprocessing, InSAR time series analysis using PSI and SBAS methods, and outlier and collapse hotspot detection.</p>
Full article ">Figure 3
<p>The ground deformation rate map; (<b>a</b>) PSI velocity map in Abshenasan Blvd.; (<b>b</b>) SBAS velocity map in Abshenasan Blvd.; (<b>c</b>) PSI velocity map in Shariati St. (<b>d</b>) SBAS velocity map in Shariati St.; (<b>e</b>) PSI velocity map in Karimkhan St.; (<b>f</b>) SBAS velocity map in Karimkhan St.</p>
Full article ">Figure 4
<p>Comparison between point displacement rates estimated based on the PSI and SBAS methods, the y-axis represents the PS displacement rate in the PSI algorithm and the x-axis represents the estimated displacement rate in the corresponding pixel based on the SBAS algorithm; (<b>a</b>) scatter plot in Abshenasan Blvd.; (<b>b</b>) scatter plot in Karimkhan St.; (<b>c</b>) scatter plot in Shariati St.</p>
Full article ">Figure 5
<p>The process of outlier detection in a buffer for the first and second interval by Z-score and IQR algorithms; (<b>a</b>,<b>c</b>,<b>e</b>) detected as outlier; (<b>b</b>,<b>d</b>,<b>f</b>) selected for further steps.</p>
Full article ">Figure 6
<p>The performance of PSI and SBAS algorithms in identifying the point prone to collapse (<b>a</b>) time series of LOS displacement in PSI algorithm in Abshenasan Blvd.; (<b>b</b>) time series of LOS displacement in SBAS algorithm in Abshenasan Blvd.; (<b>c</b>) time series of LOS displacement in PSI algorithm in Shariati St.; (<b>d</b>) time series of LOS displacement in SBAS algorithm in Shariati St.; (<b>e</b>) time series of LOS displacement in PSI algorithm in Karimkhan St.; (<b>f</b>) time series of LOS displacement in SBAS algorithm in Karimkhan St.</p>
Full article ">Figure 7
<p>PSs and points detected in a 50 m buffer. The red point indicates the identified collapse-prone point and the red circle is the buffer area (<b>a</b>) in Abshenasan Blvd. by the PSI algorithm; (<b>b</b>) in Abshenasan Blvd. by the SBAS algorithm; (<b>c</b>) in Shariati St. by the PSI algorithm; (<b>d</b>) in Shariati St. by the SBAS algorithm; (<b>e</b>) in Karimkhan St. by the PSI algorithm; (<b>f</b>) in Karimkhan St. by the SBAS algorithm.</p>
Full article ">Figure 8
<p>The comparison between the ground surface deformation rate obtained from SBAS InSAR; (<b>a</b>) S1 velocity map in Moghadam Sq; (<b>b</b>) TSX velocity map in Moghadam Sq.</p>
Full article ">Figure 9
<p>Comparison between velocities extracted from both S1 and TSX dataset in Moghadam Sq.</p>
Full article ">Figure 10
<p>The performance of SBAS algorithms in identifying the point prone to collapse in Moghadam Sq; (<b>a</b>) time series of LOS displacement in S1 dataset; (<b>b</b>–<b>d</b>) time series of LOS displacement in TSX dataset.</p>
Full article ">Figure 11
<p>Points detected in a 50 m buffer in Moghadam Sq. by SBAS algorithm. The red point indicates the identified collapse-prone point and the red circle is the buffer area; (<b>a</b>) S1 dataset; (<b>b</b>) TSX dataset.</p>
Full article ">
8 pages, 3869 KiB  
Proceeding Paper
Urban Road Collapse Risk Assessment Study Based on InSAR Spatiotemporal Data
by Juncai Jiang, Wenfeng Bai, Yizhao Wang, Fei Wang, Qinglun He, Long Chen, Yuming Qiao, Zhi Wang and Haitao Luo
Proceedings 2024, 110(1), 24; https://doi.org/10.3390/proceedings2024110024 - 12 Dec 2024
Viewed by 640
Abstract
Urban road collapse is a common disaster in modern cities, posing a severe threat to the safety of life and property of urban residents. Effective risk assessment is crucial for preventing collapses. This study proposes a novel method for assessing the risk of [...] Read more.
Urban road collapse is a common disaster in modern cities, posing a severe threat to the safety of life and property of urban residents. Effective risk assessment is crucial for preventing collapses. This study proposes a novel method for assessing the risk of urban road collapse by integrating Interferometric Synthetic Aperture Radar (InSAR) technology with Long Short-Term Memory (LSTM) networks. We conduct experimental analysis using collapse incidents and corresponding subsidence maps from five areas in Guangzhou. The results demonstrate that the proposed model successfully establishes a mapping relationship between the temporal features of InSAR data and the risk of road collapse, enabling a risk assessment of specific areas solely based on InSAR data. This model overcomes the difficulty of obtaining numerous assessment indicators in traditional risk assessment methods, providing strong support for the prevention and control of regional road collapse risks. Full article
(This article belongs to the Proceedings of The 31st International Conference on Geoinformatics)
Show Figures

Figure 1

Figure 1
<p>The location of Guangzhou city.</p>
Full article ">Figure 2
<p>Distribution of collapse accidents in 5 urban areas of Guangzhou.</p>
Full article ">Figure 3
<p>InSAR imagery in the study area.</p>
Full article ">Figure 4
<p>LSTM structure.</p>
Full article ">Figure 5
<p>Loss function training.</p>
Full article ">Figure 6
<p>Risk map.</p>
Full article ">
23 pages, 28195 KiB  
Article
Slow-Moving Landslide Hazard Assessment Using LS-Unilab Deep Learning Model with Highlighted InSAR Deformation Signal
by Xiangyang Li, Peifeng Ma, Song Xu, Hong Zhang, Chao Wang, Yukun Fan and Yixian Tang
Remote Sens. 2024, 16(24), 4641; https://doi.org/10.3390/rs16244641 - 11 Dec 2024
Viewed by 789
Abstract
Slow-moving landslides are often precursors of catastrophic failure, posing a major threat to human life and property safety. Interferometric synthetic aperture radar (InSAR) has become a crucial tool for investigating slow-moving landslides hazard because of its high-precision detection capability for slow surface deformation. [...] Read more.
Slow-moving landslides are often precursors of catastrophic failure, posing a major threat to human life and property safety. Interferometric synthetic aperture radar (InSAR) has become a crucial tool for investigating slow-moving landslides hazard because of its high-precision detection capability for slow surface deformation. However, landslides usually occur in alpine canyon areas and vegetation coverage areas where InSAR measurements are still limited by temporal and spatial decorrelation and atmospheric influences. In addition, there are several difficulties in monitoring the multiscale characterization of landslides from the InSAR results. To address this issue, this paper proposes a novel method for slow-moving landslide hazard assessment in low-coherence regions. A window-based atmosphere correction method is designed to highlight the surface deformation signals of InSAR results in low-coherence regions and reduce false alarms in landslide hazard assessment. Then, the deformation annual velocity rate map, coherence map and DEM are used to construct the InSAR sample set. A landslide hazard assessment model named Landslide-SE-Unilab is subsequently proposed. The global–local relationship aggregation structure is designed to capture the spatial relationship between local pixel-level deformation features and global landslides, which can reduce the number of missed assessments and false assessments of small-scale landslides. Additionally, a squeeze-and-excitation network is embedded to adjust the weight relationship between the features of each channel in order to enhance the performance of network evaluation. The method was evaluated in Kangding city and the Jinsha River Valley in the Hengduan Mountains, where a total of 778 potential landslides with slow deformation were identified. The effectiveness and accuracy of this approach for low-coherence landslide hazard assessment are demonstrated through comparisons with optical images and previous research findings, as well as evaluations via time-series deformation results. Full article
Show Figures

Figure 1

Figure 1
<p>Flowchart of the proposed technique.</p>
Full article ">Figure 2
<p>Flowchart of SBAS-InSAR with a window-based atmospheric correction.</p>
Full article ">Figure 3
<p>Flowchart of the sample production process.</p>
Full article ">Figure 4
<p>The LS-Unilab model. The deformation annual velocity rate map, coherence map, and DEM are selected for the model input.</p>
Full article ">Figure 5
<p>Study area and fault distribution. The black lines represent faults (source: <a href="https://docs.gmt-china.org/latest/dataset-CN/CN-faults/" target="_blank">https://docs.gmt-china.org/latest/dataset-CN/CN-faults/</a>, accessed on 16 May 2024). The red dots denote the earthquake locations since 2008 (source: <a href="https://data.earthquake.cn/" target="_blank">https://data.earthquake.cn/</a>, accessed on 16 May 2024), and the black boxes represent the Sentinel-1 data coverage used in this work. The background is the SRTM1 DEM (source: <a href="http://step.esa.int/auxdata/dem/SRTMGL1/" target="_blank">http://step.esa.int/auxdata/dem/SRTMGL1/</a>, accessed on 16 May 2024).</p>
Full article ">Figure 6
<p>Field photographs of the landslides along the Jinsha River.</p>
Full article ">Figure 7
<p>Annual velocity rate of path 26 from Sentinel-1 images from Jan 2022 to Sep 2023 and statistical results, where regions A–D are selected for detailed analysis. (<b>a</b>) The uncorrected results; (<b>b</b>) the elevation correction results; (<b>c</b>) the window based atmospheric correction results; and (<b>d</b>) the statistical results of (<b>a</b>,<b>c</b>).</p>
Full article ">Figure 8
<p>Annual deformation velocity of Kangding city (<b>a</b>) and the Jinsha River Gorge (<b>b</b>) from Sentinel-1 images from January 2022 to September 2023, where regions I–VI are selected for detailed analysis. (<b>c</b>) Zoomed-in view of area IV in (<b>b</b>), where the locations of P1–P6 correspond to the field photographs in <a href="#remotesensing-16-04641-f006" class="html-fig">Figure 6</a>. (<b>d</b>,<b>e</b>) (corresponding to areas (4) and (3) in <a href="#remotesensing-16-04641-f009" class="html-fig">Figure 9</a>) Corresponded to regions A and B in black circle of (<b>a</b>); (<b>f</b>,<b>g</b>) (corresponded to areas (2) and (1) in <a href="#remotesensing-16-04641-f009" class="html-fig">Figure 9</a>) Corresponded to regions V and VI in (<b>b</b>).</p>
Full article ">Figure 9
<p>Assessment results of slow-moving landslides, where regions in circles are selected for detailed analysis.</p>
Full article ">Figure 10
<p>Landslide hazard assessment results and statistical results for Kangding city (<b>a</b>,<b>c</b>) and the Jinsha River Gorge (<b>b</b>,<b>d</b>). The red triangles represent the locations of slow-moving landslides.</p>
Full article ">Figure 11
<p>Validation region in Kangding City; (<b>a</b>) Annual deformation rate map; (<b>b</b>) base image of the Sentinel-2 optical image; (<b>c</b>) model identification results; (<b>d</b>) threshold separation results.</p>
Full article ">Figure 12
<p>Validation region in the Jinsha River Gorge; (<b>a</b>) annual deformation rate map; (<b>b</b>) base image of the Sentinel-2 optical image; (<b>c</b>) model identification results; (<b>d</b>) threshold separation results.</p>
Full article ">Figure 13
<p>The upper background image is a Google Earth image overlaid with deformation rates, with red rectangles indicating the landslide identification results; the lower part shows the time-series deformation results of the monitoring points. (<b>a</b>–<b>d</b>) Areas of verification points in <a href="#remotesensing-16-04641-f010" class="html-fig">Figure 10</a>.</p>
Full article ">Figure 14
<p>Compared with other research results, the deformation period time is marked. (<b>a</b>) Results obtained by Zou et al. [<a href="#B45-remotesensing-16-04641" class="html-bibr">45</a>]; (<b>c</b>,<b>d</b>) results obtained by Liu et al. [<a href="#B12-remotesensing-16-04641" class="html-bibr">12</a>,<a href="#B13-remotesensing-16-04641" class="html-bibr">13</a>]; (<b>e</b>) results obtained by Zhang et al. [<a href="#B50-remotesensing-16-04641" class="html-bibr">50</a>]; (<b>b</b>,<b>f</b>) results obtained in the present study; (<b>g</b>–<b>k</b>) the Sentinel-2 optical imagery of areas delineated by black rectangles in (<b>f</b>). Red circels are selected for detailed analysis. The legend of the original text in the figure was redrawn.</p>
Full article ">
25 pages, 41258 KiB  
Article
The Deformation Monitoring Capability of Fucheng-1 Time-Series InSAR
by Zhouhang Wu, Wenjun Zhang, Jialun Cai, Hongyao Xiang, Jing Fan and Xiaomeng Wang
Sensors 2024, 24(23), 7604; https://doi.org/10.3390/s24237604 - 28 Nov 2024
Viewed by 801
Abstract
The Fucheng-1 (FC-1) satellite has successfully transitioned from its initial operational phase and is now undergoing a detailed performance assessment for time-series deformation monitoring. This study evaluates the surface deformation monitoring capabilities of the newly launched FC-1 satellite using the interferometric synthetic aperture [...] Read more.
The Fucheng-1 (FC-1) satellite has successfully transitioned from its initial operational phase and is now undergoing a detailed performance assessment for time-series deformation monitoring. This study evaluates the surface deformation monitoring capabilities of the newly launched FC-1 satellite using the interferometric synthetic aperture radar (InSAR) technique, particularly in urban applications. By analyzing the observation data from 20 FC-1 scenes and 20 Sentinel-1 scenes, deformation velocity maps of a university in Mianyang city were obtained using persistent scatterer interferometry (PSI) and distributed scatterer interferometry (DSI) techniques. The results show that thanks to the high resolution of 3 × 3 m of the FC-1 satellite, significantly more PS points and DS points were detected than those detected by Sentinel-1, by 13.4 times and 17.9 times, respectively. The distribution of the major deformation areas detected by both satellites in the velocity maps is generally consistent. FC-1 performs better than Sentinel-1 in monitoring densely structured and vegetation-covered areas. Its deformation monitoring capability at the millimeter level was further validated through comparison with leveling measurements, with average errors and root mean square errors of 1.761 mm and 2.172 mm, respectively. Its high-resolution and high-precision interferometry capabilities make it particularly promising in the commercial remote sensing market. Full article
(This article belongs to the Special Issue Recent Advances in Synthetic Aperture Radar (SAR) Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Coverage areas of Sentinel-1 (purple) and FC-1 (brown), study area location marked by a five-pointed star, and COPDEM topographic map. (<b>b</b>) Google Maps image of the study area.</p>
Full article ">Figure 2
<p>Flow chart of DSI and PSI.</p>
Full article ">Figure 3
<p>(<b>a</b>) Spatio-temporal baseline map of FC-1 single master image. (<b>b</b>) Spatio-temporal baseline map of Sentinel-1 single master image.</p>
Full article ">Figure 4
<p>(<b>a</b>,<b>b</b>) Vertical deformation velocity maps from FC-1 using the DSI and PSI methods. (<b>c</b>,<b>d</b>) Vertical deformation velocity maps from Sentinel-1 using the DSI and PSI methods. (<b>e</b>) Drone orthophoto of the reference point.</p>
Full article ">Figure 5
<p>(<b>a</b>,<b>b</b>) Histograms of deformation velocity from FC-1 using the DSI and PSI methods. (<b>c</b>,<b>d</b>) Histograms of deformation velocity from Sentinel-1 using the DSI and PSI methods.</p>
Full article ">Figure 6
<p>(<b>a</b>) Schematic diagram of the research area on Google Earth. (<b>b</b>,<b>c</b>) Deformation rate maps of region R1 obtained by FC-1 and Sentinel-1 using the PSI method, with a drone image as the base map. (<b>d</b>–<b>g</b>) Deformation rate maps of regions R2 and R3 obtained by FC-1 and Sentinel-1 using the PSI method, with Google Earth as the base map. (<b>h</b>–<b>k</b>) Deformation rate maps of regions R4 and R5 obtained by FC-1 and Sentinel-1 using the DSI method, with Google Earth or a drone image as the base map.</p>
Full article ">Figure 7
<p>(<b>a</b>) Deformation velocity points obtained by FC-1 using the PSI method overlaid onto a drone image. (<b>b</b>) Deformation velocity points obtained by Sentinel-1 using the PSI method overlaid onto a drone image.</p>
Full article ">Figure 8
<p>(<b>a</b>,<b>d</b>) Deformation velocity maps from FC-1 and Sentinel-1 using the PSI method, with schematic maps of ZZ1 and ZZ2 locations. (<b>b</b>,<b>c</b>) PS deformation points from FC-1 overlaid onto drone oblique images of ZZ1 and ZZ2. (<b>e</b>,<b>f</b>) PS deformation points from Sentinel-1 overlaid onto drone oblique images of ZZ1 and ZZ2.</p>
Full article ">Figure 9
<p>(<b>a</b>,<b>b</b>) Deformation velocity points obtained by FC-1 using DSI and PSI methods overlaid onto Google imagery. (<b>c</b>,<b>d</b>) Deformation velocity points obtained by Sentinel-1 using DSI and PSI methods overlaid onto Google imagery.</p>
Full article ">Figure 10
<p>Diagram of road profile location.</p>
Full article ">Figure 11
<p>(<b>a</b>,<b>b</b>) Deformation velocity profile of FC-1 under the DSI and PSI methods. (<b>c</b>,<b>d</b>) Deformation velocity profile of Sentinel-1 under the DSI and PSI methods.</p>
Full article ">Figure 12
<p>(<b>a</b>,<b>b</b>) Deformation rate profiles of Sentinel-1 and FC-1 under the DSI method. (<b>c</b>) Diagram of position of vegetation section line. (<b>d</b>) UAV 3D model of vegetation area.</p>
Full article ">Figure 13
<p>(<b>a</b>,<b>b</b>) Coherence histograms and average coherence values for the PSI method with FC-1 and Sentinel-1. (<b>c</b>,<b>d</b>) Coherence histograms and average coherence values for the DSI method with FC-1 and Sentinel-1.</p>
Full article ">Figure 14
<p>(<b>a</b>,<b>b</b>) Standard deviation maps of deformation velocity for Sentinel-1 using PSI and DSI methods. (<b>c</b>,<b>d</b>) Standard deviation maps of deformation velocity for FC-1 using PSI and DSI methods.</p>
Full article ">Figure 15
<p>(<b>a</b>) Diagram of locations of four regions A, B, C and D. (<b>b</b>–<b>e</b>) Time-series settlement maps of FC-1 and Sentinel-1 using the DSI and PSI methods in four regions.</p>
Full article ">Figure 16
<p>(<b>a</b>–<b>d</b>) Spearman’s correlation matrix heatmaps of the time-series settlement amounts obtained by FC-1 and Sentinel-1 using the DSI and PSI methods in four regions.</p>
Full article ">Figure 17
<p>(<b>a</b>–<b>d</b>) Pearson’s correlation matrix plots of the time-series subsidence values between FC-1 and Sentinel-1 using the DSI and PSI methods in four regions.</p>
Full article ">Figure 18
<p>Illustrative Google Earth map showing the locations of level points.</p>
Full article ">Figure 19
<p>The subsidence measured by FC-1 using the DSI method compared to the subsidence measured by leveling.</p>
Full article ">Figure 20
<p>The subsidence measured by FC-1 and Sentinel-1 using the DSI method compared to the subsidence measured by leveling.</p>
Full article ">
17 pages, 7864 KiB  
Article
Three-Dimensional Monitoring of Zelongnong Glacier, China, with the PO-MSBAS Technique
by Xinyi Zhai, Chaoying Zhao, Bin Li, Wenpei Wang and Xiaojie Liu
Remote Sens. 2024, 16(23), 4462; https://doi.org/10.3390/rs16234462 - 28 Nov 2024
Viewed by 536
Abstract
High-precision monitoring of glacier motion provides crucial information for a thorough understanding of the dynamic characteristics and development patterns of glaciers, which serves as a scientific basis for the prevention and management of glacier-related disasters. Zelongnong Glacier, located in Tibet, China, has experienced [...] Read more.
High-precision monitoring of glacier motion provides crucial information for a thorough understanding of the dynamic characteristics and development patterns of glaciers, which serves as a scientific basis for the prevention and management of glacier-related disasters. Zelongnong Glacier, located in Tibet, China, has experienced glacier surges, collapse, and hazard chains four times in the last 70 years. On 10 September 2020, a major glacier hazard chain occurred in this region. To reveal the influencing factors of the glacier motion, we monitor the Zelongnong Glacier motions with 65 scenes of TerraSAR/PAZ images from 2022 to 2023, where the Pixel Offset Multidimensional Small Baseline Subset (PO-MSBAS) method is employed for three-dimensional time series inversion. As the registration window size directly affects the matching success rate, deformation accuracy, and signal-to-noise ratio (SNR) during the offset tracking processing, we adopt a variable window-weighted cross-correlation strategy. The strategy balances the advantages of different window sizes, effectively reducing noise while preserving certain details in the offset results. The standard deviation in stable areas is also significantly lower than that obtained using smaller window sizes in conventional methods. The results reveal that the velocity of the southern glacier tributary was larger than the one in the northern tributary. Specifically, the maximum velocity in the northern tributary reached 45.07 m/year in the horizontal direction and −7.45 m/year in the vertical direction, whereas in the southern tributary, the maximum velocity was 50.15 m/year horizontally and 50.66 m/year vertically. The southern tributary underwent two bends before merging with the mainstream, leading to a more complex motion pattern. Lastly, correlation reveals that the Zelongnong Glacier was affected by the combined influence of temperature and precipitation with a common period of around 90 days. Full article
(This article belongs to the Section Engineering Remote Sensing)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Geographical location of the study area. (<b>a</b>) Optical image of Zelongnong Glacier, with the Snout Located at 29.616N, 94.988E. (<b>b</b>) The location of Zelongnong valley and the coverage of the SAR datasets.</p>
Full article ">Figure 2
<p>Block diagram of the PO-MSBAS method.</p>
Full article ">Figure 3
<p>The illustration of the offset tracking method based on variable window-weighted cross-correlation.</p>
Full article ">Figure 4
<p>Three-dimensional velocity maps of Zelongnong Glacier. (<b>a</b>) Three-dimensional glacier velocity map, where the horizontal velocities are shown in vectors, and the vertical velocities are color-coded; (<b>b</b>) distribution of four feature points and four profiles along with their stages on the vertical velocity map. The feature points P1–P4 are represented by red pentagrams, while the profiles 1–4 are shown as red line segments. The stages are indicated by I–V.</p>
Full article ">Figure 5
<p>Three-dimensional glacier velocity and elevation along the profile, where the error bars indicate the standard deviation within 3 × 3 pixels.</p>
Full article ">Figure 6
<p>Time series of EW glacier motion from 2 June 2022 to 8 June 2023.</p>
Full article ">Figure 7
<p>Time series of NS glacier motion from 2 June 2022 to 8 June 2023.</p>
Full article ">Figure 8
<p>Time series of UD glacier motion from 2 June 2022 to 8 June 2023.</p>
Full article ">Figure 9
<p>Three-dimensional displacement time series for points P1–P4, where the positions are shown in <a href="#remotesensing-16-04462-f004" class="html-fig">Figure 4</a>b.</p>
Full article ">Figure 10
<p>The azimuth deformation was calculated using different methods: (<b>a</b>–<b>c</b>) the azimuth deformation calculated with different fixed window sizes; (<b>d</b>) the azimuth deformation calculated using the method proposed in this study.</p>
Full article ">Figure 11
<p>The standard deviation of azimuth deformation with traditional and new methods.</p>
Full article ">Figure 12
<p>The relationship between glacier velocity and temperature and precipitation. (<b>a</b>) Glacier velocity versus temperature at points P1–P2; (<b>b</b>) glacier velocity versus temperature at points P3–P4; (<b>c</b>) glacier velocity versus precipitation at points P1–P2; and (<b>d</b>) glacier velocity versus precipitation at points P3–P4.</p>
Full article ">Figure 12 Cont.
<p>The relationship between glacier velocity and temperature and precipitation. (<b>a</b>) Glacier velocity versus temperature at points P1–P2; (<b>b</b>) glacier velocity versus temperature at points P3–P4; (<b>c</b>) glacier velocity versus precipitation at points P1–P2; and (<b>d</b>) glacier velocity versus precipitation at points P3–P4.</p>
Full article ">Figure 13
<p>Wavelet energy spectrum of glacier velocity crossed with temperature and precipitation.</p>
Full article ">
Back to TopTop