Glacier Area Research Articles

Transport of nitrate-containing groundwater to coastal areas through buried tunnel valleys, Denmark

Nitrogen impact on the aquatic environment, including coastal areas, is too high in many countries worldwide, particularly in areas with intensive agriculture. Efficient mitigation initiatives demand that important pathways and the fate of nitrate in the hydrological cycle are known. In this study, we focus on groundwater nitrate contamination in two near-shore catchment areas in north-west Denmark. Groundwater in the area is mainly located in buried tunnel valleys, which are subsurface structures eroded by meltwater during Pleistocene glaciations in former glaciated areas. Groundwater samples from the aquifers inside the buried valleys reveal the presence of up to 120 mg/l nitrate down to 10 m below sea level and about 1 km down from the stream outlet towards the coast. We interpret the complex tunnel-valley infill to be responsible for the spatial heterogeneity of the groundwater geochemistry, where sandy geological windows create localised hydraulic pathways and complex redox structures. Groundwater and stream water chemistry in the study area clearly demonstrate the role of groundwater in nitrate transport within the catchment as well as the direct pathway to the coast bypassing the stream and riverine systems. Our results show that the buried tunnel valleys potentially contribute to submarine groundwater discharge and therefore could be responsible for a hitherto unaccounted input of nitrogen to the marine environment.

Composition and constraints of lithium isotopes in cryoconite from various remote glacier areas of the Tibetan Plateau

Lithium isotope is one of the most promising indicators for the study of continental silicate weathering, and lithium concentrations and its isotopic compositions in earth surface can provide a better understanding of the geochemical behavior and isotopic fractionation during weathering and erosion. This work focused on the composition and distribution of Li isotope in cryoconite deposited on various glacier areas in a large range of the Tibetan Plateau and surroundings, as well as its implications for cryoconite dust provenances. Results showed that δ7Li in cryoconite varied within the same order of magnitude (−2.14 ‰-7.74 ‰), which is characterized by geographic distribution of higher δ7Li value of cryoconite in northern glaciers (e.g. Yuzhufeng Glacier), and lower δ7Li value in southern glaciers. In comparison with other global materials, the cryoconite dust shows a lighter δ7Li isotopic composition due to constraints of climatic conditions and land surface weathering intensity. Compared with dust materials in the surrounding Asian dust sources (e.g. large deserts and Gobi), we find that, the primary sources of Li isotope in cryoconite of the northern locations were from both local dust/soils of the TP surface and the surrounding large deserts. Moreover, the products of anthropogenic activities (e.g. coal-burning) may also influence the isotopic composition of the cryoconite dust, and Li isotope may serve as potential tracers of anthropogenic source activities. Therefore, this work provides a complete view of the composition and distribution of Lithium isotopes in cryoconite from various glacier areas of the Tibetan Plateau, and the research significance of its transport processes and source constraints of Li isotopes in cryoconite is proposed.

Comparison of Terrestrial Water Storage Changes in the Tibetan Plateau and Its Surroundings Derived from Gravity Recovery and Climate Experiment (GRACE) Solutions of Different Processing Centers

The GRACE twin satellite gravity mission from 2002 to 2017 has considerably improved investigations on global and regional hydrological changes. However, there are different GRACE solutions and products available which may yield different results for certain regions despite applying the same postprocessing and time span. This is especially the case for the Tibetan Plateau (TP) with its special hydrological conditions represented by localized but strong signals that can overlap or merge with signals inside the plateau, which can falsify the determination of terrestrial water storage (TWS) changes in the TP area. To investigate the effect of GRACE solution selection on inverted TWS changes, we analyze quantitatively the secular and monthly changes for 14 glacier areas and 10 water basins in and around the TP area that have been calculated from 16 different available GRACE solutions. Our analysis provides expectable results. While trend results from different spherical harmonic (SH) GRACE solutions match well, there are significant differences to and between mascon GRACE solutions. This is related to the different processing concepts of mascon solutions and their forced handling in our comparisons. SH solution time series match each other when mass changes are strong with a large amplitude and regular periodicity. However, for regions where small TWS changes are associated with small amplitudes, trends, and/or unstable signal periods, SH solutions can also yield different results. Such behavior is known from a time series analysis. Interestingly though, we find that the COST-G and ITSG SH GRACE solutions are closest to the average of all solutions. Therefore, these solutions appear to be preferable for TWS investigations in regions with highly variable hydrological conditions, such as in the Tibetan Plateau and its surroundings. This also indicates that combined solutions such as COST-G provide a promising pathway for an improved TWS analysis, which should be further elaborated.

Impacts of anthropogenic activities and glacial processes on the distribution of chemical elements in Billefjord, Svalbard, Arctic

The Arctic region is undergoing rapid and extensive transformations due to global climate change. This study investigated the spatial distribution of 31 chemical elements in eight locations in Billefjord, Svalbard, Arctic, with varying degrees of anthropogenic and glacial influences. The west coast of Billefjord has experienced a greater historical anthropogenic impact, while the east coast has larger glaciers and shows less visible evidence of direct human impact. Over 450 topsoil samples collected in the west (abandoned mining town Pyramiden, and glacial valleys of Elsa, Ferdinand, Sven) and east coast of the fjord (glacial valleys of Ebba, Pollock, Ragnar and nearby the Nordenskiöld glacier). These samples were extracted and analyzed by ICP-OES. The results revealed complex distributions of elements among the locations. Nordenskiöld glacier area, along with other locations in the eastern part of the Billefjord, had significantly higher levels of most elements (20 out of 31; As, B, Ca, Cd, Co, Cr, Cu, K, Li, Mg, Mo, Sb, Se, Sn, Sr, Ti, Tl, U, V, Zr). In contrast, Ferdinand Valley and other locations on the western side of the fjord had the lowest mean concentrations of most elements (18 out of 31; B, Ca, Cu, Cd, K, Li, P, Mg, Mo, Sb, Se, Sn, Sr, Ti, Tl, U, V, Zr). These findings highlight the significant influence of glacial processes on the elemental composition of soils within the region. The meltwater flow originating from glaciers in the sampled valleys contributes to the local element load, while the loss of glacier mass is associated with decreased element concentrations within these valleys. These results underscore the complexity of element distribution in the study area and emphasize the necessity for continuous monitoring efforts in this unique and environmentally sensitive region.

Implications of Accuracy of Global Glacier Inventories in Hydrological Modeling: A Case Study of the Western Himalayan Mountain Range

Alpine glaciers are a fundamental component of the cryosphere and are significantly sensitive to climate change. One such region is the Hindukush Karakoram Himalaya (HKH) and Tibetan Plateau (TP) region, which contains more than 40,000 glaciers. There are more than 12 glacier inventories available covering parts of (or the entire) HKH region, but these show significant uncertainties regarding the extent of glaciers. Researchers have used different glacier inventories without assessing their accuracy. This study, therefore, assessed the implications of the accuracy of global glacier inventories in hydrological modeling and future water resource planning. The accuracy assessment of most commonly used two global glacier inventories (Global Land Ice Monitoring from Space-GLIMS v 2.0 and Randolph Glacier Inventory-RGI v 6.0) has been carried out for three sub-basins of the Upper Indus Basin—the Swat, the Chitral, and the Kabul River basins (combined, this is referred to as the Great Kabul River Basin)—with a total basin area of 94,552.86 km2. Glacier outlines have been compared with various Landsat 7 ETM+, Landsat 8, high-resolution Google Earth images, and manually digitized debris-covered glacier outlines during different years. The total glacier area for the Great Kabul River Basin derived from RGI and GLIMS is estimated to be 2120.35 km2 and 1789.94 km2, respectively, which was a difference of 16.9%. Despite being sub-basins of the Great Kabul River Basin, the Swat, and the Chitral River basins were different by 54.74% and 19.71%, respectively, between the two inventories, with a greater glacierized area provided by RGI, whereas the Kabul River basin was different by 54.72%, with greater glacierized area provided by GLIMS. The results and analysis show that GLIMS underestimates glacier outlines in the Swat and the Chitral basins and overestimates glacier extents in the Kabul River basin. The underestimation is mainly due to the non-representation of debris-covered glaciers. The overestimation in GLIMS data is due to the digitization of seasonal snow as part of the glaciers. The use of underestimated GLIMS outlines may result in 5–10% underestimation of glacier-melt contribution to flows in the Swat River basin, while an underestimation of 7% to 15% is expected in the Chitral River Basin, all compared to RGI v 6.0 outlines. The overestimation of glacier-melt contribution to flows in the Kabul River basin is insignificant (1% to 2%) using GLIMS data. In summary, the use of the GLIMS inventory will lead to underestimated flows and show that the Great Kabul River Basin (particularly the Chitral River Basin) is less sensitive to climate change effects. Thus, the current study recommends the use of RGI v 6.0 (best glacier inventory) to revisit the existing biased hydro-climate studies and to improve future hydro-climate studies with the concomitant rectification of the MODIS snow coverage data. The use of the best glacier inventory will provide the best estimates of flow sensitivity to climate change and will result in well-informed decision-making, precise and accurate policies, and sustainable water resource management in the study area. The methodology adopted in the current study may also be used in nearby areas with similar hydro-climate conditions, as well as for the most recently released RGI v 7.0 data.

Divergent glacier area and elevation changes across the Tibetan Plateau in the early 21st century

With accelerated warming, mountain glaciers in most parts of the world have been in a state of continuous retreat in recent decades. Assessing glacier change and analyzing its influencing factors are essential for developing climate change mitigation and adaptation measures for a given region. This study provides a spatially explicit assessment and quantification of glacier changes in the early 21st century on the Tibetan Plateau (TP) at individual glacier and basin scales. We established a one-to-one correspondence between the Second Chinese Glacier Inventory (CGI-2, collected from 2004 to 2011) and a dataset of glacier inventory in Western China during 2017–2018 (CGI-2018). The majority of TP’s glaciers decreased in size with a mean area retreat rate during the investigated period of 4.1%/decade. In addition, a mean change of the median elevation of the glaciers of 6.7 m/decade was detected. Approximately 2.5% of the total number of glaciers mapped in CGI-2 disappeared, while 681 of them divided to 1758 glaciers as they retreated. The observed variations follow local trends and have different regional characteristics. Generally, the glaciers with the lowest retreat rates are found in the Karakorum and Kunlun Mountains, while those with high retreat rates are concentrated along the Gangdis and Tangula ranges. The observed changes in glaciers are mainly attributed to a significant increase in temperature. Other factors including glacier size, debris cover, orientation and mean elevation also contribute to the heterogeneity of glacier variability. This study provides for the first time a detailed spatially explicit analysis of the glacial changes on the TP in the early 21st century, substantially improving the understanding of glacier response patterns and supporting more sustainable utilization of regional water resources in the TP in the context of climate warming in the 21st century.

Cryosphere changes and their impacts on regional water resources in the Chinese Altai Mountains from 2000 to 2021

The cryosphere has an important impact on regional water resources and ecosystems in the Chinese Altai Mountains and its piedmont zone. Using the latest remote-sensing datasets of cryosphere changes and combining with in-situ observation data from glacier monitoring stations and snow cover surveys, the main cryosphere elements including glaciers, snow cover, and permafrost are investigated with emphasis on their changes since 2000 and the current situation. Their water resource effects are also discussed. The results indicate that although the glaciers in the region have experienced continuous and intensive melting, mass loss has slowed because both glacier area shrinkage and thickness reduction were larger during 2000–2010 than during 2010–2021. Snow-cover water equivalent (w.e.) has increased due to obvious increases in snow depth, although snow-cover area has decreased slightly. Permafrost has been degrading. Overall, cryosphere contributions to the regional water resource are approximately 40.9% since 2000, among which snow-cover melting is the largest, contributing 37.1% to water resources in the Irtysh River Basin and significantly more in the mountainous sub-basins with increased snowfall. Glacier melting contributes 2.9%∼3.4%, lower than earlier estimations of 3.4%∼3.6% for the late 20th century. Permafrost thaw caused by active layer thickening contributes approximately 0.59%. Meteorological data shows a warming and wetting trend, but summer temperature has a much lower increase rate and a slowing increase trend after 2013. Moreover, snowfall frequency has increased. In the future, glacier water resource contribution will continue to decrease, but the water resource effects of snow-cover melting and permafrost degradation would increase.

Debris cover effect on the evolution of Northern Caucasus glaciers in the 21st century

More than 13% of the area of the Caucasus glaciers is covered by debris affecting glacier mass balance. Using the Caucasus as example, we introduce a new model configuration that incorporates a physically-based subroutine for the evolution of supraglacial debris into the Global Glacier Evolution Model (GloGEMflow), enabling its application at a regional level. Temporal evolution of debris cover is coupled to glacier dynamics allowing the thickest debris to accumulate in the areas with low velocity. The future evolution of glaciers in the Northern Caucasus is assessed for five Shared Socioeconomic Pathways (SSP) and significance of explicitly incorporating debris-cover formulation in regional glacier modeling is evaluated. Under the more aggressive scenarios, glaciers are projected to disappear almost entirely except on Mount Elbrus, which reaches 5,642 m above sea level, by 2,100. Under the SSP1-1.9 scenario, glacier ice volume stabilizes by 2040. This finding stresses the importance of meeting the Paris Climate Agreement goals and limiting climatic warming to 1.5 °C. We compare evolution of glaciers in the Kuban (more humid western Caucasus) and Terek (drier central and eastern Caucasus) basins. In the Kuban basin, ice loss is projected to proceed at nearly double the rate of that in the Terek basin during the first half of the 21st century. While explicit inclusion of debris cover in modeling leads to a less pronounced projected ice loss, the maximum differences in glacier length, area, and volume occur before 2,100, especially for large valley glaciers diminishing towards the end of the century. These projections show that on average, fraction of debris-covered ice will increase while debris cover will become thinner towards the end of the 21st particularly under the more aggressive scenarios. Overall, the explicit consideration of debris cover has a minor effect on the projected regional glacier mass loss but it improves the representation of changes in glacier geometry locally.

Automatic Extraction of the Calving Front of Pine Island Glacier Based on Neural Network

Calving front location plays a crucial role in studying ice–ocean interaction, mapping glacier area change, and constraining ice dynamic models. However, relying solely on visual interpretation to extract annual changes in the calving front of ice shelves is a time-consuming process. In this study, a comparative analysis was conducted on the segmentation obtained from fully convolutional networks (FCN), U-Net, and U2-Net models, revealing that U2-Net exhibited the most effective classification. Notably, U2-Net outperformed the other two models by more than 30 percent in terms of the F1 parameter. Therefore, this paper introduces an automated approach that utilizes the U2-Net model to extract the calving front of ice shelves based on a Landsat image, achieving an extraction accuracy of 58 m. To assess the model’s performance on additional ice shelves in the polar region, the calving front of the Totten and Filchner ice shelves were also extracted for the past decade. The findings demonstrated that the ice velocity of the Filchner ice shelf exceeded that of the Totten ice shelf. Between February 2014 and March 2015, the majority of the calving fronts along the Filchner Ice Shelf showed an advancing trend, with the fastest-moving front measuring 3532 ± 58 m/yr.

Monitoring glacier calving using underwater sound

Abstract. Climate shifts are particularly conspicuous in glaciated areas. Satellite and terrestrial observations show significant increases in the melting and breakup of tidewater glaciers and their influence on sea level rise and ocean mixing. Increasing melt rates are creating an urgency to better understand the link between atmospheric and oceanic conditions and glacier frontal ablation through iceberg calving and melting. Elucidating this link requires a combination of short- and long-timescale measurements of terminus activity. Recent work has demonstrated the potential of using underwater sound to quantify the time and scale of calving events to yield integrated estimates of ice mass loss (Glowacki and Deane, 2020). Here, we present estimates of subaerial calving flux using underwater sound recorded at Hansbreen, Svalbard, in September 2013 combined with an algorithm for the automatic detection of calving events. The method is compared with ice calving volumes estimated from geodetic measurements of the movement of the glacier terminus and an analysis of satellite images. The total volume of above-water calving during the 26 d of acoustical observation is estimated to be 1.7±0.7×107 m3, whereas the subaerial calving flux estimated by traditional methods is 7±2×106 m3. The results suggest that passive cryoacoustics is a viable technique for long-term monitoring of mass loss from marine-terminating glaciers.

Supervised Methods for Modeling Spatiotemporal Glacier Variations by Quantification of the Area and Terminus of Mountain Glaciers Using Remote Sensing

Glaciers are important indictors of climate change as changes in glaciers physical features such as their area is in response to measurable evidence of fluctuating climate factors such as temperature, precipitation, and CO2. Although a general retreat of mountain glacier systems has been identified in relation to centennial trends toward warmer temperatures, there is the potential to extract a great deal more information regarding regional variations in climate from the mapping of the time history of the terminus position or surface area of the glaciers. The remote nature of glaciers renders direct measurement impractical on anything other than a local scale. Considering the sheer number of mountain glaciers around the globe, ground measurements of terminus position are only available for a small percentage of glaciers and ground measurements of glacier area are rare. In this project, changes in the terminal point and area of Franz Josef and Gorner glaciers were quantified in response to climate factors using satellite imagery taken by Landsat at regular intervals. Two supervised learning methods including a parametric method (multiple regression) and a nonparametric method (generalized additive model) were implemented to identify climate factors that impact glacier changes. Local temperature, CO2, and precipitation were identified as significant factors for predicting changes in both Franz Josef and Gorner glaciers. Spatiotemporal quantification of glacier change is an essential task to model glacier variations in response to global and local climate factors. This work provided valuable insights on quantification of surface area of glaciers using satellite imagery with potential implementation of a generic approach.

Micrometeorological Analysis and Glacier Ablation Simulation in East Kunlun

Worldwide, there are great challenges for meteorological monitoring and glacier ablation monitoring in high-altitude mountain areas. It is often difficult to capture fine-scale climate and glacial changes in high-altitude mountainous areas due to the harsh natural environment and the extreme lack of observational sites. Based on high-altitude meteorological stations erected on the eastern shore of Aqikkule Lake (AQK) and at the terminus of Shenshechuan Glacier (SSG), as well as on mass balance data from SSG, the characteristics and correlation of temperature, solar radiation, relative humidity, precipitation, wind speed and direction of the two regions, and the mass balance in the ablation area of SSG from 30 May 2022 to 18 May 2023 were analyzed, and the average melting depth of SSG was simulated. The results indicate the following: (1) The average annual temperature of AQK and the terminus of SSG is −3.7 °C and −7.7 °C, respectively, and the vertical lapse rate of temperature in the summer half of the year is greater than that in the winter half of the year. Precipitation timing has a great influence on daily temperature differences. (2) Precipitation in both places is concentrated in summer; the glaciers in this area are of the summer recharge type, and precipitation has a significant reducing effect on the solar incident radiation and increases the relative humidity in this region. (3) AQK and SSG both have local circulation development, in the area of AQK all year round due to the lake effect, while the terminus of SSG only has the development of valley winds in the summer, being controlled in the winter by the westerly wind belt. (4) The average mass balance value of the ablation area of SSG was −1786 mm as measured by the range poles method. The average annual ablation depth of SSG simulated by using the empirical formula was 587–597 mm, which is not large compared with other glacier areas in the Tibetan Plateau, and it has the characteristics of typical continental-type glaciers.

Spatial distribution and controls of snowmelt runoff in a sublimation-dominated environment in the semiarid Andes of Chile

Abstract. Sublimation is the main ablation component of snow in the upper areas of the semiarid Andes (∼ 26 to ∼ 32∘ S and ∼ 69 to ∼ 71∘ W). This region has elevations up to 6000 m, is characterized by scarce precipitation, high solar radiation receipt, and low air humidity, and has been affected by a severe drought since 2010. In this study, we suggest that most of the snowmelt runoff originates from specific areas with topographic and meteorological features that allow large snow accumulation and limited mass removal. To test this hypothesis, we quantify the spatial distribution of snowmelt runoff and sublimation in a catchment of the semiarid Andes using a process-based snow model that is forced with field data. Model simulations over a 2-year period reproduce point-scale records of snow depth (SD) and snow water equivalent (SWE) and are also in good agreement with an independent SWE reconstruction product as well as satellite snow cover area and indices of winter snow absence and summer snow persistence. We estimate that 50 % of snowmelt runoff is produced by 21 %–29 % of the catchment area, which we define as “snowmelt hotspots”. Snowmelt hotspots are located at mid-to-lower elevations of the catchment on wind-sheltered, low-angle slopes. Our findings show that sublimation is not only the main ablation component: it also plays an important role shaping the spatial variability in total annual snowmelt. Snowmelt hotspots might be connected with other hydrological features of arid and semiarid mountain regions, such as areas of groundwater recharge, rock glaciers, and mountain peatlands. We recommend more detailed snow and hydrological monitoring of these sites, especially in the current and projected scenarios of scarce precipitation.

Glacier changes since the Last Glacial Maximum on two slopes of Mt Noijin Kang‐Sang, Southern Tibetan Plateau

ABSTRACTIn recent decades, most glaciers have been melting, thinning and retreating globally in response to continuously increasing temperatures. We simulated ice thicknesses and volumes on the east and west slopes of Mt Noijin Kang‐Sang, Southern Tibetan Plateau (TP), since the Last Glacial Maximum (LGM), using a glacial flowline model. The simulated average ice thicknesses during the LGM and Lateglacial (LG), Early Holocene, Neoglacial and the Little Ice Age (LIA) periods were 1.4, 1.3, 1.1 and 1.2 times greater than those of modern glaciers in the Gangbu (eastern slope) and western valleys, respectively. In the Gangbu Valley, areas from the LGM‐LG to LIA periods were 1.7, 1.5, 1.4 and 1.2 times greater than the modern glacier area, and the volume expansion indexes were 2.2, 1.9, 1.5 and 1.4. In the western valleys, the area expansion indexes were 2.2, 1.9, 1.5 and 1.4 times greater than the modern glacier, and volumes were 5.4, 4.4, 3.4 and 2.9 times greater, respectively. Glaciers in the western valleys retreated more extensively than those on the eastern slopes of Gangbu Valley after the LGM‐LG, probably due to the smaller glaciers on the leeward western slopes with lower elevations. Equilibrium‐line altitudes dropped ~425 m in Gangbu Valley and ~385 m in the western valleys during the LGM‐LG, corresponding to a 3.8–4.1 °C temperature decrease in this region, which was consistent with cooling on the central Tibetan Plateau, but lower than cooling at the south‐eastern edge of the Tibetan Plateau.

Prediction of glaciated area fraction over the Sikkim Himalayan Region, India: a comparative study using response surface method, random forest, and artificial neural network.

Glacier area fraction at high altitude mountains is a serious worry in today's time triggered by climate change. The current information on this natural resource is very important for the survival of humanity as it affects the water, food, and energy security of people dependent on it. Due to its problematic accessibility and tough environmental condition, ground monitoring is quite challenging. This study investigates the impact of environmental parameters and pollutants on glacier area fraction over the Eastern Himalaya region and its prediction through random forest (RF), multilayer perceptron (MLP), radial basis function analysis (RBFN), and response surface methodology (RSM) models. The data are obtained from the Goddard Earth Sciences Data and Information Services Center (GES DISC), NASA's data archive portal ( https://giovanni.gsfc.nasa.gov ). The collinearity of independent variables reveals that all selected input parameters are highly correlated with R2 value > 0.9. The RSM and RF model provided valuable insight of the predictor's significance in addition to their capability to predict the response. The model performance was evaluated in terms of R2 value and the error matrices. The model's R2 value was found to be 0.843, 0.839, 0.838, and 0.743 for MLP, RBFN, RF, and RSM respectively. Although, the neural network model R2 values are the highest, but the most reliable and suitable model is RF as the error matrices for this model are much lower than others. This study encourages the investigation of the hybridization of these models for more accurate prediction.

Surging Glaciers in High Mountain Asia between 1986 and 2021

High Mountain Asia (HMA) is one of the concentrated areas of surging glaciers in the world. The dynamic movement of surging glaciers not only reshapes the periglacial landscape but also has the potential to directly or indirectly trigger catastrophic events. Therefore, it is crucial to understand the distribution patterns, periodicities, and occurrence mechanisms of surging glaciers. Based on Landsat TM/ETM+/OLI remote sensing images from 1986 to 2021, a total of 244 surging glaciers were identified in HMA in this study, covering an area of 11,724 km2 and accounting for 12.01% of the total area of glaciers in this region. There are 185 surging glaciers identified within the Karakoram Range and Pamirs, which constitute the primary mountainous regions in HMA. From 1986 to 2021, these surging glaciers advanced at least 2802 times and exhibited different temporal and spatial patterns. A total of 36 glaciers in HMA experienced 2 or more surges during this period, with the highest number observed in the Pamirs (19), followed by the Karakorum (13), with the other regions having fewer occurrences. Obvious differences exist in the surge phase and the quiescent phase of glaciers in different regions of HMA. The surge phase of surging glaciers in the Karakoram Range and Pamirs is generally short, mostly in the range of 2~6 years. The quiescent phase lasts for 5~19 years and the overall surge cycle ranges from 9 to 24 years. The complex nature of glacier surges in HMA suggests that multiple mechanisms may be at play, necessitating further research.

Impact of far-field glacially-induced stresses on fault stability in the eastern Paris Basin

Stable Continental Regions (SCR) feature very low seismicity levels, but historical seismicity reveals that these areas can experience large damaging earthquakes of magnitude up to Mw ≈ 7, on faults that were previously undocumented or thought to be stable. Glacial Isostatic Adjustment (GIA) is commonly considered as a factor of SCR deformation and stress perturbations that can trigger fault failure in formerly glaciated areas. Studies suggest that fault reactivation due to GIA may not only be restricted to glaciated areas, but could also occur in their periphery, several hundreds of kilometers away from the former glaciations. In this study, we model GIA associated to the Alpine, Massif Central, Celtic and Scandinavian glaciations from the last glacial cycle, in order to estimate the impact of this process on fault stability in the eastern Paris Basin, an area featuring very low strain rates and seismic activity, but of particular interest due to the proposed construction of a radioactive waste disposal facility. Our computations predict stress perturbations far (hundreds of kilometers away) from the ice loads, of similar amplitude to those in areas where fault reactivation has been observed (1–8 MPa at LGM, 0.8–2 MPa at present-day in our study area). We also show that the interaction between several GIA systems can impact the timing of maximum fault destabilization in a different way than common models considering only a single ice load. The computed far-field glacially-induced stress perturbations remain small compared to ambient crustal stresses, but they are large enough to potentially destabilize faults close to Andersonian geometries, with known faults in the eastern Paris Basin meeting these criteria, mainly in the Marne and Saint-Martin-de-Bossenay fault systems. Yet, there is no known evidence of Quaternary fault activity in the Paris Basin, and discrepancies between the localization of potentially unstable faults due to GIA and seismically active areas in northern France suggests that other processes must be at play.

Closing Greenland's Mass Balance: Frontal Ablation of Every Greenlandic Glacier From 2000 to 2020

AbstractIn Greenland, 87% of the glacierized area terminates in the ocean, but mass lost at the ice‐ocean interface, or frontal ablation, has not yet been fully quantified. Using measurements and models we calculate frontal ablation of Greenland's 213 outlet and 537 peripheral glaciers and find a total frontal ablation of 481.8 ± 24.0 for 2000–2010 and 510.2 ± 18.6 Gt a−1 for 2010–2020. Ice discharge accounted for ∼90% of frontal ablation during both periods, while mass loss due to terminus retreat comprised the remainder. Only 16 glaciers were responsible for the majority (>50%) of frontal ablation from 2010 to 2020. These estimates, along with the climatic‐basal balance, allow for a more complete accounting of Greenland Ice Sheet and peripheral glacier mass balance. In total, Greenland accounted for ∼90% of Northern Hemisphere frontal ablation for 2000–2010 and 2010–2020.

Historical and projected evolutions of glaciers in response to climate change in High Mountain Asia

Glacier changes are regarded as the conspicuous icon of climate change in High Mountain Asia (HMA) alpine environments. Multi-temporal glacier coverage is essential for mass balance estimations and understanding glacial changes in response to climate variability. However, consistent multi-temporal glacier area datasets across the HMA are limited due to challenges posed by seasonal snow and cloud cover in remote sensing satellite images. In this study, a new method is proposed to estimate glacier and nonseasonal snow (GNS) areas by combining the normalized difference snow index (NDSI) and bright temperature bands from Landsat remote sensing images in 1990s, 2000s, 2010s and 2020s. Meanwhile, the integrated ice dynamic Open Global Glacier Model (OGGM) is applied to simulate nine typical glaciers and project their changes using 15 general circulation models (GCMs) under four shared socioeconomic pathways (SSPs). The results showed that total GNS areas in HMA were estimated at 71.26 × 103 km2 in the 1990s, and decreased to 53.17 × 103 km2 in the 2020s. Massive proportions of GNS were mainly distributed in the southeast and northwest region. Tipping points of GNS indicating a transition from mass balance gain to loss during four periods were detected in the Upper Tarim basin, Inner Tibetan Plateau, Upper Yangtze river basin, Upper Brahmaputra river basin, and Upper Salween river basin. The reduction of total GNS area was attributed to significantly rising temperature and decreasing solid precipitation. The projections indicate that glacier shrinkages will be the dominant trend in the future, and most typical glaciers will experience a mass loss of 60% under SSP245 (SSP2-RCP4.5) and over 90% under SSP585 (SSP5-RCP8.5) by the end of the century. The response of glaciers to climate change is also found to be influenced by factors including local topography and elevation ranges. The findings of this study provides a better understanding of glacial dynamics in response to climate change and highlights the importance of making water resources management strategies to mitigate the influence of glacier retreat.

Unlocking archival maps of the Hornsund fjord area for monitoring glaciers of the Sørkapp Land peninsula, Svalbard

Abstract. Archival maps are an important source of information on the state of glaciers in polar zones and are very often basic research data for analysing changes in glacier mass, extent, and geometry. However, basing a quantitative analysis on archival maps requires that they be standardised and precisely matched against modern-day cartographic materials. This can be achieved effectively using techniques and tools from the field of geographic information systems (i.e. GIS). The objective of this research was to accurately register archival topographic maps of the area surrounding the Hornsund fjord (southern Spitsbergen) published by the Polish Academy of Sciences and to evaluate their potential for use in studying changes in the geometry of glaciers in the north-western part of the Sørkapp Land peninsula in the following periods: 1961–1990, 1990–2010, and 1961–2010. The area occupied by the investigated glaciers in the north-western Sørkapp Land decreased in the years 1961–2010 by 45.6 km2, i.e. by slightly over 16 %. The rate of glacier area change varied over time and amounted to 0.85 km2 yr−1 in the period 1961–1990 and sped up to 1.05 km2 yr−1 after 1990. This process was accompanied by glacier surface lowering by about 90–100 m for the largest land-terminating glaciers on the peninsula and by up to more than 120 m for tidewater glaciers (above the line marking their 1984 extents). The dataset is now available from the Zenodo web portal: https://doi.org/10.5281/zenodo.4573129 (Dudek and Pętlicki, 2021).