Search

• UAVs used for methane flux quantification at a landfill site. • Methane emissions measured: 150.7 kg h−1, with a 1σ uncertainty range of 83.0 kg h−1 to 209.5 kg h−1. • Refined UAV-based mass balance flux method. • Study addresses UAV onboard wind measurements. Accurate quantification of methane emissions from landfills is crucial for improving greenhouse gas inventories and mitigating climate change impacts. Existing methodologies, such as theoretical gas production models and labour-intensive measurement approaches, present limitations including large uncertainties and high operational costs. This study adds to a growing body of research and applications which aim to bridge this gap. To this end, we present a case study using Unmanned Aerial Vehicles (UAVs) equipped with methane and wind instrumentation for a survey of a landfill site in Bury, Manchester, UK, in summer 2022, in order to evaluate and reflect the challenges of the UAV-based mass balance method for quantification of methane emissions from a large heterogeneous source such as landfill. This study offers guidance on defining an appropriate methane background concentration, geospatial interpolation of sampled date, survey sampling strategy, and more importantly, addresses the challenges surrounding UAV wind measurements and spatial characterisation of emission plumes. For the period of the case study, we quantified methane flux for the landfill site to be 150.7 kg h−1 with a 1 standard deviation uncertainty range of 83.0 kg h−1 to 209.5 kg h−1. [ABSTRACT FROM AUTHOR]

Glaciers in the Alps are losing mass at unprecedented and accelerating rates. Monitoring of glacier mass change as well as relevant atmospheric parameters plays an important role in improving understanding of local and downstream impacts. We present subseasonal, seasonal, and annual glaciological mass balance data and meteorological observations from Mullwitzkees and Venedigerkees, two glacier monitoring sites in the Hohe Tauern range of the Austrian Alps. Ablation stake networks were established on Mullwitzkees in 2006/07 and on Venedigerkees in 2011/12. Monitoring is ongoing. In addition to stake readings at subseasonal intervals, accumulation measurements (snow pits and probing) are carried out seasonally. The glaciological data set consists of subseasonal floating date measurements as well as fixed date seasonal and annual values. Fixed date glacier wide mass balance was derived from annual point mass balance values. Automatic weather stations measuring standard meteorological parameters were installed near Mullwitzkees and Venedigerkees in 2020 and 2019, respectively. Meteorological data is provided in 10 minute intervals. Uncertainties for individual point mass balance measurements were computed following the approach of the Swiss glacier monitoring service (GLAMOS), taking into account estimated density and reading errors. The subseasonal mass balance records highlight shorter term variability in mass loss and the linkage with meteorological conditions. The most negative annual point mass balance recorded in the period of record was -5.8 ± 0.66 m w.e. at an elevation of 2536 on Venedigerkees. 2022 stands out as the most negative mass balance year to date in both time series, particularly at higher elevations. The cumulative specific mass balance (glacier wide) over the period of record was -14.68 m w.e. at Mullwitzkees and -8.79 m w.e. at Venedigerkees. Data is available in PANGAEA publication series and the associated datasets. The main publication series are updated annually. The Mullwitzkees mass balance datasets can be found at: doi: 10.1594/PANGAEA.965660 and doi: 10.1594/PANGAEA.965719. The Venedigerkees data can be found at doi: 10.1594/PANGAEA.965648 and doi: 10.1594/PANGAEA.965729. [ABSTRACT FROM AUTHOR]

In recent decades, glaciers in the southeastern Tibetan Plateau (SETP) have been rapidly melting and showing a large scale of glacier mass loss. Due to the lack of large-scale, high-resolution, and high-precision observations, knowledge on the spatial distribution of the glacier mass balance and the response to climate change is limited in this region. We propose a TanDEM-X bi-static InSAR (Interferometric Synthetic Aperture Radar) algorithm with a non-local mean filter method and difference strategy, to improve the precision of glacier surface elevation change detection. Moreover, we improved the glacier mass balance estimation algorithm with a correction method for multi-source system errors and an uncertainty evaluation method based on error propagation theory to reduce the uncertainty of estimations. We used 13 pairs of TanDEM-X bi-static InSAR images to obtain the glacier mass balance data for the entire SETP. The total area of glaciers monitored was 5821 km2 and the total number of glaciers monitored was 2321; the glacier surface elevation change rate was −0.505 ± 0.005 m/yr, and the glacier mass balance estimation was −454.5 ± 13.1 mm w.eq. during 2000–2014. Additionally, we analyzed the spatial distribution of the glacier mass balance within the SETP using the sub-watershed analysis method. The results showed that the mass loss rate had a decreasing trend from the southeast to the northwest. Furthermore, the temperature change and the glacier mass loss rate showed a positive correlation from the southeast to the northwest in this region. This study greatly advances our understanding of the regularities of glacier dynamics in this region, and can provide scientific support for major national goals such as the rational utilization of surrounding water resources and construction of important transportation projects. [ABSTRACT FROM AUTHOR]

Recent satellite measurements of glacier mass balances show mountain glaciers all over the world had generally negative mass balances in the first decades of the 21st century. Mean summer temperatures all over the world rose from the 1961–1990 period to the 1991–2020 period, implying increasingly negative mass balances. We studied archived annual balances for 38 northern hemisphere glaciers to assess changes within the 1961–2020 period. We used a modified double-mass curve to visualize mass balance changes occurring around 1990. Mean balances in 1961–1990 were already small negative for many of the studied glaciers and became even more negative in 1991–2020 for glaciers in the Alps, at high latitudes and in western North America. The largest mass balance changes were for some glaciers in the Alps. We are unable to explain the lack of change in mean balance for one glacier in High Mountain Asia. We found complex changes for eight glaciers in Scandinavia, even including one glacier with a positive balance. We explain these changes by visualizing the deviations in winter and summer balances from their respective 1961–1990 mean values. High winter balances in the 1990s for Scandinavia partly obscured the emerging trend of increasingly negative summer balances, which we expect to continue in the future. [ABSTRACT FROM AUTHOR]

Recent studies have raised concerns regarding the reconstruction of glacier mass balance using tree-ring data. The method relies on a stable relationship between both variables and summer (June–August) or melt season (May–September) temperature. However, with recent anthropogenic climate change the stability of this relationship is challenged by both, a divergence between tree-ring and temperature, as well as mass balance and temperature data. Establishing to what extent this divergence influences the reconstruction of mass balance using tree-ring data is important to assess the future viability and applicability of the method. In this paper we analyze the relationship between the Torneträsk tree-ring and Storglaciären mass balance records, their response to climate change, and investigate changes in the relationship. We provide evidence for a sensitivity loss in the Torneträsk record and quantify its impact on the reconstruction of summer mass balance of Storglaciären. We find that by including years post 1980, the amplitude of reconstructed variability is reduced, but it remains possible to explain the variance of the record significantly. Our results suggest that for glaciers without an extensive mass balance record the applicability of the method is challenged. [ABSTRACT FROM AUTHOR]

Drifting and blowing snow are important features in polar and high mountain regions. They control the surface mass balance in windy conditions and influence sublimation of snow and ice surfaces. Despite their importance, model representations in weather and climate assessments have high uncertainties because the associated physical processes are complex and highly variable in space and time. This contribution investigates the saltation system, which is the lower boundary condition for drifting and blowing snow models. Using a combination of (previous) measurements and new physics-based modeling with large-eddy simulation (LES), we show that the prevailing parameterizations that describe the saltation system in atmospheric models are based on contradictory assumptions: while some scaling laws are typical of a saltation system dominated by aerodynamic entrainment, others represent a saltation system controlled by splash. We show that both regimes can exist, depending on the friction velocity. Contrary to sand saltation, aerodynamic entrainment of surface particles is not negligible. It is important at low wind speeds, leading to a saltation height and near-surface particle velocity which increase with the friction velocity. In a splash-dominated saltation regime at higher friction velocities, the saltation height and near-surface particle velocity become invariant with the friction velocity and closer to what is observed with sand. These findings are accompanied by a detailed description of the theoretical, experimental and numerical arguments behind snow saltation parameterizations. This work offers a comprehensive understanding of the snow saltation system and its scaling laws, useful for both modelers and experimentalists. [ABSTRACT FROM AUTHOR]

Kangerlussuaq Glacier (KG) is a major contributor to central‐eastern Greenland mass loss, but existing estimates of its mass balance over the last century are inconsistent, and specific drivers of change remain poorly understood. We present a novel approach that combines numerical modeling and a 1933–2021 climate forcing to reconstruct its mass balance over the past century. The model's final state aligns remarkably well with present‐day observations. The model reveals a total ice mass loss of 285 Gt over the past century, equivalent to 0.68 mm global sea level rise, 51% of which occurred since 2003 alone. Dynamic thinning from ice front retreat is responsible for 88% of mass change since 1933, with short‐term ice front variations having minimal impact on centennial mass loss. Compared to earlier studies, our findings suggest that KG lost 59% (or 301 Gt) less mass over the century than previously thought. Plain Language Summary: Kangerlussuaq Glacier (KG) is a primary source of mass loss from the central‐eastern sector of the Greenland Ice Sheet. Despite various efforts, there is still disagreement on how much mass KG has lost over the past century and the reason for the mass loss. We used computer models and climate data to determine that the glacier has lost 285 billion metric tons of ice over the last 100 years. That is enough to raise global sea levels by about two‐thirds of a millimeter. Over half of this ice loss has happened in the last 20 years, and the rates of mass loss in the past two decades were significantly higher than for most of the last century. While the ice front of the glacier can retreat without losing mass, if the snowfall is significant enough, we found that most ice loss was related to the glacier retreating. This new estimate of ice loss is 59% lower than previous studies suggested, showing that the glacier's decline is less severe than we thought. Key Points: Dynamic thinning due to the retreat of the ice front explains 88% of the mass loss since 1933Short‐term variability in the ice front has a lesser impact on the mass balance on centennial time scalesKangerlussuaq Glacier lost 59% (or 301 Gt) less mass over the century than previously thought [ABSTRACT FROM AUTHOR]

The mass loss of glaciers outside the polar ice sheets has been accelerating during the past several decades and has been contributing to global sea-level rise. However, many of the mechanisms of this mass loss process are not well understood, especially the calving dynamics of marine-terminating glaciers, in part due to a lack of high-resolution calving front observations. Svalbard is an ideal site to study the climate sensitivity of glaciers as it is a region that has been undergoing amplified climate variability in both space and time compared to the global mean. Here we present a new high-resolution calving front dataset of 149 marine-terminating glaciers in Svalbard, comprising 124 919 glacier calving front positions during the period 1985–2023 (10.5281/zenodo.10407266, Li et al., 2023). This dataset was generated using a novel automated deep-learning framework and multiple optical and SAR satellite images from Landsat, Terra-ASTER, Sentinel-2, and Sentinel-1 satellite missions. The overall calving front mapping uncertainty across Svalbard is 31 m. The newly derived calving front dataset agrees well with recent decadal calving front observations between 2000 and 2020 (Kochtitzky and Copland, 2022) and an annual calving front dataset between 2008 and 2022 (Moholdt et al., 2022). The calving fronts between our product and the latter deviate by 32 ± 65 m on average. The R2 of the glacier calving front change rates between these two products is 0.98, indicating an excellent match. Using this new calving front dataset, we identified widespread calving front retreats during the past four decades, across most regions in Svalbard except for a handful of glaciers draining the ice caps Vestfonna and Austfonna on Nordaustlandet. In addition, we identified complex patterns of glacier surging events overlaid with seasonal calving cycles. These data and findings provide insights into understanding glacier calving mechanisms and drivers. This new dataset can help improve estimates of glacier frontal ablation as a component of the integrated mass balance of marine-terminating glaciers. [ABSTRACT FROM AUTHOR]

Many scientific and societal questions that draw on ice sheet modeling necessitate sampling a wide range of potential climatic changes and realizations of internal climate variability. For example, coastal planning literature demonstrates a demand for probabilistic sea level projections with quantified uncertainty. Further, robust attribution of past and future ice sheet change to specific processes or forcings requires a full understanding of the space of possible ice sheet behaviors. The wide sampling required to address such questions is computationally infeasible with sophisticated numerical climate models at the resolution required to accurately force ice sheet models. Stochastic generation of climate forcing of ice sheets offers a complementary alternative. Here, we describe a method to construct a stochastic generator for ice sheet surface mass balance varying in time and space. We demonstrate the method with an application to Greenland Ice Sheet surface mass balance for 1980–2012. We account for spatial correlations among glacier catchments using sparse covariance techniques, and we apply an elevation-dependent downscaling to recover gridded surface mass balance fields suitable for forcing an ice sheet model while including feedback from changing ice sheet surface elevation. The efficiency gained in the stochastic method supports large-ensemble simulations of ice sheet change in a new stochastic ice sheet model. We provide open source Python workflows to support use of our stochastic approach for a broad range of applications. [ABSTRACT FROM AUTHOR]

Between 15 and 19 March 2022, East Antarctica experienced an exceptional heat wave with widespread 30°–40°C temperature anomalies across the ice sheet. In Part I, we assessed the meteorological drivers that generated an intense atmospheric river (AR) that caused these record-shattering temperature anomalies. Here, we continue our large collaborative study by analyzing the widespread and diverse impacts driven by the AR landfall. These impacts included widespread rain and surface melt that was recorded along coastal areas, but this was outweighed by widespread high snowfall accumulations resulting in a largely positive surface mass balance contribution to the East Antarctic region. An analysis of the surface energy budget indicated that widespread downward longwave radiation anomalies caused by large cloud-liquid water contents along with some scattered solar radiation produced intense surface warming. Isotope measurements of the moisture were highly elevated, likely imprinting a strong signal for past climate reconstructions. The AR event attenuated cosmic ray measurements at Concordia, something previously never observed. Last, an extratropical cyclone west of the AR landfall likely triggered the final collapse of the critically unstable Conger Ice Shelf while further reducing an already record low sea ice extent. Significance Statement: Using our diverse collective expertise, we explored the impacts from the March 2022 heat wave and atmospheric river across East Antarctica. One key takeaway is that the Antarctic cryosphere is highly sensitive to meteorological extremes originating from the midlatitudes and subtropics. Despite the large positive temperature anomalies driven from strong downward longwave radiation, this event led to huge amounts of snowfall across the Antarctic interior desert. The isotopes in this snow of warm airmass origin will likely be detectable in future ice cores and potentially distort past climate reconstructions. Even measurements of space activity were affected. Also, the swells generated from this storm helped to trigger the final collapse of an already critically unstable Conger Ice Shelf while further degrading sea ice coverage. [ABSTRACT FROM AUTHOR]

Maritime glaciers in the southeastern Tibetan Plateau (TP) have experienced important changes in mass and dynamics over the past decades, challenging the regional water supply and glacier-related hazards. However, knowledge about long-term variations in the surface velocity and mass balance of maritime glaciers remains incomplete due to the lack of representative observations in the southeastern TP. In this study, offset tracking is employed to measure spatiotemporal variation in the surface velocity of the Hailuogou Glacier (HLG) in Mount Gongga of the southeastern TP using Sentinel-1A imagery, while the time series of the HLG mass balance is reconstructed since 1950 by a physically based energy–mass balance model. Our satellite-based results find that HLG surface velocity shows significant spatial heterogeneity with a double-peak pattern along the flow line, and sustained slowdown below the icefall zone has been observed during the past nearly 40 years, although the icefall zone and the area above it have become relatively active. Our modeling indicates a persistent increase in mass loss over the last seven decades with an average rate of −0.58 m water equivalent (w.e.) year−1, which has accelerated in the past two decades. Sustained slowdown on the glacier is concomitant with pronounced negative mass balance, thereby enhancing glacier wastage in recent decades. The long-term trend in HLG mass loss is mainly driven by an increase in positive air temperature that decreases surface albedo and solid precipitation ratio and increases longwave incoming radiation, besides the influence of supraglacial debris cover. Large-scale atmospheric circulation patterns in the Eurasian region provide important implications for regional-to-local climate variability, unsustainably intensifying the trend of the negative mass balance of the HLG in the southeastern TP in the past two decades. [ABSTRACT FROM AUTHOR]

Wind-driven snow redistribution affects the glacier mass balance by eroding or depositing mass from or to different parts of the glacier's surface. High-resolution observations are used to test the ability of large-eddy simulations as a tool for distributed mass balance modeling. We present a case study of observed and simulated snow redistribution over Hintereisferner glacier (Ötztal Alps, Austria) between 6 and 9 February 2021. Observations consist of three high-resolution digital elevation models (Δx=1 m) derived from terrestrial laser scans taken shortly before, directly after, and 15 h after snowfall. The scans are complemented by datasets from three on-site weather stations. After the snowfall event, we observed a snowpack decrease of 0.08 m on average over the glacier. The decrease in the snow depth can be attributed to post-snowfall compaction and the wind-driven redistribution of snow. Simulations were performed with the Weather Research and Forecasting (WRF) model at Δx=48 m with a newly implemented snow drift module. The spatial patterns of the simulated snow redistribution agree well with the observed generalized patterns. Snow redistribution contributed -0.026 m to the surface elevation decrease over the glacier surface on 8 February, resulting in a mass loss of -3.9 kg m -2 , which is on the same order of magnitude as the observations. With the single case study we cannot yet extrapolate the impact of post-snowfall events on the seasonal glacier mass balance, but the study shows that the snow drift module in WRF is a powerful tool to improve knowledge on wind-driven snow redistribution patterns over glaciers. [ABSTRACT FROM AUTHOR]

Debris-covered glaciers exist in many mountain ranges and play an important role in the regional water cycle. However, modelling the surface mass balance, runoff contribution and future evolution of debris-covered glaciers is fraught with uncertainty as accurate observations on small-scale variations in debris thickness and sub-debris ice melt rates are only available for a few locations worldwide. Here we describe a customised low-cost unoccupied aerial vehicle (UAV) for high-resolution thermal imaging of mountain glaciers and present a complete open-source pipeline that facilitates the generation of accurate surface temperature and debris thickness maps from radiometric images. First, a radiometric orthophoto is computed from individual radiometric UAV images using structure-from-motion and multi-view-stereo techniques. User-specific calibration and correction procedures can then be applied to the radiometric orthophoto to account for atmospheric and environmental influences that affect the radiometric measurement. The thermal orthophoto reveals distinct spatial variations in surface temperature across the surveyed debris-covered area. Finally, a high-resolution debris thickness map is derived from the corrected thermal orthophoto using an empirical or inverse surface energy balance model that relates surface temperature to debris thickness and is calibrated against in situ measurements. Our results from a small-scale experiment on the Kanderfirn (also known as Kander Neve) in the Swiss Alps show that the surface temperature and thickness of a relatively thin debris layer (ca. 0–15 cm) can be mapped with high accuracy using an empirical or physical model. On snow and ice surfaces, the mean deviation of the mapped surface temperature from the melting point (∼ 0 ∘ C) was 0.6 ± 2.0 ∘ C. The root-mean-square error of the modelled debris thickness was 1.3 cm. Through the detailed mapping, typical small-scale debris features and debris thickness patterns become visible, which are not spatially resolved by the thermal infrared sensors of current-generation satellites. The presented approach paves the way for comprehensive high-resolution supraglacial debris thickness mapping and opens up new opportunities for more accurate monitoring and modelling of debris-covered glaciers. [ABSTRACT FROM AUTHOR]

A 54-year-long series of continuous instrumental measurements of mass balance and its main components has already been accumulated at the Djankuat Glacier, which is representative of the Caucasus and the most studied glacier in Russia. The anomalies of these indicators in 2017/2018–2020/2021 were evaluated against an analysis of meteorological reasons that predetermined them. Each of the four balance years under consideration represents a particular anomaly of varying severity. As for conditions of mass income, three years saw accumulation higher than average, and in one year (2018/2019) it approached the norm. As for summer ablation conditions, similarly, in one season (2019) the melting differed from the average only slightly, but in the other three it was much higher. Consequently, in one year (2020/2021) the state of the glacier was close to normal, in another (2017/2018) the budget situation was much more favorable for Djankuat, and in the other two the final losses significantly exceeded the average annual mass loss rate. At the same time, in 2019/2020, an absolute record of ablation since the beginning of monitoring in 1967/1968 was recorded (4360 mm w.e.). Nevertheless, although negative mass balance values continue to be recorded annually, signs of an inevitable slowdown in the rate of glacier degradation in the Caucasus have appeared in the last 4-year-long period: the continued growth of winter snow accumulation overlaps the ongoing intensification of summer melting. The growth of debris cover in terms of area and thickness also affects this mass loss slowdown to some extent. This inhibits ablation, exerting a heat-insulating effect. Because of this, the congruence of mass balance parameters vs. altitude curves is distorted. Also, a tendency toward increasing annual glacier mass turnover was revealed for the last half-century. This fact gradually increases the energy of glaciation and indirectly indicates a weakening of continentality in the climate of the Caucasian highlands. [ABSTRACT FROM AUTHOR]

Glaciers are crucial sources of freshwater in particular for the arid lowlands surrounding High Mountain Asia. To better constrain glacio-hydrological models, annual, or even better, seasonal information about glacier mass changes is highly beneficial. In this study, we evaluate the suitability of very-high-resolution Pléiades digital elevation models (DEMs) to measure glacier mass balance at annual and seasonal scales in two regions of High Mountain Asia (Muztagh Ata in Eastern Pamirs and parts of western Nyainqêntanglha, south-central Tibetan Plateau), where recent estimates have shown contrasting glacier behaviour. The average annual mass balance in Muztagh Ata between 2019 and 2022 was -0.07 ± 0.20 m w.e. a -1 , suggesting the continuation of a recent phase of slight mass loss following a prolonged period of balanced mass budgets previously observed. The mean annual mass balance in western Nyainqêntanglha was highly negative for the same period (-0.60 ± 0.15 m w.e. a -1), suggesting increased mass loss rates compared to the approximately previous 5 decades. The 2022 winter (+0.13 ± 0.24 m w.e.) and summer (-0.35 ± 0.15 m w.e.) mass budgets in Muztagh Ata and western Nyainqêntanglha (-0.03 ± 0.27 m w.e. in winter; -0.63 ± 0.07 m w.e. in summer) suggest winter- and summer-accumulation-type regimes, respectively. We support our findings by implementing the Sentinel-1-based Glacier Index to identify the firn and wet-snow areas on glaciers and characterize the accumulation type. The good match between the geodetic and Glacier Index results supports the potential of very-high-resolution Pléiades data to monitor mass balance at short timescales and improves our understanding of glacier accumulation regimes across High Mountain Asia. [ABSTRACT FROM AUTHOR]

The largest regional drivers of current surface elevation increases in the Antarctic Ice Sheet are associated with ice flow reconfiguration in previously active ice streams, highlighting the important role of ice dynamics in mass balance calculations. Here, we investigate controls on the evolution of the flow configuration of the Vanderford and Totten glaciers – key outlet glaciers of the Aurora Subglacial Basin (ASB) – the most rapidly thinning region of the East Antarctic Ice Sheet (EAIS). We synthesise factors that influence the ice flow in this region and use an ice sheet model to investigate the sensitivity of the catchment divide location to changes in surface elevation due to thinning at the Vanderford Glacier (VG) associated with ongoing retreat and thickening at the Totten Glacier (TG) associated with an intensification of the east–west snowfall gradient. The present-day catchment divide between the Totten and Vanderford glaciers is not constrained by the geology or topography but is determined by the large-scale ice sheet geometry and its long-term evolution in response to climate forcing. Furthermore, the catchment divide migrates under relatively small changes in surface elevation, leading to ice flow and basal water piracy from the Totten to the Vanderford Glacier. Our findings show that ice flow reconfigurations occur not only in regions of West Antarctica like the Siple Coast but also in the east, motivating further investigations of past, and the potential for future, ice flow reconfigurations around the whole Antarctic coastline. Modelling of ice flow and basal water piracy may require coupled ice sheet thermomechanical and subglacial hydrology models constrained by field observations of subglacial conditions. Our results have implications for ice sheet mass budget studies that integrate over catchments and the validity of the zero flow assumption when selecting sites for ice core records of past climate. [ABSTRACT FROM AUTHOR]

The mass loss of glaciers outside the polar ice sheets has been accelerating during the past several decades and has been contributing to global sea-level rise. However, many of the mechanisms of this mass loss process are not well understood, especially the calving dynamics of marine-terminating glaciers, in part due to a lack of high-resolution calving front observations. Svalbard is an ideal site to study the climate sensitivity of glaciers as it is a region that has been undergoing amplified climate variability in both space and time compared to the global mean. Here we present a new high-resolution calving front dataset of 149 marine-terminating glaciers in Svalbard, comprising 124919 glacier calving front positions during the period of 1985–2023 (https://doi.org/10.5281/zenodo.8399899) (Li et al., 2023). This dataset was generated using a novel automated deep learning framework and multiple optical and SAR satellite images from Landsat, Terra-ASTER, Sentinel-2, and Sentinel-1 satellite missions. The overall calving front mapping uncertainty across Svalbard is 46 ± 21 m. The newly derived calving front dataset agrees well with recent decadal calving front observations between 2000 and 2020 (Kochtitzky and Copland, 2022) and an annual calving front dataset between 2008 and 2022 (Moholdt et al., 2022). The R² of the glacier calving front change rates between our product and the latter is 0.98, indicating an excellent match. Using this new calving front dataset, we identified widespread calving front retreats during the past three decades, across most regions in Svalbard except for a handful of glaciers draining the ice caps Vestfonna and Austfonna on Nordaustlandet. In addition, we identified complex patterns of glacier surging events overlaid with seasonal calving cycles. These data and findings provide insights into understanding glacier calving mechanisms and drivers. This new dataset can help improve estimates of glacier frontal ablation as a component of the integrated mass balance of marine-terminating glaciers. [ABSTRACT FROM AUTHOR]

Glaciers and ice caps in Norway are presently undergoing mass loss, areal reduction, and frontal retreat, mainly a result of increased summer ablation due to rising summer temperatures over Scandinavia, especially after 2000 CE. In this paper, the glacier mass-balance response of 10 Norwegian glaciers with continuous mass balance observations (>10 years) to climate projections from 1971–2000 to 2071–2100 have been estimated. Projected changes in mean summer temperature and mean winter precipitation from 1971–2000 to 2071–2100, applying the RCP8.5 emission scenario for five different regions in Norway; 'Sogn og Fjordane' and 'Hordaland', now Vestland County in western Norway, 'Oppland', now part of Innlandet County in eastern Norway, and Nordland County and Finnmark County, both in Northern Norway), range between +3.5°C and +5.0°C, and between +5% and +25%, respectively. These climate projections have been converted (by linear regression with overlapping observational mass-balance data) into specific surface glacier mass balance [winter balance (B w), summer balance (B s), and annual balance (B a) for 10 glaciers in Norway with mass-balance series [Ålfotbreen, Nigardsbreen (part of Jostedalsbreen), Austdalsbreen (part of Jostedalsbreen), Rembesdalskåka (part of Hardangerjøkulen), Blomstølskardsbreen (part of Søre Folgefonna), Storbrean, Hellstugubrean, Gråsubrean, Engabreen (part of Vestre Svartisen, Langfjordjøkelen (data: http://glacier.nve.no/glacier/viewer/ci/no/) yielding a total, cumulative surface glacier mass loss from 2000 to 2100 CE in the range of -85.2 ± 4 to -197.3 ± 10 m water equivalents. The estimated changes in equilibrium-line altitudes (ELAs), in the range of 230 ± 10 to 630 ± 30 m, indicate that the mean ELA may reach the upper part of 7 of the 10 glaciers included in this study [Ålfotbreen, Austdalsbreen, Rembesdalskåka, Blomstølskardsbreen, Gråsubrean, Engabreen and Langfjordjøkelen] by the end of the 21st century. The projected glacier mass loss and ELA rise, and thus changes in glacier length, area and volume, will most likely have profound consequences for future glacier hydrology (runoff), hydropower production, wildlife, ecosystems, glacier hazards, and tourism. [ABSTRACT FROM AUTHOR]

Despite their extreme elevation, glaciers on the Tibetan Plateau are losing mass in response to atmospheric warming, the pattern of which purportedly reflects regional contrasts in climate. Here we examine the evolution of glaciers along ~500 km of the Tanggula Shan, Central-Eastern Tibetan Plateau. Using remotely sensed datasets, we quantified changes in glacier mass, area and surface velocity, and compared these results to time series of meteorological observations, in order to disentangle drivers of glacier mass loss since the 1960s. Glacier mass loss has increased (from −0.21 ± 0.12 m w.e. a−1 in 1960s–2000s, to −0.52 ± 0.18 m w.e. a−1 in 2000s–2015/18) in association with pervasive positive temperature anomalies (up to 1.85°C), which are pronounced at the end of the now lengthened ablation season. However, glacier mass budget perturbations do not mirror the magnitude of temperature anomalies in sub-regions, thus additional factors have heightened glacier recession. We show how proglacial lake expansion and glacier surging have compounded glacier recession over decadal/multi-decadal time periods, and exert similar influence on glacier mass budgets as temperature changes. Our results demonstrate the importance of ice loss mechanisms not often incorporated into broad-scale glacier projections, which need to be better considered to refine future glacier runoff estimates. [ABSTRACT FROM AUTHOR]

In the Tien Shan, few modelling studies exist that examine in detail how individual ice bodies are responding to climate change. Nonetheless, earlier research demonstrated that the glacier response to climate change in this mountain range is heterogeneous. Here, we use several measurements and reconstructions of the ice thickness, surface elevation, surface mass balance, and ice temperature to model in depth six different ice bodies in the Kyrgyz Tien Shan: five valley glaciers and one ice cap. The selected ice masses are located in different sub-regions of the Tien Shan with different climatic and topographic settings, and they are all characterized by detailed recent glaciological measurements. A three-dimensional higher-order thermomechanical ice-flow model is calibrated and applied to simulate the evolution of the ice masses since the end of the Little Ice Age (1850) and to make a prognosis of the future evolution up to 2100 under different Coupled Model Intercomparison Project Phase 6 (CMIP6) shared socioeconomic pathway (SSP) climate scenarios. The results reveal a strong retreat of most of the ice masses under all climate scenarios, albeit with notable variations in both timing and magnitude. These can be related to the specific climate regime of each of the ice bodies and their geometry. Under a moderate warming scenario, the ice masses characterized by a limited elevation range undergo complete disappearance, whereas the glaciers with a larger elevation range manage to preserve some ice at the highest altitudes. Additionally, our findings indicate that glaciers that primarily receive precipitation during the late spring and summer months exhibit a more rapid retreat in response to climate change, while the glaciers experiencing higher precipitation levels or more winter precipitation remain for a longer duration. Projections concerning the overall glacier runoff reveal that the maximum water discharge from the ice masses is expected to occur around or prior to the middle of the 21st century and that the magnitude of this peak is contingent upon the climate scenario, with a higher warming scenario resulting in a higher peak. [ABSTRACT FROM AUTHOR]

Due to the scarcity of observational data and the intricate precipitation–runoff relationship, individually applying physically based hydrological models and machine learning (ML) techniques presents challenges in accurately predicting floods within data-scarce glacial river basins. To address this challenge, this study introduces an innovative hybrid model that synergistically harnesses the strengths of multi-source remote sensing data, a physically based hydrological model (i.e., Spatial Processes in Hydrology (SPHY)), and ML techniques. This novel approach employs MODIS snow cover data and remote sensing-derived glacier mass balance data to calibrate the SPHY model. The SPHY model primarily generates baseflow, rain runoff, snowmelt runoff, and glacier melt runoff. These outputs are then utilized as extra inputs for the ML models, which consist of Random Forest (RF), Gradient Boosting (GDBT), Long Short-Term Memory (LSTM), Deep Neural Network (DNN), Support Vector Machine (SVM) and Transformer (TF). These ML models reconstruct the intricate relationship between inputs and streamflow. The performance of these six hybrid models and SPHY model is comprehensively explored in the Manas River basin in Central Asia. The findings underscore that the SPHY-RF model performs better in simulating and predicting daily streamflow and flood events than the SPHY model and the other five hybrid models. Compared to the SPHY model, SPHY-RF significantly reduces RMSE (55.6%) and PBIAS (62.5%) for streamflow, as well as reduces RMSE (65.8%) and PBIAS (73.51%) for floods. By utilizing bootstrap sampling, the 95% uncertainty interval for SPHY-RF is established, effectively covering 87.65% of flood events. Significantly, the SPHY-RF model substantially improves the simulation of streamflow and flood events that the SPHY model struggles to capture, indicating its potential to enhance the accuracy of flood prediction within data-scarce glacial river basins. This study offers a framework for robust flood simulation and forecasting within glacial river basins, offering opportunities to explore extreme hydrological events in a warming climate. [ABSTRACT FROM AUTHOR]

During the second half of the 20th century, the Antarctic Peninsula region has undergone a long and sustained warming period, followed by a shorter but also sustained cooling period, and then a very recent return to warming conditions. All of these have profoundly impacted the glaciers peripheral to the Antarctic Peninsula. This paper focuses on the analysis of the surface mass balance monitoring of such glaciers by the glaciological method, complemented by the analysis of mass-balance estimates by geodetic methods, as well as frontal ablation estimates. We aim to summarize the current knowledge and outline the main challenges faced by investigating the mass balance of such peripheral glaciers and their current contribution to sea-level rise. [ABSTRACT FROM AUTHOR]

Alpine glaciers are sensitive indicators of regional climate change, which can affect regional ecological stability and social development. Variations in glacier mass balance (GMB) are an important parameter in studying glacier change. In this study, data from the Ice, Cloud, and Land Elevation Satellite-1 (ICESat-1), the Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2), and CryoSat-2 (Ku-band) were combined, and high-resolution ALOS DEM was employed to denoise. After that, the polynomial fitting method was used to analyze the characteristics of glacier surface elevation (GSE) variations from 2003–2020 in the Tomur Peak Region of the Central Asian Tianshan Mountains and the regional GMB was calculated. Research results showed that: (1) From 2003–2020, the GSE of the Tomur Peak Region had an overall −8.95 ± 4.48 m variation, the average rate of which was −0.53 ± 0.26 m/yr (/yr is /year). Overall, elevations of most glaciers in the Tomur Peak Region had downward trends, with a rate of change of −0.5 to 0 m/yr. The fastest rate of elevation decline in the Koxkar Glacier Tongue was −1.5 m/yr. The elevation of some altimetric points in the Eastern Tomur Peak Region showed a rising state, with a maximum rate of variation of 1.0 m/yr. (2) From 2003–2020, the average GMB in the Tomur Peak Region was −1.51 ± 0.04 Gt/yr. In the region of elevation below 4000 m, small glaciers dominated, with a GMB of −0.61 ± 0.04 Gt/yr. With increasing elevation, the melting rate of glaciers gradually slowed down, but overall, the mass balance remained in a state of decline. (3) Climate was the main driving factor of GMB change in the study area. From 2003–2020, in the Tomur Peak Region, the average annual temperature continued to increase at a rate of 0.04 ± 0.02 °C/yr, and this was the main influencing factor for the negative GMB in the Tomur Peak Region. In the same period, the annual precipitation showed a rising trend with a linear variation rate of 0.12 ± 0.06 mm/yr, and the rising precipitation was the influencing factor for the gradually slowing change in the GMB in the study area. [ABSTRACT FROM AUTHOR]

Debris-covered glaciers are a common feature of the mountain cryosphere in the southeastern Tibetan Plateau. A better understanding of these glaciers is necessary to reduce the uncertainties in regional water resource variability and to anticipate potential cryospheric risks. In this study, we quantified the seasonal thinning and surface mass balance patterns of two neighbouring debris-covered glaciers (23K Glacier and 24K Glacier) in the southeastern Tibetan Plateau with four unpiloted aerial vehicle surveys and in situ measurements. We observed that the thinning of 23K Glacier was ∼2 –7 times greater than that of the 24K Glacier for annual and cold periods. The surface velocity of the 24K Glacier is higher than that of the 23K Glacier (∼5 –6 times) for all periods. In contrast with the thinning patterns, the surface mass balance patterns of the two glaciers closely agree across the different periods. We found that the surface mass balance distribution strongly correlated with the spatial distribution of debris thickness for both glaciers. Ice cliffs and supraglacial ponds are prevalent on the surface of these glaciers (∼4.4 – 7.2±0.5 %) and enhance melt overall (enhancement factor: ∼2.5) but do not control the surface mass balance pattern of either glacier. This comparison study of two neighbouring glaciers confirms the significance of both glacier dynamics and debris thickness in controlling thinning and melt for the different debris-covered glaciers of the southeastern Tibetan Plateau in the context of climate change. [ABSTRACT FROM AUTHOR]

Earth system models (ESMs) can help to improve the understanding of climate-induced cryospheric–hydrological impacts in complex mountain regions, such as High Mountain Asia (HMA). Coarse ESM grids, however, have difficulties in representing cryospheric–hydrological processes that vary over short distances in complex mountainous environments. Variable-resolution (VR) ESMs can help to overcome these limitations through targeted grid refinement. This study investigates the ability of the VR Community Earth System Model (VR-CESM) to simulate cryospheric–hydrological variables such as the glacier surface mass balance (SMB) over HMA. To this end, a new VR grid is generated, with a regional grid refinement up to 7 km over HMA. Two coupled atmosphere–land simulations are run for the period 1979–1998. The second simulation is performed with an updated glacier cover dataset and includes snow and glacier model modifications. Comparisons are made to gridded outputs derived from a globally uniform 1 ∘ CESM grid, observation-, reanalysis-, and satellite-based datasets, and a glacier model forced by a regional climate model (RCM). Climatological biases are generally reduced compared to the coarse-resolution CESM grid, but the glacier SMB is too negative relative to observation-based glaciological and geodetic mass balances, as well as the RCM-forced glacier model output. In the second simulation, the SMB is improved but is still underestimated due to cloud cover and temperature biases, missing model physics, and incomplete land–atmosphere coupling. The outcomes suggest that VR-CESM could be a useful tool to simulate cryospheric–hydrological variables and to study climate change in mountainous environments, but further developments are needed to better simulate the SMB of mountain glaciers. [ABSTRACT FROM AUTHOR]

Long-term climate records which help decipher past climate variability and its impact are scarce in the tough terrain of the Himalayan region. Therefore, in order to fill the climate data gap and understand the glacier climate linkage, we developed a 231 year long (1785–2015 CE) March–June temperature record using ring-width chronology of Himalayan fir (Abies pindrow (Royle ex D.Don) Royle) for the Din Gad valley, Dokriani glacier region, Uttarkashi, Uttarakhand, in the Western Himalaya. The Din Gad, originating from the Dokriani glacier, is a meltwater river contributing to Bhagirathi catchment in the headwaters of the socio-economically vital Ganga River. The 21-year running mean of the temperature record showed 1978–1998 CE as the coldest period followed by 1925–1945 CE, and 1890–1910 CE as the warmest period followed by 1946–1966 CE over the entire time series. The reconstruction matches well with tree-ring based temperature records available from the Garhwal Himalaya. It also shows similarity to tree-ring based temperature reconstructions from the Western Himalaya, Nepal, Tibetan Plateau and Bhutan, thus displaying a regional scale climate signal. The low frequency fluctuation patterns of the March–June temperature also matches with Asia and Northern hemisphere temperature records. Reconstructed March–June temperature record showed a statistically negligible warming temperature trend during 1901–1989 CE in the 20th century. It, however, captured a warming spike from 1990s CE which continues rising into the 21st century, which is also evident in the Northern hemisphere temperature record. Moreover, temperature rise is not anomalous in the past 231 years and well within range of the rest of the series. The present temperature record exclusively from the glacier region revealed a strong linkage with the benchmark Dokriani glacier's winter mass balance (November–April) revealing mass loss (gain) episodes occurred in warm (cool) phases. This first such record from the glacier valleys in Ganga headwaters would be of great value at providing insight into past climate variability and glacier behaviour with respect to climate change in long term perspective, and thus would provide valuable information for water resource management in light of climate change. [ABSTRACT FROM AUTHOR]

Despite projections of extreme reduction in glacier volume in the Qilian Mountains by the end of the century, comprehensive studies of regional glacier-wide mass and energy balance characteristics remain lacking. This study undertook a comparative analysis of the surface energy and mass balance characteristics of Laohugou glacier No. 12 (LHG glacier) in the Shule River Basin (western Qilian Mountains) and Bailanghe glacier No. 21 (BLH glacier) in the Heihe River Basin (middle Qilian Mountains) based on in situ measurements from September 2020 to August 2021. During the cold season (September–April), precipitation was greater on the BLH glacier than on the LHG glacier. This resulted in a more positive mass balance on the BLH glacier during the cold season, and less melting in May–June owing to the higher incoming shortwave radiation. During the ablation season (May–August), snowfall was greater on the LHG glacier owing to its higher elevation, while melting was also greater owing to the anomalously low cloud fraction during summer 2021. The annual mass balance was notably more negative on the LHG glacier than on the BLH glacier at the same elevation below 5000 m, whereas the annual glacier-wide mass balance was just slightly more negative on the LHG glacier than on the BLH glacier because most of the area in the LHG glacier is at higher elevation. The equilibrium line altitude varied between 4900 and 5100 on the glaciers of Qilian Mountains during recent two decades, which signified that only 5.8–25.5% of the total glacier area is within the accumulation zone in the Heihe River Basin; thus, further global warming will place the regional glaciers in a very vulnerable position. [ABSTRACT FROM AUTHOR]

A distributed mass-balance model is used over a 10-year period for the re-analysis of a glaciological mass-balance time series obtained from Brewster Glacier, New Zealand. Mass-balance modelling reveals glaciological mass balance has been overestimated, with an average mass loss of −516 mm w.e. a−1 not captured by observations at the end of the ablation season, which represents 35% of the annual mass balance. While the average length of the accumulation season (199 days) remains longer than the ablation season (166 days), melting is not uncommon in the core part of the accumulation season, with 2–32% of total snowfall being melted. Refreezing of meltwater is also important, with 10% of surface and subsurface melt being refrozen in the present climate. Net radiation, driven primarily by net shortwave radiation, is the main contributor to melt energy, with melt variability mainly influenced by the turbulent heat fluxes, net longwave radiation and the heat flux from precipitation in the ablation season. Snowfalls in summer are an important moderator of melt, highlighting the critical role of the ice-albedo feedback and phase of precipitation on seasonal mass balance. A complete homogenisation of the long-term glaciological mass balance for Brewster Glacier is still required. [ABSTRACT FROM AUTHOR]

Despite their importance for regional water resource planning and as indicators of climate change, records of in situ glacier mass balance remain short and spatially sparse in the Himalaya. Here, we present an updated series of in situ mass-balance measurements from Rikha Samba Glacier, Nepal, between 2011 and 2021. The updated in situ mass balance is −0.39 ± 0.32 m w.e. for this period. We use an energy-mass balance model to extend the annual mass-balance series back to 1974. The model is forced using daily meteorological variables from ERA5-Land reanalysis data that is linearly bias-corrected using observations from an automatic weather station situated near the glacier terminus. The modeled mass balance is consistent with the in situ mass-balance series measured 2011–2021 and with previous glaciological and geodetic estimates. The model results indicate a mass balance of −0.56 ± 0.27 m w.e. a−1 over the reconstruction period of 1974–2021, which is comparable to the mass losses experienced by other Himalayan glaciers during this time. An assessment of the sensitivity of the glacier mass balance to meteorological forcing suggests that a change in temperature of ±1 K has a stronger effect on the calculated mass balance compared to a ±20% change in either precipitation, or relative humidity, or solar radiation. [ABSTRACT FROM AUTHOR]

Decadal-scale, high-resolution geodetic measurements of glacier thinning have transformed our understanding of glacier response to climate change. Annual glacier mass balance can be estimated using remote-sensing proxies like snow-line altitude. These methods require field data for calibration, which are not available for most glaciers. Here we propose a method that combines multiple remotely-sensed proxies to obtain robust estimates of the annual glacier-wide balance using only remotely-sensed decadal-scale geodetic mass balance for calibration. The method is tested on Chhota Shigri, Argentière and Saint-Sorlin glaciers in the Himalaya and the Alps between 2001 and 2020, using four remotely-sensed proxies – the snow-line altitude, the minimum summer albedo over the glacier and two statistics of normalised difference snow index over the off-glacier area around the ablation zone. The reconstructed mass balance compares favourably with the corresponding glaciological field data (correlation coefficient 0.81 − 0.90, p < 0.001; root mean squared error 0.38 − 0.43 m w.e. a−1). The method presented may be useful to study interannual variability in mass balance on glaciers where no field data are available. [ABSTRACT FROM AUTHOR]

Sediment mass‐balance analysis provides key constraints on stratigraphic architecture and its controls. We use the data‐rich Middle Jurassic Brent Delta sediment routing system in the proto‐Viking Graben, Northern North Sea, to estimate sediment budgets and mass‐balance between source areas and depositional sinks. Published studies are synthesised to provide an age‐constrained sequence stratigraphic framework, consisting of four previously defined genetic sequences (J22, J24, J26, J32). Genetic sequence J32 (3.9 Myr) records transverse progradation of basin‐margin deltas, sourced from the Shetland Platform to the west and Norwegian Landmass to the east. Genetic sequences J24 (1.1 Myr) and J26 (0.9 Myr) record the rapid progradation and subsequent aggradation of the Brent Delta along the basin axis, sourced from the uplifted Mid‐North Sea High to the south, and the western and eastern source regions. Genetic sequence J32 (2.2 Myr) records the retreat of the Brent Delta. Sediment budgets for the four genetic sequences are estimated using palaeogeographical reconstructions, isopach maps, and sedimentological analysis of core and well‐log data. The estimated net‐depositional sediment budget for the mapped Brent Delta system is 2.0–2.8 Mt/year. Temporal variations in net‐depositional sediment budget were driven by changes in tectonic boundary conditions, such as the onset of uplift before the deposition of genetic sequence J24. Over the same time period, the Shetland Platform, Norwegian Landmass and Mid‐North Sea High source regions are estimated to have supplied 2.3–5.6, 5.0–14.1, and 2.8–9.4 Mt/year of sediment, respectively, using the BQART sediment load model and independent geometrical reconstruction of eroded volumes, which are constrained by isostatic uplift estimates based on the geochemistry of syn‐depositional volcanic rocks. The net‐depositional sediment budget in the sink is an order‐of‐magnitude smaller than the total sediment budget supplied by the source regions (13.9–23 Mt/year). This discrepancy suggests that along‐shore transport by wave‐generated currents into the coeval Faroe‐Shetland Basin and/or down‐dip transport by gravity flows into the coeval western Møre Basin played a key role in redistributing sediments away from the Brent Delta system. [ABSTRACT FROM AUTHOR]

The Tibetan Plateau (TP) has the largest number of high-altitude glaciers on Earth. As a source of major rivers in Asia, this region provides fresh water to more than one billion people. Any terrestrial water storage (TWS) changes there have major societal effects in large parts of the continent. Due to the recent acceleration in global warming, part of the water environment in TP has become drastically unbalanced, with an increased risk of water disasters. We quantified secular and monthly glacier-mass-balance and TWS changes in water basins from April 2002 to December 2021 through the Gravity Recovery and Climate Experiment and its Follow-on satellite mission (GRACE/GRACE-FO). Adequate data postprocessing with destriping filters and gap filling and two regularization methods implemented in the spectral and space domain were applied. The largest glacier-mass losses were found in the Nyainqentanglha Mountains and Eastern Himalayas, with rates of −4.92 ± 1.38 Gt a−1 and −4.34 ± 1.48 Gt a−1, respectively. The Tien Shan region showed strong losses in its eastern and central parts. Furthermore, we found small glacier-mass increases in the Karakoram and West Kunlun. Most of the glacier mass change can be explained by snowfall changes and, in some areas, by summer rainfall created by the Indian monsoon. Major water basins in the north and south of the TP exhibited partly significant negative TWS changes. In turn, the endorheic region and the Qaidam basin in the TP, as well as the near Three Rivers source region, showed distinctly positive TWS signals related to net precipitation increase. However, the Salween River source region and the Yarlung Zangbo River basin showed decreasing trends. We suggest that our new and improved TWS-change results can be used for the maintenance of water resources and the prevention of water disasters not only in the TP, but also in surrounding Asian countries. They may also help in global change studies. [ABSTRACT FROM AUTHOR]

A large number of different-sized lakes exist in the inland area of the Tibetan Plateau (TP), which are examples of the important connection between the atmosphere and hydrosphere through the analysis of lake surface convergence and evaporation processes. The evaporation level changes that occur in middle–large-sized lakes (surface area > 50 km2) in the area directly influence the regional mass and energy balance values, atmospheric boundary layer heat and humidity structures, and weather processes occurring in the lower-reach areas. The studies conducted in the literature at present, concerning lake evaporation processes, generally overlook the differences in lake heat storage behavior due to the reduced amount of data in the literature concerning lake bathymetry. According to the in situ bathymetric data obtained for 68 middle–large-sized lakes in the inner basin of the TP, in this study, we calculated their heat storage (G) change values by using the different vertical-depth water-temperature-change integral method, and we established a regression equation for the heat storage and lake surface net radiation values for 68 lakes. The evaporation rates of 134 middle–large-sized lakes larger than 50 km2 in the inland are of the TP were calculated by obtaining the G regression result and adopting it into the Penman model, as well as estimating the evaporation losses of theses 134 lakes from 2002 to 2018. The result shows that the annual average evaporation rate for these lakes is 927.39 mm/year, with an insignificant upward trend (0.10 mm/year). This method achieved good accuracy compared with the Bowen ratio method, which estimates the evaporation rate during the ice-free season, with a high correlation coefficient (R) value of 0.95 and least root mean square error (RMSE) value of 61 mm. The annual mean evaporation rate can be divided into the southern and northern lake groups along a 34°N line with a difference of 314.41 mm/year. The annual average evaporation volume of these lakes was 25.02 km3 and showed an upward trend of 0.35 km3/year. Among them, the annual average evaporation volume contribution ratio of level-1 lakes (50 km2 ≤ lake's area < 100 km2, 61 lakes) was 14.04%, showing an upward trend, and the contribution of level-3 lakes (lake's area ≥ 500 km2, 10 lakes) was 41.50%, showing a downward trend. There were no obvious changes in the level-2 lakes (100 km2 ≤ lake's area < 500 km2, 63 lakes), which maintained at the same level in approximately 44.46%. Air temperature is the most important factor affecting the evaporation rate of lakes, while the lake surface area is the main factor affecting lake evaporation volume. Our study, considering the actual lake heat storage value, provides a useful reference for further improving lake water budget balance values and watershed hydrologic features in the inland closed lakes located in the TP. [ABSTRACT FROM AUTHOR]

Surface topography across the marginal zone of the Greenland Ice Sheet is constantly evolving in response to changes in weather, season, climate and ice dynamics. Yet current Digital Elevation Models (DEMs) for the ice sheet are usually based on data from a multi-year period, thus obscuring these changes over time. We here present four 500-meter resolution annual (2019-2022) summer DEMs of the Greenland ice sheet marginal zone (PRODEMs). The PRODEMs cover 50km wide band from the ice edge, and they capture all outlet glaciers of the Greenland ice sheet. Each PRODEM is based on data fusion of CryoSat-2 radar altimetry and ICESat-2 laser altimetry using a regionally-varying Kriging method. They are validated using leave-one-out cross-validation, showcasing their ability to correctly represent surface elevations within the associated spatially varying prediction uncertainties, which have a median value of 1.4m. We observe a general lowering of surface elevations compared to ArcticDEM, but the spatial pattern of change is highly complex and with annual changes superimposed. The PRODEMs will enable studies of the marginal ice sheet elevation changes in great detail, temporally as well as spatially. With their high spatio-temporal resolution, the PRODEMs will be of value to a wide range of researchers and users studying ice sheet dynamics and monitoring how the ice sheet responds to changing environmental conditions. Incorporating the PRODEM surface elevations in estimates of the solid ice discharge from Greenland outlet glaciers will e.g. improve our assessment of the mass balance of the Greenland ice sheet and its interannual variability. PRODEMs from summers 2019 through 2022 are available at https://doi.org/10.22008/FK2/52WWHG (Winstrup, 2023), and we plan to annually update the product henceforth. [ABSTRACT FROM AUTHOR]

To get a better overview of atmosphere-driven mass changes at Urumqi Glacier No.1, Chinese Tien Shan, the surface energy budget and mass balance is modeled by linking the COupled Snowpack and Ice surface energy and MAss balance model (COSIMA) with in-situ measured meteorological records for the ablation period in 2018. The COSIMA is calibrated by manual optimization and the modeled results agree well with the in-situ surface temperature, snow height and seasonal mass balance. Our results reveal that Urumqi Glacier No.1 experienced a significant mass loss, with an average value of − 0.77 m w.e. over the ablation period 2018. The surface energy budget components can be classified into two categories: radiation (shortwave and longwave) and turbulent fluxes. Surface melt and solid precipitation were dominated components of mass balance. The COSIMA can reproduce the glaciological mass balance compared with other models. Sensitivity analysis showed that the mass balance was more sensitive to the temperature than precipitation, and mass loss caused by temperature increase of 1 K needed to be compensated by at least 40% precipitation increase. Air temperature during the ablation period was more important than annual precipitation in controlling mass balance changes on Urumqi Glacier No. 1. These findings will enhance our understanding of the mechanisms underlying mass balance processes of ablation period and their contribution to the acceleration of glacier retreat in Tien Shan. [ABSTRACT FROM AUTHOR]

The Arctic is undergoing increased warming compared to the global mean, which has major implications for freshwater runoff into the oceans from seasonal snow and glaciers. Here, we present high-resolution (2.5 km) simulations of glacier mass balance, runoff, and snow conditions on Svalbard from 1991–2022, one of the fastest warming regions in the world. The simulations are created using the CryoGrid community model forced by Copernicus Arctic Regional ReAnalysis (CARRA) (1991–2021) and AROME-ARCTIC forecasts (2016–2022). Updates to the water percolation and runoff schemes are implemented in the CryoGrid model for the simulations. In situ observations available for Svalbard, including automatic weather station data, stake measurements, and discharge observations, are used to carefully evaluate the quality of the simulations and model forcing. We find a slightly negative climatic mass balance (CMB) over the simulation period of - 0.08 mw.e.yr-1 but with no statistically significant trend. The most negative annual CMB is found for Nordenskiöldland (- 0.73 mw.e.yr-1), with a significant negative trend of - 0.27 mw.e. per decade for the region. Although there is no trend in the annual CMB, we do find a significant increasing trend in the runoff from glaciers of 0.14 mw.e. per decade. The average runoff was found to be 0.8 mw.e.yr-1. We also find a significant negative trend in the refreezing of - 0.13 mw.e. per decade. Using AROME-ARCTIC forcing, we find that 2021/22 has the most negative CMB and highest runoff over the 1991–2022 simulation period investigated in this study. We find the simulated climatic mass balance and runoff using CARRA and AROME-ARCTIC forcing are similar and differ by only 0.1 mw.e.yr-1 in climatic mass balance and by 0.2 mw.e.yr-1 in glacier runoff when averaged over all of Svalbard. There is, however, a clear difference over Nordenskiöldland, where AROME-ARCTIC simulates significantly higher mass balance and significantly lower runoff. This indicates that AROME-ARCTIC may provide similar high-quality predictions of the total mass balance of Svalbard as CARRA, but regional uncertainties should be taken into consideration. The simulations produced for this study are made publicly available at a daily and monthly resolution, and these high-resolution simulations may be re-used in a wide range of applications including studies on glacial runoff, ocean currents, and ecosystems. [ABSTRACT FROM AUTHOR]

Two dynamic versions of the simplified bioaccessibility extraction test (SBET) were developed—an off-line procedure and an on-line procedure coupled directly to ICP-MS. Batch, on-line, and off-line procedures were applied to simulated PM10 samples prepared by loading NIST SRM 2711A Montana II Soil and BGS RM 102 Ironstone Soil onto 45-mm TX40 filters widely used in air quality monitoring. Three real PM10 samples were also extracted. A polycarbonate filter holder was used as an extraction unit for the dynamic procedures. Arsenic, Cd, Cr, Cu, Fe, Mn, Ni, Pb, and Zn were determined in the extracts using an Agilent 7700 × ICP-MS instrument. The residual simulated PM10 samples following application of the SBET were subjected to microwave-assisted aqua regia digestion and a mass balance calculation performed with respect to digestion of a separate test portion of the SRM. Leachates were collected as subfractions for the off-line analysis or continuously introduced to the nebuliser of the ICP-MS for the on-line analysis. The mass balance was generally acceptable for all versions of the SBET. Recoveries obtained with the dynamic methods were closer to pseudototal values than those obtained in batch mode. Off-line analysis performed better than on-line analysis, except for Pb. Recoveries of bioaccessible Pb relative to the certified value in NIST SRM 2711A Montana II Soil (1110 ± 49 mg kg−1) were 99, 106, and 105% for the batch, off-line, and on-line methods, respectively. The study demonstrates that dynamic SBET can be used to measure bioaccessibility of potentially toxic elements in PM10 samples. [ABSTRACT FROM AUTHOR]

Ice sheet surface mass balance describes the net snow accumulation at the ice sheet surface. On the Antarctic ice sheet, winds redistribute snow, resulting in a surface mass balance that is variable in both space and time. Representing wind-driven snow redistribution processes in models is critical for local assessments of surface mass balance, repeat altimetry studies, and interpretation of ice core accumulation records. To this end, we have adapted Alpine3D, an existing distributed snow modeling framework, to downscale Antarctic surface mass balance to horizontal resolutions up to 1 km. In particular, we have introduced a new two-dimensional advection-based wind-driven snow redistribution module that is driven by an offline coupling between WindNinja, a wind downscaling model, and Alpine3D. We then show that large accumulation variability can be at least partially explained by terrain-induced wind speed variations which subsequently redistribute snow around rolling topography. By comparing Alpine3D to airborne-derived snow accumulation measurements within a testing domain over Pine Island Glacier in West Antarctica, we demonstrate that our Alpine3D downscaling approach improves surface mass balance estimates when compared to the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), a global atmospheric reanalysis which we use as atmospheric forcing. In particular, when compared to MERRA-2, Alpine3D reduces simulated surface mass balance root mean squared error by 23.4 mmw.e.yr-1 (13 %) and increases variance explained by 24 %. Despite these improvements, our results demonstrate that considerable uncertainty stems from the employed saltation model, confounding simulations of surface mass balance variability. [ABSTRACT FROM AUTHOR]

Mass loss at the Greenland Ice Sheet is influenced by atmospheric processes controlling its surface mass balance, and by submarine melt and calving where glaciers terminate in fjords. There, an ice mélange – a composite matrix of calved ice bergs and sea ice – may provide a buttressing force on a glacier terminus and control terminus dynamics. Kangerlussuaq Glacier is a major outlet of the Greenland Ice Sheet, for which recent major retreat events in 2004/2005 and 2016–2018 coincided with the absence of an ice mélange in Kangerlussuaq Fjord. To better understand the response of Kangerlussuaq Glacier to climatic and oceanic drivers, a 2D flowline model is employed. Results indicate that an ice mélange buttressing force exerts a major control on calving frequency and rapid retreat. When an ice mélange forms in Kangerlussuaq Fjord, it provides stabilising forces and conditions favourable for winter terminus re-advance. When it fails to form during consecutive years, model results indicate that Kangerlussuaq Glacier is primed to retreat into the large overdeepenings in Kangerlussuaq Fjord, and to terminus positions more than 30 km farther inland, implying that excessive mass loss from Kangerlussuaq Glacier by the year 2065 cannot be excluded. [ABSTRACT FROM AUTHOR]

Glaciers and ice caps (GICs) are important contributors of meltwater runoff and to global sea level rise. However, knowledge of GIC mass changes is largely restricted to the last few decades. Here we show the extent of 5327 Greenland GICs during Little Ice Age (LIA) termination (1900) and reveal that they have fragmented into 5467 glaciers in 2001, losing at least 587 km3 from their ablation areas, equating to 499 Gt at a rate of 4.34 Gt yr−1. We estimate that the long‐term mean mass balance in glacier ablation areas has been at least −0.18 to −0.22 m w.e. yr−1 and note the rate between 2000 and 2019 has been three times that. Glaciers with ice‐marginal lakes formed since the LIA termination have had the fastest changing mass balance. Considerable spatial variability in glacier changes suggest compounding regional and local factors present challenges for understanding glacier evolution. Plain Language Summary: Glaciers and ice caps of Greenland peripheral to the ice sheet are important contributors of meltwater to the oceans and to global sea‐level rise. In this study we map the extent of 5467 glaciers during the Little Ice Age (LIA) termination c. 1900 and calculate that they have lost at least 587 km3. The rate of mass change of these glaciers between 2000 and 2019 was three times more negative than the long‐term average (of 4.34 Gt yr−1) since the LIA. Lake‐terminating glaciers now lose mass the fastest compared with land‐ or marine‐terminating glaciers. Considerable spatial variability in glacier responses suggests local factors are important and makes glacier evolution complex. Key Points: Total volume loss of at least 587 km3 since the Little Ice Age (LIA) termination, equating to 499 Gt and to 1.38 mm sea level equivalentGlacier mass balance from 2000 to 2019 is three times more negative than since the LIA but five times more negative in the North regionLake‐terminating glaciers have experienced the greatest change in rate of mass loss [ABSTRACT FROM AUTHOR]

The main objectives of this review were to: (1) review different methods/techniques to assess gaseous N-losses from manure (2) review N-gaps, attributed to dinitrogen loss as the difference between directly measured N compounds summed as total N loss and indirectly measured N loss through a mass balance in livestock manure systems, and (3) provide approaches to close the N-gap. In literature, N-gaps run up to 80% of total N loss, this undermines N emission assessments and leaves a huge part of the emission unexplained. However, studies that measure N-gaps are scarcely available or are limited in their evaluation, hence more study is needed. Three approaches are introduced to research N-gaps: (1) measure N 2 through a suggested Gas Flow Soil Core (GFSC) technique and compare the sum of all measured N losses with the indirect method, (2) assume N 2 loss as being the N-gap and (3) include N 2 as an estimate based on ratios from literature. In a hypothetical example for poultry manure, assumed values for measurement error of 50% and variance due to physical differences between the experimental units of 50% led to a total standard deviation of 131% in the N-gap. Variance of N-gap was reduced with 80% point when assuming 16 vessels compared to single vessel. Using literature-based-ratios to estimate losses of N compounds led to variation of N-gap from 0.06% initial N overestimation to 26% of initial N underestimation. Future research should address this variance and apply methods to measure N 2 to close N-gaps. [Display omitted] • Considerable N-gaps occur in livestock manure systems being up to 80% of total N loss. • Main reason for N-gap is the lack of technique and accurate estimates to determine N 2. • Measuring, estimating and equaling N 2 to N-gap are three approaches to research N-gap. • N 2 loss for laying hen manure was estimated 32% of total N loss (N mass balance). • Ratios to estimate N 2 led to 0.06% overestimation and 26% underestimation of N-gap. [ABSTRACT FROM AUTHOR]

The Greenland Ice Sheet's (GrIS) firn layer buffers the ice sheet's contribution to sea level rise by storing meltwater in its pore space. However, available pore space and meltwater retention capability is lost due to ablation of the firn layer and refreezing of meltwater as near-surface ice slabs in the firn. Understanding how firn properties respond to climate is important for constraining the GrIS's future contribution to sea level rise in a warming climate. Observations of firn density provide detailed information about firn properties, but they are spatially and temporally limited. Here we use two firn models, the physics-based SNOWPACK model and the Community Firn Model configured with a semi-empirical densification equation (CFM-GSFC), to quantify firn properties across the GrIS from 1980 through 2020. We use an identical forcing (Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2) atmospheric reanalysis) for SNOWPACK and the CFM-GSFC in order to isolate firn model differences. To evaluate the models, we compare simulated firn properties, including firn air content (FAC), to measurements from the Surface Mass Balance and Snow on Sea Ice Working Group (SUMup) dataset of snow and firn density. Both models perform well (mean absolute percentage errors of 14 % in SNOWPACK and 16 % in the CFM-GSFC), though their performance is hindered by the spatial resolution of the atmospheric forcing. In the ice-sheet-wide simulations, the 1980–1995 average spatially integrated FAC (i.e., air volume in the firn) for the upper 100 m is 34 645 km 3 from SNOWPACK and 28 581 km 3 from the CFM-GSFC. The discrepancy in the magnitude of the modeled FAC stems from differences in densification with depth and variations in the sensitivity of the models to atmospheric forcing. In more recent years (2005–2020), both models simulate substantial depletion of pore space. During this period, the spatially integrated FAC across the entire GrIS decreases by 3.2 % (-66.6 km 3 yr -1) in SNOWPACK and 1.5 % (-17.4 km 3 yr -1) in the CFM-GSFC. These differing magnitudes demonstrate how model differences propagate throughout the FAC record. Over the full modeled record (1980–2020), SNOWPACK simulates a loss of pore space equivalent to 3 mm of sea level rise buffering, while the CFM-GSFC simulates a loss of 1 mm. The greatest depletion in FAC is along the margins and especially along the western margin where observations and models show the formation of near-surface, low-permeability ice slabs that may inhibit meltwater storage. [ABSTRACT FROM AUTHOR]

Shortwave radiation penetrating beneath an ice-sheet surface can cause internal melting and the formation of a near-surface porous layer known as the weathering crust, a dynamic hydrological system that provides home to impurities and microbial life. We develop a mathematical model, incorporating thermodynamics and population dynamics, for the evolution of such layers. The model accounts for conservation of mass and energy, for internal and surface-absorbed radiation, and for logistic growth of a microbial species mediated by nutrients that are sourced from the melting ice. It also accounts for potential melt–albedo and microbe–albedo feedbacks, through the dependence of the absorption coefficient on the porosity or microbial concentration. We investigate one-dimensional steadily melting solutions of the model, which give rise to predictions for the weathering crust depth, water content, melt rate, and microbial abundance, depending on a number of parameters. In particular, we examine how these quantities depend on the forcing energy fluxes, finding that the relative amounts of shortwave (surface-penetrating) radiation and other heat fluxes are particularly important in determining the structure of the weathering crust. The results explain why weathering crusts form and disappear under different forcing conditions and suggest a range of possible changes in behaviour in response to climate change. [ABSTRACT FROM AUTHOR]

Many scientific and societal questions that draw on ice sheet modelling could be best addressed by sampling a wide range of potential climatic changes and realizations of internal climate variability. For example, coastal planning literature demonstrates a demand for probabilistic sea-level projections with quantified uncertainty. Further, robust attribution of past and future ice sheet change to specific processes or forcings requires a full understanding of the space of possible ice sheet behaviors. The wide sampling required to address such questions is computationally infeasible with sophisticated numerical climate models at the resolution required to accurately force ice sheet models. Stochastic generation of climate forcing of ice sheets offers a complementary alternative. We construct a stochastic generator of Greenland Ice Sheet surface mass balance in time and space. We find that low-order autoregressive models are sufficient to accurately reproduce the interannual variability in process-model simulations of recent Greenland surface mass balance at the glacier-catchment scale. We account for spatial correlations among glacier catchments using sparse covariance techniques, and we apply an elevation-dependent downscaling to recover gridded surface mass balance fields suitable for forcing an ice sheet model while including feedbacks to ice sheet surface elevation. The efficiency gained in the stochastic method supports large ensemble simulations of ice sheet change in a new stochastic ice sheet model. We provide open source Python workflows to support use of our stochastic approach for a broad range of applications. [ABSTRACT FROM AUTHOR]

Spaceborne radar altimetry over ice sheets has exclusively focused on assessing volume changes through changes in surface elevation. For use in mass balance calculations, these measurements are supplemented with surface property information derived from regional climate models with limited large‐scale observational validation. Simultaneously, the strength at which a radar signal is reflected from the surface contains information on these same near‐surface properties. Here we show that a quantitative interpretation of European (ESA) CryoSat‐2 and French/Indian (CNES/ISRO) SARAL surface echo powers yields timeseries of pan‐Greenland wavelength‐scale roughness and surface density. Individually, the CryoSat‐2 and SARAL results strive toward providing an observational lens to enhance our current understanding of the surface processes affecting ice sheet mass balance and validate their representation by computational models. Taken together, they highlight how diversity in currently operational and future radar satellite altimetry missions can shed light on near‐surface vertical heterogeneity. Plain Language Summary: For three decades, satellite‐based radar altimeters have been measuring areas of thickening and thinning across the Greenland ice sheet (GrIS). From these measurements we can calculate a mass balance, which is an important component in quantifying Greenland's contribution to global mean sea level rise. However, this calculation requires knowing ice sheet surface density. Even though predictable with specialized computer programs, observations are needed to ensure the computer results are accurate. In this study, we use three different radar data sets from two different satellites, to measure the surface density and roughness of the GrIS based on the strength of the radar altimetry surface echoes. As expected, surface roughness increases toward the ice sheet margin and in the vicinity of faster flowing ice. Surface density varies through time and space but appears strongly influenced by melt events (e.g., summer 2012). By comparing measurements between satellites, we can start piecing together how density changes vertically in the near surface; reflecting the processes operating at the time (i.e., snow fall, melting, melt percolation). These observations are the first step toward providing a new data set for refining the computer programs used to project Greenland evolution into the future. Key Points: In‐depth knowledge of near‐surface density is crucial for altimetry ice sheet mass balance estimatesWe provide the first, monthly, pan‐Greenland observations of surface roughness and density between 2013 and 2018Dual‐frequency results reflect deposition, melting, and percolation processes through their influence on vertical density heterogeneity [ABSTRACT FROM AUTHOR]

The Jostedalsbreen ice cap is mainland Europe's largest ice cap and accommodates 20% (458 km2 in 2019) of the total glacier area of mainland Norway. Jostedalsbreen and its meltwater contribute to global sea-level rise and to local water management, hydropower and tourism economies and livelihoods. In this study, we construct a digital terrain model (DTM) of the ice cap from 1966 aerial photographs, which by comparing to an airborne LiDAR DTM from 2020, we compute changes in surface elevation and geodetic mass balances. The area mapped in both surveys cover about 3/4 of the ice cap area and 49 of 82 glaciers. The measured glacier area has decreased from 363.4 km2 in 1966 to 332.9 km2 in 2019, i.e. a change of −30 km2 or −8.4% (−0.16% a−1), which is in line with the percentage reduction in area for Jostedalsbreen as a whole. The mean geodetic mass balance over the 49 glaciers was −0.15 ± 0.01 m w.e. a−1, however, large variability is evident between glaciers, e.g. Nigardsbreen (−0.05 m w.e. a−1), Austdalsbreen (−0.28 m w.e. a−1) and Tunsbergdalsbreen (−0.36 m w.e. a−1) confirming differences also found by the glaciological records for Nigardsbreen and Austdalsbreen. [ABSTRACT FROM AUTHOR]

The Tibetan Plateau (TP) has experienced amplified warming in recent decades, causing glaciers to melt and affecting river runoff. It is well established that the southern and northern areas of the TP have responded to climate changes differently, with the north dominated by a westerly climate and the south by the Indian monsoon. While there are more glaciers in the TP than in any other region outside the polar areas, most of these glaciers are tiny, and only a limited number of them have been monitored to study mass balance and downward runoff. This study used the mass balance measured at two glaciers along with in situ and satellite data to drive a hydrological model called the Alpine Runoff Predictor that includes glacier melt to simulate glacial melting and the accompanying hydrological processes of the two glacierized basins, analyze their contributions to the river runoffs, and investigate their responses to local climate changes. The results show that the glacier meltwater in both river basins showed an increasing trend, with values of 0.001 × 108 m3 a−1 in the Kyanjing River basin and 0.0095 × 108 m3 a−1 in the Tuole River basin. However, their multi-year average contributions to the runoff were 12.5% and 5.6%, respectively. In contrast to the Tuole River basin, where runoff is increasing (0.0617 × 108 m3 a−1), the Kyanjing River basin has decreasing runoff (−0.0216 × 108 m3 a−1) as a result of decreasing precipitation. This result highlights the dominant role played by precipitation changes in the two basins under study, which are characterized by small glacier meltwater contributions. [ABSTRACT FROM AUTHOR]

In many settings, data acquisition generates outliers that can obscure inference. Therefore, practitioners often either identify and remove outliers or accommodate outliers using robust models. However, identifying and removing outliers is often an ad hoc process that affects inference, and robust methods are often too simple for some applications. In our motivating application, scientists drill snow cores and measure snow density to infer densification rates that aid in estimating snow water accumulation rates and glacier mass balances. Advanced measurement techniques can measure density at high resolution over depth but are sensitive to core imperfections, making them prone to outliers. Outlier accommodation is challenging in this setting because the distribution of outliers evolves over depth and the data demonstrate natural heteroscedasticity. To address these challenges, we present a two-component mixture model using a physically motivated snow density model and an outlier model, both of which evolve over depth. The physical component of the mixture model has a mean function with normally distributed depth-dependent heteroscedastic errors. The outlier component is specified using a semiparametric prior density process constructed through a normalized process convolution of log-normal random variables. We demonstrate that this model outperforms alternatives and can be used for various inferential tasks. [ABSTRACT FROM AUTHOR]

Accurate monitoring of surface temperature and melting on the Greenland Ice Sheet (GrIS) is important for tracking the ice sheet's mass balance as well as global and Arctic climate change. Using a moderate-resolution-imaging-spectroradiometer (MODIS)-derived land-surface-temperature (LST) data product with a resolution of 1 km from 2000 to 2020, the temporal and spatial variations of annual and seasonal 'clear-sky' surface temperature were evaluated. We also monitored summer surface melting and studied the relationship between the mass balance of the ice sheet and changes in surface temperature and melting. The results show that the mean annual LST during the study period is −24.86 ± 5.46 °C, with the highest of −22.48 ± 5.61 °C in 2010 and the lowest temperature of −26.49 ± 5.30 °C in 2015. With the change of season, the spatial variation of the ice-sheet surface temperature changes greatly. 2012 and 2019 experienced the warmest summers (−5.92 ± 4.01 °C and −6.51 ± 3.93 °C), with extreme cumulative melting detected on the ice-sheet surface (89.9% and 89.7%, respectively), and 2002 also experienced a greater extent of melting. But short period of melt in 2002 and 2019 (30.6% and 31.4%, respectively), accounted for a larger proportion, with neither the duration nor intensity of the melt reaching that of 2012. There is a strong correlation between the GrIS surface temperature and its mass balance. By fitting the relationship between surface temperature and mass balance, it was found that 93.83% (6.17%) of the ice-sheet response to surface-temperature change was via surface-mass balance (discharge and basal-mass balance). [ABSTRACT FROM AUTHOR]

Global glacier mass balance decreased rapidly over the last two decades, exceeding mass loss from the Greenland and Antarctic Ice Sheets. In Greenland, peripheral glaciers and ice caps (GICs) cover only ~5% of Greenland's area but contributed ~20% of the island's ice mass loss between 2000 and 2018. Although Greenland GIC mass loss due to surface meltwater runoff has been estimated using atmospheric models, mass lost to changes in ice discharge into oceans (i.e., dynamic mass loss) remains unquantified. We use the flux gate method to estimate discharge from Greenland's 585 marine-terminating peripheral glaciers between 1985 and 2018, and compute dynamic mass loss as the discharge anomaly relative to the 1985–98 period. Greenland GICs discharged between 2.94 ± 0.23 and 4.03 ± 0.23 Gt a−1 from 1985 to 1998, depending on the gap-filling method, and abruptly increased to 5.10 ± 0.21 Gt a−1 from 1999 to 2018. The resultant ~1–2 Gt a−1 dynamic mass loss was driven by synchronous widespread acceleration around Greenland. The mass loss came predominantly from the southeast region, which contains 39% of the glaciers. Although changes in discharge over time were small relative to surface mass-balance changes, our speed and discharge time series suggest these glaciers may quickly accelerate in response to changes in climate. [ABSTRACT FROM AUTHOR]