Multi-Tool (LA-ICPMS, EMPA and XRD) Investigation on Heavy Minerals from Selected Holocene Peat-Bog Deposits from the Upper Vistula River Valley, Poland

: Peat sediments represent important environmental and climatic archives, as well as recording information on the processes affecting the formation of these deposits; combined these data can be used for paleoreconstruction of peat- bogs. In this paper we characterize heavy mineral-rich sandy layers from two peat- bog sites in Mizerów and Strumień ( Poland ). In both cases, the most common identified mineral suite is: epidote, staurolite, tourmaline (dravite and sch ö rl), garnet, sp inel, Al 2 SiO 5 polymorphs (sillimanite, kyanite, andalusite), amphibole (mainly hornblende), pyroxene (e.g., richterite, diopside), perovskite, topaz, cordierite, apatite, monazite, chromite, ilmenite, chlorite, iron oxides, rutile and siderite. This minera l suite is characteristic of a metamorphic aureole surrounding a magmatic body. Pyrite is likely authigenic in origin. Apatite and monazite were employed for U - Pb and CHIME dating, respectively. Based on the U - Pb age information composition and textural features of selected minerals, different provenance areas were indicated: the Tatra Massif, the Bohemian Massif, and the Silesian Basin area. Transport of the investigated mineral phases was linked to development of both the Odra (praOdra) and the Vistula valleys.


Introduction
Peat-bogs are the important components of temperate landscapes, which form under particular environmental conditions. Moreover, they are associated with specific vegetation types. Peat-bogs are regarded as archives of past vegetation changes associated with both the regional climate and the evolution of the local environment. This type of biogenic reservoir may be exposed to external factors, including periodic hydrological changes, fires, etc. When peat-bogs are frequently flooded, allochthonous sandy material may be delivered into the bog system. Another common process that can supply detritus is aeolian activity [1].
Analysis of heavy minerals in Quaternary sediments is used to determine local changes in transport direction e.g., [2,3], mode of sediment transport, as well as the depositional processes and the environment of accumulation [4][5][6][7][8]. For accurate interpretation of the heavy minerals present in Quaternary sediments, it is important to determine their resistance to chemical weathering and mechanical abrasion [9,10] due to the variable preservation of mineral grains [11][12][13][14].
For this study we used quantitative X-ray diffraction analysis, scanning electron microscope (SEM) imaging, electron microprobe analysis (EMPA), as well as apatite U-Pb dating. This methodology provides a multi-tool approach for understanding the composition of detrital heavy minerals in peat-bogs and is used to constrain their provenance area. Moreover, we present the first mineralogical investigation of peats located close to the "Moravian Gate" (Figure 1), which is located between the Carpathian Mountains range (the Tatra Massif) in the east and the Bohemian Massif (the Sudety Mountains) in the west. The selected area of study has never been investigated in terms of geochemistry, as well as paleobotany. The results of this study provide information on factors affecting the accumulation of peat, determining local changes in the type of the allochthonous material transport in the Upper Vistula Valley.  [15] with marked sampling areas. All magmatic rock localities (the teschenite-picrite association in the Silesian Basin) after Włodyka [16].

Geology
The Strumień and Mizerów peat-bogs are located in the Oświęcim Basin in the Upper Vistula River Valley [17] (Pszczyna County, Poland). The studied area belongs to the Fore-Carpathian Depression. The basement contains Carboniferous organic-rich sedimentary rocks, overlain by Miocene sediments (clays intercalated with sands). The youngest sediments comprise Quaternary fluvioglacial sands and gravels, river-profluvium sands and gravels, limnic-glacial sediments, peat (Holocene in age), which are locally overlain by sandy loess and loess-like sediments in the northern part of the study area ( Figure 1; [18]). Quaternary glaciations (Odranian, Sanian 1, Sanian 2) [19] changed the flow of rivers, and also influenced the local depositional environments. In the study area, an alluvial fan occurs on the Vistula River. The Vistula River starts in the Carpathian Foothills, where it is up to 3-4 m deep and is up to 10 km wide. Its total fan length is 23 m [20,21]. Palynological analysis and radiocarbon dating suggest that the alluvial fan on the Vistula River and other Carpathian-derived rivers formed mainly in the Pleniglacial of the last glaciation, while the older proluvial remnants are preserved as residues [21]. Fossil peats (e.g., Drogomyśl, Chybie, Pierścieniec) occurring within the alluvial fan of the Vistula River represent interstadials yielding a wide range of ages from >45 ka BP to the early Holocene [22,23].

Peat Samples
Peat sampling was undertaken between August and October 2015 in two bogs (the Strumień and Mizerów sites). Sampling was performed with a Instorf-Eijkelkamp Ø 80 auger (Eijkelkamp Soil & Water, The Netherlands). The type of sediment and peat units were distinguished using Tolpa's classification [25,26]. The degree of peat decomposition was defined using the von Post scale [27].
The Strumień peat-bog has an area c. 0.9 ha, and is covered by wetland vegetation and sparse Alnus trees. The sampling profile was 1.20 m thick. In the upper part of the profile, a clay layer (0.00-0.20 m thick; Table S1) is present, with identified plant macro-remains (Carex canescens, Carex nigra, Carex rostrata, Juncus sp.; Table S2). Layers of sand were identified at 0.55-0.60 m and 0.95-1.00 m depth ( Figure 2). In the Strumień peat-bog profile, mainly Magnocaricioni/Carici and Magnocaricioni/Carici-Phragmitetii peats were deposited in rush communities. In the interval 0.8-0.9 m, as a result of the change of vegetation in the bog, Alnioni/Alneti peat was deposited. Gyttja (mud formed from the partial decay of peat) was identified in the base interval of 1.00-1.20 m ( Figure 2).
The Mizerów peat-bog covers an area of 136 ha, and is overgrown by wetland vegetation, and is partly surrounded by both forest and farmland. Clay-rich gyttja was noted in the upper part of the profile (0.00-0.60 m), with a zonation of plant macro-remains (Table S3). In the 0.00-0.30 m interval, the identified species (Table S4) include: Alnus glutinosa, Carex sp., Drepanocladus sp., Phramites australis, Sambucus sp., and Urtica sp. In the interval between 0.30 and 0.50 m, Asteraceae sp., Nuphar pumila, Drepanocladus sp., Sambucus sp. were distinguished. In the lower part of the gyttja horizon are plant macro-remains, Nuphar pumila, Poaceae sp. Potamogeton natans, Schoenoplectus sp. In the interval between 0.60 and 0.80 m, the layers of the gyttja pass downwards into peat up to a depth of 1.90 m. In the profile, three layers of sand were noted at 1.00-1.05 m and 1.50-1.55 m depth, as well as in the floor part of core ( Figure 2). In the Mizerów peat-bog profile two types of peat were deposited: Magnocaricioni/Carici-Phragmitetii and Magnocaricioni/Phragmitetii peat, whose formation is related to rush vegetation.

Heavy Mineral Separation
High-density accessory minerals were separated from the sandy layers by using standard density-separation techniques, including: crushing, hydrofracturing, washing, Wilfley shaking table, Frantz magnetic separator and handpicking. The separation was undertaken at the Institute of Geological Sciences, Polish Academy of Sciences (Cracow, Poland).

Scanning Electron Microscope (SEM) and Electron Microprobe Analysis (EMPA)
The morphology and chemical homogeneity of the surface of selected crystals of heavy accessory minerals were investigated using a scanning FET Philips 30 electron microscope (Thermo Fisher Scientific, Waltham, MA, USA; 15 kV and 1 nA) equipped with an EDS (EDAX) detector at the University of Silesia Faculty of Earth Sciences (Sosnowiec, Poland).
Apatite mineral chemistry analyses (major and minor elements) were carried out at the Inter-Institution Laboratory of Microanalyses of Minerals and Synthetic Substances (University of Warsaw, Warsaw, Poland) using a CAMECA SX-100 electron microprobe, (CAMECA, Gennevilliers Cedex, France; 15 kV, 20 nA).
An electron microprobe was used for age determination using the Chemical Th-U-total Pb Isochron Method (CHIME) dating method of Th-and/or U-bearing minerals by Suzuki and Adachi [28,29], modified by Montel et al. [30]. For monazite mineral chemistry analyses, the operating conditions were 20 kV acceleration voltage, beam current 50 nA, beam diameter 2 μm, counting times were 200 s (2 × 100 s), and 400 s (2 × 200 s) and 600 s (2 × 300 s) for peak and background positions (in parentheses), respectively.

X-Ray Diffraction (XRD)
XRD analyses were performed on heavy minerals and powdered samples using a PANalytical X'Pert Pro MPD diffractometer (Malvern-Panalytical, Almelo, The Netherlands) powered by a PW3040/60 X-ray generator (Philips) and fitted with a 1D silicon strip detector (X'Celerator). The measurements were made using Cu Kα-radiation, with a wavelength of 0.1542 nm, an acceleration voltage of 45 kV, a current of 30 mA, and with 0.01θ step sizes between the angles of 2.5° and 65° 2θ, and a 300 s measurement time per step. The data obtained were processed using HighScore+ software and the ICSD database. All XRD analyses were performed at the University of Silesia, Faculty of Earth Sciences (in Katowice, Sosnowiec, Poland).

U-Pb Laser Ablation-Inductively Coupled Plasma Mass Spectrometry (LA-ICP-MS)
Apatite U-Pb data were acquired at the Trinity College Department of Geology (Dublin, Ireland). Twenty-eight isotopes ( 31 P, 35 Cl, 43  were acquired using a 47 μm laser spot, a 4 Hz laser repetition rate, and a fluence of 3.31 J/cm 2 . A c. 1 cm sized crystal of Madagascar apatite which has yielded a isotope dilution-thermal ionisation mass spectrometry (ID-TIMS) concordia age of 473.5 ± 0.7 Ma [31,32], was used as the primary apatite reference material in this study. McClure Mountain apatite from syenite (the rock from which the 40 Ar/ 39 Ar hornblende standard MMhb is derived), was used as a secondary standard. The McClure Mountain syenite has crystallization age (weighted mean 207 Pb/ 235 U age of 523.51 ± 2.09 Ma). The NIST 612 glass standard was used as the apatite trace-element concentration reference material, and a crushed aliquot of Durango apatite that has been characterised by solution quadrupole-ICP-MS analyses [33], was used as the apatite trace-element secondary standard.
The raw isotope data were reduced using the "Vizual Age" data reduction scheme (DRS) of Petrus and Kamber [34], within the freeware IOLITE package of Paton et al. [35]. User-defined time intervals are established for the baseline correction procedure, to calculate session-wide baselinecorrected values for each isotope. The time-resolved fractionation response of individual standard analyses is then characterised using a user-specified down-hole fractionation correction model. The data reduction scheme then fits this appropriate session-wide "model" U-Th-Pb fractionation curve to the time-resolved standard data and the unknowns. Sample-standard bracketing is applied after the correction of down-hole fractionation to account for long-term drift in isotopic or elemental ratios by normalizing all ratios to those of the U-Pb reference standards. Common Pb in the apatite standards was corrected using the 207 Pb-based correction method using a modified version of the VizualAge DRS that accounts for the presence of variable common Pb in the primary standard materials [36]. During the analytical session, McClure Mountain apatite ( 207 Pb/ 235 U ID-TIMS age of 523.51 ± 2.09 Ma; [37] yielded a U-Pb Tera-Wasserburg concordia lower intercept age of 517.6 ± 4.2 Ma with an MSWD = 1.4. The lower intercept was anchored using a 207 Pb/ 206 Pb value of value of 0.88198 derived from an apatite ID-TIMS total U-Pb isochron [37].
The investigated crystals of Al2SiO5 from both sites do not show any etching textures. The grains are spherical, and in most cases form grains up to 200 μm in length (Figure 3a-c). Sillimanite is present as aggregates of small and elongated crystals (aspect ratio 20:1; Figure 3c). Kyanite and andalusite were not recognized by SEM and EDS, but their presence was confirmed by XRD (Table S5). Al2SiO5 polymorphs were found only in the Miz 1 and Miz 2 samples.

Staurolite
Staurolite commonly occurs in all samples and forms spherical grains, which are up to 100 μm in diameter (Figure 3d,e). A few grains have surfaces dominated by mechanically produced features, such as breakage blocks and conchoidal fractures. EDS does not reveal any differences in the chemistry of the staurolite minerals. This mineral is characterized in all samples from the Mizerów site.

Garnet Group Minerals
The typical size of garnet grains is no bigger than 200 μm in diameter. Several garnet grains have smaller-scale (down to ~50 μm) surface roughness or scalloping (Figure 3f,g). Generally, the garnets do not have any faceted or conchoidal surfaces. EDS investigation show that the whole garnet population consists of spessartine, pyrope and almandine species, as well as hydrogarnets. The garnet group minerals were noted in all samples from Mizerów and the Str 1 sample.

Amphibole
Amphiboles were found in all studied samples. The most common type is hornblende, with a content up to ~25% in samples Miz 1 and Miz 2. The Miz 1 sample is characterised by ca. 15% actinolite. In the sample Miz 3, potassium-rich richterite is observed. In most cases, the amphibole grains are well rounded. Sometimes, the grains reveal faceted or conchoidal surfaces. A few kaersutite grains has been found in all the studied samples. Their presence was confirmed by EDS. In most cases, the kaersutite grains are up to 80 μm in size and had conchoidal edges.

Tourmaline
XRD investigations show that dravite dominates over schörl in all studied samples. Only in one sample (the Miz 3 sample), the schörl content reached 5 vol. %. The tourmaline grain population may be divided into two morphological types: prismatic grains with somewhat rounded edges, and rounded prisms (Figure 3l,m). The size ranges from 50 μm up to 220 μm. Tourmaline grains, especially the prismatic ones, are characterised by aspect ratios of ~3:1.

Epidote Group Minerals
The grain habit of epidote in ythe investigated samples is prismatic, and sometimes exhibit slender, prismatic forms ( Figure 3n). Only in a few cases was the epidote well-rounded. The grains vary in size from 80 μm to 150 μm and were found in the Miz 2 and Miz 3 sample.

Apatite
Two different morphological types of apatites were found only in the Strumień peat-bog. All grains were handpicked from the heavy mineral concentrate. The first apatite type is composed of elongated crystals or its fragments. Locally, these apatite grains show strong abrasion, but in most cases the grains reveal sharp edges, and well-defined crystal faces ( Figure 4A). Their size varies from c. 10 μm to 80 μm. The second apatite type is characterised by spherical grains. They are up to 50-80 μm in diameter ( Figure 4C). Both types were used for REE and U-Pb age determinations. The first type of apatite is characterized by the REE patterns typical of igneous apatites e.g. [38][39][40], with LREE predominating over the HREE (Table S6; Figure 4B). The first sample from the Strumień site is characterized by strong REE fractionation (avg. (Ce/Yb)N = 21.57), positive Eu anomalies (avg. Eu/Eu* = 1.08) and slightly negative Ce anomalies (avg. Ce/Ce* = 0.98). Only the first type gave a meaningful result from U-Pb LA-ICP-MS dating, yielding an anchored lower intercept age of 116 ± 12 Ma (MSWD = 3.1; Table S7; Figure 4A). Dating of the second type of apatites in the samples failed to give meaningful data (Table S7, Figure 4C). However, the REE fractionation in these apatites is characterised by strong fractionation (avg. (Ce/Yb)N = 56.47), and slightly negative Eu anomalies (avg. Eu/Eu* = 0.81) and positive Ce anomalies (avg. Ce/Ce* = 1.01; Table S6, Figure 4D).

Zircon
In the zircon population, numerous different crystal morphologies were identified (Figure 6a,b). Up to 20% of all inspected zircon crystals reveal prismatic shapes and habit (Figure 6c-e). The rest of the grains range from slightly rounded (c.a. 15%), prismatic and subhedral shapes to sub-rounded types (c.a. 10%). Sporadically, spherical forms were also noted. Zircon crystals are characteristic for all samples from the Mizerów site.

Authigenic Heavy Minerals
The authigenic heavy mineral suite is limited to a few framboidal pyrite aggregates, which were found on the margin with the peat (Figure 6f). The pyrite, both framboids and cube-like-octahedral forms, are the most common heavy mineral in the surrounding peat units.

Transport Type
During the Quaternary, the most significant environmental and geographic changes in the regional drainage basins took place with the development of the contemporary fluvial network of the upper Vistula River Valley [21]. Based on pollen analysis of the Strumień site peat, the subatlantic period has been recognized. However, data obtained by Pazdur [41] in the Pszczynka valley indicate the accumulation of peat in this area for the entire duration of the Holocene [42]. The presence of the locally accumulated layers of sand in the studied profiles is likely related to a floodplain environment dominated by wetlands and riverside marshes. These environments are suitable for peat bogs and phytogenic deposit formation [43], with the layers of sandy material in the peat indicating periods of flooding. Isolated and numerous ephemeral floods, as well as periods of sustained rainfall were characteristic in this area during Holocene and caused significant changes in the river valleys [44]. There were no significant climate fluctuations in both the subboreal and subatlantic periods [45]. During these intervals, palynological data show an increase in the presence of vegetation as a result of human activity. The pollen analysis for the Strumień site showed dense but treeless vegetation, which may also inhibit aeolian deposition. Horizons with high concentrations of heavy minerals accumulated mostly through alluvial transport but aeolian transport cannot be excluded as a different process leading to the formation of heavy mineral-rich sands. Alluvial reworking of aeolian deposits that concentrate heavy minerals can produce rounding of grains, such as that observed in most samples (Figures 3 and 6). In addition, the lack of muscovite (Sample Miz 1, Miz 3, and Str 2), and chlorite (Sample Miz 1 and Miz 2), defined by Chlebowski and Lindner [46] as minerals susceptible to deflation and aeolian transport, is an indication that aeolian processes are not predominant. In the literature, an increase in tourmaline content is associated with aeolian processes [8]. However, in the results obtained only in the sample Str 1 (depth 0.50-0.55 m), tourmaline content reaches a value above 20%. There are key differences between the studied heavy mineral-rich sediments from both the Strumień and Mizerów sites (Figures 3 and 6). Actinolite, garnet, ilmenite, andalusite, chlorite and zircon dominate in the Mizerów site, while tourmaline, hornblende, rutile, muscovite and magnetite are characteristic for the Strumień site. Generally, at both peat-bogs studied, the characteristic features of the heavy mineral suite can be used for catchment area reconstruction. This can be done in case of sediments, where less resistant primary components (e.g., feldspar, olivine), have been eliminated by weathering processes.
Staurolite and Al2SiO5 polymorphs, when present, appear to be the most suitable for assessing the intensity of weathering in the peat sediments. This, in turn, may allow for an assessment of the other minerals. For instance, strongly weathered staurolite in a sample likely indicates that minerals that are more susceptible to weathering, such as topaz or corundum, were eliminated. The presence of staurolite, often associated with andalusite and kyanite, in regionally metamorphosed rocks, and provides information on sand provenance. This situation is characteristic for the Mizerów site. Moreover, staurolite, kyanite, andalusite, zircon, tourmaline and hornblende, especially from the Strumień site, appear to have been resistant to acidic leaching.
The distribution of zircon is variable in the studied sites. The presence of rounded and subrounded zircon may indicate a relatively strong reworking, with most grains probably coming from igneous or metamorphic rocks. Provenance areas include from the Bohemian Massif or the Tatra Massif. No zircon crystals were found in the Strumień site.
The highest mean abundance of rutile was found in the Strumień site (up to 22 vol. %), while the Mizerów site sample is characterised by lower rutile content (up to 5 vol. %; Table S5). The rutile population is represented by grains with differences in size and shape. Some of the grains show typical bladed or prismatic characters, with well-developed cleavage traces or less common leaching edges. The presence of rounded rutile indicates reworked sediments, most probably coming from metamorphic source rocks.
The epidote population, in almost all investigated grains, lacks corrosion textures, with grains retaining typical bladed or prismatic habits, with well-developed cleavage planes. In some cases, early stages of epidote corrosion are manifested by oriented pits, but as dissolution proceeds, epidote develops large-scale 'hacksaw' terminations. This situation has been described in many sediments [47][48][49]. Epidote can form as a rock-forming mineral, typically as products of regional metamorphism, and in some felsic igneous rocks as well as in hydrothermally altered rocks. Epidote group minerals are common products in calcareous sedimentary rocks in contact metamorphic aureoles. They can also be found in pegmatites and granites as primary magmatic phases, as well as secondary phases (e.g., plagioclase-saussuritization) e.g. [50,51].
Phosphates are represented by (Ce)-monazite and fluorapatite. Both species are common accessory minerals in igneous and metamorphic rocks of different chemistry and provenance. Phosphates, especially monazite and apatite, are suitable for geochronology. The monazite population seems to have a magmatic origin as it exhibits internal magmatic zoning. The main monazite CHIME age ( Figure 5), may be interpreted as the Variscan magmatic event, which is widely observable in numerous nearby plutons [52][53][54], and references therein. The obtained U-Pb apatite LA-ICP-MS age, chemical homogeneity of the grains and REE profiles implies a magmatic source ( Figures 4A,B). The origin of the second apatite population is less clear (Figures 4C,D). The REE patterns are more characteristic of metamorphic or hydrothermal apatite cf [55,56].

Heavy Mineral Provenance Area Characterization
Two different source rock areas, the Tatra Mountains and the Sudety Mountains, are indicated for both investigated peat-bog sites. The Tatra Mountains are an example of a complex crystalline massif, composed of polygenic granitoid intrusions which are partly covered by its metamorphic envelope in the western region [52 and references therein]. The metamorphic cover is composed of gneisses, schists, amphibolites, and locally also eclogites [53,57] and references therein. This massif, uplifted during Alpine orogenesis [58], could be the source of detrital minerals in this study, which are eroded and transported away from the Carpathians by rivers. The Vistula and the Danube rivers are the largest rivers that transport material from the Tatra Mountains to the north and the Carpathians to the south of the sampling sites respectively.
The detrital mineral suite, produced by weathering and erosion of rocks from the Tatra Mountains, is dominated by magmatic components such as zircon, tourmaline, apatite. However, because the metamorphic envelope is mostly located on southern part of the crystalline core in Slovakia, all detrital metamorphic grains sourced from the Tatras (kyanite, sillimanite, garnet, staurolite, amphibole (mostly hornblende), would likely be transported to the south.
The Bohemian Massif is one of the largest continuously outcropping fragments of the extensive Variscan basement that crops out beneath the Variscan platform and younger sediments. The Bohemian Massif is a complex geological area with numerous granitoid intrusion types and metamorphic units. Due to its size and tectonic character, several regions differ in lithology and age and can be readily distinguished. They are as follows: the Moldanubian, Teplá-Barrandian, Saxo-Thuringian (subdivided by the younger Elbe Fault Zone into the Krušné Hory mountain zone and the Lugicum-West Sudetic Zone) and Moravo-Silesian units e.g. [59]. The Moravo-Silesian Unit is the nearest unit and is located to the west from both investigated peat sites. The eastern region is the Bardzkie Mountains, which are a part of the Kłodzko-Złotostocki Massif (the Central Sudetes). The whole massif is a potential source for a diverse suite of heavy minerals, but in the most eastern part, metamorphic assemblages dominate.
The third possible supply area for the studied material is the teschenite-picrite association (TPA) rocks. Generally, these rocks are exposed between the Bohemian Massif (the Sudetes) and the Tatra Massif. This early Cretaceous alkaline volcanic region, of which the TPA of volcanic rocks are associated, is unique to the western part of the Outer Western Carpathians. This magmatic province is 15-25 km wide and extends in a NE direction for over 100 km from Hranice in Moravia, the Czech Republic, to Cieszyn and Bielsko-Biała in Poland (the Silesian Basin; Figure 1). The TPA contains a wide range of intrusive rocks, which belong to two main groups of alkaline (teschenite, picrite, syenite, lamprophyre) and sub-alkaline (typically dolerite) igneous rocks. The most common heavy minerals related to the TPA are: diopside, clinopyroxene, amphibole (kaersutite), mica (annitesiderophyllite), chlorite, haematite, apatite and titanite. The timing of volcanic activity of the TPA in the Silesian Basin took place from 128 to 103 Ma, and most likely peaked between 128-120 Ma [60][61][62][63].

Development of Relief and River System in Neogene and Pleistocene
Heavy minerals are one of the best indicators of the origin of terrigenous detrital material, and allow for the determination of provenance areas for sediments whose primary components, especially those less resistant to weathering, have been eliminated by weathering processes. In the studied area, the most crucial factor for determining the mineral transport directions was the regional river valley evolution. The transport and accumulation of detrital material occurred in several phases within the Moravian Gate and its surrounding areas. In the initial phase of river formation in the Pleistocene, the main direction of fluvial outflow was dominated by a river (pre-Odra) [64] on the Fore-Sudetic Block, which was flowing from the present Silesian Uplands and the Eastern Sudetes. The river flowed broadly north-east, meandering over the Fore-Sudetic block and sporadically overflowing its flood plain.
In Pliocene and Early Pleistocene, the hydrographic network was modified. The newly created Odra Valley became the drainage area for the Carpathian basin. Numerous modifications in the drainage basin caused by glacial transgressions lead to changes in the flow directions of all the Carpathian rivers. The modern northern Vistula tributaries (Pszczynka, Gostynka) also experienced a modification of its drainage network (Figure 7). The studied area of the Oświęcim Basin is located within the region of the South Polish Glaciations (Elsterian), where two episodes of the Saalian I stage [64] are fingerprinted. The climate changes (mostly cooling) during the Dryas stage was responsible for the river debris increasing in the Carpathian area [64]. During the next stage (Early Pleistocene?) which was characterized by intense landscape weathering, the regional river system was formed, and has continued in this configuration with minor changes up to the present day. The main river channel (pre-Odra) in the western part of the study area shifted north, probably beyond the line of the present-day Odra channel, while a system of local river valleys flowing from the Sudetes developed on the Fore-Sudetic block, and rivers drained the Sudetes Mountains [63,[65][66][67][68]. River network modifications of the whole Vistula River Valley occurred at the Vistulian/Holocene boundary. Subsequently annual flood-thaw cycles were replaced by more common rain floods. At this time, braided rivers changed into meandering fluvial systems [64,69]. Starkel [70][71][72][73] indicated that climate changes in Holocene supported cyclic inundation of the Upper Vistula River Valley, which could be responsible for redeposition and/or mixing of all previously deposited sediments. The present-day Vistula River represents mostly a braided river, overcharged with sediments [69], mainly as a result of human activity in the catchment areas.

Conclusions
This paper presents a detailed characterisation of heavy minerals from peats in the Carpathian Foredeep in Poland.
(1) The studied heavy mineral suite reveals up to three different source areas, which were confirmed by mineral composition and age.
(2) These different source areas are: (a) the Tatra Massif, (b) the Bohemian Massif, and (c) the Silesian Basin area. However, some parts of the mineral assemblages may point to the South Polish Glaciations (Elsterian) sediments as the source. Originally, this material could have been deposited in the Oświęcim Basin by glaciers.
(3) The studied heavy mineral suite as well as sedimentology indicates an alluvial origin related to both the Odra (praOdra) and the Vistula valley development. Any aeolian component for the provenance of the material is insignificant.
(4) Numerous mineral grains reveal characteristic dissolution features caused by long-term stagnation in a water-rich environment.

Supplementary Materials:
The following are available online at www.mdpi.com/xxx/s1, Table S1: Simplified lithological profile including peat classifications from the Strumień site [25,27], Table S2: Plant macro-remains identified from the Strumień site, Table S3: Simplified lithological profile including peat classifications from the Mizerów site [25,27], Table S4: Plant macro-remains identified from the Mizerów site, Table S5: Mineral composition based on X-ray diffraction of investigated heavy fraction from the sandy layers from the Strumień and Mizerów sites, Table S6: Representative LA-ICPMS apatite REE analyses from this study, Table S7: Representative LA-ICPMS U-Pb apatite data for samples from this study, Table S8: Representative EMPA monazite data and their chemical formulae based on 4O 2− .
Author Contributions: Krzysztof Szopa, Sylwia Skreczko, and Artur Szymczyk wrote the paper. Sylwia Skreczko and Krzysztof Szopa collected the studied material. David Chew and Tomasz Krzykawski obtained U-Pb and powder X-ray diffraction data, respectively. Krzysztof Szopa obtained EMPA and calculated the CHIME age. Sylwia Skreczko and Artur Szymczyk recognized the flora. All authors have read and agreed to the published version of the manuscript.
Funding: This work was financially supported by the Institute of Earth Sciences at the Faculty of Natural Sciences, University of Silesia in Katowice, Poland. DC acknowledges support from research grant 13/RC/2092 from Science Foundation Ireland which is co-funded under the European Regional Development Fund and by PIPCO RSG and its member companies.