Search

The glacier surface is composed not only of ice or snow but also of a heterogeneous mixture of various materials. The presence of light-absorbing impurities darkens the glacier surface, reducing local reflectance and thereby accelerating the glacier melting process. However, our understanding of the spatial distribution of these impurities remains limited, and there is a lack of studies on quantifying the dirty degree of glacier surfaces. During the Sentinel satellite overpass on 21 August 2023, we used an ASD FieldSpec3 spectrometer to measure the reflectance spectra of glacier surfaces with varying degrees of dirtiness on the Qiyi glacier, Qinghai–Tibet Plateau. Using Multiple Endmember Spectral Mixture Analysis (MESMA), the Sentinel imagery was decomposed to generate fraction images of five primary ice surface materials as follows: coarse-grained snow, slightly dirty ice, moderately dirty ice, extremely dirty ice, and debris. Using unmanned aerial vehicle (UAV) imagery with a 0.05 m resolution, the primary ice surface was delineated and utilized as reference data to validate the fraction images. The findings revealed a strong correlation between the fraction images and the reference data (R2 ≥ 0.66, RMSE ≤ 0.21). Based on pixel-based classification from the UAV imagery, approximately 80% of the glacier surface is covered by slightly dirty ice (19.2%), moderately dirty ice (33.3%), extremely dirty ice (26.3%), and debris (1.2%), which significantly contributes to its darkening. Our study demonstrates the effectiveness of using Sentinel imagery in conjunction with MESMA to map the degree of glacier surface dirtiness accurately. [ABSTRACT FROM AUTHOR]

The glaciers in southeastern Tibet Plateau (SETP) influenced by oceanic climate are sensitive to global warming, and there remains a notable deficiency in accurate multitemporal change analyses of these glaciers. We conduct glacier inventories in the Yigong Zangbo River Basin (YZRB) in SETP for the years 1988, 2015, and 2023 utilizing Landsat and Sentinel-2 imagery, and analyze the glacier spatiotemporal variation incorporating the existing glacier inventory data. Since the 1970s until 2023, the glaciers significantly retreated at a rate of 0.76 ± 0.11%·a−1, with the area decreasing from 2583.09 ± 88.80 km2 to 1635.89 ± 71.74 km2, and the ice volume reducing from 221.7017 ± 7.9618 km3 to 152.7429 ± 6.1747 km3. The most significant retreat occurred in glaciers smaller than 1 km2. Additionally, glaciers on southern aspects retreated slower than the northern counterparts. The glaciers in the western YZRB witnessed a significantly greater shrinkage rate than those in the eastern section, with the most pronounced changes occurring in Aso Longbu River Basin. Furthermore, severe glacier mass deficits were observed from 2000 to 2019, averaging a loss rate of 0.57 ± 0.06 m w.e. a−1. The continuous rise in air temperature has primarily induced a general widespread glacier change in the YZRB. However, diverse topography led to spatial variability in glacier changes with discrepancies as large as several times. The features of individual glaciers, such as glacier size, debris cover, and the development of ice-contact glacial lakes enhanced the local complexity of glacier change and elusive response behaviors to climate warming led by the different topographic conditions. [ABSTRACT FROM AUTHOR]

Tibet's Qilian Mountains (QM) include critical water conservation areas and important ecological barriers, which help maintain downstream inland river oasis stability. For decades, contemporaneous climatic and cryospheric variation has severely impacted QM's hydrological processes, challenging local water resource management and sustainable development. However, due to prevalent data and methodological limitations, QM research has primarily focused on runoff change at a basin scale. Spatial distributions and temporal changes in runoff subsequently remain unclear. Based on multi‐source data and the literature, we estimated that QM's mountain outlets generate approximately 15.671 km3 in total annual runoff, exhibiting a spatially decreasing pattern from northeast to southwest. Moreover, runoff distribution and trend variation at seasonal and annual scales depend upon the river replenishment source type. Beginning in the 1950s and 1960s, eastern rain‐fed rivers experienced a downward trend while those dominated by meltwater or simultaneously fed by multiple sources in its central and western regions experienced an upward trend. As an integrated product of mixed multi‐factor effects, runoff is regulated by temperature, precipitation, and cryospheric meltwater. Moreover, the main controlling runoff factors varied seasonally under different water source concentrations. Annually, precipitation was the main driver for runoff change in the eastern region while, correspondingly, temperature was in the western region where glaciers and the snow line boundary predominant. Besides, this study highlighted that the existing literature has significant limitations in understanding interactions among different cryospheric components and hydrologic process mechanisms when exploring the reasons for runoff variations, which needs further exploration in the future. Plain Language Summary: As a critical water conservation area and ecological security barrier in western China, the Qilian Mountains (QM) is a core region within the Silk Road Economic Belt. It is regarded as the "water tower" to inland rivers in the Hexi Corridor, the eastern margin of the Qaidam Basin desert, and upstream areas of the Yellow River, nurturing more than 5 million people. However, spatiotemporal patterns and driving factors of runoff variations under climate change conditions remain unclear. Based on multi‐source data and published literature, this study estimated that 15.671 km3 of total annual runoff was generated on the QM, which generally decreased from northeast to southwest. Although seasonal distributions of runoff in major rivers were similar, their annual runoff trends showed a substantial difference between the eastern (i.e., downward trend) and western (i.e., upward trend) regions. From east to center to west, river replenishment types changed from rainfall dominant to mixed replenishment/groundwater dominant to meltwater dominant, causing different sensitivities to temperature and precipitation among rivers. Currently, the main approaches exploring physical mechanisms of QM runoff are qualitative analyses. Comprehensive quantitative research that encompasses various cryospheric components is needed to enhance scientific understanding of hydrological processes in this region. Key Points: Average annual total river runoff on the Qilian Mountains was estimated at approximately 15.671 km3, decreasing from northeast to southwestA downward trend in annual runoff was observed in the eastern region, while an upward trend was observed in its central and western regionsRiver replenishment types changed from rainfall dominant (eastern) to mixed replenishment (central) to meltwater dominant (western) [ABSTRACT FROM AUTHOR]

Lake-terminating glaciers have some peculiar behaviors compared to land-terminating glaciers, but in-depth observation is still limited regarding their formation, which is crucial for understanding the glacier–lake interaction. Here, the long-term evolutions of Tanymas Lake and the Fedchenko Glacier were investigated based on Landsat images, Google Earth imagery, KH-9 images, glacier surface elevation and velocity change datasets, and meteorological records. The results indicate that Tanymas Lake is both an ice-contact proglacial lake and an ice-dammed lake. It covered an area of 1.10 km2 in September 2022, and it is one of the largest glacial lakes in Pamir and even in HMA. The initial basin of Tanymas Lake is a moraine depression in Tanymas Pass, and the blocked dam is the Tanymas-5 Glacier and its terminal moraine. Tanymas Lake was in an embryonic stage before August 2005, in a formation and expansion stage from August 2005 to September 2018, and in a new expansion stage after September 2018. In this process, the Tanymas terminus of the Fedchenko Glacier also transformed from a land terminus to a partial lake terminus, and then to a complete lake terminus. The formation of Tanymas Lake is associated with the accumulation of glacial meltwater and the blockage of drainage, while the slow expansion of Tanymas Lake is related to the cold climate and slight glacier mass loss of Central Pamir. In the coming decades, with the accelerated mass loss of the Tanymas terminus of the Fedchenko Glacier, the area, depth, and water storage of Tanymas Lake will continue to increase, accompanied by the growing GLOF risk. [ABSTRACT FROM AUTHOR]

Synthetic Aperture Radar images have recently been utilized in glacier surface flow velocity research due to their continuously improving imaging technology, which increases the resolution and scope of research. In this study, we employed the offset tracking and multidimensional small baseline subset (MSBAS) technique to extract the surface flow velocity of the Siachen Glacier from 253 Sentinel-1 images. From 2017 to 2021, the Siachen Glacier had an average flow velocity of 38.25 m a−1, with the highest flow velocity of 353.35 m a−1 located in the upper part of a tributary due to the steep slope and narrow valley. The inter-annual flow velocity fluctuations show visible seasonal patterns, with the highest flow velocity observed between May and July and the lowest between December and January. Mass balance calculated by the geodetic method based on AST14DEM indicates that the Siachen Glacier experienced a positive mass change (0.07 ± 0.23 m w.e. a−1) between 2008 and 2021. However, there was significant spatial heterogeneity revealed in the distribution, with surface elevation changes showing a decrease in the glacier tongue while thickness increased in two other western tributaries of the Siachen Glacier. The non-surface parallel flow component is correlated with the strain rate and mass balance process, and correlation analysis indicates a positive agreement between these two variables. Therefore, using glacier flow velocities obtained from the SAR approach, we can evaluate the health of the glacier and obtain crucial factors for the glacier's dynamic model. Two western tributaries of the Siachen Glacier experienced mass gain in the past two decades, necessitating close monitoring of flow velocity changes in the future to detect potential glacier surges. [ABSTRACT FROM AUTHOR]

Under global warming, many glaciers worldwide are receding. However, recent studies have suggested the extension of the Karakoram Anomaly, a region of anomalous glacier mass gain, into the western Kunlun and eastern Pamir mountains. However, the eastern limit of this anomaly in the Kunlun Mountains is unclear. This study, using changes in glacier area and surface elevation, estimates the eastern limit of the Kunlun-Pamir-Karakoram anomaly at ~85°E. Over the past 50 years, glaciers west of 85°E in the Kunlun Mountains decreased in area from 8401 to 7945 km2 at a rate of −0.12 ± 0.07% a−1, showed a reduction in the rate of retreat through time and have recently gained mass, with surface elevation changes of 0.15 ± 0.35 m a−1 over the period of 2000–2013. Glaciers east of 85°E have experienced greater rates of area change (−61 ± 12 km2 and −0.43 ± 0.13% a−1) over the past 50 years, accelerated area loss in recent years and elevation change rate of −0.51 ± 0.18 m a−1 between 2000 and 2013. These patterns of elevation and area change are consistent with regional increases in summer temperature in the eastern Kunlun Mountains and slight cooling in the western Kunlun Mountains. [ABSTRACT FROM AUTHOR]

Palaeoglacier modelling is an important approach for reconstructing the palaeoclimate. The timing of glaciations in the central part of the Gangdise Mountains has been constrained previously, but the palaeoclimate remains unclear. In this paper, the palaeo-temperature and precipitation of the early marine isotope stage (MIS) 2, the Last Glacial Maximum (LGM), and the early Holocene were reconstructed using coupled mass balance and ice flow models. The results show that a series of temperature changes (ΔT) and precipitation factors (Fp) resulted in optimum palaeoglacial extents. The modelled palaeoglaciers during the early MIS 2, the LGM, and the early Holocene cover areas of ~18.1 km2, ~17.4 km2, and ~16.3 km2, respectively, with ice volumes of ~2.18 km3, ~1.99 km3, and ~1.95 km3, respectively. Previous studies on ice cores, pollen samples, and lake sediments were referenced to narrow the range of palaeo-temperatures and precipitations. The reconstructed temperatures during the early MIS 2, LGM, and early Holocene were constrained to 2.4–2.9 °C, 2.15–3.05 °C, and 0.95–1.5 °C lower than today, respectively. Their precipitation levels were 60–80%, 50–80%, and 100–150% of the present-day level, respectively. [ABSTRACT FROM AUTHOR]

In the context of global warming, glaciers in the Asian High Mountains (AHMs) are shrinking at an accelerating rate. Projecting their future change is helpful for understanding the hydrological and climatic effects related to glacier retreat. Here, we projected glacier change in the AHMs from 1979 to 2100 under shared socioeconomic pathway (SSP) scenarios from the perspective of temperature, equilibrium-line altitude (ELA), and accumulation area. The annual mean temperature in the AHMs increased by 1.26°C from 1979 to 2014, corresponding to an increase of 210 m in the mean ELA and a decrease of 1.7×104 km2 in the glacier accumulation area. Under the SSP2-4.5 (SSP5-8.5) scenario, the annual mean temperature in the AHMs would increase by 2.84°C (3.38°C) in 2040–2060 relative to that in 1850–1900, leading to the mean ELA reaching an elevation of 5661 m (5777 m). The accumulation area in the AHMs decreased by 46.3% from 1995 to 2014 and was projected to decrease by 60.1% in 2040–2060. Moreover, the annual mean temperature in the AHMs was projected to increase by 3.76°C (6.44°C) in 2080–2100 relative to that in 1850–1900, corresponding to the ELA reaching an elevation of 5821 m (6245 m) and the accumulation area decreasing to 1.8×104 km2 (0.5×104 km2). These data suggest that the conditions for glacier development will disappear in most of the AHMs, except for extreme high-altitude regions in the Tianshan, Pamir, and Himalaya Mountains. Under the SSP2-4.5 (SSP5-8.5) scenario, when the global mean temperature increases 1.5°C (2°C) above pre-industrial levels, the annual mean temperature will increase by 2.12°C (2.86°C) and the accumulation area will decrease by 15% (48%) in the AHMs compared with that in 1995–2015. Therefore, a 1.5°C increase in global warming would keep 40% more of the glacial accumulation area (1.5×104 km2) in the AHMs compared to a 2°C increase in global warming. [ABSTRACT FROM AUTHOR]

Glacier snow line altitude (SLA) at the end of the ablation season is an indicator of the equilibrium line altitude (ELA), which is a key parameter for calculating and assessing glacier mass balance. Here, we present an automated algorithm to classify bare ice and snow cover on glaciers using Landsat series images and calculate the minimum annual glacier snow cover ratio (SCR) and maximum SLA for reference glaciers during the 1985–2020 period in Google Earth Engine. The calculated SCR and SLA values are verified using the observed glacier accumulation area ratio (AAR) and ELA. We select 14 reference glaciers from High Mountain Asia (HMA), the Caucasus, the Alps, and Western Canada, which represent four mountainous regions with extensive glacial development in the northern hemisphere. The SLA accuracy is ~73%, with a mean uncertainty of ±24 m, for 13 of the reference glaciers. Eight of these glaciers yield R2 > 0.5, and the other five glaciers yield R2 > 0.3 for their respective SCR–AAR relationships. Furthermore, 10 of these glaciers yield R2 > 0.5 and the other three glaciers yield R2 > 0.3 for their respective SLA–ELA relationships, which indicate that the calculated SLA from this algorithm provides a good fit to the ELA observations. However, Careser Glacier yields a poor fit between the SLA calculations and ELA observations owing to tremendous surface area changes during the analyzed time series; this indicates that glacier surface shape changes due to intense ablation will lead to a misclassification of the glacier surface, resulting in deviations between the SLA and ELA. Furthermore, cloud cover, shadows, and the Otsu method limitation will further affect the SLA calculation. The post-2000 SLA values are better than those obtained before 2000 because merging the Landsat series images reduces the temporal resolution, which allows the date of the calculated SLA to be closer to the date of the observed ELA. From a regional perspective, the glaciers in the Caucasus, HMA and the Alps yield better results than those in Western Canada. This algorithm can be applied to large regions, such as HMA, to obtain snow line estimates where manual approaches are exhaustive and/or unfeasible. Furthermore, new optical data, such as that from Sentinel-2, can be incorporated to further improve the algorithm results. [ABSTRACT FROM AUTHOR]

Monitoring the variations in terrestrial water storage (TWS) is crucial for understanding the regional hydrological processes, which helps to allocate and manage basin-scale water resources efficiently. In this study, the impacts of climate change, glacier mass loss, and human activities on the variations in TWS of the Qaidam Basin over the period of 2002−2020 were investigated by using Gravity Recovery and Climate Experiment (GRACE) and GRACE Follow-On (GRACE-FO) data, and other hydrological and meteorological data. The results indicate that TWS anomalies (TWSA) derived from five GRACE solutions experienced significant increasing trends over the study period, with the change rates ranging from 4.85 to 6.90 mm/year (1.37 to 1.95 km3/year). The GRACE TWSA averaged from different GRACE solutions exhibited an increase at a rate of 5.83 ± 0.12 mm/year (1.65 ± 0.03 km3/year). Trends in individual components of TWS indicate that the increase in soil moisture (7.65 mm/year) contributed the most to the variations in TWS. Through comprehensive analysis, it was found that the temporal variations in TWS of the Qaidam Basin were dominated by the variations in precipitation, and the spatial variations in TWS of the Qaidam Basin were mostly driven by the increase in glacier meltwater due to climate warming, particularly in the Narin Gol Basin. In addition, the water consumption associated with human activities had relatively fewer impacts. [ABSTRACT FROM AUTHOR]

Physical erosion and chemical weathering rates beneath glaciers are expected to increase in a warming climate with enhanced melting but are poorly constrained. We present a global dataset of cations in meltwaters of 77 glaciers, including new data from 19 Asian glaciers. Our study shows that contemporary cation denudation rates (CDRs) beneath glaciers (2174 ± 977 Σ*meq+ m−2 year−1) are ~3 times higher than two decades ago, up to 10 times higher than ice sheet catchments (~150-2000 Σ*meq+ m−2 year−1), up to 50 times higher than whole ice sheet means (~30-45 Σ*meq+ m−2 year−1) and ~4 times higher than major non-glacial riverine means (~500 Σ*meq+ m−2 year−1). Glacial CDRs are positively correlated with air temperature, suggesting glacial chemical weathering yields are likely to increase in future. Our findings highlight that chemical weathering beneath glaciers is more intense than many other terrestrial systems and may become increasingly important for regional biogeochemical cycles. Global glacial chemical denudation is one of the largest contributors to global elemental cycles and, amplified by climate warming, will significantly impact nutrient loads in downstream ecosystems. [ABSTRACT FROM AUTHOR]

Palaeoflood records have provided the extended perspective required to understand linkages between extreme floods and climate variability. During our palaeo‐hydrological investigations, two units of palaeoflood slackwater deposits (SWD) were found on the river banks in the Maqu‐Oula reach along the mainstream of the Yellow River on the northeast (NE) Tibetan Plateau. A series of samples were taken from the sediment profiles at Taiwuruo (TWR) and Daerqiong (DEQ) sites for laboratory analysis. Grain‐size distribution and micro‐morphological characteristics of quartz grains indicate that these well‐sorted SWD units were derived from the suspended sediment load in floodwater. They have recorded two episodes of overbank flooding in the Yellow River source area. Using the optically stimulated luminescence (OSL) dating technique, the first episode of overbank flooding was dated to 13.6–13.0 ka and correlated with the Bølling–Allerød warm interval of the last deglaciation. The second episode occurred at 11.0–10.0 ka, that is, the early Holocene. In a warm‐humid climate, these overbank flooding episodes may have resulted from the accelerated melting of mountain glaciers and snow over the headwaters of the alpine and subalpine regions on the NE Tibetan Plateau. These results are important in facilitating an in‐depth understanding of the hydro‐climatic response of the Tibetan Plateau to global climate warming occurring presently. [ABSTRACT FROM AUTHOR]

The Tibetan Plateau is a high-elevation area in Asia and contains the largest volumes of glaciers outside the polar regions. Reconstruction of the glacier surface temperature history in this area is crucial for better understanding of the process of climate change in the Tibetan Plateau. The Tikhonov regularization method has been used on borehole temperatures measured at Malan Glacier, located on the north Tibetan Plateau, to reconstruct the surface temperature history in the twentieth century. We found that the glacier surface temperature, which rose significantly after the 1960s, increased about 1.1°C over the last century. The warmest period occurred in the 1980s to the 1990s and the highest temperature variation could be 1.5°C to 1.6°C. The results were also compared with those of nearby instrumental temperatures by Wudaoliang meteorological station and the stable oxygen isotope ( δ 18 O) from the Malan ice core. [ABSTRACT FROM AUTHOR]

Glacier is the dominant cold habitat in terrestrial environments, providing a model ecosystem to explore extremophilic strategies and study early lives on Earth. The dominant form of life in glaciers is bacteria. However, little is known about past evolutionary processes that bacteria underwent during adaptation to the cryosphere and the connection of their genomic traits to environmental stressors. Aiming to test the hypothesis that bacterial genomic content and dynamics are driven by glacial environmental stressors, we compared genomes of 21 psychrophilic Cryobacterium strains, including 14 that we isolated from three Tibetan ice cores, to their mesophilic counterparts from the same family Microbacteriaceae of Actinobacteria. The results show that psychrophilic Cryobacterium underwent more dynamic changes in genome content, and their genomes have a significantly higher number of genes involved in stress response, motility, and chemotaxis than their mesophilic counterparts (P < 0.05). The phylogenetic birth-and-death model imposed on the phylogenomic tree indicates a vast surge in recent common ancestor of psychrophilic Cryobacterium (gained the greatest number of genes by 1,168) after the division of the mesophilic strain Cryobacterium mesophilum. The expansion in genome content brought in key genes primarily of the categories "cofactors, vitamins, prosthetic groups, pigments," "monosaccharides metabolism," and "membrane transport." The amino acid substitution rates of psychrophilic Cryobacterium strains are two orders of magnitude lower than those in mesophilic strains. However, no significantly higher number of cold shock genes was found in psychrophilic Cryobacterium strains, indicating that multi-copy is not a key factor for cold adaptation in the family Microbacteriaceae, although cold shock genes are indispensable for psychrophiles. Extensive gene acquisition and low amino acid substitution rate might be the strategies of psychrophilic Cryobacterium to resist low temperature, oligotrophy, and high UV radiation on glaciers. The exploration of genome evolution and survival strategies of psychrophilic Cryobacterium deepens our understanding of bacterial cold adaptation. [ABSTRACT FROM AUTHOR]

In this study, we collected 60 snowpit samples in nine glaciers from the northern to the southern Tibetan Plateau (TP), to study the levels and spatial distributions of levoglucosan and dissolved organic carbon (DOC). The lowest concentration of levoglucosan was found in the Yuzhufeng (YZF) glacier with a mean value of 0.24 ± 0.08 ng mL − 1 , while the highest concentration of levoglucosan was detected in the Gurenhekou (GRHK) glacier with a mean value of 11.72 ± 15.61 ng mL − 1 . However, the average DOC concentration in TP glaciers were comparable, without significant regional differences. The levoglucosan/DOC ratio ranged from 0.02 to 6.03% in the Tibetan Plateau glaciers. This ratios and the correlations between levoglucosan and DOC suggested that biomass burning products contributed only marginally to DOC levels in the TP glaciers. Moreover, the analysis of air mass backward trajectories suggested that levoglucosan and DOC in TP glaciers should be transported from the northwestern TP, internal TP, Central Asia, South and East Asia regions. [ABSTRACT FROM AUTHOR]

Accurately estimating the changes in glacier mass balance and water storage in lakes and reservoirs is critical to studying the water cycle in the inland river basin in northwest China. We used high-resolution satellite images to analyze the changes in water surface area of lakes and reservoirs in Tarim Basin, and used the ICESat-GLAS altimeter data to estimate their water level changes and the glacier mass balance change, over the period 2003-2009. The results showed the average glacier thinning in the entire basin was at a rate of 0.34 ± 0.25 m w.e./year equivalent height of water, which means that the glacier mass balance occurred −6.8 ± 1.2 km water equivalent over the study period. However, the mean water level of nearby lakes decreased by 0.41 ± 0.2 m even with the influx of glacial melt water, indicating that the lake level declines were caused by the withdrawal of lake water for human activities. [ABSTRACT FROM AUTHOR]

Surface temperatures on mountain glaciers are a direct reflection of the heat budget at the glacier surface, which controls the degree of ablation. An understanding of surface temperatures is therefore crucial in simulating the evolution of mountain glaciers and in understanding a glacier's response to climate change. In situ observations indicate that when the exitance angle equals 75°, relative to the nadir angle, the emissivity variation at λ = 12.5 μm in bare ice is as much as 0.24. However, research involving the remote-sensing-retrieval of mountain glacier surface temperatures does not generally take into account the topography (i.e., the exitance angle), and the validation of remote sensing approaches using surficial climate measurements on mountain glaciers is lacking. Here we utilize the ETM + multispectral dataset and resampled water vapor content data from MOD05, combined with the emissivity model of ice and snow, taking the exitance angle into account (Hori et al., 2013), and employ the single-channel algorithm developed by Jiménez-Muñoz et al. (2009), to calculate the remote-sensing-retrieved surface temperature distribution across Qiyi glacier, China, from July 2012 to September 2013. The remote-sensing-retrieved temperatures are validated by the 10-min mean measurements of surface temperature from two automatic weather stations at different altitudes on the glacier, and the results show good agreement at the two sites, with an average bias of 0.96 K and an average root mean square error of 1.5 K. Of note, in summer the influence of thermal emissions from surrounding topography means that the error associated with retrieved temperatures for the glacier tongue and margins might exceed 1.5 K. [ABSTRACT FROM AUTHOR]

Firn line altitude is an indicator of equilibrium line altitude (ELA) and plays an important role in predicting climate behavior and global change. Landsat 30-m resolution data from clear sky at the end of 21 summer melt seasons (July and August, 1990-2011) were used to obtain the firn lines and extract the glacier boundary at Qiyi Glacier in the Qilian Mountains of China. High-resolution Advanced Spaceborne Thermal Emission and Reflection Radiometer 30-m resolution Global Digital Elevation Model (ASTER GDEM) was employed to obtain the firn line altitudes. The firn line altitudes on the Qiyi Glacier increased from 4540 m a.s.l. to 5000 m a.s.l. during the study period, a rise of 460 m. Changes in firn line altitudes had good positive correlations with air temperature data and were in good agreement with measured and simulated ELA changes. The calculated firn zone area decreased (increased) with the increase (decrease) of firn line altitude. Linear regression was used to fit the variational trend, indicating that the firn zone area decreased by -1.19 km2, with a change rate 0.057 km² a-1, during the study period. The firn zone area ratio decreased from 76.2% to 18.1% during the same time. [ABSTRACT FROM AUTHOR]

Albedo, which impacts the surface heat budget of a glacier, is a dominant factor governing snow and ice melt variability. We analyzed variations in albedo on the Dongkemadi Glacier (92°4.94'E, 33°05.88'N) using MODIS daily snow albedo products. The results show that the summer albedo over the accumulation zone and the ablation zone all presented a decreasing trend from 2002 to 2012, and both precipitation and summer temperature were critical factors that impacted variations in the albedo of the glacier. We linked the albedo variation with the measured mass balance of the glacier and found that 2006 and 2010 were high-melt years at the Dongkemadi Glacier, with lower albedo and a more negative mass balance. The area around the equilibrium line altitude (ELA) on the glacier was very sensitive to albedo feedback. The albedo change rate reached approximately -34% when the summer average temperature increased about 1 °C and the summer precipitation declined 94.5 mm. We built a linear model between the annual lowest albedo and the annual mass balance of the Xiao Dongkemadi Glacier, one branch of the Dongkemadi Glacier. This model was used to evaluate the annual mass balance of another branch of the Dongkemadi Glacier, the Da Dongkemadi Glacier. The mass balance of the Da Dongkemadi Glacier was significantly negative in 2006 (-768.7 mm w.e.) and 2010 (-696.9 mm w.e.). This work indicates that the glacier surface annual lowest albedo may be a reasonable proxy for evaluating the inter-annual variability of the mass balance of the glacier. [ABSTRACT FROM AUTHOR]

Dyadobacter tibetensis Y620-1 is the type strain of the species Dyadobacter tibetensis, isolated from ice at a depth of 59 m from a high altitude glacier in China (5670 m above sea level). It is psychrotolerant with growth temperature ranges of 4 to 35°C. Here we describe the features of this organism, together with the draft genome sequence and annotation. The 5,313,963 bp long genome contains 4,828 protein-coding genes and 39 RNA genes. To the best of our knowledge, this is the first Dyadobacter strain that was isolated from glacial ice. This study provides genetic information of this organism to identify the genes linked to its specific mechanisms for adaption to extreme glacial environment. [ABSTRACT FROM AUTHOR]

Understanding changes in glacier mass balance is important because it is indicative of changes in climate and the hydrologic cycle. The latter also has particular influence on people living near glaciers and/or glacier-fed rivers. The Western Himalayas remain one of the regions where recent changes in glacier mass balance are not well-known. The temporal and spatial changes in firn line altitudes are an indicator of equilibrium line altitudes and thus reflect changes in glacier mass balance. Here, we use Himalayan Landsat TM/ETM+ data in July and August (the late summer melt season) to quantify changes in firn line altitudes from 1998 to 2009. We produced reflectance maps through radiometric calibration and atmospheric correction and use topographic correction to remove or reduce terrain or shadow effects. The real ‘surface albedo’ is obtained by narrowband-to-broadband (NTB) albedo conversion from the combined solar radiation. The firn line altitude was then extracted by combining the ‘surface albedo’ with pre-registered digital elevation model. The individual firn line altitude varies by region. The Western Himalayas display the largest range of firn line variability, where the firn line altitudes vary from 4840m a.s.l. to 5770m a.s.l. The individual glacier mean firn line altitude from 1998 to 2009 rose from 5072±77m a.s.l. to 5640±74m a.s.l. in the Western Himalayas. The mean firn line altitude increased from 1998 to 2009. The lowest mean recorded firn line altitude recorded was 5237±166m a.s.l. in 1998, whereas the highest was 5397±135m a.s.l. in 2000. We also observed a difference between the changes in fine line altitudes of northern and southern slopes of the western Himalayans, as the northern slope glaciers display a greater increase in firn line altitudes than the southern slope glaciers. In the southern slope, changes in firn line altitudes correlate with NCDC-NOAA temperature and precipitation data. This sustained increase of firn line altitudes and associated loss of glacier mass imply a persistent loss of stored freshwater in the Western Himalaya. [ABSTRACT FROM AUTHOR]

Twenty snow samples were collected from the Qiyi glacier in Qilian Mountains, the Yuzhufeng glacier in eastern Kunlun Mountains, the Xiaodongkemadi glacier in Tanggula Mountains, and the Gurenhekou glacier in Nyainqêntanglha Range over the Tibetan Plateau. The concentration and distribution features of sixteen priority Polycyclic Aromatic Hydrocarbons (PAHs) were determined by gas chromatography equipped with a mass spectrometry detector (GC-MS). The sources of these PAHs were explored as well. Our results indicated that the average concentrations of PAHs in snow were in the range of 20.45-60.57 ng/L. Maximum PAHs levels were found in the YZF glacier and minimum in the XDKMD glacier. However, no apparent regional distribution pattern of PAHs was found in the glaciers over the Tibetan Plateau. Moreover, the 2-4 ring low molecular weight PAHs predominated in snow samples and the concentrations of phenanthrene was the highest. Integrated factor analysis and isomer pair ratios suggested that PAHs of glaciers over the Tibetan Plateau were derived from low temperature combustion of coal and biomass, and partially from the exhaust gas of locomotives. Air mass back trajectory indicated that organic compounds detected in snowpit of these four glaciers, in the period of time they represented, mainly came from Central Asia and the arid area of Northwest China by westerly wind circulation. [ABSTRACT FROM AUTHOR]

Based on observations of the equilibrium line altitude (ELA) of the Qiyi Glacier in the Qilian Mountain, we established a statistical model between ELA and its major influencing factors, warm season air temperature (air temperature averages for September, July and August) and cold season precipitation (total precipitation in the period January through March). Warm season air temperature was the leading climatic factor influencing ELA variations. The glacier ELA ascends (descends) 172 m when warm season air temperature increases (decreases) by 1°C, and ascends (descends) 62 m when cold season precipitation decreases (increases) by 10%. In the period 1958-2008, the glacier ELA showed a general increasing trend, ascending 230 m and reaching its highest altitude in 2006 at 5131 m a.s.l., close to the glacier summit. If future climate is similar to that in the period 2001-2008, the Qiyi Glacier will not stabilize until it retreats by 2.08 km. [ABSTRACT FROM AUTHOR]

Based on the field observations on Qiyi Glacier during the warm season of 2007, using a digital elevation model (DEM, 15 m resolution), we developed a distributed surface energy- and mass-balance model with an hourly resolution. The model described the effect of topography on shortwave solar radiation, and used a new parameterization for glacier albedo. The model was applied to Qiyi Glacier in the Qilian Mountain, China, for the period 20 : 00 30 June to 12 : 00 10 October 2007, to simulate the firn-line changes, the temporal and spatial variations of mass balance, and the glacial meltwater runoff. The results indicated that the patterns of altitudinal profile of glacier mass-balance were affected mainly by the altitudinal profile of albedo, and the status of the glacier mass balance was influenced directly by the values of albedo. The parameter sensitivity test showed that the model was sensitive to the air temperature lapse rate and precipitation gradient, and also sensitive to the threshold temperature for solid/liquid precipitation. Furthermore, the climate sensitivity test showed that the mass balance was more sensitive to air temperature than precipitation, and the response of mass balance to air temperature change was nonlinear while the response to precipitation change linearly. The negative mass balance trend of the glacier can not be reversed when precipitation increases by 20% and meanwhile air temperature rises by 1°C. [ABSTRACT FROM AUTHOR]

Previous studies found extremely high d-excess in both ice core and glacial melt water in Dasuopu glacier, Xixiabangma, middle of Himalayas. These values are much higher than the global average and those measured in southwest monsoon precipitation. The d-excess variation in over one year at Nyalam station will clarify this phenomenon. Studies show that the high d-excess is related to the seasonal variation of moisture transport to this region. The d-excess values are low during the southwest monsoon active periods, when moisture originated from the humid ocean surface. The d-excess values are higher in non-monsoon months, when moisture is derived from westerly transport. Winter and spring precipitation accounts for a substantial portion of the annual precipitation, resulting in higher d-excess in the yearly precipitation in the middle of Himalayas than other parts of the southern Tibetan Plateau. This finding reveals that the precipitation In the middle of Himalayas is not purely from southwest monsoon, but a large portion from the westerly transport, which is very important for ice core study in this area. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

An introduction is presented in which the editor discusses various reports within the issue on topics including mapping permafrost distribution in China, impact of vegetation on the active layer and permafrost in northwestern China, and changes in summer surface albedo on a glacier.

Black carbon (BC), by the combustion of fossil fuels and biomass, has profound effects on climate change and glacier retreat in industrial eras. In the present study, we report refractory BC (rBC) in an ice core spanning 1850–2014, retrieved from the Hariqin Glacier of the Tanggula Mountains in the central Tibetan Plateau, measured using a single particle soot photometer (SP2). The rBC concentration shows a three-fold increase since the 1950s. The mean rBC concentration was 0.71 ± 0.52 ng mL−1 during 1850s–1940s and 2.11 ± 1.60 ng mL−1 during 1950s–2010s. The substantial increase in rBC since the 1950s is consistent with rBC ice core records from the Tibetan Plateau and Eastern Europe. According to the predominant atmospheric circulation patterns over the glacier and timing of changes in regional emissions, the post-1950 amplification of rBC concentration in the central Tibetan Plateau most likely reflects increases in emissions in Eastern Europe, former USSR, the Middle East, and South Asia. Despite the low-level background rBC concentrations in the ice cores from the Tibetan Plateau, the present study highlights a remarkable increase in anthropogenic BC emissions in recent decades and the consequent influence on glaciers in the Tibetan Plateau. Image 1 • BC deposition (1850–2014) in a Central Tibetan glacier has tripled since the 1950s. • BC emission from Indian monsoon domain mainly contributes to the increase. • BC deposition in Tibetan glaciers is at the lowest level in Northern Hemisphere. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

As a receptor of atmospheric deposition, glaciers are considered an ideal archive in the study of climate change and geochemical cycles. The deposition of methanesulfonic acid (MSA) in the glaciers provides good opportunities to study the biogeochemical cycle of sulfur in the cryosphere. In the present work, snow samples were collected from six High Asia glaciers along a north-to-south transect to determine the spatial distribution of MSA and elucidate its potential sources. The median MSA concentration in the Urumqi Glacier No.1 of Tien Shan was 138.8 ng mL−1, which was distinctly higher than those observed in the Tibetan Plateau (TP) glaciers and polar regions. The levels of MSA in the interior TP glaciers were higher than those observed in the margins of northeastern and southeastern TP. Good correlations between MSA and K+ (r = 0.86, n = 30, α = 0.01), Mg2+ (r = 0.86), and NH 4 + (r = 0.73) were observed in continental glaciers. Principal component analysis indicated that MSA may have terrigenous material inputs. At Yulong Snow Mountain, MSA was correlated with Na+ (r = 0.76, n = 8, α = 0.1), a sea-salt tracer ion, suggesting that MSA may be derived from marine environments. According to dimethyl sulfide (DMS) production and NH 3 emissions in High Asia, we deduced that the high concentrations of MSA in continental glaciers are possibly related to the sources of hypersaline soil environments and animal husbandry in nomadic areas. This work is useful for further studies on regional sulfur cycling and the impacts of human activities on climate change. • The highest level of methanesulfonic acid (MSA) was in Urumqi Glacier No.1 (TS). • MSA was correlated with ions from natural dust and biomass in continental glaciers. • MSA may be derived from marine environments in Yulong Snow Mountain. • The sources of MSA might be hypersaline soil environments and animal husbandry. [ABSTRACT FROM AUTHOR]

Large-scale spatiotemporal variations in land surface temperature (LST) can be derived from satellite observations, particularly for remote high-elevated mountain glaciers where in situ observations are scarce. Our previous LST retrieval study on Qiyi Glacier of China employed a generalized single-channel algorithm and Landsat thermal infrared data, but a clear LST summertime overestimation was observed, particularly on the glacier tongue. Here we extend our previous study by calculating the thermal radiance on Qiyi Glacier that is emitted from the surrounding terrain, and incorporating this thermal radiation term into two generalized single-channel LST algorithms to analyze its influence on the errors of glacier surface temperatures retrieved from remote sensing data. In addition, we also discuss other possible causes of this overestimation in LST retrieval. The key findings of our study are as follows. (1) The surrounding terrain provides nonnegligible summertime thermal radiation to the glacier in the Landsat7 band 6 spectrum range (10.4–12.5 μm) that is 2.6 times larger than the atmospheric down-welling thermal radiance on the glacier tongue. The retrieved LSTs yield an overestimation of up to 0.28 K when the effect of thermal radiation from the surrounding slopes is not considered. (2) The average errors due to the surrounding topographic thermal radiation and the cubic convolution resampling method (CC method) employed in the Landsat7 thermal infrared data processing are only 0.07 K and − 0.01 K, respectively. These are much smaller than the average LST errors (1.79 K), indicating that the neither the surrounding topographic thermal radiation nor the CC method are primary factors in the retrieved LST overestimation. (3) The sole dependence of the generalized single-channel algorithm on the atmospheric water vapor content may therefore be the most likely cause of this overestimation, as high-altitude regions generally possess low atmospheric water contents. Radiative transfer models (e.g., MODTRAN) may provide a solution to this problem by simulating the influence of other atmospheric constituents on LST retrieval under conditions of low water vapor, which can then indicate any additional input variables that should be incorporated into the LST algorithm. • Thermal radiation from surrounding slopes is incorporated into the SC algorithm. • An error of up to 0.28 K when thermal radiation of surrounding slopes is neglected. • Surrounding terrain and CC method are not the major reasons for LST overestimation. [ABSTRACT FROM AUTHOR]

Glaciers have been recognized as biomes, dominated by microbial life. Many novel species have been isolated from glacier ecosystems, and their physiological features are well characterized. However, genomic features of bacteria isolated from the deep ice core are poorly understood. In this study, we performed a comparative genomic analysis to uncover the genomic features of strain Dyadobacter tibetensis Y620-1 isolated from a 59 m depth of the ice core drilled from a Tibetan Plateau glacier. Strain D. tibetensis Y620-1 had the smallest genome among the 12 cultured Dyadobacter strains, relatively low GC content, and was placed at the root position of the phylogenomic tree. The gene family based on a nonmetric multidimensional scaling (NMDS) plot revealed a clear separation of strain D. tibetensis Y620-1 from the reference strains. The genome of the deep ice core isolated strain contained the highest percentage of new genes. The definitive difference is that all genes required for the serine-glyoxylate cycle in one-carbon metabolism were only found in strain D. tibetensis Y620-1, but not in any of the reference strains. The placement of strain D. tibetensis Y620-1 in the root of the phylogenomic tree suggests that these new genes and functions are of ancient origin. All of these genomic features may contribute to the survival of D. tibetensis Y620-1 in the glacier. [ABSTRACT FROM AUTHOR]

Glaciers in the Qilian Mountains are important sources of fresh-water for sustainable development in the Hexi Corridor in the arid northwest China. Over the last few decades, glaciers have generally shrunk across the globe due to climate warming. In order to understand the current state of glaciers in the Qilian Mountains, we compiled a new inventory of glaciers in the region using Landsat Operational Land Imager (OLI) images acquired in 2015, and identified 2748 glaciers that covered an area of 1539.30 ± 49.50 km2 with an ice volume of 81.69 ± 7.40 km3, among which the Shule River basin occupied the largest portion of glaciers (24.8% in number, 32.3% in area, and 35.6% in ice volume). In comparison to previous inventories, glacier area was found to shrink by 396.89 km2 (20.5%) in total, and 109 glaciers with an area of 8.94 km2 disappeared over the period from the 1960s to 2015. This situation was primarily caused by the increase in air temperature, and also related with the size of glacier and some local topographic parameters. In addition, the change of glaciers in the Qilian Mountains showed a distinct spatial pattern, i.e., their shrinking rate was large in the east and small in the west. [ABSTRACT FROM AUTHOR]

[Display omitted] • Three glacial events (1002 ± 190 CE, 1623 ± 149 CE, and 1874 ± 116 CE) were dated in the Drongtso Lungpa Co Valley. • Temperatures were 2.8–0.75 °C lower than today during these glacial events. • Precipitations were up to 70 % higher than today during these glacial events. • Glaciers of southeastern TP have experienced frequent variations over the last ∼ 1000 years, most marked during ∼ 950–1200 CE, ∼1200–1700 CE, and ∼ 1750–1900 CE. • Glacial events in the Drongtso Lungpa Co Valley occurred simultaneously with those in the SE Tibetan Plateau. The temperate glaciers of the southeastern Tibetan Plateau (TP) are highly sensitive to climate change. In this study, 11 exposure ages were obtained from the Drongtso Lungpa Co Valley in the southeastern TP to constrain the timing of glaciations during the late Holocene. The 10Be exposure dating method was employed to reconstruct the palaeoclimate during glaciations, using coupled mass balance and ice flow models. The results uncovered three typical glacial events, constrained at 1002 ± 190 CE, 1623 ± 149 CE, and 1874 ± 116 CE. The modelling results showed that the temperature during the glacial event at 1002 ± 190 CE was 2.8–1.8 °C lower than today, while the precipitation amounts were 20–70 % higher than the present-day level. The modelled temperature during the 1623 ± 149 CE event was 1.5–0.6 °C lower than today, while the modelled precipitation amounts were 10–50 % higher than today. The modelled temperature during the 1874 ± 116 CE event was 1.5–0.75 °C lower than today, while the modelled precipitation amounts were up to 30 % higher than the present-day level. A comparison with the ages of glacial events of the southern TP indicated that the timing of glacial variations in the Drongtso Lungpa Co Valley kept pace with other places in the southern TP. These other places experienced glacial events during ∼ 950–1200 CE, ∼1200–1700 CE, and ∼ 1750–1900 CE. [ABSTRACT FROM AUTHOR]

Delineating glacier change during the Little Ice Age (LIA) is of great importance when attempting to understand regional climatic changes and can also help to improve the understanding of any predictions of future glacial changes. However, such knowledge is still lacking for some critical regions of the Tibetan Plateau (TP). In this study, we mapped 4188 contemporary glaciers and reconstructed 1216 LIA areas of glacial coverage in the Gangdise Mountains to the north of the Himalaya using Google Earth satellite imagery. We estimated their paleoglacial areas and equilibrium line altitudes (ELAs) based on the toe-to-headwall altitude ratio (THAR) method. Results show that most glaciers are small (<1 km 2 ), with slope/hanging glaciers the most common (2844 out of 4188 glaciers), while valley glaciers have the greatest coverage (1009.0 km 2 out of a total area of 1723.7 km 2 ). Contemporary glaciers have retreated significantly since the LIA, with reductions in length of between 5.5% and 94.7% (mean glacier length 696 m; mean reduction in length 41.7%) and reductions in glacier area of between 4.1% and 94.9% (mean glacier area 0.42 km 2 ; mean reduction in area 44.8%). These reductions have occurred under different local climatic and topographic conditions. The contemporary ELA ranges from 5516 to 6337 m asl; the LIA ELA ranged from 5476 to 6329 m asl. Contemporary and LIA ELA values rise from southeast to northwest. As a general rule, the rise in the ELA value decreases from the eastern to the central Gangdise Mountains and then increases westward, with a mean ELA rise of 45 m. Multiple regression models suggest that 46.8% of the glacier area loss can be explained by glacier elevation, area, and slope. However, only 15.5% of the rise in ELA values can be explained by glacial geometric, topographic, or locational parameters. The spatial pattern of modern ELA values in this region appears inversely related to precipitation, which decreases from southeast to northwest, implying that precipitation is one of the key controls of ELAs. This is also consistent with results from elsewhere in High Asia. In contrast to the Gangdise Mountains' eastern and western sectors, glaciers in the central sector have undergone less change, i.e., in terms of reductions in length, area loss, and rises in ELA. Topography can of course also influence glacial change by creating shielding and/or rainshadow effects and by affecting local temperatures. [ABSTRACT FROM AUTHOR]