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.
@en

About 20 % of all snow accumulation in Antarctica occurs on the ice shelves. There, ice rises control the spatial surface mass balance (SMB) distribution by inducing snowfall variability and wind erosion due to their topography. Moreover these ice rises buttress the ice flow and represent ideal drilling locations for ice cores. In this study we assess the connection between snowfall variability and wind erosion to provide a better understanding of how ice rises impact SMB variability, how well this is captured in the regional atmospheric climate model RACMO2 and the implications of this SMB variability for ice rises as an ice core drilling site. By combining ground-penetrating radar (GPR) profiles from two ice rises in Dronning Maud Land with ice core dating, we reconstruct spatial and temporal SMB variations from 1983 to 2018 and compare the observed SMB with output from RACMO2 and SnowModel. Our results show snowfall-driven differences of up to 1.5 times higher SMB on the windward side of both ice rises than on the leeward side as well as a local erosion-driven minimum at the ice divide of the ice rises. RACMO2 captures the snowfall-driven differences but overestimates their magnitude, whereas the erosion on the peak can be reproduced by SnowModel with RACMO2 forcing. Observed temporal variability of the average SMBs, retrieved from the GPR data for four time intervals in the 1983-2018 range, are low at the peak of the easternmost ice rise (∼ 0.06 m w.e. yr⁻¹), while they are higher (∼ 0.09 m w.e. yr⁻¹) on the windward side of the ice rise. This implies that at the peak of the ice rise, higher snowfall, driven by orographic uplift, is balanced out by local erosion. As a consequence of this, the SMB recovered from the ice core matches the SMB from the GPR at the peak of the ice rise but not at the windward side of the ice rise, suggesting that the SMB signal is damped in the ice core.

@en

Snow and ice albedo schemes in present-day climate models often lack a sophisticated radiation penetration scheme and do not explicitly include spectral albedo variations. In this study, we evaluate a new snow albedo scheme in the Regional Atmospheric Climate Model (RACMO2) for the Greenland ice sheet, version 2.3p3, that includes these processes. The new albedo scheme uses the TwostreAm Radiative TransfEr in Snow (TARTES) model and the Spectral-To-NarrOWBand ALbedo (SNOWBAL) module, version 1.2. Additionally, the bare-ice albedo parameterization has been updated. The snow and ice broadband and narrowband albedo output of the updated version of RACMO2 is evaluated using the Programme for Monitoring of the Greenland Ice Sheet (PROMICE) and Kangerlussuaq transect (K-Transect) in situ data and Moderate Resolution Imaging Spectroradiometer (MODIS) remote-sensing observations. Generally, the modeled narrowband and broadband albedo is in very good agreement with satellite observations, leading to a negligible domain-Averaged broadband albedo bias for the interior. Some discrepancies are, however, observed close to the ice margin. Compared to the previous model version, RACMO2.3p2, the broadband albedo is considerably higher in the bare-ice zone during the ablation season, as atmospheric conditions now alter the bare-ice broadband albedo. For most other regions, however, the updated broadband albedo is lower due to spectral effects, radiation penetration or enhanced snow metamorphism.

@en

Since the early 1990s, the Greenland ice sheet (GrIS) has been losing mass at an accelerating rate, primarily due to enhanced meltwater runoff following atmospheric warming. Here, we show that a pronounced latitudinal contrast exists in the GrIS response to recent warming. The ablation area in north Greenland expanded by 46%, almost twice as much as in the south (+25%), significantly increasing the relative contribution of the north to total GrIS mass loss. This latitudinal contrast originates from a different response to the recent change in large-scale Arctic summertime atmospheric circulation, promoting southwesterly advection of warm air toward the GrIS. In the southwest, persistent high atmospheric pressure reduced cloudiness, increasing runoff through enhanced absorption of solar radiation; in contrast, increased early-summer cloudiness in north Greenland enhanced atmospheric warming through decreased longwave heat loss. This triggered a rapid snowline retreat, causing early bare ice exposure, amplifying northern runoff.

@en

We evaluate modelled Antarctic ice sheet (AIS) near-surface climate, surface mass balance (SMB) and surface energy balance (SEB) from the updated polar version of the regional atmospheric climate model, RACMO2 (1979-2016). The updated model, referred to as RACMO2.3p2, incorporates upper-air relaxation, a revised topography, tuned parameters in the cloud scheme to generate more precipitation towards the AIS interior and modified snow properties reducing drifting snow sublimation and increasing surface snowmelt.

Comparisons of RACMO2 model output with several independent observational data show that the existing biases in AIS temperature, radiative fluxes and SMB components are further reduced with respect to the previous model version. The model-integrated annual average SMB for the ice sheet including ice shelves (minus the Antarctic Peninsula, AP) now amounts to 2229ĝ€Gtĝ€yĝ'1, with an interannual variability of 109ĝ€Gtĝ€yĝ'1. The largest improvement is found in modelled surface snowmelt, which now compares well with satellite and weather station observations. For the high-resolution ( ĝ1/4 ĝ€5.5ĝ€km) AP simulation, results remain comparable to earlier studies.

The updated model provides a new, high-resolution data set of the contemporary near-surface climate and SMB of the AIS; this model version will be used for future climate scenario projections in a forthcoming study.

@en

We evaluate modelled Greenland ice sheet (GrIS) near-surface climate, surface energy balance (SEB) and surface mass balance (SMB) from the updated regional climate model RACMO2 (1958-2016). The new model version, referred to as RACMO2.3p2, incorporates updated glacier outlines, topography and ice albedo fields. Parameters in the cloud scheme governing the conversion of cloud condensate into precipitation have been tuned to correct inland snowfall underestimation: snow properties are modified to reduce drifting snow and melt production in the ice sheet percolation zone. The ice albedo prescribed in the updated model is lower at the ice sheet margins, increasing ice melt locally. RACMO2.3p2 shows good agreement compared to in situ meteorological data and point SEB/SMB measurements, and better resolves the spatial patterns and temporal variability of SMB compared with the previous model version, notably in the north-east, south-east and along the K-transect in south-western Greenland. This new model version provides updated, high-resolution gridded fields of the GrIS present-day climate and SMB, and will be used for projections of the GrIS climate and SMB in response to a future climate scenario in a forthcoming study.

@en

The Canadian Arctic Archipelago comprises multiple small glaciers and ice caps, mostly concentrated on Ellesmere and Baffin Islands in the northern (NCAA, Northern Canadian Arctic Archipelago) and southern parts (SCAA, Southern Canadian Arctic Archipelago) of the archipelago, respectively. Because
these glaciers are small and show complex geometries, current regional climate models, using 5- to 20-km horizontal resolution, do not properly resolve surface mass balance patterns. Here we present a 58-year (1958–2015) reconstruction of daily surface mass balance of the Canadian Arctic Archipelago, statistically
downscaled to 1 km from the output of the regional climate model RACMO2.3 at 11 km. By correcting for biases in elevation and ice albedo, the downscaling method significantly improves runoff estimates over narrow outlet glaciers and isolated ice fields. Since the last two decades, NCAA and SCAA glaciers
have experienced warmer conditions (+1.1 ∘ C) resulting in continued mass loss of 28.2±11.5 and 22.0±4.5 Gt/year, respectively, more than doubling (11.9 Gt/year) and doubling (11.9 Gt/year) the pre-1996 average. While the interior of NCAA ice caps can still buffer most of the additional melt, the lack
of a perennial firn area over low-lying SCAA glaciers has caused uninterrupted mass loss since the 1980s. In the absence of significant refreezing capacity, this indicates inevitable disappearance of these highly sensitive glaciers.
@en

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.@en