Directional wavelet analysis on the sphere: Application to gravity and topography of the terrestrial planets

[1] The spectral relations (admittance and correlation) between gravity and topography are often used to obtain information on the density structure, flexural support, and heat flow of planetary lithospheres. Mapping spatial variations in these quantities requires spatiospectral analysis techniques. Here we describe the application of a directional, continuous spherical wavelet transform using a wavelet basis constructed from the superposition of azimuthally adjacent complex Morlet wavelets, in a manner similar to the “fan” wavelet developed in the plane. The method is applied to gravity and topography of the Earth, Venus, Mars, and the Moon. The wavelet coefficients are used to compute isotropic and directional wavelet autospectra and cross spectra, which are then combined to form the admittance and correlation functions. The resulting maps offer insights into lithospheric structure of the terrestrial planets. In particular we show that the Earth and Venus have uniformly low positive admittance and high correlation, whereas Mars and the Moon display hemispherical contrasts with large negative and anisotropic coefficients coinciding with lowlands. As has long been known, the two largest impact basins in the inner solar system, the South Pole–Aitken basin on the Moon and the Hellas basin on Mars, display low positive admittance and high correlation, indicating isostatic compensation. In contrast, most other impact basins, particularly the Martian and lunar mascons, show negative coefficients at low wavelet degrees suggesting flexural support by a strong lithosphere. These results imply that, although simple isotropic flexural models can account for most observations, future models may need to incorporate anisotropy as an additional parameter.

[1]  Robert B. Smith,et al.  Strength and rheology of the western U.S. Cordillera , 1995 .

[2]  D. McKenzie,et al.  Measurements of the elastic thickness under ancient lunar terrain , 2005 .

[3]  Bradford H. Hager,et al.  Localization of the gravity field and the signature of glacial rebound , 1997, Nature.

[4]  F. Simons,et al.  Localized spectral analysis on the sphere , 2005 .

[5]  C. Amante,et al.  ETOPO1 arc-minute global relief model : procedures, data sources and analysis , 2009 .

[6]  Jason D. McEwen,et al.  Detection of the integrated Sachs–Wolfe effect and corresponding dark energy constraints made with directional spherical wavelets , 2007 .

[7]  Frederik J. Simons,et al.  Isostatic response of the Australian lithosphere: Estimation of effective elastic thickness and anisotropy using multitaper spectral analysis , 2000 .

[8]  Matthew E. Pritchard,et al.  The Constitution and Structure of the Lunar Interior , 2006 .

[9]  Robert B. Smith,et al.  Flexural rigidity of the Basin and Range‐Colorado Plateau‐Rocky Mountain transition from coherence analysis of gravity and topography , 1994 .

[10]  J. Head,et al.  Lunar mascon basins - Lava filling, tectonics, and evolution of the lithosphere , 1980 .

[11]  Sean C. Solomon,et al.  Localized gravity/topography admittance and correlation spectra on Mars: Implications for regional and global evolution , 2002 .

[12]  Mark A. Wieczorek,et al.  Spatiospectral Concentration on a Sphere , 2004, SIAM Rev..

[13]  J. Kirby,et al.  A reassessment of spectral Te estimation in continental interiors: The case of North America , 2009 .

[14]  Frederik J. Simons,et al.  Spatiospectral localization of isostatic coherence anisotropy in Australia and its relation to seismic anisotropy: Implications for lithospheric deformation , 2003 .

[15]  J. D. McEwen,et al.  A high-significance detection of non-Gaussianity in the Wilkinson Microwave Anisotropy Probe 1-yr data using directional spherical wavelets , 2004, astro-ph/0406604.

[16]  Jean-Claude Mareschal,et al.  Wavelet analysis of the coherence between Bouguer gravity and topography: application to the elastic thickness anisotropy in the Canadian Shield , 2007 .

[17]  M. Gurnis,et al.  Mantle upwellings above slab graveyards linked to the global geoid lows , 2010 .

[18]  Syed Twareque Ali,et al.  Two-Dimensional Wavelets and their Relatives , 2004 .

[19]  Belgium,et al.  Correspondence principle between spherical and euclidean wavelets , 2005, astro-ph/0502486.

[20]  P. Vandergheynst,et al.  Wavelets on the 2-sphere: A group-theoretical approach , 1999 .

[21]  J. Stewart,et al.  Wavelet transform mapping of effective elastic thickness and plate loading: Validation using synthetic data and application to the study of , 2001 .

[22]  J. Kirby,et al.  Mapping the mechanical anisotropy of the lithosphere using a 2D wavelet coherence, and its application to Australia , 2006 .

[23]  Donald W. Forsyth,et al.  Subsurface loading and estimates of the flexural rigidity of continental lithosphere , 1985 .

[24]  Maria T. Zuber,et al.  Thickness of the Martian crust: Improved constraints from geoid-to-topography ratios , 2004 .

[25]  O. Aharonson,et al.  Mega-impact formation of the Mars hemispheric dichotomy , 2008, Nature.

[26]  Global Surveyor Topography and Gravity Internal Structure and Early Thermal Evolution of Mars from Mars , 2011 .

[27]  G. Balmino,et al.  The gravity field of Mars: results from Mars Global Surveyor. , 1999, Science.

[28]  P. Lognonné,et al.  Constraints on the Martian lithosphere from gravity and topography data and implication for volcanism and possible mantle heterogeneities , 2005 .

[29]  J. D. McEwen,et al.  Detection of the ISW effect and corresponding dark energy constraints , 2006 .

[30]  J. F. Kirby,et al.  Which wavelet best reproduces the Fourier power spectrum? , 2005, Comput. Geosci..

[31]  Suzanne E. Smrekar,et al.  Global mapping of crustal and lithospheric thickness on Venus , 2006 .

[32]  R. J. Willemann,et al.  Role of membrane stresses in the support of planetary topography , 1981 .

[33]  D. McKenzie The Relationship between Topography and Gravity on Earth and Venus , 1994 .

[34]  K. Lambeck,et al.  Gravity fields of the terrestrial planets: Long‐wavelength anomalies and tectonics , 1980 .

[35]  I. Velicogna,et al.  On the recovery of effective elastic thickness using spectral methods: Examples from synthetic data and from the Fennoscandian Shield , 2003 .

[36]  Qinghui Liu,et al.  An improved lunar gravity field model from SELENE and historical tracking data: Revealing the farside gravity features , 2010 .

[37]  Sean C. Solomon,et al.  Heterogeneities in the thickness of the elastic lithosphere of Mars - Constraints on heat flow and internal dynamics , 1990 .

[38]  David N. Barnett,et al.  Elastic Thickness Estimates for Venus Using Line of Sight Accelerations from Magellan Cycle 5 , 2000 .

[39]  Sean C. Solomon,et al.  Localization of gravity and topography: constraints on the tectonics and mantle dynamics of Venus , 1997 .

[40]  Mark A. Richards,et al.  A conceptual model for the relationship between coronae and large-scale mantle dynamics on Venus , 2003 .

[41]  J W Head,et al.  Internal structure and early thermal evolution of Mars from Mars Global Surveyor topography and gravity. , 2000, Science.

[42]  D. McKenzie Estimating Te in the presence of internal loads , 2002 .

[43]  David E. Smith,et al.  The global topography of Mars and implications for surface evolution. , 1999, Science.

[44]  L. Dorman,et al.  Experimental isostasy: 2. An isostatic model for the U.S.A. derived from gravity and topographic data , 1970 .

[45]  L. Dorman,et al.  Experimental isostasy: 3. Inversion of the isostatic Green function and lateral density changes , 1972 .

[46]  Jafar Arkani-Hamed,et al.  Lunar mascons revisited , 1997 .

[47]  Walter H. F. Smith,et al.  Global Sea Floor Topography from Satellite Altimetry and Ship Depth Soundings , 1997 .

[48]  David A. Yuen,et al.  Continuous wavelet-like filter for a spherical surface and its application to localized admittance function on Mars , 2003 .

[49]  R. Parker,et al.  Isostatic compensation on a continental scale: local versus regional mechanisms , 1977 .

[50]  D. Yuan,et al.  The relationship between Martian gravity and topography , 2002 .

[51]  J. Oberst,et al.  Lunar Global Shape and Polar Topography Derived from Kaguya-LALT Laser Altimetry , 2009, Science.

[52]  Michael P. Hobson,et al.  Fast Directional Continuous Spherical Wavelet Transform Algorithms , 2005, IEEE Transactions on Signal Processing.

[53]  H. Uno,et al.  Mechanical controls on the deformation of continents at convergent margins , 2007 .

[54]  M. Zuber,et al.  The Borealis basin and the origin of the martian crustal dichotomy , 2008, Nature.

[55]  Walter H. F. Smith,et al.  Marine gravity anomaly from Geosat and ERS 1 satellite altimetry , 1997 .

[56]  Strength of Martian lithosphere beneath large volcanoes , 2000 .

[57]  K. Gorski,et al.  HEALPix: A Framework for High-Resolution Discretization and Fast Analysis of Data Distributed on the Sphere , 2004, astro-ph/0409513.

[58]  R. Parker,et al.  Isostasy in Australia and the Evolution of the Compensation Mechanism , 1978, Science.

[59]  D. McKenzie,et al.  Estimates of the effective elastic thickness of the continental lithosphere from Bouguer and free air gravity anomalies , 1997 .

[60]  L. Dorman,et al.  Experimental isostasy: 1. Theory of the determination of the Earth's isostatic response to a concentrated load , 1970 .

[61]  Roger J. Phillips,et al.  Potential anomalies on a sphere: Applications to the thickness of the lunar crust , 1998 .

[62]  G. Houseman,et al.  ELASTIC THICKNESS ESTIMATES FOR CORONAE ASSOCIATED WITH CHASMATA ON VENUS. , 2005 .

[63]  D. Sandwell,et al.  Lithospheric flexure on Venus , 1994 .

[64]  A. Kalnins,et al.  Thin elastic shells , 1967 .

[65]  J. D. McEwen,et al.  A high-significance detection of non-Gaussianity in the WMAP 5-yr data using directional spherical wavelets , 2006, 0803.2157.

[66]  F. Nimmo,et al.  Implications of an impact origin for the martian hemispheric dichotomy , 2008, Nature.