Multiband 2D InSAR deformation models with error estimates from natural neighbour interpolation: Case study in the Latrobe Valley, Australia

Abstract

A long term objective for Australia is to develop a continent-wide 2D deformation model based on satellite interferometric synthetic aperture radar (InSAR). Nation-wide InSAR processing has been demonstrated in European countries and Japan. However, there are some key differences to the Australian situation for a 2D model (i.e., vertical and east–west surface movement rates) where SAR data from ascending and descending satellite orbits are needed. The primary difference is that while ESA’s Sentinel-1 ascending and descending satellite observations underpin 2D models in other parts of the world, the Australian continent is only covered by descending Sentinel-1 SAR data in most areas. Combined with Australia’s large size of 7.63 million km2, this means a different approach is needed to realise 2D deformation models. In this article, we propose the novel approach of (1) combining supplementary multiband SAR (ascending X-band and L-band) with the descending C-band Sentinel-1 data, and (2) using natural neighbour (NN) interpolation as a computationally efficient interpolation method for combining the multiband InSAR data on common grids prior to 2D model computation. This is an alternative to the computationally intensive Ordinary Kriging (OK) approach. A disadvantage of NN compared to OK is that it does not compute an error estimate, which we address by deriving an error estimate for the NN interpolation so that errors are provided for each 2D rate to aid with their correct interpretation. To demonstrate our new approach, we compute experimental 2D deformation models that propagate the full error budget from InSAR, a GNSS connection to ITRF2014 and NN interpolation into an error grid in the Latrobe Valley in southeast Australia. The significance of our approach is that it provides a suitable combination of methods that may be used to develop 2D continental deformation models for Australia’s unique circumstances. Our results demonstrate the accuracy of NN and OK interpolation are at the same level (within error), and the error estimates from both interpolation techniques are also comparable: X band (TerraSAR-X) ∼ 1 mm/yr; C band (Sentinel-1) ∼ 2.5 mm/yr, and L band (ALOS PALSAR1) 4 mm/yr. For the first time, we reveal 2D deformation in the Latrobe Valley region showing horizontal east–west deformation in the area surrounding an above-ground mine and reaching 20 mm/yr, with vertical subsidence rates reaching −40 mm/yr.