Introduction
Cosmogenic nuclides are rare nuclides made in near-surface rocks and minerals by cosmic rays. The concentration of a cosmogenic nuclide in a surface is directly proportional to the time the surface was most recently uncovered by receding ice. As such, measuring cosmogenic nuclide concentrations is a common way of studying glacier chronologies (Schaefer and others, Reference Schaefer2022). By measuring cosmogenic nuclides at different elevations above glaciers, we can constrain both the past thickness and timing and pattern of thinning (Ackert and others, Reference Ackert1999; Stone and others, Reference Stone2003), typically following the Last Glacial Maximum (LGM). These geologic constraints are used to validate the results of numerical ice sheet models investigating deglaciation (e.g., Whitehouse and others, Reference Whitehouse, Bentley and Le Brocq2012; Pittard and others, Reference Pittard, Whitehouse, Bentley and Small2022), informing model parameter selection and ultimately reducing uncertainties when these same models are used to simulate the future response of ice sheets to a changing climate.
Concentrations of cosmogenic nuclides are converted to exposure ages using production rates and, for radioactive nuclides, their half-lives. Most exposure dating studies use 10Be or combine it with 26Al (with half-lives of 1.4 and 0.7 Myr, respectively) (Balco, Reference Balco2011). Exposure dating studies rely on the assumption that concentrations accumulated in a single phase of exposure. Cosmogenic nuclides are predominantly produced in the upper few metres of rock, and we rely on erosion during glaciations to ‘reset’ surfaces. Preserved beneath cold-based (nonerosive) ice, long-lived nuclides like 10Be can persist for multiple glacial-interglacial cycles, breaking the assumption of one period of exposure. Another cosmogenic nuclide, in situ 14C, has a much shorter half-life (5700 ± 30 yr), making concentrations and resulting exposure ages less sensitive to this nuclide ‘inheritance’. The half-life is so short that 14C accumulated prior to the LGM will have decayed away, regardless of how much erosion took place. In situ 14C exposure ages are therefore essentially free of inheritance, providing unambiguous evidence for the timing of glacier thinning or retreat.
Another useful aspect of in situ 14C is the potential to constrain the maximum extent of LGM ice. A balance between production and decay is reached after about 5.5 times the half-life of a radioactive cosmogenic nuclide, at which point a surface is ‘saturated’. This means a surface is saturated with in situ 14C after ≈30 ka of exposure, assuming minimal erosion. When we measure a concentration equivalent to saturation, we know that the sample has been exposed for at least 30 ka, and thus was not covered during the LGM. Hence, surfaces saturated with in situ 14C provide unambiguous evidence for the extent of ice during the LGM. In summary, in situ 14C is useful for studying deglaciation because (i) concentrations are essentially uninfluenced by previous periods of exposure, providing exposure ages that are more likely to reflect the true age of deglaciation when compared with those from long-lived nuclides and (ii) measurements can provide upper constraints on the extent of ice at the LGM, a type of constraint that cannot be provided by measuring long-lived nuclides. Both of these aspects of in situ 14C mean measuring it is particularly useful for benchmarking the results of numerical ice sheet models.
A growing number of laboratories capable of extracting in situ 14C and automation of the extraction process have led to the nuclide being measured at an enhanced rate over the last decade (Fig. 1a). These measurements are advancing our knowledge of the most recent deglaciation in Antarctica, especially where inferences from long-lived nuclides are limited. How in situ 14C is measured, where in Antarctica it has been measured and what these measurements have shown us about deglaciation, are described below. Potential research questions that can be addressed using this nuclide are also outlined, including estimating glacial erosion rates by combining measurements of in situ 14C and 10Be, assessing the outputs of numerical ice sheet models with single measurements of in situ 14C, and investigating and quantifying a proposed Holocene thinning (beneath present) episode with subglacial measurements of in situ 14C. Much of our knowledge of the past of the Antarctic ice sheets is based on periods when the ice sheets were larger than today, such as during the LGM, because evidence for past ice extent is preserved in rock and sediments above, adjacent to and offshore of present ice margins. Knowledge of contracted configurations of the Antarctic ice sheets, gained through subglacial measurements of in situ 14C, is key to understanding the future of the ice sheets given that they are predicted to continue losing mass (DeConto and others, Reference DeConto2021).
How is in situ 14C measured?
The utility of in situ 14C exposure dating has long been known (e.g., Lal, Reference Lal1987, Reference Lal1988; Lal and Jull, Reference Lal and Jull2001) but it was not until the 2010s that improved reproducibility, increased reliability in extraction systems and an accompanying reduction in blank levels, helped make measuring it more routine. Building on methods for measuring in situ 14C in extraterrestrial samples (Goel and Kohman, Reference Goel and Kohman1962; Suess and Wänke, Reference Suess and Wänke1962), Lifton and others (Reference Lifton, Jull and Quade2001) developed the methods for extracting carbon from quartz used in laboratories today. While methods differ with laboratory, the key steps are similar: carbon is liberated through the heating of quartz under vacuum, oxidised to form CO2, then purified using liquid nitrogen. Some extraction lines use a tube furnace and fuse quartz in a lithium metaborate flux (Lifton and others, Reference Lifton, Goehring, Wilson, Kubley and Caffee2015; Lamp and others, Reference Lamp2019; Goehring and others, Reference Goehring, Wilson and Nichols2019a), while others use an electron bombardment or resistance furnace to release in situ 14C by diffusion through the crystal lattice (Fülöp and others, Reference Fülöp, Wacker and Dunai2015, Reference Fülöp2019; Lupker and others, Reference Lupker2019). Samples are sent for AMS measurement as CO2 (Hippe and others, Reference Hippe2013; Lupker and others, Reference Lupker2019) or after dilution and graphitisation (e.g., Lifton and others, Reference Lifton, Goehring, Wilson, Kubley and Caffee2015). Isotope ratios are used to determine in situ 14C concentrations following Hippe and Lifton (Reference Hippe and Lifton2014). In situ 14C concentrations, combined with sample density, thickness, elevation, latitude and longitude and topographic shielding, are then used to calculate exposure ages, usually using an online exposure age calculator such as the online calculators formerly known as the CRONUS-Earth online calculators (Balco and others, Reference Balco, Stone, Lifton and Dunai2008).
The rise in the number of studies applying in situ 14C (Fig. 1a) is fuelled by a number of factors, among which most notable are a growing number of extraction lines, automation of extraction, decreasing blank levels and the widespread adoption of data reduction and production rate calibration methods. Most importantly, an increasing number of laboratories are capable of extracting carbon from quartz. Automation of the extraction process has increased sample throughput, particularly at Tulane University (Goehring and others, Reference Goehring, Wilson and Nichols2019a). A gradual reduction in 14C in process blanks has improved the detection limit. Repeat measurements of the in situ 14C concentration of the interlaboratory comparison material CRONUS-A (Jull and others, Reference Jull, Scott and Bierman2015) have been used to characterise the reproducibility of in situ 14C measurements (approximately 6%; Nichols and others, Reference Nichols2019) and calibrate the production rate used by the online exposure age calculators (Balco and others, Reference Balco, Stone, Lifton and Dunai2008) and the Informal Cosmogenic-Nuclide Exposure-age Database (ICE-D, ice-d.org, Balco, Reference Balco2020). Standardisation of data reduction (Hippe and Lifton, Reference Hippe and Lifton2014) and the identification of a source of contamination from a commonly used method of quartz isolation (Nichols and Goehring, Reference Nichols and Goehring2019) have also contributed to the now relatively routine measurement and application of in situ 14C.
Advances based on in situ 14C
Measurements of in situ 14C are reported from all sectors of Antarctica but are focused in the Ross, Weddell and Amundsen sea embayments, with a dearth of measurements in East Antarctica and few on the Antarctic Peninsula (Fig. 1b). Post-LGM exposure ages constrain deglaciation at most sites, and saturated measurements constrain the limit of LGM ice in the Shackleton Range (Nichols and others, Reference Nichols2019), close to the West Antarctic Ice Sheet (WAIS) Divide (Spector and others, Reference Spector, Stone and Goehring2019) and adjacent to Prydz Bay (Berg and others, Reference Berg2016). Samples saturated with in situ 14C are also observed on blue ice moraines in Queen Maud Land (Akçar and others, Reference Akçar2020). CRONUS-A, a sandstone sample sourced from 1679 m asl in Arena Valley in the Dry Valleys, Antarctica (Jull and others, Reference Jull, Scott and Bierman2015; Fig. 1b), is saturated with 14C and has been measured at least 75 times.
The most obvious places to measure in situ 14C for exposure dating studies are those yielding solely or primarily pre-LGM exposure ages from long-lived nuclides, and thus inferences on the extent of LGM ice are limited. This is the case at the Lassiter Coast in the Weddell Sea Embayment (Fig. 1b), where the majority of 10Be exposure ages of deposits, presumably from the LGM or most recent deglaciation, exceed 100 ka (Fig. 2a). Taken at face value, one could infer that ice has not been thicker here for hundreds of thousands of years, certainly not during the LGM. However, in situ 14C measurements made at the same site, with many of the same samples, yield Holocene exposure ages (Fig. 2a), showing that (i) ice was at least 380 m thicker than present at the LGM, (ii) deglaciation occurred relatively rapidly and (iii) this region was covered by cold-based ice that preserved 10Be that accumulated during previous periods of exposure.
A similar pattern of pre-LGM 10Be exposure ages and post-LGM in situ 14C exposure ages is observed at other sites in the Weddell Sea Embayment. Limited LGM thickening inferred from predominantly pre-LGM 10Be exposure ages in the Shackleton Range (Hein and others, Reference Hein, Fogwill, Sugden and Xu2011) and Pensacola Mountains (Balco and others, Reference Balco2016; Bentley and others, Reference Bentley2017) was used to benchmark ice sheet models for some time (e.g., Whitehouse and others, Reference Whitehouse2017; Kingslake and others, Reference Kingslake2018). These interpretations led to relatively little post-LGM ice volume change in the Weddell Sea Embayment (when compared with previous reconstructions, see Bentley and Anderson (Reference Bentley and Anderson1998)) becoming the predominant reconstruction among the palaeo community (Hillenbrand and others, Reference Hillenbrand2014). Subsequent measurements of in situ 14C yielded post-LGM exposure ages at both locations, showing that, rather than limited thickening, ice was at least 310 and 800 m thicker than present at the LGM (Nichols and others, Reference Nichols2019). Other locations with multiple samples yielding pre-LGM exposure ages from long-lived nuclides and post-LGM in situ 14C exposure ages are the Flower Hills and Meyer Hills in the Ellsworth Mountains (Fogwill and others, Reference Fogwill2014) and the Darwin–Hatherton Glacier System in the Ross Sea Embayment (Hillebrand and others, Reference Hillebrand2021).
In situ 14C can also be useful at sites yielding post-LGM 10Be exposure ages. For example, in the Amundsen Sea Embayment, Johnson and others (Reference Johnson2020) use measurements of in situ 14C to identify a smaller degree of inheritance in their exposure ages. Here, 10Be exposure ages (n = 9) indicate deglaciation happened about 17 ka, while in situ 14C exposure ages (n = 8) show it occurred about 6 ka, a difference of 11 ka, which is significant for establishing an accurate deglacial chronology of the region. Samples with post-LGM 10Be exposure ages and younger in situ 14C exposure ages are also observed at sites in northern Victoria Land (Balco and others, Reference Balco, Todd, Goehring, Moening-Swanson and Nichols2019; Goehring and others, Reference Goehring, Balco, Todd, Moening-Swanson and Nichols2019b), with an additional sample in the Flower Hills (Fogwill and others, Reference Fogwill2014). Evidently, even when 10Be exposure ages at a site postdate the LGM and thus we know the degree of ice thickness change, there could still be a detectable amount of inheritance skewing our understanding of the timing of deglaciation.
Glacier chronologies are constrained solely with in situ 14C measurements (without accompanying long-lived nuclides) at some locations, such as the Whitmore Mountains close to WAIS Divide (Spector and others, Reference Spector, Stone and Goehring2019) and some sites in northern Victoria Land (Goehring and others, Reference Goehring, Balco, Todd, Moening-Swanson and Nichols2019b). Additionally, concordant in situ 14C and 10Be exposure ages are observed at many sites in Antarctica (White and others, Reference White, Fülöp, Bishop, Mackintosh and Cook2011; Balco and others, Reference Balco, Todd, Goehring, Moening-Swanson and Nichols2019; Goehring and others, Reference Goehring, Balco, Todd, Moening-Swanson and Nichols2019b; Hillebrand and others, Reference Hillebrand2021).
Future research priorities
Measuring in situ 14C
While we have learnt much about deglaciation in Antarctica from in situ 14C in recent years, we have also learnt much about measuring the nuclide itself, and some questions remain unanswered. Some studies observe in situ 14C concentrations in excess of theoretical limits (Balco and others, Reference Balco2016; Akçar and others, Reference Akçar2020), while another observes measurement reproducibility lower than that expected from measurement uncertainties alone (Nichols and others, Reference Nichols2019). When sample contamination can be ruled out, mass movement and supraglacial transport could explain elevated concentrations (Balco and others, Reference Balco2016), while unrecognised measurement error could explain the limited reproducibility. Further work dedicated to method development is needed to isolate what is (i) limiting measurement reproducibility and (ii) contributing toward concentrations exceeding theoretical limits. Most studies measure in situ 14C in quartz, but the nuclide is produced in other materials such as calcium carbonate (Handwerger and others, Reference Handwerger, Cerling and Bruhn1999) and olivine (Pigati and others, Reference Pigati, Lifton, Jull and Quade2010). Establishing methods for the extraction of carbon from these materials would expand the number of locations we can study with in situ 14C beyond only those rich in quartz.
Applying in situ 14C
Further exposure dating studies using in situ 14C would be useful in areas unstudied with cosmogenic nuclides or those yielding solely or primarily pre-LGM exposure ages from long-lived nuclides (e.g., Hodgson and others, Reference Hodgson2012). By filling spatial gaps in our knowledge of deglaciation in Antarctica, these measurements of in situ 14C would provide more geologic constraints for the benchmarking of numerical ice sheet models. Additionally, there are a few applications beyond traditional exposure dating yet to be used (or to their full potential) that could improve our knowledge of the history and glaciology of the Antarctic ice sheets.
How much did glaciers erode into bedrock during the Holocene? This question can be answered by combining measurements of in situ 14C and 10Be in recently exposed proglacial bedrock. Making direct measurements of subglacial erosion is complicated by the difficulty of accessing the beds of glaciers. Using cosmogenic nuclides to estimate past erosion rates provides knowledge of glacial processes over a longer period than contemporary point measurements, extending our knowledge of glacial erosion rates beyond the observational record. The relationship between the in situ 14C and 10Be concentration of a proglacial bedrock sample is related to the depth to which a glacier eroded into bedrock in the Holocene, allowing the estimation of Holocene glacial erosion rates (Rand and Goehring, Reference Rand and Goehring2019). Because this method requires proglacial bedrock, it may be limited to smaller glaciers such as those on the Antarctic Peninsula or at high elevations on the continent.
How closely do numerical ice sheet model outputs reflect the timing of both the advance and retreat phases of deglaciation contained in geologic archives? Measurements of in situ 14C, rather than long-lived nuclides, can answer this question. Many numerical ice sheet models are benchmarked against exposure age datasets recording only deglaciation. By assuming samples were saturated prior to LGM burial, individual measurements of in situ 14C can be used to assess the timing of both advance and retreat phases of model outputs (Spector and others, Reference Spector, Stone and Goehring2019), reducing uncertainties when these same models are used to simulate the future of ice sheets. More generally, targeting exposed surfaces high above modern ice elevations could help provide more upper constraints on LGM ice thicknesses to help validate numerical ice sheet model outputs.
To what extent, and where, did parts of the Antarctic ice sheets readvance in the Holocene? This is perhaps the most important question that in situ 14C can answer. A number of studies, both through geologic observations (Siegert and others, Reference Siegert, Ross, Corr, Kingslake and Hindmarsh2013; Wolstencroft and others, Reference Wolstencroft2015; Greenwood and others, Reference Greenwood, Simkins, Halberstadt, Prothro and Anderson2018; King and others, Reference King, Watson and White2022) and modelling (Kingslake and others, Reference Kingslake2018), infer that some parts of the Antarctic ice sheets were smaller than present in the Holocene and subsequently grew to their present configuration. Through measuring carbon in subglacial sediments, two studies (Venturelli and others, Reference Venturelli2020 and Neuhaus and others, Reference Neuhaus2021) report the first direct evidence of a Holocene grounding line readvance in the Ross sector. Further direct evidence for a Holocene readvance can be obtained through in situ 14C measurements in subglacial bedrock, because significant concentrations in subglacial bedrock unambiguously requires Holocene exposure, either complete or through relatively thin ice (Johnson and others, Reference Johnson2022). Constraining the scale of this readvance, both in ice thickness change and geographic extent, could shed light on the processes causing the mass loss and subsequent gain (e.g., ocean forcings, glacioisostatic adjustment), information that can then be used with numerical models to replicate this ice sheet behaviour. Given that current Antarctic ice sheet mass loss is predicted to continue (DeConto and others, Reference DeConto2021), knowing the processes that helped recover ice mass loss in a climate relatively similar to that of today is key to understanding the reversibility of current and future Antarctic ice mass loss. While previous studies have investigated long term changes in the Greenland Ice Sheet by measuring long-lived nuclides in subglacial material (Schaefer and others, Reference Schaefer2016; Christ and others, Reference Christ2021), there are no published subglacial measurements of in situ 14C from beneath any ice sheet. If above background in situ 14C indicative of a Holocene readvance is measured in samples collected from beneath the Antarctic ice sheets, multiple studies will be required to confirm if this ice sheet behaviour is widespread or localised.
Conclusions
To summarise, the cosmogenic nuclide in situ 14C has been measured at an enhanced rate over the last decade, fuelled by the automation of the extraction process and an increasing number of laboratories now capable of extracting it. Measurements of in situ 14C have been used in exposure dating studies to shed light on deglaciation in all sectors of Antarctica, but especially in the Weddell Sea Embayment. Some studies observe in situ 14C concentrations exceeding theoretical limits and also measurement reproducibility lower than expected, which can hopefully be addressed with dedicated work on understanding the extraction process and geomorphic scatter. While there are many locations in Antarctica where traditional in situ 14C exposure dating studies would be useful, there are also a number of other applications of the nuclide that hold much potential, including using subglacial measurements to constrain episodes of thinning and rethickening in the Holocene.