An On-Demand Web Tool for the Unsupervised Retrieval of Earth’s Surface Deformation from SAR Data: The P-SBAS Service within the ESA G-POD Environment

Elastic interaction between Mauna Loa and Kīlauea evidenced by independent component analysis

The contrasting dynamics between Mauna Loa and Kīlauea have been studied over the last 100 years from multiple viewpoints. The fact that dynamic changes of one volcano trigger a dynamic response of the other volcano indicates a connection may exist. Petrological works show a direct relationship between the magmatic systems of these two volcanoes is not possible. We analysed DInSAR data and GPS measurements of ground deformation patterns associated with the activity of Mauna Loa and Kīlauea volcanoes. The DInSAR SBAS dataset spans the interval between 2003 and 2010, and was acquired along ascending and descending orbits of the ENVISAT (ESA) satellite under different look angles. Of the 10 tracks that cover the Big Island (Hawai‘i), 4 cover both volcanic edifices. Using GPS measurements, we computed the areal strain on 15 triplets of stations for Kīlauea volcano and 11 for Mauna Loa volcano. DInSAR data was analysed by applying Independent Component Analysis (ICA) to decompose the time-varying ground deformation pattern of both volcanoes. The results revealed anticorrelated ground deformation behaviour of the main calderas of Mauna Loa and Kīlauea, meaning that the opposite response is seen in the ground deformation of one volcano with respect to the other. At the same time, Kīlauea exhibits a more complex pattern, with an additional component, which appears not to be correlated with the dynamics of Mauna Loa. The GPS areal strain time series support these findings. To corroborate and help interpret the results, we performed inverse modelling of the observed ground deformation pattern using analytical source models. The results indicate that the ground deformation of Mauna Loa is associated with a dike-shaped source located at 6.2 km depth. In comparison, the anticorrelated ground deformation of Kīlauea is associated with a volumetric source at 1.2 km depth. This excludes a hydraulic connection as a possible mechanism to explain the anticorrelated behaviour; instead, we postulate a stress-transfer mechanism. To support this hypothesis, we performed a 3D numerical modelling of stress and strain fields in the study area, determining the elastic interaction of each source over the others. The most relevant finding is that the Mauna Loa shallow plumbing system can affect the shallowest magmatic reservoir of Kīlauea, while the opposite scenario is unlikely. Conversely, the second independent component observed at Kīlauea is associated to a sill-shaped source located at a depth of 3.5 km, which is less affected by this interaction process.

FLATSIM: The ForM@Ter LArge-Scale Multi-Temporal Sentinel-1 InterferoMetry Service

The purpose of the ForM@Ter LArge-scale multi-Temporal Sentinel-1 InterferoMetry service (FLATSIM) is the massive processing of Sentinel-1 data using multi-temporal interferometric synthetic aperture radar (InSAR) over large areas, i.e., greater than 250,000 km2. It provides the French ForM@ter scientific community with automatically processed products using a state of the art processing chain based on a small baseline subset approach, namely the New Small Baseline (NSBAS). The service results from a collaboration between the scientific team that develops and maintains the NSBAS processing chain and the French Spatial Agency (CNES) that mirrors the Sentinel-1 data. The proximity to Sentinel-1 data, the NSBAS workflow, and the specific optimizations to make NSBAS processing massively parallel for the CNES high performance computing infrastructure ensures the efficiency of the chain, especially in terms of input and output, which is the key for the success of such a service. The FLATSIM service is made of a production module, a delivery module and a user access module. Products include interferograms, surface line of sight velocity, phase delay time series and auxiliary data. Numerous quality indicators are provided for an in-depth analysis of the quality and limits of the results. The first national call in 2020 for region of interest ended up with 8 regions spread over the world with scientific interests, including seismology, tectonics, volcano-tectonics, and hydrological cycle. To illustrate the FLATSIM capabilities, an analysis is shown here on two processed regions, the Afar region in Ethiopa, and the eastern border of the Tibetan Plateau.

Sentinel-1 Big Data Processing with P-SBAS InSAR in the Geohazards Exploitation Platform: An Experiment on Coastal Land Subsidence and Landslides in Italy

The growing volume of synthetic aperture radar (SAR) imagery acquired by satellite constellations creates novel opportunities and opens new challenges for interferometric SAR (InSAR) applications to observe Earth’s surface processes and geohazards. In this paper, the Parallel Small BAseline Subset (P-SBAS) advanced InSAR processing chain running on the Geohazards Exploitation Platform (GEP) is trialed to process two unprecedentedly big stacks of Copernicus Sentinel-1 C-band SAR images acquired in 2014–2020 over a coastal study area in southern Italy, including 296 and 283 scenes in ascending and descending mode, respectively. Each stack was processed in the GEP in less than 3 days, from input SAR data retrieval via repositories, up to generation of the output P-SBAS datasets of coherent targets and their displacement histories. Use-cases of long-term monitoring of land subsidence at the Capo Colonna promontory (up −2.3 cm/year vertical and −1.0 cm/year east–west rate), slow-moving landslides and erosion landforms, and deformation at modern coastal protection infrastructure in the city of Crotone are used to: (i) showcase the type and precision of deformation products outputting from P-SBAS processing of big data, and the derivable key information to support value-adding and geological interpretation; and (ii) discuss potential and challenges of big data processing using cloud/grid infrastructure.

InSAR Greece with Parallelized Persistent Scatterer Interferometry: A National Ground Motion Service for Big Copernicus Sentinel-1 Data

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.

State of the Art and Recent Advancements in the Modelling of Land Subsidence Induced by Groundwater Withdrawal

Land subsidence is probably one of the most evident environmental effects of groundwater pumping. Globally, freshwater demand is the leading cause of this phenomenon. Land subsidence induced by aquifer system drainage can reach total values of up to 14.5 m. The spatial extension of this phenomenon is usually extensive and is often difficult to define clearly. Aquifer compaction contributes to many socio-economic effects and high infrastructure-related damage costs. Currently, many methods are used to analyze aquifer compaction. These include the fundamental relationship between groundwater head and groundwater flow direction, water pressure and aquifer matrix compressibility. Such solutions enable satisfactory modelling results. However, further research is needed to allow more efficient modelling of aquifer compaction. Recently, satellite radar interferometry (InSAR) has contributed to significant progress in monitoring and determining the spatio-temporal land subsidence distributions worldwide. Therefore, implementation of this approach can pave the way to the development of more efficient aquifer compaction models. This paper presents (1) a comprehensive review of models used to predict land surface displacements caused by aquifer drainage, as well as (2) recent advances, and (3) a summary of InSAR implementation in recent years to support the aquifer compaction modelling process.

Spatio-Temporal Assessment of Land Deformation as a Factor Contributing to Relative Sea Level Rise in Coastal Urban and Natural Protected Areas Using Multi-Source Earth Observation Data

The rise in sea level is expected to considerably aggravate the impact of coastal hazards in the coming years. Low-lying coastal urban centers, populated deltas, and coastal protected areas are key societal hotspots of coastal vulnerability in terms of relative sea level change. Land deformation on a local scale can significantly affect estimations, so it is necessary to understand the rhythm and spatial distribution of potential land subsidence/uplift in coastal areas. The present study deals with the determination of the relative vertical rates of the land deformation and the sea-surface height by using multi-source Earth observation—synthetic aperture radar (SAR), global navigation satellite system (GNSS), tide gauge, and altimetry data. To this end, the multi-temporal SAR interferometry (MT-InSAR) technique was used in order to exploit the most recent Copernicus Sentinel-1 data. The products were set to a reference frame by using GNSS measurements and were combined with a re-analysis model assimilating satellite altimetry data, obtained by the Copernicus Marine Service. Additional GNSS and tide gauge observations have been used for validation purposes. The proposed methodological approach has been implemented in three pilot cases: the city of Alexandroupolis in the Evros Delta region, the coastal zone of Thermaic Gulf, and the coastal area of Killini, Araxos (Patras Gulf) in the northwestern Peloponnese, which are Greek coastal areas with special characteristics. The present research provides localized relative sea-level estimations for the three case studies. Their variation is high, ranging from values close to zero, i.e., from 5–10 cm and 30 cm in 50 years for urban areas to values of 50–60 cm in 50 years for rural areas, close to the coast. The results of this research work can contribute to the effective management of coastal areas in the framework of adaptation and mitigation strategies attributed to climate change. Scaling up the proposed methodology to a continental level is required in order to overcome the existing lack of proper assessment of the relevant hazard in Europe.

The Status of Earth Observation Techniques in Monitoring High Mountain Environments at the Example of Pasterze Glacier, Austria: Data, Methods, Accuracies, Processes, and Scales

Earth observation offers a variety of techniques for monitoring and characterizing geomorphic processes in high mountain environments. Terrestrial laserscanning and unmanned aerial vehicles provide very high resolution data with high accuracy. Automatic cameras have become a valuable source of information—mostly in a qualitative manner—in recent years. The availability of satellite data with very high revisiting time has gained momentum through the European Space Agency’s Sentinel missions, offering new application potential for Earth observation. This paper reviews the status of recent techniques such as terrestrial laserscanning, remote sensed imagery, and synthetic aperture radar in monitoring high mountain environments with a particular focus on the impact of new platforms such as Sentinel-1 and -2 as well as unmanned aerial vehicles. The study area comprises the high mountain glacial environment at the Pasterze Glacier, Austria. The area is characterized by a highly dynamic geomorphological evolution and by being subject to intensive scientific research as well as long-term monitoring. We primarily evaluate landform classification and process characterization for: (i) the proglacial lake; (ii) icebergs; (iii) the glacier river; (iv) valley-bottom processes; (v) slope processes; and (vi) rock wall processes. We focus on assessing the potential of every single method both in spatial and temporal resolution in characterizing different geomorphic processes. Examples of the individual techniques are evaluated qualitatively and quantitatively in the context of: (i) morphometric analysis; (ii) applicability in high alpine regions; and (iii) comparability of the methods among themselves. The final frame of this article includes considerations on scale dependent process detectability and characterization potentials of these Earth observation methods, along with strengths and limitations in applying these methods in high alpine regions.

An Open-Source Web Platform to Share Multisource, Multisensor Geospatial Data and Measurements of Ground Deformation in Mountain Areas

Nowadays, the increasing demand to collect, manage and share archives of data supporting geo-hydrological processes investigations requires the development of spatial data infrastructure able to store geospatial data and ground deformation measurements, also considering multisource and heterogeneous data. We exploited the GeoNetwork open-source software to simultaneously organize in-situ measurements and radar sensor observations, collected in the framework of the HAMMER project study areas, all located in high mountain regions distributed in the Alpines, Apennines, Pyrenees and Andes mountain chains, mainly focusing on active landslides. Taking advantage of this free and internationally recognized platform based on standard protocols, we present a valuable instrument to manage data and metadata, both in-situ surface measurements, typically acquired at local scale for short periods (e.g., during emergency), and satellite observations, usually exploited for regional scale analysis of surface displacement. Using a dedicated web-interface, all the results derived by instrumental acquisitions and by processing of remote sensing images can be queried, analyzed and downloaded from both expert users and stakeholders. This leads to a useful instrument able to share various information within the scientific community, including the opportunity of reprocessing the raw data for other purposes and in other contexts.

Wide-Area InSAR Survey of Surface Deformation in Urban Areas and Geothermal Fields in the Eastern Trans-Mexican Volcanic Belt, Mexico

This paper provides the first wide-area Interferometric Synthetic Aperture Radar (InSAR) survey of the whole eastern Trans-Mexican Volcanic Belt (42,200 km2). The aims are to identify ground deformation hotspots within major urbanized areas and rural valleys, establish baselines in geothermal exploration sites, and analyze deformation at geothermal exploitation sites and its relationship with energy production. The whole 2003–2010 ENVISAT C-band SAR archive available over the region was processed with the Small BAseline Subset (SBAS) InSAR method to retrieve over 840,000 coherent targets and estimate their ground displacement rates and time series. Land subsidence hotspots due to aquifer drawdown are found within the city of Puebla (up to −53 mm/year vertical rates, groundwater pumping for industrial use), Tlaxcala and Apizaco (−17 mm/year, industrial and public), the valley of Tecamachalco (−22 mm/year, agricultural), Tulancingo (−55 mm/year, public, industrial and agricultural), and in the eastern Mexico City metropolitan area (−44 mm/year, agricultural). The baseline for the Acoculco caldera complex shows widespread ground stability. Conversely, localized subsidence patterns of −5 to −10 mm/year exist around Las Derrumbadas and Cerro Pinto in the Serdán-Oriental basin, due to intense groundwater pumping for agriculture. A well-defined land subsidence area with −11 mm/year maximum rates is found at Los Humeros volcanic complex within Los Potreros collapse, correlating well with energy production infrastructure location and historical steam production rates. Field surveys carried out in Acoculco and Los Humeros in 2018 provide supporting evidence for the identification of hydrothermal manifestations, and understanding of the landscape and surface deformation patterns within the geothermal fields.

COSMO-SkyMed SAR for Detection and Monitoring of Archaeological and Cultural Heritage Sites

Synthetic aperture radar (SAR) imagery has long been used in archaeology since the earliest space radar missions in the 1980s. In the current scenario of SAR missions, the Italian Space Agency (ASI)’s COnstellation of small Satellites for Mediterranean basin Observation (COSMO-SkyMed) has peculiar properties that make this mission of potential use by archaeologists and heritage practitioners: high to very high spatial resolution, site revisit of up to one day, and conspicuous image archives over cultural heritage sites across the globe. While recent literature and the number of research projects using COSMO-SkyMed data for science and applied research suggest a growing interest in these data, it is felt that COSMO-SkyMed still needs to be further disseminated across the archaeological remote sensing community. This paper therefore offers a portfolio of use-cases that were developed in the last two years in the Scientific Research Unit of ASI, where COSMO-SkyMed data were analysed to study and monitor cultural landscapes and heritage sites. SAR-based applications in archaeological and cultural heritage sites in Peru, Syria, Italy, and Iraq, provide evidence on how subsurface and buried features can be detected by interpreting SAR backscatter, its spatial and temporal changes, and interferometric coherence, and how SAR-derived digital elevation models (DEM) can be used to survey surface archaeological features. The use-cases also showcase how high temporal revisit SAR time series can support environmental monitoring of land surface processes, and condition assessment of archaeological heritage and landscape disturbance due to anthropogenic impact (e.g., agriculture, mining, looting). For the first time, this paper provides an overview of the capabilities of COSMO-SkyMed imagery in StripMap Himage and Spotlight-2 mode to support archaeological studies, with the aim to encourage remote sensing scientists and archaeologists to search for and exploit these data for their investigations and research activities. Furthermore, some considerations are made with regard to the perspectives opened by the upcoming launch of ASI’s COSMO-SkyMed Second Generation constellation.

Multiscale Analysis of DInSAR Measurements for Multi-Source Investigation at Uturuncu Volcano (Bolivia)

Uturuncu volcano (southwestern Bolivia) is localized within one of the largest updoming volcanic zones, the Altiplano Puna Volcanic Complex (APVC). In several geodetic studies the observed uplift phenomenon is analyzed and modeled by considering a deep source, related to the Altiplano Puna Magma Body (APMB). In this framework, we perform a multiscale analysis on the 2003–2010 ENVISAT satellite data to investigate the existence of a multi-source scenario for this region. The proposed analysis is based on Cross-correlation and Multiridge method, pointing out the spatial and temporal multiscale properties of the deformation field. In particular, we analyze the vertical component of ground deformation during two time interval: within the 2005–2008 time interval an inflating source at 18.7 km depth beneath the central zone of the APVC is retrieved; this result is in good agreement with those proposed by several authors for the APMB. Between August 2006 and February 2007, we identify a further inflating source at 4.5 km depth, beneath Uturuncu volcano; the existence of this latter, located just below the 2009–2010 seismic swarm, is supported by petrological, geochemical, and geophysical evidence, indicating as a possible interpretative scenario the action of shallow, temporarily trapped fluids.

Multiridge Method for Studying Ground-Deformation Sources: Application to Volcanic Environments

Volcanic phenomena are currently monitored by the detection of physical and chemical observations. Generally, the ground deformation field is the most relevant shallow expression of the geometric and physical parameters variations in the magmatic reservoir. In this study, we propose a novel method for the direct estimation of the geometric parameters of sources responsible for volcanic ground deformation detected via the DInSAR technique. Starting with the biharmonic properties of the deformation field, we define an approach based on the Multiridge and ScalFun methods to achieve relevant information about both the positions and shapes of active sources, such as the Mogi source. Our methodology is definitely different from the methods currently used for modeling ground-deformation sources, mainly based on forward or inverse techniques. In fact, (i) it does not require any assumptions about the source type, and (ii) it is not influenced by the distribution of medium elastic parameters or (iii) the presence of high-frequency noise in the dataset. For synthetic cases, we accurately estimate the depth to the source within a 3% error. Finally, we study the real case of the Okmok volcano ground-deformation field and achieve results compatible with those in previous works.