Abstract
Keywords
Introduction
Aeolian sediments and related landforms reflect variations in atmospheric circulation patterns on continental (Lancaster, 1990; Renssen et al., 2007) to more regional or local scale (Björck and Clemmensen, 2004; Nielsen et al., 2016a, 2016b). Compared to many temperate and subarctic parts of the world like the European Sand Belt (Isarin et al., 1997; Kasse, 1997, 2002; Koster, 1988; Tolksdorf and Kaiser, 2012), however, only scattered areas of windblown deposits are described from Fennoscandia (Lancaster et al., 2016).
Windblown deposits and related landforms occur in Finnish Lapland (Clarke and Käyhkö, 1997; Käyhkö et al., 1999; Matthews and Seppälä, 2014; Seppälä, 1971, 1981, 1995, 2004; Vliet-Lanoë et al., 1993), Swedish Lapland (Seppälä, 1972) and in central Sweden (Agrell and Hultman, 1971; Alexanderson et al., 2016; Alexanderson and Bernhardson, 2016, 2019; Alexanderson and Fabel, 2015; Bergqvist, 1981; Bernhardson et al., 2019; Bernhardson and Alexanderson, 2017, 2018; Högbom, 1913, 1923; Hörner, 1927; Stevens et al., 2022). Along the southwestern coast of Sweden, several investigations have been carried out (Björck and Clemmensen, 2004; de Jong et al., 2006, 2007; Kylander et al., 2013, 2016, 2018; Sjöström et al., 2022).
Many sites with windblown material occur in SE Norway (Klemsdal, 1969). Among these, the aeolian dune fields at Starmoen and in the Jømna valley east of Elverum were investigated more thoroughly (Alexanderson and Henriksen, 2015; Klemsdal, 2010). However, the origin of the dune fields at Starmoen and in the Jømna valley as primarily aeolian has recently been challenged and reinterpreted by Hansen et al. (2020) as the result of catastrophic megafloods from the drainage of ice-dammed lake Nedre Glåmsjø (Longva and Bakkejord, 1990; Longva and Thoresen, 1991). Several aeolian sites exist along the coast of Norway (Klemsdal, 1969). Among these studies have been carried out in Finnmark (Samuelsen, 1935; Sjögren, 2009), Andøya (Nielsen et al., 2016a), Langøya in Vesterålen (Nielsen et al., 2016b), Jæren (Prøsch-Danielsen and Selsing, 2009; Selsing and Mejdahl, 1994) and at Lista (Høyland, 1974).
First described by Leonhard (1824), loess is terrestrial clastic sediment. Formed by the accumulation of windblown dust, loess is predominantly composed of silt-sized grain fractions (0.02–0.05 mm) with smaller amounts of incorporated fine sand and clay (Muhs, 2013; Pye, 1995; Smalley et al., 2011; Smalley and Marković, 2019; Újvári et al., 2016; Vandenberghe, 2013; Vandenberghe et al., 2018). For loess to form, three fundamental time-dependent requirements are necessary to fulfil: (1) a sustained source of dust, (2) available wind energy to transport the dust, and (3) a suitable accumulation site/trap (Albani et al., 2015; Muhs, 2013; Pye, 1995; Seppälä, 2004). Sources of loess may be local, neighbouring or distant, and grain-size distribution appears to shift towards finer grains with increasing distance from the source area (Muhs, 2013; Pye, 1995; Vandenberghe et al., 2018). Loess is usually homogenous, loosely cemented by calcium carbonate (CaCO3), and often yellowish brown or buff by colour. In addition to a characteristic CaCO3 content of 5–35%, loess typically consists of 50–70% quartz (SiO2) and 10–20% feldspars and micas (Albani et al., 2015; Editors of Encyclopaedia Britannica, 2010; Muhs, 2013; Pye, 1995). The composition of loess is more complicated than that of aeolian sand as the mineralogy reflects the geology of the source rocks (Seppälä, 2004; Stevens et al., 2007). Postdepositional processes also make it challenging to find a commonly accepted definition of loess (Smalley et al., 2011; Sprafke and Obreht, 2016). During the Quaternary, loess primarily formed in periglacial environments, but it also occurs in semi-arid parts of some lowland deserts and on the margins of high mountain ranges (Frechen, 2011; Pye, 1995).
Except for windblown dust in Iceland and some coastal areas, the occurrence of aeolian silt is not usual, and soils consisting of
Inland dunes are ‘

Overview map of southern and central Scandinavia showing study area and sites taken into account when discussing regional implications of the occurrence of early Holocene loess in Folldal.
Mention of loess in the scientific literature from Norway is near non-existent, but Björlykke (1916, 1918) discuss the possible occurrence of loess in silty soils at Romerike (‘
The shift from glacial to interglacial conditions occurred at the Younger Dryas-Preboreal transition in NW Europe. This shift dates to 11,530 +40/−60 cal. yr BP at Kråkenes on the outmost coast of western Norway (Figure 1) (Gulliksen et al., 1998). During the first 500 years after the Younger Dryas-Preboreal transition, an abrupt and distinct rise in July temperature of 8–10°C occurred at Kråkenes (Birks et al., 2000; Birks and Ammann, 2000). However, summer (July) temperatures close to present values first took place at about 10,000 cal. yr BP (Birks and Ammann, 2000; Nesje et al., 1991; Nesje and Kvamme, 1991.
In central Scandinavia, the final transition from glacial to interglacial conditions occurred during the late Preboreal/early Boreal (ca. 10,500–9500 cal. yr BP). During this period, the large ice-dammed lakes S-SE of the main watershed drained (Longva and Bakkejord, 1990; Longva and Thoresen, 1991; Stroeven et al., 2016; and references therein), and the pine-tree limit as a proxy of July/tetraterm (June-September) temperature obtained the Holocene maximum altitude from about 9500 to 9000 cal. yr BP (Aas and Faarlund, 1988; Dahl and Nesje, 1996; Kullman, 1980, 1995; Moe, 1979; Paus and Haugland, 2017; and references therein). However, the timing and causes of annual to seasonally prevailing wind patterns (Bernhardson et al., 2019; Bernhardson and Alexanderson, 2017, 2018; Sundborg, 1955) and hydroclimate (Digerfeldt, 1988; Digerfeldt et al., 2013) related to this significant shift in Scandinavia are not fully known or understood.
The main objectives of this study are to verify and date the possible occurrence and deposition of loess in Folldal, east-central southern Norway (Figure 1), and to discuss its origin and local and regional (synoptic scale) climatic implications. In addition, the results are discussed against relevant data reflecting hydroclimate and wind patterns in southern and central Sweden.
Study area
Folldal is a ca. 70 km long and generally broad valley from Hjerkinn on Dovrefjell to Alvdal in Østerdalen, Innlandet County, east-central southern Norway (Figure 1). The W-NW to E-SE orientated valley has a thick bedrock cover of superficial sediments dominated by till and large deposits of glaciofluvial, fluvial and glaciolacustrine sediments (Sollid and Carlson, 1980; Sørbel and Tolgensbakk, 2006; Thoresen and Follestad, 1999). The lack of over-deepened bedrock basins, common in Norway because of glacial erosion, is striking in Folldal. Water-filled kettle holes, however, exist in glaciofluvial terraces in Folldal/Grimsdalen (Sollid and Carlson, 1980) (Figure 2). Among these, Butjønna, Kroktjønna and Røtjønna at Stormoen are calcareous lakes based on the occurrence of charophytes (Figure 2) (Klepsland, 2007; Larsen et al., 2014; Miljødirektoratet Naturbase Faktaark, 2013). Next to Kroktjønna and Røtjønna at Stormoen, there are several peat-filled kettle holes. A core from a kettle hole named Stormoen Rus 15-01 is studied in detail (Figures 2 and 3a). There are no observations of colluvial activity or inlets in the catchment of the kettle hole.

Overview map of study area with key sites and place names. Geographical extent of calcareous relict valley sandur in upper Grimsdalen is indicated. Overflow gaps (OGs) with altitude for meltwater drainage to upper Grimsdalen during the deglaciation and related calcareous bedrock bands (modified from Bjerkgård et al., 2002; Nilsen and Wolff, 1989) influencing on the relict valley sandur are shown. To visualise the impact of calcareous material in the study area the reported observations of the calcicolous plant mountain aven (Dryas octopetala L.) are shown (Artsdatabanken, 2020). A-A′ represents the schematic transect in the conceptual model from source area to loess accumulation in Figure 8.

(a) Photo showing the relict marl pond (kettle hole) at Stormoen which is now filled up with peat. The Russian core Rus 15-01 was retrieved from central parts of the kettle hole. (b) Photo showing the shallow extant marl pond Lomslåtjønne seen from south. (c) Photo showing the shallow extant marl pond Myrlitjønne (our name) seen from south. (d) Photo showing details from the shallow bottom of the extant marl pond Myrlitjønne. The tracks seen on the photo are suggested to be after a moose (Alces alces L.). All photos: S.O. Dahl.
From the watershed towards Gudbrandsdalen at an altitude of 1160–1112 m to the west, Grimsdalen is a ca. 50 km long tributary valley which coalesces with Folldal at Stormoen (alt. ca. 660 m). The valley has a general W-SW to E-NE orientation (Figure 2). The only exception is a few km in the middle part with a NW to SE direction to pass Streitkampen (mountain) (alt. 1214 m). The upper ca. 14 km is narrower and more V-shaped, while the remaining valley is generally wide and without over-deepened bedrock basins. Exposed bedrock is typical in the uppermost V-shaped part and between the lower and middle Grimsdalen. The remaining valley has a high bedrock coverage of superficial sediments. Up to 50 m high terraces dominated by glaciofluvial and glaciolacustrine sediments surround the valley bottom along river Grimsa in the lower part. In the middle, significant accumulations of glaciofluvial sediments dominate along the southern valley side. Upper Grimsdalen has a ca. 11 km long stretch with a relatively flat valley bottom (altitudinal range from ca. 960–880 m) with terraces along Grimsa. These terraces consist of glaciofluvial, glaciolacustrine and fluvial sediment, and high-lying parts are suggested remnants after a relict glaciofluvial valley sandur. The occasionally active low-lying terraces along the river represent a modern flood plain (modified from Sollid and Carlson, 1980). Sediment input through meltwater overflow gaps from W-SW and N-NW during the last deglaciation likely contributed to the relict sandur in upper Grimsdalen. Just E-NE of Streitkampen, and between Folldal and Grimsdalen, two small and shallow ponds named Lomslåtjønne and Myrlitjønne (our name) have been investigated as part of this study (see Figures 2 and 3b–d). There are no observations of colluvial activity or inlets in the catchments of the two ponds.
The bedrock in Folldal and Grimsdalen is strongly affected by local allochtons (nappes) from the Caledonian orogen (500–400 Ma). Most of the area south of Folldal and the lower and middle parts of Grimsdalen consists of metasedimentary rocks from the middlemost allochthon with sandstone, quartzite, some conglomerate, and others in places. To the north of central Folldal and lower and middle Grimsdalen, and on both sides of upper Grimsdalen, the uppermost allochthon dominates. This allochthon consists of an ophiolite complex with metasedimentary and metavolcanic rocks, including phyllite, mica schist, marble, graphite schist, calcareous mica schist, calc-silicate schist, and others (Norges Geologiske Undersøkelse, 2021). Based on Nilsen and Wolff (1989) and Bjerkgård et al. (2002), Figure 2 shows the extent of calcareous phyllite and phyllite/calcareous schist in the catchment of Grimsdalen and central/upper Folldal. Note the calcareous rocks surrounding upper Grimsdalen, which are particularly interesting for this study.
Located within the ‘
Methods
The core Stormoen Rus 15-01 was retrieved using a Russian corer and was stored dark and chilly at 4°C. The core was first visually inspected, photographed and logged in the sediment laboratory. After this, geochemical analysis was performed using an Itrax XRF (X-ray fluorescence) core scanner (Croudace et al., 2006). First, the measurement was run using a Chrome (Cr) X-ray tube every 200 μm (
Organic- and inorganic content, as well as dry bulk density (DBD), was measured every 0.5 cm (1 cm3) throughout the core by the loss-on-ignition (LOI) method (Dean, 1974; Heiri et al., 2001). The samples were first weighed for wet bulk density and dried at 105°C for 12 h before being weighed for dry bulk density (g/cm3) and water content (%). Next, the samples were ignited at 550°C for 1 h, cooled in a desiccator, and weighed. LOI550 was calculated based on the percentage change between the ignition residue and dry weight. Afterwards, the inorganic carbonate (LOI950) was determined every 5 cm by heating the samples a second time to 950°C for 1 h and multiplying the ash residue by a factor of 1.36 (Bengtsson and Enell, 1986; Heiri et al., 2001).
Grain size analyses were examined continuously between 439 and 404 cm in Stormoen Rus 15-01, as well as four samples from Myrlitjønne and Lomslåtjønne, using a Malvern Mastersizer 3000 laser diffraction particle size analyser with a Hydro LV unit (Sperazza et al., 2004). Most samples (408–439 cm,
The mineralogy in Stormoen Rus 15-01 was identified and compared with deposits from Lomslåtjønne and Myrlitjønne by visual description, grain size analyses and content of LOI950 inorganic carbonate. Finally, the geochemical composition was analysed by use of qualitative x-ray diffraction (XRD) (Epp, 2016. Selected samples for XRD analyses were initially dried and ground to a fine powder using a mortar and pestle. A
Nine samples of organic material from Stormoen Rus 15-01 were submitted for 14C-dating. First, slices of sediment (0.5–2 cm thick) from the core were sieved through 125 μm using distilled water. Next, terrestrial plant remains were identified and isolated under the microscope before being dried overnight (50°C) and stored dark and chilly in sealed vials. Subsequently, samples were dated by accelerator mass spectrometry (AMS) radiocarbon dating at the Poznań Radiocarbon Laboratory in Poland before using the IntCal20-calibration curve to calibrate the dates to calendar years (Reimer et al., 2020) using Calib 8.2 (Stuiver et al., 2021). Finally, the age-depth modelling of the dates followed the routine in Clam 2.4 (Blaauw, 2010) using the open-source software R.
The age-depth model was calculated based on 6 of the radiocarbon-dated samples using smooth-spline interpolation (smoothing 0.2), calibration curve ‘IntCal.20C14’ for the Northern Hemisphere and 0.95 probability. The radiocarbon dates at 404, 416 and 420 cm depth are older than underlying ages and are therefore marked as outliers in the model. Outliers in aeolian erosional and depositional environments like this are expected because of shifting strength and direction of the wind.
Bacon (Blaauw and Christen, 2011) was evaluated for the age-depth modelling. There are, however, features like marked shifts in the accumulation rates in the investigated core Rus 15-01 that represent substantial challenges for this type of Bayesian analysis. You need to define the expected mean accumulation rate and depth resolution as priors in Bacon, and changing these values will effectively decide the dates which are treated as outliers. Because of this it is not possible to construct complete and reasonable age-depth models for Rus 15-01 in Bacon. Clam, however, uses a less advanced statistical approach, and it allows us to use the assumption of gradually changes in the sedimentation rate (through a smooth spline). Hence, and even if Clam may underestimate when there are large gaps between dates (Trachsel and Telford, 2017), this is the preferred choice for the age-depth modelling.
A field investigation of the relict valley sandur and modern flood plain in upper Grimsdalen is an element of this study. As part of this investigation, a compilation of superficial sediments and available data linked to the extent of calcareous rocks and related flora for central/upper Folldal and the entire catchment of Grimsdalen has been performed (Figure 2).
Results
Lithostratigraphy and chronology of Stormoen Rus 15-01
A core from the early deglaciation (Stormoen Rus 15-01) was retrieved by a Russian corer from a peat-covered kettle hole at Stormoen (Figures 2 and 3a). The core has been visually described and analysed for various sediment parameters and geochemical composition (Table 1). The chronology depends on 9 AMS radiocarbon-dated samples, eight using terrestrial plant macrofossils and one bulk sample (Table 2). The lithostratigraphy and key sediment parameters of Stormoen Rus 15-01 are shown in Figure 4, and grain size distributions of individual samples are displayed in Figure 5. An example of qualitative XRD analysis is presented in Figure 6, and age-depth modelling, including sedimentation rates, appears in Figure 7.
Description of the loess/marl localities Lomslåtjønne, Myrlitjønne and Stormoen Rus 15-01.
Distance from upstream start of calcareous relict sandur along valley bottom towards E-NE in upper Grimsdalen.
Measured at 70 cm depth next to margin – most likely more in central parts of ponds.
Measured on surface sediments close to margin. Most likely higher values deeper in the stratigraphy and in central parts of lake (Pentecost, 2009, and references therein).
Number of samples.
Suggested to be extant/relict clastic marl ponds according to definition by Pentecost (2009). See text for further discussion.
AMS radiocarbon dates from Stormoen Rus 15-01 from Stormoen, Folldal.
Calibrated by use of IntCal20/Calib 8.2 (Stuiver et al., 2021).
Calibation is rounded to nearest 10 year for samples with standard deviation greater than 50 years.

Lithostratigraphy of Stormoen Rus 15-01 with organic and inorganic content based on loss-onignition (LOI) 550°C and 950°C, respectively, dry bulk density (DBD), magnetic susceptibility (MS), XRF analysis represented by Ca/Inc+coh and Ti/Inc+coh adjusted for density changes and water content of the sediments (Kylander et al., 2011) and grain size analysis (the interval between 404 and 407 were digested in a 30% H2O2 aqueous solution prior to the analysis because of high organic content). Radiocarbon dates and units based on sediment parameters are shown. Note the yellowish-brownish colour of units C-F and partly unit G suggested to be representative of loess. The lithostratigraphy is further discussed in the text.

Continuous grain size analysis with 1 cm resolution between 439 and 404 cm depth in Stormoen Rus 15-01. All samples between 439 and 405 cm depth were run without pre-treatment, whereas the uppermost 4 samples (407–404 cm depth) were digested in a 30% H2O2 aqueous solution prior to the analysis because of high organic content. Unit A consists of very fine to fine sand and represents the early deglaciation, unit B consists of very coarse silt and mark the transition from early deglaciation to an environment dominated by loess deposition, units C, D, E and F consist of coarse silt (~22.4 μm) and are dominated by loess deposition, and unit G consists of coarse silt which is influenced by increased organic production.

Qualitative XRD analysis showing the identification of calcite (Ca) as part of CaCO3 from Stormoen Rus 15-01 and surface sediments from Myrlitjønne.

Age-depth model based on 6 AMS radiocarbon dates (marked with blue) using terrestrial plant macrofossils from Stormoen Rus 15-01. Two inverted dates (marked with red) based on terrestrial plant macrofossils are not used in the age-depth modelling. One date on bulk sediment is not shown (see Table 2 for overview of radiocarbon dates). The samples are dated at Poznań Radiocarbon Laboratory. The IntCal20 calibration curve was used for calibrating the samples to calendar years (Reimer et al., 2020) using Calib 8.2 (Stuiver et al., 2021). The dates were then applied for age-depth modelling based on Clam 2.4 (Blaauw, 2010) using the open-source software R. Sedimentation rates for the different units as defined in Figure 4 show a slow sedimentation rate (0.035 cm/yr) during the final deglaciation (unit A), and an increase (0.062 cm/yr) during the transition (unit B) towards an environment dominated by loess deposition. A high sedimentation rate of loess (0.161 cm/yr) is recorded in unit C. This is followed by a marked decrease (0.030 cm/yr) in unit D, before the sedimentation rate of loess again increases (0.055 cm/yr) in unit E. The very high sedimentation rate (1.00 cm/yr) in unit F is attributed to a marked increase in organic production. See text for further discussion.
Unit A, 444.5–437.0 cm depth. As reflected by LOI550, the organic content is very low (1–3%), DBD is between 1.2 and 1.0 g/cm3 and MS (10−6) is from 1.1 to 0.4. Calcium (Ca) adjusted for density changes and water content (Ca/inc+coh) is very low, as is the inorganic content based on LOI950. Titanium (Ti) adjusted for density changes and water content (Ti/inc+coh) has mean values of about 0.3 but is weakly declining towards the top of the unit. The colour of unit A is dark brown to more blackish towards the top, and grain sizes are from fine to very fine sand. Based on the age-depth modelling, unit A lasts from ca. 10,500 to 10,390 cal yr BP (~110 years). The mean sedimentation rate of the unit is 0.036 cm/yr (Figures 4, 5 and 7).
Unit B, 437.0–435.5 cm depth. LOI550 values are low (2–4%), DBD is about 1.1 g/cm3, and MS (10-6) drops from about 0.4 to near zero. Calcium (Ca/inc+coh) increases from near zero to about 12. The calcium increase is reflected in the inorganic content, with values of LOI950 rising from close to zero to ca. 20%. Titanium (Ti/inc+coh) drops from about 0.3 to 0.1, and the mean grain size drops from very fine sand to very coarse/coarse silt at the top. The colour of unit B is black. Age-depth modelling indicates a duration of the unit from ca. 10,390 to 10,365 cal yr BP (~ 25 years), and it has a mean sedimentation rate of 0.055 cm/yr (Figures 4, 5 and 7).
Unit C, 435.5–422.0 cm depth. Values of LOI550 are rising upwards from about 3 to ca. 4–8%, DBD drops from ca. 0.8 to 0.4 g/cm3, and MS (10−6) is close to zero. Calcium (Ca/inc+coh) increases from 12 to stabilising between 14 and 15, but with a maximum at the upper end of the unit. The same general pattern for this unit is observed for the inorganic content as reflected by values of LOI950, whereas titanium (Ti/inc+coh) drops to values close to zero. Grain size distribution remains stable within coarse silt. The colour of unit C is yellowish brown with some sporadic speckles/thin layers, which are more brownish. Based on age-depth modelling, unit C lasts from about 10,365–10,190 cal yr BP (~ 175 years), with a mean sedimentation rate of 0.158 cm/yr (Figures 4, 5 and 7).
Unit D, 422.0–417.0 cm depth. LOI550 is between ca. 6 and 15%, DBD drops from ca. 0.6 to 0.4 g/cm3 and MS (10−6) is close to zero. Calcium (Ca/inc+coh) drops from a value of ca. 15 to a level between 10 and 11, as also recorded in LOI950. Titanium (Ti/inc+coh) is close to zero, whereas mean grain size is stable within coarse silt. The colour of unit D is yellowish brown, but with some thin brownish layers towards the top. Age-depth modelling indicates that unit D lasts from ca. 10,190–10,020 cal yr BP (~170 years) with a mean sedimentation rate of 0.030 cm/yr (Figures 4, 5 and 7).
Unit E, 417.0–414.5 cm depth. Values of LOI550 vary between ca. 8.6 and 22% but is lowest in the upper part. DBD varies from 0.3 to 0.5 g/cm3 but is highest in the upper end, and MS (10−6) is close to zero. Calcium (Ca/inc+coh) increases from ca. 11 to 14, and a similar pattern is observed in values of LOI950, whereas Titanium (Ti/inc+coh) is close to zero. Mean grain size has only minor fluctuations within coarse silt. Unit E is yellowish brown, but brownish layers become more abundant towards the top. The age-depth modelling indicates a duration of unit E from ca. 10,020 to 9950 cal yr BP (~70 years) with a mean sedimentation rate of 0.036 cm/yr (Figures 4, 5 and 7).
Unit F, 414.5–409.0 cm depth. LOI550 rises from ca. 7.5% to 13.6%, DBD drops from ca. 0.5 to 0.45 g/cm3, and MS (10−6) is close to zero. Calcium (Ca/inc+coh) gradually drops from ca. 14 to 7. This decrease is also observed in inorganic carbon based on LOI950, whereas Titanium (Ti/inc+coh) remains close to zero. Mean grain size remains stable within coarse silt, and the unit’s colour transitions gradually from yellowish brown to more brownish at the top. Based on age-depth modelling, the unit lasts from ca. 9955 to 9850 cal yr BP (~105 years), with a mean sedimentation rate of 0.055 cm/yr (Figures 4, 5 and 7).
Unit G, 409.0–402.0 cm depth. LOI550 values rise quickly from 13.6% to 85.7%, DBD values drop from ca. 0.45 to 0.1 g/cm3, and MS (10−6) is close to zero. Calcium (Ca/inc+coh), inorganic carbon based on LOI950 and Titanium (Ti/inc+coh) are all close to zero. Grain size analyses are not available (not possible to perform). The colour of unit G is black. Based on age-depth modelling, the unit lasts from ca. 9850 to 9780 cal yr BP (~70 years) with a mean sedimentation rate of 0.095 cm/yr (Figures 4, 5 and 7).
Unit H, from 402.0 cm depth and upwards. LOI550 values are stable above 85%, DBD is close to 0.1 g/cm3, and MS (10−6) is close to zero. Calcium (Ca/inc+coh), inorganic carbon based on LOI950 and Titanium (Ti/inc+coh) are all close to zero. Grain size analyses are not available (not possible to perform). The colour of unit G is black. Age-depth modelling is unavailable, but the unit starts at ca. 9780 cal. yr BP (Figure 4).
Qualitative XRD analysis, inorganic carbonate and grain size
Samples from units C and E in Stormoen Rus 15-01, sub-surface sediments from Lomslåtjønne and surface sediments from Myrlitjønne were analysed using qualitative XRD, and the amount (%) of inorganic carbonate (CaCO3) in the same samples were obtained using LOI950. Table 1 shows the results of the analyses and thickness estimates of calcareous sediments. Calcite (CaCO3) is identified at all three sites using XRD (see Figure 6). Based on LOI950, units C and E in Stormoen Rus 15-01 have maximum values of inorganic carbonate of 50.72% and 51.68%, respectively. Along the margin of Lomslåtjønne and beneath 55 cm of peat, analysis shows a maximum content of inorganic carbonate at 70 cm depth of 44.94%. From surface sediments near the margin of Myrlitjønne, the amount of inorganic carbonate based on LOI950 is 52.17–51.99%.
The thickness estimate from Stormoen Rus 15-01 is reliable. In contrast, figures from Lomslåtjønne and Myrlitjønne are only minimum estimates based on observations from the surrounding bog and the margin of the small ponds, respectively.
The investigated kettle hole at Stormoen (core Rus 15-01) has a peat cover (Figure 3a). It differs markedly from the small, open, shallow ponds Lomslåtjønne (Figure 3b) and Myrlitjønne (Figure 3c and d), which, seen from above, has a brownish-yellowish and whitish-yellowish colour, respectively. The calcareous sediments at all three sites are well sorted, contain significant amounts of quartz, and have mean grain size within coarse silt (~22.4 µm) (Figures 5 and 6). Based on Nilsen and Wolff (1989), the bedrock at Lomslåtjønne, Myrlitjønne and Stormoen Rus 15-01 consists of greenschist and amphibolite, and according to Bjerkgård et al. (2002), the bedrock at Lomslåtjønne and Myrlitjønne is tuffitic greenschist. The bedrock at the Stormoen site is hard to map, as it has a thick cover of glaciofluvial deposits dominated by feldspar and quartz originating south of Folldal.
No survey for charophytes to evaluate if the ponds are calcareous exists for Lomslåtjønne and Myrlitjønne. However, charophytes like delicate stonewort (
Relict valley sandur and modern floodplain in upper Grimsdalen
Calcareous ponds/peat-filled kettle holes in non-calcareous superficial sediments/rocks need a source. The geographical distribution led to an investigation of a relict valley sandur/modern floodplain in upper Grimsdalen towards W-SW (Figure 2). The valley sandur was built out in an ice-dammed lake with erosion basis at an altitude of ca. 880 m. It covered the central and upper parts of Grimsdalen, and based on optically stimulated luminescence (OSL) dates adjusted for mean water content and depth from Svartdalsbekken (Figure 2), the formation of the valley sandur took place between ca. 17.2 and 13.8 ka (modified from Bøe et al., 2007). Quartz/feldspar-rich sediments deposited from S-SE dominate the downstream central part of the ice-dammed lake. In contrast, calcareous material of both local and more distant origin strongly influences the relict valley sandur/modern floodplain in upper Grimsdalen. Distant input of calcareous material originates from meltwater streams through overflow gaps from W-SW and N-NW during the deglaciation (Figure 2) Mountain aven (
Discussion
Occurrence, origin and timing of loess deposition in Folldal
Yellowish brown clastic sediments in Folldal, consisting of homogenous coarse silt (~22.4 µm) and some clay and fine sand, are very well sorted and contain high amounts of carbonate and quartz. (Figures 4–6). As the geographical distribution can only be explained by aeolian processes (Figures 2 and 8), these sediments are suggested to be loess (

Conceptual model for loess accumulation from the source area in upper Grimsdalen to the accumulation area (sink) in Folldal.
The present distribution of loess in Folldal is very scattered. Primarily the windblown dust occurs in sediment traps like peat-covered kettle holes with a bottom below the local groundwater table or in small and shallow ponds. The present distribution and the need for a calcareous source area strongly indicate a W-SW origin for the loess. Hence, the provenance of the loess is suggested to be a relict valley sandur/existing flood plain in upper Grimsdalen just north of the Rondane mountain massif (Figures 2 and 8). About 25–30 km downwind from the suggested source area, an about 35 cm long core section of loess from a kettle hole at Stormoen in Folldal (Stormoen Rus 15-01) is investigated. Based on age-depth modelling using 9 AMS radiocarbon dates (Table 2, Figure 7), loess deposition took place between ca. 10,390 and 9780 cal. yr BP (ca. 610 years). The highest loess accumulation appears to have occurred in two main periods from ca. 10,390 to 10,190 (200 years) and from ca. 10,020 to 9950 (70 years) cal. yr BP. A reduced deposition rate separates the two periods from 10,190 to 10,020 cal. yr BP (170 years) (Figures 7 and 9).

Deposition of loess during the early Holocene in Folldal. Note the somewhat reduced accumulation of loess in between the two deposition maxima of loess during Folldal Loess Event (FLE) 1 and 2.
Considering a suitable source area and available winds, the response time to erode, transport and accumulate loess in Folldal is suggested to be short, season to 1 year.
The mean grain size, however, is in the finer range of average loess (Muhs, 2013; Pye, 1995; Vandenberghe et al., 2018). The ca. 28.9–20.1 km long transportation and approximately 250–300 m airborne uphill move across Streitkampen from the suggested source area to the accumulation site (trap) at Stormoen may explain why the windblown dust is in the finer range of average loess (Table 1, Figure 8).
Marl lakes/ponds
Marl (mergel in German) lakes/ponds can be defined based on the composition of the lacustrine sediments and the amount of dissolved carbonate in the water. Morgan and Britton (1977) defined marl lakes as those lakes whose waters, on average, contained ⩾100 mg/L dissolved calcium carbonate (CaCO3), whereas Kelts and Hsu (1978) defined them as lakes whose lacustrine sediments contained >60 % CaCO3. Pentecost (2009), however, argued for a change in the definition of marl lakes as lakes with calcium carbonate as the main component (>50 % CaCO3). Despite Schnurrenberger et al. (2003) and others advocating to replace the term
The amount of CaCO3 in Lomslåtjønne based on LOI950 analysis on a sediment sample beneath (70 cm depth) the present peaty margin of the pond is 44.94% (Table 1). In contrast, the maximum amount of CaCO3 based on LOI950 recorded in surface samples of lacustrine sediments near the margin of Myrlitjønne and beneath the peat at Stormoen Rus 15-01, are 52.17% and 51.68%, respectively. The figure from Lomslåtjønne, however, is suggested to underestimate the actual values with depth and in central parts of the pond because of more interbedding of organic material along the margins (Pentecost, 2009, and references therein). Hence, according to Pentecost (2009), Lomslåtjønne and Myrlitjønne are suggested to be extant marl ponds, and the kettle hole at Stormoen Rus 15-01 is a relict marl pond.
Based on the occurrence of charophytes like delicate stonewort (
The composition of the local bedrock at the investigated sites cannot explain the high content of CaCO3 in the lacustrine sediments. Hence, these sediments must have an allochthonous origin and are likely to be clastic marl due to the input of suspended carbonate (Pentecost, 2005: 302). The accumulation of allochthonous carbonate may result from aeolian, colluvial or (glacio-) fluvial activity (Pentecost, 2009; Schnurrenberger et al., 2003). However, only input of suspended carbonate as loess by aeolian processes can explain the geographical distribution of the extant small and shallow marl ponds Lomslåtjønne and Myrlitjønne, the relict peat-filled marl pond Stormoen Rus 15-01 and the calcareous ponds Kroktjønna and Røtjønna (Figures 2 and 8).
Climate during loess deposition in Folldal
Strong wind is most likely to occur in Folldal during autumn and winter. As reflected in the present spatial distribution of leeward accumulation of snow, the predominant snow-bearing wind direction in Rondane is from SW (Bjørbæk, 1993). Exposed towards W-SW, the broad watershed (wind gap) area between Grimsdalen and Gudbrandsdalen may have played an essential role in directing and strengthening wintertime W-SW winds towards Folldal. Hence, tracking the loess in Folldal to its source area indicates a provenance and erosion/transport by strong W-SW winds from the calcareous relict sandur/existing flood plain in upper Grimsdalen. The occurrence of plant remains (e.g. leaves) from mountain aven (
The present dry wintertime climate with slight snow was likely even more pronounced during loess deposition in Folldal. Gradually becoming more distinct upwards, however, thin organic laminae recorded in the loess sequence at Stormoen in Folldal may indicate relatively warm but likely very dry (semi-arid) summers with access to water as the limiting factor for organic production. Loess accumulation ended ca. 9780 cal. yr BP. After that, a rapid increase (a few decades?) in organic production which quickly reached stable LOI values of about 85% (Figure 4), took place. The probable result of a shift to a warmer and wetter climate, this transition is suggested to represent the initiation of the local Holocene climatic optimum. An AMS radiocarbon date on needles from scots pine (
Predominantly W-NW winds linked to the deposition of aeolian sand and related landforms along a S-N transect in south-central Sweden have recently been reviewed by Bernhardson et al. (2019). Primarily based on OSL dating, most dune fields stabilised ca. 10.5–9.0 ka despite the local deglaciation along the transect taking place over almost 4000 years (Stroeven et al., 2016). Hence, regional environmental conditions with an unstable and dry early Holocene climate likely explain the prolonged dune activity followed by a nearly synchronous stabilisation. The lowest Holocene lake-level stands recorded between 10,500 and 9500 cal. yr BP in lake Bysjön (Digerfeldt, 1988) and lake Igelsjön (Digerfeldt et al., 2013) in south-central Sweden (Figure 1) are supporting dry conditions in Folldal. Probably due to a low precipitation/evaporation (P/E) ratio, the low groundwater level caused delayed plant colonisation of the dune fields before a vegetation cover made them more resilient to climate change.
OSL dates reflect the timing for the last aeolian sand transport prior to deposition/stabilisation. However, Bernhardson et al. (2019) found a peak in OSL ages around 10.0 ± 0.5 ka at sites along the S-N transect in south-central Sweden. At Bonäsheden, the largest continuous dune field in Sweden (Figure 1), Bernhardson and Alexanderson (2017) identified an apparent shift in the prevailing wind pattern from NW to W c. 10 ka. Local topography appears to have had little influence on the dune-forming regional winds. However, NW katabatic winds from the shrinking Scandinavian Ice Sheet likely impacted dune formation soon after the local deglaciation. In a later investigation of small dune fields in central Sweden, Bernhardson and Alexanderson (2018) could neither confirm nor reject the hypothesis. However, the loess deposition in Folldal by W-SW winds centred around 10,085 ± 305 cal yr BP (Figure 9) supports the hypothesis by Bernhardson and Alexanderson (2017). Strong W-SW wintertime winds are in southern Norway strongly linked to a positive mode of the North Atlantic Oscillation (NAO) (Nesje et al., 2000, and references therein).
Impact on permafrost and last remnants of the Scandinavian Ice sheet
Central parts of the Scandinavian Ice Sheet, including the Folldal area, are suggested to have been cold based during the Last Glacial Maximum (LGM), Younger Dryas and the following deglaciation (see Figure 4 in Stroeven et al., 2016, and references therein). Cold-based ice implies that the underlying soft sediments had sub-zero temperatures until the last ice remnants melted away. The present climate in Folldal is continental, with frigid and dry winters (Aune, 1993; Førland, 1993). However, the local loess deposition indicates a periglacial environment with even drier and windier conditions than at present. As a result, most permafrost in the predominantly soft-sediment-covered Folldal area (Sollid and Carlson, 1980) is likely to have survived until Late Preboreal/Early Boreal.
Summer (July) temperatures close to present values first took place in western Norway at about 10,000 cal. yr BP (Birks and Ammann, 2000; Nesje et al., 1991; Nesje and Kvamme, 1991), and coincides with loess deposition in Folldal centred about 10,085 ± 305 cal. yr BP (Figure 8). During and immediately after loess deposition in Folldal, a transition with widespread permafrost melting probably occurred before establishing an early Holocene steady-state landscape. Melting of the permafrost implies that most landforms with incorporated ice or frost-sensitive sediments experienced slope instability (Etzelmüller et al., 2020; Miesen et al., 2021) and most likely first stabilised about 10,000 cal yr BP. This steady-state landscape is, however, still influenced by stochastic episodes of intense snow melt or heavy spring/summer precipitation, which may cause major local and regional river floods and related colluvial processes (Bøe et al., 2006; Nesje et al., 2001).
The dry regional wintertime climate, as indicated by the loess deposition in Folldal between 10,390 and 9780 cal. yr BP probably had a marked impact on the winter balance (accumulation season) of the remaining Scandinavian Ice Sheet. Hence, the final drainage of the sizeable ice-dammed lake Nedre Glåmsjø in southeastern Norway, about 10,200 ± 200 cal yr BP (Høgaas and Longva, 2016), may have been influenced by this period of precipitation starvation.
Conclusions
Based on the presented data and discussion, the following main conclusions of local and regional importance are suggested:
In Folldal, east central southern Norway (Figure 1), windblown yellowish-brown dust fulfilling physical and geochemical criteria of loess, including a high carbonate content, is identified (Figure 4). The presence of
The two extant marl ponds Lomslåtjønne and Myrlitjønne, and the relict marl pond Stormoen Rus 15-01 are investigated. The marl ponds are located on low-carbonate bedrock, and an allochthonous origin for the marl is suggested. The most likely mechanism explaining the marl ponds is the input of suspended clastic carbonate as loess by aeolian processes (Figures 2 and 8, Table 1).
The scattered geographical distribution of loess in Folldal and the need for a calcareous source area suggest an upper Grimsdalen provenance from a relict valley sandur/existing flood plain (Figures 2 and 8).
The mean grain size (~ 22.4 µm) (Figures 5 and 6) of the aeolian silt is in the finer range of average loess. The ca. 25–30 km long transport distance and approximately 250–300 m airborne uphill move across Streitkampen from the suggested source area to the sediment trap at Stormoen Rus 15-01 may explain the finer grain size (Figure 8).
Having a suitable source area and available winds, the response time to erode, transport and accumulate loess in Folldal is likely to be seasonal to 1 year.
Based on age-depth modelling using 9 AMS radiocarbon dates from Stormoen Rus 15-01 (Table 2, Figures 4 and 7), loess deposition took place between ca. 10,390 and 9780 cal. yr BP (ca. 610 years). The highest loess accumulation appears to have occurred in two main periods from ca. 10,390 to 10,190 (200 years) and from ca. 10,020 to 9950 (70 years) cal. yr BP. A reduced deposition rate of loess separates the two periods from 10,190 to 10,020 cal. yr BP (170 years) (Figures 7 and 9).
Loess deposition in Folldal took place during late autumn/wintertime. Together with minimal summertime organic production in kettle holes beneath the groundwater table, this is suggested to reflect very dry to semi-arid climatic conditions with access to water as the limiting factor. The dry climate is in accordance with the lowest Holocene lake-level stands recorded between 10,500 and 9500 cal. yr BP in south-central Sweden (Figure 1) (Digerfeldt, 1988; Digerfeldt et al., 2013).
Loess deposition in Folldal by W-SW winds about 10,085 ± 305 cal. yr BP agrees with an apparent shift from prevailing NW to W winds at Bonäsheden in central Sweden ca. 10 ka (Figure 1) (Bernhardson and Alexanderson, 2017).
The loess deposition in Folldal is related to prevailing W-SW winds during late autumn and wintertime. This wind direction is closely associated with a positive mode of the North Atlantic Oscillation (NAO) in southern Norway.
The Scandinavian Ice Sheet was probably cold-based in Folldal during the last deglaciation (Stroeven et al., 2016), implying subglacial frozen ground. During and just after loess deposition in Folldal, widespread permafrost melting probably occurred, followed by an immediate increase in organic production as early initiation of a warmer and wetter Holocene climatic optimum (Figure 4). The ice remnants damming the sizeable ice-dammed lake Nedre Glåmsjø may have been affected by precipitation starvation during loess deposition prior to the final drainage 10,200 ± 200 cal yr BP.
