The redistribution of past and present ice and ocean loading on Earth's surface causes solid Earth deformation and geoid changes, known as glacial isostatic adjustment. The deformation is controlled by elastic and viscous material parameters, which are inhomogeneous in the Earth. We present a new viscoelastic solid Earth deformation model in ASPECT (Advanced Solver for Problems in Earth's ConvecTion): a modern, massively parallel, open-source finite element code originally designed to simulate convection in the Earth's mantle. We show the performance of solid Earth deformation in ASPECT and compare solutions to TABOO, a semianalytical code, and Abaqus, a commercial finite element code. The maximum deformation and deformation rates using ASPECT agree within 2.6% for the average percentage difference with TABOO and Abaqus on glacial cycle (∼100 kyr) and contemporary ice melt (∼100 years) timescales. This gives confidence in the performance of our new solid Earth deformation model. We also demonstrate the computational efficiency of using adaptively refined meshes, which is a great advantage for solid Earth deformation modeling. Furthermore, we demonstrate the model performance in the presence of lateral viscosity variations in the upper mantle and report on parallel scalability of the code. This benchmarked code can now be used to investigate regional solid Earth deformation rates from ice age and contemporary ice melt. This is especially interesting for low-viscosity regions in the upper mantle beneath Antarctica and Greenland, where it is not fully understood how ice age and contemporary ice melting contribute to geodetic measurements of solid Earth deformation.

The rate at which global mean sea level (GMSL) rose during the 20th century is uncertain, with little consensus between various reconstructions that indicate rates of rise ranging from 1.3 to 2 mm•y^-1. Here we present a 20th-century GMSL reconstruction computed using an area-weighting technique for averaging tide gauge records that both incorporates up-to-date observations of vertical land motion (VLM) and corrections for local geoid changes resulting from ice melting and terrestrial freshwater storage and allows for the identification of possible differences compared with earlier attempts. Our reconstructed GMSL trend of 1.1 ± 0.3 mm•y^-1 (1s) before 1990 falls below previous estimates, whereas our estimate of 3.1 ± 1.4 mm•y^-1 from 1993 to 2012 is consistent with independent estimates from satellite altimetry, leading to overall acceleration larger than previously suggested. This feature is geographically dominated by the Indian Ocean-Southern Pacific region, marking a transition from lower-than-average rates before 1990 toward unprecedented high rates in recent decades. We demonstrate that VLM corrections, area weighting, and our use of a common reference datum for tide gauges may explain the lower rates compared with earlier GMSL estimates in approximately equal proportion. The trends and multidecadal variability of our GMSL curve also compare well to the sum of individual contributions obtained from historical outputs of the Coupled Model Intercomparison Project Phase 5. This, in turn, increases our confidence in process-based projections presented in the Fifth Assessment Report of the Intergovernmental Panel on Climate Change.