Advances in synthetic aperture radar (SAR) interferometry have enabled the seamless monitoring of the Earth’s crust deformation. The dense archive of the Sentinel-1 Copernicus mission provides unprecedented spatial and temporal coverage; however, time-series analysis of such big data volumes requires high computational efficiency. We present a parallelized-PSI (P-PSI), a novel, parallelized, and end-to-end processing chain for the fully automated assessment of line-of-sight ground velocities through persistent scatterer interferometry (PSI), tailored to scale to the vast multitemporal archive of Sentinel-1 data. P-PSI is designed to transparently access different and complementary Sentinel-1 repositories, and download the appropriate datasets for PSI. To make it efficient for large-scale applications, we re-engineered and parallelized interferogram creation and multitemporal interferometric processing, and introduced distributed implementations to best use computing cores and provide resourceful storage management. We propose a new algorithm to further enhance the processing efficiency, which establishes a non-uniform patch grid considering land use, based on the expected number of persistent scatterers. P-PSI achieves an overall speed-up by a factor of five for a full Sentinel-1 frame for processing in a 20-core server. The processing chain is tested on a large-scale project to calculate and monitor deformation patterns over the entire extent of the Greek territory—our own Interferometric SAR (InSAR) Greece project. Time-series InSAR analysis was performed on volumes of about 12 TB input data corresponding to more than 760 Single Look Complex Sentinel-1A and B images mostly covering mainland Greece in the period of 2015–2019. InSAR Greece provides detailed ground motion information on more than 12 million distinct locations, providing completely new insights into the impact of geophysical and anthropogenic activities at this geographic scale. This new information is critical to enhancing our understanding of the underlying mechanisms, providing valuable input into risk assessment models. We showcase this through the identification of various characteristic geohazard locations in Greece and discuss their criticality. The selected geohazard locations, among a thousand, cover a wide range of catastrophic events including landslides, land subsidence, and structural failures of various scales, ranging from a few hundredths of square meters up to the basin scale. The study enriches the large catalog of geophysical related phenomena maintained by the GeObservatory portal of the Center of Earth Observation Research and Satellite Remote Sensing BEYOND of the National Observatory of Athens for the opening of new knowledge to the wider scientific community.
The exploitation of remote sensing techniques has substantially improved pre- and post- disaster landslide management over the last decade. A variety of landslide susceptibility methods exists, with capabilities and limitations related to scale and spatial accuracy issues, as well as data availability. The Interferometric Synthetic Aperture Radar (InSAR) capabilities have significantly contributed to the detection, monitoring, and mapping of landslide phenomena. The present study aims to point out the contribution of InSAR data in landslide detection and to evaluate two different scale landslide models by comparing a heuristic to a statistical method for the rainfall-induced landslide hazard assessment. Aiming to include areas with both high and low landslide occurrence frequencies, the study area covers a large part of the Aetolia–Acarnania and Evritania prefectures, Central and Western Greece. The landslide susceptibility product provided from the weights of evidence (WoE) method proved more accurate, benefitting from the expert opinion and the landslide inventory. On the other hand, the Norwegian Geological Institute (NGI) methodology has the edge on its immediate implementation, with minimum data requirements. Finally, it was proved that using sequential SAR image acquisitions gives the benefit of an updated landslide inventory, resulting in the generation of, on request, updated landslide susceptibility maps.