Iterative atmospheric phase screen compensation for near-real-time ground-based InSAR measurements over a mountainous slope

In this article, an atmospheric phase screen (APS) compensation algorithm for a near real-time ground-based interferometry synthetic aperture radar (GB-InSAR) over a mountainous area is investigated. A novel APS compensation scheme is proposed to compensate the fluctuated APS caused by a spatial 3-D inhomogeneous refractivity index distribution without any a priori knowledge of moving location. The proposed method simultaneously addresses to identify moving pixels by a criterion of absolute velocity estimated by the coherent pixels technique (CPT). The proposed method consists mainly of three steps: 1) the stratified APS compensation; 2) identification of moving pixel candidate; and 3) the residual APS [remained APS after 1)] compensation by Kriging interpolation. The steps mentioned above are iteratively applied in order to increase the accuracy of the whole process. In this framework, we develop the 2-D quadratic polynomial model of the refractivity index with respect to slant range and topographic height for modeling the stratified APS. Furthermore, a prediction of the residual APS is achieved by applying the intrinsic random function of order k (IRF-k) Kriging interpolation, taking into account the nonstationarity of the residual APS. We evaluate the proposed method using zero-baseline GB-differential InSAR (GB-DInSAR) data over a mountainous area located in Minami-Aso, Kumamoto, Japan, through the near real-time post-landslide measurement campaign

Atmospheric artifacts correction with a covariance-weighted linear model over mountainous regions

© 2018 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in any current or future media, including reprinting/republishing this material for advertising or promotional purposes,creating new collective works, for resale or redistribution to servers or lists, or reuse of any copyrighted component of this work in other works.Mitigating the atmospheric phase delay is one of the largest challenges faced by the differential synthetic aperture radar (SAR) interferometry community. Recently, many publications have studied correcting the stratified tropospheric phase delay by assuming a linear model between them and the topography. However, most of these studies have not considered the effect of turbulent atmospheric artifacts when adjusting the linear model to data. In this paper, we present an improved technique that minimizes the influence of the turbulent atmosphere in the model adjustment. In the proposed algorithm, the model is adjusted to the phase differences of pixels instead of using the unwrapped phase of each pixel. In addition, the different phase differences are weighted as a function of its atmospheric phase screen covariance estimated from an empirical variogram to reduce, in the model adjustment, the impact of pixel pairs with a significant turbulent atmosphere. The good performance of the proposed method has been validated with both the simulated and real Sentinel-1A SAR data in the mountainous area of Tenerife island, Spain.Peer ReviewedPostprint (author's final draft

Atmospheric artifacts correction for InSAR using empirical model and numerical weather prediction models

lnSAR has been proved its unprecedented ability and merits of monitoring ground deformation on large scale with centimeter to millimeter scale accuracy. However, several factors affect the reliability and accuracy of its applications. Among them, atmospheric artifacts due to spatial and temporal variations of atmosphere state often pose noise to interferograms. Therefore, atmospheric artifacts m itigalion remains one of the biggest challenges to be addressed in the In SAR community. State-of-the-art research works have revealed atmospheric artifacts can be partially compensated with empirical models, temporal-spatial filtering approach in lnSAR time series, pointwise GPS zenith path delay and numerical weather prediction models. In this thesis, firstly, we further develop a covariance weighted linear empirical model correction method. Secondly, a realistic LOS direction integration approach based on global reanalysis data is employed and comprehensively compared with the conventional method that integrates along zenith direction. Finally, the realistic integration method is applied to local WRF numerical forecast model data. l'vbreover, detailed comparisons between different global reanalysis data and local WRF model are assessed. In terms of empirical models correcting methods, many publications have studied correcting stratified tropospheric phase delay by assuming a linear model between them and topography. However, most of these studies ha\19 not considered the effect of turbulent atmospheric artefacts when adjusting the linear model to data. In this thesis, an improved technique that minimizes the influence of turbulent atmosphere in the model adjustment has been presented. In the proposed algorithm, the model is adjusted to the phase differences of pixels instead of using the unwrapped phase of each pixel. In addition, the different phase differences are weighted as a function of its APS covariance estimated from an empirical variogram to reduce in the model adjustment the impact of pixel pairs with significant turbulent atmosphere. The performance of the proposed method has been validated with both simulated and real Sentinel-1 SAR data in Tenerife island, Spain. Considering methods using meteorological observations to mitigate APS, an accurate realistic com puling strategy utilizing global atmospheric reanalysis data has been implemented. With the approach, the realistic LOS path along satellite and the monitored points is considered, rather than converting from zenith path delay. Com pared with zenith delay based method, the biggest advantage is that it can avoid errors caused by anisotropic atmospheric behaviour. The accurate integration method is validated with Sentinel-1 data in three test sites: Tenerife island, Spain, Almeria, Spain and Crete island, Greece. Compared to conventional zenith method, the realistic integration method shows great improvement. A variety of global reanalysis data are available from different weather forecasting organizations, such as ERA-Interim, ERAS, MERRA2. In this study, the realistic integration mitigation method is assessed on these different reanalysis data. The results show that these data are feasible to mitigate APS to some extent in most cases. The assessment also demonstrates that the ERAS performs the best statistically, compared to other global reanalysis data. l'vbreover, as local numerical weather forecast models have the ability to predict high spatial resolution atmospheric parameters, by using which, it has the potential to achieve APS mitigation. In this thesis, the realistic integration method is also employed on the local WRF model data in Tenerife and Almeria test s ites. However, it turns out that the WRF model performs worse than the original global reanalysis data.Las técnicas lnSAR han demostrado su capacidad sin precedentes y méritos para el monitoreo de la deformaci6n del suelo a gran escala con una precisión centimétrica o incluso milimétrica. Sin embargo, varios factores afectan la fiabilidad y precisión de sus aplicaciones. Entre ellos, los artefactos atmosféricos debidos a variaciones espaciales y temporales del estado de la atm6sfera a menudo añaden ruido a los interferogramas. Por lo tanto, la mitigación de los artefactos atmosféricos sigue siendo uno de los mayores desafíos a abordar en la comunidad lnSAR. Los trabajos de investigaci6n de vanguardia han revelado que los artefactos atmosféricos se pueden compensar parcialmente con modelos empíricos, enfoque de filtrado temporal-espacial en series temporales lnSAR, retardo puntual del camino cenital con GPS y modelos numéricos de predicción meteorológica. En esta tesis, en primer lugar, desarrollamos un método de corrección de modelo empírico lineal ponderado por covarianza. En segundo lugar, se emplea un enfoque realista de integracion de dirección LOS basado en datos de reanálisis global y se compara exhaustivamente con el método convencional que se integra a lo largo de la dirección cenital. Finalmente, el método de integraci6n realista se aplica a los datos del modelo de pronóstico numérico WRF local. Ademas, se evalúan las comparaciones detalladas entre diferentes datos de reanálisis global y el modelo WRF local. En términos de métodos de corrección con modelos empíricos, muchas publicaciones han estudiado la corrección del retraso estratificado de la fase troposférica asumiendo un modelo lineal entre ellos y la topografía. Sin embargo, la mayoría de estos estudios no han considerado el efecto de los artefactos atmosféricos turbulentos al ajustar el modelo lineal a los datos. En esta tesis, se ha presentado una técnica mejorada que minimiza la influencia de la atm6sfera turbulenta en el ajuste del modelo. En el algoritmo propuesto, el modelo se ajusta a las diferencias de fase de los pixeles en lugar de utilizar la fase sin desenrollar de cada pixel. Además, las diferentes diferencias de fase se ponderan en función de su covarianza APS estimada a partir de un variograma empírico para reducir en el ajuste del modelo el impacto de los pares de pixeles con una atm6sfera turbulenta significativa. El rendimiento del método propuesto ha sido validado con datos SAR Sentinel-1 simulados y reales en la isla de Tenerife, España. Teniendo en cuenta los métodos que utilizan observaciones meteorológicas para mitigar APS, se ha implementado una estrategia de computación realista y precisa que utiliza datos de reanálisis atmosférico global. Con el enfoque, se considera el camino realista de LOS a lo largo del satélite y los puntos monitoreados, en lugar de convertirlos desde el retardo de la ruta cenital. En comparación con el método basado en la demora cenital, la mayor ventaja es que puede evitar errores causados por el comportamiento atmosférico anisotrópico. El método de integración preciso se valida con los datos de Sentinel-1 en tres sitios de prueba: la isla de Tenerife, España, Almería, España y la isla de Creta, Grecia. En comparación con el método cenital convencional, el método de integración realista muestra una gran mejora

Atmospheric artifacts correction for InSAR using empirical model and numerical weather prediction models

lnSAR has been proved its unprecedented ability and merits of monitoring ground deformation on large scale with centimeter to millimeter scale accuracy. However, several factors affect the reliability and accuracy of its applications. Among them, atmospheric artifacts due to spatial and temporal variations of atmosphere state often pose noise to interferograms. Therefore, atmospheric artifacts m itigalion remains one of the biggest challenges to be addressed in the In SAR community. State-of-the-art research works have revealed atmospheric artifacts can be partially compensated with empirical models, temporal-spatial filtering approach in lnSAR time series, pointwise GPS zenith path delay and numerical weather prediction models. In this thesis, firstly, we further develop a covariance weighted linear empirical model correction method. Secondly, a realistic LOS direction integration approach based on global reanalysis data is employed and comprehensively compared with the conventional method that integrates along zenith direction. Finally, the realistic integration method is applied to local WRF numerical forecast model data. l'vbreover, detailed comparisons between different global reanalysis data and local WRF model are assessed. In terms of empirical models correcting methods, many publications have studied correcting stratified tropospheric phase delay by assuming a linear model between them and topography. However, most of these studies ha\19 not considered the effect of turbulent atmospheric artefacts when adjusting the linear model to data. In this thesis, an improved technique that minimizes the influence of turbulent atmosphere in the model adjustment has been presented. In the proposed algorithm, the model is adjusted to the phase differences of pixels instead of using the unwrapped phase of each pixel. In addition, the different phase differences are weighted as a function of its APS covariance estimated from an empirical variogram to reduce in the model adjustment the impact of pixel pairs with significant turbulent atmosphere. The performance of the proposed method has been validated with both simulated and real Sentinel-1 SAR data in Tenerife island, Spain. Considering methods using meteorological observations to mitigate APS, an accurate realistic com puling strategy utilizing global atmospheric reanalysis data has been implemented. With the approach, the realistic LOS path along satellite and the monitored points is considered, rather than converting from zenith path delay. Com pared with zenith delay based method, the biggest advantage is that it can avoid errors caused by anisotropic atmospheric behaviour. The accurate integration method is validated with Sentinel-1 data in three test sites: Tenerife island, Spain, Almeria, Spain and Crete island, Greece. Compared to conventional zenith method, the realistic integration method shows great improvement. A variety of global reanalysis data are available from different weather forecasting organizations, such as ERA-Interim, ERAS, MERRA2. In this study, the realistic integration mitigation method is assessed on these different reanalysis data. The results show that these data are feasible to mitigate APS to some extent in most cases. The assessment also demonstrates that the ERAS performs the best statistically, compared to other global reanalysis data. l'vbreover, as local numerical weather forecast models have the ability to predict high spatial resolution atmospheric parameters, by using which, it has the potential to achieve APS mitigation. In this thesis, the realistic integration method is also employed on the local WRF model data in Tenerife and Almeria test s ites. However, it turns out that the WRF model performs worse than the original global reanalysis data.Las técnicas lnSAR han demostrado su capacidad sin precedentes y méritos para el monitoreo de la deformaci6n del suelo a gran escala con una precisión centimétrica o incluso milimétrica. Sin embargo, varios factores afectan la fiabilidad y precisión de sus aplicaciones. Entre ellos, los artefactos atmosféricos debidos a variaciones espaciales y temporales del estado de la atm6sfera a menudo añaden ruido a los interferogramas. Por lo tanto, la mitigación de los artefactos atmosféricos sigue siendo uno de los mayores desafíos a abordar en la comunidad lnSAR. Los trabajos de investigaci6n de vanguardia han revelado que los artefactos atmosféricos se pueden compensar parcialmente con modelos empíricos, enfoque de filtrado temporal-espacial en series temporales lnSAR, retardo puntual del camino cenital con GPS y modelos numéricos de predicción meteorológica. En esta tesis, en primer lugar, desarrollamos un método de corrección de modelo empírico lineal ponderado por covarianza. En segundo lugar, se emplea un enfoque realista de integracion de dirección LOS basado en datos de reanálisis global y se compara exhaustivamente con el método convencional que se integra a lo largo de la dirección cenital. Finalmente, el método de integraci6n realista se aplica a los datos del modelo de pronóstico numérico WRF local. Ademas, se evalúan las comparaciones detalladas entre diferentes datos de reanálisis global y el modelo WRF local. En términos de métodos de corrección con modelos empíricos, muchas publicaciones han estudiado la corrección del retraso estratificado de la fase troposférica asumiendo un modelo lineal entre ellos y la topografía. Sin embargo, la mayoría de estos estudios no han considerado el efecto de los artefactos atmosféricos turbulentos al ajustar el modelo lineal a los datos. En esta tesis, se ha presentado una técnica mejorada que minimiza la influencia de la atm6sfera turbulenta en el ajuste del modelo. En el algoritmo propuesto, el modelo se ajusta a las diferencias de fase de los pixeles en lugar de utilizar la fase sin desenrollar de cada pixel. Además, las diferentes diferencias de fase se ponderan en función de su covarianza APS estimada a partir de un variograma empírico para reducir en el ajuste del modelo el impacto de los pares de pixeles con una atm6sfera turbulenta significativa. El rendimiento del método propuesto ha sido validado con datos SAR Sentinel-1 simulados y reales en la isla de Tenerife, España. Teniendo en cuenta los métodos que utilizan observaciones meteorológicas para mitigar APS, se ha implementado una estrategia de computación realista y precisa que utiliza datos de reanálisis atmosférico global. Con el enfoque, se considera el camino realista de LOS a lo largo del satélite y los puntos monitoreados, en lugar de convertirlos desde el retardo de la ruta cenital. En comparación con el método basado en la demora cenital, la mayor ventaja es que puede evitar errores causados por el comportamiento atmosférico anisotrópico. El método de integración preciso se valida con los datos de Sentinel-1 en tres sitios de prueba: la isla de Tenerife, España, Almería, España y la isla de Creta, Grecia. En comparación con el método cenital convencional, el método de integración realista muestra una gran mejora.Postprint (published version

Multi-technique approach to rockfall monitoring in the Montserrat massif (Catalonia, NE Spain)

Montserrat Mountain is located near Barcelona in Catalonia, in the northeast of Spain, and its massif is formed by conglomerate interleaved by siltstone/sandstone with steep slopes very prone to rockfalls. The increasing number of visitors in the monastery area, reaching 2.4 million per year, has highlighted the risk derived from rockfalls for this building area and also for the terrestrial accesses, both roads and the rack railway. A risk mitigation plan has been launched, and its first phase during 2014-2016 has been focused largely on testing several monitoring techniques for their later implementation. The results of the pilot tests, performed as a development from previous sparse experiences and data, are presented together with the first insights obtained. These tests combine four monitoring techniques under different conditions of continuity in space and time domains, which are: displacement monitoring with Ground-based Synthetic Aperture Radar and characterization at slope scale, with an extremely non-uniform atmospheric phase screen due to the stepped topography and atmosphere stratification; Terrestrial Laser Scanner surveys quantifying the frequency of small or even previously unnoticed rockfalls, and monitoring rock block centimetre scale displacements; the monitoring of rock joints implemented through a wireless sensor network with an ad hoc design of ZigBee loggers developed by ICGC; and, finally, monitoring singular rock needles with Total Station.Peer ReviewedPostprint (author's final draft

Ground-based polarimetric SAR interferometry for the monitoring of terrain displacement phenomena-part II: applications

©2014 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in any current or future media, including reprinting/republishing this material for advertising or promotional purposes, creating new collective works, for resale or redistribution to servers or lists, or reuse of any copyrighted component of this work in other works.Urban subsidence and landslides are among the greatest hazards for people and infrastructure safety and they require an especial attention to reduce their associated risks. In this framework, ground-based synthetic aperture radar (SAR) interferometry (GB-InSAR) represents a cost-effective solution for the precise monitoring of displacements. This work presents the application of GB-InSAR techniques, particularly with the RiskSAR sensor and its processing chain developed by the Remote Sensing Laboratory (RSLab) of the Universitat Politecnica de Catalunya (UPC), for the monitoring of two different types of ground displacement. An example of urban subsidence monitoring over the village of Sallent, northeastern of Spain, and an example of landslide monitoring in El Forn de Canillo, located in the Andorran Pyrenees, are presented. In this framework, the key processing particularities for each case are deeply analyzed and discussed. The linear displacement maps and time series for both scenarios are showed and compared with in-field data. For the study, fully polarimetric data acquired at X-band with a zero-baseline configuration are employed in both scenarios. The displacement results obtained demonstrate the capabilities of GB-SAR sensors for the precise monitoring of ground displacement phenomena.Peer ReviewedPostprint (author's final draft

Temporal Coherence Estimators for GBSAR

Many Ground-Based Synthetic Aperture Radar (GBSAR) applications demand preliminary analysis to select areas with high-quality signal. That is, areas in which the phase can be processed to extract the desired information. The interferometric coherence and the amplitude dispersion index are important tools widely used in the literature to assess the quality of GBSAR images. So far, no direct relation has been found between the two. Indeed, they are parameters of different natures: amplitude dispersion index is calculated with only amplitude values, while coherence provides information also on the signal phase. The purpose of this article is to find a relation between the two parameters. Indeed, the amplitude dispersion index provides some practical advantages if compared to coherence estimators, especially to perform fast preliminary analysis. In this article, a theoretical relation between amplitude dispersion index and coherence is retrieved. GBSAR measurements acquired in different scenarios, at different working frequencies are presented and used to validate such a relation

Post-failure evolution analysis of a rainfall-triggered landslide by multi-temporal interferometry SAR approaches integrated with geotechnical analysis

Persistent Scatterers Interferometry (PSI) represents one of the most powerful techniques for Earth's surface deformation processes' monitoring, especially for long-term evolution phenomena. In this work, a dataset of 34 TerraSAR-X StripMap images (October 2013–October 2014) has been processed by two PSI techniques - Coherent Pixel Technique-Temporal Sublook Coherence (CPT-TSC) and Small Baseline Subset (SBAS) - in order to study the evolution of a slow-moving landslide which occurred on February 23, 2012 in the Papanice hamlet (Crotone municipality, southern Italy) and induced by a significant rainfall event (185 mm in three days). The mass movement caused structural damage (buildings' collapse), and destruction of utility lines (gas, water and electricity) and roads. The results showed analogous displacement rates (30–40 mm/yr along the Line of Sight – LOS-of the satellite) with respect to the pre-failure phase (2008–2010) analyzed in previous works. Both approaches allowed detect the landslide-affected area, however the higher density of targets identified by means of CPT-TSC enabled to analyze in detail the slope behavior in order to design possible mitigation interventions. For this aim, a slope stability analysis has been carried out, considering the comparison between groundwater oscillations and time-series of displacement. Hence, the crucial role of the interaction between rainfall and groundwater level has been inferred for the landslide triggering. In conclusion, we showed that the integration of geotechnical and remote sensing approaches can be seen as the best practice to support stakeholders to design remedial works.Peer ReviewedPostprint (author's final draft