Elastic tidal response of a laterally heterogeneous planet: a complete perturbation formulation

S U M M A R Y Constraining laterally varying structures in planetary interiors is important for understanding both the composition and the internal dynamics of a planet. Recognizing that seismic imaging technique is currently only viable for studying the Earth’s interior structures, methods that can be supported by advanced space geodetic techniques may become alternatives to ‘image’ the interiors of other planets. The method of tidal tomography is one possibility, and it relies on high precision measurement of the response of a planet to its body tide. However, it is essential to develop an efficient analytical tool that computes the dependence of tidal response to 3-D interior structures. In this paper, we present a complete formulation of such an analytical tool, which calculates to high accuracy the tidal response of a terrestrial planet with lateral heterogeneities in its elastic and density structures. We treat the lateral heterogeneities as small perturbations and derive the governing equations based on the perturbation theory. In a spherical harmonic representation, equations at each order of perturbation are reduced into multiple matrix equations at harmonics that are allowed by mode couplings, and the total response equals the sum of all those single-harmonic responses, which can be solved semi-analytically. We test our perturbation method by applying it to the Moon with a harmonic degree-1 mantle structure for which the perturbation solutions of the tidal response are compared with those from a fully numerical method. The remarkable agreement between results from these two methods validates the perturbation method. As an example, we then use the perturbation method to evaluate the impact of lunar crustal thickness variations on tidal response of the Moon. We find that lunar crust produces much smaller degree-3 tidal responses than a relatively weak degree-1 structure in the deep lunar mantle. Our calculations show that degree-3 tidal response measurements may hold key constraints on possible degree-1 mantle structure of the Moon, as suggested from previous modelling results.

[1]  J. Tromp,et al.  A normal mode treatment of semi-diurnal body tides on an aspherical, rotating and anelastic Earth , 2015 .

[2]  J. Wahr,et al.  A perturbation method and its application: elastic tidal response of a laterally heterogeneous planet , 2014 .

[3]  David E. Smith,et al.  Lunar interior properties from the GRAIL mission , 2014 .

[4]  David E. Smith,et al.  High‒degree gravity models from GRAIL primary mission data , 2013 .

[5]  David E. Smith,et al.  The JPL lunar gravity field to spherical harmonic degree 660 from the GRAIL Primary Mission , 2013 .

[6]  A. Love,et al.  Some Problems of Geodynamics , 2013 .

[7]  David E. Smith,et al.  Gravity Field of the Moon from the Gravity Recovery and Interior Laboratory (GRAIL) Mission , 2013, Science.

[8]  Sami W. Asmar,et al.  The Crust of the Moon as Seen by GRAIL , 2012, Science.

[9]  J. Wahr,et al.  Computations of the viscoelastic response of a 3-D compressible Earth to surface loading: an application to Glacial Isostatic Adjustment in Antarctica and Canada , 2012 .

[10]  J. Wahr,et al.  Can tidal tomography be used to unravel the long‐wavelength structure of the lunar interior? , 2012 .

[11]  S. Zhong,et al.  Correlation of deep moonquakes and mare basalts: Implications for lunar mantle structure and evolution , 2012 .

[12]  Doris Breuer,et al.  Asymmetric thermal evolution of the Moon , 2012 .

[13]  Haje Korth,et al.  The Global Magnetic Field of Mercury from MESSENGER Orbital Observations , 2011, Science.

[14]  Renee C. Weber,et al.  Seismic Detection of the Lunar Core , 2011, Science.

[15]  R. Pappalardo,et al.  Modeling stresses on satellites due to nonsynchronous rotation and orbital eccentricity using gravitational potential theory , 2009 .

[16]  J. Mitrovica,et al.  Body tides on a 3-D elastic earth: Toward a tidal tomography , 2009 .

[17]  G. Neukum,et al.  Cassini Observes the Active South Pole of Enceladus , 2006, Science.

[18]  Rosaly M. C. Lopes,et al.  Cassini Encounters Enceladus: Background and the Discovery of a South Polar Hot Spot , 2006, Science.

[19]  M. Watkins,et al.  The gravity recovery and climate experiment: Mission overview and early results , 2004 .

[20]  A. Paulson,et al.  Three-dimensional finite-element modelling of Earth's viscoelastic deformation: effects of lateral variations in lithospheric thickness , 2003 .

[21]  Paul G. Lucey,et al.  Iron abundances on the lunar surface as measured by the Lunar Prospector gamma‐ray and neutron spectrometers , 2002 .

[22]  M. Ritzwoller,et al.  Monte-Carlo inversion for a global shear-velocity model of the crust and upper mantle , 2002 .

[23]  Maria T. Zuber,et al.  The crust and mantle of Mars , 2001, Nature.

[24]  M. Zuber,et al.  A dynamic origin for the global asymmetry of lunar mare basalts , 2000 .

[25]  J. Tromp,et al.  Surface loading of a viscoelastic earth—I. General theory , 1999 .

[26]  J. Tromp,et al.  Surface loading of a viscoelastic earth—II. Spherical models , 1999 .

[27]  Carle M. Pieters,et al.  Mineralogy of the lunar crust: Results from Clementine , 1999 .

[28]  J. Tromp,et al.  Theoretical Global Seismology , 1998 .

[29]  S. Maurice,et al.  Global elemental maps of the moon: the Lunar Prospector gamma-Ray spectrometer. , 1998, Science.

[30]  Wei-jia Su,et al.  Simultaneous inversion for 3-D variations in shear and bulk velocity in the mantle , 1997 .

[31]  David E. Smith,et al.  Topography of the Moon from the Clementine lidar , 1997 .

[32]  David E. Smith,et al.  The Shape and Internal Structure of the Moon from the Clementine Mission , 1994, Science.

[33]  B. Hager,et al.  Long-wavelength variations in Earth’s geoid: physical models and dynamical implications , 1989, Philosophical Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences.

[34]  W. Farrell Deformation of the Earth by surface loads , 1972 .

[35]  I. M. Longman A Green's function for determining the deformation of the Earth under surface mass loads: 2. Computations and numerical results , 1963 .

[36]  I. M. Longman A Green's function for determining the deformation of the Earth under surface mass loads: 1. Theory , 1962 .

[37]  V. S. Scott,et al.  The Lunar Orbiter Laser Altimeter Investigation on the Lunar Reconnaissance Orbiter Mission , 2010 .

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

[39]  Yosio Nakamura,et al.  Farside deep moonquakes and deep interior of the Moon , 2005 .