Articles | Volume 13, issue 11
https://doi.org/10.5194/tc-13-2935-2019
https://doi.org/10.5194/tc-13-2935-2019
Research article
 | Highlight paper
 | 
08 Nov 2019
Research article | Highlight paper |  | 08 Nov 2019

New Last Glacial Maximum ice thickness constraints for the Weddell Sea Embayment, Antarctica

Keir A. Nichols, Brent M. Goehring, Greg Balco, Joanne S. Johnson, Andrew S. Hein, and Claire Todd
Abstract

We describe new Last Glacial Maximum (LGM) ice thickness constraints for three locations spanning the Weddell Sea Embayment (WSE) of Antarctica. Samples collected from the Shackleton Range, Pensacola Mountains, and the Lassiter Coast constrain the LGM thickness of the Slessor Glacier, Foundation Ice Stream, and grounded ice proximal to the modern Ronne Ice Shelf edge on the Antarctic Peninsula, respectively. Previous attempts to reconstruct LGM-to-present ice thickness changes around the WSE used measurements of long-lived cosmogenic nuclides, primarily 10Be. An absence of post-LGM apparent exposure ages at many sites led to LGM thickness reconstructions that were spatially highly variable and inconsistent with flow line modelling. Estimates for the contribution of the ice sheet occupying the WSE at the LGM to global sea level since deglaciation vary by an order of magnitude, from 1.4 to 14.1 m of sea level equivalent. Here we use a short-lived cosmogenic nuclide, in situ-produced 14C, which is less susceptible to inheritance problems than 10Be and other long-lived nuclides. We use in situ 14C to evaluate the possibility that sites with no post-LGM exposure ages are biased by cosmogenic nuclide inheritance due to surface preservation by cold-based ice and non-deposition of LGM-aged drift. Our measurements show that the Slessor Glacier was between 310 and up to 655 m thicker than present at the LGM. The Foundation Ice Stream was at least 800 m thicker, and ice on the Lassiter Coast was at least 385 m thicker than present at the LGM. With evidence for LGM thickening at all of our study sites, our in situ 14C measurements indicate that the long-lived nuclide measurements of previous studies were influenced by cosmogenic nuclide inheritance. Our inferred LGM configuration, which is primarily based on minimum ice thickness constraints and thus does not constrain an upper limit, indicates a relatively modest contribution to sea level rise since the LGM of < 4.6 m, and possibly as little as < 1.5 m.

1 Introduction

We describe new constraints on Last Glacial Maximum (LGM, ca. 26 to 15 ka; Peltier and Fairbanks, 2006) ice thickness changes from three locations within the Weddell Sea Embayment (WSE) of Antarctica (Fig. 1). The WSE drains approximately one-fifth of the total area of the Antarctic ice sheets (AISs) (Joughin et al., 2006) and is thus an important contributor to LGM-to-present and, potentially, future sea level change. Previous attempts to reconstruct LGM-to-present ice thickness changes around the WSE used measurements of long-lived cosmogenic nuclides, primarily 10Be (half-life 1.387±0.012 Ma; Chmeleff et al., 2010; Korschinek et al., 2010) and 26Al (half-life 705±17 ka; Norris et al., 1983), sourced from bedrock and erratic cobbles proximal to modern glacier surfaces. Through measuring the cosmogenic nuclide concentration of samples of glacial deposits and bedrock, one can constrain the magnitude and timing of past changes in the thickness of adjacent ice masses. However, an absence of post-LGM apparent exposure ages at many sites around the WSE led to LGM thickness reconstructions that were spatially highly variable and inconsistent with flow line modelling (e.g. Whitehouse et al., 2017). Consequently, estimates based on ice models constrained by field evidence (Le Brocq et al., 2011) and by relative sea level records and earth viscosity models (Bassett et al., 2007) for the contribution of the sector to global sea level since deglaciation began vary by an order of magnitude, from 1.4 to 14.1 m, respectively. The lack of geological evidence for LGM thickening is also manifest in a misfit between present-day geodetic uplift rate measurements in southern Palmer Land and predicted uplift rates from a glacial isostatic adjustment (GIA) model (Wolstencroft et al., 2015). Constraining the previous vertical extent of ice provides inputs to numerical models investigating both the response of the ice sheet to past and potential future changes in climate and sea level (e.g. Briggs et al., 2014; Pollard et al., 2016, 2017; Whitehouse et al., 2017) as well as the response of the solid earth to past ice load changes to quantify present-day ice-mass loss (e.g. Wolstencroft et al., 2015). Furthermore, quantifying the LGM dimensions of the WSE sector of the AIS is required to further constrain the offset between estimates for post-LGM sea level rise and estimates of the total amount of ice melted since the LGM. The former is sourced from sea level index points, and the latter is sourced from our knowledge of the dimensions of ice masses at the LGM (Simms et al., 2019). Currently, the “missing ice” accounts for between 15.6±9.6 and 18.1±9.6 m of global sea level rise since the LGM (Simms et al., 2019).

https://www.the-cryosphere.net/13/2935/2019/tc-13-2935-2019-f01

Figure 1The Weddell Sea Embayment, including all locations referred to within the text. SH, WH, and TH are the Schmidt, Williams, and Thomas Hills, respectively. FH, P/M, and MH are the Flower Hills, Patriot and Marble Hills, and the Meyer Hills, respectively. Black is exposed rock. Red boxes show extent of satellite images in Fig. 4. Exposed rock and coastline sourced from the SCAR Antarctic Digital Database. Bathymetry sourced from the International Bathymetric Chart of the Southern Ocean V1.0 (IBSCO; Arndt et al., 2013). Surface topography (shading) is sourced from the Reference Elevation Model of Antarctica (REMA; Howat et al., 2019). PS1423-2 is a marine sediment core from Crawford et al. (1996).

Although the use of cosmogenic nuclide geochronology to study the AIS is clearly proven (e.g. Stone et al., 2003; Ackert et al., 2007), applications in the WSE are challenging. Many studies, despite making multiple cosmogenic nuclide measurements from relatively large numbers of samples, observed no or few post-LGM exposure ages (Hein et al., 2011, 2014; Balco et al., 2016; Bentley et al., 2017). With no evidence for LGM ice cover, it was not clear whether sites were covered by ice at the LGM, or whether sites were covered but the ice left no fresh deposits on top of those yielding pre-LGM ages. It is therefore currently unknown whether ice was thicker than present during the LGM at the Schmidt Hills in the Pensacola Mountains and in the Shackleton Range (Figs. 1 and 2). Results from the Schmidt Hills (Fig. 2) indicating no LGM thickening of the Foundation Ice Stream (FIS) are particularly problematic, as thickening of 500 m from the Williams Hills, 50 km upstream of the Schmidt Hills, produces a LGM surface slope that is steeper than glaciological models permit and is also steeper than present-day ice surface slopes (Balco et al., 2016). Cold-based ice and an associated lack of subglacial erosion is the likely cause of the complex 10Be data sets, evidenced by numerous studies in the WSE that report 10Be and 26Al ratios significantly below those predicted for continuous exposure which is indicative of significant periods of non-erosive burial (e.g. Bentley et al., 2006; Sugden et al., 2017). Cold-based ice preserves surfaces (e.g. Stroeven et al., 2002; Sugden et al., 2005; Gjermundsen et al., 2015), allowing nuclide concentrations to persist within surfaces from previous periods of exposure to the present, a phenomenon known as inheritance. Long-lived nuclides are particularly susceptible to inheritance due to their long half-lives which, when protected from erosion beneath cold-based ice, require long periods of burial to reduce concentrations to below measurable levels. When covered by cold-based ice during glaciations, concentrations of long-lived nuclides record exposure during multiple separate ice-free periods rather than just the most recent one. Inheritance thus hinders interpretations of cosmogenic nuclide measurements.

https://www.the-cryosphere.net/13/2935/2019/tc-13-2935-2019-f02

Figure 2Current terrestrial ice thickness constraints inferred from measurements of long-lived nuclides around the WSE. Acronyms are as in Fig. 1. Constraints for the SH, WH, and TH are sourced from Balco et al. (2016) and Bentley et al. (2017). MB is Mount Bragg (Bentley et al., 2017). Thickness estimate for the Dufek Massif (DM) is sourced from Hodgson et al. (2012). Constraints for the P/M are sourced from Hein et al. (2016). For the MH and FH, the LGM thickness constraints are sourced from Fogwill et al. (2014). The thickness constraints sourced from Fogwill et al. (2014) were interpreted using modern ice surface elevations for the Rutford Ice Stream and Union Glacier measured using the Reference Elevation Model of Antarctica (REMA; Howat et al., 2019). Thickness constraints for the Shackleton Range are sourced from Hein et al. (2011, 2014). The range of LGM thicknesses for the Behrendt Mountains are sourced from multiple locations (Bentley et al., 2006).

https://www.the-cryosphere.net/13/2935/2019/tc-13-2935-2019-f03

Figure 3(a) Hypothetical ice surface elevation change at a nunatak partially covered by ice at the LGM. (b) Resulting in situ 14C concentration, assuming no surface erosion, in samples collected at 100 m intervals along an elevation transect on the surface of the nunatak. Thin black lines indicate isochrons of exposure duration. Thick black line with grey shading represents the saturation concentration and associated error envelope. Error envelope represents typical analytical uncertainty. “True exposure” refers to the resulting 14C concentration associated with the ice surface change history on the left plot. “Apparent exposure” is the resulting concentration that includes an inherited component, which is a residual 14C inventory remaining from the hypothetical samples which were saturated prior to ice cover.

Download

We resolve conflicting LGM thickening estimates based on 10Be measurements by using measurements of in situ-produced 14C, a cosmogenic nuclide that is, owing to a short half-life of 5730 years, largely insensitive to inheritance. We present the in situ 14C analysis of transects of erratic and bedrock samples from the Shackleton Range, Lassiter Coast, and Pensacola Mountains (Fig. 1). Our results constrain the LGM thickness of the Slessor Glacier to between 310 and up to 655 m. We show that ice was at least 385 m thicker than present during the LGM at the Lassiter Coast, proximal to the modern Ronne Ice Shelf edge. Our data also constrain the LGM thickness of the FIS to at least 800 m at the Schmidt Hills. Replicate measurements made from four samples revealed higher-than-expected variability of in situ 14C measurements, which is discussed in Sect. 4.1. Our thickness estimates are comparable to those of Hein et al. (2016) in the Ellsworth Mountains, as well as those of Balco et al. (2016) and Bentley et al. (2017) in the Williams and Thomas Hills. Although our results show that locations around the WSE were buried by hundreds of metres of ice, this is less than called for by some reconstructions. Our inferred LGM configuration, which is primarily based on minimum ice thickness constraints and thus does not constrain an upper limit, indicates a relatively modest contribution to sea level rise since the LGM of < 4.6 m, and possibly as little as < 1.5 m.

1.1 The Last Glacial Maximum in the Weddell Sea Embayment

Although it is clear that grounded ice in the WSE has been thicker in the past (Bentley and Anderson, 1998), there is little evidence as to the thickness and grounding line position of the ice sheet at the LGM, with contrasting evidence from marine sources, and those inferred from terrestrial studies (Hillenbrand et al., 2014). Terrestrial evidence for the extent of ice in the WSE during the LGM takes the form of numerous cosmogenic nuclide studies. Bentley et al. (2006) measured the 10Be and 26Al content of erratics on the southern Antarctic Peninsula. Studies report cosmogenic nuclide concentrations from the Meyer Hills, Patriot Hills, Marble Hills, and the Flower Hills, all in the Ellsworth Mountains (Bentley et al., 2010; Fogwill et al, 2014; Hein et al., 2016; Sugden et al; 2017), the Pensacola Mountains (Hodgson et al., 2012; Balco et al, 2016; Bentley et al., 2017), and the Shackleton Range (Fogwill et al., 2004; Hein et al., 2011, 2014). Figure 2 summarises the ice thickness estimates from these studies. The majority of estimates are sourced from 10Be measurements, with some accompanying 26Al measurements. Two exceptions are Fogwill et al. (2014) and Balco et al. (2016), who combined some in situ 14C measurements with 10Be measurements to constrain the thickness of the Rutford and Institute ice streams and the Foundation Ice Stream, respectively. The highest elevation post-LGM exposure ages at each site delineate the minimum vertical extent of ice at the LGM. Ice thickness estimates vary spatially around the embayment, ranging from zero to hundreds of metres of LGM thickening.

Marine geological and geophysical evidence in the southern Weddell Sea indicates a significantly expanded WSE LGM configuration, with subglacial till, subglacial bedforms, and a grounding zone wedge found towards the shelf edge (Hillenbrand et al., 2012, 2014; Larter et al., 2012; Arndt et al., 2017). As a result, there is currently a disconnect between marine evidence for a greatly expanded WSE sector and terrestrial evidence indicating little to no vertical change at the LGM in some areas. Hillenbrand et al. (2014) propose two potential LGM configurations of the WSE sector of the AIS. The first scenario, based on terrestrial evidence for vertical LGM ice thicknesses, involves a complex configuration with the grounding line of the ice sheet situated towards the shelf edge and a largely ice-free Filchner Trough and western margin of the WSE. The second scenario, based on marine evidence, places the grounding line of the ice sheet at the shelf edge across the width of the WSE. Flow line modelling of the response of the FIS, which occupied the Filchner Trough at the LGM, to the onset of glacial conditions shows that there are two plausible LGM grounding line positions for the ice stream: one situated at the shelf edge and another at the northern margin of Berkner Island (Whitehouse et al., 2017).

1.2 In situ 14C exposure dating

Cosmogenic nuclides 10Be and 26Al have half-lives that are much longer than glacial–interglacial cycles, so 10Be and 26Al concentrations produced in previous interglacials persist to the present if buried by non-erosive, cold-based ice. The short half-life of in situ 14C means that only short periods of burial are required to significantly reduce concentrations from previous periods of exposure, making in situ 14C less sensitive to inheritance than longer-lived nuclides. For example, a burial duration beneath non-erosive, cold-based ice of 11 kyr results in ca. 74 % of the original in situ 14C decaying away. Furthermore, continuously exposed, slowly eroding surfaces reach an equilibrium between production and decay of in situ 14C (“saturation”) after approximately 30 to 35 kyr. A sample that has reached saturation thus requires low erosion and continuous exposure from before the LGM, whilst a sample that yields a concentration below saturation requires ice cover during the last ca. 35 kyr. Surfaces yielding saturation concentrations therefore provide an upper limit on LGM thickening. Figure 3 shows a hypothetical ice surface elevation change history at a nunatak partially buried by cold-based, non-erosive ice during the LGM, with associated in situ 14C measurements from samples collected along an elevation transect on the surface of the nunatak. There is a transition from undersaturated to saturated samples, a discontinuity in the 14C concentrations which constrains the LGM ice thickness. The “true exposure” data points represent in situ 14C concentrations with resulting exposure ages matching the post-LGM ice-surface-lowering history. The “apparent exposure” data points were saturated at the onset of ice cover and include in situ 14C that persists to the present due to an insufficient amount of time passing for it to decay away. For the five undersaturated samples, which were buried by ice for differing durations, a range of ∼2 to ∼4 % of the 14C accumulated prior to burial will persist to the present. In terms of the effect on resulting exposure ages, the sample exposed at 10 ka yields an apparent exposure age of 11.41 ka (∼13 % increase), and the sample exposed at 2 ka yields an apparent exposure age of 2.17 ka (∼8 % increase). Without knowing the burial duration of the samples or whether or not the samples were indeed saturated upon burial by LGM ice, we do not know the exact quantity of in situ 14C inherited in the samples. The in situ 14C exposure ages are therefore maximum deglaciation ages. In the same hypothetical scenario with the same samples, ca. 98 % and 97 % of the 10Be and 26Al accumulated prior to burial will persist to the present, respectively.

We report in situ 14C concentrations measured from both erratic and bedrock samples, with primarily erratic samples from the Shackleton Range and the Pensacola Mountains, and solely bedrock from the Lassiter Coast. We assume both materials provide the same information regarding the timing of ice retreat and constraining LGM ice thicknesses. For example, we assume that both erratics and bedrock samples saturated with in situ 14C indicate that their respective sampling locations were ice-free for the last 30 to 35 kyr. With the exception of two samples, all of our erratic samples have previously been measured for their 10Be content (Hein et al., 2011, 2014; Balco et al., 2016), with the vast majority yielding ages far in excess of the LGM. It is highly likely that these erratic samples have been repeatedly covered and exposed by cold-based ice. Having been covered and uncovered in situ, the erratic samples can thus effectively be considered bedrock. There are, however, potential situations where our assumption that bedrock and erratic samples provide the same information with respect to the timing of changes in ice thickness is not met and resulting 14C concentrations misrepresent the age of deglaciation, creating scatter in the measured in situ 14C data.  Erratic samples may, for example, be sourced from mass movement onto glacier surfaces, producing spuriously high 14C concentrations (see Balco et al., 2019). Spuriously high, in excess of saturation, in situ 14C concentrations sourced from bedrock samples, however, can only result from analytical errors and thus provide an important test for the premise of the technique. Additionally, erratic cobbles may have undergone downslope movement post-deposition and may have flipped over, or may have been subjected to high erosion rates, which could produce in situ 14C concentrations with resulting exposure ages lower than the true age of deglaciation. Snow shielding of sample locations is another mechanism leading to exposure ages which underestimate the age of deglaciation and can influence both bedrock and erratic samples. Whilst not without challenges, our in situ 14C measurements provide an opportunity to unambiguously show whether sites around the WSE were covered by ice at the LGM.

1.3 Sample sites

1.3.1 Shackleton Range

The Shackleton Range is located in Coats Land in northeastern WSE, adjacent to Slessor Glacier (Figs. 1 and 4a). Slessor Glacier drains ice from the East Antarctic Ice Sheet (EAIS) into the Filchner Ice Shelf. Mt Skidmore is located approximately 25 km upstream of the modern Slessor Glacier grounding line, with the Köppen and Stratton glaciers respectively joining the Slessor Glacier to the north and south of Mt Skidmore (Fig. 4a). Proximal to sampling locations are Ice Tongue A and Ice Tongue B of the Stratton Glacier, as well as the Snow Drift Glacier (Supplement Fig. S1). We assume that samples collected from Mt Skidmore record changes in the thickness of the Slessor Glacier. However, it is possible that samples collected proximal to the smaller ice masses may have been buried by them, rather than by the Slessor Glacier, potentially complicating the interpretation of results. The modern Slessor Glacier surface is situated at ∼200 m a.s.l. proximal to Mt Skidmore, with exposed surfaces of Mt Skidmore located up to over ∼820 m a.s.l. Mt Provender is located adjacent to the Slessor Glacier grounding line and is bounded by the Stratton and Blaiklock glaciers to the north and south, respectively. Exposed rock of Mt Provender rises from the modern ice surface up to over ∼900 m a.s.l. We analysed 11 samples from the Shackleton Range (Table S1 in the Supplement), with two from Mt Provender and nine from Mt Skidmore (Fig. 4a). At Mt  Provender we analysed one erratic cobble from near the modern ice surface and one bedrock sample from ∼650 m above it (Fig. S2). Samples from Mt Skidmore include one bedrock sample and eight cobbles that form an elevation transect from near the modern ice surface to ∼300 m above it (Fig. S1). The two highest elevation samples collected from Mt Skidmore are proximal to the main trunk of the Stratton Glacier more so than the Slessor Glacier and were collected from ca. 115 and 130 m above the modern Stratton Glacier surface. The two highest elevation samples on Mt Skidmore therefore may represent a Stratton Glacier ice surface lowering more so than the Slessor Glacier and thus are presented as a separate sample group to those collected proximal to the Slessor Glacier.

https://www.the-cryosphere.net/13/2935/2019/tc-13-2935-2019-f04

Figure 4Landsat imagery of study sites. Location of each image is shown in Fig. 1. Green dots show sample locations. Arrows show ice flow directions. (a) Mt Skidmore and Mt Provender, Shackleton Range. (b) Lassiter Coast, southern Palmer Land. (c) Schmidt Hills, Pensacola Mountains. (d) Thomas Hills, Pensacola Mountains. Landsat 8 imagery courtesy of the U.S. Geological Survey. Grounding line positions sourced from the MEaSUREs program V2 (Rignot et al., 2011, 2014, 2016).

1.3.2 Lassiter Coast

The Lassiter Coast is located on the east coast of southern Palmer Land, adjacent to the present position of the Ronne Ice Shelf edge (Fig. 1). The modern ice surface is situated at 490 m a.s.l. Johnson et al. (2019) collected samples from several sites in this area (Fig. 4b) and carried out 10Be measurements; we subsequently carried out 14C measurements on these samples as part of the present study, and the 14C results are reported both here and in Johnson et al. (2019). Here we discuss results for a total of eight bedrock samples from Mt Lampert and the Bowman Peninsula collected from 20 to 385 m above the modern ice surface (Figs. 4b and S3); see Table S1 for sample data and Johnson et al. (2019) for 10Be measurements. The adjacent Johnston Glacier drains ice from central Palmer Land into the WSE (Fig. 4b). We interpret the samples together as effectively a single elevation transect that records changes in the thickness of grounded ice in the WSE immediately east of these sites after the LGM.

1.3.3 Pensacola Mountains

The Schmidt Hills are a series of nunataks adjacent to the FIS in the southeast WSE (Figs. 1 and 4c) proximal to the modern grounding line. The FIS is a major ice stream that drains ice from both the EAIS and West Antarctic Ice Sheet (WAIS) into the WSE. The surface of the FIS adjacent to the Schmidt Hills is situated ca. 200 m a.s.l., with exposed surfaces of the Schmidt Hills reaching up to 1100 m a.s.l. The Thomas Hills are another series of nunataks adjacent to the FIS, located ∼130 km upstream of the Schmidt Hills (Figs. 1, 4d). The main trunk of the FIS adjacent to the Thomas Hills is near 550 m a.s.l., with the Thomas Hills rising up to 1050 m a.s.l. The local ice margin of the FIS at the Thomas Hills is situated ∼75 m below the centre of the FIS. We analysed 17 samples from the Pensacola Mountains (Table S1), 15 from the Schmidt Hills, and 2 from the Thomas Hills. We made a further seven repeat measurements from four samples collected from the Schmidt Hills. Samples from the Schmidt Hills were collected from Mount Coulter and No Name Spur (Figs. 4c and S4) from close to the modern ice surface to approximately 800 m above it. We also analysed two samples from the Thomas Hills which were collected from Mount Warnke ca. 320 m above the FIS ice margin (Figs. 4d and S5). The highest elevation sample from the Schmidt Hills, collected from ca. 1035 m a.s.l., is the only bedrock sample analysed from the Pensacola Mountains, with the rest being erratic cobbles.

2 Methods

We used between 0.5 and 10 g of quartz from each sample for in situ 14C analysis. The methodology used for the isolation of quartz varies for samples from different sample sites because quartz was previously isolated for prior cosmogenic nuclide studies (see Hein et al., 2011; Balco et al., 2016). For samples processed at the Tulane University Cosmogenic Nuclide Laboratory (primarily those from the Lassiter Coast), quartz was isolated through crushing, sieving, magnetic separation and froth flotation (modified from Herber, 1969) of sample material. Samples were then etched for at least two periods of 24 h on both a shaker table in 5 % HF∕HNO3 and then in an ultrasonic bath in 1 % HF∕HNO3. This leaching procedure removes the organic compound laurylamine used in the froth flotation procedure (Nichols and Goehring, 2019) that could otherwise potentially contaminate our samples with modern carbon.

Carbon was extracted using the Tulane University Carbon Extraction and Graphitization System (TU-CEGS), following the method of Goehring et al. (2019). Quartz is step-heated in a lithium metaborate (LiBO2) flux and a high-purity O2 atmosphere, first at 500 C for 30 min and then at 1100 C for 3 h. Released carbon species are oxidised to form CO2 via secondary hot-quartz-bed oxidation, followed by cryogenic collection and purification. Sample yields are measured manometrically, and samples are diluted with 14C-free CO2. A small aliquot of CO2 is collected for δ13C analysis, and the remaining CO2 is graphitised using H2 reduction over an Fe catalyst. We measured 14C∕13C isotope ratios at either the Lawrence Livermore National Laboratory Center for Accelerator Mass Spectrometry (LLNL-CAMS) or Woods Hole National Ocean Sciences Accelerator Mass Spectrometry (NOSAMS) facility (Table S2). Stable carbon isotope ratios were measured at the UC Davis Stable Isotope Facility.

Apparent exposure ages were calculated using version 3 of the online calculators formerly known as the CRONUS-Earth online calculators (Balco et al., 2008). The online calculators use the production rate scaling method for neutrons, protons, and muons of Lifton et al. (2014) (also known as LSDn). We use repeat measurements of the in situ 14C concentration of the CRONUS-A interlaboratory comparison standard (Jull et al., 2015; Goehring et al., 2019) to calibrate the 14C production rate. We assume CRONUS-A is saturated with respect to in situ 14C, given that, based on geological mapping and an ash chronology, the sampling location has remained ice-free since > 11.3 Ma (Marchant et al., 1993). All reported in situ 14C measurements from CRONUS-A, made at multiple laboratories, yield concentrations equivalent to saturation based on other calibration data from elsewhere in the world (e.g. Jull et al., 2015; Fülöp et al., 2019; Goehring et al., 2019; Lamp et al., 2019). We use the CRONUS-A measurements to calibrate the 14C production rate to reduce scaling extrapolations. Repeat measurements of both CRONUS-A and other samples using the TU-CEGS show that the reproducibility of in situ 14C measurements is approximately 6 %. We therefore use a 6 % uncertainty for our measured in situ 14C concentrations when calculating exposure ages, as this exceeds the reported analytical uncertainty for all of our in situ 14C measurements. Ages are included in Table S2 for completeness but are primarily discussed in the text as either finite or infinite ages. Infinite ages are those for which the measured concentration is above the uncertainty of the saturation concentration for the elevation of a given sample.

We made seven replicate measurements from four samples from the Schmidt Hills that initially yielded saturation or near-saturation in situ 14C concentrations. We made the first four replicate measurements using the same samples to test the validity of the saturation or near-saturation initial measurements. The second set of measurements produced in situ 14C concentrations below saturation. Given the difference between the initial measurements and the replicates, we made a further three measurements from three of the same four samples.

3 Results

The vast majority of the 10Be ages reported by Hein et al. (2011, 2014) in the Shackleton Range exceed 100 ka, whilst we find finite 14C ages at both Mt Skidmore and Mt Provender (Figs. 5 and S6). At Mt Skidmore, finite ages are evident across the entire Mt Skidmore transect, including those sampled proximal to the Stratton Glacier (Fig. 5). Samples were collected from multiple ridges of Mt Skidmore and thus would not necessarily be expected to form a single age-elevation line. The uppermost sample proximal to the Slessor Glacier, collected ∼310 m above the modern ice surface, provides a lower limit for the LGM ice thickness of the ice mass. The two samples proximal to the Stratton Glacier, an erratic and bedrock sample with ∼17 m a.s.l. between them, are indistinguishable from one another within uncertainties and constrain the LGM thickness to at least 130 m thicker than present. At Mt Provender, one sample collected proximal to the local Slessor Glacier margin yields a finite age. A second sample from ∼890 m a.s.l. (∼655 m above the modern ice surface) yields an infinite age, placing an upper limit on the LGM thickness at ∼655 m larger than present. We note that the upper limit of 655 m is based on a single in situ 14C measurement and discuss this limitation further in Sect. 4.2. If quartz was available for additional samples previously collected from Mt Provender (Hein et al., 2011, 2014), then further measurements could have been made to validate this measurement. The quartz was, however, exhausted in the process of measuring long-lived nuclides. One sample from Mt Skidmore, collected from ∼284 m a.s.l., yields an infinite age. Above the saturated sample we observe seven finite-age samples which require significant periods of burial beneath ice to account for their in situ 14C concentrations. It is glaciologically impossible to have the sample at ∼284 m a.s.l. exposed for ca. 35 kyr whilst those above it were covered presumably by the Slessor and Stratton glaciers. The infinite age of the sample could be due to scatter within the 14C measurements, and the fact that the sample is an erratic does allow the possibility of an unlikely geomorphic scenario. As described in Sect. 1.2, erratic samples may be sourced from mass movement onto glacier surfaces, producing spuriously high 14C concentrations (Balco et al., 2019).

https://www.the-cryosphere.net/13/2935/2019/tc-13-2935-2019-f05

Figure 5(a) Elevation versus in situ 14C concentration of samples from the Shackleton Range. Circles are erratic cobbles; triangles are bedrock. Some error bars are smaller than their respective data points. Horizontal dashed lines show the approximate elevation of the modern ice surface at each site. Light grey lines indicate isochrons of exposure duration. Thick black line and grey shading are the saturation concentration and associated error envelope. (b) Exposure ages from this study (in situ 14C) and 10Be ages of Hein et al. (2011, 2014). Samples yielding infinite in situ 14C ages are not presented on the right-hand plot.

Download

On the Lassiter Coast, Johnson et al. (2019) report 10Be ages which, with the exception of three measurements, all exceed ∼100 ka, whilst all of the in situ 14C ages are finite and fall within the Holocene (Figs. 6 and S7). The associated in situ 14C concentrations are similar over the range of sample elevations (Fig. 6). The uppermost sample, collected ∼385 m above the modern ice surface, provides a lower limit on the thickness of LGM ice at the Lassiter Coast. The small range of ages across a ca. 300 m transect indicate that ice thinning occurred rapidly at this study site (Johnson et al., 2019).

https://www.the-cryosphere.net/13/2935/2019/tc-13-2935-2019-f06

Figure 6(a) Elevation versus in situ 14C concentration of samples collected from the Lassiter Coast. All samples are bedrock. (b) In situ 14C exposure ages with 10Be ages of Johnson et al. (2019).

Download

At the Schmidt Hills, 10Be ages from Balco et al. (2016) and Bentley et al. (2017) range from ∼140 ka to 3 Ma (Fig. 7). We observe finite ages at low elevations and finite, close to infinite, and infinite ages at higher elevations (Figs. 7 and S8). Given that higher elevations cannot be covered by ice unless lower elevations were also covered, we remeasured the apparently infinite- and near-infinite-age samples (∼500 to ∼920 m a.s.l., or ∼270 to ∼690 m above the modern ice surface). The replicate results (Fig. 7) show high variability, greater than that observed in previous repeat measurements of CRONUS-A and other samples (Goehring et al., 2019). There is no apparent analytical reason for the initial measurements yielding infinite or near-infinite ages and then yielding differing concentrations with repeat measurements. Samples yielding only finite ages (those that were not measured multiple times) are observed up to ∼420 m a.s.l., or ∼190 m above the modern ice surface. In addition, the bedrock sample collected from ca. 1035 m a.s.l. yields a finite age, indicative of a LGM thickness at least ∼800 m larger than present for the FIS at the Schmidt Hills. The agreement between the bedrock age and the finite measurements from lower elevations means we conclude that, at the Schmidt Hills, the FIS was 800 m thicker than present at the LGM. This conclusion and the repeat measurements with a high degree of scatter are discussed in Sect. 4.1 and 4.2.

https://www.the-cryosphere.net/13/2935/2019/tc-13-2935-2019-f07

Figure 7(a) Elevation versus in situ 14C concentration of samples collected from the Schmidt Hills. All samples are erratics with the exception of the highest elevation sample, shown with a triangle. Samples with replicate measurements are displayed with differing symbols. (b) Schmidt Hills exposure ages from this study (in situ 14C) and those of Balco et al. (2016) and Bentley et al. (2017) (10Be). Measurements yielding infinite in situ 14C ages are not presented on the right-hand plot.

Download

The two samples collected from the Thomas Hills yield finite ages within ∼0.2 ka of one another (Figs. 8 and S9). Results thus indicate that the FIS was at least ∼320 m thicker than present at the LGM at the Thomas Hills. The apparent in situ 14C ages, at ∼10 ka, are consistent with a cluster of 10Be ages between 7 and 9 ka in the Thomas Hills reported by Balco et al. (2016) from 225 m above the modern FIS surface, as well as a 10Be age of 4.2 ka reported by Bentley et al. (2017) collected 125 m above the modern ice surface. Considering the evidence for significant LGM thickening of the FIS from our in situ 14C results from the Thomas Hills, as well as 10Be ages of Balco et al. (2016) and Bentley et al. (2017) from both the Williams and Thomas Hills, we infer that it is likely that the FIS reached up to 800 m above its present thickness at the LGM at the Schmidt Hills. We discuss this inference further in the following section.

https://www.the-cryosphere.net/13/2935/2019/tc-13-2935-2019-f08

Figure 8(a) Elevation versus in situ 14C concentration of samples collected from the Thomas Hills. All samples are erratics. Note that both plots contain in situ 14C data for two samples within close agreement, such that the points overlap. (b) Thomas Hills exposure ages from this study (in situ 14C) and those of Balco et al. (2016) and Bentley et al. (2017) (10Be).

Download

4 Discussion

4.1 Assessment of 14C elevation transects

The premise of our study is that one can clearly infer if a site was ice-covered at the LGM by determining whether the in situ 14C concentration of samples from that site are at or below saturation. In this section we assess the success of the approach. To assess the validity of this method, we can, for example, identify where the in situ 14C data record ice thinning, with saturated samples or the oldest exposure ages at the highest elevations and a trend of decreasing in situ 14C age toward modern ice surfaces. Consistency between in situ 14C data and other nuclide concentrations (e.g. 10Be) could also help validate the in situ 14C measurements. We also look at factors beyond the in situ 14C concentrations, such as the glaciological link between study sites, which may add clarity where the in situ 14C measurements show a high degree of scatter.

At some sites our results are consistent with the premise, as well as internally consistent. At the Lassiter Coast, ages decrease toward the present ice sheet surface. Though limited by the number of samples, two measurements from Mt Provender align with the premise of our study, in that we find a finite age located at a low elevation with an infinite age above it. In the Thomas Hills we see consistency between the finite 14C ages and previously published 10Be ages (Balco et al., 2016; Bentley et al., 2017). Fogwill et al. (2014) also observe consistency between 14C and 10Be ages, which constrain the LGM thickness and dynamics of the Rutford Ice Stream. However, we observe apparently finite ages above apparently infinite ages at the Schmidt Hills, a scenario that is glaciologically impossible if our assumptions are correct that samples are indeed glacial erratics that have either been deposited previously and repeatedly covered by cold-based ice or delivered to their sampling location during the last glaciation and were sourced subglacially. The scatter observed in the repeat measurements (Fig. 7) is greater than that of repeat measurements made of CRONUS-A and other samples made in our laboratory (Goehring et al., 2019). Three samples from the Schmidt Hills (006-COU, 008-NNS, and 046-NNS, collected from ∼920, ∼710, and ∼500 m a.s.l., respectively) were previously measured for their in situ 14C content and were published by Balco et al. (2016). All three of the samples previously measured by Balco et al. (2016) yielded higher concentrations (two of which were above saturation with the third at saturation) than their new measurements presented in this study. Furthermore, two of the three samples (006-COU and 046-NNS) were measured multiple times (in this study) and display the high scatter under discussion. Balco et al. (2016) proposed unrecognised measurement error as the cause of the spuriously high in situ 14C concentrations. Why the replicate measurements from samples from the Schmidt Hills display a high degree of scatter remains to be determined.

The most likely reason for 14C measurement error is contamination by modern 14C, which would result in a spuriously high concentration. In contrast, a spuriously low concentration is less likely, and we are not aware of any documented instances of this. In our laboratory we have found that it is relatively easy to contaminate a sample with modern carbon through the use of organic compounds in the froth flotation mineral separation procedure (Nichols and Goehring, 2019). However, froth flotation was not used to isolate the quartz of any of the samples for which replicate measurements were made. On multiple occasions we have observed spuriously high 14C concentrations, far in excess of saturation concentrations, from quartz separates of fine grain sizes (ca. 60 µm) that were not isolated using froth flotation. We do not yet know the reason for the fine grain sizes yielding elevated 14C concentrations, but one hypothesis is that the finite-age replicate measurements were unintentionally made using quartz separates with a coarser average grain size than the initial infinite measurements. We believe the above observations indicate that the increased scatter may be the result of measurement difficulties, perhaps lithology- or grain size-specific.

Regardless of the cause of the high degree of scatter observed in the replicate measurements, we need to discuss possible explanations for apparently infinite ages at lower elevations than apparently finite ages to isolate which measurements (infinite vs. finite replicates) are the most valid to base interpretations on. At the Schmidt Hills, the hypothesis that infinite ages situated below finite ages are spurious and due to measurement errors is consistent with the glaciological relationship amongst the Schmidt, Thomas, and Williams Hills (see Sect. 4.2) and is also consistent with the finite bedrock age sourced from a higher elevation. The bedrock age is a robust constraint because the sample cannot have been subjected to geomorphic scenarios that could cause the resulting age to misrepresent the timing of deglaciation. The hypothesis that the infinite ages are correct produces a steep LGM surface slope and is not consistent with thickness estimates from the Williams and Thomas Hills. We elaborate on this point in Sect. 4.2.

As described in Sect. 1.2, it is theoretically possible for in situ 14C saturated erratic samples to occur at lower elevations than finite ages in rare situations if the former were transported by LGM ice. Balco et al. (2019) observed an apparently saturated sample beneath finite-age samples. Supported by field observations, Balco et al. (2019) propose that the saturated sample was sourced from a rockfall upstream and transported to the study site as supraglacial debris, explaining the elevated in situ 14C concentration. Whilst this could explain the low-elevation saturated sample at Mt Skidmore, as well as infinite measurements situated beneath finite measurements at the Schmidt Hills, it does not explain the poor reproducibility of the Schmidt Hills measurements.

We conclude that the basic concept works, as shown at the Lassiter Coast and the Shackleton Range, as well as in other aforementioned studies. In the following section we discuss the implications for LGM ice sheet reconstructions. However, it is clear that more investigation into laboratory issues and geological and geomorphic factors is required to identify the cause or causes of apparently site- or lithology-specific excess scatter in situ 14C measurements.

https://www.the-cryosphere.net/13/2935/2019/tc-13-2935-2019-f09

Figure 9Terrestrial ice thickness constraints inferred from measurements of cosmogenic nuclides around the WSE. Constraints for the Lassiter Coast (LC), Shackleton Range, and the Schmidt Hills are sourced from this study. All other ice thickness values and locations are the same as in Fig. 2.

4.2 LGM ice thicknesses in the Weddell Sea Embayment

Our LGM thickness estimates are summarised in Fig. 9. The new in situ 14C concentrations indicate that the vast majority, if not the entirety, of Mt Skidmore and presumably much of Mt Provender were covered by ice at the LGM (Figs. 5 and 10). The highest elevation samples on Mt Skidmore proximal to the Slessor Glacier yield infinite ages and indicate that the ice stream was at least 300 m thicker at the LGM than at present. This assumes the samples were not influenced by expansion of local ice masses from the southeast (Fig. S1). If so, and assuming the surface gradient of Slessor Glacier during the LGM was similar to today, this would suggest the Slessor Glacier was ∼300 m thicker at Mt Provender at the LGM. With no high-elevation infinite ages found on Mt Skidmore, our thickness estimates for the Slessor Glacier are likely conservative estimates. Finite ages are observed across the entire Mt Skidmore transect, and there is only a single exposed peak (between Ice Tongue A and Ice Tongue B of the Stratton Glacier, Fig. S1) that is at a higher elevation than our sampling locations (ca. 25 m higher). Presumably, given the evidence for the expansion of the Slessor and Stratton glaciers, this small peak was covered by these or local ice masses at the LGM. Our data therefore indicate that, regardless of the source, the Mt Skidmore site was covered by ice during the LGM, whilst the top of Mt Provender remained exposed. Whilst the upper limit of LGM ice at Mt Provender is based on a single sample, we believe this sample is a reliable indicator of LGM ice thickness for the following reasons. The sample is sourced from bedrock and therefore cannot have been subjected to geomorphic scenarios causing the exposure age to misrepresent the timing of ice retreat. Furthermore, froth flotation, which introduces modern carbon to sample material (Nichols and Goehring, 2019), was not used to isolate quartz for this sample. Our thickness constraints (∼300–655 m) supersede those of previous exposure dating studies that found no evidence from long-lived isotopes for a thicker Slessor Glacier at the LGM (Hein et al, 2011, 2014). Our LGM thickness constraints for the Slessor Glacier are consistent with our other sites as well as those of previous authors for a significantly thicker FIS at the LGM (Balco et al., 2016; Bentley et al., 2017).

https://www.the-cryosphere.net/13/2935/2019/tc-13-2935-2019-f10

Figure 10(a) Exposure-age results projected onto an elevation profile along flow line of the Slessor Glacier. (b) Flow line location is shown in the map. Infinite 14C measurements are offset in regard to their distance along flow line to improve readability. The 10Be data included are those from Hein et al. (2011, 2014) which yield exposure ages below 12 ka (LSDn scaling, http://antarctica.ice-d.org, last access: 2 November 2019). Elevation data for ice surfaces and map shading are sourced from the Reference Elevation Model of Antarctica (REMA; Howat et al., 2018). Grounding line positions sourced from the MEaSUREs program V2 (Rignot et al., 2011, 2014, 2016). Minimum LGM surface is the modern-day surface profile with the elevation increased above present using our minimum LGM thickness estimates.

The new in situ 14C results from the Lassiter Coast show that bedrock surfaces 385 m above the modern ice surface were covered by ice at the LGM. As with results in the Pensacola Mountains, with only a lower limit for the LGM thickness of 360 m, there could have been thicker ice on the Lassiter Coast at the LGM. The in situ 14C measurements contrast with 10Be measurements that were likely influenced by cold-based ice cover, resulting in nuclide inheritance (Johnson et al., 2019). The finite in situ 14C ages of samples collected from 628 to 875 m a.s.l., with ages between 6.0±0.7 ka and 7.5±0.9 ka, are consistent with a minimum age of grounded ice retreat from a marine sediment core close to the modern ice shelf edge of 5.3±0.3 cal kyr BP (Hedges et al., 1995; Crawford et al., 1996; Fig. 1). The fact that significant thinning occurred in the Holocene may help explain the misfit between GIA models and GPS measurements in Palmer Land (Wolstencroft et al., 2015). A thicker ice load at the LGM than that used by current ice models or present ice load estimates that persist into the Holocene are two potential solutions postulated by Wolstencroft et al. (2015) to explain the misfit. Further work is needed to take our new ice history into account and to investigate if a minimum of 385 m of ice at the LGM and subsequent rapid thinning at ∼7 ka at the Lassiter Coast can help account for the offset.

Our in situ 14C data indicate that the FIS was at least 800 m thicker than present at the Schmidt Hills at the LGM, which contrasts with previous studies which found no evidence for the LGM thickness of the FIS at the Schmidt Hills (Balco et al., 2016; Bentley et al., 2017). We base our LGM thickness estimate on the aforementioned finite-age repeat measurements and the finite-age bedrock sample rather than on the poorly reproduced infinite-age measurements. There is robust evidence for a FIS that was at least 500 m thicker than present at the LGM at the Williams Hills, located only 50 km upstream of the Schmidt Hills (Figs. 1 and 11; Balco et al., 2016; Bentley et al., 2017). Given the evidence for a significantly thicker FIS proximal to the Schmidt Hills, we argue that the repeat measurements and the bedrock measurement indicative of the FIS being 800 m thicker are glaciologically most likely and thus base our LGM ice thickness estimates on them. Using the infinite measurements and accompanying constraint at the Schmidt Hills for the LGM thickness of 320 m thicker than present produces a steep surface slope from the nearby Williams Hills (Fig. 11), though less so than the surface slope produced when no LGM thickening is inferred at the Schmidt Hills based on 10Be measurements (Balco et al., 2016). The two measurements from the Thomas Hills provide a lower limit for the LGM thickness, but the possibility remains that there was more thickening than the ca. 320 m in situ 14C constraint. Fig. 11 tentatively indicates that the FIS may have been ∼900 m thicker when using the modern surface profile of the FIS increased in elevation up to the height of the finite ages from the Schmidt Hills and post-12 ka 10Be ages from Balco et al. (2016) and Bentley et al. (2017) from the Williams Hills. This is a tentative interpretation because, if thickening is sea level controlled, there would be progressively less thinning expected upstream.

https://www.the-cryosphere.net/13/2935/2019/tc-13-2935-2019-f11

Figure 11(a) Exposure-age results projected onto an elevation profile along flow line of the Foundation Ice Stream. (b) Flow line location is shown in the map. Infinite 14C measurements are offset in regard to their distance along flow line to improve readability. The 10Be data included are those from Balco et al. (2016) and Bentley et al. (2017) which yield exposure ages below 12 ka (LSDn scaling, antarctica.ice-d.org). Elevation data for ice surfaces and map shading are sourced from the Reference Elevation Model of Antarctica (REMA; Howat et al., 2019). Local ice margins are highly simplified. Grounding line positions are sourced from the MEaSUREs program V2 (Rignot et al., 2011, 2014, 2016). Minimum LGM surface is the modern-day surface profile with the elevation increased above present using our minimum LGM thickness estimates.

Our LGM ice thickness constraints are consistent with evidence for significantly thicker ice at the LGM in the Ellsworth Mountains (Hein et al., 2016, Fig. 2) and also likely consistent with measurements in Bentley et al. (2006). The post-LGM exposure ages of Hein et al. (2016) constrain LGM thicknesses to between 475, 373, and 247 m larger than present at three study sites in the Ellsworth Mountains. A pulse of up to 410 m of thinning appears similar both in scale and timing to the rapid ice surface lowering of 385 m recorded at the Lassiter Coast. Furthermore, measurements of long-lived nuclides by Bentley et al. (2006) show that there has been at least 300 m of thinning since the LGM in the Behrendt Mountains.

4.3 Grounding line position and flow line modelling comparison

Whitehouse et al. (2017) use their flow line model to reproduce the modern FIS ice surface profile and investigate the response of the ice stream to the onset of glacial and interglacial conditions. The following results from Whitehouse et al. (2017) are from their experiments in which the FIS is routed to the east of Berkner Island, which is believed to have been done during the LGM based on modelling studies (Le Brocq et al., 2011; Whitehouse et al., 2012) and aforementioned marine geological evidence for the former presence of grounded ice (Sect. 1.1). Under glacial conditions the FIS thickens by ∼300 to ∼500 m adjacent to the Thomas Hills, ∼200 to ∼400 m adjacent to the Williams Hills, ∼150 to ∼350 m adjacent to the Schmidt Hills, and ∼100 to ∼300 m proximal to the Shackleton Range. The lower value for each location is sourced from flow line experiments during which the grounding line of the FIS reaches a stable position at the northern margin of Berkner Island, with the higher value sourced from a scenario during which the grounded ice stream stabilises at the shelf edge. Our in situ 14C LGM thickness constraints at each study location in the Pensacola Mountains and Shackleton Range exceed the upper estimates of the FIS flow line model of Whitehouse et al. (2017) under glacial conditions. The flow line model shows that the FIS, a major contributor to the total WSE ice flux, is able to reach a stable position at the shelf edge when tuned using LGM thickness constraints lower than those presented here. Therefore, our thickness estimates add strength to the hypothesis that grounded ice occupying the WSE during the LGM reached a stable position located at the shelf edge (Bentley and Anderson, 1998; Hillenbrand et al., 2014).

4.4 Sea level contribution

To estimate the contribution to post-LGM sea level rise of the WSE we use a highly simplified scenario in which a range of minimum LGM thickness change estimates are distributed evenly across the WSE using an area for the sector defined by Hillenbrand et al. (2014). Distributing the lowest of our minimum LGM thickness constraints, 310 m for the Slessor Glacier, over the entire WSE produces a minimum sea level equivalent (SLE) of 2.2 m. When using the highest of our minimum thickness estimates, 800 m for the FIS, the minimum SLE increases to 5.8 m. Using the average minimum LGM thickness constraint for our three study sites (580 m) produces a minimum SLE for the sector of 4.2 m. This scenario lacks any glaciological basis and is unrealistic, with no variation in ice thickness with location and no consideration of ice dynamics, isostasy, or bathymetry. Hence, further work is required to produce a realistic SLE for the WSE using our in situ 14C thickness constraints.

We compare our in situ 14C LGM thickness estimates with the predicted LGM thickness change of three published ice sheet models at each of our study sites to evaluate our minimum SLE estimates. We also quantify the WSE-sourced SLE for each model output. By comparing our data with the predicted LGM thickness from the model outputs, we can see which models predict LGM thickness changes in excess of and below our in situ 14C thickness constraints. We compare our data with the predicted LGM thickness change at each of our sites from the ice sheet modelling of Le Brocq et al. (2011), Whitehouse et al. (2012), and Golledge et al. (2014), which predict a SLE for the WSE of ca. 3.0, 1.5, and 4.6 m, respectively.

From Fig. 12 it is apparent that the model output of Golledge et al. (2014) exceeds the thickness constraints at each of our sites. With a SLE of ca. 4.6 m for the WSE, this places a more robust upper limit on the minimum SLE contribution of the WSE using our data, showing that our upper minimum SLE estimate of 5.8 m is likely an overestimation due to the limitations outlined above. The only site where our minimum LGM ice thickness constraint exceeds any of the predicted LGM thickness changes from the model outputs is at the Lassiter Coast, where a LGM thickness of 385 m larger than present exceeds the model output-based thickness estimate of both Le Brocq et al. (2011) and Whitehouse et al. (2012). The Lassiter Coast data indicate that the lower limit for the SLE for the WSE is between 3.0 and 4.6 m, whilst evidence from all other sites suggests it was < 1.5 m.

https://www.the-cryosphere.net/13/2935/2019/tc-13-2935-2019-f12

Figure 12Predicted LGM ice thickness change from three ice sheet model outputs at each of our study sites and their associated sea level equivalent (SLE) for the Weddell Sea Embayment (WSE) sector of the Antarctic ice sheet. (a) is Mt Skidmore, (b) Mt Provender, (c) the Lassiter Coast, (d) Schmidt Hills, and (e) Thomas Hills. Vertical blue lines show the interpreted LGM thickness change at each site based on our in situ 14C data. For (a) and (b), the two vertical blue lines show the range of thickness estimates for the two sites, with the upper limit constrained by the highest elevation saturated sample at Mt Provender. “G2014” refers to Golledge et al. (2014), “LB2011” refers to Le Brocq et al. (2011), and “W2012” refers to Whitehouse et al. (2012). Errors are not provided for the model outputs. The average error of published SLEs associated with model outputs for the entire ice sheet is 1.45 m (see Simms et al., 2019). We therefore use an error of 0.3 m for the three model SLEs, which is 22 % of the average error (22 % is the proportion of the AIS that the WSE drains; see Joughin et al., 2006).

Download

Based on the above, we conclude that our minimum LGM thickness constraints indicate that the WSE contributed < 4.6 m, and possibly as little as < 1.5 m, toward postglacial sea level rise. This is a range of minimum contributions to sea level rise and not a minimum–maximum range, as the values are informed using only minimum thickness constraints. Because this is an estimate for the lower limit of the SLE for the WSE, we cannot rule out a larger contribution.

A SLE value of < 4.6 m places our estimate between those modelled by Bentley et al. (2010) (1.4 to 2 m) and Bassett et al. (2007) (13.1 to 14.1 m). Using the estimate based on all sites with the exception of the Lassiter Coast data, the minimum SLE estimate of < 1.5 m is consistent with the lower end of published SLEs for the sector. Our exposure ages indicate the Weddell Sea sector contributed to sea level during the early to mid-Holocene, though they do not preclude a significant contribution earlier than this. Our estimates imply that the sector provided a modest contribution to global sea level. Whitehouse et al. (2017) estimate the sea level contribution of the FIS to between ∼0.05 and ∼0.13 m. Given that our 14C thickness constraints for the FIS, including those in the Shackleton Range, exceed all of those used by Whitehouse et al. (2017) to tune their flow line model, we propose that the sea level contribution for the FIS was greater than their upper estimate of ∼0.13 m.

5 Conclusions

We present LGM ice thickness constraints for three locations within the WSE of Antarctica. In situ 14C measurements constrain the LGM thickness of the Foundation Ice Stream to at least ca. 800 m thicker than present in the Schmidt Hills and at least 320 m thicker than present in the Thomas Hills, both in the Pensacola Mountains. The Slessor Glacier was at least 310 m and up to 655 m thicker than present at the LGM. Finally, LGM ice was at least 385 m thicker than present at the Lassiter Coast. Our thickness constraints resolve a significant disconnect between previous terrestrial evidence for minimal LGM thickening in some locations from long-lived nuclides and marine evidence for a significantly laterally expanded ice sheet with the grounding line located at the shelf edge. Our in situ 14C measurements made from samples at the Schmidt Hills exhibit higher-than-expected scatter in replicate measurements. Identifying the source of excess scatter will take further work. In terms of the contribution of the ice sheet sector to global sea level rise since the LGM, we estimate, primarily based on minimum estimates which do not constrain the upper limit of ice thickness changes, that the WSE contributed < 4.6 m, and possibly < 1.5 m.

Data availability

All sample data, including photographs when available, are available in the Informal Cosmogenic-Nuclide Exposure-Age Database (ICE-D) (http://antarctica.ice-d.org/pub/120; Balco, 2019).

Supplement

The supplement related to this article is available online at: https://doi.org/10.5194/tc-13-2935-2019-supplement.

Author contributions

The study was conceived by BMG and GB. Sample material was originally collected and provided by ASH, JSJ, GB, and CT. Sample preparation and analysis undertaken by KAN and BMG. The manuscript was written by KAN with help from all authors.

Competing interests

The authors declare that they have no conflict of interest.

Acknowledgements

Keir A. Nichols and Brent M. Goehring acknowledge support from NSF-OPP grant 1542936. Greg Balco acknowledges support from NSF-OPP 1542976 and from the Ann and Gordon Getty Foundation. Geospatial support for this work was provided by the Polar Geospatial Center under NSF-OPP awards 1043681 and 1559691. This work forms part of the British Antarctic Survey Polar Science for Planet Earth programme, funded by the Natural Environment Research Council.

Financial support

This research has been supported by the National Science Foundation (grant nos. 1542936 and 1542976).

Review statement

This paper was edited by Arjen Stroeven and reviewed by Jane L. Andersen and one anonymous referee.

References

Ackert, R. P., Mukhopadhyay, S., Parizek, B. R., and Borns, H. W.: Ice elevation near the West Antarctic Ice Sheet divide during the Last Glaciation, Geophys. Res. Lett., 34, 1–6, https://doi.org/10.1029/2007GL031412, 2007. 

Arndt, J. E., Schenke, H. W., Jakobsson, M., Nitsche, F. O., Buys, G., Goleby, B., Rebesco, M., Bohoyo, F., Hong, J., Black, J., Greku, R., Udintsev, G., Barrios, F., Reynoso-Peralta, W., Taisei, M., and Wigley, R.: The international bathymetric chart of the Southern Ocean (IBCSO) version 1.0-A new bathymetric compilation covering circum-Antarctic waters, Geophys. Res. Lett., 40, 3111–3117, https://doi.org/10.1002/grl.50413, 2013. 

Arndt, J. E., Hillenbrand, C.-D., Grobe, H., Kuhn, G., and Wacker, L.: Evidence for a dynamic grounding line in outer Filchner Trough, Antarctica, until the early Holocene, Geology, 45, 1035–1038, https://doi.org/10.1130/G39398.1, 2017. 

Balco, G.: Informal cosmogenic-nuclide exposure-age database (ICE-D) Antarctica, available at: http://antarctica.ice-d.org/pub/120, 2019. 

Balco, G., Stone, J. O., Lifton, N. A., and 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, 174–195, https://doi.org/10.1016/j.quageo.2007.12.001, 2008. 

Balco, G., Todd, C., Huybers, K., Campbell, S., Vermeulen, M., Hegland, M., Goehring, B. M., and Hillebrand, T. R.: Cosmogenic-nuclide exposure ages from the Pensacola Mountains adjacent to the foundation ice stream, Antarctica, Am. J. Sci., 316, 542–577, https://doi.org/10.2475/06.2016.02, 2016. 

Balco, G., Todd, C., Goehring, B. M., Moening-Swanson, I., and Nichols, K.: Glacial geology and cosmogenic-nuclide exposure ages from the Tucker Glacier – Whitehall Glacier confluence, northern Victoria Land, Antarctica, Am. J. Sci., 319, 255–286, https://doi.org/10.2475/04.2019.01, 2019. 

Bassett, S. E., Milne, G. A., Bentley, M. J., and Huybrechts, P.: Modelling Antarctic sea-level data to explore the possibility of a dominant Antarctic contribution to meltwater pulse IA, Quaternary Sci. Rev., 26, 2113–2127, https://doi.org/10.1016/j.quascirev.2007.06.011, 2007. 

Bentley, M. J. and Anderson, J. B.: Glacial and marine geological evidence for the ice sheet configuration in the Weddell Sea–Antarctic Peninsula region during the Last Glacial Maximum, Antarct. Sci., 10, 309–325, https://doi.org/10.1017/s0954102098000388, 1998. 

Bentley, M. J., Fogwill, C. J., Kubik, P. W., and Sugden, D. E.: Geomorphological evidence and cosmogenic 10Be∕26Al exposure ages for the Last Glacial Maximum and deglaciation of the Antarctic Peninsula Ice Sheet, Geol. Soc. Am. Bull., 118, 1149–1159, https://doi.org/10.1130/B25735.1, 2006. 

Bentley, M. J., Fogwill, C. J., Brocq, A. M. Le, Hubbard, A. L., Sugden, D. E., Dunai, T. J., and Freeman, S. P. H. T.: Deglacial history of the West Antarctic Ice Sheet in the Weddell Sea embayment: Constraints on past ice volume change, Geology, 38, 411–414, https://doi.org/10.1130/G30754.1, 2010. 

Bentley, M. J., Hein, A. S., Sugden, D. E., Whitehouse, P. L., Shanks, R., Xu, S., and Freeman, S. P. H. T.: Deglacial history of the Pensacola Mountains, Antarctica from glacial geomorphology and cosmogenic nuclide surface exposure dating, Quaternary Sci. Rev., 158, 58–76, https://doi.org/10.1016/j.quascirev.2016.09.028, 2017. 

Briggs, R. D., Pollard, D., and Tarasov, L.: A data-constrained large ensemble analysis of Antarctic evolution since the Eemian, Quaternary Sci. Rev., 103, 91–115, https://doi.org/10.1016/j.quascirev.2014.09.003, 2014. 

Chmeleff, J., Von Blanckenburg, F., Kossert, K., and Jakob, D.: Determination of the 10Be half-life by multicollector ICP-MS and liquid scintillation counting, Nucl. Instrum. Meth.-B, 268, 192–199, https://doi.org/10.1016/j.nimb.2009.09.012, 2010. 

Crawford, K., Kuhn, G., and Hambrey, M.: Changes in the character of glaciomarine sedimentation in the southwestern Weddell Sea, Antarctica: evidence from the core PS1423-2, Ann. Glaciol., 22, 200–204, 10.3189/1996AoG22-1-200-204, 1996. 

Fogwill, C. J., Bentley, M. J., Sugden, D. E., Kerr, A. R., and Kubik, P. W.: 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, https://doi.org/10.1130/G19795.1, 2004. 

Fogwill, C. J., Turney, C. S. M., Golledge, N. R., Rood, D. H., Hippe, K., Wacker, L., Wieler, R., Rainsley, E. B., and Jones, R. S.: Drivers of abrupt Holocene shifts in West Antarctic ice stream direction determined from combined ice sheet modelling and geologic signatures, Antarct. Sci., 26, 674–686, https://doi.org/10.1017/S0954102014000613, 2014. 

Fülöp, R. H., Fink, D., Yang, B., Codilean, A. T., Smith, A., Wacker, L., Levchenko, V., and Dunai, T. J.: The ANSTO – University of Wollongong in-situ 14C extraction laboratory, Nucl. Instrum. Meth.-B, 438, 204–213, https://doi.org/10.1016/j.nimb.2018.04.018, 2019. 

Gjermundsen, E. F., Briner, J. P., Akçar, N., Foros, J., Kubik, P. W., Salvigsen, O., and Hormes, A.: Minimal erosion of Arctic alpine topography during late Quaternary glaciation, Nat. Geosci., 8, 789–792, https://doi.org/10.1038/ngeo2524, 2015. 

Goehring, B. M., Wilson, J., and Nichols, K.: A fully automated system for the extraction of in situ cosmogenic carbon-14 in the Tulane University cosmogenic nuclide laboratory, Nucl. Instrum. Meth.-B, 455, 284–292, https://doi.org/10.1016/j.nimb.2019.02.006, 2019. 

Golledge, N. R., Menviel, L., Carter, L., Fogwill, C. J., England, M. H., Cortese, G., and Levy, R. H.: Antarctic contribution to meltwater pulse 1A from reduced Southern Ocean overturning, Nat. Commun., 5, 1–10, https://doi.org/10.1038/ncomms6107, 2014. 

Hedges, R. E. M., Housley, R. A., Ramsey, C. B., and Van Klinken, G. J.: Radiocarbon Dates from the Oxford AMS System: Archaeometry Datelist 19, Archaeometry, 37, 195–214, 1995. 

Hein, A. S., Fogwill, C. J., Sugden, D. E., and Xu, S.: Glacial/interglacial ice-stream stability in the Weddell Sea embayment, Antarctica, Earth Planet. Sc. Lett., 307, 211–221, https://doi.org/10.1016/j.epsl.2011.04.037, 2011. 

Hein, A. S., Fogwill, C. J., Sugden, D. E., and Xu, S.: Geological scatter of cosmogenic-nuclide exposure ages in the Shackleton Range, Antarctica: Implications for glacial history, Quat. Geochronol., 19, 52–66, https://doi.org/10.1016/j.quageo.2013.03.008, 2014. 

Hein, A. S., Marrero, S. M., Woodward, J., Dunning, S. A., Winter, K., Westoby, M. J., Freeman, S. P. H. T., Shanks, R. P., and Sugden, D. E.: Mid-Holocene pulse of thinning in the Weddell Sea sector of the West Antarctic ice sheet, Nat. Commun., 7, 1–8, https://doi.org/10.1038/ncomms12511, 2016. 

Herber, L. J.: Separation of feldspar from quartz by flotation, Am. Mineral., 54, 1212–1215, https://doi.org/10.4144/rpsj1954.25.192, 1969. 

Hillenbrand, C. D., Melles, M., Kuhn, G., and Larter, R. D.: Marine geological constraints for the grounding-line position of the Antarctic Ice Sheet on the southern Weddell Sea shelf at the Last Glacial Maximum, Quaternary Sci. Rev., 32, 25–47, https://doi.org/10.1016/j.quascirev.2011.11.017, 2012. 

Hillenbrand, C. D., Bentley, M. J., Stolldorf, T. D., Hein, A. S., Kuhn, G., Graham, A. G. C., Fogwill, C. J., Kristoffersen, Y., Smith, J. A., Anderson, J. B., Larter, R. D., Melles, M., Hodgson, D. A., Mulvaney, R., and Sugden, D. E.: Reconstruction of changes in the Weddell Sea sector of the Antarctic Ice Sheet since the Last Glacial Maximum, Quaternary Sci. Rev., 100, 111–136, https://doi.org/10.1016/j.quascirev.2013.07.020, 2014. 

Hodgson, D. A., Bentley, M. J., Schnabel, C., Cziferszky, A., Fretwell, P., Convey, P., and Xu, S.: Glacial geomorphology and cosmogenic 10Be and 26Al exposure ages in the northern Dufek Massif, Weddell Sea embayment, Antarctica, Antarct. Sci., 24, 377–394, https://doi.org/10.1017/S0954102012000016, 2012. 

Howat, I. M., Porter, C., Smith, B. E., Noh, M.-J., and Morin, P.: The Reference Elevation Model of Antarctica, The Cryosphere, 13, 665–674, https://doi.org/10.5194/tc-13-665-2019, 2019. 

Johnson, J. S., Nichols, K. A., Goehring, B. M., Balco, G., and Schaefer, J. M.: Abrupt mid-Holocene ice loss in the western Weddell Sea Embayment of Antarctica, Earth Planet. Sc. Lett., 518, 127–135, https://doi.org/10.1016/j.epsl.2019.05.002, 2019. 

Joughin, I., Bamber, J. L., Scambos, T., Tulaczyk, S., Fahnestock, M., and MacAyeal, D. R.: Integrating satellite observations with modelling: Basal shear stress of the Filcher-Ronne ice streams, Antarctica, Philos. T. R. Soc. A, 364, 1795–1814, https://doi.org/10.1098/rsta.2006.1799, 2006. 

Jull, A. J. T., Scott, E. M., and Bierman, P.: The CRONUS-Earth inter-comparison for cosmogenic isotope analysis, Quat. Geochronol., 26, 3–10, https://doi.org/10.1016/j.quageo.2013.09.003, 2015. 

Korschinek, G., Bergmaier, A., Faestermann, T., Gerstmann, U. C., Knie, K., Rugel, G., Wallner, A., Dillmann, I., Dollinger, G., Lierse von Gostomski, C., Kossert, K., Maiti, M., Poutivtsev, M., and Remmert, A.: A new value for the half-life of 10Be by Heavy-Ion Elastic Recoil Detection and liquid scintillation counting, Nucl. Instrum. Meth.-B, 268, 187–191, https://doi.org/10.1016/j.nimb.2009.09.020, 2010. 

Lamp, J. L., Young, N. E., Koffman, T., Schimmelpfennig, I., Tuna, T., Bard, E., and Schaefer, J. M.: Update on the cosmogenic in situ 14C laboratory at the Lamont-Doherty Earth Observatory, Nucl. Instrum. Meth.-B, 465, 157–162, https://doi.org/10.1016/j.nimb.2019.05.064, 2019. 

Larter, R. D., Graham, A. G. C., Hillenbrand, C.-D., Smith, J. A., and Gales, J. A.: Late Quaternary grounded ice extent in the Filchner Trough, Weddell Sea, Antarctica: New marine geophysical evidence, Quaternary Sci. Rev., 53, 111–122, https://doi.org/10.1016/j.quascirev.2012.08.006, 2012. 

Le Brocq, A. M., Bentley, M. J., Hubbard, A., Fogwill, C. J., Sugden, D. E., and Whitehouse, P. L.: Reconstructing the Last Glacial Maximum ice sheet in the Weddell Sea embayment, Antarctica, using numerical modelling constrained by field evidence, Quaternary Sci. Rev., 30, 2422–2432, https://doi.org/10.1016/j.quascirev.2011.05.009, 2011. 

Lifton, N., Sato, T., and Dunai, T. J.: Scaling in situ cosmogenic nuclide production rates using analytical approximations to atmospheric cosmic-ray fluxes, Earth Planet. Sc. Lett., 386, 149–160, https://doi.org/10.1016/j.epsl.2013.10.052, 2014. 

Marchant, D. R., Denton, G. H., and Swisher, C. C.: Miocene-Pliocene-Pleistocene Glacial History of Arena Valley, Quartermain Mountains, Antarctica, Geogr. Ann. Ser. A, 75, 269–301, https://doi.org/10.2307/521204, 1993. 

Nichols, K. A. and Goehring, B. M.: Isolation of quartz for cosmogenic in situ 14C analysis, Geochronology, 1, 43–52, https://doi.org/10.5194/gchron-1-43-2019, 2019. 

Norris, T. L., Gancarz, A. J., Rokop, D. J., and Thomas, K. W.: Half-Life of 26Al, J. Geophys. Res., 88, B331–B333, 1983. 

Peltier, W. R. and Fairbanks, R. G.: Global glacial ice volume and Last Glacial Maximum duration from an extended Barbados sea level record, Quaternary Sci. Rev., 25, 3322–3337, https://doi.org/10.1016/j.quascirev.2006.04.010, 2006. 

Pollard, D., Chang, W., Haran, M., Applegate, P., and DeConto, R.: Large ensemble modeling of the last deglacial retreat of the West Antarctic Ice Sheet: comparison of simple and advanced statistical techniques, Geosci. Model Dev., 9, 1697–1723, https://doi.org/10.5194/gmd-9-1697-2016, 2016. 

Pollard, D., Gomez, N., and Deconto, R. M.: Variations of the Antarctic Ice Sheet in a Coupled Ice Sheet-Earth-Sea Level Model: Sensitivity to Viscoelastic Earth Properties, J. Geophys. Res.-Earth, 122, 2124–2138, https://doi.org/10.1002/2017JF004371, 2017. 

Rignot, E., Mouginot, J., and Scheuchl, B.: Antarctic grounding line mapping from differential satellite radar interferometry, Geophys. Res. Lett., 38, 1–6, https://doi.org/10.1029/2011GL047109, 2011. 

Rignot, E., Mouginot, J., Morlighem, M., Seroussi, H., and Scheuchl, B.: Widespread, rapid grounding line retreat of Pine Island, Thwaites, Smith, and Kohler glaciers, West Antarctica, from 1992 to 2011, Geophys. Res. Lett., 41, 3502–3509, https://doi.org/10.1002/2014GL060140, 2014. 

Rignot, E., Mouginot, J., and Scheuchl, B.: MEaSUREs Antarctic Grounding Line from Differential Satellite Radar Interferometry, Version 2, [All subsets], NASA National Snow and Ice Data Center Distributed Active Archive Center, https://doi.org/10.5067/IKBWW4RYHF1Q, 2016.  

Simms, A. R., Lisiecki, L., Gebbie, G., Whitehouse, P. L., and Clark, J. F.: Balancing the last glacial maximum (LGM) sea-level budget, Quaternary Sci. Rev., 205, 143–153, https://doi.org/10.1016/j.quascirev.2018.12.018, 2019. 

Stone, J. O., Balco, G. A., Sugden, D. E., Caffee, M. W., Sass, L. C., Cowdery, S. G., and Siddoway, C.: Holocene deglaciation of Marie Byrd Land, West Antarctica, Science, 299, 99–102, https://doi.org/10.1126/science.1077998, 2003. 

Stroeven, A. P., Fabel, D., Harbor, J., Hättestrand, C., and Kleman, J.: Quantifying the erosional impact of the Fennoscandian ice sheet in the Torneträsk-Narvik corridor, northern Sweden, based on cosmogenic radionuclide data, Geogr. Ann. A, 84, 275–287, https://doi.org/10.1111/j.0435-3676.2002.00182.x, 2002. 

Sugden, D. E., Balco, G., Sass, L. C., Cowdery, S. G., and Stone, J. O.: Selective glacial erosion and weathering zones in the coastal mountains of Marie Byrd Land, Antarctica, Geomorphology, 67, 317–334, https://doi.org/10.1016/j.geomorph.2004.10.007, 2005. 

Sugden, D. E., Hein, A. S., Woodward, J., Marrero, S. M., Rodes, A., Dunning, S. A., Stuart, F. M., Freeman, S. P. H. T., Winter, K., and Westoby, M. J.: The million-year evolution of the glacial trimline in the southernmost Ellsworth Mountains, Antarctica, Earth Planet. Sc. Lett., 469, 42–52, https://doi.org/10.1016/j.epsl.2017.04.006, 2017. 

Whitehouse, P. L., Bentley, M. J., and Le Brocq, A. M.: A deglacial model for Antarctica: Geological constraints and glaciological modelling as a basis for a new model of Antarctic glacial isostatic adjustment, Quaternary Sci. Rev., 32, 1–24, https://doi.org/10.1016/j.quascirev.2011.11.016, 2012. 

Whitehouse, P. L., Bentley, M. J., Vieli, A., Jamieson, S. S. R., Hein, A. S., and Sugden, D. E.: Controls on Last Glacial Maximum ice extent in the Weddell Sea embayment, Antarctica, J. Geophys. Res.-Earth, 122, 371–397, https://doi.org/10.1002/2016JF004121, 2017. 

Wolstencroft, M., King, M. A., Whitehouse, P. L., Bentley, M. J., Nield, G. A., King, E. C., McMillan, M., Shepherd, A., Barletta, V., Bordoni, A., Riva, R. E. M., Didova, O., and Gunter, B. C.: Uplift rates from a new high-density GPS network in Palmer Land indicate significant late Holocene ice loss in the southwestern Weddell Sea, Geophys. J. Int., 203, 737–754, https://doi.org/10.1093/gji/ggv327, 2015. 

Download
Short summary
We studied the history of ice masses at three locations in the Weddell Sea Embayment, Antarctica. We measured rare isotopes in material sourced from mountains overlooking the Slessor Glacier, Foundation Ice Stream, and smaller glaciers on the Lassiter Coast. We show that ice masses were between 385 and 800 m thicker during the last glacial cycle than they are at present. The ice masses were both hundreds of metres thicker and remained thicker closer to the present than was previously thought.