A long-standing hypothesis for synchronous global ice-sheet evolution on orbital timescales invokes an interhemispheric sea level forcing, whereby sea-level rise due to ice loss in the Northern Hemisphere (NH) in response to insolation and greenhouse gas forcing causes grounding line retreat of marine-based sectors of the Antarctic Ice Sheet (AIS) 1-3 . Recent evidence indicates that the AIS experienced substantial millennial-scale variability during and after the last deglaciation 4-7 , further suggesting a possible sea-level forcing. Global sea-level change from ice-sheet mass loss is strongly nonuniform 8 , however, suggesting that the response of AIS grounding lines to NH sea-level forcing is likely more complicated than previously considered 1,2,6 . Here we show, using a coupled ice sheet -global sea-level model, that a large or rapid NH sea-level forcing during
A coupled ice sheet-Earth-sea level model is used to estimate the modern elevations of shoreline features that were formed at high sea level stands during the warm mid Pliocene~3 million years ago. Knowledge of global mean sea level during this period is important as an indicator of possible future ice sheet retreat and sea level rise. However, local shoreline elevations can deviate from the eustatic mean by various geologic processes over the last 3 million years, including glacial isostatic adjustment of the solid Earth and gravitational field due to both Pliocene ice-cover changes and more recent glacial cycles. Our coupled model includes glacial isostatic adjustment processes and simulates Antarctic ice sheet, global sea level, and solid Earth variations in the warmest mid-Pliocene and over the last 40,000 years. Global maps of estimated modern elevations of Pliocene shoreline markers are produced for a standard radial profile of Earth viscosity and lithospheric thickness. Results are compared to an earlier study with an uncoupled Earth-sea level model and a different methodology (Raymo et al., Nature Geoscience, 2011, https://doi.org:/ 10.1038/ngeo1118). As in that study, Pliocene shoreline elevations diverge significantly from the eustatic value in widespread regions, especially in the vicinity of present and former ice sheets. In some other regions, elevations are close to eustatic. The results emphasize that care should be taken in interpreting elevations of paleo-shoreline markers.Although previous work has addressed each of these processes independently, an all-inclusive model would require substantial development and would involve large uncertainties, especially associated with dynamic POLLARD ET AL.
Abstract. Retreat and advance of ice sheets perturb the gravitational field, solid surface and rotation of the Earth, leading to spatially variable sea-level changes over a range of timescales O(100−6 years), which in turn feed back onto ice-sheet dynamics. Coupled ice-sheet–sea-level models have been developed to capture the interactive processes between ice sheets, sea level and the solid Earth, but it is computationally challenging to capture short-term interactions O(100−2 years) precisely within longer O(103−6 years) simulations. The standard forward sea-level modelling algorithm assigns a uniform temporal resolution in the sea-level model, causing a quadratic increase in total CPU time with the total number of input ice history steps, which increases with either the length or temporal resolution of the simulation. In this study, we introduce a new “time window” algorithm for 1D pseudo-spectral sea-level models based on the normal mode method that enables users to define the temporal resolution at which the ice loading history is captured during different time intervals before the current simulation time. Utilizing the time window, we assign a fine temporal resolution O(100−2 years) for the period of ongoing and recent history of surface ice and ocean loading changes and a coarser temporal resolution O(103−6 years) for earlier periods in the simulation. This reduces the total CPU time and memory required per model time step while maintaining the precision of the model results. We explore the sensitivity of sea-level model results to the model temporal resolution and show how this sensitivity feeds back onto ice-sheet dynamics in coupled modelling. We apply the new algorithm to simulate sea-level changes in response to global ice-sheet evolution over two glacial cycles and the rapid collapse of marine sectors of the West Antarctic Ice Sheet in the coming centuries and provide appropriate time window profiles for each application. The time window algorithm reduces the total CPU time by ∼ 50 % in each of these examples and changes the trend of the total CPU time increase from quadratic to linear. This improvement would increase with longer simulations than those considered here. Our algorithm also allows for coupling time intervals of annual temporal scale for coupled ice-sheet–sea-level modelling of regions such as West Antarctica that are characterized by rapid solid Earth response to ice changes due to the thin lithosphere and low mantle viscosities.
Abstract. Accurate glacial isostatic adjustment (GIA) modelling in the cryosphere is required for interpreting satellite, geophysical and geological records and for assessing the feedbacks of Earth deformation and sea-level change on marine ice-sheet grounding lines. GIA modelling in areas of active ice loss in West Antarctica is particularly challenging because the ice is underlain by laterally varying mantle viscosities that are up to several orders of magnitude lower than the global average, leading to a faster and more localised response of the solid Earth to ongoing and future ice-sheet retreat and necessitating GIA models that incorporate 3-D viscoelastic Earth structure. Improvements to GIA models allow for computation of the viscoelastic response of the Earth to surface ice loading at sub-kilometre resolution, and ice-sheet models and observational products now provide the inputs to GIA models at comparably unprecedented detail. However, the resolution required to accurately capture GIA in models remains poorly understood, and high-resolution calculations come at heavy computational expense. We adopt a 3-D GIA model with a range of Earth structure models based on recent seismic tomography and geodetic data to perform a comprehensive analysis of the influence of grid resolution on predictions of GIA in the Amundsen Sea Embayment (ASE) in West Antarctica. Through idealised sensitivity testing down to sub-kilometre resolution with spatially isolated ice loading changes, we find that a grid resolution of ∼ 13 of the radius of the load or higher is required to accurately capture the elastic response of the Earth. However, when we consider more realistic, spatially coherent ice loss scenarios based on modern observational records and future ice-sheet model projections and adopt a viscoelastic Earth, we find that predicted deformation and sea-level change along the grounding line converge to within 5 % with grid resolutions of 7.5 km or higher, and to within 2 % for grid resolutions of 3.75 km and higher, even when the input ice model is on a 1 km grid. Furthermore, we show that low mantle viscosities beneath the ASE lead to viscous deformation that contributes to the instrumental record on decadal timescales and equals or dominates over elastic effects by the end of the 21st century. Our findings suggest that for the range of resolutions of 1.9–15 km that we considered, the error due to adopting a coarser grid in this region is negligible compared to the effect of neglecting viscous effects and the uncertainty in the adopted mantle viscosity structure.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
customersupport@researchsolutions.com
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
This site is protected by reCAPTCHA and the Google Privacy Policy and Terms of Service apply.
Copyright © 2025 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.