A Robust Estimation of the 3‐D Intraplate Deformation of the North American Plate From GPS

Glacial isostatic adjustment (GIA) is the main cause of deformation in intraplate North America. Here we use up to 3,271 Global Positioning System station velocities to image this 3-D deformation across the entire plate. We apply a new robust strain rate estimation algorithm (median estimation of local deformation), which does not require the assumption that part of the plate is unaffected by GIA, an assumption we show to be false. Our results show extension in the area underneath the Laurentide ice sheet, contrasted with a semiannular belt of horizontal contraction of up to ~4 × 10 9 yr 1 around the former ice sheet. This contractional belt is kinematically linked to an ~1–2 mm yr 1 far-field horizontal motion directed toward the ice sheet. Our results, together with a new robustly imaged vertical velocity field, are consistent with GIA as themain cause of deformation, although the contractional strain rates around the former ice sheet and the far-field horizontal velocities are significantly higher than those predicted by the ICE6G_C(VM5a) model. This finding suggests that our results, including the location of reversal in sign of horizontal motion relative to the ice sheet, will be useful to reevaluate the mantle viscosity structure used by GIA models. Besides the GIA-attributed deformation, we find almost no other region with significant strain accumulation. A plate-scale spatial correlation between strain rate and seismicity is absent, suggesting that GIA is a limited driver of contemporary seismogenesis and that intraplate seismicity must be attributable to factors other than secular strain accumulation. Plain Language Summary The theory of plate tectonics says that tectonic plates move as rigid blocks along the Earth’s surface and that the Earth’s crust should only deform at the boundary between plates. However, the recent explosion in the number of high-precision Global Positioning System stations allowed us to capture some subtle deformation patterns inside the North American plate that only became apparent by a very careful analysis of the relative motions between thousands of stations. We found that most of the plate is moving at 1–2 mm/year towards central Canada. Consequently, around most of Canada there is a zone where the crust is contracting. Within Canada, the crust is extending outward and is moving upward rapidly. These patterns can be explained by the process of the crust and mantle still rebounding from a time when it was covered by a thick ice sheet about 16,000 years ago. The fact that this causes the land to move towards the former ice sheet is an unexpected result that will be useful in understanding the relaxation properties of the underlying mantle. Moreover, we found that earthquakes inside the North American plate do not occur where we see the crust deform, which leaves these events still enigmatic.

[1]  Z. Altamimi,et al.  ITRF2014 plate motion model , 2017 .

[2]  B. Hager,et al.  Geoid Anomalies in a Dynamic Earth , 1984 .

[3]  W. Peltier,et al.  Constraining Models of Postglacial Rebound Using Space Geodesy , 2008 .

[4]  Angelyn W. Moore,et al.  The Antarctica component of postglacial rebound model ICE-6G_C (VM5a) based on GPS positioning, exposure age dating of ice thicknesses, and relative sea level histories , 2014 .

[5]  J. Rupp,et al.  Preliminary Results from a GPS Geodetic Network in the Southern Illinois Basin , 2002 .

[6]  Peter Rousseeuw,et al.  Computing location depth and regression depth in higher dimensions , 1998, Stat. Comput..

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

[8]  Peter Molnar,et al.  Earthquake recurrence intervals and plate tectonics , 1979 .

[9]  K. Tiampo,et al.  GPS coordinate time series measurements in Ontario and Quebec, Canada , 2017, Journal of Geodesy.

[10]  Guoquan Wang,et al.  GPS-derived ground deformation (2005–2014) within the Gulf of Mexico region referred to a stable Gulf of Mexico reference frame , 2015 .

[11]  Jian Wang,et al.  Bedrock displacements in Greenland manifest ice mass variations, climate cycles and climate change , 2012, Proceedings of the National Academy of Sciences.

[12]  Xavier Collilieux,et al.  ITRF2008 plate motion model , 2011 .

[13]  T. Dixon,et al.  Accelerating uplift in the North Atlantic region as an indicator of ice loss , 2010 .

[14]  E. Larour,et al.  Mass transport waves amplified by intense Greenland melt and detected in solid Earth deformation , 2017 .

[15]  Andrew V. Newman,et al.  Should Memphis build for California's earthquakes? , 2003 .

[16]  Richard B. Langley,et al.  Canadian High Arctic Ionospheric Network (CHAIN) , 2009 .

[17]  S. Stein,et al.  2000 years of migrating earthquakes in North China: How earthquakes in midcontinents differ from those at plate boundaries , 2011 .

[18]  James L. Davis,et al.  Investigation of glacial isostatic adjustment in the northeast U.S. using GPS measurements , 2002 .

[19]  G. Blewitt,et al.  MIDAS robust trend estimator for accurate GPS station velocities without step detection , 2016, Journal of geophysical research. Solid earth.

[20]  S. F. Obermeier,et al.  Geologic Evidence for Recurrent Moderate to Large Earthquakes Near Charleston, South Carolina , 1985, Science.

[21]  T. Dixon,et al.  Subsidence along the Atlantic Coast of North America: Insights from GPS and late Holocene relative sea level data , 2016 .

[22]  J. Nocquet,et al.  Deformation of the North American plate interior from a decade of continuous GPS measurements , 2006 .

[23]  J. Nocquet,et al.  Seismology: Tectonic strain in plate interiors? , 2005, Nature.

[24]  Patrick Wu,et al.  Effects of lateral variations in lithospheric thickness and mantle viscosity on glacially induced surface motion in Laurentia , 2005 .

[25]  M. Hamburger,et al.  Geodetic Observations of Active Intraplate Crustal Deformation in the Wabash Valley Seismic Zone and the Southern Illinois Basin , 2010 .

[26]  R. Snay,et al.  Modeling 3‐D crustal velocities in the United States and Canada , 2016 .

[27]  K. Lambeck,et al.  Viscosity structure of Earth's mantle inferred from rotational variations due to GIA process and recent melting events , 2015 .

[28]  Newman,et al.  Slow deformation and lower seismic hazard at the new madrid seismic zone , 1999, Science.

[29]  F. Tupin,et al.  Time series analysis of Mexico City subsidence constrained by radar interferometry , 2009 .

[30]  W. J. Morgan,et al.  Horizontal motions due to post‐glacial rebound , 1990 .

[31]  S. Stein,et al.  Mid-continental earthquakes: Spatiotemporal occurrences, causes, and hazards , 2016 .

[32]  P. Basham,et al.  Seismicity and Seismotectonics of Canada’s Eastern Margin and Craton , 1989 .

[33]  M. Tamisiea,et al.  GRACE Gravity Data Constrain Ancient Ice Geometries and Continental Dynamics over Laurentia , 2007, Science.

[34]  J. Mitrovica,et al.  Radial profile of mantle viscosity: Results from the joint inversion of convection and postglacial , 1997 .

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

[36]  C. Kreemer,et al.  Pacific plate deformation from horizontal thermal contraction , 2014 .

[37]  D. Rudolph,et al.  Analysis of long‐term land subsidence near Mexico City: Field investigations and predictive modeling , 1999 .

[38]  J. Wahr,et al.  The use of GPS horizontals for loading studies, with applications to northern California and southeast Greenland , 2013 .

[39]  M. A. Goudarzi,et al.  Present-Day 3D Velocity Field of Eastern North America Based on Continuous GPS Observations , 2016, Pure and Applied Geophysics.

[40]  D. Argus,et al.  Tests of the rigid-plate hypothesis and bounds on intraplate deformation using geodetic data from very long baseline interferometry , 1996 .

[41]  Enzo Boschi,et al.  Glacial isostasy, sea-level and mantle rheology , 1991 .

[42]  A CRUSTAL DEFORMATION STUDY OF THE CHARLEVOIX SEISMIC ZONE IN QUEBEC , 2010 .

[43]  Z. Martinec,et al.  Glacial isostasy and plate motion , 2008 .

[44]  Ilya Zaliapin,et al.  Clustering analysis of seismicity and aftershock identification. , 2007, Physical review letters.

[45]  A. Crone,et al.  Style and Timing of Holocene Surface Faulting on the Meers Fault, Southwestern Oklahoma , 1990 .

[46]  Archie Paulson,et al.  FAST TRACK PAPER: Inference of mantle viscosity from GRACE and relative sea level data , 2007 .

[47]  Geoffrey Blewitt,et al.  GPS Imaging of vertical land motion in California and Nevada: Implications for Sierra Nevada uplift , 2016, Journal of geophysical research. Solid earth.

[48]  J. Tromp,et al.  Influence of lithospheric thickness variations on 3‐D crustal velocities due to glacial isostatic adjustment , 2005 .

[49]  A. Frankel,et al.  Significant Motions between GPS Sites in the New Madrid Region: Implications for Seismic Hazard , 2012 .

[50]  M. Bevis,et al.  Vertical and horizontal surface displacements near Jakobshavn Isbræ driven by melt‐induced and dynamic ice loss , 2012 .

[51]  T. James,et al.  Absolute gravity calibration of GPS velocities and glacial isostatic adjustment in mid‐continent North America , 2011 .

[52]  Charles S. Mueller,et al.  2014 Update of the United States National Seismic Hazard Maps , 2012 .

[53]  B. Hager,et al.  Inversion for mantle viscosity profiles constrained by dynamic topography and the geoid, and their estimated errors , 2000 .

[54]  Arthur Frankel Comment on Should Memphis Build for California's Earthquakes? , 2003 .

[55]  Hansheng Wang,et al.  Optimal locations for GPS measurements in North America and northern Europe for constraining Glacial Isostatic Adjustment , 2010 .

[56]  Isabella Velicogna,et al.  Regional acceleration in ice mass loss from Greenland and Antarctica using GRACE time‐variable gravity data , 2014 .

[57]  J. Gomberg,et al.  Comment: Caution urged in revising earthquake hazard estimates in New Madrid Seismic Zone , 1999 .

[58]  I. Shapiro,et al.  A spectral formalism for computing three‐dimensional deformations due to surface loads: 2. Present‐day glacial isostatic adjustment , 1994 .

[59]  R. Wheeler,et al.  A Seismotectonic Model for the 300-Kilometer-Long Eastern Tennessee Seismic Zone , 1994, Science.

[60]  Arthur L. Lerner-Lam,et al.  Geodetic constraints on the rigidity and relative motion of Eurasia and North America , 2000 .

[61]  S. Owen,et al.  The angular velocities of the plates and the velocity of Earth's centre from space geodesy , 2010 .

[62]  D. Yuen,et al.  Deep Continental Roots: The Effects of Lateral Variations of Viscosity on Post-Glacial Rebound , 1991 .

[63]  C. Kreemer,et al.  The Global Moment Rate Distribution within Plate Boundary Zones , 2013 .

[64]  L. Bollinger,et al.  Evidence for the release of long‐term tectonic strain stored in continental interiors through intraplate earthquakes , 2016 .

[65]  Timothy H. Dixon,et al.  Space geodetic imaging of rapid ground subsidence in Mexico City , 2008 .

[66]  P. Rousseeuw,et al.  High-dimensional computation of the deepest location , 2000 .

[67]  M. Zoback,et al.  Intraplate earthquakes, regional stress and fault mechanics in the Central and Eastern U.S. and Southeastern Canada , 2012 .

[68]  K. Kelson,et al.  Multiple late Holocene earthquakes along the Reelfoot fault, central New Madrid seismic zone , 1996 .

[69]  John Townend,et al.  State of stress in central and eastern North American seismic zones , 2010 .

[70]  J. Wahr,et al.  Elastic uplift in southeast Greenland due to rapid ice mass loss , 2007 .

[71]  Jean-Mathieu Nocquet,et al.  Geodetic constraints on glacial isostatic adjustment in Europe , 2005 .

[72]  W. Peltier,et al.  Space geodesy constrains ice age terminal deglaciation: The global ICE‐6G_C (VM5a) model , 2015 .

[73]  G. Spada,et al.  POSTGLACIAL REBOUND AND LATERAL VISCOSITY VARIATIONS : A SEMI-ANALYTICAL APPROACH BASED ON A SPHERICAL MODEL WITH MAXWELL RHEOLOGY , 1997 .

[74]  Michael R. Craymer,et al.  Observation of glacial isostatic adjustment in “stable” North America with GPS , 2007 .

[75]  Guoquan Wang,et al.  Introduction to GPS geodetic infrastructure for land subsidence monitoring in Houston, Texas, USA , 2015 .

[76]  Ian J. Ferguson,et al.  Investigating Canada's Lithosphere and earthquake hazards with portable arrays , 2005 .

[77]  Wei Jiang,et al.  Glacial rebound and plate spreading: Results from the first countrywide GPS observations in Iceland , 2009 .

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

[79]  Boudewijn Ambrosius,et al.  Recent Geodetic Results in the Azores Triple Junction Region , 2004 .

[80]  Steven N. Ward A multidisciplinary approach to seismic hazard in southern California , 1994 .

[81]  P. Rousseeuw,et al.  Alternatives to the Median Absolute Deviation , 1993 .

[82]  I. Shapiro,et al.  Constraining proposed combinations of ice history and earth rheology using VLBI determined baseline length rates in North America , 1993 .

[83]  Soroosh Sorooshian,et al.  SuomiNet: A Real-Time National GPS Network for Atmospheric Research and Education. , 2000 .

[84]  S. Stein,et al.  Estimation of intraplate strain accumulation in the New Madrid seismic zone from repeat GPS surveys , 1998 .

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

[86]  Longwei Xiang,et al.  Increased water storage in North America and Scandinavia from GRACE gravity data , 2012, Nature Geoscience.

[87]  C. A. R. Hoare,et al.  Algorithm 64: Quicksort , 1961, Commun. ACM.

[88]  Guoquan Wang,et al.  Current Land Subsidence and Groundwater Level Changes in the Houston Metropolitan Area (2005–2012) , 2015 .

[89]  Michael G. Sideris,et al.  Use of GRACE determined secular gravity rates for glacial isostatic adjustment studies in North-America , 2008 .

[90]  C. Kreemer,et al.  A comparison of strain rates and seismicity for Fennoscandia: depth dependency of deformation from glacial isostatic adjustment , 2015 .

[91]  R. Cox,et al.  Paleoseismic Evidence for Multiple Mw≥6 Earthquakes in the Eastern Tennessee Seismic Zone during the Late Quaternary , 2017 .

[92]  S. Stein,et al.  Mid-Continent Earthquakes as a Complex System , 2009 .

[93]  W. Peltier,et al.  Rheological stratification of the lithosphere: A direct inference based upon the geodetically observed pattern of the glacial isostatic adjustment of the North American continent , 2008 .

[94]  T. Dixon,et al.  Focused study of interweaving hazards across the Caribbean , 2012 .

[95]  Guy Masters,et al.  An inversion for radial viscosity structure using seismic tomography , 1992 .

[96]  É. Calais,et al.  Strain accumulation in the New Madrid and Wabash Valley seismic zones from 14 years of continuous GPS observation , 2014 .

[97]  Timothy H. Dixon,et al.  Tectonic control of subsidence and southward displacement of southeast Louisiana with respect to stable North America , 2006 .

[98]  A. Barr,et al.  Measuring water accumulation rates using GRACE data in areas experiencing glacial isostatic adjustment: The Nelson River basin , 2013 .

[99]  S. Stein,et al.  A new paradigm for large earthquakes in stable continental plate interiors , 2016 .

[100]  W. Ellsworth Injection-Induced Earthquakes , 2013, Science.

[101]  Ingo Sasgen,et al.  Towards the inversion of GRACE gravity fields for present-day ice-mass changes and glacial-isostatic adjustment in North America and Greenland , 2012 .

[102]  Jeffrey T. Freymueller,et al.  A viscoelastic and afterslip postseismic deformation model for the 1964 Alaska earthquake , 2009 .

[103]  Michael Bevis,et al.  On computing the geoelastic response to a disk load , 2016 .

[104]  Geoffrey Blewitt,et al.  Terrestrial reference frame NA12 for crustal deformation studies in North America , 2013 .

[105]  M. Ellis,et al.  Space geodetic evidence for rapid strain rates in the New Madrid seismic zone of central USA , 2005, Nature.

[106]  J. Mitrovica,et al.  A new inference of mantle viscosity based upon joint inversion of convection and glacial isostatic adjustment data , 2004 .

[107]  A. Moore,et al.  Geodynamics and tectonics The Antarctica component of postglacial rebound model ICE-6 G _ C ( VM 5 a ) based on GPS positioning , exposure age dating of ice thicknesses , and relative sea level histories , 2014 .

[108]  T. James,et al.  GPS crustal strain, postglacial rebound, and seismic hazard in eastern North America: The Saint Lawrence valley example , 2005 .

[109]  S. Stein,et al.  How rigid is the stable interior of the North American Plate , 1996 .

[110]  Eugene S. Schweig,et al.  THE -ENIGMA OF THE NEW MADRID EARTHQUAKES OF 1811-1812 , 1996 .

[111]  K. Tiampo,et al.  Analysis of GPS Measurements in Eastern Canada Using Principal Component Analysis , 2009, Pure and Applied Geophysics.

[112]  W. Holt,et al.  Earthquake strain rates and instantaneous relative motions within central and eastern Asia , 1995 .

[113]  J. Zumberge,et al.  Precise point positioning for the efficient and robust analysis of GPS data from large networks , 1997 .

[114]  K. Simon,et al.  A data-driven model for constraint of present-day glacial isostatic adjustment in North America , 2017 .

[115]  Robert H. Lafferty,et al.  The Earthquake Potential of the New Madrid Seismic Zone , 2002 .

[116]  M. Zoback,et al.  Rapid Intraplate Strain Accumulation in the New Madrid Seismic Zone , 1992, Science.

[117]  G. Spada,et al.  Lateral viscosity variations and post‐glacial rebound: Effects on present‐day VLBI baseline deformations , 1997 .

[118]  Patrick Wu,et al.  Water storage changes in North America retrieved from GRACE gravity and GPS data , 2015 .

[119]  S. Hough,et al.  The New Madrid Seismic Zone: Not Dead Yet , 2014, Science.

[120]  Geoffrey Blewitt,et al.  Effect of annual signals on geodetic velocity , 2002 .

[121]  Ganapati Panda,et al.  Identification of Glacial Isostatic Adjustment in Eastern Canada Using S Transform Filtering of GPS Observations , 2011, Pure and Applied Geophysics.

[122]  P. Talwani,et al.  Paleoseismic Evidence for Recurrence of Earthquakes near Charleston, South Carolina , 1985, Science.

[123]  S. Zhao Lithosphere thickness and mantle viscosity estimated from joint inversion of GPS and GRACE-derived radial deformation and gravity rates in North America , 2013 .

[124]  R. Madole Stratigraphic evidence of Holocene faulting in the mid-continent: The Meers fault, southwestern Oklahoma , 1988 .

[125]  Z. Altamimi,et al.  ITRF2005 : A new release of the International Terrestrial Reference Frame based on time series of station positions and Earth Orientation Parameters , 2007 .

[126]  Yan Jiang,et al.  City subsidence observed with persistent scatterer InSAR , 2010 .

[127]  C. A. R. Hoare,et al.  Algorithm 65: find , 1961, Commun. ACM.

[128]  Richard D. Miller,et al.  Shallow seismic reflection profile of the Meers Fault, Comanche County, Oklahoma , 1987 .

[129]  H. Dragert,et al.  Crustal motion and deformation monitoring of the canadian landmass , 2019 .

[130]  Robert J. Renka,et al.  Algorithm 772: STRIPACK: Delaunay triangulation and Voronoi diagram on the surface of a sphere , 1997, TOMS.

[131]  Ronald G. Blom,et al.  Post‐glacial sediment load and subsidence in coastal Louisiana , 2007 .

[132]  W. Prescott,et al.  Crustal deformation rates in central and eastern U.S. inferred from GPS , 2001 .