Khumbu Glacier Research Articles

Heterogeneous decadal glacier downwasting at the Mt. Everest (Qomolangma) from 2000 to ~ 2012 based on multi-baseline bistatic SAR interferometry

Remote sensing based geodetic observations can be used as alternatives data to map glacier height changes because the harsh environment complicates in situ observations. In this study, we analysed five pairs of X-band bistatic TerraSAR-X/TanDEM-X images from 2011 and 2014 for Mt. Everest (Qomolangma). Glacier height changes were derived using the D-InSAR method respect to SRTM DEM of 2000. An iterative D-InSAR method using multi-baseline bistatic SAR interferograms was employed, which increased reliability and accuracy of glacier height changing observations. From 2000 to ~2012, the geodetic glacier mass balance for the Mt. Everest and the surrounding region was −0.38±0.04m w.e. a−1. The spatial pattern of the glacier mass loss was heterogeneous. The regional heterogeneity may possibly reflect debris-covering rates, terminating type, temperature rising rates and glacier flow rates. Comparison to long-period geodetic glacier mass balance data provided by previous studies since 1970 revealed that our results showed more rapid increases in the glacier mass loss rates after 2000 in the area around Khumbu Glacier in the southern slope of the Mt. Everest, whereas glacier mass loss rates kept stable in the Rongbuk Catchment at its northern slope.

Everest's thinning glaciers: implications for tourism and mountaineering

Glacier mass loss in the Everest region of Nepal is accelerating in response to a warming climate, which is a trend observed across the central and eastern Himalaya. Thinning glaciers and the development of supraglacial (surface) ponds and large glacial lakes will increasingly restrict access to glacier surfaces and will affect popular trekking routes and mountaineering activities in the region. Through quantifying glacier accessibility and supraglacial pond expansion, we estimate that the Kongma La Pass trail across the Khumbu Glacier is likely to be impassable by 2020 due to supraglacial pond expansion and glacier thinning, and will require significant re‐routing. An estimated 197 649 227 m3 of ice melted over the period 1984–2015 on the Khumbu Glacier. Additionally, expert opinion from Everest mountaineers suggest that rockfall activity is likely to increase in the high‐mountain environment as snow and ice melts from mountain slopes, requiring changes to climbing routes on the world's highest peaks. Similarly, route difficulty will be affected by changing monsoon precipitation patterns, which determines windows of opportunity for ascents, and the distribution and quantity of snowfall. We conclude that increased collaboration between the scientific, local, and mountaineering communities offers mutual benefits for data collection and dissemination, and we identify key areas that should be investigated further.

1980年ヒマラヤ山脈エベレスト・クンブ氷河ウエスタンクームで掘削された 浅層アイスコア解析

1980年ヒマラヤ山脈エベレスト・クンブ氷河ウエスタンクームで掘削された 浅層アイスコア解析

Supraglacial Ponds Regulate Runoff From Himalayan Debris‐Covered Glaciers

AbstractMeltwater and runoff from glaciers in High Mountain Asia is a vital freshwater resource for one‐fifth of the Earth's population. Between 13% and 36% of the region's glacierized areas exhibit surface debris cover and associated supraglacial ponds whose hydrological buffering roles remain unconstrained. We present a high‐resolution meltwater hydrograph from the extensively debris‐covered Khumbu Glacier, Nepal, spanning a 7 month period in 2014. Supraglacial ponds and accompanying debris cover modulate proglacial discharge by acting as transient and evolving reservoirs. Diurnally, the supraglacial pond system may store >23% of observed mean daily discharge, with mean recession constants ranging from 31 to 108 h. Given projections of increased debris cover and supraglacial pond extent across High Mountain Asia, we conclude that runoff regimes may become progressively buffered by the presence of supraglacial reservoirs. Incorporation of these processes is critical to improve predictions of the region's freshwater resource availability and cascading environmental effects downstream.

Quantifying ice cliff evolution with multi-temporal point clouds on the debris-covered Khumbu Glacier, Nepal

ABSTRACTMeasurements of glacier ice cliff evolution are sparse, but where they do exist, they indicate that such areas of exposed ice contribute a disproportionate amount of melt to the glacier ablation budget. We used Structure from Motion photogrammetry with Multi-View Stereo to derive 3-D point clouds for nine ice cliffs on Khumbu Glacier, Nepal (in November 2015, May 2016 and October 2016). By differencing these clouds, we could quantify the magnitude, seasonality and spatial variability of ice cliff retreat. Mean retreat rates of 0.30–1.49 cm d−1 were observed during the winter interval (November 2015–May 2016) and 0.74–5.18 cm d−1 were observed during the summer (May 2016–October 2016). Four ice cliffs, which all featured supraglacial ponds, persisted over the full study period. In contrast, ice cliffs without a pond or with a steep back-slope degraded over the same period. The rate of thermo-erosional undercutting was over double that of subaerial retreat. Overall, 3-D topographic differencing allowed an improved process-based understanding of cliff evolution and cliff-pond coupling, which will become increasingly important for monitoring and modelling the evolution of thinning debris-covered glaciers.

Evaluating morphological estimates of the aerodynamic roughness of debris covered glacier ice

AbstractAerodynamic roughness length (z0), the height above the ground surface at which the extrapolated horizontal wind velocity profile drops to zero, is one of the most poorly parameterised elements of the glacier surface energy balance equation. Microtopographic methods for estimating z0 have become prominent in the literature in recent years, but are rarely validated against independent measures and are yet to be comprehensively analysed for scale or data resolution dependency. Here, we present the results of a field investigation conducted on the debris covered Khumbu Glacier during the post‐monsoon season of 2015. We focus on two sites. The first is characterised by gravels and cobbles supported by a fine sandy matrix. The second comprises cobbles and boulders separated by voids. Vertical profiles of wind speed recorded by a tower comprising five cup anemometers and deployed over both sites enable us to derive measurements of aerodynamic roughness that reflect their observed surface characteristics (0.0184 m and 0.0243 m, respectively). At the second site, z0 also varied through time following snowfall (0.0055 m) and during its subsequent melt (0.0129 m), showing the importance of fine resolution topography for near‐surface airflow. To compare the wind profile data with microtopographic methods, we conducted structure from motion multi‐view stereo (SfM‐MVS) surveys across each patch and calculated z0 using three previously published approaches. The fully three‐dimensional cloud‐based approach is shown to be most stable across different scales and these z0 values are most correct in relative order when compared with the wind tower data. Popular profile‐based methods perform less well providing highly variable values across different scales and when using data of differing resolution. These findings hold relevance for all studies using microtopographic methods to estimate aerodynamic roughness lengths, including those in non‐glacial settings. Copyright © 2017 John Wiley & Sons, Ltd.

Using structure‐from‐motion to create glacier DEMs and orthoimagery from historical terrestrial and oblique aerial imagery

AbstractIncreased resolution and availability of remote sensing products, and advancements in small‐scale aerial drone systems, allows observations of glacial changes at unprecedented levels of detail. Software developments, such as structure‐from‐motion (SfM), now allow users an easy and efficient method to generate three‐dimensional (3D) models and orthoimages from aerial or terrestrial datasets. While these advancements show promise for current and future glacier monitoring, many regions still suffer a lack of observations from earlier time periods. We report on the use of SfM to extract spatial information from various historic imagery sources. We focus on three geographic regions, the European Alps, high Arctic Norway and the Nepal Himalayas. We used terrestrial field photographs from 1896, high oblique aerial photographs from 1936 and aerial handheld photographs from 1978 to generate digital elevation models (DEMs) and orthophotos of the Rhone glacier, Brøggerhalvøya and the lower Khumbu glacier, respectively. Our analysis shows that applying SfM to historic imagery can generate high quality models using only ground control points. Limited camera/orientation information was largely reproduced using self‐calibrated model data. Using these data, we calculated mean ground sampling distances across each site which demonstrates the high potential resolution of resulting models. Vertical errors for our models are ±5.4 m, ±5.2 m and ±3.3 m. Differencing shows similar patterns of thinning at lower Rhone (European Alps) and Brøggerhalvøya (Norway) glaciers, which have mean thinning rates of 0.31 m a−1 (1896–2010) to 0.86 m a−1 (1936–2010) respectively. On these clean ice glaciers thinning is highest in the terminus region and decreasing up‐glacier. In contrast to these glaciers, uneven topography, exposed ice‐cliffs and debris cover on the Khumbu glacier create a highly variable spatial distribution of thinning. The mean thinning rate for the Khumbu study area was found to be 0.54 ± 0.9 m a−1 (1978–2015). Copyright © 2017 John Wiley & Sons, Ltd.

The dynamics of supraglacial ponds in the Everest region, central Himalaya

The dynamics of supraglacial pond development in the Everest region are not well constrained at a glacier scale, despite their known importance for meltwater storage, promoting ablation, and transmitting thermal energy englacially during drainage events. Here, we use fine-resolution (~0.5–2m) satellite imagery to reveal the spatiotemporal dynamics of 9340 supraglacial ponds across nine glaciers in the Everest region, ~2000–2015. Six of our nine study glaciers displayed a net increase in ponded area over their observation periods. However, large inter- and intra-annual changes in ponded area were observed of up to 17% (Khumbu Glacier), and 52% (Ama Dablam) respectively. Additionally, two of the fastest expanding lakes (Spillway and Rongbuk) partially drained over our study period. The Khumbu Glacier is developing a chain of connected ponds in the lower ablation area, which is indicative of a trajectory towards large lake development. We show that use of medium-resolution imagery (e.g. 30m Landsat) is likely to lead to large classification omissions of supraglacial ponds, on the order of 15–88% of ponded area, and 77–99% of the total number of ponds. Fine-resolution imagery is therefore required if the full spectrum of ponds that exist on the surface of debris-covered glaciers are to be analysed.

Assessment of the evolution in velocity of two debris‐covered valley glaciers in nepal and new zealand

Feature tracking of orthorectified pairs of Advanced Spaceborne Thermal Emission and Reflection Radiometer satellite images is used to calculate velocities for the Tasman Glacier, New Zealand (2002–2014) and the Khumbu Glacier, Nepal (2001–2008). Velocities in the middle and upper ablation zones of both glaciers show a long‐term decrease of ∼10–20%, while the terminus of Khumbu Glacier has remained near stagnation throughout the study period. In contrast, there has been a recent acceleration of the lower terminus of Tasman Glacier, from ∼5 m a–1 in 2002 to 40 m a–1 in 2014. Both of these glaciers have an extensive supraglacial debris cover across their lower ablation regions, with the Khumbu Glacier terminating on land and the Tasman Glacier terminating in a proglacial lake. The rapid recent increase in velocity of the terminus of Tasman Glacier is closely correlated with the increase in size of its proglacial lake. These results indicate the complex dynamic changes that mountain valley glaciers may undergo in response to long‐term negative mass balance.

Modelling the feedbacks between mass balance, ice flow and debris transport to predict the response to climate change of debris-covered glaciers in the Himalaya

Many Himalayan glaciers are characterised in their lower reaches by a rock debris layer. This debris insulates the glacier surface from atmospheric warming and complicates the response to climate change compared to glaciers with clean-ice surfaces. Debris-covered glaciers can persist well below the altitude that would be sustainable for clean-ice glaciers, resulting in much longer timescales of mass loss and meltwater production. The properties and evolution of supraglacial debris present a considerable challenge to understanding future glacier change. Existing approaches to predicting variations in glacier volume and meltwater production rely on numerical models that represent the processes governing glaciers with clean-ice surfaces, and yield conflicting results. We developed a numerical model that couples the flow of ice and debris and includes important feedbacks between debris accumulation and glacier mass balance. To investigate the impact of debris transport on the response of a glacier to recent and future climate change, we applied this model to a large debris-covered Himalayan glacier—Khumbu Glacier in Nepal. Our results demonstrate that supraglacial debris prolongs the response of the glacier to warming and causes lowering of the glacier surface in situ, concealing the magnitude of mass loss when compared with estimates based on glacierised area. Since the Little Ice Age, Khumbu Glacier has lost 34% of its volume while its area has reduced by only 6%. We predict a decrease in glacier volume of 8–10% by AD2100, accompanied by dynamic and physical detachment of the debris-covered tongue from the active glacier within the next 150 yr. This detachment will accelerate rates of glacier decay, and similar changes are likely for other debris-covered glaciers in the Himalaya.

Temporal Changes in Elevation of the Debris-Covered Ablation Area of Khumbu Glacier in the Nepal Himalaya since 1978

We evaluated elevation changes at four sites on debris-covered ablation area of Khumbu Glacier, Nepal Himalaya, since 1978. In 2004, we carried out a ground survey by differential GPS in the upper- and lowermost areas of the ablation area. The amount of surface lowering was calculated by comparing digital elevation models (DEMs) with 30-m grid size, as generated from survey data corrected in 1978, 1995, and in the present study. Because we could not access the middle parts of the debris-covered area due to surface roughness, for this area we used an ASTER-DEM calibrated by the ground survey data. The amount of surface lowering during the period 1978–2004 was insignificant near the terminus. A remarkable acceleration of surface lowering was found in the middle part of the debris-covered ablation area, where the glacier surface is highly undulating. In the uppermost area, surface lowering has continued at a steady rate. Surface flow speeds have decreased since 1956, revealing that the recent decrease in ice flux from the upper accumulation area would have accelerated the rate of surface lowering of the debris-covered area of Khumbu Glacier during the period 1995–2004.

Sedimentological, geomorphological and dynamic context of debris-mantled glaciers, Mount Everest (Sagarmatha) region, Nepal

This paper presents the sediment, landform and dynamic context of four avalanche-fed valley glaciers (Khumbu, Imja, Lhotse and Chukhung) in the Mount Everest (Sagarmatha) region of Nepal. All four glaciers have a mantle of debris dominated by sandy boulder-gravel that suppresses melting to an increasing degree towards the snout, leading to a progressive reduction in the overall slope of their longitudinal profile. Prominent lateral–terminal moraine complexes, also comprising sandy bouldergravel, enclose the glaciers. These terminal moraines originally grew by accretion of multiple sedimentary facies of basal glacial and supraglacial origin, probably by folding and thrusting when the glaciers were more dynamic during the Little Ice Age. The four glaciers are in various stages of recession, and demonstrate a range of scenarios from down-wasting of the glacier tongue, through morainedammed lake development, to post-moraine-dam breaching. Khumbu Glacier is at the earliest stage of supraglacial pond formation and shows no sign yet of developing a major lake, although one is likely to develop behind its >250 m high composite terminal moraine. Imja Glacier terminates in a substantial body of water behind a partially ice-cored moraine dam (as determined from geophysical surveys), but morphologically appears unlikely to be an immediate threat. Chukhung Glacier already has a breached moraine and a connected debris fan, and therefore no longer poses a threat. Lhotse Glacier has an inclined, free-draining tongue that precludes hazardous lake development. From the data assembled, a conceptual model, applicable to other Himalayan glaciers, is proposed to explain the development of large, lateral-terminal moraine complexes and associated potentially hazardous moraine dams. – 2008 Elsevier Ltd. All rights reserved .

Englacial drainage systems formed by hydrologically driven crevasse propagation

AbstractRecent work has shown that surface-to-bed drainage systems re-form annually on parts of the Greenland ice sheet and some High Arctic glaciers, leading to speed-up events soon after the onset of summer melt. Surface observations and geophysical data indicate that such systems form by hydrologically driven fracture propagation (herein referred to as ‘hydrofracturing’), although little is known about their characteristics. Using speleological techniques, we have explored and surveyed englacial drainage systems formed by hydrofracturing in glaciers in Svalbard, Nepal and Alaska. In Hansbreen, Svalbard, vertical shafts were followed through ∼60 m of cold ice and ∼10 m of temperate basal ice to a subglacial conduit. Deep hydrofracturing occurred at this site due to a combination of extensional ice flow and abundant surface meltwater at a glacier confluence. The englacial drainage systems in Khumbu Glacier, Nepal, and Matanuska Glacier, Alaska, USA, formed in areas of longitudinal compression and transverse extension and consist of vertical slots that plunge down-glacier at angles of 55° or less. The occurrence of englacial drainages initiated by hydrofracturing in diverse glaciological regimes suggests that it is a very widespread process, and that surface-to-bed drainage can occur wherever high meltwater supply coincides with ice subjected to sufficiently large tensile stresses.

A cut-and-closure origin for englacial conduits in uncrevassed regions of polythermal glaciers

AbstractOn uncrevassed regions of polythermal glaciers, englacial conduits can form by incision of supraglacial stream channels followed by roof closure. The origin and evolution of examples in Longyearbreen, Svalbard, and Khumbu Glacier, Nepal, were determined by speleological survey. The development of perennial incised channels requires that incision is significantly faster than glacier surface ablation, and thus will be favoured by high meltwater discharges in combination with cool climatic conditions or thick debris cover. Incised canyons can become blocked by drifted winter snow, refrozen meltwater, ice rafting from non-local sources (allochthonous breccias) and roof collapses (autochthonous breccias). Conduit closure can also occur in response to ice creep, particularly at depth. Following isolation from the surface, englacial conduits continue to evolve by vadose incision down to local base level. In the case of Longyearbreen, incision allowed the channel to reach the glacier bed, but on Khumbu Glacier deep incision is prevented because an effectively impermeable terminal moraine provides a high base level for the glacier drainage system. During our period of observations, deeper parts of the Longyearbreen conduit became blocked by a combination of ice accumulation and creep, causing the stream course to be re-routed to higher levels. On that glacier, incision, blockage and upward re-routing are cyclic. We conclude that ‘cut and closure’ is the dominant mechanism of englacial conduit formation on uncrevassed regions of polythermal glaciers.

Identification of glacier motion and potentially dangerous glacial lakes in the Mt. Everest region/Nepal using spaceborne imagery

Abstract. Failures of glacial lake dams can cause outburst floods and represents a serious hazard. The potential danger of outburst floods depends on various factors like the lake's area and volume, glacier change, morphometry of the glacier and its surrounding moraines and valley, and glacier velocity. Remote sensing offers an efficient tool for displacement calculations and risk assessment of the identification of potentially dangerous glacial lakes (PDGLs) and is especially helpful for remote mountainous areas. Not all important parameters can, however, be obtained using spaceborne imagery. Additional interpretation by an expert is required. ASTER data has a suitable accuracy to calculate surface velocity. Ikonos data offers more detail but requires more effort for rectification. All investigated debris-covered glacier tongues show areas with no or very slow movement rates. From 1962 to 2003 the number and area of glacial lakes increased, dominated by the occurrence and almost linear areal expansion of the moraine-dammed lakes, like the Imja Lake. Although the Imja Lake will probably still grow in the near future, the risk of an outburst flood (GLOF) is considered not higher than for other glacial lakes in the area. Potentially dangerous lakes and areas of lake development are identified. There is a high probability of further lake development at Khumbu Glacier, but a low one at Lhotse Glacier.

Sedimentological, geomorphological and dynamic context of debris-mantled glaciers, Mount Everest (Sagarmatha) region, Nepal

This paper presents the sediment, landform and dynamic context of four avalanche-fed valley glaciers (Khumbu, Imja, Lhotse and Chukhung) in the Mount Everest (Sagarmatha) region of Nepal. All four glaciers have a mantle of debris dominated by sandy boulder-gravel that suppresses melting to an increasing degree towards the snout, leading to a progressive reduction in the overall slope of their longitudinal profile. Prominent lateral–terminal moraine complexes, also comprising sandy boulder-gravel, enclose the glaciers. These terminal moraines originally grew by accretion of multiple sedimentary facies of basal glacial and supraglacial origin, probably by folding and thrusting when the glaciers were more dynamic during the Little Ice Age. The four glaciers are in various stages of recession, and demonstrate a range of scenarios from down-wasting of the glacier tongue, through moraine-dammed lake development, to post-moraine-dam breaching. Khumbu Glacier is at the earliest stage of supraglacial pond formation and shows no sign yet of developing a major lake, although one is likely to develop behind its >250 m high composite terminal moraine. Imja Glacier terminates in a substantial body of water behind a partially ice-cored moraine dam (as determined from geophysical surveys), but morphologically appears unlikely to be an immediate threat. Chukhung Glacier already has a breached moraine and a connected debris fan, and therefore no longer poses a threat. Lhotse Glacier has an inclined, free-draining tongue that precludes hazardous lake development. From the data assembled, a conceptual model, applicable to other Himalayan glaciers, is proposed to explain the development of large, lateral-terminal moraine complexes and associated potentially hazardous moraine dams.

The potential of satellite radar interferometry and feature tracking for monitoring flow rates of Himalayan glaciers

The stagnation of Himalayan glaciers will only be exacerbated by a warming climate, increasing hazards associated with glacier lake development. Frontal recession is a poor indicator of negative balance in this region where most glaciers are debris covered, but reduced glacier flow rates provide a valuable warning of stagnation. In-situ measurements of velocity are understandably rare and remote sensing measurements have hitherto been limited by severe topography and lack of suitable data. In this paper, we compare two techniques to measure flow rates in Everest region glaciers. Firstly, satellite radar interferometry (SRI) has recently become more viable because of improved digital elevation data and precision orbit models, which now allow snapshots of velocities to be retrieved from ERS imagery in the mid 1990s without ground control information. Secondly, feature tracking between Synthetic Aperture Radar (SAR) images is found to be successful, even on the small glaciers of this region, because of the patterning associated with debris cover. This technique extends through the 1990s and beyond, the time frame over which velocities can be measured. Using the examples of the Khumbu glacier flowing south into Nepal and the Kangshung glacier flowing East into China from the Everest massif, we compare SRI and SAR feature tracking in ERS images, consider the errors in these techniques, and examine the potential to monitor flow rates across the region. The techniques are found to be highly complementary and appropriate for monitoring flow rates in this region, suggesting that routine monitoring of Himalayan glaciers is possible. The Kangshung glacier appears to sustain relatively high flow rates of up to 10 cm/day which, along with a convex topographic profile, indicates a healthy glacier not exhibiting significant negative mass balance. In comparison, much of the Khumbu glacier appears to have stagnated, the surface velocity at Gorak Shep is only 20% of a field-based measurement in the 1970's, and the concave surface profile confirms the downwasting associated with a debris-covered glacier in recession.

Characteristics of Khumbu Glacier, Nepal Himalaya: recent change in the debris-covered area

AbstractUsing satellite data, the longitudinal distribution of the ablation rate and flow velocity were estimated for the ablation area, where glacier ice is covered with supraglacial debris. The ablation rate, small around Everest Base Camp (EBC) just below the equilibrium line, increased down-glacier for about 3 km, then decreased gradually toward the apparent terminus, located about 10 km from EBC. The velocity decreased almost linearly from EBC to the terminus. The results allowed estimation of the recent change of ice thickness using the continuity equation. The glacier has thinned recently in the ablation area. The rate of thinning was large near EBC, where the surface is either bare ice or covered with very thin debris cover, but relatively small and rather uniform at lower sites, where the supraglacial debris layer was thick and supraglacial lakes and ice cliffs predominate. The general pattern is compatible with field observations.

Characteristics of atmospheric particles (airborne and in snowpack)at Khumbu glacier in the Nepalese Himalayas (EV-K2-CNR project)

Atmospheric aerosol samplings taken at 5000 and 5650 m asl in a remote region of Nepalese Himalayas were determined by scanning electron microscopy for size distribution and by X-ray spectrometry for elemental composition. While the particles at both sites showed the same size distribution and elemental composition, many also had a high incidence of Si and Al, a fact suggesting their crustal origin. A stratigraphic study of a section of a high altitude glacier was then made to gain information concerning prevailing air mass trajectories and removal episodes. Periodicities in the chemical composition of melted ice and in residual microparticles are also reported.

Pedogenic gradients for iron and aluminum accumulation and phosphorus depletion in arctic and alpine soils as a function of time and climate

Pedogenically significant chemical extracts of iron, aluminum, and phosphorus are related mainly to time in chronosequences from Baffin Island (Canadian Arctic), the alpine Sierra Nevada and Wind River Range (western United States), the alpine Khumbu Glacier area (Himalaya, Nepal), and the alpine Southern Alps (New Zealand). Based on the accumulation index for Fe and Al, and the depletion index for P in each region, the following ranking of chemical pedogenic development is obtained: Southern Alps > Khumbu Glacier area ≅ Wind River Range > Sierra Nevada > coastal Baffin Island > inland Baffin Island. The ranking mainly follows regional climate, with the greatest accumulation and depletion to the warmest and wettest environment, and the least accumulation and depletion in the coldest and driest environment.