Rise of Great Lakes Surface Water, Sinking of the Upper Midwest of the United States, and Viscous Collapse of the Forebulge of the Former Laurentide Ice Sheet

Great Lakes water levels rose 0.7–1.5 m from 2013 to 2019, increasing surface water volume by 285 km3. Solid Earth's elastic response to the increased mass load is nearly known: The Great Lakes floor fell 8–23 mm, and the adjacent land fell 3–14 mm. Correcting GPS measurements for this predicted elastic loading (1) straightens position‐time series, making the evolution of position more nearly a constant velocity and (2) reduces estimates of subsidence rate in Wisconsin, Michigan, and southern Ontario by 0.5–2 mm/yr, improving constraints on postglacial rebound. GPS records Wisconsin and Michigan to have subsided at 1–4 mm/yr. We find this sinking to be produced primarily by viscous collapse of the former Laurentide ice sheet forebulge and secondarily by elastic Great Lakes loading. We infer water on land in the Great Lakes watershed to be total water change observed by GRACE minus Great Lakes surface water smeared by a Gaussian distribution. Water stored on land each year reaches a maximum in March, 6 months before Great Lakes water levels peak in September. The seasonal oscillation of water on land in the Great Lakes basin, 100 km3 (0.20 m water thickness), is twice that in a hydrology model. In the seasons, groundwater in the Great Lakes watershed increases by 60 km3 (0.12 m) each autumn and winter and decreases by roughly an equivalent amount each spring and summer. In the long term, groundwater volume remained constant from 2004 to 2012 but increased by 50 km3 (0.10 m) from 2013 to 2019.

[1]  Pascal Willis,et al.  GipsyX/RTGx, a new tool set for space geodetic operations and research , 2020, 2004.13124.

[2]  Junhua Li,et al.  Surface deformation observed by InSAR shows connections with water storage change in Southern Ontario , 2020 .

[3]  G. Blewitt,et al.  Improved GPS position time series spanning up to 25 years for over 18,000 stations resulting from IGS Repro3 products, JPL's GipsyX software, and more advanced modeling techniques , 2019 .

[4]  H. McMillan,et al.  A Decade of Water Storage Changes Across the Contiguous United States From GPS and Satellite Gravity , 2019, Geophysical Research Letters.

[5]  E. Small,et al.  Subsurface Water Dominates Sierra Nevada Seasonal Hydrologic Storage , 2019, Geophysical Research Letters.

[6]  Grzegorz Michalak,et al.  The GFZ GRACE RL06 Monthly Gravity Field Time Series: Processing Details and Quality Assessment , 2019, Remote. Sens..

[7]  L. Rivera,et al.  LoadDef: A Python‐Based Toolkit to Model Elastic Deformation Caused by Surface Mass Loading on Spherically Symmetric Bodies , 2019, Earth and space science.

[8]  C. Shum,et al.  Monitoring Water Level Variations over the Great Lakes Using Contemporary Satellite Geodetic Observations , 2018 .

[9]  R. Kopp,et al.  Postglacial relative sea-level histories along the eastern Canadian coastline , 2018, Quaternary Science Reviews.

[10]  G. Blewitt,et al.  Harnessing the GPS Data Explosion for Interdisciplinary Science , 2018, Eos.

[11]  Geoffrey Blewitt,et al.  A Robust Estimation of the 3‐D Intraplate Deformation of the North American Plate From GPS , 2018 .

[12]  Donald F. Argus,et al.  Comment on “An Assessment of the ICE‐6G_C (VM5a) Glacial Isostatic Adjustment Model” by Purcell et al. , 2018 .

[13]  R. Reedy,et al.  Global models underestimate large decadal declining and rising water storage trends relative to GRACE satellite data , 2018, Proceedings of the National Academy of Sciences.

[14]  A. Moore,et al.  Sustained Water Loss in California's Mountain Ranges During Severe Drought From 2012 to 2015 Inferred From GPS , 2017 .

[15]  A. Gronewold,et al.  Resolving Hydrometeorological Data Discontinuities along an International Border , 2017 .

[16]  Reply to Comment by W. R. Peltier, D. F. Argus, and R. Drummond on “An Assessment of the ICE6G_C (VM5a) Glacial Isostatic Adjustment Model” , 2017 .

[17]  R. König,et al.  A new high-resolution model of non-tidal atmosphere and ocean mass variability for de-aliasing of satellite gravity observations: AOD1B RL06 , 2017 .

[18]  W. Peltier,et al.  Space-geodetic and water level gauge constraints on continental uplift and tilting over North America: regional convergence of the ICE-6G_C (VM5a/VM6) models , 2017 .

[19]  B. Tapley,et al.  Long‐term Caspian Sea level change , 2017 .

[20]  Matt A. King,et al.  Uncertainty in geocenter estimates in the context of ITRF2014 , 2017 .

[21]  K. Lambeck,et al.  The North American Late Wisconsin ice sheet and mantle viscosity from glacial rebound analyses , 2017 .

[22]  Ralf Schmid,et al.  IGS14/igs14.atx: a new framework for the IGS products , 2016 .

[23]  Srinivas Bettadpur,et al.  High‐resolution CSR GRACE RL05 mascons , 2016 .

[24]  M. Watkins,et al.  Quantifying and reducing leakage errors in the JPL RL05M GRACE mascon solution , 2016 .

[25]  Z. Altamimi,et al.  ITRF2014: A new release of the International Terrestrial Reference Frame modeling nonlinear station motions , 2016 .

[26]  Song S. Qian,et al.  Hydrological drivers of record‐setting water level rise on Earth's largest lake system , 2016 .

[27]  J. Ray,et al.  The IGS contribution to ITRF2014 , 2016, Journal of Geodesy.

[28]  T. Dixon,et al.  A three-dimensional surface velocity field for the Mississippi Delta: Implications for coastal restoration and flood potential , 2015 .

[29]  M. Watkins,et al.  Improved methods for observing Earth's time variable mass distribution with GRACE using spherical cap mascons , 2015 .

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

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

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

[33]  Brian F. Thomas,et al.  River basin flood potential inferred using GRACE gravity observations at several months lead time , 2014 .

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

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

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

[37]  Maria Grazia D'Urso,et al.  On a generalized Love's problem , 2013, Comput. Geosci..

[38]  R. Dill,et al.  Numerical simulations of global‐scale high‐resolution hydrological crustal deformations , 2013 .

[39]  J. Camp,et al.  Antarctica, Greenland and Gulf of Alaska land-ice evolution from an iterated GRACE global mascon solution , 2013, Journal of Glaciology.

[40]  Craig A. Stow,et al.  Coasts, water levels, and climate change: A Great Lakes perspective , 2013, Climatic Change.

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

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

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

[44]  W. van der Wal,et al.  Detectability of groundwater storage change within the Great Lakes Water Basin using GRACE , 2012 .

[45]  Duncan Carr Agnew,et al.  SPOTL: Some Programs for Ocean-Tide Loading , 2012 .

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

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

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

[49]  D. Scavia,et al.  Great Lakes Restoration & the Threat of Global Warming , 2008 .

[50]  Donald F. Argus,et al.  Defining the translational velocity of the reference frame of Earth , 2007 .

[51]  H. Schuh,et al.  Troposphere mapping functions for GPS and very long baseline interferometry from European Centre for Medium‐Range Weather Forecasts operational analysis data , 2006 .

[52]  M. Craymer,et al.  Present-day tilting of the Great Lakes region based on water level gauges , 2005 .

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

[54]  Katharine Hayhoe,et al.  Confronting climate change in the Great Lakes region : impacts on our communities and ecosystems , 2003 .

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

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

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

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

[59]  D. Meredith Temperature effects on Great Lakes water balance studies , 1975 .

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