Collapse of the West Antarctic Ice Sheet after local destabilization of the Amundsen Basin

Edited by Michael Oppenheimer, Princeton University, Princeton, NJ, and accepted by the Editorial Board October 6, 2015 (received for review June 25, 2015)
November 2, 2015
112 (46) 14191-14196

Significance

The Antarctic Ice Sheet is losing mass at an accelerating rate, and playing a more important role in terms of global sea-level rise. The Amundsen Sea sector of West Antarctica has most likely been destabilized. Although previous numerical modeling studies examined the short-term future evolution of this region, here we take the next step and simulate the long-term evolution of the whole West Antarctic Ice Sheet. Our results show that if the Amundsen Sea sector is destabilized, then the entire marine ice sheet will discharge into the ocean, causing a global sea-level rise of about 3 m. We thus might be witnessing the beginning of a period of self-sustained ice discharge from West Antarctica that requires long-term global adaptation of coastal protection.

Abstract

The future evolution of the Antarctic Ice Sheet represents the largest uncertainty in sea-level projections of this and upcoming centuries. Recently, satellite observations and high-resolution simulations have suggested the initiation of an ice-sheet instability in the Amundsen Sea sector of West Antarctica, caused by the last decades’ enhanced basal ice-shelf melting. Whether this localized destabilization will yield a full discharge of marine ice from West Antarctica, associated with a global sea-level rise of more than 3 m, or whether the ice loss is limited by ice dynamics and topographic features, is unclear. Here we show that in the Parallel Ice Sheet Model, a local destabilization causes a complete disintegration of the marine ice in West Antarctica. In our simulations, at 5-km horizontal resolution, the region disequilibrates after 60 y of currently observed melt rates. Thereafter, the marine ice-sheet instability fully unfolds and is not halted by topographic features. In fact, the ice loss in Amundsen Sea sector shifts the catchment's ice divide toward the Filchner–Ronne and Ross ice shelves, which initiates grounding-line retreat there. Our simulations suggest that if a destabilization of Amundsen Sea sector has indeed been initiated, Antarctica will irrevocably contribute at least 3 m to global sea-level rise during the coming centuries to millennia.
The Antarctic Ice Sheet is losing mass at an accelerating rate, and thus increasingly contributes to global sea-level rise (1, 2). To what extent this is caused by anthropogenic warming of the atmosphere is unclear. The Amundsen Sea sector of West Antarctica accounts for the largest part of the ice loss, manifesting in a speed-up, thinning, and retreat of tributaries such as Pine Island Glacier and Thwaites Glacier during the last 4 decades (35). Large portions of the West Antarctic Ice Sheet (WAIS) are grounded on bedrock below sea level (6, 7). These so-called marine parts of the ice sheet hold an ice volume that would elevate global mean sea level by about 3.3 m (8). The bed below this marine ice is generally down-sloping in the inland direction. A grounding line (the line that separates the grounded ice sheet from the floating ice shelf) that is located on such bed has been shown to be potentially unstable (9, 10). The associated marine ice-sheet instability (i.e., self-sustained ice-sheet retreat) can be hindered by the buttressing of ice shelves, which has been investigated in different modeling frameworks (1114).
The recently inferred destabilization of the Amundsen Sea sector coincides with an increase in the volume and temperature of relatively warm circumpolar deep water that reaches into the ice-shelf cavities in the Amundsen Sea (15, 16). The resulting substantially enhanced sub-ice-shelf melting facilitated the retreat of the grounding line of several tributaries onto the potentially unstable down-sloping bed section. Complex, high-resolution regional modeling studies (17, 18) confined to the Pine Island and Thwaites catchments, respectively, concluded that the retreat of the ice is likely to continue in the near future. The findings are supported by observations that indicate ongoing destabilization (3), including the possibility of a developing marine ice-sheet instability (4). These state-of-the-art regional simulations need to be complemented by simulations of the long-term evolution of the region to deduce implications for the interior region of the WAIS.
Here we provide such simulations and investigate the response of the WAIS to a destabilization of the Amundsen Sea sector in ice-dynamical simulations with the Parallel Ice Sheet Model (PISM). Previous numerical simulations with PISM showed an ocean-triggered destabilization of a basin of marine ice in the Wilkes basin of East Antarctica (19), but no instability in the catchment of the Filchner–Ronne Ice Shelf (20), even for strong projected ocean warming (21). Simulations provided here cover the whole WAIS and allow for an examination of the long-term evolution of the WAIS and its associated sea-level contribution, respectively.

Results and Discussion

Our thermomechanically coupled, 3D, ice-sheet model (22, 23) (Methods) applies a superposition of the shallow-ice and the shallow-shelf approximation (24) of the stress balance that ensures a smooth transition between the different flow regimes of an ice sheet, ranging from bed-frozen inland ice to well-lubricated, fast-flowing ice streams and freely floating ice shelves. Stress transmission across the grounding line allows for the effect of ice-shelf buttressing. The grounding line is free to evolve, and a reversible grounding-line motion that is comparable to results from full-Stokes simulations was demonstrated in the Marine Ice Sheet Model Intercomparison Project for plan-view models for horizontal resolutions of 5 km and finer (25).

Initialization of the Simulations and Model Performance.

The result of this study is an if–then statement, saying that if the Amundsen Sea Sector is destabilized, then the entire marine part of West Antarctica will be discharged into the ocean. To this end, we spin up our model, using present-day initial ice and bed topography under present-day constant atmosphere conditions and 20th-century sub-ice-shelf melting. The model domain comprises the entire WAIS, the Antarctic Peninsula, and part of East Antarctica (Fig. 1). We obtain an ice sheet in equilibrium with stable ice margins comparable to present-day observations (Fig. 2), with a total ice volume that differs by −1.8% from observations. The largest deviations from the present-day grounding line position occur in the Siple Coast region, (i.e., the Ross Ice Shelf), where our modeled position is landward-biased by several 10 km (Figs. 1B and 2). Difficulties in reproducing the observed grounding-line position with our model arise particularly from the region's flat bed topography beneath the ice streams and from limitations in our parameterization of basal friction (Methods). In contrast, the modeled grounding line is in good agreement with observations for the Filchner–Ronne Ice Shelf (with the exception of some narrow troughs) and in the Amundsen Sea sector (Fig. 1 B and C), which is the relevant region for our perturbation experiments. Although we do not resolve the West Antarctic ice velocity field in full detail, the overall observed flow pattern of the WAIS, including the location of finely branched ice streams, is captured in our equilibrium simulation (compare Fig. 2 A and B). In the Amundsen Sea sector, both the Pine Island and Thwaites tributaries are captured by the model, although they are less separated from each other than observed. The same holds for the Siple Coast ice streams that feed the Ross Ice Shelf from the WAIS direction. The ice streams that originate in East Antarctica and terminate in the Ross Ice Shelf, being strongly confined by the Transantarctic Mountains, are captured individually by the model. Observed and modeled flow patterns match best for the Filchner–Ronne catchment, where most of the East and West Antarctic tributaries are clearly resolved in our simulation.
Fig. 1.
Collapse of the WAIS in response to 60 y of currently observed sub-ice-shelf melting in the Amundsen Sea. (A) Simulated initial steady-state grounding line and calving front before perturbation (dark blue contours) and at the end of perturbation (red contours) underlaid by bed topography (blue shading). Grounding-line positions during the collapse in 1,000-y time steps (black contours). Final steady-state grounded ice body in transparent gray with gray contours of surface elevation. Cross sections for ice and bed topography of B and C (yellow lines). (B and C) Transects are displayed from the Amundsen Sea Sector (Left) into the Ross Ice Shelf (B) and Filchner–Ronne Ice Shelf (C) (Right). Colors as in A, with additional profiles from present-day observations (49) in green. (Inset) Region outside the model domain (hatched) and area of A (blue rectangle).
Fig. 2.
(A) Observed present-day (50) and (B) modeled surface speed in equilibrium for the whole model domain (region outside model domain hatched). Observed (49) and modeled grounding line and calving front, respectively (black contours), and bathymetry (light blue shading). Gray contour in A delineates present-day catchment basin of the Amundsen Sea sector (51).

Destabilization of Amundsen Sea Sector.

The field of sub-ice-shelf melting we apply during model spin-up is obtained from a 20th-century simulation of the Finite Element Sea Ice-Ocean Model (Methods) and underestimates present-day sub-ice-shelf melt rates in the Amundsen Sea by an order of magnitude (21). To account for the currently observed enhanced ice-shelf melt that is considered to have destabilized the ice in the Amundsen Sea sector, we perturb our equilibrium state by scaling the regionally underestimated melt-rate field, making the applied melt rates compatible with present-day observations. The perturbation is confined to the Amundsen Sea and is applied as a pulse for periods of different duration (20–200 y), after which the original ocean conditions with melting strongly below the presently observed values are applied.
With the onset of our perturbation of the steady-state ice sheet, the simulated net mass loss from the Amundsen Sea sector increases from approximately zero to an average rate of ∼0.2 mm sea-level equivalent per year (sle/y; black curve in Fig. 3A, see lowermost blue curve for comparison with unperturbed control run). Measurements during the last 2 decades range from a minimum of 0.05 mm sle/y in the mid-1990s to a maximum of 0.32 mm sle/y in 2009 (26). Using the measured rates to extrapolate cumulative net ice loss (dashed lines that span gray area in Fig. 3A), we show that our modeled sea-level contribution from the perturbed region lies well within the range of observations. Simulated average effective melt rates below the Amundsen ice shelves during the perturbation period range from 30 to 33 m/y. Maximum melt rates of 75–145 m/y, in agreement with recent observations for Pine Island and Thwaites glaciers, occur close to the grounding line. For comparison, average observed rates are 16–33 m/y (27, 28), and maximum rates are observed to be above 100 m/y at the grounding line (29). The modeled ice sheet responds to the applied enhanced melting by inland-spreading ice thinning and grounding-line retreat (compare blue with red profiles in Fig. 1 B and C), consistent with the observed present-day situation in the Amundsen catchment basins (3).
Fig. 3.
Time evolution of the sea-level contribution after the beginning of the perturbation of the equilibrium ice sheet for (A) the first 250 y, covering the perturbation phase (black), and (B) the full simulation. Each colored line represents an individual simulation characterized by a different perturbation duration (end of perturbation after 0, 20, 40, 50, 60, 80, 100, 150, and 200 y, respectively, marked by vertical bars). Stable and unstable simulations in blue and red, respectively. Gray shading shows the range of cumulative sea-level contribution spanned when extrapolating rates of net ice loss from Amundsen Sea sector that were measured for five periods during the last two decades [dashed lines (26)]. Horizontal dotted gray line in B gives the calculated sea level equivalent for a complete drainage of the Amundsen catchment basins (3).

Ice Sheet Response after Destabilization.

After cessation of the perturbation, the original 20th-century melt rates are applied and ice-discharge rates slightly level off (Fig. 3A) during a phase of local ice adjustment. For simulations with a perturbation duration of at least 60 y, an abrupt increase in ice loss then marks the onset of the self-sustained ice loss (Fig. 3B). In these destabilized simulations, the major part of West Antarctica's marine ice disintegrates: The major part of the catchment of the Pine Island and Thwaites glaciers (gray contour in Fig. 2A) is drained in less than 3,000 y (Fig. 4B). The grounding line retreats rapidly inland, following the negative bed gradient into the deep Byrd Subglacial Basin (Figs. 1 and 4 and Movie S1). Further retreat into the interior of the WAIS comes to a halt at the steep bed of the Ellsworth Subglacial Highlands but continues laterally into the marine parts of Ellsworth Land and Marie Byrd Land (Fig. 4C).
Fig. 4.
Snapshots of modeled surface speed at (A) 0, (B) 3,000, (C) 6,000, and (D) 15,000 y after onset of the 60-y perturbation (zoom into region depicted by rectangle in Fig. 2B, same region as in Fig. 1). Areas of very small ice velocities mark the ice-divide location that shifts from Amundsen into the Ross and Filchner–Ronne basins during ice-sheet retreat (compare A, B, and C). Comparison of B and C reveals grounding-line destabilization from the landward direction in the catchment basin of the Ross Ice Shelf.
Today, these areas act as reservoirs that feed the Filchner–Ronne and Ross ice shelves, respectively. In the course of the destabilization, the widespread thinning causes the ice divides that mark the catchment areas of Pine Island and Thwaites glaciers at present to migrate into the drainage basins of Filchner–Ronne and Ross ice shelves, respectively (indicated by blue regions of small surface speed in Fig. 4 AC). The thinning signal reaches well beyond the ice divides and eventually arrives at the grounding lines of these basins. This grounding-line perturbation from the landward direction induces ice-sheet destabilization, also from the Filchner–Ronne and Ross directions (black contours in Fig. 1A, also compare Fig. 4 B and C). At some point, the rate of grounding-line retreat is similar for both landward and ocean directions (Fig. 1B). The underlying mechanism has been introduced previously, using a 2D conceptual framework (30), but here it also manifests in a complex, real-world domain that includes the important effect of buttressing.
The remnants of the ice sheet equilibrate after ∼13 ky, and total sea-level contribution by that time has reached about 3 m (Fig. 3B). In the new equilibrium state, an area of around 0.96· 106 km2 of originally grounded ice several kilometers in thickness is covered by an extensive ice shelf that is less than 500 m thick on average. The trans-Antarctic seaways between the previously unconnected Amundsen, Weddell, and Ross seas are opened (Fig. 4D), which is a situation that probably prevailed the last time during the Pleistocene (31).
Consistent with previous simulations with models of higher resolution (17, 18), our results show an earlier onset of strong ice loss with increasing forcing duration (Fig. 3B). Independent of the associated mass loss during perturbation, the eventual total sea-level rise contribution of about 3 m is the same for all destabilized simulation runs (Fig. 5, red circles). A difference of ∼1,500 Gt (4.2 mm sle) of ice loss during perturbation between the experiments with 50-y and 60-y long forcing, respectively, separates the stable from the unstable regime in our simulations (Fig. 5, gray shading). Our model's response to basal perturbations has been shown to be faster by a factor of around two compared with full-Stokes simulations (32). Hence, the speed of the WAIS collapse in our experiments might be overestimated.
Fig. 5.
Long-term sea-level contribution as a function of net ice loss from the Amundsen Sea sector during perturbation. Gray shading marks the critical threshold of ice loss between stabilization (blue) and destabilization (red) of the Amundsen Sea sector, followed by disintegration of the WAIS, corresponding to a sea-level rise contribution of almost 3 m.
The constant present-day atmosphere conditions applied during ice-sheet relaxation neglect the effect of rising temperatures and an associated increase in snowfall that are expected in a warming future climate (33). Although increased ice-sheet mass gain through snowfall may counteract imbalance, it can also lead to an increase in ice discharge resulting from the difference in surface elevation change between grounded and floating regions (34). The ice softening through higher surface temperatures tends to intensify ice flow and associated discharge. The question of to which extent these effects might force or inhibit a marine ice-sheet instability remains unanswered by our study. Our approach to apply subshelf melt rates from ocean model simulations unidirectionally to the ice model does not allow for ice–ocean interaction. Although fully coupled ice-sheet climate models are important (35, 36), they represent a compromise between model resolution and coupling degree. Here we account for the main effect of changing ice-shelf thickness on melt rates by applying a pressure-dependent melt adaptation in our model (Methods). Our ice-sheet model applies physically motivated approximations to the stress balance that determines ice-sheet dynamics. While accounting for the relevant stress components, the model allows computationally efficient simulations of the large region covered by the WAIS on a relatively high horizontal resolution over millennial timescales. Other ice-sheet models approximate the stress balance to lesser degree (37), or even solve the full-Stokes stress balance (38), and incorporate elaborated numerical methods such as grid refinement (39), enabling a more accurate simulation of ice-sheet evolution. The application of these models demands much more computational capacity, and to our knowledge, a study such as the one presented here would not be feasible to date.

Conclusions

Using a state-of-the-art continental ice-sheet model, we complement earlier simulations of higher resolution by investigating the long-term response of the WAIS to a destabilization of Amundsen Sea sector. Starting from a simulated equilibrium ice sheet close to observations, we apply a local perturbation to the Amundsen region that is constrained by present-day measurements of net ice loss and observations of sub-ice-shelf melt rates. After destabilization of the Amundsen Sea sector, we find that the marine ice-sheet instability unfolds and is not halted by topographic features. We find that ice thinning spreads far inland and shifts the ice divide of the Amundsen catchment toward the Filchner–Ronne and Ross ice shelves. Consistent with an earlier conceptual study (30), this initiates grounding-line retreat there. Thus, our simulations serve as an example for a scenario in which a destabilization of Amundsen Sea sector will not remain limited to its catchment basins but will lead to a complete disintegration of West Antarctica's marine ice on a millennial timescale. The currently observed retreat in West Antarctica hence might mark the beginning of a millennial period of self-sustained ice discharge from West Antarctica and require long-term global adaptation of coastal protection, such as the building or rebuilding or raising of dykes, the construction of seawalls, or the realization of land fills in the hinterland. For more details on the effects and adaptation strategies to sea-level rise, see the Intergovernmental Panel on Climate Change Working Group II’s latest assessment report (40, 41).
We emphasize that our findings are based on a single realization of an ice-sheet model that applies approximations to the ice dynamics. Thus, further modeling studies using different approaches (including perturbed-physics ensembles; the investigation of different scenarios of surface accumulation and temperature forcing; coupling between atmosphere, ice, and ocean models; or the application of models of higher numerical order) are much needed for comparison, setting our findings into context.

Methods

Ice-Sheet Model.

We use PISM (22, 23) to carry out regional simulations of the WAIS at 5-km horizontal resolution. The model applies a superposition of the shallow-ice approximation and the shallow-shelf approximation. This hybrid scheme ensures a smooth transition between different ice-sheet flow regimes and allows for stress transmission across the grounding line. We adjust for unmodeled anisotropy, as well as for other uncertainties relating to unobserved variations in viscosity and basal resistance in our sheet-shelf system, by using an enhancement factor of 4.5 for the shallow-ice approximation and 0.8 for the shallow-shelf approximation, respectively. A linear interpolation of the freely evolving grounding line, and accordingly interpolated basal friction, enable realistic grounding-line motion also at medium to low horizontal resolution (32). Basal friction is calculated using a nonlinear Weertman-type sliding law (24) with a sliding exponent of 1/3 in combination with a Mohr–Coulomb model for plastic till (42, 43) that accounts for the effect of evolving ice thickness and the associated change in overburden pressure on the basal till (44). The till friction angle is parameterized with bed elevation, as in ref. 42, and the basal pore water pressure is limited to a maximum fraction of 0.93 of the overburden pressure. This friction scheme ensures a continuous transition from quasi-nonslip regimes in elevated regions to the marine areas where basal resistance is low. We use a kinematic first-order calving law (45, 46), with a prescribed proportionality factor of 2·106 m⋅y and a minimum ice thickness at the calving front of 200 m. The physically motivated calving law takes into account the eigenvalues of the horizontal strain-rate tensor, allowing for geographically confined ice shelves and dynamic calving-front positions. A physical stress boundary condition is applied at the calving front to close the equations of the shallow-shelf approximation (23).

Boundary Conditions.

The landward boundary of our model domain is chosen along present-day East Antarctic ice divides, and hence far away from the marine parts of the WAIS. Along this fixed boundary, ice velocities are set to zero, whereas the coastal ice margin is free to evolve.
We spin up the model for 50,000 y under present-day time-constant atmosphere and 20th-century ocean boundary conditions, respectively, to obtain an ice sheet in equilibrium that serves as initial state for our perturbation experiments. Relative drift in total ice volume is less than 10−3 % y−1 after 50,000 y.
Subshelf melt rates are applied as an ocean boundary condition in our simulations, taken from the Finite Element Sea Ice-Ocean Model (FESOM) (21). FESOM is a hydrostatic, primitive-equation ocean model that is coupled to a dynamic-thermodynamic sea ice model. The ocean model includes an ice-shelf component that uses a three-equation system for the computation of temperature and salinity in the boundary layer between ice and ocean and the melt rate at the ice-shelf base. The global unstructured grid has a horizontal resolution of 10 km along the Antarctic coast, 7 km under the larger ice shelves in the Ross and Weddell seas, and 4 km under the small ice shelves in the Amundsen and Bellingshausen seas. Ice-shelf draft, cavity geometry, and ocean bathymetry used in FESOM consider data from many of the most recent surveys of the Antarctic continental shelf (47). The ocean model is initialized using temperature and salinity fields from National Oceanic and Atmospheric Administration’s World Ocean Atlas 2001 and is driven by atmospheric conditions from the HadCM3 global climate model (referred to as the HadCM3-20C experiment in ref. 21). Our equilibrium simulation is forced by time-averaged FESOM data for the period 1970–1999, representing 20th-century conditions. Compared with observations, sub-ice-shelf melting from FESOM is reasonably captured for Filchner–Ronne and Ross ice shelves, but significantly underestimated beneath the ice shelves in the Amundsen Sea (figure 10 in ref. 21). For this region, average FESOM melt rates are around 3 m/y, whereas present-day average melt rates are observed to be one order of magnitude larger (15, 29). We thus scale the field of melt rates below the Amundsen ice shelves toward values that are consistent with current observations to represent the observed enhanced melting below these ice shelves, perturbing the equilibrium ice sheet locally. The perturbation is applied for lengths of 20, 40, 50, 60, 80, 100, 150, and 200 y (Fig. 3). After the end of the perturbation pulse, the original (but weaker than observed) 20th-century melt rates are applied, and simulations are integrated for several thousand years until a stable grounding-line position is reached.
In FESOM simulations, the melt rates are computed for fixed present-day ice-shelf geometry and grounding-line position. We hence filled in missing ocean data in regions covered at present by land ice by a simple diffusive algorithm into regions below sea level, leading to melt rates further inland that are similar to the respective ice-shelf cavity melt rates. To account for local ice-shelf thickness change, we apply a pressure adaption of the ocean melt rates in PISM (see supplementary material in ref. 20). Ocean melting is also applied at the grid boxes below sea level, with grounded ice adjacent to the ocean through a buoyancy-dependent basal melt interpolation scheme. Prescribed fields of present-day surface mass balance and temperature obtained from the ALBMAPv1 dataset (48) act as atmospheric boundary conditions in all simulations. Initial bed and ice topography are obtained from the recent Bedmap2 data set (49).

Acknowledgments

We thank two anonymous reviewers for their valuable comments. We are grateful to Ralph Timmermann for providing the Finite Element Sea Ice-Ocean Model (FESOM) data. The research leading to these results has received funding from the European Union Seventh Framework Programme FP7/2007-2013 under Grant Agreement 603864. The study was financially supported by the German Environmental Foundation (DBU). Model development for the Parallel Ice Sheet Model (PISM) at the University of Alaska, Fairbanks, was supported by NASA Grants NNX13AM16G and NNX13AK27G.

Supporting Information

Supporting Information (PDF)
Supporting Information
pnas.1512482112.sm01.avi

References

1
M Mcmillan, et al., Increased ice losses from Antarctica detected by CryoSat-2. Geophys Res Lett 41, 3899–3905 (2014).
2
E Rignot, I Velicogna, M Van den Broeke, A Monaghan, J Lenaerts, Acceleration of the contribution of the Greenland and Antarctic ice sheets to sea level rise. Geophys Res Lett 38, L05503 (2011).
3
E Rignot, J Mouginot, Widespread, rapid grounding line retreat of Pine Island, Thwaites, Smith, and Kohler glaciers, West Antarctica, from 1992 to 2011. Geophys Res Lett 41, 3502–3509 (2014).
4
J Mouginot, E Rignot, B Scheuchl, Sustained increase in ice discharge from the Amundsen Sea Embayment, West Antarctica, from 1973 to 2013. Geophys Res 41, 1576–1584 (2014).
5
A Shepherd, et al., A reconciled estimate of ice sheet mass balance. 338(6111):1183–1189. (2012).
6
DG Vaughan, et al., New boundary conditions for the West Antarctic ice sheet: Subglacial topography beneath Pine Island Glacier. Geophys Res Lett 33, L09501 (2006).
7
JW Holt, et al., New boundary conditions for the West Antarctic Ice Sheet: Subglacial topography of the Thwaites and Smith glacier catchments. Geophys Res Lett 33, L09502 (2006).
8
JL Bamber, REM Riva, BL Vermeersen, AM LeBrocq, Reassessment of the potential sea-level rise from a collapse of the West Antarctic Ice Sheet. Science 324, 901–903 (2009).
9
J Weertman, Stability of the junction of an ice sheet and an ice shelf. J Glaciol 13, 3–11 (1974).
10
C Schoof, Ice sheet grounding line dynamics: Steady states, stability, and hysteresis. J Geophys Res 112, F03S28 (2007).
11
D Docquier, D Pollard, F Pattyn, Thwaites Glacier grounding-line retreat: Influence of width and buttressing parameterizations. J Glaciol 60, 305–313 (2014).
12
AP Wright, et al., Sensitivity of the Weddell Sea sector ice streams to sub-shelf melting and surface accumulation. Cryosph 8, 2119–2134 (2014).
13
GH Gudmundsson, Ice-shelf buttressing and the stability of marine ice sheets. Cryosph 7, 647–655 (2013).
14
D Goldberg, DM Holland, C Schoof, Grounding line movement and ice shelf buttressing in marine ice sheets. J Geophys Res 114, F04026 (2009).
15
A Jenkins, et al., Observations beneath Pine Island Glacier in West Antarctica and implications for its retreat. Nat Geosci 3, 468–472 (2010).
16
AJ Payne, A Vieli, AP Shepherd, DJ Wingham, E Rignot, Recent dramatic thinning of largest West Antarctic ice stream triggered by oceans. Geophys Res Lett 31, L021284 (2004).
17
L Favier, et al., Retreat of Pine Island Glacier controlled by marine ice-sheet instability. Nat Clim Chang 4, 117–121 (2014).
18
I Joughin, BE Smith, B Medley, Marine ice sheet collapse potentially under way for the Thwaites Glacier Basin, West Antarctica. Science 344, 735–738 (2014).
19
M Mengel, A Levermann, Ice plug prevents irreversible discharge from East Antarctica. Nat Clim Chang 4, 451–455 (2014).
20
M Mengel, J Feldmann, A Levermann, Linear sea-level response to abrupt ocean warming of major West Antarctic ice basin. Nat Clim Chang, October 5, 2015).
21
R Timmermann, HH Hellmer, Southern Ocean warming and increased ice shelf basal melting in the twenty-first and twenty-second centuries based on coupled ice-ocean finite-element modelling. Ocean Dyn 63, 1011–1026 (2013).
22
E Bueler, J Brown, Shallow shelf approximation as a “sliding law” in a thermomechanically coupled ice sheet model. J Geophys Res Solid Earth 114, F03008 (2009).
23
R Winkelmann, et al., The Potsdam Parallel Ice Sheet Model (PISM-PIK) – Part 1: Model description. Cryosph 5, 715–726 (2011).
24
R Greve, H Blatter Dynamics of Ice Sheets and Glaciers (Springer, Berlin, 2009).
25
F Pattyn, et al., Grounding-line migration in plan-view marine ice-sheet models: Results of the ice2sea MISMIP3d intercomparison. J Glaciol 59, 410–422 (2013).
26
B Medley, et al., Constraining the recent mass balance of Pine Island and Thwaites glaciers, West Antarctica, with airborne observations of snow accumulation. Cryosphere 8, 1375–1392 (2014).
27
E Rignot, S Jacobs, J Mouginot, B Scheuchl, Ice-shelf melting around Antarctica. Science 341, 266–270 (2013).
28
SS Jacobs, A Jenkins, CF Giulivi, P Dutrieux, Stronger ocean circulation and increased melting under Pine Island Glacier ice shelf. Nat Geosci 4, 519–523 (2011).
29
P Dutrieux, et al., Pine Island glacier ice shelf melt distributed at kilometre scales. Cryosph 7, 1543–1555 (2013).
30
J Feldmann, A Levermann, Interaction of marine ice-sheet instabilities in two drainage basins: Simple scaling of geometry and transition time. Cryosph 9, 631–645 (2015).
31
DK Barnes, CD Hillenbrand, Faunal evidence for a late quaternary trans-Antarctic seaway. Glob Change Biol 16, 3297–3303 (2010).
32
J Feldmann, T Albrecht, C Khroulev, F Pattyn, A Levermann, Resolution-dependent performance of grounding line motion in a shallow model compared with a full-Stokes model according to the MISMIP3d intercomparison. J Glaciol 60, 353–360 (2014).
33
K Frieler, et al., Consistent evidence of increasing Antarctic accumulation with warming. Nat Clim Chang 5, 348–352 (2015).
34
R Winkelmann, A Levermann, MA Martin, K Frieler, Increased future ice discharge from Antarctica owing to higher snowfall. Nature 492, 239–242 (2012).
35
DN Goldberg, et al., Investigation of land ice-ocean interaction with a fully coupled ice-ocean model: 1. Model description and behavior. J Geophys Res Earth Surf 117, F02037 (2012).
36
D Holland, D Holland, On the rocks: The challenges of predicting sea level rise. Eos 96. Avilable at https://eos.org/features/on-the-rocks-the-challenges-of-predicting-sea-level-rise. Accessed October 19, 2015. (2015).
37
E Larour, H Seroussi, M Morlighem, E Rignot, Continental scale, high order, high spatial resolution, ice sheet modeling using the Ice Sheet System Model (ISSM). J Geophys Res 117, F01022 (2012).
38
L Favier, O Gagliardini, G Durand, T Zwinger, A three-dimensional full Stokes model of the grounding line dynamics: Effect of a pinning point beneath the ice shelf. Cryosphere 6, 101–112 (2012).
39
SL Cornford, et al., Century-scale simulations of the response of the West Antarctic Ice Sheet to a warming climate. Cryosph 9, 1579–1600 (2015).
40
CB Field, et al. IPCC, 2014: Climate Change 2014: Impacts, Adaptation, and Vulnerability. Part A: Global and Sectoral Aspects. Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (Cambridge University Press, Cambridge, 2014).
41
VR Barros, et al. IPCC, 2014: Climate Change 2014: Impacts, Adaptation, and Vulnerability. Part B: Regional Aspects. Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (Cambridge University Press, Cambridge, 2014).
42
MA Martin, et al., The Potsdam Parallel Ice Sheet Model (PISM-PIK) – Part 2: Dynamic equilibrium simulation of the Antarctic ice sheet. Cryosph 5, 727–740 (2011).
43
KM Cuffey, WSB Paterson The Physics of Glaciers (Academic Press, Burlington, MA, 2010).
44
GKC Clarke, Subglacial processes. Annu Rev Earth Planet Sci 33, 247–276 (2005).
45
A Levermann, et al., Kinematic first-order calving law implies potential for abrupt ice-shelf retreat. Cryosph 6, 273–286 (2012).
46
T Albrecht, M Martin, M Haseloff, R Winkelmann, A Levermann, Parameterization for subgrid-scale motion of ice-shelf calving fronts. Cryosph 5, 35–44 (2011).
47
R Timmermann, et al., A consistent dataset of Antarctic ice sheet topography, cavity geometry, and global bathymetry. Earth Syst Sci Data 2, 261–273 (2010).
48
A Le Brocq, AJ Payne, A Vieli, An improved Antarctic dataset for high resolution numerical ice sheet models (ALBMAPv1). Earth Syst Sci Data 2, 247–260 (2010).
49
P Fretwell, et al., Bedmap2: Improved ice bed, surface and thickness datasets for Antarctica. Cryosph 7, 375–393 (2013).
50
E Rignot, J Mouginot, B Scheuchl, Ice flow of the Antarctic ice sheet. Science 333, 1427–1430 (2011).
51
HJ Zwally, MB Giovinetto, MA Beckley, JL Saba, Antarctic and Greenland Drainage Systems, GSFC Cryospheric Sciences Laboratory. Available at icesat4.gsfc.nasa.gov/cryo_data/ant_grn_drainage_systems.php. Accessed March 1, 2015. (2012).

Information & Authors

Information

Published in

Go to Proceedings of the National Academy of Sciences
Proceedings of the National Academy of Sciences
Vol. 112 | No. 46
November 17, 2015
PubMed: 26578762

Classifications

Submission history

Published online: November 2, 2015
Published in issue: November 17, 2015

Keywords

  1. West Antarctic Ice Sheet
  2. sea-level rise
  3. tipping point
  4. instability
  5. marine ice-sheet instability

Acknowledgments

We thank two anonymous reviewers for their valuable comments. We are grateful to Ralph Timmermann for providing the Finite Element Sea Ice-Ocean Model (FESOM) data. The research leading to these results has received funding from the European Union Seventh Framework Programme FP7/2007-2013 under Grant Agreement 603864. The study was financially supported by the German Environmental Foundation (DBU). Model development for the Parallel Ice Sheet Model (PISM) at the University of Alaska, Fairbanks, was supported by NASA Grants NNX13AM16G and NNX13AK27G.

Notes

This article is a PNAS Direct Submission. M.O. is a guest editor invited by the Editorial Board.

Authors

Affiliations

Johannes Feldmann
Potsdam Institute for Climate Impact Research, 14473 Potsdam, Germany;
Institute of Physics, Potsdam University, 14476 Potsdam, Germany
Anders Levermann1 [email protected]
Potsdam Institute for Climate Impact Research, 14473 Potsdam, Germany;
Institute of Physics, Potsdam University, 14476 Potsdam, Germany

Notes

1
To whom correspondence should be addressed. Email: [email protected].
Author contributions: A.L. designed research; J.F. performed research; J.F. and A.L. analyzed data; and J.F. and A.L. wrote the paper.

Competing Interests

The authors declare no conflict of interest.

Metrics & Citations

Metrics

Note: The article usage is presented with a three- to four-day delay and will update daily once available. Due to ths delay, usage data will not appear immediately following publication. Citation information is sourced from Crossref Cited-by service.


Citation statements




Altmetrics

Citations

If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.

Cited by

    Loading...

    View Options

    View options

    PDF format

    Download this article as a PDF file

    DOWNLOAD PDF

    Login options

    Check if you have access through your login credentials or your institution to get full access on this article.

    Personal login Institutional Login

    Recommend to a librarian

    Recommend PNAS to a Librarian

    Purchase options

    Purchase this article to access the full text.

    Single Article Purchase

    Collapse of the West Antarctic Ice Sheet after local destabilization of the Amundsen Basin
    Proceedings of the National Academy of Sciences
    • Vol. 112
    • No. 46
    • pp. 14103-E6410

    Media

    Figures

    Tables

    Other

    Share

    Share

    Share article link

    Share on social media