Weichselian Research Articles

Postglacial relative sea level change in Norway

We present the first comprehensive postglacial relative sea-level (RSL) database for the Norwegian coast from Oslo to the Kola Peninsula. The database spans the last 20 kyrs and is composed of 413 index points and 610 limiting data points derived from raised beaches, archeological data, glaciomarine terraces, and sedimentary indicators from isolation basins, salt marshes, and peat bogs. The data are quality controlled, assigned standardized indicative meanings, and recalibrated to current standards. We use an ensemble of Bayesian statistical models, trained on sea-level index points and weighted by their fit to both index and sea-level limiting points, to assess the spatiotemporal patterns of Norwegian RSL change. Continuous RSL fall driven by isostatic rebound in response to Eurasian ice sheet collapse dominates the RSL signal at every inland location in Norway. A first transgression (episode of RSL rise), which occurred in southwest Norway during the Younger Dryas (14–11.7 ka), increased RSL by as much as 15 m in some locations. A second transgression, named the Tapes transgression, occurred during the early-mid Holocene between 10 and 5.5 ka. The spatiotemporal model ensemble constrains the timing, amplitude, and spatial distribution of the Younger Dryas and Tapes transgressions. Based on our modeling results, we speculate that the Tapes Transgression was the result of global mean sea level rise temporarily outpacing isostasy-driven RSL fall, while the Younger Dryas Transgression was likely the result of local ice sheet readvance combined with low viscosity asthenosphere and weak lithosphere in the region. We also describe the effects of peripheral bulge migration on Norwegian RSL, which caused increased RSL in the early Holocene and a delayed Tapes transgression. We show that postglacial RSL data in Norway contain complex spatiotemporal patterns of nearfield RSL change that can best be estimated by combining a high-quality data compilation with glacial isostatic adjustment modeling via a robust statistical model.

Norway spruce postglacial recolonization of Fennoscandia

Contrasting theories exist regarding how Norway spruce (Picea abies) recolonized Fennoscandia after the last glaciation and both early Holocene establishments from western microrefugia and late Holocene colonization from the east have been postulated. Here, we show that Norway spruce was present in southern Fennoscandia as early as 14.7 ± 0.1 cal. kyr BP and that the millennia-old clonal spruce trees present today in central Sweden likely arrived with an early Holocene migration from the east. Our findings are based on ancient sedimentary DNA from multiple European sites (N = 15) combined with nuclear and mitochondrial DNA analysis of ancient clonal (N = 135) and contemporary spruce forest trees (N = 129) from central Sweden. Our other findings imply that Norway spruce was present shortly after deglaciation at the margins of the Scandinavian Ice Sheet, and support previously disputed finds of pollen in southern Sweden claiming spruce establishment during the Lateglacial.

The Elbe Palaeovalley: Evolution from an ice-marginal valley to a sedimentary trap (SE North Sea)

The southeastern North Sea was a sub-aerially exposed periglacial environment during the late Weichselian. At that time, the approximately 210 km long and 40 km wide SE–NW trending Elbe Palaeovalley (EPV) was the main conduit for meltwater from the southern margin of the Scandinavian ice sheet and rivers draining the North European plain. Only few studies investigated the environmental organization of the Northwest European shelves during deglacial sea level rise resulting in vague reconstructions of the drainage systems. Here, a detailed analysis of sediment echo-sounder data and core samples are presented to study the EPV morphology and stratigraphic infill. The valley base is about 65 m below the present sea-level. This SE-NW palaeo-drainage with its braided system developed during the MIS 2 sea level fall and lowstand in a periglacial environment. During deglacial sea-level rise, the EPV evolved into an estuary with tributaries, intertidal and subtidal flats. After the early Holocene (∼10 cal ka BP), ongoing sea-level rise drowned the EPV system which became the major sedimentary trap in the southern North Sea. Successive sea level rise and associated changes in the hydrodynamic regime led to three stages of shelf organization that were recorded in the valley infill successions: (1) shallow basin with restricted connections to Atlantic water; (2) the southwestern connection with the Dover Strait; (3) the open marine environment at 5.8 cal ka BP after the drowning of Dogger Bank.

Late Plio-Pleistocene evolution of the Eurasian Ice Sheets inferred from sediment input along the northeastern Atlantic continental margin

High-latitude marine sediment archives may contain information about the configuration and dynamics of former ice sheets, paleoclimate and the intensity of glacial erosion and uplift in catchment areas. Compiling information on the Late Cenozoic sediment packages along the NE Atlantic continental margin shows that large sediment volumes (c. 982 × 103 km3) were deposited during the Northern Hemisphere Glaciations (NHG), and that significant spatial and temporal variations in sediment input occurred during three distinct NHG Phases. NHG Phase I (2.7–1.5 Ma) was characterized by high sediment input to the trough mouth fan systems offshore Svalbard, suggesting strong glacial erosion and the development of large ice sheets over Svalbard. Comparatively moderate sedimentation and erosion rates are observed along the SW Barents Sea and the Norwegian margins during NHG Phase I. This indicates more restricted ice sheets over Fennoscandia compared to Svalbard, although periods of shelf edge glaciation most likely occurred. The most prominent overall margin development occurred during NHG Phase II (1.5–0.8 Ma), when recurrent large-scale, continental shelf edge, glaciations are suggested for the entire Eurasian Ice Sheets (EurIS). Compared to NHG Phase I, average sedimentation rates are three (91 cm/kyr) and two times higher (20 cm/kyr), respectively, in the Kara-Barents Sea-Svalbard Ice Sheet (KBSIS) region and the Fennoscandian/British-Irish Ice Sheet (FIS/BIIS) region. During NHG Phase III (0.8–0 Ma), sediment input decreases considerably (73%) along the marine margin of the KBSIS, while increasing significantly (62%) from the FIS/BIIS, in comparison to NHG Phase II. These estimates mark a major transition in the evolution of the EurIS, where the submergence below sea level of the Barents Sea region and the initiation of the Norwegian Channel Ice Stream are suggested to be key factors in this change. Furthermore, the submergence of the Barents Sea region affected ocean-atmosphere coupling and circulation which may have had a potential impact on global long-term climate change.

Influence of 3D Earth Structure on Glacial Isostatic Adjustment in the Russian Arctic

AbstractAnalyses of glacial isostatic adjustment (GIA) and deglacial relative sea‐level (RSL) change in the Russian Arctic deliver important insights into the Earth's viscosity structure and the deglaciation history of the Eurasian ice sheet complex. Here, we validate the 1D GIA models ICE‐6G_C (VM5a) and ICE‐7G_NA (VM7) and select new 3D GIA models in the Russian Arctic against a quality‐controlled deglacial RSL database of >500 sea‐level data points from 24 regions. Both 1D models correspond to the RSL data along the southern coast of the Barents Sea and Franz Josef Land from ∼11 ka BP to present but show notable misfits (>50 m at 10 ka BP) with the White Sea data. We find 3D model predictions of deglacial RSL resolve most of the misfits with the observed data for the White Sea while retaining comparable fits in other regions of the Russian Arctic. Our results further reveal: (a) RSL in the western Russian Arctic is sensitive to elastic lithosphere with lateral thickness variation and 3D viscosity structure in the upper mantle; and (b) RSL in the whole Russian Arctic is less sensitive to 3D viscosity structure in the lower mantle compared to the upper mantle. The 3D models reveal a compromise in the upper mantle between the background viscosity and scaling factor to best fit the RSL data, which needs to be considered in future 3D GIA studies.

Distribution and characteristics of Pleistocene ground thermal contraction polygons in Europe from satellite images

AbstractThe study of Pleistocene polygons in Europe carried out using the collection of satellite images available in Google Earth™ provides new data on the distribution of ground cracking by thermal contraction during glacial periods and sheds light on some factors controlling their formation. The distribution map shows that thermal contraction cracking affected terrain between latitude 43.5°N and the southern limit of the Fennoscandian Ice Sheet at 15 ka, with a concentration of polygons north of latitude 51°N. A clear asymmetry exists in relation to longitude, with a greater southward extension of polygons in France (43.5°N) than in central Europe (47°N). Analysis of the characteristics of polygons at the European scale reveals an association with lithology, latitude, and the age of the sediments in which they formed. The morphological evolution over time reconstructed from observations of polygons of contrasting ages indicates that a stable mature phase, characterized by small (mean size 15 m), regular polygons with numerous Y‐junctions, is reached after ca. 4 ka of thermal contraction cracking activity.

Soil-development differentiation across a glacial–interglacial cycle, Saalian upland, E Poland

Three soil profiles at the Koczery site (E Poland), studied through soil and sediment evaluation, reveal a polygenetic development of soils in areas unglaciated during the Weichselian glaciation (MIS 2) period. Pedogenesis had developed on glacial till (till) in period from the Saalian glaciation (MIS 6) to the Holocene (MIS 1). Wide-spectrum analyses demonstrate that the profiles represent a Planosol. Four stages (I–IV) of soil formation were distinguished in the study area. Stage I, comprising of a nutrient- and element-rich material, correlates to the accumulation of Saalian glacial till (MIS 6). Stage II, that took place during the Eemian Interglacial (MIS 5e), corresponds to the initial development of Cambisol as a result of carbonate leaching process. This was followed by the development of a Luvisol through the processes of illuviation of clay minerals and iron components. Stage III bears evidence of intensive frost weathering processes along with a development of sand-wedge cast that occurred under permafrost conditions of Weichselian glaciation (MIS 2). Further, the formation of an aeolian pavement along with erosion and removal of upper portions of the Eemian soil took place. An accumulation of coeval-age coversand led to a restart of the pedological clock, which resulted in the development of cryogenically-altered soils. Stage IV ushered the Holocene (MIS 1) delivering a progressive development of pedogenesis together with plant-succession modifications under temperate-climate conditions. Brunic or luvic horizons are formed within the coversand deposits. Pseudogleyic conditions developed along the lithological discontinuities. Of note, sand-wedge casts (SWCs) occupy specific locations in the profiles that provide the vertical migration of water, base components and iron. Reducing conditions that occurred along the bottom of ice-wedge casts favour the formation of gleyic horizon. SWCs that have been ultimately modified into ice-wedge and acted as sand-wedge-affiliated geochemical corridors (SWGC) for shallow groundwater forming natural tile drains.

Glaciotectonic disintegration of roches moutonnées during glacial ripping in east Sweden

ABSTRACT Roches moutonnées are typical landforms of glacial erosion developed in hard rocks, with an asymmetric profile caused by abrasion and lee-side plucking. In eastern Sweden, some roches moutonnées show extensive damage, including open fractures, disintegration into blocks, fracture caves and short boulder trains. Disintegration increases along ice-flow directions during deglaciation of the last Weichselian Fennoscandian Ice Sheet, indicating a subglacial origin: limited edge rounding can be explained by a combination of hard rock, slow abrasion rates and disintegration just prior to deglaciation. The roches moutonnées initially developed in kernels of gneissic rocks with a wide fracture spacing (large block size) and interlocking fracture pattern, and hence high overall rock mass strength. Dilated fractures and ‘fracture caves’ occur up to 15 m below the ice-bed interface. It is proposed that hydraulic jacking by overpressured water opened up the rock mass along pre-existing fractures. Jacking reduced rock mass strength, allowing glaciotectonic deformation of the roches moutonnées. Uneven hydraulic jacking led to uplift of individual fracture-bound blocks above the pre-existing smooth, abraded surface of the roches moutonnées, creating blunt, step-like edges. These edges allowed high ice pushing forces to act on large blocks: where blocks extend into the deeper rock mass, they further aided the disintegration of the rock mass. The disintegrated roches moutonnées can be regarded as transient features between intact bedrock and complete disintegration into boulders. The jacking-disintegration-transport sequence is characteristic of glacial ripping and very different from classic lee-side plucking.

Potential drivers of disparity in early Middle Pleistocene interglacial climate response over Eurasia

Poor understanding of the differential evolution of interglacial climate over various regions in Eurasia greatly limits our ability to predict the specific local impacts of future climate change. Here we demonstrate starkly opposing trends in interglacial intensities in Asia and Europe over the early Middle Pleistocene and Mid-Brunhes Transitions based on the study of various climate proxies in the Paks loess record (Hungary), a key profile in the European Loess Belt. These contrasting climate trends imply major but unexplained differences in the regional response of Quaternary climate to shifts in forcing mechanisms during two major reorganizations of Earth's climate across Eurasia. Using new rock magnetic datasets from the Paks record we suggest that the changing dominance of Mediterranean, Atlantic (Westerlies) and continental air masses under the influence of the quasi permanent high pressure centres over the Fennoscandian ice sheet played a key role in the observed differences in the evolution of early Middle Pleistocene climate. We also propose that the intensification of the early Middle Pleistocene glaciation of Central Asian mountains may have strengthened and shifted the Siberian High westward on multi-millennial/orbital timescales, which in turn forced geographically contrasting expressions of the MIS 19 to 11 interglacials in Eurasia. This is the first coherent explanation for the geographically diverse response of regional climate to the early Middle Pleistocene climate transitions and points to the clear role of global cooling and expanded mid-latitude glaciers in driving these events. • Wind patterns and precipitation are reconstructed from a European key loess sequence. • Disparity in early Middle Pleistocene interglacial climate response over Eurasia is recognized. • Potential drivers of the disparity are proposed, and a conceptual climate model was built.

Maximální rozšíření pleistocenního ledovcového štítu ve Vnějších Západních Karpatech a Východních Sudetech – současné znalosti a nový důkaz

A single find of erratic Baltic flint in the Slezské Beskydy Mts. represents a rare evidence of glaciation in the area. The find originates from fluvial deposits of Weichselian valley terrace of the Hluchová Stream, tributary of the Olza River. The erratic was reworked from Elsterian sediments deposited by the Scandinavian Ice Sheet in a very marginal position near its maximum extent in Outer Western Carpathians. Based on morphological analysis, the surface of the erratic was abraded by wave action in the source area on the Baltic shore and subsequently by glacial and (glacio)fluvial transport. Position of the find proves the glaciation of the Olza River valley by a greater ice lobe, than previously assumed.

Inventory of glacial curvilineations (GCLs) at the southern periphery of the last Scandinavian Ice Sheet

Glacial curvilineations (GCLs) are enigmatic landforms recently discovered in the southern peripheral area of the last Scandinavian Ice Sheet (SIS). They consist of multiple parallel, winding ridges several metres high and up to several kilometers long, systematically arranged in groups of common geomorphic appearance. GCLs create spectacular landscapes whose origin generates much controversy. Expanding from the first GCLs finding in central Poland, we now present an inventory of other GCLs occurrences along the SIS margin in Poland, Germany and Denmark based on analysis of high-resolution LiDAR digital elevation models. We document a total of 59 glacial curvilineation fields and focus on the geomophic characteristics of the 43 most conspicuous ones. Each of these fields consists of smaller units called GCL swarms, 137 in total. The fields differ strongly in areal coverage (between 1 and 358 km2), number of swarms (between 1 and 21) and number of individual landforms (between 2 and 255). All GCL fields occur in tunnel valleys terminating at the ice margin, typically at various topographical levels. Our findings, interpreted in the context of glacial geological processes, support the origin of GCLs as erosional remnants of antecedent landscapes carved by channelized subglacial meltwater flows.

Diversity of deformation structures in varves of the Weichselian Pyrzyce Basin (NW Poland)

Within the sediments of the Pyrzyce Ice-dammed Lake (Polish Lowlands), a number of soft-sediment deformation structures (SSDS) in varve series were distinguished. The unique morphology of the ice-dammed lake—the occurrence of level IV overload and the presence of dead ice in the central zone (Lake Miedwie) and the marginal reservoir—had a decisive impact on the formation and nature of sedimentary structures. Dead ice melting and loss of support could be the cause of plastic deformation. The highest intensity deformation in layered sediments of both plastic and rigid types was observed in the marginal (highest) part of the reservoir. Their formation is associated with the melting of rooted lumps of dead ice in sub-aquatic conditions. The exception is deformation resulting from desiccation under subaerial conditions. The genesis of the deformation observed in rhythmically-layered varves is associated with the last stage of functioning and disappearance of the reservoir during the recession of an ice sheet of the Pomeranian phase of the Weichselian glaciation. These studies detail the current state of knowledge on the processes and factors responsible for the formation of soft structures occurring in glaciolimnic basin.

Detrital zircon U–Pb age analysis of last glacial loess sources and proglacial sediment dynamics in the Northern European Plain

Loess deposits along the northern fringe of the European loess belt potentially record past changes in dust emission from areas proximal to former ice sheets. Recent chronologies from loess deposits across this region generally agree on greatly enhanced dust deposition rates when the Fennoscandian Ice Sheet reached its maximum extent during the late last glacial. However, uncertainties over the material's source and origin limit understanding of the causes of this enhanced dust activity. In particular, loess in southwestern Poland has been attributed to multiple origins, mainly involving glaciofluvial outwash plains along the Fennoscandian Ice Sheet margin and/or local sources in the mountainous areas of the Sudetes and Western Carpathians. Here we apply detrital zircon U–Pb age analyses for a large number of grains recovered from four loess samples taken from different stratigraphic units exposed at Biały Kościół in southwestern Poland, previously luminescence dated to MIS 4–2, to assess loess provenance as well as its temporal evolution during last glacial Fennoscandian Ice Sheet fluctuations. Furthermore, we analysed the detrital zircon U–Pb age spectra of five samples from potential source sediments to constrain the history of sediment recycling and mixing within the Northern European Plain prior to deflation and loess deposition. The broad range of zircon age components detected in the four loess samples suggests both Fennoscandian and closer Peri-Gondwanan proto-sources while a narrow, dominant Carboniferous age peak is consistent with sourcing from the local Strzelin Hills in the Sudetic foreland. However, the presence of both Fennoscandian and Peri-Gondwanan derived grains in samples from potential source sediments reveals that this mixture of sediment sources is widespread across the Northern European Plain, as a result of long-term glacial and fluvial reworking of cover sediments in the proglacial area throughout the Quaternary. Local rivers draining the Sudetic foreland transported this Fennoscandian-Peri-Gondwanan sediment mixture along with particles denuded from the Strzelin Hills, resulting in a nearby, temporally stable dust source for the Biały Kościół loess during MIS 4–2, while dust emission rates were substantially increased during the last glacial maximum. Given that our model for loess formation at Biały Kościół essentially involves sediment distribution via rivers prior to short distance aeolian transport, we infer that the proportion of northern ice sheet derived particles in European loess deposits is mainly controlled by the drainage pattern of major rivers in relation to Pleistocene ice margins where glaciofluvial sediment is abundant. Based on the presence of Fennoscandian derived zircon grains in European loess deposits, we constrain a southern limit of the influence of northern ice sheet dust sources along the central European highlands that currently divide drainage between the Northern European Plain and the Danube Basins.

Deglaciation of the Scandinavian Ice Sheet and a Younger Dryas ice cap in the outer Hardangerfjorden area, southwestern Norway

Understanding past responses of ice sheets to climate change provides an important long‐term context for observations of present day, and projected future, ice‐sheet change. In this work, we reconstruct the deglaciation of the marine‐terminating western margin of the Scandinavian Ice Sheet in the outer Hardangerfjorden area of southwestern Norway, following the Last Glacial Maximum (LGM) until the start of the Holocene. We base our interpretations on a combination of geomorphological mapping using high‐resolution (LiDAR) terrain models, 68 new cosmogenic nuclide 10Be exposure ages and radiocarbon‐dated lake sediment cores, supported by the stratigraphic position of the 12.1 ka Vedde Ash. We show that even the highest mountain summits in the area (~1200–1400 m a.s.l.) were ice‐covered during the LGM, thus settling debates concerning the Scandinavian Ice Sheet thickness in this region. These summits emerged as nunataqs through the ice sheet about 22–18 ka, potentially owing to upstream ice thinning caused by the break‐up and retreat of the Norwegian Channel Ice Stream. Following the break‐up of the Norwegian Channel Ice Stream, the ice margin seemingly stabilized at the outermost coast for 3500–5500 years before the mouth of Hardangerfjorden became ice free at c. 14.5 ka. Subsequently, during the Bølling and Allerød periods, the ice sheet retreated rapidly into the inner parts of Hardangerfjorden before a major ice sheet re‐advance during the Younger Dryas. We identify and reconstruct a sizeable, independent ice cap on the Ulvanosa mountain massif during the Younger Dryas (YD), a massif that earlier was mapped as covered by the Scandinavian Ice Sheet during the YD. We also document ice‐free areas that are more extensive than previously thought between Hardangerfjorden and Matersfjorden during the YD.

Sediment sequence at Muhos, western Finland – a window to the Pleistocene history of the Scandinavian Ice Sheet

In this study, we analysed a ~54‐m sediment core consisting of Quaternary sediments overlying the Neoproterozoic Muhos Formation in central western Finland, adjacent to the Gulf of Bothnia. The sediments recovered were logged, and their sedimentological characteristics defined. Two fine‐grained sediment units were subjected to biostratigraphical studies using pollen and diatom analyses. In addition, two sand‐rich units and a wooden stick were dated by the optically stimulated luminescence (OSL) and 14C‐AMS methods. The core sediments were divided into six units, where several diamicton, sand and gravel, and silt‐and‐clay‐dominated beds were studied. The results indicate that the sediment succession of the core beneath the Holocene Litorina Sea and the Ancylus Lake sediments of the Baltic Basin were deposited in glacial and lacustrine environments that existed in the Oulu River valley during the time period between the Saalian glaciation (MIS 6) and the Holocene. The stratigraphical evidence, supported by the OSL ages, suggests that the Scandinavian Ice Sheet (SIS) entered the Muhos area during the Saalian glaciation, and at least during three separate time intervals in the Weichselian stage. Stratigraphically controlled and age‐bracketed evidence shows that the first Weichselian SIS advance extended further south in the eastern part of Fennoscandia than previously estimated, and that this ice growth phase occurred during the Early Weichselian Herning Stadial (MIS 5d). The subsequent ice growth phases occurred during the Middle (MIS 4) and Late (MIS 2) Weichselian substages. The lacustrine and littoral sediments of the Muhos core were correlated with the late Eemian interglacial (MIS 5e) and two Weichselian interstadials (MIS 5c and MIS 3).

Western Siberia experienced rapid shifts in moisture source and summer water balance during the last deglaciation and early Holocene

ABSTRACTThe Russian Arctic is an extensive region, with relatively few long‐duration paleoclimate reconstructions compared to other terrestrial Arctic regions. We present a 24 000‐year reconstruction of climate in the Polar Ural Mountains using n‐alkanoic acid hydrogen isotopes from Lake Bolshoye Shchuchye. Major last deglaciation climate changes in the North Atlantic are present in this record, including transitions associated with the Bølling–Allerød, Younger Dryas and Holocene. However, the magnitude of the last deglaciation isotopic shifts at Bolshoye Shchuchye are small relative to the North Atlantic, and are dwarfed by a shift to 2H‐enriched values starting at 10.5k cal a bp at this site that is not present in most other records. The last deglaciation changes may be due to variations in local temperature, sea ice cover in the Barents and Kara seas, and plant community shifts impacting transpiration. The enrichment starting at 10.5‐k cal a bp probably records a shift towards modern climate conditions, caused by the loss of the Scandinavian Ice Sheet, increased locally sourced moisture from the Barents and Kara seas, and northward treeline migration causing enhanced transpiration. Future warming may increase summer precipitation in this region, with changes to local ecosystems and carbon cycling.

Freshwater routing in eddy-permitting simulations of the last deglacial: the impact of realistic freshwater discharge

Abstract. Freshwater, in the form of glacial runoff, is hypothesized to play a critical role in centennial- to millennial-scale climate variability, such as the Younger Dryas and Dansgaard–Oeschger events, but this relationship is not straightforward. Large-scale glacial runoff events, such as Meltwater Pulse 1a (MWP1a), are not always temporally proximal to subsequent large-scale cooling. Moreover, the typical design of hosing experiments that support this relationship tends to artificially amplify the climate response. This study explores the impact that limitations in the representation of runoff in conventional “hosing” simulations has on our understanding of this relationship by examining where coastally released freshwater is transported when it reaches the ocean. We particularly focus on the impact of (1) the injection of freshwater directly over sites of deep-water formation (DWF) rather than at runoff locations (i.e. hosing), (2) excessive freshwater injection volumes (often by a factor of 5), and (3) the use of present-day (rather than palaeo) ocean gateways. We track the routing of glaciologically constrained freshwater volumes from four different inferred injection locations in a suite of eddy-permitting glacial ocean simulations using the Massachusetts Institute of Technology General Circulation Model (MITgcm) under both open and closed Bering Strait conditions. Restricting freshwater forcing values to realistic ranges results in less spreading of freshwater across the North Atlantic and indicates that the freshwater anomalies over DWF sites depend strongly on the geographical location of meltwater input. In particular, freshwater released into the Gulf of Mexico generates a very weak freshwater signal over DWF regions as a result of entrainment by the turbulent Gulf Stream. In contrast, freshwater released into the Arctic with an open Bering Strait or from the Eurasian ice sheet is found to generate the largest salinity anomalies over DWF regions in the North Atlantic and GIN (Greenland–Iceland–Norwegian) seas region respectively. Experiments show that when the Bering Strait is open, the Mackenzie River source exhibits more than twice as much freshening of the North Atlantic deep-water formation regions as when the Bering Strait is closed. Our results illustrate that applying freshwater hosing directly into the North Atlantic with even “realistic” freshwater amounts still overestimates the amount of terrestrial runoff reaching DWF regions. Given the simulated salinity anomaly distributions and the lack of reconstructed impact on deep-water formation during the Bølling–Allerød, our results support that the majority of the North American contribution to MWP1a was not routed through the Mackenzie River.

Case Study: Road Embankment Foundation on Weak Organic Soils

The paper presents the problem of embankment foundation during reconstruction and extension of regional road 110 on the section Witomierz - Grądy. On this section the road crosses the valley of meltwater formed during the last Weichselian glaciation. It is a watershed area between the streams Stuchowska Struga and Otoczka Reska. In the Holocene a 961 ha fen peat called Wielki Smogorze has been formed. There are favorable conditions in the area for the formation of fen peats. The thickness of the peats is over 4 m, reaching almost 8.0 m in its peak. Since 1964 the Przybiernówko-Grądy II deposit (402 ha) is continuously exploited. Currently, the deposit is being used by a company called "Lasland". Extraction is conducted on the basis of a relevant concession within a designated mining area of 242.9 ha. The mining concession is valid until the end of 2030. The pits are deep and are located on both sides of the regional road no. 110. Peat is transported by narrow-gauge railroad to the nearby processing plant, where it is sieved, sorted and packed. Based on the analysis of available archival materials, the road as found today was functioning already in the middle of the 19th century. In the 1970s the road was widened to 6.0m In 2003 due to the bad condition of the surface the asphalt layer was renovated by applying a grid and new asphalt layers. During the renovation longitudinal cracks have been reported and there were problems with the compaction of the mix. The direct cause was the shallow layer of peat located just 1.1 - 1.3 m under the road. Conducted renovation did not bring expected results, so in 2019 the documentation for the reconstruction of the road was prepared. Different methods of road foundation were analyzed, from soil replacement through the use of piles. In the end the decision was made to directly settle the embankment with the use of geosynthetics. This study presents a selected solution and shows the results of calculations. Changes during the execution of the reconstruction were discussed. The applied solution allowed for simultaneous functioning of the mining plant and reconstruction of the regional road on the section of 1.2 km. geotextiles and geonets were used. The main argument for such solution was relatively low cost of reconstruction and time of execution. During construction there were problems with obtaining the parameters on the substructure for asphalt layers. Two alternative solutions were proposed, one of which required increasing the thickness of embankment reinforced with geo-mesh and an additional layer of CBGM mix substructure. The solution allowed to meet the design requirements and complete the reconstruction according to the plan without any problems.

Palaeoenvironmental changes and their chronology during the latter half of MIS 5 on the south-eastern coast of the Gulf of Finland

This paper presents an overview of the results obtained by collaborative use of integrative multidisciplinary investigation of a highly unique geological section discovered in the south-eastern part of the Gulf of Finland. The main aim of this study is to identify the Late Pleistocene palaeoenvironmental and palaeoclimatic changes in the study area during one of the most controversial periods in the Late Pleistocene history, which is the MIS 5. Our approach to achieving the goal is based on the investigation of feldspar-based infrared optically stimulated luminescence (IR-OSL) dated records of pollen-induced palaeoclimate and palaeoenvironmental proxies. Palynological analysis and IR-OSL dates of 29 samples for the first time produced a detailed, characterised and dated sequence of palynological zones M6 to M8 ranging from ca. 94 ka to ca. 70 ka and made it possible to determine the chronological sequence of climate changes in this interval. Representative pollen spectra from 69 samples provided convincing evidence of a warm interval of interglacial rank at least during the latter half of MIS 5. The data obtained in the present study from well-preserved deposits of the second half of MIS 5 provided a novel insight into the palaenvironmental conditions during this period. These results seem to contradict the widely held view according to which the last interglacial stage with the climate conditions similar or warmer than today corresponds to only a small part of marine isotope stage (MIS) 5 of the deep-sea stratigraphy, i.e. MISs 5e lasting ca. 11 ka, and the period from 100 ka to 80 ka is attributed to an Early Weichselian glaciation, which in Northern Eurasia reached its maximum at about 90 ka. However, the findings obtained during the present study are also convincingly confirmed by those of our previous long-term investigations of warm climate-related shell-bearing deposits of various origins from Northern Eurasia including climatically highly sensitive Eurasian Arctic palaeo-shelf area. These results indicated that the overwhelming majority of the warm climate-related dates (ca. 82%) is also concentrated in the latter half of MIS 5 in the time range between 110 ka and 70 ka.

Climate and ice sheet evolutions from the last glacial maximum to the pre-industrial period with an ice-sheet–climate coupled model

Abstract. The last deglaciation offers an unique opportunity to understand the climate–ice-sheet interactions in a global warming context. In this paper, to tackle this question, we use an Earth system model of intermediate complexity coupled to an ice sheet model covering the Northern Hemisphere to simulate the last deglaciation and the Holocene (26–0 ka). We use a synchronous coupling every year between the ice sheet and the rest of the climate system and we ensure a closed water cycle considering the release of freshwater flux to the ocean due to ice sheet melting. Our reference experiment displays a gradual warming in response to the forcings, with no abrupt changes. In this case, while the amplitude of the freshwater flux to the ocean induced by ice sheet retreat is realistic, it is sufficient to shut down the Atlantic meridional overturning circulation from which the model does not recover within the time period simulated. However, with reduced freshwater flux we are nonetheless able to obtain different oceanic circulation evolutions, including some abrupt transitions between shut-down and active circulation states in the course of the deglaciation. The inclusion of a parameterisation for the sinking of brines around Antarctica also produces an abrupt recovery of the Atlantic meridional overturning circulation, absent in the reference experiment. The fast oceanic circulation recoveries lead to abrupt warming phases in Greenland. Our simulated ice sheet geometry evolution is in overall good agreement with available global reconstructions, even though the abrupt sea level rise at 14.6 ka is underestimated, possibly because the climate model underestimates the millennial-scale temperature variability. In the course of the deglaciation, large-scale grounding line instabilities are simulated both for the Eurasian and North American ice sheets. The first instability occurs in the Barents–Kara seas for the Eurasian ice sheet at 14.5 ka. A second grounding line instability occurs ca. 12 ka in the proglacial lake that formed at the southern margin of the North American ice sheet. With additional asynchronously coupled experiments, we assess the sensitivity of our results to different ice sheet model choices related to surface and sub-shelf mass balance, ice deformation and grounding line representation. While the ice sheet evolutions differ within this ensemble, the global climate trajectory is only weakly affected by these choices. In our experiments, only the abrupt shifts in the oceanic circulation due to freshwater fluxes are able to produce some millennial-scale variability since no self-generating abrupt transitions are simulated without these fluxes.