The final part of the Late Pleniglacial and Late Glacial periods took place when the areas of the Central European Lowlands were subject to fundamental environmental evolution resulting from climate change and the retreat of the Scandinavian Ice Sheet (SIS). Rapid changes to terrestrial sedimentary environments are reflected in sedimentary succession called fluvio-aeolian to aeolian (Kasse, 2002 and Zieliński
Cold periods of the Late Glacial (Oldest, Older and Younger Dryas) were separated by periods of climate warming (Bølling-Allerød interstadial). Warm periods were responsible for creating appropriate climatic conditions for soil-formation processes (Manikowska 1991). During cold intervals, intensive sedimentation of aeolian sediments took place in the ESB (Manikowska, 1995; Goździk, 2007; and Zieliński, 2016). Under debate is the question about the significance of particular cold periods in aeolian activity and the formation of dune fields (Nowaczyk, 1986; Kasse, 2002). However, what is also of importance is the fact that the duration of the most well-known cold periods such as the Younger Dryas are comparable with those periods, spanning several centuries, which can be dated precisely using dosimetric (mainly OSL luminescence) techniques. This simple fact significantly contributes to issues pertaining to the interpretation of the results obtained for the sediments from these cold periods.
One of the crucial problems is the accurate dating of particular periods within the Late Pleniglacial and the Late Glacial. Between the optically stimulated luminescence (OSL) and radiocarbon methods, there are frequently significant discrepancies (Kalińska-Nartiša
One of the areas with fluvio-aeolian sediments is the Niemodlin Plateau, SW Poland (Fig. 1). Two sites with palaeosols were investigated in this area. A sedimentary succession of sandy deposits provided an opportunity to estimate the time scheme of their deposition as well as forming a comparison with the latest results from other parts of the ESB.
Fig 1
Microregions of the Niemodlin Plain and location of investigated sites, Siedliska and Przechód, with digital elevation model (DEM).

The first researcher to study dune structures in Opole, Silesia, was Pernarowski (1968). He performed extensive field research and basic laboratory tests such as granulometric analysis based on which he concluded that the nature of the dune formation in this area was highly complex. This complexity, according to Pernarowski, resulted from the large changeability of local wind-direction. These winds have repeatedly changed their direction and speed over the last several thousand years. What is interesting, however, is that despite examining 24 dune complexes in the Opole region, Pernarowski found the remains of fossil soils in only two places. It is worth emphasizing that no research was performed in this particular area after Pernarowski, which means that our particular investigation is the first research in this region conducted using modern scientific methods. Our initial research aims at referring to Pernarowski’s theory, thus the main goals of this article are: (i) to reconstruct sedimentary environments during the Late Pleniglacial and the Late Glacial periods; (ii) to determine a stratigraphic chart for sedimentary processes based on OSL and radiocarbon dating; and (iii) to compare the obtained results with research from other parts of the ESB.
The investigated sites (inland dunes) are located on the Niemodlin Plateau and form part of the Niemodlin Plain (Solon
Quaternary units are underlined by the Gozdnica Formation, a Mio-Pliocene piedmont alluvial sedimentary cover (Piwocki and Ziembińska-Tworzydło, 1997). The Niemodlin Plateau is mainly built of glacial sands and gravels deposited during the Odra Glaciation (Przybylski and Badura, 2001). The inland dune areas occur mostly in the western part of the Niemodlin Plateau close to the morphological edge (approx. 20 m high) of the Ścinawa Niemodlińska River valley. The biggest and the best-developed dunes are located in SW part of the Niemodlin Plateau.
The dunes of the Niemodlin Plateau adopt mainly parabolic, longitudinal and irregular hilly forms. The parabolic dunes are mainly U-shaped and V-shaped. These create compound hemicyclic, digitate or nested dunes (according to Pye and Tsoar, 2009). The height of these forms does not usually exceed 10−12 m. Asymmetric development of the dune arms, length and their orientation is common and clearly visible. Northern arms of dunes are longer and orientated NW−SE. There are often longitudinal dunes on the shoulder extension of the arms (up to 3−5 m high). Most longitudinal dunes seem to be transitional forms between parabolic and longitudinal (Vandenberghe 1983). The sites are located in the area of most typical for parabolic dunes.
The Siedliska profile is located in the well-developed northern parabolic dune arm in the immediate neighborhood of the nose. The shoulder height here is 3.5 m, the length is 1200 m and has a NW−SE orientation. The southern arm, like the nose, is highly reduced. The height is about 1 m, with a length of 600 m and a W−E orientation (Fig. 1).
The Przechód profile represents a compound dune, which is a combination of complex forms such as hemicyclic and nested (Pye and Tsoar 2009). The main, outer ridge of the form, over 2.6 km long, has a height of 2 m in the southern part, 10 m in the northern and nose. The dunes ‘inside’ the main form are smaller − 2−7 m high and 200−400 m long along the ridge. The profile is located in the N−E part of the dune, in the area connecting the arms of the basic dunes (Fig. 1).
Two excavations were prepared for further investigations in a forested area near the villages Siedliska and Przechód. Twenty-seven samples were taken from the Siedliska and 21 from the Przechód profile for analysis in terms of grain size and other lithological properties. Nine OSL samples were collected and three samples of charcoals (radiocarbon dating) were gathered for further measurements (Fig. 2).
Fig 2
Sedimentary succession, radiocarbon and OSL ages of deposits at Siedliska (A) and Przechód (B) sites.

To establish the origin of the sediments, the three-level lithofacies analysis (lithofacies, lithofacial association and lithofacial complex) was used based on the methodology proposed by Miall (1978, 1996) with subsequent modifications (Zieliński, 2014). For aeolian sediments, the methodology proposed by Zieliński (2016) was implemented. Measurements of directional elements, analyses of cryogenic structures and fossil soils were also performed. Soil horizon symbols used for morphological description of palaeosols are in accordance with the international standards (FAO 2006).
The procedures included the collection of samples for luminescence dating from a clean vertical section using thin-walled steel pipes from an Eijkelkamp system. At the same time, about 1 kg of sediment was collected into plastic bags from around the tubes, for gamma spectrometry.
For OSL measurements, coarse grains of quartz (125–200 μm) were extracted from the sediment samples by routine treatment with 20% hydrochloric acid (HCl) and 20% hydrogen peroxide (H2O2). The quartz grains were separated using density separation with the application of sodium polytung-state solutions leaving grains of densities between 2.62 g/cm3 and 2.75 g/cm3. The grains were sieved, before etching with concentrated hydrofluoric acid (HF), for 40 min. The grains were then placed on small steel disks (10-mm diameter) using silicon oil and placed into the luminescence reader.
All OSL measurements were performed using an automated Risø TL/OSL DA-20 reader fitted with a calibrated 90Sr/90Y beta source delivering about 7.0 Gy/min to grains at sample position and a 6-mm Hoya U-340 filter was used for OSL detection.
Equivalent doses were determined using the single-aliquot regenerative-dose (SAR) protocol (Murray and Wintle 2000). The OSL SAR protocol which was used in our measurements contained the following steps:
Irradiation with the regenerative beta dose Di
Pre-heat at a temperature of 260°C for 10 s
Blue light stimulation at the temperature of 125°C for 100 s
Irradiation with the test dose Dt (10% of the natural dose, but not less than 2 Gy)
Cut-heat at a temperature of 220°C and
Blue light stimulation at a temperature of 125°C for 100 s.
Intensities measured in Steps 3 and 6 were used for equivalent dose determination. For equivalent dose calculation, the first second of the signal was used and the background was estimated from the last 10 s. The SAR dose-response curves were best represented by a single saturating exponential. Final equivalent dose (De) values were calculated using the Central Age Model (CAM) (Galbraith
To assess the dose rates that arise from decay chains and potassium we used high-resolution γ spectrometry (HRGS) and the μDose system. Both systems were calibrated with the same reference materials, namely IAEA-RGU-1, IAEARGTh-1 and IAEA-RGK-1 obtained from International Atomic Energy Agency reference materials (IAEA 1987).
The HRGS measurements were performed using a low background, high-resolution HPGe detector with a resolution at full width and half maximum of 1.8 keV and relative efficiency of 40% at the energy of 1332 keV manufactured by Canberra. The lines 234Th, 234mPa, 214Pb, 214Bi and 210Pb were used to measure 238U decay chain, and 228Ac, 212Pb, 212Bi and 208Tl lines were used for 232Th decay chain.
For luminescence dating, coarse grains of quartz were used with standard measurement procedures, however, when determining the annual dose of sediment, we also used, in addition to the standard high resolution γ spectrometry, the new μDose system. Dose rates, from both systems, were calculated using software that is built in the μDose system (Tudyka
Basic information about investigated samples. Dose rate values obtained with different system detectors and final Bayesian estimation. Final obtained age using CAM model.
Sampling depth | Dose rate (Gy/ka) | ||||||
---|---|---|---|---|---|---|---|
Lab.code | Sample ID | (cm) | Water content (%) | Gamma spectrometry | μDose system | Bayesian averaging | OSL age (ka) |
GdTL-2638 | Przechód_12 | 120 | 7±3 | 1.06±0.08 | 0.88±0.07 | 0.98±0.05 | 11.85±0.91 |
GdTL-2639 | Przechód_36 | 360 | 7±3 | 0.90±0.07 | 0.89±0.07 | 0.89±0.04 | 12.68±1.00 |
GdTL-2640 | Przechód_45 | 450 | 7±3 | 0.81±0.06 | 0.58±0.06 | 0.74±0.04 | 14.88±1.11 |
GdTL-2641 | Przechód_60 | 600 | 7±3 | 0.74±0.05 | 0.62±0.06 | 0.70±0.04 | 18.96±1.46 |
GdTL-2642 | Siedliska_150 | 150 | 7±3 | 0.94±0.07 | 0.85±0.07 | 0.92±0.05 | 12.39±0.98 |
GdTL-2643 | Siedliska_320 | 320 | 7±3 | 0.90±0.07 | 0.91±0.07 | 0.90±0.05 | 12.47±1.04 |
GdTL-2644 | Siedliska_365 | 365 | 7±3 | 1.33±0.10 | 1.34±0.08 | 1.35±0.06 | 13.21±1.04 |
GdTL-2645 | Siedliska_425 | 425 | 7±3 | 1.11±0.09 | 1.08±0.07 | 1.10±0.06 | 14.26±1.27 |
GdTL-2646 | Siedliska_450 | 450 | 7±3 | 1.37±0.11 | 1.33±0.08 | 1.37±0.07 | 14.50±1.21 |
The dry dose rates (Adamiec and Aitken, 1998 and Guerin
Three samples were collected for radiocarbon dating. Two samples from Przechód and one from the Siedliska profile. From Przechód, about 3 kg of sediment with tiny pieces of charcoals were collected for each sample. (We used LSC technique for measurements, so we needed much more organic material in comparison to the AMS dating.) In the laboratory, this material was treated with the standard procedure acid-alkali-acid (AAA) to remove humic acids, chitin, fungal products, etc. and next converted to benzene for liquid scintillation counting (LSC). LSC measurements were performed using a Quantulus 1220 spectrometer. Detailed calibration procedure (Pawlyta
For the sample from Siedliska, it was possible to collect some small pieces of charcoals for AMS dating. The typical chemical pre-treatment based on the AAA method before radiocarbon dating was used. The AAA pretreatment involves rinsing the samples in hot hydrochloric acid (HCl) followed by sodium hydroxide (NaOH). In the radiocarbon accelerator mass spectrometry (AMS) technique, the 14C concentration is measured in graphite prepared from the carbon contained in the sample. For this purpose, a line for preparation of graphite targets in the Gliwice 14C laboratory was used. All 14C results were calibrated using the OxCal program v4.2 (Bronk Ramsey, 2009 and Southon
Calibrated 14C ages of charcoals. The age was calibrated using OxCal program v4.2, (Bronk Ramsey, 2009 and Southon
No. | Sample name | Age 14C (BP) | Range of calendar (calibrated) age 68% confidence level | Range of calendar (calibrated) age 95% confidence level |
---|---|---|---|---|
1 | Przechód _top | 11270 ± 170 | 13320BP - 12960BP | 13450BP - 12780BP |
2 | Przechód_bottom | 10780 ± 80 | 12750BP - 12640BP | 12810BP - 12560BP |
3 | Siedliska | 10915 ± 45 | 12800BP - 12730BP | 12880BP - 12700BP |
Fig 3
Depositional succession at Siedliska (A-D) and Przechód (E-H) sites: A – general view at the Siedliska site; B – S-2 unit with large-scale involutions and ice-wedge cast (arrow); C – fluvio-aeolian deposits (S-3 unit) with the fossil soil (arrow) in top of unit; D – aeolian deposits within coversand (S-4 unit) and stoss slope of dune (S-5 unit); E – top of P-1 unit, and aeolian deposits within coversand and bottom part of moving dune lee slope; F – soil within the lee slope of dune (top of P-3 unit); G – aeolian deposits within arm of dune (P-4 unit); H – climbing ripple lithofacies (arrow), deposited up the lee slope.

Fig 4
Relative probability density functions (Berger 2010) for all samples are presented on the scale of age.

The correlation of lithofacial units from both key sites and the results of absolute dating allows distinguishing of four phases of the sedimentary development process. Local processes have been compared with the results of earlier research and correlated with global climate changes. The chronology of biostratigraphic horizons is defined following Litt
During the Late Pleniglacial age, the global climate reached a maximum in terms of cooling, which was expressed by a maximum extent of terrestrial ice-sheets in the northern hemisphere (Last Glacial Maximum – LGM; Ehlers and Gibbard, 2004). In the extraglacial areas, a wide zone of permafrost developed, which affected sedimentary processes (Last Permafrost Maximum – LPM; Vandenberghe
Different types of sedimentation (fluvio-aeolian and aeolian) took place simultaneously at both key sites. At the Siedliska site, sandy-silty rhythmic succession indicates fast and cyclic changes of sedimentary conditions. Horizontally laminated sand layers represent the supercritical flow of sheet floods. Wavy laminated silty lithofacies were deposited as a result of the migration of adhesion ripples along a wet surface (Schwan, 1986 and Kasse, 1997, 2002). This sandy-silty rhytmite is interpreted as a result of fluvioaeolian deposition. At the Przechód site, inclined stratified sands with climbing ripples at the bottom are interpreted as a distal part of migrating aeolian dune (Pye and Tsoar, 2009). Climbing ripples indicate their reverse migration at the foot of the dune slope (Fedorowicz and Zieliński, 2009). The palaeocurrent measurements suggest transport of material towards the NE. This direction indicates that it could be a small moving transversal dune or even a barchan (compare Zieliński, 2016).
There are no dated sediments between 17.5 and 15.5 ka. This may be interpreted as a stratigraphic gap caused by subsequent erosion or non-depositional conditions. This gap can be correlated with the Beuningen Gravel Bed, a marker horizon widely described from the Netherlands and Belgium. This unconformity is overlaid by a deflation pavement and represents a similar time dated at 18−15 ka (Bateman and Van Huissteden, 1999; Kasse
Climate amelioration during Allerød interstadial is expressed by palaeosol horizon in both the key sites. However, palaeosols represent different types of morphology. Palaeosol from the Przechód site is two-fold with sandy layer in-between (Fig. 3F). Soil layers have different thicknesses and reversed 14C dates (Table 2). This is interpreted as being a result of colluvial processes on the distal dune slope and redeposition due to washing out of organic matter. Therefore, it is possible to interpret this as soil of colluvial type, whereas palaeosol in the Siedliska site does not reveal disturbances. The palaeosol from the Siedliska site correlate with Usselo soil type (Kasse 1997 and Vandenberghe
During the last climate cooling of the Younger Dryas, aeolian sedimentary processes prevailed in both the key sites. Large-scale inclined stratified sands are interpreted as a deposition on the distal slope of the aeolian dune. Here, the transverse transport of material played a significant role. Such conditions prevail at the distal part of parabolic dune or its branches (Mc Kee, 1966; Borówka, 1990; Pye and Tsoar, 2009; and Zieliński, 2016). Climbing ripples observed on the distal slope (Figs. 2B and 3H), within P-4 unit, confirm our interpretation. Taking this fact into account as well as the location of the outcrop at the branch of a compound dune (Fig. 1C), it is possible to surmise that deposition took place as a result of WNW and NW wind directions.
The OSL ages obtained from S-4, S-5 and P-4 units allowed for dating of investigated sediments within the Younger Dryas (Fig. 5). This was the main phase of aeolian activity and dune formation within the area being studied. The results are in good agreement with the stratigraphic scheme from Central and Western Europe (Bohncke
Fig 5
Model of environmental changes at the turn of Weichselian Pleniglacial and Late Glacial; 18O/16O ratio after Rasmussen

The studies performed provided a stratigraphic model for sedimentary succession. The results provide a stratigraphic succession. Sedimentation took place in the final part of Late Pleniglacial and Late Glacial. However, direct absolute data was not obtained from periglacial structures of the Siedliska site, where size and frequency make it possible to be compared and correlated with LPM (Vandenberghe
fluvial succession (Erkens
The beginning and finish of fluvio-aeolian sedimentation in Western Europe took place during the Late Pleniglacial. In the Netherlands fluvio-aeolian deposits are described after Van der Hammen (1971) as Older Coversand I and dated between 20−17 ka (Bateman, Van Huissteden, 1999; Kasse
The location of study sites on interfluves is presumably the reason for the stratigraphic gap between cryogenic structure formation (LPM) and the start of fluvio-aeolian deposition. Interfluves were source areas for the most of the Late Pleniglacial, and it was only during the release of sand from upper river terraces that material was transported by the wind into these areas (Goździk, 2007).
One of the most important questions is which cold phase in the Late Glacial is mostly responsible for aeolian dune formation. Results from Western Europe reveal that pure aeolian coversand (Older coversand II) started before Bølling interstadial (Kasse
Previous models of fluvio-aeolian successions provided a scheme for comparison with data presented in the articles by Kasse (2002) and Vandenberghe
The present study used multi-proxy results to establish chronology and the types of depositional processes in the study area. Based on the results, several conclusions can be drawn:
A comparison of OSL and radiocarbon data reveal high compliance of dating results. Attention paid to OSL methodology shows that the elimination of 14C/ OSL discrepancies is possible.
LPM is visible in the well-developed periglacial complex at the Siedliska site. It documents the climate minimum at the Late Pleniglacial (24−17 ka).
The morphological position of the study area (interfluve) was presumably a reason for 17−15 ka stratigraphic gap.
The first phase of sediment accumulation took place in the time framework 15−13.2 ka (from the Oldest Dryas to the Older Dryas). There are no signs of a break-in sedimentation and development of soil processes during Bølling interstadial.
Fluvio-aeolian and aeolian (dune formation) activities developed at the same time in the Oldest Dryas.
Allerød warm phase is described in well-dated palaeosol of Usselo type.
The main dune formation phase took place during the Younger Dryas. Within this period, the initial low dunes were overbuilt. Climate amelioration at the beginning of the Holocene stopped aeolian processes. There is no evidence for younger aeolian activity.
Fig 1

Fig 2

Fig 3

Fig 4

Fig 5

Basic information about investigated samples. Dose rate values obtained with different system detectors and final Bayesian estimation. Final obtained age using CAM model.
Sampling depth | Dose rate (Gy/ka) | ||||||
---|---|---|---|---|---|---|---|
Lab.code | Sample ID | (cm) | Water content (%) | Gamma spectrometry | μDose system | Bayesian averaging | OSL age (ka) |
GdTL-2638 | Przechód_12 | 120 | 7±3 | 1.06±0.08 | 0.88±0.07 | 0.98±0.05 | 11.85±0.91 |
GdTL-2639 | Przechód_36 | 360 | 7±3 | 0.90±0.07 | 0.89±0.07 | 0.89±0.04 | 12.68±1.00 |
GdTL-2640 | Przechód_45 | 450 | 7±3 | 0.81±0.06 | 0.58±0.06 | 0.74±0.04 | 14.88±1.11 |
GdTL-2641 | Przechód_60 | 600 | 7±3 | 0.74±0.05 | 0.62±0.06 | 0.70±0.04 | 18.96±1.46 |
GdTL-2642 | Siedliska_150 | 150 | 7±3 | 0.94±0.07 | 0.85±0.07 | 0.92±0.05 | 12.39±0.98 |
GdTL-2643 | Siedliska_320 | 320 | 7±3 | 0.90±0.07 | 0.91±0.07 | 0.90±0.05 | 12.47±1.04 |
GdTL-2644 | Siedliska_365 | 365 | 7±3 | 1.33±0.10 | 1.34±0.08 | 1.35±0.06 | 13.21±1.04 |
GdTL-2645 | Siedliska_425 | 425 | 7±3 | 1.11±0.09 | 1.08±0.07 | 1.10±0.06 | 14.26±1.27 |
GdTL-2646 | Siedliska_450 | 450 | 7±3 | 1.37±0.11 | 1.33±0.08 | 1.37±0.07 | 14.50±1.21 |
Calibrated 14C ages of charcoals. The age was calibrated using OxCal program v4.2, (Bronk Ramsey, 2009 and Southon et al., 2013) with calibration curve IntCal13 (Reimer et al., 2013).
No. | Sample name | Age 14C (BP) | Range of calendar (calibrated) age 68% confidence level | Range of calendar (calibrated) age 95% confidence level |
---|---|---|---|---|
1 | Przechód _top | 11270 ± 170 | 13320BP - 12960BP | 13450BP - 12780BP |
2 | Przechód_bottom | 10780 ± 80 | 12750BP - 12640BP | 12810BP - 12560BP |
3 | Siedliska | 10915 ± 45 | 12800BP - 12730BP | 12880BP - 12700BP |