Skip to main content

Main menu

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
    • Front Matter Portal
    • Journal Club
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Home
Home
  • Log in
  • My Cart

Advanced Search

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
    • Front Matter Portal
    • Journal Club
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
Research Article

Aspect controls the survival of ice cliffs on debris-covered glaciers

View ORCID ProfilePascal Buri and View ORCID ProfileFrancesca Pellicciotti
  1. aInstitute of Environmental Engineering, Hydrology and Water Resources Management Group, ETH Zürich, 8093 Zürich, Switzerland;
  2. bFaculty of Engineering and Environment, Department of Geography, Northumbria University, Newcastle upon Tyne, NE18ST, United Kingdom

See allHide authors and affiliations

PNAS April 24, 2018 115 (17) 4369-4374; first published April 9, 2018; https://doi.org/10.1073/pnas.1713892115
Pascal Buri
aInstitute of Environmental Engineering, Hydrology and Water Resources Management Group, ETH Zürich, 8093 Zürich, Switzerland;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Pascal Buri
  • For correspondence: buri@ifu.baug.ethz.ch
Francesca Pellicciotti
bFaculty of Engineering and Environment, Department of Geography, Northumbria University, Newcastle upon Tyne, NE18ST, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Francesca Pellicciotti
  1. Edited by Andrea Rinaldo, École polytechnique fédérale, Lausanne, Switzerland, and approved March 9, 2018 (received for review August 6, 2017)

  • Article
  • Figures & SI
  • Info & Metrics
  • PDF
Loading

Significance

Glaciers in High Mountain Asia (HMA) are important water sources for millions of people downstream. Ice cliffs on debris-covered glaciers act as hot spots for melt and may explain anomalously high glacier mass losses in HMA, but their temporal evolution remains unknown, hindering sound parameterizations of these features in glacier models. We simulate the evolution of cliff systems with different aspects, show that south-facing ice cliffs disappear within a few weeks in the Northern Hemisphere, and explain the processes driving this. Cliffs that persist melt 10 times faster than the surrounding glacier surfaces. These findings provide a basis for understanding the surface evolution of debris-covered glaciers with implications for their dynamics, mass balance, and hydrology.

Abstract

Supraglacial ice cliffs exist on debris-covered glaciers worldwide, but despite their importance as melt hot spots, their life cycle is little understood. Early field observations had advanced a hypothesis of survival of north-facing and disappearance of south-facing cliffs, which is central for predicting the contribution of cliffs to total glacier mass losses. Their role as windows of energy transfer suggests they may explain the anomalously high mass losses of debris-covered glaciers in High Mountain Asia (HMA) despite the insulating debris, currently at the center of a debated controversy. We use a 3D model of cliff evolution coupled to very high-resolution topographic data to demonstrate that ice cliffs facing south (in the Northern Hemisphere) disappear within a few months due to enhanced solar radiation receipts and that aspect is the key control on cliffs evolution. We reproduce continuous flattening of south-facing cliffs, a result of their vertical gradient of incoming solar radiation and sky view factor. Our results establish that only north-facing cliffs are recurrent features and thus stable contributors to the melting of debris-covered glaciers. Satellite observations and mass balance modeling confirms that few south-facing cliffs of small size exist on the glaciers of Langtang, and their contribution to the glacier volume losses is very small (∼1%). This has major implications for the mass balance of HMA debris-covered glaciers as it provides the basis for new parameterizations of cliff evolution and distribution to constrain volume losses in a region where glaciers are highly relevant as water sources for millions of people.

  • debris cover
  • ice cliffs
  • glaciers
  • High Mountain Asia
  • modeling

Many glacier tongues in High Mountain Asia (HMA) are heavily debris-covered (1, 2). Despite the insulating effect of a mantle composed by rock debris on the underlying ice (3, 4), large-scale, satellite-based studies have suggested that thinning rates of debris-covered glaciers are comparable to those of clean ice glaciers (5, 6). Although recent studies at the catchment and glacier scale do not support analogous thinning (7, 8), it has by now been established that strong local increases in glacier ablation are associated with supraglacial ponds and cliffs (9⇓⇓–12). Cliffs forming on the surface of debris-covered glaciers contribute to the glacier mass balance through enhanced melt rates but also affect glacier dynamics, and knowledge about their life cycle and distribution is important to predict future evolution of debris-covered glaciers (13). The understanding of processes acting at the scale of single cliffs has been dramatically improved recently through modeling approaches that have simulated energy fluxes and melt (11, 14) and estimated volume losses (15) of single cliffs. The rate at which cliffs can affect glacier mass balance and dynamics depends on their distribution and persistence in time, but how cliffs form, evolve, and decline is not yet understood, precluding a holistic understanding of their role on longer term mass balance patterns beyond the few observations over a melt season. A hypothesis of persisting north-facing and disappearing south-facing cliffs has been first proposed more than one decade ago (16) based on observations and conceptual assumptions on the importance of solar radiation on ice cliff melt (17, 18). The hypothesis seems to be supported by inventories of cliff distribution from satellite observations of single or selected glaciers in the Khumbu region (Nepalese Himalaya) (12, 19). Conceptual intuition, supported by sparse observational evidence, has postulated that cliff faces oriented to the south are reburied rapidly and do not persist over debris-covered glaciers, independently of glacier flow direction. No study, however, has been able so far to explain the absence of south-facing cliffs on debris-covered glaciers.

Backwasting of South-Facing Cliffs

Here, we simulate the evolution of south-facing ice cliffs to understand the effect of enhanced solar radiation compared with observed north-facing cliffs. Our aim is to establish whether south-facing supraglacial cliffs persist beyond the length of a melt season, as observed northerly facing cliffs do, or if they disappear more rapidly, and to identify the causes for their behavior. To test this, we run a 3D numerical model of cliff backwasting that was able to reproduce the evolution of north-facing cliffs (14). We force the model with hourly meteorological data from an on-glacier automatic weather station (AWS) (20) and initialize it with a digital elevation model (DEM) (21) of submeter resolution over the debris-covered Lirung Glacier (Nepalese Himalaya; Fig. 1A).

Fig. 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 1.

Observed cliffs on Lirung Glacier, Langtang Valley, Nepalese Himalaya. (A) Lirung Glacier with debris-covered tongue (orange) and accumulation area (violet). (B) Lirung Glacier surface around cliffs 1 and 2 (marked by colors indicating their observed aspect). Encircled crosses denote positions where terrestrial images (C and D) were taken; triangle shows the location of the AWS. (C) Cliff 2 photographed from the location indicated in B (Top) in May 2013, and rotated cliff system shown as a 3D elevation model (Bottom). (D) Cliff 1 (Left) as observed in a photo (taken from the location shown in B in May 2013), and rotated cliff system shown as a 3D elevation model (Right). Background images: (A) ALOS Orthoimage December 2010 and ASTER GDEM2 hillshade; (B) unmanned aerial vehicle (UAV) Orthoimage May 2013 and UAV DEM May 2013 hillshade; (C and D) Picture E. Miles May 2013 and partly rotated UAV DEM May 2013.

Initial conditions for our simulations were created by rotating north-facing ice cliff topographies as observed on Lirung Glacier (Fig. 1B) toward the south, including the surrounding glacier surface and ponds (Fig. 1C). Hence, the artificially derived south-facing cliffs were embedded into a realistic cliff-glacier topography and therefore directly comparable to the north-facing cliffs in terms of size, shape, and surrounding topography. We applied a dynamic, physically based backwasting model (14) on the two rotated cliffs over one ablation season (May to October 2013). Cliff melt is derived from distributed surface energy balance calculations and shapes the cliff surface by biweekly geometry updates. Melt at water-contact zones is enhanced to account for thermo-erosion by adjacent supraglacial ponds (10, 11). Depending on the slope and the amount of debris cells seen at the cliff margins, the cliffs can expand or shrink (because of reburial by debris).

We simulate continuous shrinkage of the south-facing cliffs, resulting in a significant reduction in extent after just a few weeks already (Fig. 2A). This is a striking difference compared with the evolution of the original north-facing cliffs (observed in the field and confirmed by our simulations) (14), shown in the background of Fig. 2A, which backwaste maintaining a self-similar geometry that allows the cliffs to persist until the end of the ablation season. The reason for the rapid shrinking of the south-facing cliffs is the progressive flattening of their surface (Fig. 2B), which allows reburial by debris. The complete reburial of the debris-free cliff areas occur after less than 3 (cliff 1) to 5 months (cliff 2; SI Appendix, Table S2). Even when the cliff is not entirely reburied, large sections of its surface disappear, reducing consistently the area available for melt (Fig. 3C). In contrast, the north-facing cliffs show stable profiles backwasting with a constant slope (cliff 2) or only minimal regrading (cliff 1; Fig. 2B).

Fig. 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 2.

Simulated outlines and elevation profiles of south-facing cliffs. (A) Cliffs 1 (Top) and 2 (Bottom) outlines simulated by the model with biweekly geometry updates (yellow to green lines). For comparison, also the observed shapes of north-facing cliffs are shown (light and dark gray polygons), rotated toward south for consistency. [The model was able to simulate the evolution of the original north-facing cliffs (14).] Dashed line indicates profile; thick circles indicate the debris–ice transitions. (B) Elevation profiles of rotated cliffs 1 (Top) and 2 (Bottom) as simulated with biweekly geometry updates (yellow to green lines). Profiles of observed north-facing cliffs are also shown (light and dark gray areas), rotated toward the south for consistency. Circles indicate debris–ice transitions of modeled (yellow to green) and observed (thick black) cliff profiles. The last of the colored lines (darkest green) indicates the last cliff profile before the cliffs disappear. None of the two cliffs survives for the duration of the ablation season, disappearing after day 85 (cliff 1) and day 141 (cliff 2). Days are counted from the start of the simulations, on May 19, 2013.

Fig. 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 3.

Modeled surface energy fluxes for cliff 2 rotated to various aspects. (A) Diurnal cycle (May–October 2013) of direct shortwave radiation receipt averaged in space over cliff 2 rotated to eight different aspects by increments of 45°. (B) Diurnal cycle (May–October 2013) of melt energy averaged for cliff 2 rotated to eight different aspects by increments of 45°. (C) Cliff persistence per aspect (angular scale in days, counted from the start of the simulations, on May 19, 2013); black lines indicate the time when more than 50% of the initially inclined area has disappeared (with ΔA indicated, providing the percentage of area that has disappeared at that time); red lines indicate the range of directions for which cliffs never reach that threshold (i.e., never lose more than 50% of their inclined area). Cliffs with aspects indicated in white (W, NW, N, NE) persisted for the entire season. In the background, the average daily sum of simulated incoming solar radiation per aspect is shown (blue to yellow).

Radiative Forcing at the Cliff Surface

To understand what controls the simulated cliffs’ evolution, we rotated cliffs 1 and 2 together with their surrounding topography by increments of 45° from the north into eight additional directions and modeled the seasonal surface energy balances. We then calculated diurnal cycles of the spatially averaged energy fluxes for the rotated cliff surfaces (Fig. 3 A and B and SI Appendix, Fig. S5) and spatial totals of energy fluxes and melt energy (Fig. 4C and SI Appendix, Figs. S7 C and D and S8).

Fig. 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 4.

Distributed energy fluxes modeled over cliff 2 rotated to various aspects. (A) Distributed incoming shortwave radiation averaged over melt season (May–October 2013) rotated to eight different aspects by increments of 45° (indicated by label at crest of each cliff), shown for hours 8 (Left), 12 (Middle), and 16 (Right) of the day, respectively. (B) Map of melt energy per pixel averaged over the melt season (May–October 2013) calculated over cliff 2 rotated to eight different aspects by increments of 45°, shown for hours 8 (Left), 12 (Middle), and 16 (Right) of the day, respectively. (C) Distributed daily sum (averaged May–October 2013) of incoming shortwave radiation. (D) Map of daily sums (averaged May–October 2013) of melt energy.

The longwave radiation component, composed of radiation emitted by the debris surfaces around the cliff and of the longwave radiation emitted by the atmosphere, shows no aspect-related differences in amount and timing (SI Appendix, Figs. S5 E and F and S8 A and B). This is not surprising as these fluxes depend on the surface (debris) and air (atmosphere) temperatures, which have no obvious dependence on aspect, and on the local topographical horizons (which are approximately constant for all directions). In contrast, a very high aspect-dependence is evident for the simulated shortwave radiation and its direct component in particular (Figs. 3 A and C and 4 A and C). Differences between directions are evident in both the timing and total amount of solar energy received. East-facing cliffs receive direct solar radiation earliest in the day, followed by south- and west-oriented slopes (Figs. 3A and 4A). The lowest amounts are received by cliffs with aspects in the range north to southwest (Fig. 3 A and C). East- and southeast-facing cliffs receive the highest direct solar radiation [up to 67% more than the original cliff and exceeding by 3.6 times the energy input at the northwest-facing cliff (SI Appendix, Table S1)], followed by south-facing ones. These cliffs do not survive the duration of the ablation season but disappear or undergo a substantial loss in area (Fig. 3C and SI Appendix, Table S2). The apparently anomalous behavior of south- and southwest-facing cliffs, which receive as little radiation as those with a prevalent northerly aspect, is likely due to the presence of cloud cover in the afternoon. During the ablation season, which coincides with the monsoon in this region, in the afternoon, when the south-facing cliffs are theoretically exposed to high solar radiation receipts, thick clouds and rain prevail with regularity and prevent high solar radiation incomes in the Langtang Valley (22, 23). This decreases the solar radiation receipt of southwesterly aspects considerably (18) and therefore dampens the all-year average of incoming shortwave energy (Fig. 3). The daily cycle and spatial patterns of melt energy closely reflect those of the solar radiation inputs, with the highest amount of melt energy for east- and southeast-facing cliffs (Figs. 3B and 4B). As a result of the energy forcing, cliffs with aspect in the range east to southwest do not survive, while cliffs facing northeast to west do (Fig. 3C and SI Appendix, Table S2).

Spatial variability in solar radiation and melt energy is high over a single cliff (coefficients of variation for direct shortwave radiation up to 238% at the west-oriented surface of cliff 2; SI Appendix, Table S3). Solar radiation is highest at the top of the cliff, and this effect is stronger at noon because of the high sun angle (Fig. 4A). The top sections of the cliffs receive also the highest amount of atmospheric longwave radiation (SI Appendix, Fig. S8B), thus amplifying the solar radiation control. This cannot be counterbalanced by the longwave flux emitted by the debris surface surrounding the slopes, which is highest at the cliff margins (SI Appendix, Fig. S8A). Total melt energy results from the interactions of these spatially variable fluxes and their temporal variability: It is highest at the top of the cliff (Fig. 4B) for most aspects and decreases toward the cliffs bottom. This energy gradient is small (with minimum differences between the energy at the top and bottom of the cliff) for cliffs with northwest and western aspects (Fig. 4 B and D and SI Appendix, Table S3). These are those that survive (Fig. 3C) because a rather uniform distribution of solar radiation and melt energy allows their backwasting and maintenance of a constant steep slope, rather than downwasting and reburial by debris. The flux of energy emitted by the surrounding debris and received by the cliff margins is not high enough to counterbalance the atmospheric fluxes of shortwave and longwave radiation at south-oriented cliffs. High receipts of solar radiation at the top of these cliffs cause a progressive flattening. The cliff flattening, controlled by the sky view factor and hence the amount of sky to which the cliff sections are exposed, is thus strongly aspect-dependent. Longitudinal profiles of progressively higher solar radiation amounts from base to top will have a much stronger vertical gradient for those aspects that receive much higher solar radiation in the morning hours (northeast to southeast). The upper part of cliff 2 shows a 20 to 30% higher sky view factor compared with the base zone (SI Appendix, Fig. S9C). The reduction in sky-openness toward the cliff bottom is the combined result of the topography in front of the cliff face and the steep slopes at the cliff bottom. The combination of a very high shortwave radiation income together with a decreasing sky view factor toward the cliff base causes the cliffs with southerly to easterly aspect to flatten progressively over time, as the upper section recedes at much higher rates than the lower parts, until they reach a slope that can be reburied by debris.

Discussion

Our model results show that south-facing supraglacial ice cliffs progressively shrink and disappear within a few weeks. We thus provide an explanation for previous observations and conceptual suggestions that (in the Northern Hemisphere) cliffs with a southern aspect are not part of the cliff population on glacier surfaces, as they do not persist on time scales relevant for glacier mass balance considerations. This narrows the knowledge gap concerning distribution and evolution of cliffs as the population of cliff systems can be reduced to northerly to westerly facing ones. We can explain this distribution with the enhanced solar radiation received by the cliffs with southern aspects. Southeast- and northwest-oriented cliffs are likely the extremes of cliff life expectancy, as exposure to solar radiation and shadowing, respectively, are highest for these aspects.

We have also further established that ice cliffs are melt hot spots that efficiently convey a large amount of atmospheric energy into the glacier ice. The daily melt rates of the two ice cliffs vary between 4.6 (for northwest orientation) and 6.3 cm (for east and south-eastern orientations; SI Appendix, Table S1) and exceed the observed daily subdebris melt of 0.5 cm on Lirung Glacier (21) for the same period by about 10 times. However, starting from very high melt rates for all cliff orientations (and for the predominantly north- and predominantly south-facing cliffs), their behavior diverges significantly over the course of the melt season: On southerly facing cliffs, the spatial distribution of the energy fluxes leads to the progressive flattening and disappearance of the cliffs (Fig. 2), while on northerly facing cliffs, the distinct interaction of cliff topography and energy distribution maintains self-consistent, persistent cliffs.

We are able to reproduce the flattening of southerly facing cliffs induced by much higher direct solar radiation, compared with northerly oriented cliffs, and the increase of the shortwave radiation-relevant sky view factor from cliff base to crest. The increasing debris view factor toward the cliff base (along a vertical gradient) and boundary zones (along a horizontal gradient from the cliff center, SI Appendix, Fig. S9D) results in a higher longwave radiation receipt from the surrounding debris at these cliff zones. This, however, is not able to counterbalance the extremely high solar radiation receipt of southerly aspects (as it is the case for cliff slopes facing north) (11, 14, 20). Importantly, we have shown that the effect of adjacent ponds (which act on cliffs through enhanced melt through thermo-erosion at the low-lying cliff–pond contact zone) is not sufficient to maintain southerly facing cliffs steep and thus allow their persistence (Fig. 2B), as they are able to do for northerly oriented cliffs (14). We show that there is a range of cliff aspects that determine their disappearance as a result of energy flux interaction and a range of aspects within which cliffs over monsoon-dominated central Himalayan glaciers will survive over the melting season: Aspects from northeast to west are associated with cliff persistence, and those from east to southwest are associated with progressive flattening and disappearance (Fig. 3C).

To test our results, we manually mapped all cliffs and ponds from UAV images in May 2014 as well as from a terrestrial photogrammetry survey carried out in October 2014 on Lirung Glacier (15). Since the UAV surveys cover only a portion of the glacier, we used a SPOT6 orthoimage from April 2014 (very close to the UAV survey of May 2014) for both Lirung Glacier and Langtang Glacier, the largest glacier in the valley (SI Appendix, Fig. S10 and Methods). For Langtang Glacier, we additionally use UAV imagery from May 2014 and October 2015 to map cliffs and lakes at very high resolution. There are no southerly facing cliffs on Lirung Glacier in the portion covered by the UAV survey in either May or October. There are a total of four south-facing cliffs on the entire Lirung Glacier in April 2014 (from the SPOT6 image), three on Langtang Glacier on the portion covered by the UAV in May 2014, and nine in total over the entire glacier in April 2014 (SI Appendix, Figs. S12 and S14 and Table S6). All south-facing cliffs are very small, covering 0.01% of the entire Lirung Glacier in April 2014, 2.93% of the portion of Langtang covered by the UAV in May 2014, and 0.07% of the entire Langtang Glacier in April 2014. While cliffs cover a total of 1.29% of the debris-covered area, only 5.14% of this total cliff area is made of southerly facing cliffs (SI Appendix, Table S6), with two orders of magnitude difference in the extension of southerly facing cliffs compared with the entire population (SI Appendix, Table S6).

We also run the cliff model on all of the cliffs on the two glaciers (Methods and SI Appendix, Section 5) to estimate their total contribution to the mass losses of the two glaciers for the period between May and October 2014. Cliffs are major contributors to total glacier mass losses (with contributions of 36.43% and 19.84% for Lirung and Langtang Glaciers, respectively, relative to the debris-covered glacier area; SI Appendix, Table S6). Southerly facing cliffs, however, contribute only to a very small percentage of these mass losses (1.2% on Langtang Glacier and 0% on Lirung Glacier, as all southerly oriented cliffs disappear with the first geometry update; SI Appendix, Table S6).

The glacier scale observational evidence and large-scale modeling confirm the main findings of our modeling experiment. Southerly facing cliffs are very few on two of the main glaciers of the Langtang Valley, both at the beginning and at the end of the ablation season (SI Appendix, Figs. S11 and S12), suggesting that indeed southerly facing cliffs do not form part of the population of stable cliffs on the glaciers of the Langtang catchment. The two glaciers differ in area, dynamics, and elevation ranges, and while Lirung has a quasi-stagnant tongue, Langtang Glacier is much larger (40.2 km2) and more active (7), suggesting that our results are largely independent of flow dynamics, at least within the range of velocities of the Langtang Valley glaciers (7, 24). It is not clear, however, how the south-facing cliffs form, for lack of a general understanding on the formation of cliffs in general. Our work has established how cliffs evolve and decay and that the solar radiation received by a cliff and the shadowing of steep cliff surfaces is the first-order control of cliff melt, evolution, and distribution. However, while radiation seems to ultimately control the evolution and disappearance of supraglacial ice cliffs, their appearance and the mechanisms controlling their formation are still largely unknown. Different hypotheses have been advanced, from subsurface developments such as collapsing of empty melt water channels close to the surface to surface changes induced by glacier dynamics or subdebris melt (16, 25, 26), but none has been demonstrated conclusively. The picture is complicated by the fact that little is known on the distribution and characteristics of debris cover worldwide and in HMA in particular. Initial observational and satellite evidence suggests that debris characteristics (thickness and spatial distribution) might vary substantially along the extreme climatic and geomorphological gradient of HMA. And yet, cliffs and ponds do appear to form on most of the region’s debris-covered glaciers, from the stagnant tongues of central Himalayan glaciers to the much more active, winter accumulation-type Karakoram glaciers. This is an important field of future investigation that will need to be addressed to understand debris-covered glaciers’ mass balance and dynamics. It can substantially benefit from the availability of new high-resolution satellite images and very high-resolution UAV surveys, as we have shown that high-resolution topographical information of both the cliff and surrounding glacier surface is crucial to understand and correctly represent cliff backwasting patterns.

Materials and Methods

We mapped two supraglacial ice cliffs on the debris-covered tongue of Lirung Glacier (Langtang Valley, Nepalese Himalaya) using a high-resolution orthoimage and DEM, which were derived from a UAV survey in May 2013 (21). No south-facing cliffs were observed on Lirung Glacier (nor on the other glaciers of the Langtang valley). Therefore, we rotated the two cliffs including their surrounding topography and ponds (within 100 m in xy-direction) by applying a 2D-matrix rotation around a common center coordinate. The rotation angle is defined as the deviation between the observed mean cliff aspect and the target direction in degrees. For our simulations, we selected two observed cliffs of different size (one relatively large and one relatively small), aspect (northeast and northwest), and bottom configuration (in contact with a supraglacial pond and with no water contact). Both were located within 100 m from an on-glacier AWS, which allowed forcing the cliff energy balance and backwasting model with local, high-resolution and accurate meteorological input.

A physically based, dynamic 3D-backwasting model (14), which has previously been tested for four cliffs (of which two are investigated in this study) on the same glacier and for the same period, allowed us to test the behavior of the south-facing cliffs generated by rotation of the two original cliffs. The model has been validated for the two original cliffs with multiple independent datasets (14), lending confidence to its use for this experiment. For this study, we further improved the model algorithm for more stable and computationally efficient simulations (SI Appendix, Section 3). The use of a high-resolution DEM for initial conditions and hourly meteorological data recorded on-glacier allow the model to calculate radiation and shading at the cliff surface with a very high level of detail (11). Simulated melt from calculation of the cliff surface energy balance (SI Appendix, Section 3.A) was accumulated for every cliff cell over biweekly intervals, after which the cliff geometry was updated accordingly (14) (SI Appendix, Section 3.D). Enhanced melt rates were applied to cliff sections in direct contact with ponded water (14), accounting for thermo-erosion (10) (SI Appendix, Section 3.J). The model algorithm also considered expansion and shrinkage of marginal cliff zones based on slope and debris-view thresholds as described in ref. 14 and in SI Appendix (Sections 3.G and H).

To provide the context of our modeling experiment, we mapped supraglacial cliffs and ponds on both the Lirung Glacier and the much larger Langtang Glacier, the largest and most remote glacier in the Langtang catchment (ref. 7 and SI Appendix, Fig. S10). We have used a UAV survey from May 2014 (orthoimage and DEM with 0.1 m and 0.2 m spatial resolution, respectively) and SPOT6 imagery from April 2014 (orthoimage and DEM with 1.5 m and 3 m spatial resolution, respectively) to delineate supraglacial ice cliffs and ponds and to derive their initial topographies. Mapping was carried out manually, based on visual interpretation using the high-resolution orthoimages and topography (slope) information (SI Appendix, Section 5.A.1) (7). We used these inventories to determine the distribution of southerly facing cliffs within the total cliff distribution on both glaciers. We then applied the 3D ice cliff ablation and backwasting model (SI Appendix, Section 3) to all of the cliffs on the two glaciers to calculate the volume losses associated with all cliffs and specifically south-facing cliffs, respectively, over one ablation season. We use a fully distributed, physically based glacio-hydrological model [Topographic Kinematic Wave Approximation and Integration ETH Zurich (TOPKAPI-ETH)] run over the same period and the same spatial domain to calculate the mass losses of the two glaciers (SI Appendix, Section 5.A.2). The models are run with meteorological input data from AWSs on-glacier and in the valley, extrapolated to each single cliff location with local lapse rates (SI Appendix, Section 5.A.2). Radiative fluxes are modeled (SI Appendix, Sections 3.A and B) and cliff geometries are updated (SI Appendix, Section 3.D) two times during the melt season.

Acknowledgments

We thank two anonymous reviewers, who provided very thorough and constructive reviews that contributed to improving the manuscript. This study was funded by the SNF (Swiss National Science Foundation) project UNCOMUN (“Understanding Contrasts in High Mountain Hydrology in Asia,” Grant No. 146761).

Footnotes

  • ↵1To whom correspondence should be addressed. Email: buri{at}ifu.baug.ethz.ch.
  • Author contributions: F.P. designed research; P.B. performed research; P.B. and F.P. analyzed data; and P.B. and F.P. wrote the paper.

  • The authors declare no conflict of interest.

  • This article is a PNAS Direct Submission.

  • This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1713892115/-/DCSupplemental.

Published under the PNAS license.

References

  1. ↵
    1. Scherler D,
    2. Bookhagen B,
    3. Strecker MR
    (2011) Spatially variable response of Himalayan glaciers to climate change affected by debris cover. Nat Geosci 4:156–159.
    OpenUrlCrossRef
  2. ↵
    1. Bolch T, et al.
    (2012) The state and fate of Himalayan glaciers. Science 336:310–314.
    OpenUrlAbstract/FREE Full Text
  3. ↵
    1. Gunnar Ø
    (1959) Ice melting under a thin layer of moraine, and the existence of ice cores in moraine ridges. Geogr Ann 41:228–230.
    OpenUrl
  4. ↵
    1. Evatt GW, et al.
    (2015) Glacial melt under a porous debris layer. J Glaciol 61:825–836.
    OpenUrl
  5. ↵
    1. Gardelle J,
    2. Berthier E,
    3. Arnaud Y
    (2012) Slight mass gain of Karakoram glaciers in the early twenty-first century. Nat Geosci 5:322–325.
    OpenUrlCrossRef
  6. ↵
    1. Kääb A,
    2. Berthier E,
    3. Nuth C,
    4. Gardelle J,
    5. Arnaud Y
    (2012) Contrasting patterns of early twenty-first-century glacier mass change in the Himalayas. Nature 488:495–498.
    OpenUrlCrossRefPubMed
  7. ↵
    1. Ragettli S,
    2. Bolch T,
    3. Pellicciotti F
    (2016) Heterogeneous glacier thinning patterns over the last 40 years in Langtang Himal, Nepal. Cryosphere 10:2075–2097.
    OpenUrl
  8. ↵
    1. Vincent C
    (2016) Reduced melt on debris-covered glaciers: Investigations from Changri Nup glacier, Nepal. Cryosphere 10:1845–1858.
    OpenUrl
  9. ↵
    1. Sakai A,
    2. Takeuchi N,
    3. Fujita K,
    4. Nakawo M
    (2000) Role of supraglacial ponds in the ablation process of a debris-covered glacier in the Nepal Himalayas. Debris-Covered Glaciers, Nakawo M, Raymond CF, Fountain A (International Association of Hydrological Sciences, Oxford, UK), Vol 265, pp 119–132.
  10. ↵
    1. Miles ES, et al.
    (2016) Refined energy-balance modelling of a supraglacial pond, Langtang Khola, Nepal. Ann Glaciol 57:29–40.
    OpenUrlCrossRef
  11. ↵
    1. Buri P,
    2. Pellicciotti F,
    3. Steiner JF,
    4. Miles ES,
    5. Immerzeel WW
    (2016a) A grid-based model of backwasting of supraglacial ice cliffs on debris-covered glaciers. Ann Glaciol 57:199–211.
    OpenUrlCrossRef
  12. ↵
    1. Thompson S,
    2. Benn DI,
    3. Mertes J,
    4. Luckman A
    (2016) Stagnation and mass loss on a Himalayan debris-covered glacier: Processes, patterns and rates. J Glaciol 62:467–485.
    OpenUrl
  13. ↵
    1. Rowan AV,
    2. Egholm DL,
    3. Quincey DJ,
    4. Glasser NF
    (2015) Modelling the feed-backs between mass balance, ice flow and debris transport to predict the response to climate change of debris-covered glaciers in the Himalaya. Earth Planet Sci Lett 430:427–438.
    OpenUrl
  14. ↵
    1. Buri P, et al.
    (2016b) A physically based 3-D model of ice cliff evolution over debris-covered glaciers. J Geophys Res Earth Surf 121:2471–2493.
    OpenUrl
  15. ↵
    1. Brun F, et al.
    (2016) Quantifying volume loss from ice cliffs on debris-covered glaciers using high-resolution terrestrial and aerial photogrammetry. J Glaciol 62:684–695.
    OpenUrl
  16. ↵
    1. Sakai A,
    2. Nakawo M,
    3. Fujita K
    (2002) Distribution characteristics and energy balance of ice cliffs on debris-covered glaciers, Nepal Himalaya. Arct Antarct Alp Res 34:12–19.
    OpenUrl
  17. ↵
    1. Inoue J,
    2. Yoshida M
    (1980) Ablation and heat exchange over the Khumbu glacier. J Jpn Soc Snow Ice 41:26–33.
    OpenUrl
  18. ↵
    1. Sakai A,
    2. Nakawo M,
    3. Fujita K
    (1998) Melt rate of ice cliffs on the Lirung glacier, Nepal Himalayas 1996. Bull Glacier Res 16:57–66.
    OpenUrl
  19. ↵
    1. Scott Watson C,
    2. Quincey DJ,
    3. Carrivick JL,
    4. Smith MW
    (2017) Ice cliff dynamics in the Everest region of the Central Himalaya. Geomorphology 278:238–251.
    OpenUrl
  20. ↵
    1. Steiner JF, et al.
    (2015) Modeling ice cliff backwasting on a debris covered glacier in the Nepalese Himalayas. J Glaciol 61:889–907.
    OpenUrl
  21. ↵
    1. Immerzeel WW, et al.
    (2014) High-resolution monitoring of Himalayan glacier dynamics using unmanned aerial vehicles. Remote Sens Environ 150:93–103.
    OpenUrlCrossRef
  22. ↵
    1. Fujita K,
    2. Sakai A,
    3. Chhetri TB
    (1997) Meteorological observation in Langtang Valley, Nepal Himalayas 1996. Bull Glacier Res 15:71–78.
    OpenUrl
  23. ↵
    1. Immerzeel WW,
    2. Petersen L,
    3. Ragettli S,
    4. Pellicciotti F
    (2014) The importance of observed gradients of air temperature and precipitation for modeling runoff from a glacierized watershed in the Nepalese Himalayas. Water Resour Res 50:2212–2226.
    OpenUrl
  24. ↵
    1. Pellicciotti F, et al.
    (2015) Mass-balance changes of the debris-covered glaciers in the Langtang Himal, Nepal, 1974–99. J Glaciol 61:373–386.
    OpenUrl
  25. ↵
    1. Kirkbride MP
    (1993) The temporal significance of transitions from melting to calving termini at glaciers in the central Southern Alps of New Zealand. Holocene 3:232–240.
    OpenUrlCrossRef
  26. ↵
    1. Benn DI, et al.
    (2012) Response of debris-covered glaciers in the Mount Everest region to recent warming, and implications for outburst flood hazards. Earth Sci Rev 114:156–174.
    OpenUrl
PreviousNext
Back to top
Article Alerts
Email Article

Thank you for your interest in spreading the word on PNAS.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Aspect controls the survival of ice cliffs on debris-covered glaciers
(Your Name) has sent you a message from PNAS
(Your Name) thought you would like to see the PNAS web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Aspect controls the survival of ice cliffs on debris-covered glaciers
Pascal Buri, Francesca Pellicciotti
Proceedings of the National Academy of Sciences Apr 2018, 115 (17) 4369-4374; DOI: 10.1073/pnas.1713892115

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Request Permissions
Share
Aspect controls the survival of ice cliffs on debris-covered glaciers
Pascal Buri, Francesca Pellicciotti
Proceedings of the National Academy of Sciences Apr 2018, 115 (17) 4369-4374; DOI: 10.1073/pnas.1713892115
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Mendeley logo Mendeley

Article Classifications

  • Physical Sciences
  • Environmental Sciences
Proceedings of the National Academy of Sciences: 115 (17)
Table of Contents

Submit

Sign up for Article Alerts

Jump to section

  • Article
    • Abstract
    • Backwasting of South-Facing Cliffs
    • Radiative Forcing at the Cliff Surface
    • Discussion
    • Materials and Methods
    • Acknowledgments
    • Footnotes
    • References
  • Figures & SI
  • Info & Metrics
  • PDF

You May Also be Interested in

Setting sun over a sun-baked dirt landscape
Core Concept: Popular integrated assessment climate policy models have key caveats
Better explicating the strengths and shortcomings of these models will help refine projections and improve transparency in the years ahead.
Image credit: Witsawat.S.
Model of the Amazon forest
News Feature: A sea in the Amazon
Did the Caribbean sweep into the western Amazon millions of years ago, shaping the region’s rich biodiversity?
Image credit: Tacio Cordeiro Bicudo (University of São Paulo, São Paulo, Brazil), Victor Sacek (University of São Paulo, São Paulo, Brazil), and Lucy Reading-Ikkanda (artist).
Syrian archaeological site
Journal Club: In Mesopotamia, early cities may have faltered before climate-driven collapse
Settlements 4,200 years ago may have suffered from overpopulation before drought and lower temperatures ultimately made them unsustainable.
Image credit: Andrea Ricci.
Steamboat Geyser eruption.
Eruption of Steamboat Geyser
Mara Reed and Michael Manga explore why Yellowstone's Steamboat Geyser resumed erupting in 2018.
Listen
Past PodcastsSubscribe
Birds nestling on tree branches
Parent–offspring conflict in songbird fledging
Some songbird parents might improve their own fitness by manipulating their offspring into leaving the nest early, at the cost of fledgling survival, a study finds.
Image credit: Gil Eckrich (photographer).

Similar Articles

Site Logo
Powered by HighWire
  • Submit Manuscript
  • Twitter
  • Facebook
  • RSS Feeds
  • Email Alerts

Articles

  • Current Issue
  • Special Feature Articles – Most Recent
  • List of Issues

PNAS Portals

  • Anthropology
  • Chemistry
  • Classics
  • Front Matter
  • Physics
  • Sustainability Science
  • Teaching Resources

Information

  • Authors
  • Editorial Board
  • Reviewers
  • Subscribers
  • Librarians
  • Press
  • Site Map
  • PNAS Updates
  • FAQs
  • Accessibility Statement
  • Rights & Permissions
  • About
  • Contact

Feedback    Privacy/Legal

Copyright © 2021 National Academy of Sciences. Online ISSN 1091-6490