2023 | 76/2 | 73–85 | 10 Figs. | www.geologia-croatica.hr Journal of the Croatian Geological Survey and the Croatian Geological Society 1. INTRODUCTION Soil pedocomplexes are unique, highly sensitive natural multi- component systems capable of forming valuable archives for nu- merous signals that are indicative of intense palaeoclimatic and palaeoenvironmental changes over time (BRONGER et al., 1998; BRONGER, 2003). Soil memory preserves diverse components which vary due to the type of evolutionary processes that have taken place during the time since deposition. Therefore, combin- ing studies of these complex physical and geo-chemical compo- nents allows detailed reconstruction of geoecological trends that drive the naturally induced soil transformation over different spa- tio-temporal scales (e.g. VASILJEVIĆ et al., 2011; 2014). Unlike modern soils, palaeosols are fossil relics which have undergone pedogenic processes in the geological past, capturing palaeocli- mate components of ancient landscapes (NETTLETON et al., 2000). The Danube watershed encompasses complex geological terrains comprised of diverse geomorphic units, such as the Alps, Bohemian Massif, Carpathians, Dinarides, and the Balkan Mountains. The catchment has been a valuable source of thick loess deposits for study, which are especially well preserved within the Middle Danube basin, with sections dating back until at least the mid-Quaternary (e.g. FENN et al., 2022; WACHA et Initial results of the colorimetric indices of the oldest exposed pedocomplex (Titel loess plateau, Serbia) Tin Lukić1, Milica G. Radaković1*, Rastko Marković2, Warren Thompson3, Tanja Micić Ponjiger1, Biljana Basarin1, Nemanja Tomić4, Vladimir M. Tomović5, Jovanka Popov Raljić6, Milivoj B. Gavrilov1, Aleksandar Valjarević7, György Sipos8, Dejan Filipović7, Cezar Morar9 and Slobodan B. Marković1,10 1 University of Novi Sad, Chair of Physical Geography, Department of Geography, Tourism and Hotel Management, Faculty of Sciences, Trg Dositeja Obradovića 3, 21000 Novi Sad, Serbia; (*corresponding author: milicar@dgt.uns.ac.rs) 2 University of Niš, Chair of Regional Geography, Department of Geography, Faculty of Sciences, Višegradska 33, 18000 Niš, Serbia 3 Radiation Physics, Technical University of Denmark, Risë Campus, Fysikvej, 311, Kgs. Lyngby, Denmark 4 University of Novi Sad, Chair of Tourism, Department of Geography, Tourism and Hotel Management, Faculty of Sciences, Trg Dositeja Obradovića 3, 21000 Novi Sad, Serbia 5 University of Novi Sad, Faculty of Technology, Bulevar cara Lazara 1, 21000 Novi Sad, Serbia 6 Singidunum University, Faculty of Tourism and Hotel Management, Danijelova 32, 11010 Belgrade, Serbia 7 University of Belgrade, Faculty of Geography, Studentski Trg 3/III, 11000 Belgrade, Serbia 8 University of Szeged, Geomorphological and Geochronological Research Group, Department of Geoinformatics, Physical and Environmental Geography, H-6722 Szeged, Egyetem u. 2-6, Hungary 9 University of Oradea, Department of Geography, Tourism and Territorial Planning, 410087 Oradea, Romania 10 Serbian Academy of Sciences and Arts, Kneza Mihaila 35, 11000 Belgrade, Serbia doi: 10.4154/gc.2023.05 Abstract In this study we present an in-depth description of the colorimetric values for the lowest section of the Dukatar Loess Palaeosol Sequence (LPS) pedocomplex S5. Formed during the Marine Isotope Stage (MIS) 13-15, it represents the oldest pedocomplex exposed at the base of the Ti- tel loess plateau (TLP), near the confluence of the Tisa and Danube rivers in Vojvodina (north- ern Serbia). The results of low-field magnetic susceptibility measurements (χlf) were compared to colour properties (obtained by conventional methods as well as instrumental measuring) and quantified Soil Development Indices (SDI). Of these measurements we found that the Redness Index (RI1) yielded the most useful results, as this index appears most sensitive to lithological changes and soil development intensity. It was also observed that a high level of correlation ex- isted between χlf, and a* chromaticity. The initial results of this study highlight the utility of colo- rimetric methods as an interdisciplinary tool when evaluating the presence of ferromagnetics, and the application of rock magnetism to the Middle and Upper Pleistocene LPS of the Middle Danube Basin. The presented approach can be used to observe the evolution of climatic and ecological conditions in the given study area, and for establishing correlations between sites ex- tending over the Eurasian LPS provinces. al., 2021; BANAK et al., 2016; GALOVIĆ & PEH, 2016; SÜMEGI et al., 2019; MARKOVIĆ et al., 2013). Therefore, Danubian river catchment loess deposits preserve some of the longest terrestrial records of palaeoenvironmental and palaeoclimatic changes in Central and Southeastern Europe, and possess great palaeoeco- logical significance (e.g. RUBINIĆ et al., 2018; GALOVIĆ, 2014, 2016; BANAK et al., 2013; MARKOVIĆ et al., 2015). In Serbia, loess-palaeosol sequences (LPS) have been documented predom- inantly in the northern part of the country (MARKOVIĆ et al., 2006; MARKOVIĆ et al., 2008, MARKOVIĆ et al., 2015; and references therein), but with a few exceptions occurring in south- ern parts of the Carpathian basin (BASARIN et al., 2014; OBREHT et al., 2014; BOESKEN et al., 2017; MARKOVIĆ et al., 2021b). Palaeosols developed on loess can be traced over ex- tensive parts of Europe and Asia (KUKLA, 1987), and it is pos- sible to correlate their chronostratigraphy based on different prox- ies and thus reconstruct pan-continental climate fluctuations. The most common proxies encompass: magnetic properties of soils and loess (RADAKOVIĆ et al., 2019), grain size or grain size in- dices (VANDENBERGHE, 2013), luminescence dating (BERGER et al., 1992; PERIĆ et al., 2019), stable isotopes (SCHATZ et al., 2011), chemical analysis (POTTER et al., 2021; Article history: Manuscript received June 13, 2022 Revised manuscript accepted May 12, 2023 Available online June 16, 2023 Keywords: loess, pedocomplex, soil color indices, Vojvodina, Titel loess plateau, Dukatar, North Serbia G eo lo gi a C ro at ic a Geologia Croatica 76/274 FENN et al., 2022), tephrochronology (TIMAR-GABOR et al., 2017), palynology (ZHANG et al., 2017), malacology (LUDWIG et al., 2021), and colour indices (SPRAFKE et al., 2020). VASILJEVIĆ et al. (2014) asserts that the aforementioned fea- tures of this widely spread sediment and pedocomplexes have proven to be of the utmost scientific, archaeological and agricul- tural significance within the Eurasian loess belt. This study fo- cuses on the colour indices proxy. Different colours of LPS were recorded by the pioneering loess researcher, Aloisius Ferdinando Marsigli in his “Danubius Pannonico Mysicus” (MARSIGLI, 1726). One and a half centu- ries later, the geologist Aleksandar Popovics (1847-1877) wrote about the transition of yellow colour of the loess sediments on the Fruška Gora Mountain in Serbia. This study recorded the appear- ance of a red colour in loess, which is today known as a feature of the regions palaeosols (POPOVICS, 1876). The main topic of our study is the colour of the Dukatar LPS, which was described primarily by JELENA MARKOVIĆ- MARJANOVIĆ (1950), who noted the strong red colour visible at the bottom of the profile. In Danubian loess stratigraphy, this layer is referred to as S5; where “S” stands for soil and 5 means that it was formed in the fifth interglacial period, counting from the present (MARKOVIĆ et al., 2015). Before the above-men- tioned proxies were applied to loess research, the colour intensity was used to make connections between multiple loess profiles in Serbia and allude to certain aspects of loess plateau geomorpho- logical evolution. Counting of the loess (yellow) or palaeosol (red) layers was utilised in around 70 papers in the former Yugoslavia in an attempt to correlate loess profiles in Serbia. As time passed, the more detailed the descriptions became, and the more differ- entiated the LPS became (MARKOVIĆ-MARJANOVIĆ, 1965, and references therein). Presently, many LPS in Serbia still do not have determined colorimetric perspectives. The work of LUKIĆ et al. (2014) was the first detailed colour investigation, describing the Orlovat LPS (Fig. 1), in the Banat region of Serbia. A strong correlation be- tween frequency-dependent magnetic susceptibility and Redness Index was discovered. The Orlovat LPS spans the last glacial cy- cle, with the lowest palaeosol belonging to Marine Isotope Stage (MIS) 5, or S1 in Danube loess stratigraphy. The second LPS in Serbia to be characterized in terms of its colorimetric properties was the profile at Stalać (OBREHT et al., 2016). This archive is a valuable record of the transition from the Mediterranean to con- tinental climate in the central Balkan region, which occurred be- tween MIS 9 and MIS 1. The most recent profile to be described in terms of its colour proxies is the Zemun LPS, on the right bank of the Danube (LAAG et al., 2021). The respective authors ex- panded the search to MIS 11, or palaeosol S4. In general, it has been found that the luminance of the LPS sediments follows gla- cial cycles; although in Zemun the S4, which is the palest of all palaeosols, the signal is not as distinguishable as in other palaeo- sols. Fortunately, there are other parameters of colour such as redness and blueness, with indices derived from them, that assist our interpretations of the LPS, and pedogenetic processes (ZYKINA et al., 2021). Interestingly, coloured soils with shades of red and brown are already recognized as touristic attractions, for example in Chamarel, Mauritius (SHETH, et al., 2010). High- resolution colour measurements provide detailed palaeoenviron- mental markers, which can even be linked to the migration of anatomically modern humans (SPRAFKE et al., 2020). This study aims to characterize the oldest pedocomplex of the Titel Loess Plateau, in terms of its magnetic and colorimetric properties, and the palaeoecological features that outline the sci- entific and educational values of the given LPS. These sites are somewhat endangered by various factors, as presently they re- main for the most part unprotected. 2. MATERIAL AND METHODS 2.1. Study area The Dukatar site (45°16’59” N, 20°14’33” E) is located 500 m downstream from the settlement of the same name and represents the most complete series of loess-palaeosol sequences of the TLP (Fig. 1). It is situated at the foot of a loess section and is available for sampling only during the low water level of the Tisa River. The sampled S5 pedocomplex part of the Dukatar sequence, is approximately 2.7 m thick and was formed during the Early Mid- dle Pleistocene (e.g. MARKOVIĆ et al., 2012). This site belongs to the highest protected zone of the Special Nature Reserve TLP, as indicated with a thick red line in Fig. 1 (SGRS, 2012). The site is situated in the contact zone of the Tisa River and TLP. After building the dam at Novi Bečej, north of the TLP, the Tisa River has a controlled regime. When the water level is high (April- May), it is not possible to reach the Dukatar LPS, while in Sep- tember and October it is accessible. The maximum recorded am- plitude of the Tisa river at the TLP is 876 cm (PAVIĆ, 2006). In the north of Dukatar there are two famous LPS: Veliki Surduk (Big Gully) (PERIĆ et al., 2019) and Feudvar (NAMIER et al., 2021). In the south, three more exist: Rogulić Surduk, Stara ciglana (Old brickyard) and Kalvarija. The list of protected LPS in Serbia is quite short: Surduk (VUJIČIĆ et al., 2011), Batajnica (VIŠNJIĆ et al., 2016), Stari Slankamen (VASILJEVIĆ et al., 2011), and Zemun (ŠARIĆ, 2008). 2.2. Sampling strategy A field survey was undertaken prior to the onset of sampling the selected profiles for palaeoclimatic and palaeoecological recon- struction using colorimetric parameters. The field survey, carried out in 2013, involved a detailed description of LPS, only then were samples collected for laboratory analysis. Efforts were com- plicated by erosion and rockfall over the surface of the profile, resulting in colluvial deposits of considerable thickness; this greatly affected the consistency of the laboratory results. Initial steps included the removal of colluvial deposits and recent veg- etation in order to expose the undisturbed surface of the loess section. The next step involved cutting a vertical trench and sam- pling of the profile. On the TLP, at the site located near the Du- katar settlement, a strongly rubified S5 pedocomplex was ex- posed and sampled. Sampling of the S5 pedocomplex on the Dukatar profile was undertaken at a high resolution of 3 cm from the bottom part of the L5 unit, across the upper part of the A ho- rizon, on into the Ck horizon in the basal part of this strongly ru- bified palaeosol. A total of 100 bulk samples were taken for pal- aeopedological research (Fig. 2A). In addition, during sampling of this pedocomplex particular care was taken to avoid contact with widespread carbonate concretions in order to obtain valid results of the laboratory analysis. 2.3. Magnetic and colorimetric measurements Measurement of the χlf of palaeosol samples from the Dukatar profile were undertaken at the laboratory of the Loess Research Group at the Department of Geography, Tourism and Hotel Man- agement, Faculty of Sciences, University of Novi Sad, using a portable Bartington susceptibility meter set to 300 Hz. Determin- G eologia C roatica 75 ing the colour of samples in wet conditions, drying of the sam- ples, and re-examination of the colour in dry conditions was car- ried out during March and April 2014 at the Laboratory for Paleoecological Reconstruction (LAPER), Department of Bio- logy and Ecology, Faculty of Sciences, University of Novi Sad. Following the approach of LUKIĆ et al. (2014), the samples were dried in the oven for at least 12 h at 39 °C. Instrumental measure- ment of loess-palaeosol sample colour involved the use of a tris- timulus colorimeter, and expressing the obtained values in the CIE L a* b* colorimetric system (Hünter Lab, 2001). This was undertaken at the Department of Food Preservation Technology, Faculty of Technology, University of Novi Sad. The calculation of Rubification (RI) and Melanization (MI) Indices was per- formed using the following the approaches provided by HARDEN (1982), HARDEN and TAYLOR (1983), VIDIĆ et al. (2004) and VIDIĆ and LOBNIK (1997). This methodology was based on the Lotus Template by TAYLOR (1988), as it possesses robust ap- proaches for quantification procedures and index calculations. The intensity of rubification indicates the formation of haematite crystals suspended in the soil matrix, which are responsible for the distinctive bright red colour of the soil. The MI indicates ac- cumulation of humus and humic substances in the soil, and the RI quantifies the degree of soil reddening relative to the parent material. MI is used to calculate the degree of darkening (VIDIĆ et al., 2004, SCHAETZL, and ANDERSON, 2005). The Redness Index (RI1), proposed by HARDEN (1982), can be seen as a proxy for semi-quantitative soil reddening, and is indicative of the proportion of haematite in the pedocomplex. The Figure 1. The upper panel shows the wider region of northern Serbia and the higher resolution part covers the Tisa and Danube confluence, with the loess pal- aeosol sequences over the TLP; the lower panel shows the geological map and protected area of the highest importance indicated with a red line. G eo lo gi a C ro at ic a Geologia Croatica 76/276 Redness Rating Index (RR) translates the increase in redness be- tween a soil or palaeosol and its parent material into a numerical value by comparing the changes in hue and chroma (BUGGLE et al., 2014). Based on the obtained parameters, RI1 (BARRON & TORRENT, 1986; ROSSEL et al., 2006) (Equation 1) and RR were calculated (TORRENT et al., 1983), (Equation 2): RI L a b b L 1 102 2 0 5 8 2 6 = + × × ( * * ) * . . , (1) RR H C V = − ×( )10 , (2) where L stands for lightness, a* for green-red colour, b* for blue- yellow colour calculated by Konica Minolta Chroma meter CR- 400, H for hue, C for chroma, and V for value in the Munsell col- our chart. As previously mentioned, the colour was determined in both dry and wet conditions. Their difference was used in the calculation, following the approach by Buggle et al. (2014). The indices described above belong to the group of “true” index val- ues, presented by the non-dimensional number. 3. RESULTS 3.1. Pedostratigraphy and Chronostratigraphy The Dukatar profile is a section where the oldest sediments on the TLP were discovered and which has a thickness of 10 m (e.g. MARKOVIĆ et al., 2012, 2015). Both the oldest loess unit (L5) and the oldest exposed pedocomplex (S5) of TLP are observed within this site. For the purposes of this study, the lowest part of the L5 loess unit and the S5 pedocomplex, with a total thickness of approximately 2.7 m, were sampled. The sampled sequence (300 cm thick) extends from the contact between the upper part of the humus accumulating horizon A (~ 20 cm thick) and the basal part of the L5 loess unit. Beneath the A horizon there is a well-developed dark brown rubified cambic Bw horizon, with moderately developed coarse polyhedral structure that is approxi- mately 110 cm thick (Fig. 2B). In the central part of the pedo- complex, S5 is a Ck horizon containing many carbonate concre- tions, and bioturbation caused by roots that are now preserved as rhizoliths. The lower Bw is a weakly rubified horizon disturbed by many hydromorphic features. In the bottom part of the pedo- complex S5 there is a horizon containing massive carbonate con- Figure 2. A) The picture of Dukatar LPS (photo by Lukić, T., September 2014) with indicated lithostratigraphy; B) Pedostratigraphy of Dukatar LPS. Figure 3. Histograms indicating A) L; B) a* chromaticity; C) b* chromaticity; D) The same parameters plotted on the depth scale of Dukatar LPS. G eologia C roatica 77 cretions overlying loessoid material of alluvial origin (also with numerous hydromorphic forms). The description presented here corresponds well to the previous palaeopedological interpreta- tions put forward by BRONGER, (1976, 2003), MARKOVIĆ et al., 2012, 2015; BASARIN et al., 2014). Chronostratigraphic mod- els of the LPSs in Serbia, presented by MARKOVIĆ et al. (2012, 2015) and BASARIN et al. (2014), provide correlation between the pedocomplex S5 at the Dukatar profile (TLP) with the ac- cepted MIS 13-15 S5. MARKOVIĆ et al. (2015) indicated a re- markable accordance between Danubian and Chinese loess re- cords, thus opening up the possibility for a transcontinental correlation of European, Central Asian and Chinese LPSs, using a standardised nomenclature and chronostratigraphic model. 3.2. Colorimetric and magnetic interpretation The histograms in Fig. 3 illustrates the distribution of results for the colorimetric measurements L, a*, and b* parameters from the Dukatar profile (S5). Based on instrumental colour measure- ments, marked variations in L colorimetric values were observed. These variations displayed the following ranges: between the bot- tom part of the L5 unit, and the upper part of the A horizon within the S5 pedocomplex (64-66 units), and the bottom part of the Bw horizon, and its contact with Ck basal horizon (66-75 units) (Fig. 3A and Fig. 3D). The lowest values (56) are associated with the Bw dark brown horizon, at 150 cm depth. Variations are particu- larly noticeable between 180 and 195 cm from the upper part of the examined section (values vary between 61 and 68 units) and result from the presence of large diameter carbonate concretions that are characteristic of this part of the pedocomplex. Again, the values of L increase towards the calcified Ck horizon in the basal part of the S5 pedocomplex. The values of a* chromaticity (Fig. 3B and Fig. 3D) display variation that correspond to the detected variations in L colori- metric values. Note that the a* values are mirrored to symmetri- cally match the L. They range from 3.4 units to 8.2 units at a depth of 15 cm to 120 cm. Subsequently, the values gradually decline to a depth of 192 cm, where the lowest value (3.6) was recorded. The a* chromaticity gradually increases to the value of 6.5 units at a depth of 252 cm, and then begins to continuously decline to a depth of 270 cm (2.4), which is the lowest value in the pedo- complex. The highest measured a* values are 8.2 and were re- corded in the upper part of the examined section, in the central part of the Bw rubified horizon. Observations of b* chromaticity (Fig. 3C and Fig. 3D) dis- play fluctuating values that vary slightly starting from the top to the middle of the examined profile. These values range from 19.5 to 23.9 units. In the transitional part between the L5 unit, and the upper part of the A horizon within the S5 pedocomplex the val- ues rise from 21.6 to 19.8 units (between the bottom part of the Bw horizon and its contact with the Ck basal horizon). As noted in the previously mentioned colorimetric parameters values, stronger fluctuations in b* chromaticity values are observed in the central part of the S5 pedocomplex. In general, these values gradually increase from the upper part of the pedocomplex to 171 cm, ranging from 21.5 units to 23.9 units, respectively. After this, b* chromaticity values start to gradually decrease towards the value of 21.2 units observed at a Ck horizon, at the depth of 192 cm. The b* chromaticity values from this point display continu- ous fluctuations that are most likely to be the result of the pres- ence of hydromorphic features and carbonate concretions that are common at the basal part of the Bw horizon and its contact with the Ck horizon. The depth plot of these parameters is presented in Fig. 3D. Results of the χlf measurements indicate that χlf variations are in good agreement with the pedostratigraphy (Fig. 4A). The his- togram Fig. 4B illustrates the results of χlf for the Dukatar profile (S5). Values of χlf range from 40×10-8 m3 kg-1 in the upper part up to 132×10-8 m3 kg-1 at a depth of 114 cm. Thereafter, the values gradually decline to 43×10-8 m3 kg-1 at a depth of 171 cm, after which they begin to progressively increase to the value of 93×10-8 m3 kg-1, observed at 229 cm from the top of the examined section (Fig. 4A). A decrease in the χlf value was identified in the deeper parts of the horizon, with a sharp transition detected at a depth of 270 cm (61×10-8 m3 kg-1) and continuing to the basal part of the profile, where the values of 18×10-8 m3 kg-1 were recorded. The values represent the absolute minimum found on the examined profile. These low χlf values are probably the result of gleization due to variation in underground and surface water levels. The re- sults display very similar variations to those recorded during pre- vious research of LPSs in Serbia; also to those on the Eastern and Figure 4. A) The low-frequency magnetic susceptibility measurements (χlf) for the Dukatar profile; B) Histogram for χlf. G eo lo gi a C ro at ic a Geologia Croatica 76/278 Western coast of the Black Sea and in Central and East Asia (China) (e.g. HELLER & EVANS, 1995; BUGGLE et al., 2009; BUGGLE et al., 2014; MARKOVIĆ et al., 2009, 2012; BASARIN et al., 2014). In Fig. 5 the four panels present the linear and polynomial fit of L and a* with magnetic properties of the S5 pedocomplex. The correlation between the χlf and L of the samples in Dukatar is very strong; high values of L correlate to the low values of χlf. This can be explained by the higher concentrations of calcium carbonate (CaCO3) along some horizons of the pedocomplex, which do not possess Fe minerals, and provides the white colour. The coeffi- cient of determination R2 is best expressed as the polynomial fit, and has the value of 0.87. In contrast, the higher a* values cor- respond to higher χlf, and the trend is steadily rising. The high values of a* present the red colour which most likely comes from ferromagnetic haematite. The R2 value is also fitted better in poly- nomial case, with the value of 0.97. 3.3. Soil Developed Indices (SDI): Rubification (RI) and melanization (MI) indices In moist conditions, the hues determined using a Munsell Soil Colour Chart for the first 66 cm of the examined sequence hues, range between 10YR3/6 and 10YR5/6. In the next 90 cm the hues vary from 7.5YR3/4 to 7.5YR4/6; for a depth from 159 cm to 288 cm hues range from 10YR3/4 and 10YR5/6. In the remaining part of the profile, the hue has the value of 2,5Y5/6. At the same depths in dry conditions, the values ranged from 2.5Y6/4 to 10YR7/4, then from 10YR6/6 to 7.5YR5/6, from 10YR4/6 to 10YR7/4 and 2.5Y7/4 respectively. Quantified soil development indices (SDI), rubification (RI) and melanization (MI) are shown in Fig. 6. The RI shows rela- tively low variability in the first 57 cm from the upper part of the sampled section. Calculated values vary from 10 to 40, with only one extreme value observed at 24 cm (10). From 60 cm to 150 cm Figure 5. Polynomial amd linear fit of the magnetic susceptibility data χlf, L, and a* for Dukatar LPS. The y and x in the regression equalions refer to the mentioned parameters. Figure 6. The Rubification index (RI), Melanization index (MI), Redness Index (RI1), and the Redness Rating index (RR) for the Dukatar LPS. G eologia C roatica 79 of the section, RI values progressively increase from 40 to 100. A pronounced decline in values was observed at a depth between 165 cm and 170 cm (90 to 40) which is consistent with measured a* chromaticity values. At a depth of 216 cm of the examined section, RI values gradually decline from 60 to 20. Melanization index (MI) in the first 63 cm of the profile is characterized by moderate variations ranging from –10 to –30. Thereafter, index values vary from –30 to –50 to a depth of 84 cm. Then a progressive decline in values occurs, with –10 being the most extreme one, detected at 171 cm from the top of the pro- file. At a depth between 171 cm and 258 cm relatively stable MI values (–30) are present. In this part of the section, only one value of –40 was noted (at a depth of 222 cm). From 258 cm, a gradual decrease in values ranging from –30 to –10 was apparent. RI and MI values, obtained through soil sample analysis using a Munsell Soil Colour Chart (observed with the naked eye) also largely cor- respond to the results from measurement of a* chromaticity and L lightness values using a tristimulus colorimeter. 3.4. Soil colour proxies: Redness (RI1) and Redness rating (RR) indices The calculated redness index values (RI1) are shown in Fig. 6. Values of RI1 that may relate to relative changes in haematite content (BARRON & TORRENT, 1986) in the examined pedo- complex, appear to correspond to the detected variations recorded in the signal of χlf (Fig. 4A). These values remain relatively low in the top 30 cm of the upper part of the section, and vary between 0.12 and 0.11, presumably indicating a slightly lower haematite content. After this depth, the possible haematite concentration gradually increases to a depth of 126 cm, with a maximum value of 0.30. Thereafter, the RI1 values decrease to a depth of 189 cm, where the minimum value was recorded (0.10). The values then gradually increase to a depth of 234 cm from the upper part (0.21), before decreasing again towards the basal part of the sequence. Consequently, it can be noted that RI1 values generally corre- spond to the measured values of L, a*, and b* colorimetric indi- cators and SDI RI and MI. Soil redness rating index (RR) values, associated with the pedostratigraphy of the section are also shown in Fig. 6. There- fore, RI1 as well as RR indices could reflect the haematite content in the examined Dukatar S5 pedocomplex in a semi-quantitative manner. The values appear relatively uniform for the first 78 cm from the upper part of the section and amount to 0, with one ex- ception, observed at 72 cm (7.5). Thereafter, the index values in- crease to a depth of 156 cm, ranging from 7.5 to 11.25. At a depth of 171 cm the values again drop to zero. The basal part of this site is characterized by the highest RR values of 12. Soil RR values correspond to a large degree to the values of L, a*, and b* colo- rimetric indicators, SDI RI and MI and χlf values for the first 156 cm from the top of the examined section. Values detected in the second half of the Dukatar section at the contact between hori- zons Bw and Ck show slightly higher deviations, especially in the basal part. This can be explained by the presence of massive carbonate concretions. 4. DISCUSSION 4.1. Magnetic susceptibility and colour as indicators of palaeoclimatic and palaeoecological processes at the Dukatar pedocomplex S5 Studies focused solely on instrumental colour determination of palaeosols are relatively scarce in Serbia. The results presented in this study represent the first palaeoclimatic/palaeoecological reconstruction based on colorimetric values of the strongly rubi- fied interglacial S5 pedocomplex, formed in the time frame be- tween MIS 13-15 at the Dukatar site. The observed changes, from yellowish-brown to reddish- brown, along with the corresponding indices, most likely indicate variations between haematite and goethite concentrations (JI et al., 2004; SCHEINOST & SCHWERTMANN, 1999). Studies undertaken in China indicate that on the loess plateaus, where soils are generally well drained, haematite and goethite concen- tration is closely associated with variations in soil temperature and precipitation/humidity. Magnetic susceptibility measurement of palaeosols from the Chinese Loess Plateau (CLP) indicate pa- laeo-precipitation variability between 590-1120 mm-1 (SARTORI et al., 2005). Comparable palaeo-precipitation data for the South- eastern and Central Pannonian Basin indicates variations of 580 to 650 mm a−1 over a similar timescale (PANAIOTU et al., 2001; BRADÁK et al., 2011). Palaeoclimatic forcing affects the production of pedogeni- cally formed haematite and goethite. KÄMPF & SCHWERT- MANN (1983) assert that pedogenic ferrihydrite conversion to haematite is initiated by dehydration in the soil, increasing arid- ity and/or a sustained increase in temperature. The same research states that unlike haematite, goethite production from ferrihydrite does not display a high degree of sensitivity to temperature change, with its formation being initiated under more humid con- ditions. Therefore, haematite prevalence in relation to goethite rises with an increase in mean annual air temperature, while goe- thite prevalence is associated with increased soil moisture. Based on these characteristics, the utilization of SDI along with mag- netic measurements for detection of iron-bearing minerals can serve as a useful tool in LPS research. Increases in the haematite/ goethite ratio (Hm/Gt) indicate arid conditions, while lower val- ues are an indicator of more humid conditions (LIU et al., 2006). BUGGLE et al. (2014) in their study on Quaternary climate change in LPSs from South-East Europe, presented results ob- tained by using diffuse reflectance spectroscopy (DRS) and colo- rimetric proxies. These were based on predetermined values from the Munsell Soil Colour Chart such as RI1 and RR. The RI1 and the RR rely upon the positive relationship between chroma, and the chroma to lightness ratio, to the total iron oxide content, as well as the power of hue to discriminate between the relative proportions of haematite and goethite (HURST, 1977; BUGGLE et al., 2014). The aforementioned colour proxies in correlation with the χlf values were used to determine the variation in the presence of ferromagnetic minerals in the LPSs of Batajnica and Stari Slankamen in Serbia, and Mircea Vodă in Romania, as well as the initial results for this study. Due to a general lack of absolute dating of LPSs, chro- nostratigraphic models based on the correlation of χlf values on a regional scale, along with deep-sea sediment records, are crucial for the reconstruction of the palaeoclimatic change dynamics (BUGGLE et al., 2009). The stratigraphic model developed for the Vojvodina loess, while based upon the Chinese loess-strati- graphic system, actually resulted from the systematically con- ducted research of MARKOVIĆ et al. (2004a, b, 2005, 2006, 2007, 2008, 2009, 2011, 2012, 2014, 2015, 2021a). This strati- graphic model, enabled loess units and pedocomplexes to be linked with the corresponding MIS for the first time. In general, cyclic fluctuations that occur in the χlf signal correlate well with variations of oxygen isotopes in deep-sea sediments. European stratigraphic members from S0 to S5 are in good agreement with G eo lo gi a C ro at ic a Geologia Croatica 76/280 stratigraphic members from the CLP (LU et al., 1999), corre- sponding to the MIS 1, 5, 7, 9 and 11 interglacials. The thick S5 pedocomplex at Dukatar correlates well with MIS 13 to 15 (MARKOVIĆ et al., 2015). According to its pa- laeopedological characteristics, this stratigraphic unit must have formed over a longer period of warm and humid climate as it dis- plays elements of more advanced pedogenetic modification than can be found in more recent soils from the region (BRONGER et al., 1998; BRONGER, 2003). The results of BUGGLE et al., (2014) also support this assertion. Samples from the Dukatar S5 pedocomplex display a general increase in the χlf value. This pedogenically conditioned increase in magnetic susceptibility indicates the formation of ferromag- netic minerals. Consequently, more investigation is required to clarify pedogenic and palaeoclimatic background of the investi- gated pedocomplex. A similar increase in S5 pedocomplex χlf values was also observed in the Batajnica (MARKOVIĆ et al., 2009) and Stari Slankamen LPSs (MARKOVIĆ et al., 2011). In Fig. 7, all of the parameters and indices from this study are plotted with depth. Sampled at a resolution of 3 cm, the ob- tained χlf values display variations with increasing depth: from the contact between the basal part of loess unit L5 and the upper part of the A horizon, through to the contact between the bottom part of the Ck horizon and loessoid material at the base of the ex- posed part of the profile. The SDIs RI and MI, as well as L, a*, and b* colorimetric indicators co vary with the χlf values. Mea- sured values of a* chromaticity could be considered as an indi- cator of the presence of pedogenetic haematite. In addition, these values show a strong correlation with the values of the RI1 (R2=0.91, not shown) at the significance level p<0.01. 4.2. Possible relationship between colorimtric indices and ferromagnetic composition 4.2.1. The haematite/goethite model problem Iron minerals such as haematite and goethite are significant con- stituents of most soils, and serve as proxy records of climatic and ecological changes. As has been previously asserted, warm and dry conditions are favourable for the formation of haematite, and lead to its separation from ferrihydrite (Fh), while humid condi- tions are more favourable for the formation of goethite, and lead to its formation by direct excretion from any Fe source (JI et al., 2004). Soils containing only, or almost exclusively goethite, have yellow colour, whereas in soils rich in haematite the colour is usu- ally reddish, due to the fact that red haematite successfully masks yellow goethite (SCHEINOST & SCHWERTMANN, 1999). Un- fortunately, the haematite/goethite model in loess (e.g. GUO et al., 2009; LUKIĆ et al., 2014), is not uniformly applicable in all LPSs across Europe. Contrary to this model, ZEEDEN et al. (2015) found haematite present in a loess unit at Krems-Wacht- berg LPS. Furthermore, at the Paks LPS the well developed MB and PD2 palaeosols contain significant amounts of pedogenic magnetite and maghemite, but no haematite (BRADAK et al., 2019). Haematite was found both in loess and palaeosols in Czech Republic (OCHES & BANERJEE, 1996). However, in Serbia, the study of LAAG et al. (2021) suggested that tracking the mag- netite/maghemite content on its own may not provide a clear enough interpretation of pedogenic processes during the inter- glacials. They advocated the use of colorimetric data, as it pro- vides direct insight into haematite and goethite content. The same authors revealed that between MIS2 and MIS11 palaeosols con- tain relatively higher amounts of haematite compared to goethite. 4.2.2. The significance of rubified fossil soils in LPSs As pointed out by numerous international studies colour varia- tions have been used as a proxy for mineral concentrations in LPSs (e.g. CHEN et al., 2002; VIDIĆ et al., 2004; LUKIĆ et al., 2014; BUGGLE et al., 2014). The reddening of palaeosols is shown to depend of presence of the amount of ferromagnetic min- erals such as maghaemite (SCHWERTMANN, 1993) and hae- matite (CHEN et al., 2002; BUGGLE et al., 2014). LUKIĆ et al. (2014) and BUGGLE et al. (2014) pointed out that a* chromatici ty values can be considered to be a proxy for relative abundance of pedogenic haematite similar to RI1 and RR. This pedogenic hae- matite may also be responsible for the darkening of fossil soils (as shown on CLP by CHEN et al. (2002)) affecting the MI and the L chromaticity values, regardless of soil organic matter con- tent. VIDIĆ and LOBNIK (1997) emphasized that the decrease in MI depends on time, nature and quality of soil organic matter. There are a few cases where the MI obtained in this study does not fully correspond with the other quantified indices at some depths, such as the slight increase in MI in weakly rubified Bw horizon. This feature may be attributed to soil components other than organic matter such as Mn-mottles. Based on these initial findings, soil colour indices represent additional sensitive pa- leoenvironmental proxies, which, therefore, should be more widely applied in similar investigations including rock magnetic Figure 7. All of the calculated parameters and indices for Dukatar LPS. Figure 8. The linear function fitting the Redness Index (RI1) and low-frequency magnetic susceptibility measurements (χlf) of Dukatar LPS. G eologia C roatica 81 properties. Significant correlation (p<0.01) is observed between RI and MI, wherein higher RI values indicate more intensive red colour, while higher MI values indicate its darkening (R2=0.76, not shown). The Dukatar pedocomplex is characterized by the significant presence of carbonates which are excreted from higher horizons and deposited in the central and basal part of the pa- laeosol S5. The light, whitish colour of the carbonates can be reg- istered through the L values and melanization index (MI) values. This can be particularly well noted in values observed at a depth between 180 and 195 cm from the upper part of the examined section. As mentioned in the previous chapter, combined utiliza- tion of parameters for determining rock magnetism and colori- metric studies has a rather good perspective (LUKIĆ et al., 2014). 4.2.3. Interpretation of the Dukatar pedocomplex colorimetric data Figure 8. illustrates the significance (p<0.01) of the positive re- lationship between χlf and RI1 (R2=0.88); indicating that during the Middle Pleistocene, throughout MIS 13-15 (when pedocom- plex S5 was formed) conditions were, most likely, favourable for haematite formation. The regression function for the investigated section is y=0.001x+0.032. BUGGLE et al. (2014) noted that the haematite/goethite value ratio (Hm/(Hm + Gt)) increased be- tween the S4 and S5 palaeosols. These results are somewhat in accordance with the presented initial palaeopedological and geo- chemical indicators which possibly suggest more intensive devel- opment and the effect of eluvial processes in older pedocomplexes (BUGGLE et al., 2013). The soil redness index (RI1) of the Du- katar S5 pedocomplex corresponds to a large extent with the de- tected variations recorded in the χlf. Therefore, RI1 can be used as a convenient alternative to magnetic susceptibility for genera- ting paleoclimatic data in the loess province of Vojvodina (north Serbia). The haematite content appears to be higher in the upper part of the examined section, with the highest values occurring be- tween 146 cm and 180 cm. These findings are confirmed by the values of the soil redness rating index (RR). This is in direct con- trast to RI1, which, after the calcified part of the profile, displays a continuous decrease in values towards the bottom part of the section. With regard to haematite prevalence, the RR index be- tween 171 cm to 288 cm in the upper part of the profile records values equal to 0. The observed differences between these proxies can be ex- plained by their nature, which uses distinct variables upon which differences in soil redness and haematite content are expressed. It should be noted that the RR index represents a value which, in a semi-quantitative manner, indicates possible haematite content in the soil. The RR index therefore represents a somewhat more conservative method, as the threshold hue used is 7.5YR or greater. In contrast, the RI1 is not limited in terms of colour hue, and therefore can indicate in a more sensitive manner lesser con- centrations of ferromagnetic minerals that cannot produce colour values of at least 7.5YR. Based upon the soil colour index values recorded at the Dukatar section, the use of the RI1 index is deemed to be a more appropriate metric, since this index is more sensitive to lithological changes and soil development intensity (e.g. MARKOVIĆ et al., 2009). During future studies of this pedo- complex, diffuse reflectance spectroscopy (DRS) should be used in order to determine haematite content in the most precise way possible, as well as haematite and goethite ratios in this pedo- complex. Based on findings from BUGGLE et al. (2014), which refer to other S5 pedocomplexes in north Serbia, possible palaeoenvi- ronmental conditions could be noted: in common with the Ba- tajnica and Stari Slankamen LPSs, the investigated pedocomplex at Dukatar was formed during the MIS 13-15. It is understood that these stages were characterized by higher summer tempera- tures and/or longer drought periods. Slightly higher temperatures during the winter months (>0 ºC) could serve as an initial expla- nation for the increased production of haematite in the S5 pedo- complex, since the presence of organic matter hinders its forma- tion. In addition, a marked increase in magnetic particle concentrations in the Dukatar S5 pedocomplex can be seen on the χlf plot. Likewise, it is worth mentioning that the S5 pedocom- plex in a chronostratigraphic context corresponds well with the S5 lithological unit on the CLP (Fig. 9). The results of χlf for the Dukatar pedocomplex and Luochuan S5 on the CLP (LU et al., 1999) are presented and highlighted on the Figure 9. Namely, this highly developed pedocomplex (corresponding to the MIS 13-15) is used as a stratigraphic marker in central China because of its dark brown-reddish colour (CHEN et al., 2002; GUO et al., 2013; HAN et al., 1998, LU et al., 1999). Based on the studied pedoge- netic characteristics, HAN et al. (1998) assert that the S5 pedo- complex was formed during a period characterized by relatively stable and drier climatic conditions in relation to the climatic con- ditions during the Holocene. According to YIN and GUO (2008), the existing climatic conditions were spurred by the strengthen- ing of summer monsoon circulation in the northern hemisphere. In addition, the measured values of χlf presented in this study cor- relate well with research results of other LPSs both in the region and on the Eurasian continent (BUGGLE et al., 2009; MARKO- VIĆ et al., 2009, 2012b, 2015). In general, higher χlf values in the S5 pedocomplex on the Dukatar profile also indicate a possible reduction in the excess moisture in the soil, which is associated with slightly higher tem- peratures during the winter months; this results in a reduction in the precipitation/evaporation ratio. The aforementioned interpre- tations suggest climatic conditions similar to those in the Medi- terranean, prevailed in the loess region of both the Lower and Middle Danube Basin during MIS 13-15 (BUGGLE et al., 2014). The prevailing influence of the Mediterranean climate on soil formation persists for around 180 km southwards of the investi- gated Dukatar site, as described at the Stalać LPS. Here, a well rubified cambisol was able to form during both MIS 9 and MIS 7, when the Middle Danube basin was already experiencing colder, continental climate conditions (OBREHT et al., 2016). 4.3. Comparative views on interglacial palaeoclimatic and palaeoecological conditions during the formation of S1 and S5 pedocomplexes in northern Serbia The same methodological approach undertaken for this study was also adopted by LUKIĆ et al. (2014) at the Orlovat LPS (Figure 1). As the Orlovat and Dukatar are the only LPSs investigated from a joint rock magnetic and colorimetric perspective in the north part of Serbia, the comparison between the sections is pre- sented here. Well defined differences in interglacial palaeoclimatic and palaeoecological conditions become apparent when studying pa- laeosol S1 at the Orlovat-brickyard, and pedocomplex S5 Dukatar presented in this study (Fig. 10). The S1 palaeosol is a strongly developed soil of approximately 400 cm thickness formed during the MIS 5 period. The aforementioned thickness indicates higher accumulation rates during the final part of MIS 5 and MIS 4 (MARKOVIĆ et al., 2014). The S1 palaeosol has significantly higher χlf values ranging from 35 to 70 × 10-8m3kg-1 in contrast G eo lo gi a C ro at ic a Geologia Croatica 76/282 to the L1 loess unit which shows small variations in χlf values (cca. 20 × 10-8m3kg-1). These values gradually increase towards the base of S1, possibly indicating pedogenetic formation of su- perparamagnetic and single-domain ferromagnetics (LUKIĆ et al., 2014). The χlf plot shows high similarity with the correspond- ing palaeosols at other sites in Vojvodina, the Carpathian region and in the CLP (MARKOVIĆ et al., 2009, 2012, 2014). The mea- sured values of a* chromaticity generally follow the χlf values, RI and MI. Warmer climatic conditions during the lower Eemian are presumably reflected in higher a* chromaticity, RI and MI values at the top of the S1 lithological member. Humid episodes, which provided sufficient moisture in the soil, occurred during the late spring and early summer. The existence of a highly dif- ferentiated Ah horizon, characterized by the accumulation of or- Figure 10. The low-frequency magnetic susceptibility measurements (χlf), a* chromaticity, Rubification index (RI) and Melanization index (MI) for A) Orlovat LPS (LUKIĆ et al., 2014); and B) Dukatar LPS (covered by this study). Figure 9. The normalized magnetic susceptibility of some of the oldest LPS in Serbia: Batajnica (MARKOVIĆ et al., 2009), TLP (MARKOVIĆ et al., 2012), Stari Slanka- men (MARKOVIĆ et al., 2011), and Orlovat (LUKIĆ et al., 2014) presented on a depth scale and compared to Lingtai/Zhaojiachuan (SUN et al., 2006) and LR04 (LISIECKI & RAYMO, 2005) benthic oxygen isotope record presenting the time scale. MIS stages as well as magnetic reversals are presented. Note the data from Dukatar S5 next to the TLP LPS. G eologia C roatica 83 ganic matter supports this conclusion (LUKIĆ et al., 2014). In general, this chernozem palaeosol indicates a weakening of pedo- genesis from interglacial conditions towards the early glacial pe- riod (MARKOVIĆ et al., 2012, 2014). The χlf record of Dukatar S5 shows variations due to the pre- sence of carbonate concretions and hydromorphic features. Based upon the data presented in Fig. 10, it can be observed that the concentration of ferromagnetic minerals is significantly higher in the Dukatar S5 pedocomplex compared to the S1 in Orlovat. Significantly higher rubification index values in the S5 pedocom- plex indicate the dominance of different pedogenetic conditions to those in which the S1 steppe palaeosol was formed. This find- ing, to a large extent, coincides with other loess-palaeosol studies in which rubification index values were presented (e.g. VIDIĆ et al., 2004; MARKOVIĆ et al., 2009, 2011; BUGGLE et al., 2014); all indicate a strongly expressed palaeoclimatic transition from sub-Mediterranean to arid continental climatic conditions in this part of Europe during the last five glacial-interglacial cycles. The variations in the RI and a* values presented in this study provide further evidence of palaeoclimatic evolution during different phases of Pleistocene. Compared to the S1 palaeosol, MI has sig- nificantly lower values due to the presence of different prevailing palaeoclimatic and palaeoecological conditions during the ob- served interglacial periods. The study of LAAG et al. (2021) as- serts that that the general climatic trend was characterized by progressively reduced precipitation and air temperature over the last 430 kyr (covering the formation of palaeosols S4 to S1), grad- ually leading to more arid climate conditions (from sections about 50 km south of the Dukatar LPS). 5. CONCLUSIONS The aim of this paper is to present the first reconstruction of pa- laeoenvironmental and climate variability for the MIS13-15 S5 Dukatar pedocomplex based on high resolution results of the SDI, colour indicators and proxies, supported by χlf. The observed dif- ference in variations of the soil redness index RI1 appear to be influenced by haematite content in this rubified pedocomplex, whereas it is significantly reduced in the soil redness rating index (RR). Based on the presented values of the soil colour index for the Dukatar section, it could be noted that the use of the RI1 in- dex is deemed more appropriate, since this index appears more sensitive to lithological changes and soil development intensity, while the RR index is considered a somewhat more conservative method. The RI1 can therefore be considered a more useful pal- aeoenvironmental proxy, which should be more widely applied in similar multiproxy investigations as a possible tracer of iron- bearing minerals. Moreover, it should be emphasized the multi- proxy approach set out in this study has proved a very favourable methodology for assessing palaeoclimatic and palaeoenviron- mental conditions from palaeosols. While the studied indices have confirmed the assertions of earlier studies, regarding the climatic transition from humid Mediterranean towards the more arid continental one, more detailed investigation is required to verify precise ferromagnetic contents, their ratio, and correlation to colorimetric perspectives, in order to clarify pedogenic and paleoclimatic context. ACKNOWLEDGEMENT The authors gratefully acknowledge the financial support of the Ministry of Science, Technological Development and Innovation of the Republic of Serbia (Grant No. 451-03-47/2023-01/200125). This research was partially funded by the Provincial Secretariat for Higher Education and Scientific Research of Autonomous Province of Vojvodina, Republic of Serbia (Grant No. 142-451- 2655/2021-01/2). The authors are grateful to Đorđije VASILJEVIĆ and Miroslav VUJIČIĆ for their help during the field work. We confirm that all the authors made an equal contribution to the study and its development. Finally, the authors are grateful to the anonymous reviewers whose comments and suggestions greatly improved the manuscript. REFERENCES BANAK, A., MANDIC, O., SPROVIERI, M., LIRER, F. & PAVELIĆ, D. (2016): Stable isotope data from loess malacofauna: Evidence for climate changes in the Pannon- ian Basin during the Late Pleistocene.– Quatern. Int., 415, 15–24. doi: 10.1016/j. quaint.2015.10.102 BANAK, A., PAVELIĆ, D., KOVAČIĆ, M. & MANDIC, O. (2013): Sedimentary char- acteristics and source of loess in Baranja (Eastern Croatia).– Aeolian Res., 11, 129–139. doi: 10.1016/j.aeolia.2013.08.002 BARRÓN, V. & TORRENT, J. (1986): Use of the Kubelka – Munk theory to study the influence of iron oxides on soil colour.– J. Soil Sci., 37/4, 499–510. doi: 10.1111/ j.1365-2389.1986.tb00382.x BASARIN, B., BUGGLE, B., HAMBACH, U., MARKOVIĆ, S.B., DHAND, K.O. H., KOVAČEVIĆ, A., ... & LUKIĆ, T. (2014): Time-scale and astronomical forcing of Serbian loess–paleosol sequences.– Global Planet. Change, 122, 89–106. doi: 10.1016/j.gloplacha.2014.08.007 BERGER, G.W., PILLANS, B.J. & PALMER, A.S. (1992): Dating loess up to 800 ka by thermoluminescence.– Geology, 20/5, 403–406. doi: 10.1130/0091-7613(1992) 020<0403:DLUTKB>2.3.CO;2 BOESKEN, J., KLASEN, N., ZEEDEN, C., OBREHT, I., MARKOVIĆ, S.B., HAM- BACH, U. & LEHMKUHL, F. (2017): New luminescence‐based geochronology framing the last two glacial cycles at the southern limit of European Pleistocene loess in Stalać (Serbia).– Geochronometria, 44/1, 150–161. doi: 10.1515/geo- chr-2015-0062 BRADÁK, B., SETO, Y., CSONKA, D., VÉGH, T., & SZEBERÉNYI, J. (2019): The hematite–goethite enhancement model of loess and an ‘irregular’case from Paks, Hungary.– J. of Quaternary Sci. 34/4-5, 299-308. BRADÁK, B., THAMÓ-BOZSÓ, E., KOVÁCS, J., MÁRTON, E., CSILLAG, G., & HORVÁTH, E. (2011): Characteristics of Pleistocene climate cycles identified in Cérna Valley loess–paleosol section (Vértesacsa, Hungary).– Quaternary Interna- tional, 234/1–2, 86–97. doi: 10.1016/j.quaint.2010.05.002 BRONGER, A. (1976): Zur quartären Klima- und Landschaftsentwicklung des karpaten- beckens auf (paläo) pedologischer Grundlage.– Habilitationschrift. Kieler Geogra- phischen Schriften, 45, Kiel, 1–269. BRONGER, A. (2003): Correlation of loess–paleosol sequences in East and Central Asia with SE Central Europe: towards a continental Quaternary pedostratigraphy and paleoclimatic history.– Quatern. Int., 106, 11–31. doi: 10.1016/S1040- 6182(02)00159-3 BRONGER, A., WINTER, R. & SEDOV, S. (1998): Weathering and clay mineral forma- tion in two Holocene soils and in buried paleosols in Tadjikistan: towards a Quater- nary paleoclimatic record in Central Asia.–Catena, 34/1–2, 19–34. doi: 10.1016/ S0341-8162(98)00079-4 BUGGLE, B., HAMBACH, U., GLASER, B., GERASIMENKO, N., MARKOVIĆ, S. B., GLASER, I. & ZÖLLER, L. (2009): Stratigraphy and spatial and temporal pa- leoclimatic trends in East European loess paleosol sequences.– Quatern. Int., 196, 86–106. doi: 10.1016/j.quaint.2008.07.013 BUGGLE, B., HAMBACH, U., KEHL, M., ZÖLLER, L., MARKOVIĆ, S.B. & GLA- SER, B. (2013): The progressive evolution of a continental climate in SE-Central European lowlands during the Middle Pleistocene recorded in loess paleosol se- quences.– Geology, 41, 771–774. doi: 10.1130/G34198.1 BUGGLE, B., HAMBACH, U., MÜLLER, K., ZÖLLER, L., MARKOVIĆ, S.B. & GLA- SER, B. (2014): Iron mineralogical proxies and Quaternary climate change in SE- European loess–paleosol sequences.– Catena, 117, 4–22. doi: 10.1016/j.cate- na.2013.06.012 CHEN, J., JI, J., BALSAM, W., CHEN, Y., LIU, L. & AN, Z. (2002): Characterization of the Chinese loess–paleosol stratigraphy by whiteness measurement.– Palaeogeogr. Palaeocl., 183/3–4, 287–297. doi: 10.1016/S0031-0182(02)00246-8 FENN, K., MILLAR, I.L., DURCAN, J. A., THOMAS, D.S., BANAK, A., MARKOVIĆ, S.B., ... & STEVENS, T. (2022): The provenance of Danubian loess.– Earth-Science Reviews, 103920. doi: 10.1016/j.earscirev.2022.103920 GALOVIĆ, L. (2014): Geochemical archive in the three loess/paleosol sections in the Eastern Croatia: Zmajevac I, Zmajevac and Erdut.– Aeolian Res., 15, 113–132. doi: 10.1016/j.aeolia.2014.07.004 GALOVIĆ, L. (2016): Sedimentological and mineralogical characteristics of the Pleisto- cene loess/paleosol sections in the Eastern Croatia.– Aeolian Res., 20, 7–23. doi: 10.1016/j.aeolia.2015.10.007 G eo lo gi a C ro at ic a Geologia Croatica 76/284 GALOVIĆ, L. & PEH, Z. (2016): Mineralogical discrimination of the pleistocene loess/ paleosol sections in Srijem and Baranja, Croatia.– Aeolian Res., 21, 151–162. doi: 10.1016/j.aeolia.2016.04.006 GUO ZT, BERGER A, YIN, QZ & QIN, L. (2009). Strong asymmetry of hemispheric climates during MIS‐13 inferred from correlating China loess and Antarctica ice records.– Clim. Past. 5, 21–31. doi: 10.5194/cp-5-21-2009 GUO, X., LIU, X., LI, P., LÜ, B., GUO, H., CHEN, Q., ... & MA, M. (2013): The mag- netic mechanism of paleosol S5 in the Baoji section of the southern Chinese Loess Plateau.– Quatern. Int., 306, 129-136. doi: 10.1016/j.quaint.2013.02.033 HAN, J., FYFE, W. S., & LONGSTAFFE, F. J. (1998): Climatic implications of the S5 paleosol complex on the southernmost Chinese Loess Plateau.– Quaternary Res., 50/1, 21–33. doi: 10.1006/qres.1998.1976 HARDEN, J.W. (1982): A quantitative index of soil development from field descriptions: Examples from a chronosequence in central California.– Geoderma, 28/1, 1–28. doi: 10.1016/0016-7061(82)90037-4 HARDEN, J.W. & TAYLOR, E.M. (1983): A quantitative comparison of soil develop- ment in four climatic regimes.– Quaternary Res., 20/3, 342–359. doi: 10.1016/0033- 5894(83)90017-0 HELLER, F. & EVANS, M.E. (1995): Loess magnetism.– Rev. of Geophysics, 33/2, 211–240. doi: 10.1029/95RG00579 HÜNTER LAB (2001): Hünter Lab versus CIE 1976 La*b* color scale. Insight on color, 1–15. HURST, V.J. (1977): Visual estimation of iron in saprolite.– Geol. Soc. Am. Bull., 88/2, 174–176. doi: 10.1130/0016-7606(1977)88%3C174:VEOIIS%3E2.0.CO;2 JI, J., CHEN, J., BALSAM, W., LU, H., SUN, Y. & XU, H. (2004): High resolution he- matite/goethite records from Chinese loess sequences for the last glacial‐intergla- cial cycle: Rapid climatic response of the East Asian Monsoon to the tropical Pa- cific.– Geophys. Res. Lett., 31/3. doi: 10.1029/2003GL018975 KÄMPF, N. & SCHWERTMANN, U. (1983): Goethite and hematite in a climosequence in southern Brazil and their application in classification of kaolinitic soils.– Geo- derma, 29/1, 27–39. doi: 10.1016/0016-7061(83)90028-9 KUKLA, G. (1987): Loess stratigraphy in central China.– Quaternary Sci. Rev., 6/3–4, 191–219. doi: 10.1016/0277-3791(87)90004-7 LAAG, C., HAMBACH, U., ZEEDEN, C., LAGROIX, F., GUYODO, Y., VERES, D., ... & MARKOVIĆ, S. (2021): A Detailed Paleoclimate Proxy Record for the Mid- dle Danube Basin Over the Last 430 kyr: A Rock Magnetic and Colorimetric Study of the Zemun Loess-Paleosol Sequence.– Front. Earth Sci., 9, 600086. doi: 10.3389/ feart.2021.600086 LIU, Q., BLOEMENDAL, J., TORRENT, J. & DENG, C. (2006): Contrasting behavior of hematite and goethite within paleosol S5 of the Luochuan profile, Chinese Loess Plateau.– Geophys. Res. Lett., 33/20. doi: 10.1029/2006GL027172 LU, H., LIU, X., ZHANG, F., AN, Z. & DODSON, J. (1999): Astronomical calibration of loess–paleosol deposits at Luochuan, central Chinese Loess Plateau.– Palaeoge- ography, Palaeoclimatology, Palaeoecology, 154/3, 237–246. doi: 10.1016/S0031- 0182(99)00113-3 LUDWIG, P., GAVRILOV, M. B., RADAKOVIĆ, M. G. & MARKOVIĆ, S. B. (2021): Malaco temperature reconstructions and numerical simulation of environmental conditions in the southeastern Carpathian Basin during the Last Glacial Maximum.– J. Quaternary Sci., 36/8, 1426–1435. doi: 10/1002/jqs.3318 LUKIĆ, T., BASARIN, B., BUGGLE, B., MARKOVIĆ, S.B., TOMOVIĆ, V.M., RALJIČ, J.P., ... & GAVRILOV, M. B. (2014): A joined rock magnetic and colori- metric perspective on the Late Pleistocene climate of Orlovat loess site (Northern Serbia).– Quatern. Int., 334, 179–188. doi: 10.1016/j.quaint.2014.03.042 MARKOVIĆ, S.B., BOKHORST, M.P., VANDENBERGHE, J., MCCOY, W.D., OCHES, E.A., HAMBACH, U., ... & MACHALETT, B. (2008): Late Pleistocene loess‐paleosol sequences in the Vojvodina region, north Serbia.– J. Quaternary Sci.: Published for the Quaternary Research Association, 23/1, 73–84. doi: 10/1002. jqs.1124 MARKOVIĆ, S.B., GAVRILOV, B.M., PERIĆ, Z. & RADAKOVIĆ, M. (2021a): Лесна палеоземљишта у Србији - Зашто треба чувати земљишне ресурсе, лекције из прошлости. [Loess-paleosol sequences in Serbia – Why should land resources be preserved, lessons from the past – In Serbian].– Sustainable use of land and water in Serbia, 197, 19. MARKOVIĆ, S.B., HAMBACH, U., CATTO, N., JOVANOVIĆ, M., BUGGLE, B., MACHALETT, B., ... & FRECHEN, M. (2009): Middle and late Pleistocene loess sequences at Batajnica, Vojvodina, Serbia.– Quatern. Int., 198/1–2, 255–266. doi: 10.1016/j.quaint.2008.12.004 MARKOVIĆ, S.B., HAMBACH, U., STEVENS, T., BASARIN, B., O’HARA-DHAND, K., GAVRILOV, M.M., ... & TEOFANOV, N. (2012a): Relating the astronomical timescale to the loess–paleosol sequences in Vojvodina, Northern Serbia.– Climate change, 65–78. doi: 10.1007/978-3-7091-0973-1_5 MARKOVIĆ, S.B., HAMBACH, U., STEVENS, T., JOVANOVIĆ, M., O’HARA- DHAND, K., BASARIN, B., ... & ZÖLLER, L. (2012b): Loess in the Vojvodina region (Northern Serbia): an essential link between European and Asian Pleistocene environments.– Neth. J. Geosci., 91/1–2, 173–188. doi: 10.1017/ S0016774600001578 MARKOVIĆ, S.B., HAMBACH, U., STEVENS, T., JOVANOVIĆ, M., O’HARA- DHAND, K., BASARIN, B., ... & ZÖLLER, L. (2013): Loess in the Vojvodina region (Northern Serbia): an essential link between European and Asian Pleistocene environments.– Netherlands Journal of Geosciences, 91/1–2, 173–188. doi: 10.1017/S0016774600001578 MARKOVIĆ, S.B., HAMBACH, U., STEVENS, T., KUKLA, G.J., HELLER, F., MC- COY, W.D., ... & ZÖLLER, L. (2011): The last million years recorded at the Stari Slankamen (Northern Serbia) loess-palaeosol sequence: revised chronostratigraphy and long-term environmental trends.– Quaternary Sci. Rev., 30/9–10, 1142–1154. doi: 10.1016/j.quascirev.2011.02.004 MARKOVIĆ, S.B., KOSTIĆ, N.S., & OCHES, E.A. (2004a): Paleosols in the Ruma loess section (Vojvodina, Serbia).– Rev. Mex. Cienc. Geol., 21/1, 79–87. MARKOVIĆ, S.B., MCCOY, W.D., OCHES, E.A., SAVIĆ, S., GAUDENYI, T., JOVANOVIĆ, M., STEVENS, T., WALTHER, R., IVANIŠEVIĆ, P., & GALOVIĆ, Z. (2005): Paleoclimate record in the Late Pleistocene loess–paleosol sequence at Petrovaradin Brickyard (Vojvodina, Serbia).– Geol. Carpath., 56, 545–552. MARKOVIC, S.B., OCHES, E.A., GAUDENYI, T., JOVANOVIC, M., HAMBACH, U., ZÖLLER, L., & SÜMEGI, P. (2004b): Palaeoclimate record in the Late Pleis- tocene loess-paleosol sequence at Miseluk (Vojvodina, Serbia).– Quaternaire, 15/4, 361–368. doi: 10/3406/quate.2004.1781 MARKOVIĆ, S.B., OCHES, E.A., MCCOY, W.D., FRECHEN, M., & GAUDENYI, T. (2007): Malacological and sedimentological evidence for “warm” glacial climate from the Irig loess sequence, Vojvodina, Serbia.– Geochemistry, Geophysics, Geo- systems, 8/9. doi: 10.1029/2006GC001565 MARKOVIĆ, S. B., OCHES, E. A., PERIĆ, Z. M., GAUDENYI, T., JOVANOVIĆ, M., SIPOS, G., ... & GAVRILOV, M. B. (2021b): The Požarevac loess–paleosol se- quence: a record of increased aridity in the south‐eastern margin of the Carpathian Basin during the last 350 ka. –J. Quaternary Sci., 36/8, 1436–1447. doi: 10.1002/ jqs.3327 MARKOVIĆ, S.B., OCHES, E., SÜMEGI, P., JOVANOVIĆ, M. & GAUDENYI, T. (2006): An introduction to the Middle and Upper Pleistocene loess–paleosol se- quence at Ruma brickyard, Vojvodina, Serbia.– Quatern. Int., 149/1, 80–86. doi: 10.1016/j.quaint.2005.11.020 MARKOVIĆ, S.B., STEVENS, T., KUKLA, G.J., HAMBACH, U., FITZSIMMONS, K.E., GIBBARD, P., ... & SVIRČEV, Z. (2015): Danube loess stratigraphy – Towards a pan-European loess stratigraphic model.– Earth-Science Reviews, 148, 228–258. doi: 10.1016/j.earscirev.2015.06.005 MARKOVIĆ, S. B., TIMAR‐GABOR, A., STEVENS, T., HAMBACH, U., POPOV, D., TOMIĆ, N., ... & GAVRILOV, M. B. (2014): Environmental dynamics and lumi- nescence chronology from the Orlovat loess–palaeosol sequence (Vojvodina, north- ern Serbia).– J. Quaternary Sci., 29/2, 189–199. doi: 10.1002/jqs.2693 MARKOVIĆ-MARJANOVIĆ, J. (1950): Contribution to the geological structure of Ti- tel loess plateau. [Prilog za geološku građu Titelskog brega – In Serbian].– Zbornik radova Geološkog instituta SANU, Beograd, 3/1, 91–121. MARKOVIĆ-MARJANOVIĆ, J. (1965): A review of the loess problems of Yugoslavia in the time of Jovan Cvijić and today. [Osvrt na poznavanje lesnih problema Jugo- slavije u doba Jovana Cvijića i danas – In Serbian].– Bulletin of the Serbian Geo- graphical society, 45, 2, 99–113. MARSIGLI, A.F. (1726): Danubius Pannonico Mysicus. P. Grosse, Chr. Alberts, P de Hoodt Haga, Herm. Uytwert and Franc Changuion Amsterdamus. NAMIER, N., GAO, X., HAO, Q., MARKOVIĆ, S. B., FU, Y., SONG, Y., ... & GUO, Z. (2021): Mineral magnetic properties of loess–paleosol couplets of northern Ser- bia over the last 1.0 Ma.– Quaternary Res., 103, 35–48. doi: 10.1017/qua.2021.41 NETTLETON, W.D., OLSON, C.G. & WYSOCKI, D.A. (2000): Paleosol classification: problems and solutions.– Catena, 41/1–3, 61–92. doi: 10.1016/S0341- 8162(00)00109-0 OBREHT, I., BUGGLE, B., CATTO, N., MARKOVIČ, S.B., BÖSEl, S., VANDEN- BERGHE, D.A., ... & JOVIĆ, G. (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.– Quatern. Int., 334, 10–19. doi: 10.1016/j.quaint.2013.05.037 OBREHT, I., ZEEDEN, C., HAMBACH, U., VERES, D., MARKOVIĆ, S.B., BÖSKEN, J., ... & LEHMKUHL, F. (2016): Tracing the influence of Mediterranean climate on Southeastern Europe during the past 350,000 years.– Sci. Rep-UK, 6/1, 1–10. doi: 10.1038/srep36334 OCHES, E.A. & BANERJEE, S.K. (1996): Rock‐magnetic proxies of climate change from loess‐paleosol sediments of the Czech Republic.– Studia Geophysica et Ge- odætica. 40, 287–300. doi: 10.1007/BF02300744 PANAIOTU, C.G., PANAIOTU, E.C., GRAMA, A., & NECULA, C., (2001): Paleocli- matic record from a loess–paleosol profile in southeastern Romania.– Physics and Chemistry of the Earth (A) 26, 893–898. doi: 10.1016/S1464-1895(01)00138-7 PAVIĆ, D. (2006): Potamological characteristics of the Tisza in Serbia and predispositions for the development of nautical tourism. [Potamološke karakteristike Tise u Srbiji i predispozicije za razvoj nautičkog turizma – In Serbian].– PhD thesis, Faculty of Science, University of Novi Sad, 343 p. G eologia C roatica 85 PERIĆ, Z., ADOLPHI, E.L., STEVENS, T., ÚJVÁRI, G., ZEEDEN, C., BUYLAERT, J.P., ... & GAVRILOV, M.B. (2019): Quartz OSL dating of late quaternary Chinese and Serbian loess: a cross Eurasian comparison of dust mass accumulation rates.– Quatern. Int., 502, 30–44. doi: 10.1016/j.quaint.2018.01.010 POPOVICS, A.B. (1876): Jelentés a Fruscagora hegysében tett geológiai gyüjtés és ku- tatásról – II Üledékes közetek.– Földtani Közlöny, Budapest, 286–301. PÖTTER, S., VERES, D., BAYKAL, Y., NETT, J.J., SCHULTE, P., HAMBACH, U., & LEHMKUHL, F. (2021): Disentangling sedimentary pathways for the Pleniglacial Lower Danube loess based on geochemical signatures.– Front. Earth Sci., 9, 150. doi: 10.3389/feart.2021.600010 RADAKOVIĆ, M.G., GAVRILOV, M.B., HAMBACH, U., SCHAETZL, R.J., TOŠIĆ, I., NINKOV, J., ... & MARKOVIĆ, S.B. (2019): Quantitative relationships between climate and magnetic susceptibility of soils on the Bačka Loess Plateau (Vojvodina, Serbia).– Quatern. Int., 502, 85–94. doi: 10.1016/j.quaint.2018.04.040 ROSSEL, R.V., MINASNY, B., ROUDIER, P., & MCBRATNEY, A.B. (2006): Colour space models for soil science.– Geoderma, 133/3–4, 320–337. doi: 10.1016/j.geo- derma.2005.07.017 RUBINIĆ, V., GALOVIĆ, L., LAZAREVIĆ, B., HUSNJAK, S. & DURN, G. (2018): Pseudogleyed loess derivates–the most common soil parent materials in the Pan- nonian region of Croatia.– Quatern. Int., 494, 248–262. doi: 10.1016/j.quaint. 2017.06.044 ŠARIĆ, J. (2008): Paleolithic and mesolithic finds from profile of the Zemun loess.– Stari- nar, 58, 9–27. SARTORI, M., EVANS, M.E., HELLER, F., TSATSKIN, A. & HAN, J. M. (2005): The last glacial/interglacial cycle at two sites in the Chinese Loess Plateau: mineral mag- netic, grain-size and 10Be measurements and estimates of palaeoprecipitation.– Pa- laeogeography, Palaeoclimatology, Palaeoecology, 222/1–2, 145–160. doi: 10.1016/j.palaeo.2005.03.013 SCHAETZL, R. & ANDERSON, S. (2005): Soil genesis and profile differentiation.– In: Soils: genesis and geomorphology. New York: Cambridge, 347–461. SCHATZ, A.K., ZECH, M., BUGGLE, B., GULYÁS, S., HAMBACH, U., MARKOVIĆ, S.B., ... & SCHOLTEN, T. (2011): The late Quaternary loess record of Tokaj, Hun- gary: Reconstructing palaeoenvironment, vegetation and climate using stable C and N isotopes and biomarkers.– Quatern. Int., 240/1–2, 52–61. doi: 10.1016/j. quaint.2010.10.009 SCHEINOST, A.C. & SCHWERTMANN, U. (1999): Color identification of iron oxides and hydroxysulfates: use and limitations.– Soil Sci. Soc. Am. J., 63/5, 1463–1471. doi: 10.2136/sssaj1999.6351463x SCHWERTMANN, U. (1993): Relations between iron oxides, soil color, and soil forma- tion.– Soil color, 31, 51–69. doi: 10.2136/sssaspecpub31.c4 SGRS (2012): Official Gazette of the Republic of Serbia, Decree on the proclamation of a special nature reserve “Titel loess plateau” [Službeni glasnik Rrepublike Srbije, Uredba o proglašenju specijalnog rezervata prirode “Titelski breg” – in Serbian], 56. SHETH, H.C., JOHNSON, C. P., & OLLIER, C.D. (2010): The “seven-coloured earth” of Chamarel, Mauritius.– J. Afr. Earth Sci., 57/1–2, 169–173. doi: 10.1016/j. jafrearsci.2009.07.009 SPRAFKE, T., SCHULTE, P., MEYER-HEINTZE, S., HÄNDEL, M., EINWÖGERER, T., SIMON, U., ... & TERHORST, B. (2020): Paleoenvironments from robust loess stratigraphy using high-resolution color and grain-size data of the last glacial Krems- Wachtberg record (NE Austria).– Quaternary sci. rev., 248, 106602. doi: 10.1016/j. quascirev.2020.106602 SÜMEGI, P., GULYÁS, S., MOLNÁR, D., SÜMEGI, B.P., TÖRŐCSIK, T., ALMOND, P.C., ... & KOLOSZÁR, L. (2019): Periodicities of paleoclimate variations in the first high-resolution non-orbitally tuned grain size record of the past 1 Ma from SW Hungary and regional, global correlations.– Aeolian Res., 40, 74–90. doi: 10.1016/j. aeolia.2019.07.002 SUN, Y., CLEMENS, S.C., AN, Z., & YU, Z. (2006): Astronomical timescale and pal- aeoclimatic implication of stacked 3.6-Myr monsoon records from the Chinese Loess Plateau.– Quaternary Sci. Rev., 25/1–2, 33–48. doi: 10.1016/j.quasci- rev.2005.07.005 TAYLOR, E.M. (1988): Instructions for the soil development index template: Lotus 1-2- 3 (No. USGS-OFR-88-233). US Geological Survey, Denver, CO (United States). TIMAR-GABOR, A., PANAIOTU, C., VEREȘ, D., NECULA, C. & CONSTANTIN, D. (2017): The lower Danube loess, new age constraints from luminescence dating, magnetic proxies and isochronous tephra markers.– In: Landform dynamics and evolution in Romania. Springer, Cham, 679–697. doi: 10.1007/978-3-319-32589- 7_29 TORRENT, J., SCHWERTMANN, U., FECHTER, H. & ALFEREZ, F. (1983): Quanti- tative relationships between soil color and hematite content.– Soil Sci., 136/6, 354–358. doi: 10.1097/00010694-198312000-00004 VANDENBERGHE, J. (2013): Grain size of fine-grained windblown sediment: A pow- erful proxy for process identification.– Earth-Science Reviews, 121, 18–30. doi: 10.1016/j.earscirev.2013.03.001 VASILJEVIĆ, D.A., MARKOVIĆ, S.B., HOSE, T.A., DING, Z., GUO, Z., LIU, X., ... & VUJIČIĆ, M. D. (2014): Loess–palaeosol sequences in China and Europe: com- mon values and geoconservation issues.– Catena, 117, 108–118. doi: 10.1016/j. catena.2013.06.005 VASILJEVIĆ, D.A., MARKOVIĆ, S.B., HOSE, T. A., SMALLEY, I., BASARIN, B., LAZIĆ, L. & JOVIĆ, G. (2011): The Introduction to Geoconservation of loess- palaeosol sequences in the Vojvodina region: Significant geoheritage of Serbia.– Quatern. Int., 240/1–2, 108–116. doi: 10.1016/j.quaint.2010.07.008 VIDIC, N. J. & LOBNLK, F. (1997): Rates of soil development of the chronosequence in the Ljubljana Basin, Slovenia.– Geoderma, 76/1–2, 35–64. doi: 10.1016/S0016- 7061(96)00098-5 VIDIĆ, N.J., SINGER, M.J. & VEROSUB, K.L. (2004): Duration dependence of mag- netic susceptibility enhancement in the Chinese loess–palaeosols of the past 620 ky.–Palaeogeogr., Palaeoclim. Paleoeco., 211/3–4, 271–288. doi: 10.1016/j.pal- aeo.2004.05.012 VIŠNIĆ, T., SPASOJEVIĆ, B. & VUJIČIĆ, M. (2016): The potential for geotourism de- velopment on the Srem Loess Plateau based on a preliminary geosite assessment model (GAM).– Geoheritage, 8/2, 173–180. doi: 10.1007/s12371-015-0149-9 VUJIČIĆ, M.D., VASILJEVIĆ, D.A., MARKOVIĆ, S.B., HOSE, T.A., LUKIĆ, T., HADŽIĆ, O. & JANIĆEVIĆ, S. (2011): Preliminary geosite assessment model (GAM) and its application on Fruška Gora Mountain, potential geotourism destina- tion of Serbia.– Acta geogr. Slov., 51/2, 361–376. doi: 10.3986/AGS51303 LISIECKI, L.E. & RAYMO, M.E. (2005): A Pliocene‐Pleistocene stack of 57 globally distributed benthic δ18O records.– Paleoceanography, 20/1. doi: 10.1029/2004PA001071 WACHA, L., LAAG, C., GRIZELJ, A., TSUKAMOTO, S., ZEEDEN, C., IVANIŠEVIĆ, D., ... & FRECHEN, M. (2021): High-resolution palaeoenvironmental reconstruc- tion at Zmajevac (Croatia) over the last three glacial/interglacial cycles.– Palaeoge- ogr., Palaeocl., 576, 110504. doi: 10.1016/j.palaeo.2021.110504 YIN, Q.Z. & GUO, Z.T. (2008): Strong summer monsoon during the cool MIS-13.– Clim. Past, 4/1, 29–34. doi: 10.5194/cp-4-29-2008 ZEEDEN, C., HAMBACH, U., HÄNDEL, M. (2015): Loess magnetic fabric of the Krems‐Wachtberg archaeological site.– Quatern. Int., 372, 188–194. ZHANG, W., LU, H., LI, C., DODSON, J., & MENG, X. (2017): Pollen preservation and its potential influence on paleoenvironmental reconstruction in Chinese loess de- posits.– Rev. Palaeobot. Palyno., 240, 1–10. doi: 10.1016/j.revpalbo.2017.01.002 ZYKINA, V.S., ZYKIN, V.S., VOLVAKH, A.O., RADAKOVIĆ, M.G., GAVRILOV, M.B. & MARKOVIĆ, S.B. (2020): Late Pleistocene loess-paleosol sequence at the Belovo section, south of Western Siberia, Russia: Preliminary results.– Quatern. Int., 620, 75–84. doi: 10.1016/j.quaint.2020.10.069