Surface Deformation Measurements Scientific Requirements & Challenges

Transcription

1 Surface Deformation Measurements Scientific Requirements & Challenges 1st Science and Application Workshop for Germany-Japan Next-Generation SAR M. Eineder, C. Minet, A. Parizzi Tokyo,

2 Tandem-L Mission: Two Operation Modes Geosphere: Repeat pass dinsar technique deformation measurements of volcanoes, earthquakes, tectonic, cities, CCS Across track PolInSAR technique selective penetration forest height, biomass L-band DEM

3 Applications Driving Geosphere Requirements Earthquakes Landslides Volcanic Activity Urban Subsidence

4 Motivation Many InSAR examples have been shown - exciting technique - but results opportunistic (selected) and incomplete Why? Temporal decorrelation, e.g. TerraSAR-X, ERS, ASAR Insufficient coverage, e.g. TerraSAR-X Infrequent revisit: ALOS-PALSAR Improper viewing geometry: 1D, layover, Non-systematic data evaluation Challenge: operational & systematic & accurate Earth surface motion monitoring with an InSAR system tailored to needs (cf. SRTM, TanDEM- X): Tandem-L German study / proposal (Moreira et al.) German Helmholtz Study JAXA Study

5 Tandem-L Requirements: Coverage Green: biomass/forest areas (different heights) Red: weekly high strain areas (UNAVO+seismic) + four global deformation images / year, Blue: weekly volcanos Yellow: weekly cities + TBD weekly landslides, CCS sites, etc.

6 Accuracy Challenge: Small Interseismic Strain Rates L x,y,z Requirement formulation: x=1 mm, L=30 km, t=1 year, 1-5 y measurement x 9 1 strain _ rate 3310 y Lt

7 Requirements (Generalized & Compressed) Displacement Accuracy Imaged area Resolution Frequency Global 1cm/50km (5a) approx. 150 Mio. km² 50 m 4/year Earthquakes 1mm/30km (5a) or better 3D-Vector approx. 25 Mio. km² m up to 4/week Volcanoes Anthropogenenic deformation cm-range 3D-Vector approx. 5 Mio. km² 50 m (1-2 m?) mm/cm-range Metropolitan aereas 20 m (1-2 m?) weekly monthly Landslides cm-range 3D-Vector Selected mountainous regions 20 m 1-4/month + DEM times series (from TanDEM-X experience) (1-2 m?): Demand for higher resolution (from TerraSAR-X exp.) Slide 7

8 Biggest Challenge: Atmospheric Water Vapor Horizontal Variogram Hanssen, Netherlands stddev [mm] Wadge, Lanyi Meyer, Warsaw Parizzi, Berlin-C Parizzi, Berlin-X Distance [km] Comparsion of reported atmospheric zenith delay InSAR phase errors 1. R. Hanssen, Radar Interferometry, Kluwer, Dordrecht, F. Meyer, R. Bamler, R. Leinweber, J. Fischer, A Comparative Analysis of Tropospheric Water Vapor Measurements from MERIS and SAR, Proc. of IEEE IGARSS EVINSAR, Proposal for an ESA Earth Exporer Mission, G. Wadge 4. Parizzi, GITEWS Project (DLR internal TN) + high vertical stratification! Countermeassures: >1 year data averaging: PSI (!) Low res tropospheric model for correction of vertical stratification High res. tropospheric model for correction of horizontal variation TBD?

9 Compensation of Stratified Atmospheric Delay Corrected PS Differential Phase Corrected Differential Phase Atmospheric Phase (from ECMWF) Stromboli Volcano, Italy (926m) TerraSAR-X / Baseline: 3.5 m P-9

10 Ionospheric Effects in L-Band IGARSS 11: 2.6 cycles (30 cm) / TECU! Countermeasures: Isolate dispersive Ionosphere using split-bandwidth techniques Use available global TEC maps Use satellite GNSS measurements

11 Ionospheric Correction: Global Models vs. Local GNS CODE In orbit: Satellite to satellite tracking See: C. Gisinger, EGU 2013, Local ionospheric corrections derived from GNSS,

12 Atmospheric Correction - TSX Corner Reflector azimuth offset [m] Correction of IERS 2010 Conventions (earth tides, pole tides, ) Tropospheric delay Ionospheric delay Contintental drift Azimut error σ = 4.3 cm Range error σ= 1.2 cm range offset [m] (bias removed)

13 Enemy of Stacking: Temporal Decorrelation Nisiros / Greece Smaller wavelength lower coherence t e t drift (Zebker, 1992) time constant scales with square of λ ( t) ( 0 ) e Our current model k t k

14 Coherence History in Indonesia (Java/Semarang) 0, 46, 92, 138, 1472 Days (4 Years) Coherence Evolution (small baselines)

15 Coherence Evolution in Indonesia (Java/Semarang) 0, 46, 92, 138, 1472 Days (4 Years) Coherence Evolution 25 PALSAR scenes, small baselines only Optimal observation and processing strategies will be developed

16 Semarang Persistent Scatterer Interferometry Fast motion, not measurable in X-Band! -60 [mm/year] 60

17 Challenge: 2/3-Dimensional Vector Wanted 31 cm uplift y z x 17 cm subsidence Co-seismic Bam / Iran interferogram (ASAR) > 2/3 Observations from different angles required

18 MAP-Estimator for 3D-Vector from n Observations y Kv Model x sin1 cos1 sin1 sin1 cos1 y z y 1 z v x y Cxx ˆˆ xˆ MAP T T 1 K C K C K C y xx C xˆˆ x T 1 K C K 1 Covariance matrix of required 3D Vector estimate, diagonal elements give accuracies in x, y, z (e.g. East, North, Height)

19 Bad 3D-Vector Conditioning with Polar Orbits y =1, (Asc, Desc) (30, 45 ) Inclination=98 Incl.=110 slant 1 Possible improvements: non-polar orbit (other problems: power, ground stations, cryosphere) left and right look direction (requires maneuvers & special satellite design) use second satellite speckle tracking to support azimuth component

20 Earthquake Shifts from InSAR & Speckle Tracking az 3 1 2N 2 Azimuth (correlation) Range (interferometry) Bam / Iran (ASAR)

21 Ground Displacement from Speckle Tracking / Radargrammetry GRSL 2012:

22 Processing Challenge: 3-4 TB/day PGS online processing L0 instrument frames complex image Concept & Product Tree L1 product 8 day interferogram L2 product n cycle 1D rolling buffer Online Quicklook offline or user 2/3 D vector with space time model 1D PSI model select L3 product L3 product long term archive temporary User product

23 Conclusions Challenges are accuracy over large areas ionospheric and tropospheric correction required 2/3D-vector ascending + descending orbits (inclusion of geophysical models) Decorrelation L-band promising Exciting mission - challenges can be met

24 Biosphere Geo-/ Lithosphere Cryo- & Hydrosphere Tandem-L Science Products Forest height Above ground biomass Vertical forest structure Plate tectonics Volcanoes Landslides Deformation Glacier flow Soil moisture Water level change Snow water equivalent Ice structure change Ocean currents Resolutio n m Revisit 16 days - seasonal m weekly m weekly weekly on demand seasonal seasonal weekly Global Digital terrain and surface model m yearly