Multi-Temporal InSAR for transport infrastructure monitoring: Recent trends and challenges

Worldwide, transport infrastructure is increasingly vulnerable to ageing-induced deterioration and climate-related hazards. Oftentimes inspection and maintenance costs far exceed available resources, and numerous assets lack any rigorous structural evaluation. Space-borne Synthetic Aperture Radar Interferometry (InSAR) is a powerful remote-sensing technology, which can provide cheaper deformation measurements for bridges and other transport infrastructure with short revisit times, while scaling from the local to the global scale. As recent studies have shown the InSAR accuracy to be comparable with traditional monitoring instruments, InSAR could offer a cost-effective tool for long-term, near-continuous deformation monitoring, with the possibility to support inspection planning and maintenance prioritisation, while maximising functionality and increasing the resilience of infrastructure networks. However, despite the high potential of InSAR for structural monitoring, some important limitations need to be considered when applying it in reality. This paper identifies and discusses the challenges of using InSAR for the purpose of structural monitoring, with a specific focus on bridges and transport networks. Examples are presented to illustrate current practical limitations of InSAR; possible solutions and promising research directions are identified. The aim of this study is to motivate future action in this area and highlight the InSAR advances needed to overcome current challenges.

A Robust InSAR Phase Unwrapping Method via Phase Gradient Estimation Network

Phase unwrapping is a critical step in synthetic aperture radar interferometry (InSAR) data processing chains. In almost all phase unwrapping methods, estimating the phase gradient according to the phase continuity assumption (PGE-PCA) is an essential step. The phase continuity assumption is not always satisfied due to the presence of noise and abrupt terrain changes; therefore, it is difficult to get the correct phase gradient. In this paper, we propose a robust least squares phase unwrapping method that works via a phase gradient estimation network based on the encoder–decoder architecture (PGENet) for InSAR. In this method, from a large number of wrapped phase images with topography features and different levels of noise, the deep convolutional neural network can learn global phase features and the phase gradient between adjacent pixels, so a more accurate and robust phase gradient can be predicted than that obtained by PGE-PCA. To get the phase unwrapping result, we use the traditional least squares solver to minimize the difference between the gradient obtained by PGENet and the gradient of the unwrapped phase. Experiments on simulated and real InSAR data demonstrated that the proposed method outperforms the other five well-established phase unwrapping methods and is robust to noise.

Rapid characterisation of the extremely large landslide threatening the Rules Reservoir (Southern Spain)

When an active landslide is first identified in an artificial reservoir, a comprehensive study has to be quickly conducted to analyse the possible hazard that it may represent to such a critical infrastructure. This paper presents the case of the El Arrecife Landslide, located in a slope of the Rules Reservoir (Southern Spain), as an example of geological and motion data integration for elaborating a preliminary hazard assessment. For this purpose, a field survey was carried out to define the kinematics of the landslide: translational in favour of a specific foliation set, and rotational at the foot of the landslide. A possible failure surface has been proposed, as well as an estimation of the volume of the landslide: 14.7 million m3. At the same time, remote sensing and geophysical techniques were applied to obtain historical displacement rates. A mean subsidence rate of up to 2 cm/yr was obtained by means of Synthetic Aperture Radar Interferometry (InSAR) and Ground Penetrating Radar (GPR) data, during the last 5 and 22 years, respectively. The Structure-from-Motion (SfM) technique provided a higher rate, up to 26 cm/yr during the last 14 years, due to compaction of a slag heap located within the foot of the landslide. All of this collected information will be valuable to optimise the planning of future monitoring surveys (i.e. Differential Global Positioning Systems, inclinometers, ground drilling and InSAR) that should be applied in order to prevent further damage on the reservoir and related infrastructures.

Assessing Urban Landslide Dynamics through Multi-Temporal InSAR Techniques and Slope Numerical Modeling

Landslides threaten more than before the urbanized areas and are a worldwide growing problem for the already affected communities and the local authorities committed to landslide risk management and mitigation. For this reason, it is essential to analyze landslide dynamics and environmental conditioning factors. Various techniques and instruments exist for landslide investigation and monitoring. Out of these, Multi-temporal Synthetic Aperture Radar Interferometry (MT-InSAR) techniques have been widely used in the last decades. Their capabilities are enhanced by the availability of the active Sentinel-1 mission, whose 6-day revisiting time enables near real-time monitoring of landslides. Interferometric results, coupled with ground measurements or other approaches such as numerical modeling, significantly improve the knowledge of the investigated surface processes. In this work, we processed the C-band SAR images of the available European Space Agency (ESA) satellite missions, using MT-InSAR methods to identify the surface deformations related to landslides affecting the Iași Municipality (Eastern Romania). The results (i.e., velocity maps) point out the most active landslides with velocities of up to 20 mm/year measured along the satellite Line of Sight (LOS). Following, we focused on the most problematic landslide that affects the Țicău neighborhood and is well-known for its significant implications that it had. To better understand its behavior and the sensitivity of the displacements to the environmental factors (i.e., rainfall), we carried out 2D numerical modeling using a finite difference code. The simulated displacement field is consistent with the InSAR displacements and reveals the most active sectors of the landslide and insights about its mechanism.

The Role of Satellite InSAR for Landslide Forecasting: Limitations and Openings

The paper explores the potential of the satellite advanced differential synthetic aperture radar interferometry (A-DInSAR) technique for the identification of impending slope failure. The advantages and limitations of satellite InSAR in monitoring pre-failure landslide behaviour are addressed in five different case histories back-analysed using data acquired by different satellite missions: Montescaglioso landslide (2013, Italy), Scillato landslide (2015, Italy), Bingham Canyon Mine landslide (2013, UT, USA), Big Sur landslide (2017, CA, USA) and Xinmo landslide (2017, China). This paper aimed at providing a contribution to improve the knowledge within the subject area of landslide forecasting using monitoring data, in particular exploring the suitability of satellite InSAR for spatial and temporal prediction of large landslides. The study confirmed that satellite InSAR can be successful in the early detection of slopes prone to collapse; its limitations due to phase aliasing and low sampling frequency are also underlined. According to the results, we propose a novel landslide predictability classification discerning five different levels of predictability by satellite InSAR. Finally, the big step forward made for landslide forecasting applications since the beginning of the first SAR systems (ERS and Envisat) is shown, highlighting that future perspectives are encouraging thanks to the expected improvement of upcoming satellite missions that could highly increase the capability to monitor landslides’ pre-failure behaviour.

On the integration of multi-temporal synthetic aperture radar interferometry products and historical surveys data for buildings structural monitoring

The management and the safeguard of existing buildings and infrastructures are actual tasks for structural engineering. Non-invasive structural monitoring techniques can provide useful information for supporting the management process and the safety evaluation, reducing at once the impact of disturbances on the structure’s functionality. This paper focuses on the exploitation of advanced multi-temporal differential synthetic aperture radar interferometry (DInSAR) products for the structural monitoring of buildings and infrastructures, subjected to different external actions. In this framework, a methodological approach is proposed, based on the integration of DInSAR measurements with historical sources, accurate 3D modelling and consistent positioning of the reflecting targets in the GIS environment. Documentary sources can prove particularly helpful in collecting technical information, to reconstruct an accurate 3D geometry of the building under monitoring, limiting in-situ surveys. The analysis of DInSAR-based displacements time series and mean deformation velocity values allows the identification of possible critical situations for buildings to be monitored. The paper presents different approaches, with increasing accuracy levels, to study the active deformative processes of the examined buildings and the related damage assessment. An insight into these interpretative approaches is given through the application of the proposed procedure to two case studies in the city of Rome (Italy), the residential building named Torri Stellari in Valco San Paolo (1951–1953) and the housing complex referred to as Corviale (1967–1983), by exploiting the whole COSMO-SkyMed data archive (both ascending and descending acquisitions), collected during the 2011–2019 time interval. Pros and cons of the various approaches are deeply discussed, together with an estimation of the required computational effort.

A High-Resolution Antarctic Grounding Zone Product from ICESat-2 Laser Altimetry

The Antarctic grounding zone, which is the transition between the fully grounded ice sheet to freely floating ice shelf, plays a critical role in ice sheet instability, mass budget calculations and ice sheet model projections. It is therefore important to continuously monitor its location and migration over time. Here we present the first ICESat-2-derived high-resolution grounding zone product of the Antarctica Ice Sheet, including three important boundaries: the inland limit of tidal flexure (Point F), inshore limit of hydrostatic equilibrium (Point H) and the break-in-slope (Point Ib). This dataset was derived from automated techniques developed in this study, using ICESat-2 laser altimetry repeat tracks between 30 March 2019 and 30 September 2020. The new grounding zone product has a near complete coverage of the Antarctica Ice Sheet with a total of 21346 Point F, 18149 Point H and 36765 Point Ib identified, including the difficult to survey grounding zones, such as the fast-flowing glaciers draining into the Amundsen Sea Embayment. The locations of newly derived ICESat-2 landward limit of tidal flexure agree well with the most recent differential synthetic aperture radar interferometry (DInSAR) observations in 2018, with the mean absolute separation and standard deviation of 0.02 and 0.02 km, respectively. By comparing the ICESat-2-derived grounding zone with the previous grounding zone products, we find up-to 15 km grounding line retreat on the Crary Ice Rise of Ross Ice Shelf and the pervasive landward grounding line migration along the Amundsen Sea Embayment during the past two decades. We also identify the presence of ice plain on the Filchner-Ronne Ice Shelf and the influence of oscillating ocean tides on the grounding zone migration. The product derived from this study is available at https://doi.org/10.5523/bris.bnqqyngt89eo26qk8keckglww (Li et al., 2021) and is archived and maintained at the National Snow and Ice Data Center.

On the Assessment GPS-Based WRFDA for InSAR Atmospheric Correction: A Case Study in Pearl River Delta Region of China

The accuracy and applications of synthetic aperture radar interferometry (InSAR) are severely suppressed by tropospheric error. Numerical Weather Models (NWMs) and GPS-derived tropospheric delays have been widely used to correct the tropospheric error considering their complete spatial coverage or high accuracy. However, few studies focus on the fusion of both NWMs and GPS for the tropospheric error correction. In this study, we used the Weather Research and Forecasting (WRF) to obtain NWMs with a higher spatial-temporal resolution of 3 km and 20 s from both ERAI (79 km and 6 h) and ERA5 (0.25° and 1 h). After that, we utilized the WRF Data Assimilation (WRFDA) system to assimilate the GPS ZTD into these enhanced NWMs and generate merged NWMs products. The tropospheric correction effectiveness from different NWMs products was evaluated in a case in the Pearl River Delta region of China. The results showed that all the NWMs products could correct the stratified component in the interferogram but could not mitigate the turbulence well, even after improving the spatial-temporal resolution. As for the trend component, the merged NWMs products showed obvious superiority over other products. From the statistics perspective, the stdev of the interferogram decreased further over 20% by the merged NWMs products than other products when using both ERAI and ERA5, indicating the significant effectiveness of GPS ZTD assimilation.

New evidence for active talus-foot rock glaciers at Øyberget, southern Norway, and their development during the Holocene

Synthetic aperture radar interferometry (InSAR) measurements demonstrate that lobate, blocky depositional landforms, located in southern Norway at an altitude of ~530 m above sea level, with an estimated mean annual air temperature of ~1.6°C, currently exhibit deformation attributed to viscous creep. Five years of InSAR measurements for six lobes demonstrate average surface velocities of 1.2–22.0 mm/year with maximum rates of 17.5–55.6 mm/year. New Schmidt-hammer exposure-age dating (SHD) of two proximal lobes reveals mid-Holocene ages (7.6 ± 1.3 ka and 6.0 ± 1.2 ka), which contrast with the early-Holocene SHD and 10Be ages obtained previously from distal lobes, and late-Holocene SHD ages presented here from two adjacent talus slopes (2.3 ± 1.0 ka and 2.4 ± 1.0 ka). Although passive transport of boulders on the surfaces of these small, slow-moving rock glaciers affected by compressive flow means that the exposure ages are close to minimum estimates of the time elapsed since lobe inception, disturbance of boulders on rock glaciers is a source of potentially serious underestimates of rock-glacier age. Rock-glacier development at Øyberget began shortly after local deglaciation around 10 ka before present and continued throughout the Holocene in response to microclimatic undercooling within the coarse blocky surface layer of the talus and rock-glacier lobes. We suggest this enhanced cooling lowers mean annual surface-layer temperature by at least ~3.6°C, which is needed at such a low altitude to sustain sporadic permafrost and avoid fast thawing as atmospheric temperatures rise. Our results point to circumstances where inferences about rock glaciers as indicators of regional climate should be interpreted with caution, and where they may be less useful in palaeoclimatic reconstruction than previously thought.