Subglacial Drainage Network Research Articles

Isotope Hydrograph Separation Reveals Rainfall on the Glaciers Will Enhance Ice Meltwater Discharge to the Himalayan Rivers

AbstractThe Indian Summer Monsoon (ISM) and meltwater from the Himalayan are the two most important sources of water in the Indian subcontinent. However, the impact of ISM on Himalayan glaciers and subsequent stream hydrology remains largely unknown. To provide new insight into the impact of rainfall on glacial hydrology, here we present hydro‐meteorological and time‐series observations of meltwater stable water isotope compositions from the snout of the Chorabari glacier in the Upper Ganga Basin, Central Himalayas across the ablation season corresponding to 2019. We observe that rainfall events (>2 mm d−1) on the glacier enhance discharge driven by ice meltwater in River Mandakini. Energy balance calculations reveal that one of the drivers behind enhanced ice meltwater contribution could be rain‐induced melting of the glacier where rainfall on the ice surface melts the glacier producing up to 13% of the total discharge at the glacier snout. Further, rainfall on glacier surface have other control on glacial processes—for example, snow metamorphism, ice flow dynamics such as short‐term acceleration in ice speed flow, and reorganization of the englacial and subglacial drainage network—that are poorly studied and needs further investigation. We conclude rainfall events on the glacier have a complex control on mountain hydrology. This study, therefore, provides an interpretative framework that calls for additional assessments of the direct and indirect impact of rainfall in glacial hydrology.

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.

Trace element, rare earth element and trace carbon compounds in Subglacial Lake Whillans, West Antarctica

Whillans Subglacial Lake (SLW) lies beneath 801 m of ice in the lower portion of the Whillans Ice Stream (WIS) in West Antarctica and is part of an extensive and active subglacial drainage network. Here, the geochemical characterization of SLW rare earth elements (REE), trace elements (TE), free amino acids (FAA), and phenolic compounds (PC) measured in lakewater and sediment porewater are reported. The results show, on average, higher values of REEs in the lakewater than in the porewater, and clear changes in all REE concentrations and select redox sensitive trace element concentrations in porewaters at a depth of ~15 cm in the 38 cm lake sediment core. This is consistent with prior results on the lake sediment redox conditions based on gas chemistry and microbiological data. Low concentrations of vanillyl phenols were measured in the SLW water column with higher concentrations in porewater samples and their concentration profiles in the sediments may also reflect changing redox conditions in the sediments. Vanillin concentrations increased with depth in the sediments as oxygenation decreases, while the concentrations of vanillic acid, the more oxidized component, were higher in the more oxygenated surface sediments. Collectively these results indicate redox changes occurring with the upper 38 cm of sediment in SLW and provide support for the existence of a seawater source, already hypothesized, in the sediments below the lowest measured depth, and of a complex and dynamic geochemical system beneath the West Antarctic Ice Sheet.Our results are the first to detail geochemical properties from an Antarctic subglacial environment using direct sampling technology. Due to their isolation from the wider environment, subglacial lakes represent one of our planets last pristine environments that provide habitats for microbial life and natural biogeochemical cycles but also impact the basal hydrology and can cause ice flow variations.

Threshold response to melt drives large-scale bed weakening in Greenland.

Ice speeds in Greenland are largely set by basal motion1, which is modulated by meltwater delivery to the ice base2-4. Evidence suggests that increasing melt rates enhance the subglacial drainage network's capacity to evacuate basal water, increasing bed friction and causing the ice to slow5-10. This limits the potential of melt forcing to increase mass loss as temperatures increase11. Here we show that melt forcing has a pronounced influence on dynamics, but factors besides melt rates primarily control its impact. Using a method to examine friction variability across the entirety of western Greenland, we show that the main impact of melt forcing is an abrupt north-to-south change in bed strength that cannot be explained by changes in melt production. The southern ablation zone is weakened by 20-40 per cent compared with regions with no melt, whereas in northern Greenland the ablation zone is strengthened. We show that the weakening is consistent with persistent basal water storage and that the threshold is linked to differences in sliding and hydropotential gradients, which exert primary control on the pressures within drainage pathways that dewater the bed. These characteristics are mainly set by whether a margin is land or marine terminating, suggesting that dynamic changes that increase mass loss are likely to occur in northern Greenland as temperatures increase. Our results point to physical representations of these findings that will improve simulated ice-sheet evolution at centennial scales.

Subglacial evolution from distributed to channelized drainage: Evidence from the Lake Murtoo area in SW Finland

AbstractDuring the last millennia of the Weichselian glaciation (c. 11 700–14 000 cal BP) and the first millennium of the Holocene interglacial (c. 10 500–11 700 cal BP), melting of the ice sheet and subsequent subglacial drainage resulted in the formation of murtoos and murtoo‐related landforms in the Fennoscandian Ice Sheet area in Finland and Sweden. Murtoos are composed of clay‐poor diamictons produced by sediment‐concentrated creep and material sorting in a subglacial environment under effective pressure close to zero, and potentially represent a transition form from non‐channelized to channelized subglacial drainage networks. Here, we explore the geomorphology and substratum characteristics of the Lake Murtoo area, one of the locations where murtoos were initially found and described in SW Finland. Our interpretation suggests an interconnected development from inefficient distributed or semi‐distributed cavities (murtoos) to efficient channelized drainage (esker) via erosional escarpments, all displayed in the same region. In addition, our excavation into a murtoo‐related erosional escarpment (MRE) revealed subglacial erosional channels alongside MREs that were filled with fine‐grained finely laminated sediments. We suggest that these channel fill sediments were deposited subglacially, because they are covered by flow till and show spatial and lithological characteristics that deviate significantly from the proglacial basin fill clay (varve) sediments in the area. The appearance of rhythmically laminated silt and clay sediments indicates sedimentation in a suspension‐load‐dominated environment in shallow subglacial channels with relatively low (but repetitiously changing) water flow velocities, possibly reflecting disconnected spill‐over routes of subglacial drainage in small tributaries that were only periodically used due to water pressure changes in the up‐ice direction.

Subglacial meltwater routes of the Fennoscandian Ice Sheet

ABSTRACT Subglacial drainage systems are crucial elements of glaciers and ice sheets because they modulate ice flow velocity. However, logistical challenges of measuring subglacial processes beneath contemporary ice and natural limitations in long-term monitoring hinder our understanding about their spatio-temporal evolution. Subglacial meltwater landforms created by palaeo-ice sheets are records of past subglacial drainage systems and offer the potential to study their large-scale development throughout deglaciation. Although collectively recording subglacial drainage, individual meltwater landforms such as eskers, meltwater channels and meltwater corridors, which comprise tunnel valleys and meltwater tracks (assemblages of landforms in broad, elongated paths with irregular surface texture), have mostly been investigated as separate entities. Using high-resolution (1–2 m) digital elevation models, we map integrated networks of subglacial meltwater landforms, herein called subglacial meltwater routes, on an ice-sheet scale in Fennoscandia. Our map provides a basis for future research on the long-term evolution of subglacial drainage networks and its effect on ice dynamics of the Fennoscandian Ice Sheet.

Seasonality of Solute Flux and Water Source Chemistry in a Coastal Glacierized Watershed Undergoing Rapid Change: Wolverine Glacier Watershed, Alaska

AbstractAs glaciers around the world rapidly lose mass, the tight coupling between glaciers and downstream ecosystems is resulting in widespread impacts on global hydrologic and biogeochemical cycling. However, a range of challenges make it difficult to conduct research in glacierized systems, and our knowledge of seasonally changing hydrologic processes and solute sources and signatures is limited. This in turn hampers our ability to make predictions on solute composition and flux. We conducted a broad water sampling campaign in order to understand the present‐day partitioning of water sources and associated solutes in Alaska's Wolverine Glacier watershed. We established a relationship between electrical conductivity and streamflow at the watershed outlet to divide the melt season into four hydroclimatic periods. Across hydroclimatic periods, we observed a shift in nonglacial source waters from snowmelt‐dominated overland and shallow subsurface flow paths to deeper groundwater flow paths. We also observed the shift from a low‐ to high‐efficiency subglacial drainage network and the associated flushing of water stored subglacially with higher solute loads. We used calcium, the dominant dissolved ion, from watershed outlet samples to estimate solute fluxes for each hydroclimatic period across two melt seasons. We found between 40% and 55% of Ca2+ export occurred during the late season rainy period. This partitioning of the melt season coupled with a characterization of the chemical makeup and magnitude of solute export provides new insight into a rapidly changing watershed and creates a framework to quantify and predict changes to solute fluxes across a melt season.

Diversity of murtoos and murtoo‐related subglacial landforms in the Finnish area of the Fennoscandian Ice Sheet

Murtoos are recently discovered triangular‐shaped subglacial landforms that form under warm‐based ice and in association with significant subglacial meltwater flow. They appear in distinct fields and commonly occur in the area that was covered by the Fennoscandian Ice Sheet during glacial periods. Murtoos potentially represent a transition form from non‐channelized to channelized subglacial drainage networks. In the present study, we analyse and classify murtoos and murtoo‐related landforms in the Finnish area of the Fennoscandian Ice Sheet based on their characteristics and appearance in LiDAR‐based digital elevations models. Combined with morphometric analyses, the observations suggest that five types of murtoos and murtoo‐related landforms are common and widespread in Finland: (i) triangle‐type murtoos (TTMs), (ii) chevron‐type murtoos (CTMs), (iii) lobate‐type murtoos (LTMs), (iv) murtoo‐related ridges and escarpments (MREs), and (v) other murtoo‐related polymorphous landforms (PMRs) that look like small mounds and ridges. The morphometric characteristics of the different types are described here in detail, and it is shown that they are spatially and geomorphologically related. In addition, we provide examples of murtoos other than the TTMs to demonstrate that different murtoo types and murtoo‐related landforms are composed of similar sediments and architectural characteristics. The diversity of murtoo landforms and the transition between distinct murtoo types indicate rapid and complicated variations in the configuration of subglacial hydrology at different spatial and temporal scales. This study emphasizes the essential role of subglacial meltwater in the shaping of glacial landscapes and the redistribution of large volumes of sediments during the deglaciation of the Fennoscandian Ice Sheet.

Parallel tunnel channels of the Huron-Erie Lobe of the Laurentide Ice Sheet: A variation of tunnel channels with implications on lobe history and subglacial meltwater dynamics

Newly available high resolution elevation data across North America and Europe have revealed previously unknown glacial landforms. This study focuses on the morphology and stratigraphy of a variation of tunnel channels with unique morphology that we have named parallel tunnel channels (PTCs). They are found in northeastern Indiana, northwestern Ohio, and southeastern Michigan, which were previously glaciated by the Huron-Erie Lobe (HEL) of the Laurentide Ice Sheet (LIS). PTCs are discontinuous, quasi-parallel pairs of depressions separated by a continuous central ridge that can extend upwards of 50 km, cross-cut multiple recessional moraines, and transition into standard tunnel channels. The PTCs range in width from 200 to 900 m, and the depressions range in width from 100 to 200 m and in depth from 3 to 10 m. Using water well logs and multiple electrical resistivity surveys, we determined that the stratigraphy beneath and surrounding the PTCs consists of a lower layer of sand and gravel (confined aquifer) and an upper layer of New Holland till, both of the Lagro Formation. A vibracore and multiple resistivity surveys reveal that depressions are filled with ~4 m of lacustrine sediment and occasionally are floored with gravel channel deposits. The depressions are erosional in origin, rather than a product of ice block melt-out and are cut into the upper surface of the New Holland till. We interpret the PTCs as subglacial, erosional landforms, similar to tunnel channels, which represent distinct pathways transporting subglacial meltwater from the interior to the ice margin. Identification of these landforms has several important implications for understanding the deglaciation of the HEL: 1) the major recessional moraines may not record the final deglaciation of the HEL (i.e., are palimpsest features); 2) a large portion of the HEL must have stagnated to preserve the subglacial landscape including the PTCs; and 3) the combined network of PTCs, tunnel channels, and proglacial spillways may represent a partially complete subglacial drainage network with distinct morphological changes as meltwater was transported from the interior to the margin of the HEL. An initial formational model is proposed for the formation of the PTCs by way of a supra- and/or subglacial lake draining through an ice-walled subglacial tunnel, but additional work is required to understand the distinct pattern of depressions observed.

Characterization of subglacial marginal channels using 3-D analysis of high-density ground-penetrating radar data

AbstractStudying subglacial drainage networks is important for understanding the potential relationship between channel dynamics and rapid glacier recession as well as the role of subglacial channels in subglacial sediment evacuation. In order to delineate the planform geometry of snout marginal subglacial channels, densely spaced ground-penetrating radar (GPR) measurements at a frequency of ~70 MHz were carried out over the snout marginal zones of two temperate glaciers in the southwestern Swiss Alps, the Haut Glacier d'Arolla and the Glacier d'Otemma. Three-dimensional (3-D) data processing and amplitude analysis of the GPR reflection along the glacier bed was used to map the channels. At the Haut Glacier d'Arolla, two relatively straight channels of several meters in width were identified. The positions of these channels correspond well with the locations of channel outlets at the glacier terminus, as well as with fractures appearing on the glacier surface one month after the GPR data acquisition. The latter are believed to represent the beginning of ice collapse above the subglacial channels. At the Glacier d'Otemma, a major subglacial conduit was detected with similar dimensions to those identified at the Haut Glacier d'Arolla, but greater sinuosity. The position of this channel was confirmed by drone-based imagery acquired after glacier margin collapse. Our results confirm that high-density 3-D GPR surveys can be used to map subglacial channels near temperate alpine glacier margins.

Rapid mineral surface weathering beneath the Greenland Ice Sheet shown by radium and uranium isotopes

Glaciated watersheds are regions of intense physical and chemical weathering. In order to gain new insight on subglacial weathering processes, we measured uranium and radium isotopes from a proglacial river draining the Greenland Ice Sheet (GrIS). Time series samples were collected from the spring to mid-summer, a time period during which subglacial drainage pathways are thought to transition from slow-inefficient (distributed) to fast-efficient (channelized) systems. The 228Ra/226Ra activity ratio of the dissolved load varied from 5.2 ± 0.9–16.9 ± 3.6, which was significantly higher than the 228Ra/226Ra ratio of a suspended sediment load sample of 2.1 ± 0.07 and crustal values of ~1. The high 228Ra/226Ra in the dissolved load relative to the source rock material is indicative of mineral surface weathering induced by rapid and continuous flushing of the subglacial drainage network during the course of the melt season and those prior. The 234U/238U ratio (δ234U) varied between 33 and 106‰ with a discharge-weighted mean of 67‰; the seasonal evolution of δ234U did not correlate to geochemical indicators of subglacial meltwater storage time. An experiment designed to measure changes in δ234U with increasing meltwater storage times found that δ234U in the dissolved phase decreased rapidly with increasing storage time. Similarly, samples collected along a transect moving downstream from the ice sheet terminus had decreasing δ234U values from 63 to 15‰ further indicating that with increased weathering, the δ234U of meltwater decreases. Coupled with the relatively low δ234U and high 228Ra/226Ra, U appears to be impacted by rapid chemical weathering of subglacial and suspended sediments. The Leverett River discharge weighted U concentration was 0.13 nM; if this system is considered representative of the broader GrIS, then the total dissolved U flux from the GrIS would be on the order of 6.4 × 104 mol/y. Using a similar set of assumptions, the dissolved 228Ra and 226Ra flux from the GrIS was ~1.1 × 1014 dpm/y and ~ 5.5 × 1013 dpm/y, respectively. These estimates suggest that the 226Ra flux to the ocean from the GrIS is globally significant and that the 228Ra flux in particular is larger than most river inputs.

A two-dimensional glacier–fjord coupled model applied to estimate submarine melt rates and front position changes of Hansbreen, Svalbard

ABSTRACTWe have developed a two-dimensional coupled glacier–fjord model, which runs automatically using Elmer/Ice and MITgcm software packages, to investigate the magnitude of submarine melting along a vertical glacier front and its potential influence on glacier calving and front position changes. We apply this model to simulate the Hansbreen glacier–Hansbukta proglacial–fjord system, Southwestern Svalbard, during the summer of 2010. The limited size of this system allows us to resolve some of the small-scale processes occurring at the ice–ocean interface in the fjord model, using a 0.5 s time step and a 1 m grid resolution near the glacier front. We use a rich set of field data spanning the period April–August 2010 to constrain, calibrate and validate the model. We adjust circulation patterns in the fjord by tuning subglacial discharge inputs that best match observed temperature while maintaining a compromise with observed salinity, suggesting a convectively driven circulation in Hansbukta. The results of our model simulations suggest that both submarine melting and crevasse hydrofracturing exert important controls on seasonal frontal ablation, with submarine melting alone not being sufficient for reproducing the observed patterns of seasonal retreat. Both submarine melt and calving rates accumulated along the entire simulation period are of the same order of magnitude, ~100 m. The model results also indicate that changes in submarine melting lag meltwater production by 4–5 weeks, which suggests that it may take up to a month for meltwater to traverse the englacial and subglacial drainage network.

Proglacial sediment dynamics from daily to seasonal scales in a glaciated Alpine catchment (Bossons glacier, Mont Blanc massif, France)

AbstractThe sediment yields of Alpine catchments are commonly determined from streamload measurements made some distance downstream from glaciers. However, this approach indiscriminately integrates erosion processes occurring in both the glacial and proglacial areas. A specific method is required to ascertain the respective inputs from (i) subglacial and supraglacial sediments, (ii) proglacial hillslopes and (iii) proglacial alluvial areas or sandurs. This issue is addressed here by combining high‐resolution monitoring (2 min) of suspended sediment concentrations at different locations within a catchment with discharge gauging and precipitation data. This methodological framework is applied to two proglacial streams draining the Bossons glacier (Mont Blanc massif, France): the Bossons and Crosette streams. For the Bossons stream, discharge and suspended load data were acquired from June to October 2013 at 1.15 and 1.5 km from the glacial terminus, respectively upstream and downstream from a small valley sandur. These hydro‐sedimentary data are compared with the Crosette stream dataset acquired at the outlet of the Bossons glacier subglacial drainage system. A fourfold analysis focusing on seasonal changes in streamload and discharge, multilinear regression modelling, evaluation of the sandur flux balance and probabilistic uncertainty assessment is used to determine the catchment sediment budget and to explain the proglacial sediment dynamics. The seasonal fluctuation of the sediment signal observed is related to the gradual closing of the subglacial drainage network and to the role of the proglacial area in the sediment cascade: the proglacial hillslopes appear to be disconnected from the main channel and the valley sandur acts as a hydrodynamic sediment buffer both daily and seasonally. Our findings show that an understanding of proglacial sediment dynamics can help in evaluating paraglacial adjustment and subglacial erosion processes. Copyright © 2017 John Wiley & Sons, Ltd.

Bedrock morphology reveals drainage network in northeast Baffin Bay

A subglacial drainage network underneath the paleo-ice sheet off West Greenland is revealed by a new compilation of high-resolution bathymetry data from Melville Bay, northeast Baffin Bay. This drainage network is an indicator for ice streaming and subglacial meltwater flow toward the outer shelf. Repeated ice sheet advances and retreats across the crystalline basement together with subglacial meltwater drainage had their impact in eroding overdeepened troughs along ice stream pathways. These overdeepenings indicate the location of a former ice sheet margin. The troughs inherit characteristics of glacial and subglacial meltwater erosion. Most of the troughs follow tectonic weakness zones such as faults and fractures in the crystalline bedrock. Many of these tectonic features correspond with the orientations of major fault axes in the Baffin Bay region. The troughs extend from the present (sub) glacial fjord systems at the Greenland coast and parallel modern outlet-glacier pathways. The fast flowing paleo-ice streams were likely accelerated from the meltwater flow as indicated by glacial landforms within and along the troughs. The ice streams flowed along narrow tributary troughs and merged to form large paleo-ice streams bedded in the major cross-shelf troughs of Melville Bay. Apart from the troughs, a rough seabed topography characterises the bedrock, and we see a sharp geomorphic transition where ice flowed onto sedimentary rock and deposits.

Carbon dating reveals a seasonal progression in the source of particulate organic carbon exported from the Greenland Ice Sheet

AbstractSurface melt from the Greenland Ice Sheet (GrIS) collects particulate organic carbon (POC) as it drains into subglacial environments and transports it downstream where it serves as a microbial substrate. We hypothesized that older POC is entrained by meltwaters as the subglacial drainage network expands upglacier over the summer. To test this, POC samples were collected from a meltwater river exiting the GrIS over an ablation season and 14C dated. Resulting values were compared with meltwater hydrochemistry and satellite observations of the catchment area. We found that POC ages increased from ~5000 to ~9000 years B.P. until peak discharge and catchment size. Afterward, significant fluctuations in POC age were observed, interpreted to result from periods of high and low subglacial hydrological pressure and sediment supply and subsequent exhaustion. These observations suggest a seasonal progression in the source of POC exported from the GrIS and provide evidence for a seasonally evolving subglacial drainage system.

Deglaciation of the Eurasian ice sheet complex

The Eurasian ice sheet complex (EISC) was the third largest ice mass during the Last Glacial Maximum with a span of over 4500 km and responsible for around 20 m of eustatic sea-level lowering. Whilst recent terrestrial and marine empirical insights have improved understanding of the chronology, pattern and rates of retreat of this vast ice sheet, a concerted attempt to model the deglaciation of the EISC honouring these new constraints is conspicuously lacking. Here, we apply a first-order, thermomechanical ice sheet model, validated against a diverse suite of empirical data, to investigate the retreat of the EISC after 23 ka BP, directly extending the work of Patton et al. (2016) who modelled the build-up to its maximum extent. Retreat of the ice sheet complex was highly asynchronous, reflecting contrasting regional sensitivities to climate forcing, oceanic influence, and internal dynamics. Most rapid retreat was experienced across the Barents Sea sector after 17.8 ka BP when this marine-based ice sheet disintegrated at a rate of ∼670 gigatonnes per year (Gt a−1) through enhanced calving and interior dynamic thinning, driven by oceanic/atmospheric warming and exacerbated by eustatic sea-level rise. From 14.9 to 12.9 ka BP the EISC lost on average 750 Gt a−1, peaking at rates >3000 Gt a−1, roughly equally partitioned between surface melt and dynamic losses, and potentially contributing up to 2.5 m to global sea-level rise during Meltwater Pulse 1A. Independent glacio-isostatic modelling constrained by an extensive inventory of relative sea-level change corroborates our ice sheet loading history of the Barents Sea sector. Subglacial conditions were predominately temperate during deglaciation, with over 6000 subglacial lakes predicted along with an extensive subglacial drainage network. Moreover, the maximum EISC and its isostatic footprint had a profound impact on the proglacial hydrological network, forming the Fleuve Manche mega-catchment which had an area of ∼2.5 × 106 km2 and drained the present day Vistula, Elbe, Rhine and Thames rivers through the Seine Estuary. During the Bølling/Allerød oscillation after c. 14.6 ka BP, two major proglacial lakes formed in the Baltic and White seas, buffering meltwater pulses from eastern Fennoscandia through to the Younger Dryas when these massive proglacial freshwater lakes flooded into the North Atlantic Ocean. Deglaciation temporarily abated during the Younger Dryas stadial at 12.9 ka BP, when remnant ice across Svalbard, Franz Josef Land, Novaya Zemlya, Fennoscandia and Scotland experienced a short-lived but dynamic re-advance. The final stage of deglaciation converged on present day ice cover around the Scandes mountains and the Barents Sea by 8.7 ka BP, although the phase-lagged isostatic recovery still continues today.

Active subglacial lakes and channelized water flow beneath the Kamb Ice Stream

Abstract. We identify two previously unknown subglacial lakes beneath the stagnated trunk of the Kamb Ice Stream (KIS). Rapid fill-drain hydrologic events over several months are inferred from surface height changes measured by CryoSat-2 altimetry and indicate that the lakes are probably connected by a subglacial drainage network, whose structure is inferred from the regional hydraulic potential and probably links the lakes. The sequential fill-drain behavior of the subglacial lakes and concurrent rapid thinning in a channel-like topographic feature near the grounding line implies that the subglacial water repeatedly flows from the region above the trunk to the KIS grounding line and out beneath the Ross Ice Shelf. Ice shelf elevation near the hypothesized outlet is observed to decrease slowly during the study period. Our finding supports a previously published conceptual model of the KIS shutdown stemming from a transition from distributed flow to well-drained channelized flow of subglacial water. However, a water-piracy hypothesis in which the KIS subglacial water system is being starved by drainage in adjacent ice streams is also supported by the fact that the degree of KIS trunk subglacial lake activity is relatively weaker than those of the upstream lakes.

Glacier naled evolution and relation to the subglacial drainage system based on water chemistry and GPR surveys (Werenskioldbreen, SW Svalbard)

ABSTRACTGlacier naledi are extrusive ice masses that appear in front of glaciers as a consequence of refreezing of meltwater seepage during the accumulation season. These structures provide a unique opportunity to understand subglacial drainage activity during the accumulation season; however, only few detailed studies have previously focused on their characteristics. Here, we investigated glacier-derived naled assemblages in the proglacial zone of the polythermal glacier Werenskioldbreen (27.4 km2) in SW Svalbard. We determined the spatial distribution of naledi using ground penetrating radar surveys. The main subglacial drainage pattern was related to a channel under the medial moraine, and three sources are linked to a distributed subglacial drainage network. The relation between atmospherically-corrected (Ca2+ + Mg2+) and (SO42−) in sub-naled waters was closely related to sulphide oxidation coupled with carbonate dissolution (r = 0.99; slope = 1.6). This is consistent with the local lithology, which is dominated by schist containing carbonates. We also found high carbonate saturation indices in pale white ice layers within the naled. We conclude that sulphide oxidation coupled with carbonate dissolution is the dominant chemical weathering process in the subglacial drainage system of Werenskioldbreen during the accumulation season.

Stable oxygen isotope variability in two contrasting glacier river catchments in Greenland

Abstract. Analysis of stable oxygen isotope (δ18O) characteristics is a useful tool to investigate water provenance in glacier river systems. In order to attain knowledge on the diversity of δ18O variations in Greenlandic rivers, we examined two contrasting glacierised catchments disconnected from the Greenland Ice Sheet (GrIS). At the Mittivakkat Gletscher river, a small river draining a local temperate glacier in southeast Greenland, diurnal oscillations in δ18O occurred with a 3 h time lag to the diurnal oscillations in run-off. The mean annual δ18O was −14.68 ± 0.18 ‰ during the peak flow period. A hydrograph separation analysis revealed that the ice melt component constituted 82 ± 5 % of the total run-off and dominated the observed variations during peak flow in August 2004. The snowmelt component peaked between 10:00 and 13:00 local time, reflecting the long travel time and an inefficient distributed subglacial drainage network in the upper part of the glacier. At the Kuannersuit Glacier river on the island Qeqertarsuaq in west Greenland, the δ18O characteristics were examined after the major 1995–1998 glacier surge event. The mean annual δ18O was −19.47 ± 0.55 ‰. Despite large spatial variations in the δ18O values of glacier ice on the newly formed glacier tongue, there were no diurnal oscillations in the bulk meltwater emanating from the glacier in the post-surge years. This is likely a consequence of a tortuous subglacial drainage system consisting of linked cavities, which formed during the surge event. Overall, a comparison of the δ18O compositions from glacial river water in Greenland shows distinct differences between water draining local glaciers and ice caps (between −23.0 and −13.7 ‰) and the GrIS (between −29.9 and −23.2 ‰). This study demonstrates that water isotope analyses can be used to obtain important information on water sources and the subglacial drainage system structure that is highly desired for understanding glacier hydrology.

Implications for post-comminution processes in subglacial suspended sediment using coupled radiogenic strontium and neodymium isotopes

Enhanced physical weathering rates in subglacial systems promote high levels of comminution, transport, and deposition of fine-grained sediment within the subglacial drainage network. The impact of shifts in sediment loads from variations in meltwater flux, and their effects on downstream ecosystems, remains poorly quantified and places a fundamental importance on our ability to characterize subglacial depositional environments. Here, for the first time, we assess the seasonal evolution of the subglacial suspended sediment using coupled radiogenic strontium (87Sr/86Sr) and neodymium (143Nd/144Nd) isotopic ratios with elemental ratios and in situ measurements. Weathering rates in fluvial and riverine systems have been traditionally assessed using radiogenic isotopic tracers: 143Nd/144Nd ratios relate to the crustal age whereas 87Sr/86Sr ratios relate to age and preferential mineral dissolution. Thus relative shifts in these ratios will allow us to characterize distinct sediment transport networks. We apply this technique to the Lemon Creek Glacier (LCG), Alaska, USA, and to the Athabasca Glacier (AG), Alberta, CA. At the LCG, the 143Nd/144Nd values range from εNd of −4.6 (0.9) to −8.7 (0.2), which suggests a poorly mixed sediment flux. However, the greatest period of variability may correlate with the drainage of a supraglacial lake and suggests caution should be exerted in time-scale 143Nd/144Nd provenance studies that may be affected by climatic disturbances. In contrast, limited variation is observed within the AG 143Nd/144Nd seasonal record. A consistent, direct relation between the Rb/Sr elemental ratio and the 87Sr/86Sr ratio proves interesting as it enables us to unravel incongruent weathering trends in the radiogenic Sr record. Correlation between the 87Sr/86Sr and total discharge suggests that the process is partially controlled by mantling of the bedrock, which can be detected using post-comminution ages. While the subglacial structure may be enabled by the subglacial till beneath the AG, our study supports the use of Sr-Nd as a new proxy in the subglacial environment.