Optimal interpolation of spatially discretized geodetic data

Abstract We present an algorithm to calculate horizontal strains (or strain rates) through interpolation of geodetically derived displacements (or velocities). This is an underdetermined inverse problem to derive smoothly distributed strains (or strain rates) using spatially discretized geodetic observations. A priori information, in the form of weighted smoothing, is critical to facilitate the solution. At a given site, the horizontal displacement (or velocity) field in its vicinity is approximated by a bilinear function and represented by rigid block translation, rotation, and strains (or their rates). The weighted displacement (or velocity) data in the neighborhood are used to estimate the field parameters through a least‐squares inversion procedure. Optimal weightings are prescribed for the neighboring data, based on their distances to the interpolation site and their spatial coverage. Nonelastic strains resulted from surface fault rupture and creep may also be excluded from the solution. We apply this method to the Southern California Earthquake Center Crustal Motion Map version 4.0 velocity field and derive the strain‐rate field in southern California. Our result shows that (1) distance‐dependent weighting can be optimally achieved by employing either a Gaussian or quadratic decay function, with the former offering a slightly sharper result than the latter. (2) Spatially dependent weighting is important to improve the interpolation, and can be done by invoking either an azimuthal weighting or a Voronoi cell areal weighting function. (3) The strain‐rate pattern in southern California is dominated by dextral shear of the San Andreas fault (SAF) system, and the secondary faults surrounding the Big Bend of the SAF strike at oblique angles with respect to the maximum shear direction, suggesting that tectonic deformation field on and off the SAF is dominated by mechanic processes of the SAF.

[1]  Akemi Noda,et al.  Physics-based GPS data inversion to estimate three-dimensional elastic and inelastic strain fields , 2010 .

[2]  D. Sandwell.  BIHARMONIC SPLINE INTERPOLATION OF GEOS-3 AND SEASAT ALTIMETER DATA , 1987 .

[3]  A. Allam,et al.  Seismic velocity structures in the southern California plate-boundary environment from double-difference tomography , 2011 .

[4]  G. Matheron.  Random Functions and their Application in Geology , 1970 .

[5]  X. Pichon,et al.  Velocity field in Asia inferred from Quaternary fault slip rates and Global Positioning System observations , 2000 .

[6]  Interpolation by Regularized Spline with Tension � , 2022 .

[7]  Mark Simons,et al.  Interseismic crustal deformation in the Taiwan plate boundary zone revealed by GPS observations, seismicity, and earthquake focal mechanisms , 2008 .

[8]  W. Holt,et al.  A procedure for obtaining the complete horizontal motions within zones of distributed deformation from the inversion of strain rate data , 1993 .

[9]  W. Holt,et al.  Deformation kinematics in the western United States determined from Quaternary fault slip rates and recent geodetic data , 1999 .

[10]  Steven Fortune,et al.  Voronoi Diagrams and Delaunay Triangulations , 2004, Handbook of Discrete and Computational Geometry, 2nd Ed..

[11]  A. John Haines,et al.  An integrated global model of present‐day plate motions and plate boundary deformation , 2003 .

[12]  Paul Wessel,et al.  Interpolation using a generalized Green's function for a spherical surface spline in tension , 2008 .

[13]  P. Gazis,et al.  A large‐scale survey of corotating interaction regions and their successors in the outer heliosphere , 2000 .

[14]  T. Kato,et al.  Continuous distribution of the horizontal strain in the Tohoku district, Japan, predicted by least-squares collocation , 1998 .

[15]  Geoffrey Blewitt,et al.  A geodetic plate motion and Global Strain Rate Model , 2014 .

[16]  Piyush Agram,et al.  Aseismic slip and seismogenic coupling along the central San Andreas Fault , 2015 .

[17]  William H. Prescott,et al.  An extension of Frank's method for obtaining crustal shear strains from survey data , 1976 .

[18]  T. Kato,et al.  Distribution of vertical crustal movement rates in the Tohoku district, Japan, predicted by least-squares collocation , 1997 .

[19]  Paul Wessel,et al.  Interpolation with Splines in Tension: A Green's Function Approach , 1998 .

[20]  S. Ward.  On the consistency of earthquake moment rates, geological fault data, and space geodetic strain: the United States , 1998 .

[21]  F. C. Frank,et al.  Deduction of earth strains from survey data , 1966 .

[22]  Yehuda Bock,et al.  Creep along the Imperial Fault, southern California, from GPS measurements , 2002 .

[23]  D. Jackson,et al.  Southern California Deformation , 1997, Science.

[24]  H. Mitásová,et al.  Interpolation by regularized spline with tension: I. Theory and implementation , 1993 .

[25]  Reinhard Dietrich,et al.  Crustal deformation across the Southern Patagonian Icefield observed by GNSS , 2016 .

[26]  Kenneth W. Hudnut,et al.  Crustal deformation measured in Southern California , 1997 .

[27]  Duncan Carr Agnew,et al.  A unified analysis of crustal motion in Southern California, 1970–2004: The SCEC crustal motion map , 2011 .

[28]  Michael Edward Hohn,et al.  An Introduction to Applied Geostatistics: by Edward H. Isaaks and R. Mohan Srivastava, 1989, Oxford University Press, New York, 561 p., ISBN 0-19-505012-6, ISBN 0-19-505013-4 (paperback), $55.00 cloth, $35.00 paper (US) , 1991 .

[29]  Y. Kagan,et al.  Implications of Geodetic Strain Rate for Future Earthquakes, with a Five-Year Forecast of M5 Earthquakes in Southern California , 2007 .

[30]  Walter H. F. Smith,et al.  New, improved version of generic mapping tools released , 1998 .

[31]  David D. Jackson,et al.  Crustal deformation across and beyond the Los Angeles basin from geodetic measurements , 1996 .

[32]  G. Blewitt,et al.  A Geodetic Strain Rate Model for the Pacific-North American Plate Boundary, western United States , 2012 .

[33]  R. Coleman,et al.  A comparison of computation methods for crustal strains from geodetic measurements , 1981 .