Click here to go to Geological & Planetary Science's Website Research


Home
Research
Courses
Publications
People
Resources
Activities & More
News

Caltech Logo

>Return to Research

Multi-scale analysis of geophysical fields such as gravity, topography, and crustal deformation
 

My Ph.D. thesis (at MIT, under the supervision of Sean Solomon and Brad Hager) focused on the development and application of a multi-resolution (wavelet-like) perspective on the analysis of geodynamics problems. This particular effort developed a crude approach (analogous to a Gabor Transform, but with a frequency dependent window) for analyzing data on a sphere. More specifically, our approach focused on the case when the original data is provided as spherical harmonic coefficients. The challenge was to maintain some sense of spectral or scale-dependent information, while being able to discriminate processes affecting different regions. My bias is that it is key to stay in the spherical domain and not suffer all the pitfalls of projecting data into a local Cartesian space. Any adopted approach should not violate any sense of conservation of original information. The approach we developed included a Nyquist condition for dealing with truncated spherical harmonic data. We applied these spatio-spectral localization techniques to studies of topographic compensation on Venus (2, 3), Mars (14), and Earth (4).

More recently, we have begun to adopt a multi-resolution perspective in other aspects of our work. We are looking at wavelet-based temporal filtering of continuous GPS data as well as a multi-scale estimation of the surface strain field as constrained by GPS.
 










 
(Top) Free-air gravity in an oblique Mercator projection. (Middle) Gravity profiles along the great circle path shown above. (Bottom) RMS amplitude spectrogram of the free air gravity. The spectrogram is estimated using the full two-dimensional gravity field. We associate the spatially and spectrally concentrated signal at degrees 5-9 centered on Hudson Bay to be the result of incomplete post-glacial rebound. For more details, see ref 4.

This multi-scale strain estimation is of interest since most approaches tend to mix the scales over which strain is estimated, thereby leading to results that can be difficult to interpret. At the same time, we want to estimate strain at the smallest scale allowable by local data distribution. A multi-scale strain estimator may also form the basis of an efficient “event” detector in regions with dense GPS networks such as in Japan, Taiwan, and the U.S..

14 Localized gravity/topography admittance and correlation spectra on Mars: Implications for regional and global evolution, McGovern P. J., S. C. Solomon, D. E. Smith,M. T. Zuber, M. Simons, M. A. Wieczorek, R. J. Phillips, G. A. Neumann, O. Aharonson, and J. W. Head, J. Geophys. Res., 107 (E12), 5136, doi: 10, 1029/2002JE001854, 2002,; Correction, J. Geophys. Res.,109, E07007, doi:10.1029/2004JE002286, July 2004, [PDF]

4 Localization of the gravity field and the signature of glacial rebound, M. Simons and B. H. Hager, Nature, 390, 500-504, 1997. [PDF]

3 Localization of gravity and topography: Constraints on the tectonics and mantle dynamics of Venus, M. Simons, S. C. Solomon, and B. H. Hager, Geophys. J. Int., 131, 24-44, 1997.

2 Global variations in the geoid/topography admittance of Venus, M. Simons, B. H. Hager, and S. C. Solomon, Science, 264, 798-803, 1994. [PDF]

Mark Simons' Paper Collection: Entire paper including figures are all made available online (within the bounds of copyright restrictions).

Last modified 2 October, 2007 | © California Institute of Technology
Click here for Seimological Laboratory Website Click here to go to Geological & Planetary Science's Website