- Award ID(s):
- 1745074
- NSF-PAR ID:
- 10180987
- Date Published:
- Journal Name:
- Journal of Climate
- Volume:
- 33
- Issue:
- 2
- ISSN:
- 0894-8755
- Page Range / eLocation ID:
- 443 to 459
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
Abstract. Accurate glacial isostatic adjustment (GIA) modelling in the cryosphere is required for interpreting satellite, geophysical and geological recordsand for assessing the feedbacks of Earth deformation and sea-level change on marine ice-sheet grounding lines. GIA modelling in areas of active ice lossin West Antarctica is particularly challenging because the ice is underlain by laterally varying mantle viscosities that are up to several orders ofmagnitude lower than the global average, leading to a faster and more localised response of the solid Earth to ongoing and future ice-sheet retreatand necessitating GIA models that incorporate 3-D viscoelastic Earth structure. Improvements to GIA models allow for computation of the viscoelasticresponse of the Earth to surface ice loading at sub-kilometre resolution, and ice-sheet models and observational products now provide the inputs toGIA models at comparably unprecedented detail. However, the resolution required to accurately capture GIA in models remains poorly understood, andhigh-resolution calculations come at heavy computational expense. We adopt a 3-D GIA model with a range of Earth structure models based on recentseismic tomography and geodetic data to perform a comprehensive analysis of the influence of grid resolution on predictions of GIA in the AmundsenSea Embayment (ASE) in West Antarctica. Through idealised sensitivity testing down to sub-kilometre resolution with spatially isolated ice loadingchanges, we find that a grid resolution of ∼ 13 of the radius of the load or higher is required to accurately capture the elasticresponse of the Earth. However, when we consider more realistic, spatially coherent ice loss scenarios based on modern observational records andfuture ice-sheet model projections and adopt a viscoelastic Earth, we find that predicted deformation and sea-level change along the grounding lineconverge to within 5 % with grid resolutions of 7.5 km or higher, and to within 2 % for grid resolutions of 3.75 km andhigher, even when the input ice model is on a 1 km grid. Furthermore, we show that low mantle viscosities beneath the ASE lead to viscousdeformation that contributes to the instrumental record on decadal timescales and equals or dominates over elastic effects by the end of the 21stcentury. Our findings suggest that for the range of resolutions of 1.9–15 km that we considered, the error due to adopting a coarser gridin this region is negligible compared to the effect of neglecting viscous effects and the uncertainty in the adopted mantle viscosity structure.more » « less
-
SUMMARY Earth structure beneath the Antarctic exerts an important control on the evolution of the ice sheet. A range of geological and geophysical data sets indicate that this structure is complex, with the western sector characterized by a lithosphere of thickness ∼50–100 km and viscosities within the upper mantle that vary by 2–3 orders of magnitude. Recent analyses of uplift rates estimated using Global Navigation Satellite System (GNSS) observations have inferred 1-D viscosity profiles below West Antarctica discretized into a small set of layers within the upper mantle using forward modelling of glacial isostatic adjustment (GIA). It remains unclear, however, what these 1-D viscosity models represent in an area with complex 3-D mantle structure, and over what geographic length-scale they are applicable. Here, we explore this issue by repeating the same modelling procedure but applied to synthetic uplift rates computed using a realistic model of 3-D viscoelastic Earth structure inferred from seismic tomographic imaging of the region, a finite volume treatment of GIA that captures this complexity, and a loading history of Antarctic ice mass changes inferred over the period 1992–2017. We find differences of up to an order of magnitude between the best-fitting 1-D inferences and regionally averaged depth profiles through the 3-D viscosity field used to generate the synthetics. Additional calculations suggest that this level of disagreement is not systematically improved if one increases the number of observation sites adopted in the analysis. Moreover, the 1-D models inferred from such a procedure are non-unique, that is a broad range of viscosity profiles fit the synthetic uplift rates equally well as a consequence, in part, of correlations between the viscosity values within each layer. While the uplift rate at each GNSS site is sensitive to a unique subspace of the complex, 3-D viscosity field, additional analyses based on rates from subsets of proximal sites showed no consistent improvement in the level of bias in the 1-D inference. We also conclude that the broad, regional-scale uplift field generated with the 3-D model is poorly represented by a prediction based on the best-fitting 1-D Earth model. Future work analysing GNSS data should be extended to include horizontal rates and move towards inversions for 3-D structure that reflect the intrinsic 3-D resolving power of the data.
-
Abstract Glacial isostatic adjustment explains topographic change in formerly and currently glaciated regions, but the role of small (∼100s km) regions of unusually low‐viscosity mantle is poorly understood. We developed viscoelastic models with low‐viscosity regions in the upper mantle, and measured the effect of these regions on solid earth uplift resulting from contemporary surface ice melt. We found viscous uplift occurring on decadal timescales above the low‐viscosity region, at rates comparable to or larger than those from elastic uplift or the viscous response to ice age melting. We find that uplift rates are sensitive to the location, dimensions, and viscosity of the low‐viscosity region, and that the largest uncertainty in uplift rates likely comes from the low‐viscosity region's horizontal extent. Rapid viscous ground uplift can impact ice dynamics if the low‐viscosity region is located close to an ice sheet margin, as for Antarctica and Greenland.
-
Abstract 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.
-
Geodetic, seismic, and geological evidence indicates that West Antarctica is underlain by low-viscosity shallow mantle. Thus, as marine-based sectors of the West Antarctic Ice Sheet (WAIS) retreated during past interglacials, or will retreat in the future, exposed bedrock will rebound rapidly and flux meltwater out into the open ocean. Previous studies have suggested that this contribution to global mean sea level (GMSL) rise is small and occurs slowly. We challenge this notion using sea level predictions that incorporate both the outflux mechanism and complex three-dimensional viscoelastic mantle structure. In the case of the last interglacial, where the GMSL contribution from WAIS collapse is often cited as ~3 to 4 meters, the outflux mechanism contributes ~1 meter of additional GMSL change within ~1 thousand years of the collapse. Using a projection of future WAIS collapse, we also demonstrate that the outflux can substantially amplify GMSL rise estimates over the next century.more » « less