Seasonal variation in total water storage in California inferred from GPS observations of vertical land motion

GPS is accurately recording vertical motion of Earth's surface in elastic response to seasonal changes in surface water storage in California. California's mountains subside up to 12 mm in the fall and winter due to the load of snow and rain and then rise an identical amount in the spring and summer when the snow melts, the rain runs off, and soil moisture evaporates. We invert the GPS observations of seasonal vertical motions to infer changes in equivalent water thickness. GPS resolves the distribution of change in total water across California's physiographic provinces at a resolution of 50 km, compared to 200 km resolution from the Gravity Recovery and Climate Experiment. The seasonal surface water thickness change is 0.6 m in the Sierra Nevada, Klamath, and southern Cascade Mountains and decreases sharply to about 0.1 m east into the Great Basin and west toward the Pacific coast. GPS provides an independent inference of change in total surface water, indicating water storage to be on average 50% larger than in the NLDAS‐Noah hydrology model, likely due to larger changes in snow and reservoir water than in the model.

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

[2]  M. Heflin,et al.  Horizontal motion in elastic response to seasonal loading of rain water in the Amazon Basin and monsoon water in Southeast Asia observed by GPS and inferred from GRACE , 2013 .

[3]  Jeffrey T. Freymueller,et al.  Repeated large Slow Slip Events at the southcentral Alaska subduction zone , 2013 .

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

[5]  M. Cheng,et al.  Deceleration in the Earth's oblateness , 2013 .

[6]  L. Jia,et al.  Load Love numbers and Green's functions for elastic Earth models PREM, iasp91, ak135, and modified models with refined crustal structure from Crust 2.0 , 2012, Comput. Geosci..

[7]  D. Argus.  Uncertainty in the velocity between the mass center and surface of Earth , 2012 .

[8]  J. Freymueller,et al.  Seasonal hydrological loading in southern Alaska observed by GPS and GRACE , 2012 .

[9]  F. Landerer,et al.  Accuracy of scaled GRACE terrestrial water storage estimates , 2012 .

[10]  J. Thepaut,et al.  The ERA‐Interim reanalysis: configuration and performance of the data assimilation system , 2011 .

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

[12]  S. Swenson,et al.  Satellites measure recent rates of groundwater depletion in California's Central Valley , 2011 .

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

[14]  Empirical Modeling of Solar Radiation Pressure Forces Affecting GPS Satellites , 2010 .

[15]  F. Nievinski,et al.  Can we measure snow depth with GPS receivers? , 2009 .

[16]  C. Faunt,et al.  Groundwater availability of the Central Valley Aquifer, California , 2009 .

[17]  E. Small,et al.  Use of GPS receivers as a soil moisture network for water cycle studies , 2008 .

[18]  D. Chambers,et al.  Estimating Geocenter Variations from a Combination of GRACE and Ocean Model Output , 2008 .

[19]  P. Stark.  Bounded-Variable Least-Squares: an Algorithm and Applications , 2008 .

[20]  Jeffrey T. Freymueller,et al.  Coseismic deformation of the 2002 Denali fault earthquake: Contributions from synthetic aperture radar range offsets , 2007 .

[21]  Jürgen Kusche,et al.  Surface mass redistribution inversion from global GPS deformation and Gravity Recovery and Climate Experiment (GRACE) gravity data , 2005 .

[22]  M. A. D'Alessio,et al.  Slicing up the San Francisco Bay Area : Block kinematics and fault slip rates from GPS-derived surface velocities , 2005 .

[23]  Michael B. Heflin,et al.  Interseismic strain accumulation and anthropogenic motion in metropolitan Los Angeles , 2005 .

[24]  J. D. Tarpley,et al.  The multi‐institution North American Land Data Assimilation System (NLDAS): Utilizing multiple GCIP products and partners in a continental distributed hydrological modeling system , 2004 .

[25]  Jeffrey P. Walker,et al.  THE GLOBAL LAND DATA ASSIMILATION SYSTEM , 2004 .

[26]  Lifeng Luo,et al.  Snow process modeling in the North American Land Data Assimilation System (NLDAS): 1. Evaluation of model‐simulated snow cover extent , 2003 .

[27]  Lifeng Luo,et al.  Snow process modeling in the north american Land Data Assimilation System (NLDAS): 2. Evaluation of model simulated snow water equivalent : GEWEX Continental-Scale International Project, Part 3 (GCIP3) , 2003 .

[28]  R. Bürgmann,et al.  Interactions between the Landers and Hector Mine, California, Earthquakes from Space Geodesy, Boundary Element Modeling, and Time-Dependent Friction , 2002 .

[29]  D. Argus,et al.  Present tectonic motion across the Coast Ranges and San Andreas fault system in central California , 2001 .

[30]  Gerald W. Bawden,et al.  Tectonic contraction across Los Angeles after removal of groundwater pumping effects , 2001, Nature.

[31]  K. Heki.  Seasonal Modulation of Interseismic Strain Buildup in Northeastern Japan Driven by Snow Loads , 2001, Science.

[32]  F. Bryan,et al.  Time variability of the Earth's gravity field: Hydrological and oceanic effects and their possible detection using GRACE , 1998 .

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

[34]  T. van Dam,et al.  Displacements of the Earth's surface due to atmospheric loading: Effects on gravity and baseline measurements , 1987 .

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

[36]  W. Farrell.  Deformation of the Earth by surface loads , 1972 .