Subways have been an important method for relieving traffic pressures in urban areas, but ground subsidence, during construction and operation, can be a serious problem as it may affect the safety of its operation and that of the surrounding buildings. Thus, conducting long-term ground deformation monitoring and modeling for subway networks are essential. Compared with traditional geodetic methods, the Persistent Scatterer Interferometric Synthetic Aperture Radar (PS-InSAR) technique offers wider coverage and denser measurements along subway lines. In this study, we mapped the surface deformation of the Guangzhou subway network with Advanced Synthetic Aperture Radar (ASAR) and Phased Array Type L-band Synthetic Aperture Radar (PALSAR) data using the Interferometric Point Target Analysis (IPTA) technique. The results indicate that newly excavated tunnels have regional subsidence with an average rate of more than 8 mm/year, as found on Lines Two, Three, Six, and GuangFo (GF). Furthermore, we determined the spatio-temporal subsidence behavior of subways with PALSAR in delta areas using Peck’s formula and the logistic time model. We estimated the tunneling-related parameters in soft soil areas, which had not been previously explored. We examined a section of line GF, as an example, to estimate the ground settlement trough development. The results showed the maximum settlement increased from −5.2 mm to −23.6 mm and its ground loss ratio ranged from 1.5–8.7% between 13 July 2008 and 19 January 2011. In addition, we found that the tunnels in line GF will become stable after a period of about 2300 days in peak subsidence areas. The results show that the proposed approach can help explain the dynamic ground subsidence along a metro line. This study can provide references for urban subway projects in delta areas, and for the risk assessment of nearby buildings and underground pipelines along metro lines.

Because of the advantages of low cost, large coverage and short revisit cycle, Landsat 8 images have been widely applied to monitor earth surface movements. However, there are few systematic studies considering the error source characteristics or the improvement of the deformation field accuracy obtained by Landsat 8 image. In this study, we utilize the 2013 Mw 7.7 Balochistan, Pakistan earthquake to analyze error spatio-temporal characteristics and elaborate how to mitigate error sources in the deformation field extracted from multi-temporal Landsat 8 images. We found that the stripe artifacts and the topographic shadowing artifacts are two major error components in the deformation field, which currently lack overall understanding and an effective mitigation strategy. For the stripe artifacts, we propose a small spatial baseline (<200 m) method to avoid the stripe artifacts effect on the deformation field. We also propose a small radiometric baseline method to reduce the topographic shadowing artifacts and radiometric decorrelation noises. Those performances and accuracy evaluation show that these two methods are effective in improving the precision of deformation field. This study provides the possibility to detect subtle ground movement with higher precision caused by earthquake, melting glaciers, landslides, etc., with Landsat 8 images. It is also a good reference for error source analysis and corrections in deformation field extracted from other optical satellite images.

This paper discusses the potential and limitations of high-resolution P-band polarimetric synthetic aperture radar (SAR) interferometry (PolInSAR) in underlying topography estimation over forest areas. Time-frequency (TF) analysis in the azimuth direction is utilized to separate the ground scattering contribution from the total PolInSAR signal, without the use of any physical model, because the P-band PolInSAR data have a significant penetration depth and sufficient observation angle interval. To achieve this goal, a one-dimensional polynomial fitting (PF) method is proposed for correcting the residual motion error (RME). The Krycklan catchment test site, which is covered with pine forest, was selected to test the performance of the digital elevation model (DEM) inversion. The results show that the PF method can correct the RMEs for the sub-look interferograms well. When compared to the existing line-fit method, the TF+PF method can provide a more accurate DEM (the accuracy is improved by 26.9%). Moreover, the performance of the DEM inversion is free from the random-volume-over-ground assumption.

Forests are the most important part of terrestrial ecosystems. In the context of China’s industrialization and urbanization, mining activities have caused huge damage to the forest ecology. In the Ulan Mulun River Basin (Ordos, China), afforestation is standard method for reclamation of coal mine degraded land. In order to understand, manage and utilize forests, it is necessary to collect local mining area’s tree information. This paper proposed an improved Faster R-CNN model to identify individual trees. There were three major improved parts in this model. First, the model applied supervised multi-policy data augmentation (DA) to address the unmanned aerial vehicle (UAV) sample label size imbalance phenomenon. Second, we proposed Dense Enhance Feature Pyramid Network (DE-FPN) to improve the detection accuracy of small sample. Third, we modified the state-of-the-art Alpha Intersection over Union (Alpha-IoU) loss function. In the regression stage, this part effectively improved the bounding box accuracy. Compared with the original model, the improved model had the faster effect and higher accuracy. The result shows that the data augmentation strategy increased AP by 1.26%, DE-FPN increased AP by 2.82%, and the improved Alpha-IoU increased AP by 2.60%. Compared with popular target detection algorithms, our improved Faster R-CNN algorithm had the highest accuracy for tree detection in mining areas. AP was 89.89%. It also had a good generalization, and it can accurately identify trees in a complex background. Our algorithm detected correct trees accounted for 91.61%. In the surrounding area of coal mines, the higher the stand density is, the smaller the remote sensing index value is. Remote sensing indices included Green Leaf Index (GLI), Red Green Blue Vegetation Index (RGBVI), Visible Atmospheric Resistance Index (VARI), and Normalized Green Red Difference Index (NGRDI). In the drone zone, the western area of Bulianta Coal Mine (Area A) had the highest stand density, which was 203.95 trees ha⁻¹. GLI mean value was 0.09, RGBVI mean value was 0.17, VARI mean value was 0.04, and NGRDI mean value was 0.04. The southern area of Bulianta Coal Mine (Area D) was 105.09 trees ha⁻¹ of stand density. Four remote sensing indices were all the highest. GLI mean value was 0.15, RGBVI mean value was 0.43, VARI mean value was 0.12, and NGRDI mean value was 0.09. This study provided a sustainable development theoretical guidance for the Ulan Mulun River Basin. It is crucial information for local ecological environment and economic development.

The TerraSAR-X add-on for Digital Elevation Measurements (TanDEM-X) bistatic system provides high-resolution and high-quality interferometric data for global topographic measurement. Since the twin TanDEM-X satellites fly in a close helix formation, they can acquire approximately simultaneous synthetic aperture radar (SAR) images, so that temporal decorrelation and atmospheric delay can be ignored. Consequently, the orbital error becomes the most significant error limiting high-resolution SAR interferometry (InSAR) applications, such as the high-precision digital elevation model (DEM) reconstruction, subway and highway deformation monitoring, landslide monitoring and sub-canopy topography inversion. For rugged mountainous areas, in particular, it is difficult to estimate and correct the orbital phase error in TanDEM-X bistatic InSAR. Based on the rigorous InSAR geometric relationship, the orbital phase error can be attributed to the baseline errors (BEs) after fixing the positions of the master SAR sensor and the targets on the ground surface. For the constraint of the targets at a study scene, the freely released TanDEM-X DEM can be used, due to its consistency with the TanDEM-X bistatic InSAR-measured height. As a result, a parameterized model for the orbital phase error estimation is proposed in this paper. In high-resolution and high-precision TanDEM-X bistatic InSAR processing, due to the limited precision of the navigation systems and the uneven baseline changes caused by the helix formation, the BEs are time-varying in most cases. The parameterized model is thus built and estimated along each range line. To validate the proposed method, two mountainous test sites located in China (i.e., Fuping in Shanxi province and Hetang in Hunan province) were selected. The obtained results show that the orbital phase errors of the bistatic interferograms over the two test sites are well estimated. Compared with the widely applied polynomial model, the residual phase corrected by the proposed method contains little undesirable topography-dependent phase error, and avoids unexpected height errors ranging about from −6 m to 3 m for the Fuping test site and from −10 m to 8 m for the Hetang test site. Furthermore, some fine details, such as ridges and valleys, can be clearly identified after the correction. In addition, the two components of the orbital phase error, i.e., the residual flat-earth phase error and the topographic phase error caused by orbital error, are separated and quantified based on the parameterized expression. These demonstrate that the proposed method can be used to accurately estimate and mitigate the orbital phase error in TanDEM-X bistatic InSAR data, which increases the feasibility of reconstructing high-resolution and high-precision DEM. The rigorous geometric constraint, the refinement of the initial baseline parameters, and the assessment for height errors based on the estimated BEs are investigated in the discussion section of this paper.

To address the challenge of retrieving sub-canopy topography using single-baseline single-polarization TanDEM-X InSAR data, we propose a novel InSAR processing framework. Our methodology begins by employing the SINC model to estimate the penetration depth (PD). Subsequently, we establish a linear relationship between PD and phase center height (PCH) to generate a wall-to-wall PCH product. To achieve this, space-borne LiDAR data are employed to capture the elevation bias between actual ground elevation and InSAR-derived elevation. Finally, the sub-canopy topography is derived by subtracting the PCH from the conventional InSAR-based DEM. Moreover, this approach enables the simultaneous estimation of forest height from single-baseline TanDEM-X data by combining the estimated PD and PCH components. The approach has been validated against Airborne Lidar Scanning data over four diverse sites encompassing different forest types, terrain conditions, and climates. The derived sub-canopy topography in the boreal and hemi-boreal forest sites (Krycklan and Remningstorp) demonstrated notable improvement in accuracy. Additionally, the winter acquisitions outperformed the summer ones in terms of inversion accuracy. The achieved RMSEs for the winter scenarios were 2.45 m and 3.83 m, respectively, representing a 50% improvement over the InSAR-based DEMs. And the forest heights are also close to the ALS measurements, with RMSEs of 2.70 m and 3.33 m, respectively. For the Yanguas site in Spain, characterized by rugged terrain, sub-canopy topography in forest areas was estimated with an accuracy of 4.27m, a 35% improvement over the original DEM. For the denser tropical forest site, only an average elevation bias could be corrected.