A model-backfeed deformation estimation method for revealing 20-year surface dynamics of the Groningen gas field using multi-platform SAR imagery

Abstract

• The model-backfeed deformation estimation method offers more reliable InSAR deformation parameter estimations. • Multiple Hypothesis Testing is employed to determine the best deformation models. • Geolocation error ellipsoids and Monte Carlo methods are used for spatial integration. • Three SAR missions are used to map 20-year Groningen gas field subsidence. The Groningen gas field, the largest natural gas field in Europe, was discovered in 1959 and started its production in 1963. The Earth surface above it experienced subsidence over the past six decades because of gas extraction activities. To accurately reveal this surface movement with satellite SAR data, our study first proposes and demonstrates a model-backfeed (MBF) InSAR deformation estimation method to improve InSAR deformation time series modeling. This method allowed us to include a priori knowledge and to iteratively optimize functional and stochastic models. Using this method and employing a spatio-temporal SAR data integration method based upon Monte Carlo and Multiple Hypothesis Testing methods, we retrieved the 20-year subsidence history of the Groningen gas field by integrating 32 ERS-1/2, 68 Envisat and 82 Radarsat-2 SAR images. The results show that the maximum cumulative surface subsidence in this gas field has been as much as 15.5 cm between 1995 and 2015. In terms of precision and accuracy, our MBF method offered a better result than the standard Multi-Temporal InSAR analysis method: the Ensemble Coherence increased by 10%–19% and Spatio-Temporal Consistency decreased by 2%–20%. In terms of accuracy, our results better concur with the external GNSS reference observations. We further show that the spatio-temporal SAR data integration method has better links with multi-platform SAR data if the uncertainties of the InSAR geolocation and temporal deformation are included. The study demonstrates that the MBF method optimized the estimation of deformation parameters and mitigated unwrapping errors.