Drainage Events Research Articles

Step-change in supraglacial pond area on Tshojo Glacier, Bhutan, and potential downstream inundation patterns due to pond drainage events

Step-change in supraglacial pond area on Tshojo Glacier, Bhutan, and potential downstream inundation patterns due to pond drainage events

Evidence of Ecosystem Tipping Point on St. Lawrence Island: Widespread Lake Drainage Events After 2018

AbstractInfluenced by climate change, numerous lakes in permafrost regions are draining, showing significant spatial variability. This study focuses on St. Lawrence Island, where over the last two decades, 771 of 3,271 lakes have drained—a rate around 40 times higher than across the entire northern permafrost region. The surge in lake drainage began in 2018, coinciding with record low sea ice extent in the Bering Sea and unprecedented bird mortalities. Using satellite imagery and machine learning methods, we analyzed drainage events to identify the climatic drivers and potential climate thresholds affecting the island's lake ecosystems. Our findings indicate that autumn peak temperatures above 6°C more than triple the drainage probability, and warming‐induced permafrost thawing may be the direct driver of lake drainage. This research highlights the vulnerability of Arctic lake ecosystems to climate change and assists in developing predictive models for permafrost response, crucial for mitigating impacts on Arctic communities.

Insights into glacial processes from micromorphology of silt-sized sediment

Abstract. Silt-rich meltwater plume deposits (MPDs) analyzed from marine sediment cores have elucidated relationships that are clearly connected, yet difficult to constrain, between subglacial hydrology, ice-marginal landforms, and grounding-zone retreat patterns for several glacial catchments. Few attempts have been made to infer details of subglacial hydrology, such as flow regime, geometry of drainage pathways, and mode(s) of sediment transport through time, from grain-scale characteristics of MPDs. Using sediment samples from MPD, till, and grounding-zone proximal diamicton collected offshore of six modern and relict glacial catchments in both hemispheres, we examine grain shape distributions and microtextures (collectively, grain micromorphology) of the silt fraction to explore whether grains are measurably altered from their subglacial sources via meltwater action. We find that 75 % of all imaged grains (n = 9400) can be described by 25 % of the full range of measured shape morphometrics, indicating grain shape homogenization through widespread and efficient abrasive processes in subglacial environments. Although silt grains from MPDs exhibit edge rounding more often than silt grains from tills, grain surface textures indicative of fluvial transport (e.g., v-shaped percussions) occur in only a modest number of grains. Furthermore, MPD grain surfaces retain several textures consistent with transport beneath glacial ice (e.g., straight or arcuate steps, (sub)linear fractures) in comparable abundances to till grains. Significant grain shape alteration in MPDs compared to their till sources is observed in sediments from glacial regions where (1) high-magnitude, potentially catastrophic meltwater drainage events are inferred from marine sediment records and (2) submarine landforms suggest supraglacial melt contributed to the subglacial hydrological budget. This implies that quantifiable grain shape alteration in MPDs could reflect a combination of high-energy flow of subglacial meltwater, persistent sediment entrainment, and/or long sediment transport distances through subglacial drainage pathways. Integrating grain micromorphology into analysis of MPDs in site-specific studies could therefore aid in distinguishing periods of persistent, well-connected subglacial discharge from periods of sluggish or disorganized drainage. In the wider context of deglacial marine sedimentary and bathymetric records, a grain micromorphological approach may bolster our ability to characterize ice response to subglacial meltwater transmission through time. This work additionally demonstrates that glacial and fluvial surface textures are retained on silt-sized quartz grains in adequate amounts for microtexture analysis, which has heretofore been conducted exclusively on the sand fraction. Therefore, grain microtextures can be examined on silt-rich glaciogenic deposits that contain little to no sand as a means to evaluate sediment transport processes.

Implications of MOLA topography on the presence of liquid water under the south polar cap of Mars

The existence of liquid water at the base of the south polar layered deposits (SPLD) on Mars has been proposed on the basis of radar data, but dry interpretations have been advanced as well. Here, we consider topography of the SPLD surface above the location of hypothesized water using individual altimetry tracks rather than gridded or interpolated data. Liquid water could create negative surface relief through drainage events and subsequent collapse, or positive surface relief through ice expansion. We find meters of positive topographic relief relative to the surrounding region, consistent with thermal expansion of the ice column from the added geothermal heat that would be required to maintain liquid water at the SPLD base. However, topographic analysis of randomly chosen areas throughout the SPLD suggest that positive topographic relief of this amplitude has a 20% chance of occurring through other processes or statistical noise. Therefore, we conclude that both wet and dry interpretations of the radar data remain viable.

A supervised, externally validated machine learning model for artifact and drainage detection in high-resolution intracranial pressure monitoring data.

In neurocritical care, data from multiple biosensors are continuously measured, but only sporadically acknowledged by the attending physicians. In contrast, machine learning (ML) tools can analyze large amounts of data continuously, taking advantage of underlying information. However, the performance of such ML-based solutions is limited by different factors, for example, by patient motion, manipulation, or, as in the case of external ventricular drains (EVDs), the drainage of CSF to control intracranial pressure (ICP). The authors aimed to develop an ML-based algorithm that automatically classifies normal signals, artifacts, and drainages in high-resolution ICP monitoring data from EVDs, making the data suitable for real-time artifact removal and for future ML applications. In their 2-center retrospective cohort study, the authors used labeled ICP data from 40 patients in the first neurocritical care unit (University Hospital Zurich) for model development. The authors created 94 descriptive features that were used to train the model. They compared histogram-based gradient boosting with extremely randomized trees after building pipelines with principal component analysis, hyperparameter optimization via grid search, and sequential feature selection. Performance was measured with nested 5-fold cross-validation and multiclass area under the receiver operating characteristic curve (AUROC). Data from 20 patients in a second, independent neurocritical care unit (Charité - Universitätsmedizin Berlin) were used for external validation with bootstrapping technique and AUROC. In cross-validation, the best-performing model achieved a mean AUROC of 0.945 (95% CI 0.92-0.969) on the development dataset. On the external validation dataset, the model performed with a mean AUROC of 0.928 (95% CI 0.908-0.946) in 100 bootstrapping validation cycles to classify normal signals, artifacts, and drainages. Here, the authors developed a well-performing supervised model with external validation that can detect normal signals, artifacts, and drainages in ICP signals from patients in neurocritical care units. For future analyses, this is a powerful tool to discard artifacts or to detect drainage events in ICP monitoring signals.

Tracking rapid and slow ice-dammed lake changes through optical satellites and local knowledge: a case study of Tystigbreen in Norway

Abstract The number of glacial lakes has grown globally concurrently with the retreat of glaciers in the last few decades, increasing the risk of potentially hazardous glacial lake outburst floods (GLOFs) and posing a threat to downstream communities. Norway has several known ice-dammed lakes that produce repeated GLOFs but as we show here, the existing GLOF database is incomplete and needs to be improved through continuous monitoring of glaciers and glacial lakes. This study examines the case of an ice cap in central Norway hosting at least four drainage-prone lakes. We reconstruct the sequence of lake drainage events through a combination of remote sensing, ground-truthing and citizen science while scrutinising the applicability of the PlanetScope imagery vs Sentinel-2 and Landsat-8 OLI products. As opposed to the Landsat imagery that often fails to resolve even the largest glacial lakes of Tystigbreen, both PlanetScope and Sentinel-2 are helpful in identifying previously unrecognised glacial lakes and undocumented drainage events. Our analysis suggests that a fusion of the two satellite products may be beneficial for automated tracking of glacial lake changes. We also demonstrate that local knowledge and systematic involvement of citizens in data collection have a potential to enrich GLOF databases.

Multi‐proxy palaeoenvironmental reconstruction of the Skagerrak from the Lateglacial to Middle Holocene

This study uses a multiproxy approach including the first use of 187Os/188Os, %C37:4 biomarkers, carbonate content, sedimentological grain size, geochemical X‐ray fluorescence and microfossil benthic foraminifera species combined with radiocarbon dating, measured on six cores from across the Skagerrak, in order to study the Lateglacial to Middle Holocene history of the area. A new chronostratigraphic framework is developed based on the appearance of specific benthic foraminifera species along with changes in carbonate/X‐ray fluorescence and grain size data. This allowed the correlation of cores based on a series of radiocarbon dated tie points. Analysing the cores together reveals several events recorded in the Skagerrak including: (i) an increased freshwater input (bracketed between 13.3 and 11.3 cal. ka BP) signified by radiogenic 187Os/188Os values, high %C37:4 values and an increase in sand content; (ii) the Glomma drainage event, signified by a sudden appearance of Valvulineria as well as higher %C37:4; and (iii) the opening of the Danish Straits and English Channel leading to the development of modern‐day conditions and circulation patterns in the Skagerrak, signified by the appearance of Hyalinea balthica and a fall in %C37:4.

Downscaling air temperatures for high-resolution niche modeling in a valley of the Amazon lowland forests: A case study on the microclima R package.

The forests of the Amazon basin are threatened by climate and land use changes. Due to the transition towards a drier climate, moisture-dependent organisms such as canopy epiphytes are particularly affected. Even if the topography in the Amazon lowland is moderate, mesoscale nocturnal katabatic flows result from cold air production related to radiative cooling. From a certain level of mass the cold air starts to flow downslope towards the valley centers leading to temperature inversions. The resulting cooling in the valleys drives localized fog formation in the valleys at night. This correlates with high epiphyte abundance and diversity in the valleys, which is much less pronounced upslope. The underlying temperature dynamics are, however, not sufficiently included in coarse-resolution reanalysis models such as ERA5-Land. Since high resolution climate data are needed e.g. for proper niche modeling of locally distributed species such as canopy epiphytes, downscaling models such as microclima have been developed and include micro- and mesoscale effects. However, it is unclear how well the elevation-related diurnal course of air temperature can be simulated. Here, we test functions for downscaling coarse-resolution temperature data to high spatial resolution data implemented in the R-package microclima for the South American tropical lowland forests. To do so we compared microclima-downscaled ERA5-Land air temperature data with meteorological station data. We found that the microclima functions only properly detect 73 temperature inversions out of 412 nocturnal cold air drainage (CAD) events during the dry season study period and only 18 out of 400 during the wet season with default settings. By modifying default values such as the emissivity threshold and time frames of possible CAD condition detection, we found 345 of 412 CAD events during the dry season and 177 out of 400 during the wet season. Despite problems with the distinction between CAD and non-CAD events the microclima algorithms show difficulties in correctly modeling the diurnal course of the temperature data and the amplitudes of elevational temperature gradients. For future studies focusing on temperature downscaling approaches, the modules implemented in the microclima package have to be adjusted for their usage in tropical lowland forest studies and beyond.

Monitoring of Supraglacial Lake Distribution and Full-Year Changes Using Multisource Time-Series Satellite Imagery

Change of supraglacial lakes (SGLs) is an important hydrological activity on the Greenland ice sheet (GrIS), and storage and drainage of SGLs occur throughout the year. However, current studies tend to split SGL changes into melt/non-melt seasons, ignoring the effect of buried lakes in the exploration of drainage, and the existing threshold-based approach to SGL extraction in a synthetic aperture radar (SAR) is influenced by the choice of the study area mask. In this study, a new method (Otsu–Canny–Otsu (OCO)), which accesses the features of SGLs on optical and SAR images objectively, is proposed for full-year SGL extraction with Google Earth Engine (GEE). The SGLs on the Petermann Glacier were monitored well by OCO throughout 2021, including buried lakes and more detailed rapid drainage events. Some SGLs’ extent varied minimally in a year (area varying by 10–25%) while some had very rapid drainage (a rapid drainage event from July 26 to 30). The SGL extraction results were influenced by factors such as the mode of polarization, the surface environment, and the depth of the lake. The OCO method can provide a more comprehensive analysis for SGL changes throughout the year.

A global insight on sensitivity of nitrate leaching to drainage in arable cropping systems

Nitrate (NO3–) leaching leads to widespread deterioration of water quality. NO3– leaching can be regulated by both biochemical and hydrological processes. However, the extent to which drainage accounts for the variation in NO3– leaching and the sensitivity of NO3– leaching to drainage remain unclear, particularly at the global scale. This raises uncertainties in prioritizing the regulation of hydrological (e.g., drainage) processes for reduced NO3– leaching. A global meta-analysis was conducted for arable cropping systems to test two hypotheses: (1) Global variations in NO3– leaching are explained more by drainage and associated influencing factors (e.g., water input); (2) NO3– leaching is more sensitive to the changes in drainage under conditions with lower risks of preferential flow. In this meta-analysis, the average NO3– leaching loss from arable cropping systems was 22.2 kg N ha−1 year−1. Water input, followed by N input, were the most influential in explaining global variations in NO3– leaching. Compared with NO3– concentration, drainage and associated influencing factors explained more of the global variations in NO3– leaching losses. NO3– leaching was more sensitive to drainage in soils with finer texture and higher soil organic carbon content, in regions with low to medium mean annual temperature (<15 ℃) and medium mean annual precipitation (600–1000 mm), where water input during the measurement period was low (<2 mm day−1), nitrogen fertilizer rates were high (>200 kg N ha−1), and where conservation tillage systems (no-tillage or non-inversion tillage) were adopted. Therefore, more attention should be given to avoiding excessive water application in soils characterized by low risk of preferential flow. Preferential flow reduces the chance for drainage water to interact with NO3– in the soil matrix, particularly in areas with high N inputs. While this study has successfully tested our hypotheses, it further highlights the importance of water input management for mitigating global NO3– leaching risks, particularly under environments and changing climate conditions more conducive to drainage events.

Tracking lake drainage events and drained lake basin vegetation dynamics across the Arctic

Widespread lake drainage can lead to large-scale drying in Arctic lake-rich areas, affecting hydrology, ecosystems and permafrost carbon dynamics. To date, the spatio-temporal distribution, driving factors, and post-drainage dynamics of lake drainage events across the Arctic remain unclear. Using satellite remote sensing and surface water products, we identify over 35,000 (~0.6% of all lakes) lake drainage events in the northern permafrost zone between 1984 and 2020, with approximately half being relatively understudied non-thermokarst lakes. Smaller, thermokarst, and discontinuous permafrost area lakes are more susceptible to drainage compared to their larger, non-thermokarst, and continuous permafrost area counterparts. Over time, discontinuous permafrost areas contribute more drained lakes annually than continuous permafrost areas. Following drainage, vegetation rapidly colonizes drained lake basins, with thermokarst drained lake basins showing significantly higher vegetation growth rates and greenness levels than their non-thermokarst counterparts. Under warming, drained lake basins are likely to become more prevalent and serve as greening hotspots, playing an important role in shaping Arctic ecosystems.

Tracing the origin of Lago Mare biota: Ostracods and mollusks from the late Neogene of the Slavonian mountains in the southern Pannonian Basin (NE Croatia)

Lake Pannon was a huge central European, long-lived endorheic lake occupying the Pannonian Basin System during the late Neogene. The benthic fauna of this originally sparsely populated brackish lake went through a spectacular adaptive radiation leading to a great number of autochthonous species. Although largely endemic, some of these species are excellent stratigraphic markers of the Lago Mare interval in the Mediterranean Basin. This situation apparently reflects a short-termed unidirectional migration-wave in the latest Miocene. The latter interval marks the final phase of a major environmental perturbation followed by a massive evaporite deposition termed the Messinian Salinity Crisis. The Section Bozara, located in the southern Pannonian Basin at the southern slopes of Mt. Papuk, bears a well-preserved benthic fauna representative of Lake Pannon. It is embedded into a 25-m-thick succession of alternating mud and sand packages. We detected therein 25 ostracod and 18 mollusk taxa enabling, along with a sedimentary facies analysis, an integrated evaluation of the depositional setting, biostratigraphic position, and paleogeographic distribution pattern. Accordingly, the depositional environment shows a shallowing upward trend from deep-water sublittoral and distal prodelta settings to deltaic high-energy littoral conditions. The biostratigraphic markers, such as the bivalve Rhombocongeria rhomboidea (Hörnes, 1867) and the ostracod Caspiocypris pontica (Sokač, 1972), constrained its stratigraphic position to the Portaferrian interval (8.0–4.5 Ma). Of 16 ostracod taxa determined at the species level, 10 are shared with the Eastern Paratethys, whereas only 3 are present in the Mediterranean Lago Mare deposits as well. In contrast, among 12 corresponding mollusk taxa, only 4 are known from the Eastern Paratethys, all being absent in the Lago Mare interval. Such a paleobiogeographic pattern suggests that the time window of the outflow event from Lake Pannon into the Eastern Paratethys was not synchronous with, but must have preceded, the hypothetic Eastern Paratethys drainage event generating the Lago Mare depositional conditions in the Mediterranean Basin.

The Subglacial Lake That Wasn't There: Improved Interpretation From Seismic Data Reveals a Sediment Bedform at Isunnguata Sermia

AbstractRadio Echo Sounding (RES) surveys conducted in May 2010 and April 2011 revealed a 2 km2 flat area with increased bed reflectivity at the base of Isunnguata Sermia at the western margin of the Greenland Ice Sheet. This flat reflector was located within a localized subglacial hydraulic potential (hydropotential) minimum, as part of a complex and elongated trough system. By analogy with comparable features in Antarctica, the initial interpretation of such a feature was a potential subglacial lake. In September 2013 a co‐located seismic survey revealed a 1,750 m by 540 and 37 m thick stratified lens‐shaped bedform at the base of a subglacial trough system. Amplitude Versus Angle (AVA) analysis yields a derived reflection coefficient R = 0.09 ± 0.14 indicative of consolidated sediments possibly overlain by dilatant till. The bed and flank on the northern side of the trough consist of unconsolidated, possibly water‐bearing sediments with R = −0.10 ± 0.08, whereas on the southern side it consists of more consolidated material. We interpret the trough as a key component of the wider subglacial drainage network, for which the sediments on its northern side act as a localized water‐storage reservoir. Given the observation of seasonally forming and rapidly draining supraglacial meltwater lakes in this area, we interpret the lens‐shaped bedform as deposited by episodically ponding meltwater within the subglacial trough system. Our results highlight the importance of transient subglacial hydrological and sedimentological processes such as drainage events for the interaction of ice sheets and their substrates, to understand ice dynamics in a warming climate.

Monitoring Thermokarst Lake Drainage Dynamics in Northeast Siberian Coastal Tundra

Thermokarst lakes in permafrost regions are highly dynamic due to drainage events triggered by climate warming. This study focused on mapping lake drainage events across the Northeast Siberian coastal tundra from 2000 to 2020 and identifying influential factors. An object-based lake analysis method was developed to detect 238 drained lakes using a well-established surface water dynamics product. The LandTrendr change detection algorithm, combined with continuous Landsat satellite imagery, precisely dated lake drainage years with 83.2% accuracy validated against manual interpretation. Spatial analysis revealed the clustering of drained lakes along rivers and in subsidence-prone Yedoma regions. The statistical analysis showed significant warming aligned with broader trends but no evident temporal pattern in lake drainage events. Our machine learning model identified lake area, soil temperature, summer evaporation, and summer precipitation as the top predictors of lake drainage. As these climatic parameters increase or surpass specific thresholds, the likelihood of lake drainage notably increases. Overall, this study enhanced the understanding of thermokarst lake drainage patterns and environmental controls in vulnerable permafrost regions. Spatial and temporal dynamics of lake drainage events were governed by complex climatic, topographic, and permafrost interactions. Integrating remote sensing with field studies and modeling will help project lake stability and greenhouse gas emissions under climate change.

Drainage Mechanisms of Submerged Sinkholes

Subsurface salt layer dissolution along the western shores of the Dead Sea is considered to be the primary cause for extensive large sinkhole formation in the past 40 years. Many of these sinkholes are arranged in clusters and are filled with water from nearby springs. The Mineral Beach resort was built in an area with a thermal spring with water emerging at around 40 °C at the Shalem sinkhole cluster. Unfortunately, the same spring was responsible for the destruction of the resort as it supplied water undersaturated with respect to halite, which promoted dissolution and sinkhole formation. The sinkholes in the Shalem cluster drain out in sudden catastrophic events and then slowly fill again. The drainage mechanisms of this phenomenon are studied in the Shalem-2 sinkhole cluster using leveling data collectors and ground-based LiDAR surveys over a period of 5 years, including thirty-five drainage events. Drainage volume and fluxes calculated using water level and topographic data obtained by LiDAR scans suggest that the formation of additional sinkholes beneath the pond’s bottom triggers drainage events. The subsequent flux shows that the evolution of the newly formed sinkholes either improves the hydraulic connection or temporarily seals the connection between the surface pond and deeper caverns/aquifers. The drainage event ends when either the hydraulic connection is sealed or when the level of water in the pond drops to the level of the newly formed sinkhole. The large volumes of drained water and drainage fluxes imply the existence of a well-developed active underground draining system.

Ice-dammed lakes and deglaciation history of the Scandinavian Ice Sheet in central Jämtland, Sweden

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.

Assessment of the interactions between soil–biosphere–atmosphere (ISBA) land surface model soil hydrology, using four closed-form soil water relationships and several lysimeters

Abstract. Soil water drainage is the main source of groundwater recharge and river flow. It is therefore a key process for water resource management. In this study, we evaluate the soil hydrology and the soil water drainage, simulated by the interactions between soil–biosphere–atmosphere (ISBA) land surface model currently used for hydrological applications from the watershed scale to the global scale, where parameters are generally not calibrated. This evaluation is done using seven lysimeters from two long-term model approach sites measuring hourly water dynamics between 2009 and 2019 in northeastern France. These 2 m depth lysimeters are filled with different soil types and are either maintained as bare soil or covered with vegetation. Four closed-form equations describing soil water retention and hydraulic conductivity functions are tested, namely the commonly used equations from Brooks and Corey (1966) and van Genuchten (1980), a combination of the van Genuchten (1980) soil water retention function with the Brooks and Corey (1966) unsaturated hydraulic conductivity function, and, for the very first time in a land surface model (LSM), a modified version of the van Genuchten (1980) equations, with a new hydraulic conductivity curve proposed by Iden et al. (2015). The results indicate good performance by ISBA with the different closure equations in terms of soil volumetric water content and water mass. The drained flow at the bottom of the lysimeter is well simulated, using Brooks and Corey (1966), while some weaknesses appear with van Genuchten (1980) due to the abrupt shape near the saturation of its hydraulic conductivity function. The mixed form or the new van Genuchten (1980) hydraulic conductivity function from Iden et al. (2015) allows the solving of this problem and even improves the simulation of the drainage dynamic, especially for intense drainage events. The study also highlights the importance of the vertical heterogeneity of the soil hydrodynamic parameters to correctly simulate the drainage dynamic, in addition to the primary influence of the parameters characterizing the shape of the soil water retention function.

The morphostratigraphic imprint of the Baltic Ice Lake drainage event in southern Finland

Digital elevation models, based on laser scanning imageries (LiDAR-DEM) and aided by ground penetrating radar (GPR) data, were used to study glaciofluvial Gilbert-type ice-contact deltas in the Younger Dryas Salpausselkä end-moraine zone in southern Finland. The geomorphological data analysed were used to reconstruct the water-level drop of the late glacial Baltic Ice Lake to the early Holocene Yoldia Sea and tie these changes to a wider stratigraphic context. The results indicate that the sudden drainage event at around 11 650 cal. yrs BP left its imprint not only on the varved sediments but also on ice-contact glaciofluvial deltas in the Second Salpausselkä zone throughout southern Finland. This morphostratigraphic boundary can be placed at locations where the ice-contact deltas occur at two different levels; the higher-level deltas formed during the Baltic Ice Lake B III water-level stage and the lower-level deltas during the Yoldia Sea Y I water-level stage. This morphostratigraphically defined boundary in southern Finland marks the Pleistocene/Holocene chronostratigraphic boundary in southern Finland and shows the corresponding positions of the Scandinavian Ice Sheet’s Finnish Lake District Ice Lobe and the Baltic Sea Ice Lobe.

Episodic Subglacial Drainage Cascades Below the Northeast Greenland Ice Stream

AbstractSubglacial hydrology can exert an important control on ice flow by affecting friction at the ice‐bedrock interface. Here, we report on a series of subglacial drainage events along the Northeast Greenland Ice Stream (NEGIS), initiating as far inland as 500 km from the margin of Zachariae Isstrøm. The drainage events exhibit local transient uplift, followed by prolonged subsidence, measured by differential satellite synthetic aperture radar interferometry (DInSAR). In downstream regions, drainage events are associated with temporary acceleration in ice flow. The high spatiotemporal resolution of the DInSAR measurements allows for a detailed mapping of the drainage propagation pathway. We show that multiple drainage cascades have occurred along the same pathway over the years 2020–2022. Finally, the propagation speed of subglacial water flow is found to vary greatly along NEGIS, suggesting that fundamental differences could exist in the subglacial environment.

Post-Little Ice Age glacial lake evolution in Svalbard: inventory of lake changes and lake types

AbstractThe rapid formation of glacial lakes is one of the most conspicuous landscape changes caused by atmospheric warming in glacierised regions. However, relatively little is known about the history and current state of glacial lakes in the High Arctic. This study aims to address this issue by providing the first inventory of glacial lakes in Svalbard, focusing in particular on the post-Little Ice Age evolution of glacial lakes and their typology. To do so, we used aerial photographs and satellite imagery together with archival topographic data from 1936 to 2020. The inventory comprises the development of 566 glacial lakes (146 km2) that were still in direct contact with glaciers during the period 2008–2012. The results show a consistent increase in the total area of glacial lakes from the 1930s to 2020 and suggest an apparent link between climatic and geological factors, and the formation of specific lake dam types: moraine, ice, or bedrock. We also detected 134 glacial lake drainage events that have occurred since the 1930s. This study shows that Svalbard has one of the highest rates of glacial lake development in the world, which is an indicator of the overall dynamics of landscape change in the archipelago in response to climate change.