Next Article in Journal
Reconstructing Paleoflood Occurrence and Magnitude from Lake Sediments
Previous Article in Journal
Palaeogeographical Reconstruction of Ancient Diolkos Slipway by Using Beachrocks as Proxies, West Corinth Isthmus, Greece
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Late Holocene Stable Isotope and Carbon Accumulation Record from Teringi Bog in Southern Estonia

1
Department of Geology and Environmental Geosciences, Northern Illinois University, 312 Davis Hall, Normal Road, DeKalb, IL 60115, USA
2
Department of Geological Sciences, University of Alaska Anchorage, 3211 Providence Drive, Anchorage, AK 99508, USA
3
Institute of Ecology, Tallinn University, Uus-Sadama 5, 10120 Tallinn, Estonia
4
Department of Earth and Environmental Sciences, University of Michigan, Earth and Environmental Sciences, Ann Arbor, MI 48109, USA
5
Triton College, Chicago, IL 60171, USA
6
Waubonsee Community College, Sugar Grove, IL 60554, USA
7
Huntley High School, Huntley, IL 60142, USA
*
Author to whom correspondence should be addressed.
Quaternary 2022, 5(1), 8; https://doi.org/10.3390/quat5010008
Submission received: 17 August 2021 / Revised: 20 December 2021 / Accepted: 4 January 2022 / Published: 1 February 2022

Abstract

:
Radiocarbon-dated peat cores collected from an ombrotrophic bog in southern Estonia record shifting environmental conditions and carbon accumulation rates in northern Europe during the late Holocene. Modern observations indicate that the water balance of the peatland is highly influenced by changes in relative humidity, followed by temperature and precipitation. The modern δ18O and δ2H values of surface water suggest that the groundwater is an integration of several months of precipitation. There also appears to be little or no direct influence of surface evaporation on the water within the bog, suggesting that water loss is preferentially through transpiration and sub-surface flow. Bulk peat δ13C values exhibit a trend of higher values through the late Holocene, suggesting a pattern of overall increased surface wetness. The δ15N values were low from ~4130 to 3645 cal yr BP, suggesting drier conditions, followed by intermediate values until ~2995 cal yr BP. The δ15N values decrease again from ~2995 to 2470 cal yr BP, suggesting a return to drier conditions, followed by intermediate values until ~955 cal yr BP. The δ15N values were high, suggesting wetter conditions from ~955 to 250 cal yr BP, followed by intermediate values through the modern. Carbon accumulation rates were low to intermediate from ~4200 to 2470 cal yr BP, followed by intermediate-to-high values until ~1645 cal yr BP. Carbon accumulation rates were then low until ~585 cal yr BP, followed by intermediate values through the modern. The geochemical data, combined with observed changes in peat composition and regional proxies of temperature and water table fluctuations through the late Holocene, suggest that carbon accumulation rates were relatively low under dry and warm conditions, whereas accumulation was generally higher (up to ~80 g C m−2 yr−1) when the climate was wetter and/or colder. These findings further suggest that future environmental changes affecting the regional water balance and temperature will impact the potential for northern peatlands to capture and store carbon.

1. Introduction

Peatlands are an integral part of the global carbon cycle, and they are sensitive to changing environmental conditions. High northern latitude peatlands cover approximately 4 million km2 and store up to 550 Pg (1 Pg = 1015 g) of carbon [1,2,3], equivalent to almost one-third of the world’s soil carbon [4], and more than half of today’s atmospheric CO2 [5]. Nevertheless, it is unclear how this internationally important carbon reservoir responds to specific climate forcings and how these responses will drive ecosystem-climate feedbacks [6,7,8,9]. Previous studies highlight the impacts of climate change on peatland carbon dynamics at a large scale [10]; however, the local- and regional-scale responses of peatland hydrology and carbon accumulation to environmental changes need to be better evaluated to improve dynamic modeling of peat accumulation through time [11].
The ecological and biological characteristics of peatlands are largely governed by hydrological conditions. For example, the trophic status of a peatland is classified principally on hydrologic regime, and the height of the water table is the primary delineating factor separating physical and chemical processes in these systems [12]; the depth of the water table separates the acrotelm with oxic conditions (above) from the catotelm with anoxic conditions (below). Recent studies suggest that there is an intermediate zone, the mesotelm [13], where the fluctuating water table causes variations in redox conditions that likely increase carbon turn over [14]. As such, any process that affects water table conditions has the potential to drive environmental responses in these highly sensitive ecosystems. In ombrotrophic systems, these processes should be dominated by changes in humidity, temperature, precipitation amounts, and residence times of water flowing through peat. The oxygen isotopic composition of precipitation (δ18Op) of Estonian surface water tracks these environmental changes [15,16,17,18,19], and in turn, provides insight into the modern hydrologic processes operating in bog systems. From a paleo-perspective, carbon and nitrogen isotopes (δ13C and δ15N) serve as proxies of surface wetness [15] and connectedness to the water table [16].
Ultimately, the carbon balance of peatlands is related to productivity, respiration, and decomposition. Whether a peatland will be a source or a sink of greenhouse gases is also closely related to the hydrologic dynamics described above because the position of the water table determines the boundary between aerobic and anaerobic microbial respiration. Both biotic production and peat decomposition rates in peatlands are controlled largely by water availability and temperature [1]. Most newly formed biomass is subsequently decomposed by aerobic decay and then respired. As plant litter transitions into the anaerobic catotelm, some is decomposed to CO2 and CH4, while the rest is disconnected from atmospheric exchanges and stored as peat. Much of the CO2 and CH4 is lost at the surface through diffusion, bubbling, or stem and root aerenchyma [12]. Therefore, depending on hydrologic conditions, peatlands serve as both sinks and sources of CO2 and CH4 [3,17], and the fate of the stored fraction of this carbon has important consequences for global climate feedbacks [10,18,19,20].
Here, we present a ~4200-year-old stable isotope (δ13C, δ15N) and carbon accumulation record from Teringi Bog, an ombrotrophic system in southern Estonia, to reconstruct the timing, magnitude, and direction of centennial- to millennial-scale shifts in hydroclimate and carbon accumulation rates during the late Holocene. Located in the Baltic region, it is sensitive to North Atlantic-driven ocean-atmospheric processes, and therefore also has the potential to provide a local perspective of hydroclimate variability in response to regional circulation dynamics. We compare these new records to insolation [21] and regional proxy data of water table, temperature and precipitation changes to improve the paleo-perspective of climate and carbon dynamics from continental northern Europe.

2. Study Site

Teringi Bog (57°58′33″ N, 25°33′30″ E) is a rainfall-dominated peatland system in southern Estonia (Figure 1), which is located in the boreo-nemoral vegetation zone [22]. The region was ice covered during the latest Pleistocene as part of the Scandinavian glaciation, and it has remained ice free for the last ~14,000 years BP [23]. Teringi, like most peatlands in southern Estonia, sits atop remnant glacial till and sands, although the raised bog conditions have separated most of it from the mineral sediments beneath.
The biodiversity at Teringi is similar to other peatlands in southern Estonia [24,25]. In brief summary, the bog is surrounded by mixed forests of Picea (spruce), Betula (birch) and Pinus (pine), as well as cultivated land and meadows. Within the bog, Pinus is the most common tree species, followed by Betula pubescens and Betula nana on the edges (Figure 2). The Pinus forms vary by shape and size based on connectedness to underlying mineral soils and peat thickness. These pines include f. litwinowii as the most common throughout the bog, followed by f. uliginosa in the areas with more mineral-rich soils, f. wilkommi in areas of thicker peat, and f. pumila in areas of nutrient-poor conditions and thick moss. There are several species of Salix, Sorbus and Populus on the edges of the bog, while Picea, Frangula alnus and Juniperus are less common. The dwarf shrub layer in the bog contains common peatland plants: Caluna vulgaris, Ledum palustre, Andromeda polifolia, Chamaedaphne calyculata and Vaccinium uliginosum. The grass and herb layer consists of Eriophorum vaginatum, Trichophorum cespitosum, Rhynchospora alba, Rubus chamaemorus and Drosera spp. The moss layer contains mostly Polytrichum strictum.
The undecomposed top layer of peat at Teringi is entirely dominated by Sphagnum sect. Acutifolia (e.g., Sphagnum fuscum) and sect. Cuspidata (e.g., Sphagnum angustifolium), which are distributed in relation to the microtopography of the bog. Species comprising the Sphagnum sect. Acutifolia live in hummocks, several centimeters above the water table; species of sect. Cuspidata live much closer to the water surface in lawns or hollows.
Estonia has a temperate climate that is transitional between oceanic and continental, with changes in North Atlantic Ocean-driven circulation dynamics affecting the environmental conditions of the region [26,27]. Stronger westerly zonal winds deliver warmer ocean air to the interior, causing warmer temperatures, while weaker circulation causes colder conditions [28]. The majority of moisture is sourced from the Atlantic Ocean while the Mediterranean Sea and Black Sea provide small secondary sources [29]. The Baltic Sea plays a negligible role in the development of precipitation [30]. Average annual precipitation amounts at Teringi (measured at Valga) are ~685 mm. Average January temperature at Teringi is ~ −5.5 °C while July is ~17.3 °C (Figure 3).

3. Materials and Methods

3.1. Modern Hydrology

A water level and temperature logger (Solinst 3001 LT Levelogger Junior) was installed in summer 2014 to monitor hydrological changes, and measurements were collected every hour from 14 August 2014 to 1 July 2021 (Figure 4). The well housing for the logger was constructed of PVC pipe approximately 6 cm in diameter and 1.1 m in length. The bottom ~40 cm of the well contained mesh-covered holes, to allow for equilibration with groundwater and to prevent clogging of the PVC pipe by surrounding peat. The pressure logger was positioned at ~80 cm depth in the well and both ends of the well were capped. In order to determine changes in peatland water-table depth, the pressure data from the logger was corrected using atmospheric pressure data collected from a barometric logger placed ~2 m above ground at a location ~15 km north of the peatlands (Karksi, Estonia, see Figure 1A).
A remote, automated weather station (HOBO U30-NRC and various sensors) was established near the peatland site (also at Karksi, Estonia) in May, 2015 in order to record weather variables including precipitation (tipping bucket), relative humidity (RH) and temperature. Sensors scanned every 1 min and logged every 1 h (Figure 4).

3.2. Modern Water Sampling

A total of 148 event-based rainfall, lake, wetland, and river samples for δ18O and δ2H analysis were collected in air-tight 15 mL vials, with no headspace, and analyzed as part of a multi-year monitoring effort [31]. A total of 26 surface-water and well-water samples were collected across a west–east surface transect through Teringi during the summers of 2014, 2015 and 2016 (Figure 5). An additional 6 surface-water and well-water samples were collected during December, 2014. All Estonia water isotope samples were measured for δ18O and δ2H ratios using a Piccaro water isotope analyzer in the Stable Isotope Laboratory at the University of Alaska Anchorage. All water samples were screened to verify that organic content was low and not a source of absorption interference spectra [32]. Each sample was then analyzed six times and reanalyzed if the standard deviation of the six replicates was greater than 0.3‰ for δ18O and/or 3‰ for δ2H, or if the internal standard for the run differed from the accepted value by greater than ±0.2‰ or 2‰, for δ18O and δ2H, respectively. All isotope results are reported in delta notation (δ) relative to standards according to the equation δ [‰] = (Rsample − Rstandard /Rstandard) ∗ 1000‰ [33].

3.3. Core Sample Collection

A peat record that was 528 cm in length was collected from the center of Teringi Bog (57°58′33.56″ N, 25°33′30.90″ E) by extracting continuous 50 cm sections from the same location using a Russian peat corer. The corer was unable to penetrate through any material below 528 cm. Cores were extruded in the field into split PVC tubes and wrapped in plastic. The uppermost ~50 cm of peat were collected separately by hand as a monolith, wrapped in plastic and stored in PVC tubing. All peat samples were kept in refrigerated storage after collection.

3.4. Chronology

Age control for the composite peat record was established using accelerator mass spectrometry (AMS) 14C analysis of discrete peat (Sphagnum) stem material (Table 1). The 14C samples were pretreated at Northern Illinois University using standard acid-base-acid wash techniques and processed and measured at the W.M. Keck Carbon Cycle AMS Laboratory at the University of California, Irvine. The 14C ages were calibrated using the IntCal20 calibration curve [34]. An age-depth model was created using a cubic spline interpolator with Bayesian statistical methods [35], and the peat core time-series data presented are based on the weighted mean ages modeled using Bacon software in R (Figure 6).

3.5. Species Identification

The core was sampled every 10 cm for species identification, and bulk Sphagnum sp. remains were placed in a small Petri dish and submerged in distilled water. Approximately 20 to 50 leaves were mounted on glass microscope slides, and individual species were determined by leaf morphology and placement of pores in the leaf structure under 40×–100× magnification following established methods [36] (Figure 7).

3.6. Carbon and Nitrogen Concentrations and Stable Isotopes

A total of 53 samples were collected down-core from sections 1 cm in thickness. Samples that contained both stems and branches were collected consistently down-core following standard methods [37]. The samples were dried overnight at 70 °C in a drying oven and homogenized using a combination of mortar and pestle and razor blade. The dried and homogenized material was weighed (approximately 3 mg for N and 0.3 mg for C) and loaded into tin capsules. Isotopic analyses and total carbon and nitrogen concentrations were conducted at the Stable Isotope Laboratory at the Environmental and Natural Resources Institute at University of Alaska Anchorage. Carbon (δ13C) and nitrogen (δ15N) isotope values, as well as C and N concentrations were measured using a coupled Costech Instruments elemental analyzer and ThermoFinnigan Electron Delta Plus XP continuous flow isotope ratio mass spectrometer (IRMS). All isotope values were calibrated with internal standards. Analytical precision was 0.08‰ for carbon and 0.1‰ for nitrogen.
Bulk density and loss on ignition (LOI) was also measured on 1 cm3 samples collected every 4 cm down-core at 550 °C to estimate the percent organic matter [38,39,40]. Organic carbon content was determined by multiplying the amount of organic matter by 0.423, the typical value used for Sphagnum-dominated systems [41]. Flux of organic carbon was determined for both the IRMS and LOI data by multiplying the percent of organic carbon of dry bulk density (g/cm3) by the accumulation rate (cm/yr) (Figure 8 and Figure 9).
Time-series analyses were performed on all environmental datasets. Significant regime shifts in all time-series data were identified using Sequential Regime Shift Detection software (v. 6.2) [42], with a 95% confidence level, a cutoff length of 5, and a Huber’s tuning constant of 1 (Figure 8 and Figure 9). A Principal Component Analysis (PCA) of the sedimentological, stable isotope (δ13C and δ15N) and bulk geochemical data was done using the R FactoMineR package [43], after scaling the data by variance (Figure 10). The PCA also includes late Holocene temperature data from Fennoscandia [44] and water table data from Männikjärve Bog in Estonia [45].

4. Results

4.1. Weather Station Data

Weather station data reveal variations consistent with seasonal changes (Figure 4). The seasonal patterns, as well as the maximum and minimum values of the air temperature data from the automated weather station at Karksi are similar to the 30 year averages of reprocessed data for the region (Figure 3). Temperatures are highest from late May through September, peaking in late July and early August. The coldest temperatures are from January through March. Precipitation events occur throughout the year, with generally wetter conditions in late summer and fall. The measured precipitation values in winter are well below the 30 year averages from Valga, but this is likely because the battery powered weather station at Karksi is incapable of recording accurate winter precipitation amounts. Relative humidity values are highest in late fall through winter, lowest in spring, and intermediate through summer.
Peatland water temperatures varied seasonally, and as expected, the magnitude of changes was lower than that of air temperatures (Figure 4). The highest water temperatures are generally from June through September, peaking in late July and early August. The coldest temperatures are generally from January to April. There is an apparent abrupt increase in temperature at the bog during the spring melt season and this lag relative to changes in air temperature might be exaggerated as the ice covered datalogger probably records the onset of thawed conditions. The water-level data revealed sharp drawdowns of water levels from May through September of most years. Water levels generally remain low during the early-to-mid summer months, followed by fluctuating but relatively high values through late summer, fall, winter and early spring.

4.2. Surface Water and Precipitation Stable Isotopes

Surface water and precipitation isotope values vary seasonally and by location (Figure 5). Cold season (November–April) precipitation and surface water samples have relatively low δ18O values, generally in the range of ~−23 to −9‰ VSMOW, while warm season (May–October) values are higher (~−10 to 0‰ VSMOW). For comparison, the Teringi Bog, river, and other open-basin lakes all plot near the global meteoric water line (GMWL) with little to no influence of evaporation. The bog samples collected between May and August from both the well and surface water generally range from ~−10.5‰ to −9.1‰ VSMOW. Surface water samples collected across the bog in December range from ~−18.9‰ to −9.0‰ VSMOW while the well water at the coring site was ~−9.5‰ VSMOW.

4.3. Chronology

The initial 14C, calibrated, and Bayesian age-depth model results are presented in Table 1 and Figure 6. The record spans the last ~4200 years with no hiatuses. The youngest sample (140 14C BP) ranged from ~60 to 120 cal yr BP with a median probability of 110 cal yr BP, and a weighted mean modeled age of 270 cal yr BP. The other 5 samples have median probability values that are within the >90% probability distributions and the weighted mean modeled ages are within the 2-sigma calibrated age ranges.

4.4. Species Identification

Acutifolia was the primary Sphagnum section throughout the core (Figure 7 and Table 2), and the dominant species were S. fuscum and S. capillifolium [46]. S. fuscum, S. capillifolium and S. angustifolium make up the bulk of the record from ~4200 to 3960 cal yr BP. S. fuscum, S. capillifolium and S. medium are dominant from ~3960 to 2920 cal yr BP, followed by S. fuscum, S. capillifolium and S. russowii from ~2920 to 1700 cal yr BP, and then S. fuscum, S. capillifolium and S. angustifolium from ~1700 to 1130 cal yr BP. S. fuscum and S. capillifolium are dominant from ~1130 to 960 cal yr BP, followed by S. fuscum and S. rubellum from ~960 to 700 cal yr BP. From 700 cal yr BP to the modern, species are much more diverse and include S. fuscum and S. capillifolium, S. rubellum, S. balticum, S. fallax, and S. pulchrum.

4.5. Peat Stable Isotope Records

The δ13C values range from ~−27.8‰ to −24.5 ‰ (VPDB), which is typical of C3 forming peat plants [47]. The δ13C values were intermediate from ~4130 to 2130 cal yr BP, followed by lower values until ~1810 cal yr BP and a return to intermediate values through the modern. The δ13C values show an overall increasing trend from the beginning of the late Holocene to the present.
The δ15N values ranged from ~−6.6‰ to 1‰ (versus air), with most values varying between ~−4‰ to 0‰ (Figure 8 and Figure 9). The δ15N values were low from ~4130 to 3645 cal yr BP, followed by intermediate values until ~2995 cal yr BP, and a shift to lower values until ~2470 cal yr BP. The δ15N values were then intermediate until ~955 cal yr BP, followed by higher values until ~250 cal yr BP, and then a return to intermediate levels through the modern. The lowest δ15N value was in the upper-most (2 cm) sample of the core.

4.6. Carbon, Nitrogen and Carbon Accumulation Rates

The bulk density values were intermediate from ~4200 to 3350 cal yr BP, followed by lower values until ~1860 cal yr BP, and then a return to high values until ~1790 cal yr BP. Values then decreased until ~1675 cal yr BP, followed by high values until ~1500 cal yr BP, and a return to lower values until ~1295 cal yr BP. Values were then intermediate until ~585 cal yr BP, followed by a shift to higher values until ~340 cal yr BP. Bulk density returned to intermediate values through the modern.
The mass ratios of C and N (C:N) were low from ~4130 to 3170 cal yr BP, followed by higher values until ~1630 cal yr BP and a return to low values through the modern. The %C values were intermediate from ~4130 to 1630 cal yr BP, followed by a slight increase until ~170 cal yr BP and then a shift to lower values through the modern. The %N values were intermediate from ~4130 to 3840 cal yr BP, followed by higher values until ~3170 cal yr BP. There was then a return to lower % N values until ~1630 cal yr BP, followed by intermediate values through the modern.
The carbon accumulation (flux) rates for the LOI and IRMS methods produced similar results. The percentage of organic matter in all the LOI samples was >97%, and we thus assume that any mineral sediments are negligible. The carbon accumulation rates were low from ~4200 to 2470 cal yr BP, followed by intermediate-to-high values until ~1645 cal yr BP. Carbon accumulation rates were then lower until ~585 cal yr BP, followed by a return to intermediate values through the modern.

4.7. Principal Component Analysis

The results of the PCA (Table 3 and Figure 10) highlight varying associations between past changes in water table depth, temperature and proxies for environmental changes at Teringi Bog. For example, the %N and C:N data correlate strongly with PC1 as shown with high positive (0.94) and negative loadings (−0.92). The %C values have only a 0.46 loading with PC1 and do not correlate strongly with either PC2 or PC3. Organic carbon flux, peat accumulation rates and bulk density correlate with PC2 with loadings of 0.88, 0.66 and 0.65. Water table depths also have a moderately strong correlation with PC2 with a loading of −0.62. July temperatures, peat accumulation rates and õ13C correlate with PC3 with loadings of 0.50, 0.45 and −0.69. The õ15N values weakly correlate with PC1, PC2 and PC3 with loadings of 0.56, 0.44 and -0.53.

5. Discussion

5.1. Surface Water Source and Seasonality

The δ18Op of mid-to-high latitude rainfall is primarily controlled by air mass distillation and atmospheric temperature [48,49]. This results in a positive correlation between winter δ18Op and temperature in Europe [50]. Additionally, Punning et al. [30] identified a temperature gradient of isotopic composition of δ18Osurfacewater of 0.32‰/°C for Estonia. Moreover, there is a strong negative correlation between winter (December–March) precipitation amount and δ18Op [50], indicating that seasonality is a dominant driver of δ18Op variability recorded in Estonian surface waters. Variations in snow pack thickness [51,52] should be negligible for precipitation-fed systems like Teringi that lack other water sources [53]. Most of the surface water δ18O values for samples collected from Teringi vary only by ~2‰ throughout the year (from −11 to −9‰), which are intermediate between warm and cold season precipitation values recorded at the nearby Tartu GNIP station (Figure 3). This suggests that the bog water is an integration of several months of precipitation and that the residence time of water in the bog is probably at least several months to years. There appears to be little or no direct influence of surface evaporation on the water within the bog, suggesting that water loss is preferentially through transpiration and sub-surface flow.
The local weather station and water-logger data suggest that the timing of water level fluctuations corresponds most closely to changes in RH, followed by temperature and rainfall changes (Figure 4). For example, the water-logger data revealed sharp drawdowns of water levels from May through September of most years when RH values were low to intermediate. Water levels generally remained low during the early-to-mid summer months when temperatures were high and there were fewer rainfall events. Water levels then fluctuated but remained relatively high through late summer, fall, winter and early spring. These modern observations of environmental conditions suggest that changes in either hydroclimate or temperatures had the potential to impact the carbon dynamics and greenhouse gas exchanges at Teringi during the late Holocene.

5.2. Changes in Sphagnum

Species of Sphagnum shift under changing environmental conditions. For example, S. fuscum and S. capillifolium tend to be more productive under stressful conditions [54], and S. rubellum are more productive at higher temperatures [55]. Hummock species like S. fuscum and S. medium are generally more resistant to drought than lawn or hollow species like S. balticum and S. russowi that typically represent wetter periods [45,56,57]. Lab experiments indicate that S. rubellum is generally unaffected by modest water table changes [55]. S. angustifolium is better than species like S. fallax at growing under drier and warmer conditions [58].
We therefore broadly interpret wetter intervals as times when S. balticum and/or S. russowi increase, warmer and/or drier intervals when S. medium, S. angustifolium and/or S. rubellum have higher proportions, and intermediate conditions when S. fuscum and/or S. capillifolium dominate. It should be noted that the typical time needed for hollow species to replace hummocks in Estonia is ~180 years, and 320 years are required for hummock species to replace hollows [45]. Acknowledging that there is a lag in the plant response time that has not been accounted for here, the composition of the core suggests that the climate was overall warm and dry from ~4200 to 3960 cal yr BP, followed by intermediate-to-dry conditions until ~2920 cal yr BP. Conditions were generally wetter from ~2920 to 1700 cal yr BP, followed by a shift to somewhat warmer and drier conditions until ~1130 cal yr BP. The climate conditions were milder from ~1130 to 960 cal yr BP, followed by higher temperatures until ~700 cal yr BP. The remaining section of the core with abundant species variations through the modern likely represents overall wetter conditions and intermediate-to-cold temperatures.

5.3. Carbon and Nitrogen Isotopes

Carbon isotopes in peatland systems are affected by multiple environmental processes, but in certain conditions they can be an additional indicator of bog surface wetness [15]. Sphagnum mosses lack guard cells surrounding their stomata, so they cannot regulate their uptake of atmospheric CO2, and the δ13C value reflects the relative wetness or dryness of the plant’s environment. The CO2 enters the Sphagnum leaf through the stomata into the large and dead hyaline cells, and it must then pass through a membrane between the hyaline and photosynthetic cells before it can be utilized for photosynthesis [59]. When the hyaline cells are filled with air, CO2 can more easily diffuse through the membrane, allowing the plant to preferentially uptake the lighter isotope, which leads to lower δ13C values. When the hyaline cells are filled with water, diffusion is more difficult, and the plant cannot discriminate against the heavy isotope as easily, leading to higher δ13C values. For Sphagnum mosses, once the CO2 enters the photosynthetic cells, the process of fixing CO2 into usable compounds produces a fractionation of approximately −26.5‰ [59,60]
In ombrotrophic bogs, Sphagnum takes in nitrogen directly from the atmosphere and water [16], and measurements of nitrogen isotopes (δ15N) on peat have the potential to identify changes in past hydroclimate conditions. Hummock species are more influenced by rainwater/surface bog water than deeper bog water, and rain water has lower δ15N values than bog water. Therefore, hummock species should have lower δ15N values than hollow species [16,61,62]. During drier conditions, a lower water table will cause the plant to utilize more rain water, decreasing the δ15N values [16].
The δ13C and δ15N values from the Teringi record provide complementing perspectives of likely shifts in hydrologic conditions during the late Holocene. The δ13C values taken from bulk peat at Teringi show a slight step-wise increase through time, suggesting that conditions have become increasingly wetter. The δ15N values in the core are all generally within range (−4‰ to 0‰) of typical Sphagnum mosses [16] and are more variable than the δ13C values. The δ15N values were low from ~4130 to 3645 cal yr BP, suggesting drier conditions, followed by intermediate values until ~2995 cal yr BP. There was another shift to likely drier conditions from ~2995 to ~2470 cal yr BP. The δ15N values suggest intermediate hydroclimate conditions until ~955 cal yr BP, followed by wetter conditions until ~250 cal yr BP, and then values were intermediate again through the modern. Here, again, the actual response time for plant species to climate forcing is uncertain, but any shifts in δ15N values probably lag any hydroclimate forcing by at least 2 centuries [45]. The results of the PCA also indicate that while there is a moderately strong relationship between water table fluctuations in Estonia and δ13C and δ15N values at Teringi, multiple environmental variables likely also drove variations in these peatland systems through time (Figure 10 and Table 3). Likewise, the lack of correlation between %C and %N with the available proxies for temperature and wetness suggests that these systems within the bog respond to complex interactions of multiple environmental variables.

5.4. Carbon Accumulation Rates and Hydrologic Conditions

Generally speaking, carbon accumulation rates for peatland systems vary largely as a function of bulk density and age-model calculations [41]. For Teringi, there are multiple possible calibrated ages for the uppermost radiocarbon sample; however, recalculating flux based on the maximum and minimum possible values (Table 1) resulted in no significant differences in the timing of shifts in the carbon accumulation record compared to the weighted mean modeled ages. Humification of peat (process of decomposing litter into humic substances) can also lead to higher bulk density and carbon flux rates; however, the C:N values in the Teringi core exhibit only slight shifts through time, suggesting that there were no prolonged periods of organic matter decay during the late Holocene [63].
Bulk density is generally higher in sections of the Teringi core where S. fuscum is more dominant and lower in sections where the composition is more diverse. Compared to most of the core, the bulk density values are high from ~1860 to 1790 cal yr BP, from ~1670 to 1500 cal yr BP, and again from ~585 to 340 cal yr BP. The bulk density values are lowest from ~3350 to 1860 cal yr BP, low from ~1790 to 1675 cal yr BP, and low again from ~1500 to 1295 cal yr BP. Otherwise, the remaining sections of core had intermediate values. Hummock species become more dominant under drier conditions while hollow species are more common when it is wetter [45], and hummock species have higher growth rates than hollow species [64]. Combined with the δ15N and C:N data, these observations suggest that bulk density values in the Teringi Bog record are probably controlled more by changes in the species that are present than by decomposition. The results of the PCA suggest that bulk density and carbon accumulation values are generally low in sections of the core that represent dry periods and high during wetter periods while temperature played only a limited role (Figure 10).

5.5. Climate Dynamics during the Late Holocene (4200 Cal Yr BP to Present)

Broad-scale patterns of shifting mean-state climate conditions during the late Holocene in Estonia can be inferred from the Teringi Bog core data combined with marine and terrestrial paleoclimate records. The onset of Neoglaciation and declining Northern Hemisphere high-latitude summer insolation [65] likely led to an overall colder climate (Figure 9). Pollen records from Lakes Raigastvere in Estonia [66] and Tsoulbmajavri in Finland [44] indeed indicate a trend toward colder conditions in the northern Baltic region during the late Holocene. Lake sediment stable isotope records suggest that precipitation patterns and circulation dynamics were variable, but that there was also a trend of overall wetter conditions through the late Holocene [31,53]. Thus, the overall trend of higher δ13C values recorded at Teringi Bog suggests that a wetter climate seems to coincide with a broader pattern of colder and wetter conditions in the Baltic region.
The basal radiocarbon age from Teringi (~4200 cal yr BP) is notable when compared to the timing and pattern of Holocene Northern Hemisphere climate change. The start of the late Holocene is defined as ~4200 cal yr BP [67] when there was a major shift in circulation dynamics and hydrologic balance across northern Europe at the onset of Neoglaciation and declining Northern Hemisphere high-latitude summer insolation [68]. It is unclear what lies below ~530 cm depth at Teringi Bog, as the material below was impenetrable. This same impenetrable layer was present in a transect across the bog [46]. Interestingly, Männikjärve Bog in northeast Estonia has a basal radiocarbon age of ~4360 cal yr BP that similarly marks the initiation of ombrotrophic conditions [45]. We therefore suggest that the environmental conditions became conducive for Sphagnum growth in the center of both bogs sometime before ~4200 cal yr BP, possibly as a result of wetter and cooler conditions in the region compared to the middle Holocene, and that these conditions remained relatively cool and wet through the late Holocene.
Multi-centennial- to centennial-scale periods of hydroclimate variability in northern Europe likely affected peatland carbon dynamics. For example, changes in water table depth based on testate amoebae from Männikjärve Bog identify drier conditions in Estonia while carbon accumulation rates were low at Teringi Bog from ~4200 to 2470 cal yr BP (Figure 9). Carbon accumulation rates were then generally the highest at Teringi Bog (up to ~80 g C m−2 yr−1) from ~2470 to 1645 cal yr BP when the bulk of the proxy records suggest that it was wet and temperatures were intermediate. These carbon accumulation values are relatively high compared to much of Europe [41], and while it is presently unclear exactly why these values are higher, it might be because Sphagnum dominates throughout the core and it is remarkably well-preserved. The Teringi δ13C and δ15N data, combined with other paleoclimate records from the region [69,70,71,72,73], suggest that hydroclimate conditions were likely intermediate to dry for several centuries prior to, and including the Medieval Climate Anomaly (MCA; ~1000 to 700 cal yr BP), compared to earlier in the late Holocene (Figure 9). Similarly, carbon accumulation rates also appear to be low at Teringi from ~1645 to ~585 cal yr BP. The Little Ice Age (LIA) in northern Europe that followed is defined as a relatively cold period from ~450 to 100 cal yr BP [74]. During the LIA, there was a shift to intermediate carbon accumulation rates compared to several centuries prior in the Teringi record that continued through present day (Figure 9). This suggests that colder and wetter LIA conditions were generally favorable for peat growth and carbon accumulation.

5.6. Future Implications

Future shifts in hydroclimate, water table dynamics, and/or temperature will impact the potential for Teringi Bog and other regional peatlands to capture and store carbon. The Baltic region will likely receive more precipitation and become 4–5 °C warmer during the next century [75,76,77,78,79]. However, there is a high degree of uncertainty in future model projections: down-scaled regional climate models (RCMs) forecast either no substantial change, or up to a ~20% increase in precipitation for the Baltic States [80,81,82]. The uncertainties of these models stem largely from a lack of climate reconstructions extending beyond the instrumental period [83], and the current model ensembles simply do not accurately represent the full range of precipitation variability that exists across the Baltic States [84]. In a range of RCM simulations, the rate of evapotranspiration in regions south of the Baltic Sea will likely increase in the future [85], which should lead to overall drier conditions. Likewise, a range of long-term climate projections indicate that the Baltic region will likely experience a decrease in relative humidity as temperatures increase, even if there is an increase in precipitation amounts [86]. It is also possible that human impacts and changes in biodiversity might outpace any climatic influences, which could significantly alter the carbon dynamics of Baltic peatlands. Nevertheless, periods of the late Holocene when there were warmer and drier conditions, combined with lower carbon accumulation rates, might serve as analogues for future environmental conditions in Estonian peatlands.

6. Conclusions

The water balance and past carbon dynamics of Teringi Bog appear to be driven mostly by variations in humidity and precipitation followed by temperature changes. Analyses of stable isotopes (δ13C and δ15N) measured on peat cores taken from the bog, combined with multiple proxy records from the region, suggest that there was a trend toward wetter conditions in Estonia during the late Holocene. Carbon accumulation rates were low from ~4200 to 2470 cal yr BP when it was generally warm and dry. Carbon accumulation rates were then higher at Teringi Bog from ~2470 to 1645 cal yr BP when conditions were wetter and temperatures were intermediate. Several centuries leading up to and including the MCA had overall low carbon accumulation rates. The LIA was also a period of increased carbon accumulation when it was colder, suggesting that evaporation rates and/or decomposition was lower during that interval. At present, it is unclear if ombrotrophic peatlands will become sinks or sources of carbon emissions in the coming decades, and more detailed local-scale studies are needed in order to understand the fate of greenhouse gases in northern European wetlands. Nevertheless, the results presented here suggest that decreased precipitation amounts, lower relative humidity and/or higher temperatures will likely lead to a negative carbon balance for Teringi Bog.

Author Contributions

Conceptualization, methodology, supervision, investigation, formal analysis, project administration, writing—original draft preparation, and funding acquisition, N.D.S.; conceptualization, methodology, supervision, investigation, formal analysis, and writing—review and editing, E.S.K.; methodology, formal analysis, investigation, and writing—review and editing, K.H.; supervision, investigation, formal analysis, and writing—review and editing, J.T. and J.D.; investigation and writing—review and editing, M.E. and B.P.; investigation, M.B., M.E., C.F., A.F., R.G. and B.S.; All authors have read and agreed to the published version of the manuscript.

Funding

This project was funded by the U.S. National Science Foundation, International Research Experiences for Students (IRES) grant (OISE-1827135). Additional funding for this project was provided by the Research and Artistry grant program and Goldich funds from Northern Illinois University.

Data Availability Statement

The peat core data for this project are available at https://www.ncdc.noaa.gov/paleo/study/34995 (accessed on 1 August 2021).

Acknowledgments

We thank Nicole LaDue, Luule Lauk, Hannes Saag, Sten Sarap, Evelin Stansell, Sheldon Turner, Tiit Vaasma, David Voorhees and the Paavel family. The Institute of Ecology at Tallinn University graciously hosted NDS for a sabbatical while this paper was written. This manuscript benefited greatly from comments by Martin Küttim and 2 anonymous reviewers.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gorham, E. Northern Peatlands: Role in the Carbon Cycle and Probable Responses to Climatic Warming. Ecol. Appl. 1991, 1, 182–195. [Google Scholar] [CrossRef]
  2. Turunen, J.; Tomppo, E.; Tolonen, K.; Reinikainen, A. Estimating carbon accumulation rates of undrained mires in Finland–application to boreal and subarctic regions. Holocene 2002, 12, 69–80. [Google Scholar] [CrossRef]
  3. Yu, Z.; Loisel, J.; Brosseau, D.P.; Beilman, D.W.; Hunt, S.J. Global peatland dynamics since the Last Glacial Maximum. Geophys. Res. Lett. 2010, 37, L13402. [Google Scholar] [CrossRef]
  4. Jobbágy, E.G.; Jackson, R.B. The vertical distribution of soil organic carbon and its relation to climate and vegetation. Ecol. Appl. 2000, 10, 423–436. [Google Scholar] [CrossRef]
  5. Battle, M.; Bender, M.; Tans, P.P.; White, J.; Ellis, J.; Conway, T.; Francey, R. Global carbon sinks and their variability inferred from atmospheric O2 and δ13C. Science 2000, 287, 2467–2470. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Grosse, G.; Harden, J.; Turetsky, M.; McGuire, A.D.; Camill, P.; Tarnocai, C.; Frolking, S.; Schuur, E.A.; Jorgenson, T.; Marchenko, S. Vulnerability of high-latitude soil organic carbon in North America to disturbance. J. Geophys. Res. Biogeosci. 2011, 116, G00K06. [Google Scholar] [CrossRef] [Green Version]
  7. Turetsky, M.R.; Kane, E.S.; Harden, J.W.; Ottmar, R.D.; Manies, K.L.; Hoy, E.; Kasischke, E.S. Recent acceleration of biomass burning and carbon losses in Alaskan forests and peatlands. Nat. Geosci. 2011, 4, 27–31. [Google Scholar] [CrossRef]
  8. Davidson, E.A.; Janssens, I.A. Temperature sensitivity of soil carbon decomposition and feedbacks to climate change. Nature 2006, 440, 165–173. [Google Scholar] [CrossRef] [PubMed]
  9. McGuire, A.; Hayes, D.J.; Kicklighter, D.W.; Manizza, M.; Zhuang, Q.; Chen, M.; Follows, M.J.; Gurney, K.; Mcclelland, J.W.; Melillo, J. An analysis of the carbon balance of the Arctic Basin from 1997 to 2006. Tellus B Chem. Phys. Meteorol. 2010, 62, 455–474. [Google Scholar] [CrossRef] [Green Version]
  10. Jones, M.C.; Yu, Z. Rapid deglacial and early Holocene expansion of peatlands in Alaska. Proc. Natl. Acad. Sci. USA 2010, 107, 7347–7352. [Google Scholar] [CrossRef] [Green Version]
  11. Frolking, S.; Roulet, N.T.; Tuittila, E.; Bubier, J.L.; Quillet, A.; Talbot, J.; Richard, P. A new model of Holocene peatland net primary production, decomposition, water balance, and peat accumulation. Earth Syst. Dyn. 2010, 1, 1–21. [Google Scholar] [CrossRef] [Green Version]
  12. Rydin, H.; Jeglum, J.K. The Biology of Peatlands, 2nd ed.; Oxford University Press: Oxford, UK, 2013. [Google Scholar]
  13. Clymo, R.; Bryant, C. Diffusion and mass flow of dissolved carbon dioxide, methane, and dissolved organic carbon in a 7-m deep raised peat bog. Geochim. Cosmochim. Acta 2008, 72, 2048–2066. [Google Scholar] [CrossRef]
  14. Tfaily, M.M.; Wilson, R.M.; Cooper, W.T.; Kostka, J.E.; Hanson, P.; Chanton, J.P. Vertical Stratification of Peat Pore Water Dissolved Organic Matter Composition in a Peat Bog in Northern Minnesota. J. Geophys. Res. Biogeosci. 2018, 123, 479–494. [Google Scholar] [CrossRef]
  15. Loisel, J.; Garneau, M.; Hélie, J.-F. Sphagnum δ13C values as indicators of palaeohydrological changes in a peat bog. Holocene 2010, 20, 285–291. [Google Scholar] [CrossRef]
  16. 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]
  17. Asada, T.; Warner, B.G.; Aravena, R. Nitrogen isotope signature variability in plant species from open peatland. Aquat. Bot. 2005, 82, 297–307. [Google Scholar] [CrossRef]
  18. Laine, J.; Vasander, H. Ecology and vegetation gradients of peatlands. Peatl. Finl. 1996, 10, 20. [Google Scholar]
  19. Hugelius, G.; Loisel, J.; Chadburn, S.; Jackson, R.B.; Jones, M.; MacDonald, G.; Marushchak, M.; Olefeldt, D.; Packalen, M.; Siewert, M.B.; et al. Large stocks of peatland carbon and nitrogen are vulnerable to permafrost thaw. Proc. Natl. Acad. Sci. USA 2020, 117, 20438–20446. [Google Scholar] [CrossRef] [PubMed]
  20. Ives, S.L.; Sullivan, P.F.; Dial, R.; Berg, E.E.; Welker, J.M. CO2 exchange along a hydrologic gradient in the Kenai Lowlands, AK: Feedback implications of wetland drying and vegetation succession. Ecohydrology 2013, 6, 38–50. [Google Scholar] [CrossRef]
  21. Punning, J.M.; Toots, M.; Vaikmäe, R. Oxygen-18 in Estonian Natural Waters. Isot. Isot. Environ. Health Stud. 1987, 23, 232–234. [Google Scholar] [CrossRef]
  22. Klein, E.S.; Booth, R.K.; Yu, Z.; Mark, B.G.; Stansell, N.D. Hydrology-mediated differential response of carbon accumulation to late Holocene climate change at two peatlands in Southcentral Alaska. Quat. Sci. Rev. 2013, 64, 61–75. [Google Scholar] [CrossRef]
  23. Berger, A.; Loutre, M.F. Insolation values for the climate of the last 10 million years. Quat. Sci. Rev. 1991, 10, 297–317. [Google Scholar] [CrossRef]
  24. Paal, J. Rare and threatened plant communities of Estonia. Biodivers. Conserv. 1998, 7, 1027–1049. [Google Scholar] [CrossRef]
  25. Kalm, V. Pleistocene chronostratigraphy in Estonia, southeastern sector of the Scandinavian glaciation. Quat. Sci. Rev. 2006, 25, 960–975. [Google Scholar] [CrossRef]
  26. Masing, V.; Paal, J.; Kuresoo, A. Biodiversity of Estonian wetlands. Biodivers. Wetl. Assess. Funct. Conserv. 2000, 1, 259–279. [Google Scholar]
  27. Paal, J. Estonian mires. Mires—Sib. Tierra Del Fuego. Stapfia85 Zugleich Kat. Der OO. Landesmuseen Neue Ser. 2005, 35, 117–146. [Google Scholar]
  28. Jaagus, J. The impact of climate change on the snow cover pattern in Estonia. Clim. Change 1997, 36, 65–77. [Google Scholar] [CrossRef]
  29. Jaagus, J. Climatic changes in Estonia during the second half of the 20th century in relationship with changes in large-scale atmospheric circulation. Theor. Appl. Climatol. 2006, 83, 77–88. [Google Scholar] [CrossRef]
  30. Bladé, I.; Liebmann, B.; Fortuny, D.; Oldenborgh, G. Observed and simulated impacts of the summer NAO in Europe: Implications for projected drying in the Mediterranean region. Clim. Dyn. 2012, 39, 709–727. [Google Scholar] [CrossRef]
  31. Mätlik, O.; Post, P. Synoptic Weather Types That Have Caused Heavy Precipitation in Estonia in the Period 1961-2005. Est. J. Eng. 2008, 14, 195–208. [Google Scholar] [CrossRef] [Green Version]
  32. Wassenaar, L.I.; Coplen, T.B.; Aggarwal, P.K. Approaches for Achieving Long-Term Accuracy and Precision of δ18O and δ2H for Waters Analyzed using Laser Absorption Spectrometers. Environ. Sci. Technol. 2014, 48, 1123–1131. [Google Scholar] [CrossRef] [PubMed]
  33. Coplen, T.B. Guidelines and recommended terms for expression of stable-isotope-ratio and gas-ratio measurement results. Rapid Commun. Mass Spectrom. 2011, 25, 2538–2560. [Google Scholar] [CrossRef]
  34. Reimer, P.J.; Austin, W.E.; Bard, E.; Bayliss, A.; Blackwell, P.G.; Ramsey, C.B.; Butzin, M.; Cheng, H.; Edwards, R.L.; Friedrich, M. The IntCal20 Northern Hemisphere radiocarbon age calibration curve (0–55 cal kBP). Radiocarbon 2020, 62, 725–757. [Google Scholar] [CrossRef]
  35. Blaauw, M.; Christen, J.A. Flexible paleoclimate age-depth models using an autoregressive gamma process. Bayesian Anal. 2011, 6, 457–474. [Google Scholar] [CrossRef]
  36. Laine, J.; Harju, P.; Timonen, T.; Laine, A.; Tuittila, E.-S.; Minkkinen, K.; Vasander, H. The Intricate Beauty of Sphagnum Mosses—A Finnish Guide to Identification; University of Helsinki, Department of Forest Ecology: Helsinki, Finland, 2009. [Google Scholar]
  37. Loader, N.; McCarroll, D.; van der Knaap, W.O.; Robertson, I.; Gagen, M. Characterizing carbon isotopic variability in Sphagnum. Holocene 2007, 17, 403–410. [Google Scholar] [CrossRef]
  38. Dean, W.E. Determination of carbonate and organic matter in calcareous sediments and sedimentary rocks by loss on ignition: Comparison with other methods. J. Sediment. Petrol. 1974, 44, 242–248. [Google Scholar]
  39. Heiri, O.; Lotter, A.F.; Lemcke, G. Loss on ignition as a method for estimating organic and carbonate content in sediments: Reproducibility and comparability of results. J. Paleolimnol. 2001, 25, 101–110. [Google Scholar] [CrossRef]
  40. Boyle, J. A comparison of two methods for estimating the organic matter content of sediments. J. Paleolimnol. 2004, 31, 125–127. [Google Scholar]
  41. Loisel, J.; Yu, Z.; Beilman, D.W.; Camill, P.; Alm, J.; Amesbury, M.J.; Anderson, D.; Andersson, S.; Bochicchio, C.; Barber, K.; et al. A database and synthesis of northern peatland soil properties and Holocene carbon and nitrogen accumulation. Holocene 2014, 24, 1028–1042. [Google Scholar] [CrossRef]
  42. Rodionov, S.N. A sequential algorithm for testing climate regime shifts. Geophys. Res. Lett. 2004, 31, L09204. [Google Scholar] [CrossRef] [Green Version]
  43. Lê, S.; Josse, J.; Husson, F. FactoMineR: An R package for multivariate analysis. J. Stat. Softw. 2008, 25, 1–18. [Google Scholar] [CrossRef] [Green Version]
  44. Seppä, H.; Birks, H.J.B. July mean temperature and annual precipitation trends during the Holocene in the Fennoscandian tree-line area: Pollen-based climate reconstructions. Holocene 2001, 11, 527–539. [Google Scholar] [CrossRef] [Green Version]
  45. Sillasoo, U.; Mauquoy, D.; Blundell, A.; Charman, D.A.N.; Blaauw, M.; Daniell, J.R.G.; Toms, P.; Newberry, J.; Chambers, F.M.; Karofeld, E. Peat multi-proxy data from Männikjärve bog as indicators of late Holocene climate changes in Estonia. Boreas 2007, 36, 20–37. [Google Scholar] [CrossRef]
  46. Hill, K. Late Holocene Climate Controls on Carbon Dynamics at Teringi Bog, Estonia; Northern Illinois University: DeKalb, IL, USA, 2016. [Google Scholar]
  47. Ménot, G.; Burns, S.J. Carbon isotopes in ombrogenic peat bog plants as climatic indicators: Calibration from an altitudinal transect in Switzerland. Org. Geochem. 2001, 32, 233–245. [Google Scholar] [CrossRef]
  48. Dansgaard, W. Stable isotopes in precipitation. Tellus 1964, 16, 436–468. [Google Scholar] [CrossRef]
  49. Rozanski, K.; Araquás-Araquás, L.; Gonfiantini, R. Relation Between Long-Term Trends of Oxygen-18 Isotope Composition of Precipitation and Climate. Science 1992, 258, 981–985. [Google Scholar] [CrossRef]
  50. Baldini, L.M.; McDermott, F.; Foley, A.M.; Baldini, J.U.L. Spatial variability in the European winter precipitation δ18O-NAO relationship: Implications for reconstructing NAO-mode climate variability in the Holocene. Geophys. Res. Lett. 2008, 35, L04709. [Google Scholar] [CrossRef]
  51. St Amour, N.A.; Hammarlund, D.; Edwards, T.W.; Wolfe, B.B. New insights into Holocene atmospheric circulation dynamics in central Scandinavia inferred from oxygen-isotope records of lake-sediment cellulose. Boreas 2010, 39, 770–782. [Google Scholar] [CrossRef]
  52. Hammarlund, D.; Björck, S.; Buchardt, B.; Israelson, C.; Thomsen, C.T. Rapid hydrological changes during the Holocene revealed by stable isotope records of lacustrine carbonates from Lake Igelsjön, southern Sweden. Quat. Sci. Rev. 2003, 22, 353–370. [Google Scholar] [CrossRef]
  53. Heikkilä, M.; Edwards, T.W.D.; Seppä, H.; Sonninen, E. Sediment isotope tracers from Lake Saarikko, Finland, and implications for Holocene hydroclimatology. Quat. Sci. Rev. 2010, 29, 2146–2160. [Google Scholar] [CrossRef]
  54. Luken, J.O. Zonation of Sphagnum Mosses: Interactions among Shoot Growth, Growth Form, and Water Balance. Bryologist 1985, 88, 374–379. [Google Scholar] [CrossRef]
  55. Robroek, B.J.M.; Limpens, J.; Breeuwer, A.; Schouten, M.G.C. Effects of water level and temperature on performance of four Sphagnum mosses. Plant Ecol. 2007, 190, 97–107. [Google Scholar] [CrossRef]
  56. Hájek, T.; Beckett, R.P. Effect of water content components on desiccation and recovery in Sphagnum mosses. Ann. Bot. 2008, 101, 165–173. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Payette, S. Late-Holocene development of subarctic ombrotrophic peatlands: Allogenic and autogenic succession. Ecology 1988, 69, 516–531. [Google Scholar] [CrossRef]
  58. Rastogi, A.; Antala, M.; Gąbka, M.; Rosadziński, S.; Stróżecki, M.; Brestic, M.; Juszczak, R. Impact of warming and reduced precipitation on morphology and chlorophyll concentration in peat mosses (Sphagnum angustifolium and S. fallax). Sci. Rep. 2020, 10, 8592. [Google Scholar] [CrossRef]
  59. Moschen, R.; Kühl, N.; Rehberger, I.; Lücke, A. Stable carbon and oxygen isotopes in sub-fossil Sphagnum: Assessment of their applicability for palaeoclimatology. Chem. Geol. 2009, 259, 262–272. [Google Scholar] [CrossRef]
  60. Ménot-Combes, G.; Burns, S.J.; Leuenberger, M. Variations of 18O/16O in plants from temperate peat bogs (Switzerland): Implications for paleoclimatic studies. Earth Planet. Sci. Lett. 2002, 202, 419–434. [Google Scholar] [CrossRef]
  61. Melillo, J.M.; Aber, J.D.; Linkins, A.E.; Ricca, A.; Fry, B.; Nadelhoffer, K.J. Carbon and nitrogen dynamics along the decay continuum: Plant litter to soil organic matter. Plant. Soil 1989, 115, 189–198. [Google Scholar] [CrossRef]
  62. Malmer, N.; Svensson, G.; Wallén, B. Mass balance and nitrogen accumulation in hummocks on a South Swedish bog during the late Holocene. Ecography 1997, 20, 535–549. [Google Scholar] [CrossRef]
  63. Kuhry, P.; Vitt, D.H. Fossil Carbon/Nitrogen Ratios as a Measure of Peat Decomposition. Ecology 1996, 77, 271–275. [Google Scholar] [CrossRef]
  64. Ohlson, M.; Dahlberg, B. Rate of peat increment in hummock and lawn communities on Swedish mires during the last 150 years. Oikos 1991, 61, 369–378. [Google Scholar] [CrossRef]
  65. Berger, A.L. Long-term variations of caloric insolation resulting from the Earth’s orbital elements. Quat. Res. 1978, 9, 139–167. [Google Scholar] [CrossRef]
  66. Seppä, H.; Poska, A. Holocene annual mean temperature changes in Estonia and their relationship to solar insolation and atmospheric circulation patterns. Quat. Res. 2004, 61, 22–31. [Google Scholar] [CrossRef]
  67. Walker, M.J.C.; Berkelhammer, M.; Björck, S.; Cwynar, L.C.; Fisher, D.A.; Long, A.J.; Lowe, J.J.; Newnham, R.M.; Rasmussen, S.O.; Weiss, H. Formal subdivision of the Holocene Series/Epoch: A Discussion Paper by a Working Group of INTIMATE (Integration of ice-core, marine and terrestrial records) and the Subcommission on Quaternary Stratigraphy (International Commission on Stratigraphy). J. Quat. Sci. 2012, 27, 649–659. [Google Scholar] [CrossRef]
  68. Morley, A.; Rosenthal, Y.; deMenocal, P. Ocean-atmosphere climate shift during the mid-to-late Holocene transition. Earth Planet. Sci. Lett. 2014, 388, 18–26. [Google Scholar] [CrossRef] [Green Version]
  69. Cowling, S.A.; Sykes, M.T.; Bradshaw, R.H.W. Palaeovegetation-model comparisons, climate change and tree succession in Scandinavia over the past 1500 years. J. Ecol. 2001, 89, 227–236. [Google Scholar] [CrossRef] [Green Version]
  70. Poska, A.; Saarse, L.; Koppel, K.; Nielsen, A.B.; Avel, E.; Vassiljev, J.; Väli, V. The Verijärv area, South Estonia over the last millennium: A high resolution quantitative land-cover reconstruction based on pollen and historical data. Rev. Palaeobot. Palynol. 2014, 207, 5–17. [Google Scholar] [CrossRef]
  71. Helama, S.; Meriläinen, J.; Tuomenvirta, H. Multicentennial megadrought in northern Europe coincided with a global El Niño–Southern Oscillation drought pattern during the Medieval Climate Anomaly. Geology 2009, 37, 175–178. [Google Scholar] [CrossRef] [Green Version]
  72. Luoto, T.P.; Helama, S. Palaeoclimatological and palaeolimnological records from fossil midges and tree-rings: The role of the North Atlantic Oscillation in eastern Finland through the Medieval Climate Anomaly and Little Ice Age. Quat. Sci. Rev. 2010, 29, 2411–2423. [Google Scholar] [CrossRef]
  73. Tiljander, M.I.A.; Saarnisto, M.; Ojala, A.E.K.; Saarinen, T. A 3000-year palaeoenvironmental record from annually laminated sediment of Lake Korttajarvi, central Finland. Boreas 2003, 32, 566–577. [Google Scholar] [CrossRef]
  74. Tarand, A.; Nordli, P.Ø. The Tallinn Temperature Series Reconstructed Back Half a Millennium by Use of Proxy Data. In The Iceberg in the Mist: Northern Research in Pursuit of a “Little Ice Age”; Ogilvie, A.E.J., Jónsson, T., Eds.; Springer: Dordrecht, The Netherlands, 2001; pp. 189–199. [Google Scholar]
  75. Stouffer, R.J.; Manabe, S. Response of a Coupled Ocean–Atmosphere Model to Increasing Atmospheric Carbon Dioxide: Sensitivity to the Rate of Increase. J. Clim. 1999, 12, 2224–2237. [Google Scholar] [CrossRef] [Green Version]
  76. BACC. Assessment of Climate Change for the Baltic Sea Basin; Germany, 2008; Available online: https://link.springer.com/book/10.1007/978-3-540-72786-6 (accessed on 1 August 2021).
  77. EEA. European Environment Agency: National Adaptation Policy Processes in European Countries 2014. 2014; Available online: https://www.eea.europa.eu/publications/national-adaptation-policy-processes (accessed on 1 August 2021).
  78. IPCC. Climate Change 2014: Impacts, Adaptation, and Vulnerability. Part. B: Regional Aspects. Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Barros, V.R.C.B., Field, D.J., Dokken, M.D., Mastrandrea, K.J., Mach., T.E., Bilir, M., Chatterjee, K.L., Ebi, Y.O., Estrada, R.C., Genova, B., et al., Eds.; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2014. [Google Scholar]
  79. Kont, A.; Jaagus, J.; Aunap, R. Climate change scenarios and the effect of sea-level rise for Estonia. Glob. Planet. Change 2003, 36, 1–15. [Google Scholar] [CrossRef]
  80. Kauker, F.; Meier, H.E.M. Modeling decadal variability of the Baltic Sea: 1. Reconstructing atmospheric surface data for the period 1902–1998. J. Geophys. Res. Ocean. 2003, 108, 3267. [Google Scholar] [CrossRef] [Green Version]
  81. Meier, H.E.M.; Kauker, F. Modeling decadal variability of the Baltic Sea: 2. Role of freshwater inflow and large-scale atmospheric circulation for salinity. J. Geophys. Res. Ocean. 2003, 108, 3368. [Google Scholar] [CrossRef] [Green Version]
  82. Meier, H.E.M. Baltic Sea climate in the late twenty-first century: A dynamical downscaling approach using two global models and two emission scenarios. Clim. Dyn. 2006, 27, 39–68. [Google Scholar] [CrossRef]
  83. Harrison, S.P.; Bartlein, P.J.; Izumi, K.; Li, G.; Annan, J.; Hargreaves, J.; Braconnot, P.; Kageyama, M. Evaluation of CMIP5 palaeo-simulations to improve climate projections. Nat. Clim. Change 2015, 5, 735–743. [Google Scholar] [CrossRef]
  84. Rutgersson, A.; Jaagus, J.; Schenk, F.; Stendel, M. Observed changes and variability of atmospheric parameters in the Baltic Sea region during the last 200 years. Clim. Res. 2014, 61, 177–190. [Google Scholar] [CrossRef] [Green Version]
  85. Graham, L.P. Climate Change Effects on River Flow to the Baltic Sea. AMBIO A J. Hum. Environ. 2004, 33, 235–241. [Google Scholar] [CrossRef]
  86. Collins, M.; Knutti, R.; Arblaster, J.; Dufresne, J.-L.; Fichefet, T.; Friedlingstein, P.; Gao, X.; Gutowski, W.J.; Johns, T.; Krinner, G. Long-term climate change: Projections, commitments and irreversibility. In Climate Change 2013-The Physical Science Basis: Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK, 2013; pp. 1029–1136. [Google Scholar]
Figure 1. Location map of the study site panels (A) and orthophoto of Teringi Bog panels (B). The yellow star indicates the sampling well and coring location (57°58′33.56″ N, 25°33′30.90″ E).
Figure 1. Location map of the study site panels (A) and orthophoto of Teringi Bog panels (B). The yellow star indicates the sampling well and coring location (57°58′33.56″ N, 25°33′30.90″ E).
Quaternary 05 00008 g001
Figure 2. Typical vegetation at Teringi Bog and the surrounding margins. Panel (A) represents the boundary between the open bog and forest. Panels (BD) are at the coring location in the center of the bog (see Figure 1). Panel (E) represents trees growing near the center of the bog. Panel (F) is the forest on the western edge of Teringi.
Figure 2. Typical vegetation at Teringi Bog and the surrounding margins. Panel (A) represents the boundary between the open bog and forest. Panels (BD) are at the coring location in the center of the bog (see Figure 1). Panel (E) represents trees growing near the center of the bog. Panel (F) is the forest on the western edge of Teringi.
Quaternary 05 00008 g002
Figure 3. Panel (A): Estonia Global Network of Isotopes in Precipitation (GNIP) precipitation-weighted δ18Op data for Tartu shown with average monthly values for 2014 through 2018, with the annual mean shown with a gray solid line. Panel (B): Monthly precipitation and temperature values at Valga, Estonia for the period spanning 1945 to 2020 (Estonian Weather Service).
Figure 3. Panel (A): Estonia Global Network of Isotopes in Precipitation (GNIP) precipitation-weighted δ18Op data for Tartu shown with average monthly values for 2014 through 2018, with the annual mean shown with a gray solid line. Panel (B): Monthly precipitation and temperature values at Valga, Estonia for the period spanning 1945 to 2020 (Estonian Weather Service).
Quaternary 05 00008 g003
Figure 4. Automatic weather station from Karksi, Estonia (B,D,E) and levelogger data from Teringi Bog (A,C). The peat water temperature and water level were measured continuously from August, 2014 to July, 2021. The continuous automatic weather station measurements are from May 2015 to July 2021.
Figure 4. Automatic weather station from Karksi, Estonia (B,D,E) and levelogger data from Teringi Bog (A,C). The peat water temperature and water level were measured continuously from August, 2014 to July, 2021. The continuous automatic weather station measurements are from May 2015 to July 2021.
Quaternary 05 00008 g004
Figure 5. Modern water isotopes from Estonia plotted versus the global meteoric water line (GMWL). The Teringi Bog surface and well water samples plot near the GMWL, reflecting the isotopic composition of source precipitation and the water is a mix of both warm and cold season precipitation. Cold season precipitation samples have relatively low δ18O and δ2H values, while warm season values are higher. The regional evaporation line (REL), determined from closed-basin lakes in the area, is shown for comparison and plots well below the Teringi Bog surface and well water samples.
Figure 5. Modern water isotopes from Estonia plotted versus the global meteoric water line (GMWL). The Teringi Bog surface and well water samples plot near the GMWL, reflecting the isotopic composition of source precipitation and the water is a mix of both warm and cold season precipitation. Cold season precipitation samples have relatively low δ18O and δ2H values, while warm season values are higher. The regional evaporation line (REL), determined from closed-basin lakes in the area, is shown for comparison and plots well below the Teringi Bog surface and well water samples.
Quaternary 05 00008 g005
Figure 6. Age-depth model for Teringi Bog calculated with Bacon. The horizontal bars and grey lines represent the modeled calibration range of each 14C sample (see Table 1).
Figure 6. Age-depth model for Teringi Bog calculated with Bacon. The horizontal bars and grey lines represent the modeled calibration range of each 14C sample (see Table 1).
Quaternary 05 00008 g006
Figure 7. Sample leaf morphologies used for modern species identification (Hill, 2016). The top left panel is S. fuscum, top middle is S. capillifolium, top right is S. rubellum, bottom left is likely S. balticum, bottom middle is S. fallax and bottom right is S. pulchrum.
Figure 7. Sample leaf morphologies used for modern species identification (Hill, 2016). The top left panel is S. fuscum, top middle is S. capillifolium, top right is S. rubellum, bottom left is likely S. balticum, bottom middle is S. fallax and bottom right is S. pulchrum.
Quaternary 05 00008 g007
Figure 8. Peat core data. Solid black lines represent the mean values computed by regime-shift analysis.
Figure 8. Peat core data. Solid black lines represent the mean values computed by regime-shift analysis.
Quaternary 05 00008 g008
Figure 9. Teringi Bog data plotted versus regional paleoclimate records mentioned in text. Summer insolation (A) [21] decreased through the late Holocene as there was an overall trend to lower temperatures recorded in lake sediment records (B) [44]. Fluctuations in water table depth recorded at Männikjärve Bog in Estonia (C) [45] generally correspond to changes in hydroclimate recorded in the carbon (D) and nitrogen (E) isotope records from Teringi Bog. Note the inverted vertical axes for δ13C and δ15N. The carbon accumulation rates at Teringi Bog (F) were generally higher under wetter conditions during the late Holocene.
Figure 9. Teringi Bog data plotted versus regional paleoclimate records mentioned in text. Summer insolation (A) [21] decreased through the late Holocene as there was an overall trend to lower temperatures recorded in lake sediment records (B) [44]. Fluctuations in water table depth recorded at Männikjärve Bog in Estonia (C) [45] generally correspond to changes in hydroclimate recorded in the carbon (D) and nitrogen (E) isotope records from Teringi Bog. Note the inverted vertical axes for δ13C and δ15N. The carbon accumulation rates at Teringi Bog (F) were generally higher under wetter conditions during the late Holocene.
Quaternary 05 00008 g009
Figure 10. PCA of Teringi Bog core data and regional paleoclimate records of water table fluctuations [45] and July temperatures recorded in the Baltic region [44].
Figure 10. PCA of Teringi Bog core data and regional paleoclimate records of water table fluctuations [45] and July temperatures recorded in the Baltic region [44].
Quaternary 05 00008 g010
Table 1. Radiometric ages used in study. All samples were measured on Sphagnum stems.
Table 1. Radiometric ages used in study. All samples were measured on Sphagnum stems.
14C AgeAnalytical Uncertainty ± 1σMedian Calibrated Age (cal yr BP)Minimum Calibrated Age (cal yr BP) 2σMaximum Calibrated Age (cal yr BP) 2σWeighted Mean Modeled Age (cal yr BP)Minimum Modeled Age (cal yr BP) 2σMaximum Modeled Age (cal yr BP) 2σ
–651 –65–68–62
140151095811817263271
154515140513691419140113511488
202515196319252000196819012043
238020239223462465244223522654
307515329132323358329732183363
381020419541454252420040974347
Table 2. Sphagnum species identification at each sampled depth. A ‘1’ indicates that a species was the most common at that interval. A ‘2’ indicates that the species was present, but not the most common.
Table 2. Sphagnum species identification at each sampled depth. A ‘1’ indicates that a species was the most common at that interval. A ‘2’ indicates that the species was present, but not the most common.
Section AcutifoliaSection CuspidataSection SubsecundaSection Sphagnum
Depth (cm)Age (cal yr BP)S. fuscumS. capilli-foliumS. rubellumS. russowiS. pulchrumS. fallaxS. angusti-foliumS.balticumS. lind-bergiispp.S. medium
32182 22
422662 1
52352 12 2
624391 2
725251
8261322 2
9269922 2 2
1027862 1
112870 1
12295721
1321042 1
142113012 2
15212201 2
162131122 1
17213971 2
182146012 2 2
19215211 2
20215811
212164222 1
222170111 2
232176211 2
242182321
252188521
262194612 2
272201112
282207812
292214412
302221012
312227612 1
32223431
332241112 2
34224871 2 2
352257512 2
362266312
372275021
382283721
392292321 2
402301112
41230991 2
422318712
432327512
44233701
452346912
46235671
47236651 1
48237641
49238621 1
502395912 1
51240551
52241541
Table 3. Correlation coefficients for PCA (p values are all <0.001).
Table 3. Correlation coefficients for PCA (p values are all <0.001).
VariableDim. 1Dim. 2Dim. 3
%N0.938−0.0160.048
C:N−0.9180.047−0.072
õ15N0.5570.439−0.531
Organic C flux0.1410.8840.365
Peat accum. Rate−0.3470.6570.446
Bulk density0.4330.6540.111
Water table depth0.316−0.6240.255
õ13C−0.2640.274−0.692
July temp0.072−0.2460.504
%C0.461−0.168−0.13
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Stansell, N.D.; Klein, E.S.; Hill, K.; Terasmaa, J.; Dodd, J.; Boes, M.; Eensalu, M.; Fortney, C.; Fritts, A.; Garcia, R.; et al. A Late Holocene Stable Isotope and Carbon Accumulation Record from Teringi Bog in Southern Estonia. Quaternary 2022, 5, 8. https://doi.org/10.3390/quat5010008

AMA Style

Stansell ND, Klein ES, Hill K, Terasmaa J, Dodd J, Boes M, Eensalu M, Fortney C, Fritts A, Garcia R, et al. A Late Holocene Stable Isotope and Carbon Accumulation Record from Teringi Bog in Southern Estonia. Quaternary. 2022; 5(1):8. https://doi.org/10.3390/quat5010008

Chicago/Turabian Style

Stansell, Nathan D., Eric S. Klein, Kristyn Hill, Jaanus Terasmaa, Justin Dodd, Maxwell Boes, Mariliis Eensalu, Carolyn Fortney, Annabella Fritts, Roxana Garcia, and et al. 2022. "A Late Holocene Stable Isotope and Carbon Accumulation Record from Teringi Bog in Southern Estonia" Quaternary 5, no. 1: 8. https://doi.org/10.3390/quat5010008

APA Style

Stansell, N. D., Klein, E. S., Hill, K., Terasmaa, J., Dodd, J., Boes, M., Eensalu, M., Fortney, C., Fritts, A., Garcia, R., Price, B., & Swanson, B. (2022). A Late Holocene Stable Isotope and Carbon Accumulation Record from Teringi Bog in Southern Estonia. Quaternary, 5(1), 8. https://doi.org/10.3390/quat5010008

Article Metrics

Back to TopTop