Glacier Evolution Research Articles

Be‐10 Dating of Ice‐Marginal Moraines in the Khumbu Valley, Nepal, Central Himalaya, Reveals the Response of Monsoon‐Influenced Glaciers to Holocene Climate Change

AbstractThe dynamic response of large mountain glaciers to climatic forcing operates over timescales of several centuries and therefore understanding how these glaciers change requires observations of their behavior through the Holocene. We used Be‐10 exposure‐age dating and geomorphological mapping to constrain the evolution of glaciers in the Khumbu Valley in the Everest region of Nepal. Khumbu and Lobuche Glaciers are surrounded by high‐relief lateral and terminal moraines from which seven glacial stages were identified and dated to 7.4 ± 0.2, 5.0 ± 0.3, 3.9 ± 0.1, 2.8 ± 0.2, 1.3 ± 0.1, 0.9 ± 0.02, and 0.6 ± 0.16 ka. These stages correlate to each of the seven latest Holocene regional glacial stages identified across the monsoon‐influenced Himalaya, demonstrating that a coherent record of high elevation terrestrial palaeoclimate change can be extracted from dynamic mountain landscapes. The time‐constrained moraine complex represents a catchment‐wide denudation rate of 0.8–1.4 mm a−1 over the last 8 kyr. The geometry of the ablation area of Khumbu Glacier changed around 4 ka from a broad, shallow ice tongue to become narrower and thicker as restricted by the topographic barrier of the terminal moraine complex.

Volume, evolution, and sedimentation of future glacier lakes in Switzerland over the 21st century

Abstract. Ongoing climate change and associated glacier retreat is causing rapid environmental change, including shifts in high-alpine landscapes. Glacier lakes, which can form in topographical depressions left behind by glacier retreat, are prominent features within such landscapes. Whilst model-based estimates for the number and area of future glacier lakes exist for various mountain regions across the world, the exact morphology and temporal evolution remain largely unassessed. Here, we leverage a recently released, measurement-based estimate for the subglacial topography of all glaciers in the Swiss Alps to provide an estimate about the number, size, time of emergence, and sediment infill of future glacier lakes. The topographical information is based on 2450 km of measured ice thickness profiles, whilst the temporal evolution of glaciers is obtained from a glacier evolution model forced with an ensemble of climate projections. We estimate that up to 683 potential lakes with an area >5000 m2 and a depth >5 m could emerge across the Swiss Alps if glaciers were to disappear completely, with the potential to hold a total water volume of up to 1.16 [1.05, 1.32] km3 (numbers and 95 % confidence interval). For a middle-of-the-road climate scenario, we estimate that about 10 % (0.12 [0.04, 0.18] km3) and 48 % (0.56 [0.26, 0.67] km3) of this volume could be realized by 2050 and 2100, respectively. In a first-order assessment, we also estimate that ca. 45 % of the newly emerging glacier lakes (260 out of 570) will be transient features, i.e. will disappear again before the end of the century owing to refilling with sediments released by glacial erosion and proglacial sediment transport.

Multi-centennial mass balance of perennial ice deposits in Alpine caves mirrors the evolution of glaciers during the Late Holocene

Mid-latitude alpine caves preserve a record of past solid precipitation during winter, locally spanning several centuries to millennia. Dating organic macro-remains trapped in ice layers allows the determination of timing and duration of past periods of positive and negative ice mass balance. We present here the largest comparative study of ice cave sites yet published, using Bayesian age-modelling on a database comprising 107 radiocarbon dates, spread over eight caves in the Austrian Alps. We show that periods of positive mass balance coincide with past glacier advances. We find organic and macro-remain rich layers dated to the Medieval Climate Anomaly (between 850 and 1200 CE) marking widespread ice retreat. We demonstrate positive ice mass balance at all studied sites for the Little Ice Age, coinciding with the largest glacier advances in the Holocene between 1400 and 1850 CE. At the sites with records spanning over 2000 years, positive mass balance is also observed during the periods from 300 BCE to 100 CE and 600–800 CE. These subterranean ice deposits show widespread evidence of accelerated negative mass balances in recent years and their record is under imminent threat of disappearing.

Validation of ERA5-Land temperature and relative humidity on four Peruvian glaciers using on-glacier observations

Weather and climate conditions drive the evolution of tropical glaciers which play an important role as water reservoirs for Peruvian inhabitants in the arid coast and semi-arid Andean region. The scarcity of long-term high-quality observations over Peruvian glaciers has motivated the extensive use of reanalysis data to describe the climatic evolution of these glaciers. However, the representativeness and uncertainties of these reanalysis products over these glaciers are still poorly constrained. This study evaluates the ability of the ERA5-Land reanalysis (ERA5L) to reproduce hourly and monthly 2 m air temperature and relative humidity (T2m and Rh2m, respectively) over several Peruvian glaciers. We compared the ERA5L with data from four on-glacier automatic weather stations (AWS), whose hourly time series were completed with nearby stations, for the period January 2017 to December 2019. Results indicates a better performance of the reanalysis for T2m (r >0.80) than for Rh2m (∼0.4< r <∼0.6) in all four glaciers. Concerning the observations, both parameters show a daily cycle influenced by the presence of the glacier. This influence is more prominent during the dry months when the so-called glacier damping and cooling effects are stronger. On a monthly time scale, the ERA5L validation for both parameters are better in wet outer tropical sites (RMSE between ±0.2°C for T2m and between 3%–7% for Rh2m) rather than in dry outer tropical sites (RMSE between ±0.2°C for T2m and between 3%–7% for Rh2m). Among all sites considered in the study, the Rh2m bias is the highest in the Cavalca glacier (correlation of 0.81; RMSE 13%, MAE 11% and bias 8.3%) and the lowest in Artesonraju glacier (correlation of 0.96; RMSE 3%; MAE 2.3% and bias — 0.8%). Based on certain considerations outlined in this paper, it is appropriate to use ERA5L to characterize T2m and Rh2m conditions on Peruvian glaciers, particularly in the wet outer tropics.

Snow Multidata Mapping and Modeling (S3M) 5.1: a distributed cryospheric model with dry and wet snow, data assimilation, glacier mass balance, and debris-driven melt

Abstract. By shifting winter precipitation into summer freshet, the cryosphere supports life across the world. The sensitivity of this mechanism to climate and the role played by the cryosphere in the Earth's energy budget have motivated the development of a broad spectrum of predictive models. Such models represent seasonal snow and glaciers with various complexities and generally are not integrated with hydrologic models describing the fate of meltwater through the hydrologic budget. We present Snow Multidata Mapping and Modeling (S3M) v5.1, a spatially explicit and hydrology-oriented cryospheric model that simulates seasonal snow and glacier evolution through time and that can be natively coupled with distributed hydrologic models. Model physics include precipitation-phase partitioning, snow and glacier mass balances, snow rheology and hydraulics, a hybrid temperature-index and radiation-driven melt parametrization, and a data-assimilation protocol. Comparatively novel aspects of S3M are an explicit representation of the spatial patterns of snow liquid-water content, the implementation of the Δh parametrization for distributed ice-thickness change, and the inclusion of a distributed debris-driven melt factor. Focusing on its operational implementation in the northwestern Italian Alps, we show that S3M provides robust predictions of the snow and glacier mass balances at multiple scales, thus delivering the necessary information to support real-world hydrologic operations. S3M is well suited for both operational flood forecasting and basic research, including future scenarios of the fate of the cryosphere and water supply in a warming climate. The model is open source, and the paper comprises a user manual as well as resources to prepare input data and set up computational environments and libraries.

Existence of Glacier Anomaly in the Interior and Northern Tibetan Plateau between 2000 and 2012

There was sufficient evidence to indicate a nearly balanced glacier mass change (termed glacier anomaly) for Karakoram Mts. since the 1970s, in contrast to worldwide glacier mass losses caused by climate warming. Recently, this anomalous phenomenon was detected over the neighboring western Kunlun and Pamir Mts. However, the southeastern limit of this glacier anomaly remains uncertain, owing to the paucity of glacier mass balance observations across the interior and northern Tibetan Plateau (INTP). In this study, we presented a decadal glacier mass balance estimation in the INTP by differencing the SRTM DEM with the topographic data produced from TanDEM-X bistatic InSAR images. From 2000 to 2012, decade-average glacier mass balances of between −0.339 ± 0.040 and 0.237 ± 0.078 m w.e. yr−1 were detected over 22 glacierized areas. Significantly, we found a gradient and switch of glacier mass loss over the southeastern portion to glacier mass gain over the northwestern portion. This varying spatial pattern illustrates that glacier anomaly has existed over the northwestern or even central zone of the INTP since the early 21st century. This study provides important evidence for the model simulation of both glacier evolution and atmospheric circulations in investigating the prevailing mechanism of the regional anomalous phenomenon.

Recent evolution of glaciers in the Cocuy‐Güican Mountains (Colombian Andes) and the hydrological implications

AbstractThis study extends knowledge of the evolution of glacier shrinkage in the Cocuy‐Güican Mountains since the maximum glacier extent of the Little Ice Age (LIA). Mass balance data for the Ritacuba Glacier since 2009 were acquired and compared with available data for the Conejeras Glacier (Los Nevados National Park). This study also investigated the hydrological significance of Colombian glaciers, which is still largely unknown because the available information is very limited. Glaciers in the Cocuy‐Güican Mountains covered 13.2 km2 in 2019 compared with 127.8 km2 during the maximum LIA, representing a shrinkage of 89.7%. Analysis of glacier cover observations made in 1955, 1994, 2010, and 2019 revealed that the rate of ice loss was greatest from 1994 to 2010 (0.59 km2 yr−1), and was then almost halved from 2010 to 2019 (0.34 km2 yr−1). This slowing of glacier retreat is consistent with the moderate negative mass balance measured for 2009–2019, and an accumulated loss of 1766 mm w.e. (mm water equivalent). The progressive confinement of glaciers to higher elevations, ice accumulation in topographic locations providing shelter from solar radiation, and an absence of recent marked climatic anomalies could explain why the Cocuy‐Güican glaciers have temporally reached a near equilibrium state conditions. This contrasts with the Conejeras Glacier, where 47,000 mm w.e. has been lost in the same period. The available data on runoff and isotope tracers of streamflow and precipitation suggest that precipitation rather than glacier melt water exerts primary control over the hydrological variability at high elevation sites.

The Remotely and Directly Obtained Results of Glaciological Studies on King George Island: A Review

Climate warming has become indisputable, and it is now crucial to increase our understanding of both the mechanisms and consequences of climate change. The Antarctic region is subjected to substantial changes, the trends of which have been recognized for several decades. In the South Shetland Islands, the most visible effect of climate change is progressive deglaciation. The following review focuses on past glaciological studies conducted on King George Island (KGI). The results of collected cryosphere element observations are discussed herein in a comprehensive manner. Our analysis showed that there is a lack of temporal as well as spatial continuity for studies on the basic mass balance parameters on the entire KGI ice dome and only Bellingshausen Dome has a relatively long history of data collection. The methodologies of past work, which have improved over time, are also discussed. When studying the glacier front fluctuations, the authors most frequently use a 1956 aerial photography as reference ice coverage. This was the case for seven papers, while other sources are seldomly mentioned. In other papers as many as 41 other sources were used, and therefore comparison to photos taken up to 60 years later can give misleading trends, as small glaciers may have both advanced and retreated in that time. In the case of glacial velocities there is also an apparent lack of consistency, as different glaciers were indicated as the fastest on KGI. Only Lange, Anna, Crystal, Eldred, and eastern part of Usher glaciers were determined by more than one author as the fastest. Additionally, there are gaps in the KGI Ground Penetrating Radar (GPR) survey area, which includes three ice domes: the Warszawa Icefield, the Krakow Icefield, and eastern part of King George Island. Ideas for further work on the topic are also suggested, allowing for easier access to data and thus contributing to a better understanding of glacier development mechanisms.

Modelling the mass budget and future evolution of Tunabreen, central Spitsbergen

Abstract. The 26 km long tidewater glacier Tunabreen is the most frequently surging glacier in Svalbard, with four documented surges in the past 100 years. We model the evolution of this glacier with a minimal glacier model (MGM), in which ice mechanics, calving, and surging are parameterized. The model geometry consists of a flow band to which three tributaries supply mass. The calving rate is set to the mean observed value for the period 2012–2019 and kept constant. For the past 120 years, a smooth equilibrium line altitude (ELA) history is reconstructed by finding the best possible match between observed and simulated glacier length. There is a modest correlation between this reconstructed ELA history and an ELA history based on meteorological observations from Longyearbyen. Runs with and without surging show that the effect of surging on the long-term glacier evolution is limited. Due to the low surface slope and associated strong height–mass-balance feedback, Tunabreen is very sensitive to changes in the ELA. For a constant future ELA equal to the reconstructed value for 2020, the glacier front will retreat by 8 km during the coming 100 years. For an increase in the ELA of 2 m a−1, the retreat is projected to be 13 km, and Tunabreen becomes a land-terminating glacier around 2100. The calving parameter is an important quantity: increasing its value by 50 % has about the same effect as a 35 m increase in the ELA, with the corresponding equilibrium glacier length being 17.5 km (as compared to 25.8 km in the reference state). Response times vary from 150 to 400 years, depending on the forcing and on the state of the glacier (tidewater or land-terminating).

EXTRACTION OF ICE FLOW VELOCITY BY COMBINATION OF DINSAR AND OFFSET TRACKING METHODS FOR PINE GLACIER, WEST ANTARCTICA

Abstract. West Antarctica is the main contributor to global sea level rise at present and in the coming decades, since it occupies 80% of mass loss in Antarctica. In this paper, a combination of Differential Interferometric Synthetic Aperture Radar (DInSAR) and offset tracking technology is used to extract the ice flow velocity of the Pine Island Glacier (PIG), a typical glacier in West Antarctica. Due to the large deformation gradients in PIG, DInSAR technology is used to extract ice flow velocity in bare rock and mountains, then offset tracking technology is used to extract ice flow velocity in areas where glaciers collapse frequently. Finally, the above two results are mosaiced into a new image of the interannual ice flow velocity of PIG in 2017. Through qualitative and quantitative evaluation, it is found that the ice flow velocity extracted by the combination has high accuracy in both high and low velocity areas. In summary, we concluded that the combination of DInSAR and Offset tracking can obtain reliable ice velocity products in glaciers that change rapidly. This combination is of scientific significance for monitoring the movement and evolution of glaciers in the West Antarctica.

Changes in equilibrium-line altitude and implications for glacier evolution in the Asian high mountains in the 21st century

Changes in equilibrium-line altitude and implications for glacier evolution in the Asian high mountains in the 21st century

Ice‐Dynamical Glacier Evolution Modeling—A Review

AbstractGlaciers play a crucial role in the Earth System: they are important water suppliers to lower‐lying areas during hot and dry periods, and they are major contributors to the observed present‐day sea‐level rise. Glaciers can also act as a source of natural hazards and have a major touristic value. Given their societal importance, there is large scientific interest in better understanding and accurately simulating the temporal evolution of glaciers, both in the past and in the future. Here, we give an overview of the state of the art of simulating the evolution of individual glaciers over decadal to centennial time scales with ice‐dynamical models. We hereby highlight recent advances in the field and emphasize how these go hand‐in‐hand with an increasing availability of on‐site and remotely sensed observations. We also focus on the gap between simplified studies that use parameterizations, typically used for regional and global projections, and detailed assessments for individual glaciers, and explain how recent advances now allow including ice dynamics when modeling glaciers at larger spatial scales. Finally, we provide concrete recommendations concerning the steps and factors to be considered when modeling the evolution of glaciers. We suggest paying particular attention to the model initialization, analyzing how related uncertainties in model input influence the modeled glacier evolution and strongly recommend evaluating the simulated glacier evolution against independent data.

Modelling supraglacial debris-cover evolution from the single-glacier to the regional scale: an application to High Mountain Asia

Abstract. Currently, about 12 %–13 % of High Mountain Asia’s glacier area is debris-covered, which alters its surface mass balance. However, in regional-scale modelling approaches, debris-covered glaciers are typically treated as clean-ice glaciers, leading to a bias when modelling their future evolution. Here, we present a new approach for modelling debris area and thickness evolution, applicable from single glaciers to the global scale. We derive a parameterization and implement it as a module into the Global Glacier Evolution Model (GloGEMflow), a combined mass-balance ice-flow model. The module is initialized with both glacier-specific observations of the debris' spatial distribution and estimates of debris thickness. These data sets account for the fact that debris can either enhance or reduce surface melt depending on thickness. Our model approach also enables representing the spatiotemporal evolution of debris extent and thickness. We calibrate and evaluate the module on a selected subset of glaciers and apply GloGEMflow using different climate scenarios to project the future evolution of all glaciers in High Mountain Asia until 2100. Explicitly accounting for debris cover has only a minor effect on the projected mass loss, which is in line with previous projections. Despite this small effect, we argue that the improved process representation is of added value when aiming at capturing intra-glacier scales, i.e. spatial mass-balance distribution. Depending on the climate scenario, the mean debris-cover fraction is expected to increase, while mean debris thickness is projected to show only minor changes, although large local thickening is expected. To isolate the influence of explicitly accounting for supraglacial debris cover, we re-compute glacier evolution without the debris-cover module. We show that glacier geometry, area, volume, and flow velocity evolve differently, especially at the level of individual glaciers. This highlights the importance of accounting for debris cover and its spatiotemporal evolution when projecting future glacier changes.

The Relative Impacts of Initialization and Climate Forcing in Coupled Ice Sheet‐Ocean Modeling: Application to Pope, Smith, and Kohler Glaciers

AbstractCoupled ice sheet‐ocean models are beginning to be used to study the response of ice sheets to ocean warming. Initializing an ice‐ocean model is challenging and can introduce nonphysical transients, and the extent to which such transients can affect model projections is unclear. We use a synchronously‐coupled ice‐ocean model to investigate evolution of Pope, Smith and Kohler Glaciers, West Antarctica, over the next half‐century. Two methods of initialization are used: In one, the ice‐sheet model is constrained with observed velocities in its initial state; in another, the model is constrained with both velocities and grounded thinning rates over a 4‐year period. Each method is applied to two basal sliding laws. For each resulting initialization, two climate scenarios are considered: one where ocean conditions during the initialization period persist indefinitely, and one where the ocean is in a permanent “warm” state. At first, model runs initialized with thinning data exhibit volume loss rates much closer to observed values than those initialized with velocity only, but after 1–2 decades, the forcing primarily determines rates of volume loss and grounding line retreat. Such behavior is seen for both basal sliding laws, although volume loss rates differ quantitatively. Under the “warm” scenario, a grounding line retreat of ∼30 km is simulated for Smith and Kohler, although variation in total retreat due to initialization is nearly as large as that due to forcing. Furthermore it is questionable whether retreat will continue due to narrowing of submarine troughs and limiting of heat transport by bathymetric obstacles.

LGM glaciers in the SE Mediterranean? First evidence from glacial landforms and 36Cl dating on Mount Lebanon

The hydrological conditions during the Last Glacial Maximum in the Eastern Mediterranean have long been debated. In particular, very little is known about the development of high-altitude glaciers in this region. In the present study, morphological and sedimentological evidence, such as U-shape valleys, glacial cirques, and morainic ridges, are used to identify past glaciers on the western slopes of Mount Lebanon, which is by far the highest topography along the eastern Mediterranean coast (3088 m a.s.l.). The geometry of three paleo-glaciers is reconstructed: one 8-km-long glacial tongue stretching from the highest peak down to an elevation of about 2000 m a.s.l., and two small glacial cirques located at about 2300 m a.s.l. The age of maximum glacier advance is constrained from a 10-m-deep 36Cl depth-profile within the diamicton of one terminal moraine, and yields an early LGM timing (between 25.4 ± 3.1 ka and 31.2 ± 4.5 ka), which agrees with the chronology of the last maximum glacier extension in the Eastern Mediterranean. The three glaciers together form a consistent group, compatible with an Equilibrium Line Altitude (ELA) at 2400 m estimated using the Accumulation Area Ratio method. A degree-day model is further used to infer the temperature and precipitation conditions required to depress the LGM ELA to 2400 m. Our results suggest that the observed glacier extent is compatible with early LGM paleo-temperature records of the area only under the condition of decreased precipitation relative to modern times. These findings thus support an early LGM climate of the Levant region that was drier than today.

Modern sand supply of the Tengger Desert and temporal variations in sand provenance driven by northern Hemisphere glaciation

Knowledge of surface processes in source-to-sink system of sand is important for the development and utilization of natural resources in desert. However, the surface processes within Chinese deserts, especially their response to long-term global cooling since the Late Pliocene, are still poorly understood. Here, we use geochemical analyses of aeolian sand (63–500 μm) samples from a drill core in the central Tengger Desert of China, together with widely-distributed surface samples, to elucidate the present and past surface processes. The main sand provenances for the southwestern and northeastern parts of the Tengger Desert are the northeastern Tibetan Plateau and the Alxa Block, respectively. After the inception of Northern Hemisphere glaciation at ∼2.7 Ma, the expansion of ice sheets in high northern latitudes was the major trigger of an increased sand supply from the Alxa Block to the desert hinterland, via the strengthening of the Siberian High. After the Mid-Pleistocene Transition (MPT), during ∼1.2–0.8 Ma, the development of mountain glaciers in the Tibetan Plateau was the major trigger of an enhanced fluvial sand supply from the northeastern Tibetan Plateau to the hinterland. Further, we propose that these two triggers may subsequently have caused substantial variations in the sand supply to most of the deserts surrounding the northern Tibetan Plateau rather than only the Tengger, due to the similarity of their sand supply mechanisms. And the enhanced fluvial sand supply from the northern Tibetan Plateau may be of significance to the continuous existence of a sand sea landscape in the surrounding basins during the post-MPT interval. Overall, our study suggests that glaciation in the high latitudes and the high mountains of the Northern Hemisphere have a significant impact on the two sand supply agents, i.e., wind and water, of Chinese deserts respectively.

Recent Evolution of Glaciers in the Manaslu Region of Nepal From Satellite Imagery and UAV Data (1970–2019)

Glacierized mountain ranges such as the Himalaya comprise a variety of glacier types, including clean and debris-covered glaciers. Monitoring their behaviour over time requires an assessment of changes in area and elevation along with surface features and geomorphology. In this paper we quantify the surface evolution of glacier systems in the Manaslu region of Nepal over the last five decades using 2013/2019 multi-sensor imagery and elevation data constructed from 1970 declassified Corona imagery and 1970 declassified Corona imagery. We investigate area changes, glacier thickness, geodetic glacier mass balance and surface velocity changes at regional scales and focus on the Ponkar Glacier and Thulagi Glacier and Lake for an in-depth assessment of surface geomorphology and surface feature dynamics (ponds, vegetation and ice cliffs). The time series of surface elevation changes for the lower ablation area of Ponkar Glacier is extended using 2019 UAV-based imagery and field-based ablation rates measured over the period 2016–2019. Glaciers in the Manaslu region experienced a mean area loss of −0.26 ± 0.0001% a−1 between 1970 and 2019. The mean surface lowering was −0.20 ± 0.02 ma−1 over the period 1970 to 2013, corresponding to a regional geodetic mass balance of −0.17 ± 0.03 m w. e.a−1. Overall, debris-covered glaciers had slightly higher thinning rates compared to clean ice glaciers; lake-terminating glaciers had double thinning rates compared to land-terminating glaciers. Individual glacier mass balance was negatively controlled by glacier slope and mean glacier elevation. During the period 1970 to 2013, Ponkar Glacier had a geodetic mass balance of −0.06 ± 0.01 m w. e.a−1, inversely correlated with parts of the central trunk thickening. Between 2013 and 2019 there was a nine-fold increase in the thinning rates over the lower parts of the glacier tongue relative to the period 1970–2013. Ice-surface morphology changes between 1970 and 2019 on Ponkar Glacier include a decrease in ogives and open crevasses, an increase in ice cliffs and ponds and the expansion of the supraglacial debris and ice-surface vegetation. These changes point to reduced ice-dynamic activity and are commensurate with the observed recession and negative glacier mass balance over the last five decades.

Evolution of geodetic mass balance over the largest lake-terminating glacier in the Tibetan Plateau with a revised radar penetration depth based on multi-source high-resolution satellite data

The southeast Tibetan Plateau has become the region with the fastest rate of mass loss within High Mountain Asia (HMA), and the mass loss from lake-terminating glaciers in the region is the most significant. However, the evolution of the mass change of lake-terminating glaciers is still unclear, and accurate penetration depth estimates of C-band and X-band radar in glacierized areas are unknown, which greatly hinders the monitoring and understanding of the glacier evolution. Hence, in this study, taking the largest lake-terminating glacier in HMA (the Yanong Glacier) as the study object, we first evaluated the X-band penetration depth using Pléiades and TanDEM-X data, and further combined the penetration depth difference between the C-band and X-band radar (1.14 ± 0.14 m) to estimate the C-band penetration depth based on a newly proposed area-weighted strategy, according to surface categories. The results indicate a region-wide average penetration depth of 1.84 ± 0.59 m and 2.98 ± 0.61 m for the X-band and C-band radar, respectively, demonstrating either an underestimation of 60% or an overestimation of 202% for the previous C-band penetration corrections. On this basis, we determined the multi-temporal glacier mass balance for the Yanong Glacier, using KH-9, SRTM DEM, SPOT-7, and multi-orbit TanDEM-X data. The results show that the Yanong Glacier has been in a state of serious mass loss (at a rate of −0.73 ± 0.13 m w.e./a) during 1974–2015, and the rate of mass loss has doubled since 2000 (i.e., −0.48 ± 0.20 (1974–2000), −0.95 ± 0.20 (2000−2012), and − 1.02 ± 0.53 m w.e./a (2012–2015)). Moreover, the inter-annual mass change has shown a highly volatile and accelerating trend during 2012–2015 (i.e., −0.47 ± 0.85 (2011−2012), −0.87 ± 0.66 (2012−2013), −1.26 ± 0.72 (2013–2014), and − 1.58 ± 0.44 m w.e./a (2014–2015)). By further analyzing the ERA5 reanalysis data and the changes in the proglacial lake and glacier dynamics, we can conclude that, qualitatively: 1) the long-term acceleration of mass loss at the Yanong Glacier has mainly been controlled by climate change, but has also been affected by the glacier dynamics to a limited extent; and 2) the inter-annual acceleration of mass loss has, to some extent, been affected by subglacial/englacial melting. The findings of this study will not only provide an accurate penetration depth correction for future relevant studies, but will also significantly improve the understanding of the evolution of lake-terminating glaciers in HMA and other high mountain areas of the world.

Evolution and origin of global glaciers and their impacts on the environment

Glacier is one of the most important freshwater resources. However, the situation of glaciers is changed nowadays because of global warming. This study aimed at two types of glaciers, and analyzed their evolution, origin and impacts on the environment. As for typical melting glaciers, the Vatnajokull glacier and Urumqi Glacier No.1, Tianshan Mountains gradually melt over time. Their causes are related to atmospheric River events, black carbon pollutants, ice recycling, global warming, and human activities. The impacts of glacier retreat mainly include sea-level rise and can alleviate global warming and cause surface deformation. However, there are some “weird” glaciers because they are expanded in the past few years rather than shrinking. Those glaciers are highly concentrated in the Karakoram, Pamir Mountains, Svalbard, the Canadian Arctic islands, Alaska, and Iceland, which are called surged-type glaciers. Hydrological control, thermal regime, deformable bed hypothesis, and critical mass are reasons for glacier expansion. There is no clear and unified opinion on the advantages or disadvantages of glacier expansion.

Coincident evolution of glaciers and ice-marginal proglacial lakes across the Southern Alps, New Zealand: Past, present and future

Global glacier mass loss is causing expansion of proglacial landscapes and producing meltwater that can become impounded as lakes within natural topographic depressions or ‘overdeepenings’. It is important to understand the evolution of these proglacial landscapes for water resources, natural hazards and ecosystem services. In this study we (i) overview contemporary loss of glacier ice across the Southern Alps of New Zealand, (ii) analyse ice-marginal lake development since the 1980s, (iii) utilise modelled glacier ice thickness to suggest the position and size of future lakes, and (iv) employ a large-scale glacier evolution model to suggest the timing of future lake formation and future lake expansion rate. In recent decades, hundreds of Southern Alps glaciers have been lost and those remaining have fragmented both by separation of tributaries and by detachment of ablation zones. Glaciers with ice-contact margins in proglacial lakes (n > 0.1 km2 = 20 in 2020) have experienced the greatest terminus retreat and typically twice as negative mass balance compared to similar-sized land-terminating glaciers. Our analysis indicates a positive relationship between mean glacier mass balance and rate of lake growth (R2 = 0.34) and also with length of an ice-contact lake boundary (R2 = 0.44). We project sustained and relatively homogenous glacier volume loss for east-draining basins but in contrast a heterogeneous pattern of volume loss for west-draining basins. Our model results show that ice-marginal lakes will increase in combined size by ~150% towards 2050 and then decrease to 2100 as glaciers disconnect from them. Overall, our findings should inform (i) glacier evolution models into which ice-marginal lake effects need incorporating, (ii) studies of rapid landscape evolution and especially of meltwater and sediment delivery, and (iii) considerations of future meltwater supply and water quality.