Glacier Surface Mass Balance Research Articles

Strong contrast in mass and energy balance between a coastal mountain glacier and the Greenland ice sheet

ABSTRACTWe show a strong difference in surface mass and energy balance of a mountain glacier and two sites on the ice sheet at 64°N in West Greenland using stake and automated weather station observations. Net surface mass balance is on average 2.2 m w.e. less negative at the coast compared with the ice sheet in the same elevation. We find a larger energy turnover at the ice sheet margin on Qamanarssup Sermia than measured on the coastal mountain glacier Qassigiannguit with both energy input and output being of larger absolute value. More cloudiness and a thicker snow cover at the relatively humid coastal glacier result in smaller gains in net-shortwave radiation and smaller losses in net-longwave radiation and a less negative mass balance. Lower wind speeds at the coastal glacier result in weaker turbulent heat exchange between atmosphere and ice surface. On annual average, 17 W m−2more energy is available for melt at the ice-sheet margin compared with the coastal glacier in the same elevation.

Surface mass balance of Davies Dome and Whisky Glacier on James Ross Island, north-eastern Antarctic Peninsula, based on different volume-mass conversion approaches

This study presents surface mass balance of two small glaciers on James Ross Island calculated using constant and zonally-variable conversion factors. The density of 500 and 900 kg·m–3 adopted for snow in the accumulation area and ice in the ablation area, respectively, provides lower mass balance values that better fit to the glaciological records from glaciers on Vega Island and South Shetland Islands. The difference between the cumulative surface mass balance values based on constant (1.23 ± 0.44 m w.e.) and zonally-variable density (0.57 ± 0.67 m w.e.) is higher for Whisky Glacier where a total mass gain was observed over the period 2009–2015. The cumulative surface mass balance values are 0.46 ± 0.36 and 0.11 ± 0.37 m w.e. for Davies Dome, which experienced lower mass gain over the same period. The conversion approach does not affect much the spatial distribution of surface mass balance on glaciers, equilibrium line altitude and accumulation-area ratio. The pattern of the surface mass balance is almost identical in the ablation zone and very similar in the accumulation zone, where the constant conversion factor yields higher surface mass balance values. The equilibrium line altitude and accumulation-area ratio determined for the investigated glaciers differ by less than 2m and 0.01, respectively. The annual changes of equilibrium line altitude and the mean values determined over the period 2009–2015 for Whisky Glacier (311 ± 16 m a.s.l.) and Davies Dome (393 ± 18 m a.s.l.) coincide with the values reported from Bahía del Diablo Glacier on Vega Island but differ from the glaciological records on South Shetland Islands.

Quantification of annual glacier surface mass balance for the Chhota Shigri Glacier, Western Himalayas, India using an Equilibrium-Line Altitude (ELA) based approach

ABSTRACTIn line with the increasing scientific interest on the Himalayan glaciers, this study focuses on estimating a long-term annual surface mass balance time series of the Chhota Shigri glacier, a ‘benchmark’ glacier in the western Himalayas. The approach used here is based on the fact that the annual glacier-wide surface mass balance can be deduced from the equilibrium-line altitude (ELA). Depending on the distribution and availability of multiple cloud free remotely sensed images during ablation period, a multi-temporal approach has been used to estimate ELA. When compared with field-based ELA, the results indicate that the multi-temporal approach resulted in better estimates of ELA than the conventional single image approach. Likewise, the annual surface mass balances quantified from this study closely match with field estimates over the common period (2003–2014) and even better than some estimates from earlier studies based on other proxies (meteorological data or glacier surface albedo). A sensitivity analysis shows that the annual surface mass balance quantified from the ELA-based approach is not very sensitive to changes in the mass balance gradient and average mass balance. Hence, the approach has been further applied to reconstruct the long-term annual surface mass balance series of the Chhota Shigri Glacier over the period 1989–2017. Our results show a good agreement between the reconstructed surface mass balance and estimates of other long-term studies. Therefore, this study indicates the great potential for this approach for quantifying the annual surface mass balance for glaciers with no ground data lying in same climatic zone.

Multi-Criteria Evaluation of Snowpack Simulations in Complex Alpine Terrain Using Satellite and In Situ Observations

This work presents an extensive evaluation of the Crocus snowpack model over a rugged and highly glacierized mountain catchment (Arve valley, Western Alps, France) from 1989 to 2015. The simulations were compared and evaluated using in-situ point snow depth measurements, in-situ seasonal and annual glacier surface mass balance, snow covered area evolution based on optical satellite imagery at 250 m resolution (MODIS sensor), and the annual equilibrium-line altitude of glaciers, derived from satellite images (Landsat, SPOT, and ASTER). The snowpack simulations were obtained using the Crocus snowpack model driven by the same, originally semi-distributed, meteorological forcing (SAFRAN) reanalysis using the native semi-distributed configuration, but also a fully distributed configuration. The semi-distributed approach addresses land surface simulations for discrete topographic classes characterized by elevation range, aspect, and slope. The distributed approach operates on a 250-m grid, enabling inclusion of terrain shadowing effects, based on the same original meteorological dataset. Despite the fact that the two simulations use the same snowpack model, being potentially subjected to same potential deviation from the parametrization of certain physical processes, the results showed that both approaches accurately reproduced the snowpack distribution over the study period. Slightly (although statistically significantly) better results were obtained by using the distributed approach. The evaluation of the snow cover area with MODIS sensor has shown, on average, a reduction of the Root Mean Squared Error (RMSE) from 15.2% with the semi-distributed approach to 12.6% with the distributed one. Similarly, surface glacier mass balance RMSE decreased from 1.475 m of water equivalent (W.E.) for the semi-distributed simulation to 1.375 m W.E. for the distribution. The improvement, observed with a much higher computational time, does not justify the recommendation of this approach for all applications; however, for simulations that require a precise representation of snowpack distribution, the distributed approach is suggested.

Ice flow modelling to constrain the surface mass balance and ice discharge of San Rafael Glacier, Northern Patagonia Icefield

ABSTRACTWe simulate the ice dynamics of the San Rafael Glacier (SRG) in the Northern Patagonia Icefield (46.7°S, 73.5°W), using glacier geometry obtained by airborne gravity measurements. The full-Stokes ice flow model (Elmer/Ice) is initialized using an inverse method to infer the basal friction coefficient from a satellite-derived surface velocity mosaic. The high surface velocities (7.6 km a−1) near the glacier front are explained by low basal shear stresses (<25 kPa). The modelling results suggest that 98% of the surface velocities are due to basal sliding in the fast-flowing glacier tongue (>1 km a−1). We force the model using different surface mass-balance scenarios taken or adapted from previous studies and geodetic elevation changes between 2000 and 2012. Our results suggest that previous estimates of average surface mass balance over the entire glacier (Ḃ) were likely too high, mainly due to an overestimation in the accumulation area. We propose that most of SRG imbalance is due to the large ice discharge (−0.83 ± 0.08 Gt a−1) and a slightly positiveḂ(0.08 ± 0.06 Gt a−1). The committed mass-loss estimate over the next century is −0.34 ± 0.03 Gt a−1. This study demonstrates that surface mass-balance estimates and glacier wastage projections can be improved using a physically based ice flow model.

Multi-decadal mass balance series of three Kyrgyz glaciers inferred from modelling constrained with repeated snow line observations

Abstract. Glacier surface mass balance observations in the Tien Shan and Pamir are relatively sparse and often discontinuous. Nevertheless, glaciers are one of the most important components of the high-mountain cryosphere in the region as they strongly influence water availability in the arid, continental and intensely populated downstream areas. This study provides reliable and continuous surface mass balance series for selected glaciers located in the Tien Shan and Pamir-Alay. By cross-validating the results of three independent methods, we reconstructed the mass balance of the three benchmark glaciers, Abramov, Golubin and Glacier no. 354 for the past 2 decades. By applying different approaches, it was possible to compensate for the limitations and shortcomings of each individual method. This study proposes the use of transient snow line observations throughout the melt season obtained from satellite optical imagery and terrestrial automatic cameras. By combining modelling with remotely acquired information on summer snow depletion, it was possible to infer glacier mass changes for unmeasured years. The model is initialized with daily temperature and precipitation data collected at automatic weather stations in the vicinity of the glacier or with adjusted data from climate reanalysis products. Multi-annual mass changes based on high-resolution digital elevation models and in situ glaciological surveys were used to validate the results for the investigated glaciers. Substantial surface mass loss was confirmed for the three studied glaciers by all three methods, ranging from −0.30 ± 0.19 to −0.41 ± 0.33 m w.e. yr−1 over the 2004–2016 period. Our results indicate that integration of snow line observations into mass balance modelling significantly narrows the uncertainty ranges of the estimates. Hence, this highlights the potential of the methodology for application to unmonitored glaciers at larger scales for which no direct measurements are available.

Surface mass balance of small glaciers on James Ross Island, north-eastern Antarctic Peninsula, during 2009–2015

ABSTRACTTwo small glaciers on James Ross Island, the north-eastern Antarctic Peninsula, experienced surface mass gain between 2009 and 2015 as revealed by field measurements. A positive cumulative surface mass balance of 0.57 ± 0.67 and 0.11 ± 0.37 m w.e. was observed during the 2009–2015 period on Whisky Glacier and Davies Dome, respectively. The results indicate a change from surface mass loss that prevailed in the region during the first decade of the 21st century to predominantly positive surface mass balance after 2009/10. The spatial pattern of annual surface mass-balance distribution implies snow redistribution by wind on both glaciers. The mean equilibrium line altitudes for Whisky Glacier (311 ± 16 m a.s.l.) and Davies Dome (393 ± 18 m a.s.l.) are in accordance with the regional data indicating 200–300 m higher equilibrium line on James Ross and Vega Islands compared with the South Shetland Islands. The mean accumulation-area ratio of 0.68 ± 0.09 and 0.44 ± 0.09 determined for Whisky Glacier and Davies Dome, respectively, is similar to the value reported for Vega Island and within the range of typical values for high-latitude glaciers.

Relative performance of empirical and physical models in assessing the seasonal and annual glacier surface mass balance of Saint-Sorlin Glacier (French Alps)

Abstract. This study focuses on simulations of the seasonal and annual surface mass balance (SMB) of Saint-Sorlin Glacier (French Alps) for the period 1996–2015 using the detailed SURFEX/ISBA-Crocus snowpack model. The model is forced by SAFRAN meteorological reanalysis data, adjusted with automatic weather station (AWS) measurements to ensure that simulations of all the energy balance components, in particular turbulent fluxes, are accurately represented with respect to the measured energy balance. Results indicate good model performance for the simulation of summer SMB when using meteorological forcing adjusted with in situ measurements. Model performance however strongly decreases without in situ meteorological measurements. The sensitivity of the model to meteorological forcing indicates a strong sensitivity to wind speed, higher than the sensitivity to ice albedo. Compared to an empirical approach, the model exhibited better performance for simulations of snow and firn melting in the accumulation area and similar performance in the ablation area when forced with meteorological data adjusted with nearby AWS measurements. When such measurements were not available close to the glacier, the empirical model performed better. Our results suggest that simulations of the evolution of future mass balance using an energy balance model require very accurate meteorological data. Given the uncertainties in the temporal evolution of the relevant meteorological variables and glacier surface properties in the future, empirical approaches based on temperature and precipitation could be more appropriate for simulations of glaciers in the future.

Monitoring glacier albedo as a proxy to derive summer and annual surface mass balances from optical remote-sensing data

Abstract. Less than 0.25 % of the 250 000 glaciers inventoried in the Randolph Glacier Inventory (RGI V.5) are currently monitored with in situ measurements of surface mass balance. Increasing this archive is very challenging, especially using time-consuming methods based on in situ measurements, and complementary methods are required to quantify the surface mass balance of unmonitored glaciers. The current study relies on the so-called albedo method, based on the analysis of albedo maps retrieved from optical satellite imagery acquired since 2000 by the MODIS sensor, on board the TERRA satellite. Recent studies revealed substantial relationships between summer minimum glacier-wide surface albedo and annual surface mass balance, because this minimum surface albedo is directly related to the accumulation–area ratio and the equilibrium-line altitude. On the basis of 30 glaciers located in the French Alps where annual surface mass balance data are available, our study conducted on the period 2000–2015 confirms the robustness and reliability of the relationship between the summer minimum surface albedo and the annual surface mass balance. For the ablation season, the integrated summer surface albedo is significantly correlated with the summer surface mass balance of the six glaciers seasonally monitored. These results are promising to monitor both annual and summer glacier-wide surface mass balances of individual glaciers at a regional scale using optical satellite images. A sensitivity study on the computed cloud masks revealed a high confidence in the retrieved albedo maps, restricting the number of omission errors. Albedo retrieval artifacts have been detected for topographically incised glaciers, highlighting limitations in the shadow correction algorithm, although inter-annual comparisons are not affected by systematic errors.

Model-based calculations of surface mass balance of mountain glaciers for the purpose of water consumption planning: focus on Djankuat Glacier (Central Caucasus)

Mountain glaciers act as regulators of run-off in the summer period, which is very crucial for economy especially in dynamically developing regions with rapidly growing population, such as Central Asia or the Northern Caucasus in Russia. In overall, glaciers stabilize water consumption in comparatively arid areas and provide conditions for sustainable development of the economy in mountainous regions and in the surrounding territories. A proper prediction of the glacial run-off is required to elaborate strategies of the regional development. This goal can be achieved by implementation of mathematical modeling methods into planning methodologies. In the paper, we consider one of the first steps in glacier dynamical modeling – surface mass balance simulation. We focus on the Djankuat Glacier in the Central Caucasus, where regular observations have been conducted during the last fifty years providing an exceptional opportunity to calibrate and to validate a mathematical model.

Relationship between glacier melting and atmospheric circulation in the southeast Siberia

The interaction between climate and cryosphere is a key issue in recent years. Changes in surface mass balance of mountain glaciers closely correspond to differential changes in atmospheric circulation. Mountain glaciers in southeast Siberia located on East Sayan, Baikalsky and Kodar ridges have been continuously shrinking since the end of the Little Ice Age. In this study we used daily synoptic weather maps (Irkutsk Center of Hydrometeorology and Environmental Monitoring), 500 hPa, 700 hPa and 850 hPa geopotential height and air temperature data of NCEP/NCAR reanalysis to assess relationships between atmospheric circulation patterns and the sum of positive temperature (SPT), a predictor of summer ice/snow ablation. Results show that increased SPT (ablation) is generally associated with anticyclones and anticyclonic pressure fields (with cloudless weather conditions) and warm atmospheric fronts. Decreased SPT (ablation) is strongly correlated with cyclones and cyclonic type pressure fields, cold atmospheric fronts and air advections. Significant correlations have been found between ablation and cyclonic/anticyclonic activity. Revealed decreasing trends in the SPT in three glaciarized ridges at the beginning of the 21st century led to changes of air temperature and snow/ice melt climates.

The response of surface mass and energy balance of a continental glacier to climate variability, western Qilian Mountains, China

To understand how a continental glacier responds to climate change, it is imperative to quantify the surface energy fluxes and identify factors controlling glacier mass balance using surface energy balance (SEB) model. Light absorbing impurities (LAIs) at the glacial surface can greatly decrease surface albedo and increase glacial melt. An automatic weather station was set up and generated a unique 6-year meteorological dataset for the ablation zone of Laohugou Glacier No. 12. Based on these data, the surface energy budget was calculated and an experiment on the glacial melt process was carried out. The effect of reduced albedo on glacial melting was analyzed. Owing to continuous accumulation of LAIs, the ablation zone had been darkening since 2010. The mean value of surface albedo in melt period (June through September) dropped from 0.52 to 0.43, and the minimum of daily mean value was as small as 0.1. From the records of 2010–2015, keeping the clean ice albedo fixed in the range of 0.3–0.4, LAIs caused an increase of +7.1 to +16 W m−2 of net shortwave radiation and an removal of 1101–2663 mm water equivalent. Calculation with the SEB model showed equivalent increases in glacial melt were obtained by increasing air temperature by 1.3 and 3.2 K, respectively.

The importance of accurate glacier albedo for estimates of surface mass balance on Vatnajökull: evaluating the surface energy budget in a regional climate model with automatic weather station observations

Abstract. A simulation of the surface climate of Vatnajökull ice cap, Iceland, carried out with the regional climate model HIRHAM5 for the period 1980–2014, is used to estimate the evolution of the glacier surface mass balance (SMB). This simulation uses a new snow albedo parameterization that allows albedo to exponentially decay with time and is surface temperature dependent. The albedo scheme utilizes a new background map of the ice albedo created from observed MODIS data. The simulation is evaluated against observed daily values of weather parameters from five automatic weather stations (AWSs) from the period 2001–2014, as well as in situ SMB measurements from the period 1995–2014. The model agrees well with observations at the AWS sites, albeit with a general underestimation of the net radiation. This is due to an underestimation of the incoming radiation and a general overestimation of the albedo. The average modelled albedo is overestimated in the ablation zone, which we attribute to an overestimation of the thickness of the snow layer and not taking the surface darkening from dirt and volcanic ash deposition during dust storms and volcanic eruptions into account. A comparison with the specific summer, winter, and net mass balance for the whole of Vatnajökull (1995–2014) shows a good overall fit during the summer, with a small mass balance underestimation of 0.04 m w.e. on average, whereas the winter mass balance is overestimated by on average 0.5 m w.e. due to too large precipitation at the highest areas of the ice cap. A simple correction of the accumulation at the highest points of the glacier reduces this to 0.15 m w.e. Here, we use HIRHAM5 to simulate the evolution of the SMB of Vatnajökull for the period 1981–2014 and show that the model provides a reasonable representation of the SMB for this period. However, a major source of uncertainty in the representation of the SMB is the representation of the albedo, and processes currently not accounted for in RCMs, such as dust storms, are an important source of uncertainty in estimates of snow melt rate.

Annual and Seasonal Glacier-Wide Surface Mass Balance Quantified from Changes in Glacier Surface State: A Review on Existing Methods Using Optical Satellite Imagery

Glaciers are one of the terrestrial essential climate variables (ECVs) as they respond very sensitively to climate change. A key driver of their response is the glacier surface mass balance that is typically derived from field measurements. It deserves to be quantified over long time scales to better understand the accumulation and ablation processes at the glacier surface and their relationships with inter-annual changes in meteorological conditions and long-term climate changes. Glaciers with in situ monitoring of surface mass balance are scarce at the global scale, and satellite remote sensing provides a powerful tool to increase the number of monitored glaciers. In this study, we present a review of three optical remote sensing methods developed to quantify seasonal and annual glacier surface mass balances. These methodologies rely on the multitemporal monitoring of the end-of-summer snow line for the equilibrium-line altitude (ELA) method, the annual cycle of glacier surface albedo for the albedo method and the mapping of the regional snow cover at the seasonal scale for the snow-map method. Together with a presentation of each method, an application is illustrated. The ELA method shows promising results to quantify annual surface mass balance and to reconstruct multi-decadal time series. The other two methods currently need a calibration on the basis of existing in situ data; however, a generalization of these methods (without calibration) could be achieved. The two latter methods show satisfying results at the annual and seasonal scales, particularly for the summer surface mass balance in the case of the albedo method and for the winter surface mass balance in the case of the snow-map method. The limits of each method (e.g., cloud coverage, debris-covered glaciers, monsoon-regime and cold glaciers), their complementarities and the future challenges (e.g., automating of the satellite images processing, generalization of the methods needing calibration) are also discussed.

Surface mass balance on Glacier No. 31 in the Suntar–Khayata Range, eastern Siberia, from 1951 to 2014

This study presents a 64-year (1951–2014) reconstruction of the surface mass balance of Glacier No. 31, located in the Suntar-Khayata Range of the eastern Siberia, where the ablation zone is characterized by the extensive dark ice surface. We use a temperature index-based glacier mass-balance model, which computes all major components of glacier mass budget and is forced by daily air temperature and precipitation from a nearby meteorological station. The glacier shows a mean annual mass balance of–0.35 m w.e.a–1 during the past 64 years, with an acceleration of–0.50 m w.e. a–1 during the recent years. A cumulative mass loss of the glacier is ~22.3 m w.e. over the study period, about 56% of which is observed during 1991–2014. In addition to the contribution of temperature rise and precipitation decrease to recent mass loss of the glacier, an experimental analysis, in which the clean and dark ice surfaces are respectively assumed to cover the entire ablation zone, indicates that dark ice surface, caused by insoluble impurities consisting of mineral dusts, cryoconite granules, and ice algae, plays a crucial role in the changing mass balance through enhancing melt rates in the ablation zone of the glacier.

Glacier surface temperatures in the Canadian High Arctic, 2000–15

ABSTRACTCanada's Queen Elizabeth Islands (QEI) contain ~14% of the world's glacier and ice-cap area. Sparse in-situ measurements indicate that interannual variability in glacier surface mass balance in this region is driven primarily by variations in summer melt, and that the annual surface mass balance of four index glaciers has become increasingly negative since 2007. Here, we use a 16-a record of satellite-derived mean summer (June–August) land surface temperatures (LST) from NASA's Moderate Resolution Imaging Spectroradiometer to investigate large-scale spatial and temporal variability in the duration and intensity of summer melt across glaciated surfaces in the QEI from 2000 to 2015. During this period, QEI mean summer glacier surface temperatures increased at an average rate of 0.06 ± 0.04°C a−1, for a total of nearly 1°C. Most of this increase occurred between 2005 and 2012, when mean summer near-surface (2 m) and upper-air (700 hPa) temperatures were 1.0–1.2°C higher than the 1948–2015 mean. There is a strong correlation between the glacier LST and 700 hPa air temperature records (r> 0.8). The period 2005–12, when mean summer LSTs were anomalously high, was likely the warmest period in the region since at least 1948.

A mass-flux perspective of the tidewater glacier cycle

AbstractI explore the tidewater glacier cycle with a 1-D, depth- and width-integrated flow model that includes a mass-flux calving parameterization. The parameterization is developed from mass continuity arguments and relates the calving rate to the terminus velocity and the terminus balance velocity. The model demonstrates variable sensitivity to climate. From an advanced, stable configuration, a small warming of the climate triggers a rapid retreat that causes large-scale drawdown and is enhanced by positive glacier-dynamic feedbacks. Eventually, the terminus retreats out of deep water and the terminus velocity decreases, resulting in reduced drawdown and the potential for restabilization. Terminus readvance can be initiated by cooling the climate. Terminus advance into deep water is difficult to sustain, however, due to negative feedbacks between glacier dynamics and surface mass balance. Despite uncertainty in the precise form of the parameterization, the model provides a simple explanation of the tidewater glacier cycle and can be used to evaluate the response of tidewater glaciers to climate variability. It also highlights the importance of improving parameterizations of calving rates and of incorporating sediment dynamics into tidewater glacier models.

Changes in ice dynamics, elevation and mass discharge of Dinsmoor–Bombardier–Edgeworth glacier system, Antarctic Peninsula

The northern Antarctic Peninsula is one of the fastest changing regions on Earth. The disintegration of the Larsen-A Ice Shelf in 1995 caused tributary glaciers to adjust by speeding up, surface lowering, and overall increased ice-mass discharge. In this study, we investigate the temporal variation of these changes at the Dinsmoor–Bombardier–Edgeworth glacier system by analyzing dense time series from various spaceborne and airborne Earth observation missions. Precollapse ice shelf conditions and subsequent adjustments through 2014 were covered. Our results show a response of the glacier system some months after the breakup, reaching maximum surface velocities at the glacier front of up to 8.8 m/d in 1999 and a subsequent decrease to ∼1.5 m/d in 2014. Using a dense time series of interferometrically derived TanDEM-X digital elevation models and photogrammetric data, an exponential function was fitted for the decrease in surface elevation. Elevation changes in areas below 1000 m a.s.l. amounted to at least 130±15 m between 1995 and 2014, with change rates of ∼3.15 m/a between 2003 and 2008. Current change rates (2010–2014) are in the range of 1.7 m/a. Mass imbalances were computed with different scenarios of boundary conditions. The most plausible results amount to −40.7±3.9 Gt. The contribution to sea level rise was estimated to be 18.8±1.8 Gt, corresponding to a 0.052±0.005 mm sea level equivalent, for the period 1995–2014. Our analysis and scenario considerations revealed that major uncertainties still exist due to insufficiently accurate ice-thickness information. The second largest uncertainty in the computations was the glacier surface mass balance, which is still poorly known. Our time series analysis facilitates an improved comparison with GRACE data and as input to modeling of glacio-isostatic uplift in this region. The study contributed to a better understanding of how glacier systems adjust to ice shelf disintegration.

Northern Hemisphere Glacier and Ice Cap Surface Mass Balance and Contribution to Sea Level Rise

Abstract Mass changes and mass contribution to sea level rise from glaciers and ice caps (GIC) are key components of the earth’s changing sea level. GIC surface mass balance (SMB) magnitudes and individual and regional mean conditions and trends (1979–2009) were simulated for all GIC having areas greater or equal to 0.5 km2 in the Northern Hemisphere north of 25°N latitude (excluding the Greenland Ice Sheet). Recent datasets, including the Randolph Glacier Inventory (RGI; v. 2.0), the NOAA Global Land One-km Base Elevation Project (GLOBE), and the NASA Modern-Era Retrospective Analysis for Research and Applications (MERRA) products, together with recent SnowModel developments, allowed relatively high-resolution (1-km horizontal grid; 3-h time step) simulations of GIC surface air temperature, precipitation, sublimation, evaporation, surface runoff, and SMB. Simulated SMB outputs were calibrated against 1422 direct glaciological annual SMB observations of 78 GIC. The overall GIC mean annual and mean summer air temperature, runoff, and SMB loss increased during the simulation period. The cumulative GIC SMB was negative for all regions. The SMB contribution to sea level rise was largest from Alaska and smallest from the Caucasus. On average, the contribution to sea level rise was 0.51 ± 0.16 mm sea level equivalent (SLE) yr−1 for 1979–2009 and ~40% higher 0.71 ± 0.15 mm SLE yr−1 for the last decade, 1999–2009.

All‐sky radiation over a glacier surface in the Southern Alps of New Zealand: characterizing cloud effects on incoming shortwave, longwave and net radiation

ABSTRACTClouds are important features of many high‐altitude and glaciated areas, yet detecting their presence and specifying their effects on incoming shortwave (SW↓), longwave (LW↓) and net all‐wave radiation (Rnet) remains challenging in these environments. These limitations hamper efforts to understand atmospheric controls on glacier surface mass balance (SMB) in the Southern Alps of New Zealand, as both cloud and airmass forcing accompanies key synoptic controls on SMB. Multi‐year datasets of four‐component broadband radiation from two sites at Brewster Glacier, Southern Alps of New Zealand, are used here to develop cloud metrics to account for the effects of clouds on SW↓, LW↓ and Rnet. On average 23% of top‐of‐atmosphere shortwave radiation (SWTOA) is attenuated by the clear‐sky atmosphere, while clouds attenuate a further 31%, resulting in <50% of SWTOA being received at the surface. The transmission of shortwave radiation by clouds (trc) during overcast conditions is found to vary with season and airmass characteristics. A simple parameterization is developed to account for lower trc observed during periods of higher water vapour pressure. Cloud metrics derived at the site show overcast conditions are frequent (45% of period) and strongly dependent on wind direction, highlighting the dominant role of orography in cloud formation and enhancement in the Southern Alps. The effect of clouds on Rnet exhibits a distinct seasonal variation; during summer when albedo and trc are lower, clouds decrease Rnet by 20–40 W m−2, while during autumn, winter and spring, clouds enhance Rnet by approximately 20 W m−2. Idealized modelling shows that these patterns are strongly dependent on albedo and extend across the elevation range of glaciers in the Southern Alps. Thus, overcast conditions appear to aid the extension of ablation into spring and autumn by increasing the energy available for snow and ice melt.