A GIS-based multi-proxy analysis of the evolution of subglacial dynamics of the Quebec–Labrador ice dome,northeastern Quebec,Canada

The glacial records of the inner-core regions of the Laurentide Ice Sheet (LIS) document complex yet coherent patterns reflecting ice-sheet change (e.g. ice-divide migration), providing unique insights into past glacial conditions. This study develops a conceptual model of subglacial dynamics evolution within a major ice-dispersal centre of the LIS in northeastern Quebec, Canada using a GIS-based analysis of the surficial geologic record. Multiple proxies of subglacial conditions (subglacial streamlined landforms, lake density and lake area over thin drift/bedrock) were analysed through grid-overlay techniques and then classified based on different proxy variables ranging from highly mobile warm-based to immobile cold-based conditions. An additional proxy (till blanket) was used to identify areas of thick till deposition, but with few proxies (few lake or landform metrics). Based on local ice-flow reconstructions, the most ‘relict’ glacial terrain zone (GTZ1) has warm-based conditions over 66% of its area and is remarkably well preserved, suggesting laterally extensive warm-based conditions during the oldest identified ice-flow phase. This relict glacial terrain is partially overprinted by two subsequent ice-flow phases in spatially restricted zones in the northeast (73% warm-based), east-central (41% warm-based), and northwest (33% warm-based) of the study area. A zone of more sluggish conditions (only 3% warm-based) was identified in the highlands at the centre of the study area, characterized by thin till cover, few landforms, yet with large patches of relatively abundant small lakes, indicative of areal scouring. No clear evidence of sustained cold-based conditions (i.e. high chemical index of alteration values or high ¹⁰Be abundances) was found in the study area. These results suggest that warm-based conditions (active erosion and/or deposition) were uniformly widespread during the earliest ice-flow phase, later becoming more spatially restricted with broader sluggish ice conditions. These spatially restricted regions of warm-based subglacial regimes were likely controlled by surrounding and down-flow ice streaming. © 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd     

VARIABILITY IN THE CHEMICAL COMPOSITION OF IN SITU SUBGLACIAL MELTWATERS

Meltwaters collected from boreholes drilled to the base of the Haut Glacier d'Arolla, Switzerland have chemical compositions that can be classified into three main groups. The first group is dilute, whereas the second group is similar to, though generally less concentrated in major ions, than contemporaneous bulk glacial runoff. The third group is more concentrated than any observed bulk runoff, including periods of flow recession. Waters of the first group are believed to represent supraglacial meltwater and ice melted during drilling. Limited solutes may be derived from interactions with debris in the borehole. The spatial pattern of borehole water levels and borehole water column stratification, combined with the chemical composition of the different groups, suggest that the second group represent samples of subglacial waters that exchange with channel water on a diurnal basis, and that the third group represent samples of water draining through a ‘distributed’ subglacial hydraulic system. High NO⁻₃ concentrations in the third group suggest that snowmelt may provide a significant proportion of the waters and that the residence time of the waters at the bed in this particular section of the distributed system is of the order of a few months. The high NO⁻₃ concentrations also suggest that some snowmelt is routed along different subglacial flowpaths to those used by icemelt. The average SO²⁻₄: (HCO⁻₃ + SO²⁻₄) ratio of the third group of meltwaters is 0.3, suggesting that sulphide oxidation and carbonate dissolution (which gives rise to a ratio of 0.5) cannot provide all the HCO⁻₃ to solution. Hence, carbonate hydrolysis may be occurring before sulphide oxidation, or there may be subglacial sources of CO₂, perhaps arising from microbial oxidation of organic C in bedrock, air bubbles in glacier ice or pockets of air trapped in subglacial cavities. The channel marginal zone is identified as an area that may influence the composition of bulk meltwater during periods of recession flow and low diurnal discharge regimes. © 1997 by John Wiley & Sons, Ltd.     

VELOCITY–DISCHARGE RELATIONSHIPS DERIVED FROM DYE TRACER EXPERIMENTS IN GLACIAL MELTWATERS: IMPLICATIONS FOR SUBGLACIAL FLOW CONDITIONS

Repeated dye tracer tests were undertaken from individual moulins at Haut Glacier d'Arolla, Switzerland, over a number of diurnal discharge cycles during the summers of 1989–1991. It was hoped to use the concepts of at-a-station hydraulic geometry to infer flow conditions in subglacial channels from the form of the velocity–discharge relationships derived from these tests. The results, however, displayed both clockwise and anticlockwise velocity–discharge hysteresis, in addition to the simple power function relationship assumed in the hydraulic geometry approach. Clockwise hysteresis seems to indicate that a moulin drains into a small tributary channel rather than directly into an arterial channel, and that discharges in the two channels vary out of phase with each other. Anticlockwise hysteresis is accompanied by strong diurnal variations in the value of dispersivity derived from the dye breakthrough curve, and is best explained by hydraulic damming of moulins or sub/englacial passageways. Despite the complex velocity–discharge relationships observed, some indication of subglacial flow conditions may be obtained if tributary channels comprise only a small fraction of the drainage path and power function velocity–discharge relationships are derived from dye injections conducted during periods when the supraglacial discharge entering the moulin and the bulk discharge vary in phase. Analyses based on this premise suggest that both open and closed channel flow occur beneath Haut Glacier d'Arolla, and that flow conditions are highly variable at and between sites.     

Subglacial processes from drumlins in Clew Bay,western Ireland

Interpretation of subglacial processes and environments can be usefully informed by the stratigraphy and structures of sediments preserved within drumlins, and the down‐ice variability of these sediments. Drumlins in Clew Bay, western Ireland, were formed by westward late Pleistocene ice flow onto the Atlantic shelf, but the depositional processes and environments of these sediments remain uncertain. This study describes and interprets the drumlin stratigraphy and sediment properties and structures from two drumlins on the south side of Clew Bay. Drumlin sediments are dominated by massive to stratified diamicton deposited subglacially by deformation of flows of varying rheology. Folds, shears and clastic dikes within the sediments suggest active subglacial water flow related to variations in ice flow and ice‐bed coupling. The distal ends of the drumlins are dominated by stratified diamicton and gravel units with sandy interbeds. These sediments reflect deposition in a leeside subglacial cavity formed coevally with bedform development. A model is proposed that can account for sediment stratigraphic patterns and drumlin formation. Copyright © 2015 John Wiley & Sons, Ltd.     

Time‐transgressive tunnel‐valley infill revealed by a three‐dimensional sedimentary model,Hamburg, north‐west Germany

Deep, elongated incisions, often referred to as tunnel valleys, are among the most characteristic landforms of formerly glaciated terrains. It is commonly thought that tunnel valleys were formed by meltwater flowing underneath large ice sheets. The sedimentary infill of these features is often highly intricate and therefore difficult to predict. This study intends to improve the comprehension of the sedimentology and to establish a conceptual model of tunnel‐valley infill, which can be used as a predictive tool. To this end, the densely sampled, Pleistocene tunnel valleys in Hamburg (north‐west Germany) were investigated using a dataset of 1057 deep wells containing lithological and geophysical data. The stratigraphic correlations and the resulting three‐dimensional lithological model were used to assess the spatial lithological distributions and sedimentary architecture. The sedimentary succession filling the Hamburg area tunnel valleys can be subdivided into three distinct units, which are distinguished by their inferred depositional proximity to the ice margin. The overall trend of the succession shows a progressive decrease in transport energy and glacial influence through time. The rate of glacial recession appears to have been an important control on the sedimentary architecture of the tunnel‐valley fill. During periods of stagnation, thick ice‐proximal deposits accumulated at the ice margin, while during rapid recession, only a thin veneer of such coarse‐grained sediments was deposited. Ice‐distal and non‐glaciogenic deposits (i.e. lacustrine, marine and terrestrial) fill the remaining part of the incision. The infill architecture suggests formation and subsequent infill of the tunnel valleys at the outer margin of the Elsterian ice sheet during its punctuated northwards recession. The proposed model shows how the history of ice‐sheet recession determines the position of coarse‐grained depocentres, while the post‐glacial history controls the deposition of fines through a progressive infill of remnant depressions.     

Surface melt-driven seasonal behaviour (englacial and subglacial) from a soft-bedded temperate glacier recorded by in situ wireless probes

We investigate the spatial and temporal englacial and subglacial processes associated with a temperate glacier resting on a deformable bed using the unique Glacsweb wireless in situ probes (embedded in the ice and the till) combined with other techniques [including ground penetrating radar (GPR) and borehole analysis]. During the melt season (spring, summer and autumn), high surface melt leads to high water pressures in the englacial and subglacial environment. Winter is characterized by no surface melting on most days (‘base’) apart from a series of positive degree days. Once winter begins, a diurnal water pressure cycle is established in the ice and at the ice/sediment interface, with direct meltwater inputs from the positive degree days and a secondary slower englacial pathway with a five day lag. This direct surface melt also drives water pressure changes in the till. Till deformation occurred throughout the year, with the winter rate approximately 60% that of the melt season. We were able to show the bed comprised patches of till with different strengths, and were able to estimate their size, relative percentage and temporal stability. We show that the melt season is characterized by a high pressure distributed system, and winter by a low pressure channelized system. We contrast this with studies from Greenland (overlying rigid bedrock), where the opposite was found. We argue our results are typical of soft bedded glaciers with low englacial water content, and suggest this type of glacier can rapidly respond to surface-driven melt. Based on theoretical and field results we suggest that the subglacial hydrology comprises a melt season distributed system dominated by wide anastomosing broad flat channels and thin water sheets, which may become more channelized in winter, and more responsive to changes in meltwater inputs. © 2019 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

HYDROCHEMISTRY AS AN INDICATOR OF SUBGLACIAL DRAINAGE SYSTEM STRUCTURE: A COMPARISON OF ALPINE AND SUB-POLAR ENVIRONMENTS

The anion compositions (SO²₄, HCO⁻₃ and Cl⁻) of runoff from the Haut Glacier d'Arolla, Switzerland and Austre Brøggerbreen, Svalbard are compared to assess whether or not variations in water chemistry with discharge are consistent with current understanding of the subglacial drainage structure of warm- and polythermal-based glaciers. These glacial catchments have very different bedrocks and the subglacial drainage structures are also believed to be different, yet the range of anion concentrations show considerable overlap for SO²⁻₄ and HCO⁻₃. Concentrations of Cl⁻ are higher at Austre Brøggerbreen because of the maritime location of the glacier. Correcting SO²⁻₄ for the snowpack component reveals that the variation in non-snowpack SO²⁻₄ with discharge and with HCO⁻₃ is similar to that observed at the Haut Glacier d'Arolla. Hence, if we assume that the provenance of the non-snowpack SO²⁻₄ is the same in both glacial drainage systems, a distributed drainage system also contributes to runoff at Austre Brøggerbreen. We have no independent means of testing the assumption at present. The lower concentrations of non-snowpack SO²⁻₄ at Austre Brøggerbreen may suggest that a smaller proportion of runoff originates from a distributed drainage system than at the Haut Glacier d'Arolla.     

OUTBURST FLOODS FROM GLACIER-DAMMED LAKES: THE EFFECT OF MODE OF LAKE DRAINAGE ON FLOOD MAGNITUDE

Published accounts of outburst floods from glacier-dammed lakes show that a significant number of such floods are associated not with drainage through a tunnel incised into the basal ice—the process generally assumed—but rather with ice-marginal drainage, mechanical failure of part of the ice dam, or both. Non-tunnel floods are strongly correlated with formation of an ice dam by a glacier advancing from a tributary drainage into either a main river valley or a pre-existing body of water (lake or fiord). For a given lake volume, non-tunnel floods tend to have significantly higher peak discharges than tunnel-drainage floods. Statistical analysis of data for floods associated with subglacial tunnels yields the following empirical relation between lake volume V and peak discharge Q_p : Q_p = 46V^0.66 (r² = 0.70), when Q_p is expressed in metres per second and V in millions of cubic metres. This updates the so-called Clague–Mathews relation. For non-tunnel floods, the analogous relation is Q_p = 1100V^0.44 (r² = 0.58). The latter relation is close to one found by Costa (1988) for failure of constructed earthen dams. This closeness is probably not coincidental but rather reflects similarities in modes of dam failure and lake drainage. We develop a simple physical model of the breach-widening process for non-tunnel floods, assuming that (1) the rate of breach widening is controlled by melting of the ice, (2) outflow from the lake is regulated by the hydraulic condition of critical flow where water enters the breach, and (3) the effect of lake temperature may be dealt with as done by Clarke (1982). Calculations based on the model simulate quite well outbursts from Lake George, Alaska. Dimensional analysis leads to two approximations of the form Q_p ∝ V^qf(h_i, θ₀), where q = 0.5 to 0.6, h_i is initial lake depth, θ₀ is lake temperature, and the form of f (h_i, θ₀) depends on the relative importance of viscous dissipation and the lake's thermal energy in determining the rate of breach opening. These expressions, along with the regression relations, should prove useful for assessing the probable magnitude of breach-type outburst floods.     

青藏高原东南部贡嘎山海螺沟冰川冰下沉积物形成与变形

冰岩界面的冰川动力学是冰川系统的重要组成部分, 海螺沟冰川地处温暖湿润的海洋环境, 冰川运动速度较快, 冰川底部接近压融点, 是研究冰下过程的较理想地点. 在海螺沟冰川大型磨光面上浅显侵蚀坑内发现了碎屑物质. 对碎屑物质理化特征研究表明: 粒度特征、地球化学与石英砂SEM 分析表明沉积在冰岩界面上的物质来自于冰川底部的底碛层, 而不是冰上环境的产物. 偏光显微镜下观察到的冰下沉积物呈现出一系列塑性变形(微旋转、褶皱)和脆性变形(线性结构、支撑结构、断层)微观结构和构造. 两种变形结构的存在是碎屑物质在形成过程中其含水量波动情况的反映. 冰下碎屑物质是冰下融出、滞碛作用的共同产物. 在整个冰下碎屑物质形成与变形过程中, 由于冰下水系季节性变化带来的冰岩界面上冰川融水含量的波动起了决定性作用.