Search

Since the launch of Sentinel-1 mission, automated processing systems have been developed for near real-time monitoring of ground deformation signals. Here, we perform a regional analysis of 5 years over 64 volcanic centres located along the East African Rift System (EARS). We show that the correction of atmospheric signals for the arid and low-elevation EARS volcanoes is less important than for other volcanic environments. We find that the amplitude of the cumulative displacements exceeds three times the temporal noise of the time series (3σ) for 16 of the 64 volcanoes, which includes previously reported deformation signals, and two new ones at Paka and Silali volcanoes. From a 5-year times series, uncertainties in rates of deformation are ∼0.1 cm/yr, whereas uncertainties associated with the choice of reference pixel are typically 0.3–0.6 cm/yr. We fit the time series using simple functional forms and classify seven of the volcano time series as ‘linear’, six as ‘sigmoidal’ and three as ‘hybrid’, enabling us to discriminate between steady deformation and short-term pulses of deformation. This study provides a framework for routine volcano monitoring using InSAR on a continental scale. Here, we focus on Sentinel-1 data from the EARS, but the framework could be expanded to include other satellite systems or global coverage.

Several methods allow accurate measurement of terrain surface motions. Global navigation satellite systems (GNSSes) and interferometry with synthetic aperture radar (InSAR) stand out in terms of measurement accuracy among them. In principle, both methods make it possible to evaluate a three-dimensional vector of the motion of points on the terrain surface. In this work, we dealt with the evaluation of motions in the up–down (U–D) and east–west direction (E–W) over underground gas storage (UGS) from InSAR. One crucial step in breaking down PSInSAR line of sight (LOS) measurements to U–D and E–W components is getting time series derived from individual tracks to the same time frame. This is usually performed by interpolation, but we used an innovative approach: we analyzed individual time series using the Lomb–Scargle periodogram (LSP), which is suitable for periodic noisy and irregularly sampled data; we selected the most significant period, created LSP models, and used them instead of the original time series. Then, it was possible to derive time series values for any arbitrary time step. To validate the results, we installed one GNSS receiver in the Tvrdonice UGS test area to perform independent measurements. The results show a good agreement in the evaluation of motions by both methods. The correlation coefficient between horizontal components from both PSInSAR and GNSS was 0.95 in the case of the E–W component, with an RMSE of 1.75 mm; for U–D they were 0.78 and 2.35 mm, respectively. In addition to comparing the motions in the U–D and E–W directions, we also created a comparison by converting GNSS measurements to a line of sight of the Sentinel-1 satellite to evaluate the conformity of InSAR and GNSS measurements. Based on descending track, the correlation coefficient between LOS from both methods is, on average, 0.97, with an RMSE of 2.70 mm.

Remote-sensing-driven urban change detection has been studied in many ways for decades for a wide field of applications, such as understanding socio-economic impacts, identifying new settlements, or analyzing trends of urban sprawl. Such kinds of analyses are usually carried out manually by selecting high-quality samples that binds them to small-scale scenarios, either temporarily limited or with low spatial or temporal resolution. We propose a fully automated method that uses a large amount of available remote sensing observations for a selected period without the need to manually select samples. This enables continuous urban monitoring in a fully automated process. Furthermore, we combine multispectral optical and synthetic aperture radar (SAR) data from two eras as two mission pairs with synthetic labeling to train a neural network for detecting urban changes and activities. As pairs, we consider European Remote Sensing (ERS-1/2) and Landsat 5 Thematic Mapper (TM) for 1991–2011 and Sentinel 1 and 2 for 2017–2021. For every era, we use three different urban sites—Limassol, Rotterdam, and Liège—with at least 500km2 each, and deep observation time series with hundreds and up to over a thousand of samples. These sites were selected to represent different challenges in training a common neural network due to atmospheric effects, different geographies, and observation coverage. We train one model for each of the two eras using synthetic but noisy labels, which are created automatically by combining state-of-the-art methods, without the availability of existing ground truth data. To combine the benefit of both remote sensing types, the network models are ensembles of optical- and SAR-specialized sub-networks. We study the sensitivity of urban and impervious changes and the contribution of optical and SAR data to the overall solution. Our implementation and trained models are available publicly to enable others to utilize fully automated continuous urban monitoring.

The Sentinel-1 satellite system continuously observes European countries at a relatively high revisit frequency of six days per orbital track. Given the Sentinel-1 configuration, most areas in Czechia are observed every 1–2 days by different tracks in a moderate resolution. This is attractive for various types of analyses by various research groups. The starting point for interferometric (InSAR) processing is an original data provided in a Single Look Complex (SLC) level. This work represents advantages of storing data augmented to a specifically corrected level of data, SLC-C. The presented database contains Czech nationwide Sentinel-1 data stored in burst units that have been pre-processed to the state of a consistent well-coregistered dataset of SLC-C. These are resampled SLC data with their phase values reduced by a topographic phase signature, ready for fast interferometric analyses (an interferogram is generated by a complex conjugate between two stored SLC-C files). The data can be used directly into multitemporal interferometry techniques, e.g., Persistent Scatterers (PS) or Small Baseline (SB) techniques applied here. A further development of the nationwide system utilising SLC-C data would lead into a dynamic state where every new pre-processed burst triggers a processing update to detect unexpected changes from InSAR time series and therefore provides a signal for early warning against a potential dangerous displacement, e.g., a landslide, instability of an engineering structure or a formation of a sinkhole. An update of the processing chain would also allow use of cross-polarised Sentinel-1 data, needed for polarimetric analyses. The current system is running at a national supercomputing centre IT4Innovations in interconnection to the Czech Copernicus Collaborative Ground Segment (CESNET), providing fast on-demand InSAR results over Czech territories. A full nationwide PS processing using data over Czechia was performed in 2017, discovering several areas of land deformation. Its downsampled version and basic findings are demonstrated within the article.

Space-borne Synthetic Aperture Radar (SAR) Interferometry (InSAR) is now a key geophysical tool for surface deformation studies. The European Commission’s Sentinel-1 Constellation began acquiring data systematically in late 2014. The data, which are free and open access, have global coverage at moderate resolution with a 6 or 12-day revisit, enabling researchers to investigate large-scale surface deformation systematically through time. However, full exploitation of the potential of Sentinel-1 requires specific processing approaches as well as the efficient use of modern computing and data storage facilities. Here we present Looking Into Continents from Space with Synthetic Aperture Radar (LiCSAR), an operational system built for large-scale interferometric processing of Sentinel-1 data. LiCSAR is designed to automatically produce geocoded wrapped and unwrapped interferograms and coherence estimates, for large regions, at 0.001° resolution (WGS-84 coordinate system). The products are continuously updated at a frequency depending on prioritised regions (monthly, weekly or live update strategy). The products are open and freely accessible and downloadable through an online portal. We describe the algorithms, processing, and storage solutions implemented in LiCSAR, and show several case studies that use LiCSAR products to measure tectonic and volcanic deformation. We aim to accelerate the uptake of InSAR data by researchers as well as non-expert users by mass producing interferograms and derived products.

Underground gas storage facilities are an important element of the natural gas supply system. They compensate for seasonal fluctuations in natural gas consumption. Their expected lifetime is in tens of years. Continuous monitoring of underground gas storage is therefore very important to ensure its longevity. Periodic injection and withdrawal of natural gas can cause, among other things, vertical movements of the terrain surface. Radar interferometry is a commonly used method for tracking changes in the terrain height. It can register even relatively small height changes (mm/year). The primary aim of our research was to verify whether terrain behavior above a relatively deep underground gas storage can be monitored by this method and to assess the possibility of detecting the occurrence of anomalous terrain behavior in an underground gas storage area such as reactivation of faults in the area. The results show a high correlation between periodic injection and withdrawal of natural gas into/from the underground reservoir and periodic changes in terrain height above it (the amplitude of the height changes is in centimeters), which may allow the detection of anomalous phenomena. We documented special behavior of storage structures in the Vienna Basin: the areas adjacent to the underground gas storages show exactly the opposite phase of vertical movements, i.e., while the terrain above the underground reservoirs rises as natural gas is injected, the adjacent areas subside, and vice versa. Based on the analysis of geological conditions, we tend to conclude that this behavior is conditioned by the tectonic fault structure of the studied area.

The Benínar Dam, located in Southeastern Spain, is an earth-fill dam that has experienced filtration issues since its construction in 1985. Despite the installation of various monitoring systems, the data collected are sparse and inadequate for the dam’s lifetime. The present research integrates Multi-Temporal Interferometric Synthetic Aperture Radar (MT-InSAR) and dam modeling to validate the monitoring of this dam, opening the way to enhanced integrated monitoring systems. MT-InSAR was proved to be a reliable and continuous monitoring system for dam deformation, surpassing previously installed systems in terms of precision. MT-InSAR allowed the almost-continuous monitoring of this dam since 1992, combining ERS, Envisat, and Sentinel-1A/B data. Line-of-sight (LOS) velocities of settlement in the crest of the dam evolved from maximums of −6 mm/year (1992–2000), −4 mm/year (2002–2010), and −2 mm/year (2015–2021) with median values of −2.6 and −3.0 mm/year in the first periods (ERS and Envisat) and −1.3 mm/year in the Sentinel 1-A/B period. These results are consistent with the maximum admissible modeled deformation from construction, confirming that settlement was more intense in the dam’s early stages and decreased over time. MT-InSAR was also used to integrate the monitoring of the dam basin, including critical slopes, quarries, and infrastructures, such as roads, tracks, and spillways. This study allows us to conclude that MT-InSAR and dam modeling are important elements for the integrated monitoring systems of embankment dams. This conclusion supports the complete integration of MT-InSAR and 3D modeling into the monitoring systems of embankment dams, as they are a key complement to traditional geotechnical monitoring and can overcome the main limitations of topographical monitoring.

Worldwide, cetaceans are impacted by human activities, and those populations that occur in shallow-nearshore habitats are particularly vulnerable. We present the results of the first long-term study on the responses of a coastal population of endangered Irrawaddy dolphins to widespread habitat changes. We particularly investigated their responses in terms of distribution and abundance. Boat-based, line-transect surveys were conducted during 12 discrete survey periods in 7 survey years spanning a 15-year period (totaling 78 days and 4,630 km of effort) in Balikpapan Bay, East Kalimantan, Indonesia. Irrawaddy dolphins were sighted on 136 occasions. Through DISTANCE analysis, a decrease in population density in the inner Bay area was observed from 0.45 dolphins/km2 in 2000–2001 (CV = 24%) to 0.34 and 0.32 dolphins/km2 in 2008 and 2015 (CV = 31% and 25%). A shift in distribution was noted between the periods 2000–2002 and 2008–2015 with significantly lower occurrence in the lower Bay segment compared to upper Bay segments. No sightings were made in the outer Bay area in later years, which coincided with increased shipping traffic in these areas. A peak in stranding events in 2016 and 2018 followed extremely high phenol levels within Bay waters in 2015 and a large-scale oil spill in 2018. The mean annual mortality rates of 0.67 Irrawaddy dolphins/year is unsustainable based on the lower potential biological removal (PBR) values for best abundance estimates of 2015 (Ndistance = 45 and Nmark–recapture = 73). Other threats to local dolphins include unsustainable fishing, underwater noise caused by construction, particularly piling activities. The research helped to identify Balikpapan Bay as an Important Marine Mammal Area by the IUCN MMPA Taskforce. Serious concerns remain for the concrete plans to move Indonesia’s capital city to the area north of the Bay, in terms of increased shipping traffic and harbor construction in the upper Bay segments that represent primary dolphin habitat. We recommend that protected areas be assigned for marine mammals and artisanal fisheries and shipping traffic and piling activities be excluded from these areas. We also recommend a legislated requirement of a mitigation protocol compulsory for piling and seismic activities within Indonesia.

For the past five years, the 2-satellite Sentinel-1 constellation has provided abundant and useful Synthetic Aperture Radar (SAR) data, which have the potential to reveal global ground surface deformation at high spatial and temporal resolutions. However, for most users, fully exploiting the large amount of associated data is challenging, especially over wide areas. To help address this challenge, we have developed LiCSBAS, an open-source SAR interferometry (InSAR) time series analysis package that integrates with the automated Sentinel-1 InSAR processor (LiCSAR). LiCSBAS utilizes freely available LiCSAR products, and users can save processing time and disk space while obtaining the results of InSAR time series analysis. In the LiCSBAS processing scheme, interferograms with many unwrapping errors are automatically identified by loop closure and removed. Reliable time series and velocities are derived with the aid of masking using several noise indices. The easy implementation of atmospheric corrections to reduce noise is achieved with the Generic Atmospheric Correction Online Service for InSAR (GACOS). Using case studies in southern Tohoku and the Echigo Plain, Japan, we demonstrate that LiCSBAS applied to LiCSAR products can detect both large-scale (>100 km) and localized (~km) relative displacements with an accuracy of