Spherical harmonic analysis and synthesis using FFT: Application to temporal gravity variation

FFT and complex algebra-based methods of spherical harmonic analysis and synthesis are presented. Two computer programs in FORTRAN are developed based on the methods. Both general and special cases are discussed. Special cases involve the analyses of gravity changes of the hydrological origin and the atmospheric origin. Functionals of the Earth's gravity field such as gravity anomaly and geoidal height can also be computed via synthesis. Thermal-corrected sea level anomaly from TOPEX/Poseidon and atmospheric pressure from ECMWF are used to compute changes of geopotential coefficients due to oceanic and atmospheric mass redistributions. Interesting phenomena in the changes of geopotential coefficients have been identified. The two computer programs can facilitate analyses and syntheses of gravity products from satellite missions such as GRACE.

[1]  Walter H. F. Smith,et al.  New version of the generic mapping tools , 1995 .

[2]  Gabriele Steidl,et al.  Fast and stable algorithms for discrete spherical Fourier transforms , 1998 .

[3]  Fernando Sansò,et al.  International Summer School of Theoretical Geodesy on “Theory of Satellite Geodesy and Gravity Field Determination” , 1990 .

[4]  M. Watkins,et al.  GRACE Measurements of Mass Variability in the Earth System , 2004, Science.

[5]  B. Chao,et al.  Detection of a Large-Scale Mass Redistribution in the Terrestrial System Since 1998 , 2002, Science.

[6]  O. Colombo.  Numerical Methods for Harmonic Analysis on the Sphere , 1981 .

[7]  B. Tapley,et al.  Oceanic effects on polar motion determined from an ocean model and satellite altimetry: 1993-2001 , 2004 .

[8]  N. K. Pavlis,et al.  The Development of the Joint NASA GSFC and the National Imagery and Mapping Agency (NIMA) Geopotential Model EGM96 , 1998 .

[9]  Reiji Suda,et al.  A fast spherical harmonics transform algorithm , 2002, Math. Comput..

[10]  M. K. Paul.  Recurrence relations for integrals of Associated Legendre functions , 1978 .

[11]  Martin J. Mohlenkamp.  A fast transform for spherical harmonics , 1997 .

[12]  Sean S. B. Moore,et al.  FFTs for the 2-Sphere-Improvements and Variations , 1996 .

[13]  F. Bryan,et al.  Time variability of the Earth's gravity field: Hydrological and oceanic effects and their possible detection using GRACE , 1998 .

[14]  Gary A Dilts,et al.  Computation of spherical harmonic expansion coefficients via FFT's , 1985 .

[15]  I. Fukumori,et al.  Recent Earth Oblateness Variations: Unraveling Climate and Postglacial Rebound Effects , 2002, Science.

[16]  Theodossios Engelis,et al.  Spherical harmonic expansion of the Levitus Sea surface topography , 1987 .

[17]  C. Hwang.  Gravity recovery using COSMIC GPS data: application of orbital perturbation theory , 2001 .

[18]  P. A. Cross,et al.  Lecture notes in Earth sciences: Vol. 12. S. Turner (Editor), Applied Geodesy VIII, Springer, Berlin, F.R.G., 1987, 393pp, DM78.00, ISBN 3 540 182195 , 1989 .

[19]  J. Wahr,et al.  The viscoelastic relaxation of a realistically stratified earth, and a further analysis of postglacial rebound , 1995 .

[20]  D. Rockmore,et al.  Computational Harmonic Analysis for Tensor Fields on the Two-Sphere , 2000 .