Abstract
Recent studies suggest that Antarctica has the potential to contribute up to ~15 m of sea-level rise over the next few centuries. The evolution of the Antarctic Ice Sheet is driven by a combination of climate forcing and non-climatic feedbacks. In this review we focus on feedbacks between the Antarctic Ice Sheet and the solid Earth, and the role of these feedbacks in shaping the response of the ice sheet to past and future climate changes. The growth and decay of the Antarctic Ice Sheet reshapes the solid Earth via isostasy and erosion. In turn, the shape of the bed exerts a fundamental control on ice dynamics as well as the position of the grounding line—the location where ice starts to float. A complicating issue is the fact that Antarctica is situated on a region of the Earth that displays large spatial variations in rheological properties. These properties affect the timescale and strength of feedbacks between ice-sheet change and solid Earth deformation, and hence must be accounted for when considering the future evolution of the ice sheet.
Similar content being viewed by others
Introduction
The solid Earth, along with the oceans and the atmosphere, exerts a strong influence on the dynamics of the Antarctic Ice Sheet (AIS). The pre-glacial topography of the Antarctic continent determined the location and style of glacial inception ~34 Ma ago1, whereas today the shape of the bed and the properties of the ice-bed interface exert a first-order control on the contemporary pattern of ice flow2. In the intervening period the bed of the AIS has been continuously reshaped by erosion and sedimentation, periodically flexed by glacial isostasy, and abruptly altered by tectonic and volcanic activity, with the latter two processes also playing a role in determining the thermal conditions at the bed (Fig. 1). Basal heat flux affects ice flow via its influence on subglacial hydrology and ice rheology3, but its spatial variation is currently poorly quantified4. Parts of Antarctica are underlain by active volcanoes, notably West Antarctica (Fig. 1). This is of interest because in other regions, e.g., Iceland5, volcanism has been shown to increase during periods of deglaciation. However, little is known of this phenomenon in Antarctica6 and a detailed review of this effect is therefore not possible. We instead focus on those processes that are better known and that control long-wavelength changes to the shape of the bed beneath Antarctica – glacial isostasy and erosion – and feedbacks between these processes and ice sheet evolution.
Interactions between ice sheets and the solid Earth have long been studied within the field of Glacial Isostatic Adjustment (GIA)—defined here as the response of the solid Earth and the global gravity field to changes in the distribution of ice and water on Earth’s surface. The first numerical models of GIA were developed in the 1970s7 but they have received renewed attention over the last decade, reflecting the important role they play in the interpretation of satellite measurements of contemporary ice-sheet change8. A number of recent studies have also sought to understand the strength of feedbacks acting in the opposite direction, that is, the impact of solid Earth deformation on ice dynamics and the potential for this deformation to delay or prevent unstable ice-sheet retreat9,10. Such feedbacks were first considered in the 1980s11 but have only recently begun to be implemented into fully-coupled models12, where the evolving shape of the solid Earth and the depth of the ocean adjacent to an ice sheet grounded below sea level act as fundamental boundary conditions on the dynamics of the ice sheet.
The strength of any feedbacks between glacial isostasy and ice dynamics depends on the rate at which the solid Earth responds to ice-sheet change, which, in turn, depends on the rheological properties of the mantle. Seismic evidence13 indicates that there are significant spatial variations in mantle properties beneath Antarctica, which suggests mantle viscosities, and hence relaxation timescales, may vary by up to several orders of magnitude from the global mean. Indeed, in the northern Antarctic Peninsula, geodetic evidence14 suggests that contemporary ice loss is triggering a viscous response orders of magnitude more rapidly than was previously assumed possible in Antarctica – on decadal rather than millennial timescales. Feedbacks on ice dynamics are likely to be enhanced in such regions15, and the ongoing GIA signal is likely to be dominated by recent ice-sheet change (few millennia or less), as opposed to the ice loss that followed the Last Glacial Maximum (LGM)16. Here, we define the GIA signal to be the ongoing response of the solid Earth, the gravity field, and/or relative sea level to past ice-sheet change.
The growth and decay of the AIS is driven by a combination of climate forcing and non-climatic feedbacks, but modelling studies that seek to understand the controls on AIS change often neglect to consider how the ice sheet alters its own boundary conditions. Over the timescale of the last deglaciation, GIA model output17,18 suggests that the response to surface load change can alter bed slopes across West Antarctica by 0.25–0.4 m/km (values may be an underestimate due to assumptions of strong mantle rheology) and that water-depth change around the margin of the ice sheet can deviate from eustatic by >100 m. Coupled ice sheet-sea level models can now be used to quantify the impact of such changes on ice sheet evolution15,19. Over timescales of millions of years, changes to the shape of the ice sheet bed will additionally reflect processes associated with erosion and deposition, the isostatic response to sediment redistribution, and mantle convection. Feedbacks between ice sheet evolution and long-term landscape evolution have been hypothesised20, but they have not been modelled within a coupled framework. The impacts of long-wavelength isostasy-driven changes to ice sheet boundary conditions are reviewed here, but smaller-scale subglacial controls, such as the material properties of the bed and variations in subglacial hydrology and geomorphology, are discussed elsewhere2.
In this review we proceed by briefly summarising the current state of knowledge on ice-sheet change and glacial isostasy across Antarctica, before discussing the impact of spatial variations in Earth rheology on glacial isostasy and the impact of GIA on ice dynamics. Feedbacks between ice sheet and landscape evolution are discussed, and we conclude by identifying a number of future research priorities that link ice history, Earth rheology, and the overarching issue of global sea-level change.
Antarctic ice-sheet change and glacial isostasy
Forward model predictions of GIA-related solid Earth deformation, gravity-field change, and polar motion rely on: a reconstruction of the spatial and temporal pattern of past ice loading, a viscoelastic Earth model that describes the time-dependent response of the solid Earth to surface loading changes, and the iterative consideration of physical feedbacks associated with polar motion and coastline position as they are altered by the deformation of the solid Earth and the gravity field21,22. While the third component is well defined by physical theory, Antarctic ice-sheet reconstructions and solid Earth rheology are subject to substantial uncertainty and debate.
GIA modelling can be used to infer past large-scale ice-sheet change via comparison of model output with a range of constraint data relating to sea-level change and solid Earth deformation18,23. On a global scale, such data-model comparisons have been used to estimate the timing of ice-sheet growth and retreat over glacial cycles24,25, but the combined dependence of model output on Earth rheology and ice history leads to non-uniqueness26, and it remains challenging to determine how ice was partitioned between the different ice sheets27.
Reconstructions of the LGM AIS suggest an ice sheet larger than present, equivalent to ~5–22 m global mean sea level28. This range reflects differences in methodology and the interpretation of data constraints. The most common reconstruction approaches consider either numerical models of ice-sheet dynamics19,29 or local geological, geomorphological, and geodetic data constraints on past ice extent30, or a hybrid31,32. Data constraints are sparse as they rely on sampling of spatially limited bedrock exposures33,34 or interpretation of ice core records35, while numerical ice-sheet models are restricted by a limited representation of reality. Continent-wide ice-sheet reconstructions can be produced by simulating the response of dynamic models to past climate changes19,31,36,37,38, but at present there is not enough data to tightly constrain such models. Consequently, there is significant variation in the predicted magnitude and spatial pattern of isostatic deformation across Antarctica due to past ice-sheet change39 (Fig. 2).
It is important to note that some of the differences in Fig. 2 reflect the timescale of ice-sheet change considered to contribute to the deformation signal; traditional forward models of GIA do not consider ice-sheet change during the last millennium (Fig. 2a-f), but this is accounted for in solutions derived via coupled modelling (Fig. 2g, h) or data inversion (Fig. 2i-l). Altogether, four approaches to determining the present-day GIA signal are represented in Fig. 2, reflecting a number of recent methodological advances. We do not seek to quantitatively assess the accuracy of each GIA solution, partly due to the difficulty of defining validation data sets, but we briefly note the advantages and limitations of each approach in Table 1.
The effect of GIA on geodetic and gravimetric measurements, and the resulting impact on estimates of ice-mass change, has driven much of the recent interest in Antarctic GIA. Contemporary ice-sheet change can be determined from analysis of Earth’s time-varying gravity field, as measured by the Gravity Recovery and Climate Experiment (GRACE) satellites between 2002 and 201740 and GRACE Follow-On (from May 2018). However, the time-varying gravity field also contains a gravity-change signal associated with past ice-sheet change; GIA modelling is used to determine the ‘GIA correction’ that should be applied to remove the signal associated with past change8. Uncertainties associated with this process remain the dominant source of error in gravimetry-based estimates of the Antarctic contribution to contemporary sea-level rise41. However, newer GIA model predictions17,30,32 yield improved inter-model agreement, as do inverse GIA solutions derived by combining ice elevation and GRACE data42,43,44. Recent estimates for the magnitude of the net Antarctic GIA signal vary over the range ~40–80 Gt/yr45, but basin-level differences remain substantial and in some cases different studies do not even agree on the sign of the mass change in each basin43, which hampers the advance of glaciological insights into the processes governing present change.
Differences in estimates of the present-day GIA signal are particularly acute in the region of the Amundsen Sea embayment, where forward models have substantially less GIA signal than inverse solutions. This may be due to the fact that the GIA signal in this low mantle viscosity region predominantly reflects significant decadal-to-centennial ice load change that is not accounted for in most forward models46,47. For the same reasons, forward models underestimate the signal in the northern Antarctic Peninsula, but then so do inverse solutions because the GIA signal in this region has a shorter spatial wavelength than inverse solutions can resolve (~300 km). Comparison of the new generation of GIA models with GPS-derived uplift rates—which must be corrected for elastic deformation associated with contemporaneous ice-mass change48,49—demonstrates that important differences remain, notably in West Antarctica and especially in the Amundsen Sea and northern Antarctic Peninsula regions39.
Earth structure and rheology beneath Antarctica
The pattern and magnitude of the solid Earth response to ice-sheet growth and decay is strongly dependent on the rheology of the interior of the solid Earth, and the choice of rheological model has a large effect on the modelled GIA response (see Box 1). For example, ice history is often tuned in tandem with Earth rheology, so adoption of a different rheological model can lead to large differences in the assumed ice history. Earth rheology beneath Antarctica is spatially variable and, as detailed in the next section, this has implications for the behaviour of the ice sheet.
A number of rheological models have been adopted to explain the solid Earth response to surface load change across Antarctica, including linear Maxwell viscoelastic rheology, power law rheology, and Burgers rheology (see Box 1). Regardless of the choice of rheological model, the need for future modelling efforts to consider spatial variations in Earth structure beneath Antarctica is now clear. Seismic studies of the upper mantle beneath Antarctica suggest large lateral variations in material properties13,50,51,52,53, with greater heterogeneity than observed in areas of Northern Hemisphere continental glaciation54. East Antarctica shows high upper mantle seismic velocities characteristic of cold cratonic regions worldwide, whereas West Antarctica shows upper mantle structure consistent with much warmer tectonically active zones (Fig. 3a-c). The range of seismic velocities observed beneath Antarctica provides strong evidence for lateral variations in mantle viscosity, but absolute values are currently poorly known. Preliminary GIA studies that explore the effect of including lateral variations in mantle viscosity reveal significant differences in predicted patterns and rates of deformation across Antarctica55,56,57,58, motivating the need for better constraints on absolute mantle viscosity in this region.
Upper mantle viscosity variations of several orders of magnitude beneath West Antarctica have been quantified by modelling observed GPS uplift in specific settings. Nield et al.14 identified rapid viscoelastic deformation occurring in the northern Antarctic Peninsula due to a well-observed change in glacial loading resulting from the 2002 breakup of the Larsen B Ice Shelf. Comparison of GPS time series and modelled uplift in the region suggests upper mantle viscosities of between 6 × 1017 and 2 × 1018 Pa s in the northern Antarctic Peninsula. Similarly, Zhao et al.59 used the changes in ice load resulting from thinning of Fleming Glacier over recent decades to estimate upper mantle viscosity, 500 km further south along the southern Antarctic Peninsula, to be at least 2 × 1019 Pa s. Other estimates of upper mantle viscosity in West Antarctica include 1–3 × 1020 Pa s beneath the southwestern Weddell Sea60 and 4 × 1018 Pa s for the Amundsen Sea Coast47. These represent variations of two orders of magnitude. In contrast, estimates of spatially-averaged upper mantle viscosity beneath the whole of Antarctica range from 2 × 1020 to 1 × 1021 Pa s (ref. 16,17), which are similar to estimates of upper mantle viscosity in Fennoscandia61. The emerging picture suggests cratonic East Antarctica is characterised by higher upper mantle viscosity than West Antarctica, with exceptionally low upper mantle viscosity, on the order of 1018 to 1019 Pa s, beneath some regions of West Antarctica.
Existing constraints on mantle viscosity across Antarctica draw on our ability to measure the solid Earth response to known surface-load change. In regions where this is not possible, the three-dimensional shear velocity structure of the upper mantle beneath Antarctica can be used to estimate the lateral and depth variation of viscosity. Although there is no direct physical correspondence between shear velocity and viscosity, their variation in the upper mantle is largely controlled by temperature62,63. Ivins and Sammis64 formulated an approach for converting shear velocity anomalies to viscosity anomalies relative to a global reference 1-D viscosity model, using the observed scaling of shear velocity and density anomalies to infer the temperature scaling, and then using olivine diffusion creep (linear) flow laws to constrain the viscosity perturbation. Wu et al.65 proposed a similar method that uses experimentally determined temperature derivatives of shear velocity, including the effect of anelasticity. An alternative method is to estimate the temperature structure of the mantle from the shear velocity structure, and then use a composite rheology that computes the strain associated with both diffusion and dislocation creep mechanisms56. Since the dislocation creep rheology is non-linear, the calculated effective viscosity of the mantle will depend on the stress field. One key assumption in any of these methods is that the observed seismic velocity anomalies result entirely from temperature anomalies, whereas it is well known that a portion of the velocity perturbations result from (spatially variable) compositional anomalies66. Wu et al.65 adjust their viscosity structure by the percentage of the velocity anomalies thought to be caused by thermal variations, which in their study was estimated to be between 65% and 100%. An alternative approach would involve correcting for (poorly known) spatial variations in mantle composition.
Using the approach outlined in Wu et al.65, the estimated mantle viscosity beneath Antarctica, at depths of 100 and 250 km, varies by 2–3 orders of magnitude across the continent (Fig. 4), with extremely low viscosity predicted beneath the Antarctic Peninsula, the Amundsen Sea coast, Marie Byrd Land, and the Transantarctic Mountain Front. The implications of these viscosity variations, and in particular the anomalously low viscosities beneath parts of West Antarctica, are explored in the next section.
Feedbacks between GIA processes and ice dynamics
It has long been recognised that the evolution of ice sheets is influenced by the geometry and deformation of the underlying solid Earth67, and that the stability and dynamics of marine ice sheets (ice sheets which rest on ground below sea level) are sensitive to the depth of water at their grounding lines, i.e., the point where they begin to float68,69,70. Marine-terminating ice sheets such as Antarctica lose most of their mass following the flow of grounded ice across their grounding lines into floating ice shelves (Fig. 5). The ice flux across the grounding line is very sensitive to the thickness of ice there, and the thickness is in turn proportional to the depth of water such that a small increase in water depth at the grounding line leads to a large reduction in grounded ice70. Marine ice sheets are widely thought to be prone to runaway retreat when resting upon beds that slope down towards the interior of the ice sheet, i.e., reverse bed slopes70, as is the case for West Antarctica and much of East Antarctica. This canonical argument is based on the idea that when the grounding line retreats into deeper water, the ice loss at the grounding line increases, leading to further retreat (Fig. 5b). However, this reasoning does not account for the fact that as ice mass loss occurs, the solid Earth rebounds, and the gravitational attraction of the ice on the surrounding water is diminished, lowering the local sea surface (Fig. 5c). The bathymetry shallows and bed slopes are altered near the margins of the ice sheet where ice mass loss is occurring. Within a modelling framework this local shallowing of water reduces the loss of ice across the grounding line, acting to stabilise the modelled ice sheet71, slowing, and in some cases halting, migration of the modelled grounding line along reverse bed slopes12. Ongoing viscous uplift of the bed following a halt in ice sheet retreat can also initiate readvance of the grounding line in marine areas (e.g., Fig. 6d, ref. 15).
Water depth changes can also influence the degree to which ice shelves are able to stabilise the ice sheet. Ice shelves play an important role in the stability and evolution of marine ice sheets by providing resistance to the flow of grounded ice across the grounding line. The spreading of ice shelves is inhibited by friction along their sides and base, particularly where the ice shelf becomes locally grounded on bumps or pinning points in the bathymetry, forming ice rises72. A local decrease in water depth can enhance grounding of the ice shelf at ice rises, stabilising the ice sheet, while an increase in water depth can lead to ungrounding at the ice rise, enhancing flow across the grounding line of the ice sheet. Modelling studies of the AIS have shown that accelerated retreat and thinning and eventually large-scale collapse of marine sectors of the ice sheet can occur when the surrounding ice shelves break up73, but that if ice shelves are able to re-ground, this can have a stabilising effect on the ice sheet74. Furthermore, changes in bathymetry may have implications for ocean circulation and heat transfer under the ice shelves. The role of these and other ice shelf processes in controlling the dynamics of the grounded portion of the ice sheet are discussed in more detail in a companion paper by Smith et al. (manuscript submitted).
Initial coupled modelling studies, using simplified flowline ice-sheet modelling and bedrock geometry, demonstrated the stabilising feedback of sea-level changes on marine ice sheets12,71. Recently, a series of more realistic coupled models have been developed that capture Antarctic ice sheet and ice shelf dynamics, global sea level and solid Earth deformation, and the interactions between these systems9,10,15,19,36. It has been shown that GIA-related sea-level and solid Earth changes, including changes to the slope of the underlying bed, alter the stress field of the ice sheet in a way that acts to dampen and slow past19,75 and future9,10,76 ice-sheet growth and retreat in Antarctica (Fig. 6a, c). An important process that is also accounted for in these coupled models is the feedback between isostatically-driven ice surface elevation change and surface mass balance75,77,78.
The strength of the feedback between GIA processes and ice dynamics depends on the rheology of the solid Earth. Models representing the palaeo79,80,81 and long-term future evolution of ice sheets over millennia73,82 generally account for ice-load-driven Earth deformation by adopting simplified treatments such as the Elastic Lithosphere Relaxing Asthenosphere (ELRA) model, which treats the asthenosphere as a time-lagged relaxation towards equilibrium11,73,83, or a model with flow in a viscous half-space below an elastic plate82,84. The Earth rheology models adopted in the newly-developed fully-coupled models described above19,36,75 capture more realistically the full multi-normal-mode response of the Earth to both ice and water loading, thus enabling the computation of gravitationally self-consistent variations of the sea surface, and Earth rotational effects. By solving the sea-level equation7, these studies also account for migrating shorelines, including migration into regions previously occupied by marine-based ice. Ice model simulations15 over the last deglaciation incorporating an ELRA bed model (black line) and the full sea-level coupling (blue and red lines) are shown in Fig. 6a. Note that the bed topography at the start of each simulation shown in Fig. 6a will be different; this is to ensure that the final modelled topography is close to the modern observed topography in each case. Due to the complexity and spatial variations of bedrock geometry, ice dynamics and climate-ice interactions, the importance of GIA processes on ice sheet evolution cannot be quantified universally and must be considered on a case-by-case basis for different regions, time periods, and climate forcings.
When considering short-term change, ice-sheet models designed to simulate decadal to centennial-scale transient ice dynamics in response to future warming85,86 often do not include bedrock deformation and sea-level changes. Simulated Antarctic ice volume changes under moderate future climate warming with fixed bed and sea surface are compared with results from a coupled model9 in which these surfaces are allowed to vary (Fig. 6c). The impact of sea-level changes on ice dynamics has previously been considered negligible on these timescales, in particular under strong warming, because the viscous response of the Earth is only considered relevant on millennial timescales and longer. However, as discussed above, the Earth structure underneath the AIS is highly variable, and viscosities may be as low as 1018 Pa s beneath parts of West Antarctica, leading to substantial (i.e., metres to tens of metres of) viscoelastic uplift occurring on centennial or even decadal timescales14,47, with consequent implications for ice sheet evolution.
Recent coupled modelling studies have begun to quantify the sensitivity of predicted ice dynamics, GIA and crustal deformation in Antarctica, as well as the AIS’s contribution to past and future global sea-level change, to the adopted radial viscosity structure of the solid Earth9,10,15. Note that, in a coupled modelling context, altering the Earth structure influences GIA predictions both directly by altering the timing and geometry of the Earth’s response to surface loading, and indirectly by changing the ice loading itself. The blue and red lines in Fig. 6a-d compare coupled model simulations9,15 over the last 40 ka, during a Pliocene warm period, and into the future under moderate and strong climate warming scenarios (see figure caption for details) adopting two different radially-varying models of Earth viscosity and lithospheric thickness within the coupled model. One is a relatively high viscosity Earth model with a thick lithosphere (HV), similar to models adopted in global GIA studies, and the other is representative of the Earth structure beneath much of West Antarctica, having a thinned lithosphere and a zone of low viscosity down to 200 km in the upper mantle (LVZ). For each of the time periods considered in Fig. 6, the timing and extent of ice-sheet retreat is sensitive to the adopted Earth rheology, and the predicted contribution of Antarctica to sea-level rise is lowered in the LVZ simulation as compared to the HV simulation. In the LVZ case the solid Earth responds faster to ice loading changes and deformation is localised to near the edges of the ice sheet where the ice loss occurs. The resulting sea-level fall at the grounding line can more effectively act to slow ice-sheet retreat as compared to the simulations with the HV Earth model.
A comparison of the future simulations shown in Fig. 6c, d highlights that the influence of Earth structure on ice sheet evolution depends on both the strength of the climate forcing and the physics adopted in the ice-sheet model. For a moderate climate warming, uplift of the LVZ Earth model preserves much of West Antarctica as compared to the simulation with the HV Earth model (Fig. 6c). While, for the simulation where strong RCP 8.5 climate warming is applied and new rapid-retreat-promoting ice physics are added (hydrofracturing and cliff failure73), West Antarctica collapses early on regardless of the choice of Earth model, and differences in the ice sheet evolution between the LVZ and HV simulations occur mostly in East Antarctica (Fig. 6d).
Given the sensitivity of the ice-Earth-sea level system in Antarctica to differences in radially varying Earth structure (Fig. 6), and the known variability in Earth structure beneath Antarctica (Fig. 3), no single, radial Earth model is able to accurately represent all of Antarctica and hence consideration of lateral variations in Earth structure is motivated. Incorporating lateral variations in Earth structure into a coupled ice sheet–sea level model represents a large jump in computational cost. Gomez et al.87 developed the first coupled ice sheet–sea level model that incorporates 3-D variations in Earth structure and applied it to model Antarctic evolution over the last deglaciation. They show that substantial localised differences can arise in ice cover, sea level and crustal movement, which introduces substantial uncertainty into the GIA corrections that should be applied to contemporary geodetic observations.
Feedbacks between ice-sheet change and landscape evolution
The previous section describes the impact of ice-load-driven changes to the elevation of the solid Earth surface and gravity field on ice dynamics. Many ice-sheet modelling studies account for ice-driven isostatic adjustment88, but over multiple glacial cycles several other processes also alter the underlying topography. Thermal subsidence following tectonic extension has lowered central West Antarctica by several hundred metres over the last ~34 Ma (ref. 89) and changes in dynamic topography over the last 3 Ma are postulated to have altered the stability of some sectors of East Antarctica90. However, since the inception of the AIS, glacially-driven erosion and sedimentation, and the accompanying isostatic response91,92,93,94, have been the main drivers of topographic change across Antarctica. This topographic change has, in turn, played a role in controlling the sensitivity of the ice sheet to climate forcing. This section focuses on the extent to which feedbacks between the ice sheet and the solid Earth have shaped the long-term evolution of both.
Prior to 34 Ma, Antarctica was characterised by a fluvial landscape95 and higher mean topography96. Dated offshore sediment packages are testament to the volume of material that has been removed from Antarctica since 34 Ma (ref. 97). Backstripping techniques can be used to unpick the history of progradation and isostatic subsidence, and hence reconstruct the palaeotopography of the continental shelf98. It is more challenging to reconstruct onshore palaeotopography due to the difficulty of determining the volume of material that has been removed, and the source of that material. By using assumed drainage pathways89 or process-based erosion modelling tuned to match offshore sediment volumes92, and accounting for the isostatic response to ice and ocean loading, as well as sediment erosion and deposition, bounds can be placed on the pre-glaciation topography of Antarctica96,99. Current palaeotopography reconstructions have been used to suggest that the West Antarctic Ice Sheet could have formed much earlier than previously thought, during the Eocene-Oligocene transition100. However, when two palaeotopography end members are used to define boundary conditions for ice-sheet growth under ‘cold’ mid-Miocene conditions, the difference in modelled ice volume is equivalent to 20 m sea level101, with the shallower topography leading to the growth of a larger ice sheet. It is clear that the configuration of the ice sheet is highly sensitive to the underlying topography, but the question remains as to how Antarctic topography has evolved, and the degree to which this evolution was coupled to ice-sheet change via glacial erosion and the resulting isostatic response.
To understand the role of internal, i.e., non-climatic, processes in controlling long-term ice-sheet change it is necessary to consider feedbacks between ice dynamics, erosion, deposition, isostasy, and water depth change. Extending the coupled ice sheet - sea level modelling approach described above to account for landscape evolution processes would allow the testing of hypotheses that seek to explain the dramatic changes that took place as Antarctica evolved from a terrestrial to a marine ice sheet against a background of long-term cooling. Specifically, between 34 and 14 Ma the volume of the AIS fluctuated significantly102, then at 14 Ma there was a switch to cold, polar conditions103, which resulted in the establishment of a cold-based, less erosive ice sheet104. A range of hypotheses have been proposed to explain this transition to the current ‘icehouse’ world105,106, but it remains to be tested whether landscape evolution played a role, either through a rapid change in topographic boundary conditions or by causing the system to pass some internal threshold as large portions of the ice sheet became marine grounded.
Ice sheet-driven processes that will have had an impact on ice-sheet dynamics include progradation of the continental shelf and long wavelength erosion and deposition2. The resulting isostatic response of onshore uplift and offshore subsidence will have focused erosion on the inner shelf92, leading to the development of a reverse bed slope that is too deep to support ice-sheet advance98,107 and is susceptible to unstable grounding line retreat68. These processes are conceptually simple, but they will have been subject to control by spatially variable sea-level change. As the mean topography of Antarctica decreased through erosion, the ice sheet will have become more sensitive to sea-level change108 and ocean forcing109. Applying appropriate boundary conditions to understand long-term ice-sheet change is challenging due to uncertainties associated with palaeotopography and the growth of the Northern Hemisphere ice sheets, but initial studies that have explored the impact of spatially variable sea-level change on the pre-Pleistocene ice sheet highlight the importance of accounting for peripheral bulge growth when considering ice110 and ocean111 dynamics, and the damping effects of ice sheet - sea level feedbacks in regions of weak mantle viscosity15. The impact of spatial variability in Earth rheology on coupled ice sheet - landscape evolution has not yet been investigated.
Regions of the AIS that are currently located on over-deepened reverse bed slopes have been suggested to be fundamentally unstable86. The details of how the ice sheet and the underlying topography reached this state are not known, but it is clear that feedbacks will have played a role as the fluvially- and tectonically-controlled topography of 34 Ma was eroded and warped into its current configuration94, and it has been suggested that future landscape evolution will render the ice sheet even more unstable20. Models that seek to explore the past and future long-term evolution of the AIS should account for evolution of the underlying topography due to both sediment redistribution and isostasy, ideally within a coupled framework. The development of such models will improve understanding of sediment transport pathways, permitting stronger conclusions to be drawn from provenance studies112. Better quantification of bedrock erosion rates and offshore sediment packages are important targets for improving our understanding of feedbacks between ice sheet and landscape change.
Outstanding problems and future outlook
The previous section highlights a number of advances that are required to further our understanding of the long-term evolution of the AIS and the role of changing topography. These include the development of more sophisticated numerical models that consider feedbacks between ice dynamics, glacial erosion, isostasy, and global sea-level change. Such processes will also play a role in controlling contemporary and future ice dynamics but, as highlighted above, quantifying present change is exacerbated by the need to interpret geodetic observations in terms of the response to both past and present ice-sheet change. In this section we discuss several areas that should be prioritised as we seek to better quantify the GIA signal across Antarctica and hence understand the processes responsible for contemporary ice-sheet change.
Characterisation of absolute mantle viscosity is preliminary. Beneath West Antarctica low viscosity mantle approaches isostatic equilibrium more quickly compared with global-average timescales. The consequence of this is that the present-day rate of adjustment will depend heavily on recent, and relatively localised, surface loading changes87. Characteristic relaxation times of a Maxwell fluid may be approximated by dividing the viscosity by the shear modulus (~7 × 1010 Pa for the upper mantle). The relaxation time of mantle material with a viscosity of 1019 Pa s is therefore only ~5 years, or 1–2 orders of magnitude faster than global averages. Actual relaxation times will be somewhat larger due to the presence of higher-viscosity layers within the rheological profile, but this calculation illustrates that, for regions underlain by low viscosity mantle, detailed knowledge of glacial load changes over the past decades to centuries is needed to quantify the GIA signal in regions of low viscosity113. For viscosities of 1018 Pa s or lower the viscous relaxation almost occurs contemporaneously with the surface load change and, by implication, the elastic deformation, with the viscous component being around an order of magnitude greater than the elastic component14,47. Such a rapid response is at odds with the traditional idea that the present-day GIA signal is related to post-LGM ice loss, although, given the much longer relaxation time of the lower mantle, some component of the present-day GIA signal is likely still associated with large-scale, post-LGM ice-sheet change.
Across East Antarctica, spatial variations in Earth rheology are currently poorly constrained, as they are across all offshore regions, due to significant uncertainty in Earth structure resulting from the absence of seismic stations. Across all Antarctica, uncertainties also exist in relation to the rheological law that should be used to describe mantle deformation. Although traditional GIA models do not parameterise viscosity to be directly related to the physical properties of the mantle, physically-based approximations to mantle creep processes are being implemented in new models56, and these require quantification of mantle temperature, grain size, and water content, which have varying measurement uncertainties. Mantle temperatures, and hence viscosities, may be inferred from seismic velocity perturbations, but different approaches yield different results. When using a power-law rheology, the stress in the mantle prior to the change in surface load also plays a role in defining the viscosity, but this is often taken as zero114. Such stresses are expected to be largely a function of long-term mantle convection processes but could be more complex in the shallow mantle and asthenosphere due to ongoing isostatic relaxation or earthquake-related deformation.
Constraints on ice loading history are extremely sparse across all time periods and locations (Fig. 7). However, this problem is particularly acute for late Holocene loading changes across West Antarctica and post-LGM changes in East Antarctica. Very few data record changes to the AIS from the late Holocene up until the commencement of the satellite record115, although there is evidence of a dynamic West Antarctic Ice Sheet116,117,118,119 and large changes in net accumulation during this period120,121. Progress can be made by studying the internal structure of the ice sheet to determine past flow patterns122,123, but at present, global and Antarctic-focused GIA models generally assume no ice-load change over the last 1–2 ka (see refs. 14,17,116,124 for some recent regional exceptions). Combined with low mantle viscosities in much of West Antarctica, this means that the predicted upper-mantle component of present-day deformation is likely erroneous. For vast sections of East Antarctica almost no data exist on past ice extent and retreat history125 (Fig. 7). Furthermore, there are few data to constrain the spatial extent and history of the post-LGM margin due to very limited bathymetric sampling. These limitations provide the motivation for new fieldwork and the reanalysis of historic and high resolution palaeo datasets59,126.
Data with which to validate or constrain the GIA signal across Antarctica are also sparse (Fig. 7). Less than 20 records of past sea-level change exist across the continent (e.g., through dating raised beach terraces), and these are presently limited to the last 12–15 ka (ref. 31). It is possible to extend data to earlier time-periods through the study of sediment from submerged offshore basins127. Opportunities exist for new absolute gravity128, InSAR129 and, especially, GNSS measurements at rock outcrops with large sections of East Antarctica sparsely or un-observed. Continuous measurements allow both increased precision but also, in regions of low mantle viscosity, probing of mantle properties by time-varying surface loading14,59. The largely unexploited horizontal deformation field promises new insights into local or regional-scale deformation patterns of both elastic and viscoelastic processes59,130, although a robust approach to remove tectonic plate motion and post-seismic deformation prevents continent-wide analyses at present131,132. A methodology is not yet available to precisely measure bedrock displacement under the present ice sheet or offshore, and this precludes model validation for these regions.
An important supplement to GNSS measurements of deformation at discrete bedrock outcrops is provided by inverse estimates of the spatially continuous deformation field42,43; such estimates are useful in separating competing conventional ‘forward model’ predictions of the GIA signal124. While possessing their own inter-solution variation due largely to differences in altimeter snow-densification corrections, inverse estimates tend to differ more from the spatial patterns of deformation predicted by current forward models (Fig. 2). In regions where mantle viscosities are thought to be ~1018 Pa s or lower, viscoelastic deformation in response to contemporary surface load change will perturb the deformation and gravity fields due to the short response time of the mantle, but the effect of this signal on inverse solutions is yet to be addressed.
Quantification of model prediction uncertainty is immature with most attempts limited to sampling a generally small set of Earth models. Ice history uncertainty is rarely taken into account, in part due to the sparse information available from which to sample probabilistically133 or lack of rigorous measurement uncertainties to propagate formally134. This prevents the robust propagation of GIA uncertainties into other quantities that make use of GIA model predictions, for example GRACE-derived ice mass balance estimates.
A practical problem for those wishing to employ viscoelastic models is the absence of open source software that includes state-of-the-art model physics. While open source software are available and widely used for purely elastic135,136 or viscoelastic137,138 solutions, the viscoelastic models do not solve the full sea-level equation139, which requires consideration of gravitationally self-consistent meltwater redistribution on a rotating, spherical Earth with polar-motion feedbacks and migrating coastlines, and they do not have the capacity to account for 3-D Earth structure, compressibility and a selection of rheologies (e.g., transient, linear and power-law). Given that there is just one real Earth, a standardised viscoelastic software framework that allows the consistent treatment of GIA and post-seismic viscoelastic deformation would have distinct advantages; it would pave the way for new insights into Earth’s rheology, provide a framework to advance Earth’s viscoelastic predictability, and help resolve ongoing debates regarding model robustness140,141,142.
Summary
Recent modelling advances have highlighted the importance of understanding the role of the solid Earth in moderating the response of the global ice sheets to past and future climate change9,10,19,36. It has long been understood that the growth and decay of an ice sheet will alter the shape of the solid Earth and result in spatially variable sea-level change7. We are beginning to understand how those changes affect the dynamics of ice sheets71, but further work is needed to quantify feedbacks between ice sheet evolution and landscape evolution before we can fully explain the long-term evolution of the Antarctic Ice Sheet.
Coupled models that consider interactions between ice sheets and the solid Earth via dynamic modelling of both systems represent a new type of model that can provide insight into the extent, timing, and mechanisms of past ice sheet evolution and sea-level change. Such models also provide a new tool for understanding the link between past climate and global sea-level change. Coupled models have been used to demonstrate that the solid Earth response to ice loss, and the accompanying fall in sea level, will have a stabilising effect on grounding line dynamics71. This mechanism has the potential to slow or even halt grounding line retreat along a reverse bed slope, depending on the rate of solid Earth rebound. Such negative feedbacks have been modelled in association with ice loss across West Antarctica9,10, but crucially the rheology of the underlying mantle is poorly constrained, and so the strength of the feedback cannot currently be quantified.
Interdisciplinary approaches to determining mantle rheology reveal large differences between East and West Antarctica52,56, with low viscosity mantle inferred to lie beneath West Antarctica. The short mantle relaxation time associated with such regions means that estimates of the contemporary GIA signal derived by forward modelling will be biased if they do not account for ice-sheet change over the last few millennia, with implications for estimates of current ice mass balance8. Improved quantification of late Holocene ice history, as well as tighter constraints on mantle rheology, are needed to reduce uncertainty on the contemporary GIA signal across Antarctica. Low viscosity mantle will also enhance the strength of the stabilising effect of GIA on grounding line dynamics, highlighting the importance of considering such feedbacks when modelling the future evolution of the Antarctic Ice Sheet. A number of studies have successfully incorporated 3-D Earth structure into GIA models beneath modern ice sheets56,143,144, and it is now apparent that coupled modelling and inclusion of 3-D Earth structure should both be considered when modelling solid Earth-cryosphere feedbacks15. Such models are being developed87 but further interdisciplinary work, combining modelling and observational approaches, is needed to calibrate such models and better understand controls on the evolution of the Antarctic Ice Sheet.
References
Rose, K. C. et al. Early East Antarctic Ice Sheet growth recorded in the landscape of the Gamburtsev Subglacial Mountains. Earth Planet Sci. Lett. 375, 1–12 (2013).
Colleoni, F. et al. Spatio-temporal variability of processes across Antarctic ice-bed-ocean interfaces. Nat. Commun. 9, 2289 (2018).
Seroussi, H., Ivins, E. R., Wiens, D. A. & Bondzio, J. Influence of a West Antarctic mantle plume on ice sheet basal conditions. J. Geophys. Res. Solid Earth 122, 7127–7155 (2017).
Martos, Y. M. et al. Heat flux distribution of Antarctica unveiled. Geophys. Res. Lett. 44, 11417–11426 (2017).
Maclennan, J., Jull, M., McKenzie, D., Slater, L. & Gronvold, K. The link between volcanism and deglaciation in Iceland. Geochem. Geophy. Geosy. 3, 1062 (2002).
Smellie, J. L. & Edwards, B. R. Glaciovolcanism on Earth and Mars. (Cambridge University Press, Cambridge, UK 2016).
Farrell, W. E. & Clark, J. A. On postglacial sea level. Geophys. J. R. Astron. Soc. 46, 647–667 (1976).
King, M. A. et al. Lower satellite-gravimetry estimates of Antarctic sea-level contribution. Nature 491, 586–589 (2012).
Gomez, N., Pollard, D. & Holland, D. Sea-level feedback lowers projections of future Antarctic Ice-Sheet mass loss. Nat. Commun. 6, 8798 (2015).
Konrad, H., Sasgen, I., Pollard, D. & Klemann, V. Potential of the solid-Earth response for limiting long-term West Antarctic Ice Sheet retreat in a warming climate. Earth Planet. Sci. Lett. 432, 254–264 (2015).
Lingle, C. S. & Clark, J. A. A Numerical-Model of Interactions between a Marine Ice-Sheet and the Solid Earth - Application to a West Antarctic Ice Stream. J. Geophys. Res. Oceans 90, 1100–1114 (1985).
Gomez, N., Pollard, D., Mitrovica, J. X., Huybers, P. & Clark, P. U. Evolution of a coupled marine ice sheet-sea level model. J. Geophys. Res. Earth 117, F01013 (2012).
Heeszel, D. S. et al. Upper mantle structure of central and West Antarctica from array analysis of Rayleigh wave phase velocities. J. Geophys. Res. Solid Earth 121, 1758–1775 (2016).
Nield, G. A. et al. Rapid bedrock uplift in the Antarctic Peninsula explained by viscoelastic response to recent ice unloading. Earth Planet. Sci. Lett. 397, 32–41 (2014).
Pollard, D., Gomez, N. & DeConto, R. M. Variations of the Antarctic Ice Sheet in a coupled ice sheet-Earth-sea level model: sensitivity to viscoelastic Earth properties. J. Geophys. Res. Earth Surf. 122, 2124–2138 (2017).
Ivins, E. R. & James, T. S. Antarctic glacial isostatic adjustment: a new assessment. Antarct. Sci. 17, 541–553 (2005).
Whitehouse, P. L., Bentley, M. J., Milne, G. A., King, M. A. & Thomas, I. D. A new glacial isostatic adjustment model for Antarctica: calibrated and tested using observations of relative sea-level change and present-day uplift rates. Geophys J. Int 190, 1464–1482 (2012).
Peltier, W. R., Argus, D. F. & Drummond, R. Space geodesy constrains ice age terminal deglaciation: the global ICE-6G_C (VM5a) model. J. Geophys. Res. Solid Earth 120, 450–487 (2015).
Gomez, N., Pollard, D. & Mitrovica, J. X. A 3-D coupled ice sheet-sea level model applied to Antarctica through the last 40 ky. Earth Planet. Sci. Lett. 384, 88–99 (2013).
Jamieson, S. S. R. et al. The glacial geomorphology of the Antarctic ice sheet bed. Antarct. Sci. 26, 724–741 (2014).
Mitrovica, J. X., Milne, G. A. & Davis, J. L. Glacial isostatic adjustment on a rotating earth. Geophys. J. Int. 147, 562–578 (2001).
Martinec, Z. & Hagedoorn, J. The rotational feedback on linear-momentum balance in glacial isostatic adjustment. Geophys. J. Int. 199, 1823–1846 (2014).
Lambeck, K., Smither, C. & Johnston, P. Sea-level change, glacial rebound and mantle viscosity for northern Europe. Geophys. J. Int. 134, 102–144 (1998).
Peltier, W. R. Global glacial isostasy and the surface of the ice-age earth: The ICE-5G (VM2) model and GRACE. Annu. Rev. Earth Planet. Sci. 32, 111–149 (2004).
Lambeck, K., Rouby, H., Purcell, A., Sun, Y. Y. & Sambridge, M. Sea level and global ice volumes from the Last Glacial Maximum to the Holocene. Proc. Natl Acad. Sci. USA 111, 15296–15303 (2014).
Whitehouse, P. L. Glacial isostatic adjustment modelling: historical perspectives, recent advances, and future directions. Earth Surf. Dynam 6, 401–429 (2018).
Caron, L., Metivier, L., Greff-Lefftz, M., Fleitout, L. & Rouby, H. Inverting Glacial Isostatic Adjustment signal using Bayesian framework and two linearly relaxing rheologies. Geophys. J. Int. 209, 1126–1147 (2017).
Clark, P. U. & Tarasov, L. Closing the sea level budget at the Last Glacial Maximum. Proc. Natl Acad. Sci. USA 111, 15861–15862 (2014).
Golledge, N. R. et al. Glaciology and geological signature of the Last Glacial Maximum Antarctic ice sheet. Quat. Sci. Rev. 78, 225–247 (2013).
Ivins, E. R. et al. Antarctic contribution to sea-level rise observed by GRACE with improved GIA correction. J. Geophys. Res.: Solid Earth 118, 3126–3141 (2013).
Whitehouse, P. L., Bentley, M. J. & Le Brocq, A. M. A deglacial model for Antarctica: geological constraints and glaciological modelling as a basis for a new model of Antarctic glacial isostatic adjustment. Quat. Sci. Rev. 32, 1–24 (2012).
Argus, D. F., Peltier, W. R., Drummond, R. & Moore, A. W. 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. Geophys. J. Int. 198, 537–563 (2014).
Bentley, M. J. et al. Deglacial history of the West Antarctic Ice Sheet in the Weddell Sea embayment: constraints on past ice volume change. Geology 38, 411–414 (2010).
Clark, P. U. Deglacial history of the West Antarctic Ice Sheet in the Weddell Sea embayment: constraints on past ice volume change: COMMENT. Geology 39, E239–E239 (2011).
Parrenin, F. et al. 1-D-ice flow modelling at EPICA Dome C and Dome Fuji, East Antarctica. Clim. Past. 3, 243–259 (2007).
de Boer, B., Stocchi, P. & van de Wal, R. S. W. A fully coupled 3-D ice-sheet-sea-level model: algorithm and applications. Geosci. Model Dev. 7, 2141–2156 (2014).
Ritz, C., Rommelaere, V. & Dumas, C. Modeling the evolution of Antarctic ice sheet over the last 420,000 years: Implications for altitude changes in the Vostok region. J. Geophys Res-Atmos. 106, 31943–31964 (2001).
Huybrechts, P. Sea-level changes at the LGM from ice-dynamic reconstructions of the Greenland and Antarctic ice sheets during the glacial cycles. Quat. Sci. Rev. 21, 203–231 (2002).
Martín-Español, A. et al. An assessment of forward and inverse GIA solutions for Antarctica. J. Geophys. Res. Solid Earth 121, 6947–6965 (2016).
Velicogna, I. & Wahr, J. Measurements of time-variable gravity show mass loss in Antarctica. Science 311, 1754–1756 (2006).
Shepherd, A. et al. A reconciled estimate of ice-sheet mass balance. Science 338, 1183–1189 (2012).
Gunter, B. C. et al. Empirical estimation of present-day Antarctic glacial isostatic adjustment and ice mass change. Cryosphere 8, 743–760 (2014).
Martin-Espanol, A. et al. Spatial and temporal Antarctic Ice Sheet mass trends, glacio-isostatic adjustment, and surface processes from a joint inversion of satellite altimeter, gravity, and GPS data. J. Geophys. Res. Earth Surf. 121, 182–200 (2016).
Sasgen, I. et al. Joint inversion estimate of regional glacial isostatic adjustment in Antarctica considering a lateral varying Earth structure (ESA STSE Project REGINA). Geophys. J. Int. 211, 1534–1553 (2017).
Shepherd, A. et al. Mass balance of the Antarctic Ice Sheet from 1992 to 2017. Nature 558, https://doi.org/10.1038/s41586-018-0179-y (2018).
Groh, A. et al. An investigation of Glacial Isostatic Adjustment over the Amundsen Sea sector, West Antarctica. Glob. Planet. Change 98-99, 45–53 (2012).
Barletta, V. R. et al. Observed rapid bedrock uplift in Amundsen Sea Embayment promotes ice-sheet stability. Science 360, 1335 (2018).
Bevis, M. et al. Geodetic measurements of vertical crustal velocity in West Antarctica and the implications for ice mass balance. Geochem. Geophy. Geosy. 10, https://doi.org/10.1029/2009GC002642 (2009).
Thomas, I. D. et al. Widespread low rates of Antarctic glacial isostatic adjustment revealed by GPS observations. Geophys. Res. Lett. 38, L22302 (2011).
Ritzwoller, M. H., Shapiro, N. M., Levshin, A. L. & Leahy, G. M. Crustal and upper mantle structure beneath Antarctica and surrounding oceans. J. Geophys. Res. Sol. Earth 106, 30645–30670 (2001).
Morelli, A. & Danesi, S. Seismological imaging of the Antarctic continental lithosphere: a review. Glob. Planet Change 42, 155–165 (2004).
Lloyd, A. J. Seismic Tomography of Antarctica and the Southern Oceans: Regional and Continental Models from the Upper Mantle to the Transition Zone. (PhD thesis, Washington University in St Louis, USA 2018).
Shen, W. et al. The crust and upper mantle structure of central and West Antarctica from Bayesian inversion of Rayleigh wave and receiver functions. J. Geophys. Res. https://doi.org/10.1029/2017JB015346 (2018).
Schaeffer, A. J. & Lebedev, S. Global shear speed structure of the upper mantle and transition zone. Geophys. J. Int. 194, 417–449 (2013).
A., G., Wahr, J. & Zhong, S. J. Computations of the viscoelastic response of a 3-D compressible Earth to surface loading: an application to Glacial Isostatic Adjustment in Antarctica and Canada. Geophys. J. Int. 192, 557–572 (2013).
van der Wal, W., Whitehouse, P. L. & Schrama, E. J. O. Effect of GIA models with 3D composite mantle viscosity on GRACE mass balance estimates for Antarctica. Earth Planet. Sci. Lett. 414, 134–143 (2015).
Kaufmann, G., Wu, P. & Ivins, E. R. Lateral viscosity variations beneath Antarctica and their implications on regional rebound motions and seismotectonics. J. Geodyn. 39, 165–181 (2005).
Nield, G. A. et al. The impact of lateral variations in lithospheric thickness on glacial isostatic adjustment in West Antarctica. Geophys. J. Int. https://doi.org/10.1093/gji/ggy158 (2018).
Zhao, C. et al. Rapid ice unloading in the Fleming Glacier region, southern Antarctic Peninsula, and its effect on bedrock uplift rates. Earth Planet. Sci. Lett. 473, 164–176 (2017).
Wolstencroft, M. et al. Uplift rates from a new high-density GPS network in Palmer Land indicate significant late Holocene ice loss in the southwestern Weddell Sea. Geophys. J. Int. 203, 737–754 (2015).
Steffen, H. & Wu, P. Glacial isostatic adjustment in Fennoscandia-A review of data and modeling. J. Geodyn. 52, 169–204 (2011).
Hirth, G. & Kohlstedt, D. in Inside the Subduction Factory 138 Geophysical Monograph Series (ed. Eiler. J.) 83–105 (AGU, Washington D.C., USA 2003).
Faul, U. H. & Jackson, I. The seismological signature of temperature and grain size variations in the upper mantle. Earth Planet. Sci. Lett. 234, 119–134 (2005).
Ivins, E. R. & Sammis, C. G. On Lateral Viscosity Contrast in the Mantle and the Rheology of Low-Frequency Geodynamics. Geophys. J. Int. 123, 305–322 (1995).
Wu, P., Wang, H. S. & Steffen, H. The role of thermal effect on mantle seismic anomalies under Laurentia and Fennoscandia from observations of Glacial Isostatic Adjustment. Geophys. J. Int. 192, 7–17 (2013).
Lee, C. T. A. Compositional variation of density and seismic velocities in natural peridotites at STP conditions: Implications for seismic imaging of compositional heterogeneities in the upper mantle. J. Geophys. Res. Solid Earth 108, https://doi.org/10.1029/2003jb002413 (2003).
Le Meur, E. & Huybrechts, P. A model computation of the temporal changes of surface gravity and geoidal signal induced by the evolving Greenland ice sheet. Geophys. J. Int. 145, 835–849 (2001).
Weertman, J. Stability of the junction of an ice sheet and an ice shelf. J. Glaciol. 13, 3–13 (1974).
Mercer, J. H. West Antarctic Ice Sheet and CO2 Greenhouse Effect - Threat of Disaster. Nature 271, 321–325 (1978).
Schoof, C. Ice sheet grounding line dynamics: Steady states, stability, and hysteresis. J. Geophys. Res. Earth Surf. 112, F03S28 (2007).
Gomez, N., Mitrovica, J. X., Huybers, P. & Clark, P. U. Sea level as a stabilizing factor for marine-ice-sheet grounding lines. Nat. Geosci. 3, 850–853 (2010).
Matsuoka, K. et al. Antarctic ice rises and rumples: Their properties and significance for ice-sheet dynamics and evolution. Earth-Sci. Rev. 150, 724–745 (2015).
DeConto, R. M. & Pollard, D. Contribution of Antarctica to past and future sea-level rise. Nature 531, 591–597 (2016).
Kingslake, J. et al. Extensive retreat and re-advance of the West Antarctic Ice Sheet during the Holocene. Nature 558, 430 (2018).
Konrad, H. et al. The Deformational Response of a Viscoelastic Solid Earth Model Coupled to a Thermomechanical Ice Sheet Model. Surv. Geophys. 35, 1441–1458 (2014).
Adhikari, S. et al. Future Antarctic bed topography and its implications for ice sheet dynamics. Solid Earth 5, 569–584 (2014).
van den Berg, J., van de Wal, R. S. W., Milne, G. A. & Oerlemans, J. Effect of isostasy on dynamical ice sheet modeling: a case study for Eurasia. J. Geophys. Res. 113, B05412 (2008).
Abe-Ouchi, A. et al. Insolation-driven 100,000-year glacial cycles and hysteresis of ice-sheet volume. Nature 500, 190–193 (2013).
Le Meur, E. & Huybrechts, P. A comparison of different ways of dealing with isostasy: examples from modelling the Antarctic ice sheet during the last glacial cycle. Ann. Glaciol. 23, 309–317 (1996).
Golledge, N. R. et al. Antarctic contribution to meltwater pulse 1A from reduced Southern Ocean overturning. Nat. Commun. 5, 5107 (2014).
Pollard, D., Chang, W., Haran, M., Applegate, P. & DeConto, R. Large ensemble modeling of the last deglacial retreat of the West Antarctic Ice Sheet: comparison of simple and advanced statistical techniques. Geosci. Model Dev. 9, 1697–1723 (2016).
Golledge, N. R. et al. The multi-millennial Antarctic commitment to future sea-level rise. Nature 526, 421 (2015).
Pollard, D. & DeConto, R. M. Description of a hybrid ice sheet-shelf model, and application to Antarctica. Geosci. Model Dev. 5, 1273–1295 (2012).
Bueler, E., Lingle, C., & Brown, J. Fast computation of a viscoelastic deformable Earth model for ice-sheet simulations. Annals of Glaciology, 46, 97–105 (2007).
Cornford, S. L. et al. Century-scale simulations of the response of the West Antarctic Ice Sheet to a warming climate. Cryosphere 9, 1579–1600 (2015).
Joughin, I., Smith, B. E. & Medley, B. Marine Ice Sheet Collapse Potentially Under Way for the Thwaites Glacier Basin, West Antarctica. Science 344, 735–738 (2014).
Gomez, N., Latychev, K. & Pollard, D. A coupled ice sheet-sea level model incorporating 3D Earth structure: Variations in Antarctica during the last deglacial retreat. J. Clim. 31, 4041–4054 (2018).
Pollard, D. & DeConto, R. M. Modelling West Antarctic ice sheet growth and collapse through the past five million years. Nature 458, 329–332 (2009).
Wilson, D. S. & Luyendyk, B. P. West Antarctic paleotopography estimated at the Eocene-Oligocene climate transition. Geophys. Res. Lett. 36, L16302 (2009).
Austermann, J. et al. The impact of dynamic topography change on Antarctic ice sheet stability during the mid-Pliocene warm period. Geology 43, 927–930 (2015).
Stern, T. A., Baxter, A. K. & Barrett, P. J. Isostatic rebound due to glacial erosion within the Transantarctic Mountains. Geology 33, 221–224 (2005).
Jamieson, S. S. R., Sugden, D. E. & Hulton, N. R. J. The evolution of the subglacial landscape of Antarctica. Earth Planet. Sci. Lett. 293, 1–27 (2010).
Paxman, G. J. G. et al. Erosion-driven uplift in the Gamburtsev Subglacial Mountains of East Antarctica. Earth Planet. Sci. Lett. 452, 1–14 (2016).
Paxman, G. J. G. et al. Uplift and tilting of the Shackleton Range in East Antarctica driven by glacial erosion and normal faulting. J. Geophys Res-Solid Earth 122, 2390–2408 (2017).
Sugden, D. & Denton, G. Cenozoic landscape evolution of the Convoy Range to Mackay Glacier area, Transantarctic Mountains: Onshore to offshore synthesis. Geol. Soc. Am. Bull. 116, 840–857 (2004).
Wilson, D. S. et al. Antarctic topography at the Eocene-Oligocene boundary. Palaeogeogr. Palaeocl 335, 24–34 (2012).
Lindeque, A., Gohl, K., Wobbe, F. & Uenzelmann-Neben, G. Preglacial to glacial sediment thickness grids for the Southern Pacific Margin of West Antarctica. Geochem. Geophy. Geosy. 17, 4276–4285 (2016).
De Santis, L., Prato, S., Brancolini, G., Lovo, M. & Torelli, L. The Eastern Ross Sea continental shelf during the Cenozoic, implications for the West Antarctic ice sheet development. Glob. Planet Change 23, 173–196 (1999).
Paxman, G. J. G. et al. Bedrock Erosion Surfaces Record Former East Antarctic Ice Sheet Extent. Geophys. Res. Lett. 45, https://doi.org/10.1029/2018GL077268 (2018).
Wilson, D. S., Pollard, D., DeConto, R. M., Jamieson, S. S. R. & Luyendyk, B. P. Initiation of the West Antarctic Ice Sheet and estimates of total Antarctic ice volume in the earliest Oligocene. Geophys. Res. Lett. 40, 4305–4309 (2013).
Gasson, E., DeConto, R. M., Pollard, D. & Levy, R. H. Dynamic Antarctic ice sheet during the early to mid-Miocene. Proc. Natl Acad. Sci. USA 113, 3459–3464 (2016).
Naish, T. R. et al. Orbitally induced oscillations in the East Antarctic ice sheet at the Oligocene/Miocene boundary. Nature 413, 719–723 (2001).
Shevenell, A. E., Kennett, J. P. & Lea, D. W. Middle Miocene Southern Ocean cooling and Antarctic cryosphere expansion. Science 305, 1766–1770 (2004).
Lewis, A. R., Marchant, D. R., Ashworth, A. C., Hemming, S. R. & Machlus, M. L. Major middle Miocene global climate change: Evidence from East Antarctica and the Transantarctic Mountains. Geol. Soc. Am. Bull. 119, 1449–1461 (2007).
Holbourn, A., Kuhnt, W., Schulz, M. & Erlenkeuser, H. Impacts of orbital forcing and atmospheric carbon dioxide on Miocene ice-sheet expansion. Nature 438, 483–487 (2005).
Colleoni, F. et al. Past continental shelf evolution increased Antarctic ice sheet sensitivity to climatic conditions. Sci. Rep. Uk 8, 11323 (2018).
Taylor, J. et al. Topographic controls on post-Oligocene changes in ice-sheet dynamics, Prydz Bay region, East Antarctica. Geology 32, 197–200 (2004).
Denton, G. H. & Hughes, T. J. Global Ice-Sheet System Interlocked by Sea-Level. Quat. Res. 26, 3–26 (1986).
Bart, P. J. & Iwai, M. The overdeepening hypothesis: How erosional modification of the marine-scape during the early Pliocene altered glacial dynamics on the Antarctic Peninsula’s Pacific margin. Palaeogeogr. Palaeocl. 335, 42–51 (2012).
Stocchi, P. et al. Relative sea-level rise around East Antarctica during Oligocene glaciation. Nat. Geosci. 6, 380–384 (2013).
Rugenstein, M., Stocchi, P., von der Heydt, A., Dijkstra, H. & Brinkhuis, H. Emplacement of Antarctic ice sheet mass affects circumpolar ocean flow. Glob. Planet. Change 118, 16–24 (2014).
Pierce, E. L. et al. Evidence for a dynamic East Antarctic ice sheet during the mid-Miocene climate transition. Earth Planet. Sci. Lett. 478, 1–13 (2017).
Ivins, E. R., Raymond, C. A. & James, T. S. The influence of 5000 year-old and younger glacial mass variability on present-day crustal rebound in the Antarctic Peninsula. Earth Planets Space 52, 1023–1029 (2000).
Barnhoorn, A., van der Wal, W., Vermeersen, B. L. A. & Drury, M. R. Lateral, radial, and temporal variations in upper mantle viscosity and rheology under Scandinavia. Geochem. Geophy. Geosy. 12, Q01007 (2011).
Bentley, M. J. et al. A community-based geological reconstruction of Antarctic Ice Sheet deglaciation since the Last Glacial Maximum. Quat. Sci. Rev. 100, 1–9 (2014).
Bradley, S. L., Hindmarsh, R. C. A., Whitehouse, P. L., Bentley, M. J. & King, M. A. Low post-glacial rebound rates in the Weddell Sea due to Late Holocene ice-sheet readvance. Earth Planet. Sci. Lett. 413, 79–89 (2015).
Conway, H. et al. Switch of flow direction in an Antarctic ice stream. Nature 419, 465–467 (2002).
Vaughan, D. G., Corr, H. F. J., Smith, A. M., Pritchard, H. D. & Shepherd, A. Flow-switching and water piracy between Rutford Ice Stream and Carlson Inlet, West Antarctica. J. Glaciol. 54, 41–48 (2008).
Jenkins, A. et al. Observations beneath Pine Island Glacier in West Antarctica and implications for its retreat. Nat. Geosci. 3, 468 (2010).
Thomas, E. R., Bracegirdle, T. J., Turner, J. & Wolff, E. W. A 308 year record of climate variability in West Antarctica. Geophys. Res. Lett. 40, 2013GL057782 (2013).
Thomas, E. R., Marshall, G. J. & McConnell, J. R. A doubling in snow accumulation in the western Antarctic Peninsula since 1850. Geophys. Res. Lett. 35, L01706 (2008).
Siegert, M., Ross, N., Corr, H., Kingslake, J. & Hindmarsh, R. Late Holocene ice-flow reconfiguration in the Weddell Sea sector of West Antarctica. Quat. Sci. Rev. 78, 98–107 (2013).
Kingslake, J., Martin, C., Arthern, R. J., Corr, H. F. J. & King, E. C. Ice-flow reorganization in West Antarctica 2.5kyr ago dated using radar-derived englacial flow velocities. Geophys. Res. Lett. 43, 9103–9112 (2016).
Nield, G. A., Whitehouse, P. L., King, M. A. & Clarke, P. J. Glacial Isostatic Adjustment in response to changing Late Holocene behaviour of ice streams on the Siple Coast, West Antarctica. Geophys. J. Int. 205, 1–21 (2016).
Mackintosh, A. N. et al. Retreat history of the East Antarctic Ice Sheet since the Last Glacial Maximum. Quat. Sci. Rev. 100, 10–30 (2014).
Fieber, K. D. et al. Rigorous 3D change determination in Antarctic Peninsula glaciers from stereo WorldView-2 and archival aerial imagery. Remote Sens Environ. 205, 18–31 (2018).
Hodgson, D. A. et al. Rapid early Holocene sea-level rise in Prydz Bay, East Antarctica. Glob. Planet Change 139, 128–140 (2016).
Makinen, J., Amalvict, M., Shibuya, K. & Fukuda, Y. Absolute gravimetry in Antarctica: Status and prospects. J. Geodyn. 43, 339–357 (2007).
Auriac, A. et al. Iceland rising: Solid Earth response to ice retreat inferred from satellite radar interferometry and visocelastic modeling. J. Geophys. Res. Solid Earth 118, 1331–1344 (2013).
Wahr, J. et al. The use of GPS horizontals for loading studies, with applications to northern California and southeast Greenland. J. Geophys. Res. 118, 1–12 (2013).
King, M. A., Whitehouse, P. L. & van der Wal, W. Incomplete separability of Antarctic plate rotation from glacial isostatic adjustment deformation within geodetic observations. Geophys J. Int. 204, 324–330 (2016).
King, M. A. & Santamaría-Gómez, A. Ongoing deformation of Antarctica following recent Great Earthquakes. Geophys. Res. Lett. 43, 1918–1927 (2016).
Briggs, R., Pollard, D. & Tarasov, L. A glacial systems model configured for large ensemble analysis of Antarctic deglaciation. Cryosphere 7, 1949–1970 (2013).
Briggs, R. D. & Tarasov, L. How to evaluate model-derived deglaciation chronologies: a case study using Antarctica. Quat. Sci. Rev. 63, 109–127 (2013).
Melini, D., Gegout, P., King, M., Marzeion, B. & Spada, G. On the Rebound: Modeling Earth’s Ever-Changing Shape Eos 96, https://doi.org/10.1029/2015EO033387 (2015).
Agnew, D. C. NLOADF: A program for computing ocean-tide loading. J. Geophys. Res. 102, 5109–5110 (1997).
Spada, G. & Stocchi, P. SELEN: A Fortran 90 program for solving the “sea-level equation”. Comput. Geosci.-Uk 33, 538–562 (2007).
Moresi, L. et al. CitcomS v3.3.1 [software]. Computational Infrastructure for Geodynamics. https://geodynamics.org/cig/software/citcoms/ (2014).
Kendall, R. A., Mitrovica, J. X. & Milne, G. A. On post-glacial sea level - II. Numerical formulation and comparative results on spherically symmetric models. Geophys. J. Int. 161, 679–706 (2005).
Purcell, A., Tregoning, P. & Dehecq, A. An assessment of the ICE6G_C(VM5a) glacial isostatic adjustmentmodel. J. Geophys. Res-Solid Earth 121, 3939–3950 (2016).
Peltier, W. R., Argus, D. F. & Drummond, R. Comment on the paper by Purcell et al (2016) entitled “An Assessment of the ICE-6G_C (VM5a) glacial isostatic adjustment model”. J. Geophys. Res. Solid Earth 122, https://doi.org/10.1002/2016JB013844 (2017).
Purcell, A., Tregoning, P. & Dehecq, A. 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”. J. Geophys. Res. Sol. Earth 123, 2029–2032 (2018).
Hay, C. C. et al. Sea-level fingerprints in a region of complex Earth structure: the case of WAIS. J. Clim. 30, 1881–1892 (2017).
Milne, G. A. et al. The influence of lateral Earth structure on glacial isostatic adjustment in Greenland. Geophys J. Int 214, 1252–1266 (2018).
Briggs, R. D., Pollard, D. & Tarasov, L. A data-constrained large ensemble analysis of Antarctic evolution since the Eemian. Quat. Sci. Rev. 103, 91–115 (2014).
Sasgen, I. et al. Antarctic ice-mass balance 2003 to 2012: regional reanalysis of GRACE satellite gravimetry measurements with improved estimate of glacial-isostatic adjustment based on GPS uplift rates. Cryosphere 7, 1499–1512 (2013).
Kustowski, B., Ekstrom, G. & Dziewonski, A. M. Anisotropic shear-wave velocity structure of the Earth’s mantle: A global model. J. Geophys Res-Sol. Ea 113, B06306 (2008).
Fretwell, P. et al. Bedmap2: improved ice bed, surface and thickness datasets for Antarctica. Cryosphere 7, 375–393 (2013).
Peltier, W. R. Impulse Response of a Maxwell Earth. Rev. Geophys. 12, 649–669 (1974).
Spada, G. et al. A benchmark study for glacial isostatic adjustment codes. Geophys. J. Int. 185, 106–132 (2011).
Ranalli, G. & Fischer, B. Diffusion Creep, Dislocation Creep, and Mantle Rheology. Phys. Earth Planet Int. 34, 77–84 (1984).
Karato, S. & Wu, P. Rheology of the Upper Mantle - a Synthesis. Science 260, 771–778 (1993).
van der Wal, W. et al. Glacial isostatic adjustment model with composite 3-D Earth rheology for Fennoscandia. Geophys. J. Int. 194, 61–77 (2013).
Mainprice, D., Barruol, G. & Ben Ismail, W. in Earth’s Deep Interior: Mineral Physics and Tomography from the Atomic to the Global Scale (eds. Karato, S. Liebermann, R.C., Masters, G. & Stixrude, L.) 237–264 (American Geophysical Union, Washington D.C., USA 2000).
Beghein, C., Trampert, J. & van Heijst, H. J. Radial anisotropy in seismic reference models of the mantle. J. Geophys. Res. Solid Earth 111, B02303 (2006).
Gasperini, P., Dal Forno, G. & Boschi, E. Linear or non-linear rheology in the Earth’s mantle: the prevalence of power-law creep in the postglacial isostatic readjustment of Laurentia. Geophys. J. Int. 157, 1297–1302 (2004).
van der Wal, W., Wu, P., Wang, H. S. & Sideris, M. G. Sea levels and uplift rate from composite rheology in glacial isostatic adjustment modeling. J. Geodyn. 50, 38–48 (2010).
Burgmann, R. & Dresen, G. Rheology of the lower crust and upper mantle: evidence from rock mechanics, geodesy, and field observations. Annu. Rev. Earth Planet Sci. 36, 531–567 (2008).
Pollitz, F. F. Transient rheology of the upper mantle beneath central Alaska inferred from the crustal velocity field following the 2002 Denali earthquake. J. Geophys. Res. Solid Earth 110, B08407 (2005).
Yuen, D. A., Sabadini, R. C. A., Gasperini, P. & Boschi, E. On transient rheology and glacial isostasy. J. Geophys. Res. 91, 1420–1438 (1986).
Acknowledgements
We thank the following people for providing model output shown in Fig. 2: Erik Ivins, Dick Peltier, Lev Tarasov, Geruo A, Wouter van der Wal, Ingo Sasgen, Brian Gunter, and Alba Martín-Español. P.L.W. was supported by a UK Natural Environment Research Council (NERC) Independent Research Fellowship (NE/K009958/1). N.G. was supported by National Sciences and Engineering Research Council (NSERC), the Canada Research Chairs Program and McGill University. M.A.K. was supported by the Australian Research Council Special Research Initiative for Antarctic Gateway Partnership (Project ID SR140300001). This article is a contribution to the SCAR SERCE program.
Author information
Authors and Affiliations
Contributions
P.L.W. drafted the structure of the article, wrote the introductory and summary text, and the section on long-term landscape evolution. N.G. wrote the section on feedbacks between GIA and ice dynamics. M.A.K. wrote the sections on current knowledge of glacial isostasy and outstanding problems. D.A.W. wrote the section on rheology of the solid Earth. All authors commented on the final version of the article.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Journal peer review information: Nature Communications thanks the anonymous reviewers for their contribution to the peer review of this work. Peer reviewer reports are available.
Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Whitehouse, P.L., Gomez, N., King, M.A. et al. Solid Earth change and the evolution of the Antarctic Ice Sheet. Nat Commun 10, 503 (2019). https://doi.org/10.1038/s41467-018-08068-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-018-08068-y
- Springer Nature Limited
This article is cited by
-
Ocean cavity regime shift reversed West Antarctic grounding line retreat in the late Holocene
Nature Communications (2024)
-
A thicker Antarctic ice stream during the mid-Pliocene warm period
Communications Earth & Environment (2023)
-
Overshooting the critical threshold for the Greenland ice sheet
Nature (2023)
-
A topographic hinge-zone divides coastal and inland ice dynamic regimes in East Antarctica
Communications Earth & Environment (2023)
-
Mass Balances of the Antarctic and Greenland Ice Sheets Monitored from Space
Surveys in Geophysics (2023)