Introduction
The Kola Peninsula, next to the White and Barents Seas, was subjected to several marine transgressions and glacier expansions during successive late Pleistocene warm and cold stages. In the early Marine Isotope Stage (MIS) 5, huge areas in NW Russia, including valleys and plains of the Kola Peninsula, were flooded by transgressive seas associated with a number of global hydrological, geologic, and climatic processes. The Boreal marine transgression that is identified by many researchers (Lavrova, Reference Lavrova1960; Gudina and Yevzerov, Reference Gudina and Yevzerov1973; Faustova, Reference Faustova and Velichko1984; Molodkov and Raukas, Reference Molodkov and Raukas1998; Ikonen and Ekman, Reference Ikonen and Ekman2001; Molodkov and Bolikhovskaya, Reference Molodkov and Bolikhovskaya2002, Reference Molodkov and Bolikhovskaya2009; Lunkka et al., Reference Lunkka, Kaparulina, Putkinen and Saarnisto2018; Korsakova, Reference Korsakova2021) was the largest in the late Pleistocene history of the Kola Peninsula. It is believed that the Kola region and adjacent White Sea depression were covered by the Scandinavian Ice Sheet (SIS) during the last cold stage (late Valdaian, late Weichselian, MIS 2) (Lavrova, Reference Lavrova1960; Ekman and Iljin, Reference Ekman and Iljin1991; Kleman et al., Reference Kleman, Hättestrand, Borgström and Stroeven1997; Yevzerov, Reference Yevzerov, Matishov, Denisov, Dzhenyuk and Tarasov2001; Svendsen et al., Reference Svendsen, Alexanderson, Astakhov, Demidov, Dowdeswell, Funder and Gataullin2004; Demidov et al., Reference Demidov, Houmark-Nielsen, Kjær and Larsen2006; Hättestrand et al., Reference Hättestrand, Kolka and Stroeven2007; Hughes at al., Reference Hughes, Gyllencreutz, Lohne, Mangerud and Svendsen2016; Stroeven et al., Reference Stroeven, Hättestrand, Kleman, Heyman, Fabel, Fredin and Goodfellow2016; Lunkka et al., Reference Lunkka, Kaparulina, Putkinen and Saarnisto2018; Boyes et al., Reference Boyes, Linch, Pearce, Kolka and Nash2021). These two late Pleistocene events are clearly traced by palaeoenvironmentally controlled sedimentary processes and are relatively well represented in the depositional successions. Details of the Kola region's sedimentary history during the late MIS 5, MIS 4, and MIS 3 periods are more elusive, although the general outlines of the late Pleistocene stratigraphy are established. Known litho-chrono-biostratigraphic data, as well as sedimentary analyses and numerical dating studies and their interpretation, are summarized in Gudina and Yevzerov (Reference Gudina and Yevzerov1973), Ikonen and Ekman (Reference Ikonen and Ekman2001), Molodkov and Bolikhovskaya (Reference Molodkov and Bolikhovskaya2002, Reference Molodkov and Bolikhovskaya2009), Svendsen et al. (Reference Svendsen, Alexanderson, Astakhov, Demidov, Dowdeswell, Funder and Gataullin2004), Korsakova (Reference Korsakova2019, Reference Korsakova2021), Hughes et al. (Reference Hughes, Gyllencreutz, Lohne, Mangerud and Svendsen2016), Korsakova et al. (Reference Korsakova, Kolka and Semenova2016), Lunkka et al. (Reference Lunkka, Kaparulina, Putkinen and Saarnisto2018), Zaretskaya et al. (Reference Zaretskaya, Korsakova and Panin2019, Reference Zaretskaya, Korsakova, Molodkov, Ruchkin, Baranov, Rybalko, Lugovoy and Merkuliev2022), and Boyes et al. (Reference Boyes, Linch, Pearce, Kolka and Nash2021, Reference Boyes, Linch, Pearce and Nash2023).
However, different opinions still remain concerning the extent and age of the late Pleistocene glacial advances, as well as the number and timing of marine transgressions. Mollusc-bearing clayey mud, silt, and sand containing marine diatoms are known in the fluvial and marine terraces from several localities in coastal Kola Peninsula (Grave et al., Reference Grave, Gunova, Devyatova, Lavrova, Lebedeva, Samsonova, Cheremisinova, Grave and Koshechkin1969; Gudina and Yevzerov, Reference Gudina and Yevzerov1973; Ikonen and Ekman, Reference Ikonen and Ekman2001; Molodkov and Yevzerov, Reference Molodkov and Yevzerov2004; Korsakova et al., Reference Korsakova, Molodkov and Kolka2004, Reference Korsakova, Kolka and Semenova2016; Korsakova, Reference Korsakova2021). Early work by Gudina and Yevzerov (Reference Gudina and Yevzerov1973) attributed these marine beds to two middle Valdaian (middle Weichselian, MIS 3) transgressions—the relatively warm-water Ponoy transgression and the colder-water Strelna one, represented by so-called Ponoy Beds and Strelna Beds, respectively. Later, the Ponoy transgression was correlated with MIS 5e, and the Strelna transgression with the MIS 3 interglacial period (Yevzerov, Reference Yevzerov1993). On the other hand, based on biostratigraphic data (Lavrova, Reference Lavrova1960; Nikonov and Vostrukhina, Reference Nikonov and Vostrukhina1964; Armand and Lebedeva, Reference Armand, Lebedeva, Grave and Armand1966; Grave et al., Reference Grave, Gunova, Devyatova, Lavrova, Lebedeva, Samsonova, Cheremisinova, Grave and Koshechkin1969; Gudina and Yevzerov, Reference Gudina and Yevzerov1973; Strelkov et al., Reference Strelkov, Yevzerov, Koshechkin, Rubinraut, Afanasiev, Lebedeva and Kagan1976; Ikonen and Ekman, Reference Ikonen and Ekman2001), uranium-thorium (230Th/U), electron spin resonance (ESR), feldspar-based infrared optically stimulated luminescence (IR-OSL) dates, and amino acid ratios (Arslanov et al., Reference Arslanov, Yevzerov, Tertychnyy, Gerasimov, Lokshin, Velichko and Faustova1981; Miller and Mangerud, Reference Miller and Mangerud1985; Ikonen and Ekman, Reference Ikonen and Ekman2001; Molodkov and Yevzerov, Reference Molodkov and Yevzerov2004; Korsakova et al., Reference Korsakova, Molodkov and Kolka2004, Reference Korsakova, Kolka and Semenova2016; Zaretskaya et al., Reference Zaretskaya, Korsakova and Panin2019, Reference Zaretskaya, Korsakova, Molodkov, Ruchkin, Baranov, Rybalko, Lugovoy and Merkuliev2022), it has been concluded that the Ponoy Beds and Strelna Beds were deposited on the Kola Peninsula by two phases of a single marine transgression during the last interglacial, that is, the Mikulinian warming (MIS 5); a third marine unit formed during the MIS 3 warming.
It is considered that at least two late Pleistocene glacial events are represented on the Kola Peninsula, suggesting two separate ice advances (Nikonov, Reference Nikonov and Grichuk1966; Armand et al., Reference Armand, Armand, Grave, Yevzerov, Lebedeva, Grave and Koshechkin1969; Grave et al., Reference Grave, Gunova, Devyatova, Lavrova, Lebedeva, Samsonova, Cheremisinova, Grave and Koshechkin1969; Apukhtin, Reference Apukhtin1978; Hättestrand et al., Reference Hättestrand, Kolka and Stroeven2007; Yevzerov, Reference Yevzerov2016; Boyes et al., Reference Boyes, Linch, Pearce, Kolka and Nash2021, Reference Boyes, Linch, Pearce and Nash2023; Korsakova, Reference Korsakova2021) during MIS 5d–MIS 4 and MIS 2. Situated in a coalescing zone of the several ice sheets, Kola Peninsula was generally accepted to be covered by the SIS and Barents Sea Ice Sheet (BSIS) (Grosswald, Reference Grosswald2001; Grosswald and Hughes, Reference Grosswald and Hughes2002; Svendsen et al., Reference Svendsen, Alexanderson, Astakhov, Demidov, Dowdeswell, Funder and Gataullin2004; Kjær et al., Reference Kjær, Larsen, Funder, Demidov, Jensen, Hеkansson and Murray2006; Lambeck et al., Reference Lambeck, Purcell, Funder, Kjær, Larsen and Möller2006, Reference Lambeck, Purcell, Zhao and Svenson2010; Larsen et al., Reference Larsen, Kjær, Demidov, Funder, Grosfjeld, Houmark-Nielsen, Linge and Lyså2006; Stroeven et al, Reference Stroeven, Hättestrand, Kleman, Heyman, Fabel, Fredin and Goodfellow2016; Korsakova, Reference Korsakova2021). However, almost all researchers involved in the glacial geology of European Russia agree that neither the late Valdaian (late Weichselian, MIS 2) BSIS nor the Barents-Kara Ice Sheet advanced into the Kola Peninsula (Ekman and Iljin, Reference Ekman and Iljin1991; Yevzerov, Reference Yevzerov, Matishov, Denisov, Dzhenyuk and Tarasov2001; Svendsen et al., Reference Svendsen, Alexanderson, Astakhov, Demidov, Dowdeswell, Funder and Gataullin2004; Kjær et al., Reference Kjær, Larsen, Funder, Demidov, Jensen, Hеkansson and Murray2006; Larsen et al., Reference Larsen, Kjær, Demidov, Funder, Grosfjeld, Houmark-Nielsen, Linge and Lyså2006; Pasanen et al., Reference Pasanen, Lunkka and Putkinen2010; Stroeven et al., Reference Stroeven, Hättestrand, Kleman, Heyman, Fabel, Fredin and Goodfellow2016; Korsakova et al., Reference Korsakova, Vashkov, Nosova, Palacio, Hughes, Garcia-Ruiz and Andres2021), although the BSIS may have reached the Kola Peninsula in MIS 4 (Svendsen et al., Reference Svendsen, Alexanderson, Astakhov, Demidov, Dowdeswell, Funder and Gataullin2004; Lambeck et al., Reference Lambeck, Purcell, Funder, Kjær, Larsen and Möller2006; Larsen et al., Reference Larsen, Kjær, Demidov, Funder, Grosfjeld, Houmark-Nielsen, Linge and Lyså2006; Korsakova et al., Reference Korsakova, Kolka, Zozulya, Lavrushin, Khoreva and Chistyakova2007; Korsakova, Reference Korsakova2021).
Later glacial advances and related glacioaqueous erosion considerably destroyed unconsolidated deposits formed before these events. However, they have been preserved quite completely on the southern Kola Peninsula, where they outcrop along riverbanks or in coastal cliffs. Here, lithologic and chronological investigations of sedimentary successions were undertaken to obtain data about the ages, textures, and structure of deposits that would allow the history of the main changes of sedimentary settings to be determined. The primary objectives of this paper are to discuss and evaluate the obtained and available published results, identify the constituent units of the late Pleistocene depositional sequences in the southern Kola Peninsula, and associate them with the main marine and glacial events, as well as to determine the chronology of the main depositional stages. In this contribution, we present new lithostratigraphic observations and chronological data from the series of the key sections in the southern Kola Peninsula that can shed light on the late Pleistocene depositional settings and palaeogeographic events in the European north.
Study Area
The investigated sites are located on the western and southeastern Tersky Coast of the White Sea (Fig. 1). In the southern Kola Peninsula, bedrock is composed of Archaean and Lower and Upper Proterozoic metamorphic gneiss, amphibolite, migmatite, intrusive igneous granodiorite, tonalite, trondhjemite, and metasedimentary arcosic and polimict sandstone and siltstone (Mitrofanov, Reference Mitrofanov2001). The position of the river valleys tends to be predetermined by the neotectonic pattern of the northeastern Fennoscandian Shield, including the Kola Peninsula (Koshechkin, Reference Koshechkin1979; Strelkov et al., Reference Strelkov, Yevzerov, Koshechkin, Rubinraut, Afanasiev, Lebedeva and Kagan1976).
White Sea and location of the key sections: 1, Varzuga; 2, Chavanga; 3, Kamenka; 4, Chapoma 1; 5, Chapoma 2; 6, Bolshaya Kumzhevaya. Dotted coloured lines show glaciation limits during MIS 4 at 60 ka (magenta line; after Svendsen et al., Reference Svendsen, Alexanderson, Astakhov, Demidov, Dowdeswell, Funder and Gataullin2004) and at 70–65 ka (orange line; after Larsen et al., Reference Larsen, Kjær, Demidov, Funder, Grosfjeld, Houmark-Nielsen, Linge and Lyså2006) and during the last glacial maximum (LGM) (black line; after Demidov et al., Reference Demidov, Houmark-Nielsen, Kjær and Larsen2006). The insert shows the study area (blue rectangle) on the map of Europe. Digital terrain model by Farr et al. (Reference Farr, Rosen, Caro, Crippen, Duren, Hensley and Kobrick2007).
Quaternary deposits in the study area are mainly composed of glacial, glaciofluvial, glaciolacustrine, glaciomarine, marine, and aeolian sediments (Astafiev et al., Reference Astafiev, Bogdanov, Voinova, Zhuravlev, Nogina, Paramonova and Peshkova2012). The main geomorphological features on the Tersky Coast are primary shaped by the last glaciation, with late Valdaian (late Weichselian, MIS 2) tills and various glacial deposits dominating the landscape on the Kola Peninsula. Landforms such as fluting moraine, hummocky moraine, end moraine complexes, glacifluvial and glaciolacustrine toplands, eskers, and glacifluvial deltas have been identified through mapping efforts (Hättestrand and Clark, Reference Hättestrand and Clark2006; Boyes et al., Reference Boyes, Linch, Pearce, Kolka and Nash2021, Reference Boyes, Linch, Pearce and Nash2023). Coastal landforms in the late glacial and Holocene are largely composed of marine sediments and are found within a narrow coastal zone. The upper marine limit, which decreases from west to east, has been reconstructed to be approximately 35–40 m above sea level (m asl) and 10–12 m asl in the lower Varzuga and Chapoma Rivers, respectively (Kolka and Korsakova, Reference Kolka, Korsakova and Lisitsin2017). Additionally, Holocene aeolian landforms occur in the coastal zone. Some late Pleistocene glacial, glacioaqueous, and marine deposits, though localized, are buried beneath the late Valdaian (late Weichselian, MIS 2) and Holocene sediments.
Key sections with different types of deposits are exposed in the lower reach of the Varzuga, Chavanga, Kamenka, Chapoma, and Bolshaya Kumzhevaya Rivers, which all enter into the White Sea.
Methods
Field investigations
Fieldwork through the Tersky Coast of the White Sea included sedimentologic, lithostratigraphic, and other traditional geologic investigations of the sedimentary sequences in the outcrops, sampling for numerical dating, and unmanned aerial vehicle survey. Sedimentary sequences were sketched in the field and logged with a vertical resolution of 1:50. Sediments were classified according to the size of their components, as follows: (1) boulders: more than 100 mm; (2) pebbles: 10–100 mm; (3) gravel: 1–10 mm; (4) sand: 1.0–0.01 mm; and (5) silt and clay: <0.01 mm. The colour of deposits was identified according to the Munsell Soil Classification System. The samples for IR-OSL dating were collected according to a standard procedure by driving ca. 30-cm-long opaque plastic tubes into a cleaned face of the outcrop, then capping and sealing them to retain moisture. In total, six river and coastal sections were studied (Fig. 1). Elevation above sea level (asl) of the sections and lithostratigraphic units was measured using the 5× magnifying CB-17-630 Berger Hand Level and the Differential Global Positioning System total station.
IR-OSL dating
IR-OSL ages in this paper were obtained at the Research Laboratory for Quaternary Geochronology (RLQG) at Tallinn University of Technology.
The IR-OSL method uses potassium feldspar found in the Earth's crust to date sediments (Huntley et al., Reference Huntley, Godfrey-Smith and Thewalt1985; Godfrey-Smith et al., Reference Godfrey-Smith, Huntley and Chen1988). It relies on trapped charges (electrons) within the crystal structure of feldspar minerals, released when the grains are exposed to infrared radiation in the lab. The brightness of the luminescence signal corresponds to the natural ionizing radiation dose (palaeodose), correlating with the sedimentary age.
To determine a sample's age, accumulated palaeodose (P) and total natural radiation dose rate (D Σ) are considered. The age-correlating palaeodose has been accumulated in the mineral grains from the time the sedimentary material was last exposed to daylight, while the total dose rate includes contributions from cosmic radiation, sedimentary matrix, and the decay of 40K within the mineral. Potassium feldspar grains (100‒150 μm) were extracted from sediment using a specific procedure (Molodkov and Bitinas, Reference Molodkov and Bitinas2006). IR-OSL measurements were conducted at room temperature using a thermoluminescence (TL)/IR-OSL reader.
The dose–response curves were used to determine palaeodose through the multiple-aliquot additive-dose protocol. External beta and gamma contributions to the total dose rate (D Σ) were estimated in the lab using specific factors (Adamiec and Aitken, Reference Adamiec and Aitken1998). The in situ water content was considered. The contribution of cosmic dose rates to total dose rates was calculated using formulae provided by Prescott and Hutton (Reference Prescott and Hutton1994). The internal beta dose from decay of 40K and 87Rb was obtained based on recommended concentration estimates from Huntley and Baril (Reference Huntley and Baril1997) and Huntley and Hancock (Reference Huntley and Hancock2001). The K-feldspar samples were not preheated or thermally treated but stored for about 1 month at room temperature to reduce postirradiational phosphorescence (Molodkov and Bitinas, Reference Molodkov and Bitinas2006).
230Th/U dating
The 230Th/U dating was performed at the Laboratory for Polar and Marine Research, Earth Science Institute, St.-Petersburg University (Russia) on shell samples (Arctica islandica) from the marine deposits. The 230Th/U dating of shells includes the determination of uranium and thorium isotopes either in the whole carbonate skeleton or in its individual parts. The content of uranium and thorium isotopes in the external and internal fractions was quantified using radiochemical techniques (Maximov et al., Reference Maksimov, Sharin, Kuznetsov, Okunev, Grigoriev and Petrov2016). Often-observed variance in the distribution of uranium and thorium isotopes between the inner and outer layers of fossil shells can be attributed to their disparate postmortem uptake and migration into mollusc shells from the surrounding environment (Kaufman et al., Reference Kaufman, Chaleb, Wehmiller and Hillaire-Marsel1996; Labonne and Hillaire-Marcel, Reference Labonne and Hillaire-Marcel2000; Maksimov et al., Reference Maksimov, Sharin, Kuznetsov, Okunev, Grigoriev and Petrov2016; Ayling et al., Reference Ayling, Eggins, McCulloch, Chappell, Grün and Mortimer2017). According to the known equation (Kaufman and Broecker, Reference Kaufman and Broecker1965), the 230Th/U age and its error for the corresponding fractions were calculated from the values of activity ratios of thorium and uranium isotopes.
Results
Varzuga section
The Varzuga section is located on the right bank of the Varzuga River, 20 km upstream from its mouth, on the westernmost Tersky Coast (1 in Fig. 1). Here, strong river lateral erosion impacts an elongated hill with an elevation of 53–54 m asl at its peak. Its pitted surface is 15–25 m above adjacent areas. Due to the active slumping, the available sequence varies from year to year. Having been discovered in 1898 (Rippas, Reference Rippas1899), this outcrop, elongated for 2 km along an abrupt slope of the river valley, was subsequently and repeatedly investigated (Apukhtin, Reference Apukhtin1957; Lavrova, Reference Lavrova1960; Armand and Lebedeva, Reference Armand, Lebedeva, Grave and Armand1966; Nikonov, Reference Nikonov and Grichuk1966; Gudina and Yevzerov, Reference Gudina and Yevzerov1973; Yevzerov et al., Reference Yevzerov, Yelovicheva, Lebedeva, Rayamyae and Evzerov1981; Korsakova et al., Reference Korsakova, Molodkov and Kolka2004, Reference Korsakova, Molodkov, Yelovicheva and Kolka2019; Molodkov and Yevzerov, Reference Molodkov and Yevzerov2004; Lunkka et al., Reference Lunkka, Kaparulina, Putkinen and Saarnisto2018; Zaretskaya et al., Reference Zaretskaya, Korsakova, Molodkov, Ruchkin, Baranov, Rybalko, Lugovoy and Merkuliev2022).
The late Pleistocene sedimentary succession studied by us lies in the middle part of the Varzuga outcrop (66°23.77′N, 36°39.00′E) at 16–31 m asl, where corresponding deposits are superposed on middle Pleistocene marine and glacial sediments (Korsakova et al., Reference Korsakova, Molodkov, Yelovicheva and Kolka2019). The sedimentary log is shown in Figure 2 (section 1). The lower part of the studied sequence is composed of dark grayish-brown clay and silty clay, which are dated to 245 ± 27/21 ka by the 230Th/U method (section 1 in Fig. 2, Table 1). The late Pleistocene sequence begins with dark gray-brownish compact gravelly clayey silt, silt, and sand that occurred at 16.0–19.9 m asl; an IR-OSL date of 92.8 ± 14.4 ka was obtained from a sandy lens located 1.1 m above the lower contact of this unit. At 19.9–26.5 m asl, this depositional unit is overlain with an erosional contact by an alternation of yellow, olive-yellow, or gray, laminated, solid, fine-grained sand and sandy silt, which are horizontally, ripple-, wavy-, or cross-bedded, deformed upward with microlags, folds, streaked-out ripples, and fine silty partings similar to flames. The two IR-OSL dates of 69.5 ± 5.5 and 66.5 ± 7.9 ka were obtained from this unit at 22.6 and 25.5 m asl, respectively (Zaretskaya et al., Reference Zaretskaya, Korsakova, Molodkov, Ruchkin, Baranov, Rybalko, Lugovoy and Merkuliev2022; 1 in Fig. 2, Table 1). At 26.5–31.0 m asl in the studied sequence, there is a solid matrix-supported massive diamicton represented by reddish-brown clayey silt with gravel, pebbles, and boulders; its lower contact is sharp and even. The top of the Varzuga section is represented by light brownish-gray medium and fine sand with intercalations of gravel and pebbles.
Sediment logs from the studied outcrops. The locations of the sections are shown in Fig. 1.
Chavanga section
The Chavanga section (66°9.03′N, 37°47.13′E) is located on the southern Tersky Coast at an altitude of 28–37 m asl on the right bank of the Chavanga River, 4.1 km from its mouth (2 in Fig. 1).
In the sedimentary succession (2 in Fig. 2), gray compact silty mud was identified near the river water level at 28.5–29.0 m asl. At 29.0–30.9 (32.0) m asl, it is sharply overlain by coarse sand with gravel and well-rounded flattened pebbles. This unit is sharply overlain with very uneven contact by fine- to medium-grained yellowish sand including small lenses of pinkish-gray sand with fine horizontal lamination and ripple marks: there are pockets and interbeds of fine, medium, and coarse sand at the interval of ca. 30.9 (32.0) to 32.0 (32.4) m asl. The sand at 31.2–31.3 m asl is IR-OSL dated to 53.5 ± 4.1 ka (RLQG 2610-069) (Zaretskaya et al., Reference Zaretskaya, Korsakova, Molodkov, Ruchkin, Baranov, Rybalko, Lugovoy and Merkuliev2022; 2 in Fig. 2, Table 1). At ca. 32.0 (32.4) to 35.9 m asl, yellowish-brown sand with indistinct lamination and spots of grayish sand and lenses of sandy silt, probably glaciotectonically deformed, occurs; IR-OSL dates show these deposits accumulated no earlier than 18.4–17.8 ka (Zaretskaya et al., Reference Zaretskaya, Korsakova, Molodkov, Ruchkin, Baranov, Rybalko, Lugovoy and Merkuliev2022; 2 in Fig. 2). At the top of the sequence, compact poorly sorted sandy and clayey silt with gravel, pebbles, and small boulders occurs at 35.9–36.2 m asl and is overlain with trough-shaped contact by poorly sorted sand with gravel at the interval of 36.2–37.0 m asl.
Kamenka section
The Kamenka section (66°5.83′N, 38°16.65′E) is located at an altitude of 44–46 m asl on the abrupt right bank of the Kamenka River valley, 4.7 km from its mouth (3 in Fig. 1).
The basal part of the sediment succession (3 in Fig. 2) consists of dark grayish-brown silty clay with boulders up to 30–40 cm in diameter, abundant pebbles, and gravel; local fissility appears as indistinct horizontal thin partings, that is, the diamict splits into numerous lens-shaped plates, a few millimetres to one or few centimetres thick. Light brownish-gray fine sand and silt with sharp uneven contact at 40.8–40.9 m asl lies immediately on the clayey sediments; ferruginous rust stains and dark-gray smudges are common for the sand. The IR-OSL method has yielded the age of fine sand and silt to be 40.1 ± 4.6 ka (RLQG 2611-079) (3 in Fig. 2, Table 1). Overlying the fine sand is a lens-shaped unit composed of brown very fine sand with abundant gravel. A pocket of washed fine sand IR-OSL dated to 37.2 ± 4.3 ka (RLQG 2612-079) (3 in Fig. 2, Table 1) occurs at 41.8–42.6 m asl. At 42.6–46.5 m asl, dark-brown silty clay with abundant boulders, pebbles, and gravel was exposed; sandy lenses and streaks are common in this unit. The upper silty clay is moist and massive with gritty texture; included gravel and pebbles are well-rounded and evenly distributed, and some are eroded. Hair-like inclusions similar to aquatic plant remains are visible in these sediments. At 44.0 m asl, a streak of reddish-brown, poorly sorted, predominantly medium sand occurred and is dated to 58.2 ± 3.9 ka (RLQG 2662-070). A pocket of greenish-black silty clay with broken shells and well-rounded small boulders, pebbles, and gravel overlies the lens-shaped sandy texture. A reddish-dark brown sand–gravel–pebble streak of 3–15 cm depth occurs at 44.6–45.1 m asl. These sediments are IR-OSL dated to 54.4 ± 3.7 ka (RLQG 2663-070) (3 in Fig. 2, Table 1). The IR-OSL-dated sandy streaks are suggested to be inclusions in the clayey unit, which occurs at 42.6–46.5 m asl (3 in Fig. 2). Above the streak up to the top, dark-brown silty clay with abundant boulders, pebbles, and gravel is represented in the investigated sequence.
Chapoma 1 section
The Chapoma 1 section (previously published as Chapoma) is located on the abrupt left slope of the Chapoma River valley, 3.5 km from its mouth (4 in Fig. 1), at an altitude up to 29–30 m asl.
In the Chapoma 1 sequence (66°6.93′N, 38°51.92′E), the lowermost weak red compact silty clay with abundant boulders, pebbles, and gravel of mainly local rocks and scarce erratic pebbles is exposed above the river shoreline, located here at ca. 10.0 m asl, and above it up to 11.05 m asl (4 in Fig. 2). The upper sharp contact, identified according to colour and changes in grain composition, dips steeply to the west. On this diamict, moisture-laden greenish-gray compact silt and silty clay (mud) with numerous shells and shell fragments, scarce pebbles, and gravel occurs at 11.5–13.6 m asl. Collected during latest study of the section from the mud, the A. islandica shell has been 230Th/U dated to 136 ± 15/13 ka (LU-479b) (4 in Fig. 2, Table 1). Up the section, this mud is covered with silt and fine sand including abundant pebbles. Upper contact of this unit is uneven, with pockets of greenish-gray silty clay and yellowish-brown, poorly sorted, predominantly coarse sand with well-rounded pebbles; isolated weathered pebbles are represented here. Yellowish-brown, compact, poorly sorted sand with abundant boulders up to 25 cm in diameter, pebbles, and gravel occurs at 13.6–14.75 m asl. These mixed pebbly–gravelly–sandy sediments are overlain with dark grayish-brown compact silt and clayey silt with abundant rubble, differently rounded pebbles, and scarce shell fragments with a sharp and uneven contact. Above 14.75 m asl, dark grayish-brown silt gradually passes into grayish-brown and greenish-gray, weakly laminated, compact silt and sand with scarce gravel and pebble. At 14.75–17.5 m asl, thin (2–3 mm) splitting sandy laminae and scarce shell fragments are included in the silty sediments. Compact structureless clayey silt, silt, and sand with solitary rounded pebbles of local rocks and weathered far-travelled carbonate rocks occur at 17.5–18.3 m asl. At 18.3–20.0 m asl, compact silt and clayey silt with gravel and pebbles assume a special structure, breaking up into separate angular fragments; scarce fragments of shells are again visible in these deposits. The upper contact of these deposits is very uneven and sharp. At 20.0–23.5 m asl, brownish-gray, brownish-yellow, and pale-olive sands overlie the clayey silts with disconformity. Individual fragments of shells and flattened pebbles are visible in the lower part of this sand, which is IR-OSL dated to 29.1 ± 3.4 ka (RLQG 2599-069) (4 in Fig. 2, Table 1). In the sand, a deformed parting of compact dark grayish-brown silt and clayey silt occurs at 21.2–21.5 m asl. In the lower 1.5 m of this member, sands with unclear wavy bedding include pockets and angular-shaped fragments of this silt and clayey silt, as well as the broken shell fragments. A silty-clayey diapir-like fold with a disturbed fine-grained sandy core and small pockets including angular silty-clay fragments and silt is in the northern part of the outcrop. At 22.3–23.5 m asl, strongly deformed sand, silt, and clayey silt, including lenses and pockets with gravel and pebbles, as well as pockets of ripple bedding sand with individual very small shell fragments are presented. Two different IR-OSL dates were obtained from this deformed part of the Chapoma 1 section, namely, 23.0 ± 2.7 ka (RLQG 2600-069) from the sand at 22.5 m asl and 76.7 ± 8.9 ka (RLQG 2601-069) from the sand at 23.3 m asl (4 in Fig. 2, Table 1). Upward, dark olive-gray homogenized diamicton represented by compact clayey silt with boulders, pebbles, and gravel lies on the deformed sandy deposits with sharp uneven contact at 23.5–27.5 m asl. The top of the Chapoma 1 outcrop is covered by turf up to 29.0 m asl.
Chapoma 2 section
The Chapoma 2 section is located 0.05 km from the Chapoma River mouth at the junction of the abrupt left slope of the river valley and an eroded sea cliff up to 23.8 m asl (Zaretskaya Mount) (5 in Fig 1). This section was previously published as Faktoria (Zaretskaya et al., Reference Zaretskaya, Korsakova, Molodkov, Ruchkin, Baranov, Rybalko, Lugovoy and Merkuliev2022).
The Chapoma 2 sedimentary sequence was observed in two trenches spaced 40 m apart (5S and 5N in Fig. 2). Basal sedimentologic units were studied in the trench located in the southern part of the outcrop (66°6.22′N, 38°53.18′E) facing the White Sea shore (5S in Fig. 2). At ca. 6.0 m asl, a basal shingle and boulders are covered with compact very dark gray clayey silt, which gradually changed to dark-gray clayey silt and silt up to 8.9 m asl. These deposits contain scarce shell detritus, gravel, and pebbles. In the uppermost part of this unit, 0.10- to 0.12-m-length lenses and pockets with gray wavy-laminated fine sand occur. The upper contact, detected by particle-size distribution and colour, is distinct and uneven, with fragments of clayey silt intercalated into overlying sediments. Above this erosional contact marked with well-rounded pebbles, loose, indistinctly horizontally laminated, reddish-brown, poorly sorted, predominantly coarse sand with red arkosic sandstone cobbles, well-rounded gravel, and very small shell fragments were observed. In the southern trench, the apparent thickness of these sediments is ca. 1 m; upward, talus covers the coastal slope. The sand sample taken from basal part was IR-OSL dated to 95.3 ± 8.1 ka (RLQG 2661-070) (Zaretskaya et al., Reference Zaretskaya, Korsakova, Molodkov, Ruchkin, Baranov, Rybalko, Lugovoy and Merkuliev2022; 5S in Fig. 2, Table 1).
Other sedimentary units of the Chapoma 2 sequence were studied in the northern part of the outcrop (66°6.27′N, 38°53.15′E) facing the river valley (5N in Fig. 2). Here, dark grayish-brown fine sand and silt with shell fragments and small pebbles is exposed in the basal part of the sequence; festooned flame-like structures, seams, and pockets formed by gray fine sand or reddish-gray medium sand occur here. A sand sample at 7.5 m asl was IR-OSL dated to 72.4 ± 5.6 ka (RLQG 2604-069) (Zaretskaya et al., Reference Zaretskaya, Korsakova, Molodkov, Ruchkin, Baranov, Rybalko, Lugovoy and Merkuliev2022; 5N in Fig. 2, Table 1). These lowermost deformed sandy and silty deposits were observed here at 6.0–8.5 m asl in the form of a lens, which seems to be leaning against the clayey silt, silt, and sand exposed in the Chapoma 2 south section (5S in Fig. 2). Along the stretch of outcrop, these deformed sediments in section 5N replace the clayey silt, silt, and sand visible in section 5S. Upper contact of this member is not clearly visible due to thick talus; its location was roughly identified according to implicit data. Upward, the thick sandy unit, identified at 8.5–14.15 m asl, is composed of gray and reddish-gray fine and medium sand with shell fragments and displays horizontal and wavy lamination throughout the entire thickness. In its upper part, lenses and interbeds of gravel and coarse sand with shell fragments, sparse flattened arkosic sandstone cobbles, and pockets of horizontal, weakly laminated, fine-grained sand are visible at 14.0–14.15 m asl. Above them at 14.15–16.8, predominantly wavy or oblique laminated sand was IR-OSL dated to 36.2 ± 2.8 ka (RLQG 2603-069) and 32.5 ± 3.8 ka (RLQG 2602-069) (Zaretskaya et al., Reference Zaretskaya, Korsakova, Molodkov, Ruchkin, Baranov, Rybalko, Lugovoy and Merkuliev2022; 5N in Fig.2, Table 1). At 16.8–17.5 m asl, compact dusky red with thin interbeds of dark-gray fine sand occurs; the lower contact is clear, uneven, and identifiable by colour and grain-size composition. At 17.5–23.8 m asl, a compact dark reddish-gray diamicton, represented by silt with gravel and pebbles and having a platy structure, turning upwards into massive clayey silt, is observed; an interbed of pale-brown, poorly sorted, predominantly fine sand with 0.2-m-thick included well-rounded gravel and pebbles occurs in the middle part of this unit. The diamicton overlaps the underlying sands, with a gradual transition marked with a pocket of fine sand.
Bolshaya Kumzhevaya section
The Bolshaya Kumzhevaya section (66°13.60′N, 39°40.48′E) is located on the steep left bank of the Bolshaya Kumzhevaya River valley, 3.1 km from the mouth of the river at an altitude of 16–27 m asl (6 in Fig. 1). The natural outcrop exposes a wide terrace composed of marine and glacial deposits and is located 200 m from its bench, which has a steepness of about 45° and a relative height of about 10 m.
In the basal part of the Bolshaya Kumzhevaya sequence (6 in Fig. 2), water-bearing grayish-green fine sand and silt with gravel and rounded pebbles occur at 16.0–18.0 m asl; the complanate large pebbles and scarce boulders are oriented subhorizontally. Upper contact is sharp, identified by colour and grain composition. At 18.0–18.5 m asl, subbasal grayish-brown 0.1-m-thick silt with gravel and pebbles is overlain by water-bearing black and dark-gray, 0.1-m-thick, poorly sorted, predominantly medium sand with cobbles located under pockets of coarse sand and gravel, which gradually change upward into dark grayish-brown and brown, 0.3-m-thick medium sand with pockets of coarse sand, scarce gravel, and cobbles. The sand at 18.25 m asl was IR-OSL dated to 95.9 ± 6.3 ka (RLQG 2660-070) (6 in Fig. 2, Table 1). The thick sandy unit, identified at 18.5–23.5 m asl, is composed of three beds. The lower one is represented by light yellowish-brown 0.7-m-thick thin sand with horizontally weak lamination and differently directed oblique lamination underlying light brownish-gray, weakly laminated, thin sand and silt of 0.4 m thick; a sand sample at 19.6 m asl was IR-OSL dated to 57.6 ± 6.7 ka (RLQG 2607-069) (6 in Fig. 2, Table 1). Upward, this sand changes to 0.7-m-thick massive grayish-brown thin sand with scarce gravel and small pebbles, which is overlaid by 1.7-m-thick, horizontal, weakly laminated, thin sand IR-OSL dated to 55.5 ± 5.1 ka (RLQG 2606-069) at 21.55 m asl (6 in Fig. 2, Table 1). The uppermost bed is composed of horizontal or wavy weakly laminated thin sand with thin interbeds and pockets of light brownish-gray silt and medium sand. The upper erosional contact of this sandy layer is clear, uneven, and identifiable by colour and grain-size composition. Upward, reddish-gray diamicton represented by compact silty clay splitting into indistinct horizontal thin plates and clayey silt with boulders, pebbles, and gravel lies on the sandy member with sharp uneven contact at 23.5–27.0 m asl. A 0.25-m-thick parting of pinkish-gray, wavy laminated, fine sand and silt, IR-OSL dated to 15.4 ± 1.8 ka (RLQG 2605-069) (6 in Fig. 2, Table 1), with erosional sharp contacts, occurs at 23.9 m asl.
Discussion
The reliability of the dating results
To validate the reliability of the IR-OSL dates presented in this paper (Table 1), we provide a series of comparative results obtained earlier using two methods: IR-OSL and ESR (Fig. 3, Table 2). These results mainly stem from parallel dating using ESR and IR-OSL techniques on shells and feldspar grains, respectively, extracted from the same sediment sample. Our comparative results showcase the potential of both combined and independent palaeodosimetric dating methods applied to two different minerals— biogenic carbonate and feldspar—to chronologically organize middle to late Pleistocene palaeoenvironmental events.
Comparison of parallel dating results obtained by two main independent dating methods used in the present study: electron spin resonance (ESR) and infrared optically stimulated luminescence (IR-OSL). The figure indicates the ages (in ka).
Results of parallel dating of the same marine sediment sample using mollusc-based electron spin resonance (ESR) and feldspar-based infrared optically stimulated luminescence (IR-OSL) methods.
a Sampling locations: samples 2, 3, and 10: Kola Peninsula; samples 1, 8, and 11: the easternmost coastal area of the Barents Sea; samples 4–7 and 9: the Kara Sea coastal area.
The results obtained are particularly significant due to frequent references in the literature to challenges such as incomplete bleaching of sediment grains. Such incomplete bleaching may result from factors like transport in turbid waters with high suspended particle concentrations or reworking of initial subaqueous sediments. These factors, or similar ones, would be expected for marine sediments, potentially leading to age overestimations. Additionally, it is assumed that alkali feldspars exhibit specific athermal (anomalous) fading of their luminescence, which could hinder the reliable luminescence dating of Quaternary deposits based on feldspar extractions, potentially causing unpredictable age underestimations. For instance, Wallinga et al. (Reference Wallinga, Murray, Duller and Tornqvist2001) reported a nearly 50% underestimation of feldspar IR-OSL ages compared with quartz OSL ages measured on the same samples, presumably due to anomalous fading.
The reliability of the palaeodosimetric dating techniques employed in this study is also demonstrated through the series of comparative results obtained using a series of independent numerical dating methods, such as K-feldspar-based IR-OSL, mollusc shell–based ESR, quartz-based optically stimulated afterglow, U-Th, and 14C, applied to the same sedimentary samples (Molodkov, Reference Molodkov2020), as well as relative dating methods (Molodkov, Reference Molodkov2012). A notable consistency observed in dating results from feldspar and quartz grains, as well as mollusc shells taken from the same sedimentary context, indicates no noticeable athermal fading of feldspar dosimetric traps over geologic timescales, as detailed in previous studies (Vasilchenko et al., Reference Vasilchenko, Molodkov and Jaek2005; Molodkov et al., Reference Molodkov, Jaek and Vasilchenko2007; Jaek et al., Reference Jaek, Molodkov and Vasilchenko2008, Reference Jaek, Molodkov and Vasilchenko2010), nor does it indicate incomplete bleaching of sediment grains, even under complex conditions expected in marine environments.
The results obtained using the 230Th/U age radioisotope method on subfossil thick-walled, well-preserved mollusc shells are somewhat questionable due to the potential open-system behaviour of subfossil shells with respect to uranium (see Kaufman et al., Reference Kaufman, Chaleb, Wehmiller and Hillaire-Marsel1996). In the postsedimentation period, the shells can be subject to both secondary uranium enrichment, as a result of their contact with marine or ground water, and leaching of uranium and thorium compounds. In addition, it is possible for mineral suspension to penetrate into shells, leading to thorium contamination, which is detected by the presence of the 232Th isotope (Choukri et al., Reference Choukri, Hakam, Reyss and Plaziat2007; Rowe et al., Reference Rowe, Turner, Andrews, Leeder, Calsteren and Thomas2015). Thus, the reliability of 230Th/U age, depending on the local conditions of shell fossilization, is not an accurate radioisotope chronometer and can rather be used to approximate the age of the enclosing sediments (Labonne and Hillaire-Marcel, Reference Labonne and Hillaire-Marcel2000; Choukri et al., Reference Choukri, Hakam, Reyss and Plaziat2007).
Interpretation of sedimentary successions and depositional environments
In accordance with the regional subdivisions of the late Pleistocene derived from previous reconstructions (Molodkov and Bolikhovskaya, Reference Molodkov and Bolikhovskaya2002; Korsakova et al., Reference Korsakova, Molodkov and Kolka2004, Reference Korsakova, Kolka and Semenova2016, Reference Korsakova, Molodkov, Yelovicheva and Kolka2019; Korsakova, Reference Korsakova2019, Reference Korsakova2021; Zaretskaya et al., Reference Zaretskaya, Korsakova and Panin2019, Reference Zaretskaya, Korsakova, Molodkov, Ruchkin, Baranov, Rybalko, Lugovoy and Merkuliev2022), the stratigraphic sequence is subdivided into four informal lithostratigraphic formations (units and subunits). These units are U1 (MIS 5 warming) and corresponding marine (m) Ponoy (U1a) and Strelna (U1b) beds; then U2 (severe cooling correlated to MIS 4) and U3 (a time interval with complex and variable environment correlated to MIS 3) (see, e.g., Kind, Reference Kind1974; Bolikhovskaya and Molodkov, Reference Bolikhovskaya, Molodkov, Johansson and Sarala2007), represented mainly by glaciomarine (gm) and marine (m) formations, as well as glacial (g) and glaciofluvial (f) deposits; and U4 (MIS 2 cooling) represented by deposits interpreted as glacial (g), glaciofluvial (f), glaciolacustrine (l), and glaciomarine (gm). In the studied sequences, Unit H includes mainly the Holocene marine, aeolian, and other terrestrial surficial deposits; middle Pleistocene deposits known in the Kola Peninsula correlate with Units MP 1 (MIS 9) and MP 2 (nonsegmented glacial–interglacial beds, MIS 8–6) (Korsakova et al., Reference Korsakova, Molodkov, Yelovicheva and Kolka2019).
Varzuga sequence
In the Kola region, sedimentary sequences exposed in the Varzuga outcrop (Figs. 2 and 4) represent the most complete and longest record of late Pleistocene history over the MIS 5–4 stages. Earliest late Pleistocene dark brownish-gray mud with dominating clayey silt includes shells and shells fragment (Unit 1a), and overlying laminated, cross-bedded, sandy and silty sediments (Unit 1b) were described in the Varzuga outcrop. Previous investigations (Korsakova et al., Reference Korsakova, Molodkov and Kolka2004, Reference Zaretskaya, Korsakova and Panin2019; Molodkov and Yevzerov, Reference Molodkov and Yevzerov2004) showed that along the stretch of the outcrop downstream of the Varzuga River, the Unit 1b sandy and silty sediments are partly replaced by clayey silt containing mollusc shells (Unit 1a) at 25–31 m asl (Fig. 4, sections a1 and a2). The weighted mean of several ESR age determinations on different shells from the Unit 1a clayey silt was estimated at 103.0 ± 4.2 ka (RLQG 309-12-042), which agrees with the age of 104.0 ± 8.3 ka (RLQG 1405-031) obtained by IR-OSL on potassium feldspar extracted from the enclosing deposits (Molodkov and Yevzerov, Reference Molodkov and Yevzerov2004). The Unit 1a sediments are also suggested in Varzuga section b3 (Fig. 4) in accordance with its stratigraphic position in the sequence presented by Lunkka et al. (Reference Lunkka, Kaparulina, Putkinen and Saarnisto2018). The lowermost part of the sandy Unit 1b was later dated by IR-OSL at 95.9 ± 11.1 ka (Korsakova et al., Reference Korsakova, Molodkov, Yelovicheva and Kolka2019; Fig. 4, section a1) and by OSL at 88 ± 7 ka (Lunkka et. al., 2018; Fig. 4, section b2).
Late Pleistocene sedimentologic successions from the Varzuga sequence: sections 1a and 1b derived from Korsakova et al. (Reference Korsakova, Kolka and Semenova2016); sections b1, b2, and b3 derived from Lunkka et al. (Reference Lunkka, Kaparulina, Putkinen and Saarnisto2018); section c presented in this paper. Date samples and dating results: 1, infrared optically stimulated luminescence (IR-OSL) (this paper); 2, 230Th/U (this paper); 3, electron spin resonance (ESR) (Korsakova et al., Reference Korsakova, Molodkov and Kolka2004, Reference Korsakova, Molodkov, Yelovicheva and Kolka2019; Molodkov and Yevzerov, Reference Molodkov and Yevzerov2004); 4, IR-OSL (Molodkov and Yevzerov, Reference Molodkov and Yevzerov2004; Korsakova et al., Reference Korsakova, Molodkov, Yelovicheva and Kolka2019); 5, OSL (Lunkka et al., Reference Lunkka, Kaparulina, Putkinen and Saarnisto2018); 6, 230Th/U (Arslanov et al., Reference Arslanov, Yevzerov, Tertychnyy, Gerasimov, Lokshin, Velichko and Faustova1981). The main subdivisions are designated as units correlating with the Marine Isotope Stages (MIS) according to presented data: MP (MIS 9), MP 2 (MIS 8–6), U1a and U1b (MIS 5 with corresponding Ponoy and Strelna beds, respectively), U2 (MIS 4), U3 (MIS 3), U4 (MIS 2), and H (MIS 1).
The juxtaposing of lithologic and chronological data from several sediment successions studied by many researchers (Fig. 4, Table 1) provide evidence that the late Pleistocene depositional sequences in the Varzuga sections show four formations (Units 1a, 1b, 2, and 4 in Fig. 4); three sharp contacts suggested erosional events throughout the depositional history. The Unit 1a sediment succession with clayey and silty deposits including warm-water molluscan fauna, as well as foraminifera and diatoms (Gudina and Yevzerov, Reference Gudina and Yevzerov1973), that is, Ponoy beds in the regional stratigraphy (Korsakova, Reference Korsakova2019), is locally present in the Varzuga sections and has been identified in a few points at ca. 25–31 and 20–24 m asl (Fig. 4, sections a2 and b3, respectively), indicating its lens-shaped bedding. Results of multimethod chronological investigations (Molodkov and Yevzerov, Reference Molodkov and Yevzerov2004) show that Unit 1a was formed here at 104–103 ka (Fig. 4, section a2, Table 1), which corresponds to further development of the interglacial warm-water Boreal marine transgression in the White Sea depression. Unit 1b is correlated to one or more interglacial marine events in the Kola Peninsula known in the regional stratigraphic scheme as the Belomorian transgression (Lavrova, Reference Lavrova1960; Korsakova et al., Reference Korsakova, Kolka and Semenova2016, Reference Korsakova, Molodkov, Yelovicheva and Kolka2019), which is suggested to be the final phase of the Boreal marine transgression (about 110–70 ka). In the Varzuga sequence, sandy and silty sediments of Unit 1b dated to 95.9, 92.8, and 88.0 ka (Fig. 4, sections a1 and b2, Table 1) sharply overlie silty sediments of Unit 1a (Ponoy beds) in section b3; their composition suggests the ca. 1.0-m-thick transgression succession at 24–25 m asl and a thicker regression at 25–32.5 m asl. In sections a1, b1, b2, and c (Fig. 4), the Unit 1b deposits sharply overlie the middle Pleistocene sediments with erosion. Although there is the opinion (Lunkka et. al., 2018) that the basal part of this unit (i.e., lower 1.0-m-thick transgression succession of Unit 1a in section b3 in Fig. 4) represents subglacially deformed sediments or mass movement diamicton conformably overlain by laminated silt (i.e., several meters thick regression succession in section b3 in Fig. 4) accumulated in a basin with glaciolacustrine conditions, we suggest their marine origin because of lithology and lithostratigraphic position and correlation to other sites on the Tersky Coast. The MIS 4 deposits were identified because of the new geochronological data from the uppermost part of Unit 2 in section c, which was dated to 69.5 ± 5.5 (RLQG 2613-079) and 66.5 ± 7.9 ka (RLQG 2614-079), as well as from basal part of Unit 2 in section b2 dated to 70 ka (Lunkka et al, Reference Lunkka, Kaparulina, Putkinen and Saarnisto2018; Fig. 4, Table 1). Unit 2, which is composed of sandy deposits, should have consistently accumulated during a cold environment, which, according to geochronological data, is comparable to the cold MIS 4. The origin of these sandy deposits can only be assumed with a high degree of uncertainty. It is most likely that the Unit 2 sediment accumulated in unstable coastal and shallow-water environments, possibly in a marine setting, for example, during the Mezen transgression, which developed in this region during MIS 4 (Zaretskaya et al., Reference Zaretskaya, Korsakova, Molodkov, Ruchkin, Baranov, Rybalko, Lugovoy and Merkuliev2022). The late Valdaian (late Weichselian, MIS 2) clayey silt of Unit 4 overlaps sandy Unit 2 with sharp contact and erosion (Fig. 4), which demonstrates a hiatus in the middle–late Valdaian (middle–late Weichselian) interglacial–glacial sediment succession. In the Varzuga sections, the observed sedimentary sequences show listric folds, thrust planes, and other folded structures in the sandy deposits underlying the Unit 4 sediments; these structures are associated with the glaciotectonic impact of the SIS (Lunkka et. al., Reference Lunkka, Kaparulina, Putkinen and Saarnisto2018; Korsakova et al., Reference Korsakova, Molodkov, Yelovicheva and Kolka2019). The uppermost Unit H sediments were attributed to the end of the late glacial and Holocene, suggesting their lacustrine and aeolian origin (Lunkka et. al., Reference Lunkka, Kaparulina, Putkinen and Saarnisto2018). In the lower Varzuga River area, the vast aeolian landscapes are presented now in the coastal zone.
Chavanga sequence
Units 1a–2–3–4 sedimentologic succession was identified in the Chavanga outcrop (Figs. 2 and 5). Previously investigated (Gudina and Yevzerov, Reference Gudina and Yevzerov1973; Korsakova et al., Reference Korsakova, Molodkov and Kolka2004), section a (Fig. 5) demonstrates the Unit 1a deposits with fragments and valves of subfossil shells of arctic, arctic–boreal, and boreal mollusc species ESR dated to 99.0 ± 7.6 ka (Fig. 5, Table 1). Unit 1a was attributed to Ponoy beds accumulated in the end of initial phase of the Ponoy transgression, which corresponds to the first phase of the Boreal transgression. These sediments are sharply overlaid with erosion by the Unit 2 pebbles and gravels in a sandy and silty matrix (subunit 1, section a in Fig. 5) covered by laminated sands with rounded pebbles (subunit 2, section a in Fig. 5). Here, subunit 2 of Unit 2 (section a in Fig. 5) is similar to Unit 2 in the Varzuga sequence, specifically section c (Fig. 4), because of its lithology and chronology, suggesting similar settings during their accumulation in MIS 4 under marine environments. The 53.5 ± 4.1 ka date (section b in Fig. 5, Table 1) could indicate a longer accumulation stage in the marine environment, embracing part of MIS 4 and the onset of MIS 3. Subunit 2 of Unit 2, dated to 63.65 ± 8.0 ka (section a in Fig. 5, Table 1), and laminated sand of Unit 3 (section b in Fig. 5) have a marine origin and accumulated in unstable coastal and shallow-water environments. Based on lithologic composition, it can be interpreted that subunit 1 of Unit 2 could have accumulated in MIS 4 in a glaciofluvial setting. The late Valdaian (late Weichselian, MIS 2) Unit 4 is subdivided into three subunits (section b in Fig. 5). Dated to 18.1–17.8 ka, sand of subunits 1 with folded structures and thrust planes can be identified as deformed till that was deposited during the last glacial/ice sheet advance. Subunits 2 and 3 of Unit 4 can be interpreted as supraglacial allochthonous sediments.
Late Pleistocene sedimentologic successions from the Chavanga sequence: section a derived from Korsakova et al. (Reference Korsakova, Molodkov and Kolka2004); section b presented in this paper. See Fig. 4 legend for symbols.
Kamenka sequence
A reduced Units 2–3–4 sedimentary sequence is represented in the Kamenka outcrop (Figs. 2 and 6). This sedimentary succession was previously investigated (Gudina and Yevzerov, Reference Gudina and Yevzerov1973; Korsakova et al., Reference Korsakova, Molodkov and Kolka2004; Astafiev et al., Reference Astafiev, Bogdanov, Voinova, Zhuravlev, Nogina, Paramonova and Peshkova2012). The silty clay and sand of marine units sandwiched here with sharp contacts between the diamicton indicated by pebbles and gravel in a clayey, silty, or sandy matrix (section a in Fig. 6) were recognized. As in other sequences on the Tersky Coast, Unit 2 is correlated to MIS 4. Here, Unit 2 consists of a glacial diamicton, which is superposed by marine sediments overlapping it with erosion. The glacial origin of the Unit 2 gravelly silty clay was recognized according to its lithologic composition and lithostratigraphic position below the silty clay, silt, and sand with mollusc shells and marine diatoms (Gudina and Yevzerov, Reference Gudina and Yevzerov1973). Fragments of subfossil shells of arctic–boreal and arctic mollusc species, that is, Mya sp., Chlamys islandica, Astarte crenata, and Astarte crebricostata, as well as marine diatoms Isthmia sp., Isthmia nervosa, Rhabdonema sp., and Paralia sulcata (Grave et al., Reference Grave, Gunova, Devyatova, Lavrova, Lebedeva, Samsonova, Cheremisinova, Grave and Koshechkin1969; Korsakova, Reference Korsakova2021), were found in silty clay, silt, and sand corresponding to Unit 3 (section a in Fig. 6). The Unit 2 diamicton, demonstrating structures in the form of indistinct horizontal thin partings into thin microplate and underlying marine sediment with uneven sharp contact, appears to be a basal till and can be correlated with the MIS 4 glacial event on the Kola Peninsula. Accumulation of the Unit 3 marine sediments, dated to 37.2 ± 4.4 ka (RLQG 2612-079) and 40.1 ± 4.6 ka (RLQG 2611-07) (section b in Fig. 6, Table 1), suggests a cold-water basin with a changing coastline during the time interval 58–37 ka (Fig. 6, Table 1) correlated to MIS 3. No thermophilic diatom species have been found here; the mollusc shells from the lower part of this unit were ESR dated to 58.7 ± 4.4 (RLQG 344-073) and 52.0 ± 4.3 (RLQG 416-119) ka (Korsakova et al., Reference Korsakova, Molodkov and Kolka2004; Astafiev et al., Reference Astafiev, Bogdanov, Voinova, Zhuravlev, Nogina, Paramonova and Peshkova2012; Fig. 6, Table 1). Spores from these sediments show that a forest–tundra with sparse growth of Betula–Pinus trees covered the adjacent land (Grave et al., Reference Grave, Gunova, Devyatova, Lavrova, Lebedeva, Samsonova, Cheremisinova, Grave and Koshechkin1969). The uppermost part of the Kamenka sequence, that is, Unit 4, appears to be waterfront glacial sediments represented by glaciomarine deformed till with pockets of allochthonous inclusions originated most probably from underlying marine deposits, IR-OSL dated to 58.2 ± 3.9 ka (RLQG 2662-070) and 54.4 ± 3.7 ka (RLQG 2663-070) (Fig. 6, Table 1). The top position of Unit 4 in the sequence indicates its late Valdaian (late Weichselian, MIS 2) age.
Late Pleistocene sedimentologic successions from the Kamenka sequence: section a derived from Korsakova et al. (Reference Korsakova, Molodkov and Kolka2004); section b presented in this paper. See Fig. 4 legend for symbols.
Chapoma 1 sequence
In the Chapoma 1 (Figs. 2 and 7), Units 1a–1b–2–4 sediment succession was suggested. In the twentieth century, this sequence was studied by Lavrova (Reference Lavrova1960), Malyasova (Reference Malyasova and Pokrovskaya1960), Cheremisinova (Reference Cheremisinova1962), Grave et al. (Reference Grave, Gunova, Devyatova, Lavrova, Lebedeva, Samsonova, Cheremisinova, Grave and Koshechkin1969), Gudina and Yevzerov (Reference Gudina and Yevzerov1973), and Arslanov et al. (Reference Arslanov, Yevzerov, Tertychnyy, Gerasimov, Lokshin, Velichko and Faustova1981). Corresponding litho-, bio-, and chronostratigraphic data were summarized in Korsakova (Reference Korsakova2019, Reference Korsakova2021). The data previously obtained demonstrate that compact silty clay with gravel, pebbles, and boulders and scarce marine foraminifera (Grave et al., Reference Grave, Gunova, Devyatova, Lavrova, Lebedeva, Samsonova, Cheremisinova, Grave and Koshechkin1969) comprises the lowermost part of the sequence (Unit MP 2 in Fig. 7). Upwards with a sharp contact, silt and silty clay (Chapoma 1 section, Unit 1a in Fig. 7) with abundant subfossil mollusc shells of A. islandica and Astarte elliptica ESR dated at 126.2 ± 10.1 ka (RLQG 267-129) (Molodkov and Bolikhovskaya, Reference Molodkov and Bolikhovskaya2002), 134.5 ± 9.5 ka (RLQG 395-039A) and 137.0 ± 9.6 ka (RLQG 395-039B) (Molodkov, Reference Molodkov2012) were identified. The lower part of this unit was 230Th/U dated at 136 ± 15.13 ka and the upper part at 86.0 ± 3.9 ka (Chapoma 1 section in Fig. 7, Table 1). The latter is a very preliminary date, because this was the first attempt to date shells by 230Th/U method in the Kola Peninsula.
Late Pleistocene sedimentary successions from the Chapoma 1, Chapoma 2, and Bolshaya Kumzhevaya sequences. See Fig. 4 legend for symbols.
According to available data, the Unit 1a compact silt and silty clay with numerous shells and shell fragments, scarce pebbles, and gravel started to accumulate at ca. 138–128 ka and can be attributed to the initial phase of the interglacial Boreal transgression. The abundant fauna and diatoms suggest hydrobiological conditions in that marine basin to be more favourable than those of the present-day White Sea (Cheremisinova, Reference Cheremisinova1962; Grave et al., Reference Grave, Gunova, Devyatova, Lavrova, Lebedeva, Samsonova, Cheremisinova, Grave and Koshechkin1969; Gudina and Yevzerov, Reference Gudina and Yevzerov1973; Korsakova, Reference Korsakova2021), likely due to the influence of the warm North Atlantic Current (Zatsepin and Filatova, Reference Zatsepin and Filatova1961; Funder et al., Reference Funder, Demidov and Yelovicheva2002; Grøsfjeld et al., Reference Grøsfjeld, Funder, Seidenkratz and Glaister2006). The Unit 1b sediments with shells superpose the Unit 1a with a sharp erosional contact, suggesting some gap in the marine sedimentation. Successively overlaid, pebbly and gravelly sand (subunit 1 of Unit 1b in the Chapoma 1 sequence in Fig. 7), a thin layer of silt (subunit 2 of Unit 1b in the Chapoma 1 sequence in Fig. 7), and weakly laminated silty-sandy deposits (subunit 3 of Unit 1b in the Chapoma 1 sequence in Fig. 7) appear to be the attributable to the latter half of the Boreal transgression (Molodkov and Yevzerov, Reference Molodkov and Yevzerov2004), which can be correlated with MIS 5c–a. In the Chapoma 1 sequence, Unit 2 clayey silt and silt with scarce shell fragments can be interpreted as the sedimentary succession associated with the shallowing marine basin that existed during the final phase of MIS 5a and onset of MIS 4. Subdivided into subunits, Unit 4 is similar to Unit 4 in the Chavanga sequence (section b in Fig. 5). In the Chapoma 1 sequence, subunit 1 of Unit 4 was dated by three samples at 29.1 ± 3.4 ka (RLQG 2599-069), 23.0 ± 2.7 ka (RLQG 2600-069), and 76.7 ± 8.9 ka (RLQG 2601-069) (Fig. 7, Table 1). The last, in addition to reliable IR-OSL dating, is confirmed by the results of TL analysis performed on all three samples of subunit 1 of Unit 4. Figure 8 shows that the TL curves of feldspars from RLQG 2599-069 (29.1 ka) and RLQG 2600-069 (23.0 ka) exhibit a striking similarity, whereas the curve for RLQG 2601-069 (76.7 ka) displays noticeable differences. This observation suggests that the source areas for samples dated at 29.1 ka and 23.0 ka share a relatively consistent origin, unlike the third sample. It appears that the material corresponding to this sample, dated at 76.6 ka (MIS 5a), originated from a distinct source, most likely of marine origin, and was redeposited through various physical processes. In this case, it is probable that the redeposition occurred due to the action of glacier ice. Subunit 2 of Unit 4 (Chapoma 1 section in Fig. 7) is suggested to be a basal till. In previous studies, glaciofluvial coarse-grained sediments were identified in the uppermost part of Unit 4 in the Chapoma 1 sequence (Gudina and Yevzerov, Reference Gudina and Yevzerov1973; Korsakova, Reference Korsakova2021).
Thermoluminescence (TL) glow curve shapes for natural (left) and irradiated to 200 Gy (right) potassium feldspar grains extracted from sedimentary samples RLQG 2599-069 (29.1 ± 3.4 ka), RLQG 2600-069 (23.0 ± 2.7 ka), and RLQG 2601-069 (76.1 ± 5.9 ka) taken from the Chapoma 1 section at 20.5, 22.3, and 23.4 m asl, respectively.
The TL glow curve shapes of the natural and irradiated feldspar grains are practically indistinguishable for the two lower samples dated back to 29.1 ka and 23.0 ka, but the two lower samples diverge significantly from the upper sample dated back to 76.1 ka. This suggests a similarity in the chemical and physical environments during the crystallization of the feldspars in the two lower samples throughout their geologic history and a contrast with the conditions for the upper sample. Consequently, it is inferred that the sedimentary material in the upper layer originated from a different source area than the sediments in the lower layers.
Chapoma 2 sequence
In the Chapoma 2 outcrop (Figs. 2 and 7), lithologic and chronological data suggest the late Pleistocene Units 1a–1b–2–3–4 sediment succession. Clayey silt and silt of Unit 1a in the Chapoma 2 sequence are attributable to the initial phase of the Boreal transgression, as their lowermost lithostratigraphic position under marine sandy member (Unit 1b) was dated with the IR-OSL method to 95.3 ± 8.1 ka (RLQG 2661-070) in the Chapoma 2 south section and 72.4 ± 5.6 ka (RLQG 2604-069) in the Chapoma 2 north section (5N in Fig. 7), with the south and north sections being located about 40 m apart. In the north section at the Chapoma 2 site, a late phase of the Boreal transgression (Belomorian transgression) is marked by silt and sand IR-OSL dated to 72.4 ± 5.6 ka (Unit 1b in the Chapoma 2 north section, Fig. 7). Here, Unit 2 sand is correlated to MIS 4 according to its lithostratigraphic position in the sequence; sedimentologic features of this sand, which is locally laminated, included flattened pebbles and shell fragments and interbedded with gravel, demonstrate its accumulation in the shallow-water marine environments of the Mezen transgression (Zaretskaya et al., Reference Zaretskaya, Korsakova, Molodkov, Ruchkin, Baranov, Rybalko, Lugovoy and Merkuliev2022). The overlying wavy or oblique laminated sand of subunit 1 of Unit 3, IR-OSL dated to 36.2 ± 2.8 ka (RLQG 2603-069) and 32.5 ± 3.8 ka (RLQG 2602-069), and laminated silt of subunit 2 of Unit 3 suggest it was deposited during MIS 3 in a marine basin. The uppermost Unit 4 deposits are interpreted to be glaciomarine sediments accumulated at the glacier–sea contact during the late Valdaian (late Weichselian, MIS 2) glaciation.
Bolshaya Kumzhevaya sequence
In the Bolshaya Kumzhevaya outcrop, the Units 1b–2(?)–3–4 succession was identified based on the lithology and IR-OSL dating data (Fig. 7). This sedimentologic succession was investigated for the first time in 2003 (Korsakova et al., Reference Korsakova, Molodkov and Kolka2004); some differences were found only in the upper part of the excavated sequence in 2018.
The lowermost part of the Bolshaya Kumzhevaya sequence (Unit 1b in Fig. 7) consists of silty and sandy sediment that can be correlated with the second phase of the marine Boreal transgression. Judging by an IR-OSL date of 95.9 ± 6.3 ka (RLQG 2660-070) (Fig. 7, Table 1), only a minor sedimentary component of this marine event has been preserved in the section. Forming three lithologic units at 18.5–23.5 m asl and IR-OSL dated to 57.6, 55.5, and 44.4 ka (Fig. 7, Table 1), the sandy deposits accumulated in the shallow-water marine basin. The Unit 2(?) sediments suggest accumulation in an unstable littoral environment at the end of MIS 4 or at the onset of MIS 3. The Unit 3 deposits also accumulated in marine environments in early MIS 3; IR-OSL data demonstrate their MIS 3 age. No glacial or glacioaqueous sediments are obviously present in the Units 1a–3 succession here. Glacial clayey silt with sandy parting dated to 15.4 ka comprises Unit 4 and correlates with MIS 2.
Wider implications: the late Pleistocene palaeoenvironment on the Tersky Coast
The multilayered nature of the sequences provides evidence of changes in palaeoenvironmental conditions and associated depositional processes. Lithologic and chronological data on the sediments under consideration show that the sequences, as expected, consist mainly of marine and occasionally glaciomarine deposits (Fig. 9). Marine units are correlated not only to the interglacial epoch but also with periods of significant cooling (Fig. 10). Based on the available research material, the sedimentologic successions from the sections on the southern Kola Peninsula (Figs. 2, 4–7, 9, and 10) are indicative of the Boreal marine transgression, encompassing at least two phases of marine environments (Units 1a and 1b in the Varzuga sections 1(a2) and 1(b3) and in sections 2 [Chavanga], 4 [Chapoma 1], 5N and 5S [Chapoma 2 north and Chapoma 2 south], and 6 [Bolshaya Kumzhevaya] in Fig. 9). The interglacial marine environments also occurred in MIS 3 (Unit 3 in sections 2(b) [Chavanga], 3 [Kamenka], 5N [Chapoma 2 north], 6 [Bolshaya Kumzhevaya] in Fig. 9), covering in the studied sections a shorter time span than the entirety of MIS 3 (Fig. 10), in contrast to the data on the eastern regions of the Eurasian North, where MIS 3 marine sediments are ESR dated on mollusc shells to the interval of ca. 59–25 ka (Molodkov, Reference Molodkov2020). Furthermore, the presence of glacial expansion into the coastal areas of the White Sea during early MIS 4 (Unit 2 in the sections 1(b3) and 1(c) [Varzuga], 2(b) [Chavanga], 3 [Kamenka], and 4 [Chapoma 1] in Figs. 9 and 10) and the MIS 3 marine transgressions is observed. According to previous research in the areas adjacent to the Kola Peninsula (Svendsen et al., Reference Svendsen, Alexanderson, Astakhov, Demidov, Dowdeswell, Funder and Gataullin2004; Lambeck et al., Reference Lambeck, Purcell, Funder, Kjær, Larsen and Möller2006; Larsen et al., Reference Larsen, Kjær, Demidov, Funder, Grosfjeld, Houmark-Nielsen, Linge and Lyså2006), it is suggested that the BSIS or a local ice sheet may have reached the coastal Kola Peninsula in MIS 4 (see orange line in Fig. 1). Indirect evidence of the BSIS expansion is in the form of pebbles and boulders of exotic rocks that are not known in the bedrock on the Kola region. On the eastern Tersky Coast and adjacent Gorlo Strait shelf, boulders and pebbles of massive and organogenic-detrital carbonate rocks containing abundant fauna of brachiopods, bryozoans, and echinoderms were found in contemporaneous marine deposits (Korsakova et al., Reference Korsakova, Kolka, Zozulya, Lavrushin, Khoreva and Chistyakova2007; Korsakova, Reference Korsakova2021). Occurrence of short-lived marine environments and glacial expansion during the late Valdaian (late Weichselian, MIS 2) glaciation (Unit 4 in the studied sections in Fig. 9) were also demonstrated.
Late Pleistocene lithostratigraphic units and correlations of the studied sequences. The numbers 1–6 indicate the sections under consideration; their locations are shown in Fig. 1. Lithostratigraphic units (U) correspond to those in Figs. 4–7; the genesis of sediments is shown by letters: g, glacial; f, glaciofluvial; l, glaciolacustrine; gm, glaciomarine.
The Varzuga section provides a comprehensive late Pleistocene history from the MIS 5–4 stages. Age determinations reveal marine sediments with ESR and IR-OSL ages ranging from 104.0 ± 8.3 ka to 95.9 ± 11.1 ka, indicating the latter half of the MIS 5 marine environment. The unit in the lowermost part of the Chapoma 1 section indicates an initial phase of the Boreal transgression starting around 138–128 ka, while the following unit in Chapoma 2 section shows its late phase, dated at 72.4 ± 5.6 ka. The lower part of the Bolshaya Kumzhevaya section dates back to 95.9 ± 6.3 ka and is correlated with the middle of the Boreal transgression.
The presence of Unit 2 glacial deposits in the Kamenka sequence (section 3 in Fig. 9) and glacioaqueous deposits in the Chavanga and Chapoma 1 sequences (sections 2 and 4 in Fig. 9) indicates that during the MIS 4 period, the ice sheet covered only the eastern part of the Kola Peninsula. During this time, marine environments predominated in the White Sea depression, with marine and glacial–marine sediments identified and dated to 69‒32 ka in the studied sediment successions (sections: 1(c) [Varzuga] 2 [Chavanga], 4 [Chapoma 1], 5N [Chapoma 2 north], 6 [Bolshaya Kumzhevaya] in Fig. 9). Figure 10 illustrates that these sediments correlate not only with the warm MIS 3 but also with the cold MIS 4.
The presence of marine deposits with ages ranging from approximately 37 ka to 59 ka in the Chavanga, Kamenka, Chapoma 1, and Chapoma 2 sections suggests an accumulation stage in the marine environment, encompassing the very end of MIS 4 and almost the entirety of MIS 3. The research material from the studied sections provides evidence of the Boreal marine transgression during MIS 5, a short glacial expansion into coastal areas of the White Sea during early MIS 4, and glaciation in the late MIS 3 and MIS 2 (Fig. 10).
Conclusions
The study of the multilayered sedimentary sequence on the southern coast of the Kola Peninsula provided valuable insights into the dynamics of the palaeoenvironment and accumulation processes for most of the late Pleistocene period.
It was found that the sedimentary sequences are associated with the Boreal marine transgression during MIS 5, characterized by at least two phases of marine environments (Units 1a and 1b). Additionally, the expansion of glaciers into the coastal areas of the White Sea during MIS 4 (Unit 2) and the alternating transgression and regression phases of the sea during MIS 3 (Unit 3) was observed, along with the formation of short-lived marine environments and glacier expansion during the late Valdaian (late Weichselian, MIS 2) glaciation (Unit 4).
Marine environments are correlated not only with the interglacial stages, but also with glacial stages. During MIS 4, the BSIS or a local ice sheet reached the coastal Kola Peninsula and a marine basin continuously existed in the White Sea depression during MIS 5–3 and early MIS 2. The Tersky Coast of the White Sea was glaciated in the late MIS 3–MIS 2.
The results of this study emphasize the significance and complexity of palaeoenvironmental processes on the southern Kola Peninsula and provide new insights for understanding climate changes and glacial events in the late Pleistocene. Overall, our study underscores the importance of an integrated approach to studying sedimentology and palaeoenvironmental conditions for advancing our knowledge of past palaeoenvironmental dynamics. These findings contribute to our understanding of late Pleistocene geologic history and highlight the complex alternation between marine transgressions and glacial expansions and how it affected environmental conditions in the northern European region.
Data availability statement
Data presented in this manuscript are available on request.
Funding
This research was carried out within the framework of the research projects of the Geological Institute of the Kola Science Centre, Russian Academy of Sciences FMEZ-2024-0007 and Institute of Geography, Russian Academy of Sciences FMWS-2024-0003 and partly (dating samples) supported by RSF grant 22-17-00081.
Competing interests
The authors state they have no known competing financial interests or personal relationships that could have appeared to influence the work described in this paper.
