Fourier-domain Green's function for an elastic semi-infinite solid under gravity, with applications to earthquake and volcano deformation

SUMMARY We present an analytic solution in the Fourier domain for an elastic deformation in a semi-infinite solid due to an arbitrary surface traction. We generalize the so-called Boussinesq's and Cerruti's problems to include a restoring buoyancy boundary condition at the surface. Buoyancy due to a large density contrast at the Earth's surface is an approximation to the full effect of gravity that neglects the perturbation of the gravitational potential and the change in density in the interior. Using the perturbation method, and assuming that the effect of gravity is small compared to the elastic deformation, we derive an approximation in the space domain to the Boussinesq's problem that accounts for a buoyancy boundary condition at the surface. The Fourier- and space-domain solutions are shown to be in good agreement. Numerous problems of elastostatic or quasi-static time-dependent deformation relevant to faulting in the Earth's interior (including inelastic deformation) can be modelled using equivalent body forces and surface tractions. Solving the governing equations with the elastic Green's function in the space domain can be impractical as the body force can be distributed over a large volume. We present a computationally efficient method to evaluate the elastic deformation in a 3-D half space due to the presence of an arbitrary distribution of internal forces and tractions at the surface of the half space. We first evaluate the elastic deformation in a periodic Cartesian volume in the Fourier domain, then use the analytic solutions to the generalized Boussinesq's and Cerruti's problems to satisfy the prescribed mixed boundary condition at the surface. We show some applications for magmatic intrusions and faulting. This approach can be used to solve elastostatic problems involving spatially heterogeneous elastic properties (by employing a homogenization method) and time-dependent problems such as non-linear viscoelastic relaxation, poroelastic rebound and non-steady fault creep under the assumption of spatially homogeneous elastic properties.

[1]  Sylvain Barbot,et al.  Three-dimensional models of elastostatic deformation in heterogeneous media, with applications to the Eastern California Shear Zone , 2009 .

[2]  Y. Fialko.  Probing the mechanical properties of seismically active crust with space geodesy: Study of the coseismic deformation due to the 1992 Mw7.3 Landers (southern California) earthquake , 2004 .

[3]  Paul Segall,et al.  Imaging of aseismic fault slip transients recorded by dense geodetic networks , 2003 .

[4]  Paul Segall,et al.  A transient subduction zone slip episode in southwest Japan observed by the nationwide GPS array , 2003 .

[5]  Yuri Fialko,et al.  Interseismic strain accumulation and the earthquake potential on the southern San Andreas fault system , 2006, Nature.

[6]  Rongjiang Wang,et al.  PSGRN/PSCMP - a new code for calculating co- and post-seismic deformation, geoid and gravity changes based on the viscoelastic-gravitational dislocation theory , 2006, Comput. Geosci..

[7]  Fred F. Pollitz,et al.  Gravitational viscoelastic postseismic relaxation on a layered spherical Earth , 1997 .

[8]  Y. Okada.  Internal deformation due to shear and tensile faults in a half-space , 1992, Bulletin of the Seismological Society of America.

[9]  Sylvain Barbot,et al.  Seismic and geodetic evidence for extensive, long-lived fault damage zones , 2009 .

[10]  K. Mogi.  Relations between the Eruptions of Various Volcanoes and the Deformations of the Ground Surfaces around them , 1958 .

[11]  A. Love.  A treatise on the mathematical theory of elasticity , 1892 .

[12]  Mark Simons,et al.  Deformation due to a pressurized horizontal circular crack in an elastic half-space, with applications to volcano geodesy , 2001 .

[13]  WangRongjiang,et al.  Computation of deformation induced by earthquakes in a multi-layered elastic crust , 2003 .

[14]  Yehuda Bock,et al.  Parkfield earthquake: Stress-driven creep on a fault with spatially variable rate-and-state friction parameters , 2009 .

[15]  D. Sandwell,et al.  A three-dimensional semianalytic viscoelastic model for time-dependent analyses of the earthquake cycle , 2004 .

[16]  Raymond D. Mindlin,et al.  Thermoelastic Stress in the Semi-Infinite Solid , 1950 .

[17]  Rongjiang Wang,et al.  Erratum to: "Computation of deformation induced by earthquakes in a multi-layered elastic crust - FORTRAN programs EDGRN/EDCMP": [Computers & Geosciences, 29(2) (2003) 195-207] , 2006, Comput. Geosci..

[18]  Y. Fialko.  Fracture and Frictional Mechanics: Theory , 2007 .

[19]  Huajian Gao,et al.  Quasi‐static dislocations in three dimensional inhomogeneous media , 1997 .

[20]  Brendan J. Meade,et al.  Algorithms for the calculation of exact displacements, strains, and stresses for triangular dislocation elements in a uniform elastic half space , 2007, Comput. Geosci..

[21]  R. D. Mindlin.  Note on the Galerkin and Papkovitch stress functions , 1936 .

[22]  G. Bonnet,et al.  Green’s operator for a periodic medium with traction-free boundary conditions and computation of the effective properties of thin plates , 2008 .

[23]  Leon Knopoff,et al.  Body Force Equivalents for Seismic Dislocations , 1964 .

[24]  S. Nemat-Nasser,et al.  Micromechanics: Overall Properties of Heterogeneous Materials , 1993 .

[25]  Yehuda Bock,et al.  Frictional Afterslip Following the 2005 Nias-Simeulue Earthquake, Sumatra , 2006, Science.

[26]  D. Agnew,et al.  The complete (3‐D) surface displacement field in the epicentral area of the 1999 MW7.1 Hector Mine Earthquake, California, from space geodetic observations , 2001 .

[27]  Roland Bürgmann,et al.  Evidence of power-law flow in the Mojave desert mantle , 2004, Nature.

[28]  D. Sandwell,et al.  Effect of a compliant fault zone on the inferred earthquake slip distribution , 2008 .

[29]  Y. Fialko,et al.  Mechanics of active magmatic intraplating in the Rio Grande Rift near Socorro, New Mexico , 2010 .

[30]  D. Wolf.  Viscoelastodynamics of a stratified, compressible planet: incremental field equations and short- and long-time asymptotes , 1991 .

[31]  Sylvain Barbot,et al.  A unified continuum representation of post-seismic relaxation mechanisms: semi-analytic models of afterslip, poroelastic rebound and viscoelastic flow , 2010 .

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

[33]  H. M. Westergaard.  General solution of the problem of elastostatics of an $n$-dimensional homogeneous isotropic solid in an $n$-dimensional space , 1935 .

[34]  I. Okumura.  ON THE GENERALIZATION OF CERRUTI'S PROBLEM IN AN ELASTIC HALF-SPACE , 1995 .

[35]  Sylvain Barbot,et al.  Space geodetic investigation of the coseismic and postseismic deformation due to the 2003 Mw7.2 Altai earthquake: Implications for the local lithospheric rheology , 2008 .

[36]  John B. Rundle,et al.  Viscoelastic‐gravitational deformation by a rectangular thrust fault in a layered Earth , 1982 .

[37]  J. Steketee.  ON VOLTERRA'S DISLOCATIONS IN A SEMI-INFINITE ELASTIC MEDIUM , 1958 .

[38]  S. Orszag,et al.  Advanced Mathematical Methods For Scientists And Engineers , 1979 .

[39]  Paul Johnston,et al.  Dependence of horizontal stress magnitude on load dimension in glacial rebound models , 2002 .

[40]  James H. Dieterich,et al.  Deformation from Inflation of a Dipping Finite Prolate Spheroid in an Elastic Half‐Space as a Model for Volcanic Stressing , 1988 .

[41]  Cnlumhut Kutitrrattn,et al.  FORCE AT A POINT IN THE INTERIOR OF A SEMI-INFINITE SOLID , 2009 .

[42]  Raymond D. Mindlin,et al.  Nuclei of Strain in the Semi‐Infinite Solid , 1950 .

[43]  L. Rivera,et al.  Coseismic Deformation from the 1999 Mw 7.1 Hector Mine, California, Earthquake as Inferred from InSAR and GPS Observations , 2002 .

[44]  F. Pollitz,et al.  Mantle Flow Beneath a Continental Strike-Slip Fault: Postseismic Deformation After the 1999 Hector Mine Earthquake , 2001, Science.