BOG PINE AND DECIDUOUS TREES CHRONOLOGIES RELATED TO PEAT SEQUENCES STRATIGRAPHY OF THE PODEMSZCZYZNA PEATLAND (SANDOMIERZ BASIN, SOUTHEASTERN POLAND)

ABSTRACT The Podemszczyzna peatland (Sandomierz Basin, SE Poland) is a place of peat exploitation for balneological purposes. The thickness of organic sediments (minerogenic peat) reaches 4.0 m, while the beginning of peat accumulation was dated using the radiocarbon method (14C) at 13,517–13,156 cal BP. During the peat exploitation numerous fragments of subfossil wood (of various species) were excavated and, based on dendrochronological analyzes and 14C dating (wiggle-matching), two short floating chronologies were elaborated: bog pine chronology (147 years long) and deciduous trees (oak, elm) chronology (139 years long). 14C dating has shown that the bog pine chronology (ca. 9980–9830 mod. cal BP) is the oldest pine chronology found in the Polish peatlands so far. It was synchronous with the Preboreal decline of fluvial activity and peat formation, whereas dying off of trees was connected with distinct rise of fluvial activity. Floating chronology of deciduous trees is much younger and encompasses time interval of ca. 680–545 cal BP. The trees’ encroachment on the peatland was related to the terrestrialization of the depositional fen, recorded in the loss on ignition curve in the form of mineral sediment delivery to the bog, as well as it is marked in the pollen record.

In the Polish territory dendrochronological studies of bog pine and bog oak trees have been so far conducted for the Rucianka peatland (Polish Lowland) ( Figure 1A-site 12) (Barniak et al. 2014). Bog pine chronology was also studied in the Imszar raised bog (Polish Lowland) ( Figure 1A-site 13) (Margielewski et al. 2022a) and in the Puścizna Wielka raised bog (Polish Carpathians, southern Poland-see Figure 1A-site 14) . The newest site, where subfossil trunks were excavated during the peat exploitation, is the Podemszczyzna peatland situated in Southeastern Poland ( Figure 1A-site 15). In this case, bog pine and deciduous trees dendrochronology, along with peat analysis using pollen, non-pollen palynomorphs, Cladocera, geochemical and lithological analyses, provided the basis for the examination of the influence of climate on the phases of growth and dying off of trees in the peatland area.

THE STUDY AREA
Podemszczyzna peatland is located in the village of Podemszczyzna, in the SE part of the Tarnogród Plateau, which is part of the Sandomierz Basin, in close proximity to the Roztocze Upland ( Figure 1A) (Kondracki 2001). A fen, with an area of 25 ha, was formed in a depression developed within the Świdnica river valley which contemporary flows in the vicinity ( Figure 1B). Nowadays, peat is being exploited here (with the use of peat excavator) for balneological purposes of the nearby Horyniec health resort ( Figure 1C). Two artificial ponds were formed in post-exploitation depressions (Figure 1 C-D). During the peat extraction, several dozen pieces of wood were excavated, including numerous subfossil tree trunks which were used for dendrochronological analysis (Figure 2).

Coring and Sampling
Coring was performed in the central parts of the peatland, in the area between two postexploitation ponds ( Figure 1B) (GPS coordinates 50º12.768'N, 23º17.861'E). Drilling was performed using an Instorf sampler, 10 cm in diameter. Two complete logs of peat were taken for multiproxy analysis.
Samples for dendrochronological analysis were collected from subfossil tree trunks extracted during the peat exploitation, using a gas-engine chainsaw ( Figure 2C and D). In total, 40 samples of wood (wood slices) of various species were collected.

Stratigraphy of the Podemszczyzna Peatland 1559
Dendrochronological Analysis Tree species were identified anatomically (see Schweingruber 1988Schweingruber , 1990Godet 2008). Annual tree rings were precisely measured (with the accuracy of 0.01 mm), using the Dendrolab 1.0 equipment (Zielski and Krąpiec 2004). Processing of the measured tree-ring sequences was carried out using the TREE-RINGS software (Krawczyk and Krąpiec 1995), and the TSAP computer programme (Rinn 2005).
The crossdates of samples were statistically validated, using t-coefficient (Baillie and Pilcher 1973), and "Gl coefficient" (Gleichlaeufigkeit; Eckstein and Bauch 1969). The accuracy of the measurements and the quality of the cross-dating were verified by the COFECHA software (Holmes 1999).
The relative age of wood samples, enabling the development of floating chronologies, have been determined by the wiggle-matching method with a use of the OxCal computer program, v. 4.3 (Bronk Ramsey 2009) (Figure 2A and B).

Radiocarbon Dating
Conventional radiocarbon dating of organic material (peat, wood) using the LSC (liquid scintillation counting) method was conducted in the Laboratory of Absolute Dating in Kraków, Poland, according to standard procedure (Skripkin and Kovalyukh 1998). AMS dating was made in the Laboratory of Absolute Dating in Kraków, Poland (graphite preparation) and in the University of Georgia, USA (measurements) (Cherkinsky et al. 2010).
Ten 14 C dates of organic deposits and seven dates of wood (enabling wiggle-matching analysis) were obtained (Table 1). Calibrated radiocarbon dates (denoted by "cal BP" abbreviation) with 95.4% probability, were determined using the OxCal computer program v. 4.3 (Bronk Ramsey 2009 and the IntCal20 radiocarbon calibration dataset (Reimer et al. 2020).
The age-depth model was elaborated using the OxCal computer program (v. 4.4.3), and their P_sequence function with the following parameters: k 0 = 0.8, log 10 (k/k 0 ) = 1 and interpolation = 0.5 cm (Bronk Ramsey 2009), based on IntCal 20 dataset (Reimer et al. 2020). Sections of log, where the changes in accumulation rate are expected, were marked using Boundary command (Bronk Ramsey 2009). As a result, a modeled age (mod. cal BP) for each 0.5 cm sequence of deposits was obtained.

Stratigraphy of the Podemszczyzna Peatland 1561
Pollen and Non-Pollen Palynomorphs (NPPs) Analyses Samples of 1 cm 3 in volume were prepared using the standard procedure (Berglund and Ralska-Jasiewiczowa 1986;Moore et al. 1991). Pollen analysis of 60 samples (taken for each 5-10 cm) was performed ( Figure 3B). Pollen analysis was carried out with an Olympus BX43 light microscope with a magnification of 600×. More than 1000 terrestrial pollen grains were counted in each sample. NPPs were analyzed along with the pollen ( Figure 1C), according to van Geel (1978van Geel ( , 2001, Bell (1983), Jankovská and Komárek (2000) procedures.
Calculations and presentation of pollen and NPPs data were performed using TILIA Graph software (Grimm 1991). The pollen diagram of terrestrial taxa ( Figure 3B) was divided into local pollen assemblage zones (LPAZ) and local pollen assemblage subzones (LPASZ) ( Figure 1B). Aquatic pollen and non-pollen palynomorphs are shown in a separate diagram ( Figure 3C)  Samples (of 1 cm 3 in volume) taken for each 10 cm of the profile, were prepared according to procedure proposed by Frey (1986Frey ( , 1987. The extracted Cladocera and Rotifer remains, were stored in 10 mL of water with glycerine and safranine. Taxa were identified and counted at 200 or 400× magnification under a Nikon 50i microscope. All skeletal parts were counted: headshields, shells, postabdomens and postabdominal claws. The results of analyses are presented in diagram elaborated using TILIA Graph computer program (Grimm 1991) ( Figure 3D).

Geochemistry
Geochemical analysis was performed for 45 samples taken for each 10 cm of the profile. Samples were prepared according to the procedure proposed by Pansu and Gautheyrou (2006). The content of macroelements (Ca, Mg, Na, K, Fe) and trace elements (Mn, Cu, Zn, Ni, Pb) were measured using Atomic Absorption Spectrometry (AAS) with a Thermo Scientific iCE 3500 apparatus. On the geochemical diagram developed in TILIA Graph computer program (Grimm 1991), were presented both: the variability of the content of the main elements in the profile (Mg, Na, Mn, Cu, Pb), and geochemical environment indicators: Na/K, Ca/Mg, NaKMg/Ca, Fe/Mn, Fe/Ca Cu/Zn ( Figure 3E). Based on the contents of 11 geochemical elements, the Principal Component Analysis (PCA) was carried out using Statistica 13 and its results were presented on the diagram ( Figure 3E).

Sedimentary Sequence
In the bottom parts of the analyzed sediment profile, in the interval of 4.5-4.0 m, sandy silt occurs. The accumulation of organic sediments begins at a depth of approx. 4 m. It is represented by strongly decomposed peat (decomposition degree R> 65%), without mineral deposit admixture. The sludge is characterized by high losses on ignition (90%) ( Figure 3A). Higher in the sediment profile, in the interval of 2.5-1.6 m, there is a sedge peat (containing Carex sp.), with admixture of common reed (Phragmites australis) macroremnants. In the interval of 1.6-0.9 m, strongly decomposed sedge peat (R> 65%) appears again, and above it (in the interval of 0.9-0.2 m) sedge-reed peat (Magnocaricioni) was accumulated. In the upper parts of the peat complex (the last 20 cm), occurring peat is again heavily decomposed (R> 65%) and contains a large admixture of mineral sediment. The deposit is characterized by a high ash content (only 25% loss on ignition) ( Figure 3A).

Dendrochronology
Among 40 samples taken during field works, 31 samples of subfossil wood were used for dendrochronological analysis (8 samples of pine trees, and 23 samples of deciduous trees). The other samples could not be analyzed due to the poor condition of the wood or too few annual rings. The average age of analyzed trees was 83 years for pines, and 87 for deciduous trees, whereas the age of individual trees ranged from 30-50 to a maximum of 137 years (Figure 2A, B).
Cross-analysis (crossdates) of all samples, analyzed separately for Pinus and deciduous trees, allowed to identify onest showing the greatest similarity. Based on the best correlated dendrochrochonological sequences, the mean dendrochronological curve of Pinus sylvestris lasting 147 years was determined ( Figure 2B). The chronology is confirmed by high mean values of t-coefficient (6.6) and Gl being 70.5% (P<0.01). The second, younger chronology was compiled on the basis of deciduous trees samples (Oak, Elm) and is 139 years long (t=8.3; Gl=77%) ( Figure 2A).
Relative age of the both floating chronologies was obtained using the wiggle-matching technique. The tree rings selected for the radiocarbon dating and wiggle-matching analyses and results are presented in Figure 2A and B, and in Table 1. Based on wiggle-matching analysis, the time range of the floating pine chronology was determined at ca. 9980-9830 mod. cal BP ( Figure 2B). In turn, the time span of the younger floating chronology of deciduous trees is ca. 680-545 mod. cal BP (Figure 2A).

Local Vegetation History
Local vegetation history was compiled for individual chronozones based on the characteristics of the Local Pollen Assemblages Zones (LPAZ) and Local Pollen Assemblages Subzones (LPASZ) ( Figure 3B). Absolute age of pollen zones limits was determined on the basis of the modeled age (mod. cal BP).
Older Dryas Stadial (Pinus LPAZ) (before 13,386 mod. cal BP) During this period, the dominance of pine forests and the local, very small share of alder (Alnus sp.) is visible. The share of herbaceous plants is marginal, limited mainly to grasses (Poaceae), and sedges (Cyperaceae) and mugwort (Artemisia). There is also some pollen of aquatic plants (Nuphar).
Allerød Interstadial (Pinus-Betula LPAZ) (13,386-11,584 mod. cal BP). The Allerød Interstadial here is clearly bipartite ( Figure 3B). In the older part (AL1) (Pinus-Betula LPASZ) the stand is dominated by Pinus sylvestris and Betula, sometimes with an admixture of Picea and thermophilous species: Quercus, Carpinus and Corylus. Salix and Alnus occur locally. Artemisia and Chenopodiaceae have a large share among the herbaceous vegetation. The presence of grass pollen (Poaceae) and sedges (Cyperaceae) is also noted. In the younger part of Allerød Interstadial (AL2) (LPASZ: Pinus-Betula-NAP), pine (Pinus), and birch (Betula) are still dominant, but the share of other woody species, including those with higher thermal requirements (Corylus, Ulmus, Quercus, Carpinus) is significant. Picea, Fagus and (traces of) Abies alba also appear. A similar share of pollen of thermophilous plants in the Late Glacial deposits, accompanied by pollen of Fagus and Abies, was stated by K. Mamakowa (1962) in the Allerød sedimentary sequences of the Sandomierz Basin peatlands (Podbukowina, Świlcza, Obary peatlands), linking their occurrence with a presence of on-site refugia, rather than with reworking of pollen grains from older sediments or with a long transport (Mamakowa 1962). The share of Poaceae and Cyperaceae communities is increasing, and there are also Artemisia and Chenopodiaceae.
Younger Dryas Stadial. Pinus LPASZ (11,584-11,061 mod. cal BP). The forest communities were dominated by Pinus (which is characterized by expansion here) with an admixture of Picea. The cooling of the climate was clearly visible in the form of a drastic decrease in the share of most trees, including Betula, Alnus, and especially thermophilous species (Corylus, Ulmus, Quercus).

Stratigraphy of the Podemszczyzna Peatland 1565
Preboreal Phase. Pinus LPASZ (11,061-9448 mod. cal BP). In general, during the Preboreal Phase, two episodes in the initial and final period of this chronozone are clearly visible, related to the reduction of the pine's range and the short-term expansion of Betula, Alnus and thermophilous Corylus, Ulmus and Quercus. These two "warm" episodes are separated by a clear, sharp decline of a percentage of the above-described taxa, probably related to the cooling of the Preboreal Oscillation (Björck et al. 1997). In the youngest part of the Preboreal Phase, Tilia and Fagus appear, and there is also a greater proportion of ferns (Filicales monolete and Thelypteris palustris) ( Figure 3C). Herbaceous plants are relatively poorly represented in this period, mainly by Poaceae and Cyperaceae.
Boreal Phase. Pinus-Betula-Corylus LPASZ (9448-6871 mod. cal BP). The decrease in the proportion of pine pollen is accompanied by an increase (sometimes rapid) in the proportion of Betula, Alnus, Corylus, Ulmus, Fraxinus, Tilia, Quercus, as well as a small share of Carpinus and Fagus-pollen of these last two taxa was already found in Boreal sedimentary sequence of the Sandomierz Basin peatlands (Mamakowa 1962 Characteristic for the palynological profile of the peat bog are frequent fluctuations in the proportion of pine pollen in the diagram with the general dominance of pollen of this species. The local decrease in the share of pine is always accompanied by an increase in the share of the Betula pollen, as well as pollen of the thermophilous taxa: Corylus, Ulmus, Alnus, and Quercus ( Figure 3B).

Rotifera and Cladocera
In the Podemszczyzna sedimentary sequence, presence of 7 taxa belonging to two systematic groups has been confirmed: rotifers (2 taxa) and cladocerans (5 taxa) ( Figure 3D). In samples containing Cladocera it was, however, impossible to identify them to a lower systematic degree. Their chitinous remains were very damaged and very thin.
Among rotifers taxa two species were recognized: first belonging to Bdelloidea: Habrotrocha angusticollis and second belonging to Monogononta: Keratella cochlearis. Shells of H. angusticollis were found in significant densities in the sediment samples corresponding to the end part of the Preboreal Phase (PB) ( Figure 3D). The highest densities of these species were identified at the depths of 160 cm-41 ind. cm 3 ; 170 cm-35 ind. cm 3 , 210 cm-15 ind. cm 3 , and 220 cm-22 ind. cm 3 . K. cochlerais was found in a trace amount in the samples of Atlantic and Subboreal Phases ( Figure 3D).
In Cladocera group all identified taxa belonged to the Chydoridae family, which is a littoral taxa group. Due to a small density of Cladocera taxa, which ranges between 1 to 4 ind. cm 3 , the dominant taxa could not be distinguished. Cladocera remains were found in the Atlantic Phase (170-100 cm interval), in the Boreal Phase (190 cm), some single pieces in the Preboreal Phase (220 cm) and in the Allerød Interstadial (380-390 cm) ( Figure 3D).
On the basis of taxonomic analysis and environmental preferences, two groups of organisms were distinguished: the first one related to the acidic environment and wet area with low level of water (H. angusticollis, Alonella excisa, and Chydorus sphaericus) and the second one related to the open space of stagnant water Keratella cochlearis ( Figure 3D).

Geochemical Analysis
In the sedimentary profile of the Podemszczyzna peatland eight geochemical zones (GZI-GZVIII) were distinguished that differ in content of geochemical elements and indicators of the sedimentation environment ( Figure 3E).
GZ I (4. 13,440 mod. cal BP). First geochemical zone is related to the deposition of sandy clayey silt underlying the peat sequence. The horizon is marked by slightly increasing content of organic matter (from 3.5 to 18.0%) and occurrence of almost all lithophilic and trace elements (except Cu). In relation to the whole profile, values of Fe/ Mn and NaKMg/Ca ratios reach here their maximum ( geochemistry variation is directly proportional to the content of organic matter (0.76) and inversely proportional to concentrations of almost all major and trace elements (except Ca and Mg). Less important second component (PCA 2) which explain ca. 14.3% of the total geochemical variance is related mainly to content of Ca (0.85) ( Figure 3E).

DISCUSSION
Radiocarbon dating and palynological analysis indicate that the sedimentary basin was formed in the Older Dryas Stadial of the Late Glacial ( Figure 3A-B). The presence of pollen of aquatic plants (Nuphar, Sparganium) indicates that at that time in the sedimentary basin there was a permanent water reservoir and sedimentation was of a minerogenic character ( Figure 3B). This is confirmed by geochemical indicators suggesting also significant erosion in the vicinity of the reservoir at that time ( Figure 2E). The surroundings of the reservoir were dominated by pine forest with an admixture of Betula and Alnus ( Figure 3B). The peat accumulation (strongly decomposed sedge peat) began during the Allerød Interstadial climate warming ( Figure 3A-B). Delivery of mineral sediment to the peatland is visible as the decline of the loss on ignition curve in the middle Allerød depositional sequence ( Figure 3A) and it was caused by flooding of the river waters over the fen, as geochemical indicators suggest (phase GZIII-see Figure 3E). This, in turn, may be related to the cooling and growth of climate humidity by about 13,200-12,800 cal BP known as the Gerzensee Oscillation . In the Podemszczyzna profile, the sediment record of the Allerød Interstadial is bipartite, but the usual vegetation pattern-birch dominating in the older Allerød phase and pine-in the younger phase (see Latałowa 2003)-is in this case different. Two birch-pine phases are separated here by a phase with a predominance of pine ( Figure 3B). In the sedimentation sequence of Allerød there is a significant share of thermophilous plants (Corylus, Ulmus, Quercus, Carpinus), which usually accompany a simultaneous decrease in the share of pine ( Figure 3B). The significant share of pollen of thermophilous plants in the Allerød sedimentary sequence, which is also accompanied by Fagus and Abies, occurs within the aforementioned peatlands of the Sandomierz Basin (Mamakowa 1962). Also in the Polish Carpathians, in the Late Glacial sedimentary sequences of peatlands (Bølling Interstadial), there is sometimes a significant share of pollen of thermophilous plant taxa, invariably accompanied by Fagus and Abies (Margielewski et al. 2022b). In these cases, it would be difficult to explain the presence of such abundant share of pollen in deposits with their distant transport or redeposition. These sitesx probably constituted refugia during the Late Glacial, what was suggested by K. Mamakowa (1962) (see also Margielewski et al. 2022b) Climate warming in the Preboreal Phase was marked by a decrease in pine pollen concentration and an increase in a share of Betula, Alnus and thermophilous Corylus, Ulmus, Quercus ( Figure 3B). In the middle part of the Preboreal Phase, however, another regression of thermophilous species is visible, the content of Picea pollen is also increasing ( Figure 3B). This regression of termophilous plants dated at 10,744 ± 118-10,159 ± 114 mod. cal BP, is probably related to the cooling of the climate during the Preboreal Cold Oscillation, recorded at that time in lake sediments from Poland (Fiłoc et al. 2016). Geochemical indicators suggest higher water level in the fen at that time ( Figure 3E). In the younger phase of Preboreal, a short-term regression of pine is again noted and the sudden development of Salix, Betula, Alnus along with thermophilous Corylus, Ulmus Tilia and Quercus occurs ( Figure 3B). During this period, pine trees encroachment into the fen is dated at 9980-9830 mod. cal BP ( Figure 2B and 3B). The pine woodlands in the vicinity of Stratigraphy of the Podemszczyzna Peatland 1569 the fen thinned out, which is indicated by the significant share of Corylus and ferns (Filicales monolete, Thelypteris palustris) ( Figure 3B and C). The encroachment of pines into the fen (its surface was about 2.2 m below the modern ground level) coincided with a strong climate overdrying in Central Europe (ca. 10,200-9600 cal BP), causing, among other phenomena, a decline of fluvial activity of the European rivers ( Figure 4) (Starkel et al. 2013). This period was, however, not stable: the reduction of annual growth of pine indicates a local deterioration of ecological conditions ca. 9955-9945 and 9860-9845 mod. cal BP: trees dying off registered in the peat bog was associated with the latter period ( Figure 2B).
Later pine trees dying off recorded in the fen was caused by the deterioration of climatic and ecological conditions: the pollen curves show a short-term regression of thermophilous species ( Figure 3B). There is a sharp increase in the share of Cladocera and Rotifera in the sediments, Figure 4 Correlation of subfossil bog pine and bog deciduous trees (oak, elm) dendrochronology of the Podemszczyzna fen, in relation to bog pine and bog oak time range during the Holocene in other European peatlands (location of analyzed sites and authors of dendrochronological analysis-see Figure 1A), and in relation to palaeoclimatic changes in the Northern Hemisphere (Bond et al. 2001(Bond et al. , 2008Starkel et al. 2013;Wanner et al. 2015). Holocene chronostratigraphy after Mangerud et al. (1974), Starkel et al. (2013), and Walker et al. (2019).
previously unheard of, which indicates an increase in the water level in the bog ( Figure 3D). In the fen sediments, just above the horizon with pine trunks, a delivery of mineral sediment to the reservoir is visible, marked as decline on the loss on ignition curve, and related to a short-term flooding of the reservoir with water ( Figure 3A). The geochemical indices that show an increase in chemical denudation are also changing rapidly ( Figure 3E). The cooling of the climate at that time coincided with this moistening (Figure 4). In the period after 9600 cal BP, a distinct rise of fluvial activity, landslide intensification and debris flow in the European areas were recorded (Mayewski et al. 2004;Pánek et al. 2013;Starkel et al. 2013;Margielewski 2018).
Bog pine floating chronology from the Podemszczyzna fen is currently the oldest bog pine chronology in Polish and European territory (see Barniak et al. 2014;Margielewski et al. 2022a) (Figure 4). It should be mentioned that previously some short pine chronologies for the Late Glacial (Bølling-Younger Dryas) have been elaborated in the pre-Alps, German and Polish territories, but they are, however, not related to bog trees (Schaub et al. 2008;Kaiser et al. 2012;Krąpiec et al. 2020).
Since the middle of the Atlantic Phase, pollen profile shows a significant decrease in the proportion of pine: the Betula and Alnus curves and the thermophilous species are already continuous here (Corylus, Ulmus, Tilia, Quercus), which indicates a permanent reconstruction of plant communities around the reservoir ( Figure 3B). At the onset of the Subatlantic Phase, a mineral cover on the peat begins to form due to the delivery of mineral sediment to the peatland ( Figure 3A). There is a level of trunks of deciduous trees within it, suggesting that the fen area was settled at ca. 680-545 mod. cal BP (ca. 1270-1305 mod. cal AD) by Quercus and Ulmus, which grew on the surface of the bog stabilized by the mineral cover ( Figure 2A). Forest dieback which subsequently took place here could be related to the drastic changes in the climate (cooling, dampness) at the beginning of the Little Ice Age (see Grove 1988).

CONCLUSIONS
On the basis of dendrochronological analysis and wiggle-matching dating of subfossil tree trunks extracted during peat exploitation of the Podemszczyzna fen (SE Poland), which was formed in the Older Dryas Stadial, two short floating chronologies were elaborated: bog Pine chronology and deciduous trees chronology. Older, bog pine floating chronology is 147 years long and indicates that pine trees used to grow on the Podemszczyzna peatland area during the period spanning ca. 9980-9830 mod. cal BP. It is the oldest bog pine chronology reported in the European peatlands so far. The floating chronology of deciduous trees (Oak, Elm) is 139 years long and comprises a time range of 680-545 mod. cal BP. In turn, it is the youngest chronology of deciduous trees growing in peatlands in Europe.
Dendrochronology and peat multi-proxy analyzes indicate that the encroachment of trees into the fen area took place during climate warming, when the groundwater level was lowered and the fen was desiccated. Trees dying off was, in turn, caused by groundwater level rising which occurred during climate humidity growth in the Holocene.
Stratigraphy of the Podemszczyzna Peatland 1571