Search

Peatlands are one of the largest terrestrial carbon sinks on the planet, yet little is known about the carbon accumulation rates (CARs) of mountainous peatlands. The long-term variability in the size of the associated carbon sink and its drivers remain largely unconstrained, especially when the long-term anthropogenic impact is also considered. Here, we present a composite CAR record of nine peatlands from central–eastern Europe (Romania and Serbia) detailing variability in the rates of carbon accumulation during the Holocene. We show examples of extremely high long-term rates of carbon accumulation (LORCA>120 gCm-2yr-1), indicating that mountain peatlands constitute an efficient regional carbon sink at times. By comparing our data to modelled palaeoclimatic indices and to measures of anthropogenic impact we disentangle the drivers of peat carbon accumulation in the area. Variability in early- and mid-Holocene CARs is linked to hydroclimatic controls, with high CARs occurring during the early Holocene and lower CARs associated with the transition to cooler and moister mid-Holocene conditions. By contrast, after 4000 years (calibrated) before present (years BP), the trends in CARs indicate a divergence from hydroclimate proxies, suggesting that other processes became the dominant drivers of peat CARs. We propose that enhanced erosion following tree cover reduction as well as increased rates of long-distance atmospheric dust fallout might have played a role, as both processes would result in enhanced mineral and nutrient supply to bog surfaces, stimulating peatland productivity. Surprisingly though, for the last 1000 years, reconstructed temperature is significantly correlated with CARs, with rising temperatures linked to higher CARs. Under future climate conditions, which are predicted to be warmer in the region, we predict that peat growth may expand but that this is entirely dependent upon the scale of human impact directly affecting the sensitive hydrological budget of these peatlands.

Reconstructions of dust flux have been used to produce valuable global records of changes in atmospheric circulation and aridity. These studies have highlighted the importance of atmospheric dust in marine and terrestrial biogeochemistry and nutrient cycling. By investigating a 10 800-year-long paleoclimate archive from the Eastern Carpathians (Romania) we present the first peat record of changing dust deposition over the Holocene for the Carpathian–Balkan region. Using qualitative (X-ray fluorescence (XRF) core scanning) and quantitative inductively coupled plasma optical emission spectrometer(ICP-OES) measurements of lithogenic (K, Si, Ti) elements, we identify 10 periods of major dust deposition between 9500–9200, 8400–8100, 7720–7250, 6350–5950, 5450–5050, 4130–3770, 3450–2850, 2000–1450, 800–620, and 60 cal yr BP to present. In addition, we used testate amoeba assemblages preserved within the peat to infer local palaeohydroclimatic conditions. Our record highlights several discrepancies between eastern and western European dust depositional records and the impact of highly complex hydrological regimes in the Carpathian region. Since 6100 cal yr BP, we find that the geochemical indicators of dust flux have become uncoupled from the local hydrology. This coincides with the appearance of millennial-scale cycles in the dust input and changes in geochemical composition of dust. We suggest that this is indicative of a shift in dust provenance from local–regional (likely loess-related) to distal (Saharan) sources, which coincide with the end of the African Humid Period and the onset of Saharan desertification.

Here we present a speleothem isotope record (POM2) from Ascunsă Cave (Romania) that provides new data on past climate changes in the Carpathian–Balkan region from 8.2 ka until the present. This paper describes an approach to constrain the effect of temperature changes on calcite δ18O values in stalagmite POM2 over the course of the middle Holocene (6–4 ka), and across the 8.2 and 3.2 ka rapid climate change events. Independent pollen temperature reconstructions are used to this purpose. The approach combines the temperature-dependent isotope fractionation of rain water during condensation and fractionation resulting from calcite precipitation at the given cave temperature. The only prior assumptions are that pollen-derived average annual temperature reflects average cave temperature, and that pollen-derived coldest and warmest month temperatures reflect the range of condensation temperatures of rain above the cave site. This approach constrains a range of values between which speleothem δ18O changes should be found if controlled only by surface temperature variations at the cave site. Deviations of the change in δ18Ocspel values from the calculated temperature-constrained range of change are interpreted towards large-scale variability of climate–hydrology. Following this approach, we show that an additional ∼0.6‰ enrichment of δ18Oc in the POM2 stalagmite was caused by changing hydrological patterns in SW Romania across the middle Holocene, most likely comprising local evaporation from the soil and an increase in Mediterranean moisture δ18O. Further, by extending the calculations to other speleothem records from around the entire Mediterranean basin, it appears that all eastern Mediterranean speleothems recorded a similar isotopic enrichment due to changing hydrology, whereas all changes recorded in speleothems from the western Mediterranean are fully explained by temperature variation alone. This highlights a different hydrological evolution between the two sides of the Mediterranean. Our results also demonstrate that during the 8.2 ka event, POM2 stable isotope data essentially fit the temperature-constrained isotopic variability. In the case of the 3.2 ka event, an additional climate-related hydrological factor is more evident. This implies a different rainfall pattern in the Southern Carpathian region during this event at the end of the Bronze Age.

An accurate and coherent chronological framework is essential for the interpretation of climatic and environmental records obtained from deep polar ice cores. Until now, one common ice core age scale had been developed based on an inverse dating method (Datice), combining glaciological modelling with absolute and stratigraphic markers between 4 ice cores covering the last 50 ka (thousands of years before present) (Lemieux-Dudon et al., 2010). In this paper, together with the companion paper of Veres et al. (2013), we present an extension of this work back to 800 ka for the NGRIP, TALDICE, EDML, Vostok and EDC ice cores using an improved version of the Datice tool. The AICC2012 (Antarctic Ice Core Chronology 2012) chronology includes numerous new gas and ice stratigraphic links as well as improved evaluation of background and associated variance scenarios. This paper concentrates on the long timescales between 120–800 ka. In this framework, new measurements of δ18Oatm over Marine Isotope Stage (MIS) 11–12 on EDC and a complete δ18Oatm record of the TALDICE ice cores permit us to derive additional orbital gas age constraints. The coherency of the different orbitally deduced ages (from δ18Oatm, δO2/N2 and air content) has been verified before implementation in AICC2012. The new chronology is now independent of other archives and shows only small differences, most of the time within the original uncertainty range calculated by Datice, when compared with the previous ice core reference age scale EDC3, the Dome F chronology, or using a comparison between speleothems and methane. For instance, the largest deviation between AICC2012 and EDC3 (5.4 ka) is obtained around MIS 12. Despite significant modifications of the chronological constraints around MIS 5, now independent of speleothem records in AICC2012, the date of Termination II is very close to the EDC3 one.

The deep polar ice cores provide reference records commonly employed in global correlation of past climate events. However, temporal divergences reaching up to several thousand years (ka) exist between ice cores over the last climatic cycle. In this context, we are hereby introducing the Antarctic Ice Core Chronology 2012 (AICC2012), a new and coherent timescale developed for four Antarctic ice cores, namely Vostok, EPICA Dome C (EDC), EPICA Dronning Maud Land (EDML) and Talos Dome (TALDICE), alongside the Greenlandic NGRIP record. The AICC2012 timescale has been constructed using the Bayesian tool Datice (Lemieux-Dudon et al., 2010) that combines glaciological inputs and data constraints, including a wide range of relative and absolute gas and ice stratigraphic markers. We focus here on the last 120 ka, whereas the companion paper by Bazin et al. (2013) focuses on the interval 120–800 ka. Compared to previous timescales, AICC2012 presents an improved timing for the last glacial inception, respecting the glaciological constraints of all analyzed records. Moreover, with the addition of numerous new stratigraphic markers and improved calculation of the lock-in depth (LID) based on δ15N data employed as the Datice background scenario, the AICC2012 presents a slightly improved timing for the bipolar sequence of events over Marine Isotope Stage 3 associated with the seesaw mechanism, with maximum differences of about 600 yr with respect to the previous Datice-derived chronology of Lemieux-Dudon et al. (2010), hereafter denoted LD2010. Our improved scenario confirms the regional differences for the millennial scale variability over the last glacial period: while the EDC isotopic record (events of triangular shape) displays peaks roughly at the same time as the NGRIP abrupt isotopic increases, the EDML isotopic record (events characterized by broader peaks or even extended periods of high isotope values) reached the isotopic maximum several centuries before. It is expected that the future contribution of both other long ice core records and other types of chronological constraints to the Datice tool will lead to further refinements in the ice core chronologies beyond the AICC2012 chronology. For the time being however, we recommend that AICC2012 be used as the preferred chronology for the Vostok, EDC, EDML and TALDICE ice core records, both over the last glacial cycle (this study), and beyond (following Bazin et al., 2013). The ages for NGRIP in AICC2012 are virtually identical to those of GICC05 for the last 60.2 ka, whereas the ages beyond are independent of those in GICC05modelext (as in the construction of AICC2012, the GICC05modelext was included only via the background scenarios and not as age markers). As such, where issues of phasing between Antarctic records included in AICC2012 and NGRIP are involved, the NGRIP ages in AICC2012 should therefore be taken to avoid introducing false offsets. However for issues involving only Greenland ice cores, there is not yet a strong basis to recommend superseding GICC05modelext as the recommended age scale for Greenland ice cores.

We compare a variety of methods for estimating the gas/ice depth offset (Δdepth) at EPICA Dome C (EDC, East Antarctica). (1) Purely based on modelling efforts, Δdepth can be estimated combining a firn densification with an ice flow model. (2) The diffusive column height can be estimated from δ15N and converted to Δdepth using an ice flow model and assumptions about past average firn density and thickness of the convective zone. (3) Ice and gas synchronisation of the EDC ice core to the GRIP, EDML and TALDICE ice cores shifts the ice/gas offset problem into higher accumulation ice cores where it can be more accurately evaluated. (4) Finally, the bipolar seesaw hypothesis allows us to synchronise the ice isotopic record with the gas CH4 record, the later being taken as a proxy of Greenland temperature. The general agreement of method 4 with methods 2 and 3 confirms that the bipolar seesaw antiphase happened during the last 140 kyr. Applying method 4 to the deeper section of the EDC core confirms that the ice flow is complex and can help to improve our reconstruction of the thinning function and thus, of the EDC age scale. We confirm that method 1 overestimates the glacial Δdepth at EDC and we suggest that it is due to an overestimation of the glacial lock-in depth (LID) by the firn densification model. In contrast, we find that method 1 very likely underestimates Δdepth during Termination II, due either to an underestimated thinning function or to an underestimated LID. Finally, method 2 gives estimates within a few metres of methods 3 and 4 during the last deglacial warming, suggesting that the convective zone at Dome C cannot have been very large at this time, if it existed at all.

This study aims at refining the synchronisation between the EPICA Dome C (EDC) and Vostok ice cores in the time interval 0–145 kyr BP by using the volcanic signatures. 102 common volcanic events were identified by using continuous electrical conductivity (ECM), di-electrical profiling (DEP) and sulfate measurements while trying to minimize the distortion of the glaciological chronologies. This is an update and a continuation of previous works performed over the 0–45 kyr interval that provided 56 tie points to the ice core chronologies (Udisti et al., 2004). This synchronisation will serve to establish Antarctic Ice Core Chronology 2012, the next synchronised Antarctic dating. A change of slope in the EDC-depth/Vostok-depth diagram is probably related to a change of accumulation regime as well as to a change of ice thickness upstream of the Lake Vostok, but we did not invoke any significant temporal change of surface accumulation at EDC relative to Vostok. No significant phase difference is detected between the EDC and Vostok isotopic records, but depth shifts between the Vostok 3G and 5G ice cores prevent from looking at this problem accurately. Three possible candidates for the Toba volcanic super-eruption ~73 kyr ago are suggested in the Vostok and EDC volcanic records. Neither the ECM, DEP nor the sulfate fingerprints for these 3 events are significantly larger than many others in the records.