410 likes | 645 Views
Institute of Geophysics and Tectonics. Robust corrections for topographically-correlated atmospheric noise in InSAR data from large deforming regions. By David Bekaert Andy Hooper, Tim Wright and Richard Walters. Why a tropospheric correction for InSAR?. Tectonic.
E N D
Institute of Geophysics and Tectonics Robust corrections for topographically-correlated atmospheric noise in InSAR data from large deforming regions By David Bekaert Andy Hooper, Tim Wright and Richard Walters
Why a tropospheric correction for InSAR? Tectonic • To extract smaller deformation signals 100 km cm -10 13.5 Over 9 months
Why a tropospheric correction for InSAR? Troposphere Tectonic • To extract smaller deformation signals • Tropospheric delays can reach up to 15 cm • With the tropospheric delay a superposition of • Short wavelength turbulent component • Topography correlated component • Long wavelength component 100 km cm -10 13.5 Over 9 months 1 interferogram (ti –tj)
Tropospheric corrections for an interferogram • Auxiliary information (e.g.): Limitations • GPS • Weather models • Spectrometer data • Station distribution • Accuracy and resolution • Cloud cover and temporal sampling
Tropospheric corrections for an interferogram • Auxiliary information (e.g.): Limitations • GPS • Weather models • Spectrometer data Interferometric phase • Linear estimation (non-deforming region or band filtering) • Station distribution • Accuracy and resolution • Cloud cover and temporal sampling Assumes a laterally uniform troposphere isolines
A laterally uniform troposphere Interferogram Linear est • A linear correction can work in small regions Tropo GPS InSAR and GPS data property of IGN isolines
A spatially varying troposphere Topography est: Spectrometer & Linear • However • Spatial variation of troposphere A linear correction can work in small regions + + - + isolines
A spatially varying troposphere Allowing for spatial variation -9.75 rad 9.97 Interferogram (Δɸ) Why not estimate a linear function locally?
A spatially varying troposphere -9.75 rad 9.97 Interferogram (Δɸ) Why not estimate a linear function locally? Does not work as: Const is also spatially-varying and cannot be estimated from original phase!
A spatially varying troposphere -9.75 rad 9.97 Interferogram (Δɸ) Why not estimate a linear function locally? Does not work as: Const is also spatially-varying and cannot be estimated from original phase! We propose a power-law relationship that can be estimated locally
Allowing for spatial variation With h0 the lowest height at which the relative tropospheric delays ~0 • 7-14 km from balloon sounding Sounding data provided by the University of Wyoming
Allowing for spatial variation Allowing for spatial variation With h0 the lowest height at which the relative tropospheric delays ~0 • 7-14 km from balloon sounding With α a power-law describing the decay of the tropospheric delay • 1.3-2 from balloon sounding data Sounding data provided by the University of Wyoming
Power-law example -9.75 rad 9.97 Interferogram (Δɸ)
Power-law example -9.75 rad 9.97 (Y. Lin et al., 2010, G3) for a linear approach Band filtered: phase (Δɸband) & topography (h0-h)αband Interferogram (Δɸ)
Power-law example (Y. Lin et al., 2010, G3) for a linear approach Band filtered: phase (Δɸband) & topography (h0-h)αband
Power-law example (Y. Lin et al., 2010, G3) for a linear approach Band filtered: phase (Δɸband) & topography (h0-h)αband Anti-correlated! For each window: estimate Kspatial
Power-law example (Y. Lin et al., 2010, G3) for a linear approach Band filtered: phase (Δɸband) & topography (h0-h)αband Anti-correlated! For each window: estimate Kspatial
Power-law example rad/mα -1.1e-69.8e-5 Original phase (Δɸ) Tropo variability (Kspatial) Band filtered: phase (Δɸband) & topography (h0-h)αband
Power-law example 1/mα rad/mα 4.7e4 2.4e5 -1.1e-69.8e-5 -9.75 rad 9.97 Power-law est (Δɸtropo) Topography (h0-h)α Original phase (Δɸ) Tropo variability (Kspatial) Band filtered: phase (Δɸband) & topography (h-h0)αband
Power-law example Allowing for spatial variation -9.75 rad 9.97 -9.75 rad 9.97 -9.75 rad 9.97 Spectrometer est (Δɸtropo) Original phase (Δɸ) Power-law est (Δɸtropo)
Case study regions Regions: • El Hierro (Canary Island) • GPS • Weather model • Uniform correction • Non-uniform correction • Guerrero (Mexico) • MERIS spectrometer • Weather model • Uniform correction • Non-uniform correction
El Hierro Interferograms (original) -11.2 rad 10.7
El Hierro Interferograms (original) -11.2 rad 10.7 WRF (weather model)
El Hierro Interferograms (original) -11.2 rad 10.7 WRF (weather model)
El Hierro Interferograms (original) -11.2 rad 10.7 WRF (weather model) Linear (uniform)
El Hierro Interferograms (original) -11.2 rad 10.7 WRF (weather model) Linear (uniform) Power-law (spatial var)
El Hierro quantification • ERA-I run at 75 km resolution • WRF run at 3 km resolution
Mexico -9.75 rad 9.97 MERIS MERIS
Mexico -9.75 rad 9.97 MERIS MERIS Clouds
Mexico -9.75 rad 9.97 (Weather model) MERIS ERA-I MERIS ERA-I
Mexico -9.75 rad 9.97 (Weather model) MERIS ERA-I MERIS ERA-I Misfit near coast
Mexico -9.75 rad 9.97 (Weather model) Linear Linear MERIS ERA-I MERIS ERA-I
Mexico -9.75 rad 9.97 (Weather model) Linear Linear MERIS ERA-I MERIS ERA-I
Mexico -9.75 rad 9.97 (Weather model) Linear Power-law Power-law Linear MERIS ERA-I MERIS ERA-I
Mexico -9.75 rad 9.97 (Weather model) Linear Power-law Power-law Linear MERIS ERA-I MERIS ERA-I
Mexico techniques compared: profile AA’ MERIS ERA-I A’ A Linear Power-law
Mexico quantification MERIS accuracy (Z. Li et al., 2006)
Summary/Conclusions • Fixing a reference at the ‘relative’ top of the troposphere allows us to deal with spatially-varying tropospheric delays. • Band filtering can be used to separate tectonic and tropospheric components of the delay in a single interferogram • A simple power-law relationship does a reasonable job of modelling the topographically-correlated part of the tropospheric delay. • Results compare well with weather models, GPS and spectrometer correction methods. • Unlike a linear correction, it is capable of capturing long-wavelength spatial variation of the troposphere. Toolbox with presented techniques will be made available to the community