We provide a new analysis of glacial isostatic adjustment (GIA) with the goal of assembling the model uncertainty statistics required for rigorously extracting trends in surface mass from the Gravity Recovery and Climate Experiment (GRACE) mission. Such statistics are essential for deciphering sea level, ocean mass, and hydrological changes because the latter signals can be relatively small (≤2 mm/yr water height equivalent) over very large regions, such as major ocean basins and watersheds. With abundant new >7 year continuous measurements of vertical land motion (VLM) reported by Global Positioning System stations on bedrock and new relative sea level records, our new statistical evaluation of GIA uncertainties incorporates Bayesian methodologies. A unique aspect of the method is that both the ice history and 1‐D Earth structure vary through a total of 128,000 forward models. We find that best fit models poorly capture the statistical inferences needed to correctly invert for lower mantle viscosity and that GIA uncertainty exceeds the uncertainty ascribed to trends from 14 years of GRACE data in polar regions.

No abstract

Geodetic investigations of crustal motions in the Amundsen Sea sector of West Antarctica and models of ice-sheet evolution in the past 10,000 years have recently highlighted the stabilizing role of solid-Earth uplift on polar ice sheets. One critical aspect, however, that has not been assessed is the impact of short-wavelength uplift generated by the solid-Earth response to unloading over short time scales close to ice-sheet grounding lines (areas where the ice becomes afloat). Here, we present a new global simulation of Antarctic evolution at high spatiotemporal resolution that captures all solid Earth processes that affect ice sheets and show a projected negative feedback in grounding line migration of 38% for Thwaites Glacier 350 years in the future, or 26.8% reduction in corresponding sea-level contribution.

Ice sheets and continental hydrology changes on decadal time scales are the dominant drivers of decadal scale polar motion.

[1] Glacier volume is known for less than 0.1% of the world's glaciers, but this information is needed to quantify the impacts of glacier changes on global sea level and regional water resources. Observations indicate a power-law relation between glacier area and volume, with an exponent g ≈ 1.36. Through numerical simulations of 3D, high-order glacier mechanics, we demonstrate how different topographic and climatic settings, glacier flow dynamics, and the degree of disequilibrium with climate systematically affect the volumearea relation. We recommend more accurate scaling relations through characterization of individual glacier shape, slope and size. An ensemble of 280 randomly-generated valley glaciers spanning a spectrum of plausible glaciological conditions yields a steady-state exponent g = 1.46. This declines to 1.38 for glaciers that are 100 years into a sustained retreat, which corresponds exceptionally well with the observed value for present-day glaciers.

No abstract

Ocean dynamics, land motion, and changes in Earth's gravitational and rotational fields cause local sea level change to deviate from the rate of global mean sea level rise. Here we use observations and simulations of spatial structure in sea level change to estimate the likelihood that these processes cause sea level trends in the longest and highest‐quality tide gauge records to be systematically biased relative to the true global mean rate. The analyzed records have an average twentieth century rate of approximately 1.6 mm/yr, but based on the locations of these gauges, we show that the simple average underestimates the twentieth century global mean rate by 0.1 ± 0.2 mm/yr. Given the distribution of potential sampling biases, we find that <1% probability that observed trends from the longest and highest‐quality tide gauge records are consistent with global mean rates less than 1.4 mm/yr.

Due to the lack of measurements of ice velocity, mass balance, glacier geometry and other baseline data, model-based studies of glacial systems in the Nepal Himalaya are very limited. Here a numerical ice-flow model has been developed for glacier AX010 in order to study its relation to local climate and investigate the possible causes of its general retreat since the end of the Little Ice Age. First, an attempt is made to simulate the historical front variations, considering each climatic parameter separately. Good agreement between the observations and model projections can be obtained under the assumption that variations in glacier front position are a response to changes in temperature alone. The same assumption is made about future changes to explore the 21st-century evolution of the glacier. Under a no-change scenario, the glacier will retreat by another ∽600m by AD 2100, whereas it is projected to vanish completely during this century for all trends with a temperature rise larger than +2.5˚C by AD 2100 with respect to the 1980–99 mean. With constant precipitation at the 1980–99 mean, the model predicts that the glacier will cease to exist at AD 2083, 2056 or 2049 if the temperature rises linearly by 3˚C, 4.5˚C or 6˚C respectively by the end of this century. With an additional range of precipitation changes between –30% and +30%, the life expectancy of glacier AX010 varies by 18, 6 and 2 years for the respective temperature rises. Thus the role of precipitation becomes minimal for the higher trends of temperature rise.