Next Article in Journal
Climate Change Made Major Contributions to Soil Water Storage Decline in the Southwestern US during 2003–2014
Previous Article in Journal
Assessing the Fresh–Saline Groundwater Distribution in the Nile Delta Aquifer Using a 3D Variable-Density Groundwater Flow Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Reconstruction of Thermal Conditions in the Subboreal Inferred from Isotopic Studies of Groundwater and Calcareous Tufa from the Spring Mire Cupola in Wardzyń (Central Poland)

1
Institute of Hydrogeology and Engineering Geology, University of Warsaw, Żwirki i Wigury St. 93, 02-089 Warszawa, Poland
2
Institute of Geological Sciences, Polish Academy of Sciences, Twarda St. 51/55, 00-818 Warszawa, Poland
3
Laboratory of Geology, University of Łódź, Kopcińskiego St. 31, 90-142 Łódź, Poland
*
Author to whom correspondence should be addressed.
Water 2019, 11(9), 1945; https://doi.org/10.3390/w11091945
Submission received: 14 August 2019 / Revised: 4 September 2019 / Accepted: 12 September 2019 / Published: 19 September 2019
(This article belongs to the Section Hydrology)

Abstract

:
Studies with application of stable isotopes of oxygen and carbon have been performed on calcareous tufa, groundwater and dissolved inorganic carbon (DIC) from the spring mire cupola in Wardzyń. This study was focused on the verification of the a priori hypothesis that the analysed calcareous tufa is a chemical deposit and on the attempt to supplement an earlier scenario of environmental changes in the Subboreal with oscillations of water temperature. The constructed model of chemical and isotope balance, and δ13C determinations in DIC, allowed for calculating ratios of stable isotopes of carbon in particular speciations and in gaseous CO2. The obtained results coupled with δ13C values in calcite indicate that this mineral precipitated from the solution chemically (without the contribution of living organisms). Additionally, it was possible to reconstruct the temperature range at which the calcareous tufa was formed. The reconstructed scenario of changes in the thermal conditions was refined based on δ18O determinations in groundwater and calcite. Accordingly, the oldest calcareous tufa, with an age of about 5500 cal years BP, was formed in cool climate conditions (with average annual temperatures by about 3 °C lower than presently). The formation of younger series of the calcareous tufa took place between 4400–2900 cal years BP and represents a much warmer period with two distinct cooler episodes at 3900 and 3000 cal years BP, respectively. The course of the obtained temperature curves correlates well with the GISP2 curve and curves obtained for other sites in Northern, and Central Europe.

1. Introduction

Climate changes over the Holocene are a field of interest for researchers representing various branches of science from the beginning of the 20th century. These studies are of particular importance as they help to better understand the mechanism of recent climate change and can provide important premises for further changes forecasting. Interpretation of changes in climate parameters, in particular air temperature, is possible thanks to the analysis of marine and terrestrial sediments [1] and ice cores extracted from glaciers [2,3]. There is a number of indicators of climate change in deposits such as pollen (e.g. [4]) or macrofossils found in peat (e.g. [5]). An important carrier of information about climate changes are stable carbon isotopes δ13C and δ18O in endogenous calcite (tufa deposits), which were used among others to reconstruct Holocene climatic changes in southern Poland [6] and in the Bohemian Karst [7]. The isotopic composition of fine-grained calcite from lake sediments were used to reflecting variations in the precipitation seasonality in the Baltic region as evidence of changes in North Atlantic ocean-atmospheric circulation [8]. However, apart from information on the variability of climatic conditions, it is necessarily for palaeoclimatic reconstruction to refer registered changes to the time axis. The above conditions were fulfilled for the sediment series of mire spring cupola in Wardzyń.
The spring mire cupola in Wardzyń represents an accumulation of hydrogenic and chemical deposits, genetically related with a set of groundwater springs located in the marginal zone of the Wolbórka valley. Previous studies of the cupola [9] have allowed for distinguishing four levels of calcareous tufa, predominantly composed of calcite separated by peat layers. Radiocarbon dates of the peat indicate that the carbonate horizons were formed from about 5500 cal years BP to about 2900 cal years BP, and the lithological features and mineralogical composition of the deposits record changes of environmental conditions in the Subboreal. The succession of strata and the mineral composition of the deposits has allowed for linking the calcareous tufa formation with dry episodes, with the peat layers representing intervals with elevated annual precipitation [9]. This study applies determinations of stable isotopes of carbon and oxygen in calcite, coupled with 18O/16O ratios determined in groundwater and 13C/12C ratios determined in DIC. Calcite is the most commonly used carbonate in paleoenvironmental reconstructions. It precipitates in the aquatic environment as a result of carbonate imbalance. Calcite precipitation is accompanied by the fractionation of stable carbon and oxygen isotopes. The main factor determining the isotope fractionation (especially in the case of oxygen) is the water temperature. On this basis, by determining stable isotopes in calcite, one can reconstruct how the temperature of the water changed during the precipitation of this mineral from the aqueous solution.
The δ13C composition of various carbon pools is a tracer for the origin of carbon and contributes to understanding the global carbon cycle in the context of the climate change and the global worming effects. The δ18O of carbonates preserves a record of the 18O/16O composition of the water from which the calcite is precipitated and allows to determine the temperature of precipitation process; this led to formulation of the δ18OCaCO3 paleotemperature scale [10].
Terrestrial or freshwater carbonates, including spelothem, travertine, calcrete and secondary mineralization in soli and near surface fractures, precipitate from groundwaters at or near the point of discharge. The geochemical processes that initiate calcite precipitation can include depressurization and CO2 degassing (e.g. spelothem and travertine), evaporation (e.g. soil calcrete), freezing of bicarbonate waters, sulfate reduction and other biotic reactions [11]. Their C and O isotopic composition is a reflection of the geochemical and climatic (i.e. mainly temperature) conditions during precipitation, and can be used to interpret their origin.
In fact, there are two temperatures recorded by the δ18O of freshwater calcite: (i) the water temperature at which precipitation took place—i.e. the temperature of equilibrium fractionation between water and calcite, and (ii) the climatic temperature (i.e. mean annual air temperature) recorded by the δ18O of water according to the local T-δ18O relationship [11].
The main focus of these investigations was to verify the a priori hypothesis that the analysed calcareous tufa is a chemical deposit and to supplement an earlier scenario of environmental changes in the Subboreal with oscillations of water temperature determined in the spring zone.

2. Materials and Methods

2.1. Study Area

The spring mire cupola in Wardzyń is located along the margin of the moraine plateau (Romanowskie Hills) within a depression assigned to the polygenic valley of upper Wolbórka (Figure 1). The cupola covers an area of about 1.5 ha. At present, 20 permanent springs function within it; their total discharge varies within 10.5–16.5 l/s. The spring deposits infill an erosional incision, which is an elongation of a small dry valley visible in the plateau morphology. The cupola is composed of peats with 1 to a maximum of 4 layers of calcareous tufa [9]. These deposits with a total thickness reaching 3.5 m lie directly on medium-and multi-grained fluvial sands (Figure 2). The oldest calcareous tufa (series A) was probably formed at the end of the Atlantic period. Determination of its precise date is not possible due to the lack of organic deposits in its base, however the radiocarbon age of peat covering this calcareous tufa indicates its formation before 4600 ± 60 years BP. The tufa is composed of calcite encrusting plant remains. Sediments of series A occur only within the erosional incision. Above lies a layer of strongly decomposed peat with high contribution of calcium carbonate and detritic material (mainly quartz). Locally these deposits pass laterally into mineral deposits developed as silts. Younger series of the calcareous tufa occur in a narrow, elongated zone, whose orientation corresponds to the axis of the erosional incision of the substratum. The largest thicknesses (locally exceeding 1 m) are attained by these deposits in the southern and central part of the analysed cupola. Three series of calcareous tufa (B, C and D) can be distinguished here, being separated by peat layers. In the north they represent usually over ten centimetres thick interbeds of calcium carbonate within peats. In the marginal zones of calcareous tufa occurrence and in the topmost part of the cupola the deposits represent rhythmites. The tufas are covered by a layer of reed-wood peat, strongly decomposed in its topmost part.
Groundwater springs in the vicinity of the cupola are hydrodynamically related with a shallow aquifer system comprising Cenozoic and the topmost part of Mesozoic formations. The Cenozoic formation was formed as a result of processes related with the advance and retreat of the Scandinavian ice-sheets, and fluvial and lacustrine accumulation in the interglacials. The variable origin of Quaternary deposits and the variable energy of the sedimentary environment have caused that this formation is highly heterogeneous. The aquifer horizons attain diverse thicknesses and may locally be discontinuous. The filtration parameters of the deposits also vary in space. The richest aquifer horizons are linked with fluvioglacial sediments, and fluvial and lacustrine sands, which according to the stratigraphic scheme by Marks et al. [12] represent the complex of Middle Polish Glaciations. These deposits form an aquifer attaining a thickness of up to 35 m. The medium is siliceous and the material includes also numerous grains of carbonate rocks (mainly Palaeozoic limestones). Locally the horizon is bipartite, and the subdividing layer includes glacial tills and ice-dammed deposits. The aquifer sediments are characterised by good permeability. The values of hydraulic conductivity determined on the basis of pumping tests are within 2.7x10−5 to 3.0x10−4 m/s [13]. The lower part of the Quaternary formation represented by deposits of South Polish Glaciations is dominated by poorly permeable deposits (ice-dammed clays and glacial tills). These deposits together with the Neogene weathering cover form the aquitard. Below occur aquifer deposits of the Mesozoic formation assigned to the Łódź Trough. Wardzyń is located in the southern part of the small Łódź Trough–a syncline with a SSW-NNW-oriented axis. Faults with an orientation perpendicular to the syncline axis occur within it, resulting in a block structure of the trough. Aquifer is linked with Upper Cretaceous limestones, marls and opokas. These deposits are usually characterised by good permeability, and the hydraulic conductivities within the sub-Cenozoic basement are between 2.0x10−5 and 1.8x10−4 m/s [13].

2.2. Sampling and Analytical Methods

Water samples were collected from wells (samples 1–10) during their continuous exploitation. The water sample from the Wardzyń spring (11) was collected directly from the outflow from the peat zone, sample (12) was collected from a negative pressure probe with a 0.3 m long active teflon bit. Water temperature, pH, redox potential, and SEC were measured directly in the field; the content of free CO2 was determined by sodium hydroxide titration against phenylophtaleine. Water samples for chemical analyses were collected in single PE bottles, filtered through a 0.45 μm pore size membrane filter (MCE), and preserved in HNO3, if necessary. Analysis of all water samples (Ca2+, Mg2+, Na+, K+) involved flame absorption spectrophotometry SpectrAA 300 (Varian, Palo Alto, CA, USA) (detection limit is 0. 05 mg/L) and UV-vis spectrophotometry S.330 (Marcel, Poland) in the case of NH4+, NO2 (d.l. is 0. 001 mg/L) and NO3 (d.l. is 0,01 mg/L). SO42− was determined by the gravimetric method (PN-C-04566-09), Cl by Mohr’s method (PN-ISO 9297), and alkalinity according to the PN-ISO 9963-1 method (detection limit is 0.1 mg/L). The chemical analyses were performed in the Laboratory of the Department of Geology at the Faculty of Geographical Sciences, University of Łódź and the Laboratory of Computer and Analytical Techniques at the Faculty of Biology and Environmental Protection, University of Łódź.
Ten samples of groundwater from springs and shallow dug wells within the study area were taken for determination of the isotopic composition of oxygen and hydrogen. All waters were collected in 30 ml glass bottles to overflowing and tightly capped to avoid air bellows inside. The 18O/16O ratios in waters were determined by the CO2–H2O equilibration technique with application of on-line continuous flow Gas Bench II system supported by Combi PAL autosampler (CTC Analytics, Zwingen, Switzerland) and connected to MAT 253 (Thermo Fisher, Germany) isotope ratio mass spectrometer (IRMS). The equilibration process was performed for 18 h at 32 °C. The 2H/1H ratios in water samples were determined with application of H-Device automated system equipped with Combi PAL autosampler and connected to MAT 253 IRMS. For mass spectrometric analysis of the 2H/1H ratio the hydrogen gas was obtained via direct reduction of water on hot chromium in temperature 850 °C. The results were reported using δ notation with respect to the VSMOW international standard. Normalization of measured data was performed according to three international standards, namely: VSMOW (δ2H = 0.0‰, δ18O = 0.0‰), VSLAP (δ2H = −427.5‰,δ18O = −55.50‰), and GISP (δ2 = −189.5‰, δ18O = −24.76‰). The precision of the δ18O and δ2H measurements were ±0.1‰ and ±1.0‰ respectively.
Twelve samples of water were collected from the same springs and shallow wells for determination of stable carbon isotopic composition of DIC (dissolved inorganic carbon). All waters were collected in 100 ml glass bottles to overflowing, tightly capped to avoid air bellows inside, and additionally capped with parafilm to avoid evaporation and gas exchange during transport. The 13C/12C ratios of DIC in water samples were determined on the evolved CO2 gas with application of the on-line continuous flow Gas Bench II system equipped with Combi PAL autosampler and connected to MAT 253 IRMS. CO2 gas was evolved from DIC using >98% anhydrous H3PO4 acid and left for 18h equilibration at 70 °C prior to the isotopic analysis.
Thirteen samples of carbonates (i.e. >98% calcite) were collected from carbonate domes within the study area for determination of their carbon and oxygen isotope composition. The 13C/12C and 18O/16O ratios of calcites were determined on the released CO2 gas with application of the same on-line continuous flow Gas Bench II system connected to MAT 253 IRMS. CO2 gas was released into the vials headspace from solid calcites reacting with concentrated (i.e. >98%), anhydrous H3PO4 acid at 70 °C for 1h. The results of measurements of δ13C of DIC, and δ13C and δ18O of the carbonates were reported with respect to VPDB international standard. Normalization of measured data was performed according to two international standards, namely: NBS-18 (carbonatite: δ13C = −5.01‰, δ18O = −23.01‰) and NBS-19 (limestone: δ13C = +1.95‰, δ18O = −2.2‰). The precision of the δ13C and δ18O analysis were ± 0.1‰. The precision of the isotopic measurements was calculated based on long-term measurements of international reference materials.
The isotopic analysis was performed in the Institute of Geological Sciences of the Polish Academy of Sciences, Warsaw, Poland.
Six peat samples, 50 g each, were collected from two cores for the radiocarbon analysis. Determinations were made in the Absolute Age Laboratory at Skała. The scintillation technique with the use of a low-background spectrometer was applied, in which benzene obtained from the material containing organic carbon is analysed. In earlier studies, the radiocarbon age was used. In the current study the age was calibrated against the IntCal 13 curve with application of OxCal software, version 4.3 [14].

2.3. Mathematical Model

Determination of the temperature of calcium precipitation based on stable isotopes of carbon required calculation of equilibrium values of δ13C in relation to the main speciations of carbon dioxide in the solution and gaseous carbon dioxide. Therefore, calculations were based on chemical and isotopic equilibrium, taking into account solution temperature and mass balance. The first stage of the calculations included determination of molar concentrations of CO32−, HCO3 and CO2(aq) in the solution based on the chemical equilibrium model. PHREEQC-2 v. 2.0.53 software [15] was used in the calculations. The input data were results of physical-chemical parameters of the waters. Molar concentrations mCO32−, mHCO3− and mCO2(aq) calculated by the chemical mass balance model were next applied for calculating δ13C in relation to these speciations. It was assumed that DIC concentration in the solution, expressed in moles per litre, includes the sum of moles of stable isotopes of inorganic carbon per litre of solution:
D I C = m C DIC 13 + m C DIC 12
This assumption allows to express the [13C/12C]DIC ratio as:
m C DIC 13 D I C m C DIC 13 = [ C 13 C 12 ] DIC
and therefore:
m C DIC 13 = [ C 13 C 12 ] DIC D I C 1 + [ C 13 C 12 ] DIC
Analogous equations may be expressed in relation to the main speciation forms of carbon dioxide in the solution, e.g. for HCO3:
m C 13 HCO 3 = [ C 13 C 12 ] HCO 3 m HCO 3 1 + [ C 13 C 12 ] HCO 3
Because in Equation (4) the [13C/12C]HCO3 ratio is the unknown, it is easier to express it by the isotope fractionation ratio in relation to gaseous carbon dioxide α13CHCO3-CO2(g):
m C 13 HCO 3 = α 2 [ C 13 C 12 ] CO 2 ( g ) m HCO 3 1 + α 2 [ C 13 C 12 ] CO 2 ( g )
In Equation (5), the α13CHCO3-CO2(g) ratio was expressed as α2 for simplification. Similarly, the fractionation coefficients of carbon stable isotopes among other speciations and the gaseous carbon dioxide are expressed below as α1 and α3:
α 1 = α C 13 CO 3 CO 2 ( g ) α 2 = α C 13 HCO 3 CO 2 ( g ) α 3 = α C 13 CO 2 ( aq ) CO 2 ( g )
Coefficients α1, α2 and α3 were determined based on the solution temperature T [K] and the following equations [16,17,18]:
10 3 ln α 1 = 0.87 · 10 6 T 2 3.4 10 3 ln α 2 = 9.552 · 10 3 T 1 24.10 10 3 ln α 3 = 0.373 · 10 3 T 1 + 0.19
Such approach allows on one hand taking into account the thermal conditions of the solution in the calculations, and on the other minimises the number of unknowns in the system of equations, which, in order to calculate the isotope ratio [13C/12C]CO2(g) in gaseous carbon dioxide, should be expressed as:
{ m C 13 CO 3 = α 1 [ C 13 C 12 ] CO 2 ( g ) · m CO 3 1 + α 1 [ C 13 C 12 ] CO 2 ( g ) m C 13 HCO 3 = α 2 [ C 13 C 12 ] CO 2 ( g ) · m HCO 3 1 + α 2 [ C 13 C 12 ] CO 2 ( g ) m C 13 CO 2 ( aq ) = α 3 [ C 13 C 12 ] CO 2 ( g ) · m CO 2 ( aq ) 1 + α 3 [ C 13 C 12 ] CO 2 ( g ) m C 13 DIC = m C 13 CO 3 + m C 13 HCO 3 + m C 13 CO 2 ( aq )
After substituting the right sides of the first three equations to Equation (4), and transforming it, Equation (9) was obtained, in which, for simplification, unknown [13C/12C]CO2(g) was expressed as x:
α 1 α 2 α 3 ( C 13 DIC ( m CO 3 + m HCO 3 + m CO 2 ( aq ) ) ) x 3 + ( α 1 α 2 ( C 13 DIC ( m CO 3 + m HCO 3 ) ) + α 2 α 3 ( C 13 DIC ( m HCO 3 + m CO 2 ( aq ) ) ) + α 1 α 3 ( C 13 DIC ( m CO 3 + m CO 2 ( aq ) ) ) ) x 2 + ( α 1 ( C 13 DIC m CO 3 ) + α 2 ( C 13 DIC m HCO 3 ) + α 3 ( C 13 DIC m CO 2 ( aq ) ) ) x + C 13 DIC = 0
It is a cubic equation in the form of:
A x 3 + B x 2 + C x + D = 0
The coefficients of this equation were determined based on earlier calculated values of fractionation coefficients and molar concentrations for particular speciations. The concentration of isotope 13CDIC was assumed based on laboratory determinations. After calculating coefficients A, B, C and D, the real elements of this equation were determined and thus x was obtained, i.e. [13C/12C]CO2(g). Next, based on the following relationships it was possible to calculate δ13C for particular speciations:
δ 13 C CO 2 ( g ) = ( [ C 13 C 12 ] CO 2 ( g ) [ C 13 C 12 ] VPDB 1 ) · 1000
And
δ 13 C CO 3 = α 1 δ 13 C CO 2 ( g ) + α 1 · 1000 1000 δ 13 C HCO 3 = α 2 δ 13 C CO 2 ( g ) + α 2 · 1000 1000 δ 13 C CO 2 ( aq ) = α 3 δ 13 C CO 2 ( g ) + α 3 · 1000 1000
The calculated values of δ 13 C CO 2 ( g ) and δ 13 C CaCO 3 determinations in calcite obtained from samples of calcareous tufa were used to interpret the thermal conditions of tufa formation. This was achieved by using equation [19]:
10 3 ln α CO 2 ( g ) CaCO 3 = 2.988 · 10 6 T 2 + 7.6663 · 10 3 T 1 2.4612
At the assumption that calcite precipitation during tufa formation took place in isotope equilibrium conditions, Equation (13) may be used to reconstruct the water temperature. This requires determining the fractionation coefficient from equation:
α CO 2 ( g ) CaCO 3 = δ 13 C CO 2 ( g ) + 1000 δ 13 C CaCO 3 + 1000
and transformation of Equation (13) to:
( 2.4612 + 10 3 ln α CO 2 ( g ) CaCO 3 ) T 2 7.6663 · 10 3 T + 2.988 · 10 6 = 0
In the value range of α CO 2 ( g ) CaCO 3 typical for the discussed environment only one element of the equation may be used to reconstruct water temperature; it may be determined from equation:
T [ ° C ] = 7.6663 · 10 3 5.8772 · 10 7 2.9416 · 10 7 1.1952 · 10 10 ln α CO 2 ( g ) CaCO 3 4.9224 + 2 · 10 3 ln α CO 2 ( g ) CaCO 3 273.15
Additionally, a thermometer based on stable oxygen isotopes was used for palaeotemperature reconstruction. In this case the isotope equilibrium may be expressed as [20]:
10 3 ln α CaCO 3 H 2 O = 2.78 · 10 6 T 2 2.89
Similarly as in the case of stable isotopes of carbon, the value of the fractionation coefficient was calculated from equation:
α CaCO 3 H 2 O = δ 18 O CaCO 3 + 1000 δ 18 O H 2 O + 1000
As δ18OCaCO3 were accepted determinations for calcareous tufa, whereas the δ18OH2O value was accepted as a fixed value corresponding to the average value of the ratio of stable isotopes of oxygen in groundwater supplying the springs in the vicinity of the cupola. In Equation (18), δ18O values for water and calcite should be expressed in relation to the same standard. The temperature of calcite precipitation was calculated according to equation:
T [ ° C ] = 2.78 · 10 6 10 3 ln α + 2.89 273.15

3. Results

The achieved results indicate that waters of the Cretaceous and Quaternary formations have a composition typical of hypergenic zones They represent typical freshwater with mineralisation not exceeding 550 mg/L (averagely 418.9 mg/L), representing the HCO3-Ca and subordinately HCO3-SO4-Ca types (Table 1). Differences in the ionic composition can be observed between the Upper Cretaceous and the Quaternary. Quaternary waters contain more sulphates and (to a lesser degree) magnesium.
δ18O and δ2H determinations in the analysed groundwater (Table 2) are within the range of −11.43–−9.35‰ (averagely −10.29‰) and −78.89–−72.27‰ (averagely −76.56‰). The obtained range of δ18O‰ values corresponds to the range typical of precipitation water registered in stations located on continental areas. At average annual atmospheric temperature in the range of 8–10 °C, δ18O values in precipitation water are typically within −9‰–−11‰, as confirmed also by data from Poland [21,22]. The obtained results indicate that the groundwaters studied are infiltration waters, and the main factor influencing the fractionation of stable isotopes of oxygen is atmospheric temperature. The infiltration origin of the water is also indicated by the fact that the results of δ18O coupled with values of δ2H are close to the global meteoric water line (GMWL) (Figure 3). Spatial diversity of the isotopic composition results from the infiltration of meltwater and variable in time intensity of evapotranspiration [23].
δ13CDIC values (Table 2) are within the range from −12.46 to −15.35‰, with an average value of −13.51‰. Calculation results (Table 2) indicate that the dominating speciation in the groundwaters studied is HCO3. At the same time, this component is characterised by a higher contribution of 13C compared to DIC. The remaining speciations, particularly CO2(aq), are according to model calculations significantly impoverished in the heavier isotope. With regard to this, particularly interesting are calculation results for gaseous carbon dioxide. δ13CCO2(g) attains values within −20.89–−22.82‰ and thus very close to the average value in the soil carbon dioxide (−23‰) of temperate latitudes, i.e. in areas with a dominating contribution of plants, whose photosynthesis follows the Calvin cycle [11].
δ18O for calcite in samples of calcareous tufa attained values between −7.18–−7.69‰, and δ13C—between −8.38–−9.65‰ (Figure 4). These data were used in an attempt to reconstruct the changes of water temperature in the outflow zone. The thermometer based on oxygen indicated that the oldest tufa series (A) was formed at conditions much cooler than presently. According to the indication of this thermometer, water temperature in the outflow zone was about 6 °C. At present, groundwater in the outflow zone is characterised by low temperature oscillations. Between 2011–2016, regardless the time of the year at which the measurements were made, the registered temperatures were between 9.4–9.9 °C (averagely 9.6 °C), which is a value close to the average annual atmospheric temperature in this area (9.0 °C between 2007–2016, synoptic meteorological station in Łódź). The carbon thermometer points to an even lower temperature of tufa formation from series A (4 °C) and, despite the fact that this result should be treated with caution, it evidently confirms distinct cooling of atmospheric temperature, and thus also of groundwater, at the end of the Atlantic period. Younger series of calcareous tufa represent deposits of much warmer intervals, although according to the indication of the oxygen thermometer, still lower by 1–2 °C from the present day values. At the same time, according to the indication of the oxygen thermometer, the amplitude of temperature oscillations in the Subboreal (i.e. during the formation of calcareous tufa from series B, C and D) did not exceed 2 °C. A completely different scenario was obtained from the carbon thermometer. Despite the fact that this thermometer also indicates that distinct climate warming took place after deposition of series A tufa, the reconstructed temperature range is much wider, and the temperature oscillations have a different course in time.
Analysis of variability in time required dating the calcareous tufa samples, in which determinations of δ18O and δ13C were made in calcite. Radiocarbon age determinations in peat (Table 3) were used in this case. After calibration of the sample ages using the IntCal 13 curve, a sedimentary rate model was constructed (Figure 5). This allowed for recalculating the depth of tufa samples on the calibrated radiocarbon age.

4. Discussion

Previous studies of deposits within the spring mire cupola at Wardzyń [9] have allowed for correlating the sedimentary succession with variably humidity conditions in the Subboreal. Analysis of the sedimentary succession and mineral composition of the tufas have allowed to relate them with periods of relatively low precipitation, whereas the peats seem to represent more humid intervals. It was assumed that calcite dominating in the tufa composition was formed by balancing of groundwater flowing to the surface with atmospheric carbon dioxide, and the tufas represent a chemical deposit that was formed within the spring basin. The presented results of stable isotope analysis of carbon confirm this hypothesis. The chemical and isotopic composition, and the temperature of contemporary groundwater have enabled calculating the theoretical values of δ13C in calcite that could be precipitated by balancing these waters with atmospheric carbon dioxide. The obtained theoretical values of δ13CCaCO3 in relation to the analysed 12 groundwater samples are between −8.39 and −10.30‰. This range is very similar to that obtained for empirical measurements reformed in relation to calcite in calcareous tufas of the cupola (−8.38–−9.65‰). Similarity of these results evidences that calcite was formed as a chemical deposit, and isotope fractionation mostly depended on the temperature. Results of tufa analysis were used to interpret the temperature of the environment in the Subboreal, with application of the chemical and isotope equilibrium model. This model, however, has certain restrictions. The basic limit is the lack of knowledge on the chemical composition and stable isotope ratios in groundwater during the Subboreal. In some cases it is possible to reconstruct isotope ratios in water based on the isotope composition of organic matter in sediments. A good example is the isotope determination of oxygen in cellulose. The results of experimental work indicate that the fractionation coefficient of oxygen isotopes between cellulose formed as a result of photosynthesis of aquatic plants and water is constant and does not depend on temperature, plant species or type of photosynthesis [24]. This allows reconstruction of the isotope composition of water, based on δ18O determinations in cellulose present in lake sediments. However, in the case of the studied cuppola, organic matter present in the calcareous tufa contains an admixture of allochtonic material and this technique could not be used. For this reason it was necessarily assumed that this composition did not differ significantly from present day waters, and average values calculated on the basis of δ18OH2O and δ13CDIC determinations in present day waters were used in the calculations. This means that the temperatures obtained both from the carbon and oxygen thermometers should not be treated as absolute values. However, due to the low variability of ratios of stable isotopes in DIC in present day groundwater, they allow for indicative tracing of changes in the environmental temperatures in the groundwater outflow zone during the Subboreal.
Temperatures reconstructed based on carbon and oxygen isotope thermometers indicate that the oldest calcareous tufa series was formed in cold climate conditions, whereas the younger series developed when the atmospheric temperature did not differ much from present day values. In this part of the succession, however, temperatures reconstructed with the use of both thermometers show a variable course of changes in time. This most probably results from the sensitivity of the fractionation process of stable isotopes of carbon to biotic factors. If both temperature and humidity changed during the formation of the cupola deposits, it may be assumed that these changes were accompanied by transformation of the vegetation cover. Potential transformation of the soil layer could thus also have influence on the fractionation of carbon isotopes both in soil carbon dioxide and in groundwater. Variability of isotope ratios could also be influenced by variable biotic conditions in the spring basin. Here, both macrophytes and phytoplankton developing within the spring basin could have had influence on isotope fractionation [25]. In such conditions, isotope fractionation in calcite was not only the function of temperature, and thus the results of the applied model remain dubious and may be treated only as indicative. Therefore, results from the oxygen thermometer were considered as more reliable and the attempt to evaluate the thermal changes of the environment in time was based on them.
By calibrating the age of the tufa samples, it was possible to compare temperatures reconstructed from the oxygen thermometer with temperature curves obtained in other sites. δ18O values for the tufa from the spring mire cupola in Wardzyń obtained by Dobrowolski et al. [5] were also used in the comparison; temperature values were calculated with the same methodology as for the determinations presented in this report (Figure 6). The obtained curves were compared with the GISP2 Ice Core Temperature and Accumulation Data curve [2,3]. The results indicate a similar trend of thermal changes in the environment with time. The cooling and warming episodes recorded in Greenland (op.cit.) are clearly reflected in the temperature curves obtained for the spring mire cupola in Wardzyń, with two clearly marked cool episodes at about 3000 and about 3900 cal years BP, separated by a slightly warmer period. A similar trend was also observed in Northern Europe. Based on thermometers using pollen analysis, chironomids and stable isotopes in lacustrine deposits, Seppa et al. [1] (2009) distinguished cool episodes between 3000–3900 cal years BP. Similar results were obtained by Stansell et. al. [7] based on δ18O determinations; they also evidenced that the largest temperature oscillations in the Subboreal took place in the summer months. A very similar course of temperature curves for Central and Eastern Europe was presented by Davies et al. [4]. They elaborated averaged temperature curves based on a qualitative reconstruction of pollen analysis from over 500 sites in Europe. The curve of average annual temperature for Central and Eastern Europe has a very similar trend to that obtained from the Wardzyń deposits, with two cool episodes at about 3900 and about 2900 cal years BP, with the cooling particularly occurring in the summer months. A similar course have also temperature curves obtained on the basis of stable oxygen isotopes in calcareous tufa from the Czech Republic and southern Poland [6,8]. Temperature changes reconstructed based on the isotope analysis of deposits from Wardzyń correspond well also with other sites from Poland. The analysis of probability density for various types of Holocene deposits [26] indicates that between 4850–3700 cal years BP, gradual climate cooling and aridification took place. Next, a warm phase began, which ended about 2850 cal years BP. A very similar trend is visible on the temperature curve from Wardzyń (Figure 6), on which the cooler episodes at about 3900 and about 2900 cal years BP mark the limits of an interval with distinctly higher temperatures during the Late Subboreal (SB2). Transition from cool climate to much warmer conditions at about 3900 cal years BP is also confirmed by isotope analysis of calcite from lacustrine deposits. Temperature curves obtained for the Gościąż Lake [27] indicate that the change of thermal conditions was rapid and took place after 3900 cal years BP.

5. Conclusions

The model presented in this paper, based on chemical and isotope equilibrium in relation to stable isotopes of carbon has allowed for calculating δ13C with regard to the main speciation forms of carbon dioxide in the solution and to gaseous carbon dioxide remaining in equilibrium with the solution. The obtained results are very similar to δ13C values in soil carbon dioxide in areas covered by vegetation, in which photosynthesis follows the Calvin cycle. Because the plant communities in the study area represent a similar type of vegetation, this lends credence to the obtained results. Moreover, these results indicate that dissolution of carbonates in the system studied is similar to that from open carbonate systems. This means that carbon isotope fractionation is largely controlled by temperature. Based on this assumption, the theoretical values of δ13C in calcite were calculated; their levels were similar to those obtained for calcite from samples of calcareous tufa. This confirms the hypothesis of the chemical nature of the calcite (formation without the contribution of living organisms) in effect of balancing the groundwater solution with atmospheric carbon dioxide in the outflow zone (spring basin). This justifies the application of stable isotopes of oxygen and carbon in palaeotemperature reconstructions during the Subboreal (i.e. the interval during which the calcareous tufa was formed). Due to the complexity of the carbon fractionation process and the related interpretative issues, results based on δ13C were treated as indicative, the main importance being attributed to δ18O.
The calculations indicate that the oldest calcareous tufa (series A) was formed about 5500 cal years BP in a dry and significantly cooler climate compared to present day conditions. The formation of subsequent series of calcareous tufa took place between 4400 and 2900 cal years BP, at thermal conditions only slightly lower compared to present day values, with two distinct cooling episodes at about 3900 and 3000 cal years BP. The course of the obtained temperature curves correlates well with the GISP2 curve obtained from the ice core in Greenland and other sites from Northern, and Central and Eastern Europe.
Results of stable isotopes of oxygen and carbon have enabled to supplement the results of previous studies of the spring mire cupola in Wardzyń [9]. Previous analyses of the mineral composition, especially of the sulphur minerals (gypsum and pyrite) accompanying calcite, have allowed to link the calcareous tufa layers with intervals with relatively low precipitation. Analyses based on stable isotopes of carbon and oxygen indicate that the tufas were formed both in cool and warm intervals. This means that the main factor influencing the accumulation of calcite in the spring basin was humidity. At low precipitation, springs in the vicinity of the cupola became less efficient, and the longer time of groundwater contact with atmospheric air favoured balancing of the solution with atmospheric carbon dioxide and in consequence precipitation of calcite. In humid intervals the springs became more effective and conditions favouring peat accumulation developed in the spring basin. In effect, a cupola was formed, which due to the succession of organic and chemical layers is at present a valuable object for palaeoenvironmental reconstructions in the Polish Lowlands for the later parts of the Holocene (particularly the Subboreal).

Author Contributions

Conceptualization: T.G. and J.J.M.; Methodology: T.G. and M.Z.; Investigation: T.G. and M.Z.; Resources: A.R.; Data curation: M.Z. and A.R.; Writing—original draft preparation: T.G.; Writing—review and editing: M.Z.; Visualization: T.G. and A.R.; Supervision: J.J.M.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Seppa, H.; Bjune, A.E.; Telford, R.J.; Birks, H.J.B.; Veski, S. Last Nine-Thousand Years of Temperature Variability in Northern Europe. Clim. Past 2009, 5, 523–535. [Google Scholar] [CrossRef]
  2. Cuffey, K.M.; Clow, G.D. Temperature, Accumulation, and Ice Sheet Elevation in Central Greenland through the Last Deglacial Transition. J. Geophys. Res. 1997, 102, 26383–26396. [Google Scholar] [CrossRef]
  3. Alley, R.B. The Younger Dryas Cold Interval as Viewed from Central Greenland. Quat. Sci. Rev. 2000, 19, 213–226. [Google Scholar] [CrossRef]
  4. Davies, B.A.S.; Brewerb, S.; Stevensona, A.C.; Guiotc, J. The Temperature of Europe during the Holocene Reconstructed from Pollen Data. Quat. Sci. Rev. 2003, 22, 1701–1716. [Google Scholar] [CrossRef]
  5. Dobrowolski, R.; Ziułkiewicz, M.; Okupny, D.; Forysiak, J.; Bałaga, K.; Alexandrowicz, W.P.; Buczek, A.; Hałas, S. Origin and Neoholocene Evolution of Spring-Fed Fens in Wardzyń, Łódź Upland, Central Poland. Golog. Q. 2017, 61, 413–434. [Google Scholar]
  6. Pazdur, A.; Pazdur, M.F.; Starkel, L.; Szulc, J. Stable Isotopes of Holocene Calcareous Tufa in Southern Poland as Paleoclimatic Indicators. Quat. Res. 1988, 30, 177–189. [Google Scholar] [CrossRef]
  7. Stansell, N.D.; Klein, E.S.; Finkenbinder, M.S.; Fortney, C.S.; Dodd, J.P.; Terasmaa, J.; Nelson, D.B. A Stable Isotope Record of Holocene Precipitation Dynamics in the Baltic Region from Lake Nuudsaku, Estonia. Quat. Sci. Rev. 2017, 175, 73–84. [Google Scholar] [CrossRef]
  8. Zak, K.; Lozek, V.; Kadlec, J.; Hladikova, J.; Cilek, V. Climate-Induced Changes in Holocene Calcareous Tufa Formations, Bohemian Karst. Quat. Int. 2002, 91, 137–152. [Google Scholar] [CrossRef]
  9. Gruszczyński, T.; Małecki, J.J.; Ziułkiewicz, M. Reconstruction of Paleohydrological Conditions in the Late Holocene Based on the Study of Calcareous Tufa in the Spring Mire of the Wolbórka River Drainage Basin (Central Poland). J. Elem. 2017, 22, 1049–1065. [Google Scholar]
  10. Epstein, S.E.; Buchsbaum, R.; Lowenstam, H.A.; Urey, H.C. Revised Carbonate-Water Isotopic Temperature Scale. Geol. Soc. Am. Bull. 1953, 64, 1315–1326. [Google Scholar] [CrossRef]
  11. Clark, I.; Fritz, P. Environmental Isotopes in Hydrogeology; CRC Press: Boca Raton, FL, USA, 1997. [Google Scholar]
  12. Marks, L.; Dzierżek, J.; Janiszewski, R.; Kaczorowski, J.; Lindner, L.; Majecka, A.; Makos, M.; Szymanek, M.; Tołoczko-Pasek, A.; Woronko, B. Quaternary Stratigraphy and Palaeogeography of Poland. Acta Geol. Pol. 2016, 66, 403–427. [Google Scholar] [CrossRef]
  13. Ziułkiewicz, M. Pionowa Strefowość Hydrochemiczna Wód Podziemnych Na Obszarze Aglomeracji Łódzkiej (Vertical Hydrochemical Zonation of Groundwater on the Łódź Metropolitan Area); Łódzkie Towarzystwo Naukowe: Łódź, Poland, 2003. [Google Scholar]
  14. Bronk Ramsey, C. Bayesian Analysis of Radiocarbon Dates. Radiocarbon 2009, 51, 337–360. [Google Scholar] [CrossRef] [Green Version]
  15. Parkhurst, D.L.; Appelo, C.A.J. Description of Input and Examples for PHREEQC Version 3—A Computer Program for Speciation, Batch Reaction, One Dimensional Transport and Inverse Geochemical Calculations; U.S. Geological Survey Techniques and Methods: Reston, VA, USA, 2013; Volume 43. [Google Scholar]
  16. Deines, P.; Langmuir, D.; Harmon, R.S. Stable Carbon Isotope Ratios and the Existence of Gas Phase in the Evolution of Carbonate Groundwaters. Geochim. Cosmochim. Acta 1974, 38, 1147–1164. [Google Scholar] [CrossRef]
  17. Mook, W.G.; Bommerson, J.C.; Staverman, W.H. Carbon Isotope Fractionation between Dissolved Bicarbonate and Gaseous Carbon Dioxide. Earth Planet. Sci. Lett. 1974, 22, 169–176. [Google Scholar] [CrossRef]
  18. Vogel, J.C.; Grootes, P.M.; Mook, W.G. Isotope Fractionation between Gaseous and Dissolved Carbon Dioxide. Z. Phys. 1970, 230, 255–258. [Google Scholar]
  19. Bottinga, Y. Calculation of Fractionation Factors for Carbon and Oxygen in the System Calcite—Carbon Dioxide—Water. J. Phys. Chem. 1968, 72, 800–808. [Google Scholar] [CrossRef]
  20. Friedman, I.; O’Neil, J.R. Compilation of Stable Isotope Fractionation Factors of Geochemical Interest. US Geol. Surv. Bull. 1977, 440, 1–12. [Google Scholar]
  21. D’Obyrn, K.; Grabczak, J.; Zuber, A. Maps of Isotopic Composition of the Holocene Meteoric Waters in Poland (in Polish with English Summary). In Współczesne Problemy Hydrogeologii; Adam Mickiewicz University: Poznań, Poland, 1997; Volume 8, pp. 331–333. [Google Scholar]
  22. Różański, K.; Araguas-Araguas, L.; Gonfiantini, R. Isotopic Patterns in Modern Global Precipitation. In Continental Isotope Indicators of Climate; American Geophysical Union Monograph: Warsaw, Poland, 1993. [Google Scholar]
  23. Krogulec, E.; Zabłocki, S.; Sawicka, K. Changes in Groundwater Regime during Vegetation Period in Groundwater Dependent Ecosystems. Acta Geol. Pol. 2016, 66, 525–540. [Google Scholar] [CrossRef]
  24. Sternberg, L.S.L.; Anderson, W.T.; Morrison, K. Separating Soil and Leaf 180 Isotopic Signals in Plant Stem Cellulose. Geochem. Cosmochim. Acta 2003, 67, 2561–2566. [Google Scholar] [CrossRef]
  25. Leng, M.J.; Marshall, J.D. Palaeoclimate Interpretation of Stable Isotope Data from Lake Sediment Archives. Quat. Sci. Rev. 2004, 23, 811–831. [Google Scholar] [CrossRef]
  26. Starkel, L.; Michczyńska, D.J.; Krąpiec, M.; Margielewski, W.; Nalepka, D.; Pazdur, A. Progress in the Holocene Chrono-Climatostratigraphy of Polish Territory. Geochronometria 2013, 40, 1–21. [Google Scholar] [CrossRef]
  27. Różański, K.; Klisch, M.; Wachniew, P.; Gorczyca, Z.; Goslar, T.; Edwards, T.W.D.; Shemesh, A. Oxygen-Isotope Geothermometers in Lacustrine Sediments: New Insights through Combined δ18O Analyses of Aquatic Cellulose, Authigenic Calcite and Biogenic Silica in Lake Gościąż, Central Poland. Geochim. Cosmochim. Acta 2010, 74, 2957–2969. [Google Scholar] [CrossRef]
Figure 1. Location sketch-map.
Figure 1. Location sketch-map.
Water 11 01945 g001
Figure 2. Schematic cross-section through the deposits of the spring mire cupola in Wardzyń.
Figure 2. Schematic cross-section through the deposits of the spring mire cupola in Wardzyń.
Water 11 01945 g002
Figure 3. δ18O and δ2H determinations of the studied groundwaters in relation to the global meteoric water line (GMWL).
Figure 3. δ18O and δ2H determinations of the studied groundwaters in relation to the global meteoric water line (GMWL).
Water 11 01945 g003
Figure 4. Determinations of stable isotopes of oxygen and carbon in calcite and temperatures obtained on their basis.
Figure 4. Determinations of stable isotopes of oxygen and carbon in calcite and temperatures obtained on their basis.
Water 11 01945 g004
Figure 5. Sedimentary rate model constructed with the use of OxCal software [14] based on the IntCal 13 curve.
Figure 5. Sedimentary rate model constructed with the use of OxCal software [14] based on the IntCal 13 curve.
Water 11 01945 g005
Figure 6. Comparison of the GISP2 curve with temperatures based on stable oxygen isotopes from calcareous tufa samples from Wardzyń. The figure presents also temperatures obtained by Dobrowolski et al. [5] based on δ18O from the neighbouring core in the same cupola.
Figure 6. Comparison of the GISP2 curve with temperatures based on stable oxygen isotopes from calcareous tufa samples from Wardzyń. The figure presents also temperatures obtained by Dobrowolski et al. [5] based on δ18O from the neighbouring core in the same cupola.
Water 11 01945 g006
Table 1. Chemical composition of the analysed groundwater.
Table 1. Chemical composition of the analysed groundwater.
Point NumberStratygraphyTemperaturepHEhSECCa2+Mg2+Na+K+NH4+ClHCO3SO42−NO3NO2TDS
(°C)-(mV)(μS/cm)(mg/L)
1K29.97.78−4740272.35.15.61.50.513.2311.10.6<0.010.08400.0
2K29.37.38−8838264.53.63.30.90.199.1231.813.6<0.010.065327.1
3Q9.97.36−4659491.87.59.83.30.30223298.955.7<0.010.08490.4
4Q9.97.478940769.533.40.7<0.00110.118333.521.790.073325.1
5Q9.87.3−96608117.26.110.31.6<0.00138.9244112.7<0.010.026530.8
6Q9.77.24−97496100.85.58.11.90.23438.924441.3<0.010.094440.8
7Q9.97.2−3663111.35.711.25.10.02527256.2111.30.60.134528.6
8K210.97.2−7954595.7851.62.145.7390.46.8<0.010.08515.4
9Q9.97.49−2835763.33.25.42.10.04410.6170.839.51.620.092296.7
10K29.67.4−6948788.374.81.40.2935.2359.96.8<0.010.027473.7
11Q10.41.6621747271.12.46.62.2<0.00115.697.647.957.180.08300.7
12Q9.57.0512845085.54.140.4<0.0015.6280.617.1<0.01<0.001397.3
Mean9.97.38−9.9488.685.95.16.51.90.3116.1255.740.620.300.08418.9
Standard deviation0.40.2096.992.917.41.82.61.20.612.377.936.222.90.086.3
Table 2. Determinations of stable isotopes of oxygen and carbon in groundwater with calculations of chemical and isotope equilibrium with regard to the main speciations of carbon dioxide.
Table 2. Determinations of stable isotopes of oxygen and carbon in groundwater with calculations of chemical and isotope equilibrium with regard to the main speciations of carbon dioxide.
Sampling PointsMeasured ValuesCalculated Values
Temperature (K)δ18O (‰)δ18O (‰)δ13CDIC (‰)DIC (mol/L)mCO3 (mol/L)mHCO3 (mol/L)mCO2(aq) (mol/L) δ 13 C C O 3 ( ) δ 13 C H C O 3 ( ) δ 13 C C O 2 ( aq ) ( ) δ 13 C C O 2 ( g ) ( )
1283.05−11.43−75.71−13.465.10 × 10−32.77 × 10−54.86 × 10−32.13 × 10−4−15.17−13.02−23.59−22.49
2282.45−10.40−76.20−13.533.80 × 10−37.15 × 10−63.41 × 10−33.84 × 10−4−14.62−12.44−23.09−21.99
3283.05−10.48−74.57−12.924.90 × 10−31.06 × 10−54.39 × 10−35.00 × 10−4−13.98−11.82−22.41−21.31
4283.05−10.53−74.34−13.063.00 × 10−37.41 × 10−62.75 × 10−32.48 × 10−4−14.34−12.19−22.77−21.67
5282.95−10.48−72.27−12.624.00 × 10−38.21 × 10−63.53 × 10−34.59 × 10−4−13.55−11.39−22.00−20.89
6282.85−9.96−78.57−13.194.00 × 10−36.66 × 10−63.47 × 10−35.24 × 10−4−13.96−11.79−22.41−21.30
7283.05−10.22−78.48−13.494.20 × 10−36.48 × 10−63.60 × 10−35.90 × 10−4−14.15−11.99−22.58−21.48
8284.05−9.35−77.90−14.096.40 × 10−39.92 × 10−65.50 × 10−38.87 × 10−4−14.76−12.64−23.10−22.00
9283.05−10.36−78.65−12.462.80 × 10−36.94 × 10−62.57 × 10−32.22 × 10−4−13.78−11.62−22.21−21.11
10282.75−9.67−78.89−14.175.90 × 10−31.39 × 10−55.33 × 10−35.60 × 10−4−15.33−13.16−23.77−22.67
11283.55--−13.821.60 × 10−36.46 × 10−61.51 × 10−38.68 × 10−4−15.36−13.22−23.74−22.64
12282.65--−15.354.60 × 10−34.23 × 10−63.71 × 10−38.82 × 10−4−15.48−13.30−23.93−22.82
Table 3. Radiocarbon determinations and dates calibrated with IntCal13 curve.
Table 3. Radiocarbon determinations and dates calibrated with IntCal13 curve.
Profile NumberDepth [m]Radiocarbon Determination (Years BP)Calibrated Date [BP]
68.2% Probability95.4% Probability
230.632810 ± 502975–2850 (68.2%)3059–2789 (95.4%)
231.834690 ± 905577–5535 (13.2%)
5479–5319 (55.0%)
5605–5272 (90%)
5183–5120 (3.1%)
5112–5066 (2.4%)
221.182780 ± 502946–2843 (59%)
2820–2800 (9.2%)
2998–2769 (95.4%)
221.483130 ± 503440–3433 (2.7%)
3401–3325 (44.3%)
3295–3253 (21.1%)
3450–3220 (95.4%)
221.763340 ± 703679–3670 (2.8%)
3641–3479 (65.4%)
3820–3794 (1.9%)
3761–3752 (0.5%)
3725–3440 (90.5%)
3433–3400 (2.5%)
222.334600 ± 605459–5375 (31.4%)
5331–5280 (21.8%)
5164–5135 (7.8%)
5105–5076 (7.3%)
5569–5559 (0.8%)
5471–5212 (68.4%)
5193–5050 (26.2%)

Share and Cite

MDPI and ACS Style

Gruszczyński, T.; Małecki, J.J.; Romanova, A.; Ziułkiewicz, M. Reconstruction of Thermal Conditions in the Subboreal Inferred from Isotopic Studies of Groundwater and Calcareous Tufa from the Spring Mire Cupola in Wardzyń (Central Poland). Water 2019, 11, 1945. https://doi.org/10.3390/w11091945

AMA Style

Gruszczyński T, Małecki JJ, Romanova A, Ziułkiewicz M. Reconstruction of Thermal Conditions in the Subboreal Inferred from Isotopic Studies of Groundwater and Calcareous Tufa from the Spring Mire Cupola in Wardzyń (Central Poland). Water. 2019; 11(9):1945. https://doi.org/10.3390/w11091945

Chicago/Turabian Style

Gruszczyński, Tomasz, Jerzy J. Małecki, Anastasiia Romanova, and Maciej Ziułkiewicz. 2019. "Reconstruction of Thermal Conditions in the Subboreal Inferred from Isotopic Studies of Groundwater and Calcareous Tufa from the Spring Mire Cupola in Wardzyń (Central Poland)" Water 11, no. 9: 1945. https://doi.org/10.3390/w11091945

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop