ORIGINAL RESEARCH article
Disentangling Sedimentary Pathways for the Pleniglacial Lower Danube Loess Based on Geochemical Signatures
- 1Department of Geography, Chair of Physical Geography and Geoecology, RWTH Aachen University, Aachen, Germany
- 2Romanian Academy, Institute of Speleology, Cluj-Napoca, Romania
- 3Department of Earth Sciences, Uppsala University, Uppsala, Sweden
- 4BayCEER and Chair of Geomorphology, University of Bayreuth, Bayreuth, Germany
The source of aeolian sediments such as loess has been investigated since decades. Reliable knowledge on potential dust sources is crucial to understand past climatic and environmental conditions accompanying the dispersal of early modern humans (EMH) into Europe. Provenance studies are usually performed on small sample sets and most established methods are expensive and time-consuming. Here, we present the results of high-resolution geochemical analyses performed on five loess-palaeosol sequences from the Lower Danube Basin (LDB), a region, despite its importance as a trajectory for EMH, largely underrepresented in loess provenance studies. We compare our results with geochemical data of loess-palaeosol sequences from Austria, Hungary, Serbia, and Ukraine. Based on published literature, we thus evaluate five plausible sedimentary pathways for the LDB loess: 1) the Danube alluvium (DA) pathway, which constrains the transport and re-deposition of detrital material by the Danube and its tributaries; 2) the Carpathian Bending (CB) pathway, where sediment is mainly transported from the Cretaceous to Neogene flysch of the Eastern Carpathian Bending; 3) the Eastern Carpathian (EC) pathway, in which sediment is eroded from the flysch of the Outer Eastern Carpathians, transported by rivers, and deflated by northwesterly to westerly winds; 4) the glaciofluvial (GF) pathway, where dust is deflated from glacial outwash plains in nowadays Ukraine, and 5) the Black Sea (BS) pathway, where dust originates from the exposed shelf of the Black Sea. Based on geochemical data, we consider the DA pathway to be the major sediment trajectory for loess in the LDB. Especially the sequences located close to the Danube and the Dobrogea show similarities to sites in Central and Northeast Hungary as well as Northern Serbia. For the northeastern part of the LDB, we demonstrate that dust input is mainly sourced from primary material from the Eastern Carpathians. Mineralogical estimations and geochemical data render the CB pathway as an additional substantial source of detrital material for the loess of this area. We consider the influence of the GF pathway in the LDB as negligible, whereas some minor influences of the BS pathway cannot be ruled out based on geochemical data.
Since the earliest days of loess research, the origin and source of aeolian sediments such as loess have been a matter of debate. Richthofen (1882), a pioneer in loess research, named non-vegetated deserts and desiccated lake beds as potential source areas for loess in China. Throughout the decades, the definition of loess as well as its sources and origin have been subjects of debate in numerous studies (see Smalley et al., 2001; Marković et al., 2016; Sprafke and Obreht, 2016 and references therein). However, most studies agree that the formation of dust (i.e., silt-sized particles) essentially involves a series of sedimentary processes prior to aeolian transport and loess formation (Smalley and Smalley, 1983; Pécsi, 1990; Pye, 1995; Muhs and Bettis, 2003; Badura et al., 2013; Sprafke and Obreht, 2016). The knowledge of the sedimentary pathways, thus the transport chains from the initial production of silt-sized particles to its final deposition and loess formation, is crucial to understand past environmental and climatic systems in varying geographical settings (Muhs, 2018). The reconstruction of atmospheric circulation patterns and climatic gradients is crucial to understand past environmental constraints on the spread of early modern humans (Obreht et al., 2017; Staubwasser et al., 2018), especially in areas which are believed to be important trajectories for their migration into Europe (Anghelinu and Niţă, 2014; Chu, 2018; Fitzsimmons et al., 2020).
The steps involved in dust production can be generalised and summarised in the pathways of desert and glacial loess genesis (Tsoar and Pye, 1987; Wright, 2001; Muhs, 2013; Muhs et al., 2014; Lancaster, 2020 and references therein). The main differences of these two pathways are the processes of dust formation, thus the production of silt-sized particles. Silt-sized particles of desert loess are mainly generated by physical weathering in high altitude areas and desert basins. This mode of genesis is constrained to the margins of drylands in low‐ to mid‐latitudes as well as extreme continental regions (Wright, 2001; Lancaster, 2020). For glacial loess, however, grinding by mountain glaciers and continental ice sheets (CIS) as well as frost shattering of periglacial mountainous areas are the main dust formation processes (Smalley and Smalley, 1983; Smalley et al., 2014). Furthermore, various authors highlighted the general significance of rivers as agents of reworking and storage of silt-sized particles, respectively (Smalley and Leach, 1978; Smalley et al., 2009; Badura et al., 2013; Lehmkuhl et al., 2018a, 2018b). The distinction between dust formation by CIS and processes in the high mountain areas was recently postulated by Li et al. (2020), who differentiated between the continental glacier-river transport (CR) mode and the mountain provenance-river transport (MR) mode of loess genesis.
The two glacial modes of loess genesis can be assumed as the main sedimentary pathways of European loess (Figure 1; Li et al., 2020). Loess deposits in northern Central Europe, roughly extending over southern England, France, Germany, Poland, Belarus, and Ukraine, are believed to partially originate from glacially ground material derived from the Fennoscandian and British-Irish ice sheets (Antoine et al., 2009; Skurzyński et al., 2019, 2020; Stevens et al., 2020). In Central and Southeastern Europe, silt particles forming loess deposits were subjected to a multi-step transport, essentially involving river transport of silt particles generated by intensive physical weathering and/or glacial activity in the adjacent mountain ranges (Smalley and Leach, 1978; Buggle et al., 2008; Smalley et al., 2009; Újvári et al., 2012).
FIGURE 1. (A) Distribution of aeolian sediments (modified according to Haase et al. (2007). Continental glacier provenance-river transport (CR) and mountain provenance-river transport (MR) modes of loess genesis (modified according to Li et al., 2020) for Europe are displayed. Danube (Da.), Dniester (Dn.) and Dnieper (Dp.) rivers are shown for orientation. Extent of the dry continental shelfs during the last glacial maximum (LGM) is modified according to Willmes (2015). (B) Location of the Lower Danube Basin in southeast Europe. Black triangles show investigated and reference loess-palaeosol sequences (LPS): 1: Urluia; 2: Rasova; 3: Vlasca; 4: Sageata; 5: Balta Alba Kurgan; 6: Langenlois; 7: Ságvár; 8: Bodrogkeresztúr; 9: Orlovat; 10: Petrovaradin; 11: Irig; 12: Kurortne. Elevations >300 m above sea level (a.s.l.) are shown as the upper topographic boundary for loess and archaeological open-air findings in Southeast Europe (Hauck et al., 2018). (C): Detailed view of the Lower Danube Basin with investigated LPS.
However, global models of loess genesis simplify transport processes and atmospheric circulation patterns, which hampers their (supra-) regional application. Regions can be influenced by several dust sources, which correspond to both modes of loess genesis. In Eastern Europe for example, the boundary between the two modes is roughly parallel to the Eastern Carpathians (Smalley et al., 2009; Li et al., 2020; Pańczyk et al., 2020; Skurzyński et al., 2020). Hence, the Lower Danube Basin (LDB) is located near the boundary between these two modes and thus it represents a transitional depositional zone, where a mixing of several potential source areas occurs (Buggle et al., 2008).
Therefore, the aim of this study is to evaluate potential source areas for the loess of the LDB and to refine our understanding of the modes of loess genesis in the region. We hypothesise five possible regional sedimentary pathways that may account for loess formation in the LDB and asses their importance alongside results of high-resolution geochemical analyses of five loess-palaeosol sequences (LPS) distributed throughout the basin. First, geochemical compositions are compared in order to evaluate signatures of each loess-palaeosol sequence that may reflect commonalities and differences in terminal dust sources. The continuous sampling in equidistant intervals allows us to reliably assess shifts in sediment provenance over time. Second, the elemental compositions and respective ratios are compared to loess records from reference regions upstream the Danube in Lower Austria, Central and Northeastern Hungary, and the Vojvodina region in Northern Serbia as well as loess deposits from the shore of the Black Sea in Ukraine. After establishing commonalities and differences of the different loess regions, the results are compared with findings from published provenance studies from the reference regions.
Regional Setting and Study Sites
The Lower Danube Basin
The LDB is the last depositional segment of the Danube River catchment, spanning between the Iron Gates and the Danube Delta (Figures 1, 2). The Lower Danube drainage basin sensu stricto covers large parts of the Transylvanian Basin, most of the Romanian Carpathians as well as the north facing slopes of the Balkan Mountains (Figure 2). The LDB sensuJipa (2014), who defined it as a loess sedimentary basin, only covers the lowlands of the Danube Plain (Romanian/Wallachian Plain and Bulgarian Plain, hereafter combined as Danube Plain) as well as the plateau of the Dobrogea. The basin is bounded by the Carpathians to the north and west, the Balkans to the west and south, and opens toward the Black Sea in the east. The Moldavian Plateau adjoins in the northeast (Figure 2). The basin is mainly covered by detrital Late Neogene to Pleistocene/Holocene deposits of predominantly marine, fluvial, and aeolian origin (Krézsek and Olariu, 2021). The aeolian sediments consist of loess, sandy loess, and Pleistocene to Holocene sand dunes, particularly in the western part of the basin, and essentially loess and loess derivates in the central and eastern part (Lehmkuhl et al., 2021).
FIGURE 2. (A) Perspective view of the LDB, looking from the southeast. Extents are indicated in Figure 1. Shaded areas indicate landscape units within the LDB: I: Wallachian Plain; II: Bulgarian Plain; III: Dobrogea; IV: Moldavian Plateau; V: Subcarpathians; VI: Pryedbalkans; VII: Transylvanian Basin; VIII: Carpathians; IX: Balkans; X: Danube Delta. Hypothesised sedimentary pathways are shown as black arrows (DA: Danube Alluvium pathway; CB: Carpathian Bending pathway; BS: Black Sea pathway; GF: Glaciofluvial pathway). Note that the scale bar is only accurate in the front part of the figure, due to the perspective view (B) Stratigraphic sketches for the investigated sections with correlation scheme.
The loess deposits in the LDB are among the thickest in Europe and laterally continuous over large areas (Haase et al., 2007; Fitzsimmons et al., 2012; Jipa, 2014; Marković et al., 2015). The landscape of the basin is diverse and allows for its division into seven subunits. The Carpathian foothills are intersected by a multitude of Danube tributaries such as Jiu, Olt, Argeş, Ialomiţa, and Buzău (Figure 1C). The Bulgarian Plain is characterized by a hilly landscape with short, steep tributaries discharging into the Danube that incises into mainly Mesozoic sediments underlying the Late Neogene to Quaternary mantles of the plain (Figure 2; Krézsek and Olariu, 2021). The Wallachian Plain is dominated by a flat topography, gently inclined toward the south. The left-side tributaries of the Danube, draining the Carpathians and their foothills, are sediment laden and show braided river characteristics. Therefore, they intersect the widespread Late Neogene to Quaternary deposits in broad, alluvial valleys, leaving vertical bluffs which often expose thick loess sequences. The topography of the Dobrogea plateau shows many ephemeral valleys, draining the area towards the Danube and the Black Sea. The Moldavian Plateau, which adjoins to the northeast, is dissected by the broad valleys of the Siret and the Prut, large tributaries of the Danube which are draining the Eastern Carpathian flysch. Additionally, the plateau is deeply dissected by a dense dendritic network of tributary systems.
The investigated loess sections are located along a south-north transect of the eastern part of the LDB in different geomorphological settings (Figure 2).
1) The Urluia (URL) LPS is located in an abandoned limestone quarry on the limestone plateau of the Dobrogea (Fitzsimmons et al., 2013; Fitzsimmons and Hambach, 2014; Obreht et al., 2017). The outcrop preserves five loess-palaeosol couplets reaching to the Middle Pleistocene. In this study, we focus on the uppermost 8.5 m, spanning from Marine Isotope Stage (MIS) 3 to 1. The Campanian Ignimbrite (CI) tephra, which is an important isochron in the area (Fitzsimmons et al., 2013; Veres et al., 2013), can be found in a stratigraphic depth of approximately 6 m at the studied profile. At the sampling site, tephra particles are syn-depositionally mixed into the loess and are exposed in an approx. 10 cm thick characteristic bed (Obreht et al., 2017). The tephra layer is dated in its source region of the Phlegrean Fields in Italy to an age of roughly 39–40 ka (De Vivo et al., 2001; Giaccio et al., 2017). The URL LPS was additionally dated by luminescence dating. The results were incorporated into a correlative age model (Obreht et al., 2017). This age model dates the uppermost 8.5 m of the section to a time span reaching back to ca. 55 ka.
2) The Rasova (RAS) sequence is exposed in a small valley on the right bank of the Dunarea Veche branch of the Danube bordering the Dobrogea to the west. The section is 7.2 m thick and overlies weathered Cretaceous limestone (Zeeden et al., 2018). The CI tephra is found at a depth of around 6 m and was confirmed by geochemical evidence and by luminescence dating (Anechitei-Deacu et al., 2014). The tephra is preserved as a thick bed of pure volcanic ash, with glass shards also incorporated into the overlying loess. Furthermore, the chronology of RAS relies on luminescence dating and a correlative age model (Zeeden et al., 2018) linked to the NGRIP δ18O record (North Greenland Ice Core Project members, 2004). The resulting age model dates the base of the sequence to ca. 47 ka.
3) The Vlasca (VLA) LPS is exposed at the left bank of the Braţul Borcea branch of the Danube River bordering the Wallachian Plain toward the east (Obreht et al., 2017). The loess bluff has a length of more than 2 km and a height of approx. 30–40 m above the current Danube level. The investigated part of the outcrop, covering MIS 3 to 1, is 12.5 m thick and the CI tephra is situated at a depth of 9.1 m. The tephra is preserved as a lateral continuous bed, several centimetres in thickness. The section was dated by a correlative approach to URL and Greenland ice core data (Obreht et al., 2017), which led to an age model with a maximum age of ca. 44 ka in a depth of 10.72 m. Furthermore, the underlying 1.8 m of the here investigated section were correlated to the published proxy data from URL covering the entire MIS 3 (Supplementary Material 3).
4) The section Sageata (SAG) is situated on the left bank of the Buzău River and is 5.5 m thick. The LPS consists of sandy loess sediments and overlies river pond loams, which developed in the floodplain of the Buzău River. A more than 50 cm thick layer of pure ash, assigned to the CI tephra, is intercalated in the sequence. The section was not directly dated yet but correlated to the other LPS records using the CI tephra and the Holocene topsoil as tie points (Supplementary Material 3).
5) The Balta Alba Kurgan (BAK) site is situated 1.5 km west of the alkaline Balta Alba lake in the northern Danube Plain. The section is exposed in a road cut opening the edge of a loess plateau. Within the upper 2 m of the section, remains of a Neolithic to Bronze Age tumulus can be found. Below the tumulus, the sequence consists of loess with several intercalated interstadial palaeosols and weakly pedogenically overprinted horizons. The CI tephra lies at a depth of 8.06 m and is preserved as a horizontal bed of patchy tephra lenses. The investigated section is 10.9 m thick. The chronological framework for this section is based on an integrative age model, using luminescence dating, radiocarbon dating as well as (palaeo-) magnetic correlation (Scheidt et al., 2021). The resulting model shows a maximum age of ca. 53 ka at a depth of 9.8 m.
All investigated LPS were correlated using proxy data and their respective age models revealing a set of records covering the timeframe from MIS 3 to MIS 1.
Possible Sedimentary Pathways of the Lower Danube Loess
The loess deposits of the LDB are located in the transitional zone between the East European loess zone and the Caspian-Black Sea loess zone (Stephens et al., 2002). The provenance of the Danube alluvium (DA) in the Lower Danube Basin is diverse, including sediments derived from the upper reaches of the river in Central Europe but also more locally derived material carried by its main tributaries within the basin. The DA is regarded as the main source area for loess deposits within the LDB (Buggle et al., 2008; Jipa, 2014). The Carpathian Bending (CB), narrowly defined here as the flysch area, drained by the Buzău and Ialomiţa rivers consists mainly of Cretaceous and Neogene flysch deposits, as well as unconsolidated but uplifted Pliocene to Quaternary marine to brackish sediments evolving into lacustrine deposits and locally interfingered with fluvial and alluvial fan deposits (e.g, Krézsek and Olariu, 2021 and references therein). The sediment yields of these rivers are among the largest in Europe and their carrying capacities are strongly reduced when entering the basin. Their river flow regime changes from an intra-mountainous, hence erosive, to a depositional fluvial system. These rivers built up large alluvial fan systems and even nowadays show braided river characteristics (Râdoane et al., 2003; Floroiu, 2011). These sedimentary cascades are also influenced by complicated local tectonics (Matenco et al., 2016; Krézsek and Olariu, 2021). The strong tectonic subsidence in the Lower Danube area, that comprises several km thick Quaternary marine to brackish, lacustrine, and fluvial deposits (Matenco and Andriessen, 2013; Matenco et al., 2016; Krézsek and Olariu, 2021) coupled with uplifting in the Carpathian Bending area (Necea et al., 2013; Krézsek and Olariu, 2021) considerably influence the riverine morphocharacteristics and sediment yield. Thus, the exposed large valley plains act as deflation areas within and in the foreland of the Carpathian Bending area (Jipa, 2014). The third pathway, which similarly to the DA and CB pathways corresponding to the MR mode (Li et al., 2020), is the drainage of the Eastern Carpathians (EC). This mode is mainly fed by eroded Cretaceous and Neogene flysch sediments, as well as widespread unconsolidated fluvial, aeolian, and shallow marine deposits of Late Miocene to Quaternary age covering the Eastern European Platform. These areas are drained by the Danube River system with its tributaries Siret and Prut, as well as the Dniester River that discharges into the Black Sea. This sediment pathway was identified as an important dust source for loess deposits in Poland and (southwestern) Ukraine (Nawrocki et al., 2018; Pańczyk et al., 2020). The glaciofluvial outwash plains of Central and Northern Ukraine pose a major dust source for Ukrainian loess-palaeosol sequences (Nawrocki et al., 2006, 2019; Buggle et al., 2008; Pańczyk et al., 2020). This pathway, hereafter called glaciofluvial (GF) pathway, can be related to the CR mode (Li et al., 2020). Another potential sediment pathway is the dust entrainment from exposed sediments on the shelf of the Black Sea (BS) during low sea level stands. This pathway cannot be directly linked to any loess genesis mode sensuLi et al. (2020), since the catchment of the Black Sea between southern Ukraine and Bulgaria covers riverine input from several mode domains, including the DA. For example, the sediment yield of the Danube corresponds to the MR mode, whereas the sediment yield of the Bug, Dnieper, and Don is derived from the CR mode. Since the BS pathway is assumed to be an important dust source area for loess deposits in Bulgaria (Avramov et al., 2006; Jordanova et al., 2007), it is considered as a potential pathway for the LDB loess.
Sampling and Laboratory Analyses
Methodological Background of Loess Provenance Studies
The range of methodological approaches used in loess provenance research is broad. Conceptual research assessed dust sources based on the geographical and geomorphological situation of loess deposits (Smalley and Leach, 1978; Leger, 1990; Smalley et al., 2009; Badura et al., 2013; Lehmkuhl et al., 2016, 2018b) or prevailing atmospheric regimes favoring dust transport and deposition (Antoine et al., 2001, 2009; Vandenberghe et al., 2006). Spatial variability in grain size distributions can be used as a proxy for the distance from the source area (Frazee et al., 1970; Jipa, 2014; Schaetzl et al., 2018). More complex statistical analyses of grain size data, e.g. based on end member models, have been used to unmix grain size distributions and to constrain transport processes, pathways, as well as potential source areas (Prins and Vriend, 2007; Bokhorst et al., 2011; Nottebaum et al., 2015). The imbrication of grains in sediments, recorded by the anisotropy of magnetic susceptibility (AMS) allows the reconstruction of palaeowind directions and thus dust transport pathways (Lagroix and Banerjee, 2002; Nawrocki et al., 2006, 2018; Bradák, 2009; Bradák et al., 2020; Zeeden and Hambach, 2021). Heavy mineral assemblage data can be used to fingerprint loess deposits and compare them to potential source sediments. However, using heavy mineral distribution assemblages as a provenance proxy may bear the drawback of mixing multiple distinct sources (Mikulčić Pavlaković et al., 2011; Römer et al., 2016; Meng et al., 2019; Wolf et al., 2019). Recently, the advent of single-grain provenance techniques such as detrital zircon U-Pb dating allowed disentangling multiple dust sources and potentially sign specific proto-source terranes (Stevens et al., 2010; Újvári et al., 2012; Fenn et al., 2018; Pańczyk et al., 2020). However, the aeolian system is subjected to well-known gravitative distance sorting processes, potentially leading to an over- or underrepresentation of proximal sources, depending on the sample location (Schaetzl et al., 2018).
Loess geochemical provenance proxies circumvent misinterpretations of proximal and distal dust sources by including sedimentary components representative for all grain size classes and mineralogical components (Buggle et al., 2008; Skurzyński et al., 2020). Despite the aforementioned mixing of sources, distinct geochemical differences between loess sequences provide information about the detrital background of the sediment and therefore allow the interpretation of dust sources (Muhs, 2018). Usually, the geochemical compositions of loess deposits are compared to those of potential source areas to determine geochemical commonalities (Buggle et al., 2008; Obreht et al., 2016). Despite methodological limitations (Buggle et al., 2008; Schatz et al., 2015a), geochemical signatures provide a tool which is neither time-consuming nor expensive. Its application has enabled for example to distinguish dust sources of the Chinese Loess Plateau and loess deposits in southern China (Hao et al., 2010) and to exclude local sources for aeolian silt mantles on the Californian Channel Islands (Muhs et al., 2007, 2008). Campodonico et al. (2019) confirmed the dominance of Andean detritus as source material for Pampean loess deposits of Central Argentina, based on bulk geochemical evidence.
Due to its complex formation, loess represents a geological and lithological mixture of its source areas. The different mineralogical components of loess deposits are prone to a multitude of methodological issues and sensible to, e.g., instrumentation or pretreatment of samples (Goff et al., 2020). Differences between different geochemical analysis techniques, such as inductively coupled plasma mass spectrometry (ICP-MS) and x-ray fluorescence analyses (XRF) can account to up to 10% of the measured contents (Miyazaki et al., 2016). To exclude such biases (Demir et al., 2006), we only use XRF data from samples pretreated and measured at the Physical Geography Laboratory of the RWTH Aachen University (Germany).
Sampling of all sections was conducted in a continuous sampling line. Before sampling, several decimetres of sediment were removed from the vertical face of the loess wall to avoid contamination from freshly weathered exposed loess or downslope relocated material. Samples were taken contiguously in 2 cm increments for URL, RAS, VLA, and BAK and in 10 cm increments for SAG. The reference sections discussed in our regional comparison were sampled in 3–10 cm resolution (Table with sampling resolution and number of samples is provided in Supplementary Table S1). Sampling was conducted with cleaned tools and samples were packed into sterile plastic bags.
All samples were analysed by means of an energy dispersive x-ray fluorescence (EDPXRF, hereafter XRF) analysis using a Spectro Xepos. The samples were sieved to the silt fraction (<63 µm) and dried at 105 °C for 12 h. A subsample of 8 g for each sample was mixed with 2 g Fluxana Cereos wax and homogenised in a shaker. The admixture was pressed to a pellet with a pressure of 19.2 MPa for 120 s. The measurements were conducted by means of a pre-calibrated method. Each sample was measured in duplicate and the pellets were rotated by 90° between measurements to avoid matrix effects. The results of every sample were averaged. Conspicuous samples with striking differences between both measurements were remeasured in duplicate to avoid analytical artifacts. Major elements are given in oxide form (%; Na2O, Al2O3, SiO2, K2O, TiO2, Fe2O3), whereas trace elements are expressed in mg/kg (V, Cr, Ni, Rb, Sr, Y, Zr, Ba, Th). A similar analytical procedure was applied for all other studies used in this comparison (Obreht et al., 2015; Bösken et al., 2018, 2019; Marković et al., 2020; Tecsa et al., 2020).
The results were visualised in bi- and depth plots to investigate differences and similarities in space and time. To visualise the relations between soluble and insoluble elements and therefore interpret post-depositional alterations of sediments due to weathering, the Al2O3-CaO*+Na2O-K2O (molar proportions) ternary (A–CN–K) diagram (Nesbitt and Young, 1984) was plotted for the investigated LPS to display commonalities and differences between the sites. CaO* reflects the content of Ca in silicates. The ratio of the total Ca content (CaO) and Na2O in silicates is assumed to never exceed a value of 1. Therefore, when CaO content is higher than the Na2O, the siliceous Ca content (CaO*) is set to the content of Na2O (McLennan, 1993; Buggle et al., 2011).
Grain Size Analyses
For grain size analyses, the samples were dried at 35 °C, sieved to the fine earth fraction (<2 mm), and two subsets of each sample (between 0.1 and 0.3 g) were pretreated with 0.7 ml H2O2 (30%) at 70 °C for 12 h. This process was repeated until a bleaching of the material was visible (Allen and Thornley, 2004), but not longer than three days. To keep the particles dispersed during the analysis, samples were treated with 1.25 ml Na4P2O7 (0.1 mol*L−1) in an overhead shaker for 12 h. The grain size was determined with a Beckman Coulter LS 13320 laser diffractometer using the Mie-theory (Fluid RI: 1.33, Sample RI: 1.55, Imaginary RI: 0.1; Özer et al., 2010; Nottebaum et al., 2015; Schulte et al., 2016). The results of grain size analyses are visualised as heatmaps. This visualisation allows interpretations of shifts in grain size distributions with depth (Schulte et al., 2016; Schulte and Lehmkuhl, 2018). From each loess sequence, an interval of 30 cm apparently unaffected by post-depositional processes such as weathering or pedogenesis was chosen. These intervals were determined using the grain size distributions and magnetic susceptibility. Over each of these intervals, the average of the mean grain size was calculated as representative for undisturbed loess of the respective section.
Environmental Magnetic Analyses and Correlative Age Model
Samples for environmental magnetic analyses were dried, ground using a porcelain mortar, and packed into plastic boxes. The samples were fixed with cotton wool to avoid movement of particles during transport and measurement. The volumetric magnetic susceptibility was measured at two frequencies (300 and 3,000 Hz) in a static field of 300 mA/m, using a Magnon International VSFM at Bayreuth University (Germany). The results were corrected for drift and the weak diamagnetism of the sample boxes and normalised for density. Environmental magnetic data was used for correlation of the investigated LPS to create a tentative regional age model (Supplementary Material 3, Figure S1).
A correlative age model (Supplementary Material 3, Figure S1) was established based on published age models for the investigated LPS URL and VLA (Obreht et al., 2017), RAS (Zeeden et al., 2018), and BAK (Scheidt et al., 2021). SAG was correlated based on tephrochronological evidence and trends in magnetic susceptibility. The LPS were correlated on the MIS scale spanning from MIS 3 to 1. To ensure comparability to the reference regions, only (parts of) sections covering the same or comparable time frames were used (Marković et al., 2007, 2020; Obreht et al., 2015; Bösken et al., 2018, 2019; Tecsa et al., 2020; Pötter et al., 2021).
Within the geochemical data set, the tephra layers show the largest fluctuations in their geochemical composition compared to the surrounding sediments, especially for RAS and SAG (see below). However, to ensure the full integrity of the presented records, the data of the samples affected by the tephra are presented together with other samples, but not included into the discussion. The same applies for samples where volcanic glass shards are incorporated into the sediment in considerable amounts, influencing the geochemical signature noticeably. In the data visualisations, extreme outliers caused by tephra particles are displayed together and additionally marked.
The geochemical compositions of the investigated LPS exhibit broad variations (Supplementary Table S2). Mean values for the Na2O content vary between 0.92% for URL and 1.34% in SAG. The latter shows the highest values within the CI tephra layer (max.: 4.84%). Mean Al2O3 contents spread between 10.28% for URL and 12.01% for SAG. BAK shows the highest variability of SiO2 contents, spanning from 39.84 to 66.72%, with an average of 50.41%. The other sections, with exception of SAG (mean: 52.78%) exhibit lower contents (URL: 43.46%; RAS: 45.61%; VLA: 48.73%). The mean content of K2O of the southern sites of URL (1.68%), RAS (1.75%), and VLA (1.74%) are generally lower than those of the northern sites of BAK (1.89%) and SAG (2.71%). The high values in SAG again occur in the thick CI tephra bed. TiO2 values for SAG again show the broadest spread (min.: 0.48%; max.: 0.84%), whereas the mean value (0.68%) is the lowest of the investigated LPS (URL: 0.78%; RAS: 0.83%; VLA: 0.79%; BAK: 0.74%). The range of mean Fe2O3 contents of 0.16% for all sections is narrow (VLA: 4.05%; RAS: 4.21%). The trace element compositions as well as the data for the reference section can be found in Supplementary Table S2.
The bi-plot visualisations of element contents and ratios show distinct patterns when comparing the investigated LPS. The ratios of major elements (Figures 3, 4) overlap between the sites. This is especially the case in the Al2O3/K2O ratio (Figure 3). In this ratio as well as the Al2O3/SiO2 ratio, URL and RAS show very similar trends. Especially URL is characterised by the lowest values for both elements, and the point clouds of both sites are roughly parallel to the diagonal of the plot. VLA shows higher SiO2 contents compared to URL and RAS, especially in samples with lower Al2O3 values. The overlap of these three sites is generally quite large, whereas BAK and SAG are characterised by higher SiO2, K2O, and Al2O3 contents. The TiO2 values are depleted for these two sites. For SAG, the CI tephra layer is visible in all major elements with either extremely elevated or depleted values. The TiO2/Ni plot, adapted from Újvári et al. (2008), shows different trends within the five sites. BAK, SAG, and URL show elongated point clouds with a trend parallel to a positive correlation between Ti and Ni. VLA, however, shows a more circular point cloud with a trend roughly perpendicular to BAK and URL. The pattern of RAS shows both, a circular and a linear component in its pattern.
FIGURE 3. Bi-plots of major element composition for all investigated sections and reference regions. Symbols represent the here investigated loess-palaeosol sequences (Urluia (URL): green; Rasova (RAS): blue; Vlasca (VLA): black; Sageata (SAG): gray; Balta Alba Kurgan (BAK): red. Extreme outliers from the CI tephra layers and mixed loess-tephra samples are indicated by the letter T. The dashed lines show the geochemical composition of the reference regions: Lower Austria = orange; Lake Balaton (Central Hungary) = purple; Bükk Mountains (Northeast Hungary) = blue; Vojvodina = green; Black Sea Shore (Ukraine) = yellow (for location of the reference sections, see Figure 1). Note that extreme outliers for the reference regions are excluded.
FIGURE 4. Bi-plots of major element ratios for all investigated sections and reference regions. Legend according to Figure 3. Note that extreme outliers for the reference regions are excluded.
The Na2O/Al2O3 vs. K2O/Al2O3 diagram shows a large overlap between URL, RAS, VLA, and most of the BAK record (Figure 4). The rest of BAK as well as the majority of SAG are shifted toward higher K2O/Al2O3 and lower Na2O/Al2O3 ratios. In RAS and SAG, the CI tephra layer is clearly visible in strongly elevated values. The distinction between the investigated LPS is more evident in the other major element ratios. Especially the SiO2/TiO2 vs. K2O/Al2O3 diagram shows a clear differentiation between most sites. RAS and URL again are overlapping, whereas VLA is characterised by higher SiO2/TiO2 ratios. BAK shows a narrow range, with similar K2O/Al2O3 and elevated SiO2/TiO2 ratios. The range of the SAG data, even when the tephra samples are not included, is quite high for both ratios. Whereas the SiO2/TiO2 is comparable to BAK, with higher maximum values, the K2O/Al2O3 ratio shows a spread of 0.17 to more than 0.25. The Al2O3/TiO2 and to lesser extent the Fe2O3/TiO2 ratio show a similar pattern.
Trace element compositions vary strongly between the five LPS (Figure 5). Whereas URL and VLA overlap in large parts, RAS shows elevated Zr/Ni and Zr/V ratios. The samples from the tephra layer in SAG have, similarly to RAS, extremely high values in these ratios. The loess layers are characterised by depleted values. The Cr/V vs. Y/Ni diagram shows a trend subparallel to the x-axis. Here, SAG exhibits the lowest Y/Ni values, but the CI tephra samples again show very high values for VLA, BAK, and SAG compared to the other LPS. With the samples covering the CI tephra excluded, the Cr/V ratios for these three LPS vary between 0.5 and 1.5. RAS and URL show a large spread, from values around 1 to more than 4. The Ba/Sr vs. Rb/Sr diagram depicts an almost linear trend between these two ratios for most of the sample. BAK is slightly shifted toward lower Rb/Sr ratios. SAG, which shows the largest spread in both ratios, also has an almost horizontal trend.
FIGURE 5. Bi-plots of trace element composition for all investigated sections and reference regions. Legend according to Figure 3. Note that extreme outliers for the reference regions are excluded.
The A–CN–K diagram shows point clouds subparallel to the CN joint for all sections (Figure 6). Whereas this is the case for all samples from URL, VLA, and BAK, RAS and SAG have a strong vertical component in their pattern. The point clouds of the respective sections have relative offsets to each other. VLA shows the largest spread on the Al2O3 joint, whereas RAS is shifted toward the CN joint, and BAK and especially SAG are shifted toward the K2O apex. The point clouds of VLA and BAK have very similar starting points trending toward the composition of the upper continental crust (UCC).
FIGURE 6. A–CN–K plot for the investigated sections (URL: Urluia; RAS: Rasova; VLA: Vlasca; SAG: Sageata, BAK: Balta Alba Kurgan). Compositions of minerals as well as grain size, weathering and sorting trends are modified according to Buggle et al. (2008). UCC = Upper Continental Crust. Extreme outliers from the CI tephra layers and mixed loess-tephra samples are indicated by the letter T.
Grain Size Analyses
The results of the granulometric analyses for the investigated LPS are depicted in Figure 7. All five sites are dominated by strong modes of coarse silt (20–63 µm). However, the five sites show distinguished patterns of grain size distributions. The base of the investigated section at URL is characterised by a distinct mode of coarse silt, with a slightly fining tendency up-section (Figure 7A). Between 7 and 6 m below surface (b.s.), the contents of finer particles, particularly fine silt, increase. Around the CI tephra layer (approx. 6 m b.s.), a thin but distinct peak in finer particles occurs. The mode of coarse silt is less well expressed but still visible. Above the tephra, the sediment becomes coarser until 5 m depth and the mode of coarse silt as well as fine sand contents increase. Between 4 and 3 m b.s., two coarser layers with increased amounts of sand occur. Above 3 m, a decreasing trend in grain size can be seen toward the top. Generally, the grain size distributions are well sorted.
FIGURE 7. Grain size distributions for the investigated sections, sorted from south to north (A: Urluia; B: Rasova; C: Vlasca; D: Sageata, E: Balta Alba Kurgan). F: Grain size distributions for unweathered loess from the investigated sections [URL: Urluia (3–3.3 m); RAS: Rasova (2–2.3 m); VLA: Vlasca (5.1–5.4 m); SAG: Sageata (2–2.3 m); BAK: Balta Alba Kurgan (3.5–3.8 m)].
RAS generally shows a coarser distribution than URL, as the main mode is between coarse silt and fine sand (Figure 7B). The basal part is characterised by a strong mode of these two grain size classes. The CI tephra is clearly visible as a strong distinct peak of coarse grain sizes, with a strong mode in fine sand and increased middle sand contents. Above the tephra, the strong mode continues until a depth of around 5.5 m. Within the 50 cm above, two distinct coarse peaks are visible. Above that, the grain size distribution is well sorted, with a fining trend toward the top. Despite this trend, no major variations are visible.
The mode of coarse silt is also visible in VLA, where the basal part shows high percentages of this class (Figure 7C). Between 11.5 m b.s and the CI tephra at a depth of 9.1 m, increased proportions of medium and fine silt occur. The coarse silt mode dominates with increased grain size between 7 and 6 m b.s.. Above that, the distributions remain rather consistent until a depth of around 2 m, where the fine silt and especially the clay contents increase toward the top of the profile. The loess packages with a prominent mode in coarse silt are very well sorted.
The grain size distribution of SAG is very complex (Figure 7D). The base of the section, which overlies alluvial loams, shows no clear sorting. Within the thick tephra bed, a prominent mode within the fine sand fraction is visible. This continues above the ash layer, until a depth of around 3.1 m. Between three and 2.5 m b.s., the sediment fines up and shows a mode between coarse silt and fine sand. Above, the sediment is coarser and less well sorted. Between 1.7 and 1 m, the sediment is relatively well sorted and the mode ranges in the coarse silt fraction. Above 1 m b.s., middle clay contents increase.
The mode of the grain size distribution of BAK is from the base to approx. 7 m b.s rather constant, ranging in the coarse silt grain size (Figure 7E). Within this package, however, there are four layers of increased clay content as well as a coarser interval. Above that, a slight coarsening of the sediments texture is noticeable, with a more pronounced coarse silt mode. The sediment is generally very well sorted. Until 4 m b.s., four coarser layers are intercalated. Between 4 and 3 m, a coarsening trend is visible, whereas between 3 and 2 m, the clay content increases. The uppermost 2 m of the profile exhibit a mode of coarse silt to fine sand, with varying clay contents. However, since this is the anthropogenic infill of the kurgan, it is not included in the provenance discussion.
The mean grain size distributions for unaltered loess intervals are shown in Figure 7F. The distributions of URL and BAK show almost identical unimodal courses, with a prominent peak in the coarse silt fraction. Both sites also exhibits shoulders in coarse clay, as well as fine sand. The peak in VLA is shifted slightly toward coarser grain sizes, whereas the clay and fine silt contents are lower than in URL and BAK. The mode for RAS lies on the boundary between coarse silt and fine sand, and the content of fine sand is higher than at the other sites. The most complex distribution is found in SAG, where the mode is less expressed and not as narrow as for the other sites. The content of fine sand is considerably higher than, e.g., in BAK. SAG also shows the highest contents of clay, especially in the coarse and middle clay fraction.
The Geochemical Imprint of the Campanian Ignimbrite Tephra
Tephra layers alter the chemical composition of the sediment in which they are incorporated (e.g., Obreht et al., 2016). The CI tephra, which is used here as a chronostratigraphic marker horizon (Constantin et al., 2012; Fitzsimmons et al., 2013; Veres et al., 2013; Obreht et al., 2017), leaves a clear geochemical signature compared to the embedding loess. Due to its trachytic characteristics the CI tephra can be identified using geochemical proxies of bulk sediment samples even when present as cryptotephra (see Obreht et al., 2016). The influence of the CI tephra on the geochemical composition of the embedding loess varies strongly between the investigated sections. In BAK and URL, where the CI tephra is only preserved in patchy lenses or mixed with loess, respectively, the geochemical composition of the embedding loess samples seems less affected (e.g., Figures 8, 9). Contrary, for VLA the Zr/V and Zr/Ni ratio for example are increased since the ash layer is preserved as a laterally continuous stratum of several cm thickness (Figure 10, Supplementary Figure S3). The strongest and most prominent peaks in Zr/V can be seen in RAS (30 cm thick tephra) and SAG, where the tephra bed consists of pure pyroclastic material. Especially in SAG, where the CI tephra is more than 50 cm thick, parts of the layer consist of pure ash beds, which strongly influence the geochemical signal of the analysed samples. Since volcanic glass shards are usually incorporated into subsequently deposited dust (Obreht et al., 2017), the geochemical effects of the tephra can be found even some cm to dm above the actual layer forming a gradual transition to the overlying loess. Geochemical proxy data, especially those based on diagnostic elements, can be easily used to detect volcanic material in loess sequences. In the here investigated LPS, the tephra layers show typical characteristics imprinted in the geochemical data. The Al2O3/TiO2 ratio in RAS, e.g., shows a rapid increase at the lower boundary of the tephra and a gentle decline toward the above-lying loess, indicating incorporation of volcanic glass shards into the sediment. In other localities, such as BAK where the tephra is preserved as a patchy layer, the geochemical influences up-profile of the tephra are reduced (Figure 10). However, these large variations originating from the volcanic ash, both in samples of pure ash beds as well as varying amounts of glass shards in the sediment, cannot be interpreted as true loess provenance signals but must be kept in mind when discussing transport pathways of contemporaneous dust deposition.
FIGURE 8. Depth variations of the SiO2/TiO2 ratios for the investigated loess-palaeosol sequences, sorted from south to north (URL: Urluia; RAS: Rasova; VLA: Vlasca; SAG: Sageata, BAK: Balta Alba Kurgan). The box-whisker plots on top show the composition of the sections from the reference regions: I: Lower Austria; II: Vojvodina; III: Balaton (Central Hungary); IV: Bükk Mountains (Northeast Hungary) and V: Black Sea Shore (Ukraine). The interquartile area (between 25 and 75% quartile) are shown as shaded boxes. The correlative age model as well as the CI tephra are shown for temporal comparison.
FIGURE 9. Depth variations of the K2O/Al2O3 ratios for the investigated loess-palaeosol sequences, sorted from south to north (URL: Urluia; RAS: Rasova; VLA: Vlasca; SAG: Sageata, BAK: Balta Alba Kurgan). The box-whisker plots on top show the composition of the sections from the reference regions: I: Lower Austria; II: Vojvodina; III: Balaton (Central Hungary); IV: Bükk Mountains (Northeast Hungary) and V: Black Sea Shore (Ukraine). The interquartile area (between 25 and 75% quartile) are shown as shaded boxes. The correlative age model as well as the CI tephra are shown for temporal comparison.
FIGURE 10. Depth variations of the Zr/V ratios for the investigated loess-palaeosol sequences, sorted from south to north (URL: Urluia; RAS: Rasova; VLA: Vlasca; SAG: Sageata, BAK: Balta Alba Kurgan). The box-whisker plots on top show the composition of the sections from the reference regions: I: Lower Austria; II: Vojvodina; III: Balaton (Central Hungary); IV: Bükk Mountains (Northeast Hungary) and V: Black Sea Shore (Ukraine). The interquartile area (between 25 and 75% quartile) are shown as shaded boxes. The correlative age model as well as the CI tephra are shown for temporal comparison.
Geochemical Signature and Mineralogical Composition of the Lower Danube Loess
Some mineralogical approximations can be made based on loess geochemical data (Muhs, 2018). Generally high SiO2 contents show the predominance of quartz in the mineralogical composition of all five LPS. The varying contents, both between and within the LPS, are mainly governed by differences in the mineralogical compositions and grain size distributions (Figure 8, Supplementary Figure S2). The clear positive relationship between Al2O3 and SiO2 (Figure 3) indicates the presence of aluminosilicates such as mica or feldspars in all LPS investigated here. This is supported by the positive relation between Al2O3 and K2O, which suggests that the aluminum content is mainly bound to phyllosilicates or K-feldspars. The presence of Ti-bearing phyllosilicates is also supported by the weak positive correlation between Ti and Ni (Újvári et al., 2008), since Ni tends to be enriched in these minerals. This is, to varying degrees, visible in all sections (Figure 3).
The A–CN–K ternary diagram developed by Nesbitt and Young (1984) has been widely applied to LPS from varying environments to investigate influences of chemical weathering and mineralogical compositions on the geochemical composition of sediments (e.g., Yang et al., 2004; Schatz et al., 2015b; Vlaminck et al., 2018; Skurzyński et al., 2019). All investigated LPS show a trend of plagioclase weathering (Figure 6). This trend indicates a rather stable ratio of K-feldspar and albite for most of the samples. The (sub-) parallel point clouds of especially URL is characteristic for loess deposits of the Dobrogea, such as the Mircea Voda (Buggle et al., 2008) or the Costineşti LPS (Tugulan et al., 2016). The shifts of BAK and especially SAG toward the right edge of the plot indicates either increased feldspar weathering or different mineralogical composition of source areas compared to the other sites (Vlaminck et al., 2018). In the latter case, it can be interpreted to reflect more K-feldspar rich parent material (Újvári et al., 2008). The abundance of K-bearing minerals, such as K-feldspars, muscovite, or illite is indicated for all sections by the shift toward the K-apex in the A–CN–K ternary diagram (Figure 6) and a strong correlation between the Rb/Sr and Ba/Sr ratio (Figure 5; Újvári et al., 2008). High Ba/Sr ratios also show the presence of plagioclase, as indicated by the A–CN–K plot. However, the shifts of SAG and BAK toward higher Rb/Sr and Ba/Sr values, as well as their extreme values, can be interpreted as Sr depletion due to leaching processes during soil formation (Újvári et al., 2008).
The geochemical and mineralogical compositions of loess and palaeosols is not solely determined by the composition of the parent material. In subaerial deposits such as loess, post-depositional processes including weathering or pedogenesis frequently alter the material’s geochemical composition (Vlaminck et al., 2018; Bösken et al., 2019; Skurzyński et al., 2020). As seen in the A–CN–K plot (Figure 6), the majority of the samples show a discernible plagioclase weathering trend. Especially soluble or easily weathered compounds and minerals are prone to depletion by leaching. The Ba/Sr and the Rb/Sr ratios, e.g., are used in loess research as weathering indices, since Sr is easily leached during carbonate dissolution while Ba and Rb are rather immobile elements, which are relatively enriched during leaching (Chen et al., 1999; Buggle et al., 2011). When discussing loess provenance based on geochemistry, palaeosol samples are often excluded to overcome issues of geochemical alteration during pedogenesis (Buggle et al., 2008; Újvári et al., 2008; Schatz et al., 2015a). However, a clear distinction between loess, pedogenically overprinted loess and fully developed palaeosols is not always straightforward in semi-arid environments. Therefore, we included all samples and plotted provenance proxies as depth functions to complement the bi-plot visualisations. In these plots (Figures 8–13), data from the reference regions are displayed as box-whisker plots on top of the graphs. The interquartile ranges are indicated as shaded boxes of the respective color in the depth plots. The effects of solution and weathering can be seen in a variety of element contents and ratios. The Na2O/Al2O3 ratio is depleted in the recent topsoil at URL, VLA, and BAK, and increased and/or constant in the loess intervals (Supplementary Figure S4). Post-depositional alterations by weathering or pedogenesis are also expressed by the (neo-)formation of (aluminous) clay minerals or iron oxides. As the former may lead to a relative enhancement of aluminum oxides (Figure 11, Supplementary Figure S2), the latter leads to an increase in iron contents in the affected sediment (Figure 12). Since the formation if these minerals is usually linked to pedogenesis, the elevated Al and Fe contents can be related to the clay content (Supplementary Figure S2). The Al2O3/TiO2 ratio plotted as a depth function (Figure 11) as well as the direct comparison between the Al content and the grain size (Supplementary Figure S2) reflects this issue and must therefore be carefully discussed. Shifts in these ratios, such as those observed below the CI tephra in VLA, cannot be interpreted as provenance signal as they are of post-depositional nature.
FIGURE 11. Depth variations of the Al2O3/TiO2 ratios for the investigated loess-palaeosol sequences, sorted from south to north (URL: Urluia; RAS: Rasova; VLA: Vlasca; SAG: Sageata, BAK: Balta Alba Kurgan). The box-whisker plots on top show the composition of the sections from the reference regions: I: Lower Austria; II: Vojvodina; III: Balaton (Central Hungary); IV: Bükk Mountains (Northeast Hungary) and V: Black Sea Shore (Ukraine). The interquartile area (between 25 and 75% quartile) are shown as shaded boxes. The correlative age model as well as the CI tephra are shown for temporal comparison.
FIGURE 12. Depth variations of the Fe2O3/TiO2 ratios for the investigated loess-palaeosol sequences, sorted from south to north (URL: Urluia; RAS: Rasova; VLA: Vlasca; SAG: Sageata, BAK: Balta Alba Kurgan). The box-whisker plots on top show the composition of the sections from the reference regions: I: Lower Austria; II: Vojvodina; III: Balaton (Central Hungary); IV: Bükk Mountains (Northeast Hungary) and V: Black Sea Shore (Ukraine). The interquartile area (between 25 and 75% quartile) are shown as shaded boxes. The correlative age model as well as the CI tephra are shown for temporal comparison.
Since trace elements such as Ni and V are compatible elements within igneous systems, they can be used as indicators for the provenance of rocks and sediments (Újvári et al., 2014; Obreht et al., 2016). Ratios based on these trace elements reveal striking similarities between URL and VLA (Figures 3, 5). RAS shows large overlaps with the two sites in, e.g., the Th/Ni or Th/V ratios (Figure 5). However, the Zr based provenance ratios vary substantially between those sites. These differences can be explained by the enrichment of zircon minerals in RAS; in detrital sediments Zr is almost entirely bound to zircons. These heavy minerals can be enriched in loess by selective deflation of less dense particles. This so-called hydraulic sorting (Malusà et al., 2016) causes enrichment of zircons which in turn leads to elevated Zr contents in the geochemical screening of such samples. Therefore, the elevated Zr/Ni and Zr/V ratios can be interpreted as an indicator of sorting processes (Figure 5; Újvári et al., 2008).
The hydraulic sorting of minerals by water or air not only affects the trace element geochemistry, such as the Zr content of loess deposits, but is also expressed in other geochemical proxies (Ohta, 2004). The position and the shape of the point cloud of a sample set in the A–CN–K plot can be indicative for such sorting processes. Due to the abundance of Al in aluminous clay minerals, finer, often weathered sediments tend to plot close to the Al2O3 apex of the diagram, whereas coarser sediments show a contrary trend (Nesbitt and Young, 1996; Nesbitt et al., 1996). When the grain size differences within a LPS are solely driven by (plagioclase) weathering, point clouds plot (sub-) parallel to the CaO*–Na2O join (see Vlaminck et al., 2018; Bösken et al., 2019; Skurzyński et al., 2019). However, when hydraulic sorting skews the grain size distribution, point clouds in the A–CN–K plot tend to develop (sub-)vertical shapes, trending toward the Al2O3 apex (Ohta, 2004; Buggle et al., 2008). For BAK and URL, the point clouds sub-parallel to the CaO*–Na2O join indicate no major influence of hydraulic sorting processes. VLA shows a similar pattern, but the point cloud shows a slightly stronger vertical component (Figure 6). The strongest evidence for sorting processes can be found in RAS and SAG, with strong vertical components. The direct vicinity of RAS and SAG to the recent courses of the Danube and the Buzău might be the main reason for these sorting processes. The high intensity of these processes in RAS can be explained by enhanced fluvial activity of small Danube tributaries acting as dust sources, as this section is exposed in the recent dry valley of such a tributary and its base shows a sharp fluvial erosional contact (Zeeden et al., 2018). The pre-aeolian setting is therefore fluvial as it is in SAG. This section lies on top of alluvial loams, which were presumably deposited in river ponds in the floodplain of the palaeo-Buzău. Therefore, the general high fluvial activity strongly influenced the formation of these loess deposits during the Late Pleistocene. The alluvial deposits near the sites of RAS and SAG, and to a lesser extent of VLA, acted as local sources for coarse material deposited in direct vicinity while finer particles were transported over longer distances. The grain size distributions also reflect these processes, since the three sites show generally coarser distributions than URL and BAK (Figure 7).
Generally, the investigated LPS show striking similarities in some of the presented provenance proxies. Especially the sites from the Dobrogea, URL and RAS, show significant common variabilities, indicating geochemically similar parent material as expressed in the Cr/V ratios (Figures 5, 13). This ratio serves as a good indicator for the contributions of detrital material derived from mafic to ultramafic rocks since Cr is almost entirely bound to mafic strata (Zimmermann and Bahlburg, 2003; Mikes et al., 2006; Újvári et al., 2008). For the other sites, the values of this ratio are in a rather narrow spectrum without large variations. However, data from URL and RAS exhibit large scatter of this ratio with maximum values of >4. This feature is in accordance with elevated Cr contents compared to the UCC established in loess sequences from the southern Dobrogea (Tugulan et al., 2016), which indicates the influence of distinct sources for the loess deposits of the area (Figure 5, Supplementary Figure S7).
FIGURE 13. Depth variations of the Cr/V ratios for the investigated loess-palaeosol sequences, sorted from south to north (URL: Urluia; RAS: Rasova; VLA: Vlasca; SAG: Sageata, BAK: Balta Alba Kurgan). The box-whisker plots on top show the composition of the sections from the reference regions: I: Lower Austria; II: Vojvodina; III: Balaton (Central Hungary); IV: Bükk Mountains (Northeast Hungary) and V: Black Sea Shore (Ukraine). The interquartile area (between 25 and 75% quartile) are shown as shaded boxes. The correlative age model as well as the CI tephra are shown for temporal comparison.
VLA shows similarities to URL in several of the presented geochemical provenance proxies. The correspondence with URL is evident for example in the Na2O/Al2O3 vs. K2O/Al2O3 and the Al2O3/TiO2 vs. SiO2/TiO2 plot (Figure 4), suggesting similar sourced material contributed to both records. RAS, however, shows some overlap with URL, but also exhibits a distinct pattern, e.g., in the Al2O3/TiO2 ratio, which is depleted compared to the other sites. Besides the general commonalities identified within the studied LPS, significant variations are also discernible. It is noteworthy that the variations in the K2O/Al2O3 vs. SiO2/TiO2 diagram tend to plot closer to the data from BAK, with which VLA also has large overlaps. SAG, however, shows some commonalities to BAK, but is distinguishable from the other sites, indicating a substantial contribution from a different dust source area. Therefore, we propose similar dominant source areas for the URL, RAS, and VLA sequences, with a more diverse input from multiple sources for the latter. For SAG, different and especially local sources, such as the proximal input from the Carpathian Bending can be assumed while BAK represents a more mixed provenance signal of different sources. While the use of multiple insightful bulk geochemical proxies has enabled us to sufficiently differentiate loess deposits in the Lower Danube Basin, further reference data is needed to assign the respective dust sources.
Estimations of Sedimentary Pathways of the Lower Danube Loess
As stated previously, at this stage of research we compare our results with geochemical data obtained using the same instrumentation and pretreatment from other loess deposits in Central and Southeastern Europe (i.e., the Middle Danube Basin (MDB) and Ukraine), where previous multi-proxy provenance studies have constrained potential dust sources (Supplementary Table S1). This comparison potentially allows to indirectly assign potential sources to the loess deposits in the Lower Danube Basin, a region yet largely underrepresented in loess provenance studies. The discussed provenance studies were chosen according to their proximity to the reference sections from Austria, Hungary, the Vojvodina, and Ukraine. However, local and regional differences in deflation and accumulation cannot be excluded. Therefore, only a tentative identification of potential source areas can be postulated in addition to the discussed geochemical characterisation. An overview of the discussed and additional provenance (Supplementary Table S3) and palaeowind direction studies (Supplementary Table S4) can be found in the Supplementary Material.
The Danube Alluvium (DA) Sedimentary Pathway
The large overlap of especially URL, RAS, and VLA with the geochemical compositions of some of the reference sites in Lower Austria, Central Hungary, and the Vojvodina, indicates similar geochemical signatures of the dust sources to the respective loess deposits (Figures 8, 9, Supplementary Figure S4). The alluvial sediments of the MDB are recognized as the main source for the loess deposits in the area (Újvári et al., 2012; Obreht et al., 2019). This is in accordance with findings from Buggle et al. (2008), who demonstrated geochemical similarities between the Mircea Voda site in the LDB and the Batajnica/Stari Slankamen composite sequence of the Vojvodina. The authors interpreted these similarities as an indicator for the importance of the Danube alluvium (DA) as the main sediment supplier and dust source in the LDB, corresponding to the MR mode of loess genesis (Li et al., 2020). However, the geochemical composition of the alluvial sediment changes laterally, as major tributaries in the MDB contribute to the sediment load of the Danube (Nenadić et al., 2016; Vandenberghe et al., 2018; Mitrinović et al., 2020). The influence of other local to sub-regional source areas for the loess of the Middle Danube Basin is evident in the range in the geochemical data of the reference sections from the MDB. The importance of this influence increases with distance from the course of the Danube (Basarin et al., 2011; Obreht et al., 2014, 2016, 2019). Accordingly, the loess deposits of Northeast Hungary, which cover the upper to middle reaches of the Tisa catchment (Schatz et al., 2015a; Bösken et al., 2019), show different geochemical signatures as expressed in the K2O/Al2O3 vs. SiO2/TiO2 diagram (Figure 4). These patterns are also evident for other LPS of the area, which show considerable geochemical differences to deposits in direct proximity to the Danube (Schatz et al., 2015a). Enhanced contributions from other local sources, such as the Bohemian Massif for Austrian loess (Újvári et al., 2013; Újvári and Klötzli, 2015), exposed local sandstones in Hungary (Újvári et al., 2008, 2012; Sebe, 2013) or from the Deliblatska Peščara sand field in Serbia (Obreht et al., 2015) were detected using various methods. However, for most loess deposits in the Carpathian Basin, the Danube alluvium appears to be the main dust source. The large overlap between the geochemical compositions of especially VLA, RAS, and URL with sections from the Middle Danube Basin indicates that the Danube alluvium in the LDB area is geochemically similar to the upper reaches of the river and acted as the dominant regional source area for dust deflation and ultimately, loess formation. Geochemical differences are either due to local sources within the LDB, but can also be attributed to palaeoclimatic gradients (Buggle et al., 2008).
The Eastern Carpathian (EC) Sedimentary Pathway
The SiO2/TiO2 ratio also shows a good agreement between BAK and SAG and the reference section of Kurortne (Black Sea Shore) in southwest Ukraine (Tecsa et al., 2020), which is visible in the K2O/Al2O3 vs. SiO2/TiO2 diagram as well (Figures 4, 8). In general, two sediment pathways were reconstructed for Ukrainian loess covers. On the one hand, the loess deposits of Central Ukraine are linked to the glaciofluvial (GF) sedimentary pathway (Buggle et al., 2008). These deposits, inter alia located in the catchment of the Dnieper River (Veres et al., 2018), consist of aeolian dust entrained from the glaciofluvial outwash plains by strong northerly winds (Buggle et al., 2008; Nawrocki et al., 2019). On the other hand, the loess deposits of southern Ukraine rather correspond to the Eastern Carpathian (EC) sedimentary pathway sourced from the sediments of the Dniester River and its tributaries. Northwesterly palaeowinds have been reconstructed for dust transport in this region (Nawrocki et al., 2006, 2018; Pańczyk et al., 2020). Although the GF pathway cannot be ruled out completely for Kurortne, its location on the southern Ukrainian shore of the Black Sea (Figure 1; Tecsa et al., 2020) close to the Dniester estuary suggests the EC pathway to account for loess deposits in this area. Therefore, the striking overlap with especially BAK suggests similar source areas for both sections and agrees with the reconstructed strong and consistent N-NW palaeowinds in the area. This is supported by residual forms of loess dunes, so-called gredas (Rozycki, 1967; Leger, 1990; Mason et al., 1999), which can be found in the entire eastern foreland of the Carpathians, and which are usually oriented NW-SE (Figure 14; Rozycki, 1967; cited in Leger, 1990). Although detrital material derived from the Carpathians has a strong influence on the sediment yield of the Danube, the sediments of tributaries draining the (Outer) Eastern Carpathians can be distinguished by geochemistry (Buggle et al., 2008) and detrital zircon age assemblages (Ducea et al., 2018). Therefore, the differences between BAK and loess deposits further south in the LDB such as VLA and URL can be linked to the dominance of the EC sedimentary pathway for the BAK area.
FIGURE 14. A: potential sedimentary pathways for loess deposits in the Lower Danube Basin (DA: Danube alluvium; CB: Carpathian Bending; BS: Black Sea; EC: Eastern Carpathians; GF: glaciofluvial). Modes of loess genesis according to Li et al. (2020). Residual forms of gredas and longitudinal dunes in the Middle and Lower Danube Basin are adapted from (Rozycki 1967; cited in Leger, 1990). The division of the Carpathians is given to demonstrate different (proto-) source areas such as the flysch terrains of the Outer Eastern Carpathians, the volcanic chains of the Inner Eastern Carpathians or the crystalline and metamorphic bed rocks of the Southern Carpathians. The study sites are shown as triangles. B: surface trends for zircons (pink) and hornblende minerals (green). The isolines represent surface concentrations (in %) and are derived from statistical computation of mineralogical analyses of hundreds of surface samples (Andar and Codarcea, 1981; cited in; Jipa, 2014). The arrows indicate reconstructed transport directions (pink: towards Danube; green: towards Carpathians), based on the heavy mineralogical evidence.
The Carpathian Bending (CB) Sedimentary Pathway
Another potential mode of origin for loess at SAG and BAK is the Carpathian Bending (CB) sedimentary pathway, with primary material being eroded from the bending sector of the Eastern Carpathians and transported by rivers like Buzău or Ialomiţa. This pathway corresponds to the MR mode of loess genesis (Li et al., 2020). There are no direct provenance studies on loess deposits from this region, but the floodplains of the rivers draining the Southern Carpathians to the west of the CB were detected as potential source areas based on heavy mineral composition data (Jipa, 2014). Decreasing zircon concentrations toward the Danube plain can be linked to the decreasing influence of sediment input from the Carpathians where detrital material is enriched in zircons. The opposing trend in hornblende contents was interpreted as the transport direction of dust deflated from the Danube alluvium, as the hornblende minerals are enriched in the deposits of the Danube (Figure 14B; Andar and Codarcea, 1981; cited in Jipa, 2014). Figure 14B shows that BAK and SAG are located in a transitional area, where the influence of the Eastern Carpathians is strong, but also some amount of hornblende can be found within the loess deposits, pointing to the Danube alluvium as a potential source area. The EC pathway is supported by high zirconium concentrations for both sites but especially in SAG (Supplementary Table S2). As stated above, the inferred zircon enrichment is indeed an indicator for sorting processes at this site. However, the fact that according to the A–CN–K plot (Figure 6) not all samples of the SAG section are affected by sorting, the generally high values could be a sign for the input of detrital material from the Carpathians for these loess deposits. The overlap with data from northeast Hungary in proxies such as the Al2O3/TiO2 ratio (Figures 4, 11) or the elevated Zr contents (Figure 5 and Supplementary Table S2) support the influence of detrital material from the Carpathians, as northeast Hungary is strongly controlled by the Tisa and its tributaries (Thamó-Bozsó et al., 2014; Schatz et al., 2015a). Therefore, these areas experience a strong input from the Inner Eastern Carpathians and their extensive volcaniclastic forelands, which are also drained by the upper reaches of some main tributaries of the Danube in the LDB such as the Olt River.
The Black Sea (BS) and Glaciofluvial (GF) Sedimentary Pathways
The main characteristics of the BS and GF sedimentary pathways cannot be illustrated in detail due to the lack of reference data. However, the importance of these pathways has been discussed for loess deposits of the LDB and adjacent areas. Especially the BS pathway was linked to the formation of loess deposits in Bulgaria (Avramov et al., 2006; Jordanova et al., 2007). The complexity of the sedimentological composition of the Black Sea shelf, caused by multiple large rivers discharging into the basin, hinders a clear link to this pathway based on geochemistry. Large amounts of the sediments building the Black Sea shelf were, e.g., deposited by the palaeo-Danube, indicating similar geochemical and mineralogical compositions to the DA pathway (Ducea et al., 2018). The GF pathway was reconstructed for LPS in Poland and Ukraine (Buggle et al., 2008; Skurzyński et al., 2019, 2020; Pańczyk et al., 2020) and was mainly driven by enhanced aeolian activities of katabatic winds influenced by fluctuations of the Eurasian CIS (Nawrocki et al., 2019). Pańczyk et al. (2020) constrained a divide of dust provinces based on detrital zircon ages, which divides the GF from the EC mode in Poland and Ukraine. In combination with the reconstruction of atmospheric circulations for the last glacial (Meyer and Kottmeier, 1989; Obreht et al., 2017), a major input of entrained dust from the GF pathway seems unlikely, especially for the southern part of the LDB.
Spatial and Temporal Provenance Variations of Lower Danube Loess and Palaeoclimatic Implications
The provenance of Lower Danube loess deposits varies spatially. While the sites from the Dobrogea plateau and the Danube plain seem to be linked mainly to sediment input via the Danube alluvium, sites in the northern part of the LDB are influenced by a mixture of potential source areas, including the DA, the EC, and the CB sedimentary pathways.
The hypothesized sedimentary pathways for the Lower Danube loess highlight the importance of proximal dust sources, as previously stated for the Middle Danube Basin (Obreht et al., 2019). However, far traveled dust originating from, e.g., the Sahara Desert can contribute especially to the fine silt fraction of aeolian deposits (Stuut et al., 2009). Although the importance of these distal sources for Southeastern European dust records was discussed for interglacial periods (Longman et al., 2017), the proposed contribution to glacial loess is of minor importance (see Varga et al., 2013).
Changes of dust source areas over time are potentially linked to palaeoclimatic fluctuations driving dust source expansion and/or changes in dust transport pathways (Zhang et al., 2018; Nawrocki et al., 2019; Pańczyk et al., 2020; Skurzyński et al., 2020). For MIS 3, URL, VLA, and BAK show some variation in the zirconium-based provenance proxies (Figure 10, Supplementary Figure S3). Especially in VLA, the Zr/Ni and Zr/V ratios show a decreasing trend from the bottom of the investigated section towards the CI tephra, whereas above it the values remain rather constant. The same trend but to lesser extent is visible in BAK. Since the grain size and especially the median grain size do not show large variations during this period (Figure 7), zircon enrichment due to sorting seems unlikely. The decreasing ratios are rather an indicator for the weakening of the influence of the Eastern Carpathians as a source area during this period. The weakened input of material derived from the Eastern Carpathians can be explained by shifts of atmospheric patterns during this time. During the middle MIS 3, westerlies reached Southeastern Europe on a NW-SE trajectory (see also Staubwasser et al., 2018). During the late MIS 3, the Siberian High pressure system intensified, leading to a shift of westerlies southwards (Obreht et al., 2017). Local maxima in the Zr based ratios in VLA suggest temporally increased input of Carpathian detrital material, indicating short-term fluctuations in source material. In URL, the Cr/V record is clearly divided during MIS 3 (Figure 13). Below the CI tephra, some oscillations occur indicating the input of detrital material potentially originating from more mafic sources. Above the tephra, these oscillations decrease and are within the interquartile range of the reference regions upstream. In RAS, the Cr/V proxy shows large variations throughout the whole profile, with exceptions during mid-MIS 2 and MIS 1, indicating an almost constant influence of mafic sources for these intervals. The pattern during MIS 2 is similar to URL, but the values are generally higher. Accumulation effects can explain the higher values and the large variations through time. Cr rich material was eroded in the Dobrogea and transported via small rivers toward the Danube Plain. The sediment was accumulated near RAS, deflated and deposited as aeolian dust (Zeeden et al., 2018).
The early MIS 2 was generally characterized by stable source areas for the investigated LPS. URL, RAS and VLA are clearly dominated by the DA sedimentary pathway during this period, whereas BAK shows large overlaps with the EC pathway (e.g., Figure 8). In the middle of MIS 2, however, the geochemical evidence shows changes in all sequences, e.g., in the K2O/Al2O3 ratio (Figure 9). URL, RAS, and BAK show a noticeable decline in this proxy, whereas VLA and SAG show increasing values. These changes could be explained by the shift toward prevailing easterlies in large parts of Europe during the LGM (Schaffernicht et al., 2020). This would imply that the Black Sea (BS) sedimentary pathway, or long-range transport of material deflated even further to the east, becomes progressively more important throughout MIS 2, possibly contributing dust into the aeolian system of the LDB more significantly than hitherto considered. This coincides with sea level low stands of the Black Sea during the LGM (Constantinescu et al., 2015). However, since no reference data from sediments of the Black Sea shelf is available for this study, this hypothesis remains tentative. The same applies for long-range dust transport from the Caspian Basin further east, which was postulated for southwestern Ukraine especially for glacial periods (Stephens et al., 2002).
Methodological Limitations and Potential for Further Research
The general limitations of this methodological approach, such as the averaging of input from multiple dust sources, must be considered when discussing geochemical source estimations. The admixture of material from different proto-sources in a large sedimentary system such as the Danube Basin, can complicate the interpretation of distinct dust sources. When using bulk sediment (<63 µm) geochemistry as a provenance tool, signatures of multiple sources may be averaged out and assigning distinct sources might not be feasible (Bird et al., 2015). To identify specific dust (proto-)sources, it is essential that future research complements this method with further analyses, preferably single grain techniques that allow for more conclusive discrimination between potential source sediments of similar geochemical composition (Stevens et al., 2010; Fenn et al., 2018). Nevertheless, the geochemical composition of bulk sediments can be used to either rule out sources or to assess the main variability in loess provenance based on the varying contents of diagnostic elements and identifiable geochemical similarities or differences (Buggle et al., 2008; Újvári et al., 2008), as presented here for a transect of records from the Lower Danube loess.
Summary and Conclusions
Our high-resolution geochemical study on several loess-palaeosol sequences from the Lower Danube Basin is among the first multi-record analyses contributing to the understanding of loess provenance and major sedimentary pathways in different sectors of the basin. Based on the presented data, we subdivide the basin into two major domains of loess geochemical signatures. The southern part of the basin including the southern Wallachian Plain and the Dobrogea with the LPS of Urluia, Rasova, and Vlasca can be distinguished from the northern Wallachian Plain with the LPS of Sageata and Balta Alba Kurgan. The boundary between these two domains appears to be transitional rather than a sharp contact.
The comparison with reference sections from Austria, Hungary, the Vojvodina, and Ukraine reveals that the Urluia, Rasova, and Vlasca records show striking geochemical similarities with loess records from the Middle Danube Basin. Since these similarities can be seen in many provenance proxies, a similar dominant geochemical signature of source areas in the Middle and Lower Danube Basin can be assumed. Furthermore, the geochemical composition of the Balta Alba Kurgan LPS rather resembles the one of Kurortne in Ukraine, pointing to similar source areas for these regions. Lastly, Sageata, situated in the floodplain of the Buzău River, shows rather unique geochemical patterns compared to the other sections, despite some overlaps with Balta Alba Kurgan.
The link of these geochemical findings with published provenance studies from the reference regions allows a tentative identification of potential source areas for the Lower Danube loess. It confirms the previous perception that the alluvial sediments of the Danube and its major tributaries act as a main storage and as source areas for dust deflation in the basin. The comparison to loess deposits in the middle and upper reaches of the Danube supports the Danube alluvium pathway hypothesis, especially for the southern sites of Urluia, Rasova, and Vlasca. Further to the north, other sources gain more importance. For Sageata, detrital material eroded from the Carpathian Bending area and redeposited by sediment-laden rivers such as the proximal Buzău River plays a major role as a local dust source. For Balta Alba Kurgan, the Eastern Carpathian sedimentary pathway is of particular importance. Here, material is deflated from the riverbeds draining the northern slopes of the Carpathians. Therefore, the dominant mode of loess genesis here is the mountain provenance-river transport loess sensuLi et al. (2020). The input of material deflated from the glaciofluvial outwash plains of Ukraine can neither be confirmed nor ruled out. However, we consider this sedimentary pathway, corresponding to the continental glacier provenance-river transport mode sensuLi et al. (2020) to be of minor importance, especially for the southern Lower Danube Basin. Our study confirms the general assumptions of prevailing northwesterly to westerly winds during the Pleniglacial in the area. These findings should be complemented using single grain techniques such as detrital zircon U-Pb ages to illuminate the sedimentary pathways accounting for dust deposition in the Lower Danube Basin in greater detail.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
SP, PS, and UH designed the study. SP compiled the figures and wrote the first draft of the manuscript. SP, DV, PS, and UH defined the methodology. PS supervised the laboratory work and generated the discussed data. SP, DV, JN, UH, PS, and FL among others conducted fieldwork. All authors discussed the data, visualisations, and text.
The investigations were carried out in the frame of the CRC 806 “Our way to Europe,” subproject B1 “The Eastern Trajectory”: “Last Glacial Palaeogeography and Archaeology of the Eastern Mediterranean and of the Balkan Peninsula,” funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation)—Projektnummer 57444011-SFB 806.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We thank Marianne Dohms and her team for the laboratory framework. Alina Blume and Bruno Boemke helped vectorising parts of Figures 1, 14. We thank Dagmar Haase for providing the shapefiles of loess distribution for Figure 1. We thank Jacek Skurzyński for his input on the methodological discussion. DV’ work on the tephrochronology of investigated sequences was supported through a fellowship provided by Alexander von Humboldt Foundation. We thank Annett Junginger for the editorial work and MQ and RS for their helpful comments, which improved the manuscript substantially.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/feart.2021.600010/full#supplementary-material.
Andar, P., and Codarcea, V. (1981). Analiza statisticaă a variabilităaţt¸ ii mineralelor grele din loessurile şs¸ i depozitele loessoide din Dobrogea şs¸ i Câmpia Româna Occidentală (Statistical analysis of the variability of heavy minerals in the loess and loess deposits from Dobrogea and the Western Romanian Plain). Dari Seama Institutul Geol. Geofiz. LXVI, 15–41.
Anechitei-Deacu, V., Timar-Gabor, A., Fitzsimmons, K., Veres, D., and Hambach, U. (2014). Multi-method luminescence investigations on quartz grains of different sizes extracted from a loess section in Southeast Romania interbedding the Campanian Ignimbrite ash layer. Geochronometria 41, 1–14. doi:10.2478/s13386-013-0143-4
Antoine, P., Rousseau, D.-D., Moine, O., Kunesch, S., Hatté, C., Lang, A., et al. (2009). Rapid and cyclic aeolian deposition during the Last Glacial in European loess: a high-resolution record from Nussloch, Germany. Quat. Sci. Rev. 28, 2955–2973. doi:10.1016/j.quascirev.2009.08.001
Antoine, P., Rousseau, D.-D., Zöller, L., Lang, A., Munaut, A.-V., Hatté, C., et al. (2001). High-resolution record of the last Interglacial–glacial cycle in the Nussloch loess–palaeosol sequences, Upper Rhine Area, Germany. Quat. Int. 76–77, 211–229. doi:10.1016/S1040-6182(00)00104-X
Avramov, V. I., Jordanova, D., Hoffmann, V., and Roesler, W. (2006). The role of dust source area and pedogenesis in three loess-palaeosol sections from North Bulgaria: a mineral magnetic study. Stud. Geophys. Geod. 50, 259–282. doi:10.1007/s11200-006-0015-y
Basarin, B., Vandenberghe, D. A. G., Marković, S. B., Catto, N., Hambach, U., Vasiliniuc, S., et al. (2011). The Belotinac section (Southern Serbia) at the southern limit of the European loess belt: initial results. Quat. Int. 240, 128–138. doi:10.1016/j.quaint.2011.02.022
Bird, A., Stevens, T., Rittner, M., Vermeesch, P., Carter, A., Andò, S., et al. (2015). Quaternary dust source variation across the Chinese Loess Plateau. Palaeogeogr. Palaeoclimatol. Palaeoecol. 435, 254–264. doi:10.1016/j.palaeo.2015.06.024
Bokhorst, M. P., Vandenberghe, J., Sümegi, P., Łanczont, M., Gerasimenko, N. P., Matviishina, Z. N., et al. (2011). Atmospheric circulation patterns in central and eastern Europe during the Weichselian Pleniglacial inferred from loess grain-size records. Quat. Int. 234, 62–74. doi:10.1016/j.quaint.2010.07.018
Bösken, J., Obreht, I., Zeeden, C., Klasen, N., Hambach, U., Sümegi, P., et al. (2019). High-resolution paleoclimatic proxy data from the MIS3/2 transition recorded in northeastern Hungarian loess. Quat. Int. 502, 95–107. doi:10.1016/j.quaint.2017.12.008
Bösken, J., Sümegi, P., Zeeden, C., Klasen, N., Gulyás, S., and Lehmkuhl, F. (2018). Investigating the last glacial Gravettian site ‘Ságvár Lyukas Hill’ (Hungary) and its paleoenvironmental and geochronological context using a multi-proxy approach. Palaeogeogr. Palaeoclimatol. Palaeoecol. 509, 77–90. doi:10.1016/j.palaeo.2017.08.010
Bradák, B. (2009). Application of anisotropy of magnetic susceptibility (AMS) for the determination of paleo-wind directions and paleo-environment during the accumulation period of Bag Tephra, Hungary. Quat. Int. 198, 77–84. doi:10.1016/j.quaint.2007.11.005
Bradák, B., Seto, Y., Chadima, M., Kovács, J., Tanos, P., Újvári, G., et al. (2020). Magnetic fabric of loess and its significance in Pleistocene environment reconstructions. Earth-sci. Rev. 210, 103385. doi:10.1016/j.earscirev.2020.103385
Buggle, B., Glaser, B., Hambach, U., Gerasimenko, N., and Marković, S. (2011). An evaluation of geochemical weathering indices in loess–paleosol studies. Quat. Int. 240, 12–21. doi:10.1016/j.quaint.2010.07.019
Buggle, B., Glaser, B., Zöller, L., Hambach, U., Marković, S., Glaser, I., et al. (2008). Geochemical characterization and origin of Southeastern and Eastern European loesses (Serbia, Romania, Ukraine). Quat. Sci. Rev. 27, 1058–1075. doi:10.1016/j.quascirev.2008.01.018
Campodonico, V. A., Rouzaut, S., and Pasquini, A. I. (2019). Geochemistry of a Late Quaternary loess-paleosol sequence in central Argentina: implications for weathering, sedimentary recycling and provenance. Geoderma 351, 235–249. doi:10.1016/j.geoderma.2019.04.024
Chen, J., An, Z., and Head, J. (1999). Variation of Rb/Sr ratios in the loess-paleosol sequences of Central China during the last 130,000 Years and their implications for monsoon paleoclimatology. Quat. Res. 51, 215–219. doi:10.1006/qres.1999.2038
Chu, W. (2018). The Danube corridor hypothesis and the Carpathian Basin: geological, environmental and archaeological approaches to characterizing Aurignacian dynamics. J. World Prehistory 31, 117–178. doi:10.1007/s10963-018-9115-1
Constantin, D., Timar-Gabor, A., Veres, D., Begy, R., and Cosma, C. (2012). SAR-OSL dating of different grain-sized quartz from a sedimentary section in southern Romania interbedding the Campanian Ignimbrite/Y5 ash layer. Quat. Geochronol. 10, 81–86. doi:10.1016/j.quageo.2012.01.012
Constantinescu, A. M., Toucanne, S., Dennielou, B., Jorry, S. J., Mulder, T., and Lericolais, G. (2015). Evolution of the Danube deep-sea fan since the last glacial maximum: new insights into Black Sea water-level fluctuations. Mar. Geol. 367, 50–68. doi:10.1016/j.margeo.2015.05.007
De Vivo, B., Rolandi, G., Gans, P. B., Calvert, A., Bohrson, W. A., Spera, F. J., et al. (2001). New constraints on the pyroclastic eruptive history of the Campanian volcanic Plain (Italy). Mineral. Petrol. 73, 47–65. doi:10.1007/s007100170010
Demir, F., Budak, G., Baydaş, E., and Şahin, Y. (2006). Standard deviations of the error effects in preparing pellet samples for WDXRF spectroscopy. Nucl. Instrum. Methods Phys. Res. Sect. B Beam Interact. Mater. 243, 423–428. doi:10.1016/j.nimb.2005.09.019
Ducea, M. N., Giosan, L., Carter, A., Balica, C., Stoica, A. M., Roban, R. D., et al. (2018). U-PB detrital zircon geochronology of the lower Danube and its tributaries: implications for the geology of the Carpathians. Geochem. Geophys. Geosystems 19, 3208–3223. doi:10.1029/2018GC007659
Fenn, K., Stevens, T., Bird, A., Limonta, M., Rittner, M., Vermeesch, P., et al. (2018). Insights into the provenance of the Chinese Loess Plateau from joint zircon U-Pb and garnet geochemical analysis of last glacial loess. Quat. Res. 89, 645–659. doi:10.1017/qua.2017.86
Fitzsimmons, K. E., Doboş, A., Probst, M., and Iovita, R. (2020). Thinking outside the box at open-air archeological contexts: examples from loess landscapes in southeast Romania. Front. Earth Sci. 8, 561207. doi:10.3389/feart.2020.561207
Fitzsimmons, K. E., Hambach, U., Veres, D., and Iovita, R. (2013). The Campanian Ignimbrite eruption: new data on volcanic ash dispersal and its potential impact on human evolution. PLoS ONE 8, e65839. doi:10.1371/journal.pone.0065839
Fitzsimmons, K. E., Marković, S. B., and Hambach, U. (2012). Pleistocene environmental dynamics recorded in the loess of the middle and lower Danube basin. Quat. Sci. Rev. 41, 104–118. doi:10.1016/j.quascirev.2012.03.002
Giaccio, B., Hajdas, I., Isaia, R., Deino, A., and Nomade, S. (2017). High-precision 14C and 40Ar/39Ar dating of the Campanian Ignimbrite (Y-5) reconciles the time-scales of climatic-cultural processes at 40 ka. Sci. Rep. 7, 45940. doi:10.1038/srep45940
Goff, K., Schaetzl, R. J., Chakraborty, S., Weindorf, D. C., Kasmerchak, C., and Bettis, E. A. (2020). Impact of sample preparation methods for characterizing the geochemistry of soils and sediments by portable X‐ray fluorescence. Soil Sci. Soc. Am. J. 84, 131–143. doi:10.1002/saj2.20004
Haase, D., Fink, J., Haase, G., Ruske, R., Pécsi, M., Richter, H., et al. (2007). Loess in Europe—its spatial distribution based on a European Loess Map, scale 1:2,500,000. Quat. Sci. Rev. 26, 1301–1312. doi:10.1016/j.quascirev.2007.02.003
Hao, Q., Guo, Z., Qiao, Y., Xu, B., and Oldfield, F. (2010). Geochemical evidence for the provenance of middle Pleistocene loess deposits in southern China. Quat. Sci. Rev. 29, 3317–3326. doi:10.1016/j.quascirev.2010.08.004
Hauck, T. C., Lehmkuhl, F., Zeeden, C., Bosken, J., Thiemann, A., and Richter, J. (2018). The Aurignacian way of life: Contextualizing early modern human adaptation in the Carpathian Basin. Quat. Int. 485, 150–166. doi:10.1016/j.quaint.2017.10.020
Jordanova, D., Hus, J., and Geeraerts, R. (2007). Palaeoclimatic implications of the magnetic record from loess/palaeosol sequence Viatovo (NE Bulgaria): palaeoclimatic implications of the magnetic record. Geophys. J. Int. 171, 1036–1047. doi:10.1111/j.1365-246X.2007.03576.x
Krézsek, C., and Olariu, C. (2021). Filling of sedimentary basins and the birth of large rivers: the lower Danube network in the Dacian Basin, Romania. Glob. Planet. Change 197, 103391. doi:10.1016/j.gloplacha.2020.103391
Lehmkuhl, F., Bösken, J., Hošek, J., Sprafke, T., Marković, S. B., Obreht, I., et al. (2018a). Loess distribution and related quaternary sediments in the Carpathian Basin. J. Maps 14, 661–670. doi:10.1080/17445647.2018.1526720
Lehmkuhl, F., Nett, J. J., Pötter, S., Schulte, P., Sprafke, T., Jary, Z., et al. (2021). Loess landscapes of Europe - mapping, geomorphology, and zonal differentiation. Earth-sci. Rev. 215, 103496. doi:10.1016/j.earscirev.2020.103496
Lehmkuhl, F., Zens, J., Krauß, L., Schulte, P., and Kels, H. (2016). Loess-paleosol sequences at the northern European loess belt in Germany: distribution, geomorphology and stratigraphy. Quat. Sci. Rev. 153, 11–30. doi:10.1016/j.quascirev.2016.10.008
Longman, J., Veres, D., Ersek, V., Salzmann, U., Hubay, K., Bormann, M., et al. (2017). Periodic input of dust over the Eastern Carpathians during the Holocene linked with Saharan desertification and human impact. Clim. Past 13, 897–917. doi:10.5194/cp-13-897-2017
Marković, S. B., Oches, E. A., McCoy, W. D., Frechen, M., and Gaudenyi, T. (2007). Malacological and sedimentological evidence for “warm” glacial climate from the Irig loess sequence, Vojvodina, Serbia. Geochem. Geophys. Geosystems 8, Q09008. doi:10.1029/2006GC001565
Marković, S. B., Stevens, T., Kukla, G. J., Hambach, U., Fitzsimmons, K. E., Gibbard, P., et al. (2015). Danube loess stratigraphy — towards a pan-European loess stratigraphic model. Earth-sci. Rev. 148, 228–258. doi:10.1016/j.earscirev.2015.06.005
Marković, S. B., Vandenberghe, J., Stevens, T., Mihailović, D., Gavrilov, M. B., Radaković, M. G., et al. (2020). Geomorphological evolution of the Petrovaradin fortress paleolithic site (novi sad, Serbia). Quat. Res. 14, 33. 10.1017/qua.2020.88
Matenco, L., and Andriessen, P. (2013). Quantifying the mass transfer from mountain ranges to deposition in sedimentary basins: source to sink studies in the Danube Basin–Black Sea system. Glob. Planet. Change 103, 1–18. doi:10.1016/j.gloplacha.2013.01.003
Matenco, L., Munteanu, I., ter Borgh, M., Stanica, A., Tilita, M., Lericolais, G., et al. (2016). The interplay between tectonics, sediment dynamics and gateways evolution in the Danube system from the Pannonian Basin to the western Black Sea. Sci. Total Environ. 543, 807–827. doi:10.1016/j.scitotenv.2015.10.081
Meng, X., Liu, L., Zhao, W., He, T., Chen, J., and Ji, J. (2019). Distant Taklimakan Desert as an important source of aeolian deposits on the Chinese Loess Plateau as evidenced by carbonate minerals. Geophys. Res. Lett. 46, 4854–4862. doi:10.1029/2018GL081551
Meyer, H.-H., and Kottmeier, C. (1989). Die atmosphärische Zirkulation in Europa im Hochglazial der Weichsel-Eiszeit—abgeleitet von Paläowind-Indikatoren und Modellsimulationen. Eampg Quat. Sci. J. 39, 10–18. doi:10.3285/eg.39.1.02
Mikulčić Pavlaković, S., Crnjaković, M., Tibljaš, D., Šoufek, M., Wacha, L., Frechen, M., et al. (2011). Mineralogical and geochemical characteristics of Quaternary sediments from the Island of Susak (Northern Adriatic, Croatia). Quat. Int. 234, 32–49. doi:10.1016/j.quaint.2010.02.005
Mitrinović, D., Zarić, J., Anđelković, O., Sipos, G., Polomčić, D., and Dimkić, M. (2020). Dating of alluvial sediments from borehole at the lower course of the Sava river and indications of the connection between their genesis and climate changes in the Pleistocene. Quat. Int. 547, 75–85. doi:10.1016/j.quaint.2019.09.028
Miyazaki, T., Kimura, J.-I., and Katakuse, M. (2016). Geochemical records from loess deposits in Japan over the last 210 kyr: lithogenic source changes and paleoclimatic indications. Geochem. Geophys. Geosystems 17, 2745–2761. doi:10.1002/2016GC006322
Muhs, D. R., and Bettis, A. E. (2003). “Quaternary loess-Paleosol sequences as examples of climate-driven sedimentary extremes,” in Extreme depositional environments: mega end members in geologic time. Editors M. A. Chan, and A. W. Archer (Boulder: Geological Society of America). doi:10.1130/0-8137-2370-1.53
Muhs, D. R., Budahn, J., Reheis, M., Beann, J., Skipp, G., and Fisher, E. (2007). Airborne dust transport to the eastern pacific ocean off southern California: evidence from san clemente island. J. Geophys. Res. Atmospheres 112, 2006JD007577. doi:10.1029/2006JD007577
Muhs, D. R., Budahn, J. R., Johnson, D. L., Reheis, M., Beann, J., Skipp, G., et al. (2008). Geochemical evidence for airborne dust additions to soils in Channel Islands National Park, California. Geol. Soc. Am. Bull. 120, 106–126. doi:10.1130/B26218.1
Muhs, D. R., Cattle, S. R., Crouvi, O., Rousseau, D.-D., Sun, J., and Zárate, M. A. (2014). “Loess records,” in Mineral Dust. Editors P. Knippertz, and J.-B. W. Stuut (Dordrecht: Springer Netherlands), 411–441. doi:10.1007/978-94-017-8978-3_16
Nawrocki, J., Bogucki, A. B., Gozhik, P., Łanczont, M., Pańczyk, M., Standzikowski, K., et al. (2019). Fluctuations of the Fennoscandian Ice Sheet recorded in the anisotropy of magnetic susceptibility of periglacial loess from Ukraine. Boreas 48, 940–952. doi:10.1111/bor.12400
Nawrocki, J., Gozhik, P., Łanczont, M., Pańczyk, M., Komar, M., Bogucki, A., et al. (2018). Palaeowind directions and sources of detrital material archived in the Roxolany loess section (southern Ukraine). Palaeogeogr. Palaeoclimatol. Palaeoecol. 496, 121–135. doi:10.1016/j.palaeo.2018.01.028
Nawrocki, J., Polechońska, O., Boguckij, A., and Łanczont, M. (2006). Palaeowind directions recorded in the youngest loess in Poland and western Ukraine as derived from anisotropy of magnetic susceptibility measurements. Boreas 35, 266–271. doi:10.1080/03009480600584907
Necea, D., Fielitz, W., Kadereit, A., Andriessen, P. A. M., and Dinu, C. (2013). Middle Pleistocene to Holocene fluvial terrace development and uplift-driven valley incision in the SE Carpathians, Romania. Tectonophysics 602, 332–354. doi:10.1016/j.tecto.2013.02.039
Nenadić, D., Gaudenyi, T., Bogićević, K., and Tošović, R. (2016). Stratigraphic and lithologic characteristics of Pleistocene fluvial deposits in the Danube and Sava riparian area near Belgrade (Serbia). Stratigr. Geol. Correl. 24, 427–437. doi:10.1134/S0869593816040055
Nesbitt, H. W., Young, G. M., McLennan, S. M., and Keays, R. R. (1996). Effects of chemical weathering and sorting on the petrogenesis of siliciclastic sediments, with implications for provenance studies. J. Geol. 104, 525–542. doi:10.1086/629850
Nesbitt, H. W., and Young, G. M. (1984). Prediction of some weathering trends of plutonic and volcanic rocks based on thermodynamic and kinetic considerations. Geochim. Cosmochim. Acta 48, 1523–1534. doi:10.1016/0016-7037(84)90408-3
Nesbitt, H., and Young, G. (1996). Petrogenesis of sediments in the absence of chemical weathering: effects of abrasion and sorting on bulk composition and mineralogy. Sedimentology 43, 341–358. doi:10.1046/j.1365-3091.1996.d01-12.x
Nottebaum, V., Stauch, G., Hartmann, K., Zhang, J., and Lehmkuhl, F. (2015). Unmixed loess grain size populations along the northern Qilian Shan (China): relationships between geomorphologic, sedimentologic and climatic controls. Quat. Int. 372, 151–166. doi:10.1016/j.quaint.2014.12.071
Obreht, I., Buggle, B., Catto, N., Marković, S. B., Bösel, S., Vandenberghe, D. A. G., et al. (2014). The Late Pleistocene Belotinac section (southern Serbia) at the southern limit of the European loess belt: environmental and climate reconstruction using grain size and stable C and N isotopes. Quat. Int. 334, 10–19. doi:10.1016/j.quaint.2013.05.037
Obreht, I., Hambach, U., Veres, D., Zeeden, C., Bösken, J., Stevens, T., et al. (2017). Shift of large-scale atmospheric systems over Europe during late MIS 3 and implications for Modern Human dispersal. Sci. Rep. 7, 5848. doi:10.1038/s41598-017-06285-x
Obreht, I., Zeeden, C., Hambach, U., Veres, D., Marković, S. B., Bösken, J., et al. (2016). Tracing the influence of Mediterranean climate on Southeastern Europe during the past 350,000 years. Sci. Rep. 6, 36334. doi:10.1038/srep36334
Obreht, I., Zeeden, C., Hambach, U., Veres, D., Marković, S. B., and Lehmkuhl, F. (2019). A critical reevaluation of palaeoclimate proxy records from loess in the Carpathian Basin. Earth-sci. Rev. 190, 498–520. doi:10.1016/j.earscirev.2019.01.020
Obreht, I., Zeeden, C., Schulte, P., Hambach, U., Eckmeier, E., Timar-Gabor, A., et al. (2015). Aeolian dynamics at the Orlovat loess–paleosol sequence, northern Serbia, based on detailed textural and geochemical evidence. Aeolian Res. 18, 69–81. doi:10.1016/j.aeolia.2015.06.004
Ohta, T. (2004). Geochemistry of Jurassic to earliest Cretaceous deposits in the Nagato Basin, SW Japan: implication of factor analysis to sorting effects and provenance signatures. Sediment. Geol. 171, 159–180. doi:10.1016/j.sedgeo.2004.05.014
Özer, M., Orhan, M., and Işik, N. S. (2010). Effect of particle optical properties on size distribution of soils obtained by laser diffraction. Environ. Eng. Geosci. 16, 163–173. doi:10.1016/s0012-821x(01)00564-7
Pańczyk, M., Nawrocki, J., Bogucki, A. B., Gozhik, P., and Łanczont, M. (2020). Possible sources and transport pathways of loess deposited in Poland and Ukraine from detrital zircon U-Pb age spectra. Aeolian Res. 45, 100598. doi:10.1016/j.aeolia.2020.100598
Pötter, S., Schmitz, A., Lücke, A., Schulte, P., Obreht, I., Zech, M., et al. (2021). Middle to Late Pleistocene environments based on stable organic carbon and nitrogen isotopes of loess‐palaeosol sequences from the Carpathian Basin. Boreas 50, 184–204. doi:10.1111/bor.12470
Prins, M. A., and Vriend, M. (2007). Glacial and interglacial eolian dust dispersal patterns across the Chinese Loess Plateau inferred from decomposed loess grain-size records. Geochem. Geophys. Geosystems 8, Q07Q05. doi:10.1029/2006GC001563
Römer, W., Lehmkuhl, F., and Sirocko, F. (2016). Late Pleistocene aeolian dust provenances and wind direction changes reconstructed by heavy mineral analysis of the sediments of the Dehner dry maar (Eifel, Germany). Glob. Planet. Change 147, 25–39. doi:10.1016/j.gloplacha.2016.10.012
Rozycki, S. Z. (1967). Le sens des vents portant la poussière de loess à la lumière de I’analyse des formes d’accumulation du loess en Bulgarie et en Europe Centrale. Rev. Géomorphologie Dyn. 1–9, 33.
Schaetzl, R. J., Bettis, E. A., Crouvi, O., Fitzsimmons, K. E., Grimley, D. A., Hambach, U., et al. (2018). Approaches and challenges to the study of loess—introduction to the LoessFest special issue. Quat. Res. 89, 563–618. doi:10.1017/qua.2018.15
Schaffernicht, E. J., Ludwig, P., and Shao, Y. (2020). Linkage between dust cycle and loess of the last glacial maximum in Europe. Atmos. Chem. Phys. 20, 4969–4986. 10.5194/acp-20-4969-2020
Schatz, A.-K., Qi, Y., Siebel, W., Wu, J., and Zöller, L. (2015a). Tracking potential source areas of Central European loess: examples from Tokaj (HU), Nussloch (D) and Grub (AT). Open Geosci. 7, 678–720. doi:10.1515/geo-2015-0048
Schatz, A.-K., Scholten, T., and Kühn, P. (2015b). Paleoclimate and weathering of the Tokaj (Hungary) loess–paleosol sequence. Palaeogeogr. Palaeoclimatol. Palaeoecol. 426, 170–182. doi:10.1016/j.palaeo.2015.03.016
Scheidt, S., Berg, S., Hambach, U., Klasen, N., Pötter, S., Stolz, A., et al. (2021). Chronological assessment of the Balta Alba Kurgan loess-paleosol section (Romania)—a comparative study on different dating methods for a robust and precise age model. Front. Earth Sci. 14, 312. 10.3389/feart.2020.598448
Schulte, P., Lehmkuhl, F., Steininger, F., Loibl, D., Lockot, G., Protze, J., et al. (2016). Influence of HCl pretreatment and organo-mineral complexes on laser diffraction measurement of loess–paleosol-sequences. CATENA 137, 392–405. doi:10.1016/j.catena.2015.10.015
Schulte, P., and Lehmkuhl, F. (2018). The difference of two laser diffraction patterns as an indicator for post-depositional grain size reduction in loess-paleosol sequences. Palaeogeogr. Palaeoclimatol. Palaeoecol. 509, 126–136. doi:10.1016/j.palaeo.2017.02.022
Skurzyński, J., Jary, Z., Kenis, P., Kubik, R., Moska, P., Raczyk, J., et al. (2020). Geochemistry and mineralogy of the Late Pleistocene loess-palaeosol sequence in Złota (near Sandomierz, Poland): implications for weathering, sedimentary recycling and provenance. Geoderma 375, 114459. doi:10.1016/j.geoderma.2020.114459
Skurzyński, J., Jary, Z., Raczyk, J., Moska, P., Korabiewski, B., Ryzner, K., et al. (2019). Geochemical characterization of the late Pleistocene loess-palaeosol sequence in Tyszowce (Sokal Plateau-Ridge, SE Poland). Quat. Int. 502, 108–118. doi:10.1016/j.quaint.2018.04.023
Smalley, I., Jefferson, I. F., Dijkstra, T. A., and Derbyshire, E. (2001). Some major events in the development of the scientific study of loess. Earth-Sci. Rev. 54, 5–18. doi:10.1016/s0012-8252(01)00038-1
Smalley, I., and Leach, J. A. (1978). The origin and distribution of the loess in the Danube basin and associated regions of East-Central Europe—a review. Sediment. Geol. 21, 1–26. doi:10.1016/0037-0738(78)90031-3
Smalley, I., O’Hara-Dhand, K., Wint, J., Machalett, B., Jary, Z., and Jefferson, I. (2009). Rivers and loess: the significance of long river transportation in the complex event-sequence approach to loess deposit formation. Quat. Int. 198, 7–18. doi:10.1016/j.quaint.2008.06.009
Staubwasser, M., Drăguşin, V., Onac, B. P., Assonov, S., Ersek, V., Hoffmann, D. L., et al. (2018). Impact of climate change on the transition of Neanderthals to modern humans in Europe. Proc. Natl. Acad. Sci. 115, 9116–9121. doi:10.1073/pnas.1808647115
Stephens, M., Krzyszkowsku, D., Ivchenko, A., and Majewski, M. (2002). Palaeoclimate and pedosedimentary reconstruction of a Middle to Late Pleistocene loess-palaeosol sequence, Prymorske, SW Ukraine. Stud. Quat. 19, 3–17.
Stevens, T., Palk, C., Carter, A., Lu, H., and Clift, P. D. (2010). Assessing the provenance of loess and desert sediments in northern China using U-Pb dating and morphology of detrital zircons. Geol. Soc. Am. Bull. 122, 1331–1344. doi:10.1130/B30102.1
Stevens, T., Sechi, D., Bradák, B., Orbe, R., Baykal, Y., Cossu, G., et al. (2020). Abrupt last glacial dust fall over southeast England associated with dynamics of the British-Irish ice sheet. Quat. Sci. Rev. 250, 106641. doi:10.1016/j.quascirev.2020.106641
Tecsa, V., Gerasimenko, N., Veres, D., Hambach, U., Lehmkuhl, F., Schulte, P., et al. (2020). Revisiting the chronostratigraphy of Late Pleistocene loess-paleosol sequences in southwestern Ukraine: OSL dating of Kurortne section. Quat. Int. 542, 65–79. doi:10.1016/j.quaint.2020.03.001
Thamó-Bozsó, E., Kovács, L. Ó., Magyari, Á., and Marsi, I. (2014). Tracing the origin of loess in Hungary with the help of heavy mineral composition data. Quat. Int. 319, 11–21. doi:10.1016/j.quaint.2013.04.030
Tugulan, L. C., Duliu, O. G., Bojar, A.-V., Dumitras, D., Zinicovskaia, I., Culicov, O. A., et al. (2016). On the geochemistry of the Late Quaternary loess deposits of Dobrogea (Romania). Quat. Int. 399, 100–110. doi:10.1016/j.quaint.2015.06.062
Újvári, G., Klötzli, U., Kiraly, F., and Ntaflos, T. (2013). Towards identifying the origin of metamorphic components in Austrian loess: insights from detrital rutile chemistry, thermometry and U–Pb geochronology. Quat. Sci. Rev. 75, 132–142. doi:10.1016/j.quascirev.2013.06.002
Újvári, G., and Klötzli, U. (2015). U–Pb ages and Hf isotopic composition of zircons in Austrian last glacial loess: constraints on heavy mineral sources and sediment transport pathways. Int. J. Earth Sci. 104, 1365–1385. doi:10.1007/s00531-014-1139-x
Újvári, G., Varga, A., Ramos, F. C., Kovács, J., Németh, T., and Stevens, T. (2012). Evaluating the use of clay mineralogy, Sr–Nd isotopes and zircon U–Pb ages in tracking dust provenance: an example from loess of the Carpathian Basin. Chem. Geol. 304–305, 83–96. doi:10.1016/j.chemgeo.2012.02.007
Újvári, G., Varga, A., Raucsik, B., and Kovács, J. (2014). The Paks loess-paleosol sequence: a record of chemical weathering and provenance for the last 800ka in the mid-Carpathian Basin. Quat. Int. 319, 22–37. doi:10.1016/j.quaint.2012.04.004
Vandenberghe, J., Kasse, C., Popov, D., Markovic, S. B., Vandenberghe, D., Bohncke, S., et al. (2018). Specifying the external impact on fluvial lowland evolution: the last glacial tisza (Tisa) catchment in Hungary and Serbia. Quaternary 1, 14. doi:10.3390/quat1020014
Vandenberghe, J., Renssen, H., van Huissteden, K., Nugteren, G., Konert, M., Lu, H., et al. (2006). Penetration of atlantic westerly winds into Central and East Asia. Quat. Sci. Rev. 25, 2380–2389. doi:10.1016/j.quascirev.2006.02.017
Varga, G., Kovács, J., and Újvári, G. (2013). Analysis of saharan dust intrusions into the Carpathian Basin (central Europe) over the period of 1979–2011. Glob. Planet. Change 100, 333–342. doi:10.1016/j.gloplacha.2012.11.007
Veres, D., Lane, C. S., Timar-Gabor, A., Hambach, U., Constantin, D., Szakács, A., et al. (2013). The Campanian Ignimbrite/Y5 tephra layer–A regional stratigraphic marker for Isotope Stage 3 deposits in the Lower Danube region, Romania. Quat. Int. 293, 22–33.
Veres, D., Tecsa, V., Gerasimenko, N., Zeeden, C., Hambach, U., and Timar-Gabor, A. (2018). Short-term soil formation events in last glacial east European loess, evidence from multi-method luminescence dating. Quat. Sci. Rev. 200, 34–51. doi:10.1016/j.quascirev.2018.09.037
Vlaminck, S., Kehl, M., Rolf, C., Franz, S. O., Lauer, T., Lehndorff, E., et al. (2018). Late Pleistocene dust dynamics and pedogenesis in Southern Eurasia – detailed insights from the loess profile Toshan (NE Iran). Quat. Sci. Rev. 180, 75–95. doi:10.1016/j.quascirev.2017.11.010
Willmes, C. (2015). LGM sealevel change (HiRes), CRC806 database, collaborative research centre 806.Available at:http://crc806db.uni-koeln.de/dataset/show/lgm-sealevel-change-hires1436532921 [Accessed March 9, 2016].
Wolf, D., Ryborz, K., Kolb, T., Zapata, R. C., Vizcaino, J. S., Zöller, L., et al. (2019). Origins and genesis of loess deposits in central Spain, as indicated by heavy mineral compositions and grain-size variability. Sedimentology 66, 1139–1161. doi:10.1111/sed.12539
Wright, J. S. (2001). “Desert” loess versus “glacial” loess: quartz silt formation, source areas and sediment pathways in the formation of loess deposits. Geomorphology 36, 231–256. doi:10.1016/S0169-555X(00)00060-X
Yang, S. Y., Li, C. X., Yang, D. Y., and Li, X. S. (2004). Chemical weathering of the loess deposits in the lower Changjiang Valley, China, and paleoclimatic implications. Quat. Int. 117, 27–34. doi:10.1016/S1040-6182(03)00113-7
Zeeden, C., and Hambach, U. (2021). Magnetic susceptibility properties of loess from the Willendorf archaeological site: implications for the syn/post-depositional interpretation of magnetic fabric. Front. Earth Sci. 8 .10.3389/feart.2020.599491
Zeeden, C., Hambach, U., Veres, D., Fitzsimmons, K., Obreht, I., Bösken, J., et al. (2018). Millennial scale climate oscillations recorded in the Lower Danube loess over the last glacial period. Palaeogeogr. Palaeoclimatol. Palaeoecol. 509, 164–181. doi:10.1016/j.palaeo.2016.12.029
Zhang, H., Lu, H., Stevens, T., Feng, H., Fu, Y., Geng, J., et al. (2018). Expansion of dust provenance and aridification of asia since ∼7.2 ma revealed by detrital zircon U‐Pb dating. Geophys. Res. Lett. 45. doi:10.1029/2018GL079888
Keywords: loess provenance, aeolian sediments, Southeastern Europe, geochemistry, sedimentary pathways
Citation: Pötter S, Veres D, Baykal Y, Nett JJ, Schulte P, Hambach U and Lehmkuhl F (2021) Disentangling Sedimentary Pathways for the Pleniglacial Lower Danube Loess Based on Geochemical Signatures. Front. Earth Sci. 9:600010. doi: 10.3389/feart.2021.600010
Received: 28 August 2020; Accepted: 18 February 2021;
Published: 20 April 2021.
Edited by:Annett Junginger, Department of Geosciences, University of Tubingen, Germany
Reviewed by:Mingrui Qiang, South China Normal University, China
Roberto Adrián Scasso, University of Buenos Aires, Argentina
Copyright © 2021 Pötter, Veres, Baykal, Nett, Schulte, Hambach and Lehmkuhl. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Stephan Pötter, email@example.com