The repeat period of SAR data and its side-looking characteristics make InSAR time series analysis useful for water level monitoring applications. The standard approach determines corresponding scatterers by focusing the study area on the multipath radar reflections that include the water level. This paper introduces an alternative approach to identifying such signals using two metrics: cosine similarity and temporal differential coherence. The results show that temporal differential coherence can detect phase variations similar to water level by constantly returning high values even when there is an offset, while cosine similarity yields low scores. Within an urban environment, this approach finds point scatterers corresponding to water level changes in or near water, such as permanent floating objects, bridges, and buildings adjacent to water, where the highest differential coherence value was acquired from a permanent floating restaurant in open water.

@en

The growing availability of SAR data offers a real-time deformation monitoring opportunity, but data utilization can be inefficient. Our study introduces a mathematical framework using recursive least-squares and the wrapped phase, allowing efficient updates when new data arrives. This method also incorporates prior knowledge about signal smoothness for non-linear displacement estimation. Compared to the batch solution, our recursive approach achieves parameter estimation without storing past measurements while respecting signal smoothness constraints.

@en

We introduce SARXarray and STMtools, two python libraries designed to support the exploitation of modern Interferometric Synthetic Aperture Radar (InSAR) data by enabling handling of larger-than memory datasets and the incorporation and fusion with relevant contextual information.The libraries are developed upon two innovative and well-established open-source Python libraries:Xarray [5] and Dask [6]. They are implemented as Xarray extensions. SARXarray is designed to manipulate and operate on larger-than-memory coregistered raster stacks such as Single-look Complex images (SLC) or interferograms, performing scatterer selection and producing STM objects. STMtools leverage the Space-Time Matrix (STM) concept [3,4] and provides functionalities to process STM and perform enrichment/data fusion with other data sources.Both libraries are built on the Xarray library, providing support for a wide range of data formats, and utilize Dask for parallel computation, making them scalable for distributed computation infrastructures. By enabling InSAR data analysis incorporating contextual information, the two libraries enhance the potential to uncover underlying mechanisms driving deformation phenomena.

@en

Ghana's coastline has been facing erosion and sedimentation phenomena for several decades, resulting in a serious threat to life and property considering that major urban settlements are located on the coast. In this region, there has been a lack of emphasis on comprehensive, large-scale investigations into coastal changes: prior research has predominantly centered on site-specific assessments. These studies have revealed alarming erosion rates, with reports indicating that nearly ten meters are lost annually. The use of high-resolution remotely sensed data can be a consistent support in regions where physical or economic obstacles interfere with collecting in situ information. In particular, the use of continuous all-weather SAR data may facilitate the evaluation of erosion and sedimentation phenomena in coastal areas. In this paper, we apply SAR data over a time period between 2017 and 2021. Sentinel-1 data are pre-processed using the Google Earth Engine platform, and a dedicated algorithm is then applied to identify and quantify erosion and sedimentation processes. Optical images are used as a reference for detecting the location of two areas where consistent sedimentation and erosion phenomena occurred in the considered four years. The results demonstrate that SAR backscattering variations over time offer a reliable method for monitoring coastal changes. This approach enables the identification of the type of phenomena occurring - sedimentation or erosion -, and allows for the quantification of their intensity and dimensions over time. The method can be worldwide applied once the appropriate thresholds are evaluated and help in predictive studies and environmental planning.

@en

InSAR enables the estimation of spatio-temporal displacements, relative to a reference point and a reference epoch, here defined as the mother image. When dealing with time series, there are several options to treat the mother image in computing and plotting the temporal phase differences, producing distinctly different results, in terms of the estimated displacement parameters and their precision. Here we review the three approaches mostly encountered in literature, discuss the implications of the different approaches, and recommend the 'embracing mother' approach for standard InSAR analyses and visualizations.

@en

We present the preliminary results of an InSAR analysis of peatland surface motion covering a large spatial and temporal extent. This work is the first large scale analysis of the Dutch Green Heart region, and is made possible using a novel distributed scatter (DS) InSAR processing method. This method is designed to handle breakages in the observed interferometric phase time series which occur due to temporal decorrelation, which we designate with the term loss-of-lock.

@en

La plupart des méthodes exploitant les séries temporelles InSAR (interférométrie radar) présentent des caractéristiques similaires qui permettent de les classer en deux catégories de techniques. La première catégorie de techniques est basée sur des réflecteurs distribués (DS) pour la surveillance des déformations. Une façon courante de réduire la décorrélation des signaux consiste à sélectionner des interférogrammes avec des lignes de base spatiales et temporelles courtes (technique des small baselines subset, SBAS). La seconde catégorie est celle de l’approche PS (permanent/persistant scatterers), qui utilise des réflecteurs individuels dominant le signal à l’intérieur d’une cellule de résolution pour suivre la déformation dans le temps. Une technique avancée récente nous permet de combiner à la fois les PS et les DS pour surmonter la rareté des points identifiés pour l’estimation (technique PSDS combinant réflecteurs persistants et distribués). Ce chapitre décrit les deux principales familles de techniques InSAR de séries temporelles (SBAS et PSDS). @en

Due to its sensitivity to water vapor, high resolution, and global availability, interferometric satellite radar (InSAR) has a large but unexploited potential for the improvement of regional NWP models. A relatively straightforward approach is to exploit the exact instantaneous character of the InSAR data in data assimilation to improve the timing of NWP model realizations. Here we show the potential impact of InSAR data on the NWP model timing and subsequently on improved model performance. By time-shifting the model to find the best match with the InSAR data we show that we can achieve a model error reduction (one-sigma) of up to 40% in cases where weather fronts are present, while other cases show more modest improvements. Most model performance gain due to time-shifts can therefore be achieved in cases where weather fronts are present over the study area. The model-timing errors related to the maximum model error reduction for these cases are in the order of (Formula presented.) 30 min.

@en

Deformation phenomena on Earth are inherently three dimensional. With SAR interferometry (InSAR), in many practical situations the maximum number of observations is two (ascending and descending), resulting in an infinite number of possible displacement estimates. Here we propose a practical solution to this underdeterminancy problem in the form of the strapdown approach. With the strapdown approach, it is possible to obtain “3D-global/2D-local” solutions, by using minimal and largely undisputed contextual information, on the expected driving mechanisms and/or spatial geometry. It is a generic method that defines a local reference system with transversal, longitudinal, and normal (TLN) axes, with displacement occurring in the transversal-normal plane only. Since the orientation of the local frame is based on the physics of the problem at hand, the strapdown approach gives physically more relevant estimates compared to conventional approaches. Moreover, using an a-priori uncertainty approximation on the orientation of the local frame it is possible to assess the precision of the final estimates. As a result, appropriate cartographic visualization using a vector map with confidence ellipses enables an improved interpretation of the results.@en

The estimation of displacement vectors for (objects on) the Earth's surface using satellite InSAR requires geometric transformations of the observables based on orbital viewing geometries. Usually, there are insufficient viewing geometries available for full 3-D reconstruction, leading to nonunique solutions. Currently, there is no standardized approach to deal with this problem, resulting in products that are based on haphazard and/or oversimplified assumptions with biased estimates and reduced interpretability. Here, we show that a clear definition of - and subsequent adherence to - enabling conditions guarantees the validity and quality of displacement vector estimates leading to standardized interferometric products with improved interpretability. We introduce the concept of the null line as a key metric for InSAR geometry and bias estimation, assess its impact and orientation for all positions on Earth, and propose a novel reference system that is inherently unbiased. We evaluate current operational practice, leading to a taxonomy of frequently encountered misconceptions and to recommendations for InSAR product generation and interpretation. We also propose new subscript notation to uniquely distinguish different projection and decomposition products. Our propositions contribute to further standardization of InSAR product definition, improved map annotation, and robust interpretability.

@en

The observed phase in time series of interferometric synthetic aperture radar (InSAR) products is a superposition of various components. Differential topography, line-of-sight displacements, and differential atmospheric delays are the main contributions and need to be disentangled to derive accurate digital elevation model (DEM), deformation, or atmospherical products from InSAR. However, isolating the atmospheric component has been proven difficult as it is spatiotemporally highly dynamic and a superposition of two atmospheric states. Here, we propose an approach to parameterize the stochastic properties of the single-epoch atmospheric delay field as a way to define the atmospheric signal. We found that the atmospheric signal of a time series of interferograms can be characterized by structure functions, which can be used to isolate the single-epoch structure functions. Due to the scaling properties of the atmospheric signal, it is then possible to construct a parametric function per SAR acquisition, using two isotropic and three anisotropic parameters. In particular, the isotropic parameters for the short-distance variation and long-distance variation in atmospheric delay can be used to characterize the atmospheric signal. For a test set of 151 Sentinel-1 acquisitions, this results in an atmospheric energy range of about 10 for short-distance scales and about 50 for long-distance scales. Our parameterization demonstrates that we can describe the spatiotemporal variability of InSAR atmospheric delays, which provides a measure for atmospheric noise for individual epochs in deformation time series based on distance and azimuth.

@en

We present SPAMS: Simple Parameterization for the Motion of Soils, a model to describe the motion of deformable soils in the Vadose zone, mainly peat and clay, herein called shallow soft soils. The SPAMS model estimates the reversible and irreversible vertical component of surface displacement to within sub-centimetre RMSE, using only four parameters: three scaling factors and an integration time. Requiring only meteorological data as an input, its lightweight nature and simple implementation make it a powerful tool when used as a first approximation in inverse problems like those encountered in remote sensing. It has been validated against in-situ data from five test sites in The Netherlands with different Holocene soil strata.

@en

In the view of climate change, understanding and managing effects on coastal areas and adjacent cities is essential. Permanent Laser Scanning (PLS) is a successful technique to not only observe notably sandy coasts incidentally or once every year, but (nearly) continuously over extended periods of time. The collected point cloud observations form a 4D point cloud data set representing the evolution of the coast provide the opportunity to assess change processes at high level of detail. For an exemplary location in Noordwijk, The Netherlands, three years of hourly point clouds were acquired on a 1 km long section of a typical Dutch urban sandy beach. Often, the so-called level of detection is used to assess point cloud differences from two epochs. To explicitly incorporate the temporal dimension of the height estimates from the point cloud data set, we revisit statistical testing theory. We apply multiple hypothesis testing on elevation time series in order to identify different coastal processes, like aeolian sand transport or bulldozer works. We then estimate the minimal detectable bias for different alternative hypotheses, to quantify the minimal elevation change that can be estimated from the PLS observations over a certain period of time. Additionally, we analyse potential error sources and influences on the elevation estimations and provide orders of magnitudes and possible ways to deal with them. Finally we conclude that elevation time series from a long term PLS data set are a suitable input to identify aeolian sand transport with the help of multiple hypothesis testing. In our example case, slopes of 0.032 m/day and sudden changes of 0.031 m can be identified with statistical power of 80% and with 95% significance in 24-h time series on the upper beach. In the intertidal area the presented method allows to classify daily elevation time series over one month according to the dominating model (sudden change or linear trend) in either eroding or accreting behaviour.

@en

We introduce the term loss-of-lock to describe a specific form of coherence loss that results in the breakage of a synthetic aperture radar interferometric (InSAR) time series. Loss-of-lock creates a specific pattern in the coherence matrix of a multilooked distributed scatterer (DS) by which it may be detected. Along with identification, we introduce a new DS processing methodology that is designed to mitigate the effects of loss-of-lock by introducing contextual data to assist in the time-series processing. This methodology is of particular relevance to regions that suffer from severe temporal decorrelation, such as northern peatlands. We apply our new method to two subsiding cultivated peatland regions in The Netherlands which previously proved impossible to monitor using DS InSAR techniques. Our results show a very good agreement with in situ validation data as well as spatial correlation between regions and the natural terrain.

@en

Multi-epoch interferometric synthetic aperture radar (InSAR) is a highly effective technique for monitoring deformation in urban areas. However, interpreting InSAR deformation can be challenging due to various factors, including inherent geometric imaging distortion, the intricate structure and deformation properties of targets in urban scenes, and the multiple scattering of microwave signals between objects in urban scenes. This paper discusses the challenges involved in interpreting time-series InSAR deformation: (1) Precisely identifying the location of deformation signals and linking them to their corresponding objects, i.e., determining where the deformation signal occurs, (2) understanding the mechanisms and factors that cause the detected deformation signals, i.e., determining what the deformation signal represents, (3) establishing the connection among the detected deformation signals, the deformation events, and the scattering mechanisms. We suggest a parametric framework to improve the accurate interpretation of InSAR deformation. This framework includes several factors, including kinematic characteristics (deformation rate, cumulative deformation, deformation gradient, and deformation model), geometric parameters (position, size, structure, orientation, and roughness), semantic information (land cover type, terrain morphology, texture, and auxiliary information on natural and anthropogenic disturbance) and physical properties (scattering mechanism, penetrability, extensibility, conductivity, and thermal conductivity). Our approach aims to enhance the representation of coherent points for a better understanding of InSAR deformation. This paper offers a comprehensive overview of the advancements achieved in extracting parameters of InSAR coherent points and interpreting deformation based on geometric parameters, semantic information, and physical properties. High-precision 3D positioning is crucial for InSAR fine monitoring in urban areas. It helps determine the source of deformation signals and facilitates the analysis of deformation mechanisms. Semantic information, such as 3D models, high-resolution optical images, laser point cloud data, and land use data, can aid in interpreting InSAR deformation. By combining InSAR deformation data with a deep learning approach, there is an opportunity to interpret deformations effectively. In urban environments, the scattering mechanism of ground objects is complex. Multiple scattering signals can provide effective observations of deformation and information about the target’s size. However, combining the scattering mechanism of synthetic aperture radar signals to carry out parameter inversion and deformation mechanism interpretation of urban target terrain remains a challenge. The framework, which considers the geometric parameters, semantic information, and physical attributes of InSAR coherent points, will be crucial for deformation interpretation and mechanism cognition. This framework will enable fine deformation monitoring, intelligent recognition, and application in future urban scenes.

@en

InSAR enables the estimation of displacements of (objects on) the earth's surface. To provide reliable estimates, both a stochastic and mathematical model are required. However, the intrinsic problem of InSAR is that both are unknown. Here we derive the Variance-Covariance Matrix (VCM) for double differenced phase observations for an arc, i.e., the phase difference between two points relative to a reference epoch. Using the Normalized Amplitude Dispersion we subdivide the time series in multiple partitions. The method results in a more realistic stochastic model, and consequently more realistic and reliable displacement parameters. The stochastic model also allows to make statements on the precision and reliability of the estimated parameters.

@en

In the planning stage of new infrastructure or when designing renovation of existing infrastructure, information about existing slope movements or settlements is essential to make informed design decisions. Interferometric
Synthetic Aperture Radar (InSAR) techniques can be of value to identify these risks in an early stage of a project. InSAR can offer insight into the surface movements of an area from historic archives using satellite-based SAR
data. Furthermore, InSAR observations can help identify zones with displacements larger than the average of an area, and be used to plan future soil investigation more effectively. Thanks to their high temporal and spatial resolution, InSAR observations can also complement in situ conventional monitoring during the construction and operational stage. Despite these possibilities, the use of InSAR is not yet standard practice in geotechnical projects
and no formal guidelines are currently available to inform engineers, planners and infrastructure stakeholder on the use of InSAR-based monitoring within geotechnical design codes. Here we provide an operational framework for the practical integration of InSAR monitoring into current geotechnical design codes, such as Eurocode-7, for all project stages. The proposed framework is then demonstrated for the planning stage of a highway renovation project, focusing on an area potentially subjected to landslides where no conventional monitoring
data was available at this stage. We concluded that the proposed framework is a practical and operational
tool that can be used by planners and engineers in the whole lifecycle of an infrastructure project.@en

Landslides are a major geohazard in hilly and mountainous environments. In-situ inspection of downslope motion is costly, sometimes dangerous and, requires prior knowledge of the existence of a landslide. Remote sensing from space is a way to detect and characterize landslides systematically at large scale. Interferometric Synthetic Aperture Radar (InSAR) has shown to be a valuable resource of deformation information, but it requires expert knowledge and considerable computational efforts. Moreover, the successful application of InSAR for landslides requires a favorable acquisition geometry relative to the landslide deformation pattern. Consequently, there is a need for a widely applicable tool to assess the potential of InSAR at a particular location a priori. Here we present a novel, generic approach to assess the potential of InSAR-based deformation tracking, providing a standardised and automated method applicable on any slope. We define the detection potential as the sensitivity of InSAR to detect downslope displacement combined with the presence of coherently scattering surfaces. We show that deformation can be detected on at least 91% of the global landslide-prone slopes, and provide an open source Google Earth Engine tool for the quick assessment of the availability of potential coherent scatterers. This tool enables any person interested in applying InSAR to routinely assess the potential for monitoring landslide deformation in their region of interest.@en

Large parts of the Netherlands have experienced subsidence due to anthropogenic soil deformation (peat oxidation and clay shrinkage and swell) since the Middle Ages. As of the start of gas production from the Groningen reservoir in 1968, the northeastern Netherlands has seen additional subsidence (>30 cm) originating from deep reservoir compaction. Due the land being situated close to sea-level this poses a significant societal problem. Geodetic surface observations (InSAR) contain components of both shallow and deep subsidence. This study is part of a larger project on subsidence forecasting, which aims to disentangle deep and shallow (soil) drivers of subsidence by assimilating geodetic time series in geomechanical models of the subsurface. This requires the models to be highly efficient. In this study, we focus on modelling the deep subsurface (reservoir and overburden). We perform a sensitivity analysis to investigate the level of geometrical and mechanical complexity of the reservoir and overburden that can be resolved from the data, to reach optimal model efficiency. We employ a semi-analytical mechanical model for the Groningen subsurface using the PSGRN/PSCMP code by Wang et al. (2006) of laterally uniform viscoelastic layers. We use the Visvalingam–Whyatt algorithm to simplify the existing geological model to versions with decreasing levels of complexity. The surface deformation resulting from the different model versions are used in a convergence test to define the required model complexity. We find that we can achieve orders of magnitude improvement in model run time efficiency, depending on the data uncertainty. Wang, R., Lorenzo Martín, F., Roth, F. (2006): PSGRN/PSCMP - a new code for calculating co- and post-seismic deformation, geoid and gravity changes based on the viscoelastic-gravitational dislocation theory. Computers and Geosciences, 32, 4, 527-541.@en

Synthetic aperture radar (SAR) missions with short repeat times enable opportunities for near real-time deformation monitoring. Traditional multitemporal interferometric SAR (MT-InSAR) is able to monitor long-term and periodic deformation with high precision by time-series analysis. However, as time series lengthen, it is time-consuming to update the current results by reprocessing the whole dataset. Additionally, the number of coherent scatterers varies over time due to disappearing and emerging scatterers due to inevitable changes in surface scattering, and potential deformation anomalies require changes in the prevailing deformation model. Here, we propose a novel method to analyze InSAR time series recursively and detect both significant changes in scattering as well as deformation anomalies based on the new acquisitions. Sequential change detection is developed to identify temporary coherent scatterers (TCSs) using amplitude time series. Based on the predicted phase residuals, scatterers with abnormal deformation displacements are identified by a generalized ratio test, while the parameters of stable scatterers are updated using Kalman filtering. The quality of the anomaly detection is assessed based on the detectability power and the minimum detectable deformation. This facilitates (near) real-time data processing and decreases the false alarm likelihood. Experimental results show that the technique can be used for the real-time evaluation of deformation risks.

@en