GPS Imaging of vertical land motion in California and Nevada: Implications for Sierra Nevada uplift

We introduce Global Positioning System (GPS) Imaging, a new technique for robust estimation of the vertical velocity field of the Earth's surface, and apply it to the Sierra Nevada Mountain range in the western United States. Starting with vertical position time series from Global Positioning System (GPS) stations, we first estimate vertical velocities using the MIDAS robust trend estimator, which is insensitive to undocumented steps, outliers, seasonality, and heteroscedasticity. Using the Delaunay triangulation of station locations, we then apply a weighted median spatial filter to remove velocity outliers and enhance signals common to multiple stations. Finally, we interpolate the data using weighted median estimation on a grid. The resulting velocity field is temporally and spatially robust and edges in the field remain sharp. Results from data spanning 5–20 years show that the Sierra Nevada is the most rapid and extensive uplift feature in the western United States, rising up to 2 mm/yr along most of the range. The uplift is juxtaposed against domains of subsidence attributable to groundwater withdrawal in California's Central Valley. The uplift boundary is consistently stationary, although uplift is faster over the 2011–2016 period of drought. Uplift patterns are consistent with groundwater extraction and concomitant elastic bedrock uplift, plus slower background tectonic uplift. A discontinuity in the velocity field across the southeastern edge of the Sierra Nevada reveals a contrast in lithospheric strength, suggesting a relationship between late Cenozoic uplift of the southern Sierra Nevada and evolution of the southern Walker Lane.

[1]  Robert W. King,et al.  Active tectonics of northwestern U.S. inferred from GPS‐derived surface velocities , 2013 .

[2]  T. Holzer.  Elastic expansion of the lithosphere caused by groundwater depletion , 1979 .

[3]  Jeffrey T. Freymueller,et al.  Seasonal and long-term vertical deformation in the Nepal Himalaya constrained by GPS and GRACE measurements , 2012 .

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

[5]  D. Sandwell,et al.  Vertical crustal displacement due to interseismic deformation along the San Andreas fault: Constraints from tide gauges , 2014 .

[6]  Z. Altamimi,et al.  ITRF2008: an improved solution of the international terrestrial reference frame , 2011 .

[7]  M. Lisowski,et al.  Newly discovered subsidence at Lassen Peak, southern Cascade Range, California, from InSAR and GPS , 2004 .

[8]  G. Matheron.  Principles of geostatistics , 1963 .

[9]  Kenneth D. Smith,et al.  Evidence for Moho‐lower crustal transition depth diking and rifting of the Sierra Nevada microplate , 2016 .

[10]  R. Weldon,et al.  Interseismic uplift rates for western Oregon and along‐strike variation in locking on the Cascadia subduction zone , 2009 .

[11]  J. Langbein.  Deformation of the Long Valley Caldera, California: inferences from measurements from 1988 to 2001 , 2003 .

[12]  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 .

[13]  Hans-Peter Plag,et al.  Contemporary uplift of the Sierra Nevada, western United States, from GPS and InSAR measurements , 2012 .

[14]  D. Shelly,et al.  Renewed inflation of Long Valley Caldera, California (2011 to 2014) , 2015 .

[15]  M. Marcos,et al.  Vertical land motion as a key to understanding sea level change and variability , 2016 .

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

[17]  C. Andronicos,et al.  Inference of the viscosity structure and mantle conditions beneath the Central Nevada Seismic Belt from combined postseismic and lake unloading studies , 2016 .

[18]  Yuri Fialko,et al.  Evidence of fluid-filled upper crust from observations of postseismic deformation due to the 1992 Mw7.3 Landers earthquake , 2004 .

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

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

[21]  Wilhelm Burger,et al.  Digital Image Processing - An Algorithmic Introduction using Java , 2008, Texts in Computer Science.

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

[23]  G. Blewitt,et al.  Assessing the impact of vertical land motion on twentieth century global mean sea level estimates , 2016 .

[24]  Olaf Berke.  Modified median polish kriging and its application to the Wolfcamp–Aquifer data , 2001 .

[25]  Zhong Lu,et al.  Magmatic activity beneath the quiescent Three Sisters volcanic center, central Oregon Cascade Range, USA , 2002 .

[26]  T. Herring,et al.  Distribution of present‐day vertical deformation across the Southern Alps, New Zealand, from 10 years of GPS data , 2010 .

[27]  J. W. Hendley,et al.  Eruptions of Lassen Peak, California, 1914 to 1917 , 1999 .

[28]  G. L. Farmer,et al.  Tectonics of Pliocene removal of lithosphere of the Sierra Nevada, California , 2004 .

[29]  G. Blewitt,et al.  Uplift and seismicity driven by groundwater depletion in central California , 2014, Nature.

[30]  S. Hreinsdóttir,et al.  Contemporary vertical velocity of the central Basin and Range and uplift of the southern Sierra Nevada , 2008 .

[31]  G. Blewitt,et al.  Geodetic constraints on contemporary deformation in the northern Walker Lane: 3. Central Nevada seismic belt postseismic relaxation , 2009 .

[32]  Susan Owen,et al.  Fault slip rates and interseismic deformation in the western Transverse Ranges, California , 2013 .

[33]  C. Conrad.  The solid Earth’s influence on sea level , 2013 .

[34]  Yosuke Aoki,et al.  Vertical deformation of the Japanese islands, 1996–1999 , 2003 .

[35]  R. Bürgmann,et al.  Steady subsidence of Medicine Lake volcano, northern California, revealed by repeated leveling surveys , 2002 .

[36]  Zhong Lu,et al.  Systematic assessment of atmospheric uncertainties for InSAR data at volcanic arcs using large-scale atmospheric models: Application to the Cascade volcanoes, United States , 2015 .

[37]  K. Furlong,et al.  Ephemeral crustal thickening at a triple junction:The Mendocino crustal conveyor , 1999 .

[38]  Zhen Liu,et al.  Constraints on time-dependent volcanic source models at Long Valley Caldera from 1996 to 2009 using InSAR and geodetic measurements , 2011 .

[39]  A. Newman,et al.  Constraints on continued episodic inflation at Long Valley Caldera, based on seismic and geodetic observations , 2009 .

[40]  Zhong Lu,et al.  Investigating long-term subsidence at Medicine Lake Volcano, CA, using multitemporal InSAR , 2014 .

[41]  William C. Hammond,et al.  Contemporary tectonic deformation of the Basin and Range province, western United States: 10 years of observation with the Global Positioning System , 2004 .

[42]  Geodetic observation of contemporary deformation in the northern Walker Lane: 1. Semipermanent GPS strategy , 2009 .

[43]  J. Lowenstern,et al.  Volcano Hazards Assessment for Medicine Lake Volcano, Northern California , 2007 .

[44]  W. R. Peltier,et al.  Postglacial variations in the level of the sea: Implications for climate dynamics and solid‐Earth geophysics , 1998 .

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

[46]  Y. Fialko,et al.  On the effects of thermally weakened ductile shear zones on postseismic deformation , 2013 .

[47]  George E. Hilley,et al.  Resolving vertical tectonics in the San Francisco Bay Area from permanent scatterer InSAR and GPS analysis , 2006 .

[48]  D. Hill.  Unrest in Long Valley Caldera, California, 1978–2004 , 2006, Geological Society, London, Special Publications.

[49]  J. Johansson,et al.  Space-Geodetic Constraints on Glacial Isostatic Adjustment in Fennoscandia , 2001, Science.

[50]  E. Humphreys,et al.  Forces acting on the Sierra Nevada block and implications for the strength of the San Andreas fault system and the dynamics of continental deformation in the western United States , 2008 .

[51]  Falk Amelung,et al.  Postseismic Mantle Relaxation in the Central Nevada Seismic Belt , 2005, Science.

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

[53]  Carl Tape,et al.  Multiscale estimation of GPS velocity fields , 2008 .

[54]  Jan P. Weiss,et al.  Single receiver phase ambiguity resolution with GPS data , 2010 .

[55]  E. Hauksson,et al.  Internal deformation of the southern Sierra Nevada microplate associated with foundering lower lithosphere, California , 2014 .

[56]  Meghan S. Miller,et al.  Present‐day motion of the Sierra Nevada block and some tectonic implications for the Basin and Range province, North American Cordillera , 2000 .

[57]  S. Wesnousky.  The San Andreas and Walker Lane fault systems, western North America: transpression, transtension, cumulative slip and the structural evolution of a major transform plate boundary , 2005 .

[58]  M. Tamisiea,et al.  The moving boundaries of sea level change: understanding the origins of geographic variability , 2011 .

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

[60]  S. Hreinsdóttir,et al.  Increasing long-wavelength relief across the southeastern flank of the Sierra Nevada, California , 2009 .

[61]  Yehuda Bock,et al.  Southern California permanent GPS geodetic array: Spatial filtering of daily positions for estimating coseismic and postseismic displacements induced by the 1992 Landers earthquake , 1997 .

[62]  D. Argus,et al.  Current Sierra Nevada-North America motion from very long baseline interferometry:Implications for the kinematics of the western United States , 1991 .

[63]  Georg Dresen,et al.  Rheology of the Lower Crust and Upper Mantle: Evidence from Rock Mechanics, Geodesy, and Field Observations , 2008 .

[64]  P. Cashman,et al.  Late Cenozoic structure and evolution of the Great Basin-Sierra Nevada transition , 2009 .

[65]  F. Pollitz,et al.  Dislocation models of interseismic deformation in the western United States , 2008 .

[66]  F. Pollitz,et al.  Temporal evolution of continental lithospheric strength in actively deforming regions , 2008 .

[67]  Christopher H. Scholz,et al.  Vertical deformation of the Japanese islands , 2002 .

[68]  G. Milne.  Space-geodetic constraints on glacial isostatic adjustment , 2001 .

[69]  R. Simpson,et al.  Northward migration of the Cascadia forearc in the northwestern U.S. and implications for subduction deformation , 2001 .

[70]  Duncan Carr Agnew,et al.  Ongoing drought-induced uplift in the western United States , 2014, Science.

[71]  D. Stockli,et al.  Evolution and Strain Reorganization within Late Neogene Structural Stepovers Linking the Central Walker Lane and Northern Eastern California Shear Zone, Western Great Basin , 2008 .

[72]  T. Herring,et al.  Far‐reaching transient motions after Mojave earthquakes require broad mantle flow beneath a strong crust , 2007 .

[73]  Felix W. Landerer,et al.  Seasonal variation in total water storage in California inferred from GPS observations of vertical land motion , 2014 .

[74]  H. Theil.  A Rank-Invariant Method of Linear and Polynomial Regression Analysis , 1992 .

[75]  Jeffrey T. Freymueller,et al.  GPS-derived strain in northwestern California: Termination of the San Andreas fault system and convergence of the Sierra Nevada–Great Valley block contribute to southern Cascadia forearc contraction , 2006 .

[76]  N. Cressie.  The origins of kriging , 1990 .

[77]  T. Dixon,et al.  Influence of the earthquake cycle and lithospheric rheology on the dynamics of the Eastern California Shear Zone , 2001 .

[78]  D G Krige,et al.  A statistical approach to some mine valuation and allied problems on the Witwatersrand , 2015 .

[79]  P. Sen.  Estimates of the Regression Coefficient Based on Kendall's Tau , 1968 .

[80]  D. Dzurisin,et al.  Crustal subsidence, seismicity, and structure near Medicine Lake Volcano, California , 1991 .

[81]  Mark Simons,et al.  Finite source modelling of magmatic unrest in Socorro, New Mexico, and Long Valley, California , 2001 .

[82]  T. Parsons,et al.  Vertical deformation associated with normal fault systems evolved over coseismic, postseismic, and multiseismic periods , 2016 .

[83]  Michael Bevis,et al.  Love's problem , 2004 .

[84]  Robert W. King,et al.  Fault locking, block rotation and crustal deformation in the Pacific Northwest , 2007 .

[85]  R. Rau,et al.  Present-day kinematics of active mountain building in Taiwan from GPS observations during 1995-2005 , 2011 .

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

[87]  W. Holt,et al.  Dynamics of the pacific-north american plate boundary in the western united states , 2000, Science.

[88]  Zhong Lu,et al.  Time-scale and mechanism of subsidence at Lassen Volcanic Center, CA, from InSAR , 2016 .

[89]  Michael P. Poland,et al.  Constraints on the mechanism of long-term, steady subsidence at Medicine Lake volcano, northern California, from GPS, leveling, and InSAR , 2006 .

[90]  B. Hager,et al.  Postseismic relaxation across the Central Nevada Seismic Belt , 2003 .

[91]  A. Grayver,et al.  Journal of Geophysical Research : Solid Earth Stochastic Inversion of Geomagnetic Observatory Data Including Rigorous Treatment of the Ocean Induction Effect With Implications for Transition ZoneWater Content and Thermal Structure , 2018 .

[92]  Felix W. Landerer,et al.  GPS as an independent measurement to estimate terrestrial water storage variations in Washington and Oregon , 2015 .

[93]  C. Busby.  Birth of a plate boundary at ca. 12 Ma in the Ancestral Cascades arc, Walker Lane belt of California and Nevada , 2013 .