Characterization and removal of shadow bias from Optical Image Correlation; application to the 2013 Balouchistan earthquake
Abstract
In recent years, optical image correlation (OIC) has become a powerful geodetic tool (in compliment to InSAR) for the retrieval of near-field earthquake ground deformations using optical satellite data. However, because optical satellite images are acquired throughout the year, the changing illumination conditions at the time of each acquisition may be expected to subtly influence the shadow content of each pixel, which in turn may bias the correlation process, resulting in contaminated displacement maps. We use the Landsat8 archive of satellite images acquired several years after the 2013 Balouchistan earthquake (Mw 7.8; left-lateral strike-slip with thrust) to investigate the influence of seasonal variations in sun illumination on the resulting deformation field retrieved with OIC. We find the displacement field varies seasonally as a function of the difference in pre- and post-image sun elevation. The amplitude of the seasonal artifacts correlate strongly with topography (slope angle), and are greatest in the north-south displacement component, which is subject to the largest variation in shadow movement. By correlating different combinations of images, we develop functions for each pixel in the correlation describing how the seasonal displacement varies as a function of the pre- and post-image sun elevation. We then remove seasonal artifacts from various correlations derived from Landsat8 combinations spanning the co- and post-seismic periods. Finally, we invert our illumination-corrected correlations using the NSBAS processing chain, which takes advantage of the large data redundancy to produce a more robust deformation time-series. These two approaches allow us to significantly increase the precision of OIC, such that we can retrieve small (<1 m) post-seismic deformations for the Balouchistan earthquake using medium-resolution (15 m) Landsat8 imagery. We see evidence for extensional afterslip in the upper crust in the week following the earthquake (possibly related to overshoot during the mainshock). In the following months, the signal changes to compression across the fault, as would be expected from the mainshock source mechanism.
- Publication:
-
AGU Fall Meeting Abstracts
- Pub Date:
- December 2017
- Bibcode:
- 2017AGUFM.G43A0909H
- Keywords:
-
- 1207 Transient deformation;
- GEODESY AND GRAVITY;
- 1211 Non-tectonic deformation;
- GEODESY AND GRAVITY;
- 8159 Rheology: crust and lithosphere;
- TECTONOPHYSICS;
- 8163 Rheology and friction of fault zones;
- TECTONOPHYSICS