Full text

Turn on search term navigation

ARTICLE

Received 28 Apr 2015 | Accepted 15 Oct 2015 | Published 26 Nov 2015

DOI: 10.1038/ncomms9910 OPEN

Rapid Holocene thinning of an East Antarctic outlet glacier driven by marine ice sheet instability

R.S. Jones1,2, A.N. Mackintosh1,2, K.P. Norton2, N.R. Golledge1,3, C.J. Fogwill4, P.W. Kubik5,M. Christl5 & S.L. Greenwood6

Outlet glaciers grounded on a bed that deepens inland and extends below sea level are potentially vulnerable to marine ice sheet instability. This instability, which may lead to runaway ice loss, has been simulated in models, but its consequences have not been directly observed in geological records. Here we provide new surface-exposure ages from an outlet of the East Antarctic Ice Sheet that reveal rapid glacier thinning occurred approximately 7,000 years ago, in the absence of large environmental changes. Glacier thinning persisted for more than two and a half centuries, resulting in hundreds of metres of ice loss. Numerical simulations indicate that ice surface drawdown accelerated when the otherwise steadily retreating glacier encountered a bedrock trough. Together, the geological reconstruction and numerical simulations suggest that centennial-scale glacier thinning arose from unstable grounding line retreat. Capturing these instability processes in ice sheet models is important for predicting Antarcticas future contribution to sea level change.

1 Antarctic Research Centre, Victoria University of Wellington, Kelburn, Wellington 6012, New Zealand. 2 School of Geography, Environment and Earth Sciences, Victoria University of Wellington, Kelburn, Wellington 6012, New Zealand. 3 GNS Science, Avalon, Lower Hutt 5011, New Zealand. 4 Climate Change Research Centre, University of New South Wales, Sydney 2052, Australia. 5 Laboratory of Ion Beam Physics, Department of Physics, ETH Zurich, Otto Stern Weg 5, 8093 Zurich, Switzerland. 6 Department of Geological Sciences, Stockholm University, 10691 Stockholm, Sweden. Correspondence and requests for materials should be addressed to R.S.J. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 6:8910 | DOI: 10.1038/ncomms9910 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms9910

The East Antarctic Ice Sheet (EAIS) contains more than two-thirds of the total volume of Antarctic ice grounded below sea level1, which is potentially vulnerable to marine

ice sheet instability26. This instability occurs when initial grounding line retreat into deeper water leads to thicker ice at the grounding line that is closer to oatation, causing increased ice ux and glacier thinning, a positive feedback that results in further grounding line retreat and drawdown of ice5,7,8. Satellite observations suggest that parts of the EAIS are currently experiencing dynamic thinning, and outlet glacier retreat9,10. This thinning is most likely initiated at the bed and terminus11,12, and may ultimately lead to irreversible mass loss if local topographic conditions favour retreat beyond a stability threshold5,6. Current understanding of self-sustaining glacier retreat is based on the short-lived tidewater cycles of marine-terminating glaciers in temperate and sub-polar settings7,8,13,14. Geological data can be used to extend the record of ice sheet observations in Antarctica, providing rates, durations and magnitudes of ice surface lowering that may have resulted from instability processes in the period preceding the satellite era.

Here we report a new terrestrial record of ice sheet thinning from Mackay Glacier, an outlet of the EAIS (Fig. 1). Mackay Glacier is suitably located to investigate the effects of marine ice sheet instability because: (1) a tectonically and glacially formed overdeepening occurs immediately offshore of the glacier15;(2) Mackay Glacier is laterally conned, and can be simulated using a one-dimensional owline model; and (3) a relatively high concentration of offshore and onshore chronologies provide regional context, recording when Mackay Glacier became independent from grounded ice in the western Ross Sea during the last deglaciation16. Our record identies a period of abrupt glacier change during the mid-Holocene. Using empirical data to constrain a high-resolution numerical model, we link changes in glacier surface elevation to the most likely driver of ice sheet retreat. In our model simulations, retreat of the grounding line across a landward-deepening bed is accompanied onshore by centennial-scale glacier thinning of the same magnitude as that shown by the terrestrial surface-exposure dates. On this basis, we argue that our surface-exposure data provide direct evidence of rapid topographically induced ice sheet retreat in Antarctica.

ResultsChronology of ice surface lowering. Dating of glacial erratic clasts in altitudinal transects has the potential to record the past surface lowering of a glacier from a previously thicker conguration17. We targeted two nunataks in the lower reaches of Mackay Glacier to determine its thinning history from the Last Glacial Maximum (LGM) to the present day. Mt Suess/Gondola Ridge lies adjacent to the main ow path of Mackay Glacier, B714-km upstream from the present-day grounding line and 1525 km from the terminus. Present-day ice ow reaches velocities of 180 m per year18 alongside this exposed basaltic and granitic

bedrock. The steep, ice-abraded anks of Mt Suess and glacially smoothed and streamlined topography of Gondola Ridge provide evidence of thicker ice in the past. Sandstone and quartzite glacial erratics of the Beacon Supergroup are found resting on both weathered and striated bedrock, and extend up to 890 m above sea

40 W

20 W

20 E

40 E

50 S

50 S 50 S

East Antarctic Ice Sheet

West Antarctic Ice Sheet

Mackay

Glacier

4,082

0 2,000Ice surface elevation (m a.s.l.)

50 S

7,400 4,000 1,000 Bathymetry/bed elevation (m a.s.l.)

140 W 160 E

160 W 180

140 E

167 E

168 E

165 E

163 E

161 E

78 S

Explorers Cove

Marble Point

McMurdo Sound

Scott Coast

77 S

Mackay Gl.

Granite

Harbour

Cape Bird

77 S

Pre 9 kaPre ~7.8/7.5 ka ~7.5 ka~6.3 ka Present-day

Supporting evidence Surface-exposure dating transectsLineations

Grounding-zone wedge

Kilometres

166 E

164 E

162 E

161 E

160 E

100 km

Figure 1 | Study context. (a) Ice shelves (outlined in red) buttress parts of Antarctic ice sheets that are grounded below sea level. Mackay Glacier and the location of b is marked by a black box. (b) Reconstructed grounding line (solid lines) and ice shelf (dashed lines) retreat downstream of Mackay Glacier, supported by local chronological evidence and offshore geomorphology (Supplementary Table 1). Arrows denote inferred ow of past grounded ice. (c) Surface velocity18 of present-day Mackay Glacier catchment with starred sample transect locations at Mt Suess/Gondola Ridge (upper, mid-lower and lower; GRU, GRM and GRL) and Low Ridge (LR), adjacent to the main ow path. Inset shows the owline used in model simulations and surface contours.

77 S

GRU

GRM

GRU

RM GRL

G GRL

GRL

77 S

2 NATURE COMMUNICATIONS | 6:8910 | DOI: 10.1038/ncomms9910 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms9910 ARTICLE

level (m a.s.l.) on Mt Suess (B340 m above the proximal modern ice surface). Low Ridge presently anks the glacier terminus, 14.5 km downstream of the grounding line. Here numerous sandstone cobbles are found on smooth, striated granitic bedrock. We collected quartz-rich cobbles resting on glacially eroded bedrock, along four separate altitudinal transects (Supplementary Fig. 1). The Mt Suess/Gondola Ridge (upper) transect included 16 sampled erratics at 824587 m a.s.l. (26024 m above the modern ice surface), while an additional 18 samples were collected as part of two further transects at the downstream end of Gondola Ridge (mid-lower and lower). At Low Ridge, 10 samples were collected with an elevation range of 26462 m a.s.l. (2042 m above the modern ice surface).

Cosmogenic nuclide surface-exposure dating allowed us to establish a chronology of ice surface lowering (see Methods). Across all 4 transects, 44 exposure ages ranged from 22.32 thousand years to 24575 years BP (Supplementary Data 1; Fig. 2). Only four sample outliers were identied in the complete data set, by their stratigraphic positions within and between transects (Methods). The remaining 40 samples provide a direct constraint on outlet glacier thinning from the LGM to preindustrial times, with 29 of these samples revealing an episode of rapid Holocene thinning. The high coherence and sample density of our cosmogenic exposure ages with few outliers allows us to rene the timing of former ice surface lowering at Mackay Glacier using Bayesian age-elevation modelling (Methods).

The resulting chronology provides a near-complete record of outlet glacier thinning, with the Mt Suess/Gondola Ridge (upper) transect revealing ice surface lowering from B22 ka BP into the late Holocene (Fig. 3ac). Undated glacial erratics occur on Mt Suess up to 340 m above the modern ice surface (B890 m a.s.l.), indicating that Mackay Glacier may have been this thick during the LGM, or at another time in the recent geological past. Two dated samples from a surface elevation B260 m above the modern glacier date to the LGM (B22.32 and 18.33.5 ka BP). Subsequent ice surface lowering below this elevation appears to have been initially gradual, with just B30 m of thinning occurring before B8 ka BP. However, this part of the record lacks sufcient sample coverage to clarify surface elevation changes that

might have occurred between the LGM and early Holocene. The most notable ice surface elevation changes occurred in the mid-Holocene, when 480% of the total LGM-to-present thinning is recorded. This event is best preserved in the Mt Suess/Gondola Ridge (upper) transect, where rapid thinning is recorded from 230 to 50 m above the modern ice surface (180 m total lowering). The onset of rapid thinning occurred by B6.8 ka BP (see Methods;

Supplementary Fig. 2). The end of this thinning event is likely recorded by ve exposure ages at Gondola Ridge (lower), which is shown after B6 ka BP as a return to gradual rates of thinning.

Data from low ridge show that this more gradual thinning during the mid-late Holocene (a further B25 m of ice surface lowering)

persisted until B250 years ago, at least in the lower reaches of the glacier.

To estimate the rate and duration of rapid thinning at Mackay Glacier, we carried out linear regression analysis on data from the Mt Suess/Gondola Ridge (upper) and Low Ridge transects. While the peak rate of surface lowering during this episode and the corresponding duration are not possible to determine, this analysis provides the most probable time-averaged estimates for the whole episode of rapid thinning. Regression was applied randomly to exposure ages through a Monte Carlo simulation (see Methods; Supplementary Fig. 3). First, we assessed the rate of thinning, assuming a linear rate implied by the samples. Previous authors have used raw exposure ages and associated uncertainties for regression analysis to derive thinning rates19. Applied to Low Ridge, this regression analysis indicates thinning at 8.2358.8 cm per year (2 s). Due to the higher sample density at Mt Suess/

Gondola Ridge (upper, n 6), we were able to use the Bayesian-

modelled ages to estimate thinning rates, which we consider to provide the best estimates of the surface lowering history. Here regression analysis indicates thinning of 33.180.2 cm per year (2 s) and, as the rate likely varied during this period, we suggest that the upper end of this range better represents the peak of thinning (Fig. 3d). Second, the regression analysis also provides an estimate for the duration of this rapid thinning episode (Methods; Fig. 3e). At Mt Suess/Gondola Ridge (upper), the duration of this event is constrained by a dominant, normally distributed peak at B420 years (251731 years, 2 s), which also corresponds to its median and mean values. At Low Ridge, a wide and skewed distribution of possible durations are estimated with a best t of B300 years; however, the median value is consistent within 1 s of

Mt Suess/Gondola Ridge (upper) at B400 years. We consider these to be minimum estimates of duration as the onset of rapid thinning is not recorded at Mt Suess/Gondola Ridge (upper), while neither the start nor end is recorded at Low Ridge.

In summary, Bayesian age-elevation modelling and regression analysis of 10Be exposure ages demonstrate that an episode of accelerated ice surface lowering occurred at Mackay Glacier during the Holocene. This thinning was rapid both in the context of its LGM-to-present deglaciation history and of rapidly changing outlet glaciers observed in Antarctica today, such as Totten and Pine Island Glaciers whose margins are thinning at B40150 and B504200 cm per year, respectively9. At Mackay

Glacier, this episode of surface lowering persisted for at least 251 years (95% condence), providing a window into ice sheet behaviour that extends far beyond the period of modern satellite observations.

Absence of large climatic or oceanic changes in the Holocene. Substantial, rapid and prolonged ice surface lowering at Mackay Glacier occurred during a period when proxy data indicate an absence of large climatic or oceanic changes, especially in the context of LGM-to-present environmental changes (Fig. 3f). The Holocene climatic optimum in Antarctica was reached at B1012 ka BP (ref. 20), and ice core records show that

300

250

Elevation above modern ice (m)

200

150

100

0 0 5,000 10,000 15,000 20,000 25,000

Exposure age (years)

Figure 2 | Surface-exposure ages from altitudinal transects at Mackay Glacier. Displayed are the raw, unmodelled ages (1 s) with outliers included and identied (dashed circles). Surface lowering is recorded from B22 ka to B200 years ago with an episode of rapid thinning evident in all four transects, Mt Suess/Gondola Ridge upper (blue), Gondola Ridge mid-lower (purple) and lower (green), and Low Ridge (red).

NATURE COMMUNICATIONS | 6:8910 | DOI: 10.1038/ncomms9910 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms9910

300

400

100 cm per year Rates of rapid thinning

Elevation above modern ice (m) Frequency

250

50 cm per year

300

Elevation above

200

modern ice (m)

10 cm per year

1 cm per year

150

200

100

LR:

Mt Suess/Gondola Ridge (upper)

8359 cm per year

GRU: 3380 cm per year

10,000

20,000

30,000

5,000

5,500

6,000

6,500

7,000

7,500

Exposure age (years)

500

Duration of rapid thinning

400

Elevation above

modern ice (m)

GRU: 420 years

LR: 300 years

300

200

100

Gondola Ridge (lower)

2,000

4,000

6,000

8,000

10,000

500

1,000

1,500

Exposure age (years)

Duration (years)

250

Sea surface temperature (C)

200

Elevation above

modern ice (m)

Atmospheric temperature

anomaly (C)

Ice volume equiv. sea level

(m above present)

150

100

150

Low ridge

Climate forcing

2,000

4,000

6,000

8,000

10,000

5,000

10,000

15,000

20,000

25,000

Exposure age (years)

Years BP

Figure 3 | Age-modelled chronology of surface lowering at Mackay Glacier and regional climate forcing. Bayesian age-elevation modelled surface-exposure ages (black, 1 s) are shown for three of four transects (ac), linearly interpolated with 1 and 2 s uncertainty. Initial raw exposure ages (Fig. 2)

are in grey. Rapid thinning is recorded at B6,800-6,000 years before present. (d) Rapid surface lowering is estimated at Mt Suess/Gondola Ridge (upper; GRU) and Low Ridge (LR), assuming linear thinning rates between B6.8 and 6.0 ka (Supplementary Fig. 3). GRU provides a tighter estimate of33.1-80.2 cm per year. The best t (solid lines) and 95% condence bounds (dashed lines) are shown. At Low Ridge, the full extent of rapid thinning is uncertain; however, the end of this episode and a change to more gradual thinning may be recorded at Gondola Ridge (lower; green line). (e) Rapid thinning lasted for B420 years (251-731 years, 2 s), based on the higher-quality GRU chronology. As the onset of rapid thinning here may have been from the LGM surface elevation, we consider this a minimum estimate for the full duration of the episode. Median estimates are shown by solid lines, while 1 and 2 s are denoted by a shaded area (GRU) and dashed lines (LR). (f) This episode does not correspond to signicant increases in regional sea surface temperature26,

atmospheric temperature21 or sea level23, irrespective of exposure age uncertainty (grey shaded area; Supplementary Fig. 2).

4 NATURE COMMUNICATIONS | 6:8910 | DOI: 10.1038/ncomms9910 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms9910 ARTICLE

atmospheric temperatures cooled slightly from this peak to reach conditions similar to preindustrial climate by B7,000 years ago21,22. Rates of eustatic sea level rise reduced signicantly after

B8 ka BP (ref. 23), and the small amount of continuing sea level rise was most likely buffered at the Antarctic coastline by regional isostatic uplift24. This probably resulted in a lowering relative sea level in the region of Mackay Glacier during this time25. Although we lack proximal information about oceanic conditions, it also seems unlikely that very large changes in ocean temperature occurred at the time of rapid surface lowering, as far-eld ocean temperature reconstructions show little change during this period26; sea surface temperatures in the southwest Pacic Ocean also peaked before 10 ka BP, during the early Holocene climatic optimum. Therefore, while small-scale uctuations in any of these variables could have inuenced the precise timing of initial retreat8, external forcing from the atmosphere or ocean alone cannot explain the rapid nature of recorded ice surface lowering at Mackay Glacier.

Deglaciation in the western Ross Sea. Independent geological evidence constrains the timing and extent of past ice retreat in the western Ross Sea (Fig. 1b; Supplementary Table 1; Supplementary Note 1). Large-scale grounding line and ice shelf retreat to just north of Ross Island occurred by B10 ka BP (ref. 27), or perhaps slightly earlier (Supplementary Note 1). At B9.4 ka BP, Explorers Cove was still occupied by grounded ice28 and Hall et al.25 suggest grounded ice in McMurdo Sound had yet to retreat. A relative sea level curve produced from radiocarbon dating of raised beach deposits indicates that nal unloading of grounded ice adjacent to the Scott Coast occurred at B7.5 ka BP (ref. 25), which is supported by a radiocarbon dated bivalve in a sediment core collected off Cape Bird29. The remaining ice shelf disappeared between B7.5 and B6.3 ka BP, at which point Granite Harbour30, Gneiss Point and Marble Point28,31 had become open water. An ice shelf may have, however, still existed immediately south at Explorers Cove28.

The geomorphology on the present-day seabed helps infer the style and rate of past grounded ice retreat in and around McMurdo Sound during the Holocene15. Immediately, downstream of Mackay Glaciers terminus is the deep Mackay Sea Valley (below 800 m a.s.l.), which shallows and widens into

a trough and then the southern Drygalski Basin, north of McMurdo Sound (Supplementary Fig. 4). A series of grounding-zone wedges (GZWs) preserved in the outer part of the trough imply that the initial retreat of the Mackay Glacier grounding line was staggered.

Our surface-exposure chronology from Mackay Glacier indicates that rapid, uninterrupted thinning occurred over at least two and a half centuries. This most likely happened after the glacier retreated from the GZWs of the outer trough to the inner parts of Mackay Sea Valley, during a time when the large buttressing effect provided to Mackay Glacier by grounded ice, and possibly an ice shelf, in the Ross Sea was removed. Offshore chronologies indicate that these conditions were met by B7.5 and B6.3 ka BP, respectively, which is consistent with the timing of rapid thinning recorded in our onshore chronology (Supplementary Fig. 2).

Dynamic glacier thinning requires a reduction of resistive stresses at the bed, grounding line, or near the terminus7,12,32,33. Regional geological evidence indicates that the most probable processes that drove rapid thinning of Mackay Glacier were retreat over a reverse bed slope into Mackay Sea Valley and/or removal of a local ice shelf, providing a testable hypothesis.

Bathymetric proles of Mackay Sea Valley and trough reveal an overdeepening, which could have resulted in Mackay Glacier retreating as a consequence of marine ice sheet instability2. The

loss of buttressing from ungrounding at the terminus during periods of rapid grounding line retreat produces a dynamic adjustment where thinning is propagated upstream7,33,34. The conned lower reaches of Mackay Glacier during the mid-Holocene were most likely susceptible to such upstream propagation; today, surface velocities of 4150 m per year18 and probable abundant basal sliding35 help to quickly redistribute ice mass7. Abundant roche moutones and striae on Cuff Cape near the current terminus indicate that these warm-based conditions also existed in the recent past36.

Geomorphological evidence from the region suggests that Mackay Glacier may have also been inuenced by the loss of an adjacent ice shelf. Ice shelf presence may have acted to slow or halt retreat at Mackay Glacier during its earliest stages, even if the grounding line was positioned on a reverse bed slope5,11,37. Some GZWs are located on a reverse bed slope within the Mackay Glacier trough (Supplementary Fig. 4), possibly suggesting that a supporting ice shelf was present at this time. The ultimate loss of this ice shelf would have led to increased ice uxes over the Mackay Glacier grounding line38, possibly enabling a period of accelerated retreat and the rapid thinning recorded onshore.

Both of these scenarios involve rapid thinning of Mackay Glacier being driven by perturbations near its terminus, resulting in glacier retreat through an overdeepened trough.

Numerical modelling of glacier retreat and surface lowering. In an attempt to understand whether progressive retreat of Mackay Glacier through its overdeepened trough resulted in the rapid thinning recorded, we simulated the glacier using a one-dimensional owline model (see Methods). First, we simulated the modern ice surface at the sampling transects with a grounding line located upstream of the large overdeepening, and an advanced glacier that extended beyond the most distal GZW immediately offshore from Mackay Glacier and which honoured our empirically derived estimates of ice thickening along its longitudinal prole.

To assess whether the change in surface elevation at our thinning transects could be explained simply by recession of the grounding line, we forced the grounding line from the outer trough, through the zone of GZWs and across the reverse bed slope of the inner trough. This was achieved using a range of individual ocean temperature and calving rate scenarios as the forcing. We also investigated the effects of pulses of sea level rise on the glacier, for example, if the rate of eustatic sea level rise exceeded that of isostatic rebound. In all experiments, we were able to simulate retreat from the outermost GZW to upstream of the overdeepened Mackay Sea Valley in B300400 years, in good agreement with our chronology (Fig. 4; Supplementary Fig. 5a).

In all retreat scenarios, the grounding line uctuates in the outer trough, retreats more rapidly over the reverse bed slope and then re-establishes stability on the normal slope of the bedrock ridge upstream of the overdeepening (Fig. 4a). A higher frequency of grounding line positions were simulated in the outer trough, which in most cases match observed GZWs15. In the overdeepening, retreat into deeper water allows for thicker ice at the grounding line that is closer to oatation, leading to greater ice uxes and therefore further retreat (Fig. 4b). While grounding line retreat is generally faster for warmer ocean temperatures and higher relative sea level, the pattern of upstream thinning is similar in all forcing scenarios (Fig. 4a; Supplementary Fig. 5a); accelerated upstream thinning occurs as the grounding line retreats over the reverse bed slope. The increased rate of surface lowering is initially simulated at Low Ridge. Subsequently, a more pronounced thinning response occurs at Mt Suess/Gondola Ridge (upper). This spatial pattern of ice-dynamic thinning ts the

NATURE COMMUNICATIONS | 6:8910 | DOI: 10.1038/ncomms9910 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms9910

300

900 GRU

800 LR

Elevation (m above ice)

GRU

Surface elevation (m a.s.l.)

250

800

600

200

700

400

Ocean temperature

150

100

600

200

1.80 C 1.75 C 1.70 C 1.65 C

2,000

Chronology Modelled

50 5 10 15 20

500 160 180 200 220 240

0 160 180 200 220 240

Elevation (m a.s.l.)

Exposure age (ka)

Grounding line position (km)

1,000

GRU LR

Grounding line positions

(frequency)

1,000

lce flux (m3 per year) ( 104)

Model time (years)

100 200

100 120 140 160 180 200 220 240

Distance downstream (km from divide)

Figure 4 | Modelled retreat of Mackay Glacier through the overdeepened trough with corresponding accelerated ice ux and upstream thinning at our transects. (a) An evolving surface prole is shown at 10-year intervals for one model run ( 1.75 C ocean temperature forcing scenario). The modern

surface (dot-dashed line), simulated advanced prole (bold black line) and sample locations (white circles) are also shown. A histogram displays the frequency of grounding line positions for all modelled scenarios, with peaks broadly corresponding to most grounding-zone wedges (brown triangles). This modelling identies that stable glacier positions occur in the outer trough, unstable fast retreat occurs on the reverse bed slope and then restored stability occurs in the vicinity of the present-day grounding line. (b) An evolving ice ux corresponds to the surface prole. As the grounding line retreats over the reverse bed slope, thicker ice allows for greater ice uxes and therefore potential for mass loss. The largest ux of ice occurs in the deepest part of the trough, where the glacier is thickest and a steep ice surface slope facilitates higher velocities. Ice ux remains relatively large until the glacier stabilises on the normal-sloping bed, temporarily maintained by high surface velocities. The dashed black line denotes the longitudinally averaged (25 km) ice ux (5 104) over this retreat period. Simulated surface thinning at our Mt Suess/Gondola Ridge (upper; GRU) and Low Ridge (LR) transects is shown in

scatter plots (a) for a range of ocean forcing scenarios. Some initial thinning occurs in all scenarios and at both transects when the grounding line is located in the outer trough. Accelerated surface lowering occurs rst at LR simultaneous with retreat over the reverse bed slope, while the most rapid thinning is simulated upstream at GRU, matching our surface lowering chronologies.

thinning gradients of our respective transect chronologies (Fig. 3) and the idea that thinning propagated upstream following perturbations at the grounding line11,33,39. In particular, prolonged thinning that still occurs when the grounding line is located on the normal bed slope may reect the steep surface slope in this portion, which is able to facilitate high ice velocities and maintain ice loss40.

A rapid episode of dynamic thinning in our model may have occurred due to a process other than that provided by the reverse bed slope. Therefore, we also carried out an idealized experiment in which Mackay Sea Valley was replaced by a at bed (Supplementary Fig. 5b). Under this modied bed geometry, modelled retreat occurs over a longer period of B800 years without substantial jumps in the grounding line position. The e-folding response time41 of the glacier for this experiment is

58945 model years, opposed to the faster response of 290126 model years when ice retreated through the overdeepened trough. Importantly, there is no increase in the rate of surface lowering in this simulation, and it therefore fails to replicate the rapid thinning recorded in our Mt Suess/Gondola Ridge (upper) and Low Ridge transects.

In summary, simulated accelerated thinning corresponds with increased ice ux as the grounding line retreats over the reverse bed slope, irrespective of the forcing applied. Together, the modelling experiments and regional geological constraints suggest that the substantial thinning recorded at our nunatak sites most likely resulted from marine ice sheet instability2,4,5. Reduced buttressing from ice shelf disintegration5,42 is not required to explain rapid retreat and dynamic thinning but it could have helped facilitate this retreat.

6 NATURE COMMUNICATIONS | 6:8910 | DOI: 10.1038/ncomms9910 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms9910 ARTICLE

DiscussionRapid and protracted ice loss of marine-terminating glaciers that retreat into deeper water has been observed in temperate and sub-polar settings7,13,4345. Similar to Mackay Glacier, rapid retreat and ice-dynamic thinning becomes self-sustaining in these glaciers as the grounding line recedes across a reverse bed slope8,33,46. In these environments, a combination of atmospheric and oceanic warming, and rising sea level may produce a succession of negative mass balance years that is sufcient to initiate topographically driven retreat8,47. Today, accelerated ice loss occurs where warm subsurface waters inltrate glacial troughs, such as at Jakobshavn Isbr and Pine Island Bay48,49. The record at Mackay Glacier demonstrates that unstable tidewater glacier retreat via the marine instability mechanism is also possible in a more polar oceanic environment, 4600 km from the Antarctic continental shelf edge. Our modelling indicates that small changes in ocean temperature and/or sea level rise were sufcient to initiate prolonged retreat into deeper water with corresponding dynamic thinning and ice loss.

Most tidewater glacier cycles have been observed to last for up to several decades as a result of small variations in the bed topography7,13,44,45. Our chronology provides a time-averaged estimate of ice-discharge events and ice-dynamic thinning for the retreat phase of a prolonged tidewater glacier cycle. The scale of this outlet glacier retreat is similar to that observed in the twentieth and early twenty-rst centuries in well-known temperate or sub-polar tidewater environments, such as Glacier Bay, Alaska and Jakobshavn Isbr, Greenland, which have retreated B100 and 35 km, respectively, between B1850 and 2010 AD5053. However, at Mackay Glacier, we show that this retreat and associated accelerated ice surface lowering continued for at least two and a half centuries, providing a record of tidewater glacier behaviour that extends over a longer period than historic observations. The timescale of this retreat is independently supported by the modelled response time of the glacier (B300 years).

Glacial troughs similar to Mackay Sea Valley (excavated to B1 km below sea level or greater) and their associated reverse bed slopes are common features on the inner continental shelf surrounding Antarctica1,54. Ice-dynamic thinning and rapid mass loss may therefore have occurred in many parts of Antarctica, as grounding lines retreated through these troughs. Marine grounding lines retreated from the outer to inner continental shelf following the LGM, where they encountered troughs with reverse bed slopes54,55. We suggest that phases of accelerated ice sheet thinning during the Holocene in West Antarctica and on the Antarctic Peninsula17,19,5658 may also have resulted from the topographic instability that drove the retreat of Mackay Glacier.

Our new data from East Antarctica reveals that accelerated ice sheet thinning as a consequence of marine ice sheet instability can lead to hundreds of metres of ice surface lowering. We highlight that periods of rapid, topographically driven ice loss are likely typical of marine ice sheet margins, even in East Antarctica. These ndings provide condence for models that are able to simulate such ice-dynamic retreat14,59,60 and which predict centennial-scale responses in locations where overdeepened basins extend inland for many tens of kilometres from glacier or ice sheet margins14.

Methods

Sample collection. Forty-four samples were collected for surface-exposure dating from two nunataks in the lower reaches of Mackay Glacier. The upstream nunatak comprises Gondola Ridge, an elongate granitic ridge, and Mt Suess, a basaltic dome that peaks at 1,127 m a.s.l. Low Ridge is an area of exposed granitic bedrock (B205 m above the ice surface) that occurs downstream of the present-day grounding line. Geomorphological evidence shows that Mackay Glacier was larger in the geologically recent past61. Moraines at Cuff Cape, downstream of the

present-day grounding line, record small uctuations of Mackay Glacier and adjacent New Glacier termini during the last two centuries36. Before this, its glacial history is unconstrained, although regional-scale ice sheet reconstructions suggest that Mackay Glacier along with other outlets on the Victoria Land coast thickened at the LGM, due to buttressing by grounded ice in the Ross Sea62.

We carried out altitudinal sampling transects at Mt Suess/Gondola Ridge to constrain thickness at the LGM and subsequent initial thinning, and at Low Ridge to record more recent thinning near the present-day terminus (Supplementary Fig. 1). To minimize the risk of sampling clasts containing an inherited cosmogenic inventory, erratics were sampled that showed signs of glacial transport and erosion, evidenced by facetted and abraded surfaces. Erratics that rested in locations unaffected by post-depositional processes were also prioritized to provide a simple exposure signal. Such samples were generally either perched on eroded bedrock or propped up by a thin layer (13 cm) of draped glacial till. Samples consisted of cobbles (B1025 cm in length) and some small boulders (o1 m across) that were subsampled in the eld. Details of samples collected are listed in Supplementary Data 1.

Sample processing and measurement. Preparation of samples for surface-exposure dating was conducted at the sedimentology, cosmogenic nuclide and geochemistry laboratories of Victoria University of Wellington and GNS Science, New Zealand. The upper B37 cm of each sample was extracted for processing using a large circular saw. This material was then crushed with a ne jaw-crusher and sieved to retain sand-sized grains (250500 mm). Each sample was initially cleaned using a Frantz Isodynamic separator (at B0.5 A and a 10 tilt) to remove the magnetic component. Further cleaning was required to remove non-quartz minerals and etch any meteoric beryllium from the outer surface of the grains. This was achieved for each sample (up to 150 g) with a 1-day leach in hydrochloric acid and then three 2-day leaches in a 5% mixture of weak hydrouoric acid and nitric acid, warmed and rotated on hotdog rollers. A nal 1-h hydrouoric acid (7 M) leach and concentrated aqua regia cleaning was carried out to remove any remaining meteoric beryllium before sample dissolution and (in situ) beryllium extraction.

Beryllium was extracted following established geochemical procedures63.

A 9Be spike (B0.150.18 mg per sample) was added to the quartz samples, which were then dissolved in concentrated hydrouoric acid. Anion exchange columns were rst used to remove iron, and then cation exchange columns were used to isolate beryllium from iron, aluminium, titanium, sodium and magnesium. Procedural blanks yielded mean 10Be/9Be ratios of 1.38 10 15 with a s.d. of

5.91 10 16. Analytical measurements of 10Be/9Be were undertaken at ETH

Zurich mass spectrometry facilities using both Tandem and Tandy accelerator mass spectrometers64. Low measureable yields of a small number samples resulted in less precise exposure ages with B1733% uncertainties (samples GR34, GR37,

GR47, GR51, GR54, GR62b, GR64 and GR83). All samples were measured relative to the ETH Zurich in house standard S2007N (nominal 10Be/9Be ratio of28.10.8 10 12), which in turn was calibrated relative to the ICN 01-5-1

standard (10Be/9Be ratio of 27.090.3 10 12) (ref. 65), and were corrected with

procedural blanks.

Surface-exposure age calculation and outliers. Surface-exposure ages were calculated from the measured concentrations of 10Be, corrected for topographic shielding, sample thickness, quartz density (2.7 g cm3) and an Antarctic atmospheric pressure gradient, using CRONUS-Earth online calculator66 (Supplementary Data 1). The production rate of 10Be is currently unconstrained in Antarctica,and production rates from both a global data set66 (production due to spallation of4.490.39 to 4.960.43 atoms g per year) and New Zealand (NZ) calibration site67 (3.840.08 atoms per gram per year) have been applied to Antarctic chronologies in the past17,19. We prioritize the global production rate in this study as the high-precision NZ production rate (44 S) would produce ages with articially low uncertainties at Mackay Glacier (77 S), given the uncertain temporal and spatial scaling. The choice of production rate does not signicantly affect the reconstructed thinning rate at Mackay Glacier, and instead primarily inuences its absolute timing (Supplementary Fig. 2).

Four sample outliers were identied in the data set. On Gondola Ridge (mid-lower), two samples were erroneously old or young (GR38 and GR40), representing micro-inheritance and some post-depositional effects, respectively. At Gondola Ridge (lower), the oldest sample (13,0803,700 years BP, GR83) is not consistent with the thinning trends recorded in adjacent transects; if the surface elevation here was o50 m above present before B13 ka, then an unfeasible ice surface prole would exist with sites upstream and downstream, which record glacier surfaces 4200 m higher at this time. For this reason, sample GR83 most likely contains a small amount of inheritance. A single sample on Low Ridge (CC93) is not consistent with exposure ages above and below its elevation. This minor outlier has probably been affected by post-depositional rotation or spalling.

Bayesian age-elevation modelling. We applied Bayesian modelling68 to rene the timing of surface lowering recorded in our high-density chronologies, adapted here for altitudinal surface-exposure transects. This approach narrows the possible age uncertainty ranges based on the respective likely age distribution and elevation,

NATURE COMMUNICATIONS | 6:8910 | DOI: 10.1038/ncomms9910 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms9910

with the assumption that older and higher samples would be exposed by glacier thinning before the lower samples. At elevations where internal uncertainties overlapped, we calculated a weighted mean and standard error. We used external exposure age uncertainties that represent production rate uncertainty from spallation and muons, as we compare the output age ranges to other chronologies66. Although the Gondola Ridge (mid-lower) chronology is consistent with other transects, it was not possible to apply Bayesian age-elevation modelling at this site, as the ages show slightly larger stratigraphic scatter68. However, at Mt Suess/Gondola Ridge (upper), Gondola Ridge (lower) and Low Ridge, this statistical analysis of high-density sample transects allowed the uncertainty of raw exposure ages to be reduced (Fig. 3ac; Supplementary Data 1).

Regression analysis of rapid thinning episode. To estimate the rate and duration of rapid thinning at Mackay Glacier, linear regression analysis was carried out on data from the Mt Suess/Gondola Ridge (upper) and Low Ridge transects for the period of B6.86.0 ka. Error-weighted least-squares regression was applied randomly to normally distributed exposure ages (2 s) through a 4,000-iteration Monte Carlo simulation. Rates and durations of rapid thinning were estimated from the distribution of feasible, positive-sloping linear regressions, with uncertainty generally reective of the number of samples contributing to each transect and their respective uncertainties. Linear rates of thinning were calculated using the raw exposure ages at Low Ridge and the Bayesian-modelled ages at Mt Suess/ Gondola Ridge (upper). The duration of this thinning episode was additionally determined, based on the start and end points of the modelled regressions.

Glacier owline model. We employed a one-dimensional nite difference owline model to investigate whether retreat of Mackay Glacier as a consequence of grounding line retreat through an overdeepening could result in the magnitude of thinning observed at our exposure dating transects. The owline model has previously been applied to Transantarctic Mountain outlet glaciers and is fully described in Golledge and Levy69. Ice thickness change is determined using the mass conservation equation, as the balance between ice ux and net accumulation, which is then integrated through time. Zero ice ux is prescribed at the top of the domain. The shallow-ice approximation is used in combination with a longitudinal averaging scheme to calculate basal shear stress at 1-km horizontal resolution. As the past ice shelf extent is not well constrained in time or space, we only implement the numerical scheme upstream of the grounding line and do not simulate an ice shelf. Basal sliding occurs where basal temperatures are close to melting point, while ow through creep is controlled by a temperature-dependent deformation rate with an enhancement coefcient. Here we used a spatially variable basal sliding coefcient to account for differences in basal traction between the higher elevation Transantarctic Mountains portion, where Mackay Glacier is currently grounded on bedrock, and the downstream marine portion, where soft sediments occur. Enhanced sliding in the marine portion is suggested by the occurrence of elongate glacial lineations in the Mackay Glacier trough (Supplementary Fig. 4), indicative of former fast ice ow from subglacial deformation over softer sediment.

Net mass balance is determined using a positive degree-day scheme, with mean surface air temperature, annual temperature range and precipitation rate initially prescribed. Surface melt is calculated if mean surface temperatures are above the freezing point. Although ice shelves are not directly incorporated into the model, negative mass balance is imposed in areas where the glacier bed is below sea level. Sub-shelf melting is calculated using an ocean temperature-based melt scheme when the ice thickness at the grounding line is less than oatation thickness69, controlled by a given seasonal ocean temperature cycle. Otherwise, mass loss at the grounding line occurs from tidewater calving using a scalable calving rate coefcient as a factor of a calving constant69.

The model was initialized with bed topography and ice thickness obtained from multiple radar data proles. Glacier surface elevations were obtained from airborne altimetry measurements Investigating the Cryospheric Evolution of the Central Antarctic Plate (ICECAP), while the bed topography was interpolated from BEDMAP2 data1 that was corrected along the owline with radio-echoed ice thickness measurements70. The model domain extends over 244 km from the ice divide near Taylor Dome to the most distal offshore GZW (Fig. 1c; Supplementary Fig. 4). The owline used for Mackay Glacier followed the sinuous yet laterally constrained ow path through the Transantarctic Mountains, and was determined from present-day surface velocities, surface ow stripes, offshore lineations and trough bathymetry. To evaluate the thinning response resulting from migration of the grounding line without effects from lateral variations39, we applied a uniform width of 8 km along the length of the model domain, which approximates the glacier width in the lower reaches of the modern Mackay Glacier and the offshore trough.

Model experiments used climate and glaciological parameterisations from eld measurements and commonly used physical values (Supplementary Table 2). First, modern and advanced steady-state ice surface proles were simulated by adjusting the precipitation at sea level, precipitation lapse rate, calving coefcient, ocean temperature, and basal sliding and creep enhancement factors. A more extensive advanced glacier was simulated with the grounding line beyond the end of the domain and the surface elevation o100 m above the Mt Suess/Gondola Ridge transect, where an LGM ice surface is recorded prior to rapid thinning. This advanced state was simulated by decreasing both the calving coefcient and ocean

temperature while maintaining the other tuned parameters of the present-day simulation.

Transient simulations were then carried out by forcing retreat of the grounding line from the outer trough. Glacier retreat was achieved for a range of scenarios by enhancing the calving coefcient and summer ocean temperature. No retreat scenarios simulated ice thicker than oatation with a grounding line positioned within the trough in the model domain, and therefore all negative changes in mass balance were achieved through basal melt at the grounding line. Further model experiments used increased sea level to force retreat of the glacier; however, a similar pattern of retreat and thinning was observed irrespective of environmental forcing. In all scenarios, the grounding line initially uctuates in the outer trough, but then accelerates over the reverse bed slope, corresponding to increased ice ux and rapid glacier thinning. In some cases (for example, at B210 km from the divide), temporary grounding line stability that is implied by the presence of GZWs is not simulated in any retreat scenarios. This could be because of reduced trough width and complex lateral drag in this location (Supplementary Fig. 4), a stability feedback provided by the positive relief of GZW deposition, and/or enhanced stability from the presence of an ice shelf at that time, none of which were incorporated within the model.

References

1. Fretwell, P. et al. Bedmap2: improved ice bed, surface and thickness datasets for Antarctica. Cryosphere 7, 375393 (2013).

2. Weertman, J. Stability of the junction of an ice sheet and an ice shelf. J. Glaciol. 13, 311 (1974).

3. Mercer, J. H. West Antarctic ice sheet and CO2 greenhouse effect: a threat of disaster. Nature 271, 321325 (1978).

4. Thomas, R. H. & Bentley, C. R. A model for Holocene retreat of the West Antarctic ice sheet. Quat. Res. 10, 150170 (1978).

5. Schoof, C. Ice sheet grounding line dynamics: Steady states, stability, and hysteresis. J. Geophys. Res. Earth Surf. 112, F03S28 (2007).

6. Mengel, M. L. A. Ice plug prevents irreversible discharge from East Antarctica. Nat. Clim. Change 4, 17586798 (2014).

7. Howat, I. M., Joughin, I. & Scambos, T. A. Rapid changes in ice discharge from Greenland outlet glaciers. Science 315, 15591561 (2007).

8. Vieli, A., Funk, M. & Blatter, H. Flow dynamics of tidewater glaciers: a numerical modelling approach. J. Glaciol. 47, 595606 (2001).

9. Pritchard, H. D., Arthern, R. J., Vaughan, D. G. & Edwards, L. A. Extensive dynamic thinning on the margins of the Greenland and Antarctic ice sheets. Nature 461, 971975 (2009).

10. Miles, B., Stokes, C., Vieli, A. & Cox, N. Rapid, climate-driven changes in outlet glaciers on the Pacic coast of East Antarctica. Nature 500, 563566 (2013).

11. Payne, A. J., Vieli, A., Shepherd, A. P., Wingham, D. J. & Rignot, E. Recent dramatic thinning of largest West Antarctic ice stream triggered by oceans. Geophys. Res. Lett. 31, L23401 (2004).

12. Dupont, T. & Alley, R. Assessment of the importance of ice-shelf buttressing to ice-sheet ow. Geophys. Res. Lett. 32, L04503 (2005).

13. Arendt, A. A., Echelmeyer, K. A., Harrison, W. D., Lingle, C. S. & Valentine, V. B. Rapid wastage of Alaska glaciers and their contribution to rising sea level. Science 297, 382386 (2002).

14. Nick, F. M. et al. Future sea-level rise from Greenlands main outlet glaciers in a warming climate. Nature 497, 235238 (2013).

15. Greenwood, S. L., Gyllencreutz, R., Jakobsson, M. & Anderson, J. B. Ice-ow switching and East/West Antarctic Ice Sheet roles in glaciation of the western Ross Sea. Geol. Soc. Am. Bull 124, 17361749 (2012).

16. Anderson, J. B. et al. Ross Sea paleo-ice sheet drainage and deglacial history during and since the LGM. Quat. Sci. Rev. 100, 3154 (2014).

17. Stone, J. O. et al. Holocene deglaciation of Marie Byrd land, west Antarctica. Science 299, 99102 (2003).

18. Rignot, E., Mouginot, J. & Scheuchl, B. Ice ow of the Antarctic ice sheet. Science 333, 14271430 (2011).

19. Johnson, J. S. et al. Rapid thinning of Pine Island Glacier in the early Holocene. Science 343, 9991001 (2014).

20. Masson-Delmotte, V. et al. A comparison of the present and last interglacial periods in six Antarctic ice cores. Clim. Past 7, 397423 (2011).

21. Jouzel, J. et al. Orbital and millennial Antarctic climate variability over the past 800,000 years. Science 317, 793796 (2007).

22. Stenni, B. et al. Expression of the bipolar see-saw in Antarctic climate records during the last deglaciation. Nat. Geosci. 4, 4649 (2011).

23. Lambeck, K., Rouby, H., Purcell, A., Sun, Y. & Sambridge, M. Sea level and global ice volumes from the Last Glacial Maximum to the Holocene. Proc. Natl Acad. Sci. USA 111, 1529615303 (2014).

24. Argus, D. F., Peltier, W., 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. Intl 198, 537563 (2014).

25. Hall, B. L., Baroni, C. & Denton, G. H. Holocene relative sea-level history of the Southern Victoria Land Coast, Antarctica. Glob. Planet. Change 42, 241263 (2004).

8 NATURE COMMUNICATIONS | 6:8910 | DOI: 10.1038/ncomms9910 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms9910 ARTICLE

26. Pahnke, K. & Sachs, J. Sea surface temperatures of southern midlatitudes 0160 kyr B.P. Paleoceanography 21, PA2003 (2006).

27. McKay, R. et al. Retreat history of the Ross Ice Sheet (Shelf) since the Last Glacial Maximum from deep-basin sediment cores around Ross Island. Palaeogeogr. Palaeoclimatol. Palaeoecol. 260, 245261 (2008).

28. Hall, B. L. & Denton, G. H. Extent and chronology of the Ross Sea ice sheet and the Wilson Piedmont Glacier along the Scott Coast at and since the last glacial maximum. Geogr. Ann. A 82, 337363 (2000).

29. Licht, K. J., Jennings, A. E., Andrews, J. T. & Williams, K. M. Chronologyof late Wisconsin ice retreat from the western Ross Sea, Antarctica. Geology 24, 223226 (1996).

30. Domack, E. W., Jacobson, E. A., Shipp, S. & Anderson, J. B. Late Pleistocene Holocene retreat of the West Antarctic Ice-Sheet system in the Ross Sea:part 2sedimentologic and stratigraphic signature. Geol. Soc. Am. Bull 111, 15171536 (1999).

31. Hall, B. L. & Denton, G. H. New relative sea-level curves for the southern Scott Coast, Antarctica: evidence for Holocene deglaciation of the western Ross Sea.J. Quat. Sci. 14, 641650 (1999).32. Stearns, L. A., Smith, B. E. & Hamilton, G. S. Increased ow speed on a large East Antarctic outlet glacier caused by subglacial oods. Nat. Geosci. 1, 827831 (2008).

33. Nick, F. M., Vieli, A., Howat, I. M. & Joughin, I. Large-scale changes in Greenland outlet glacier dynamics triggered at the terminus. Nat. Geosci. 2, 110114 (2009).

34. Vieli, A. & Nick, F. M. Understanding and modelling rapid dynamic changes of tidewater outlet glaciers: issues and implications. Surv. Geophys. 32, 437458 (2011).

35. Pattyn, F. Antarctic subglacial conditions inferred from a hybrid ice sheet/ice stream model. Earth. Planet. Sci. Lett. 295, 451461 (2010).

36. Mller, P. Subrecent moraine ridge formation on Cuff Cape, Victoria Land, Antarctica. Geogr. Ann. A 77, 8394 (1995).

37. Scambos, T. A., Bohlander, J., Shuman, C. A. & Skvarca, P. Glacier acceleration and thinning after ice shelf collapse in the Larsen B embayment, Antarctica. Geophys. Res. Lett. 31, L18402 (2004).

38. Dupont, T. & Alley, R. Role of small ice shelves in sea-level rise. Geophys. Res. Lett. 33, L09503 (2006).

39. Jamieson, S. S. et al. Ice-stream stability on a reverse bed slope. Nat. Geosci. 5, 799802 (2012).

40. Joughin, I., Abdalati, W. & Fahnestock, M. Large uctuations in speed on Greenlands Jakobshavn Isbrae glacier. Nature 432, 608610 (2004).

41. Oerlemans, J. Climate sensitivity of Franz Josef Glacier, New Zealand, as revealed by numerical modeling. Arct. Alp. Res. 29, 233239 (1997).

42. Pattyn, F., Huyghe, A., De Brabander, S. & De Smedt, B. Role of transition zones in marine ice sheet dynamics. J. Geophys. Res. Earth Surf. 111, F02004 (2006).

43. Meier, M. F. & Post, A. Fast tidewater glaciers. J. Geophys. Res. 92, 9051 (1987).44. Mouginot, J. & Rignot, E. Ice motion of the Patagonian Iceelds of South America: 1984-2014. Geophys. Res. Lett. 42, 14411449 (2015).

45. Mercer, J. The response of fjord glaciers to changes in the rn limit. J. Glaciol. 3, 850858 (1961).

46. ONeel, S., Pfeffer, W. T., Krimmel, R. & Meier, M. Evolving force balance at Columbia Glacier, Alaska, during its rapid retreat. J. Geophys. Res. Earth Surf. 110, F03012 (2005).

47. Joughin, I. et al. Ice-front variation and tidewater behavior on Helheim and Kangerdlugssuaq Glaciers, Greenland. J. Geophys. Res. 113, F01004 (2008).48. Holland, D. M., Thomas, R. H., De Young, B., Ribergaard, M. H. & Lyberth, B. Acceleration of Jakobshavn Isbrae triggered by warm subsurface ocean waters. Nat. Geosci. 1, 659664 (2008).

49. Jacobs, S. S., Jenkins, A., Giulivi, C. F. & Dutrieux, P. Stronger ocean circulation and increased melting under Pine Island Glacier ice shelf. Nat. Geosci. 4, 519523 (2011).

50. Connor, C., Streveler, G., Post, A., Monteith, D. & Howell, W. The Neoglacial landscape and human history of Glacier Bay, Glacier Bay National Park and Preserve, southeast Alaska, USA. Holocene 19, 381393 (2009).

51. Larsen, C. F., Motyka, R. J., Arendt, A. A., Echelmeyer, K. A. & Geissler, P. E. Glacier changes in southeast Alaska and northwest British Columbia and contribution to sea level rise. J. Geophys. Res. 112, F01007 (2007).

52. Csatho, B., Schenk, T., Van Der Veen, C. J. & Krabill, W. B. Intermittent thinning of Jakobshavn Isbrae, West Greenland, since the little ice age. J. Glaciol. 54, 131144 (2008).

53. Joughin, I. et al. Continued evolution of Jakobshavn Isbrae following its rapid speedup. J. Geophys. Res. 113, F04006 (2008).

54. Livingstone, S. J. et al. Antarctic palaeo-ice streams. Earth Sci. Rev. 111, 90128 (2012).

55. Bentley, M. J. et al. A community-based geological reconstruction of Antarctic Ice Sheet deglaciation since the Last Glacial Maximum. Quat. Sci. Rev. 100, 19 (2014).

56. 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, 411414 (2010).

57. Balco, G. & Schaefer, J. M. Exposure-age record of Holocene ice sheet and ice shelf change in the northeast Antarctic Peninsula. Quat. Sci. Rev. 59, 101111 (2013).

58. Fogwill, C. J. et al. Drivers of abrupt Holocene shifts in West Antarctic ice stream direction determined from combined ice sheet modelling and geologic signatures. Antarctic Sci. 26, 674686 (2014).

59. Favier, L. et al. Retreat of Pine Island Glacier controlled by marine ice-sheet instability. Nat. Clim. Change 4, 117121 (2014).

60. Joughin, I., Smith, B. E. & Medley, B. Marine Ice Sheet Collapse Potentially Under Way for the Thwaites Glacier Basin, West Antarctica. Science 344, 735738 (2014).

61. 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, 840857 (2004).

62. Denton, G. H. & Hughes, T. J. Reconstruction of the Ross ice drainage system, Antarctica, at the last glacial maximum. Geogr. Ann. A 82, 143166 (2000).

63. Von Blanckenburg, F., Belshaw, N. & ONions, R. Separation of 9Be and cosmogenic 10Be from environmental materials and SIMS isotope dilution analysis. Chem. Geol. 129, 9399 (1996).

64. Christl, M. et al. The ETH Zurich AMS facilities: performance parameters and reference materials. Nucl. Instr. Methods Phys. Res. B 294, 2938 (2013).

65. Nishiizumi, K. et al. Absolute calibration of 10Be AMS standards. Nucl. Instr. Methods Phys. Res. B 258, 403413 (2007).

66. Balco, G., Stone, J. O., Lifton, N. A. & Dunai, T. J. A complete and easily accessible means of calculating surface exposure ages or erosion rates from

10Be and 26Al measurements. Quat. Geochronol. 3, 174195 (2008).67. Putnam, A. et al. In situ cosmogenic 10Be production-rate calibration from the Southern Alps, New Zealand. Quat. Geochronol. 5, 392409 (2010).

68. Ramsey, C. B. Deposition models for chronological records. Quat. Sci. Rev. 27, 4260 (2008).

69. Golledge, N. & Levy, R. Geometry and dynamics of an East Antarctic Ice Sheet outlet glacier, under past and present climates. J. Geophys. Res. Earth Surf. 116, F03025 (2011).

70. Calkin, P. E. Subglacial geomorphology surrounding the ice-free valleys of southern Victoria Land, Antarctica. J. Glaciol. 13, 415429 (1974).

Acknowledgements

We acknowledge eld support from Antarctica New Zealand. R.S.J. was supported by a VUW Victoria Doctoral Scholarship. A.N.M. and N.R.G. were supported by the New Zealand Government (MBIE) Past Antarctic Climate and Future Implications program. Sample preparation and measurement were supported by the Antarctic Research Centres Endowed Development Fund and School of Geography, Environment and Earth Sciences (VUW). C.J.F. is supported by the Australian Research Council. S.L.G. is supported by a Swedish Research Council grant. We thank E.G.C. Smith for assistance with regression analysis, D.A. Young for airborne radar data (ICECAP), R.M. McKay for critical discussions and three anonymous reviewers who helped to improve the manuscript.

Author contributions

R.S.J., A.N.M., K.P.N., N.R.G. and C.J.F. designed the research, and R.S.J. carried out the

10Be dating, data analysis and glacier modelling. R.S.J. was assisted by C.J.F. and K.P.N. (eldwork), K.P.N. (10Be chemistry) and N.R.G. (glacier modelling). P.W.K. and M.C. carried out AMS measurements, and S.L.G. helped interpret the marine geology. R.S.J. and A.N.M. drafted the manuscript and all authors contributed to the nal paper.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/ http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing nancial interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Jones, R. S. et al. Rapid Holocene thinning of an East Antarctic outlet glacier driven by marine ice sheet instability. Nat. Commun. 6:8910 doi: 10.1038/ ncomms9910 (2015).

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecommons.org/licenses/by/4.0/

NATURE COMMUNICATIONS | 6:8910 | DOI: 10.1038/ncomms9910 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 9

Word count: 9417

Show less

Abstract

Translate

Outlet glaciers grounded on a bed that deepens inland and extends below sea level are potentially vulnerable to 'marine ice sheet instability'. This instability, which may lead to runaway ice loss, has been simulated in models, but its consequences have not been directly observed in geological records. Here we provide new surface-exposure ages from an outlet of the East Antarctic Ice Sheet that reveal rapid glacier thinning occurred approximately 7,000 years ago, in the absence of large environmental changes. Glacier thinning persisted for more than two and a half centuries, resulting in hundreds of metres of ice loss. Numerical simulations indicate that ice surface drawdown accelerated when the otherwise steadily retreating glacier encountered a bedrock trough. Together, the geological reconstruction and numerical simulations suggest that centennial-scale glacier thinning arose from unstable grounding line retreat. Capturing these instability processes in ice sheet models is important for predicting Antarctica's future contribution to sea level change.

Details

Title

Rapid Holocene thinning of an East Antarctic outlet glacier driven by marine ice sheet instability

Author

Jones, R S; Mackintosh, A N; Norton, K P; Golledge, N R; Fogwill, C J; Kubik, P W; Christl, M; Greenwood, S L

Pages

8910

Publication year

2015

Publication date

Nov 2015

Publisher

Nature Publishing Group

e-ISSN

20411723

Source type

Scholarly Journal

Language of publication

English

DOI

https://doi.org/10.1038/ncomms9910

ProQuest document ID

1736408174

Rapid Holocene thinning of an East Antarctic outlet glacier driven by marine ice sheet instability

Jump to:

Full text

Abstract

Details

Suggested sources