The historical Saklei Shuyinj and Chateboi Glacier Dams as triggers for lake outburst cascades in the Karambar Valley, Hindukush

Abstract   Since the mid-nineteenth century, devastating glacier lake outbursts have occurred in the Karambar Valley. The exact source areas of these floods are to date unknown. The present study uses geomorphologic field evidence and interviews of local inhabitants to reconstruct nine potential glacier dams in the Karambar Valley within a horizontal distance of only 40 km. The article focuses on the geomorphologic reconstruction of the highest glacier dams, the Chateboi and Saklei Shuyinj Glaciers. Their lake basins were connected in former times resulting in a complex interfingering of lake sediments, lake terraces and glacial deposits. The outbursts of these lakes could have triggered the drainage of one of the lower ice-dammed lakes (Sokther Rabot, Chillinji, Warghut or Karambar) and therefore initiated an outburst cascade in the upper Karambar Valley. Successive glacier dams are wide spread in the Karakoram, and cascading lakes might have also played a role in other lake outburst scenarios. In the Karambar Valley, even today the Chateboi Glacier blocks the Karambar River over a distance of 4 km and represents a permanent hazard for the villages located downstream.     

Heterogeneous water storage and thermal regime of supraglacial ponds on debris‐covered glaciers

The water storage and energy transfer roles of supraglacial ponds are poorly constrained, yet they are thought to be important components of debris‐covered glacier ablation budgets. We used an unmanned surface vessel (USV) to collect sonar depth measurements for 24 ponds to derive the first empirical relationship between their area and volume applicable to the size distribution of ponds commonly encountered on debris‐covered glaciers. Additionally, we instrumented nine ponds with thermistors and three with pressure transducers, characterizing their thermal regime and capturing three pond drainage events. The deepest and most irregularly‐shaped ponds were those associated with ice cliffs, which were connected to the surface or englacial hydrology network (maximum depth = 45.6 m), whereas hydrologically‐isolated ponds without ice cliffs were both more circular and shallower (maximum depth = 9.9 m). The englacial drainage of three ponds had the potential to melt ~100 ± 20 × 10³ kg to ~470 ± 90 × 10³ kg of glacier ice owing to the large volumes of stored water. Our observations of seasonal pond growth and drainage with their associated calculations of stored thermal energy have implications for glacier ice flow, the progressive enlargement and sudden collapse of englacial conduits, and the location of glacier ablation hot‐spots where ponds and ice cliffs interact. Additionally, the evolutionary trajectory of these ponds controls large proglacial lake formation in deglaciating environments. Copyright © 2017 John Wiley & Sons, Ltd.     

A modified landform development model for the topography of drained thermokarst lake basins in fine‐grained sediments

Permafrost degradation associated with the expansion of thermokarst lakes is commonly interrupted by catastrophic drainage. Subsequently, in tundra areas, permafrost aggradation in drained basins leads to uneven topography characterized by raised centres and wet, depressed margins. The genesis of such topography has been investigated in Old Crow Flats (OCF), a glaciolacustrine plain in the continuous permafrost of northern Yukon. The thermokarst lakes of OCF have a mean depth of only 1.5 m because excess ice is dominantly found only in the uppermost 10 m of the ground. Surface conditions were measured in three drained thermokarst lake basins, including relief, snow conditions, ground temperatures, near‐surface ground ice, and sediment stratigraphy. Four nearby lakes provided information on wave base, shore recession patterns, and bathymetry before drainage: the bottoms of these lakes were not raised in the centre. An elevation difference of up to 2 m was recorded between drained basin margins and centres but was not associated with variations in ice‐wedge density or segregated ice content. Hence basin topography was not controlled by differences in volumetric ground‐ice content between margins and centres. We propose that transport of fine sediment away from eroding lake margins during lake development is the primary mechanism for the genesis of depressed margins and raised centres in drained basins of OCF. Over time, the transport results in the deposition of more and finer sediment in the central parts of lakes, where the lake bottom has subsided below wave base, than at the shallow margins, where resuspension by wave action occurs frequently. This difference in sediment volume is revealed in the topography after drainage, when permafrost aggrades in the lake‐bottom sediment and underlying talik. Copyright © 2016 John Wiley & Sons, Ltd.     

Englacial and subglacial water flow at Skálafellsjökull,Iceland derived from ground penetrating radar,in situ Glacsweb probe and borehole water level measurements

We reconstruct englacial and subglacial drainage at Skálafellsjökull, Iceland, using ground penetrating radar (GPR) common offset surveys, borehole studies and Glacsweb probe data. We find that englacial water is not stored within the glacier (water content ~0–0.3%). Instead, the glacier is mostly impermeable and meltwater is able to pass quickly through the main body of the glacier via crevasses and moulins. Once at the glacier bed, water is stored within a thin (1 m) layer of debris‐rich basal ice (2% water content) and the till. The hydraulic potential mapped across the survey area indicates that when water pressures are high (most of the year), water flows parallel to the margin, and emerges 3 km down glacier at an outlet tongue. GPR data indicates that these flow pathways may have formed a series of braided channels. We show that this glacier has a very low water‐storage capacity, but an efficient englacial drainage network for transferring water to the glacier bed and, therefore, it has the potential to respond rapidly to changes in melt‐water inputs. © 2015 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Long-term impact of the proglacial lake Jökulsárlón on the flow velocity and stability of Breiðamerkurjökull glacier,Iceland

Proglacial lakes are becoming ubiquitous at the termini of many glaciers worldwide due to continued climate warming and glacier retreat, and such lakes have important consequences for the dynamics and future stability of these glaciers. In light of this, we quantified decadal changes in glacier velocity since 1991 using satellite remote sensing for Breiðamerkurjökull, a large lake-terminating glacier in Iceland. We investigated its frontal retreat, lake area change and ice surface elevation change, combined with bed topography data, to understand its recent rapid retreat and future stability. We observed highly spatially variable velocity change from 1991 to 2015, with a substantial increase in peak velocity observed at the terminus of the lake-terminating eastern arm from ~1.00 ± 0.36 m day⁻¹ in 1991 to 3.50 ± 0.25 m day⁻¹ in 2015, with mean velocities remaining elevated from 2008 onwards. This is in stark comparison to the predominately land-terminating arms, which saw no discernible change in their velocity over the same period. We also observed a substantial increase in the area of the main proglacial lake (Jökulsárlón) since 1982 of ~20 km², equating to an annual growth rate of 0.55 km² year⁻¹. Over the same period, the eastern arm retreated by ~3.50 km, which is significantly greater than the other arms. Such discrepancies between the different arms are due to the growth and, importantly, depth increase of Jökulsárlón, as the eastern arm has retreated into its ~300 m-deep reverse-sloping subglacial trough. We suggest that this growth in lake area, forced initially by rising air temperatures, combined with the increase in lake depth, triggered an increase in flow acceleration, leading to further rapid retreat and the initiation of a positive feedback mechanism. These findings may have important implications for how increased melt and calving forced by climate change will affect the future stability of large soft-bedded, reverse-sloped, subaqueous-terminating glaciers elsewhere. © 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd     

Dating of megaflood deposits in the Russian Altai using rock surface luminescence

Catastrophic drainage of ice-dammed lakes in the Altai Mountains has been inferred from geomorphological evidence in the Katun Valley (Russia), and is presumed to have occurred during the Pleistocene. The sedimentary features have been difficult to date directly, due to the absence of organic carbon, and the improbability that luminescence signals in sand grains would be reset during transport. However, the development of rock-surface luminescence dating provides a new opportunity to date the features: clasts have a different transport history to sand grains, and their luminescence depth profiles can be inspected for evidence of bleaching before burial. Here we investigate two sites in the Altai Mountains, and use rock-surface luminescence burial dating to constrain the age of the megaflood deposits. In the Katun Valley, we sampled granite cobbles from a frozen sediment clast emplaced as a dropstone within a massive megaflood gravel terrace. Burial ages for the clasts range from 16.7 to 21.4 ka, with a mean age of 19.8 ± 1.5 ka. This represents the depositional age of the fluvial sediments that preceded the lake outburst flood, (and hence places a maximum age on the catastrophic flood). Clasts sampled from mega-ripples in the Kurai Basin are shown to have a mid-to-late Holocene burial age, which is not consistent with the possible origin of these features during a catastrophic drainage of a glacier-dammed lake. Instead, the burial age of the Kurai Basin sediments may reflect local-scale periglacial or seismic processes along the Kurai Fault Zone.     

Meltwater discharge through the subglacial bed and its land‐forming consequences from numerical experiments in the Polish lowland during the last glaciation

Numerical experiments suggest that the last glaciation severely affected the upper lithosphere groundwater system in NW Poland: primarily its flow pattern, velocities and fluxes. We have simulated subglacial groundwater flow in two and three spatial dimensions using finite difference codes for steady‐state and transient conditions. The results show how profoundly the ice sheet modifies groundwater pressure heads beneath and some distance beyond the ice margin. All model runs show water discharge at the ice forefield driven by ice‐sheet‐thickness‐modulated, down‐ice‐decreasing hydraulic heads. In relation to non‐glacial times, the transient 3D model shows significant changes in the groundwater flow directions in a regionally extensive aquifer ca. 90 m below the ice–bed interface and up to 40 km in front of the glacier. Comparison with empirical data suggests that, depending on the model run, only between 5 and 24% of the meltwater formed at the ice sole drained through the bed as groundwater. This is consistent with field observations documenting abundant occurrence of tunnel valleys, indicating that the remaining portion of basal meltwater was evacuated through a channelized subglacial drainage system. Groundwater flow simulation suggests that in areas of very low hydraulic conductivity and adverse subglacial slopes water ponding at the ice sole was likely. In these areas the relief shows distinct palaeo‐ice lobes, indicating fast ice flow, possibly triggered by the undrained water at the ice–bed interface. Owing to the abundance of low‐permeability strata in the bed, the simulated groundwater flow depth is less than ca. 200 m. Copyright © 2009 John Wiley & Sons, Ltd.     

Complex kame belt morphology,stratigraphy and architecture

The development of glacier karst at the margins of melting ice sheets produces complex glaciofluvial sediment-landform assemblages that provide information on ice sheet downwasting processes. We present the first combined geomorphological, sedimentological and geophysical investigation of the Brampton Kame Belt, an important glaciofluvial depositional zone at the centre of the last British-Irish Ice Sheet. Ground-penetrating radar (GPR) data allow the broad scale internal architecture of ridges (eskers) and flat-topped hills (ice-walled lake plains) to be determined at four sites. In combination with sediment exposures, these provide information on lateral and vertical variations in accretion styles, depositional boundaries, and grain size changes. Building on existing work on the subject, we propose a refined model for the formation of ice-walled lake plains resulting from the evolution and collapse of major drainage axes into lakes as stable glacier karst develops during deglaciation. The internal structure of esker ridges demonstrates variations in sedimentation that can be linked to differences in ridge morphologies across the kame belt. This includes low energy flow conditions and multiple accretion phases identified within large S-N oriented esker ridges; and fluctuating water pressures, hyperconcentrated flows, and significant deformation within a fragmented SW–NE oriented esker ridge. In combination with updated geomorphological mapping, this work allows us to identify two main styles of drainage within the kame belt: (1) major drainage axes aligned broadly S-N that extend through the entire kame belt and collapsed into a chain of ice-walled lakes; and (2) a series of smaller, fragmented SW–NE aligned esker ridges that represent ice-marginal drainage as the ice sheet receded south-eastwards up the Vale of Eden. Our study demonstrates the importance of integrated geomorphological, sedimentological and geophysical investigations in order to understand complex and polyphase glaciofluvial sediment-landform assemblages. © 2019 John Wiley & Sons, Ltd. © 2019 John Wiley & Sons, Ltd.     

The hydrochemistry of Lake Vostok and the potential for life in Antarctic subglacial lakes

Our understanding of Lake Vostok, the huge subglacial lake beneath the East Antarctic Ice Sheet, has improved recently through the identification of key physical and chemical interactions between the ice sheet and the lake. The north of the lake, where the overlying ice sheet is thickest, is characterized by subglacial melting, whereas freezing of lake water occurs in the south, resulting in ～210 m of ice accretion to the underside of the ice sheet. The accreted ice contains lower concentrations of the impurities normally found in glacier ice, suggesting a net transfer of material from meltwater into the lake. The small numbers of microbes found so far within the accreted ice have DNA profiles similar to those of contemporary surface microbes. Microbiologists expect, however, that Lake Vostok, and other subglacial lakes, will harbour unique species, particularly within the deeper waters and associated sediments. The extreme environments of subglacial lakes are characterized by high pressures, low temperatures, permanent darkness, limited nutrient availability, and oxygen concentrations that are derived from the ice that provides the meltwater. Copyright © 2002 John Wiley & Sons, Ltd.     

Ice flow‐unit influence on glacier structure,debris entrainment and transport

The links between structural glaciology, glacial debris entrainment and transport have been established in a number of different glacier settings. Here we document the structural evolution of a temperate Alpine valley glacier (Vadrec del Forno, Switzerland) and demonstrate that individual flow units within the glacier have very different structural and debris characteristics. The glacier consists of a broad accumulation area with multiple basins feeding a relatively narrow tongue and is formed from six distinct flow units. Each flow unit has its own characteristic structural assemblage. Flow units that narrow rapidly down‐glacier are dominated by primary stratification that has evolved into longitudinal foliation. In contrast, wider flow units preferentially develop an axial planar foliation. Glacier structure plays a limited role in the entrainment of debris, which is more strongly influenced by ice‐marginal rockfall and avalanche inputs onto the glacier surface. However, once entrained, glacier structure controls the reorientation and redistribution of debris within the ice mass. By taking a whole‐glacier approach to describing glacier structure and debris transport, we conclude that individual flow units are unique with regard to structure and debris transfer. Copyright © 2013 John Wiley & Sons, Ltd.     

The formation of supraglacial debris covers by primary dispersal from transverse englacial debris bands

Glaciological controls on debris cover formation are investigated from the perspective of primary dispersal of supraglacial debris across a melting ice surface. This involves the migration of angled debris septa outcrops across a melting, thinning glacier ablation zone. Three measures of a glacier's ability to evacuate supraglacial debris are outlined: (1) a concentration factor describing the focusing of englacial debris into specific supraglacial mass loads; (2) the rate of migration of a septum outcrop relative to the local ice surface; and (3) a downstream velocity differential between a slower septum outcrop and the faster ice surface velocity. Measures (1) and (2) are inversely related, while measure (3) increases down‐glacier to explain why slow‐moving, thinning ice rapidly becomes debris covered. Data from Glacier d'Estelette (Italian Alps) are used to illustrate these processes, and to explore the potential for debris cover formation and growth in different glaciological environments. The transition from a ‘clean’, transport‐dominated to a debris‐covered ablation‐dominated glacier is explained by the melting out of more closely‐spaced debris septa, in combination with the geometric interactions of angled septa and ice surface in a field of reducing flow and increasing ablation. The growth and shrinkage of debris covers are most sensitive to glaciological changes at glaciers with gently‐dipping debris‐bearing foliation, but less sensitive at high‐compression glaciers whose termini are constrained by moraine dams and other forms of obstruction. These findings show that a variety of debris‐covered glacier types will show a spectrum of response characteristics to negative mass balance. Copyright © 2013 John Wiley & Sons, Ltd.     

Debris entrainment by basal freeze‐on and thrusting during the 1995–1998 surge of Kuannersuit Glacier on Disko Island,west Greenland

Kuannersuit Glacier, a valley glacier on Disko Island in west Greenland, experienced a major surge from 1995 to 1998 where the glacier advanced 10·5 km and produced a ～65 m thick stacked sequence of debris‐rich basal ice and meteoric glacier ice. The aim of this study is to describe the tectonic evolution of large englacial thrusts and the processes of basal ice formation using a multiproxy approach including structural glaciology, stable isotope composition (δ¹⁸O and δD), sedimentology and ground‐penetrating radar. We argue that the major debris layers that can be traced in the terminal zone represent englacial thrusts that were formed early during the surge. Thrust overthrow was at least 200–300 m and this lead to a 30 m thick repetition of basal ice at the ice margin. It is assumed that the englacial thrusting was initiated at the transition between warm ice from the interior and the cold snout. The basal debris‐rich ice was mainly formed after the thrusting phase. Two sub‐facies of stratified basal ice have been identified; a lower massive ice facies (S_M) composed of frozen diamict enriched with heavy stable isotopes overlain by laminated ice facies (S_L) consisting of millimetre thick lamina of alternating debris‐poor and debris‐rich ice. We interpret the stratified basal ice as a continuum formed mainly by freeze‐on processes and localized regelation. First laminated basal ice is formed and as meltwater is depleted more sediment is entrained and finally the glacier freezes to the base and massive diamict is frozen‐on. The increased ability to entrain sediments may partly be associated with higher basal freezing rates enhanced by loss of frictional heat from cessation of fast flow and conductive cooling through a thin heavily crevassed ice during the final phase of the glacier surge. The dispersed basal ice facies (D) was mainly formed by secondary processes where fine‐grained sediment is mobilized in the vein system of ice. Our results have important implications for understanding the significance of basal ice formation and englacial thrusting beneath fast‐flowing glaciers and it provides new information about the development of landforms during a glacier surge. Copyright © 2010 John Wiley & Sons, Ltd.     

用氚讨论中国西部新疆干旱区的水循环率

Tritium concentrations are used to trace water circulation in the Urumqi and Turfan basins in the Xinjiang, western China. Tritium analyses were made for 77 water samples of river waters, groundwaters, spring waters, lake waters and glacier ice collected in summers in 1992 and 1994. The tritium concentrations in the waters are in a wide range from 0 to 125 TU, most of which are considerably high compared with those of most waters in Japan, because tritium levels in precipitation in the area are over ten times as high as those in Japan. River waters originating in glacier regions contain melt glacier, the proportion of which is over 0.5 to river water. The mean resi-dence time of circulating meteoric water in the mountain regions is estimated to be about 15 years. Most groundwaters and spring waters in the flat regions are mainly derived from river waters originating in glacier regions. The groundwater of greatest tritium concentrations in wells in Urumqi City is derived from Urumqi River about 25 years ago. It takes several ten years for river water to pass the underground to many springs. Some groundwaters and spring waters have taken a long time more than 40 years to travel under the ground. Enrichment of tritium in lake water by evaporation is considered to estimate the contribution of groundwater flow to the recharge of lake. Various contributions of groundwater to lakes are inferred for the various type of salinity in closed or semi-closed lakes. The inflow rates of groundwater to salt lakes are small as against fresh water lakes.     

Englacial vs lacustrine origin of volcanic table mountains: evidence from iceland

Detailed facies analysis of hyaloclastites and associated lavas from eight table mountains and similar "hyaloclastite volcanoes" in the Icelandic rift zone contradict a rapid and continuous, "monogenetic", entirely subglacial evolution of most volcanoes studied. The majority of the exposed hyaloclastite deposits formed in large, stable lakes as indicated by widespread, up to 300-m-thick, continuous sections of deep water, shallow water and emergent facies. Salient features include extensively layered or bedded successions comprising mainly debris flow deposits, turbidites, base surge and fallout deposits consisting of texturally and compositionally variable, slightly altered hyaloclastites, as well as sheet and pillow lavas. In contrast, chaotic assemblages of coarser-grained, more poorly sorted and more strongly palagonitized hyaloclastite tuffs and breccias, as well as scoria and lava are interpreted to have formed under sub- or englacial conditions in small, chimney-like ice cavities or ice-bound lakes. Irregularly shaped and erratically arranged hyaloclastite bodies produced at variable water levels appear to have resulted mainly from rapid changes of the eruptive environment due to repeated build-up and drainage of ice-bound lakes as well as the restricted space between the ice walls. We distinguish a "deep water" facies formed during high water levels of the lake, a hydroclastic shallow water and emergent facies (leakage of the lake or growth of the volcano above the water surface). Our model implies the temporary existence of large, stable lakes in Iceland probably formed by climatically induced ice melting. The highly complex edifices of many table mountains and similar volcanoes were constructed during several eruptive periods in changing environments characterized by contrasting volcanic and sedimentary processes. Received: 10 June 1997 / Accepted: 28 July 1998     

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.     

Quantitative comparison of river inflows to a rapidly expanding lake in central Tibetan Plateau

The recent rapid expansion of inland lakes on the Tibetan Plateau (TP) are a good indicator of the consequences of climate change. Quantifying the hydrological cycle of the lake basin is fundamentally important to understand the causes of lake growth. However, the hydrological processes of the TP interior are very complex and difficult to investigate because of the lack of observations. This is especially true for estimating the lake changes when run‐off inflows are affected by small lakes located in the flow routes within drainage areas. We used an integrated hydrological model, in combination with glacier melt and lake retention models, to analyse the run‐off inflows to Lake Siling Co, the largest endorheic lake in Tibet. It includes four subdrainage basins: Zhajiazangbu, Zhagenzangbu, Alizangbu, and Boquzangbu. Lake Siling Co was characterized by considerable increases during warm season from 1981 to 2012, due to the increased run‐off from Zhajiazangbu accounting for about 51–62% of the total run‐off inflows. Moreover, the dramatic increases exhibited during cold seasons were related to the increased retention water released from the small lakes within Zhagenzangbu and Alizangbu. Of the studied subdrainage basins, Boquzangbu contributed the least during both warm and cold seasons. On average, the annual amount of evaporation from lakes within the drainage area was about 2 times greater than that of glacier melt run‐off. Our results suggest that the retention effects of lakes on river inflows should receive more attention, because understanding these effects is potentially crucial to improved understanding of lake variations in the TP.     

冰川融水对山地冰冻圈冰湖水文效应的影响

冰川融水通过热量、水、物质传输对山地冰冻圈冰湖水文效应产生影响,引起广泛关注.本文从山地冰冻圈冰湖的水量、物理化学性质、生物等方面系统总结冰川融水对冰湖水文效应的影响.冰川融水被冰湖滞留能在一定程度上延缓区域冰川水资源的亏损,但也直接导致了潜在危险性冰湖数量和危险程度增大.冰川融水对冰湖物理性质的影响主要表现在降低湖水温度、影响透明度/浊度、改变湖水密度、造成湖水热力分层现象等方面,对冰湖化学性质的影响主要表现在增加湖水中的氮素、溶解有机物、持久性有机污染物、各类离子和重金属等,进而影响冰湖生物的分布、组成、结构和功能.深入系统地开展冰川融水及其变化对冰湖水文效应研究,对冰川水文与水资源、山地冰冻圈生态环境研究具有重要意义.     

Processes at the margins of supraglacial debris cover: Quantifying dirty ice ablation and debris redistribution

Current glacier ablation models have difficulty simulating the high-melt transition zone between clean and debris-covered ice. In this zone, thin debris cover is thought to increase ablation compared to clean ice, but often this cover is patchy rather than continuous. There is a need to understand ablation and debris dynamics in this transition zone to improve the accuracy of ablation models and the predictions of future debris cover extent. To quantify the ablation of partially debris-covered ice (or ‘dirty ice’), a high-resolution, spatially continuous ablation map was created from repeat unmanned aerial systems surveys, corrected for glacier flow in a novel way using on-glacier ablation stakes. Surprisingly, ablation is similar (range ~ 5 mm w.e. per day) across a wide range of percentage debris covers (~ 30–80%) due to the opposing effects of a positive correlation between percentage debris cover and clast size, countered by a negative correlation with albedo. Once debris cover becomes continuous, ablation is significantly reduced (by 61.6% compared to a partial debris cover), and there is some evidence that the cleanest ice (<~ 15% debris cover) has a lower ablation than dirty ice (by 3.7%). High-resolution feature tracking of clast movement revealed a strong modal clast velocity where debris was continuous, indicating that debris moves by creep down moraine slopes, in turn promoting debris cover growth at the slope toe. However, not all slope margins gain debris due to the removal of clasts by supraglacial streams. Clast velocities in the dirty ice area were twice as fast as clasts within the continuously debris-covered area, as clasts moved by sliding off their boulder tables. These new quantitative insights into the interplay between debris cover characteristics and ablation can be used to improve the treatment of dirty ice in ablation models, in turn improving estimates of glacial meltwater production. © 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd     

Differences and cause analysis of changes in lakes of different supply types in the north‐western Tibetan Plateau

Most lakes of the Tibetan Plateau were experiencing quick expansion in recent decades; a detailed study on the changes in lakes of different supply types will help to understand the cause of the changes by analysing area change of 34 lakes and water level change of eight lakes in the north‐western Tibetan Plateau. All lakes are classified into three types: non‐glacier‐fed lakes, upstream lakes and glacier‐fed lakes. The glacier‐fed lakes are separated into glacier‐fed_P1 (quick expansion region) and glacier‐fed_P2 (slow expansion region). Combining the changes in precipitation, temperature and evaporation, less precipitation may be the main reason for lake shrinkage from 1976 to 1996 and quickly increasing precipitation led to the lakes' expansion from 1996 to 2000. However, after 2000, upstream lakes remained in a stable state with slight variation, non‐glacier‐fed lakes and glacier‐fed_P2 lakes exhibited a slightly increasing rate with high precipitation and high evaporation, and glacier‐fed_P1 lakes still expanded quickly. On the basis of the assumption of same precipitation and evaporation, glacial meltwater may make an important contribution (>52%) to the increase in water storage of Bangdag Co and Meima Co. The results suggest that glacial meltwater augments the increased rate of glacier‐fed_P2 lakes and plays a much more important role in the expansion of glacier‐fed_P1 lakes compared to other lakes.     

Sensitivity of subglacial Lake Vostok's flow regime on environmental parameters

Subglacial lakes provide unique habitats, but the exact nature of physical and geochemical conditions are still a matter of debate and await direct sampling of water. Due to its isolation from external atmospheric forcing other environmental parameters influence the flow characteristics within the lake. In this study we use an improved treatment of the physical processes at the ice–water boundary interface to identify and quantify the impact of (1) the geothermal heat flux, (2) the heat flux from the lake into the ice, (3) the influence of the salinity of the lake water, and (4) the ice thickness on the size of the freezing area and the freeze/melt rates. We show that the modelled basal mass imbalance (that is the produced melt water minus the re-frozen water) depends on the geothermal heating as well as the heat flux into the ice. The circulation and the temperature distribution within subglacial Lake Vostok are rather stable against variations of geothermal heat flux, heat flux into the ice sheet, salinity of the lake, and small changes of the ice thickness above the lake. However, the flow regime for any subglacial lake with less than 2000 m ice thickness above, will be substantially different from those that experience higher pressures. This is because the buoyancy–temperature relationship reverses at this depth.