The Baltic Ice Stream, a large fast-flowing sector of the last Fennoscandian Ice Sheet that occupied the present-day Baltic Sea basin, was first conceptualized in the earliest days of glacial geological research in Scandinavia. Landform and sedimentological evidence from the terrestrial margins support the concept and numerical ice-sheet models demonstrate its existence and possible evolution. However, with evidence for the Baltic Ice Stream thus far limited to the terrestrial periphery, its true form, scale, function, and role in deglaciation have proven enigmatic. Here we present geomorphological evidence directly from the Baltic seabed that confirms the existence of and sheds light on the behaviour of the Baltic Ice Stream. Based on an extensive collection of high-, moderate- and low-resolution bathymetric terrain models covering a large proportion of the Baltic Sea floor, and complemented by LiDAR-data for the Baltic islands, we have identified and mapped >20 000 individual subglacial bedforms, meltwater landforms and grounding line landforms. We reconstruct a six-stage sequence of ice flow and retreat, finding that streaming was persistent in the Baltic but that pathways were variable in extent, timing and duration: different sectors of the Baltic exhibit asynchronous streaming and out-of-phase grounding line changes. During deglaciation, grounding line readvances occurred in both the southwestern and the northern Baltic Proper, and, while abundant iceberg ploughmarks attest to calving as a significant ice loss mechanism, lobate margins suggest supply to the Baltic catchment was consistently high. Our reconstruction is limited by a fragmentary geomorphic record. Here we put forward a first hypothesis for how the Baltic Ice Stream evolved, and hope it stimulates new geomorphic, stratigraphical and core data collection to extend the landform record, provide insights into subglacial and grounding line processes, and constrain the chronology for Baltic Ice Stream flow and retreat.

Extensive glacial lakes dammed in the Scandinavian Mountains during the retreat of the last Scandinavian Ice Sheet were first hypothesised over a century ago. Here, using high-resolution LiDAR, we report >4500 relict shorelines, deltas and palaeo-channels related to ice-dammed lakes over a -30 000 km2 area of central Jämtland, west-central Sweden. Shorelines occur as flights on the valley sides, a consequence of sequential lowering of palaeo-lake levels during ice margin retreat and lower threshold outlets becoming ice-free. Based on the extent and elevation of shorelines, we identify requisite lake-damming ice-margin positions and lake drainage outlets, and we reconstruct the coupled evolution of ice-dammed lakes and the retreating ice margin. Beginning as a series of smaller ice-dammed lakes along the Swedish-Norwegian border, draining westward across the present-day water divide and into the Atlantic Ocean, the lakes successively coalesced during eastward ice margin retreat to form water bodies covering 1000s of km2 with 10s of km-long calving margins. Ultimately, the lake system coalesced into a single lake: the Central Jämtland Ice Lake, which exceeded 3500 km2 in area and 360 km3 in volume. Eventually, the damming ice-margin split in two, resulting in a large (-200 km2) catastrophic glacial lake outburst flood (GLOF) that reversed the drainage of the entire lake system from the west to an eastern outlet draining to the Baltic basin. We present new radiocarbon ages for one lake drainage event prior to the eastward outburst flood and, together with previously published deglacial ages and local varve records, we suggest that the region was possibly deglaciated within just 350 years, sometime between 10.5 and 9.2 cal ka BP. We tentatively correlate the penultimate drainage of the Central Jämtland Ice Lake to the zero-varve of the Swedish Time Scale, a drainage varve at Döviken, eastern Jämtland, raising the tantalising prospect of using the evolution of the ice-dammed lake system to tie the varve-based Swedish Time Scale to the radiocarbon timescale with a new programme of radiocarbon dating in central Jämtland. 

Here we present geomorphological evidence of a previously unrecognised ∼50 km long, ice-marginal moraine complex in southwestern Skåne, southernmost Sweden, which we name the “Lund Moraine”. This lobate moraine marks a sharp boundary between heavily streamlined and gently undulating landscapes, and closely outlines the extent of the “Lund till/diamicton”. We interpret that the moraine was formed by a northward readvance, corresponding to a last Young Baltic readvance of the Scandinavian Ice Sheet into Öresund. Consequently, we infer that the “Lund till/diamicton” was formed subglacially, in contrast to earlier interpretations of it being a waterlain diamicton. Based on previously published dates, stratigraphically below “Lund till/diamicton”, we infer that this readvance occurred sometime after c. 16 cal. ka BP. This readvance could offer an explanation to the apparent discrepancy of observations of the postglacial marine limit from outside and inside the Lund Moraine. Our observations will hopefully settle the ∼50 years long controversy concerning the extent or even existence of such a readvance into Öresund. We expect that our findings will guide further work towards disentangling the complex deglacial history of Skåne and the wider Öresund region.

Meltwater drainage beneath ice sheets is a fundamental consideration for understanding ice–bed conditions and bed-modulated ice flow, with potential impacts on terminus behavior and ice-shelf mass balance. While contemporary observations reveal the presence of basal water movement in the subglacial environment and inferred styles of drainage, the geological record of former ice sheets, including sediments and landforms on land and the seafloor, aids in understanding the spatiotemporal evolution of efficient and inefficient drainage systems and their impact on ice-sheet behavior. We highlight the past decade of advances in geological studies that focus on providing process-based information on subglacial hydrology of ice sheets, how these studies inform theory, numerical models and contemporary observations, and address the needs for future research.

The Swedish Varve Chronology is an unparalleled tool for linking the deglacial history of Sweden with associated palaeo-environmental change at an annual time scale, and it forms part of Sweden's cultural heritage. A full deglacial chronology connected to the present day does not yet exist; a notable gap is in southeasternmost Sweden, where few varved records are successfully connected to reconstruct ice-margin retreat. Deglaciation in southern Sweden covers both the climate transition to the Bølling warm period (similar to 14.7 ka BP) and the ice-margin transition from a subaqueous to terrestrial terminus. To facilitate investigations into the links between ice-margin dynamics and abrupt climate change, we revisited the varve chronologies of southern Sweden. We digitized unpublished records, reanalysed existing varve thickness records, and obtained and analysed new varve series both on land and offshore. This combined approach has enabled us to refine and extend the existing south coast chronology east and 78 km northwards. Our new Skåne-Småland chronology records 725 years of deglaciation, in addition to a younger floating chronology in parts. This chronology suggests that the glacial-lake terminating Fennoscandian Ice Sheet in southern Sweden initially retreated northwards at similar to 110-160 m a(-1) slowing to 60-70 m a(-1) near the palaeo-shoreline. Between today's mainland and the (now) island of Öland the retreat rates increase three- to fivefold. Ice-margin retreat was initially oriented towards the north (as along the south coast), but later pivoted towards the northwest, signifying a landward retreat of terrestrial 'Swedish' ice that became divorced from the Baltic Sea ice-sheet catchment. Our new 725-year-long varve thickness series reveals repeated multidecadal scale episodes of increased sedimentation. These likely signify phases of enhanced ice-sheet melting that repeat and persist throughout the deglaciation of Skåne-Småland.

Projections of ice sheet behavior hinge on how ice flow velocity evolves and the extent to which marine-based grounding lines are stable. Ice flow and grounding line retreat are variably governed by the coupling between the ice and underlying terrain. We ask to what degree catchment-scale bed characteristics determine ice flow and retreat, drawing on paleo-ice sheet landform imprints from 99 sites on continental shelves worldwide. We find that topographic setting has broadly steered ice flow and that the bed slope favors particular styles of retreat. However, we find exceptions to accepted rulesof behavior: Regional topographic highs are not always an impediment to fast ice flow, retreat may proceed in a controlled, steady manner on reverse slopes and, unexpectedly, the occurrence of ice streaming is not favored on a particular geological substrate. Furthermore, once grounding line retreat is under way, readvance is rarely observed regardless of regional bed characteristics.

Palaeo-ice sheets are important analogues for understanding contemporary ice sheets, offering a record of ice sheet behaviour that spans millennia. There are two main approaches to reconstructing palaeo-ice sheets. Empirical reconstructions use the available glacial geological and chronological evidence to estimate ice sheet extent and dynamics but lack direct consideration of ice physics. In contrast, numerically modelled simulations implement ice physics, but often lack direct quantitative comparison with empirical evidence. Despite being long identified as a fruitful scientific endeavour, few ice sheet reconstructions attempt to reconcile the empirical and model-based approaches. To achieve this goal, model-data comparison procedures are required. Here, we compare three numerically modelled simulations of the former British–Irish Ice Sheet with the following lines of evidence: (a) position and shape of former margin positions, recorded by moraines; (b) former ice-flow direction and flow-switching, recorded by flowsets of subglacial bedforms; and (c) the timing of ice-free conditions, recorded by geochronological data. These model–data comparisons provide a useful framework for quantifying the degree of fit between numerical model simulations and empirical constraints. Such tools are vital for reconciling numerical modelling and empirical evidence, the combination of which will lead to more robust palaeo-ice sheet reconstructions with greater explicative and ultimately predictive power.

Realistic characterization of subglacial hydrology necessitates knowledge of the range in form, scale, and spatiotemporal evolution of drainage networks. A relict subglacial meltwater corridor on the deglaciated Antarctic continental shelf encompasses 80 convergent and divergent channels, many of which are hundreds of meters wide and several of which lack a definable headwater source. Without significant surface-melt contributions to the bed like similarly described landforms in the Northern Hemisphere, channelized drainage capacity varies non-systematically by three orders of magnitude downstream. This signifies apparent additions and losses of basal water to the bed-channelized system that relates to bed topography. Larger magnitude grounding-line retreat events occurred while the channel system was active than once channelized drainage had ceased. Overall, this corridor demonstrates that meltwater drainage styles co-exist in time and space in response to bed topography, with prolonged impacts on grounding-line behavior.

Baltic Sea bathymetric properties are analysed here using the newly released digital bathymetric model (DBM) by the European Marine Observation and Data Network (EMODnet). The analyses include hypsometry, volume, descriptive depth statistics, and kilometre-scale seafloor ruggedness, i.e. terrain heterogeneity, for the Baltic Sea as a whole as well as for 17 sub-basins defined by the Baltic Marine Environment Protection Commission (HELCOM). We compare the new EMODnet DBM with IOWTOPO the previously most widely used DBM of the Baltic Se aproduced by the Leibniz-Institut fur Ostseeforschung Warnemiinde (IOW), which has served as the primary gridded bathymetric resource in physical and environmental studies for nearly two decades. The area of deep water exchange between the Bothnian Sea and the Northern Baltic Proper across the Aland Sea is specifically analysed in terms of depths and locations of critical bathymetric sills. The EMODnet DBM provides a bathymetric sill depth of 88 m at the northern side of the Aland Sea and 60 m at the southern side, differing from previously identified sill depths of 100 and 70 m, respectively. High-resolution multibeam bathymetry acquired from this deep water exchange path, where vigorous bottom currents interacted with the seafloor, allows us to assess what presently available DBMs are missing in terms of physical characterization of the seafloor. Our study highlights the need for continued work towards complete high-resolution mapping of the Baltic Sea seafloor.

Groundwater discharge into lakes is an important component of the fluid and nutrient budgets, and a possible route for contaminant transport. However, groundwater flow beneath lakes is difficult to investigate due to the need for drilling deep boreholes. In 2012, a 2,000 m deep borehole was drilled in Lake Vättern, the second largest lake in Sweden. A continuous temperature profile was collected from the borehole. The geothermal gradient in the upper 180 m is highly non-linear, and not controlled by variability in the measured thermal properties of the sediments and rocks. The anomalous temperature profile is best explained by fluid flow into the borehole and subsequent vertical flow of warm waters towards the lake floor. Combining the temperature profile with stratigraphic information from drilling logs and seismic data, we find that fluid flow into the borehole occurs in glacial and glaciofluvial sediments deposited on top of a large sandstone aquifer (the Visingsö Group). The warm waters flowing through the glacial and glaciofluvial sediments are likely sourced from the underlying Visingsö Group sandstones. There is no evidence for substantial vertical migration of these waters through the overlying glacial and postglacial sediments. We speculate that they escape either along lake margins where overlying sediments become thinner, or along faults that are known to exist in the deeper basin.These results highlight an important hydraulic transport pathway between recognised regional aquifers and Lake Vättern. Further work is needed to evaluate the significance of groundwater discharge on the water and nutrient budget of the lake.