Scholarly article on topic 'Exposure age and ice-sheet model constraints on Pliocene East Antarctic ice sheet dynamics'

Exposure age and ice-sheet model constraints on Pliocene East Antarctic ice sheet dynamics Academic research paper on "Earth and related environmental sciences"

Share paper
Academic journal
Nature Communications

Academic research paper on topic "Exposure age and ice-sheet model constraints on Pliocene East Antarctic ice sheet dynamics"



Received 10 Oct 2014 | Accepted 24 Mar 2015 | Published 24 Apr 2015^^BoOl:ia1o38/ncomms8016| OPEN

Exposure age and ice-sheet model constraints on Pliocene East Antarctic ice sheet dynamics

Masako Yamane1,f, Yusuke Yokoyama1,2,3, Ayako Abe-Ouchi1,2,4, Stephen Obrochta1,5, Fuyuki Saito4, Kiichi Moriwaki6 & Hiroyuki Matsuzaki7

The Late Pliocene epoch is a potential analogue for future climate in a warming world. Here we reconstruct Plio-Pleistocene East Antarctic Ice Sheet (EAIS) variability using cosmogenic nuclide exposure ages and model simulations to better understand ice sheet behaviour under such warm conditions. New and previously published exposure ages indicate interior-thickening during the Pliocene. An ice sheet model with mid-Pliocene boundary conditions also results in interior thickening and suggests that both the Wilkes Subglacial and Aurora Basins largely melted, offsetting increased ice volume. Considering contributions from West Antarctica and Greenland, this is consistent with the most recent IPCC AR5 estimate, which indicates that the Pliocene sea level likely did not exceed + 20 m on Milankovitch timescales. The inception of colder climate since ~3Myr has increased the sea ice cover and inhibited active moisture transport to Antarctica, resulting in reduced ice sheet thickness, at least in coastal areas.

1 Atmosphere and Ocean Research Institute, University of Tokyo, 5-1-5 Kashiwanoha, Chiba 277-8564, Japan. 2 Department of Earth and Planetary Sciences, University of Tokyo, 7-3-1 Hongo, Bunkyoku, Tokyo 113-0033, Japan. 3 Japan Agency for Marine-Earth Science and Technology, 2-15 Natsushima-cho, Yokosuka 237-0061, Japan. 4Japan Agency for Marine-Earth Science and Technology, 3173-25, Showa-machi, Kanazawa-ku, Yokohama, Kanagawa 236-0001, Japan. 5 Faculty of International Resource Science, Akita University 1-1 Tegata Gakuin-cho, Akita 010-8502, Japan. 6 National Institute of Polar Research, 10-3, Midoricho, Tachikawa, Tokyo 190-8518, Japan. 7The University Museum, University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan. f Present address: Japan Agency for Marine-Earth Science and Technology, 2-15 Natsushima-cho, Yokosuka 237-0061, Japan. Correspondence and requests for materials should be addressed to Y.Y. (email:

The climatic optimum of the Late Pliocene epoch (3.3-3.0 Myr) is often considered an analogue for near future conditions1'2 because of globally warmer temperatures and similar CO2 concentrations relative to today3-8. This resulted in lower continental ice volume, in addition to which there is 'high confidence' of higher global sea level during the Pliocene2. Published reconstructions of the maximum Pliocene sea level vary between + 5 and + 50 m compared with present9-13, with the current best estimate being the IPCC AR5 value of no more than + 20 m (ref. 2). The Greenland Ice Sheet and the West Antarctic Ice Sheet (WAIS), both of which are thought to have largely melted during Pliocene interglacials14-16, currently store b 13 m sea-level equivalent (SLE) ice. Modelling studies17,18 and recent observational evidence19,20 indicate that the margins of East Antarctica, especially the Wilkes Subglacial and Aurora Basins, further contributed to high Pliocene sea levels.

However, there is relatively little evidence directly constraining the response of the East Antarctic Ice Sheet (EAIS) to Pliocene climate evolution. Thick Cenozoic glaciomarine sequences are exposed at Mac Robertson Land (MRL) and the Transantarctic Mountains (Pagodroma Group and Sirius Group, respectively), suggesting that the ice margin was located several hundred kilometres inland of its present location during the Miocene or Pliocene21,22. An b 80-m diatomite section recovered in the ANDRILL core indicates higher Pliocene phytoplankton productivity in a largely ice-free Ross Sea15. During the Late Pliocene, the occurrence of diamictons in the Prydz Bay (ODP Site 1167 (ref. 23)) suggests glacial expansion.

To better constrain EAIS dynamics during the Pliocene, we present newly acquired and recently compiled in situ cosmogenic nuclide exposures ages together with model simulations of Antarctic ice thickness and distribution. We first reconstruct Pliocene EAIS configuration from the cumulative exposure time of bedrock, as reflected by the concentration of in situ produced cosmogenic nuclides (for example, 3He, 10Be, 14C, 21Ne, 26Al and 36Cl) in quartz from ice-free regions. The altitude of each sampling sites in turn provides a minimum estimate of ice sheet thickness at that location. We then compare this result with our model-based reconstruction to make inference on the ice sheet history of currently glaciated areas, such as the Wilkes Subglacial and Aurora Basins, from which geologic data are scarce. We finally use the ratio of two cosmogenic nuclides (for example, 26Al/10Be) to determine exposure history (simple or complex).


New and compiled exposure ages. New 10Be and 26Al data presented here were obtained from material collected at the central and eastern portions of the Sor Rondane Mountains, Droning Maud Land (DML) by the National Institute for Polar Research during JARE 27, 28, 31 and 32 (Fig. 1; Supplementary Table 1). The Sor Rondane Mountains form the eastern portion of a discontinuous, ~ 2,000 km east-west trending range that is located 200-400 km inland24 (Supplementary Fig. 1). There are a series of glacial tills located at varying altitudes in the Sor Rondane Mountains. These tills indicate five distinct exposure stages, the oldest (4 Myr) and most weathered of which is located at b 400 m above the present ice surface24. Exposure age and weathering extent of the tills progressively decreases with altitude to a minimum of ~ 40 ka at 3 m above the present ice sheet. This indicates that the ice sheet at this location was 400 m thicker than present before 4 Myr and gradually thinned to its current


Previously published exposure ages are compiled from five icefree areas of East Antarctica and grouped on the basis of distance from the ice sheet edge (Fig. 1). Area 1 (Supplementary Note 1) is

T-r 4 3 2 1 0

[ -2 .9

Figure 1 | Antarctic bedrock elevation. Antarctic bedrock elevation with grounded ice sheet boundary and attached, floating ice shelves. Thick black line shows grounded ice sheet boundary and thick grey lines show floating ice shelves. Thin black lines show ice sheet elevation (1,000 m contour interval70). Location of exposure ages reported here is shown with a red star and all other Area 1 sites with blue circles. Area 2 and 3 sites are shown with green squares and orange triangles, respectively. CL is Coats Land, DML is Droning Maud Land and MRL is Mac. Robertson Land, PEL is Princess Elizabeth Land and VL is Victoria Land. AB and WSB refer to the Aurora Basin and Wilkes Subglacial Basin, respectively.

b 300 km inland and includes our new data with previously published ages25,26 from the Sor Rondane Mountains in DML (Supplementary Fig. 1), as well as previously reported data from the Grove Mountains27-30 in Princess Elizabeth Land (Supplementary Fig. 2). Area 2 (Supplementary Note 2) includes exposure ages from samples and outcrops that are located b 100 km inland and encompasses the Richer Hills31 (Supplementary Fig. 3), the Allan Hills25,32 (Supplementary Fig. 4), the McMurdo Sound-Dry Valleys33-35 (Supplementary Fig. 5) in Victoria Land (VL) and the Shackleton Range36 in Coats Land (CL; Supplementary Fig. 6). Area 3 (Supplementary Note 3) includes exposure ages from samples and outcrops that are proximal to the ice sheet edge and includes Prince Charles Mountains37 in MRL (Supplementary Fig. 7) and the Deep Freeze Range38,39 (Supplementary Fig. 8) in VL.

Area 1. We report ten 10Be and six 26Al ages from the Sor Rondane Mountains. Two samples exhibit a complex exposure history (Supplementary Fig. 9; Supplementary Table 2). We also review previously published exposure ages25, from Sor Rondane Mountains, several of which were calculated with standard atmospheric thickness (see Methods). Where appropriate, we revise the ages using the scaling factor of Stone (ref. 40). Our new data together with previously reported ages from the Sor Rondane Mountains indicate that the EAIS in this area was thicker by at least 400 m than present between b2.5 and B3Myr (Fig. 2a); however, the presence of an age (near the limit of Be dating) obtained from a higher altitude raises the possibility that the ice sheet was b 600 m above present on at least one occasion before 4 Myr. Exposure histories are predominately complex during the Late Pleistocene, with simple exposure histories predominating for samples with larger exposure ages. The subsequent thinning pattern varies between the upstream and downstream portions of the ice sheet. Upstream ice thickness began decreasing during the Pliocene, declining to b 200 m above present by b 2.5 ± 0.5 Myr. Downstream ice thickness declined more rapidly, beginning at b 1.0 ±0.1 Myr, and reached the current level at b 0.4 Myr.

These results are consistent with géomorphologie evidence from glacial tills indicating that the ice sheet was 400 m thicker than present before 4Myr (ref. 24).

800 600 400 200 0

12 3 4

10Be min. exposure age (Myr)

2 200 -

2 3 4 5

10Be min. exposure age (Myr)

Half-life: 1.51 Myr Half-life: 1.36 Myr

Figure 2 | Area 1 ice sheet history. (a) Exposure ages from the S0r Rondane Mountains from this study (denoted in Fig. 1) and previous studies25,26. Circles and triangles indicate samples north (downstream) and south (upstream) of 72 °S, respectively. The black symbol indicates samples lacking 26Al measurement, and filled and open symbols refer to bedrock and erratic samples, respectively. (b) Compiled Grove Mountains exposure ages27-30. In both panels, 1s uncertainty is shown, and orange and green symbols indicate simple and complex exposure history, respectively.

For consistency, we recalculate the ages previously obtained from the Grove Mountains28'29 using a Be half-life of 1.51 Myr (Supplementary Fig. 10; Supplementary Table 3). Compiled ages indicate that the EAIS in this area was thicker by at least 150-220 m than present, at least once before 3.5 ±0.7 Myr, and then the ice sheet thinned by ~ 50-100 m at ~ 1.9 ±0.9 Myr. Complex exposure histories are found for samples with 10Be ages of less than 1.6 ± 0.2 Myr; simple exposure histories predominate for samples with larger exposure ages (Fig. 2b). These results are consistent with the distribution of non-glacial and glacial erosion features in the Grove Mountains. High-altitude locations are characterized by non-glacial weathering features, with glacial modification limited to ~100m above the present ice


Area 2. Strasky et a/.31 reported 10Be and 21Ne exposure ages from the Ricker Hills, VL. The four youngest samples exhibit complex histories, while the older three show simple histories. These data indicate that ice thickness was at least 400 m higher than present at ~0.6 ± 0.1 Myr (Fig. 3a). This is consistent with glacial erosional and depositional features up to the trimline and with wind erosion at higher altitudes31.

Glacial deposits, such as the Sirius Formation, are observed in the nunataks of the Allan Hills32. Tschudi et a/.32 reported exposure ages from Sirius Formation boulders and the underlying bedrock. Nishiizumi et a/.25 measured cosmogenic 10Be and Al in rock samples from Allan Hills, which we rescale for differing atmospheric thickness40. All samples exhibit a simple exposure history (Supplementary Fig. 11; Supplementary Table 3). Ice thickness in the Allan Hills was at least 70 m higher than the present at ~ 1.55 ± 0.15 Myr (Fig. 3b).

Several studies33-35 have reported exposure ages from glacial deposits exhibiting differing degrees of weathering (that is, Taylor Il-IVb) in the McMurdo Sound-Dry Valleys area of VL, which we compile and rescale40 here. Both simple and complex exposure histories are indicated (Supplementary Fig. 12; Supplementary Table 3). We interpret that the ice sheet in this

« 600

® 600

£ 400

12 3 4

10Be min. exposure age (Myr)

12 3 4

10Be min. exposure age (Myr)

-£ 100

-§ 40

J= 800

8 600 Π0

12 3 4

10Be min. exposure age (Myr)

12 3 4

10Be min. exposure age (Myr)

Figure 3 | Area 2 ice sheet history. Compiled exposure ages from (a) Ricker Hills31, (b) Allan Hills25,32, (c) Dry Valleys33-35 and (d) Shackleton Range36. In each panel, 1s uncertainty is shown; orange and green symbols indicate simple and complex exposure histories, respectively; black symbols indicate samples lacking 26Al measurement; and filled and open symbols refer to bedrock and erratic samples, respectively.

12 3 4

10Be min. exposure age (Myr)

CD 600 -

> 400 -

tio 200 a

(m) 800

res 600 pr

| 400 h

■§ 200 k

10Be min. exposure age (Myr)

Figure 4 | Area 3 ice sheet history. Compiled exposure ages from (a) the Fisher Massif (triangles) and Amery Oasis (circles) in the Prince Charles Mountains37, and (b) the Black Ridge (circles), Mt. Keinath (triangles), Mt. Abbott (squares) and Mt. Browning (inverted triangles) in the Deep Freeze Range38,39. Systematic uncertainties of ±1s are also indicated. Orange and green symbols indicate simple and complex exposure history, respectively, and black symbols indicate samples with only 10Be measurement. Filled and open symbols refer to bedrock and erratic samples, respectively.

area was thicker by at least 700 m than present, at least once before b 2.8 ±0.4 Myr (Fig. 3c). Thickness appears to decrease slightly within the Late Pliocene. From b2.2 ± 1.0 Myr, samples exhibit a baseline thickness of b 300 m. The most recent samples indicate a thickness similar to that of the present day. Complex exposure histories are found for samples with 10Be ages of less than 0.9± 0.2 Myr, with simple exposure histories predominating for samples with larger exposure ages.

Fogwill et al.36 reported three 10Be and two 26Al exposure ages from the horizontal bedrock of a plateau in the Shackleton Range, CL. Both samples with reported 26Al concentrations exhibit 10Be ages of b 1 and 2 Myr with a simple exposure history (Supplementary Fig. 13; Supplementary Table 3). These were originally considered to indicate that the ice sheet has been in a relatively stable condition for at least the past 1 Myr. Considering the additional sample that lacks a corresponding 26Al measurement and has a much older 10Be exposure age (~3Myr), we interpret these data to indicate that the EAIS at this location was thicker by at least 750 m than present, at least once before 3Myr, before thinning by B350m by 1 Myr (Fig. 3d). This is consistent with a meteorite exposure age that indicates that thickness did not exceed 350 m during the past 415 ka and with the inferred late Miocene or Pliocene age of striated till and bedrock at a higher altitude41.

Area 3. Fink et al.37 presented a deglacial chronology based on seven cosmogenic 10Be and 26Al exposure ages from three areas of the northern Prince Charles Mountains, MRL. All data exhibit a simple exposure history (Supplementary Fig. 14; Supplementary Table 3). The elevation of the EAIS in both the Fisher Massif and

the Amery Oasis was at least 650 m higher than present at b2.15 ± 0.15 Myr (Fig. 4a). The elevation of the ice sheet in the Fisher Massif remained relatively unchanged until b 0.5 Myr, but that of the Amery Oasis thinned to B400m between b2 and b 1 Myr.

Oberholzer et al.38 and Nicola et al.39 reported exposure ages from glacially rounded bedrock and erratics obtained from nunataks of the Deep Freeze Range (Supplementary Fig. 8). Complex exposure histories are not observed at any site until after b0.6 ± 0.1 Myr (Supplementary Fig. 15; Supplementary Table 3). At b 1.9 ±0.2 Myr, ice sheet thickness was at least 750-850 m higher than the present (Fig. 4b). The ice sheet at Mt. Keinath and Mt. Abbott thinned from b 1.9 ± 0.2 Myr. Exposure at the Black Ridge and Mt. Browning occurred in the Mid Pleistocene. The young exposure age (B11ka) of an erratic at a 950-m altitude has been interpreted to indicate potential shielding or post depositional overturning and should be interpreted with caution38. These results are consistent with weathered nunatak summits that are free of glacial deposits38,39.

Exposure age reconstruction of EAIS fluctuations. As illustrated in the above compilation, there is a wide range in reconstructed thickness across each of the sites of Areas 1-3, which indicates between 150 and 800 m increase in thickness relative to present, at least once before b 4 Myr. The variability in thickness is similar to the differing altitudes of each sampling location. For example, the maximum altitude of the Grove Mountains is b 200 m above the present height of the EAIS, making this region unable to record any additional increase in thickness. Individual ages appear to range over 1-2 Myr period over the Pliocene, with an average b 1-Myr age uncertainty. Except for the >4-Myr age, all error bars (± 1s) are overlapping such that we cannot discern any significant age differences. Nonetheless, we acknowledge that these are local sites distributed over a wide spatial area, and there is the possibility that some of these ages reflect local dynamics. Considered individually, however, each of the areas exhibits an overall decrease in thickness towards the present.

To explore the overall Plio-Pleistocene EAIS fluctuations, we consider conceptual ice sheet profiles for specific time slices. This exercise indicates that thinning likely began in inland Area 1 between 3 and 2.5 Myr, followed by a similar degree of thinning in Area 2 between 2.5 and 2 Myr. The onset of thinning in Area 3, proximal to the ice sheet edge, appears to have been delayed until b 2 Myr. Further thinning to near present levels likely occurred in all areas between 1 and 0.6 Myr.

Modelled reconstruction of EAIS fluctuations. We utilize the Ice sheet model for Integrated Earth system Studies (IcIES)42,43 (Fig. 5) to further examine Pliocene EAIS thickness and to independently assess the potential complications of constructing idealized ice sheet profiles from widely distributed exposure age data. Modern (pre-industrial) and Pliocene ice sheet experiments assume equilibrium response of the Antarctic ice sheet under their respective climatic conditions. The main climatic differences considered here are larger precipitation and higher temperature during the Pliocene, primarily due to increased CO2 concentration relative to preindustrial (405 versus 280 p.p.m., respectively)44 (see Methods). Modern initial conditions are based on the current ice-sheet configuration, while Pliocene conditions are approximated45 based on information indicating that the ice sheet retreated more often in the Pliocene than the late Pleistocene.

For the Pliocene, the model produces an b — 3-m SLE change in the WAIS volume. The ice sheet interior in East Antarctica is generally thicker, consistent with the individual exposure age sites

(for example, Figs 1 and 5), as well as with the conceptual ice sheet profiles. However, the ice sheets covering the Wilkes Subglacial and Aurora Basins thin and retreat inland by several hundred kilometres (Fig. 5c). Therefore, in total the modelled EAIS contribution to the Pliocene sea level is ~ 0 m SLE because increased interior thickness largely offsets the coastal retreat and collapse of large portions of the Wilkes Subglacial and Aurora Basins, where bed topography is relatively low in East Antarctica.


Thickening of the EAIS by up to 600 m before 3 Myr (Fig. 6) at first seems to contradict inferred relative high Pliocene sea level. On the basis of present ice volumes, complete melting of the Greenland Ice Sheet and the WAIS obtains B7 and 6 m SLE, respectively, indicating that an additional B7 m SLE would be required to reconcile a potential + 20 m Pliocene sea level2. Our model produces the expected B3-m SLE decrease in volume of the WAIS, but almost no SLE decrease in that of the EAIS, where interior thickening largely offsets the near complete melting of the Wilkes Subglacial and Aurora Basins (which are unsuitable for exposure dating due to the lack of unglaciated area; Fig. 5). This is similar to the coastal retreat and interior thickening currently observed by satellite in Greenland and Antarctica under anthropogenic warming46'47 and supports the IPCC AR5 interpretation that a value of + 20 m is a maximum on Milankovitch timescales (for example, greater than ~21 k.y., although it does not rule out shorter perturbations). The full range of reconstructed positive Pliocene sea level excursions spans 5-50 m (refs 44,48-50). Our results and the ± 10-m uncertainty of benthic d18O reconstructions indicates that the lower end is most likely.

Other major basins can be eliminated as large contributors to the Pliocene sea level on the basis of our model results, and increased thickness is observed in both our model- and nuclide-based reconstructions. Similar to the WAIS, which is inferred to

45° S

65° S

70° S

75° S

80° S

45° S

65° S

70° S

75° S

80° S


Figure 6 | Schematic diagram of ice sheet-climate interactions.

(a) The Early-Mid Pliocene (>3.0Myr) and (b) the Mid-Late Pleistocene (<1.0 Myr). Dashed line represents the Antarctic ice sheet surface during the early Pliocene.

have retreated during the Pliocene15, the Wilkes Subglacial Basin is grounded below the sea level. Geological and geochemical evidence indicate that its behaviour is also dynamic19,20, and a recent provenance study indicates that offshore Pliocene detrital sediment is sourced from Wilkes Land19. These results highlight the importance of bedrock elevation to ice sheet susceptibility to melting and retreat.

A thicker Pliocene EAIS was likely the result of enhanced moisture transport, resulting in increased upstream snowfall, as opposed to higher temperatures16, and is consistent with higher atmospheric pCO2 (refs 5,7,8), and global sea surface temperature

(H--2 to 3 °C)3-6. Higher SST and decreased sea ice would have

induced more evaporation from the ocean, increasing transport of water vapour to Antarctica that in turn created a thicker EAIS (Fig. 6a). This scenario is in line with the argument put forward by Prentice and Matthews51 that the sea-surface temperature of the Southern Ocean controls snowfall at Antarctica (snow gun hypothesis). Modern satellite observations indicate thickening of the interior EAIS in response to increased snowfall52.

The interpretation of more effective heat transport and a stronger hydrologic cycle in the Southern Ocean before 3 Myr is supported by results from Atmosphere-Ocean General Circulation Model experiments that indicate the global mean annual SST at 3 Myr was 3 °C higher than today53. A Recent model comparison study suggests that the uncertainty between different models is large in the high-latitude Southern Hemisphere54. In addition, General Circulation Model results indicate an overall 5% increase in precipitation during the Pliocene55, with an increase in the mean annual total precipitation rate (0.2-2 mm per day) in coastal Antarctic regions relative to today56.

After 3 Myr, EAIS thickness decreased as global climate cooled and Northern Hemisphere Ice Sheets expanded. We interpret this to be primarily caused by decreased SST that resulted in reduced moisture transport to Antarctica (Fig. 6b). Therefore, the net transport of humidity from the lower latitude ocean to Antarctica may have decreased, reducing snowfall28,30,57.

During the Pliocene, coastal areas of the EAIS retreated and the interior thickened, similar to modern ice sheets under anthropogenic warming. Some areas of the EAIS were up to 600 m thicker than present at least once before 3 Myr. It is plausible that increased interior thickness largely offsets the decrease from melting in Wilkes Land, primarily from the Wilkes Subglacial and Aurora Basins, resulting in little change in EAIS volume. Results presented here are consistent with the IPCC AR5 maximum estimate and suggest a B10-m Pliocene sea level increase on Milankovitch Cycle timescales.


Chemical procedures and measurements of cosmogenic nuclides. Chemical procedures to prepare samples for 10Be and 26Al measurements were on the basis of those of Kohl and Nishiizumi58, modified for processing gneiss samples. Chemical processing was performed in three steps: crushing of rock samples and purification of quartz, dissolution of quartz, and isolation and purification of beryllium and aluminium. Separation of 10Be and 26Al from rock samples was performed at the University of Tokyo.

Prepared samples were then measured using accelerator mass spectrometry (AMS) at the University of Tokyo59. The standards used were KNB5-1 for 10Be (10Be/9Be ratio 2.997 x 10 -11; (ref. 60)) and KNA4-2 for 26Al (26Al/27Al ratio 3.029 x 10 - 11; (ref. 61)). The background levels typically ranged from 1.32 x 10 - 14 to 4.66 x 10 - 14 for the 10Be/9Be measurements and from 5.0 x 10 - 15 to 6.91 x 10 - 15 for 26Al/27Al. Half-lives used were 1.51 Myr for 10Be (ref. 60) and 0.70 Myr for 26Al (ref. 61). The more recent value of 1.34 Myr 10Be half-life was not used because of the lack of raw data needed to recalculate the ages from older studies. Our conclusions are not significantly altered with a 1.34-Myr 10Be half-life. Nuclide production rates at each sampling site were calculated using the conventional production rates (5.1 ±0.3 atomsg- 1yr" "1 for 10Be and 31.1 ± 1.9 atomsg- 1yr- 1 for 26Al).

Antarctic atmospheric pressure. Because nuclide production rate scales with atmospheric thickness40, we accounted for the effect of potential differences in

Antarctic atmospheric pressure between the Pliocene and modern. (This effect is most pronounced when nuclide concentration is high.) We applied a scaling factor to nuclide concentrations reported here and, when appropriate, also to the compiled data. The scaling is on the basis of the assumption that the warm Pliocene atmospheric pressure was 2-3 hPa higher than the present (see Supplementary Table 4 and Supplementary Note 4). This results in a typically small (0.2%) change in the production rate that causes little change in the calculated exposure age. These final derived site production rates were employed to deduce exposure ages for each location.

IclES ice sheet model. The numerical experiments were performed using the Ice sheet model for IcIES, as described in refs 42,43), with several modifications explained below. Boundary conditions such as bedrock topography, surface temperature and accumulation were replaced with those provided by the Sea-level Response to Ice Sheet Evolution project62. The atmospheric lapse rate is assumed to be uniform at 8Kkm- 1. The enhancement factor and the geothermal heat flux were assumed to be uniform and, respectively, retuned to 3 and 54.6 mWm - 2 to simulate the volume of the present-day Antarctic ice sheet. In addition, the basal sliding parameterization was revised to allow submelt sliding following ref. 63. The same Weertman-type parameterization64 is adopted for basal sliding with the coefficient (14.4 x 10 - 10 m8 yr - 1N - which is eight times larger than that used in ref. 42. IcIES was run on a 40 x 40-km grid with 26 vertical layers. The time step was 0.5 year to stabilize the solution to larger basal velocities. Ice shelf processes are not included, and the ice sheet portion was only calculated with the shallow ice approximation, which is most important when the EAIS change is focused for the Pliocene situation. The lateral boundary is defined by floating conditions; ice that becomes afloat is cut off immediately65.

Two sets of experiments using IcIES were conducted, one for the 'modern' (Pre-Industrial) and the other for the 'Pliocene' situation using different climatic conditions and initial conditions, assuming the Antarctica ice sheet responds to equilibria under different climatic conditions. The climatic anomaly difference was obtained and simplified from the MIROC Pliocene experiment66, which was used to conduct the experimental set-up following the Pliocene model intercomparison project (PlioMIP44,67). The main climatic difference assumed in the Pliocene experiments of PlioMIP was in the atmospheric level of CO2, 405 versus 280 p.p.m. (Pliocene and pre-industrial level, respectively) as well as differing biome distributions, land sea mask and ice sheet configuration/topography on the basis of the PRISM3 data set44,68,69. Since the assumption in the PlioMIP experiment includes a different ice sheet shape (smaller in the Pliocene), the modelled climate of Atmosphere-Ocean General Circulation Model (PlioMIP) is strongly influenced by the assumed size and shape of ice sheet (retreated in the West Antarctica, Willkes Basin and Aurora Basin) and induces a large inhomogeneity of climatic conditions, that is, very warm and moist climate over the ice-retreated area through topographic and albedo feedback. To avoid this inhomogeneous climate change caused by these assumptions, we assume a homogeneous anomaly of temperature and precipitation rate by averaging the difference of MIROC between Pliocene and pre-industrial experiments. The temperature change due to topography change was subtracted. As a result, the climatic inputs for the ice sheet model are a temperature increase of 3 °C and a precipitation ratio increased by 50% in the Pliocene experiment compared with the pre-industrial one. The initial conditions for the modern and Pliocene ice sheets are those observed for present day (on the basis of BEDMAP2) and PRISM ice sheet from an approximation of the Pliocene ice sheet45,68. As a result, the Pliocene ice sheet extent is much larger than as initially specified, indicating that the model is not so sensitive to the specific initial conditions. Rather, the separation of Eastern and West Antarctica is more important. Each experiment was performed for 200,000 years in order to obtain steady-state solutions (Supplementary Fig. 16).


1. Salzmann, U. et a/. Challenges in quantifying Pliocene terrestrial warming revealed by data-model discord. Nat. Clim. Chang 3, 969-974 (2013).

2. Masson-Delmotte, V. et a/. Information from Paleoclimate Archives. in C/imate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. (eds Stocker, T. F. et a/.) (Cambridge University Press, 2013).

3. Dowsett, H. J. et al Micropaleontological evidence for increased meridonal heat-transport in the North-Atlantic Ocean during the Pliocene. Science 258, 1133-1135 (1992).

4. Crowley, T. J. Pliocene climates: the nature of the problem. Mar. Micropaleo 27, 3-12 (1996).

5. Raymo, M. E. et al. Mid-Pliocene warmth: stronger greenhouse and stronger conveyor. Mar. Micropaleo 27, 313-326 (1996).

6. Lunt, D. J. et al. Earth system sensitivity inferred from Pliocene modelling and data. Nat. Geosci. 3, 60-64 (2010).

7. Pagani, M. et al. High Earth-system climate sensitivity determined from Pliocene carbon dioxide concentrations. Nat. Geosci. 3, 27-30 (2010).

8. Seki, O. et al. Alkenone and boron-based Pliocene pCO2 records. Earth Planet. Sci. Lett. 292, 201-211 (2010).

9. Dowsett, H. J. & Cronin, T. M. High eustatic sea level during the middle Pliocene: evidence from the southeastern U.S. Atlantic Coastal Plain. Geology 18, 435-438 (1990).

10. Krantz, D. E. A chronology of Pliocene sea-level fluctuations: the U.S. Middle Atlantic Coastal Plain record. Quat. Sci. Rev. 10, 163-174 (1991).

11. Wardlaw, B. R. & Quinn, T. M. The record of Pliocene sea-level change at Enewetak Atoll. Quat. Sci. Rev. 10, 247-258 (1991).

12. Kaufman, D. S. & Brigham-Grette, J. Aminostratigraphic correlations and paleotemperature implications, Pliocene-Pleistocene high-sea-level deposits, northwestern Alaska. Quat. Sci. Rev. 12, 21-33 (1993).

13. Miller, K. G. et al. The Phanerozoic record of global sea-level change. Science 310, 1293-1298 (2005).

14. Lunt, D. J. et al. Late Pliocene Greenland glaciations controlled by a decline in atmospheric CO2 levels. Nature 454, 1102-1105 (2008).

15. Naish, T. et al. Obliquity-paced Pliocene West Antarctic ice sheet oscillations. Nature 458, 322-328 (2009).

16. Pollard, D. & DeConto, R. M. Modelling West Antarctic ice sheet growth and collapse through the past five million years. Nature 458, 329-332 (2009).

17. Hill, D. J. et al. Characterizing ice sheets during the Pliocene: evidence from data and models. in Deep-Time Perspectives on Climate Change: Marrying the Signal from Computer Models and Biological Proxies (eds Williams, M. et al.) (Geological Society of London, 2007).

18. Dolan, A. M. et al. Sensitivity of Pliocene ice sheets to orbital forcing. Palaeogeogr. Palaeoclimatol. Palaeoecol. 309, 98-110 (2011).

19. Cook, C. P. et al. Dynamic behaviour of the East Antarctic ice sheet during Pliocene warmth. Nat. Geosci. 6, 765-769 (2013).

20. Patterson, M. O. et al. Orbital forcing of the East Antarctic ice sheet during the Pliocene and Early Pleistocene. Nat. Geosci. 7, 841-847 (2014).

21. Webb, P.-N. et al. Cenozoic marine sedimentation and ice volume variation on the East Antarctic craton. Geology 12, 287-291 (1984).

22. McKelvey, B. C. et al. The Pagodroma Group—a Cenozoic record of the East Antarctic ice sheet in the northern Prince Charles Mountains. Antarc. Sci. 13, 455-468 (2001).

23. O'Brien, P. E. et al. Late neogene ice drainage changes in Prydz Bay, East Antarctica and the interaction of Antarctic Ice Sheet evolution and climate. Palaeogeogr. Palaeoclimatol. Palaeoecol. 245, 390-410 (2007).

24. Moriwaki, K., Hirakawa, K. & Matsuoka, N. Weathering stage of till and glacial history of the central Sor Rondane Mountains, East Antarctica. Proc. NIPR Symp. Antarc. Geosci 5, 99-111 (1991).

25. Nishiizumi, K. et al. Cosmic ray produced 10Be and 26Al in Antarctic rocks: exposure and erosion history. Earth Planet. Sci. Lett. 104, 440-454 (1991).

26. Matsuoka, N. et al. Quaternary bedrock erosion and landscape evolution in the Sor Rondane Mountains, East Antarctica: reevaluating rates and processes. Geomorphology 81, 408-420 (2006).

27. Huang, F. et al. Fluctuation history of the interior East Antarctic Ice Sheet since mid-Pliocene. Antarc. Sci. 20, 197-203 (2008).

28. Kong, P. et al. Late Miocene ice sheet elevation in the Grove Mountains, East Antarctica, inferred from cosmogenic 21Ne-10Be-26Al. Glob. Planet. Chan.g 72, 50-54 (2010).

29. Lilly, K. et al. Pleistocene dynamics of the interior East Antarctic ice sheet. Geology 38, 703-706 (2010).

30. Liu, X. et al. History of ice sheet elevation in East Antarctica: Paleoclimatic implications. Earth Planet. Sci. Lett. 290, 281-288 (2010).

31. Strasky, S. et al. Surface exposure ages imply multiple low-amplitude Pleistocene variations in East Antarctic Ice Sheet, Ricker Hills, Victoria Land. Antarc. Sci. 21, 59-69 (2009).

32. Tschudi, S. et al. Surface exposure dating of Sirius formation at Allan Hills nunatak, Antarctica: new evidence for long-term ice-sheet stability. Eclogae Geol. Helv. 96, 109-114 (2003).

33. Brown, E. T. et al. Examination of surface exposure ages of Antarctic moraines using in situ produced 10Be and 26Al. Geochim. Cosmochim. Acta 55, 2269-2283 (1991).

34. Brook, E. J. et al. Chronology of Taylor Glacier Advances in Arena Valley, Antarctica, Using in situ Cosmogenic 3He and 10Be. Quat. Res. 39, 11-23 (1993).

35. Brook, E. J. et al. Constraints on age, erosion, and uplift of Neogene glacial deposits in the Transantarctic Mountains determined from in situ cosmogenic 10Be and 26Al. Geology 23, 1063-1066 (1995).

36. Fogwill, C. J. et al. Cosmogenic nuclides 10Be and 26Al imply limited Antarctic Ice Sheet thickening and low erosion in the Shackleton Range for > 1 m.y. Geology 32, 265-268 (2004).

37. Fink, D. et al. Pleistocene deglaciation chronology of the Amery Oasis and Radok Lake, northern Prince Charles Mountains, Antarctica. Earth Planet. Sci. Lett. 243, 229-243 (2006).

38. Oberholzer, P. et al. Limited Pliocene/Pleistocene glaciations in Deep Freeze Range, northern Victoria Land, Antarctica, derived from in situ cosmogenic nuclides. Antarc. Sci. 15, 493-502 (2003).

39. Nicola, L. D. et al. Multiple cosmogenic nuclides document complex Pleistocene exposure history of glacial drifts in Terra Nova Bay (northern Victoria Land, Antarctica). Quat. Res. 71, 83-92 (2009).

40. Stone, J. O. Air pressure and cosmogenic isotope production. J. Geophys. Res. 105, 23753-23759 (2000).

41. Hofle, H.-C. & Buggisch, W. Glacial geology and petrography of erratics in the Shackleton Range, Antarctica. Polarforschung 63, 183-201 (1995).

42. Saito, F. & Abe-Ouchi, A. Modelled response of the volume and thickness of the Antarctic ice sheet to the advance of the grounded area. Ann. Glaciol. 51, 41-48 (2010).

43. Abe-Ouchi, A. et al. Insolation-driven 100,000-year glacial cycles and hysteresis of ice-sheet volume. Nature 500, 190-193 (2013).

44. Haywood, A. M. et al. Pliocene Model Intercomparison Project (PlioMIP): experimental design and boundary conditions (Experiment 1). Geosci. Model Dev. 3, 227-242 (2010).

45. Dolan, A. M. et al. Pliocene Ice Sheet Modelling Intercomparison Project (PLISMIP)—experimental design. Geosci. Model Dev. 5, 963-974 (2012).

46. Krabill, W. et al. Greenland Ice Sheet: increased coastal thinning. Geophys. Res. Lett. 31, L24402 (2004).

47. Pritchard, H. D. et al. Extensive dynamic thinning on the margins of the Greenland and Antarctic ice sheets. Nature 461, 971-975 (2009).

48. Dwyer, G. S. & Chandler, M. A. Mid-Pliocene sea level and continental ice volume based on coupled benthic Mg/Ca palaeotemperatures and oxygen isotopes. Philos. Trans. R. Soc. A 367, 157-168 (2009).

49. Naish, T. R. & Wilson, G. S. Constraints on the amplitude of Mid-Pliocene (3.6-2.4 Ma) eustatic sea-level fluctuations from the New Zealand shallow-marine sediment record. Philos. Trans. R. Soc. A 367, 169-187 (2009).

50. Sosdian, S. & Rosenthal, Y. Deep-sea temperature and ice volume changes across the Pliocene-Pleistocene climate transitions. Science 325, 306-310 (2009).

51. Prentice, M. L. & Matthews, R. K. Tertiary ice sheet dynamics: the snow gun hypothesis. J. Geophys. Res. 96, 6811-6827 (1991).

52. Davis, C. H. et al. Snowfall-driven growth in East Antarctic ice sheet mitigates recent sea-level rise. Science 308, 1898-1901 (2005).

53. Haywood, A. M. et al. Large-scale features of Pliocene climate: results from the Pliocene Model Intercomparison Project. Clim. Past 9, 191-209 (2013).

54. Hill, D. J. et al. Evaluating the dominant components of warming in Pliocene climate simulations. Clim. Past 10, 79-90 (2014).

55. Sloan, L. C., Crowley, T. J. & Pollard, D. Modeling of middle Pliocene climate with the NCAR GENESIS general circulation model. Mar. Micropaleo 27, 51-61 (1996).

56. Haywood, A. M. et al. Antarctic climate during the middle Pliocene: model sensitivity to ice sheet variation. Palaeogeogr. Palaeoclimatol. Palaeoecol. 182, 93-115 (2002).

57. Altmaier, M. et al. Glaciation history of Queen Maud Land (Antarctica) reconstructed from in-situ produced cosmogenic 10Be, 26Al and 21Ne. Polar Sci. 4, 42-61 (2010).

58. Kohl, C. P. & Nishiizumi, K. Chemical isolation of quartz for measurement of in-situ-produced cosmogenic nuclides. Geochim. Cosmochim. Acta 56, 3583-3587 (1992).

59. Matsuzaki, H. et al. Current status and future direction of MALT, The University of Tokyo. Nucl. Instrum. Methods Phys. Res. B 223-224, 92-99 (2004).

60. Nishiizumi, K. et al. Absolute calibration of 10Be AMS standards. Nucl. Instrum. Methods Phys. Res. B 258, 403-413 (2007).

61. Nishiizumi, K. Preparation of 26Al AMS Standards. Nucl. Instrum. Methods Phys. Res. B 223-224, 388-392 (2003).

62. Bindschadler, R. A. et al. Ice-sheet model sensitivities to environmental forcing and their use in projecting future sea level (the SeaRISE project). J. Glaciol. 59, 195-224 (2013).

63. Greve, R. Relation of measured basal temperatures and the spatial distribution of the geothermal heat flux for the Greenland ice sheet. Ann. Glaciol. 42, 424-432 (2005).

64. Weertman, J. The theory of glacier sliding. J. Glaciol. 5, 287-303 (1964).

65. Chan, W. L., Abe-Ouchi, A. & Ohgaito, R. Simulating the mid-Pliocene climate with the MIROC general circulation model: experimental design and initial results. Geosci. Model Dev. 4, 1035-1049 (2011).

66. Payne, A. J. A thermomechanical model of ice flow in West Antarctica. Clim. Dyn. 15, 115-125 (1999).

67. Dowsett, H. J. et al. An assessment of confidence in Pliocene global sea-surface temperature. Nat. Clim. Chang. 2, 365-371 (2012).

68. Hill, D. J. Modelling Earth's Cryosphere during Peak Pliocene Warmth (Ph.D. Thesis, University of Bristol 368, 2009).

69. Dowsett, H. J. et al. The PRISM3D paleoenvironmental reconstruction. Stratigraphy 7, 123-139 (2010).

70. Fretwell, P. et al. Bedmap2: improved ice bed, surface and thickness datasets for Antarctica. Cryosphere 7, 375-393 (2013).


We thank Yosuke Miyairi for technical support during the AMS measurements. This research was partly supported by the JSPS NEXT programme (GR031 to Y.Y.) and KAKENHI (26247085 to Y.Y.).

Author contributions

Y.Y., M.Y. and A.A.-O. designed the research; K.M. collected the rock samples; M.Y. and H.M. carried out AMS measurements; M.Y. compiled and recalculated exposure ages. All authors contributed to the interpretation of the data. A.A.-O. and F.S. performed model simulations. S.O., M.Y., Y.Y. and A.A.-O. wrote the paper together with input from all other authors.

Additional information

Supplementary Information accompanies this paper at naturecommunications

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

Reprints and permission information is available online at reprintsandpermissions/

How to cite this article: Yamane, M. et al. Exposure age and ice-sheet model constraints on Pliocene East Antarctic ice sheet dynamics. Nat. Commun. 6:7016 doi: 10.1038/ ncomms8016 (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 article's 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