On February 6, 2023 (local time), two earthquakes (Mw7.8 and Mw7.7) struck central and southern Turkey, causing extensive damage to several cities and claiming a toll of 40,000 lives. In this study, we propose a method for seismic building damage assessment and analysis by combining SAR amplitude and phase coherence change detection. We determined building damage in five severely impacted urban areas and calculated the damage ratio by measuring the urban area and the damaged area. The largest damage ratio of 18.93% is observed in Nurdagi, and the smallest ratio of 7.59% is found in Islahiye. We verified the results by comparing them with high-resolution optical images and AI recognition results from the Microsoft team. We also used pixel offset tracking (POT) technology and D-InSAR technology to obtain surface deformation using Sentinel-1A images and analyzed the relationship between surface deformation and post-earthquake urban building damage. The results show that Nurdagi has the largest urban average surface deformation of 0.48 m and Antakya has the smallest deformation of 0.09 m. We found that buildings in the areas with steeper slopes or closer to earthquake faults have higher risk of collapse. We also discussed the influence of SAR image parameters on building change recognition. Image resolution and observation geometry have a great influence on the change detection results, and the resolution can be improved by various means to raise the recognition accuracy. Our research findings can guide earthquake disaster assessment and analysis and identify influential factors of earthquake damage.
In February 2023, Mw 7.8 and Mw 7.7 earthquakes struck southeastern Turkey. Generating a coseismic 3D deformation field that can directly reflect the characteristics of surface deformation is important for revealing the movement mode of a seismogenic fault and analyzing the focal mechanism. Optical image sub-pixel correlation (SPC) only captures deformation in the horizontal direction, and SAR image pixel offset tracking (POT) obtains range deformation that is not sensitive to north–south deformation signals. Thus, neither of them can capture the complete 3D deformation alone. Combining them may be able to allow the monitoring of 3D deformation. In this study, we used Sentinel-2 optical images to obtain the horizontal deformation (east–west and north–south) and Sentinel-1 and ALOS-2 data to extract the range and azimuth offsets. The least-squares method was used to fuse the optical and SAR offsets to obtain the 3D deformation field of the 2023 Turkey earthquake sequence, which indicates that the two events were both left-lateral strike-slip earthquakes. The surface deformation caused by the two large earthquakes is mainly in the east–west direction. In the vertical direction, the two earthquakes caused a small-magnitude uplift and subsidence. The findings in this paper can be used as a reference for the study of coseismic 3D deformation.
The joint action of human activities and environmental changes contributes to the frequent occurrence of landslide, causing major hazards. Using Interferometric Synthetic Aperture Radar (InSAR) technique enables the detailed detection of surface deformation, facilitating early landslide detection. The growing availability of SAR data and the development of artificial intelligence have spurred the integration of deep learning methods with InSAR for intelligent geological identification. However, existing studies using deep learning methods to detect landslides in InSAR deformation often rely on single InSAR data, which leads to the presence of other types of geological hazards in the identification results and limits the accuracy of landslide identification. Landslides are affected by many factors, especially topographic features. To enhance the accuracy of landslide identification, this study improves the existing geological hazard detection model and proposes a multi-source data fusion network termed MSFD-Net. MSFD-Net employs a pseudo-Siamese network without weight sharing, enabling the extraction of texture features from the wrapped deformation data and topographic features from topographic data, which are then fused in higher-level feature layers. We conducted comparative experiments on different networks and ablation experiments, and the results show that the proposed method achieved the best performance. We applied our method to the middle and upper reaches of the Yellow River in eastern Qinghai Province, China, and obtained deformation rates using Sentinel-1 SAR data from 2018 to 2020 in the region, ultimately identifying 254 landslides. Quantitative evaluations reveal that most detected landslides in the study area occurred at an elevation of 2500–3700 m with slope angles of 10–30°. The proposed landslide detection algorithm holds significant promise for quickly and accurately detecting wide-area landslides, facilitating timely preventive and control measures.
In recent years, increasing available synthetic aperture radar (SAR) satellite data and gradually developing interferometric SAR (InSAR) technology have provided the possibility for wide-scale ground-deformation monitoring using InSAR. Traditionally, the InSAR data are processed by the existing time-series InSAR (TS–InSAR) technology, which has inefficient calculation and redundant results. In this study, we propose a wide-area InSAR variable-scale deformation detection strategy (hereafter referred to as the WAVS–InSAR strategy). The strategy combines stacking technology for fast ground-deformation rate calculation and advanced TS–InSAR technology for obtaining fine deformation time series. It adopts an adaptive recognition algorithm to identify the spatial distribution and area of deformation regions (regions of interest, ROI) in the wide study area and uses a novel wide-area deformation product organization structure to generate variable-scale deformation products. The Turpan–Hami basin in western China is selected as the wide study area (277,000 km2) to verify the proposed WAVS–InSAR strategy. The results are as follows: (1) There are 32 deformation regions with an area of ≥1 km2 and a deformation magnitude of greater than ±2 cm/year in the Turpan–Hami basin. The deformation area accounts for 2.4‰ of the total monitoring area. (2) A large area of ground subsidence has occurred in the farmland areas of the ROI, which is caused by groundwater overexploitation. The popularization and application of facility agriculture in the ROI have increased the demand for irrigation water. Due to the influence of the tectonic fault, the water supply of the ROI is mainly dependent on groundwater. Huge water demand has led to a continuous net deficit in aquifers, leading to land subsidence. The WAVS–InSAR strategy will be helpful for InSAR deformation monitoring at a national/regional scale and promoting the engineering application of InSAR technology.
The Sentinel-1A/B satellite launched by European Space Agency (ESA) in 2014 provides a huge amount of free Terrain Observation by Progressive Scans (TOPS) data with global coverage to the public. The TOPS data have a frame width of 250 km and have been widely used in surface deformation monitoring. However, traditional Multi-Temporal Interferometric Synthetic Aperture Radar (MT-InSAR) methods require large computer memory and time when processing full resolution data with large width and long strips. In addition, they hardly correct atmospheric delays and orbital errors accurately over a large area. In order to solve these problems, this study proposes a data partition strategy based on MT-InSAR methods. We first process the partitioned images over a large area by traditional MT-InSAR method, then stitch the deformation results into a complete deformation result by correcting the offsets of adjacent partitioned images. This strategy is validated in a flat urban area (Changzhou City in Jiangsu province, China), and a mountainous region (Qijiang in Chongqing City, China). Compared with traditional MT-InSAR methods, the precision of the results obtained by the new strategy is improved by about 5% for Changzhou city and about 15% for Qijiang because of its advantage in atmospheric delay correction. Furthermore, the proposed strategy needs much less memory and time than traditional methods. The total time needed by the traditional method is about 20 h, and by the proposed method, is about 8.7 h, when the number of parallel processing is 5 in the Changzhou city case. The time will be further reduced when the number of parallel processes increases.
The optical image sub-pixel correlation (SPC) technique is an important method for monitoring large-scale surface deformation. RapidEye images, distinguished by their short revisit period and high spatial resolution, are crucial data sources for monitoring surface deformation. However, few studies have comprehensively analyzed the error sources and correction methods of the deformation field obtained from RapidEye images. We used RapidEye images without surface deformation to analyze potential errors in the offset fields. We found that the errors in RapidEye offset fields primarily consist of decorrelation noise, orbit error, and attitude jitter distortions. To mitigate decorrelation noise, the careful selection of offset pairs coupled with spatial filtering is essential. Orbit error can be effectively mitigated by the polynomial fitting method. To address attitude jitter distortions, we introduced a linear fitting approach that incorporated the coherence of attitude jitter. To demonstrate the performance of the proposed methods, we utilized RapidEye images to extract the coseismic displacement field of the 2019 Ridgecrest earthquake sequence. The two-dimensional (2D) offset field contained deformation signals extracted from two earthquakes, with a maximum offset of 2.8 m in the E-W direction and 2.4 m in the N-S direction. A comparison with GNSS observations indicates that, after error correction, the mean relative precision of the offset field improved by 92% in the E-W direction and by 89% in the N-S direction. This robust enhancement underscores the effectiveness of the proposed error correction methods for RapidEye data. This study sheds light on large-scale surface deformation monitoring using RapidEye images.