1 / 49

Satellite Altimetry

Satellite Altimetry. Ole B. Andersen. Content. The radar altimetric observations (1) : Altimetry data Contributors to sea level Retracking Crossover adjustment From altimetry to Gravity and Geoid (2): Geodetic theory FFT for global gravity fields GRAVSOFT Applications.

kellyrichard
Download Presentation

Satellite Altimetry

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Satellite Altimetry Ole B. Andersen.

  2. Content • The radar altimetric observations (1): • Altimetry data • Contributors to sea level • Retracking • Crossover adjustment • From altimetry to Gravity and Geoid (2): • Geodetic theory • FFT for global gravity fields • GRAVSOFT • Applications. • Accuray Assesment • Applications IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 2

  3. IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 3

  4. Sampling the Sea Surface IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 4

  5. Sampling the sea surface. 1 Day 3 Days IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 5

  6. Orbit Parameters The coverage of the sea surface depends on the orbit parameters (inclination of the orbit plane and repeat period). TOPEX/JASON - 10 Days IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 6

  7. ERM – GM data. ERM Data TOPEX/JASON – (280 km) ERS/ENVISAT (80 km) Geodetic Mission GEOSAT (15 Month) Drift ERS-1 (11 Month) 2 x 168 days repeat Equally spacing GEOSAT+ERS GM data is ESSENTIAL for high resolution Gravity Field mapping. IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 7

  8. Altimetric Observations Accurate ranging to the sea surface is Based on accurate time-determination Typical ocean waveforms Registred at 2000 Hz (every 3 meters) Averaged into 20 Hz. The 20 Hz height values are Also noisy Averaged to 1 Hz values (7 km averaging). IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 8

  9. Different Surfaces – Different Retrackers Ocean Echoes Extreme Terrain Desert – Australia Inland Water (River – lake) R, P. Berry, J. Freeman IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 9

  10. Benefit of retracking The Sea Ice retracker picks up many data in i.e. polar regions that is not ocean retracket. This retracker also picks up data near coast and in currents (correct echos are similar in shape) Large distribution of patch waveforms (44%) within 5 km of the coastline, and around 10% from about 25 km from the coast can be seen. At about 50 km from coastline only 60% of waveforms in this region are ocean retracked. IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 10

  11. Altimetry observes the sea surface height (SSH) Satellite Altimetry If the satellite is accurately positioned The orbital height of the space craft minus the altimeter radar ranging to the sea surface corrected for path delays and environmental corrections Yields the sea surface height: where Nis the geoid height above the reference ellipsoid, is the ocean topography, e is the error The Sea surface height mimicks the geoid. IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 11

  12. Altimetric observations The magnitudes of the contributors ranges up to The geoid NREF +/- 100 meters Terrain effect NDTM +/- 30 centimetres Residual geoidN +/- 2 meters Mean dynamic topography MDT+/- 2 meters Time varying Dyn topography (t) +/- 5 meters. (Tides + storms + El Nino……) What We want for Global Gravity is: So we need to ”take out” the rest. IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 12

  13. Remove - Restore. • “ Take These Out” • Remove-restore technique – enhancing signal to noise. • “Remove known signals and restore their effect subsequently” • Remove a global spherical harmonic geoid model (PGM04/06) • Remove terrain effect • Remove Mean Dynamic Topography (MDT) • Compute Gravity • Restore PGM04/06 global gravity field (Pavlis) • Restore the Terrain effect • No Restoring from Mean Dynamic Topography GEOID signal +/- 100 meters IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 13

  14. The Mean Dynamic Topography (MDT) +/- 2 Meters - Gives up to 3-4 mGal effect mGal 0 -3 3 SSH ~ G + MDT -> So the MDT can be determined like MDT = MSS-N IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 14

  15. Time Varying Signal Tides contribute nearly 80% to sea level variability. Removed using Ocean tide Model (GOT 2000X) Time variable signals are averaged out in ERM data but not in GM data. Errors • eorbit is the radial orbit error • etides is the errors due to remaining tidal errors • erange is the error on the range corrections. • eretrak is the errors due to retracking • enoise is the measurement noise. IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 15

  16. Errors+time varying signals. • ERM data. Most time+error average out. • Geodetic mission data (t) is not reduced • Must limit errors to avoid ”orange skin effect” • NOTICE ERRORS ARE LONG WAVELENGTH IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 16

  17. Crossover – Sea surface Slopes Enhancing residual geoid height signal for gravity (removing time variability) • Limit long wavelength contribution (time + error signal FOR GM DATA). • Using crossover adjustment. • Motivation • Assumption: the residual geoid signal is stationary at each location so residual geoid observations • should be the same on ascending and descending tracks at crossing locations. • Timevarying Dynamic sea level + orbital related signals should not be the same, and will be removed • Using Sea surface Slopes. • Motivation • Easier computation (no need for crossover adjustmenst • Theoretically straight forward wrt gravity field computation. • Time varying signal is not reduced. • Transformation from along track to east-west north-south is problematic at the Equator and at turning lat. IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 17

  18. Crossover Adjustment • dk=hi‑hj. • d=Ax+v • where x is vector containing the unknown parameters for the track-related errors. • v is residuals that we wish to minimize • Least Squares Solution to this is • Constraint is needed cTx=0 • Problem of Null space – Rank • Bias (rank=1) – mean bias is zero • Bias+Tilt (Rank = 4) Constrain to prior surf. IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 18

  19. Before - Xover IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 19

  20. After Crossover IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 20

  21. Data are now ready for computing gravity / geoid. • Corrected the range for as many known signals as possible. • Retracking enhances amount and quality of data • Removed Long wavelength Geoid part – will be restored. • XOVER: Limited errors + time varying signal (Long wavelength). • Small long wavelength errors can still be seen in sea surface heights. IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 21

  22. Part 2. Geodesy • The radar altimetric observations (1): • Altimetry data • Contributors to sea level • Crossover adjustment • From altimetric heights to Gravity (2): • Geodetic theory • FFT for global gravity field determination • DNSC Global Gravity Field. • Applications: • Accuracy assesment • Applications IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 22

  23. The Anomalous Potential. • The anomalous potential T is the difference between the actual gravity potential W and the normal potential U • T is a harmonic function outside the masses of the Earth satisfying • (²T = 0) Laplace (outside the masses) • (²T = -4) Poisson (inside the masses ( is density)) • T Harmonic -> Expand T in spherical harmonic functions: • Pij are associated Legendre's functions of degree i and order j • C+D are surface spherical harmonic functions (what you distribute) • Geoid heights, multip by 1/γ, Gravity anomalies, multip by (i-1)/R IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 23

  24. Geoid to Gravity Geoid N and T (Bruns Formula)N can be expressed in terms of a linear functional applied on T (γ is the normal gravity)Gravity and T Deflection of the vertical (n,e)Deflection of the vertical is related to geoid slope Geoid slopes (east, west) can be obtained from altimetry by tranformingthe along-track slopes to east-west slopes. IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 24

  25. Gravity from altimetry. Three feasable ways 1) Integral formulas (Stokes + Vening Meinesz + Inverse) Requires extensive computations over the whole earth. Replace analytical integrals with grids and is combined with FFT 2) Fast Fourier Techniques. Requires gridded data (will return to that). Very fast computation. Presently the most widely used method. 3) Collocation. Requires big computers. Do not require gridding. Ongoing investigating this approach IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 25

  26. Using FFT • Flat Earth approx is valid (2-300 km from computation point, (More by Sideris, (1997)). • The geodetic relations with T are then • Where F is the 2D planar FFT transform • DNSC06. • Gravity is estimated in small boxes (3 by 10° ) and pathed together globally 2800 cells • Around 100 million 1 sec ssh observations. • DNSC approach is highly parallel. Computation time is around 1-2 weeks IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 26

  27. From height to gravity using 2D FFT • An Inverse Stokes problem • High Pass filter operation enhance high frequency. • Optimal filter was designed to handle white noise + power spectral decay obtained using • Frequency domain LSC with a Wiener Filter (Forsberg and Solheim, 1997) • Power spectral decay follows Kaulas rule (k-4) • Resolution is where wavenumber k yields (k) = 0.5 • For DNSC 12 km is used IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 27

  28. Data and software • Satellite Altimetry (Major points). • Altimetry Pathfinder: http://topex-www.jpl.nasa.gov/ • RADS/NOAA (Remko): http://rads.tudelft.nl/rads/rads.shtml • NASA, ESA (Raw data). • DNSC Global Fields (ftp.spacecenter.dk) • Marine Gravity (2 min res). (http://www.spacecenter.dk/data) • Mean Sea Surface (2 min res) http://www.spacecenter.dk/data • Bathymetry (2 min res) http://www.spacecenter.dk/data • Ocean Tide model (30 min res) http://www.spacecenter.dk/data • Software (with geoid school). • GRAVSOFT • Least Squares collocation: GEOCOL, EMPCOV, COVFIT • Interpolation: GEOGRID • FFT: GEOFOUR, SPFOUR IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 28

  29. DNSC Global Marine Gravity Grid IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 29

  30. Gravity and Earth Processes IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 30

  31. IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 31

  32. Marine Gravity Prior to Satellite Altimetry IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 32

  33. IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 33

  34. DNSC05A NOAA12 0.48 0.62 5.79 4.95 82.20 47.88 GSFC00.1 DNSC05-OC 0.50 0.68 6.14 5.21 51.18 89.91 NCU01 DNSC05-DOT36 0.49 0.79 4.79 6.10 92.13 45.19 North Atlantic Region Gravity diff With KMS02 IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 34

  35. Part 3: Applications. • Sea level Change • Ocean Tides • Land Hydrology • Sea level changes: • Global coverage – open ocean • Uniform Geocentric reference • About 12 years of T/P time series used • Spatial characteristics • Calibration needed at tide gauges IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 35

  36. Sea level change / Climate From Satellite altimetry.

  37. Part 3: Applications. • Sea level changes: • Global coverage – open ocean • Uniform Geocentric reference • About 12 years of T/P time series used • Spatial characteristics • Calibration needed at tide gauges • Sea level change is currently increasing from 2.8 to 3.0 mm / year indicating acceleration….. IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 37

  38. 12 Years Sea Level Change (1993-2004) 12 Years Sea Surface Temperature change. IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 38

  39. Global Sea level change. IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 39

  40. OCEAN TIDES From Satellite altimetry.

  41. Fascinating Ocean Tides Tidal Range in Bay of Fundy and English Channel is 15 meters IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 41

  42. Tides can be ”dangerous” -BUT TIDES CAN BE PREDICTED. IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 42

  43. Tidal Forces. Force = Difference (P1-O) is the Tide generating force = At P2 the force away from the moon is = At P3 the force is directed towards O In Addition we have the centrifugal force which must be added. IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 43

  44. Alias Periods IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 44

  45. Ocean Tides - M2 loop IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 45

  46. GRACE GRACE twin satellites launched March 2002. Status: Mission successfully 4 Years in Space. 40 Monthly Level-2 solutions Use for climate studies to control global water balance better than hydrological models. Various corrections for atmosphere+ Tides +….. applied leaving Hydrology as largest contributor to gravity field variations. IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 46

  47. Land Hydrology – The Amazon Satellite altimetry In rivers Retracking is Essential (P. Berry–De Montford) IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 47

  48. Hydrology from GRACE + Satellite Altimetry GRACE Satellite Altimetry IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 48

  49. Thank You IAG 2006 Geoid School | Copenhagen 30 maj 2006 | OA | side 49

More Related