Hostname: page-component-78c5997874-4rdpn Total loading time: 0 Render date: 2024-11-18T11:15:58.380Z Has data issue: false hasContentIssue false

Palaeo-environmental significance of evaporative calcite crusts in the Untersee Oasis, East Antarctica

Published online by Cambridge University Press:  17 May 2024

Denis Lacelle*
Affiliation:
Department of Geography, Environment and Geomatics, University of Ottawa, Ottawa, ON, Canada
Michelle Christy
Affiliation:
Department of Geography, Environment and Geomatics, University of Ottawa, Ottawa, ON, Canada
Benoit Faucher
Affiliation:
Geological Survey of Canada, Ottawa, ON, Canada
Pablo Sobron
Affiliation:
Carl Sagan Center, SETI Institute, Mountain View, CA, USA
Dale Andersen
Affiliation:
Carl Sagan Center, SETI Institute, Mountain View, CA, USA
Rights & Permissions [Opens in a new window]

Abstract

Secondary carbonate precipitates on the surface of clasts have rarely been reported from Antarctica. Here, we infer the origin, age and palaeo-environmental significance of the calcite crusts in the Untersee Oasis, East Antarctica. The distribution of calcite crusts, which are up to 1 mm thick, is limited to locations with residual snow patches, and they have some of the highest δ18O (up to +17.4‰ Vienna Standard Mean Ocean Water (VSMOW)) and δ13C (up to +14.6‰ Vienna Pee Dee Belemnite (VPDB)) compositions of any carbonate deposits in terrestrial polar environments. Their δ18O and δ13C values are substantially enriched with respect to the isotopic values expected from equilibrium precipitation from the δ18O and δ13CDIC (DIC = dissolved inorganic carbon) of snow meltwater. The formation of the calcite crusts is ascribed to the evaporation of residual snow meltwater and the low relative humidity and strong winds, favouring a kinetic isotope effect. The 14C age distribution of the calcite crusts (1550 cal yr bp to modern) provides a minimum age for ice retreat and drainage of the palaeo-lake in Aurkjosen Cirque. However, in this polar desert environment in which surface melting is limited, the calcite crusts require sufficient snow accumulation and air temperatures warm enough to generate meltwater, and their age distribution corresponds to the late Holocene warm-wet climate period.

Type
Earth Sciences
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
Copyright © The Author(s), 2024. Published by Cambridge University Press on behalf of Antarctic Science Ltd

Introduction

Secondary carbonate deposits in terrestrial polar regions are common in limestone environments. They typically occur on the surface or underside of clasts (Hanshaw & Hallet Reference Hanshaw and Hallet1978, Hillaire-Marcel et al. Reference Hillaire-Marcel, Soucy and Cailleux1979, Souchez & Lemmens Reference Souchez and Lemmens1985, Courty et al. Reference Courty, Marlin, Dever, Tremblay and Vachier1994), within fissures in bedrock outcrops (Lauriol & Clark Reference Lauriol and Clark1999) and on the surface of ice in caves and icings (Lauriol et al. Reference Lauriol, Cinq Mars and Clark1991, Clark & Lauriol Reference Clark and Lauriol1997, Omelon et al. Reference Omelon, Pollard and Marion2001, Lacelle et al. Reference Lacelle, Lauriol and Clark2006). They are less prominent in non-limestone settings, but they have still been observed in crystalline and granitic bedrock environments on the surface of clasts (Blake Reference Blake1999, Lacelle Reference Lacelle2007). If the carbonate deposits precipitate under equilibrium conditions, their δ18O composition can provide palaeo-climatic information (e.g. Clark et al. Reference Clark, Lauriol, Marschner, Sabourin, Chauret and Desrochers2004, Lacelle Reference Lacelle2007); alternatively, if they precipitated from waters in equilibrium with atmospheric CO2, they can be used as palaeo-environmental proxies. For example, the 14C ages of carbonate crusts on the surface of moraines were used to infer the timing of deglaciation in east-central Ellesmere Island (high Arctic Canada; Blake Reference Blake2005), southern Baffin Island (low Arctic Canada; Lacelle et al. Reference Lacelle, Lauriol and Clark2007) and northern Pakistan (Waragai Reference Waragai2005), and the chronologies matched well with those derived from other methods.

Reports of carbonate precipitates on clasts or in soils have been spatially limited in Antarctica. In the McMurdo Dry Valleys, paedogenic carbonates were observed in the soils of multiple valleys (Claridge & Campbell Reference Claridge and Campbell1977, Campbell & Claridge Reference Campbell, Claridge and Margesin2009, Campbell et al. Reference Campbell, Claridge, Campbell and Balks2013, Lyons et al. Reference Lyons, Foley, Carey, Diaz, Bowen and Cerling2020), in soils near Shackleton Glacier (Diaz et al. Reference Diaz, Li, Michalski, Darrah, Adams and Wall2020) and in carbonate crusts near the shore of Lake Vanda (Nakai et al. Reference Nakai, Wada, Kiyosu and Takimoto1975). The paedogenic calcites were used to inform about chemical weathering (i.e. sources of calcium and carbon in the soils). Here, we report the discovery of carbonate crusts on the surface of clasts in the polar desert of the Untersee Oasis in East Antarctica. The carbonate crusts occur on morainic material around Lake Untersee and in Aurkjosen Cirque, including within a former lake basin. This study aims to determine the origin and age of the carbonate crusts in the Untersee Oasis to provide information regarding the timing of glacial retreat, drainage of a proglacial lake and climate conditions that would allow for their formation. This objective is reached by determining the mineralogy of the carbonates, their δ18O and δ13C composition to infer their process of formation and the 14C ages to infer the timing of their growth. Residual snow patches were also sampled and analysed for major ions to determine their saturation indices with respect to carbonate minerals and δ18O to assess whether the carbonates precipitated under equilibrium or kinetic conditions.

Study area

Untersee Oasis (71.3°S, 13.5°E) is located within the Gruber Mountains of Queen Maud Land, ~150 km from the coast and ~90 km south-east of the Schirmacher Oasis (Fig. 1). The Untersee Oasis includes three main ice-free regions: 1) Lake Untersee valley, an 11 km-long and 4 km-wide north-south trending valley, 2) Aurkjosen Cirque, a ~3.5 km-long and 2 km-wide east-west trending valley and 3) Pritzker Valley (informal name), a 2.5 km-long a 0.5-km wide north-east-south-west trending valley with a small ice patch near the head of the valley.

Figure 1. Map of the Untersee Oasis in Queen Maud Land, East Antarctica. (Top) The locations of sampling sites of the calcite crusts and inferred glacial limits (from Schwab Reference Schwab1998). Ages of the calcite crusts are provided as cal yr bp. (Bottom) Recessional glacial limits in the Untersee Oasis and palaeo-lake basin in Aurkjosen Cirque (from Schwab Reference Schwab1998).

The East Antarctic Ice Sheet covered the Untersee Oasis during the Late Pleistocene. Based on the 14C ages of sub-fossilized stomach oils from snow petrel nests, thinning of the ice sheet began at c. 35–30 kyr bp (40 160–34 450 cal yr bp; Hiller et al. Reference Hiller, Wand, Kämpf and Stackebrandt1988, Reference Hiller, Hermichen and Wand1995), which led to a reconfiguration of the local ice flow. According to Schwab (Reference Schwab1998), proglacial Lake Untersee probably formed between 13 and 11 kyr bp (15 560–13 000 cal yr bp) when the Anuchin Glacier started to retreat northwards. During its retreat, a proglacial lake also occupied Aurkjosen Cirque, as evidenced by palaeo-shorelines. That lake drained and/or evaporated to dryness once Anuchin Glacier receded from that valley c. 7 kyr bp (7800 cal yr bp) along a meltwater channel into Lake Untersee (i.e. the current lake ice ablation rate is 0.4–0.6 m yr-1; Faucher et al. Reference Faucher, Lacelle, Fisher, Andersen and McKay2019). The Untersee Oasis probably came to its current configuration at c. 6–4 kyr bp.

The local geology in the Oasis consists of Precambrian norite, anorthosite and anorthosite-norite alternation of the Eliseev massif complex (Kampf & Stakebrandt Reference Kampf and Stakebrandt1985, Bormann et al. Reference Bormann, Bankwitz, Bankwitz, Damn, Hurtig and Kampf1986, Paech & Stackebrandt Reference Paech, Stackebrandt, Bormann and Fritzsche1995). There are no units of carbonate rocks in the area. The surface sediments consist mainly of till and colluvium, often covered by a thin layer of aeolian sediments (Schwab Reference Schwab1998). Vegetation and lichens are absent, and the soils consist of poorly sorted sediments with very low organic content (Shamilishvili et al. Reference Shamilishvili, Abakumov, Andersen, Frank-Kamanetskaya, Vlasov, Panova and Lessovaia2020).

The region is characterized by a polar desert climate regime. Ten years of climate data (2008–2017) collected by an automated weather station along the shore of Lake Untersee (71.34°S, 13.45°E, 612 m above sea level) shows a mean annual air temperature (MAAT) of -9.5 ± 0.7°C, thawing degree-days ranging from 7 to 51 degree-days and a mean relative humidity of 42 ± 5% (Andersen et al. Reference Andersen, McKay and Lagun2015, Faucher et al. Reference Faucher, Lacelle, Fisher, Andersen and McKay2019). Despite having a relatively warm MAAT for Antarctica, the climate in the Oasis is dominated by intense ablation, which limits surface melt features (no surface runoff) due to cooling associated with latent heat of sublimation and vaporization (e.g. van den Broeke et al. Reference Van den Broke, Jan van de Berg, van Meijgaard and Reijmer2006; Hoffman et al. Reference Hoffman, Fountain and Liston2008).

Methods

White carbonate crusts were collected in November–December of 2015 and 2017 from the surface of clasts within the morainic material surrounding Lake Untersee and in Aurkjosen Cirque; in the latter, crusts were also collected on the surface of clasts within the palaeo-lake basin (Figs 1 & 2). These crusts are all found in the immediate vicinity of residual snow patches and are different from the carbonate-poor greyish lithificates that occur only along the palaeo-shorelines of Lake Untersee (Levitan et al. Reference Levitan, Kononkova, Luksha and Roshchina2012). The carbonate crusts, up to 1 mm thick, can be found on clasts up to 1 m above the surface of the soils (probably reflecting the height of residual snow patches), and they were collected with a chisel and hammer and placed in sealed Falcon tubes.

Figure 2. Field photographs of calcite crusts in the Untersee Oasis, East Antarctica. The calcite crusts are up to 1 mm thick.

The mineral compositions of five carbonate crusts were determined in the field using Raman spectroscopy. The excitation source was a frequency-doubled, Q-switched neodymium-doped yttrium aluminium garnet (Nd:YAG) pulsed laser source, and the spectra were analysed and recorded with a Kaiser Optical System fitted with an intensified charged-coupled device camera.

The δ18O and δ13C compositions of 17 carbonate crusts were determined to infer their origin (process of formation). The measurements of 18O/16O and 13C/12C ratios were made on CO2 gas produced by reacting the powdered crusts with 100% phosphoric acid in glass septum vials for 24 h at 25°C. The evolved CO2 gas was analysed in continuous flow using a Gas Bench II interfaced with a Finnigan Mat Delta+ XP isotope mass spectrometer at the Jan Veizer Laboratory (University of Ottawa, Canada). Stable isotope data for O and C were expressed in δ-notation, where δ represents the part per thousand difference of 18O/16O and 13C/12C in the sample with respect to the Vienna Pee Dee Belemnite (VPDB) standard (Kim et al. Reference Kim, Coplen and Horita2015). To facilitate comparison with the residual snow patches, the δ18O composition of the carbonate crusts were converted to the Vienna Standard Mean Ocean Water (VSMOW) scale (δ18OVSMOW = 1.0309 δ18OVPDB + 30.92; Coplen et al. Reference Coplen, Hopple, Böhlke, Peiser, Rieder and Krouse2002). Analytical reproducibility is 0.15‰ for both isotopes.

Radiocarbon measurements of 13 carbonate crusts were made to infer their carbon source and ages. The measurements were made at the Andre Lalonde AMS Laboratory (University of Ottawa, Canada) using a 3MV tandem accelerator mass spectrometer (AMS). The 14C/12C ratios are expressed as fraction of modern carbon (F14C) and corrected for spectrometer and preparation fractionation using the AMS-measured 13C/12C ratio (Crann et al. Reference Crann, Murseli, St-Jean, Zhao, Clark and Kieser2017). Radiocarbon ages are calculated as -8033ln(F14C) and reported in 14C yr bp (bp = ce 1950; Stuiver & Polach Reference Stuiver and Polach1977), and the 14C ages were calibrated to cal yr bp using OxCal and the IntCal20 calibration curve (Bronk Ramsey Reference Bronk Ramsey2009, Reimer et al. Reference Reimer, Austin, Bard, Bayliss, Blackwell and Bronk Ramsey2020).

Samples of residual snow patches in the vicinity of the carbonate crusts were also collected in November–December 2015 and 2017 to determine their geochemical and δ18O compositions. The concentrations of the major cations (Ca2+, K+, Mg2+, Na+) and anions (SO42-, Cl-, NO3-) of the filtered snow meltwater were measured using an Agilent 4200 inductively coupled plasma atomic emission spectrometer and ion chromatography (Dionex ICS-2100), respectively. Prior to the analysis of cations, water samples were acidified to pH 2 with ultra-pure HNO3 acid. Analytical reproducibility for solute analysis is ±5%. The 18O/16O and D/H ratios of the meltwater were determined using a Los Gatos Research liquid water analyser coupled to a CTC LC-PAL autosampler for simultaneous 18O/16O and D/H ratios measurements of H2O and verified for spectral interference contamination. The results are presented using the δ-notation (δ18O and δD), where δ represents the parts per thousand differences for 18O/16O or D/H in a sample with respect to VSMOW. The analytical reproducibility values for δ18O and δD are ±0.3‰ and ±1‰, respectively.

Results

Based on Raman spectroscopy and the sharp Raman band at 1088 cm-1, the carbonate crusts in the Untersee Oasis consist of calcite minerals (Fig. 3). The δ18O and δ13C of the majority calcite crusts range from +10.8 to +17.4‰ and +9.6 to +14.6‰, respectively; however, two of the 17 calcite crusts had negative δ18O and δ13C values (Fig. 4a & Table I). The calcite crusts have median bulk 14C ages in the 1589 cal yr bp to modern (2011 cal yr ce) range (Table II). The calcite crusts on the peninsula were dated from 1589 to 545 cal yr bp, those in Aurkjosen Cirque were dated from 683 cal yr bp to modern and those above the shoreline of Lake Untersee were the youngest (119 cal yr bp to modern; Fig. 1).

Figure 3. Raman spectra of five carbonate crusts in the Untersee Oasis, East Antarctica. The intense sharp Raman band at 1088 cm-1 is indicative of the mineral calcite.

Figure 4. δ13C and δ18O compositions of the calcite crusts in the Untersee Oasis, East Antarctica. a. δ13C and δ18O compositions of the calcite crusts in the Untersee Oasis compared to others in Antarctica: lithificates along palaeo-shorelines of Lake Untersee (Levitan et al. Reference Levitan, Kononkova, Luksha and Roshchina2012), evaporative calcite crusts surrounding Lake Vanda (Nakai et al. Reference Nakai, Wada, Kiyosu and Takimoto1975) and paedogenic calcite in soils surrounding Shackleton Glacier (Diaz et al. Reference Diaz, Li, Michalski, Darrah, Adams and Wall2020). b. Comparison of δ13C and δ18O compositions of the calcite crusts in the Untersee Oasis (red box) with other terrestrial polar environments (from Lacelle Reference Lacelle2007). VPDB = Vienna Pee Dee Belemnite; VSMOW = Vienna Standard Mean Ocean Water.

Table I. δ13C and δ18O measurements of calcite crusts in the Untersee Oasis, East Antarctica. The δ18O composition of the carbonate crusts were converted to the Vienna Standard Mean Ocean Water (VSMOW) scale according to: δ18OVSMOW = 1.0309 δ18OVPDB + 30.92 (Coplen et al. Reference Coplen, Hopple, Böhlke, Peiser, Rieder and Krouse2002).

VPDB = Vienna Pee Dee Belemnite.

Table II. Radiocarbon results from calcite crusts collected in the Untersee Oasis, East Antarctica. The 14C ages were calibrated using OxCal and the IntCal20 calibration curve (Bronk Ramsey Reference Bronk Ramsey2009, Reimer et al. Reference Reimer, Austin, Bard, Bayliss, Blackwell and Bronk Ramsey2020).

F14C = fraction of modern carbon.

The snow meltwater of residual snow patches has a Ca-Na-SO4 geochemical facies with [Ca2+] = 6.0 mg l-1, [Na+] = 5.0 mg l-1 and [SO42-] = 9.1 mg l-1. This geochemical facies is similar to that of regional snow (Isaksson et al. Reference Isaksson, Karlén, Gundestrup, Mayewski, Whitlow and Twickler1996, Marsh et al. Reference Marsh, Lacelle, Faucher, Cotroneo, Jasperse, Clark and Andersen2020) and of local ponds recharged by snowmelt (Faucher et al. Reference Faucher, Lacelle, Marsh, Fisher and Andersen2021). Based on the Phreeqc hydrogeochemical program, and assuming that the meltwater is in equilibrium with atmospheric CO2, the meltwater is undersaturated with respect to calcite (SIcal = -4.5). The snow meltwater has δ18O values averaging -33.5 ± 2.4‰.

Discussion

Origin of the calcite crusts

All except two of the calcite crusts collected from the surface of clasts in the Untersee Oasis have some of the highest reported δ18O and δ13C values of any type of calcite precipitates found in terrestrial polar environments (Fig. 4b). The δ18O and δ13C values are much higher than the carbonate-poor lithificates observed along the palaeo-shorelines of Lake Untersee (δ18Oavg = +4.4‰, δ13Cavg = -3.3‰; Levitan et al. Reference Levitan, Kononkova, Luksha and Roshchina2012). However, they are in the range of the paedogenic carbonates in the McMurdo Dry Valleys (δ13C in the +1.3 to +11.0‰ range) and the crusts on the dried lakebed in the basin of Lake Vanda (δ13C values in the +7.8 to +17.6‰ range; Fig. 4b).

Under most circumstances, the δ18O and δ13C of calcites can be estimated from the δ18O and δ13CDIC (DIC = dissolved inorganic carbon) of the water from which the calcite precipitated and the temperature at which the precipitation occurred (Lacelle Reference Lacelle2007). The calcite crusts in the Untersee Oasis probably precipitated at a temperature near 0°C (maximum hourly summer air temperatures are 9°C). Under equilibrium conditions (ɛCaCO3-H2O = 33.6‰ at 0°C; Kim & O'Neil Reference Kim and O'Neil1997), the δ18O values of the calcite crusts (+10.8 to +17.4‰ range) are enriched by 10–17‰ relative to the average δ18O composition of snow meltwater (–33.5 ± 2.4‰). Although the δ13CDIC of snow meltwater was not analysed, it can be assumed that its composition in a plagioclase setting would be in equilibrium with atmospheric CO2. Over the past 2000 years, the concentration of CO2 increased from 270 to 410 ppm, and the δ13C of CO2 evolved from -6 to -8‰ (Schmitt et al. Reference Schmitt, Schneider, Elsig, Leuenberger, Lourantou and Chappellaz2012, Eggleston et al. Reference Eggleston, Schmitt, Bereiter, Schneider and Fischer2016). Under equilibrium conditions, the pH of snow meltwater for these conditions for CO2 would be 6.4–6.6, with a δ13CDIC of -2.8 to -1.4‰ (εCO2aq-CO2g = -1.2‰; εHCO3-CO2aq = 10.9‰). This would produce calcite with δ13C of +0.8 to +2.1‰ (εCaCO3-HCO3 = 3.6‰), which is much lower than the δ13C value measured from the calcite crusts in the Untersee Oasis (+9.6 to +14.6‰).

There are only a few processes capable of generating calcite precipitates with enriched δ18O and δ13C compositions over those of the parent water in polar environments: 1) kinetic freezing or 2) kinetic evaporation. According to Clark & Lauriol (Reference Clark and Lauriol1992), kinetic freezing imparts an enrichment in the order of 31.2 ± 3.1‰ for δ13C (ε13CKIE CaCO3-CO2; KIE = kinetic isotope enrichment) and of 36.7 ± 1.3‰ for δ18O (ε18OKIE CaCO3-H2O). However, kinetic freezing has only been observed to develop cryogenic calcite powders, as observed in caves (Clark & Lauriol Reference Clark and Lauriol1992) and icings (Lauriol et al. Reference Lauriol, Cinq Mars and Clark1991, Lacelle et al. Reference Lacelle, Lauriol and Clark2006). As a result, this process can be ruled out because in the Untersee Oasis the calcites developed crusts on the surface of clasts and not powders.

Kinetic evaporation is the only process capable of explaining the enrichment in both δ18O and δ13C in the calcite crusts in the Untersee Oasis. The effect of evaporation on the δ18O of the residual snow meltwater and the calcite saturation index in the Untersee Oasis was modelled using the IsoVap7 computer code (appendix A in Fisher et al. Reference Fisher, Lacelle, Pollard and Faucher2020). The open system aqueous solution evaporation model uses the Freezchem database and accounts for isotopic exchanges with the ambient moisture based on the model of Sofer & Gat (Reference Sofer and Gat1975) that includes isotope salt effects. The model with nearly pure water chemistry is similar to that described in Criss (Reference Criss1999) and has been used in Antarctic studies (e.g. Lapalme et al. Reference Lapalme, Lacelle, Pollard, Fisher, Davila and McKay2017, Faucher et al. Reference Faucher, Lacelle, Marsh, Fisher and Andersen2021). Despite the relative humidity in the air being 42 ± 5%, the relative humidity at the ground surface is consistently higher than in the air (often > 70%; D. Lacelle et al., unpublished data). This is consistent with studies from the McMurdo Dry Valleys that reported ground surface relative humidities higher than in the air (e.g. Fisher et al. Reference Fisher, Lacelle, Pollard, Davila and McKay2016, McKay et al. Reference McKay, Balaban, Abrahams and Lewis2019, Marinova et al. Reference Marinova, McKay, Heldmann, Goordial, Lacelle, Pollard and Davila2022). In all numerical scenarios (Fig. 5), evaporation at ground surface relative humidity < 70% caused strong δ18O enrichment in the residual water, with enrichment in the order of 30–40‰ when most of the water had evaporated. However, for ground surface relative humidity > 85%, the equilibrium exchange between snowmelt and atmospheric moisture approached a steady-state value that is dependent on the δ18O of the moisture and relative humidity. If we assume that the snow meltwater from which the calcite crusts formed had initial δ18O near -33‰ and the δ18O water vapour was in the -44 to -37‰ range, then at least 25% of the meltwater must evaporate to reach the values measured in the calcites (Fig. 5). During evaporation of snowmelt in equilibrium with P CO2 of -3.5, calcite saturation also progressively increases as the solutes are concentrated in the residual water. Given the chemistry of the snow, calcite saturation is expected when 85% of the water has evaporated. As such, it is more likely that ground surface relative humidity was > 70% during the evaporation of snowmelt (Fig. 5). This would allow for the δ18O of the residual snow meltwater to reach a steady-state value, whereas continued evaporation would allow the solution to reach calcite saturation. For the effect of evaporation on the δ13C composition of calcites, caliche that developed on the surface of basalts in the Arizona volcanic field had δ13C values in the 4–12‰ range, and Knauth et al. (Reference Knauth, Brilli and Klonowski2003) proposed that their elevated values could be produced from the degassing of CO2 and an unknown isotope effect during evaporation. Subsequently, Lacelle et al. (Reference Lacelle, Lauriol and Clark2007) experimentally measured the 13C enrichment during kinetic evaporation. It was found that the kinetic isotope effect between CO2 and calcite ranged between 20 and 40‰ (with the kinetic isotope effect between HCO3 and calcite ranging between 11.7 and 32.1‰). If we assume that the δ13CDIC of snow meltwater was in equilibrium with atmospheric CO213CDIC in the -2.8 to -1.4‰ range), then the elevated δ13C of the calcite crusts can be explained by kinetic evaporation once ~85% of the snow meltwater has evaporated.

Figure 5. Evolution of δ18O of evaporating water (snow meltwater chemistry) for a range of relative humidities (RHs) using the Criss (Reference Criss1999) and Sofer & Gat (Reference Sofer and Gat1975) models. The initial δ18O composition of snow was set at -33‰ (average of the measurements) and that of the atmospheric moisture above the evaporating surface was set at -44‰ (equilibrium with meltwater) and -37‰, and the RH at the ground surface ranged from 50 to 95% (values near 50% represent atmospheric humidity, those > 70% probably represent those at the ground surface; i.e. Fisher et al. Reference Fisher, Lacelle, Pollard, Davila and McKay2016). VSMOW = Vienna Standard Mean Ocean Water.

In the Untersee Oasis, the calcite crusts are present on the surface of clasts (gravels to large boulders) in the immediate vicinity of residual snow patches within morainic material (Fig. 2). There, a portion of the residual snow patches can melt during the summer and interact with the plagioclase rocks, and, as it wets the surfaces of clasts, the meltwater can precipitate calcite following evaporation. The snow patches contain little dissolved organic carbon (DOC) that could increase P CO2 and decrease the δ13CDIC in the meltwater (DOC concentration in nearby ponds recharged by snowmelt range from < 0.3 to 1 ppm; Faucher et al. Reference Faucher, Lacelle, Marsh, Fisher and Andersen2021). The Untersee Oasis, with its low relative humidity and strong winds, is dominated by intense ablation, and these conditions would favour an increase in the rate of evaporation and create a kinetic isotope effect during the precipitation of the calcite crusts. This effect is caused by the rapid degassing during evaporation of a thin film of meltwater, such that there is preferential outgassing of CO2 as Ca and HCO3 concentrate prior to and during calcite precipitation. These are different from the lithificates found solely below the palaeo-shorelines of Lake Untersee (Levitan et al. Reference Levitan, Kononkova, Luksha and Roshchina2012) and from the subglacially precipitated calcites found in formerly glaciated regions (e.g. Hillaire-Marcel et al. Reference Hillaire-Marcel, Soucy and Cailleux1979, Souchez & Lemmens Reference Souchez and Lemmens1985). The latter tend to develop on the leeside of bedrock obstacles following regelation, are often striated and tend to be in isotopic equilibrium with the δ18O and δ13CDIC of the basal meltwater from which they precipitated.

Palaeo-environmental significance

Based on the 14C ages of sub-fossilized stomach oils from snow petrel nests, thinning of the ice sheet in the Untersee Oasis region began at c. 40–35 cal yr bp (Hiller et al. Reference Hiller, Wand, Kämpf and Stackebrandt1988, Reference Hiller, Hermichen and Wand1995). This led to a reconfiguration of the local ice flow in the Untersee Oasis and the formation of proglacial Lake Untersee c. 15 560–13 000 cal yr bp (Fig. 1). During the deglaciation period in the Untersee Oasis, a proglacial lake also occupied Aurkjosen Cirque, as evidenced by palaeo-shorelines. The timing of the disappearance of this lake in Aurkjosen Cirque is unknown (drainage and/or evaporation to dryness once Anuchin Glacier receded from that valley c. 7800 cal yr bp). Our 14C ages of calcite crusts can be used to infer the minimum age of ice retreat and lake disappearance in Aurkjosen Cirque (Fig. 6). Considering that some of the bulk 14C ages of the calcite crusts were modern (Table II), we can assume that their 14C content was nearly 100% at the time of their precipitation and can provide a minimum age of carbonate precipitation. In Aurkjosen Cirque, the calcite crusts have ages ranging from 683 cal yr bp to modern. This would suggest that the palaeo-lake basin in Aurkjosen Cirque came into existence at least 750 years ago. However, in this polar desert environment, the calcite crusts require sufficient snow accumulation and temperatures warm enough to generate meltwater to allow for their formation. The calcite crusts all have median 14C ages < 1589 cal yr bp, with most being < 800 cal yr bp. Over the Holocene period, the past 1500 years corresponds to a warm-wet climate interval (Fig. 6; Buizert et al. Reference Buizert, Fudge, Roberts, Steig, Sherriff-Tadano and Ritz2021). As such, even if the ice retreated from Aurkjosen Cirque and the lake drained by the Middle Holocene, climate conditions probably were not favourable for the development of calcite crusts until the late Holocene.

Figure 6. Comparison of radiocarbon age distributions obtained from evaporative calcite crusts in the Untersee Oasis, East Antarctica, with Holocene reconstructed temperatures and ice accumulation rates (Buizert et al. Reference Buizert, Fudge, Roberts, Steig, Sherriff-Tadano and Ritz2021).

Conclusions

Based on these results, the following conclusions can be drawn regarding the origin, age and palaeo-environmental significance of the calcite crusts in the Untersee Oasis:

  • The calcite crusts are located within the vicinity of residual snow patches.

  • The calcite crusts have some of the highest δ18O (up to +17.4‰ VSMOW) and δ13C (up to +14.6‰ VPDB) compositions of any carbonate deposits in terrestrial polar environments. Their δ18O and δ13C values are substantially enriched with respect to the δ18O (average of -33.5 ± 2.4‰) and δ13CDIC of snow meltwater (inferred to be in the -2.8 to -1.4‰ range).

  • The formation of the calcite crusts is ascribed to the evaporation of residual snow meltwater and the low relative humidity and strong winds in the Untersee Oasis, which would favour a kinetic isotope effect during the precipitation of the evaporative calcite crusts.

  • The 14C age distribution of the calcite crusts (1589 cal yr bp to modern) provides a minimum age for ice retreat and drainage of the palaeo-lake in Aurkjosen Cirque. However, in this polar desert environment, the calcite crusts require suitable climate conditions for their formation (sufficient snow accumulation and temperatures warm enough to generate meltwater), and their age distribution corresponds to the late Holocene warm-wet climate period.

Acknowledgements

Fieldwork was made possible by the logistical support of the Antarctic Logistics Centre International (Cape Town, South Africa) and the Arctic and Antarctic Research Institute/Russian Antarctic Expedition. We are grateful to Colonel (I.L.) J.N. Pritzker, IL ARNG (retired), Lorne Trottier and field team members for their support during the 2015 and 2017 expeditions. We thank the two reviewers for their constructive comments.

Financial support

This work was supported by contributions from TAWANI Foundation and Trottier Family Foundation to DA, NASA's Exobiology program to DA and an NSERC Discovery Grant to DL.

Competing interests

The authors declare none.

Author contributions

DL, BF and DA designed the project and collected the samples. DL, MC, BF, PS and DA contributed to sample and data analyses and to writing/editing the manuscript.

References

Andersen, D.T., McKay, C.P. & Lagun, V. 2015. Climate conditions at perennially ice-covered Lake Untersee, East Antarctica. Journal of Applied Meteorology and Climatology, 54, 10.1175/JAMC-D-14-0251.1.CrossRefGoogle Scholar
Blake, W. 1999. Glaciated landscapes along Smith Sound, Ellesmere Island, Canada and Greenland. Annals of Glaciology, 28, 10.3189/172756499781821814.CrossRefGoogle Scholar
Blake, W. 2005. Holocene carbonate precipitates on Precambrian bedrock in the High Arctic: age and potential for palaeoclimatic information. Geografiska Annaler: Series A, Physical Geography, 87, 10.1111/j.0435-3676.2005.00251.x.Google Scholar
Bormann, P., Bankwitz, P., Bankwitz, E., Damn, V., Hurtig, E., Kampf, H., et al. 1986. Structure and development of the passive continental margin across the Princess Astrid Coast, East Antarctica. Journal of Geodynamics, 373, 347373.CrossRefGoogle Scholar
Bronk Ramsey, C. 2009. Bayesian analysis of radiocarbon dates. Radiocarbon, 51, 10.2458/azu_js_rc.v51i1.3494.CrossRefGoogle Scholar
Buizert, C., Fudge, T.J., Roberts, W.H.G., Steig, E.J., Sherriff-Tadano, S., Ritz, C., et al. 2021. Antarctic surface temperature and elevation during the Last Glacial Maximum. Science, 372, 10.1126/SCIENCE.ABD2897/SUPPL_FILE/ABD2897_DATAS1.XLSX.CrossRefGoogle Scholar
Campbell, I.B. & Claridge, G.G.C. 2009. Antarctic permafrost soils. In Margesin, R., ed., Permafrost soils. Berlin: Springer, 10.1007/978-3-540-69371-0_2.Google Scholar
Campbell, I.B., Claridge, G.G.C., Campbell, D.I. & Balks, M.R. 2013. The soil environment of the McMurdo Dry Valleys, Antarctica. Antarctic Research Series, 72, 10.1029/AR072P0297.CrossRefGoogle Scholar
Claridge, G.G. & Campbell, I.B. 1977. The salts in Antarctic soils, their distribution and relationship to soil processes. Soil Science, 123, 10.1097/00010694-197706000-00006.CrossRefGoogle Scholar
Clark, I.D. & Lauriol, B. 1992. Kinetic enrichment of stable isotopes in cryogenic calcites. Chemical Geology, 102, 10.1016/0009-2541(92)90157-Z.CrossRefGoogle Scholar
Clark, I.D. & Lauriol, B. 1997. Aufeis of the Firth River Basin, northern Yukon, Canada: insights into permafrost hydrogeology and karst. Arctic and Alpine Research, 29, 10.2307/1552053.CrossRefGoogle Scholar
Clark, I.D., Lauriol, B., Marschner, M., Sabourin, N., Chauret, Y. & Desrochers, A. 2004. Endostromatolites from permafrost karst, Yukon, Canada: paleoclimatic proxies for the Holocene hypsithermal. Canadian Journal of Earth Sciences, 41, 10.1139/e04-014.CrossRefGoogle Scholar
Coplen, T.B., Hopple, J.A., Böhlke, J.K., Peiser, H.S., Rieder, S.E., Krouse, H.R., et al. 2002. Compilation of minimum and maximum isotope ratios of selected elements in naturally occurring terrestrial materials and reagents. Retrieved from https://pubs.usgs.gov/wri/wri014222/pdf/wri01-4222.pdfGoogle Scholar
Courty, M.A., Marlin, C., Dever, L., Tremblay, P. & Vachier, P. 1994. The properties, genesis and environmental significance of calcitic pendents from the High Arctic (Spitsbergen). Geoderma, 61, 10.1016/0016-7061(94)90012-4.CrossRefGoogle Scholar
Crann, C.A., Murseli, S., St-Jean, G., Zhao, X., Clark, I.D. & Kieser, W.E. 2017. First status report on radiocarbon sample preparation techniques at the A.E. Lalonde AMS Laboratory (Ottawa, Canada). Radiocarbon, 59, 10.1017/RDC.2016.55.CrossRefGoogle Scholar
Criss, R.E. 1999. Principles of stable isotope distribution. Oxford: Oxford University Press, 10.1093/oso/9780195117752.001.0001.CrossRefGoogle Scholar
Diaz, M.A., Li, J., Michalski, G., Darrah, T.H., Adams, B.J., Wall, D.H., et al. 2020. Stable isotopes of nitrate, sulfate, and carbonate in soils from the Transantarctic Mountains, Antarctica: a record of atmospheric deposition and chemical weathering. Frontiers in Earth Science, 8, 10.3389/FEART.2020.00341/BIBTEX.CrossRefGoogle Scholar
Eggleston, S., Schmitt, J., Bereiter, B., Schneider, R. & Fischer, H. 2016. Evolution of the stable carbon isotope composition of atmospheric CO2 over the last glacial cycle. Paleoceanography, 31, 10.1002/2015PA002874.CrossRefGoogle Scholar
Faucher, B., Lacelle, D., Fisher, D.A., Andersen, D.T. & McKay, C.P. 2019. Energy and water mass balance of Lake Untersee and its perennial ice cover, East Antarctica. Antarctic Science, 31, 10.1017/S0954102019000270.CrossRefGoogle Scholar
Faucher, B., Lacelle, D., Marsh, N.B., Fisher, D.A. & Andersen, D.T. 2021. Ice-covered ponds in the Untersee Oasis (East Antarctica): distribution, chemical composition, and trajectory under a warming climate. Arctic, Antarctic, and Alpine Research, 53, 10.1080/15230430.2021.2000566.CrossRefGoogle Scholar
Fisher, D.A., Lacelle, D., Pollard, W. & Faucher, B. 2020. A model for stable isotopes of residual liquid water and ground ice in permafrost soils using arbitrary water chemistries and soil-specific empirical residual water functions. Permafrost and Periglacial Processes, 32, 10.1002/ppp.2079.Google Scholar
Fisher, D.A., Lacelle, D., Pollard, W., Davila, A.F. & McKay, C.P. 2016. Ground surface temperature and humidity, ground temperature cycles and the ice table depths in University Valley, McMurdo Dry Valleys of Antarctica. Journal of Geophysical Research– - Earth Surface, 121, 10.1002/2016JF004054.CrossRefGoogle Scholar
Hanshaw, B.B. & Hallet, B. 1978. Oxygen isotope composition of subglacially precipitated calcite: possible paleoclimatic implications. Science, 200, 10.1126/SCIENCE.200.4347.1267.CrossRefGoogle Scholar
Hillaire-Marcel, C., Soucy, J.M. & Cailleux, A. 1979. Isotopic analysis of sub-glacial concretions of the Laurentide ice sheet and the oxygen-18 content of the ice. Canadian Journal of Earth Sciences, 16, 10.1139/E79-132/ASSET/IMAGES/E79-132C4.GIF.Google Scholar
Hiller, A., Hermichen, W.D. & Wand, U. 1995. Radiocarbon-dated subfossil stomach oil deposits from petrel nesting sites: novel paleoenvironmental records from Continental Antarctica. Radiocarbon, 37, 10.1017/S0033822200030617.CrossRefGoogle Scholar
Hiller, A., Wand, U., Kämpf, H. & Stackebrandt, W. 1988. Occupation of the Antarctic continent by petrels during the past 35 000 years: inferences from a 14C study of stomach oil deposits. Polar Biology, 9, 10.1007/BF00442032/METRICS.CrossRefGoogle Scholar
Hoffman, M.J., Fountain, A.G. & Liston, G.E. 2008. Surface energy balance and melt thresholds over 11 years at Taylor Glacier, Antarctica. Journal of Geophysical Research, 113, 10.1029/2008JF001029.CrossRefGoogle Scholar
Isaksson, E., Karlén, W., Gundestrup, N., Mayewski, P., Whitlow, S. & Twickler, M. 1996. A century of accumulation and temperature changes in Dronning Maud Land, Antarctica. Journal of Geophysical Research– - Atmospheres, 101, 10.1029/95JD03232.CrossRefGoogle Scholar
Kampf, H. & Stakebrandt, W. 1985. Geological investigation in the Eliseev anorthosite massif, central Dronning Maud Land, East Antarctica. Zeitschrift fur Geologische Wissenschaften, 13, 3260.Google Scholar
Kim, S.T. O'Neil, J.R. 1997. Equilibrium and nonequilibrium oxygen isotope effects in synthetic carbonates. Geochimica et Cosmochimica Acta, 61, 10.1016/S0016-7037(97)00169-5.CrossRefGoogle Scholar
Kim, S.T., Coplen, T.B. & Horita, J. 2015. Normalization of stable isotope data for carbonate minerals: implementation of IUPAC guidelines. Geochimica et Cosmochimica Acta, 158, 10.1016/J.GCA.2015.02.011.CrossRefGoogle Scholar
Knauth, L.P., Brilli, M. & Klonowski, S. 2003. Isotope geochemistry of caliche developed on basalt. Geochimica et Cosmochimica Acta, 67, 10.1016/S0016-7037(02)01051-7.CrossRefGoogle Scholar
Lacelle, D. 2007. Environmental setting, (micro)morphologies and stable C-O isotope composition of cold climate carbonate precipitat–s - a review and evaluation of their potential as paleoclimatic proxies. Quaternary Science Reviews, 26, 10.1016/j.quascirev.2007.03.011.CrossRefGoogle Scholar
Lacelle, D., Lauriol, B. & Clark, I.D. 2006. Effect of chemical composition of water on the oxygen-18 and carbon-13 signature preserved in cryogenic carbonates, Arctic Canada: implications in paleoclimatic studies. Chemical Geology, 234, 10.1016/j.chemgeo.2006.04.001.CrossRefGoogle Scholar
Lacelle, D., Lauriol, B. & Clark, I.D. 2007. Origin, age, and paleoenvironmental significance of carbonate precipitates from a granitic environment, Akshayuk Pass, southern Baffin Island, Canada. Canadian Journal of Earth Sciences, 44, 10.1139/e06-088.CrossRefGoogle Scholar
Lapalme, C., Lacelle, D., Pollard, W., Fisher, D., Davila, A. & McKay, C.P. 2017. Distribution and origin of ground ice in University Valley, McMurdo Dry Valleys, Antarctica. Antarctic Science, 29, 10.1017/S0954102016000572.CrossRefGoogle Scholar
Lauriol, B. & Clark, I. 1999. Fissure calcretes in the arctic: a paleohydrologic indicator. Applied Geochemistry, 14, 10.1016/S0883-2927(98)00090-0.CrossRefGoogle Scholar
Lauriol, B., Cinq Mars, J. & Clark, I.D. 1991. Localisation, Genese et Fonte de Quelques Naleds du Nord du Yukon (Canada). Permafrost and Periglacial Processes, 2, 10.1002/ppp.3430020306.CrossRefGoogle Scholar
Levitan, M.A., Kononkova, N.N., Luksha, V.L. & Roshchina, I.A. 2012. Holocene lithificates at the slopes of the Untersee mountain valley, East Antarctica. Geochemistry International, 50, 10.1134/S0016702912040040/METRICS.CrossRefGoogle Scholar
Lyons, B., Foley, K., Carey, A., Diaz, M., Bowen, G. & Cerling, T. 2020. The isotopic geochemistry of CaCO3 encrustations in Taylor Valley, Antarctica: implications for their origin. Acta Geographica Slovenica, 60, 10.3986/AGS.7233.CrossRefGoogle Scholar
Marinova, M.M., McKay, C.P., Heldmann, J.L., Goordial, J., Lacelle, D., Pollard, W.H. & Davila, A.F. 2022. Climate and energy balance of the ground in University Valley, Antarctica. Antarctic Science, 34, 10.1017/S0954102022000025.CrossRefGoogle Scholar
Marsh, N.B., Lacelle, D., Faucher, B., Cotroneo, S., Jasperse, L., Clark, I.D. & Andersen, D.T. 2020. Sources of solutes and carbon cycling in perennially ice-covered Lake Untersee, Antarctica. Scientific Reports, 10, 10.1038/s41598-020-69116-6.CrossRefGoogle Scholar
McKay, C.P., Balaban, E., Abrahams, S. & Lewis, N. 2019. Dry permafrost over ice-cemented ground at Elephant Head, Ellsworth Land, Antarctica. Antarctic Science, 31, 10.1017/S0954102019000269.CrossRefGoogle Scholar
Nakai, N., Wada, H., Kiyosu, Y. & Takimoto, M. 1975. Stable isotope of water and studies on the origin and geological history salts in the Lake Vanda area, Antarctica. Geochemical Journal, 9, 10.2343/GEOCHEMJ.9.7.CrossRefGoogle Scholar
Omelon, C.R., Pollard, W.H. & Marion, G.M. 2001. Seasonal formation of ikaite (CaCO3.6H2O) in saline spring discharge at Expedition Fiord, Canadian High Arctic: assessing conditional constraints for natural crystal growth. Geochimica et Cosmochimica Acta, 65, 10.1016/S0016-7037(00)00620-7.CrossRefGoogle Scholar
Paech, H.-J. & Stackebrandt, W. 1995. Geology. In Bormann, P. & Fritzsche, D., eds, The Schirmacher Oasis, Queen Maud Land, East Antarctica and its surroundings. Cambridge: Cambridge University Press, 59159.Google Scholar
Reimer, P.J., Austin, W.E.N., Bard, E., Bayliss, A., Blackwell, P.G., Bronk Ramsey, C., et al. 2020. The IntCal20 Northern Hemisphere radiocarbon age calibration curve (0–55 cal kbp). Radiocarbon, 62, 10.1017/RDC.2020.41.CrossRefGoogle Scholar
Schmitt, J., Schneider, R., Elsig, J., Leuenberger, D., Lourantou, A., Chappellaz, J., et al. 2012. Carbon isotope constraints on the deglacial CO2 rise from ice cores. Science, 336, 10.1126/science.1217161.CrossRefGoogle Scholar
Schwab, M.J. 1998. Reconstruction of the late Quaternary climatic and environmental history of the Schirmacher Oasis and the Wohlthat Massif (East Antarctica). Bremerhaven: Alfred-Wegener-Institut für Polar- und Meeresforschung, 128 pp.Google Scholar
Shamilishvili, G., Abakumov, E.V. & Andersen, D. 2020. Biogenic-abiogenic interactions and soil formation in extreme conditions of Untersee Oasis, surroundings of Lake Untersee, central Queen Maud Land, East Antarctica. In Frank-Kamanetskaya, O., Vlasov, D.Yu., Panova, E.G. & Lessovaia, S.N., eds, Processes and phenomena on the boundary between biogenic and abiogenic nature. Berlin: Springer, 10.1007/978-3-030-21614-6_25/TABLES/4.Google Scholar
Sofer, Z. & Gat, J.R. 1975. The isotope composition of evaporating brines: effect of the isotopic activity ratio in saline solutions. Earth and Planetary Science Letters, 26, 10.1016/0012-821X(75)90085-0.CrossRefGoogle Scholar
Souchez, R.A. & Lemmens, M. 1985. Subglacial carbonate deposition: an isotopic study of a present-day case. Palaeogeography, Palaeoclimatology, Palaeoecology, 51, 10.1016/0031-0182(85)90093-8.CrossRefGoogle Scholar
Stuiver, M. & Polach, H.A. 1977. Reporting of 14C data. Radiocarbon, 19, 355363.CrossRefGoogle Scholar
Van den Broke, M., Jan van de Berg, W., van Meijgaard, E. & Reijmer, C. 2006. Identification of Antarctic ablation areas using a regional atmospheric climate model. Journal of Geophysical Research - Atmospheres, 111, 10.1029/2006JD007127.Google Scholar
Waragai, T. 2005. Holocene calcrete crust deposits on the moraine of Batura Glacier, northern Pakistan. Island Arc, 14, 10.1111/J.1440-1738.2005.00492.X.CrossRefGoogle Scholar
Figure 0

Figure 1. Map of the Untersee Oasis in Queen Maud Land, East Antarctica. (Top) The locations of sampling sites of the calcite crusts and inferred glacial limits (from Schwab 1998). Ages of the calcite crusts are provided as cal yr bp. (Bottom) Recessional glacial limits in the Untersee Oasis and palaeo-lake basin in Aurkjosen Cirque (from Schwab 1998).

Figure 1

Figure 2. Field photographs of calcite crusts in the Untersee Oasis, East Antarctica. The calcite crusts are up to 1 mm thick.

Figure 2

Figure 3. Raman spectra of five carbonate crusts in the Untersee Oasis, East Antarctica. The intense sharp Raman band at 1088 cm-1 is indicative of the mineral calcite.

Figure 3

Figure 4. δ13C and δ18O compositions of the calcite crusts in the Untersee Oasis, East Antarctica. a. δ13C and δ18O compositions of the calcite crusts in the Untersee Oasis compared to others in Antarctica: lithificates along palaeo-shorelines of Lake Untersee (Levitan et al.2012), evaporative calcite crusts surrounding Lake Vanda (Nakai et al.1975) and paedogenic calcite in soils surrounding Shackleton Glacier (Diaz et al.2020). b. Comparison of δ13C and δ18O compositions of the calcite crusts in the Untersee Oasis (red box) with other terrestrial polar environments (from Lacelle 2007). VPDB = Vienna Pee Dee Belemnite; VSMOW = Vienna Standard Mean Ocean Water.

Figure 4

Table I. δ13C and δ18O measurements of calcite crusts in the Untersee Oasis, East Antarctica. The δ18O composition of the carbonate crusts were converted to the Vienna Standard Mean Ocean Water (VSMOW) scale according to: δ18OVSMOW = 1.0309 δ18OVPDB + 30.92 (Coplen et al.2002).

Figure 5

Table II. Radiocarbon results from calcite crusts collected in the Untersee Oasis, East Antarctica. The 14C ages were calibrated using OxCal and the IntCal20 calibration curve (Bronk Ramsey 2009, Reimer et al.2020).

Figure 6

Figure 5. Evolution of δ18O of evaporating water (snow meltwater chemistry) for a range of relative humidities (RHs) using the Criss (1999) and Sofer & Gat (1975) models. The initial δ18O composition of snow was set at -33‰ (average of the measurements) and that of the atmospheric moisture above the evaporating surface was set at -44‰ (equilibrium with meltwater) and -37‰, and the RH at the ground surface ranged from 50 to 95% (values near 50% represent atmospheric humidity, those > 70% probably represent those at the ground surface; i.e. Fisher et al.2016). VSMOW = Vienna Standard Mean Ocean Water.

Figure 7

Figure 6. Comparison of radiocarbon age distributions obtained from evaporative calcite crusts in the Untersee Oasis, East Antarctica, with Holocene reconstructed temperatures and ice accumulation rates (Buizert et al.2021).