Glacier Equilibrium Line Altitude Research Articles

Differential depression of the glacier equilibrium-line altitudes in the Yarlung Zangbo Downstream Basin in the Last Glacial Maximum compared to the pre-industrial era

Differential depression of the glacier equilibrium-line altitudes in the Yarlung Zangbo Downstream Basin in the Last Glacial Maximum compared to the pre-industrial era

Climate Model Simulations of the Effects of Orbital Parameters on Glacier Equilibrium Line Altitude

AbstractThe effects of obliquity and precession on conditions favorable for Northern Hemisphere glaciation are explored using an energy balance and mass balance model of equilibrium line altitude (ELA), the height on a glacier where accumulation and ablation are in balance annually. Climate forcing for the ELA model is obtained from idealized single‐forcing orbital simulations with two atm‐ocean general circulation models, Geophysical Fluid Dynamics Laboratory (GFDL) CM2.1 and National Center for Atmospheric Research (NCAR) Community Earth System Model version 1.2. Over Scandinavia and Baffin Island, the respective locations in which the Scandinavian and Laurentide ice sheets are thought to have originated, low obliquity and perihelion at the boreal winter solstice are associated with lower ELA values, as would be expected from the orbital theory of the ice ages. Linear reconstructions of ELA variations over the past 800 kyr indicate that precession dominated ELA variations in Scandinavia and Baffin Island in the GFDL model, and in Scandinavia in the NCAR model. Obliquity and precession played equal roles in Baffin Island in the NCAR model. A decomposition of the ELA responses finds that the effects of ablation on ELA are much larger than the effects of precipitation. Overall, the findings of this study point to precession being a more important factor in glacial inception than obliquity, which contrasts with previous findings in which obliquity had a slightly larger effect on positive degree days (PDDs), a simple metric for ablation. This is likely due to differences in seasonality of melt from the ELA model and PDDs.

Equilibrium line altitudes of alpine glaciers in Alaska suggest Last Glacial Maximum summer temperature was 2–5 °C lower than during the pre-industrial

Abstract. The lack of continental ice sheets in Alaska during the Last Glacial Maximum (LGM; 26–19 ka) has long been attributed to extensive aridity in the western Arctic. More recently, climate model outputs, a few isolated paleoclimate studies, and global paleoclimate synthesis products show mild summer temperature depressions in Alaska compared to much of the high northern latitudes. This suggests the importance of limited summer temperature depressions in controlling the relatively limited glacier growth during the LGM. To explore this further, we present a new statewide map of LGM alpine glacier equilibrium line altitudes (ELAs), LGM ΔELAs (LGM ELA anomalies relative to the Little Ice Age, LIA), and ΔELA-based estimates of temperature depressions across Alaska to assess paleoclimate conditions. We reconstructed paleoglacier surfaces in ArcGIS to calculate ELAs using an accumulation area ratio (AAR) of 0.58 and an area–altitude balance ratio (AABR) of 1.56. We calculated LGM ELAs (n= 480) in glaciated massifs in the state, excluding areas in southern Alaska that were covered by the Cordilleran Ice Sheet. The data show a trend of increasing ELAs from the southwest to the northeast during both the LGM and the LIA, indicating a consistent southern Bering Sea and northernmost Pacific Ocean precipitation source. Our LGM–LIA ΔELAs from the Alaska Range, supported with limited LGM–LIA ΔELAs from the Brooks Range and the Kigluaik Mountains, average to −355 ± 176 m. This value is much greater than the global LGM average of ca. −1000 m. Using a range of atmospheric lapse rates, LGM–LIA ΔELAs in Alaska translate to summer cooling of < 2–5 ∘C. Our results are consistent with a growing number of local climate proxy reconstructions and global data assimilation syntheses that indicate mild summer temperature across Beringia during the LGM. Limited LGM summer temperature depressions could be explained by the influence of Northern Hemisphere ice sheets on atmospheric circulation.

10Be exposure ages and paleoenvironmental significance of rock glaciers in the Western Tatra Mts., Western Carpathians

Relict rock glaciers are well-preserved features of high-elevated valleys in the Western Tatra Mts., Western Carpathians, but their chronology remained poorly constrained with numerical dating methods. We present the first robust set of 10Be exposure ages for eight rock glaciers with front elevation of 1376–1819 m asl. The results suggest that the rock glaciers stabilized throughout the Weichselian Lateglacial from ∼16.5 ka to 11 ka. This timing is consistent with the period of rock glacier stabilization in European mountain regions, but it extends the age span previously determined for rock glaciers in the Tatra Mts. There are differences between north- and south-facing valleys in the elevation and time of the stabilization of the rock glaciers, which were probably caused by contrasts in potential incoming solar radiation that affected the retreat of former glaciers and the subsequent formation of the rock glaciers. The initiation altitude of ∼1610 and ∼1830 m asl determined for the oldest and youngest dated rock glaciers is generally consistent with previous estimates of the glacier equilibrium line altitude during the Greenland Stadials 2.1 and 1, respectively. However, the lower limit of the rock glaciers is well above regional paleopermafrost features from the same periods, which suggests that rock glaciers may underestimate past permafrost extents and temperature declines disputing their regional validity for paleoenvironmental reconstructions.

Eocene–Oligocene glaciation on a high central Tibetan Plateau

Abstract The nature of glaciation (bipolar vs. unipolar) during the Eocene–Oligocene transition (EOT) remains unresolved. Here, we report the occurrence of frost marks, ice-rafted debris (IRD), and glendonites from the Upper Eocene to Lower Oligocene Niubao Formation (Fm.) deposited in a proglacial lake above glaciolacustrine conglomerates and diamictite facies in the Lunpola Basin, central Tibetan Plateau (CTP). Magnetostratigraphy dates these cryospheric deposits to ca. 36.2–31.8 Ma, synchronous with a stratigraphic interval containing IRD offshore of SE Greenland and in the Barents, Chukchi, and Laptev Seas, suggesting a strong continental-oceanic coupling. Our results provide robust continental evidence for intermittent cryospheric processes in the midlatitude Northern Hemisphere during the late Eocene and EOT. The global cold snap EOT-1 influenced already glacierized high-altitude mountains, lowering equilibrium line altitudes (ELAs) of glaciers and leading to local development of ice fields, ice caps, and valley glaciers with proglacial lake systems, such as the one recorded in the Niubao Fm. The record of IRD, glendonites, and frost marks before the onset of EOT-1 points to an active cryosphere on a plateau already elevated by ca. 36.2 Ma.

Cirque morphology and palaeo-climate indications along a south-north transect in High Mountain Asia

Cirques preserve important palaeo-climate information and have been investigated in various areas around the world. High Mountain Asia (HMA) contains the largest modern glacier inventory outside of polar regions, and tens of thousands of cirques. We analyzed 2831 largely ice-free cirques from five areas along a south-north transect in HMA, including central Himalayas (C. Himalayas), western Gangdise (W. Gangdise), eastern Tian Shan (E. Tian Shan), southern Altai (S. Altai), and northern Altai (N. Altai). These cirques are large relative to those in many parts of the world, and in the Himalayas, Tian Shan, and Altai their vertical dimensions are considerable. Their size variations show weakly allometric relations. The poleward decreasing trend in cirque floor altitudes is consistent with the reducing temperature trend toward high latitudes. Although just north of C. Himalayas, W. Gangdise is in their precipitation shadow and its cirque floors are higher. Relatively shallow cirques in W. Gangdise probably reflect low subglacial erosion associated with limited precipitation received from either the Indian monsoon or Westerlies due to this orographic effect. Lower cirque floor altitudes for north-facing cirques were likely caused by their lesser receipt of solar radiation. Because most mountains rose well above equilibrium-line altitudes of glaciers and were heavily glaciated, cirque aspects show symmetric vector strengths except in the lowest area, N. Altai. Eastward aspects in N. and S. Altai show the strong effects of Westerly winds in the past, as for present glaciers. Similarly, cirque aspects in C. Himalayas are more northeast-ward than modern glacier aspects, suggesting a stronger west wind influence in the past. North-facing cirques are lower and more frequent than south-facing, showing the general importance of solar radiation receipt, which is especially dominant in E. Tian Shan.

On the Shift of Glacier Equilibrium Line Altitude (ELA) under the Changing Climate

Presently available information on the glacier equilibrium line altitude (ELA) is being collected and examined. The historical course of the world’s longest ELA series of 107 years at the Claridenfirn is reviewed together with climatic elements. Further, the changes in ELAs of 70 glaciers the world over are investigated, and a linear plane model for the speed of the ELA shift is proposed as a function of the changing rates of summer temperature and winter mass balance. The four glaciers in Europe, which diverge most from the plane, are investigated in detail. The cause of the divergence is likely due to be the change in solar global radiation. Although a precise quantification of the role of radiation is not possible at this stage for the entire world, the role of solar radiation is investigated for these glaciers. Globally viewed, ten, or 15% of the 70 investigated glaciers, are expected to lose their accumulation areas within the next ten years. Half of all studied glaciers will follow the same fate by the end of this century under the present climatic conditions. If climate change is accelerated, the disappearance of glaciers will occur sooner than presented in this study.

Little Ice Age climate in southernmost Greenland inferred from quantitative geospatial analyses of alpine glacier reconstructions

Alpine glaciers are sensitive responders to climate fluctuations, especially to changes in summer temperature, and thus their past extents are invaluable indicators of past climate. Here, we leverage geomorphological evidence of past glacier extent to reconstruct the three-dimensional surfaces of 42 land-terminating paleoglaciers in southernmost Greenland and calculate their equilibrium-line altitudes (ELAs) during the presumed coldest interval of the Little Ice Age (LIA). We compare LIA paleoglacier ELAs with observed modern glacier ELAs (i.e. snowlines) to estimate the LIA summer temperature depression relative to present. On average, we find that ELAs rose 122 ± 64 m since the LIA, which corresponds to a summer temperature change of 1.1± 0.6 ° C, assuming no change in precipitation. Because there are minimal constraints on LIA precipitation anomalies in southernmost Greenland, we also report summer temperature changes for precipitation values ranging from a substantially drier to substantially wetter than modern LIA. When compared with independent temperature evidence, our results suggest LIA precipitation was similar to or slightly less than today. The alpine glaciers studied lost (on average) 56± 16% of their area and 44± 13% of their centerline length between the LIA culmination and 2019 CE. In general, smaller and lower elevation glaciers lost larger fractions of their area. We find no clear geographic patterns in the magnitudes of ELA, area, or length changes across the 9000 km2 study area despite its heterogeneous climate. This study demonstrates a workflow for regional-scale paleoglacier reconstruction and provides the first quantitative LIA summer paleotemperature estimate for southernmost Greenland via region-wide glacial evidence. Our work also suggests the need for reconstruction of a large number of glaciers to infer past regional climatic information, as individual glaciers may respond heterogeneously to climate perturbations leading to variable changes in their ELA and length.

The West Kunlun Glacier Anomaly and Its Response to Climate Forcing during 2002–2020

Research into glacial mass change in West Kunlun (WK) has been sufficient, but most of the existing studies were based on geodetic methods, which are not suitable for specific health state analyses of each glacier. In this paper, we utilize Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) imagery, applying the continuity equation to obtain altitudinal specific mass balance (SMB) for 615 glaciers (>2 km2) during 2002–2011, 2011–2020, and 2002–2020 to research glacial health and its response to climatic forcing. The results show dissimilar glacier SMB patterns between 2002–2011 (0.10 ± 0.14 m w.e. a−1), 2011–2020 (–0.12 ± 0.14 m w.e. a−1) and 2002–2020 (−0.01 ± 0.07 m w.e. a−1). Additionally, the glacier equilibrium line altitude (ELA) in WK was 5788 m, 5744 m, and 5786 m, respectively, and the corresponding accumulation area ratios (AARs) were 0.59, 0.62, and 0.58, during 2002–2011, 2011–2020, and 2002–2020, respectively. Regarding glacier response, compared with the ordinary-least-square (OLS) model, the artificial neural network (ANN) model revealed a respectively less and more sensitive glacier SMB response to extreme negative and positive summer skin temperatures. In addition, the ANN model indicated that the glacier ELA was less sensitive when the integrated water vapor transport (IVT) change exceeded 0.7 kg m−1s−1. Moreover, compared with IVT (−121.57 m/kg m−1s−1), glacier ELA shifts were chiefly dominated by summer skin temperature (+154.66 m/℃) in the last two decades. From 2002–2011 and 2011–2020, glacier SMB was more susceptible to summer skin temperature (−0.38 m w.e./℃ and −0.16 m w.e./℃, respectively), while during 2002–2020, it was more influenced by IVT (0.45 m w.e./kg m−1s−1). In contrast with eastern WK, glaciers in western WK were healthier, although mitigation measures are still needed to safeguard glacier health and prevent possible natural hazards in this region. Finally, we believe that the inconsistent change between glacier SMB and ELAs from 2002–2020 was connected with ice rheology and that the combined effects of skin temperature and IVT can explain the WK glacier anomaly.

Controls on the altitude of Scandinavian cirques: What do they tell us about palaeoclimate?

Cirques are glacially eroded, bowl-shaped depressions, characterised by steep headwalls and flat or overdeepened floors. Given their association with past glaciers, cirques are sometimes used as proxies for palaeoclimate. However, cirques are shaped over multiple glacial cycles, and their usefulness as palaeoclimate indicators therefore remains open to question. In this paper, we map 3984 glacier-free cirques across the Scandinavian Peninsula and analyse variations in cirque floor altitude (CFA). We explore the relationships between CFAs and cirque aspect, latitude, longitude, and distance to the coast. We test the validity of using CFAs as indicators of palaeoclimate through comparison with the equilibrium-line altitudes (ELAs) of 513 modern cirque glaciers. Results indicate that both CFAs and modern cirque-glacier ELAs decrease with latitude and vary with aspect, being generally lowest on east-facing slopes. However, the clearest and strongest trend in both CFAs and modern cirque glacier ELAs is an increase in elevation with distance from the modern coast (i.e., distance ‘inland’). This likely indicates that similar climatic gradients, particularly an inland reduction in precipitation, acted to regulate former sites of glacier initiation (reflected by CFAs) and modern glacier ELAs. This would imply that CFAs are a useful proxy for palaeoclimate. However, we note that both CFAs and modern ELAs reflect the general topography of this region (with increasing elevations moving inland), and the glacial history of the area (indirectly linked to palaeoclimate) may have played a role in regulating where cirques have formed. For these reasons, we suggest that palaeoclimatic interpretations derived from CFAs should be treated with caution.

Glacial equilibrium line–based relationship for paleoclimate reconstructions (Sierra Nevada, USA)

AbstractWe use PRISM climatic data (1981–2010) and Landsat images (2012–2013) to establish an empirical relationship linking annual temperature and precipitation to the equilibrium line altitude (ELA) of glaciers in the Sierra Nevada (36–41°N, California, USA). For this, we determined the present-day ELAs of 57 glaciers and the local 0°C isotherms elevation Iso0, averaged over the 1981–2010 period. The difference, for each glacier, is Y, the normalized snowline altitude (Y = ELA – Iso0). We then empirically calibrated a logarithmic relationship between this normalized snowline altitude and mean annual precipitation using data from partially covered glaciers. Our calibration is statistically distinct from that previously established for the tropical and midlatitude Andes (Fox and Bloom [1994], Journal of Geography (Chigaku Zasshi), 103, 867–885; Condom et al. [2007], Global and Planetary Change, 59, 189–202). This new relationship for North America is an easy-to-use tool to permit paleoclimatic reconstructions from paleo-ELAs. For a specific paleoglacial site, paleotemperature can be computed knowing the paleoprecipitation range, and vice versa. We also performed a test showing that, if precipitation is well known, the uncertainty associated with paleotemperature is about 1°C (1σ).

Recent Spatiotemporal Trends in Glacier Snowline Altitude at the End of the Melt Season in the Qilian Mountains, China

Glaciers in the Qilian Mountains, China, play an important role in supplying freshwater to downstream populations, maintaining ecological balance, and supporting economic development on the Tibetan Plateau. Glacier snowline altitude (SLA) at the end of the melt season is an indicator of the Equilibrium line altitude (ELA), and can be used to estimate the mass balance and climate reconstruction. Here, we employ the height zone-area method to determine the SLA at the end of the melt season during the 1989–2018 period using Landsat, MODIS (Moderate Resolution Imaging Spectroradiometer) SLA and Shuttle Radar Topography Mission (SRTM) digital elevation model (DEM) data. The accuracy of glacier SLA obtained in 1989–2018 after adding MODIS SLA data to the years without Landsat data increased by about 78 m. The difference between the remote-sensing-derived SLA and measured equilibrium line altitude (ELA) is mostly within 50 m, suggesting that the SLA can serve as a proxy for the ELA at the end of the melt season. The SLA of Qiyi Glacier in the Qilian Mountains rose from 4690 ± 25 m to 5030 ± 25 m, with an average of 4900 ± 103 m during the 30 year study period. The western, central, eastern sections and the whole range of the Qilian Mountains exhibited an upward trend in SLA during the 30 year study period. The mean glacier SLAs were 4923 ± 137 m, 4864 ± 135 m, 4550 ± 149 m and 4779 ± 149 m for the western, central, eastern sections and the whole range, respectively. From the perspective of spatial distribution, regardless of the different orientation, grid scale and basin scale, the glacier SLA of Qilian Mountains showed an upward trend from 1989 to 2018, and the glacier SLA is in general located at a comparably higher altitude in the southern and western parts of the Qilian Mountains while it is located at a comparably lower altitude in its northern and eastern parts. In an ideal condition, climate sensitivity studies of ELA in Qilian Mountains show that if the summer mean temperature increases (decreases) by 1 °C, then ELA will increase (decrease) by about 102 m. If the annual total solid precipitation increases (decreases) by 10%, then the glacier ELA will decrease (rise) by about 6 m. The summer mean temperature is the main factor affecting the temporal trend of SLA, whereas both summer mean temperature and annual total precipitation influence the spatial change of SLA.

Testing the area–altitude balance ratio (AABR) and accumulation–area ratio (AAR) methods of calculating glacier equilibrium-line altitudes

AbstractIn this study, we compare equilibrium-line altitudes (ELAs) calculated using the area–altitude balance ratio (AABR) and the accumulation–area ratio (AAR) methods, with measured ELAs derived from direct field observations. We utilise a GIS toolbox to calculate the ELA for 64 extant glaciers by applying the AABR and AAR methods to DEMs and polygons of their geometry. The calculated ELAs (c-ELAs) are then compared to measured zero-net balance ELAs (znb-ELAs) obtained from mass-balance time series held by the WGMS for the same glaciers. The correlation between znb-ELAs and AABR (1.56)/AAR (0.58) c-ELAs is very strong, with an r2 = 0.99. The smallest median difference between znb-ELAs and c-ELAs (i.e. 65.5 m) is obtained when a globally representative AABR of 1.56 is used. When applied to palaeoglacier-climate applications, this difference translates to ~0.42°C, well within the uncertainty of palaeotemperature proxies used to determine mean summer temperature at the ELA. The more widely used mean AABR of 1.75 is shown to be statistically invalid due to the skewness of the dataset. On this basis, when calculating glacier ELAs, we recommend the use of a global AABR value of 1.56.

Climatic controls on mountain glacier basal thermal regimes dictate spatial patterns of glacial erosion

Abstract. Climate has been viewed as a primary control on the rates and patterns of glacial erosion, yet our understanding of the mechanisms by which climate influences glacial erosion is limited. We hypothesize that climate controls the patterns of glacial erosion by altering the basal thermal regime of glaciers. The basal thermal regime is a first-order control on the spatial patterns of glacial erosion. Polythermal glaciers contain both cold-based portions that protect bedrock from erosion and warm-based portions that actively erode bedrock. In this study, we model the impact of various climatic conditions on glacier basal thermal regimes and patterns of glacial erosion in mountainous regions. We couple a sliding-dependent glacial erosion model with the Parallel Ice Sheet Model (PISM) to simulate the evolution of the glacier basal thermal regime and glacial erosion in a synthetic landscape. We find that both basal thermal regimes and glacial erosion patterns are sensitive to climatic conditions, and glacial erosion patterns follow the patterns of the basal thermal regime. Cold temperature leads to limited glacial erosion at high elevations due to cold-based conditions. Increasing precipitation can overcome the impact of cold temperature on the basal thermal regime by accumulating thick ice and lowering the melting point of ice at the base of glaciers. High precipitation rates, therefore, tend to cause warm-based conditions at high elevations, resulting in intensive erosion near the peak of the mountain range. Previous studies often assessed the impact of climate on the spatial patterns of glacial erosion by integrating climatic conditions into the equilibrium line altitudes (ELAs) of glaciers, and glacial erosion is suggested to be maximal around the ELA. However, our results show that different climatic conditions produce glaciers with similar ELAs but different patterns of basal thermal regime and glacial erosion, suggesting that there might not be any direct correlation between ELAs and glacial erosion patterns.

Geomorphological record and equilibrium line altitude of glaciers during the last glacial maximum in the Rodna Mountains (eastern Carpathians) – Corrigendum

Geomorphological record and equilibrium line altitude of glaciers during the last glacial maximum in the Rodna Mountains (eastern Carpathians) – Corrigendum

Climatological analysis of freezing level height over China and its implications using homogenized in‐situ data

ABSTRACTFreezing level height (FLH) is the altitude of a location where the air temperature is 0°C. The FLH denotes the approximate position for water phase transitions from liquid (or vapour) to solid but not vapour to liquid. Based on temperature and height data on the ground surface to heights with an atmospheric pressure of 400 hPa, the monthly FLH at 118 radiosonde stations located in China are derived for the 1958–2017 period. The monthly and annual climatology during 1981–2010 and the seasonal, intra‐annual and decadal variabilities of the FLH during 1958–2017, along with its long‐term trends, were quantitatively evaluated at the regional scale and for individual stations. FLH climatology is generally related to season, zone and topography. The annual FLH anomalies decline initially, followed by a rise, with synchronous changes in the temperature in the lower and middle troposphere; the anomalies in south‐eastern China are partly associated with El Niño–Southern Oscillation events. The increasing trends of FLH are closely associated with tropospheric warming and show evident regional and seasonal differences. Stronger increasing trends are observed in northern and western China in autumn or winter (particularly in recent decades). The changes in FLH have implications on long‐term trends in equilibrium line altitude (ELA) of glaciers in Tianshan Mountains, snow cover extent in melt season and hail frequency in China especially in high‐elevations. During recent decades, the evolution of the FLH and major composition in the cryosphere has accelerated as a response to tropospheric warming. The study on FLH change can provide valuable information for the evaluation of glaciers, snow cover extent and hail frequency in China from a new perspective; which is helpful for understanding regional hydrological processes and climate changes better. Further application of FLH in cryosphere monitoring can improve water management and flood forecasting in headwater regions in China.

Adjustment of glacier geometries to future equilibrium line altitudes in the Shigar River Basin of Karakoram Range, Pakistan

Pakistan’s Karakoram region has a large variety of glacier types. Equilibrium line altitude (ELA) of alpine or valley glaciers represents mass balance. Field observations for estimation of ELA of the majority of Karakoram’s glaciers are not available due to rugged glacier-covered terrains and lack of climatic data above 5000 masl. Therefore, we applied the hypsometrically controlled accumulation area ratio (AAR) and accumulation area balance ratio (AABR) methods for ELA estimation and glacier–climate reconstructions in the Shigar River Basin of the Karakoram region. Constrained by mountain topography, larger size and type of glaciers, several ranges of ELA are calculated and implemented for several ratios. Two parameters (ratio and interval) are provided to calculate AAR-ELAs between 0.4 and 0.8 with 0.05 interval and AABR-ELAs between 0.9 and 4.4 with 0.01 interval. By providing constant AAR (rather than constant glacier area), this approach adjusted glacier geometries (area) to future ELA variations. For constant AAR of 0.4–0.45, a 90-m ELA decrease from 5769 to 5679 m of Baltoro glacier adjusted its geometries by reducing ~ 5% area. The highest decrease of 140-m ELA of the same glacier is reported for constant AAR of 0.7–0.75, revealing a significant loss of 8% geometries. The projected geometry losses for all these glaciers are highly variable, with top-heavy glaciers (Biafo and Baltoro) projected to experience the major losses in glacier-ice area. It is concluded that the quality of ELAs is highly dependent on the reconstructed three-dimensional glacier surfaces.

Spatiotemporal Variability in the Glacier Snowline Altitude across High Mountain Asia and Potential Driving Factors

The glacier snowline altitude (SLA) at the end of the melt season is an indicator of the glacier equilibrium line altitude and can be used to estimate glacier mass balance and reconstruct past climate. This study analyzes the spatiotemporal variability in glacier SLA across High Mountain Asia, including the Altai Mountains, Karakoram Mountains, Western Himalayas, Gongga Mountains, Tian Shan, and Nyainqentanglha Mountains, over the past 30 years (1989–2019) to better elucidate the state of these mountain glaciers. Remote-sensing data are processed to delineate the glacier SLA across these mountainous regions, and nearby weather station data are incorporated to determine the potential relationships between SLA and temperature/precipitation. The mean SLAs across the Altai and Karakoram mountains ranged from 2860 ± 169 m to 3200 ± 152 m and from 5120 ± 159 m to 5320 ± 240 m, respectively, with both regions experiencing an average increase of up to 137 m over the past 30 years. Furthermore, the mean glacier SLAs across the Western Himalayas and Gongga Mountains increased by 190–282 m over the past 30 years, with both regions experiencing large fluctuations. In particular, the mean glacier SLA across the Western Himalayas varied from 4910 ± 190 m in 1989 to 5380 ± 164 m in 2000, and that across the Gongga Mountains varied from 4960 ± 70 m in 1989 to 5330 ± 100 m in 2012. Correlation analyses between glacier SLA and temperature/precipitation suggest that temperature is the primary factor influencing glacier SLA across these High Mountain Asia glaciers. There is a broad increase in glacier SLA from the Altai Mountains to the Karakoram Mountains, with a decrease in glacier SLA with decreasing latitude across the Himalayas; the maximum SLA occurs near the northern slopes of the Western Himalayas. The glacier SLA is lower on the eastern side of the Tibetan Plateau and exhibits a longitudinal distribution pattern. These results are expected to provide useful information for evaluating the state of High Mountain Asia glaciers, as well as their response and feedback to climate change.

A first-order flexible ELA model based on geomorphic constraints

Alpine glaciers, with their valuable combination of highly sensitive response to climate and near-global extent, are powerful tools for investigating previous and present climate changes. They also represent critical water resources for areas around the globe, with the potential for far-reaching effects in a warming world. Advancements to understand and model glacial changes and the variables influencing them are therefore paramount. Many glacier models fall into one of two endmembers: either highly complex transient models requiring careful tuning of multiple parameters to individual glaciers, or basic empirical correlations of glacier area and length with few considerations for local and regional variations in characteristics. Here we detail a physical steady-state model for alpine glaciers relating directly to glacier mass balance (via the equilibrium line altitude) while retaining the simplicity of other morphology methods, and simultaneously including error estimates. We provide custom MATLAB functions as a user-friendly and generally-applicable method to estimate glacier equilibrium line altitudes from only a limited number of glacier bed topography and glacier width measurements. As a test of the model’s efficacy, we compare the model results for present-day glaciers in the Swiss Alps with previously published estimates of equilibrium line altitudes and intermediate model outputs.•The method estimates glacier equilibrium line altitudes from a limited set of bed topography measurements and constraints on glacier width.•The method is based on continuity equations, reducing the need for empirical coefficients tuned with measured data.•The method uses Monte Carlo sampling and bootstrapping to generate uncertainty bounds on the equilibrium line altitude estimates.

Numerical Reconstruction of Three Holocene Glacial Events in Qiangyong Valley, Southern Tibetan Plateau and Their Implication for Holocene Climate Changes

The dating of well-preserved Holocene moraines in the Qiangyong Valley, southern Tibetan Plateau (TP), offers great potential for reconstructing Holocene glacier extents and examining climate changes in the region. Guided by Holocene moraine features, this study used Geographic Information System (GIS) model tools to reconstruct paleo-glacier surfaces and glacier equilibrium line altitude (ELA) depressions for three Holocene glacial stages in the valley. The GIS-based models showed that the Qiangyong Valley contained ice volumes of 8.1 × 108, 6.2 × 108, and 4.6 × 108 m3 during the early Holocene, Neoglacial, and Little Ice Age (LIA) glacial stages, and that the ELA was decreased by ~230 ± 25, ~210 ± 25, and ~165 ± 25 m, respectively, compared to modern conditions. Furthermore, the summer temperatures were estimated to be 1.56–1.79, 1.37–1.64, and 1.29–1.32 °C cooler than present to support the three Holocene glacier extents, based on the evidence that the respective precipitation increased by 20–98, 13–109, and 0.9–11 mm relative to the present, which were derived from the lacustrine pollen data for the southern TP. By comparison, this study found that the amplitudes of the ELA-based summer temperature depressions were much larger than the pollen-based counterparts for the three glacial stages, although the two proxies both showed increasing trends in the reconstructed summer temperatures.