Abstract
The occurrence of heavy rainfall events is expected to undergo significant changes under increasing anthropogenic forcing. South-eastern Europe is reacting rapidly to such changes, therefore understanding and forecasting of precipitation variability is vital to better comprehending environmental changes in this area. Here we present a sub-decadal reconstruction of enhanced rainfall events for the past 2000 years from the Southern Carpathians, Romania using peat geochemistry. Five clear periods of enhanced rainfall are identified at 125–250, 600–900, 1050–1300, 1400–1575 and 1725–1980 CE. Significant runoff is observed during the second half of the Medieval Warm Period, whilst the Little Ice Age was characterised by significant variability. The North Atlantic Oscillation appears to be the main control on regional precipitation, but changes in solar irradiance also seem to play a significant role, together with the Siberian High. Comparison of the data presented here with model outputs confirms the ability of models to predict general trends, and major shifts, but highlights the complexity of the region’s hydrological history.
Similar content being viewed by others
Introduction
Understanding the causes and effects of wetter periods and associated hydro-meteorological events (heavy precipitation, flooding) is of great interest in climate science1, as the majority of modelling studies predict higher incidence of such events in the near future2. Heavy precipitation usually leads to flooding, the risk to human life, and an increase in economic disturbances1,3,4. However, comprehension of the spatial and temporal distribution, as well as the magnitude and driving mechanisms of high rainfall, runoff and flooding events is still incomplete. South-eastern Europe is one of the most rapidly reacting areas to current climate change2,5 and to understand the long-term hydroclimate variability in this region, high-resolution and well-dated records are needed. Within eastern Europe, projections suggest a shift of North Atlantic storm tracks northward, which will lead to more extreme weather patterns6,7, with higher risk of droughts2, and stronger cyclonic activity8. Since cyclonic activity is one of the major controls on northern hemisphere precipitation variability6,9, understanding its long-term behaviour and environmental impact is vital. However, projections are extremely uncertain, with models unable in most cases to reliably resolve storm track movements without validation against long-term observational data10.
For the most recent past, and especially the last century, data on rainfall variability in the Carpathian area, and its linkage to atmospheric pressure systems may be inferred from documentary sources11,12,13, and meteorological observations14,15. So far, in order to asses long-term hydrological changes and disentangle past rainfall variability research focused mainly on geochemical records of peat bogs and lake sediments from central-western Europe16,17,18,19,20,21,22. With some exceptions, comparable high-resolution, well-dated palaeohydrological records from eastern Europe are rare23,24,25,26, particularly records at yearly, or near-yearly resolution27.
Compared to the hydrology of central-western Europe, which is primarily influenced by Atlantic forcings, meteorological observations14,28,29 and proxy data25,26,30 show south-eastern Europe is periodically subjected to the North Atlantic, Siberian and Mediterranean atmospheric pressure systems. As a result, local to regional climate is seasonally and sometimes abruptly affected by shifts in these systems23,30,31,32. In particular, shifts in the strength of North Atlantic Oscillation (NAO) have been linked to changes in winter precipitation14,28,33, and winter temperature in the region14. Several paleoenvironmental records appear to confirm this connectivity has persisted for the past few millennia23,25,34,35. However, other records, particularly from the eastern Mediterranean, show no clear link to the NAO36,37. As such, further high-resolution palaeoclimate data from south-eastern Europe may strengthen our understanding of long-term patterns and amplitudes of change in the pressure systems mentioned above, as well as their environmental impact across the continent as a whole30.
Here we present a new high-resolution geochemical record of sediment input onto Sureanu bog (Southern Carpathians, Romania, Fig. 1) linked to precipitation variability in the catchment area, and associated runoff erosion for the last 2000 years. Such work builds on a previous, low-resolution study into the paleoclimatic history of this area, which demonstrated a clear link between variations in lithogenic element concentrations and land erosion.
Results and Discussion
When interpreting X-Ray Fluorescence Core-Scanning (XRF-CS) data, particularly with peat sediment, it is important to ensure the variations truly reflect changes in depositional geochemistry38. To ensure the XRF-CS data discussed here provide a reliable method of geochemical screening for the Sureanu peat record, comparison of ITRAX results with those from destructive Inductively Coupled Plasma Optical Emission Spectrometry (ICP-OES) has been made (Fig. 2). Visually, the XRF-CS data reflects the broad changes observed in the ICP-derived geochemical series (Fig. 2). By comparing the ICP and the XRF-CS data on the same age scale (at 20 year resolution), statistical comparisons may be attempted. Positive correlations significant at the 0.01 level (p-value < 0.001) may be observed for Fe (r2 = 0.209), Sr (r2 = 0.301) and Ti (r2 = 0.496) whilst for Zr the correlation is significant at the 0.05 level (r2 = 0.106, p-value = 0.037). These correlations, although low, are considerable when the constraints of the two methods are considered. ICP analysis provides the elemental concentration of a 1cm-wide sediment slice, and therefore encompasses 5 ITRAX data-point measurements. Further, ITRAX screening analyses only the uppermost fraction of the sample surface, and results may be impacted by variations in water content and surface roughness. Finally, the fact that analysis via ICP requires bringing the sample into solution introduces another uncertainty, that of digestion. Whilst all attempts have been made, including the addition of HF to ensure total digestion, it may not be possible in all cases. For Rb no such correlation is observed, likely a result of Rb concentrations approaching the ITRAX detection limit. However, correlation between all XRF-CS-derived datasets (Table 1) suggests records of all lithogenic elements may be used to reconstruct runoff related erosion in Sureanu peat bog’s catchment area (Fig. 2). For ease of comprehension, however, we focus here on the Ti record (Fig. 3e).
Sureanu peat bog is located at the base of a glacial cirque and is surrounded on three sides by steep slopes (Figs 1b and S125). These slopes display numerous mass-wasting related geomorphological features, including erosional channels, debris lobes, and exposed rock (Figs 1b and S1, S239). This means that mass-wasting activity and superficial runoff affecting the slopes will likely be channelled towards Iezerul Sureanu Lake and the adjacent Sureanu bog. Furthermore, any overtopping of the lake during highstands will normally result in mineral debris being washed onto the bog. Longman et al. 2017 discussed a number of mechanisms (e.g. avalanches, debris flows, runoff, or lake flooding) that may result in mass-wasting activity and the deposition of minerogenic matter (MM) onto Sureanu bog25, but the common feature is their linkage with high precipitation events. The resolution of the record means that each year is covered by one or two data points, and so should be considered indicative of all erosion occurring at the site in the previous half-year or one year. As a result, the rate at which MM is exported from the slope to the bog is inconsequential, and the record should not be considered indicative of individual events, but of an average value indicative of all erosion occurring in the preceding period of time. By whatever specific mechanism, periods of increased MM deposition within Sureanu bog may be considered representative of increased precipitation in the catchment area, and thus stronger runoff.
The comparability of all lithogenic elements as analysed via ITRAX indicate similar controls on their deposition (Fig. 2), and a dominant and geochemically distinct mineral source. This contrasts with what would be expected if trends in the lithogenic elements reflect various local and/or distal lithogenic sources, including significant input of long-distance transported dust. In such cases, shifts in relative geochemical makeup of the depositional events should be observed in different lithogenic elements, linked to, for example, changing (proximal to distal) dust sources40,41,42. Furthermore, typical average grain sizes within the MM-rich layers at Sureanu are larger (40–130 μm25) than would be expected from distal atmospheric dust deposition during the Holocene (~20 μm)43. Contribution of distal dust, as observed within Mohos peat bog in the Eastern Carpathians40 cannot be ruled out, but this is clearly overprinted by the contribution from local rocks at Sureanu. Evidence for such events inputting MM into the bog are visually documented, particularly in the uppermost 100 cm of the record, with a number of detrital layers, some containing pebbles >2 mm, observable25 (Fig. 3).
Erosional regimes are not dictated solely by precipitation, with human activity and especially deforestation potentially enhancing erosion by removing protective vegetation. However, there is evidence at Sureanu25, and other records in the Southern Carpathians44 that the onset of deforestation had begun near Sureanu earlier than the onset of regular minerogenic deposition. Additionally, the fact that steep slopes encircle the bog makes the area unsuitable for grazing (Figs 1b and S1), and its location at the upper limit of the treeline (Supplementary Fig. S2), means that purposeful deforestation for timber or pasturing has not been strong. Indeed, pollen data (Fig. 3 and Supplementary Figure 3), indicates there is no association between the observed increases in lithogenic element concentration, and decreased forest cover as inferred from arboreal taxa percentages. Further, no anthropic pollen taxa appear to increase in concentration during periods of high runoff (Supplementary Fig. S3), suggesting there is no relationship between greater human impact on the landscape, and greater erosion at this particular site. As such, we confidently infer that the periods with enriched lithogenic element concentrations are indicative of increased regional rainfall and runoff, with limited or no anthropic contribution except perhaps the last decades. There are five such major periods identifiable in the record (Fig. 2): 125–250 CE, 600–900 CE, 1050–1300 CE, 1400–1575 CE and 1725–1980 CE, hereafter referred to as Enhanced Runoff Periods (ERP).
The first ERP identified in our record for the last 2000 years has been dated between c. 125–250 CE (Fig. 2) which is in agreement with previous reconstructions indicating wetter conditions at this time45,46,47 in the Carpathian area, although other records suggest drier conditions23,48. After 600 CE, the Sureanu bog record indicates an extended period of minerogenic input onto the bog, persisting until roughly 900 CE, and representing the longest period of continual lithogenic input in the record (ERP 2). Therefore, we assume this is representative of increased rainfall and runoff at Sureanu throughout this time. The resolution of the Sureanu record allows for the identification of specific peaks in rainfall/runoff within this period, with 625 ± 64 CE and 825 ± 93 CE appearing to display particularly wet conditions. Higher frequency erosional events may be linked with the Dark Ages Cold Period (600–900 CE), characterised in south-eastern Europe by increased rainfall as documented in a wide range of high-altitude palaeohydrological reconstructions48,49. Furthermore, Carpathian fluvial archives document regular flooding during this time50,51 suggestive of increased regional runoff and moisture availability.
In the Sureanu record, two distinct phases are discernible within the MWP. The first, 900–1060 CE, is characterised by low values for all lithogenic elements (Figs 2, 3) and likely reflects a period of low precipitation. After 1050 CE, a rapid and clear increase in the lithogenic elements occurs up to 1300 CE, likely reflecting a rise in regional precipitation23,24,49 and denoted here as ERP 3 (Fig. 3a). High precipitation for this period is paralleled by an increase in Fagus25, which thrives in warm and wet conditions52, and a decrease in Alnus observed at this time in the Sureanu pollen record25. Such distinct two-phase hydrological variability during MWP has not been documented in other regional climate reconstructions, perhaps due to coarser temporal resolution of these records. Studies indicate dry conditions in western Romania53,54, but other proxy data provide contrasting reconstructions, with testate amoeba-derived records indicating wet conditions throughout much of the MWP49,55. This is in line with observations of increased erosional regimes23,35, reflecting a trend seen Europe wide18,56,57,58. Other than the shift observed at 1060 CE, the Sureanu record documents few variations (Fig. 3a), suggesting the MWP was likely moderately wet, but with relatively low variability as observed elsewhere31,54. Certainly when compared to the variability (both in terms of sharp fluctuations, and in amount of runoff and minerogenic input onto the bog) of later periods in the record (Fig. 3a), the MWP is not particularly distinct in the Sureanu archive.
The moisture availability that characterised the late MWP as reconstructed from the Sureanu record gradually decreased from 1250 CE onwards, reaching its nadir by 1300 CE, with dry conditions persisting until 1400 CE (Fig. 2). Limited runoff activity is reflected in the minerogenic matter record25, where limited lithogenic input has been documented at this time. This also explains the decrease in Fagus pollen at 1300–1400 CE at Sureanu25, as dry conditions are sub-optimal for this taxa52. This is reflective of region-wide dry conditions23, and relatively low hydroclimate variability, in contrast to western European records, where a number of studies document high rainfall in the period between MWP and the Little Ice Age (LIA)18,20,57,58.
After 1400 CE, a return to wetter conditions (ERP 4) may be inferred (Fig. 3), with the timing suggesting this shift is coeval with the onset of Little Ice Age (LIA). The LIA, a period of continent-wide cool temperatures, is characterised in the Sureanu record by distinct fluctuations in minerogenic input (Fig. 3b). Of these, two periods show the most significant variability at 1400–1575 CE and 1725–1980 CE, respectively. The resolution of our record allows for an unprecedented investigation into the extent of wet periods during and post LIA in this region. Our data indicate that such fluctuations in minerogenic input onto the bog could be linked to episodic decadal-length peaks in rainfall, separated by similarly long periods of low minerogenic input, and thus low precipitation. Variability in moisture availability during the LIA in eastern Europe has been observed previously49,59, but the available analytical resolution did not allow for observing the fine details of LIA hydroclimate variability as seen in the Sureanu record (Fig. 3b).
From 1400–1575 CE, three main peaks in minerogenic input are discernible, centred on 1465 ± 103, 1490 ± 91 and 1550 ± 85 CE (Fig. 3b). These data suggest the early LIA was characterised by higher and variable moisture availability, echoing results from other studies24,60, but that wet conditions were not characteristic of the whole early LIA. This is followed by evidence of a prolonged dry period in the Sureanu record, between 1550–1715 CE (Fig. 3b). This agrees well with other proxy reconstructions of hydroclimate variability, with very dry spells identified in eastern Carpathian peat records between 1550–1750 CE49,61, and after 1580 CE60.
This period of low moisture availability is followed by a sharp transition into the most variable section of the record (ERP 5), between 1725–1980 CE, encompassing the largest peaks in lithogenic elements (Fig. 3)25. As a result, these minerogenic-rich peaks may be considered representative of the wettest interval in the record, echoing the cold and wet summers recorded across Europe62, with the variability representative of high-frequency oscillations observed at other sites in central and eastern Europe61,63,64. This instability appears symptomatic of the spatial hydroclimate variability in the broader central-eastern Europe region (including the Carpathians) during the LIA. In northern Romania, wetter conditions were not observed until after 1780 CE, whilst in south-west Romania, regular droughts are seen in the 18th Century65. Persistent dry conditions are reconstructed for northern Romania after ~1550 CE based on testate amoeba24,60. The high precipitation period between 1780–1815 CE, appears linked to several regionally representative high river discharge events, with documentary evidence attesting to a number of major floods in the period 1784–1785 CE66, with limited model results corroborating wet conditions36 (Fig. 4).
Disagreements between central and eastern European palaeohydrological records have been observed before23, likely reflecting the north-south67 and east-west37 gradients in European hydroclimate. Particularly in the LIA, the Sureanu record displays notable correspondence with flooding and precipitation records in the Alpine region18, where evidence suggests high levels of rainfall, occurring in notable pulses17,22, as at Sureanu (Fig. 5). As such, it appears the Late Holocene Sureanu record largely reflects a similar variability in Atlantic derived rainfall, but modulated also by input from the Mediterranean as recently reconstructed for south-eastern Europe in a comparison between speleothem stable isotope data and regional pollen records68.
After 1900 CE there are two other periods of enhanced runoff in the Sureanu record (Fig. 2). The first spans the middle of the 20th Century (centred on 1935 ± 98 CE), whilst the second occurs during the century’s end, peaking by 1978 ± 65 CE (Fig. 2). The end of the first period may reflect a decrease in rainfall observed across Romania roughly 1969 CE28, reflected also in river discharge data (Fig. 4). However, the lack of similarity outside of this time to other regional drought indices65 and rainfall records may be indicative of erosion linked to increased local human activity, as ski slopes, and touristic activity became significant in the area25. What is clear from the most recent section of the Sureanu record is the increase in the prevalence of enhanced runoff in the last 500 years. Whether this variability reflects solely a higher incidence of extreme precipitation events or increasingly stronger human impact in the area requires further research.
To understand the forcing mechanisms behind the hydroclimatic variability reconstructed here, we compare the Sureanu record to a number of regionally representative climatic indices potentially linked to the variability observed (Fig. 5). It has been speculated by a number of studies that the North Atlantic Oscillation (NAO) should be the most important atmospheric pressure system influencing hydroclimatic variability in eastern Europe15,25,69. This is because during periods of NAO positive phase (NAO+), intensification of the westerlies shifts Atlantic storm tracks northwards. This means Atlantic derived moisture is distributed mainly in the north and west of the continent, whereas moisture over south-eastern Europe is derived primarily from the Mediterranean sea15,31. Conversely, when the NAO is in negative phase (NAO−), weak westerlies allow for the Atlantic storm tracks to move eastwards, resulting in increased Atlantic-derived winter precipitation over south-eastern Europe25,31. From the comparison in Fig. 5, it appears that the NAO variability controls some aspects of the Sureanu runoff record, with four periods of enhanced precipitation directly ascribable to NAO− conditions, and some short-term fluctuations, especially those within the LIA, also potentially linked to a negative NAO phase.
The primary low-frequency cycles observed in the Sureanu record (200 and 500-year, Fig. 6) do not appear to correlate with Atlantic climatic fluctuations, and so may be linked to another forcing mechanism. Further, the lack of any correlation with the Atlantic Multidecadal Oscillation, which has a cyclicity of 60–80 years70, suggests it plays a limited role in the hydroclimate system of eastern Europe as reconstructed from the Sureanu record. However, short-term oscillations as observed via the multi-taper method (MTM) could be linked to the NAO, with results displaying persistent 5-year and 4-year cycles in all elements analysed (Fig. 7). Similar cycle lengths have been observed previously in NAO-controlled varve sedimentation70, so it is possible that NAO variability is playing a distinct role in driving the small-scale, short-duration minerogenic deposition cycles identified in the Sureanu data.
Meteorological evidence suggests however that winter precipitation regimes over Romania are actually controlled by the combination of the NAO, Scandinavian Pattern (SCA), and the East Asia-West Russia (EA-WR) pattern33, whilst the impact of the changing strength of the Siberian High (SH) has also been suggested to further modulate the regional hydroclimate variability71. Interestingly, the onset of LIA in the Sureanu record (around 1400 CE), characterised by increased precipitation at this time (Fig. 3), occurs during a strengthening of the SH (Fig. 5). Disentangling the impact of the SH on most recent hydroclimate variability over south-eastern Europe, however, requires further research. In the case of the LIA and more recent minerogenic record at Sureanu, after apparent pacing of its onset by a strong SH, further peaks appear linked to NAO− conditions, whereas others (e.g., ~1740 and ~1950 CE) agree well with periods of stronger SH conditions (Fig. 5). This echoes previous work on flooding in the eastern Mediterranean, which displays a similar link between flooding and SH strength in the Late Holocene72. This evidence points to the fact that interplay between the two mechanisms is, at least periodically, and particularly in the past 500 years, modulating the signal observed at Sureanu. This observation clearly displays that a denser network of records is necessary for better disentangling the interplay between these circulation modes and their long-term environmental impact in the wider Eurasian region.
In addition to atmospheric pressure systems, periods of low solar irradiance, known as grand minima, have been linked to colder and wetter climates73,74,75. Noticeable similarities between the reconstructed runoff variability at Sureanu and solar activity may be observed (Fig. 5). The onset of the wet phase in the MWP occurred during the Oort minimum, whilst the first wet interval within the LIA occurred during the Spörer minimum, and later runoff events occur during the Dalton and Gleissburg minima (Fig. 4). This reflects well existing data from various studies in which wet conditions have been linked to solar minima, with such evidence across Europe63,76,77, and North America78. Further evidence of the impact of solar irradiance cycles on the Sureanu runoff record may be observed in the existence of 500-year and 200-year cycles in much of the lithogenic geochemical data (Fig. 5), and in the MTM power spectrum for Zr and Sr (Fig. 7). Such cycles may be linked to solar irradiance, with a 500-year cycle (a harmonic of the postulated Hallstatt cycle) observed in radiocarbon records79, and in flooding archives from central Europe18. The 200-year cycle observed in the last 1000 years of the Sureanu record (Fig. 6) is most likely representative of the Seuss cycle79,80. This oscillation has been noted previously in peat archives40,81, and in marine records from the Baltic Sea82.
Changes in solar radiation may influence trophospheric circulation through downward propagation of planetary waves83, particularly impacting atmospheric circulation pattern over the North Atlantic, and the influence it exerts on Eurasian hydroclimate variability. Therefore, the postulated influence of solar irradiance upon the atmospheric pressure systems (explaining why the two appear to fluctuate in tune) and their interaction appear to be the main driver behind past runoff variability observed in Sureanu record.
To further investigate the inferred past runoff record from Sureanu based on proxy-data, and the ability of climate models to address past regional climate fluctuations, we compare the Ti-derived record with modelled precipitation values for the past 2000 years, using the model output from PaleoView84, which utilises climate data from the TRaCE21ka experiment85,86. This experiment uses the simulated projections from Community Climate System Model version 3 (CCSM387,88), a coupled atmosphere-ocean general circulation model. This model resembles present-day regional and global temperature and hydrology data and may therefore be suitable for addressing climate variability over spatially large areas, and for long-term comparisons.
Figure 8 shows a comparison between data from three model scenarios (year-long, winter (DJF) and summer (JJA) averages) of the grid cell Sureanu is located in 45–47.5°N, 22.5–25°E, for central-eastern Europe (35–50 °N, 15–30 °E) and the Ti-derived record from Sureanu peatbog. Model data is presented in anomaly form, relative to the present day84. To allow for comparison, all data was brought onto the same 10-year interval timescale (using a 30-year Gaussian window) (Fig. 8). For all three scenarios, the Sureanu runoff data suggests a discrepancy in the amount of rainfall for the first 500 years (0–600 CE) of the compared sections when compared to model reconstructions (Fig. 8). Model data suggests this period is the wettest, whilst our reconstructed runoff data indicate low runoff, particularly between 250–600 CE. This discrepancy is clearest in the JJA and whole-year model output, whilst the DJF output for this time period is much closer to the signal observed in the XRF-CS data (Fig. 8). This lack of comparability for both JJA and whole-year model data continues throughout much of the record, although some increases in erosion observed toward the end of the record are reproduced by the model, with model-related precipitation increases observable between 1400–1600 CE and around 1800 CE. However, for most of the period these two model outputs bear little resemblance to the measured data. For the DJF model output, however, more similarities may be observed (Fig. 8). The three most recent periods of enhanced runoff at Sureanu parallel the periods of increased runoff, between 1100–1250 CE, 1350–1600 CE and 1700–1850 CE, with concurrent drops in modelled precipitation also present between these peaks. This apparent agreement indicates the Sureanu runoff record may be related primarily to winter precipitation, at least during this part of the record. Since winter precipitation over central-eastern Europe is generally linked to NAO fluctuations, this may explain why the Sureanu record, when comparing to proxy data reflecting atmospheric pressure systems, bears most relation to the NAO (Fig. 5). Considering the uncertainties of the model, for the last 500 years, we compare our dataset with annual precipitation reconstructed for the area 23°–24°E, 45°–46°N89, while for the very recent interval, spanning the past 150 years, we look at river discharge data (National Center for Atmospheric Research, 2001) (Fig. 4)90. For the wetter period centred on 1780 CE, it appears that Sureanu record agrees with reconstructed precipitation89 (Fig. 4c). In addition, the latest period of higher rainfall, at the end of the 20th Century appears to resemble peaks in discharge data.
Despite the apparent ability of the data to reflect regional hydrological trends, it is clear the model does not echo such data for much of the record, with periods of non-alignment, particularly in the earliest section of the record. Precipitation and other hydroclimatic indexes are complex, localized (dependent on topography and fine-scale processes) parameters36,91. Further, the coarse resolution of the model means the small-scale fluctuations recorded via river-discharge records, and the Sureanu record, are not reproduced. Therefore, such parameters are difficult to predict via models, and as a result the periods of good correspondence may be promising indicators of the improving nature of climate models. As such, further high-resolution records of palaeoprecipitation in the region are needed to constrain the modelling results. Approaches such as the one performed here allow for the presentation of high-resolution data and will allow for greater accuracy when developing model parameters in the future.
Conclusions
We present the first record of sub-decadal fluctuations in runoff activity in the Southern Carpathians, Romania for the last 2000 years using XRF-CS and ICP-OES derived lithogenic elemental data. Our data provide a high resolution record of palaeohydrological variability for this region, documenting five periods of enhanced runoff (125–250 CE, 600–900 CE, 1050–1300 CE, 1400–1575 CE and 1725–1980 CE). These periods generally agree well with published reconstructions but provide an unprecedented high-resolution insight into short-term fluctuations. In particular, the Medieval Warm Period was characterised by an early dry phase with limited runoff, followed by a wet, but stable phase after 1050 CE. The Little Ice Age is characterised by the most variable precipitation regime of the entire record, with the largest peaks in runoff proxies, indicative of the wettest conditions.
Our data suggests the NAO is the dominant control on Southern Carpathian hydroclimate, but that, as previously suggested, other controls, such as the Siberian High may play a role in moderating the NAO’s regional impact. Further, there appears to be a strong link between solar irradiance minima and the onset of wet periods, itself likely linked to the NAO, particularly during the past 1000 years. Data-model comparison appears to confirm the apparently strong impact of the NAO. The comparability of model-derived precipitation values for winter (DJF) to the erosion-related record suggest it primarily records input of lithogenic material following winter precipitation. The lack of similarity to summer and annual precipitation appears to confirm this. Despite an apparent link between model and observed data, particularly in the last 1000 years, there appears to be a disconnection for the period 0–1000 CE, where very few similarities may be observed. Such a conclusion indicates the necessity of further high-resolution studies into palaeohydrology in this region, in order to better constrain model predictions.
Methods
Site description
Sureanu peat bog (45°34.05′N, 23°30.28′E), is a small bog located at the foot of a glacial cirque (gradient in excess of 1 in 2) in the Southern Carpathians25. The bog is uncovered, hydrologically linked to the neighbouring glacial tarn Iezerul Sureanu (Figs 1 and S1, S2), and at the uppermost forest limit in this area92, at 1840 m above sea level. As such, the vegetation on the surrounding slopes is sparse, consisting of subalpine taxa, typically grasses, shrubs (especially Juniperus) and dwarf pine (Pinus mugo). Numerous mass-wasting related geomorphological features may be observed (Fig. 1c), including avalanche channels. Climatically, the area is temperate continental, with typical average winter temperatures between −2 °C and −7 °C, and between 8–19 °C in summer93. Typically, precipitation amounts to between 900–1800 mm per year, mainly as snowfall during winter months. The site, therefore, is under snow cover for between 100–200 days per year93. Because of the bog’s location at the foot of a glacial cirque, any debris dislodged from surrounding slopes during periods of slope erosion is deposited onto its surface, thus disturbing peat growth. The reconstruction of such deposition therefore may be used as a proxy for changes in erosion and linked to rainfall variability25.
Two separate cores were extracted from the bog centre using a Russian peat corer25. The first core (SUR-1, 603 cm long) was wrapped in clingfilm, and transported to Northumbria University, where it was stored at 3 °C prior to subsampling for Inductively Coupled Plasma-Optical Emission Spectrometry (ICP-OES) analysis. The second core (SUR-2) was transported to the University of Cologne for X-Ray Fluorescence-Core Scanning (XRF-CS) analysis. As both cores were located next to one another, the correlation between them was established via alignment of geochemical data.
Geochemical methods
XRF-CS was performed on the uppermost 200 cm of SUR-2 record, using an ITRAX core scanner equipped with an Si drift detector94 at the University of Cologne (Institute of Geology and Mineralogy), Germany. Measurements were made at 2 mm intervals, with 20 s count times at each interval. Analysis was performed using a Cr X-Ray tube set to 30 kV and 50 mA. From the wide range of elements analysed, only lithogenic elements (Fe, Rb, Sr, Ti, Zr), which display high enough concentrations above the detection limits were selected for further interpretation. For ease of interpretation, a 9-point running average has been applied to the raw XRF-CS data. To ensure various issues associated with XRF-CS analysis of peat, including high but variable organic content, variable density, and changing water content are taken into account, the raw cps values have been normalised with respect to total (incoherent + coherent) scattering40,95.
To further test whether the ITRAX data reliably reproduced the geochemical signal, ICP-OES analysis of the same lithogenic elements was carried out on samples collected at roughly 2 cm intervals. These samples were brought into solution using a mixed acid (HNO3-HCl-HF) digestion25, prior to analysis via a Perkin Elmer Optima 8000 system at Northumbria University. An internal standard (either 1ppm In or 1ppm Sc) was measured to correct for instrumental drift, whilst two certified reference materials (Montana Soil 2711 and IAEA Lake Sediment) were run alongside samples, with recoveries documented in Supplementary Table S1. Blanks, with negligible contamination for all elements were also run at regular intervals. We directly compare our data with the loss on ignition-derived record of minerogenic deposition25.
Chronology
The age model is based upon 9 radiocarbon dates, analysed on bulk peat samples from the uppermost 200 cm of core SUR-125. Analysis was carried out via accelerator mass spectrometry at HEKAL AMS Laboratory, MTA ATOMKI Institute for Nuclear Research of the Hungarian Academy of Sciences in Debrecen, and at the 14C Centre at Queen’s University Belfast. Using the IntCal13 Calibration curve96, and Bacon software, an age-depth model was developed, indicating 200 cm of sediment documents just over 2000 years of deposition (Supplementary Fig. S4).
Time series analysis
To identify non-stationary cyclicities in the proxy data, continuous Morlet wavelet transform was used97,98. The Morlet wavelet is commonly employed in analysing geophysical data and produces a two-dimensional frequency-time mapping of the data, allowing identification of quasiperiodic oscillations in time series with multiple frequencies. To allow for this analysis, the normalised lithogenic element data were interpolated to equal 1-year time steps using a Gaussian window of 3 years. To further analyse the data, the Multi-Tapered Method (MTM) was performed99,100. The MTM does not have an a priori assumption about the processes producing the time series that are analysed, and is a useful technique for analysing time series which could have both continuous and singular components101. For ease of interpretation, only datasets displaying significant cyclicities are displayed here.
Data Availability
All ITRAX-derived data is uploaded as supplementary information. Other datasets generated during and/or analysed during the current study would be provided on request.
References
Kundzewicz, Z. W. et al. Flood risk and climate change: global and regional perspectives. Hydrol. Sci. J. 59, 1–28 (2014).
IPCC. The Physical Science Basis. Working Group I Contribution to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge, United Kingdom and New York, USA (2013).
Min, S. K., Zhang, X., Zwiers, F. W. & Hegerl, G. C. Human contribution to more-intense precipitation extremes. Nature 470, 378–381 (2011).
Fischer, E. M. & Knutti, R. Observed heavy precipitation increase confirms theory and early models. Nat. Clim. Chang. 6, 986–991 (2016).
Micu, D. M., Dumitrescu, A., Cheval, S. & Bîrsan, M.-V. Climate of the Romanian Carpathians: Variability and Trends, https://doi.org/10.1007/978-3-319-02886-6 (2015).
Zappa, G., Shaffrey, L. C., Hodges, K. I., Sansom, P. G. & Stephenson, D. B. A multimodel assessment of future projections of north atlantic and european extratropical cyclones in the CMIP5 climate models. J. Clim. 26, 5846–5862 (2013).
Woollings, T., Gregory, J. M., Pinto, J. G., Reyers, M. & Brayshaw, D. J. Response of the North Atlantic storm track to climate change shaped by ocean-atmosphere coupling. Nat. Geosci. 5, 313–317 (2012).
Mizuta, R. Intensification of extratropical cyclones associated with the polar jet change in the CMIP5 global warming projections. Geophys. Res. Lett. 39, 1–6 (2012).
Ulbrich, U., Leckebusch, G. C. & Pinto, J. G. Extra-tropical cyclones in the present and future climate: A review. Theor. Appl. Climatol. 96, 117–131 (2009).
Zhang, F. & Weng, Y. Predicting hurricane intensity and associated hazards: A five-year real-time forecast experiment with assimilation of airborne doppler radar observations. Bull. Am. Meteorol. Soc. 96, 25–33 (2015).
Mihăilescu, C. Clima şi hazardurile Moldovei-evoluţia, starea, predicţia. (Licorn, 2004).
Teodoreanu, E. Little climate optimum in the Carpathian-Danubian-Pontic space. Present Environ. Sustain. Dev. 11 (2017).
Topor, N. Ani ploioşi şi secetoşi în Republica Populară Romînă. (Institutul Meteorlogic, 1964).
Bojariu, R. & Paliu, D.-M. North Atlantic Oscillation Projection on Romanian Climate Fluctuations in the Cold Season in Detecting and Modelling Regional Climate Change and Associated Impacts (eds India, M. B. & Bonillo, D. L.) 345–356 (Springer Berlin Heidelberg, 2001).
Bojariu, R. & Giorgi, F. The North Atlantic Oscillation signal in a regional climate simulation for the European region. Tellus 57A, 641–653 (2005).
Swierczynski, T. et al. A 1600 yr seasonally resolved record of decadal-scale flood variability from the Austrian Pre-Alps. Geology 40, 1047–1050 (2012).
Swierczynski, T. et al. Mid- to late Holocene flood frequency changes in the northeastern Alps as recorded in varved sediments of Lake Mondsee (Upper Austria). Quat. Sci. Rev. 80, 78–90 (2013).
Wirth, S. B., Glur, L., Gilli, A. & Anselmetti, F. S. Holocene flood frequency across the Central Alps – solar forcing and evidence for variations in North Atlantic atmospheric circulation. Quat. Sci. Rev. 80, 112–128 (2013).
Bajard, M. et al. Erosion record in Lake La Thuile sediments (Prealps, France): Evidence of montane landscape dynamics throughout the Holocene. The Holocene, https://doi.org/10.1177/0959683615609750 (2015).
Arnaud, F. et al. Erosion under climate and human pressures: An alpine lake sediment perspective. Quaternary Science Reviews 152, 1–18 (2016).
Pierre, S. et al. 6-kyr record of flood frequency and intensity in the western Mediterranean Alps – Interplay of solar and temperature forcing. Quat. Sci. Rev. 170, 121–135 (2017).
Wilhelm, B., Vogel, H., Crouzet, C., Etienne, D. & Anselmetti, F. S. Frequency and intensity of palaeofloods at the interface of Atlantic and Mediterranean climate domains. Clim. Past 12, 299–316 (2016).
Haliuc, A. et al. Palaeohydrological changes over mid and late Holocene in the Carpathian area, central-eastern Europe. Glob. Planet. Change 1–43, https://doi.org/10.1016/j.gloplacha.2017.02.010 (2017).
Diaconu, A. C. et al. How warm? How wet? Hydroclimate reconstruction of the past 7500 years in northern Carpathians. Romania. Palaeogeogr. Palaeoclimatol. Palaeoecol. 482, 1–12 (2017).
Longman, J., Ersek, V., Veres, D. & Salzmann, U. Detrital events and hydroclimate variability in the Romanian Carpathians during the Mid-to-Late Holocene. Quat. Sci. Rev. 167, 78–95 (2017).
Warken, S. F. et al. Reconstruction of late Holocene autumn/winter precipitation variability in SW Romania from a high-resolution speleothem trace element record. Earth Planet. Sci. Lett. 499, 122–133 (2018).
Kern, Z. et al. Natural proxy records of temperature- and hydroclimate variability with annual resolution from the Northern Balkan–Carpathian region for the past millennium – Review & recalibration. Quat. Int. 415, 109–125 (2016).
Tomozeiu, R., Stefan, S. & Busuioc, A. Winter precipitation variability and large-scale circulation patterns in Romania. Theor. Appl. Climatol. 81, 193–201 (2005).
Rimbu, N., Stefan, S. & Necula, C. The variability of winter high temperature extremes in Romania and its relationship with large-scale atmospheric circulation. Theor. Appl. Climatol. 121, 121–130 (2015).
Obreht, I. et al. Tracing the influence of Mediterranean climate on Southeastern Europe during the past 350,000 years. Sci. Rep. 6, 36334 (2016).
Perşoiu, A. et al. Holocene winter climate variability in Central and Eastern Europe. Sci. Rep. 7 (2017).
Panagiotopoulos, F. et al. Observed Trends and Teleconnections of the Siberian High: A Recently Declining Center of Action. J. Clim. 18, 1411–1422 (2005).
Ionita, M., Chelcea, S., Rimbu, N. & Adler, M. J. Spatial and temporal variability of winter streamflow over Romania and its relationship to large-scale atmospheric circulation. J. Hydrol. 519, 1339–1349 (2014).
Morellón, M. et al. Human–climate interactions in the central Mediterranean region during the last millennia: The laminated record of Lake Butrint (Albania). Quat. Sci. Rev. 136, 134–152 (2016).
Magyari, E. et al. Palaeolimnology of the last crater lake in the Eastern Carpathian Mountains: a multiproxy study of Holocene hydrological changes. Hydrobiologia 631 (2009).
Xoplaki, E. et al. Modelling Climate and Societal Resilience in the Eastern Mediterranean in the Last Millennium. Hum. Ecol. 1–17 (2018).
Roberts, N. et al. Palaeolimnological evidence for an east-west climate see-saw in the Mediterranean since AD 900. Glob. Planet. Change 84–85, 23–34 (2012).
Longman, J., Veres, D. & Wennrich, V. Utilisation of XRF core scanning on peat and other highly organic sediments. Quat. Int. In Press, (2018).
Urdea, P. & Reuther, A. U. Some new data concerning the Quaternary glaciation in the Romanian Carpathians. Geogr. Pannonica 13, 41–52 (2009).
Longman, J. et al. Periodic input of dust over the Eastern Carpathians during the Holocene linked with Saharan desertification and human impact. Clim. Past 13, 897–917 (2017).
Allan, M. et al. Mid- and late Holocene dust deposition in western Europe: The Misten peat bog (Hautes Fagnes – Belgium). Clim. Past 9, 2285–2298 (2013).
Kylander, M. E. et al. Potentials and problems of building detailed dust records using peat archives: An example from Store Mosse (the “Great Bog”), Sweden. Geochim. Cosmochim. Acta 190, 156–174 (2016).
Stuut, J. B. & Prins, M. A. The significance of particle size of long-range transported mineral dust. PAGES Mag. 22, 70–71 (2014).
Finsinger, W. et al. Holocene fire-regime changes near the treeline in the Retezat Mts. (Southern Carpathians, Romania). Quaternary International (2016).
Feurdean, A. et al. 12,000-Years of fire regime drivers in the lowlands of transylvania (Central-Eastern Europe): A data-model approach. Quat. Sci. Rev. 81, 48–61 (2013).
Diaconu, A.-C., Grindean, R., Panait, A. & Tanţău, I. Late Holocene palaeohydrological changes in a Sphagnum peat bog from NW Romania based on testate amoebae. Stud. UBB Geol. 60 (2016).
Hubay, K. et al. Holocene environmental changes as recorded in the geochemistry of glacial lake sediments from Retezat Mountains, South Carpathians. Quat. Int, https://doi.org/10.1016/J.QUAINT.2018.02.024 (2018).
Cristea, G. et al. Carbon isotope composition as indicator for climatic changes during the middle and late Holocene in a peat bog from Maramures Mountains (Romania). The Holocene 24, 15–23 (2013).
Feurdean, A. et al. Last Millennium hydro-climate variability in Central-Eastern Europe (Northern Carpathians. Romania). The Holocene 25, 1179–1192 (2015).
Chiriloaei, F., Rǎdoane, M., Perşoiu, I. & Popa, I. Late Holocene history of the Moldova River Valley. Romania. Catena 93, 64–77 (2012).
Rădoane, M. et al. History of Holocene fluvial activity in Romania: evidences based on absolute dating. Proc. 33rd Rom. Geomorphol. Symp. 92–96 (2017).
Feurdean, A., Tanţău, I. & Fărcaş, S. Holocene variability in the range distribution and abundance of Pinus, Picea abies, and Quercus in Romania; implications for their current status. Quat. Sci. Rev. 30, 3060–3075 (2011).
Forray, F. L. et al. A Late Holocene environmental history of a bat guano deposit from Romania: an isotopic, pollen and microcharcoal study. Quat. Sci. Rev. 127, 141–154 (2015).
Onac, B. P. et al. A 2500-yr late Holocene multi-proxy record of vegetation and hydrologic changes from a cave guano-clay sequence in SW. Romania. Quat. Res. 83, 437–448 (2015).
Panait, A. et al. Hydrological conditions and carbon accumulation rates reconstructed from a mountain raised bog in the Carpathians: A multi-proxy approach. CATENA 152, 57–68 (2017).
Magny, M. Holocene climate variability as reflected by mid-European lake-level fluctuations and its probable impact on prehistoric human settlements. Quat. Int. 113, 65–79 (2004).
Wilhelm, B. et al. 1400years of extreme precipitation patterns over the Mediterranean French Alps and possible forcing mechanisms. Quat. Res. 78, 1–12 (2012).
Arnaud, F. et al. Lake Bourget regional erosion patterns reconstruction reveals Holocene NW European Alps soil evolution and paleohydrology. Quat. Sci. Rev. 51, 81–92 (2012).
Gałka, M. & Apolinarska, K. Climate change, vegetation development, and lake level fluctuations in Lake Purwin (NE Poland) during the last 8600 cal. BP based on a high-resolution plant macrofossil record and stable isotope data (δ13C and δ18O). Quat. Int. 328, 213–225 (2014).
Schnitchen, C. et al. Reconstructing hydrological variability from testate amoebae analysis in Carpathian peatlands. J. Paleolimnol. 36, 1–17 (2006).
Feurdean, A., Klotz, S., Mosbrugger, V. & Wohlfarth, B. Pollen-based quantitative reconstructions of Holocene climate variability in NW Romania. Palaeogeogr. Palaeoclimatol. Palaeoecol. 260, 494–504 (2008).
Charman, D. J., Blundell, A., Chiverrell, R. C., Hendon, D. & Langdon, P. G. Compilation of non-annually resolved Holocene proxy climate records: stacked Holocene peatland palaeo-water table reconstructions from northern Britain. Quat. Sci. Rev. 25, 336–350 (2006).
Gałka, M. et al. Disentangling the drivers for the development of a Baltic bog during the Little Ice Age in northern Poland. Quat. Int. 328–329, 323–337 (2014).
Marcisz, K. et al. Long-term hydrological dynamics and fire history over the last 2000 years in CE Europe reconstructed from a high-resolution peat archive. Quat. Sci. Rev. 112, 138–152 (2015).
Levanič, T., Popa, I., Poljanšek, S. & Nechita, C. A 323-year long reconstruction of drought for SW Romania based on black pine (Pinus Nigra) tree-ring widths. Int. J. Biometeorol. 57, 703–714 (2013).
Brázdil, R. et al. European floods during the winter 1783/1784: Scenarios of an extreme event during the ‘Little Ice Age’. Theor. Appl. Climatol. 100, 163–189 (2010).
Magny, M. et al. Geoscientific Instrumentation Methods and Data Systems North–south palaeohydrological contrasts in the central Mediterranean during the Holocene: tentative synthesis and working hypotheses. Clim. Past 9, 2043–2071 (2013).
Drăguşin, V. et al. Constraining Holocene hydrological changes in the Carpathian–Balkan region using speleothem δ18O and pollen-based temperature reconstructions. Clim. Past 10, 1363–1380 (2014).
Krichak, S. O. & Alpert, P. Signatures of the NAO in the atmospheric circulation during wet winter months over the Mediterranean region. Theor. Appl. Climatol. 82, 27–39 (2005).
Ólafsdóttir, K. B., Geirsdóttir, Á., Miller, G. H. & Larsen, D. J. Evolution of NAO and AMO strength and cyclicity derived from a 3-ka varve-thickness record from Iceland. Quat. Sci. Rev. 69, 142–154 (2013).
Obreht, I. et al. Shift of large-scale atmospheric systems over Europe during late MIS 3 and implications for Modern Human dispersal. Sci. Rep. 7, 5848 (2017).
Benito, G., Macklin, M. G., Zielhofer, C., Jones, A. F. & Machado, M. J. Holocene flooding and climate change in the Mediterranean. Catena 130, 13–33 (2015).
Steinhilber, F. et al. 9,400 Years of Cosmic Radiation and Solar Activity From Ice Cores and Tree Rings. Proc. Natl. Acad. Sci. 109, 5967–5971 (2012).
Shindell, D. T., Schmidt, G. A., Mann, M. E., Rind, D. & Waple, A. Solar forcing of regional climate change during the Maunder Minimum. Science 294, 2149–2152 (2001).
Czymzik, M., Muscheler, R. & Brauer, A. Solar modulation of flood frequency in central Europe during spring and summer on interannual to multi-centennial timescales. Clim. Past 12, 799–805 (2016).
Benito, G. et al. Recurring flood distribution patterns related to short-term Holocene climatic variability. Sci. Rep. 5, 1–8 (2015).
Martin-Puertas, C. et al. Regional atmospheric circulation shifts induced by a grand solar minimum. Nat. Geosci. 5, 397–401 (2012).
Nichols, J. E. & Huang, Y. Hydroclimate of the northeastern United States is highly sensitive to solar forcing. Geophys. Res. Lett. 39 (2012).
Damon, P. E. & Jirikowic, J. L. The Sun as a Low-Frequency Harmonic Oscillator. Radiocarbon 34, 199–205 (1992).
Lüdecke, H.-J., Weiss, C. O. & Hempelmann, A. Paleoclimate forcing by the solar De Vries/Suess cycle. Clim. Past Discuss. 11, 279–305 (2015).
Swindles, G. T., Patterson, R. T., Roe, H. M. & Galloway, J. M. Evaluating periodicities in peat-based climate proxy records. Quat. Sci. Rev. 41, 94–103 (2012).
Yu, S.-Y. Centennial-scale cycles in middle Holocene sea level along the southeastern Swedish Baltic coast. Geol. Soc. Am. Bull. 115, 1404 (2003).
Woollings, T., Lockwood, M., Masato, G., Bell, C. & Gray, L. Enhanced signature of solar variability in Eurasian winter climate. Geophys. Res. Lett. 37, 1–6 (2010).
Fordham, D. A. et al. PaleoView: a tool for generating continuous climate projections spanning the last 21 000 years at regional and global scales. Ecography (Cop.). 40, 1348–1358 (2017).
Liu, Z. et al. Evolution and forcing mechanisms of El Niño over the past 21,000 years. Nature 515, 550–553 (2014).
Otto-Bliesner, B. L. et al. Coherent changes of southeastern equatorial and northern African rainfall during the last deglaciation. Science 346, 1223–7 (2014).
Collins, W. D. et al. The Community Climate System Model Version 3 (CCSM3). J. Clim. 19, 2122–2143 (2006).
Yeager, S. G. et al. The Low-Resolution CCSM3. J. Clim. 19, 2545–2566 (2006).
International Hydrological Programme, Bodo, B. & Global Runoff Data Center. Augmented Monthly Flow Rates of World Rivers (except former Soviet Union). NCAR Research Data Archive, 10.5065/D61G0JGZ (2001).
Pauling, A., Luterbacher, J., Casty, C. & Wanner, H. Five hundred years of gridded high-resolution precipitation reconstructions over Europe and the connection to large-scale circulation. Clim. Dyn. 26, 387–405 (2006).
Smerdon, J. E. et al. Comparing proxy and model estimates of hydroclimate variability and change over the Common Era. Clim. Past 13, 1851–1900 (2017).
Cristea, V. Fitosociologie şi vegetaţia României. (Babes-Bolyai University Press, 1993).
Fărcaş, I. & Sorocovschi, V. The climate of the Retezat Mountains in The Retezat National Park, Ecological Studies (ed. Popovici, I.) 13–20 (West Side Computers, 1992).
Croudace, I. W., Rindby, A. & Rothwell, R. G. ITRAX: description and evaluation of a new multi-function X-ray core scanner in New Techniques in Sediment Core Analysis (ed. Rothwell, R. G.) 51–63, 10.1144/GSL.SP.2006.267.01.04 (Geological Society, London, Special Publications, 2006).
Kylander, M. E. et al. A novel geochemical approach to paleorecords of dust deposition and effective humidity: 8500 years of peat accumulation at Store Mosse (the ‘Great Bog’), Sweden. Quat. Sci. Rev. 69, 69–82 (2013).
Reimer, P. et al. IntCal13 and Marine13 Radiocarbon Age Calibration Curves 0–50,000 Years cal BP. Radiocarbon 55, 1869–1887 (2013).
Grinsted, A., Moore, J. C. & Jevrejeva, S. Application of the cross wavelet transform and wavelet coherence to geophysical time series. Nonlinear Process. Geophys. 11, 561–566 (2004).
Torrence, C., Compo, G. P., Torrence, C. & Compo, G. P. A Practical Guide to Wavelet. Analysis. Bull. Am. Meteorol. Soc. 79, 61–78 (1998).
Percival, D. B. & Walden, A. T. Spectral Analysis for Physical Applications. (Cambridge University Press, 1993).
Thomson, D. J. Spectrum estimation and harmonic analysis. Proc. IEEE 70, 1055–1096 (1982).
Ghil, M. et al. Advanced spectral methods for climatic time series. Rev. Geophys. 40, 1003 (2002).
Mayewski, P. A. & Maasch, K. A. Recent warming inconsistent with natural association between temperature and atmospheric circulation over the last 2000 years. Clim. Past Discuss. 2, 327–355 (2006).
Olsen, J., Anderson, N. J. & Knudsen, M. F. Variability of the North Atlantic Oscillation over the past 5,200 years. Nat. Geosci. 5, 1–14 (2012).
Trouet, V. et al. Persistent Positive North Atlantic Oscillation Mode Dominated the Medieval Climate Anomaly. Science. 324, 78–80 (2009).
Acknowledgements
We would like to thank Northumbria University for JL’s postgraduate studentship, and the QRA for their “QRA-14CHRONO Centre Radiocarbon Dating” award which provided funding for some of the radiocarbon analyses. We are indebted to F. Gogaltan and C. Popa for help during fieldwork. This is a contribution to the project PN-II-ID-PCE-2012-4-0530 “Millennial-scale geochemical records of anthropogenic impact and natural climate change in the Romanian Carpathians”.
Author information
Authors and Affiliations
Contributions
J.L., D.V. and V.E. designed the study. J.L., D.V., V.E. and A.H. wrote the manuscript. J.L. completed the ICP-OES analysis and V.W. the ITRAX analysis. All authors reviewed the manuscript.
Corresponding author
Ethics declarations
Competing Interests
The authors declare no competing interests.
Additional information
Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Longman, J., Veres, D., Ersek, V. et al. Runoff events and related rainfall variability in the Southern Carpathians during the last 2000 years. Sci Rep 9, 5334 (2019). https://doi.org/10.1038/s41598-019-41855-1
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-019-41855-1
This article is cited by
-
A millennium-long climate history of erosive storms across the Tiber River Basin, Italy, from 725 to 2019 CE
Scientific Reports (2021)
-
Historical predictability of rainfall erosivity: a reconstruction for monitoring extremes over Northern Italy (1500–2019)
npj Climate and Atmospheric Science (2020)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.