Drained and cultivated grasslands on peat soils behave as a significant source of greenhouse gasses by oxidation. However, the lack of empirical estimates of carbon losses from peatlands with adequate spatial and temporal resolution has forced researchers to rely on process-based model approximations to make quantitative, regional- or national-scale estimations. Here we use satellite-based synthetic aperture radar interferometry to estimate the land motion per parcel with a daily resolution, discriminate a reversible and an irreversible component, and convert this to an upper bound of (Formula presented.) -equivalent emissions over the western part of the Netherlands. We find an upper bound of 21.5 (Formula presented.) -eq/ha/yr, corresponding to a total regional output of 2.3 (Formula presented.) -eq/yr, or approximately 1.3% of the entire greenhouse gas emissions of the Netherlands in 2019. The method also allows us to provide estimates for future emissions as well as evaluate the efficacy of installed subsidence mitigation measures.

Peat subsidence occurs when parts of the peat soil interact with air, usually due to water table lowering, then triggers peat consolidation, shrinkage, and oxidation, releasing substantial CO₂ emissions. Managing and mitigating these impacts requires a comprehensive understanding of the mechanisms and the spatio-temporal variations of the subsidence. Advanced space geodetic techniques, particularly InSAR, enable surface displacement monitoring. While time series InSAR analysis effectively estimates displacement, its precision, accuracy, and representativity are compromised by temporal decorrelation, noise, and dynamic soil movement, especially over pastures on peat soils. Moreover, loss-of-lock events caused by an irrecoverable loss of coherence disrupt the time series and introduce arbitrary unintelligible phase offsets. Strategies such as multilooking using contextual information have improved the reliability of the InSAR displacement estimates. However, more experience in the efficacy of InSAR-based surface dynamics assessments is required. This study estimates and analyzes surface motion in a regional peat area in Midden-Delfland, The Netherlands, using Sentinel-1 data and the SPAMS model. SPAMS incorporates precipitation and evapotranspiration information to estimate surface motion parameters, distinguishing between reversible and irreversible subsidence. The results reveal an average subsidence rate of −5.4±0.7 mm/year within the study area. Irreversible subsidence is strongly correlated with climatic conditions, with the most significant subsidence observed during a prolonged dry period in the summers of 2018 and 2022. Mitigating peatland subsidence includes preserving soil water content, especially during dry periods. Integrating InSAR and SPAMS provides a valuable tool for monitoring peat surface elevation, water management, and reducing peatland degradation.

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.

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.

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.

Phase unwrapping, also known as ambiguity resolution, is an underdetermined problem in which assumptions must be made to obtain a result in SAR interferometry (InSAR) time series analysis. This problem is particularly acute for distributed scatterer InSAR, in which noise levels can be so large that they are comparable in magnitude to the signal of investigation. In addition, deformation rates can be highly nonlinear and orders of magnitude larger than neighboring point scatterers, which may be part of a more stable object. The combination of these factors has often proven too challenging for the conventional InSAR processing methods to successfully monitor these regions. We present a methodology which allows for additional environmental information to be integrated into the phase unwrapping procedure, thereby alleviating the problems described above. We show how problematic epochs that cause errors in the temporal phase unwrapping process can be anticipated by the machine learning algorithms which can create categorical predictions about the relative ambiguity level based on the readily available meteorological data. These predictions significantly assist in the interpretation of large changes in the wrapped interferometric phase and enable the monitoring of environments not previously possible using standard minimum gradient phase unwrapping techniques.