Absolute Elevation Time Series of Peat Soils from InSAR Observations using Integrated Geodetic Reference Stations

Peat areas in the Netherlands are expected to exhibit extremely dynamic vertical motion, including both reversible and irreversible components. Yet the exact behaviour as a function of time is unknown, and is spatially very variable. This results in a poorly known estimation of greenhouse gas emissions and impact to existing infrastructure, and consequently limited ability to design and deploy mitigating or adaptive measures. In situ measurements are inefficient to capture this dynamic spatio-temporal variability. To monitor wide areas, InSAR has the necessary coverage, resolution, and high temporal sampling, but is very sensitive to noise due to temporal decorrelation of vegetated areas, which in combination with the high temporal dynamics exacerbates the success rate of the phase ambiguity estimation. Moreover, motion measured with InSAR is inherently relative to a reference point which is typically not stable in time. Thus even if InSAR would yield sufficiently coherent results, interpretation is notoriously difficult. To address these problems, an Integrated Geodetic Reference Station (IGRS) was installed in a peat area. It includes SAR reflectors mechanically coupled to a GNSS antenna. Arcs from the IGRS to natural scatterers can therefore be connected to ETRS89. The resulting arcs are thus single-differences of elevation with respect to a reference epoch. In order to reduce noise, and subsequently increase the success rate of ambiguity resolution, we define Elementary Usage Polygons (EUPs), which are expected to move homogeneously. This is achieved by constraining each EUP to have only one subsurface type and one land usage, which allows for the definition of an expected functional model of the vertical movement. Using the expected functional model in combination with the GNSS measurements on the IGRS we can constrain the unwrapping on the arcs between the IGRS and nearby EUPs in order to improve the success rate of the unwrapping. By making use of elevation model data from 2020, the elevation of the reference epoch can be constrained. The single-difference arcs showing elevation change with respect to the reference epoch are transformed into elevations with respect to NAP, allowing for an absolute elevation time series from InSAR observations in the vegetated rapid-moving peat areas of the Netherlands.