Articles | Volume 19, issue 8
Research article
04 Aug 2023
Research article |  | 04 Aug 2023

Cryogenic cave minerals recorded the 1889 CE melt event in northeastern Greenland

Anika Donner, Paul Töchterle, Christoph Spötl, Irka Hajdas, Xianglei Li, R. Lawrence Edwards, and Gina E. Moseley

The investigation of cryogenic cave minerals (CCMs) has developed in recent decades to be a particularly valuable proxy for palaeo-permafrost reconstruction. Due to difficulties, however, in obtaining reliable chronologies with the so-called “fine” form of these minerals, such studies have thus far utilised the “coarse” form. In this study, we successfully investigate the northernmost-known deposit of fine-grained CCMs, which are situated in Cove Cave (Greenlandic translation: Eqik Qaarusussuaq), a low-elevation permafrost cave in northeastern Greenland (80 N). The Cove Cave CCMs display a complex mineralogy that consists of fine-grained cryogenic cave carbonates and sulfate minerals (gypsum, eugsterite, mirabilite, and löweite). Until now, previous attempts to date fine-grained CCMs have been unsuccessful; however, here we demonstrate that precise dating is possible with both isochron-based 230Th / U dating and 14C dating if the dead carbon fraction is reliably known.

The dating result (65±17 a BP; 1885±17 CE) shows that the Cove Cave CCMs formed during the late Little Ice Age, a time interval characterised by cold temperatures and abundant permafrost in northeastern Greenland, making water infiltration into Cove Cave dependent on the water amount and latent heat. We relate the CCM formation to a combination of black carbon deposition and anomalously high temperatures, which led to widespread melting over large areas of the Greenland ice sheet in the course of a few days. We propose that the anomalous weather conditions of 1889 CE also affected northeastern Greenland, where the enhanced melting of a local ice cap resulted in water entering the cave and rapidly freezing. While calcite and gypsum likely precipitated concurrently with freezing, the origin of the other sulfate minerals might not be purely cryogenic but could be linked to the subsequent sublimation of this ice accumulation in a very dry cave environment.

1 Introduction

In recent decades, cryogenic cave carbonates (CCCs), a type of speleothem associated with the formation of cave ice, have become a valuable tool for tracking evidence of past permafrost presence, particularly in the mid-latitudes in either low-elevation temperate locations (e.g. central Europe; Richter et al., 2018; Žák et al., 2012) or high-elevation periglacial environments (Bartolomé et al., 2015; Luetscher et al., 2013; Spötl et al., 2021; Spötl and Cheng, 2014). In contrast, investigations into CCCs from high-latitude caves are rare, with the exception of a few studies from northern Yukon, Canada (Clark and Lauriol, 1992; Lauriol et al., 1988; Lauriol and Clark, 1993).

The presently accepted mechanism for CCC formation is precipitation from freezing karst water (Žák et al., 2012). There are two types of CCCs which are commonly differentiated by grain size (ca. 1 mm; Žák et al., 2018), resulting in the terms coarse- and fine-grained CCCs (CCCcoarse and CCCfine), despite the fact that the isotopic composition is the main diagnostic feature (a result of open- vs. closed-system precipitation; Luetscher et al., 2013). CCCcoarse are typically found in micro-climatically stable environments of cave interiors, where they precipitate in slowly freezing pools of water carved into cave ice deposits by drip water, making them a useful proxy for the reconstruction of palaeo-permafrost (i.e. CCCcoarse= palaeo-permafrost) and negative cave temperatures close to 0 C (Žák et al., 2018). In contrast, CCCfine precipitate from a rapidly freezing film of water on top of cave ice and are commonly found in well-ventilated caves and/or near cave entrances and thus may be related to local thermal anomalies, causing them to be unsuitable for permafrost reconstruction (Žák et al., 2012). In addition to carbonates, other minerals are also known to form cryogenically in caves, e.g. gypsum and other sulfate minerals, collectively referred to as cryogenic cave minerals (CCMs; Žák et al., 2018). In this paper, we follow the approach of Žák et al. (2018) and use the term fine-grained CCMs to describe our samples, at times distinguishing the carbonate and sulfate fractions, while also comparing them to the CCM subtypes, CCCcoarse and CCCfine, which previous studies mostly focussed on.

CCCcoarse can be radiometrically dated using standard 230Th / U techniques (e.g. Koltai et al., 2021; Luetscher et al., 2013; Spötl et al., 2021; Žák et al., 2012), although the accuracy of 230Th / U ages using standard evaluation procedures has recently been called into question by showing that the conventionally used correction factors for detrital Th contamination are not universally applicable to CCCcoarse samples (Töchterle et al., 2022). Dating of CCCfine, on the other hand, has proven to be a difficult task. 14C dating of CCCfine suffers from large uncertainties associated with estimating the initial radiocarbon activity and the dead carbon fraction (DCF; Lauriol and Clark, 1993), whereas 230Th / U dating has been hindered by poor age precision due to very high detrital thorium contamination, which is particularly challenging in young (i.e. Late Holocene) samples (Spötl, 2008; Spötl and Cheng, 2014). Due to these unsuccessful dating efforts, the potential of CCCfine as a palaeoclimate archive has yet to be realised.

In this study, we investigate fine-grained CCMs from a low-elevation permafrost cave (Cove Cave, unofficial name; Eqik Qaarusussuaq in Greenlandic) in northeastern Greenland. Cove Cave is currently the northernmost-known cave containing cryogenic mineral deposits and is located within a highly climatically sensitive region adjacent to the Greenland ice sheet (Bintanja and Krikken, 2016; Bintanja and Selten, 2014; Shepherd, 2016). This paper aims to (i) extend the existing knowledge on morphology, mineralogy, and stable isotopic composition of fine-grained high-latitude CCMs; (ii) constrain the age of fine-grained CCMs (and CCCfine) using 230Th / U and 14C dating methods; and (iii) ascertain the circumstances of CCM formation at this location.

2 Study site

Cove Cave (80.25 N, 21.93 W) is located in a small tributary valley to a steep-sided canyon on Crown Prince Christian Land in northeastern Greenland (Fig. 1; Moseley et al., 2020). The Silurian limestones and dolostones in this area host numerous solutional caves (Smith and Rasmussen, 2020), some of which were first discovered in 1960 (Davies and Krinsley, 1960), followed by more cave discoveries during three subsequent caving expeditions (Loubière, 1987; Moseley et al., 2020). The area is located ca. 35 km from the coast (eastwards) and ca. 60 km from the margin of the Greenland ice sheet (southwestwards). It is characterised by an arid climate (ca. 200 mm a−1; Schuster et al., 2021), permanently frozen ground, sparse soil and vegetation cover, and permafrost landforms (Moseley et al., 2020). Under contemporary climatic conditions, common speleothems (i.e. stalactites, stalagmites, and flowstones) cannot form in the caves of the area due to the lack of water infiltration. Weather stations or long-term weather observations are absent in the area; however, mean annual air temperatures at the nearest weather stations on the ice sheet (KPC_L PROMICE weather station, July 2008–December 2021; ca. 70 km away; Fausto et al., 2019; van As et al., 2011) and near the coast (Station Nord, 1991–2020; ca. 180 km away; Danmarks Meteorologiske Institut, 2022) are -13. and −15.1C, respectively. Compared to these two stations, a more continental climate with warmer summer temperatures is expected for the study site (Donner et al., 2020). In 1960 and 1983, a small ice cap was present on the plateau in close proximity to Cove Cave, but this has since melted (Fig. 1, light blue polygon; Davies and Krinsley, 1960; Loubière, 1987; Moseley et al., 2020, 2021). Elsewhere in the area, local plateaus are now ice-free in summer (Moseley et al., 2021), though there is ample geomorphological evidence that ice caps were present in the past (Fig. 1, blue lines; Sole et al., 2020), which might have covered the plateau above Cove Cave.

Figure 1Location of Cove Cave in a tributary valley of the larger Grottedal in northeastern Greenland (map insert from Geofabrik and OpenStreetMap Contributors, 2018). The cave is located close to the margins of former ice caps. The most recent ice cap was mapped in 1960 (light blue polygon; Davies and Krinsley, 1960), still existed in 1983 (Loubière, 1987), and disappeared before 2019 (Moseley et al., 2020). Geomorphological evidence of ice margins indicates the former presence of additional ice caps, which have not been dated (blue lines; Sole et al., 2020).

Cove Cave consists of a 103 m long, gently dipping phreatic passage, currently making it the longest cave being explored in Greenland (Fig. 2; Moseley et al., 2020). The cave hosts both CCMs and inactive flowstones (i.e. without active water supply), with the latter indicating at least one phase of warmer and wetter climate in the past. The cave entrance is located 660 m a.s.l. (above sea level), with a rock overburden of ca. 25 m. Inside the entrance area, an ice pond, ice stalagmites, and hoar frost were observed during a visit in summer (Moseley et al., 2020). Beyond the entrance area, the cave is devoid of ice. Flowstone drapes the walls of a vadose canyon, and broken angular flowstone blocks, likely shattered by freeze–thaw processes, are scattered on the floor (Moseley et al., 2020). Accumulations of CCMs lie on the shattered flowstone blocks in an area ca. 3×1 m. At this location, an air temperature of −14.7C was measured in July 2019 (Fig. 2), while the outside air temperatures reached up to 18 C (Donner et al., 2020; Moseley et al., 2020). Measurements of relative humidity only exist for other parts of the cave, where temperatures were above −10C, reaching values as low as 39 % at floor level (Barton et al., 2020). It can therefore be inferred that the relative humidity in colder parts of the cave must also be low. Cold air has less capacity to hold moisture, and incoming air loses much of its humidity by resublimation, creating hoar frost on the walls close to the cave entrance (Barton et al., 2020; Lauriol et al., 1988). Deeper inside the cave, flowstone deposits block a ca. 5 m deep vadose slot, where cave air temperatures reached −17.1C at the bottom (Barton et al., 2020; Moseley et al., 2020). The cause of these low temperatures is likely density-driven flow of very cold winter air into this descending single-entrance cave. During summer, a stable air stratification inside this cold trap prevents the advection of warm outside air into the cave (Barton et al., 2020).

Figure 2Longitudinal section of Cove Cave. Ice formations such as an ice pond, ice stalagmites, and hoar frost are only present in the first third of the cave. Beyond that, the air temperatures are much lower. CCMs were found ca. 65 m behind the entrance (adapted from Moseley et al., 2020).

3 Materials and methods

3.1 Sampling and sample preparation

Sampling took place in four different spots within the 3×1 m accumulation of CCMs in July 2019. Four samples (KC19CCC-1, KC19CCC-2, KC19CCC-3, and KC19CCC-4) were picked up with a knife, wrapped in aluminium foil, and stored in plastic containers. There were no obvious visual differences between the samples during collection. During sampling, disturbance to the overall appearance of the CCM accumulation and Cove Cave was kept to a minimum.

In the lab, the samples were homogenised and transferred to glass vials in a laminar flow hood. For all subsequent analyses, macroscopic contaminants (e.g. pieces of bedrock or insect remains) were removed under a binocular microscope. To aid the calculation of the DCF, four artificial mixtures ranging from 100 wt % CCM fraction to 100 wt % non-CCM fraction were produced from KC19CCC-4 by splitting the sample into the different mineral fractions and mixing the relative amounts. Overall, analysis options and the reproducibility of analyses were limited by the available sample amount.

3.2 Mineralogy and crystal morphology

A Bruker D8 DISCOVER X-ray diffractometer (XRD), using Bragg–Brentano geometry equipped with a Cu target and a LYNXEYE detector, was used to analyse the mineralogical composition of the samples (Cu–Kα radiation 1.5406 Å; 2θ range 8–45). Furthermore, a HORIBA Jobin Yvon LabRam HR800 spectrometer, excited by a frequency-doubled Nd:YAG laser (100 mW, 532 nm), was used for the in situ determination of the mineralogy of crystals and crystal aggregates at a resolution of ca. 5 µm.

A KEYENCE VHX-6000 digital microscope was used to examine the morphology of the CCMs. To study the fine fraction, a field-emission scanning electron microscope (SEM) operating at 10 kV accelerating voltage was utilised (DSM 982 GEMINI, ZEISS Group).

3.3 Stable isotope analysis

In total, 14 aliquots of 0.15 to 0.7 mg of the carbonate fraction were taken from the four samples, and their carbon and oxygen stable isotopic composition was analysed using a Thermo Fisher DELTA V Plus isotope ratio mass spectrometer coupled with a GasBench II (Spötl and Vennemann, 2003), yielding a long-term precision of ±0.08 ‰ for δ18O and ±0.06 ‰ (1σ) for δ13C (Spötl, 2011). For comparison, carbon and oxygen isotope data from inactive common speleothems in the study area (a few square kilometres) were included. The common speleothems were collected in several caves, including Cove Cave. All results were calibrated against international standards and reported relative to the Vienna Pee Dee Belemnite (VPDB) standard.

3.4 Radioisotope dating

3.4.1230Th / U disequilibrium dating and isochron construction

For 230Th / U disequilibrium dating of the carbonate fraction, 20 mg aliquots from each of the four samples were picked in a laminar flow hood. Dating, including chemical preparation and multi-collector inductively coupled mass spectrometry, was carried out at the Trace Metal Isotope Geochemistry Laboratory at the University of Minnesota, following Edwards et al. (1987) and Shen et al. (2012). Ages are reported in years before 1950 CE (a BP), with 2σ uncertainties. Additionally, δ234U values of inactive common speleothems from Cove Cave were included in this study.

An isochron using maximum likelihood regression (Ludwig and Titterington, 1994) was constructed in IsoplotR (Vermeesch, 2018), using the 238U-normalised activities of 234U, 230Th, and 232Th. Additionally, the initial 230Th/232Th activity ratio was ascertained and applied to the detrital Th correction of the individual ages as a reliability test of the derived isochron values in order to establish whether there was one source or multiple sources of detrital 230Th (Dorale et al., 2004). Ages are reported in years before present, with 2σ measurement uncertainties.

3.4.2 Radiocarbon dating and calculation of DCF

Eight aliquots of 16–31 mg were taken from the samples, including artificial mixtures. All eight aliquots were analysed with a MICADAS accelerator mass spectrometer (AMS) in the Laboratory of Ion Beam Physics at ETH Zürich (Synal et al., 2007). Radiocarbon values are reported as conventional radiocarbon ages before present ±1σ (Stuiver and Polach, 1977) and as fraction modern ±1σ (F14C; Reimer et al., 2004, 2020), relative to 95 % of the secondary standard HOx2 (Hajdas et al., 2021).

The DCF calculation is based on the equation given by Genty and Massault (1997), using the measured F14C of the sample (F14Csample) in relation to the atmospheric F14C at the time of formation (F14Catm), which was obtained by the corrected 230Th / U age and the IntCal20 calibration curve (Reimer et al., 2020; Eq. 1):

(1) DCF = 1 - F 14 C sample F 14 C atm × 100 % .

The radiocarbon ages were then corrected for the corresponding DCFs, resulting in a radiocarbon age ±2 SE (standard error) for 0 % DCF, calculated by linear approximation and calibrated against IntCal20 (cal BP; Reimer et al., 2020), using OxCal 4.4 (Bronk Ramsey, 2009).

4 Results

4.1 Mineralogical composition and particle morphology

All four samples show a complex mineralogy and a range of particle morphologies that consist of crystal aggregates and single crystals of variable size. The aggregates range from 50 to 400 µm in diameter, rarely exceeding 500 µm, while individual crystals range from 1 µm to a few tens of micrometres. Brownish submicrometric crystals dominate sample KC19CCC-1, which contains larger translucent and opaque white and brownish crystals, as well as crystal aggregates (Fig. 3a). Sample KC19CCC-2 (Fig. 3b) is homogeneously white, with opaque and translucent crystals/crystal aggregates. The other two samples (KC19CCC-3 and KC19CCC-4; Fig. 3c) are made up of larger white and brownish crystals/crystal aggregates and brownish submicrometric crystals, resulting in a speckled appearance.

Figure 3Optical (a–c) and SEM (d) images of the studied cryogenic samples. (a) Brownish, very fine crystals (quartz, potassium feldspar, and traces of dolomite) intermixed with brownish crystal aggregates (calcite) and white/translucent single crystals made up of calcite, dolomite, and gypsum (sample KC19CCC-1). (b) White calcite, gypsum, eugsterite, mirabilite, and löweite (sample KC19CCC-2). (c) Mixed white (calcite and sulfate minerals) and brownish (calcite) particles, with traces of quartz (samples KC19CCC-3 and KC19CCC-4). (d) Intermixed dumbbell-shaped calcite and acicular sulfate minerals, which are present in all samples.


XRD results (see Appendix A) show that the samples are mostly mixtures of calcite and the sulfate minerals gypsum (CaSO42H2O), eugsterite (Na4Ca(SO4)32H2O), löweite (Na12Mg7(SO4)1315H2O), and mirabilite (Na2SO410H2O). Sample KC19CCC-1 contains additional phases such as quartz, dolomite, and potassium feldspar, while quartz is also present in samples KC19CCC-3 and KC19CCC-4.

The CCMs contain a variety of crystal morphologies, comparable to the morphologies of CCCcoarse (Žák et al., 2018), although having much smaller sizes. Predominant morphologies are spherulitic and acicular and also aggregates of elongated prismatic crystals. Furthermore, rhombic and dumbbell-shaped morphologies were observed. Optical microscopy combined with micro-Raman spectroscopy showed that calcite is mostly represented as translucent and opaque spherulitic crystals and crystal aggregates, such as chains of different colours, whereas sulfate minerals mostly consist of white fibrous, acicular, and prismatic crystals. Raman spectroscopy and SEM imagery also showed that microcrystalline calcite and sulfate minerals are often intermixed (Fig. 3d).

4.2 Stable isotopes

The stable isotopic composition of the 14 aliquots ranges from −24.0 ‰ to −16.0 ‰ for δ18O and from 7.0 ‰ to 11.4 ‰ for δ13C (Fig. 4). The δ13C values overlap with those of CCCfine from mid-latitude caves (Žák et al., 2018). The δ18O values, however, are significantly lower than values from mid-latitude CCCfine and fall at the lower end of δ18O values of mid-latitude CCCcoarse (Žák et al., 2018).

Figure 4Stable isotope composition of CCMs and common speleothems from northeastern Greenland (this study) compared to CCCfine from the Canadian Arctic Circle (Clark and Lauriol, 1992) and CCCfine and common speleothems from mid-latitude caves (Žák et al., 2018). The composition of CCCcoarse from mid-latitude caves is also shown for comparison (Žák et al., 2018). Note the latitudinal shift in δ18O (blue arrows) that can be observed in CCCfine, CCMs, and common speleothems between the mid- and high latitudes.


4.3 Radioisotope dating

Of the four 230Th / U analyses, three yielded results (KC19CCC-2, KC19CCC-3, and KC19CCC-4), while KC19CCC-1 contained too much detrital material to be analysed. Table 1 shows that the 238U content is high (1113±2 to 1278±2 ng g−1), as is the concentration of detrital thorium (indicated by the low 230Th/232Th activity ratio between 1.30±0.04 and 3.89±0.22). Uncorrected ages show a large range from 143±11 to 7335±48 a BP. However, the appearance of the sampling site and the presence of the same types of minerals in all samples suggest that they likely formed near-synchronously. The results yielded a highly correlated isochron (R2=0.998), yet the maximum likelihood regression is heavily controlled by the sample with the highest analytical precision, which is the cleanest sample (KC19CCC-2). The high mean squared weighted deviation (MSWD) of 69 indicates that the ages are overdispersed compared to the stated analytical uncertainties and that an isochron accounting for this overdispersion is needed (Vermeesch, 2018). The resulting age of 65±17 a BP agrees with the ages of two of the three samples (KC19CCC-2 is 64±12 a BP; KC19CCC-3 is 172±153 a BP) when corrected for detrital Th, using the isochron-derived initial 230Th/232Th activity of 1.39±0.006 (Fig. 5). The isochron-corrected age of KC19CCC-4 (353±74 a BP) is, however, not in agreement with the other two samples within dating uncertainty, thus suggesting multiple sources of initial 230Th (e.g. Dorale et al., 2004). Nonetheless, this disagreement could also result from an underestimation of the uncertainty by the maximum likelihood regression. The isochron-corrected δ234Uinitial is high for all three samples (1872±3 to 2193±5).

Figure 5Maximum likelihood isochron with overdispersion, based on the three 230Th / U disequilibrium ages and calculated in IsoplotR (Vermeesch, 2018).

Table 1Results of 230Th / U disequilibrium dating.

All uncertainties are 2σ. The U decay constants are λ238=1.55125×10-10 (Jaffey et al., 1971) and λ234=2.82206×10-6 (Cheng et al., 2013). Th decay constant is λ230=9.1705×10-6 (Cheng et al., 2013). a δ234U=((234U/238U)activity-1)×1000. b Corrected for detrital Th with the initial 230Th/232Th activity of 1.39±0.006.

Download Print Version | Download XLSX

Conventional radiocarbon ages for the four untreated and four mixed aliquots range from 140±27 to 2261±74 a BP, with corresponding F14C values of 0.98±0.003 to 0.06±0.0006 (Table 2). When applying the radiocarbon method to speleothems, the DCF needs to be determined (Genty and Massault, 1997) in order to correct for 14C-free material from bedrock and/or soil (Bajo et al., 2017; Hajdas et al., 2021; Hua et al., 2012). The DCFs of all analysed aliquots range from 1.24±0.50 % to 93.67±0.06 % (1σ uncertainty). For the untreated aliquots, 96.2 % of variance in radiocarbon age can be explained by the DCF (R2=0.96). When artificial mixtures are included, this value is lower (R2=0.88), yet it indicates that the artificial mixing was reasonably accurate. The radiocarbon age extrapolated to 0 % dead carbon by linear approximation is 2±28 BP. Calibration using IntCal20 (Reimer et al., 2020) yielded three distinct peaks of calibrated ages, namely 36–73, 112–139, and 226–255 cal BP, with respective probabilities of 37.8 %, 28.3 %, and 29.4 % (Fig. 6). The peak with the highest probability (36–73 BP) is in agreement with the 230Th / U isochron age of 65±17 a BP and provides the most likely timing for CCM formation. We attribute the other two peaks to a plateau of the calibration curve in this time interval.

Table 2Results of 14C dating.

All uncertainties are 1σ.

Download Print Version | Download XLSX

Figure 6Calibrated and DCF-corrected 14C dating results and the 230Th / U isochron age plotted against temperatures from two boreholes on the Greenland ice sheet (Dahl-Jensen et al., 1998) indicate that the Cove Cave CCMs formed during the late Little Ice Age. The two older peaks of the calibrated and corrected 14C dates are likely the result of a plateau of the calibration curve at that time interval.


The dating results show that the CCMs formed during the late Little Ice Age, a period of relatively cold climate conditions and glacier advances, which in Greenland lasted between about 700 and 50 a BP (i.e. 1250–1900 CE; Kjær et al., 2022).

5 Discussion

5.1 Formation of CCMs

In carbonate-hosted caves, CCCfine typically show a simple mineralogical composition, mostly consisting of calcite or, less commonly, other carbonate minerals, such as aragonite, or mixtures thereof (Žák et al., 2018). In some caves, the occurrence of non-carbonate minerals of cryogenic origin has been reported (Dublyansky et al., 2017; Žák et al., 2018). Cryogenic gypsum and, locally, other sulfate minerals are mostly known from caves hosted in gypsum rock (Kadebskaya and Tschaikovskiy, 2015; Žák et al., 2018) but also from some limestone caves, where the sulfate is derived from pyrite oxidation (Bartolomé et al., 2023; Dublyansky et al., 2017). Mirabilite is known from several caves, where its existence has been tied to evaporative conditions (e.g. Audra and Nobécourt, 2013; Bieniok et al., 2011). Although it has been reported from cold caves (Harmon et al., 1983), and its formation has been tied to preceding cryochemical processes, mirabilite was not regarded as primarily cryogenic in origin (Žák et al., 2018), and we are unaware of reports about cryogenically formed mirabilite in caves. Eugsterite is a rare sulfate mineral which has only been reported from two caves, namely Mammoth Cave, USA (White, 2017), and Chamois Cave, France (Audra and Nobécourt, 2013), where its formation is not related to cryogenic processes. Löweite has only rarely been reported from volcanic caves (Hill and Forti, 1997).

In Cove Cave, the most likely source of dissolved sulfate, necessary for the precipitation of sulfate minerals, is finely disseminated pyrite in the dark-coloured Silurian limestones in which this cave developed (Smith and Rasmussen, 2020), whereas magnesium for löweite formation was likely sourced from dolostones in the area. In contrast, the presence of sodium, which is needed for the minerals mirabilite, löweite, and eugsterite, is more challenging to explain. One possibility is that, given the proximity to the coast, the Na is sourced from aerosols from the ocean and/or the sea ice surface, which has been recognised as a dominant source of Na in coastal Arctic ice cores (Rhodes et al., 2018). The most probable source for the traces of dolomite is the host rock, while traces of quartz and potassium feldspar are, most likely, detrital material that entered the cave via water influx or aeolian transport.

The mode of formation (during freezing or sublimation) of the hydrated sulfate minerals must also be considered, especially with respect to whether their formation was synchronous with the precipitation of the carbonate fraction. While careful examination of SEM images yielded no conclusive evidence of coeval crystal growth, the high solubility of eugsterite, mirabilite, löweite, and to a lesser extent gypsum can be used to provide insights. For instance, it is clear that the four sulfate minerals did not form prior to the influx of water and subsequent freezing; otherwise pre-existing soluble minerals would have been dissolved and/or washed away. Likewise, the high solubility of these sulfate minerals suggests that the ice in which the cryogenic carbonates were embedded did not melt but rather sublimated; otherwise, these delicate crystals would not have been preserved. Furthermore, the cave section hosting the sampling site must have remained dry and ice-free since the formation of the CCMs. This also implies that the cave ventilation pattern at the time of formation was already similar to the present mode in this sag-type cave.

Based on the observation that the calcite and sulfate minerals were found intermixed, it seems likely that the formation of the sulfate minerals was associated with the formation of the carbonate fraction, in the sense that both their formations were triggered by the same event and therefore took place near-synchronously. Whether this formation was cryogenic, or occurred during the subsequent sublimation of ice, cannot be determined with certainty. A cryogenic formation of gypsum is likely, similar to the observed formation in other ice caves (see above; e.g. Žák et al., 2018), but the other sulfate minerals may have formed during the subsequent sublimation of ice after gypsum and calcite had already formed.

5.2 Stable isotopes

The relatively high δ13C (7.0 ‰–11.4 ‰) values of the Greenland CCMs are comparable to those of CCCfine from the Canadian Arctic Circle (Clark and Lauriol, 1992); however, the δ13C of the Canadian samples reach higher values (up to 17 ‰). Nevertheless, these highly positive values reflect kinetic isotope fractionation as a result of rapid freezing and the associated degassing of carbon dioxide (Lacelle et al., 2009).

The Cove Cave CCMs show highly depleted δ18O (−24.0 ‰ to −16.0 ‰) values that are much lower than the δ18O values of mid-latitude CCCfine and overlap with the lower range of δ18O values of CCCcoarse from mid-latitude caves (Fig. 4). While CCCcoarse δ18O reflects closed-system freezing in small drip water pools in ice, the low values of the Greenland samples are largely related to the isotopically light meteoric precipitation in this high-Arctic setting. We are interpreting the low δ18O values, which can be observed in Cove Cave CCMs, Cove Cave common speleothems (−17 ‰ to −13 ‰), and Canadian CCCfine (although at a smaller magnitude; Clark and Lauriol, 1992; Fig. 4), to be the result of the high-latitude setting in which both sampling sites are located. While existing CCC data are often biased towards central European/mid-latitude sites, with the isotopic compositions of CCCcoarse and CCCfine plotting in distinct O and C ranges, our data demonstrate that such borders are not universally applicable, especially at higher latitudes.

Previous studies have shown that the δ18O values of CCCfine and common speleothems often overlap (Fig. 4; e.g. Luetscher et al., 2013; Žák et al., 2018); however, the δ18O values of Cove Cave CCMs are more depleted than those of common speleothems from the same cave. While our data set is not sufficient to draw definite conclusions, we hypothesise that the isotopic composition of the source water, from which the CCMs and common speleothems precipitated, differed. While the CCMs precipitated recently from melted snow/ice (winter signal), the common speleothems formed from rainwater (annual or summer signal) during an earlier period within the Quaternary under different climatic boundary conditions.

5.3 Formation age

In this study, fine-grained CCMs were successfully dated using both a 230Th/U isochron approach and 14C dating with DCF correction. The isochron indicated that the initial 230Th/232Th activity ratio (1.39±0.006) was elevated in comparison to the bulk-Earth-derived value (i.e. 0.8; Wedepohl, 1995), which had therefore resulted in an undercorrection for detrital Th (Table 1). For dating CCMs in this study, it is sufficient to solely use 230Th / U methods, but we also show that successful and reliable dating of fine-grained CCMs is possible using 14C, provided that the DCF can be constrained and corrected. The independent age information, necessary for the DCF calculation, however, does not have to be provided by 230Th / U dating; it could, for instance, come from 14C dating of stratigraphically coeval organic matter. Where 230Th / U dating of fine-grained CCMs is impracticable (e.g. due to low 238U concentrations and/or high detrital Th contamination), 14C dating could provide a viable alternative. In this study, the 14C dating results complement the 230Th / U isochron age well (Fig. 6). However, because the DCF-corrected 14C age hits a plateau on the calibration curve (Reimer et al., 2020), the 230Th / U method is considered a better approach in this study. On the other hand, 14C dating of fine-grained CCMs could be as reliable as, or better than, 230Th / U dating for periods of time characterised by a steeply sloped calibration curve. Ultimately, the choice of dating will be study-dependent, and regardless of the method used, the dating of fine-grained CCMs could provide additional chronological control on cave ice bodies, which is an emerging (but rapidly disappearing) palaeoclimate archive (e.g. Kern and Perşoiu, 2013; Racine et al., 2022), in particular in those settings where organic inclusions are lacking.

5.4 Extreme-event-triggered CCM formation

High δ234Uinitial values suggest that the karst system of the cave was hydrologically inactive prior to the event that caused CCM formation. During arid climate periods and/or under permafrost influence, 234U accumulates in the crystal lattice of bedrock minerals through alpha recoil and is easily mobilised when water becomes available (Fleischer, 1982). Speleothems, including CCMs, may record this first mobilisation of 234U after a period of accumulation by elevated δ234Uinitial values (e.g. Wendt et al., 2020). The high δ234Uinitial values of the Cove Cave CCMs (1872±3 to 2193±5) contrast with the comparatively low δ234Uinitial values of common speleothems from the same cave, which represent a hydrologically active karst system (<150). The high δ234Uinitial values of the Cove Cave CCM therefore indicate a prolonged period of permafrost presence prior to their formation.

Combining the results of mineralogical analyses and dating, it can be concluded that the Cove Cave CCMs formed as a result of a singular event at 65±17 a BP (1885±17 CE) and therefore during the cold climate of the late Little Ice Age (ca. 700–50 a BP; Kjær et al., 2022). A reconstruction of surface temperature, based on Greenland ice core data since 1840 CE, indicates that the mid- to late 1880s CE (65–60 a BP) were some of the coolest years on record (Box et al., 2009). It can therefore be assumed that the climate of northeastern Greenland during the time of CCM formation was colder and permafrost was more abundant than today. In order for water to enter the cave through the frozen rock and not immediately freeze and sublimate on the walls, it must carry enough latent heat. Any event that led to these conditions in and around Cove Cave must have been anomalous, which is supported by the lack of CCMs observed in other caves in the area (based on the findings of two expeditions that were part of this project).

In order to explain the formation of fine-grained CCMs in Cove Cave, we consider the following scenarios: (i) an extreme rainfall event providing enough water and thus latent heat to enter Cove Cave, (ii) migration of the hoar frost boundary and melt pond deeper into the cave due to the influx of warm and moist air, (iii) an increase in the active layer thickness of (wet) permafrost due to an anomalously warm summer, year, or interval, (iv) higher cave air temperatures and water availability due to temperate ice covering the cave, and (v) enhanced melting of the local ice cap.

There is no evidence of (i) an extreme rainfall event in the high-resolution Greenland ice cores during the time of CCM formation, although it cannot be excluded that such an event occurred locally or regionally. Though a tectonic fracture exists in the roof of Cove Cave, which could promote water infiltration, the lack of CCM formation in other caves in the area renders this scenario unlikely. Furthermore, in this arid region, with ca. 200 mm yr−1 of precipitation (Schuster et al., 2021), extreme rainfall events with high volumes of water do not occur. Scenario (ii) seems unlikely, as the hoar frost boundary is located at the high point in Cove Cave (Fig. 2). Such boundaries are only observed at the highest points in the Greenland Caves (Barton et al., 2020); hence, shifting this boundary deeper into the cave is unlikely to occur, based on the sag-type geometry of the cave, which determines the location of the cold pool interface. An unusually warm summer, year, or interval that could potentially lead to (iii) an increase in the active layer thickness would have been recorded by Greenland ice cores and observational records along the coast of southern and western Greenland. Based on a merged southwestern Greenland temperature record (Vinther et al., 2006), there is no summer, year, or interval that stands out as being particularly warm in the period of 1868–1902 CE (i.e. 65±17 a BP). Another argument against scenario (iii) is, again, the lack of CCMs in other caves in the area. While it is very likely that an ice cap covered the plateau above Cove Cave in the time interval of interest (Fig. 1), scenario (iv) fails to explain why there was only one generation of CCMs found in Cove Cave, since higher cave air temperatures and water availability would probably have lasted for a prolonged interval. The ice cap does, however, appear to be a key factor together with a tectonic fracture in the roof of the cave, as they enable a mechanism in which sufficient water can pass through the permafrost and enter the underground. Specifically, enhanced melting of the local ice cap (v) could have been triggered by anomalously high temperatures and concurrent lowering of the albedo due to black carbon deposition. Both of these factors have been recognised as the cause for ice-melting conditions all over the Greenland ice sheet, including the dry snow zone, in 1889 CE (Clausen et al., 1988; Fischer et al., 1998; Keegan et al., 2014; Neff et al., 2014). This event is often referred to as the “summer melt episode” of 1889 CE, which, according to Neff et al. (2014), lasted for a few days. The summer melt episode occurred synchronously within the dating uncertainty related to the timing of the CCM formation (1885 CE; 65±17 a BP). Closer to our study site, an ice core from the Flade Isblink ice cap recorded a warming between 1920 and 1930 CE (30–20 a BP) as increased melt percentage, but it did not record enhanced melt during 1889 CE, which might simply be the result of a non-definitive timescale that does not allow for the investigation of short-term events (Lemark, 2010). The ice core from Flade Isblink did, however, record a high concentration of black carbon in 1889 CE (Eckhardt et al., 2023). We therefore infer that the unusually high air temperatures and black carbon deposition associated with this event affected our study area as well, leading to enhanced melting of the local ice cap, with more water entering moulins, reaching the base of the ice cap, and then finding its way through conduits, possibly through tectonic fractures, into the cave. This water turned to ice in the heavily undercooled cave, resulting in the precipitation of fine-grained cryogenic cave carbonates and potentially also cryogenic gypsum during a rather rapid freezing process. Subsequently, the ice sublimated in the cold and dry microclimate of the cave, releasing enclosed cryogenic calcite and gypsum particles which accumulated on the cave floor. As discussed above, the origin of the other hydrous sulfate minerals is less clear; they were, however, found intermixed with calcite and gypsum in the same spot, hinting towards a near-synchronous formation.

The presence of CCMs in Cove Cave indicates that the conditions leading to the summer melt episode of 1889 CE also reached further northeast and thus affected the area surrounding the cave. Linking the CCM formation to this short-term extreme event further provides hints about the rate of formation of fine-grained CCMs (which is a highly understudied field), while demonstrating why fine-grained CCMs cannot be used as a permafrost proxy as the presence or absence of permafrost is generally not influenced by short-term weather conditions.

6 Conclusions

Fine-grained CCMs from Cove Cave, a low-elevation, high-latitude permafrost cave in northeastern Greenland, consist of an unusual mixture of minerals, namely cryogenic calcite, potentially cryogenic gypsum, and other hydrous sulfate minerals (eugsterite, mirabilite, and löweite), whose mode of formation is less clear but presumably at least associated with cryochemical processes and/or the subsequent sublimation of the ice body in a dry cave atmosphere. The high solubility of the sulfate minerals suggests that since their formation, and after the sublimation of ice, the microclimate at the sampling site has remained cold, dry, and ice-free – as it is today.

The Greenland CCMs extend the knowledge of the stable isotopic composition of CCMs. While their δ13C values are comparable to those of mid-latitude CCCfine, the δ18O values are lower compared to both common (inactive) speleothems from the same cave and CCCfine from the mid-latitudes. These lower δ18O values can be attributed to much lower δ18O values of meteoric precipitation in the high latitudes and, to a lesser extent, the difference in the isotopic composition of the water sources during the formation of CCMs and common speleothems.

We show that precise dating of fine-grained CCMs, and hence CCCfine, is in principle possible with 230Th / U or 14C dating. Using isochrons and a site-specific initial 230Th/232Th correction factor is paramount for accurate 230Th / U ages. In this study, 230Th / U is arguably superior to 14C dating, which is dependent on the slope of a calibration curve and independent age information to constrain the DCF. However, its application is restricted to samples of low detrital Th content (and rather high U concentration). Where these conditions cannot be met, 14C dating might be a better alternative. The possibility of dating fine-grained CCMs with either 230Th / U or 14C could aid the establishment of reliable chronologies for cave ice as a palaeoclimate archive.

We conclude that CCM formation in Cove Cave was most likely the result of a single short-term extreme event that also led to the summer melt episode of 1889 CE on the Greenland ice sheet, causing enhanced melting on the local ice cap on the plateau above the cave. The Cove Cave CCMs therefore give an indication regarding the spatial extent of melting conditions in an area outside of the Greenland ice sheet, while also providing hints about the rate of CCM formation.

Appendix A: Supplementary information on XRD analysis

XRD analyses were conducted in two laboratories during different stages of the project. The displayed X-ray diffractograms show the results of one of the laboratories at a later stage in the project. While the results between the labs were mostly comparable, mirabilite was found during an analysis of KC19CCC-2 and KC19CCC-3 at an early stage of the project. This result could not be replicated during analysis later in the project. We attribute this inconsistency to the imperfect homogenisation of the samples prior to analyses, which might have been enhanced by the small sample volumes required for XRD analysis.

Figure A1X-ray diffractogram of sample KC19CCC-1. This sample contains calcite, gypsum, quartz, dolomite, and potassium feldspar.


Figure A2X-ray diffractogram of sample KC19CCC-2. As displayed, this sample contains calcite, gypsum, eugsterite, and löweite. During XRD analysis in another laboratory, mirabilite was also found in the sample.


Figure A3X-ray diffractogram of sample KC19CCC-3, which is also representative of KC19CCC-4. As shown, this sample contains calcite, gypsum, eugsterite, löweite, and quartz. During XRD analysis in another laboratory, mirabilite was also detected.


Data availability

We used unpublished data from our project for comparison. The data will be made available publicly in the near future.

Author contributions

AD designed the methodology, interpreted the data, and wrote the paper. PT assisted with fieldwork, age calculations, and preparation of the paper. CS assisted with analyses, contributed to discussions, and helped prepare the paper. IH performed 14C dating. XL contributed with 230Th / U disequilibrium analyses. RLE provided analytical 230Th / U disequilibrium dating facilities. GEM (principal investigator) designed the study, raised the funding, organised, led, and participated in the fieldwork, and contributed to the discussions and preparation of the paper.

Competing interests

The contact author has declared that none of the authors has any competing interests.


Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


We thank Manuela Wimmer for stable isotope measurements, Clivia Hejny for XRD analyses, Kristian Pfaller for SEM images, Andreas Saxer for XRD analyses and SEM images, and Bastian Joachim-Mrosko for Raman spectroscopy. We are grateful to Hazel Barton, Chris Blakeley, Pete Hodkinson, Adam Ignézi, Robbie Shone, Paul Smith, and Andrew Sole, who participated in the 2019 Greenland Caves Expedition. Data from the Programme for Monitoring of the Greenland Ice Sheet (PROMICE) were provided by the Geological Survey of Denmark and Greenland (GEUS) at (last access: 6 October 2022). The Greenland government are thanked for permission to undertake this fieldwork (KNNO Expedition Permit C-19-32; Scientific Survey Licence VU-00150; Greenland National Museum and Archives 2019/01).

Financial support

This research has been supported by the Austrian Science Fund (grant no. Y 1162N37).

Review statement

This paper was edited by Alberto Reyes and reviewed by Connor Turvey and two anonymous referees.


Audra, P. and Nobécourt, J.-C.: Rare sulfates (mirabilite, eugsterite) in the dry microclimate of chamois cave (Alpes-de-Haute-Provence, France), in: Proceedings of the 16th International Congress of Speleology, Czech Speleological Societya, 21–28 July 2013, Prague, Czech Republic, 432–436, K26-00116, 2013. 

Bajo, P., Borsato, A., Drysdale, R., Hua, Q., Frisia, S., Zanchetta, G., Hellstrom, J., and Woodhead, J.: Stalagmite carbon isotopes and dead carbon proportion (DCP) in a near-closed-system situation: An interplay between sulphuric and carbonic acid dissolution, Geochim. Cosmochim. Ac., 210, 208–227,, 2017. 

Bartolomé, M., Sancho, C., Osácar, M. C., Moreno, A., Leunda, M., Spötl, C., Luetscher, M., López-Martínez, J., and Belmonte, A.: Characteristics of cryogenic carbonates in a Pyrenean ice cave (northern Spain), Geogaceta, 58, 107–110, 2015. 

Bartolomé, M., Cazenave, G., Luetscher, M., Spötl, C., Gázquez, F., Belmonte, Á., Turchyn, A. V., López-Moreno, J. I., and Moreno, A.: Mountain permafrost in the Central Pyrenees: insights from the Devaux ice cave, The Cryosphere, 17, 477–497,, 2023. 

Barton, H. A., Breley, G. J., Töchterle, P., and Moseley, G. E.: Cryogenic features of the permafrost ice caves of Grottedal, northeast Greenland, Cave Karst Sci., 47, 93–99, 2020. 

Bieniok, A., Zagler, G., Brendel, U., and Neubauer, F.: Speleothems in the dry cave parts of the Gamslöcher-Kolowrat Cave, Untersberg near Salzburg (Austria), Int. J. Speleol., 40, 117–124, 2011. 

Bintanja, R. and Krikken, F.: Magnitude and pattern of Arctic warming governed by the seasonality of radiative forcing, Sci. Rep.-UK, 6, 1–7,, 2016. 

Bintanja, R. and Selten, F. M.: Future increases in Arctic precipitation linked to local evaporation and sea-ice retreat, Nature, 509, 479–482,, 2014. 

Box, J. E., Yang, L., Bromwich, D. H., and Bai, L.-S.: Greenland Ice Sheet surface air temperature variability: 1840–2007, J. Climate, 22, 4029–4049,, 2009. 

Bronk Ramsey, C.: Bayesian analysis of radiocarbon dates, Radiocarbon, 51, 337–360,, 2009. 

Cheng, H., Edwards, R. L., Shen, C.-C., Polyak, V. J., Asmerom, Y., Woodhead, J., Hellstrom, J., Wang, Y., Kong, X., Spötl, C., Wang, X., and Alexander Jr., E. C.: Improvements in 230Th dating, 230Th and 234U half-life values, and U–Th isotopic measurements by multi-collector inductively coupled plasma mass spectrometry, Earth Planet Sc. Lett., 371–372, 82–91,, 2013. 

Clark, I. D. and Lauriol, B.: Kinetic enrichment of stable isotopes in cryogenic calcites, Chem. Geol., 102, 217–228,, 1992. 

Clausen, H. B., Gundestrup, N. S., Johnsen, S. J., Bindschadler, R., and Zwally, J.: Glaciological investigations in the Crête area, central Greenland: A search for a new deep-drilling site, Ann. Glaciol., 10, 10–15,, 1988.  

Dahl-Jensen, D., Mosegaard, K., Gundestrup, N., Clow, G. D., Johnsen, S. J., Hansen, A. W., and Balling, N.: Past temperatures directly from the Greenland ice sheet, Science, 282, 268–271,, 1998. 

Danmarks Meteorologiske Institut: Klimanormaler for Grønland, Danmarks Meteorologiske Institut [data set], (last access: 12 July 2022), 2022. 

Davies, W. E. and Krinsley, D. B.: Caves in Northern Greenland, Natl. Speleol. Soc. Bull., 22, 114–116, 1960. 

Donner, A., Töchterle, P., and Moseley, G. E.: Basic meteorological observations in the Centrumsø region of northeast Greenland, Cave Karst Sci., 47, 104–106, 2020. 

Dorale, J. A., Edwards, R. L., Alexander, E. C., Shen, C.-C., Richards, D. A., and Cheng, H.: Uranium-series dating of speleothems: current techniques, limits, & applications, in: Studies of Cave Sediments, edited by: Sasowsky, I. D. and Mylroie, J., Springer US, Boston, MA, 177–197,, 2004. 

Dublyansky, Y., Kadebskaya, O., Menshikova, E., Demeny, A., and Spötl, C.: Polymineral cryogenic deposits in caves of the Southern Ural, Russia, as tracers of past permafrost, in: Climate Change: The Karst Record 8th International Conference, Scientific Program and Abstracts, 21–24 May 2017, Austin, USA, p. 41, 2017. 

Eckhardt, S., Pisso, I., Evangeliou, N., Zwaaftink, C. G., Plach, A., McConnell, J. R., Sigl, M., Ruppel, M., Zdanowicz, C., Lim, S., Chellman, N., Opel, T., Meyer, H., Steffensen, J. P., Schwikowski, M., and Stohl, A.: Revised historical Northern Hemisphere black carbon emissions based on inverse modelling of ice core records, Nat. Commun., 14, 217,, 2023. 

Edwards, R. L., Chen, J. H., and Wasserburg, G. J.: 238U-234U-230Th-232Th systematics and the precise measurement of time over the past 500,000 years, Earth Planet. Sc. Lett., 81, 175–192,, 1987. 

Fausto, R. S., van As, D., and Mankoff, K. D.: Programme for monitoring of the Greenland ice sheet (PROMICE): Automatic weather station data, Version: v03, GEUS – Geological survey of Denmark and Greenland [data set],, 2019. 

Fischer, H., Werner, M., Wagenbach, D., Schwager, M., Thorsteinnson, T., Wilhelms, F., Kipfstuhl, J., and Sommer, S.: Little Ice Age clearly recorded in northern Greenland ice cores, Geophys. Res. Lett., 25, 1749–1752,, 1998. 

Fleischer, R. L.: Alpha-recoil damage and solution effects in minerals: uranium isotopic disequilibrium and radon release, Geochim. Cosmochim. Ac., 46, 2191–2201,, 1982. 

Genty, D. and Massault, M.: Bomb 14C recorded in laminated speleothems: calculation of dead carbon proportion, Radiocarbon, 39, 33–48,, 1997. 

Geofabrik and OpenStreetMap Contributors: Greenland, (last access: 13 September 2022), 2018. 

Harmon, R. S., Atkinson, T. C., and Atkinson, J. L.: The mineralogy of Castleguard Cave, Columbia Icefields, Alberta, Canada, Arct. Alp. Res., 15, 503–516, 1983. 

Hajdas, I., Ascough, P., Garnett, M. H., Fallon, S. J., Pearson, C. L., Quarta, G., Spalding, K. L., Yamaguchi, H., and Yoneda, M.: Radiocarbon dating, Nat. Rev. Meth. Prim., 1, 1–26,, 2021. 

Hill, C. and Forti, P.: Cave Minerals of the World, in: 2nd Edn., National Speleological Society, Huntsville, ISBN 13:978-1879961074, 1997. 

Hua, Q., McDonald, J., Redwood, D., Drysdale, R., Lee, S., Fallon, S., and Hellstrom, J.: Robust chronological reconstruction for young speleothems using radiocarbon, Quatern. Geochronol., 14, 67–80,, 2012. 

Jaffey, A. H., Flynn, K. F., Glendenin, L. E., Bentley, W. C., and Essling, A. M.: Precision measurement of half-lives and specific activities of 235U and 238U, Phys. Rev. C, 4, 1889–1906,, 1971. 

Kadebskaya, O. and Tschaikovskiy, I.: Staging in the hypergene transformation of sulphate and carbonate rocks (based on slide-rocks under organ tubes in Kungur ice cave), Acta Carsolog., 44, 237–250, 2015. 

Keegan, K. M., Albert, M. R., McConnell, J. R., and Baker, I.: Climate change and forest fires synergistically drive widespread melt events of the Greenland Ice Sheet, P. Natl. Acad. Sci. USA, 111, 7964–7967,, 2014. 

Kern, Z. and Perşoiu, A.: Cave ice – the imminent loss of untapped mid-latitude cryospheric palaeoenvironmental archives, Quaternary Sci. Rev., 67, 1–7,, 2013. 

Kjær, K. H., Bjørk, A. A., Kjeldsen, K. K., Hansen, E. S., Andresen, C. S., Siggaard-Andersen, M.-L., Khan, S. A., Søndergaard, A. S., Colgan, W., Schomacker, A., Woodroffe, S., Funder, S., Rouillard, A., Jensen, J. F., and Larsen, N. K.: Glacier response to the Little Ice Age during the Neoglacial cooling in Greenland, Earth-Sci. Rev., 227, 1–43,, 2022. 

Koltai, G., Spötl, C., Jarosch, A. H., and Cheng, H.: Cryogenic cave carbonates in the Dolomites (northern Italy): insights into Younger Dryas cooling and seasonal precipitation, Clim. Past, 17, 775–789,, 2021. 

Lacelle, D., Lauriol, B., and Clark, I. D.: Formation of seasonal ice bodies and associated cryogenic cave carbonates in Caverne de L'Ours, Québec, Canada: Kinetic isotope effects and pseudo-biogenic crystal structures, J. Cave Karst Stud., 71, 48–62, 2009. 

Lauriol, B. and Clark, I. D.: An approach to determine the origin and age of massive ice blockages in two arctic caves, Permafrost Periglac., 4, 77–85,, 1993. 

Lauriol, B., Carrier, L., and Thibaudeau, P.: Topoclimatic zones and ice dynamics in the caves of the Northern Yukon, Canada, Arctic, 41, 215–220, 1988. 

Lemark, A.: A study of the Flade Isblink ice cap using a simple ice flow model, Master thesis, Niels Bohr Institute Copenhagen University, Centre for Ice and Climate, Denmark, 75 pp., (last access: 2 August 2023), 2010. 

Loubière, J.-F.: Observations préliminaires sur les cavités de la région du lac Centrum (NE Groenland), Karstologia, 9, 7–16,, 1987. 

Ludwig, K. R. and Titterington, D. M.: Calculation of 230Th/U isochrons, ages, and errors. Geochim. Cosmochim. Ac., 58, 5031–5042,, 1994. 

Luetscher, M., Borreguero, M., Moseley, G. E., Spötl, C., and Edwards, R. L.: Alpine permafrost thawing during the Medieval Warm Period identified from cryogenic cave carbonates, The Cryosphere, 7, 1073–1081,, 2013. 

Moseley, G. E., Barton, H. A., Spötl, C., Töchterle, P., Smith, M. P., Bjerkenås, S. E., Blakely, C., Hodkinson, P. D., Shone, R. C., Sivertsen, H. C., and Wright, M.: Cave discoveries and speleogenetic features in northeast Greenland, Cave Karst Sci., 47, 74–87, 2020. 

Moseley, G. E., Edwards, R. L., Lord, N. S., Spötl, C., and Cheng, H.: Speleothem record of mild and wet mid-Pleistocene climate in northeast Greenland, Sci. Adv., 7, 1–13,, 2021. 

Neff, W., Compo, G. P., Ralph, F. M., and Shupe, M. D.: Continental heat anomalies and the extreme melting of the Greenland ice surface in 2012 and 1889, J. Geophys. Res.-Atmos., 119, 6520–6536,, 2014. 

Racine, T. M. F., Reimer, P. J., and Spötl, C.: Multi-centennial mass balance of perennial ice deposits in Alpine caves mirrors the evolution of glaciers during the Late Holocene, Sci. Rep.-UK, 12, 1–13,, 2022. 

Reimer, P. J., Brown, T. A., and Reimer, R. W.: Discussion: Reporting and calibration of post-bomb 14C data, Radiocarbon, 46, 1299–1304,, 2004. 

Reimer, P. J., Austin, W. E. N., Bard, E., Bayliss, A., Blackwell, P. G., Bronk Ramsey, C., Butzin, M., Cheng, H., Edwards, R. L., Friedrich, M., Grootes, P. M., Guilderson, T. P., Hajdas, I., Heaton, T. J., Hogg, A. G., Hughen, K. A., Kromer, B., Manning, S. W., Muscheler, R., Palmer, J. G., Pearson, C., van der Plicht, J., Reimer, R. W., Richards, D. A., Scott, E. M., Southon, J. R., Turney, C. S. M., Wacker, L., Adolphi, F., Büntgen, U., Capano, M., Fahrni, S. M., Fogtmann-Schulz, A., Friedrich, R., Köhler, P., Kudsk, S., Miyake, F., Olsen, J., Reinig, F., Sakamoto, M., Sookdeo, A., and Talamo, S.: The IntCal20 Northern Hemisphere radiocarbon age calibration curve (0–55 cal kBP), Radiocarbon, 62, 725–757,, 2020. 

Rhodes, R. H., Yang, X., and Wolff, E. W.: Sea ice versus storms: What controls sea salt in Arctic ice cores?, Geophys. Res. Lett., 45, 5572–5580,, 2018. 

Richter, D. K., Scholz, D., Jöns, N., Neuser, R. D., and Breitenbach, S. F.: Coarse-grained cryogenic aragonite as end-member of mineral formation in dolomite caves, Sediment. Geol., 376, 136–146,, 2018. 

Schuster, L., Maussion, F., Langhamer, L., and Moseley, G. E.: Lagrangian detection of precipitation moisture sources for an arid region in northeast Greenland: relations to the North Atlantic Oscillation, sea ice cover, and temporal trends from 1979 to 2017, Weather Clim. Dynam., 2, 1–17,, 2021.  

Shen, C.-C., Wu, C.-C., Cheng, H., Edwards, R. L., Hsieh, Y.-T., Gallet, S., Chang, C.-C., Li, T.-Y., Lam, D. D., Kano, A., Hori, M., and Spötl, C.: High-precision and high-resolution carbonate 230Th dating by MC-ICP-MS with SEM protocols, Geochim. Cosmochim. Ac., 99, 71–86,, 2012. 

Shepherd, T. G.: Effects of a warming Arctic, Science, 353, 989–990,, 2016. 

Smith, M. P. and Rasmussen, J. A.: The geology of the Centrumsø area of Kronprins Christian Land, northeast Greenland, and lithological constraints on speleogenesis, Cave Karst Sci., 47, 60–65, 2020. 

Sole, A. J., Ignéczi, Á., Smith, M. P., and Clark, C. D.: Investigations to constrain retreat of the Greenland Ice Sheet: glacial geomorphology and sampling for cosmogenic exposure dating of the Centrumsø area, Kronprins Christian Land, northeast Greenland, Cave Karst Sci., 47, 66–73, 2020. 

Spötl, C.: Kryogene Karbonate im Höhleneis der Eisriesenwelt, Höhle, 59, 26–36, 2008. 

Spötl, C.: Long-term performance of the Gasbench isotope ratio mass spectrometry system for the stable isotope analysis of carbonate microsamples, Rapid Commun. Mass. Spectrom., 25, 1683–1685,, 2011. 

Spötl, C. and Cheng, H.: Holocene climate change, permafrost and cryogenic carbonate formation: insights from a recently deglaciated, high-elevation cave in the Austrian Alps, Clim. Past, 10, 1349–1362,, 2014. 

Spötl, C. and Vennemann, T. W.: Continuous-flow isotope ratio mass spectrometric analysis of carbonate minerals, Rapid Commun. Mass Spectrom., 9, 1004–1006,, 2003. 

Spötl, C., Koltai, G., Jarosch, A. H., and Cheng, H.: Increased autumn and winter precipitation during the Last Glacial Maximum in the European Alps, Nat. Commun., 12, 1839,, 2021. 

Stuiver, M. and Polach, H. A.: Discussion reporting of 14C data, Radiocarbon, 19, 355–363,, 1977. 

Synal, H.-A., Stocker, M., and Suter, M.: MICADAS: A new compact radiocarbon AMS system, Nucl. Instrum. Meth. B, 259, 7–13,, 2007. 

Töchterle, P., Steidle, S., Edwards, R. L., Dublyansky, Y., Spötl, C., Li, X., Gunn, J., and Moseley, G. E.: 230Th/U isochron dating of cryogenic cave carbonates, Geochronology, 4, 617–627,, 2022. 

van As, D., Fausto, R. S., Ahlstrøm, A. P., Andersen, S. B., Andersen, M. L., Citterio, M., Edelvang, K., Gravesen, P., Machguth, H., Nick, F. M., Nielsen, S., and Weidick, A.: Programme for monitoring of the Greenland Ice Sheet (PROMICE): first temperature and ablation records, Geol. Surv. Den. Greenl., 23, 73–76,, 2011. 

Vermeesch, P.: IsoplotR: A free and open toolbox for geochronology, Geosci. Front., 9, 1479–1493,, 2018. 

Vinther, B. M., Andersen, K. K., Jones, P. D., Briffa, K. R., and Cappelen, J.: Extending Greenland temperature records into the late eighteenth century, J. Geophys. Res., 111, D11105,, 2006. 

Wedepohl, K. H.: The composition of the continental crust, Geochim. Cosmochim. Ac., 59, 1217–1232,, 1995. 

Wendt, K. A., Pythoud, M., Moseley, G. E., Dublyansky, Y. V., Edwards, R. L., and Spötl, C.: Paleohydrology of southwest Nevada (USA) based on groundwater 234U/238U over the past 475 k.y., Geol. Soc. Am. Bull., 132, 793–802,, 2020. 

White, W. B.: Mineralogy of Mammoth Cave, in: Mammoth Cave, edited by: Hobbs III, H. H., Olson, R. A., Winkler, E. G., and Culver, D. C., Springer, Cham, 145–162,, 2017.  

Žák, K., Richter, D. K., Filippi, M., Živor, R., Deininger, M., Mangini, A., and Scholz, D.: Coarsely crystalline cryogenic cave carbonate – a new archive to estimate the Last Glacial minimum permafrost depth in Central Europe, Clim. Past, 8, 1821–1837,, 2012. 

Žák, K., Onac, B. P., Kadebskaya, O. I., Filippi, M., Dublyansky, Y., and Luetscher, M.: Cryogenic mineral formation in caves, in: Ice Caves, edited by: Perşoiu, A. and Lauritzen, S.-E., Elsevier, Amsterdam, the Netherlands, 123–162,, 2018. 

Short summary
This study investigates the first finding of fine-grained cryogenic cave minerals in Greenland, a type of speleothem that has been notably difficult to date. We present a successful approach for determining the age of these minerals using 230Th / U disequilibrium and 14C dating. We relate the formation of the cryogenic cave minerals to a well-documented extreme weather event in 1889 CE. Additionally, we provide a detailed report on the mineralogical and isotopic composition of these minerals.