Surface Mass Balance Components Research Articles

The added value of high resolution in estimating the surface mass balance in southern Greenland

Abstract. The polar version of the regional climate model RACMO2, version 2.3p1, is used to study the effect of model resolution on the simulated climate and surface mass balance (SMB) of south Greenland for the current climate (2007–2014). The model data at resolutions of 60, 20, 6.6, and 2.2 km are intercompared and compared to SMB observations using three different data refinement methods: nearest neighbour, bilinear interpolation, and a statistical downscaling method utilising the local dependency of fields on elevation. Furthermore, it is estimated how the errors induced by model resolution compare to errors induced by the model physics and initialisation. The results affirm earlier studies that SMB components which are tightly linked to elevation, like runoff, can be refined successfully, as soon as the ablation zone is reasonably well resolved in the source dataset. Precipitation fields are also highly elevation dependent, but precipitation has no systematic correlation with elevation, which inhibits statistical downscaling to work well. If refined component-wise, 20 km resolution model simulations can reproduce the SMB ablation observations almost as well as the finer-resolution model simulations. Nonetheless, statistical downscaling and regional climate modelling are complementary; the best results are obtained when high-resolution RACMO2 data are statistically refined. Model estimates in the accumulation zone do not benefit from statistical downscaling; hence, a resolution of about 20 km is sufficient to resolve the majority of the accumulation zone of the Greenland Ice Sheet with respect to the limited measurements we have. Furthermore, we demonstrate that using RACMO2, a hydrostatic model, at 2.2 km resolution led to invalid results as topographic and synoptic vertical winds exceed 10 m s−1, which violates the hydrostatic model assumptions. Finally, additional tests show that model resolution is as important as properly resolving spatial albedo patterns, correctly initialising the firn column, and uncertainties in the modelled precipitation and turbulent exchange.

Brief communication: CESM2 climate forcing (1950–2014) yields realistic Greenland ice sheet surface mass balance

Abstract. We present a reconstruction of historical (1950–2014) surface mass balance (SMB) of the Greenland ice sheet (GrIS) using a high-resolution regional climate model (RACMO2; ∼11 km) to dynamically downscale the climate of the Community Earth System Model version 2 (CESM2; ∼111 km). After further statistical downscaling to 1 km spatial resolution, evaluation using in situ SMB measurements and remotely sensed GrIS mass change shows good agreement. Comparison with an ensemble of previously conducted RACMO2 simulations forced by climate reanalysis demonstrates that the current product realistically represents the long-term average and variability of individual SMB components and captures the recent increase in meltwater runoff that accelerated GrIS mass loss. This means that, for the first time, climate forcing from an Earth system model (CESM2), which assimilates no observations, can be used without additional corrections to reconstruct the historical GrIS SMB and its recent decline that initiated mass loss in the 1990s. This paves the way for attribution studies of future GrIS mass loss projections and contribution to sea level rise.

Present‐Day Greenland Ice Sheet Climate and Surface Mass Balance in CESM2

The response of the Greenland Ice Sheet (GrIS) to a warmer climate is uncertain on long time scales. Climate models, such as those participating in the Coupled Model Intercomparison Project phase 6 (CMIP6), are used to assess this uncertainty. The Community Earth System Model version 2.1 (CESM2) is a CMIP6 model capable of running climate simulations with either one‐way coupling (fixed ice sheet geometry) or two‐way coupling (dynamic geometry) to the GrIS. The model features prognostic snow albedo, online downscaling using elevation classes, and a firn pack to refreeze percolating melt water. Here we evaluate the representation of the GrIS surface energy balance and surface mass balance in CESM2 at 1° resolution with fixed GrIS geometry. CESM2 agrees closely with ERA‐Interim reanalysis data for key controls on GrIS SMB: surface pressure, sea ice extent, 500 hPa geopotential height, wind speed, and 700 hPa air temperature. Cloudsat‐CALIPSO data show that supercooled liquid‐containing clouds are adequately represented, whereas comparisons to Moderate Resolution Imaging Spectroradiometer and CM SAF Cloud, Albedo, and Surface Radiation data set from Advanced Very High Resolution Radiometer data second edition data suggest that CESM2 underestimates surface albedo. The seasonal cycle and spatial patterns of surface energy balance and surface mass balance components in CESM2 agree well with regional climate model RACMO2.3p2, with GrIS‐integrated melt, refreezing, and runoff bracketed by RACMO2 counterparts at 11 and 1 km. Time series of melt, runoff, and SMB show a break point around 1990, similar to RACMO2. These results suggest that GrIS SMB is realistic in CESM2, which adds confidence to coupled ice sheet‐climate experiments that aim to assess the GrIS contribution to future sea level rise.

Importance of snow as component of surface mass balance of Arctic glacier (Hansbreen, southern Spitsbergen)

Importance of snow as component of surface mass balance of Arctic glacier (Hansbreen, southern Spitsbergen)

Applying the Energy- and Water Balance Model for Incorporation of the Cryospheric Component into a Climate Model. Part III. Modeling Mass Balance on the Surface of the Antarctic Ice Sheet

The Earth system models have been actively developed during recent decades. They integrate the blocks describing the dynamics of the atmosphere, the ocean and the cryosphere. The Antarctic ice sheet (AIS) is one of the most important components of the Earth’s cryosphere. Therefore, the incorporation of the AIS model into the Earth system model will significantly improve the accuracy of reconstructions of climate parameters for paleoclimate studies and for climate predictions. This paper presents a description of the energy balance block of the AIS model. This is a continuation of our recent work on integration of the cryosphere block to the climate model. The testing of the block (the EWBM-A buffer model) was carried out using the climate data generated by the INMCM model developed in the Institute of Numerical Mathematics of RAS. The results of the calculations of the surface mass balance components of the Antarctic ice sheet are presented for the 30-year pre-industrial period and are compared with the recent results obtained in the literature.

Link between the North Atlantic Oscillation and the surface mass balance components of the Greenland Ice Sheet under preindustrial and last interglacial climates: a study with a coupled global circulation model

Abstract. The relationship between the surface mass balance (SMB) components (accumulation and melting) of the Greenland Ice Sheet (GrIS) and the North Atlantic Oscillation (NAO) is examined from numerical simulations performed with a new atmospheric stretched grid configuration of the Centre National de Recherches Météorologiques Coupled Model (CNRM-CM) version 5.2 under three periods: preindustrial climate, a warm phase (early Eemian, 130 ka BP) and a cool phase (late Eemian, 115 ka BP) of the last interglacial. The horizontal grid of the atmospheric component of CNRM-CM5.2 is stretched from the tilted pole on Baffin Bay (72∘ N, 65∘ W) in order to obtain a higher spatial resolution on Greenland. The correlation between simulated SMB anomalies averaged over Greenland and the NAO index is weak in winter and significant in summer (about 0.6 for the three periods). In summer, spatial correlations between the NAO index and SMB components display different patterns from one period to another. These differences are analyzed in terms of the respective influence of the positive and negative phases of the NAO on accumulation and melting. Accumulation in south Greenland is significantly correlated with the positive (negative) phase of the NAO in a warm (cold) climate. Under preindustrial and 115 ka BP climates, melting along the margins is more correlated with the positive phase of the NAO than with its negative phase, whereas at 130 ka BP it is more correlated with the negative phase of the NAO in north and northeast Greenland.

The SUMup dataset: compiled measurements of surface mass balance components over ice sheets and sea ice with analysis over Greenland

Abstract. Increasing atmospheric temperatures over ice cover affect surface processes, including melt, snowfall, and snow density. Here, we present the Surface Mass Balance and Snow on Sea Ice Working Group (SUMup) dataset, a standardized dataset of Arctic and Antarctic observations of surface mass balance components. The July 2018 SUMup dataset consists of three subdatasets, snow/firn density (https://doi.org/10.18739/A2JH3D23R), at least near-annually resolved snow accumulation on land ice (https://doi.org/10.18739/A2DR2P790), and snow depth on sea ice (https://doi.org/10.18739/A2WS8HK6X), to monitor change and improve estimates of surface mass balance. The measurements in this dataset were compiled from field notes, papers, technical reports, and digital files. SUMup is a compiled, community-based dataset that can be and has been used to evaluate modeling efforts and remote sensing retrievals. Active submission of new or past measurements is encouraged. Analysis of the dataset shows that Greenland Ice Sheet density measurements in the top 1 m do not show a strong relationship with annual temperature. At Summit Station, Greenland, accumulation and surface density measurements vary seasonally with lower values during summer months. The SUMup dataset is a dynamic, living dataset that will be updated and expanded for community use as new measurements are taken and new processes are discovered and quantified.

Using remotely sensed data from AIRS to estimate the vapor flux on the Greenland ice sheet: Comparisons with observations and a regional climate model

Mass loss from the Greenland ice sheet (GrIS) in recent years has been dominated by runoff from surface melt. It is currently being studied extensively, while little interest has been given to the smallest component of surface mass balance (SMB): the vapor flux. Although poorly understood, it is not negligible and could potentially play a larger role in SMB in a warming climate where temperature, relative humidity, and precipitation changes remain uncertain. Here we present an innovative approach to estimate the vapor flux using the Atmospheric Infrared Sounder (AIRS) version 6 data and a modified vapor flux model (BMF13) over the GrIS between 2003 and 2014. One modification to the BMF13 model includes a new Multiangle Imaging SpectroRadiometer surface aerodynamic roughness product, which likely produces more accurate estimates of the drag coefficient on the ice sheet. When comparing AIRS data with GC-Net and Programme for Monitoring of the Greenland Ice Sheet automatic weather station observations of skin temperature, near-surface air temperature, and humidity, they agree within 2 K, 2.68 K, and 0.34 g kg−1. Largest differences occur in the ablation zone where there is significant subgrid heterogeneity. Overall, the average vapor flux from the GrIS between 2003 and 2014 was found to be 14.6 ± 3.6 Gt yr−1. No statistically significant trends were found during the data record. This data set is compared to the Regional Atmospheric Climate Model (RACMO2.3) vapor flux, and BMF13 produced smaller vapor fluxes in the summer (~0.05 Gt d−1) and slightly more deposition in the winter (~9.4 × 10−3 Gt d−1). Annually, differences between BMF13 and RACMO2.3 were only 30 ± 15%.

A daily, 1 km resolution data set of downscaled Greenland ice sheet surface mass balance (1958–2015)

Abstract. This study presents a data set of daily, 1 km resolution Greenland ice sheet (GrIS) surface mass balance (SMB) covering the period 1958–2015. Applying corrections for elevation, bare ice albedo and accumulation bias, the high-resolution product is statistically downscaled from the native daily output of the polar regional climate model RACMO2.3 at 11 km. The data set includes all individual SMB components projected to a down-sampled version of the Greenland Ice Mapping Project (GIMP) digital elevation model and ice mask. The 1 km mask better resolves narrow ablation zones, valley glaciers, fjords and disconnected ice caps. Relative to the 11 km product, the more detailed representation of isolated glaciated areas leads to increased precipitation over the southeastern GrIS. In addition, the downscaled product shows a significant increase in runoff owing to better resolved low-lying marginal glaciated regions. The combined corrections for elevation and bare ice albedo markedly improve model agreement with a newly compiled data set of ablation measurements.

The modelled surface mass balance of the Antarctic Peninsula at 5.5 km horizontal resolution

Abstract. This study presents a high-resolution (∼ 5.5 km) estimate of surface mass balance (SMB) over the period 1979–2014 for the Antarctic Peninsula (AP), generated by the regional atmospheric climate model RACMO2.3 and a firn densification model (FDM). RACMO2.3 is used to force the FDM, which calculates processes in the snowpack, such as meltwater percolation, refreezing and runoff. We evaluate model output with 132 in situ SMB observations and discharge rates from six glacier drainage basins, and find that the model realistically simulates the strong spatial variability in precipitation, but that significant biases remain as a result of the highly complex topography of the AP. It is also clear that the observations significantly underrepresent the high-accumulation regimes, complicating a full model evaluation. The SMB map reveals large accumulation gradients, with precipitation values above 3000 mm we yr−1 in the western AP (WAP) and below 500 mm we yr−1 in the eastern AP (EAP), not resolved by coarser data sets such as ERA-Interim. The average AP ice-sheet-integrated SMB, including ice shelves (an area of 4.1 × 105 km2), is estimated at 351 Gt yr−1 with an interannual variability of 58 Gt yr−1, which is dominated by precipitation (PR) (365 ± 57 Gt yr−1). The WAP (2.4 × 105 km2) SMB (276 ± 47 Gt yr−1), where PR is large (276 ± 47 Gt yr−1), dominates over the EAP (1.7 × 105 km2) SMB (75 ± 11 Gt yr−1) and PR (84 ± 11 Gt yr−1). Total sublimation is 11 ± 2 Gt yr−1 and meltwater runoff into the ocean is 4 ± 4 Gt yr−1. There are no significant trends in any of the modelled AP SMB components, except for snowmelt that shows a significant decrease over the last 36 years (−0.36 Gt yr−2).

Evaluation of the CMIP5 models in the aim of regional modelling of the Antarctic surface mass balance

Abstract. The surface mass balance (SMB) of the Antarctic Ice Sheet cannot be reliably deduced from global climate models (GCMs), both because their spatial resolution is insufficient and because their physics are not adapted for cold and snow-covered regions. By contrast, regional climate models (RCMs) adapted for polar regions can physically and dynamically downscale SMB components over the ice sheet using large-scale forcing at their boundaries. Polar-oriented RCMs require appropriate GCM fields for forcing because the response of the cryosphere to a warming climate is dependent on its initial state and is not linear with respect to temperature increase. In this context, we evaluate the current climate in 41 climate models from the Coupled Model Intercomparison Project Phase 5 (CMIP5) data set over Antarctica by focusing on forcing fields which may have the greatest impact on SMB components simulated by RCMs. Our inter-comparison includes six reanalyses, among which ERA-Interim reanalysis is chosen as a reference over 1979–2014. Model efficiency is assessed taking into account the multi-decadal variability of the fields over the 1850–1980 period. We show that fewer than 10 CMIP5 models show reasonable biases compared to ERA-Interim, among which ACCESS1-3 is the most pertinent choice for forcing RCMs over Antarctica, followed by ACCESS1-0, CESM1-BGC, CESM1-CAM5, NorESM1-M, CCSM4 and EC-EARTH. Finally, climate change over the Southern Ocean in CMIP5 is less sensitive to the global warming signal than it is to the present-day simulated sea-ice extent and to the feedback between sea-ice decrease and air temperature increase around Antarctica.

Coast-to-interior gradient in recent northwest Greenland precipitation trends (1952–2012)

The spatial and temporal variability of precipitation on the Greenland ice sheet is an essential component of surface mass balance, which has been declining in recent years with rising temperatures. We present an analysis of precipitation trends in northwest (NW) Greenland (1952–2012) using instrumental (coastal meteorological station) and proxy records (snow pits and ice cores) to characterize the precipitation gradient from the coast to the ice sheet interior. Snow-pit-derived precipitation near the coast (1950–2000) has increased (∼7% decade−1, p < 0.01) whereas there is no significant change observed in interior snow pits. This trend holds for 1981–2012, where calculated precipitation changes decrease in magnitude with increasing distance from the coast: 13% decade−1 (2.4 mm water equivalent (w.e.) decade−2) at coastal Thule air base (AB), 8.6% decade−1 (4.7 mm w.e. decade−2) at the 2Barrel ice core site 150 km from Thule AB, −5.2% decade−1 (1.7 mm w.e. decade−2) at Camp Century located 205 km from Thule AB, and 4.4% decade−1 (1.0 mm w.e. decade−2) at B26 located 500 km from Thule AB. In general, annually averaged precipitation and annually and seasonally averaged mean air temperatures observed at Thule AB follow trends observed in composite coastal Greenland time series, with both notably indicating winter as the fastest warming season in recent periods (1981–2012). Trends (1961–2012) in seasonal precipitation differ, specifically with NW Greenland summer precipitation increasing (∼0.6 mm w.e. decade−2) in contrast with decreasing summer precipitation in the coastal composite time series (3.8 mm w.e. decade−2). Differences in precipitation trends between NW Greenland and coastal composite Greenland underscore the heterogeneity in climate influences affecting precipitation. In particular, recent (1981–2012) changes in NW Greenland annual precipitation are likely a response to a weakening North Atlantic oscillation.

Cloud and precipitation properties from ground-based remote-sensing instruments in East Antarctica

Abstract. A new comprehensive cloud–precipitation–meteorological observatory has been established at Princess Elisabeth base, located in the escarpment zone of Dronning Maud Land (DML), East Antarctica. The observatory consists of a set of ground-based remote-sensing instruments (ceilometer, infrared pyrometer and vertically profiling precipitation radar) combined with automatic weather station measurements of near-surface meteorology, radiative fluxes, and snow height. In this paper, the observatory is presented and the potential for studying the evolution of clouds and precipitating systems is illustrated by case studies. It is shown that the synergetic use of the set of instruments allows for distinguishing ice, liquid-containing clouds and precipitating clouds, including some information on their vertical extent. In addition, wind-driven blowing snow events can be distinguished from deeper precipitating systems. Cloud properties largely affect the surface radiative fluxes, with liquid-containing clouds dominating the radiative impact. A statistical analysis of all measurements (in total 14 months mainly during summer–beginning of winter) indicates that these liquid-containing clouds occur during as much as 20% of the cloudy periods. The cloud occurrence shows a strong bimodal distribution with clear-sky conditions 51% of the time and complete overcast conditions 35% of the time. Snowfall occurred during 17% of the cloudy periods with a predominance of light precipitation and only rare events with snowfall &gt;1 mm h−1 water equivalent (w.e.). Three of such intense snowfall events occurred during 2011 contributing to anomalously large annual surface mass balance (SMB). Large accumulation events (&gt;10 mm w.e. day−1) during the radar-measurement period of 26 months were always associated with snowfall, but at the same time other snowfall events did not always lead to accumulation. The multiyear deployment of a precipitation radar in Antarctica allows for assessing the contribution of the snowfall to the local SMB and comparing it to the other SMB components. During 2012, snowfall rate was 110 ± 20 mm w.e. yr−1, from which surface and drifting snow sublimation removed together 23%. Given the measured yearly SMB of 52 ± 3 mm w.e., the residual term of 33 ± 21 mm w.e. yr−1 was attributed to the wind-driven snow erosion. In general, this promising set of robust instrumentation allows for improved insight into cloud and precipitation processes in Antarctica and can be easily deployed at other Antarctic stations.

Surface mass balance model intercomparison for the Greenland ice sheet

Abstract. A number of high resolution reconstructions of the surface mass balance (SMB) of the Greenland ice sheet (GrIS) have been produced using global re-analyses data extending back to 1958. These reconstructions have been used in a variety of applications but little is known about their consistency with each other and the impact of the downscaling method on the result. Here, we compare four reconstructions for the period 1960–2008 to assess the consistency in regional, seasonal and integrated SMB components. Total SMB estimates for the GrIS are in agreement within 34% of the four model average when a common ice sheet mask is used. When models' native land/ice/sea masks are used this spread increases to 57%. Variation in the spread of components of SMB from their mean: runoff 42% (29% native masks), precipitation 20% (24% native masks), melt 38% (74% native masks), refreeze 83% (142% native masks) show, with the exception of refreeze, a similar level of agreement once a common mask is used. Previously noted differences in the models' estimates are partially explained by ice sheet mask differences. Regionally there is less agreement, suggesting spatially compensating errors improve the integrated estimates. Modelled SMB estimates are compared with in situ observations from the accumulation and ablation areas. Agreement is higher in the accumulation area than the ablation area suggesting relatively high uncertainty in the estimation of ablation processes. Since the mid-1990s each model estimates a decreasing annual SMB. A similar period of decreasing SMB is also estimated for the period 1960–1972. The earlier decrease is due to reduced precipitation with runoff remaining unchanged, however, the recent decrease is associated with increased precipitation, now more than compensated for by increased melt driven runoff. Additionally, in three of the four models the equilibrium line altitude has risen since the mid-1990s, reducing the accumulation area at a rate of approximately 60 000 km2 per decade due to increased melting. Improving process representation requires further study but the use of a single accurate ice sheet mask is a logical way to reduce uncertainty among models.

Impact of spatial resolution on the modelling of the Greenland ice sheet surface mass balance between 1990–2010, using the regional climate model MAR

Abstract. With the aim to force an ice dynamical model, the Greenland ice sheet (GrIS) surface mass balance (SMB) was modelled at different spatial resolutions (15–50 km) for the period 1990–2010, using the regional climate model MAR (Modèle Atmosphérique Régional) forced by the ERA-INTERIM reanalysis. This comparison revealed that (i) the inter-annual variability of the SMB components is consistent within the different spatial resolutions investigated, (ii) the MAR model simulates heavier precipitation on average over the GrIS with decreasing spatial resolution, and (iii) the SMB components (except precipitation) can be derived from a simulation at lower resolution with an "intelligent" interpolation. This interpolation can also be used to approximate the SMB components over another topography/ice sheet mask of the GrIS. These results are important for the forcing of an ice dynamical model needed to enable future projections of the GrIS contribution to sea level rise over the coming centuries.

Spatial and temporal variability of surface mass balance near Talos Dome, East Antarctica

Predictions concerning Antarctica's contribution to sea level change have been hampered by poor knowledge of surface mass balance. Snow accumulation is the most direct climate indicator and has important implications for paleoclimatic reconstruction from ice cores. Snow accumulation measurements (stake, core, snow radar) taken along a 500‐km transect crossing Talos Dome (East Antarctica) have been used to assess accumulation signals and the representativeness of ice core records. Stake readings show that accumulation hiatuses can occur at sites with accumulation rates below 120 kg m−2 yr−1. Differences between cores and stakes can lead to statistical misidentification of annual layers determined from seasonal signals at sites with accumulation rates below 200 kg m−2 yr−1 because of nondetection of higher and lower values. Achieving ±10% accuracy in the reconstruction of snow accumulation from single cores requires high accumulation (750 kg m−2 yr−1). Low‐accumulation sites are representative if cumulative rates computed over several years are used to reach the 750 kg m−2 yr−1 threshold. Temporal variability of accumulation over the last two centuries shows no significant increase in accumulation. Wind‐driven processes are a fundamental component of surface mass balance. Spatial variations in accumulation are well correlated with surface slope changes along the wind direction and may exceed 200 kg m−2 yr−1 within 1 km. Wind‐driven sublimation rates are less than 50 kg m−2 yr−1 in plateau areas and up to 260 kg m−2 yr−1 in slope areas and account for 20–75% of precipitation, whereas depositional features are negligible in surface mass balance.

Precipitation, sublimation, and snow drift in the Antarctic Peninsula region from a regional atmospheric model

A regional atmospheric model, with a horizontal grid spacing (Δx) of 14 km, is used to study the surface mass balance components (precipitation, sublimation, and snow drift) in the region of the Antarctic Peninsula (AP). An integration is performed for the 7‐year period 1987–1993, using a realistic forcing at the lateral model boundaries and at the sea surface. Output from this integration indicates that the precipitation reaches its maximum value on the northwestern slope of the AP, where the upward motion in the atmosphere is largest. Uplift occurs upstream of the barrier, affecting the precipitation distribution over sea. The effect of the barrier on the precipitation distribution over the Bellingshausen Sea might have important implications for the ocean circulation in this region. The mean precipitation over the grounded ice of the AP (1.20 m water eq yr−1) is 6 times larger than the mean value over all the grounded ice of Antarctica. Our estimates for the surface sublimation and wind transport of snow over the grounding line toward the sea are 9% and 6 ± 1% of the precipitation, respectively. In situ data of the wind distribution at three coastal sites located on the northern, eastern, and western sides of the AP are used to evaluate the modeled wind field, which is important for the snow drift calculations. For two of the three sites considered, the prevailing wind direction and bimodal wind distribution are correctly represented by the model. The calculated distribution of accumulation and ablation due to snow drift shows a complex pattern. The wind removes snow from the spine of the AP, where the near‐surface flow field diverges, whereas deposition occurs mainly on the eastern slopes, where the near‐surface flow field converges. An intercomparison between two 7‐year integrations at different horizontal resolution (Δx = 14 km and Δx = 55 km) shows that the precipitation on the northwestern slope is very sensitive to the model resolution: In the Δx = 14 km integration, precipitation on the northwestern slope is higher than in Δx = 55 km because of higher vertical velocities, resulting in a 35% increase in average precipitation over the grounded ice of the AP.

Greenland ice sheet surface mass balance 1991–2000: Application of Polar MM5 mesoscale model and in situ data

The Polar Pennsylvania State University–National Center for Atmospheric Research Fifth‐Generation Mesoscale Model (Polar MM5) regional climate model was run over the North Atlantic region for 1991–2000. We analyze 24‐km output over the Greenland ice sheet to evaluate spatial and temporal variability of the surface mass balance and its subcomponents. The model output is compared with 3 years of automatic weather station (AWS) data from 17 sites to identify biases. Using the in situ data, we derive simple corrections for biases in melt energy and in water vapor fluxes from the surface and from blowing snow. The simulated accumulation rate is in agreement with AWS and snow pit observations. Estimates for runoff and the surface mass balance distribution over the ice sheet are produced using modeled melt volume and a meltwater retention scheme. From the decade investigated, the magnitude of interannual variability in surface mass balance components is tentatively established. The largest variability is concentrated along the ice sheet margin, where both accumulation and ablation rates are largest. The simulated interannual fluctuations suggest a large absolute variability, ±187 km3 yr−1 for total ice sheet surface mass balance. Variability in simulated equilibrium line altitude is suggestive of a dominance of thermal variability in the south with increasing importance of accumulation variability with increasing latitude. Empirical functions for the sensitivity of surface mass balance to temperature and precipitation anomalies are presented. The precise locations and regions of maximum and minimum surface energy and mass fluxes are suggested. Using an estimate for iceberg discharge and bottom melting, the total ice sheet mass balance is estimated be −78 km3 yr−1, producing 2.2 mm of eustatic sea level rise over the 1991–2000 decade and contributing 15% to the observed (1.5 mm yr−1) global sea level rise. The more negative mass balance is attributed to including blowing‐snow sublimation loss and to regional warming in the 1990s.