Weichselian Research Articles

The organization of subglacial drainage during the demise of the Finnish Lake District Ice Lobe

Abstract. Unknown basal characteristics limit our ability to simulate the subglacial hydrology of rapidly melting contemporary ice sheets. Sediment-based landforms generated beneath Late Pleistocene ice sheets, together with detailed digital elevation models, offer a valuable means of testing basal hydrology models, which describe the flow and dynamics of water in the subglacial system. However, to date no work has evaluated how well process-based subglacial hydrology models represent the hypothesized conditions associated with glaciofluvial landform formation in the palaeo setting. Previous work comparing model output to geomorphological evidence has typically done so using models that do not resolve subglacial processes and instead express likely subglacial water pathways. Here, we explore the ability of the Glacier Drainage System model (GlaDS), a process-based subglacial hydrology model, to represent the genesis conditions associated with a specific glaciofluvial landform termed “murtoos”. Distinctive triangular landforms found throughout Finland and Sweden, murtoos are hypothesized to form 40–60 km from the former Fennoscandian Ice Sheet margin within a “semi-distributed” system at the onset of channelized drainage in small cavities where water pressure is equal to or exceeds ice overburden pressure. Concentrating within a specific ice lobe of the former Fennoscandian Ice Sheet and using digital elevation models with a simulated former ice surface geometry, we forced GlaDS with transient surface melt and explored the sensitivity of our model outcomes to parameter decisions such as the system conductivity and bed topography. Our model outputs closely match the general spacing, direction, and complexity of eskers and mapped assemblages of features related to subglacial drainage in “meltwater routes”. Many of the predictions for murtoo formation are produced by the model, including the location of water pressure equal to ice overburden, the onset of channelized drainage, the transition in drainage modes, and importantly the seasonal sequence of drainage conditions inferred from murtoo sedimentology. These conclusions are largely robust to a range of parameter decisions, and we explore seasonal and inter-annual drainage behaviour associated with murtoo zones and meltwater pathways. Our results demonstrate that examining palaeo basal topography alongside subglacial hydrology model outputs holds promise for the mutually beneficial analyses of palaeo and contemporary ice sheets to assess the controls of hydrology on ice dynamics and subglacial landform evolution.

Evidence for an active ice margin during the last deglaciation: the Vimmerby Moraine, South Swedish Uplands

The Vimmerby Moraine is the only significant ice‐marginal moraine on the eastern side of southern Sweden, but no detailed studies exist on its formation during the final deglaciation of the Fennoscandian Ice Sheet. Through ground‐penetrating radar surveys and detailed sediment logging, we provide evidence for an active, oscillating ice margin during the formation of the Vimmerby Moraine, suggesting that the deglaciation of the South Swedish Uplands was, at least in some regions, dynamic. Ground‐penetrating radar surveys enabled imaging of internal sediment and delineation of the bedrock surface. These were complemented by common mid‐point surveys and sediment logging, as well as lithofacies analysis at three exposures in agravel pit. This approach revealed multiple subglacial till units partially separated by intercalated glacifluvial deposits. The glacifluvial sediments exhibit evidence of glaciotectonism, suggesting active overriding by the last ice sheet. Further evidence of an active ice margin is provided by the ground‐penetrating radar profiles collected perpendicular to the moraine crest. These contain a series of northerly dipping reflectors, which we interpret as evidence of repeated basal freeze‐on and melt‐out of sediment slabs during ice margin oscillations, as has been observed at contemporary glacier margins in Iceland, Norway, and the Alps. The data presented here demonstrate that the last Fennoscandian Ice Sheet remained active around the time of the Vimmerby Moraine formation. This work highlights the benefits of including detailed sediment logging and near‐surface geophysical surveys in the interpretation of deglaciation dynamics.

Land‐To‐Sea Mapping of the Glacial Erosion Unconformity Reveals Evolution of the Jasmund Glacitectonic Complex East of Rügen Island (SW Baltic Sea)

AbstractGlacial movements shaped vast northern parts, offering critical insights into glacial dynamics in a changing climate. Located on the island of Rügen in NE Germany, the Jasmund Glacitectonic Complex (JGC) is a key area to study the dynamics of past glaciations. Previous reconstructions focused primarily on the onshore realm, resulting in some areas remaining unexplored. Here we use more than 140 high‐resolution marine multi‐channel seismic profiles to map the erosional unconformity surrounding the JGC for the first time. Submarine glacial features match features observed onshore, allowing a consistent land‐to‐sea reconstruction of the evolution of the JGC. Our results indicate a single SW‐directed Weichselian glacier advance, suggesting that the JGC formed through three distinct glacier lobes exerting pressure from multiple directions. The ice advance encircled the Jasmund peninsula and overthrusted Cretaceous sediments on the JGC perpendicularly and laterally.

Imprint of a dead-ice environment on lidar imagery—Spatial patterns at the last Scandinavian Ice Sheet margin (N Poland)

There is a unique glacial landscape system in the southern Żarnowiec Moraine Plateau of northern Poland. Here the terrain is characterized by a very high concentration of ridges that are often adjacent to pits, which together constitutes a special type of hummocky topography. The morphological diversity of the study-area ridges and pits was determined from high-resolution lidar imagery that was qualitatively and quantitatively processed, then analysed. Fourteen test fields of spatial dimensions 0.5 km × 0.5 km were selected within the moraine plateau for detailed studies. The results show that the ridges and pits constitute a very dense network of approximately circular, diverse-oriented forms of various sizes, most of which are linked to each other in a step-wise fashion. The pits are surrounded by ridges of various heights and widths that separate one pit from another. A characteristic feature of the studied pits is the orientation of their long axis. Form-based morphological analysis revealed that their genesis was associated with a glacial rather than a periglacial environment. It is believed that ridge orientation reflects a spatial pattern related to ice crevasses. A moraine plateau constituted a terrain obstacle that acted as a counterslope to the advancing ice sheet. The studied landforms and landscape resulted from the ice sheet's breaking up into dead-ice blocks of various sizes. The very dense system of perpendicular, parallel, and diagonal crevasses was related to ice sheet advance from a northerly direction, resulting in a flow around the obstacle and differences in ice-mass speed.

Environmental changes in the Fleuve Manche paleoriver drainage system (Western Europe) linked to North Atlantic sub-millennial climate variability across Heinrich Stadial 1: Palynological evidence from the Bay of Biscay

Marine microfossils (dinoflagellate cysts and planktonic foraminifera) and geochemical (XRF-Ti/Ca)-based climatic records from a core (MD13–3438) located off the Fleuve Manche (FM) paleo-mouth have revealed that sustained warm summer sea surface temperatures (SSTs) during sub-millennial climate changes within HS1 (∼18–14.7 ka) may have played a key role in the FM regime related to the European Ice Sheet (EIS) melting rate. In this study, we have analyzed the MD13–3438 pollen content over the HS1 at a mean resolution of ∼50 years to test whether vegetation-based air temperatures were coupled to SSTs face to this rapid climate variability. First, our results highlight two major phases of pollen sources at site MD13–3438, preventing the pollen record to be interpreted as a continuous record of the evolution of vegetation and climate occupying a single watershed across HS1. The first phase, i.e. the HS1-a interval (∼18–16.8 ka), is marked by strong occurrences of boreal pollen taxa (especially Picea-Abies). Considering their spatial distribution and the coalescence of the British and Scandinavian ice sheets into the North Sea during the Last Glacial Maximum, these taxa probably originated from the North European Plain, i.e., eastern FM tributaries (east of the Rhine River), where cool-humid conditions generally prevailed. Then, the second phase, i.e. the HS1-b interval (∼16.8–14.7 ka BP), is characterized by a deceleration of the EIS retreat and the drop of boreal pollen values at site MD13–3438 further signing a less influence of the upstream FM drainage system and thus a better characterization of pollen sources related with western FM tributaries. Superimposed to these two HS1 main phases, pollen fluctuations are concomitant with sub-millennial variability in the EIS deglaciation intensity. During the early HS1 (HS1-a), we discuss two short-term increases in the ratio between deciduous trees (Quercus-Corylus-Alnus) and herbaceous plants (Plantago-Amaranthaceae-Artemisia). These events are coeval with phases of increasing dinocyst-based SST seasonality (i.e. through summer SST amplification). We associate these events with lower contribution of the upstream FM catchment as well as, possibly, atmospheric warming and regional sea-level positive oscillations. The HS1-b is composed of three main phases that appear more influenced by the downstream FM drainage system. HS1-b1 (16.8–16.3 ka BP) corresponds to the driest and coldest conditions west of the Rhine River. HS1-b2 (16.3–15.5 ka BP) is coeval with large arrivals of iceberg from the Hudson strait in the Bay of Biscay and thus likely to a major sea-level positive oscillation associated with a phase of FM valley reworking. HS1-b3 (15.5–14.7 ka BP) corresponds to persistent arid conditions that preceded the subsequent more humid conditions recorded from 14.7 ka BP at the start of the Bölling-Alleröd.

Hydraulic damage in subglacial conduits: Evidence from rock hydrofracture and hydraulic jacking for high fluid pressures during rapid melt of the Fennoscandian Ice Sheet

When large volumes of meltwater arrive suddenly at the glacier bed, the capacity of subglacial conduits is exceeded, generating transiently high fluid pressures. Yet we currently lack reliable markers for where such conditions developed on the beds of former glaciers. Here, we describe hydraulic damage that developed at close to or above overpressure (Pw > Pi) at the bed of the last Fennoscandian Ice Sheet in Uppland, east-central Sweden. Hydraulic damage is mostly confined within sharp lateral and vertical limits at the boundaries of former large subglacial conduits. The conduits include large Röthlisberger- and Nye-channels - eskers, esker corridors, and meltwater channels – that operated at full pipeflow. Wider (0.4–1.5 km) meltwater corridors opened under sheetflow with broad, shallow Hooke-channels. Hydraulic damage found in Precambrian gneiss and granite bedrock includes hydrofracturing and hydraulic jacking, with fracture propagation, dilation, and block displacement. Dilated fractures are often plugged by massive sands and diamicts or laminated silts. The fracture fills represent brief, single and multi-phase, injections of viscous, sediment-laden fluids at the ice sheet bed. The most spectacular effects of hydraulic damage are seen in the distension, disruption, and bursting of 1–20 m high rock bumps and hills. We present a detailed taxonomy for hydraulic damage and the consequent glacitectonic damage sustained during glacial ripping.Industrial fracking provides a well understood analogue for hydrofracture in nature. In subglacial settings, pathways for the flow of pressurised fluids are determined by the tensile strengths and hydraulic transmissivity of ice, till, and rock. The easiest flow path is across the glacier bed beneath sliding ice. Abrupt increases in meltwater flow force ice-bed separation. Opening of conduits under high fluid pressures triggers further propagation of hydrofractures into ice, till, and rock. On the floors of subglacial conduits, pressurised fluids entered the hard rock bed via old, mineral coated fractures. Under Poiseuille's law, fracture flow was initially slow due to the narrow apertures, but flow accelerated rapidly as fluid wedging dilated fractures. Fluid pressures in old fractures are constrained by published measurements for near-surface fracture toughness of ≤0.7 MPa at Forsmark, east-central Sweden. Formation of new rock hydrofractures overcame gneiss tensile strengths of 8–18 MPa. Hydraulic jacking required accommodation space that developed during conduit opening. Due to buoyancy and ice bed separation, lower excess pressures of <0.3 MPa were required to lift rock sheets up to 10 m thick. Hydraulic bursting developed in large conduits under very high, but rapidly fluctuating fluid pressures as dilating rock fractures extended to the glacier bed. Distension of rock surfaces was brief due to rapid leak-off and the inherent instability of overpressured conditions in dilatant, fractured rock masses.Hydraulic damage developed in Uppland when accelerating ablation generated large subglacial meltwater fluxes in response to abrupt warming at the Younger Dryas-early Holocene transition (∼10.8 ka). The large volumes of subglacial meltwater, the sudden increases in subglacial water flow and fluid pressures, and the release and transport of 2 m b-axis boulders in esker tunnels suggest that hydraulic damage was triggered at the onset of large subglacial floods. Ice thicknesses at distances of up to 50 km behind the grounding line at the retreating subaqueous ice margin indicate that steady-state hydrostatic pressures were ∼2.9–4.5 MPa. Such pressures are consistent with (i) extensive, but brief ice-bed separation in meltwater corridors under overpressure and (ii) the development of hydraulic damage in conduits close to the ice margin. More extreme fluid pressures developed instantaneously under hydraulic shock, as rapid subglacial pipeflow was suddenly stopped or diverted. The taxonomy of hydraulic damage presented here can be applied to any hard glacier bed to identify where large volumes of subglacial meltwater reached high pressures. Similar markers are widely reported from modern and former glacier beds.

Weichselian–Holocene glacial history of the Sjuøyane archipelago, northern Svalbard

The Sjuøyane archipelago is the northernmost land area of Svalbard; thus, it provides a window to study the terrestrial glacial history and dynamics of the Svalbard–Barents Sea Ice Sheet and complement marine geological studies in the region. To reconstruct the glacial history of Sjuøyane, we describe coastal sedimentary sections in Quaternary sediments and constrain their chronology by radiocarbon and optically stimulated luminescence ages. Erratic boulders and bedrock are sampled for 10Be cosmogenic exposure dating, aiming to determine the deglaciation age and exposure history. Holocene environments are studied based on lake sediments and emerging vegetation from retreating snow patches. The sedimentary sections largely consist of shallow (glacio‐)marine and/or littoral sediments deposited during high relative sea levels. The radiocarbon and luminescence ages suggest they formed during a Middle Weichselian interstadial and after the Late Weichselian glaciation. A wave‐washed bedrock erosional notch and rounded boulders at 36±1 m a.h.t. most likely formed during this interstadial. Most of the cosmogenic 10Be ages are older than the last deglaciation, likely indicating a complex exposure history. One boulder sample suggests that the lowlands were deglaciated 14.7±1.82 ka ago, and two boulder samples with ages of 18.94±3.26 and 22.89±4.05 ka suggest that the highlands were possibly ice‐free at this time. The lake sediments from Isvatnet, Phippsøya, consist of glaciolacustrine silt and clay overlain by gyttja. The gyttja has accumulated at least since 7.0 cal. ka BP. Two radiocarbon ages from emerging vegetation suggest Neoglacial cooling since 3.8 cal. ka BP. A patchy glacial drift at the surface of Sjuøyane and well‐preserved pre‐Late Weichselian sediments suggest that the Late Weichselian glaciation was non‐erosive and/or cold‐based at this part of the north margin of the Svalbard–Barents Sea Ice Sheet.

Sub- and supra-moraine loesses in the Vistula catchment (Poland) with regard to the age and extentof the Scandinavian ice-sheets

The study focuses on the geochronology and correlation between loesses and glacial tills in the Vistula catchment in Poland, in an area covered by all Scandinavian ice-sheets and characterized by thick loess successions. The basis for this correlation was a litostratigraphic analysis of selected loess sites in the uplands, which are well-documented in the literature, and the occurrences of loess and loess-like deposits under the cover of glacial deposits in the lowland part of the Vistula catchment, characterized by a much sparser documentation. The nomenclature of loess horizons follows Maruszczak (2001) and Kukla (1987), but sometimes is significantly changed, according to modern Quaternary stratigraphical schemes. According to the analysis, accumulation of almost all loess horizons distinguished so far has been documented in the Pleistocene succession of the Vistula catchment: the oldest – lower (LNd), middle (LNs), and upper (LNg); older: lower (LSd), middle (LSs), and upper (LSg); and younger: lowest (LMn), lower (LMd), middle (LMs), and upper (LMg). In most cases loess accumulation took place in steppe-tundra conditions preceding the maximum ice-sheet development during the succeeding Scandinavian glaciations. For selected sites, the loess occurrence is presented in superposition to glacial tills and interglacial deposits. The distinguished loesses and glacial tills are tied to the stratigraphic schemes of the Pleistocene in Poland (Marks et al. 2016, 2019) and correlated with Ukrainian loess horizons (Łanczont et al. 2019).

Reconstructing the Holocene glacial history of northern Troms and western Finnmark, Arctic Norway

Here we present the first Lateglacial and Holocene glacial history from Rotsunddalen, northern Troms and western Finnmark county, northern Norway, based on both relative and numerical moraine dating using Schmidt hammer, soil chronosequencing and terrestrial cosmogenic nuclide dating. We combine these chronological data with a regional map of the glacial geomorphology to hypothesize key events in the glacial history from around 14 ka to present. Our reconstruction shows that, following deglaciation of the main ice sheet across central Troms and Finnmark, mountain glaciers were terminating on land, close to the coast, between around 12.1 and 10.6 ka. Continued recession of the main Fennoscandian Ice Sheet margin towards the SE led to the isolation of several large plateau icefields and outlet glaciers that generated moraines at around 10.2–9.2 ka, which we ascribe to the Erdalen Event, and 8.4–8.2 ka, which is broadly contemporaneous with the 8.2 ka cold event. Although the latter corresponds with the Scandinavian Finse Event, very few moraines have been dated to this time and we therefore view it as a tentative hypothesis for future work to test. During the Holocene Thermal Maximum (~6.6 to 6.3 ka) most (if not all) glaciers in the region disappeared, but then regrew during the Neoglaciation and produced large moraines dated to around 4.7 ka that lie a few hundred metres distal to the prominent Little Ice Age moraines (previously dated to AD 1810s–1870s). Given the limitations of our dating approach, the preservation of moraines dated to this period in northern Norway also warrants further investigation. We also highlight that terrestrial cosmogenic nuclide dating of the moraines is not consistent with other dating approaches and the widely established deglaciation history of the region, probably owing to cosmogenic inheritance and insufficient glacial erosion.

Ice-marginal ridge relief complex in northern Kola Peninsula (NW Russia): morphology, structure, and genetic interpretation

A widely spread ridge relief complex was formed by the last ice cover in the northeastern part of the Kola Peninsula. The complex belongs to the ice-marginal formations of the ice streams of the north-eastern Fennoscandian Ice Sheet. The applied morphometric research technique (including the analysis of the digital elevation model) has allowed identifying the ridge relief complex to make a sub-longitudinal belt. The belt comprises frontal and radial ridges, ridges and hills of the distal part of the belt, and complex-shaped ridges. These landforms were studied in eight outcrops and trenches. The frontal and some complex-shaped ridges consist of basal and ablation tills and glaciofluvial deposits deformed by the east- and north-eastward advancing glacier. The radial ridges composed of glaciofluvial sediments are the eskers. Glaciofluvial deposits with a thin cover of ablation and flow tills form the ridges on the distal part of the belt. They mark the limits of glacier expansion at an individual glacial stage. Some complex-shaped ridges are composed only of the ablation melt-out till and formed in crevasses of dead ice massifs. The ridge relief complex shows morphological similarities with the Veiki moraine in North Sweden, Pulju moraine in Finland, and ridges of the “ice-walled-lake plains” in North America. The location of the moraine ridge landforms is controlled by the topography of crystalline bedrock. During the last glaciation, the glacier retreated in several stages marked by belts of the ridge landforms. Each stage was followed by a series of oscillations corresponding to chains of frontal moraine ridges. The correlation of the ice-marginal forms in other parts of the Kola and Karelia regions allows the authors to refer them to the Neva Stage (other names are Keiva I, Syamozero), which took place in the Older Dryas Stadial.

Late Middle Pleistocene (MIS 10–6) glacial–interglacial records from loess–palaeosol and fluvial sequences from northern France: a cyclostratigraphic approach

This contribution summarizes the most informative loess–palaeosol sequences (LPS) and fluvial terrace records from the late Middle Pleistocene (LMP) of northern France demonstrating the reliability of the cyclostratigraphic approach for the interpretation of pedosedimentary sequences controlled by major glacial–interglacial climatic cycles. In this area, continental mollusc assemblages from interglacial fluvial silts and calcareous tufas are particularly rich and diverse and marker species define the malacological signatures of each interglacial optimum for Marine Isotope Stages (MIS) 11, 9 and 5e. This approach shows that the forest was less developed during MIS 7, suggesting that climatic conditions were either drier and/or cooler than during other Pleistocene interglacials. In the Somme basin, the terrace system shows that two alluvial formations were deposited between early MIS 8 and late MIS 7. In some LPS, the occurrence of two pedosedimentary sub‐cycles (IIa and IIb), separated by a relatively long (~12 ka) and cold period, corresponding to an ‘aborted glacial’ (MIS 7d), underlines the complex pattern of this unusual ‘interglacial’. Overall, during the LMP, each soil complex corresponding to interglacial and early‐glacial periods from MIS 11 to MIS 5 is broadly made up of the same soil facies but exhibits a specific succession pattern or signature. Throughout the area, LPS show a huge change in both the deposition rates and the geographical extent of typical calcareous loess at the beginning of MIS 6. This so‐called ‘Loess Revolution’ probably reflects a change in the palaeogeography of the southern North Sea and eastern Channel source areas at times marked by the coalescence of the British and Scandinavian ice sheets. MIS 6 is also characterized by the oldest evidence of permafrost development in the area. In addition, this work allows the age of the Lower/Middle Palaeolithic boundary to be confirmed, with the oldest occurrence of Levallois technology being around 300 ka.

Isotopic constraints on the late Pleistocene glacial water and sediment inputs to the central Arctic ocean

The Quaternary history of the Arctic Ocean reflects pronounced periodic changes in marine and glacial environments. This variability is an important component of the Arctic paleoclimatic system. For more insight, we use isotopic measurements on the detrital (Nd-Sr-Pb) and authigenic (Nd) phases of the Quaternary sediments from the central Arctic Ocean. We analyze core ICE4 from ∼3 km water depth off the Lomonosov Ridge. This core was previously studied for other sedimentary proxies (Dong et al., 2022). The analyzed record is estimated to cover major Late Pleistocene Weichselian/Wisconsinian glaciations occurring over the last ∼100–200 ka. Our results provide new evidence that glacially-derived waters from the glaciated margins reached the bottom of the deep Arctic basins. We infer that these waters were delivered to the deep basins by hyperpycnal flows originating from subglacial outbursts accompanying overflow meltwater and iceberg discharge during glacial/deglacial events. Voluminous outbursts during the older, presumed Early/Middle Weichselian glaciation, ca. 115-35 ka, are traced to the West Siberian margin impacted by the Eurasian Ice Sheet. Similar sedimentary layers are found in longer records from the central Arctic Ocean at older stratigraphic intervals. Sediments constrained to the Last Glacial Maximum have less pronounced glacigenic inputs with a stronger Laurentide Ice Sheet signature. The circulation patterns were likely dependent on the interplay of the timing and extent of the Laurentide and Eurasian ice sheets.

Reconstructing dynamics of the Baltic Ice Stream Complex during deglaciation of the Last Scandinavian Ice Sheet

Abstract. Landforms left behind by the last Scandinavian Ice Sheet (SIS) offer an opportunity to investigate controls governing ice sheet dynamics. Terrestrial sectors of the ice sheet have received considerable attention from landform and stratigraphic investigations. In contrast, despite its geographical importance, the Baltic Sea remains poorly constrained due to limitations in bathymetric data. Both ice-sheet-scale investigations and regional studies at the southern periphery of the SIS have considered the Baltic depression to be a preferential route for ice flux towards the southern ice margin throughout the last glaciation. During the deglaciation the Baltic depression hosted the extensive Baltic Ice Lake, which likely exerted a considerable control on ice dynamics. Here we investigate the Baltic depression using newly available bathymetric data and peripheral topographic data. These data reveal an extensive landform suite stretching from Denmark in the west to Estonia in the east and from the southern European coast to the Åland Sea, comprising an area of 0.3 million km2. We use these landforms to reconstruct aspects of the ice dynamic history of the Baltic sector of the ice sheet. Landform evidence indicates a complex retreat pattern that changes from lobate ice margins with splaying lineations to parallel mega-scale glacial lineations (MSGLs) in the deeper depressions of the Baltic Basin. Ice margin still-stands on underlying geological structures indicate the likely importance of pinning points during deglaciation, resulting in a stepped retreat signal. Over the span of the study area we identify broad changes in the ice flow direction, ranging from SE–NW to N–S and then to NW–SE. MSGLs reveal distinct corridors of fast ice flow (ice streams) with widths of 30 km and up to 95 km in places, rather than the often-interpreted Baltic-wide (300 km) accelerated ice flow zone. These smaller ice streams are interpreted as having operated close behind the ice margin during late stages of deglaciation. Where previous ice-sheet-scale investigations inferred a single ice source, our mapping identifies flow and ice margin geometries from both Swedish and northern Bothnian sources. We anticipate that our landform mapping and interpretations may be used as a framework for more detailed empirical studies by identifying targets to acquire high-resolution bathymetry and sediment cores and also for comparison with numerical ice sheet modelling.

Stadial and interstadial deposits of Late Nemunas (Late Weichselian/MIS 2) glaciation in south Lithuania and their interpretation

This paper presents sedimentological and paleobotanical studies of the deposits from Last Glaciation (Late Nemunas, Late Weichselian) in the middle reaches of the River Nemunas in southern Lithuania. Spores, pollen, and other palynomorphs in the intercalated sediments between the diamicton (glacial till) beds were analyzed, and the paleogeographical and paleoclimatic conditions of their sedimentation were assessed using the modern analog technique (MAT). The two till layers from the Last Glaciation, characterized by significant differences in their composition, testify to a difference in glacier movement paths in the region. The continuous sedimentation process between glacial tills provides evidence of global climate changes during the Last Glaciation and the restructuring of glacial lobes during deglaciation in the southeastern part of the Scandinavian Ice Sheet. The sedimentation process and vegetation composition bear witness to an open landscape of steppe-tundra, scattered wooded areas, and the emergence of thermophilic vegetation afterward, suggesting increased annual temperatures and higher humidity. The available data enable the identification of the Balbieriškis interstage period within the Late Nemunas Glaciation (Upper Nemunas Formation) between the Grūda and Baltija stages (Grūda and Baltija subformations) in the regional Quaternary stratigraphic scheme. This period correlates with the Marine Isotope Stage 2 (MIS 2).

Physical modeling of ice-sheet-induced salt movements using the example of northern Germany

Abstract. Salt structures and their surroundings can play an important role in the energy transition related to a number of storage and energy applications. Thus, it is important to assess the current and future stability of salt bodies in their specific geological settings. We investigate the influence of ice sheet loading and unloading on subsurface salt structures using physical models based on the geological setting of northern Germany, which was repeatedly glaciated by the Scandinavian Ice Sheet during the Pleistocene. Apparent spatial correlations between subsurface salt structures in northern Germany and Weichselian ice marginal positions have been observed before, and the topic is a matter of ongoing debate. Recently described geomorphological features – termed surface cracks – have been interpreted as a direct result of ice-sheet-induced salt movement resulting in surface expansion. The spatial clustering and orientation of these surface cracks has not been well understood so far, owing to only a limited number of available studies dealing with the related salt tectonic processes. Thus, we use four increasingly complex physical models to test the basic loading and unloading principle, to analyze flow patterns within the salt source layer and within salt structures, and to examine the influence of the shape and orientation of the salt structures with respect to a lobate ice margin in a three-dimensional laboratory environment. Three salt structures of the northern German basin were selected as examples that were replicated in the laboratory. Salt structures were initially grown by differential loading and buried before loading. The ice load was simulated by a weight that was temporarily placed on a portion of the surface of the models. The replicated salt structures were either completely covered by the load, partly covered by the load, or situated outside the load extent. In all scenarios, a dynamic response of the system to the load could be observed; while the load was applied, the structures outside the load margin started to rise, with a decreasing tendency with distance from the load margin, and, at the same time, the structures under the load subsided. After the load was removed, a flow reversal set in, and previously loaded structures started to rise, whereas the structures outside the former load margin began to subside. The vertical displacements during the unloading stage were not as strong as during the load stage, and thus the system did not return to its pre-glaciation status. Modeled salt domes that were located at distance from the load margin showed a comparably weak reaction. A more extreme response was shown by modeled salt pillows whose margins varied from sub-parallel to sub-perpendicular to the load margin and were partly covered by the load. Under these conditions, the structures showed a strong reaction in terms of strain and vertical displacement. The observed strain patterns at the surface were influenced by the shape of the load margin and the shape of the salt structure at depth, resulting in complex deformation patterns. These physical modeling results provide more evidence for a possible interplay between ice sheets and subsurface salt structures, highlighting the significance of three-dimensional effects in dynamic geological settings. Our results lead to a better understanding of spatial patterns of the surface cracks that were mapped at the surface above salt structures and offer further room for interpretation of the influence of salt movements on the present-day landscape.

Stratigraphy and OSL chronology of the Middle–Upper Pleistocene sedimentary sequence and vegetation history during Late MIS6–MIS5e in the Neva Lowland (St. Petersburg region, Russia)

ABSTRACTThe quarry of the Sverdlov (Etalon) Factory stands out as a unique site in the northwest of the East European Plain since it exposes an almost complete Upper Pleistocene stratigraphic succession. Previous investigations have revealed a sequence of marine sediments deposited during the whole Mikulino Interglacial [Eemian, Marine Isotope Stage (MIS)5e], but the chronology of the Late Pleistocene deposits remains incomplete and requires further update. We present a study of the Middle–Upper Pleistocene units exposed in the Sverdlov Factory quarry by using geological methods, optically stimulated luminescence (OSL) dating and pollen analysis. The lower part of the Sverdlov Factory section comprises glaciomarine/glaciolacustrine varves, which accumulated during degradation of the Late Moscow (Late Saalian, Late MIS6) glaciation and correspond to pollen zone M1 (lower peak of Picea). They are covered by interglacial marine clayey silt that includes all Mikulino regional pollen zones M2–M8. According to a Bayesian model based on K‐rich feldspar luminescence ages, the interglacial marine deposition started 133 ± 8 ka ago and ended 109 ± 7 ka ago. The marine sediments are overlain by lacustrine silt with a K‐rich feldspar age of 119 ± 7 ka, lacustrine/alluvial sand with plant detritus formed 47–41 ka ago, and lacustrine sand and silt without organic matter accumulated ~39 ka ago. The upper part of the section is composed of subglacial till formed during the Late Valdai (Late Weichselian, MIS2) under the Scandinavian Ice Sheet and glaciolacustrine varved deposits of the Baltic Ice Lake. Two ages were acquired using quartz OSL from these varves, 16.9 ± 1.6 and 15.7 ± 1.0 ka, which are assumed to slightly overestimate the true age.

Millennial-scale fluctuations of palaeo-ice margin at the southern fringe of the last Fennoscandian Ice Sheet

Abstract. The paper presents the first terrestrial record of millennial-scale palaeo-ice margin oscillations at the southern fringe of the last Fennoscandian Ice Sheet (FIS) during the last glacial cycle. The study area is located in northern Poland close to the last FIS maximum limit. The chronology and dynamics of palaeo-ice margin oscillations at the southern fringe of the FIS are based on combined luminescence and 10Be surface exposure dating. Optically stimulated luminescence (OSL) was used to date sandy deposits (fluvioglacial sediments and aeolian deposits filling fossil periglacial wedges) intercalating basal till layers. The most likely age of the tills was constrained by Bayesian modelling of the sequence of OSL ages and lithostratigraphy. 10Be surface exposure dating was used on erratic boulders left during the final retreat of the last FIS and resting on the surface of glacial landforms. Our results, which are mainly based on OSL chronology and Bayesian modelling, indicate millennial-scale oscillations of the last FIS in northern Poland between ∼19 and ∼17 ka. The last FIS retreated and re-advanced over a relatively short period of time (2–3 ka), leaving lithostratigraphic records (basal tills) of three ice re-advances over a millennial-scale cycle: 19.2±1.1, 17.8±0.5 and 16.9±0.5 ka. Despite 10Be surface exposure ages obtained for 14 erratic boulders being poorly clustered, the main mode of age distribution occurs at ∼18 ka and indicates a possible signal of the ice sheet retreat after one of the re-advances. We explore the dynamics of these oscillations and compare the proposed cycles of the southern FIS advances and retreats with existing patterns of the last deglaciation and millennial-scale fluctuations of the last FIS inferred from marine records.

The influence of glacial landscape evolution on Scandinavian ice-sheet dynamics and dimensions

Abstract. The Scandinavian topography and bathymetry have been shaped by ice through numerous glacial cycles in the Quaternary. In this study, we investigate how the changing morphology has influenced the Scandinavian ice sheet (SIS) in return. We use a higher-order ice-sheet model to simulate the SIS through a glacial period on three different topographies, representing different stages of glacial landscape evolution in the Quaternary. By forcing the three experiments with the same climate conditions, we isolate the effects of a changing landscape morphology on the evolution and dynamics of the ice sheet. We find that early Quaternary glaciations in Scandinavia were limited in extent and volume by the pre-glacial bathymetry until glacial deposits filled depressions in the North Sea and built out the Norwegian shelf. From middle–late Quaternary (∼0.5 Ma) the bathymetry was sufficiently filled to allow for a faster southward expansion of the ice sheet causing a relative increase in ice-sheet volume and extent. Furthermore, we show that the formation of The Norwegian Channel during recent glacial periods restricted southward ice-sheet expansion, only allowing for the ice sheet to advance into the southern North Sea close to glacial maxima. Finally, our experiments indicate that different stretches of The Norwegian Channel may have formed in distinct stages during glacial periods since ∼0.5 Ma. These results highlight the importance of accounting for changes in landscape morphology through time when inferring ice-sheet history from ice-volume proxies and when interpreting climate variability from past ice-sheet extents.

Refining Holocene sea‐level variations for the Lofoten and Vesterålen archipelagos, northern Norway: implications for prehistoric human–environment interactions

ABSTRACTThe Lofoten and Vesterålen archipelagos are located off the outer coast of northern Norway far from the center of the former Fennoscandian Ice Sheet and near the continental shelf edge. Existing relative sea‐level (RSL) data indicate a pronounced mid‐Holocene transgression and interesting connections with the region's prehistoric human settlement history. Here we present seven new sea‐level index points from isolation basins and five terrestrial limiting points from a coastal sedimentary sequence to refine the region's RSL history. Ingression and isolation contacts in isolation basin sediment cores are identified using sedimentary geochemical data, scanning X‐ray fluorescence profiles and phytoplankton analysis. The ages of these contacts are determined using radiocarbon‐based age models. Our index points range from 11.2 to 1.5k cal a bp and are combined with previously published data to predict the spatiotemporal evolution of sea level in this region using an ensemble of spatiotemporal empirical hierarchical models (STEHME). The new RSL curve constrains the timing of the mid‐Holocene transgression, which occurred from c. 9 to 6k cal a bp when sea level increased from −4 to 7 m above present day. From c. 6 to 5k cal a bp, RSL rapidly fell to c. 4 m above present values, and more gradually declined at an average rate of c. 0.8 m ka−1 over the last 5k cal a bp. Isobase maps derived using the STEHME show a decrease in the regional shoreline gradients since the transgression maximum from 0.25 to 0.07 m km−1. Our data also better define how RSL variations influenced the location and preservation of coastal settlement locations and harbors from the early Stone Age through historic intervals, improving understanding of regional human–environment interactions.

Younger Dryas and Early Holocene ice‐margin dynamics in northwest Russia

Younger Dryas and Early Holocene ice‐margin dynamics in northwest Russia