Professional Documents
Culture Documents
Quantification of Paleo-Aquifer Changes Using Clumped Isotopes in Caves
Quantification of Paleo-Aquifer Changes Using Clumped Isotopes in Caves
PII: S0009-2541(18)30302-4
DOI: doi:10.1016/j.chemgeo.2018.05.046
Reference: CHEMGE 18799
To appear in: Chemical Geology
Received date: 6 April 2018
Revised date: 26 May 2018
Accepted date: 31 May 2018
Please cite this article as: Fernando Gázquez, Andrea Columbu, Jo De Waele, Sebastian
F.M. Breitenbach, Ci-Rong Huang, Chuan-Chou Shen, Yanbin Lu, José-María Calaforra,
Maryline J. Mleneck Vautravers, David A. Hodell , Quantification of paleo-aquifer
changes using clumped isotopes in subaqueous carbonate speleothems. Chemge (2017),
doi:10.1016/j.chemgeo.2018.05.046
This is a PDF file of an unedited manuscript that has been accepted for publication. As
a service to our customers we are providing this early version of the manuscript. The
manuscript will undergo copyediting, typesetting, and review of the resulting proof before
it is published in its final form. Please note that during the production process errors may
be discovered which could affect the content, and all legal disclaimers that apply to the
journal pertain.
ACCEPTED MANUSCRIPT
(1) School of Earth and Environmental Sciences. University of St. Andrews. St Andrews,
KY16 9AL, Scotland, United Kingdom. [email protected]
PT
(2) Italian Institute of Speleology, Department of Biological, Geological and Environmental
Sciences, University of Bologna, Via Zamboni, 67, 40126 Bologna, Italy.
(3) Sediment & Isotope Geology, Institute of Geology, Mineralogy & Geophysics. Ruhr,
RI
Universität Bochum, Universitätsstr, Bochum, Germany
SC
(4) High-Precision Mass Spectrometry and Environmental Change Laboratory (HISPEC),
Department of Geosciences, National Taiwan University, Taipei City 106, Taiwan ROC.
(5) Department of Earth Sciences, University of Minnesota, 310 Pillsbury Drive SE,
NU
Minneapolis, MN 55455-0231, USA.
(6) Water Resources and Environmental Geology Research Group. University of Almería.
MA
*corresponding author
PT
Abstract
CE
Here we track the water-table position and temperature of the Mount San Giovanni
aquifer (Iglesiente-Sulcis mining district, SW Sardinia, Italy) during the past 600 ka by
AC
determining the ages (U-Th dating) and stable isotope compositions (18O, 13C and
47) of a variety of subaqueous carbonate speleothems (e.g. calcite spars, dogtooth
calcite crystals and calcite coatings). Clumped isotopes (47) provide quantitative
estimates of carbonate formation temperatures (and thus water temperatures) that are
independent of the oxygen isotope composition of water (18Ow). Then, the 18Ow of
the paleo-water has been reconstructed from the clumped isotope temperature (T47)
and the 18O of the carbonate (18Oc). We find that high-temperature calcite spars
formed already before 600 ka at temperatures above ~120 oC. Lower-temperature
ACCEPTED MANUSCRIPT
spars (~70 oC) precipitated at ~400 ka, and cold-water subaqueous speleothems (~10-
20 oC) formed in perched ponds at different levels of the karst systems between 410
ka and 110 ka, while coeval precipitation of subaerial flowstones occurred in the upper
levels of the shallower caves until 82 ka. We infer that the groundwater level dropped
by ~120 m from ~400 to ~250 ka, with a relatively rapid rate of ~0.8 mm/y. Considering
the tectonic stability of Sardinia during the Quaternary, this high rate derives from
climate driven geomorphological processes at the surface rather than tectonic
uplifting. The 18O values of the paleo-aquifer water range from -6.0±0.7 ‰ during MIS
PT
5c, and similar to modern cave water values (-5.1 ‰), to -7.7±0.4 ‰ during the colder
RI
MIS 8. These values indicate that the groundwater reflected the 18O signal of meteoric
water, with no significant contributions from metasomatism and metamorphic waters.
SC
The observed 18Ow variability can be explained by glacial/interglacial paleoclimate
changes affecting rainfall 18O (and thus groundwater 18O). We conclude that
NU
clumped isotope thermometry on subaqueous carbonate speleothems is a powerful
tool for tracking paleo-aquifer temperatures and 18Ow reconstructions.
MA
1. Introduction
was thought only a decade years ago, a common karstification process that generated
dozens of caves worldwide (see Klimchouk et al., 2017, for a recent review on
AC
Cunningham, 1995; Dublyansky and Spötl, 2009; Palmer and Palmer, 2012; Decker
et al., 2015).
PT
is given by the reaction:
RI
Ca2+-HCO3--rich waters can precipitate carbonate minerals (i.e. calcite or aragonite)
SC
when the equilibrium reaction displaces to the left. This is generally accompanied by
CO2-degassing, evaporation and/or water temperature lowering (Audra et al., 2009).
NU
Such conditions are prone to occur when the fluids rise to the upper level of
hydrothermal aquifers and the hydrostatic pressure decreases, as water approaches
MA
the top of the aquifer, moving away from the deeper thermal source. Indeed, the
shallower sectors of many hypogene karst systems generally show greater abundance
of secondary carbonate formations, compared to deeper areas that often display
D
mainly dissolution features (Audra et al., 2009; Gázquez and Calaforra, 2013).
E
on their position with respect to the water table at the time of mineral precipitation
Subaqueous speleothems can furthermore be subdivided into phreatic (usually formed
CE
tens of meters below the water table), epiphreatic (formed in shallower aquifer levels
that can be emerged because of occasional water table lowering) (Jennings, 1985)
AC
and vadose (formed in pools/lakes in caves well above the water table). Note that in
other studies, the term “phreatic” has been used to refer to overgrowth speleothems
(POS) deposited close to the water table in anchialine caves (e.g. Fornós et al., 2002),
which would be classified as epiphreatic speleothems following the classification used
in the present study.
coatings (Polyak et al., 2013; De Waele et al., 2017; Gázquez et al., 2017), whereas
folia (e.g. Audra et al., 2009; Leél-Őssy et al., 2011; Davis, 2012; Gázquez and
Calaforra, 2013; Polyak et al., 2013), calcite rafts and calcite raft cones (e.g. Jones,
1989; Hill and Forti, 1997; Davis, 2000; Piccini, 2000; Audra et al., 2002; Leél-Őssy et
al., 2011; Polyak et al., 2013; Gázquez and Calaforra, 2013; Piccini et al., 2015), are
typical of shallower (epi)phreatic levels or the oscillation zone. The mechanisms
involved in the formation of each type of speleothem are generally associated with
specific environmental conditions. For example, calcite raft cones form exclusively in
PT
relatively shallow (less than few meters deep) subterranean water pools, and
therefore, can be used as indicators of paleo-aquifer levels (Gázquez and Calaforra,
RI
2013; De Waele et al., 2017; Kovacs et al., 2017). However, it is worth mentioning that
SC
some types of subaqueous carbonate speleothems are not exclusive of hypogene
caves but are also found in low-temperature epigene caves (e.g. cave cones and folia;
NU
Wildberger, 1987; Plan and De Waele, 2011; Davis, 2012; Fornós et al., 2002;
D’Angeli et al., 2015a; De Waele et al., 2017; Kovacs et al., 2017).
MA
temperature at the time of mineral precipitation. This is because the oxygen isotope
E
fractionation factor () between HCO3- (previously equilibrated with H2O) and CaCO3
PT
by the expression:
PT
limited in many cases because of the lack of understanding of
equilibrium/disequilibrium processes during carbonate precipitation from thin water
RI
films in vadose cave environments (Affek et al., 2008; Daëron et al., 2011; Kluge et
al., 2014a). However, recent clumped isotope analyses of mammillary coatings in
SC
Devil’s Hole (Nevada, US) suggest that subaqueous carbonate formed under isotope
equilibrium conditions between 180 ka and 27 ka, and therefore, water temperatures
NU
could be calculated (Kluge et al., 2014b). Also, T47 obtained from recent subaqueous
cave pearls formed in different climatic settings, matches the water temperature of the
MA
pool waters from which it formed (Breitenbach et al., 2018). This suggests that, under
certain conditions, subaqueous cave carbonates can precipitate in thermodynamic
and isotopic equilibrium with water and clumped isotope thermometry can be used for
D
clumped isotope thermometry has been recently applied to infer the temperature of
PT
paleo-waters by using calcite spars (Budd et al., 2013) and hydrothermal dolomite
(Ferry et al., 2011) formed in deep environments.
CE
Here we combine clumped isotope thermometry and classic 18O and 13C
measurements in calcite to reconstruct paleo-water temperature and 18Ow in Mount
AC
San Giovanni (SW Sardinia, Italy) during the Late Quaternary (<600 ka to 82 ka). The
various crystalline habits and arrangement of speleothems in the caves of Mount San
Giovanni suggest that carbonate formed under different environmental conditions,
including changing distances from the water surface (subaqueous/subaerial
conditions), water geochemistry and temperatures. However, the exact environmental
settings in which these speleothems crystallized have not been constrained yet. In
addition to temperature, clumped isotopes on subaqueous speleothems have been
used to calculate 18Ow values of the paleo-aquifer. The significance of the different
types of speleothems and potential paleoclimatic implications are discussed,
ACCEPTED MANUSCRIPT
regarding oscillations in the water table, water geochemistry and temperature changes
with time.
2. Geological setting
Mount San Giovanni Mine is located in the Iglesiente region (SW Sardinia, Italy), a Pb-
Zn mining district active since pre-Roman times (Cidu et al., 2009) (Fig. 1). Most mines
exploited Mississippi Valley Type ore bodies hosted in Lower Cambrian limestones
PT
and dolostones (Bechstadt and Boni, 1996). The carbonate succession starts with
grey dolostones (Santa Barbara Fm.), followed by darker massive dolostones and a
RI
thick succession of whitish-grey intensely karstified limestones (San Giovanni Fm.).
SC
The San Giovanni Mine exploited a very rich E-W aligned, 2 km long and 700 m wide
ore body (Brusca and Dessau, 1968) hosted in the southern limb of the Iglesias
NU
syncline (Civita et al., 1983). This ore was excavated from the top of the mountain
(400 m a.s.l.) to around 250 m below sea level, and around 100 km of mine galleries
MA
have been excavated through Mount San Giovanni. Ore exploitations lasted until
1998, when pumping operations in a nearby (connected) mine was discontinued in
1997 and the water table rapidly rose back to its original level at around 60 m a.s.l.
D
During tunnel construction many karst voids were intercepted. A total of over 50 caves
PT
have been mapped in this mine, most of which lack a natural entrance (so called “mine
caves”; Messina et al., 2005; Pagliara et al., 2010; De Waele et al., 2013; 2017;
CE
Gázquez et al., 2013b). Despite the fact that some caves are known since a very long
time (the most interesting from a scientific point of view being Santa Barbara 1 cave,
AC
Mine levels hosting significant numbers of caves are located between 170 and 226 m
a.s.l. (Figs. 2 and 3). In the Peloggio tunnel (218 m a.s.l.) a sequence of small caves
was intercepted, including Crovassa Azzurra cave (Caddeo et al., 2011). This cave is
composed of two superposed chambers, connected by a shaft, and hosts a wide
variety of speleothems and mineralizations, like blue aragonite flowstones, and calcite
ACCEPTED MANUSCRIPT
dogtooth spars. The Idina gallery (150 m a.s.l.) intercepts a large mineralized karst
void extending from this level up to the surface (370 m a.s.l.). Between 160 and 180
m a.s.l. a section of a speleothem-rich natural cave has been saved from mining
activities and is known as Crovassa Ricchi in Argento cave. The richest mine tunnel
is known as Albert (~226 m a.s.l.), which intercepted at least eight important caves.
Massa Riccardo cave is located in its upper portion and descends into a large room,
the floor of which is located at around 205 m a.s.l. (Figs. 2 and 3).
PT
3. Materials and methods
RI
3.1. Samples
SC
In 2005 a core 8 cm in diameter and 0.5 m long (sample GA-core) was retrieved from
an aragonite/calcite flowstone located at ~2 m below the cave entrance of Crovassa
NU
Azzurra (216 m a.s.l., Figs. 3 and 4E). The core was sampled using a hand-held water-
cooled HILTI® drilling machine with a 30 cm long drill bit and extensions. The core was
MA
split lengthwise and 11 subsamples for U-Th dating were taken (Fig. 4D). In the 2012
sampling campaign, the void left by the 2005 core extraction was found filled with
dripwater, and a sample of this water was collected for isotopic analysis.
E D
In addition to the cores samples, carbonate samples from three caves of the Mount
PT
San Giovanni mine were collected in 2012 (Table 1, Fig. 3). Six samples were
collected in Crovassa Azzurra cave. A fragment of calcite coatings detached from the
CE
cave ceiling was collected at ~5 m below the cave entrance (208 m a.s.l., Fig. 3). The
inner part of this sample consists of a pristine scalenohedral calcite spar (~10 cm in
size) that is coated with a layered calcite crust (~15 cm thick) composed of whitish and
AC
brownish laminae (Fig. 4A). Two additional scalenohedral calcite spars (~15 cm in size
each) were collected from the cave floor at 205 m a.s.l.. Also, three dogtooth calcite
speleothems were taken from the cave walls at ~12 m below the cave entrance (at
201 m a.s.l. in Fig. 3 and Fig. 4B). A fragment of blue aragonite flowstone was collected
from the cave floor (208 m a.s.l.) as well as a sample of the carbonate host rock.
In Crovassa Ricchi in Argento cave, four scalenohedral calcite spars were collected in
2012 from the floor in different parts of the cave, at around 180 m a.s.l. (Figs. 3). In
Massa Riccardo cave (Fig. 3), fragments of a whitish calcite flowstone (~20 cm thick)
ACCEPTED MANUSCRIPT
that deposited on a scalenohedral calcite spar (Fig. 4D), were collected at 206 m a.s.l..
Two samples of euhedral calcite, composed of aggregates of ~3 cm in size planar
calcite crystals resembling “roses” were collected from the cave wall (Fig. 4C) at 208
m a.s.l.. A fragment of the carbonate host rock was taken from the wall of Massa
Riccardo cave.
Homogeneous carbonate powders from the cleanest part of all the specimens
described above were obtained using a Dremel® drill tool. In addition, we consider for
PT
interpretations the U-Th age of a mammillary calcite coating (cave clouds) taken in
2005 in Santa Barbara 2 cave at 65 m a.s.l. and reported by De Waele and Forti (2006)
RI
(SB2-CC in Fig. 3).
subsamples of drilled flowstone cores from Crovassa Azzurra cave were analyzed for
U-Th dating at the High-precision mass spectrometry and environment change lab
D
isotope dilution method was employed to correct mass bias and determine U-Th
contents and isotopic compositions (Cheng et al., 2000; Shen et al., 2012). U-Th
CE
al., 2012; Cheng et al., 2013). Half-lives of U-Th nuclides used for 230Th age calculation
are given in Cheng et al. (2013). Ages were corrected by assuming an initial
230Th/232Th atomic ratio of 4.4±2.2 x10-6 (see Table S1 for details). Uncertainties in the
U-Th isotopic data and uncorrected-corrected ages are given at the two sigma (2)
level or two standard deviations of the mean (2m) unless otherwise stated (Table S1).
Powdered carbonate samples (n=32, ~200 g each) were flushed with CP grade
helium in Exetainer tubes (Labco®), then acidified with 104% H3PO4 and left to react
for 1h at 70 oC. Stable carbon and oxygen isotopes of CO2 were measured using a
ThermoScientific GasBench II, equipped with a CTC autosampler coupled to a
MAT253 mass spectrometer (Spötl and Vennemann, 2003) in the Godwin Laboratory
at the University of Cambridge, UK. Analytical precision was estimated at ±0.03‰ for
13C and ±0.08 ‰ for 18O by repeated analysis of an inhouse Carrara Marble
PT
standard. Results are reported in ‰ relative to the Vienna PeeDee Belemnite (V-PDB)
reference.
RI
3.4. Clumped isotopes
SC
Seven samples of subaqueous carbonate speleothems (i.e. calcite spars, dogtooth
calcite crystals and calcite coatings), were investigated for their clumped isotope
NU
composition (47, Table 1 and supplementary Table S2). Note that no samples from
the GA-core were analyzed for 47, since it is interpreted as having formed in subaerial
MA
conditions and may have been affected by isotopic disequilibrium during precipitation
(Affek et al., 2008; Daëron et al., 2011; Kluge et al., 2014a). Powder aliquots of 120 to
140 g each were loaded in borosilicate vials and analyzed on a KIEL-IV device,
D
at the University of Cambridge, UK (for details see Breitenbach et al., 2018). Between
PT
16 and 37 aliquots per sample were analyzed over a period of several months. Acid
digestion using 104 % orthophosphoric acid was performed at 70 oC. The four
CE
carbonate standards ETH1, ETH2, ETH3 and ETH4 were used for normalization
(Meckler et al., 2014) and errors are reported as 2-standard error of the mean of the
aliquots measured. The relationship between 47 and temperature for low temperature
AC
carbonate speleothems (i.e. 47 >0.5‰) was determined using our inhouse
calibration, which is based on natural subaqueous cave carbonates formed at
temperatures between 3 and 47 oC (Breitenbach et al. 2018, 47 =
(0.0447±0.006)*106/T2+(0.149±0.07); where T in K and 47 in ‰; Eq. 2). For higher
temperatures (i.e. 47 <0.5‰), we used the 47-temperature calibration of Kele et al.
(2015) (47 = (0.044±0.005)*106/T2+(0.205±0.0504); where T in K and 47 in ‰; Eq.
3). Both calibration equations were obtained using the same parameters and
ACCEPTED MANUSCRIPT
standards, but the Breitenbach et al. (2018) calibration is used here for low-
temperature samples, so it is based on the same type of subaqueous cave carbonates.
4. Results
238U content in the dated samples varies between aragonite and calcite, although in
both polymorphs concentrations are relatively high (Table S1). In aragonite, 238U is
PT
always higher than 1000 ppb, while in calcite the average is 476 ppb (min: 8 ppb, max:
1642 ppb). The content of non-authigenic 232Th is higher for calcite and lower for
RI
aragonite (Table S1); however, the 230Th/232Th activity ratio (which together with the
error associated with the calculation of 238U/234U/230Th/232Th rules the uncertainties
SC
associated with the final ages) is generally high (>50). Consequently, all ages reported
have relatively low 2errors except sample GA-CORE-2, which is close to the U-Th
NU
dating limit).
U-Th results attest that speleothems in the Mount San Giovanni karst system formed
MA
between >600 ka and 82.3±0.3 ka BP (Table 1). The “>600 ka” notation refers to
samples out of the range of the dating method, which is currently between 600 and
D
Azzurra (GA-core) started to grow beyond this limit, with basal ages of 684±244 and
PT
>500 ka BP (GA-core-2 and 3, Table 1). Its middle section, mostly characterized by
cm-long fiber calcite, grew from 570±53 to 506±19 ka BP (GA-core 4 and 5), although
CE
an evident hiatus occurs between 587±44 and 581±51 ka BP (GA-core 11 and 12),
indicated by the presence of muddy-like brownish materials. A similar hiatus appears
AC
4.2. 18O and 13C analyses in carbonates and clumped isotopes (47)
PT
thermometry
RI
Carbonate oxygen and carbon isotope values show two distinct groups of samples
(Fig. 5). A first cluster is characterized by relatively low 18Oc values (-18.8 ‰ to -13.6
SC
‰), i.e. distinctly lower than the Cambrian carbonate host rock (-11.7 ‰ to -10.4 ‰),
and relatively high and homogeneous 13C values (-0.3 ‰ to +0.2 ‰), similar to the
NU
host rock (-0.8 ‰ to -0.2 ‰). This cluster comprises all calcite spars, except one (MS-
06-03) from Massa Riccardo cave that shows considerably higher 18O (-7.0 ‰) and
MA
The second cluster comprises the remaining samples, which show rather
D
homogeneous 18O values, ranging -5.2 ‰ to -4.1 ‰, but a wider range of 13C, from
-11.1 ‰ to -4.3 ‰. In general, aragonite flowstone shows more positive 13C values
E
PT
(>-10 ‰), while most samples (i.e. calcite flowstones and calcite “roses”) are grouped
around -10.5 ‰. The dogtooth calcite crystals from Crovassa Azzurra cave show
CE
The seven samples analyzed for clumped isotopes show 47 values ranging from
AC
0.476 ‰ to 0.708 ‰ (Table 1 and supplementary Table S2). The uncertainties of the
clumped isotope analyses are between 0.004 ‰ and 0.009 ‰ (1SE) or 0.008 to 0.018
‰ as margin of error at the 95% confidence level. All the clumped isotope uncertainties
and propagated errors for T47 and 18Ow are given as 2-standard errors of the mean
(2SE). The spar samples GA-03-1, GRA-01-1 and GRA-01-3 showed the lowest 47
value (0.476±0.018 ‰, 0.485±0.009 ‰ and 0.491±0.012 ‰, respectively), whereas
the calcite crust GA-05-03 (0.696±0.012) and the dogtooth crystals (GA-01-01 and
GA-01-02; 0.680±0.013 ‰ and 0.694±0.013 ‰, respectively) of Crovassa Azzurra
ACCEPTED MANUSCRIPT
cave and the calcite “roses” MS-01 of Massa Riccardo cave (0.708±0.008 ‰) display
higher 47. After calibration using eq. 2 for 47 values >0.5 and eq. 3 for 47 values
<0.5, the higher formation temperature corresponds to the calcite spars GA-03-01
(130±14 oC), GRA-01-1 (123±6 oC) and GRA-01-3 (119±9 oC), whereas the rest of the
samples formed at lower temperatures, ranging from 10±2 oC to 17±3 oC (Table 1).
5. Discussion
PT
The clumped isotope results suggest that the paleo-aquifer of Mount St. Giovanni
RI
experienced pronounced cooling from the earlier (>600 ka) stages, when calcite spar
precipitated under deep-phreatic conditions, to the later (<410 ka) phases of
SC
epiphreatic and subaerial speleothem formation. According to the 47 data from calcite
spars (samples GA-03-1, GRA-01-1 and GRA-01-3, Table 1), the water temperature
NU
in the caves of the Mount San Giovanni was very high (probably >120 oC) before >600
ka, but lowered until 410±14 ka and ranged from 10±2 oC to 17±3 oC till at least ca.
MA
109±0.5 ka. Given the limited number of samples analyzed and the analytical
uncertainty of the measurements, small T 47 variations (e.g. 13±3 oC to 17±3 oC) are
not interpreted here. Monitoring in Santa Barbara 1 cave carried out between 2001
D
and 2003 has shown cave temperature to range between 15 and 16°C, and these
E
temperatures are comparable to those measured from time to time in the other caves
PT
located roughly at the same altitude (Chiesi, 2005). This suggests that water
temperature at the ~210 m level has been cold (i.e. 10 to 20 oC) from at least 410±14
CE
ka until present.
Combining T47 and 18Oc from the same speleothem, the 18Ow of the paleo-water
AC
can be calculated. The 18Ow value of the speleothem-forming water at ~410 ka was -
6.4±0.7 ‰ (sample GA-05-3) during the warmer-than-today marine isotope stage 11
(MIS 11). During the similarly warmer MIS 7 (~216 ka) and MIS 5c (~109 ka) the 18Ow
were -6.0±0.7 ‰ and -6.0±0.7 ‰, respectively. These values roughly agree with the
modern water 18Ow in the cave (-5.1 ‰) and the annual weighted mean 18Ow of
rainfall in Sardinia (-4.8±1.8 ‰, n=16, at the IAEA stations at Cagliari and Capo
Caccia, between 2001 and 2002). Previous investigations of the microthermometry of
fluid inclusions coupled with 18Oc measurements on epiphreatic calcite speleothems
ACCEPTED MANUSCRIPT
(i.e. carbonates formed in pools) of unknown ages in caves of Masua mine (~15 km
North from Mount San Giovanni), produced also low temperatures (15-20 oC) and
18Ow values of -4 to -5 ‰ (De Vivo et al., 1987), in good agreement with our
observations. Interestingly, the 18Ow of the speleothem-forming water at ~264 ka,
during the colder MIS 8 (-7.7±0.4 ‰,) was ~2 ‰ lower than the modern 18Ow in the
cave. These results suggest that the 18Ow of the cave water (and thus that of the
rainfall over Sardinia) may have varied at least ~2 ‰ during the past 400 ka, including
PT
warmer (MIS 11, MIS 7 and MIS 5c) and colder (MIS 8) periods. Such differences
might be explained by the temperature effect on the 18O of rainfall that tends to be
RI
lower during colder periods (Clark and Fritz, 1997), or other paleoclimatically driven
processes affecting rainfall at glacial/interglacial timescale (Fairchild et al., 2006).
SC
The undated – although likely >600 ka, according to other U-Th dated spars – high-
NU
temperature calcite spar from Crovassa Azzurra and Crovassa Ricchi in Argento
caves (GA-03-1, GRA-01-1 and GRA-01-3) show 18Ow values between -3.1±1.5 ‰
and -4.3±1.0 ‰. These values are slightly higher than the modern 18O observations
MA
of the rainfall and cave water and also higher than the reconstructed <600 ka-old
paleo-water 18O. This may have been related to seawater ingression in the aquifer
D
evaporation because of the “boiling” that enriched the water in 18O or some
contributions of water from deep geological processes (e.g. metasomatism and
metamorphic water), which generally shows higher 18O (Clark and Fritz, 1997).
CE
It can be assumed that the 18Ow of the recharge to the aquifer is equal to the mean
AC
18Ow of the rainfall in this region and varied between -5 ‰ and -6.5 ‰ during the past
600 kyrs. Then, using mean 18Ow of 5.5±1 ‰, we apply the 18Oc-18Ow approach to
estimate the formation temperatures of all our carbonate speleothems, including those
that were not analyzed for clumped isotopes (using eq. 1 for calcite samples, and the
temperature equation for inorganic aragonite by Kim et al., 2007). Most speleothems
formed at temperatures between 11 oC and 20 oC, except for calcite spars which return
precipitation temperatures ranging from 33 oC to 110 oC (T, 18Oc-18Ow in Table 1). The
discrepancies between the 47-temperatures and the 18O 18
c- Ow-approach
ACCEPTED MANUSCRIPT
temperatures (up to ~20 oC in the case of some high-temperature spars) are assumed
to be caused by temporal variations in 18Ow of the paleo-water, as discussed above.
PT
reconstruct the evolution of the Mount San Giovanni paleo-aquifer during the Late
Quaternary. In this karst system, subaqueous speleothem precipitation occurred in
RI
three different environments: i) phreatic settings, i.e. tens to hundreds of meters below
SC
the water table; ii) epiphreatic environments, i.e. less than tens of meters below the
water table and; iii) long-standing or ephemeral cave pools, located above the local
NU
water table. However, much earlier, and before the groundwater reached calcium
carbonate saturation and speleothems precipitated, karstification
(corrosion/dissolution of the carbonate host rock) was the prevailing mechanism.
MA
The earliest stages of speleogenesis (formation of karstic voids) date back to the
Cambrian-Ordovician (De Waele et al., 2013), when the upwelling of thermal water
D
and oxidation of sulfide ore bodies during repeated cycles of marine to continental
E
environmental changes (Bechstadt and Boni, 1996) led to karstification via sulfuric
PT
acid speleogenesis (Fig. 6A). Active oxidation of polymetallic sulfide ore bodies
occurred because of enhanced circulation of meteoric water inside the system (Fig.
CE
6B); this oxygen-rich water forced intense alteration of sulfides and partial re-
mobilization of ore mineralizations (De Vivo et al. 1987; Cortecci et al. 1989). Typical
hypogene morphologies related to rising thermal fluids, including bubble trails, cupolas
AC
and megascallops (De Waele et al., 2013), probably developed during the more recent
karstification interval since the Oligo-Miocene, up to present.
From the Oligocene to the Pliocene subaerial and subaqueous conditions alternated
in the entire Iglesiente region, with periods of seawater ingression into the karst (De
Vivo et al., 1987). Subsequent stages were marked by water temperature lowering, as
the solution reached the supersaturation in quartz (SiO 2) and calcite. First, quartz
precipitated from highly saline solutions (13.3–23.3 % NaCl eq.) at high temperature
(70–112 °C) (De Vivo et al., 1987). Subsequently, the water temperature decreased
ACCEPTED MANUSCRIPT
PT
spars GA-03-1 and GA-03-2 found in Crovassa Azzurra cave. Most high-temperature
18O, with 18Oc values ranging from -18.8 ‰ to -13.6
RI
spars are relatively depleted in
‰, typical of high temperature carbonates (e.g. LaRock and Cunningham, 1995;
SC
Dublyansky and Spötl, 2009; Palmer and Palmer, 2012; Decker et al., 2015), whereas
their 13C values are within a narrow range from -0.5 ‰ to +0.2 ‰. These carbon
NU
isotope values are similar to the Cambrian carbonate host rock (-0.8 ‰ to -0.2 ‰) and
suggest that the main source of dissolved inorganic carbon during this early stage was
MA
marine bedrock, with insignificant contributions from other sources (e.g. CO2 from the
overlying soil and vegetation) (Fig. 5).
Although the exact age of the early phases of calcite spar precipitation has not been
D
determined, thermal conditions might have occurred much earlier than 600 ka BP (Fig.
E
suggested for the nearby Corona ‘e Sa Craba Cave, ca. 10 km south of Mt. San
Giovanni and at around 270 m a.s.l.. Also there, the age of this episode remains
CE
unknown, although it has been estimated to have occurred during the Lower Miocene
(ca. 20 Ma), when volcanic activity characterized this part of Sardinia (Sauro et al.,
AC
2014).
Thermal activity persisted until at least ~400 ka BP (Fig. 6D), as suggested from
relatively low 18Oc at ~208 m a.s.l. in Crovassa Azzurra cave (~70 oC at 394±13 ka
BP; GA-05-1). However, even before that time (i.e. 684±244 and 570±53 ka BP)
subaerial carbonate flowstones (the base of the GA-core) formed at low temperature
(~12-16 oC) in the upper chamber of this cave, at 216 m a.s.l., where no sign of
hydrothermal calcite precipitation has been found from ca. 600 ka to 83 ka BP. This
testifies that, at ~500-400 ka, carbonate deposition in two sectors of the same cave
occurred in radically different hydrogeochemical environments (Fig. 6D). This can be
ACCEPTED MANUSCRIPT
explained by the upper cave levels were affected by meteoric cold-water infiltration
and were relatively isolated from the lower levels by more than 10 m of carbonate rock.
The age of the cold-water calcite coatings (~13±3 oC at 410±14 ka; sample GA-05-3)
deposited over the GA-05-1 spar, and thus stratigraphically younger than the euhedral
calcite crystals, overlaps with the final stage of high-temperature hydrothermal activity
recorded by this spar. Apparently, a shift from hydrothermal to low-temperature
conditions occurred in the caves at ~400 ka at the ~208 m level. A plausible
PT
explanation for this sudden change is a rapid drop of the local water table and
replacement of thermal water at the ~208 m level (e.g. the bottoms of Crovassa
RI
Azzurra, Massa Riccardo and Santa Barbara 1 caves) by cold water that occupied
perched ponds fed by meteoric seepage (Fig. 6E). The local groundwater level may
SC
be tens of meters below the perched ponds, and although the aquifer water
temperature may still be warm, this did not significantly affect the temperature of the
NU
perched waters that likely adjusted to the average external environmental temperature
during some phases of MIS 11. The coexistence of a thermal aquifer and colder
MA
perched pools at higher levels in the same karst system and disconnected from the
main groundwater level, has been reported from hydrothermal systems elsewhere
(Gázquez et al., 2012; 2013a). The relatively low 13C values found in samples
D
younger than 400 ka (-11.1 ‰ to -4.3 ‰), compared to the older calcite spars support
E
the greater contribution of edaphic CO2, and minor input of CO2 of deep/hypogene
PT
et al., 2010) and carbonate precipitation in perched pools disconnected from a main
hydrothermal reservoir.
AC
After ~400 ka, no evidence is found that the local aquifer reoccupied the ~208 m level.
Interestingly, mammillary calcite coatings (cave clouds) deposited under ephiphreatic
conditions in the lower level of Santa Barbara 2 cave at 247±17 ka (sample SB2-CC;
De Waele and Forti, 2006), probably in a relatively warm aquifer. Cave clouds
generally form in epiphreatic, warm environments, by a slow mechanism of CO2
degassing close to the water table (some tens of meters below). Thus, at ~247 ka BP
the water table probably was around 90 m a.s.l., at the same time that the calcite
“roses” formed in a cold water perched pound in Massa Riccardo cave at 208 m a.s.l..
ACCEPTED MANUSCRIPT
In synthesis, our data suggest that the phreatic level may have dropped from 208 m
a.s.l. at ~400 ka BP, to around 90 m a.s.l. at 247±17 ka BP. The latter altitude is a
conservative estimate, considering that the dated cave cloud in Santa Barbara 2 was
not the highest one in this location (mammillary calcite cave clouds occur up to 10 m
higher) (De Waele and Forti, 2006), and these peculiar speleothems can also form
very close to the water surface (Hill and Polyak, 2010). Thus, the groundwater level
dropped at least 120 m in 147±30 kyrs, which is a rate of 0.82 +0.2/-0.14 mm/year. Similar
values are found in rapidly uplifting areas (Picotti and Pazzaglia, 2008; Columbu et al.,
PT
2015), but Sardinia is thought to be a relatively stable terrain at least for the last million
years. A minimal (<0.1 mm/year) uplift rate locally affected Sardinia (D’Angeli et al.,
RI
2015b) during the Late Quaternary, constituting ~14±3 m absolute groundwater falling
SC
over this 147±30 ka timeframe. Nevertheless, future detailed studies will investigate
other processes triggering the remaining ~100 m drop. At this point of the research,
NU
we exclude higher uplift rates affecting the SW area of Sardinia and we point toward
a prolonged downwearing of the phyllites (Fig. 1) surrounding the karst terrains at the
surface. This low permeability rock, subjected to higher surface dismantling with
MA
perturbations, etc.) might have intensified their weathering and incision, accelerating
the drainage of the aquifer and thus the groundwater drop.
E
PT
6. Conclusions
CE
Clumped isotope thermometry, along with U-Th series dating of different types of
carbonate speleothems from caves of the Mount San Giovanni mine have been used
to reconstruct the temperature, isotope composition (18Ow) and the position of the
AC
groundwater level during the Late Quaternary. Subaqueous calcite spars were already
forming before 600 ka at temperatures above ~120 oC, from a solution with 18Ow
values slightly higher than the modern infiltration water in this region and the paleo-
water 18Ow during the Late Quaternary reconstructed here.
spars formed in perched cave pools at approximately the same time (410±14 ka),
within analytical errors. At 265±6 ka, the aquifer temperature remained cold, whereas
subaqueous dogtooth calcite crystals formed between 216±1 ka and 109±0.5 ka in a
subterranean pool at lower temperatures. Subaerial flowstones continued forming
above 210 m a.s.l., for at least 300 ka (from 383 to 82 ka), stopping their growth only
when climate became dry during MIS 5b.
Our data suggest that the water level fell by ~120 m from ~400 to ~250 ka; which is a
PT
rapid rate of ~0.8 +0.2/-0.14 mm/year, never observed for the tectonically stable Sardinian
landmass. Therefore, we suggest that such a high rate has been the result of climate-
RI
driven surface weathering of less permeable bounding geological formations (i.e.
SC
phyllites) that, prior to ~400 ka, confined the karst reservoir.
The reconstructed 18Ow values of the paleo-aquifer during MIS 11, MIS 8, MIS 7 and
NU
MIS 5c differ by less than 2 ‰ from the modern cave water and meteoric water in the
region. This 18Ow results can be of interest for future quantitative paleo-temperature
MA
researches.
PT
Acknowledgements
CE
We thank James Rolfe and Ian Mather for technical assistance in the Godwin
AC
Laboratory. This study was supported by the NERC project, ID E/M003752/1 to D.A.H.
The sampling surveys were supported by the Water Resources and Environmental
Geology Research Group (University of Almería). S.F.M.B., F.G. and D.A.H.
acknowledge support from the European Union’s Horizon 2020 Research and
Innovation programme for project “QUEST” (Marie Sklodowska-Curie grant
agreement No. 691037). U-Th dating was supported by grants from Taiwan ROC
MOST (104-2119-M-002-003, 105-2119-M-002-001 to C.-C.S.) and the Research
Center for Future Earth, National Taiwan University (107L901001 to C.-C.S.). The
ACCEPTED MANUSCRIPT
authors are grateful to Professor Philippe Audra and two anonymous reviews, who
contributed to improve the original manuscript.
REFERENCES
Affek, H.P., Bar-Matthews, M., Ayalon, A., Matthews, A., Eiler, J.M., 2008.
PT
recorded by ‘clumped isotope’ thermometry. Geochim. Cosmochim. Acta 72,
RI
5351–5360.
Arienzo, M. M., Swart, P. K., Vonhof, H. B., 2013. Measurement of 18O and 2H values
SC
of fluid inclusion water in speleothems using cavity ring down spectroscopy
NU
compared with isotope ratio mass spectrometry, Rapid Comm. Mass Spectrom.
27, 2616–2624.
MA
Audra, P., Bigot, J.Y., Mocochain, L., 2002. Hypogenic caves in Provence (France).
Audra, P., Hobléa, F., Bigot, J.Y., Nobécourt, J.C., 2007. The role of condensation in
E
Audra, P., Mocochain, L., Bigot, J.Y., Nobécourt, J.C., 2009. The association between
Audra, P., Gázquez, F., Rull, F., Bigot, J.Y., Camus, H., 2015. Hypogene Sulfuric Acid
Bar-Matthews, M., Marean, C.W., Jacobs, Z., Karkanas, P., Fisher, E.C., Herries,
A.I.R., Brown, K., Williams, H.M., Bernatchez, J., Ayalon, A., Nilssen, P.J., 2010.
Bechstadt, T., Boni, M., (eds), 1996. Sedimentological, stratigraphical and ore deposits
PT
field guide of the autochtonous Cambro-Ordovician of Southwestern Sardinia.
RI
Memorie descrittive della Carta Geologica d’Italia 48. Servizio Geologico Italiano,
SC
Breitenbach, S.F.M., Mleneck-Vautravers, M.J., Grauel, A.-L., Lo, L., Bernasconi, S.M.,
NU
Müller, I.A., Rolfe, J., Gázquez, F., Greaves, M., Hodell, D.A., (2018). Coupled
nel quadro della geologia del Cambrico sardo. L’Industria Mineraria 19, 1-53.
E
Budd, D.A., Frost III, E.L., Huntington, K.W., Allwardt, P.F., 2013. Syndepositional
PT
Caddeo, G.A., De Waele, J., Frau, F., Railsback, L.B., 2011. Trace element and stable
AC
Calaforra, J.M., De Waele, J., 2011. New peculiar cave ceiling forms from Carlsbad
134, 43-48.
ACCEPTED MANUSCRIPT
Cheng, H., Edwards, R.L., Hoff, J., Gallup, C.D., Richards, D.A., Asmerom, Y., 2000.
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., Alexander, E.C., 2013.
Improvements in 230Th dating, 230Th and 234U half-life values, and U-Th isotopic
PT
Earth Planet. Sci. Lett. 371-372, 82-91.
RI
Chiesi, M., 2005. Il pre-monitoraggio dei parametri ambientali della Grotta di Santa
SC
parameters in the Santa Barbara Cave), Memorie dell’Istituto Italiano di
NU
Speleologia, 2(17), 23-34.
Cidu, R., Biddau, R., Fanfani, L., 2009. Impact of past mining activity on the quality of
MA
Cidu, R., Frau, F., Da Pelo, S., 2011. Drainage at abandoned mine sites: natural
D
Civita, M., Cocozza, T., Forti, P., Perna, G., Turi, B., 1983. Idrogeologia del bacino
CE
Clark, I.D, Fritz, P., 1997. Environmental Isotopes in Hydrogeology. Lewis Publishers.
USA.
Columbu, A., De Waele, J., Forti, P., Montagna, P., Picotti, V., Pons-Branchu, E.,
Hellstrom, J., Bajo, P., and Drysdale, R. 2015. Gypsum caves as indicators of
43, 539-542.
ACCEPTED MANUSCRIPT
Columbu, A., Drysdale, R., Capron, E., Woodhead, J., De Waele, J., Sanna, L.,
Hellstrom, J. and Bajo, P., 2017. Early last glacial intra-interstadial climate
3948–3957.
PT
Cortecci, G., Fontes, J. C., Maiorani, A., Perna, G., Pintus, E., Turi, B., 1989. Oxygen,
RI
sulfur, and strontium isotope and fluid inclusion studies of barite deposits from
SC
Deposita 24(1), 34-42.
NU
Daëron, M., Blamart, D., Peral, M., Affek, H.P., 2016. Absolute isotopic abundance
ratios and the accuracy of 47 measurements. Chem. Geol. 442, 83-96.
MA
D’Angeli, I.M., De Waele, J., Ceballo Melendres, O., Tisato, N., Sauro, F., Grau-
CE
Gonzalez, E., Bernasconi, S.M., Torriani, S., Bontognali, T.R.R., 2015a. Genesis
526-535.
D'Angeli, I. M., Sanna, L., Calzoni, C., De Waele, J., 2015b. Uplifted flank margin
Decker, D.D., Polyak, V.J., Asmerom, Y., 2015. Depth and timing of calcite spar and
“spar cave” genesis: Implications for landscape evolution studies. In: Feinberg,
ACCEPTED MANUSCRIPT
J., Gao, Y., Alexander, E.C., Jr. (Eds.), Caves and Karst Across Time. Geological
De Vivo, B., Maiorani, A., Perna, G., Turi, B., 1987. Fluid inclusion and stable isotope
studies of Calcite, Quartz and Barite from karstic caves in the Masua mine,
De Waele, J., Forti, P., 2006. A new hypogean karst form: the oxidation vent. Z.
PT
Geomorphol. 147, 107-127.
RI
De Waele, J., Forti, P., Naseddu, A., 2013. Speleogenesis of an exhumed
SC
Sardinia). Earth Surf. Proc. Land. 38(12), 1369-1379.
NU
De Waele, J., Audra, P., Madonia, G., Vattano, M., Plan, L., D’Angeli, I.M., Bigot, J.-Y.,
Nobécourt, J.-C., 2016. Sulfuric acid speleogenesis (SAS) close to the water
MA
table: Examples from southern France, Austria, and Sicily. Geomorphology 253,
52-467.
D
De Waele, J., Gázquez, F., Forti, P., Naseddu, A., 2017. Inactive Hydrothermal
E
Hypogenic Karst in SW Sardinia (Italy), In: Klimchouk, A., Palmer, A., De Waele,
PT
J., Auler, A.S., Audra, P. (Eds.), Hypogene Karst Regions and Caves of the World,
CE
Simferopol, 45-49.
ACCEPTED MANUSCRIPT
Edwards, R.L., Chen, J.H., Ku, T.-L., Wasserburg, G.J., 1987. Precise timing of the
Fairchild, I. J., Smith, C. L., Baker, A., Fuller, L., Spötl, C., Mattey, D., McDermott, F.,
PT
E.I.M.F., 2006. Modification and preservation of environmental signals in
RI
speleothems. Earth Sci. Rev. 75, 105-153.
Ferry J.M., Passey B.H., Vasconcelos C., Eiler J.M., 2011. Formation of dolomite at
SC
40-80°C in the Latemar carbonate buildup, Dolomites, Italy, from clumped
NU
isotope thermometry. Geology 39, 571–574.
Ford, D., Williams, P., 2007. Karst geomorphology and hydrology. John Wiley & Sons,
MA
Chichester, 562 p.
Fornós, J.J., Gelabert, B., Ginés, A., Ginés, J., Tuccimei, P., Vesica, P., 2002. Phreatic
D
Forti, P., 2010. Genesis and evolution of the caves in the Naica mine (Chihuahua,
Gázquez, F., Calaforra, J.M., 2013. Origin of double-tower raft cones in hypogenic
Gázquez, F., Calaforra, J.M., Forti, P., Rull, F. and Martínez- Frías, J. 2012. Gypsum-
220.
Gázquez, F., Calaforra, J.M., Stoll, H., Sanna, L., Forti, P., Lauritzen, S.E., Delgado,
A., Rull, F. and Martínez-Frías, J.M. 2013a. Isotope and trace element evolution
PT
of the Naica aquifer (Chihuahua,Mexico) over the past 60,000 yr revealed by
RI
speleothems. Quat. Res. 80, 510 -521.
Gázquez, F., Calaforra, J.M., Forti, P., De Waele, J., Sanna, L., Rull, F., Sanz, A.,
SC
2013b. Corrosion of calcite crystals by metal-rich mud in caves: study case in
NU
Crovassa Ricchi in Argento Cave (SW Sardinia, Italy). Geomorphology 198, 138–
146.
MA
Gázquez, F., Calaforra, J.M., Rodriguez-Estrella, T., Ros, A., Lamusi, J.L, Sanchez, J.,
2017. Evidence for Regional Hypogene Speleogenesis in Murcia (SE Spain) In:
D
Klimchouk A., Palmer A., De Waele, J. Auler, A.S., Audra, P. (Eds.), Hypogene
E
Karst Regions and Caves of the World, Cave and Karst Systems of the World,
PT
Springer, 183-197.
CE
Ghosh, P., Patecki, M., Rothe, M., Brand, W. A., 2006. 13C–18O bonds in carbonate
1439–1456.
Hellstrom, J., 2003, Rapid and accurate U/Th dating using parallel ion-counting multi-
Hill, C.A., Forti P., 1997. Cave minerals of the World. 2nd edition. Nat. Speleol. Soc.,
Huntsville, 463 p.
ACCEPTED MANUSCRIPT
Hill., C.A., Polyak, V.J., 2010. Karst hydrology of Grand Canyon, Arizona, USA. J.
Jones, B., 1989. Calcite Rafts, Peloids, And Micrite In Cave Deposits From Cayman
Kele, S., Breitenbach, S.F.M., Capezzuoli, E., Nele Meckler, A., Ziegler, M., Millan,
PT
I.M., Kluge, T., Deák, J., Hanselmann, K., John, C.M., Yan, H., Liu, Z.,
RI
Bernasconi, S.M. (2015) Temperature dependence of oxygen- and clumped
SC
temperature range, Geochim.Cosmochim. Acta. 168, 172-192.
NU
Kim S.-T., O’Neil J. R., Hillaire-Marcel C. & Mucci A., 2007 - Oxygen isotope
perspective. National Cave and Karst Research Institute Special Paper #1,
E
Carlsbad, 106 p.
PT
100-117.
Klimchouk, A.B., Palmer, A., De Waele, J., Auler, A.S., Audra, P., 2017. Hypogene
Karst Regions and Caves of the World, Cave and Karst Systems of the World,
Kluge, T., Affek, H.P., Zhang, Y., Dublyansky, Y., Spötl, C., Immenhauser, A., Richter,
Kluge, T., Affek, H.P., Dublyansky, Y., Spötl, C., 2014b. Devils Hole
PT
Kovacs, S.E., Reinhardt, E.G., Chatters, J.C., Rissolo, D., Schwarz, H.P., Collins,
RI
S.V., Kim, S.-T., Blank, A.N., Erreguerena, P.L., 2017. Calcite raft geochemistry
SC
as a hydrological proxy for Holocene aquifer conditions in Hoyo Negro and Ich
Balam (Sac Actun Cave System), Quintana Roo, Mexico. Quat. Sci. Rev. 175,
NU
97-111.
LaRock, E.J., Cunningham, K.I., 1995. Helictite bush formation and aquifer cooling in
MA
Wind Cave, Wind Cave National Park, South Dakota. J. Cave Karst Stud. 57(1),
43-51.
D
Leél-Őssy, S., Szanyi, G., Surányi, G., 2011. Minerals and Speleothems of the József-
E
Lundberg, J., Ford, D.C., Hill, C.A., 2000. A preliminary U-Pb date on cave spar, Big
CE
Canyon, Guadalupe Mountains, New Mexico. Journal of Cave and Karst Studies
62, 144–148.
AC
Meckler, A.-N., Ziegler, M., Breitenbach, S.F.M., Millan, I., Bernasconi, S.M., 2014.
scheme for clumped isotope measurements. Rapid Comm. Mass Spectrom. 28,
1705-1715.
Messina, M., Naseddu, A., Papinuto, S., Sanna, F., Sotgia, S., Forti, P., De Waele, J.,
69-86.
Pagliara, A., De Waele, J., Forti, P., Galli, E., Rossi, A., 2010. Speleothems and
PT
Palmer, A.N., 2011. Distinction between epigenic and hypogenic caves.
RI
Geomorphology 134, 9-22
SC
Palmer, M.V., Palmer A.N. 2012. Petrographic and isotopic evidence for late-stage
processes in sulfuric acid caves of the Guadalupe Mountains, New Mexico, USA.
NU
Int. J. Speleol. 41(2), 231–250.
Piccini, L., 2000. Il carsismo di origine idrotermale del Colle di Monsummano (Pistoia–
MA
Piccini, L., De Waele, J., Galli, E., Polyak, V.J., Bernasconi, S.M., Asmerom, Y., 2015.
E
Picotti, V., Pazzaglia, F.J., 2008. A new active tectonic model for the construction of
the Northern Apennines mountain front near Bologna (Italy). J. Geophys. Res.
AC
113, B005307.
Plan, L., De Waele, J., 2011. Folia in der Odelsteinhöhle, Steiermark. Erstnachweis im
Plan, L., Tschegg, C., De Waele, J., Spötl, C., 2012. Corrosion morphology and cave
Polyak, V.J., DuChene, H.R., Davis, D.G., Palmer, A.N., Palmer, M.V., Asmerom Y.,
2013. Incision history of Glenwood Canyon, Colorado, USA, from the uranium-
PT
Rossetti, V., Zucchini, A., 1956. Baritina della Grotta di Santa Barbara. Rendiconti del
RI
Seminario della Facoltà di Scienze Università di Cagliari 26(3-4), 240- 255.
Sauro, F., De Waele, J., Onac, B.P., Galli, E., Dublyansky, Y., Baldoni, E., Sanna, L.,
SC
2014. Hypogenic speleogenesis in quartzite: The case of Corona ‘e Sa Craba
NU
Cave (SW Sardinia, Italy). Geomorphology 211, 77–88.
Shen, C.-C., Wu, C.-C., Cheng, H., Edwards, R.L., Hsieh, Y.-T., Gallet, S., Chang, C.-
MA
C., Li, T.-Y., Lam, D.D., Kano, A., Hori, M., Spötl, C., 2012. High-precision and
Spötl, C., Vennemann, T.W., 2003. Continuous-flow isotope ratio mass spectrometric
PT
1006.
Stafford, K. W., Raymond, N., Rosales-Lagarde, L., Boston, P.J., 2008. Epigene and
AC
New Mexico and Culberson County, Texas, USA. Int. J. Speleol. 37, 83-98.
Talma, A.S., Vogel, J.C., 1992. Late Quaternary paleotemperatures derived from a
speleothem from Cango Caves, Cape Province, South Africa. Quat. Res. 37,
203–213.
Urey, H.C., 1948. Oxygen isotopes in nature and in the laboratory. Science 108, 489.
ACCEPTED MANUSCRIPT
Karst Regions and Caves of Iran. In: Klimchouk, A.B., Palmer, A.N., De Waele,
J., Auler, A.S., Audra, P. (Eds.), Hypogene Karst Regions and Caves of the World,
PT
Acte 8º Congress National of the Speleological Society, Switzerland, 137-142.
RI
SC
FIGURES AND TABLE CAPTIONS
Figure 1. Geological map of the Iglesiente-Sulcis mining district (SW Sardinia, Italy)
NU
and location of the Mount San Giovanni mine (after Caddeo et al., 2011).
Figure 2. Cross-section of the Mount S. Giovanni mine (after De Waele et al., 2013)
MA
with the locations of the studied caves, the main ore bodies and the mine galleries
and shafts indicated in grey. Sampled caves in this study are shown in red.
D
sites: GA-05 (calcite spar and calcite coatings); GA-01 (calcite dogtooth crystals); GA-
PT
core (flowstone core); GA-03 (calcite spars); GA-BA (blue aragonite flowstone); GRA-
01 (calcite spars); MS-01 (calcite roses); MS-06 (calcite spar and flowstone); GA-host
CE
and MS-host (host rock); SB2-CC. Mammillary calcite coatings (cave clouds) analyzed
by De Waele and Forti (2006). Note that the relative positions of the caves in the X
AC
the top and MS-06-2 from the base) collected in Massa Riccardo cave. E and F.
Sampling of the flowstone core in Crovassa Azzurra cave (sample GA-core) and
position of eleven subsamples for U-Th dating.
Fig. 5. Oxygen and carbon isotopes in carbonate samples from the caves of the Mount
San Giovanni mine.
Figure 6. Genetic evolution of the Mount San Giovanni caves. Estimated water table
levels and temperatures are given in connection with precipitation of different types of
PT
speleothems: A. First stages of hypogene speleogenesis at high temperature during
Cambrian; B. Intense karstification during the Carboniferous, probably related to
RI
H2SO4 produced due to oxidation of polymetallic sulfide ore bodies after increased O 2-
SC
rich meteoric water input; C. The first stages of hydrothermal calcite precipitation took
place before 600 ka BP at temperature ~120 oC; D. Precipitation of hydrothermal (~70
NU
oC) calcite spars and subaerial speleothem in the upper level of the karst coexisted
between ~500 and ~400 ka BP; E. At ~400 ka BP, the water temperature decreased
(~14 oC) and subaqueous calcite coatings formed in a perched pound covering
MA
hydrothermal spars formed during previous stages in Crovassa Azzurra cave. At the
same time, subaerial aragonite flowstones formed in the upper level of this cave; F.
D
Different types of epiphreatic speleothems (calcite dogtooth crystals and calcite roses)
precipitated at low temperature (~10 oC to 20 oC) in perched pounds between 265 and
E
109 ka BP at the ~208 m level, while mammillary calcite coatings formed in deeper
PT
levels (~65 m a.s.l.) from the still warm, aquifer water. Meanwhile, aragonite flowstone
kept growing until 82 ka BP. Note that the depth scale is illustrative and may have
CE
Table 1. Sample locations and descriptions, isotopic analyses, U-Th ages (±2SD) (see
supplementary Table S1 for details) and temperature estimates using clumped
isotopes thermometry (T47 ± 2SE). The number of aliquots analyzed (n) is given. The
actual 18Ow (± 2SE) of the paleo-aquifer has been calculated from the T47 and 18Oc
using equation 1 in the main text (see supplementary Table S2 for calculations and
propagation of analytical uncertainties). Water temperatures are also calculated using
the 18Oc-18Ow approach (T18Oc-18Ow), assuming that the paleo-18Ow was -5.5±1 ‰
(2SE) for all samples, a value similar to the modern cave water (see discussion for
ACCEPTED MANUSCRIPT
details), using eq. 1 for calcite samples and the paleotemperature equation of Kim et
al. (2007) for aragonite samples. The 18Oc-18Ow approach gives temperature
uncertainties of ±5oC (2SE). *sample analyzed by De Waele and Forti (2006).
PT
RI
SC
NU
MA
E D
PT
CE
AC
ACCEPTED MANUSCRIPT
P T 20
16
0.680±0.013 (n=16)
0.694±0.013 (n=20)
17±3
14±3
-6.0±0.7
-6.0±0.7
I
GA-01-3 201 Dogtooth calcite crystals Subaqueous (low temperature) -4.2 -6.2 - 16 - - -
GA-BA-1 208 Aragonite flowstone (white layer) Subaerial (low temperature) -4.3 -5.7 - 12 - - -
GA-BA-2
GA-03-1
GA-03-2
208
205
205
Aragonite flowstone (blue layer)
Scalenohedral calcite spar
Scalenohedral calcite spar
Subaerial (low temperature)
Subaqueous (high temperature)
Subaqueous (high temperature)
-4.1
-18.8
-15.6
-4.3
-0.3
C
0.1 R
383±10
-
-
11
109
83
-
0.476±0.018 (n=20)
-
-
130±14
-
-
-3.1±1.5
-
S
GA-core-2 216 Flowstone core (36 cm from the top) Subaerial (low temperature) -5.1 -8.9 684±244 16 - - -
GA-core-3 216 Flowstone core (31.5 cm from the top) Subaerial (low temperature) -4.8 -11.1 >600 14 - - -
GA-core-4
GA-core-11
GA-core-12
216
216
216
Flowstone core (27 cm from the top)
Flowstone core (22.5 cm from the top)
Flowstone core (21.5 cm from the top)
Subaerial (low temperature)
Subaerial (low temperature)
Subaerial (low temperature)
-4.3
-4.5
N
-4.6 U -9.1
-8.2
-6.2
570±53
587±44
581±51
12
14
15
-
-
-
-
-
-
-
-
-
A
GA-core-5 216 Flowstone core (16.6 cm from the top) Subaerial (low temperature) -4.1 -8.4 506±19 11 - - -
GA-core-6 216 Flowstone core (12 cm from the top)-aragonite Subaerial (low temperature) -4.8 -7.6 228±2 14 - - -
M
GA-core-7 216 Flowstone core (8.5 cm from the top)-aragonite Subaerial (low temperature) -4.4 -7.4 208±1 12 - - -
GA-core-8 216 Flowstone core (6 cm from the top)-aragonite Subaerial (low temperature) -5.2 -7.2 198±1 16 - - -
GA-core-9 216 Flowstone core (3.6 cm from the top)-aragonite Subaerial (low temperature) -4.9 -6.8 193.6±0.8 14 - - -
GA-core-10
GA-host
MASSA RICCARDO
216
210
Flowstone core (0.5 cm from the top)
Host rock
E D
Subaerial (low temperature)
-
-4.5
-10.4
-6.5
-0.8
82.3±0.3
-
13
-
-
-
-
-
-
-
MS-06-1
MS-06-2
MS-06-3
206
206
206
Calcite flowstone (top)
Calcite flowstone (bottom)
Scalenohedral calcite spar
P T Subaerial (low temperature)
Subaerial (low temperature)
Subaqueous (medium temperature)
-4.4
-5.1
-7.0
-10.1
-10.2
-8.4
-
-
-
16
20
29
-
-
-
-
-
-
-
-
-
MS-01-1
MS-01-2
MS-host
208
208
207
Calcite “roses”
Calcite “roses”
Host rock
C E Subaqueous (low temperature)
Subaqueous (low temperature)
-
-4.4
-5.2
-11.7
-10.4
-10.4
-0.2
-
265±6
-
16
23
-
-
0.708±0.008 (n=37)
-
-
10±2
-
-
-7.7±0.4
-
CROVASSA RICCHI IN ARGENTO
GRA-01-1
GRA-01-2
GRA-01-3
GRA-01-4
180
180
180
180
A C
Scalenohedral calcite spar
Scalenohedral calcite spar
Scalenohedral calcite spar
Scalenohedral calcite spar
Subaqueous (high temperature)
Subaqueous (high temperature)
Subaqueous (high temperature)
Subaqueous (high temperature)
-18.8
-18.2
-18.8
-18.6
-0.3
0.2
0.1
0.2
>600
-
-
-
109
105
109
109
0.485±0.009 (n=20)
-
0.491±0.012 (n=18)
-
123±6
-
119±9
-
-3.8±0.6
-
-4.3±1.0
-
SANTA BARBARA 2 CAVE
SB2-CC* 65 Mammillary calcite coatings Subaqueous (medium temperature?) - - 247±17 - - - -
Table 1.
ACCEPTED MANUSCRIPT
Highlights:
We apply clumped isotopes on subaqueous speleothems to reconstruct changes in paleo-
aquifer water temperature and 18Ow.
The 18Ow of the paleo-water is then reconstructed from the clumped isotope temperature and
the 18O of the carbonate.
Hydrothermal spars (~120°C) formed before 400 ka, while low-temperature (10 to 17°C)
speleothems did from 400 to 82 ka.
The 18O values of the paleo-aquifer ranged from -6.0±0.7 ‰ to -6.4±0.7 ‰ during MIS 5c, MIS7
PT
and MIS11, similar to modern values.
The 18Ow of the paleo-aquifer was ~2‰ lower during the colder MIS 8.
RI
SC
NU
MA
E D
PT
CE
AC
Figure 1
Figure 2
Figure 3
Figure 4
Figure 5
Figure 6