Hamilton’s principle and normal mode coupling in an aspherical planet with a fluid core

Seismic free oscillations, or normal modes, provide a convenient tool to calculate low-frequency seismograms in heterogeneous Earth models. A procedure called ‘full mode coupling’ allows the seismic response of the Earth to be computed. However, in order to be theoretically exact, such calculations must involve an infinite set of modes. In practice, only a finite subset of modes can be used, introducing an error into the seismograms. By systematically increasing the number of modes beyond the highest frequency of interest in the seismograms, we investigate the convergence of full-coupling calculations. As a rule-of-thumb, it is necessary to couple modes 1–2 mHz above the highest frequency of interest, although results depend upon the details of the Earth model. This is significantly higher than has previously been assumed. Observations of free oscillations also provide important constraints on the heterogeneous structure of the Earth. Historically, this inference problem has been addressed by the measurement and interpretation of splitting functions. These can be seen as secondary data extracted from low frequency seismograms. The measurement step necessitates the calculation of synthetic seismograms, but current implementations rely on approximations referred to as self- or group-coupling and do not use fully accurate seismograms. We therefore also investigate whether a systematic error might be present in currently published splitting functions. We find no evidence for any systematic bias, but published uncertainties must be doubled to properly account for the errors due to theoretical omissions and regularization in the measurement process. Correspondingly, uncertainties in results derived from splitting functions must also be increased. As is well known, density has only a weak signal in low-frequency seismograms. Our results suggest this signal is of similar scale to the true uncertainties associated with currently published splitting functions. Thus, it seems that great care must be taken in any attempt to robustly infer details of Earth's density structure using current splitting functions.

[1]  T. Lay.  Deep Earth Structure – Lower Mantle and D″ , 2015 .

[2]  M. D. Hoop,et al.  Variational formulation of the earth's elastic-gravitational deformations under low regularity conditions , 2017, 1702.04741.

[3]  Andrew P. Valentine,et al.  The impact of approximations and arbitrary choices on geophysical images , 2016 .

[4]  Ezio Faccioli,et al.  2d and 3D elastic wave propagation by a pseudo-spectral domain decomposition method , 1997 .

[5]  John H. Woodhouse,et al.  Iteration method to determine the eigenvalues and eigenvectors of a target multiplet including full mode coupling , 2004 .

[6]  Felix E. Browder,et al.  On the spectral theory of elliptic differential operators. I , 1961 .

[7]  D. Komatitsch,et al.  The spectral element method: An efficient tool to simulate the seismic response of 2D and 3D geological structures , 1998, Bulletin of the Seismological Society of America.

[8]  H. Simpson,et al.  On the positivity of the second variation in finite elasticity , 1987 .

[9]  J. Tromp,et al.  Synthetic free-oscillation spectra: an appraisal of various mode-coupling methods , 2015 .

[10]  C. Truesdell,et al.  The Non-Linear Field Theories Of Mechanics , 1992 .

[11]  D. Al‐Attar,et al.  Particle relabelling transformations in elastodynamics , 2016 .

[12]  John H. Woodhouse,et al.  Calculation of normal mode spectra in laterally heterogeneous earth models using an iterative direct solution method , 2012 .

[13]  G. Uhlmann,et al.  Stability estimates for the hyperbolic Dirichlet to Neumann map in anisotropic media , 1998 .

[14]  Jeffrey Park,et al.  The subspace projection method for constructing coupled-mode synthetic seismograms , 1990 .

[15]  David R. Smith,et al.  Design of electromagnetic cloaks and concentrators using form-invariant coordinate transformations of Maxwell’s equations , 2007, 0706.2452.

[16]  Freeman Gilbert,et al.  Coupled free oscillations of an aspherical, dissipative, rotating Earth: Galerkin theory , 1986 .

[17]  A. Dziewoński,et al.  Seismic Surface Waves and Free Oscillations in a Regionalized Earth Model , 1982 .

[18]  T. R. Hughes,et al.  Mathematical foundations of elasticity , 1982 .

[19]  J. Virieux.  P-SV wave propagation in heterogeneous media: Velocity‐stress finite‐difference method , 1986 .

[20]  W. Noll.  The Foundations of Mechanics and Thermodynamics: Selected Papers , 1974 .

[21]  F. Dahlen.  Elastic Dislocation Theory for a Self‐Gravitating Elastic Configuration with an Initial Static Stress Field , 1972 .

[22]  Peter Moczo,et al.  Finite-difference technique for SH-waves in 2-D media using irregular grids-application to the seismic response problem , 1989 .

[23]  George E. Backus,et al.  Moment Tensors and other Phenomenological Descriptions of Seismic Sources—I. Continuous Displacements , 1976 .

[24]  D. Komatitsch,et al.  Spectral-element simulations of global seismic wave propagation: II. Three-dimensional models, oceans, rotation and self-gravitation , 2002 .

[25]  D. L. Anderson,et al.  Preliminary reference earth model , 1981 .

[26]  Plamen Stefanov,et al.  Rigidity for metrics with the same lengths of geodesics , 1998 .

[27]  R. Geller,et al.  Inversion for laterally heterogeneous upper mantle S-wave velocity structure using iterative waveform inversion , 1993 .

[28]  S. Agmon,et al.  Estimates near the boundary for solutions of elliptic partial differential equations satisfying general boundary conditions. I , 1959 .

[29]  A. Ibrahimbegovic.  Nonlinear Solid Mechanics , 2009 .

[30]  R. Abraham,et al.  Manifolds, Tensor Analysis, and Applications , 1983 .

[31]  Inversion for laterally heterogeneous earth structure using a laterally heterogeneous starting model: preliminary results , 2007 .

[32]  大森 英樹,et al.  Infinite dimensional Lie transformations groups , 1974 .

[33]  S. Stewart,et al.  The structure of terrestrial bodies: Impact heating, corotation limits, and synestias , 2017, 1705.07858.

[34]  B. Valette,et al.  Influence of liquid core dynamics on rotational modes , 2009 .

[35]  F. Dahlen.  The Normal Modes of a Rotating, Elliptical Earth , 1968 .

[36]  F. Dahlen.  Elastic Dislocation Theory for a Self‐Gravitating Elastic Configuration with an Initial Static Stress Field ii. Energy Release , 1973 .

[37]  P. Lognonné.  Normal modes and seismograms in an anelastic rotating Earth , 1991 .

[38]  Darryl D. Holm,et al.  Geometric Mechanics and Symmetry: From Finite to Infinite Dimensions , 2009 .

[39]  F. A. Dahleo.  The Passive Influence of the Oceans upon the Rotation of the Earth , 2009 .

[40]  Matti Lassas,et al.  On nonuniqueness for Calderón’s inverse problem , 2003 .

[41]  Emmanuel Chaljub,et al.  Spectral element modelling of three-dimensional wave propagation in a self-gravitating Earth with an arbitrarily stratified outer core , 2003, physics/0308102.

[42]  G. Ekström,et al.  The relationships between large‐scale variations in shear velocity, density, and compressional velocity in the Earth's mantle , 2016 .

[43]  Place Jussieu.  Matrix methods for generally stratified media , 1976 .

[44]  David R. Smith,et al.  Controlling Electromagnetic Fields , 2006, Science.

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

[46]  F. Dahlen.  The Normal Modes of a Rotating, Elliptical Earth—II Near-Resonance Multiplet Coupling , 1969 .

[47]  J. Woodhouse,et al.  On the parametrization of equilibrium stress fields in the Earth , 2010 .

[48]  J. Tromp,et al.  Summation of the Born series for the normal modes of the Earth , 1990 .

[49]  F. A. Dahlen,et al.  The Effect of A General Aspherical Perturbation on the Free Oscillations of the Earth , 1978 .

[50]  S. Agmon.  On the eigenfunctions and on the eigenvalues of general elliptic boundary value problems , 1962 .

[51]  Jeroen Tromp,et al.  A spectral-infinite-element solution of Poisson's equation: an application to self gravity , 2017, 1706.00855.

[52]  D. Boore,et al.  Finite Difference Methods for Seismic Wave Propagation in Heterogeneous Materials , 1972 .

[53]  J. Tromp,et al.  Normal-mode and free-Air gravity constraints on lateral variations in velocity and density of Earth's mantle , 1999, Science.

[54]  Joseph S. Resovsky,et al.  Probabilistic Tomography Maps Chemical Heterogeneities Throughout the Lower Mantle , 2004, Science.

[55]  J. Marsden,et al.  Groups of diffeomorphisms and the motion of an incompressible fluid , 1970 .

[56]  Gerhard A. Holzapfel,et al.  Nonlinear Solid Mechanics: A Continuum Approach for Engineering Science , 2000 .

[57]  B. Romanowicz,et al.  Modelling of coupled normal modes of the Earth: the spectral method , 1990 .

[58]  A. Valentine,et al.  Exact free oscillation spectra, splitting functions and the resolvability of Earth's density structure , 2018 .

[59]  J. Tromp,et al.  Tidal tomography constrains Earth’s deep-mantle buoyancy , 2017, Nature.

[60]  John H. Woodhouse,et al.  Theoretical free-oscillation spectra: the importance of wide band coupling , 2001 .

[61]  James Lowry Thompson,et al.  Some existence theorems for the traction boundary value problem of linearized elastostatics , 1969 .

[62]  N. Takeuchi.  Finite boundary perturbation theory for the elastic equation of motion , 2005 .

[63]  F. Gilbert.  Excitation of the Normal Modes of the Earth by Earthquake Sources , 1971 .

[64]  G. Backus.  Converting Vector and Tensor Equations to Scalar Equations in Spherical Coordinates , 1967 .

[65]  J. Woodhouse.  The coupling and attenuation of nearly resonant multiplets in the Earth's free oscillation spectrum , 1980 .

[66]  D. Komatitsch,et al.  Introduction to the spectral element method for three-dimensional seismic wave propagation , 1999 .

[67]  Anna L. Mazzucato,et al.  Partial Uniqueness and Obstruction to Uniqueness in Inverse Problems for Anisotropic Elastic Media , 2006 .

[68]  Ralph Abraham,et al.  Foundations Of Mechanics , 2019 .

[69]  J. Woodhouse.  On Rayleigh's Principle , 1976 .

[70]  Hideki Omori.  Infinite-Dimensional Lie Groups , 1996 .

[71]  On Postglacial Sea Level , 2007 .

[72]  B. Romanowicz.  Multiplet-multiplet coupling due to lateral heterogeneity: asymptotic effects on the amplitude and frequency of the Earth's normal modes , 1987 .

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

[74]  Normal mode multiplet coupling on an aspherical, anelastic earth , 1992 .