Surface Mass Balance Research Articles

High-resolution temporal gravity field data products: Monthly mass grids and spherical harmonics from 1994 to 2021

Since April 2002, Gravity Recovery and Climate Experiment (GRACE) and GRACE-FO (FollowOn) satellite gravimetry missions have provided precious data for monitoring mass variations within the hydrosphere, cryosphere, and oceans with unprecedented accuracy and resolution. However, the long-term products of mass variations prior to GRACE-era may allow for a better understanding of spatio-temporal changes in climate-induced geophysical phenomena, e.g., terrestrial water cycle, ice sheet and glacier mass balance, sea level change and ocean bottom pressure (OBP). Here, climate-driven mass anomalies are simulated globally at 1.0° × 1.0° spatial and monthly temporal resolutions from January 1994 to January 2021 using an in-house developed hybrid Deep Learning architecture considering GRACE/-FO mascon and SLR-inferred gravimetry, ECMWF Reanalysis-5 data, and normalized time tag information as training datasets. Internally, we consider mathematical metrics such as RMSE, NSE and comparisons to previous studies, and externally, we compare our simulations to GRACE-independent datasets such as El-Nino and La-Nina indexes, Global Mean Sea Level, Earth Orientation Parameters-derived low-degree spherical harmonic coefficients, and in-situ OBP measurements for validation.

Spatial and temporal variability of environmental proxies from the top 120 m of two ice cores in Dronning Maud Land (East Antarctica)

Abstract. The Antarctic ice sheet's future contribution to sea level rise is difficult to predict, mostly because of the uncertainty and variability of the surface mass balance (SMB). Ice cores are used to locally (kilometer scale) reconstruct SMB with a very good temporal resolution (up to sub-annual), especially in coastal areas where accumulation rates are high. The number of ice core records has been increasing in recent years, revealing an important spatial variability and different trends of SMB, highlighting the crucial need for greater spatial and temporal representativeness. We present records of density, water stable isotopes (δ18O, δD, and deuterium excess), major ions concentrations (Na+, K+, Mg2+, Ca2+, MSA, Cl−, SO42-, and NO3-), and continuous electrical conductivity measurement (ECM), as well as age models and resulting surface mass balance from the top 120 m of two ice cores (FK17 and TIR18) drilled on two adjacent ice rises located in coastal Dronning Maud Land and dating back to the end of the 18th century. Both environmental proxies and SMB show contrasting behaviors, suggesting strong spatial and temporal variability at the regional scale. In terms of precipitation proxies, both ice cores show a long-term decrease in deuterium excess (d-excess) and a long-term increase in δ18O, although less pronounced. In terms of chemical proxies, the non-sea-salt sulfate (nssSO42-) concentrations of FK17 are twice those of TIR18 and display an increasing trend on the long-term, whereas there is only a small increase after 1950 in TIR18. The SO42-/Na+ ratios show a similar contrast between FK17 and TIR18 and are consistently higher than the seawater ratio, indicating a dominant impact of the nssSO42- on the SO42- signature. The mean long-term SMB is similar for FK17 and TIR18 (0.57 ± 0.05 and 0.56 ± 0.05 mi.e.yr-1, respectively), but the annual records are very different: since the 1950s, TIR18 shows a continuous decrease while FK17 has shown an increasing trend until 1995 followed by a recent decrease. The datasets presented here offer numerous development possibilities for the interpretation of the different paleo-profiles and for addressing the mechanisms behind the spatial and temporal variability observed at the regional scale (tens of kilometers) in East Antarctica. The “Mass2Ant IceCores” datasets are available on Zenodo (https://doi.org/10.5281/zenodo.7848435; Wauthy et al., 2023).

An agenda for the future of Arctic snow research: the view from Svalbard

The Arctic region is warming at over twice the mean rate of the Northern Hemisphere and nearly four times faster than the globe since 1979. The local rate of warming is even higher in the European archipelago of Svalbard. This warming is transforming the terrestrial snow cover, which modulates surface energy exchanges with the atmosphere, accounts for most of the runoff in Arctic catchments and is also a transient reservoir of atmospherically deposited compounds, including pollutants. Improved observations, understanding and modelling of changes in Arctic snow cover are needed to anticipate the effects these changes will have on the Arctic climate, atmosphere, terrestrial ecosystems and socioeconomic factors. Svalbard has been an international hub of polar research for many decades and benefits from a well-developed science infrastructure. Here, we present an agenda for the future of snow research in Svalbard, jointly developed by a multidisciplinary community of experts. We review recent trends in snow research, identify key knowledge gaps, prioritize future research efforts and recommend supportive actions to advance our knowledge of present and future snow conditions pertaining to glacier mass balance, permafrost, surface hydrology, terrestrial ecology, the cycling and fate of atmospheric contaminants, and remote sensing of snow cover. This perspective piece addresses issues relevant to the circumpolar North and could be used as a template for other national or international Arctic research plans.

Assessment and analysis of correlation between remote sensing albedo in ablation season and annual mass balance for glaciers in High Mountain Asia

Assessment and analysis of correlation between remote sensing albedo in ablation season and annual mass balance for glaciers in High Mountain Asia

Uncertainty-Aware Learning With Label Noise for Glacier Mass Balance Modeling

Uncertainty-Aware Learning With Label Noise for Glacier Mass Balance Modeling

2000–2020年喀喇昆仑山洪扎河流域冰川高程变化数据集

Under the influence of climate warming, glaciers in High Mountain Asia (HMA) have undergone significant recession in recent decades. However, stable or even positive glacier mass budgets have been found in the Karakoram Mountains in the past decades, which is termed as “Karakoram anomaly”. The Hunza Basin, located in the western Karakoram, is a region with concentrated glacier development along the Sino-Pakistan Economic Corridor. The glacier mass balance in Hunza Basin plays an important role in water resource utilization and infrastructure construction. In this study, we used the SRTM DEM (2000) and ZY3-02 stereo images (2020) to obtain the glacier height changes in the Hunza Basin by employing the method of Geodetic (GeoTIFF format with a resolution of 30 m). This study covers an area of around 7600 km2, within the spatial extent of 36°N–37°N, 74°E–76°E. The statistical analysis of the residual movement in the non-glacial region indicates an overall accuracy of about 1.1 m. The findings suggest that glaciers experienced a moderate thinning during the early 21st century, with significant spatial variations. Factors such as debris cover, crevasses, supraglacial lake, ice cliff and glacier surge have distinct influences on glacier height changes. This dataset can be used as basic data for the research on glacier mass balance in the Karakoram, as well as water resources utilization, glacier disaster warning and climate change in the Sino-Pakistan Economic Corridor.

1990–2018年干城章嘉流域冰川物质平衡序列模拟数据集

As a result of climate change, the global glacier shrinkage was observed, and glaciers of different types have experienced distinct changes in their mass budgets. Debris-covered and/or lake-terminating glaciers have been identified in the Kanchenjunga Basin, and the dataset of their mass balance on a long-term scale could provide an important basic data for studies of specific responses of glaciers of various types to climate change. Using the rectified China Meteorological Forcing Dataset, we reconstructed the daily mass balance of glaciers in the Kanchenjunga Basin from 1990 to 2018 through the application of a surface energy balance model. We used a glacier elevation change dataset obtained through geodetic methods to validate the model estimation results, and the results align closely with previously published research results on a multi-year scale. We extracted the surface velocities of glaciers from 1989 to 2018 from the ITS_LIVE dataset and assessed their precision based on the statistical characteristics of the surface movements in the glacier-free regions. This dataset is expected to provide basic data for the studies of the responses of glaciers with different types to climate change, and can also offer valuable references for the studies of glacial hydrological effects and disasters.

Glacier Surface Velocity Variations in the West Kunlun Mts. with Sentinel-1A Image Feature-Tracking (2014–2023)

Glacier velocity is a crucial parameter in understanding glacier dynamics and mass balance, especially in response to climate change. Despite numerous studies on glaciers in the West Kunlun Mts., there is still insufficient knowledge about the details of inter- and intra-annual velocity changes under global warming. This study analyzed the glacier velocity changes in the West Kunlun Mts. using Sentinel-1A satellite data. Our results revealed that: (1) The velocity of glaciers across the region shows an increasing trend from 2014 to 2023. (2) Five glaciers were found to have been surged during the study period, among which two of them were not reported before. (3) The surges in the study region were potentially controlled through a combination of hydrological and thermal mechanisms. (4) The glacier N2, Duofeng Glacier, and b2 of Kunlun Glacier exhibit higher annual velocities (32.82 m a−1) compared to surging glaciers in quiescent phases (13.22 m a−1), and were speculated as advancing or fast-flowing glaciers.

Rapid mass losses of Urumqi River Basin glaciers, eastern Tianshan Mountains revealed from multi-temporal DEMs, 1964–2021

ABSTRACT Glacier volume changes have remarkable impacts on regional water resources; however, their estimation uncertainties remain large. In this study, the surface elevations of seven glaciers at the headwaters of the Urumqi River exhibited obvious thinning at the rate of 0.32 m a−1 from 1964 to 2021. Correlations between the ice volume change (dV) derived from surface elevation conversion and the area change (dS) were established for valley, cirque, and hanging glaciers. The ice volumes of valley, cirque, and hanging glaciers in the Urumqi River Basin decreased by 0.99 km3, 0.18 km3, and 0.59 km3, respectively, and the total ice volume decreased by 1.76 km3 during 1964–2021. Ice volume losses are accelerating. The correlation between dV and dS for valley glaciers produced a large estimation difference for the ice volume change of Urumqi Glacier No. 1 (UG1) compared with the glaciological mass balance (15%), while the estimation differences of ice volume changes for the three types of glaciers in the Urumqi River Basin were within acceptable limits (10%) based on the geodetic results. This new method has vast potential in estimating ice volume changes at the basin or regional scale, or even across western China.

Annual to seasonal glacier mass balance in High Mountain Asia derived from Pléiades stereo images: examples from the Pamir and the Tibetan Plateau

Abstract. 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.

Empirical projection of global sea level in 2050 driven by Antarctic and Greenland ice mass variations

Global mean sea level rise, driven by ice mass loss in Antarctic and Greenland Ice Sheets (AIS and GrIS), is a significant consequence of global warming. Although various ice sheet models have attempted to predict the ice mass balance and subsequent sea level changes, non-trivial disagreements between models exist. In this study, we employ an empirical approach to estimate the future (2050) ice mass changes for both ice sheets, assuming their historical patterns of ice dynamics would persist in the coming decades. To achieve this, we estimate decadal-scale ice discharge variations by subtracting the surface mass balance (SMB) from the observed ice mass changes and extrapolate linear trend and acceleration components of ice discharges up to 2050. We also consider future SMB data from Coupled Model Intercomparison Project phase 6 models to estimate net ice mass balance. Our estimates suggest that from 2021 to 2050, the global sea level rise due to AIS and GrIS ranges between 6–19 mm and 15–31 mm, respectively. Additionally, we investigate regional sea level variability resulting from geoid changes induced by ice mass changes in both regions, highlighting that heterogeneous sea level changes may cause more pronounced sea level rises in lower latitude regions, where major cities are located.

Hysteresis and orbital pacing of the early Cenozoic Antarctic ice sheet

Abstract. The hysteresis behaviour of ice sheets arises because of the different thresholds for growth and decline of a continental-scale ice sheet depending on the initial conditions. In this study, the hysteresis effect of the early Cenozoic Antarctic ice sheet to different bedrock elevations is investigated with an improved ice sheet–climate coupling method that accurately captures the ice–albedo feedback. It is shown that the hysteresis effect of the early Cenozoic Antarctic ice sheet is ∼180 ppmv or between 3.5 and 5 ∘C, depending only weakly on the bedrock elevation dataset. Excluding isostatic adjustment decreases the hysteresis effect significantly towards ∼40 ppmv because the transition to a glacial state can occur at a warmer level. The rapid transition from a glacial to a deglacial state and oppositely from deglacial to glacial conditions is strongly enhanced by the ice–albedo feedback, in combination with the elevation–surface mass balance feedback. Variations in the orbital parameters show that extreme values of the orbital parameters are able to exceed the threshold in summer insolation to induce a (de)glaciation. It appears that the long-term eccentricity cycle has a large influence on the ice sheet growth and decline and is able to pace the ice sheet evolution for constant CO2 concentration close to the glaciation threshold.

From atmospheric water isotopes measurement to firn core interpretation in Adélie Land: a case study for isotope-enabled atmospheric models in Antarctica

Abstract. In a context of global warming and sea level rise acceleration, it is key to estimate the evolution of the atmospheric hydrological cycle and temperature in polar regions, which directly influence the surface mass balance of the Arctic and Antarctic ice sheets. Direct observations are available from satellite data for the last 40 years and a few weather data since the 1950s in Antarctica. One of the best ways to access longer records is to use climate proxies in firn or ice cores. The water isotopic composition in these cores is widely used to reconstruct past temperature variations. We need to progress in our understanding of the influence of the atmospheric hydrological cycle on the water isotopic composition of ice cores. First, we present a 2-year-long time series of vapor and precipitation isotopic composition measurement at Dumont d’Urville Station, in Adélie Land. We characterize diurnal variations of meteorological parameters (temperature, atmospheric water mixing ratio (hereafter humidity) and δ18O) for the different seasons and determine the evolution of key relationships (δ18O versus temperature or humidity) throughout the year: we find that the temperature vs. δ18O relationship is dependent on synoptic events dynamics in winter contrary to summer. Then, this data set is used to evaluate the atmospheric general circulation model ECHAM6-wiso (model version with embedded water stable isotopes) in a coastal region of Adélie Land where local conditions are controlled by strong katabatic winds which directly impact the isotopic signal. We show that a combination of continental (79 %) and oceanic (21 %) grid cells leads model outputs (temperature, humidity and δ18O) to nicely fit the observations, at different timescales (i.e., seasonal to synoptic). Therefore we demonstrate the added value of long-term water vapor isotopic composition records for model evaluation. Then, as a clear link is found between the isotopic composition of water vapor and precipitation, we assess how isotopic models can help interpret short firn cores. In fact, a virtual firn core built from ECHAM-wiso outputs explains much more of the variability observed in S1C1 isotopic record than a virtual firn core built from temperature only. Yet, deposition and post-deposition effects strongly affect the firn isotopic signal and probably account for most of the remaining misfits between archived firn signal and virtual firn core based on atmospheric modeling.

Insights into the vulnerability of Antarctic glaciers from the ISMIP6 ice sheet model ensemble and associated uncertainty

Abstract. The Antarctic Ice Sheet represents the largest source of uncertainty in future sea level rise projections, with a contribution to sea level by 2100 ranging from −5 to 43 cm of sea level equivalent under high carbon emission scenarios estimated by the recent Ice Sheet Model Intercomparison for CMIP6 (ISMIP6). ISMIP6 highlighted the different behaviors of the East and West Antarctic ice sheets, as well as the possible role of increased surface mass balance in offsetting the dynamic ice loss in response to changing oceanic conditions in ice shelf cavities. However, the detailed contribution of individual glaciers, as well as the partitioning of uncertainty associated with this ensemble, have not yet been investigated. Here, we analyze the ISMIP6 results for high carbon emission scenarios, focusing on key glaciers around the Antarctic Ice Sheet, and we quantify their projected dynamic mass loss, defined here as mass loss through increased ice discharge into the ocean in response to changing oceanic conditions. We highlight glaciers contributing the most to sea level rise, as well as their vulnerability to changes in oceanic conditions. We then investigate the different sources of uncertainty and their relative role in projections, for the entire continent and for key individual glaciers. We show that, in addition to Thwaites and Pine Island glaciers in West Antarctica, Totten and Moscow University glaciers in East Antarctica present comparable future dynamic mass loss and high sensitivity to ice shelf basal melt. The overall uncertainty in additional dynamic mass loss in response to changing oceanic conditions, compared to a scenario with constant oceanic conditions, is dominated by the choice of ice sheet model, accounting for 52 % of the total uncertainty of the Antarctic dynamic mass loss in 2100. Its relative role for the most dynamic glaciers varies between 14 % for MacAyeal and Whillans ice streams and 56 % for Pine Island Glacier at the end of the century. The uncertainty associated with the choice of climate model increases over time and reaches 13 % of the uncertainty by 2100 for the Antarctic Ice Sheet but varies between 4 % for Thwaites Glacier and 53 % for Whillans Ice Stream. The uncertainty associated with the ice–climate interaction, which captures different treatments of oceanic forcings such as the choice of melt parameterization, its calibration, and simulated ice shelf geometries, accounts for 22 % of the uncertainty at the ice sheet scale but reaches 36 % and 39 % for Institute Ice Stream and Thwaites Glacier, respectively, by 2100. Overall, this study helps inform future research by highlighting the sectors of the ice sheet most vulnerable to oceanic warming over the 21st century and by quantifying the main sources of uncertainty.

Linking ground ice and glacier melting to lake volume change in the Dogai Coring watershed on the Tibetan Plateau

The Dogai Coring Lake on the Tibetan Plateau, with widespread permafrost and glaciers in the watershed, has expanded rapidly between 2014 and 2020. However, it remains unclear how much melting water from permafrost and glaciers contributes to the increase in lake volume. This study quantifies the potential amount of water released by the melt of permafrost and glaciers, and assesses the fraction of water storage gains in the lake offset by the loss from the cryosphere (i.e., ground ice and glaciers). The small baseline subset interferometric synthetic aperture radar (SBAS-InSAR) method was used to derive the spatiotemporal changes in surface deformation, which then assisted the estimation of the potential water release from ground ice melting. The glacier mass balance was estimated using a geodetic method with ASTER DEMs. Results show that 51.3 % of the region experienced subsidence with rates in the range of 5–20 mm∙a−1 during 2014–2020. Ground ice melting released water at a rate of 35.1 × 106 m3∙a−1, accounting for ∼19.3 % of the increase in lake volume. The glacier mass loss rate was estimated to be − 0.69 ± 0.04 m w.e. yr−1 (113.6 ± 7.3 × 106 m3∙a−1), contributing ∼62.5 % to the lake volume increase. Our results show that the water storage loss in the cryosphere has largely (∼80 %) offset the water storage gain in the lake, suggesting that cryospheric meltwater is an important factor regulating water budget change of the lake watershed.

Brief communication: Surface energy balance differences over Greenland between ERA5 and ERA-Interim

Abstract. We compare the main atmospheric drivers of the melt season over the Greenland Ice Sheet (GrIS) in ERA5 and ERA-Interim (ERAI) in their overlapping period 1979–2018. In summer, ERA5 differs significantly from ERAI, especially in the melt regions. Small-scale ERA5 − ERAI differences near the ice sheet's margins and over steep slopes can be explained by the different resolution, while the large-scale differences indicate a different representation of physical processes in the two reanalyses: averaged over the lower parts of the GrIS, the mean near-surface air temperature is 1 K lower, while the mean downward shortwave radiation at the surface is on average 15 W m−2 higher than in ERAI. Comparison with observational weather station data shows a significant warm bias in ERAI and, for ERA5, a significant positive bias in downward shortwave radiation. Consequently, methods that previously estimated the GrIS surface mass balance from the ERAI surface energy balance need to be carefully recalibrated before converting to ERA5 forcing.

An Increase in the Antarctic Surface Mass Balance during the Past Three Centuries, Dampening Global Sea Level Rise

Abstract Antarctic surface mass balance (SMB) is a direct regulator of global sea level changes, but quantification of its long-term evolution at the ice sheet scale is challenging. Here, we combine for the first time a recently complied dataset of ice core records with the outputs of five different reanalysis products and two regional climate models to produce a reconciled 310-yr reconstruction of spatially and temporally complete SMB over the Antarctic Ice Sheet (AIS). Despite greatly variable signs and magnitudes of reconstructed SMB trends in the different regions, a significant positive trend (3.6 ± 0.8 Gt yr−1 decade−1) is observed for SMB over the entire AIS during the past 300 years, with a larger increase rate since 1801. The increased SMB cumulatively dampened global sea level rise by ∼14 mm between 1901 and 2010, which did not account for the ice sheet dynamical imbalance. The first and second modes of empirical orthogonal function analysis (EOF1 and EOF2) capture 38.0% and 24.6% of the total variability in reconstructed SMB, respectively. EOF1 consists of an east–west dipole of SMB changes over West Antarctica, primarily driven by Southern Annular Mode (SAM) variability. EOF2 represents a strong signal over the whole Antarctic Peninsula and coastal West Antarctica, which is not associated with SAM, but with El Niño–Southern Oscillation (ENSO) at the decadal scale.

Higher Antarctic ice sheet accumulation and surface melt rates revealed at 2 km resolution

Antarctic ice sheet (AIS) mass loss is predominantly driven by increased solid ice discharge, but its variability is governed by surface processes. Snowfall fluctuations control the surface mass balance (SMB) of the grounded AIS, while meltwater ponding can trigger ice shelf collapse potentially accelerating discharge. Surface processes are essential to quantify AIS mass change, but remain poorly represented in climate models typically running at 25-100 km resolution. Here we present SMB and surface melt products statistically downscaled to 2 km resolution for the contemporary climate (1979-2021) and low, moderate and high-end warming scenarios until 2100. We show that statistical downscaling modestly enhances contemporary SMB (3%), which is sufficient to reconcile modelled and satellite mass change. Furthermore, melt strongly increases (46%), notably near the grounding line, in better agreement with in-situ and satellite records. The melt increase persists by 2100 in all warming scenarios, revealing higher surface melt rates than previously estimated.

Local cooling and drying induced by Himalayan glaciers under global warming

Understanding the response of Himalayan glaciers to global warming is vital because of their role as a water source for the Asian subcontinent. However, great uncertainties still exist on the climate drivers of past and present glacier changes across scales. Here, we analyse continuous hourly climate station data from a glacierized elevation (Pyramid station, Mount Everest) since 1994 together with other ground observations and climate reanalysis. We show that a decrease in maximum air temperature and precipitation occurred during the last three decades at Pyramid in response to global warming. Reanalysis data suggest a broader occurrence of this effect in the glacierized areas of the Himalaya. We hypothesize that the counterintuitive cooling is caused by enhanced sensible heat exchange and the associated increase in glacier katabatic wind, which draws cool air downward from higher elevations. The stronger katabatic winds have also lowered the elevation of local wind convergence, thereby diminishing precipitation in glacial areas and negatively affecting glacier mass balance. This local cooling may have partially preserved glaciers from melting and could help protect the periglacial environment.

A computationally efficient statistically downscaled 100 m resolution Greenland product from the regional climate model MAR

Abstract. The Greenland Ice Sheet (GrIS) has been contributing directly to sea level rise, and this contribution is projected to accelerate over the next decades. A crucial tool for studying the evolution of surface mass loss (e.g., surface mass balance, SMB) consists of regional climate models (RCMs), which can provide current estimates and future projections of sea level rise associated with such losses. However, one of the main limitations of RCMs is the relatively coarse horizontal spatial resolution at which outputs are currently generated. Here, we report results concerning the statistical downscaling of the SMB modeled by the Modèle Atmosphérique Régional (MAR) RCM from the original spatial resolution of 6 km to 100 m building on the relationship between elevation and mass losses in Greenland. To this goal, we developed a geospatial framework that allows the parallelization of the downscaling process, a crucial aspect to increase the computational efficiency of the algorithm. Using the results obtained in the case of the SMB, surface and air temperature are assessed through the comparison of the modeled outputs with in situ and satellite measurement. The downscaled products show a considerable improvement in the case of the downscaled product with respect to the original coarse output, with the coefficient of determination (R2) increasing from 0.868 for the original MAR output to 0.935 for the SMB downscaled product. Moreover, the value of the slope and intercept of the linear regression fitting modeled and measured SMB values shifts from 0.865 for the original MAR to 1.015 for the downscaled product in the case of the slope and from the value −235 mm w.e. yr−1 (original) to −57 mm w.e. yr−1 (downscaled) in the case of the intercept, considerably improving upon results previously published in the literature.