Abstract
The Central European area has been extensively studied using qualitative reconstruction techniques focusing on the development of peatlands in the past; however, research based on quantitative techniques is still lacking, especially in relation to highlands and mountainous areas. In this study, we focused on the hydrological changes and human-induced disturbances that affected two raised bogs from the Orawa-Nowy Targ Basin (Carpathian region): Puścizna Krauszowska and Puścizna Mała. We aimed to reconstruct the development of peatlands and changes in water table under different intensities of human activities. Two peat sequences from two different bogs, both possessing absolute chronologies based on high-resolution 14C dating, were examined in terms of pollen, non-pollen palynomorphs, plant macrofossils and testate amoeba records. We detected an asynchronous decrease in the water table level on the bogs that took place between the 4th and the 7th century
Keywords
Introduction
Ombrotrophic peatlands dominated by Sphagnum mosses occurring in the temperate and boreal zone of Northern Hemisphere are important carbon sinks (Yu et al., 2009) and unique habitats for relic plants of the Last Glaciation (e.g. Bastl et al., 2009; Koczur, 2006; Marcisz et al., 2015). Peat growth in the ombrotrophic bogs is dependent on precipitation, and the vegetation is isolated from mineral-rich groundwater. Furthermore, peatlands are very sensitive to climatic changes such as shifts in the water balance (the relationship between precipitation and evapotranspiration; Wieder et al., 2006) and thus perceived as highly suitable for palaeoclimatic reconstructions (e.g. Chambers et al., 2012; Gałka et al., 2016; Lamentowicz et al., 2008; Langdon et al., 2003; Nichols et al., 2009; Väliranta et al., 2012). Anoxic conditions and low pH provide suitable conditions for the preservation of pollen, plant macrofossils, testate amoebae (TA), non-pollen palynomorphs (NPPs) and other subfossils (e.g. Chambers et al., 2012; Moore et al., 1991; Rydin and Jeglum, 2006), making them invaluable archives of environmental changes.
The last two millennia were a period of substantial climatic and demographic changes in Europe. Apart from the recent climate-related warming (Hanhijärvi et al., 2013; Mann et al., 2003), many other significant climatic events occurred during this period, such as (1) 1.6 ka event (Bond et al., 1997, 2001) which occurred simultaneously with dramatic demographic changes referred to as the Migration Period (Büntgen et al., 2011); (2) the Mediaeval Climate Anomaly (MCA), also known as the ‘Mediaeval Warm Period’ (MWP;
In this study, we focused on the Orawa-Nowy Targ Basin where the biggest raised bog complex (8255.6 ha) in Poland is located. Like many other peatland areas in Europe, the study area was strongly transformed by human activities (Łajczak, 2006). Previous palaeoecological studies conducted in the Orawa-Nowy Targ Basin based on analyses of pollen, plant macrofossils, a very limited set of identified NPPs, dendrochronology and TA traits (e.g. Kołaczek et al., 2010; Krąpiec et al., 2016; Marcisz et al., 2016; Obidowicz, 1990, 1996). These works focused on general patterns in the development of peatlands and vegetation history. However, there are no detailed in-depth palaeoecological research that are supported by quantitative reconstructions of bog’s hydrology. Herein, we present a multi-proxy study on changes in the peatland ecosystem based on pollen, a broad set of NPPs, plant macrofossils and TA-based quantitative reconstruction of water table depth during the last 2000 years. We aimed to reconstruct the development of peatlands and changes in water table, paying close attention to climatic and human impact changes.
Study area
The study area is located in the Orawa-Nowy Targ Basin (covering approximately 262 km2 mainly at an altitude of 490–680 m a.s.l.), an area bounded by the Western Beskid Mountains (Babia Góra and Gorce ranges) in the north and the Tatra Mountains in the south (Figure 1). The Orawa-Nowy Targ Basin is a tectonic depression filled with Neogene and Quaternary deposits (Watycha, 1976). The limited permeability of the loam layer in this basin’s bed is a key factor that supports peatland development in this area (Łajczak, 2009; Kowanetz, 1998). Another factor which supports high-ground wetness is climate. Moderately cool climate in the basin with a tendency to stagnate cold air during the winter supports maintenance of snow cover for long periods (Kondracki, 1998; Lorenc, 2005). In addition, between May and October precipitation exceeds evaporation (Kowanetz, 1998; Obrębska-Starklowa, 1977). The heavy cloudiness and frequent fog further contribute to the condensation of water in the basin (Staszkiewicz and Szeląg, 2003). The basin’s mean annual temperature is 5–6°C; the coldest month being January (mean values range between −3.5°C and −4.0°C; Lorenc, 2005) and the warmest being July (mean temperatures reaching 15.6°C; Kowanetz, 1998). The basin’s annual precipitation is 900–1100 mm, which is highest in July (Lorenc, 2005).

Location of Puścizna Mała and Puścizna Krauszowska sites: (a) across Central, Western and Northern Europe, dots and open circles with numbers reflect location of sites referred in the ‘Discussion’ section; (b) detailed relief map of the Orawa-Nowy Targ Basin and adjacent geographical units, black rectangles represent area shown in c and d; (c) aerial view on the PM bog (GoogleEarth, date of image acquisition 19 March 2015); and (d) aerial view on the Puścizna Krauszowska bog (GoogleEarth, date of image acquisition 19 March 2015).
As many as 420 species of vascular plants have been found in the Orawa-Nowy Targ peatlands, but only 57 of them were found growing on typical raised bog domes. Bog coniferous forests occupy vast areas adjoining the peat bogs. The area is/was a refugium of glacial relics such as Pinus mugo, Rubus chamaemorus (Koczur, 2006), and until the early 19th century
Two peatlands – Puścizna Krauszowska and Puścizna Mała – were selected for this study, given their well-documented history of human-induced disturbances (Łajczak, 2006). The peat from both bogs is still being extracted; hence, the bog surfaces are continuously shrinking. According to Łajczak (2006), the contemporary Puścizna Krauszowska bog covered 79 ha, including 29 ha of the fragmentary dome, but the present surface of the dome does not exceed 5 ha (B. Fiałkiewicz-Kozieł, personal communication, 2008). The Puścizna Mała bog extends to cover an area of more than 125 ha with a dome surface area covering more than 51 ha. Historical sources revealed that the intensity of human activities varied over time. During the early Slavonic and the early Mediaeval periods, the boggy, forested and almost uninhabited area of the Orawa-Nowy Targ Basin was a natural barrier against the war raids from the south, which means human-induced disturbances were rare. From the late Mediaeval period, the increasing number of settlements brought about gradual deforestation, and the local people were encouraged to exploit the bog in different ways for subsistence, such as pastoral purposes and/or peat extraction. The peat extraction was not significant until the 20th century
Methods
Coring and absolute chronology of the profiles
For the purpose of our study, the top sections of the peat deposits were collected. The Puścizna Mała profile was retrieved as peat monolith (135 × 20 × 20 cm) from an outcrop in the bog’s dome (613 m a.s.l., 49°28′06″N, 19°56′18″E) in June 2006 (Fiałkiewicz-Kozieł et al., 2014; Kołaczek et al., 2010). The Puścizna Krauszowska profile was collected from the centre of the dome (655 m a.s.l., 49°27′36″N, 19°47′12″E) in June 2008, using a Wardenaar sampler of dimensions measuring 100 × 10 × 10 cm3 (Fiałkiewicz-Kozieł et al., 2015; Kołaczek et al., 2010).
An absolute chronology is based on (1) the age–depth models based on numerous 14C dates (11 dates for Puścizna Krauszowska and 21 dates for Puścizna Mała) and constructed using the OxCal v. 4.2 software applying the P_Sequence function (Bronk Ramsey, 1995, 2008) and the IntCal13 atmospheric curve (Reimer et al., 2013) as the calibration set and (2) 210Pb dates obtained for the topmost sections of the profiles (Fiałkiewicz-Kozieł et al., 2014, 2015; Figure 2). In the case of Puścizna Mała model, in comparison with the one published by Fiałkiewicz-Kozieł et al. (2014), two boundaries (introduced by the command Boundary) that reflect potential changes in the rate of peat accumulation were introduced at depths of (1) 114 cm – an abrupt decrease in pollen concentration and a fall in the content of amorphous organic matter and (2) at 103 cm – an increase in Sphagnum magellanicum (see the following sections) and an increase in bulk density (Fiałkiewicz-Kozieł et al., 2014). In the case of Puścizna Krauszowska profile, the age–depth model no. 3 has been selected from scenarios proposed in the paper by Fiałkiewicz-Kozieł et al. (2015; Figure 2). To facilitate a more detailed description and discussions about results, the age is presented as a µ (mean) value of the modelled age expressed by

Age–depth models from the Puścizna Krauszowska and Puścizna Mała profiles.
The 14C dating of both profiles revealed inversions in the layers between depths of 48.5 and 21 cm in Puścizna Krauszowska and 45 and 21.5 cm in Puścizna Mała, which hampered the calculation of absolute chronology (Fiałkiewicz-Kozieł et al., 2014, 2015). These fragments of the profiles were excluded from this study.
Testate amoebae (TA)
An analysis of TA was carried out to obtain quantitative estimates of the depth to water table (DWT) changes (Charman et al., 2007). Subfossil TA were extracted from 3–5 cm3 subsamples following the procedure described by Booth et al. (2010). The material was washed through a sieve of 250-µm mesh size, and the filtrate was used for microscopic analyses. Samples were counted under an upright microscope at 200× and 400× magnifications until a minimum of 150 tests was detected (Booth et al., 2010; Payne and Mitchell, 2008). Identification criteria were drawn from available literature, since we tried to achieve the highest possible taxonomical resolution (e.g. Clarke, 2003; Mazei and Tsyganov, 2006; Meisterfeld, 2001; Ogden and Hedley, 1980). DWT was estimated from TA using a training set consisting of 123 samples developed for northern Poland by Lamentowicz and Mitchell (2005) and Lamentowicz et al. (2008), which is the only training set available from Poland.
Pollen and non-pollen palynomorphs (NPPs)
Palynological analysis was used for the reconstruction of (1) terrestrial vegetation changes, (2) human impact and (3) qualitative water table fluctuations (e.g. Cugny et al., 2010; Gaillard, 2013; Moore et al., 1991; Poska et al., 2004). Preliminary results considering vegetation changes based on palynological analysis carried out on 60 samples, with very limited data considering NPPs, were published in the paper by Kołaczek et al. (2010). Consequently, in this paper, the focus is on the palaeohydrology and bog ecosystem dynamics. In this study, the resolution of pollen analysis chosen for the Puścizna Mała profile increased twice compared with the one used by Kołaczek et al. (2010), whereas with Puścizna Krauszowska profile, the number of samples from the bottom section increased. Samples (1 cm3 in volume) were prepared using standard laboratory procedures: (1) a 10% HCl solution was used to dissolve potential carbonates in the samples, (2) samples were heated in a 10% KOH to remove humic compounds and (3) samples were subjected to a minimum 24-h treatment with HF to remove mineral fractions. Next, acetolysis was carried out (Berglund and Ralska-Jasiewiczowa, 1986). One Lycopodium tablet of known number of spores (N = 10,679; produced by the Lund University) was added to each sample for calculating microfossil concentrations (Stockmarr, 1971). Using an upright microscope with 400× and 1000× magnifications, pollen and spores were counted until the number of arboreal pollen (AP) grains in a sample reached at least 500. In extreme cases where low concentrations of pollen material were present, pollen grains were counted to get a minimum of 250 AP grains. Pollen grains and spores were identified with atlases and keys (Beug, 2004; Moore et al., 1991) and the reference pollen slide collection of the W. Szafer Institute of Botany, Polish Academy of Sciences (Kraków, Poland). NPPs were identified using photographs and descriptions reported in various publications (e.g. van Geel, 1978; van Geel and Aptroot, 2006, as well as those cited in Miola, 2012). The NPP-type nomenclature follows Miola (2012). Percentages of pollen grains, spores and NPPs were calculated according to the following formula: number of taxon’s sporomorphs/TPS × 100%; here TPS (total pollen sum) includes both AP and NAP (non-AP, excluding aquatic and wetland plants such as Ericaceae and Cyperaceae, spores and NPPs).
Plant macrofossils
In the Puścizna Krauszowska profile, plant macrofossils were analysed from 1.5-cm-thick (depth: 0–15 cm) and 3-cm-thick (depth: 15–100 cm) slices of the profile, whereas in the Puścizna Mała profile, 2-cm-thick (depth: 0–28 cm) and 1-cm-thick (depth: 28–135 cm) slices were analysed. In both cases, slices (ca. 3–5 cm3 in volume) were deemed contiguous if the material was not utilized by other analyses. Each sample was washed and wet-sieved through 0.2-mm mesh screens. The carpological and vegetative remains (leaves, rootlets and epidermis) were identified using available keys and atlases (e.g. Grosse-Brauckmann, 1974; Mauquoy and van Geel, 2007). The volume percentages of the various unidentified vegetative remains and Sphagnum taxa were estimated to the nearest 5% with a stereomicroscope. The relative proportions of the taxonomic sections of Sphagnum were estimated on the basis of the branch leaves, using an upright microscope with two 22 × 22-mm cover glasses. Then, a species-level identification of the Sphagnum was separately performed using appropriate guides for identification (Hölzer, 2010; Laine et al., 2011). The moss nomenclature follows Ochyra et al. (2003), whereas the vascular plant nomenclature follows Mirek et al. (2002). The volume proportion of the amorphous organic matter, which serves as a measure of peat decomposition, was estimated during sieving in steps of 25% and was calculated on the basis of the ratio of the volume of plant remains on the screen to the volume of the sample prior to sieving (same procedure used in Gałka et al., 2013). In the case of Puścizna Mała profile, identification of fossil S. capillifolium and S. rubellum as well as S. angustifolium and S. fallax was found to be impossible due to the similar morphology of branches and stem leaves. Nonetheless, interpretation of DWT and acidity was not a serious problem because these species occur in similar habitats (Hölzer, 2010; Laine et al., 2011).
Presentation of results
The results of each analysis were summarized in diagrams drawn using the TILIA Graph programme (Grimm, 1991). The reconstructions of DWT were plotted in the C2 software (Juggins, 2003). Due to different sampling resolution in particular analyses, diagrams were divided into two common phases of bogs’ development on the basis of water table reconstructions and the level of human impact.
Results and interpretation: palaeoecological reconstruction
Phase I: Moderately wet conditions and low human impact (Puścizna Krauszowska: 100–48.5 cm, ca. 70 bc –ad 680; Puścizna Mała: 135–45 cm, ca. ad 120–1360)
The TA-inferred water table slightly fluctuated in both profiles but was the highest in the period spanned by both profiles. The mean DWT in the Puścizna Krauszowska profile was 10.7 cm, and in Puścizna Mała, it was 6.9 cm. About

Testate amoebae (TA) and TA-inferred depth to water table (DWT) from the Puścizna Krauszowska and Puścizna Mała profiles. The grey area of each column shows percentage values magnified by 10.
In both profiles, layers with very low TA concentrations were recorded (Puścizna Krauszowska: ca.

Plant macrofossils from the Puścizna Krauszowska and Puścizna Mała profiles. The grey area of each column shows percentages magnified by 10.

Non-pollen palynomorphs (NPPs) and local pollen taxa from the Puścizna Krauszowska and Puścizna Mała profiles. The grey area of each column shows percentages magnified by 10. sapr.: saprotrophic.
Pollen analysis revealed the domination of mixed forests in the region (AP: 86–99%). The forests consisted of beech (Fagus sylvatica: 10–20% in Puścizna Krauszowska, 10–24% in Puścizna Mała), fir (Abies alba: 7–22% in Puścizna Krauszowska, 7–26% in Puścizna Mała) and spruce (Picea abies: 6–26% in Puścizna Krauszowska, 5–23% in Puścizna Mała), with an addition of hornbeam (Carpinus betulus: 3–9% in Puścizna Krauszowska, 1–10% in Puścizna Mała) and oak (Quercus: 3–10% in Puścizna Krauszowska, 3–10% in Puścizna Mała). The relative proximity of the edge of mixed forests to the coring spot is highlighted by the presence of spores of parasitic fungus Kretzschmaria deusta, which is characterized by the poor dispersal ability (cf. van Geel and Andersen, 1988). The drier parts of peatlands were occupied by pines (P. sylvestris type: 7–27% in Puścizna Krauszowska, 11–45% in Puścizna Mała) and birches (Betula: 6–14% in Puścizna Krauszowska, 7–20% in Puścizna Mała), whereas lags and riverine woodlands were occupied by alder (Alnus: 8–16% in Puścizna Krauszowska, 3–18% in Puścizna Mała). Alnus visibly retreated in the 5th (Puścizna Mała) and 7th (Puścizna Krauszowska) centuries
Phase II: Dry conditions and strong human impact (Puścizna Krauszowska: 21–0 cm, ad 1900–2008; Puścizna Mała: 21.5–0 cm, ad 1840–2006)
The mean DWT was 26.3 cm for the Puścizna Krauszowska profile and 32.1 cm for the Puścizna Mała profile. In the Puścizna Krauszowska profile, DWT increased ca.
The Puścizna Krauszowska profile revealed very dynamic changes in Sphagnum communities as subsequent dominations of S. cuspidatum, S. rubellum, S. balticum, S. capillifolium and finally again S. rubellum. This suggests a sequence of an initial transition from a wet bog and/or the presence of a hollow towards dry bog conditions and/or the formation of a hummock. Afterwards, together with the occurrence of S. balticum, the microhabitat characterized by high water table reappeared, finally followed by dry conditions on the mire and/or hummock development (cf. Hölzer, 2010). The formation of a hummock and/or decrease in the water table on the bog provided a suitable habitat for C. vulgaris, whose pollen percentage maximum (40%) coincided with the domination of S. rubellum. These possible dry and wet shifts are poorly manifested by the TA composition, which revealed dry conditions during the entire phase; however, S. rubellum increases in frequency are parallel with decreases in the water table. In the Puścizna Mała profile, S. magellanicum and S. capillifolium/S. rubellum dominated, but mosses such as Polytrichum strictum and Aulacomnium palustre became more frequent. This set of taxa points to a continuation of relatively dry conditions (Bragazza and Gerdol, 1996; Manukjanová et al., 2014). Percentages of C. vulgaris reached the highest values in the Puścizna Mała profile and rapidly decreased afterwards (ca.
Considerable landscape openness in the bog’s surroundings is reflected by low AP values (Puścizna Krauszowska: 50–87%, Puścizna Mała: 61–81%). However, in the Puścizna Krauszowska profile, AP revealed an increasing trend caused mainly by a rise in P. sylvestris type percentages (from 18% to 37%), which might have been an effect of pine spread on the mire surface and/or preference for this taxon in silviculture. In the vicinity of the bogs, cultivated fields (pollen: Cerealia type, Triticum type and Avena type), meadows and ruderal communities (pollen: Poaceae undiff., Plantago lanceolata, Artemisia, Chenopodiaceae and others) occupied considerable areas (Figure 6). The probable expansion of habitats rich in nitrogen was recorded by maxima of Urtica (cf. Ellenberg and Leuschner, 2010).

Pollen of forest and open land taxa included to total pollen sum from the Puścizna Krauszowska and Puścizna Mała profiles. The grey area of each column shows percentages magnified by 10. The values of HdV-44 Kretzschmaria deusta are magnified by 10. L.c.TA: low concentrate of testate amoebae.
Discussion
Migration and Mediaeval Periods in the Orawa-Nowy Targ peatlands across climate and palaeohydrological patterns in Central European peatlands
The Migration and further Mediaeval Periods were characterized by high climate variability in Europe (Büntgen et al., 2011). One of the most climate-affecting events during this period was a discharge of the drifting ice to the Atlantic Ocean ca.

Testate amoebae (TA)–inferred depth to water table (DWT) from the Puścizna Krauszowska and Puścizna Mała profiles on the background climatic and hydro-climatic reconstruction from Central and Western Europe during the period
The other palaeohydrological reconstructions from the period
Another period of distinct climate changes, recorded in several palaeoecological archives, was the MWP (
The last event of sun activity during the Mediaeval Period being potentially recorded in studied bogs is the Wolf solar minimum (
Bog’s functioning after the drainage
The last two centuries brought the most severe transformation of the Orawa-Nowy Targ peatlands. The peat extraction dramatically intensified, which was connected to the development of a drainage network (Łajczak, 2006). This is manifested in both studied profiles by very low water table levels. During the second half of the 19th and the first half of the 20th century
Conclusions
The study of environmental change that took place in the last two millennia in the Orawa-Nowy Targ Basin (Western Carpathians) based on the Puścizna Krauszowska and Puścizna Mała bogs revealed different palaeohydrological events recorded by biotic proxies:
Two events of drought between the 4th and the 7th century
An increase in the water table reflected by the Puścizna Mała bog (ca.
A water table decrease was a possible reaction to summer and autumn drought(s) in the period of
Footnotes
Acknowledgements
Barbara Fiałkiewicz-Kozieł is acknowledged for giving us permission and full access to the peat cores designated for analyses presented in this contribution. We would like to thank Michał Słowiński and Sandra Słowińska for the preparation of samples for the analysis of testate amoebae, and B. Nowaczyńska for carrying out the laboratory preparation of samples designated for palynological analysis. We are grateful to Bas van Geel and two anonymous reviewers for their linguistic support and critical comments on the previous versions of this manuscript.
Funding
The present research was supported by the National Science Centre (Poland) – grants 2011/01/D/ST10/02579 (PI: Barbara Fiałkiewicz-Kozieł) and UMO-2014/13/B/ST10/02091 (PI: Monika Karpińska-Kołaczek); by the Swiss Contribution to the enlarged European Union – grant PSPB-013/2010 CLIMPEAT; and statutory funds of the Institute of Geoecology and Geoinformation (Faculty of Geographical and Geological Sciences, Adam Mickiewicz University) – research topics WNGiG/2012/S/P-B/017 and subsidy for young researchers –0800000000/505/060/STAT_MN (funds for Piotr Kołaczek, 2015).
