Debris Cover Research Articles

Microbial assemblages reflect environmental heterogeneity in alpine streams.

Alpine streams are dynamic habitats harboring substantial biodiversity across small spatial extents. The diversity of alpine stream biota is largely reflective of environmental heterogeneity stemming from varying hydrological sources. Globally, alpine stream diversity is under threat as meltwater sources recede and stream conditions become increasingly homogeneous. Much attention has been devoted to macroinvertebrate diversity in alpine headwaters, yet to fully understand the breadth of climate change threats, a more thorough accounting of microbial diversity is needed. We characterized microbial diversity (specifically Bacteria and Archaea) of 13 streams in two disjunct Rocky Mountain subranges through 16S rRNA gene sequencing. Our study encompassed the spectrum of alpine stream sources (glaciers, snowfields, subterranean ice, and groundwater) and three microhabitats (ice, biofilms, and streamwater). We observed no difference in regional (γ) diversity between subranges but substantial differences in diversity among (β) stream types and microhabitats. Within-stream (α) diversity was highest in groundwater-fed springs, lowest in glacier-fed streams, and positively correlated with water temperature for both streamwater and biofilm assemblages. We identified an underappreciated alpine stream type-the icy seep-that are fed by subterranean ice, exhibit cold temperatures (summer mean <2°C), moderate bed stability, and relatively high conductivity. Icy seeps will likely be important for combatting biodiversity losses as they contain similar microbial assemblages to streams fed by surface ice yet may be buffered against climate change by insulating debris cover. Our results show that the patterns of microbial diversity support an ominous trend for alpine stream biodiversity; as meltwater sources decline, stream communities will become more diverse locally, but regional diversity will be lost. Icy seeps, however, represent a source of optimism for the future of biodiversity in these imperiled ecosystems.

Automated Glacier Extraction Index by Optimization of Red/SWIR and NIR /SWIR Ratio Index for Glacier Mapping Using Landsat Imagery

Glaciers are recognized as key indicators of climate change on account of their sensitive reaction to minute climate variations. Extracting more accurate glacier boundaries from satellite data has become increasingly popular over the past decade, particularly when glacier outlines are regarded as a basis for change assessment. Automated multispectral glacier mapping methods based on Landsat imagery are more accurate, efficient and repeatable compared with previous glacier classification methods. However, some challenges still exist in regard to shadowed areas, clouds, water, and debris cover. In this study, a new index called the automated glacier extraction index (AGEI) is proposed to reduce water and shadow classification errors and improve the mapping accuracy of debris-free glaciers using Landsat imagery. Four test areas in China were selected and the performances of four commonly used methods: Maximum-likelihood supervised classification (ML), normalized difference snow and ice index (NDSI), single-band ratios Red/SWIR, and NIR/SWIR, were compared with the AGEI. Multiple thresholds identified by inspecting the shadowed glacier areas were tested to determine an optimal threshold. The confusion matrix, sub-pixel analysis, and plot-scale validation were calculated to evaluate the accuracies of glacier maps. The overall accuracies (OAs) created by AGEI were the highest compared to the four existing automatic methods. The sub-pixel analysis revealed that AGEI was the most accurate method for classifying glacier edge mixed pixels. Plot-scale validation indicated AGEI was good at separating challenging features from glaciers and matched the actual distribution of debris-free glaciers most closely. Therefore, the AGEI with an optimal threshold can be used for mapping debris-free glaciers with high accuracy, particularly in areas with shadows and water features.

Glacier Dynamics in Changme Khangpu Basin, Sikkim Himalaya, India, between 1975 and 2016

This study provides a high resolution glacier database in the Changme Khangpu Basin (CKB) using LANDSAT 8 (2014) and Sentinel-2A image (2016), mapping of 81 glaciers that cover a 75.78 ± 1.54 km2 area. Composite maps of land surface temperature, slope and Normalized differential Snow Index have been successfully utilized in delineating near accurate debris cover boundary of glaciers. The cumulative controlling parameters of aspect, elevation, slope, and debris cover have been assessed to evaluate the nature of glacier distribution and dynamics. The local topographic settings seem to have significantly determined the glacier distribution in the CKB. Almost 20% area erstwhile under glacier cover has been lost since 1975 at an average rate of −0.453 ± 0.001 km2a−1. The recent decade (2001–2016) has witnessed a higher rate of area shrinkage (−0.665 ± 0.243 km2a−1), compared to a relatively lower rate of recession (−0.170 ± 0.536 km2a−1) between 1988 and 2001. The lower rates of glacial recession can most likely be induced regionally due to relatively cooler decadal late summer temperatures and peak in the monsoon spell. Glaciers with western and north-western aspects showed more vulnerability to area loss than the rest of the aspects. Lower altitude glaciers have receded faster than ones perched up on higher elevations. The rate of glacier area recession has been nearly twice that on clean glaciers as compared to debris-covered glaciers in the CKB.

The impact of supraglacial debris on proglacial runoff and water chemistry

Debris is known to influence the ablation, topography and hydrological systems of glaciers. This paper determines for the first time how these influences impact on bulk water routing and the proglacial runoff signal, using analyses of supraglacial and proglacial water chemistry and proglacial discharge at Miage Glacier, Italian Alps. Debris does influence the supraglacial water chemistry, but the inefficient subglacial system beneath the debris-covered zone also plays a role in increasing the ion contribution to the proglacial stream. Daily hydrographs had a lower amplitude and later discharge peak compared to clean glaciers and fewer diurnal hydrographs were found compared to similar analysis for Haut Glacier d’Arolla. We attribute these observations to the attenuating effect of the debris on ablation, smaller input streams on the debris-covered area, a less efficient subglacial system, and possible leakage into a raised sediment bed beneath the glacier. Strongly diurnal hydrographs are constrained to periods with warmer than average conditions. ‘Average’ weather conditions result in a hydrograph with reverse asymmetry. Conductivity and discharge commonly show anti-clockwise hysteresis, suggesting the more dilute, rapidly-routed melt component from the mid-glacier peaks before the discharge peak, with components from higher up-glacier and the debris-covered areas arriving later at the proglacial stream. The results of this study could lead to a greater understanding of the hydrological structure of other debris-covered glaciers, with findings highlighting the need to include the influence of the debris cover within future models of debris-covered glacier runoff.

Heterogeneous Influence of Glacier Morphology on the Mass Balance Variability in High Mountain Asia

We investigate the control of the morphological variables on the 2000–2016 glacier‐wide mass balances of 6,470 individual glaciers of High Mountain Asia. We separate the data set into 12 regions assumed to be climatically homogeneous. We find that the slope of the glacier tongue, mean glacier elevation, percentage of supraglacial debris cover, and avalanche contributing area all together explain a maximum of 48% and a minimum of 8% of the glacier‐wide mass balance variability, within a given region. The best predictors of the glacier‐wide mass balance are the slope of the glacier tongue and the mean glacier elevation for most regions, with the notable exception of the inner Tibetan Plateau. Glacier‐wide mass balances do not differ significantly between debris‐free and debris‐covered glaciers in 7 of the 12 regions analyzed. Lake‐terminating glaciers have more negative mass balances than the regional averages, the influence of lakes being stronger on small glaciers than on large glaciers.

Evolution of a debris-covered glacier in the western Himalaya during the last four decades (1971–2016): A multiparametric assessment using remote sensing and field observations

This study assesses the response of the Panchi Nala glacier, Bhaga Basin, western Himalaya to climate change by monitoring multiple glacier parameters namely, length, area, debris cover, snowline altitude (SLA), surface ice velocity (SIV), and surface elevation change (Δh). These parameters were analysed using high-resolution satellite images (Corona KH-4B and Sentinel-2 MSI) and medium-resolution ones (Landsat TM/ETM+/OLI and Terra ASTER) acquired between 1971 and 2016, along with extensive field observations. Considering only dimensional changes, the glacier appears to have low rates of recession (4.54 ± 2.56 m y−1) and deglaciation (0.11 ± 0.04% y−1) over the study period (1971–2016). However, SIV and Δh values reveal a significant slowdown (33.87%; 1992/93–2015/16), and an extensive lowering (−1.62 ± 0.06 m y−1; 2000–2013) on the glacier. The debris cover area and SLA also increased by 23.50 ± 1.74% and 287.63 m, respectively, during 1971–2016. Increase in the mean winter minimum temperature, and decrease in the winter snowfall are probably the main drivers of glacier depletion. We conclude that the significant slowdown, extensive surface lowering, development of a proglacial lake and appearance of the supraglacial lakes collectively indicate a severe glacier degeneration. This forewarns that the glacier may experience an accelerated pace of wastage in the future.

Temporal monitoring of Ghamu bar glacial lakes using remote sensing and GIS

The trends of glacial lakes formation and glacial lake outburst flooding events have been increased across Himalayan Karakorum Hindu Kush (HKH) ranges during last decade due to increase in global warming. This research is addressing the temporal monitoring of ghamu bar glacial lakes using remote sensing and GIS. Landsat images of 1990, 2000, 2010 and 2015 were used to map temporal glacial lakes using normalized difference water index (NDWI) index. The results of normalized difference water index were validated through modified normalized difference water index and field photographs. Temporal variability shows that, glacier lake area has been increase from 1990 to 2015. In 1990 total area of lake was 0.052 sq, which further increased 0.0423 in 1995 than it decreases to 0.314 in 2000 due to detached of debris cover moraine from glacier tongue and it reach 0.0846 sq.km in 2005. The area gradually increased up to 0.1296 sq.km in 2010 and it goes up to 0.157 sq.km 2015. The overall increase in area are expanding at an accelerated rate in past two decades, indicating that Darkut glacier is more vulnerable toward climate change through increase in size and volume ofghamu bar glacial lakes. There is need for vigilance in monitoring of ghamu bar glacial lake through high resolution remote sensing data and development of Geo-database enabling more details about past and future lakes behaviors toward climate change impacts.

Glacier mass balance over the central Nyainqentanglha Range during recent decades derived from remote-sensing data

ABSTRACTTo obtain information on changes in glacier mass balance in the central Nyainqentanglha Range, a comprehensive study was carried out based on digital-elevation models derived from the 1968 topographic maps, the Shuttle Radar Topography Mission DEM (2000) and TerraSAR-X/TanDEM-X (2013). Glacier area changes between 1968 and 2016 were derived from topographic maps and Landsat OLI images. This showed the area contained 715 glaciers, with an area of 1713.42 ± 51.82 km2, in 2016. Ice cover has been shrinking by 0.68 ± 0.05% a−1 since 1968. The glacier area covered by debris accounted for 11.9% of the total and decreased in the SE–NW directions. Using digital elevation model differencing and differential synthetic aperture radar interferometry, a significant mass loss of 0.46 ± 0.10 m w.e. a−1 has been recorded since 1968; mass losses accelerated from 0.42 ± 0.20 m w.e. a−1 to 0.60 ± 0.20 m w.e. a−1 between 1968–2000 and 2000–2013, with thinning noticeably greater on the debris-covered ice than the clean ice. Surface-elevation changes can be influenced by ice cliffs, as well as debris cover and land- or lake-terminating glaciers. Changes showed spatial and temporal heterogeneity and a substantial correlation with climate warming and decreased precipitation.

Early twenty-first century glacier mass losses in the Indus Basin constrained by density assumptions

Glaciers in the upper Indus supply more than half of the river water, are experiencing significant melting with a debated fate. The recent melting rate is still contained by considerable uncertainties, hindering to estimate precise glacier mass change. Here we present geodetic mass balance results for the whole Indus Basin using SRTM and ALOS 30 m elevation data, improved glacier inventory, optimized glacier surface density, and validation through in-situ differential GPS and ICESat data. Our glacier inventory and derived by improving RGI6.0 boundaries and separated into debris cover and debris free parts. The derived surface elevation changes were converted into annual mass balances using separated density assumptions (four criteria) for debris-covered ice (900 ± 60 kg m−3), debris-free ice (below 20° and 25° slopes (850 ± 60 kg m−3) and above 20° and 25° slopes (600 ± 60 kg m−3)), respectively. The resulting mass balance biased between −0.20 and 0.09 m water equivalent (w.e.) a-1 using an average (850 ± 60 kg m−3) density assumption throughout the Indus Basin. In the western Himalaya and Hindu Kush, the glacier mass losses are less affected by the average density assumption compared to the Karakoram. The western (Hunza) and central (Shigar) Karakoram glaciers show negligible mass losses of −0.02 ± 0.12 and −0.01 ± 0.13 m w.e. a−1 in contrast to the relatively more negative mass balance (−0.26 ± 0.21 m w.e. a−1) in the eastern (Shyok) Karakoram. All the sub-basins exhibit negative mass balances, with the most negative values ranging from −0.34 ± 0.31 to 0.44 ± 0.27 m w.e. a−1 in the Ravi, Chenab and Jhelum sub-basins of the Himalaya. The whole Indus Basin contributes approximately +0.014 ± 0.016 mm a−1 to the global mean sea-level equivalent.

Changes in glacier volume on Mt. Gongga, southeastern Tibetan Plateau, based on the analysis of multi-temporal DEMs from 1966 to 2015

ABSTRACTThe accelerated retreat of glaciers and the reduction of glacier ice reserves caused by climate change can significantly affect regional water resources and hydrological cycles. Changes in glacier thickness are among the key indicators that reflect this process. We analyzed changes observed in the elevation of glacier surfaces in the Gongga Mountains (GGM) using multi-temporal Digital Elevation Models (DEMs) derived from topographic maps, SRTM, ICESat and ZY-3 data. The results showed that the mean rate of change in glacier surface altitude in the GGM was ~−26.7 ± 2.03 m (0.54 ± 0.04 m a−1) between 1966 and 2015. The mean melt rates differed over different time periods, latterly showing an accelerating trend. As a general rule, glaciers appear to be losing more volume at lower than at higher elevations. Further analysis of these results suggests that supraglacial debris coverage in the GGM promotes higher rates of mass loss.

An efficient representation of glacier dynamics in a semi-distributed hydrological model to bridge glacier and river catchment scales

Glacierised river catchments have been shown to be highly sensitive to climate change, while large populations depend on the water resources originating from them. Hydrological models are used to aid water resource management, yet their treatment of glacier processes is either rudimentary in large-scale applications or linked to fully distributed glacier models that prevent larger model domains. Also, data scarcity in mountainous catchments has hampered the implementation of physically based approaches over entire river catchments. A fully integrated glacier dynamics module was developed for the hydrological model SWIM (SWIM-G) that takes full account of the spatial heterogeneity of mountainous catchments but keeps in line with the semi-distributed disaggregation of the hydrological model. The glacierised part of the catchment is disaggregated into glaciological response units that are based on subbasin, elevation zone and aspect classes. They seamlessly integrate into the hydrological response units of the hydrological model. Robust and simple approaches to ice flow, avalanching, snow accumulation and metamorphism as well as glacier ablation under consideration of aspect, debris cover and sublimation are implemented in the model, balancing process complexity and data availability. The fully integrated model is also capable of simulating a range of other hydrological processes that are common for larger mountainous catchments such as reservoirs, irrigation agriculture and runoff from a diverse soil and vegetation cover. SWIM-G is initialised and calibrated to initial glacier hypsometry, glacier mass balance and river discharge. While the model is intended to be used in medium to large river basins with data-scarce and glacierised headwaters, it is here validated in the data-scarce catchment of the Upper Aksu River, Kyrgyzstan/NW China and in the relatively data-abundant catchment of the Upper Rhone River, Switzerland.

Wildfire Alters Spatial Patterns of Available Soil Nitrogen and Understory Environments in a Valley Boreal Larch Forest

Wildfire, a primary natural disturbance in many forests, affects soil nutrient availability and spatial distributions of forest plants. However, post-fire changes in soil nutrients and spatial patterns of understory environments at fine scales are poorly understood. Here, we characterized spatial patterns of soil nitrogen availability and site characteristics at a 3-year-post-fire and an unburned site in a valley boreal larch forest. We also examined the relationship between soil nitrogen availability and site characteristics. The results showed that the burned site had higher NO3− and lower NH4+ than the control. The herb, litter and coarse wood debris cover was greater at the burned site than at the control site with higher soil pH, depth of the organic horizon (DOH) and shrub cover. Relative variability (coefficient of variation) in soil nitrogen and site characteristic variables at the control site was greater than at the burned site except for shrub and regeneration tree seedling cover. Spatial structure (quantified by semi-variograms) was lacking for soil nitrogen and site characteristic variables except for DOH, herb and shrub cover at the control site, but wildfire created a strong spatial structure for all variables. Shorter spatial autocorrelation ranges of soil nitrogen (1.6–3.5 m) and site characteristic variables (2.6–6.0 m) were detected at the burned site, indicating higher heterogeneity. The spatial scale of soil NH4+ was congruent with those of herb, shrub and regeneration tree seedling cover, indicating local coupling, while that of soil NO3− was not. The number of correlations between soil nitrogen and site characteristic variables in the burned site was greater than in the control. These results indicate that fire could not only create higher heterogeneity patches of soil resources, but also strengthen the local coupling between soil resources and understory vegetation, which may impact the establishment and growth of new individual plants.

Fire legacies in eastern ponderosa pine forests.

Disturbance legacies structure communities and ecological memory, but due to increasing changes in disturbance regimes, it is becoming more difficult to characterize disturbance legacies or determine how long they persist. We sought to quantify the characteristics and persistence of material legacies (e.g., biotic residuals of disturbance) that arise from variation in fire severity in an eastern ponderosa pine forest in North America. We compared forest stand structure and understory woody plant and bird community composition and species richness across unburned, low‐, moderate‐, and high‐severity burn patches in a 27‐year‐old mixed‐severity wildfire that had received minimal post‐fire management. We identified distinct tree densities (high: 14.3 ± 7.4 trees per ha, moderate: 22.3 ± 12.6, low: 135.3 ± 57.1, unburned: 907.9 ± 246.2) and coarse woody debris cover (high: 8.5 ± 1.6% cover per 30 m transect, moderate: 4.3 ± 0.7, low: 2.3 ± 0.6, unburned: 1.0 ± 0.4) among burn severities. Understory woody plant communities differed between high‐severity patches, moderate‐ and low‐severity patches, and unburned patches (all p < 0.05). Bird communities differed between high‐ and moderate‐severity patches, low‐severity patches, and unburned patches (all p < 0.05). Bird species richness varied across burn severities: low‐severity patches had the highest (5.29 ± 1.44) and high‐severity patches had the lowest (2.87 ± 0.72). Understory woody plant richness was highest in unburned (5.93 ± 1.10) and high‐severity (5.07 ± 1.17) patches, and it was lower in moderate‐ (3.43 ± 1.17) and low‐severity (3.43 ± 1.06) patches. We show material fire legacies persisted decades after the mixed‐severity wildfire in eastern ponderosa forest, fostering distinct structures, communities, and species in burned versus unburned patches and across fire severities. At a patch scale, eastern and western ponderosa system responses to mixed‐severity fires were consistent.

Do debris-covered glaciers demonstrate distinctive hydrological behaviour compared to clean glaciers?

Supraglacial debris is known to strongly influence the distribution of glacier surface melt. Since melt inputs drive the formation and evolution of glacial drainage systems, it should follow that the drainage systems of debris-covered glaciers will differ from those of debris-free glaciers. This would have implications for the proglacial runoff regime, subglacial erosion and glacier dynamics. This paper presents analysis of return curves from 33 successful dye injections into the extensively debris-covered Miage Glacier, Italian Alps. It demonstrates that the spatial distribution of supraglacial debris influences the structure and seasonal evolution of the glacial drainage system. Where the debris cover is continuous, melt is lower and the surface topography is chaotic, with many small supraglacial catchments. These factors result in an inefficient englacial/subglacial drainage network beneath continuous debris, which drains to the conduit system emanating from the upper ablation zone. Melt rates are high in areas of clean and dirty ice above the continuous debris. Runoff from these areas is concentrated by inter-moraine troughs into large supraglacial streams, which encourages the early-season development of an efficient englacial/subglacial conduit system downstream of this area. Drainage efficiency from the debris-covered area increases over the melt season but dye-trace transit velocity remains lower than from moulins on the upper glacier. Future runoff models should account for the influence of supraglacial debris on the hydrological system.

Mapping the Late-Holocene Glacial Geomorphology and Glacier Surface Types in the Mt. Harajoriha, Central Tian Shan

Alpine glaciers in the central Tian Shan are an important indicator of climate change and also the freshwater tower for the transboundary countries in Central Asia. Knowledge about the glacier dynamics in the late Holocene, such as the Little Ice Age, and surface zones is still limited. In this study, two headwater basins, the Xiata and the Muzart basins, in the Harajoriha Mountain Range in northwestern China were selected to investigate the glacial landforms and glacier surface types using a combination of geomorphological mapping and remote sensing analysis. Several types of glacial landforms including glacial valleys, moraine complexes, moraine ridges, and trimlines were identified and manually digitized based on the 30 m Shuttle Radar Topography Mission (SRTM) digital elevation model, 10 m Sentinel-2 satellite imagery, and high-resolution images from Google Earth. In addition, an effective automated mapping algorithm was applied to the modern glaciers from a Landsat 8 scene using its optical and thermal bands to classify glacier facies, i.e. ice, snow, and slush zone, and supraglacial debris cover. Field trip to the forelands of the Aerqialeteer Glacier in the Xiata basin allowed detailed mapping of the proglacial environment and provided field checks for the mapping with GPS tracks and panoramic photos. Through this fused approach, the mapping results showed a combination of different sets of information connecting the glacier retreat since the late Holocene to contemporary glacier properties. They will be of particular value for future chronological reconstruction of past glacial events and for understanding how surface characteristics play a role in the heterogeneity of glacier responses to climate change.

Interannual variability in glacier contribution to runoff from a high‐elevation Andean catchment: understanding the role of debris cover in glacier hydrology

AbstractWe present a field‐data rich modelling analysis to reconstruct the climatic forcing, glacier response, and runoff generation from a high‐elevation catchment in central Chile over the period 2000–2015 to provide insights into the differing contributions of debris‐covered and debris‐free glaciers under current and future changing climatic conditions. Model simulations with the physically based glacio‐hydrological model TOPKAPI‐ETH reveal a period of neutral or slightly positive mass balance between 2000 and 2010, followed by a transition to increasingly large annual mass losses, associated with a recent mega drought. Mass losses commence earlier, and are more severe, for a heavily debris‐covered glacier, most likely due to its strong dependence on snow avalanche accumulation, which has declined in recent years. Catchment runoff shows a marked decreasing trend over the study period, but with high interannual variability directly linked to winter snow accumulation, and high contribution from ice melt in dry periods and drought conditions. The study demonstrates the importance of incorporating local‐scale processes such as snow avalanche accumulation and spatially variable debris thickness, in understanding the responses of different glacier types to climate change. We highlight the increased dependency of runoff from high Andean catchments on the diminishing resource of glacier ice during dry years.

Supraglacial debris thickness variability: impact on ablation and relation to terrain properties

Abstract. Shallow ground-penetrating radar (GPR) surveys are used to characterize the small-scale spatial variability of supraglacial debris thickness on a Himalayan glacier. Debris thickness varies widely over short spatial scales. Comparison across sites and glaciers suggests that the skewness and kurtosis of the debris thickness frequency distribution decrease with increasing mean debris thickness, and we hypothesize that this is related to the degree of gravitational reworking the debris cover has undergone and is therefore a proxy for the maturity of surface debris covers. In the cases tested here, using a single mean debris thickness value instead of accounting for the observed small-scale debris thickness variability underestimates modelled midsummer sub-debris ablation rates by 11 %–30 %. While no simple relationship is found between measured debris thickness and morphometric terrain parameters, analysis of the GPR data in conjunction with high-resolution terrain models provides some insight into the processes of debris gravitational reworking. Periodic sliding failure of the debris, rather than progressive mass diffusion, appears to be the main process redistributing supraglacial debris. The incidence of sliding is controlled by slope, aspect, upstream catchment area and debris thickness via their impacts on predisposition to slope failure and meltwater availability at the debris–ice interface. Slope stability modelling suggests that the percentage of the debris-covered glacier surface area subject to debris instability can be considerable at glacier scale, indicating that up to 32 % of the debris-covered area is susceptible to developing ablation hotspots associated with patches of thinner debris.

Spatio‐temporal variability in contributions to low flows in the high Alpine Poschiavino catchment

AbstractIn order to predict streamflow accurately during extended dry periods, we need to understand the spatial variability of low flows and the extent to which it is affected by the spatial organization and drainage of catchment subsurface storage areas. This is especially true in Alpine catchments with widely varying topography, lithology, sediment deposits, and soil properties. Field measurements in the Poschiavino catchment in southern Switzerland, during a winter recession period without recharge, provided a unique opportunity to demonstrate the connections between subsurface storage areas, low flows, and their variability. We measured discharge in four nested sub‐catchments during seven field campaigns in the winter of 2013–2014. We analysed stream water electrical conductivity (EC) and water chemistry to identify the areas contributing to low‐flow discharge and estimated their contributions. Sediment cover type and thickness were mapped using a recently developed tool for geomorphology‐based storage classification of mountainous terrain, to determine the physical properties of the subsurface storage areas contributing to low‐flow discharge. Recession analyses combined with water chemistry data allowed the detection of different drainage timescales and the estimation of storage potential of the unconsolidated (Quaternary) deposits. We found substantial spatial variation in storage depletion between the sub‐catchments, ranging from 54 mm to 200 mm for the four‐month monitoring period. Variability in low‐flow contributions from different catchments and different recession behavior could be related to the differences in the estimated storage potential. For some point sources, we could quantify the contributing area and thus quantify low flows at the hillslope scale. Overall, the low‐flow variability is mostly related to the fraction of precipitation that recharges subsurface storage areas and to the properties influencing their drainage. To capture these processes, we suggest low‐flow geomorphological mapping approaches that consider not only morphometric (shape of the landscape) and geologic (properties of the bedrock) controls but also the water storage potential of debris cover and weathered rock.

Global Assessment of Supraglacial Debris‐Cover Extents

AbstractRocky debris on glacier surfaces influences ice melt rates and the response of glaciers to climate change. However, scarce data on the extent and evolution of supraglacial debris cover have so far limited its inclusion in regional to global glacier models. Here we present global data sets of supraglacial debris‐cover extents, based on Landsat 8 and Sentinel‐2 optical satellite imagery. We find that about 4.4% (~26,000 km2) of all glacier areas (excluding the Greenland ice sheet and Antarctica) are covered with debris, but that the distribution is heterogeneous. The largest debris‐covered areas are located in high‐mountain ranges, away from the poles. At a global scale, we find a negative scaling relationship between glacier size and percentage of debris. Therefore, the influence of debris cover on glacier mass balances is expected to increase in the future, as glaciers continue to shrink.

Post‐little ice age paraglacial processes and landforms in the high Iberian mountains: A review

AbstractThree Iberian mountain ranges encompassed glaciers during the Little Ice Age (LIA): the Pyrenees, Cantabrian Mountains, and Sierra Nevada. The gradual warming trend initiated during the second half of the 19th century promoted the progressive shrinking of these glaciers, which completely melted during the first half of the 20th century in the Cantabrian mountains and Sierra Nevada and reduced by 80% of their LIA extent in the Pyrenees. In these formerly glaciated environments, the transition between glacial and periglacial conditions results in an accelerated paraglacial readjustment, with very active geomorphic processes. Cirque walls generate a large amount of sediments through rock‐falls and slides. LIA moraines, devoid of vegetation and composed of highly unstable sediments, are being intensely mobilized by slope processes. Inside the moraines, the shrinking of LIA glaciers favoured the development of buried ice patches, with permafrost‐related landforms, small periglacial features generated by solifluction, and cryoturbation processes and remarkable hydrological changes. Present‐day morphodynamics is mostly related to seasonal frost, though patches of permafrost have formed in contact with the buried ice, undergoing a process of degradation because it is not balanced with present‐day climate. This is reflected in the occurrence of multiple collapses and subsidence of the debris cover where the frozen bodies sit. Next to the small glaciated environments in the highest Pyrenean massifs, there is a permafrost belt undergoing also rapid geomorphic changes. Based on the observed processes, we discuss spatio‐temporal patterns of paraglacial readjustment in Iberian mountains and compare it with other midlatitude mountain environments.