Runoff has recently become the main source of mass loss from the Greenland Ice Sheet and is an important contributor to global sea level rise. Linking runoff to surface meltwater production is complex, as meltwater can be retained within the firn by refreezing or perennial liquid water storage. To constrain these uncertainties, the outputs of two offline snow/firn models of different complexity (IMAU-FDM and SNOWPACK) are compared to assess the sensitivity of meltwater retention to the model formulation (e.g., densification, irreducible water content, vertical resolution). Results indicate that model differences are largest in areas where firn aquifers form, i.e., particularly along the south-eastern margin of the ice sheet. The IMAU-FDM simulates higher densification rates for such climatic conditions and prescribes a lower irreducible water content than SNOWPACK. As a result, the model predicts substantially lower amounts of refreezing and liquid water storage. SNOWPACK performs better for this area, confirmed both by density profiles from firn cores and radar-inferred observations. Refreezing integrated over the entire ice sheet and averaged for the period 1960 - 2014 amounts to 216 Gt a-1 (IMAU-FDM) and 242 Gt a-1 (SNOWPACK), which is 41% and 46% of the total liquid water input (snowmelt and rainfall). The mean areal extents of perennial firn aquifers for 2010 - 2014 simulated by the models are 55,700 km2 (IMAU-FDM) and 90,200 km2 (SNOWPACK). Discrepancies between modeled firn profiles and observations emphasize the importance of processes currently not accounted for in most snow/firn models, such as vertical heterogeneous percolation, ponding of water on impermeable layers, lateral (sub-)surface water flow, and the issue of ill-constrained refreezing conditions at the base of firn aquifers.

The recent development of high-resolution climate models offers a promising approach in improving the simulation of precipitation, clouds and temperature. However, higher grid spacing is also a promising feature to improve the simulation of snow cover. In particular, it provides a refined representation of topography and allows for an explicit simulation of convective precipitation processes. In this study we analyze the snow cover in a set of decade-long high-resolution climate simulation with horizontal grid spacing of 2.2 km over the greater Alpine region. Results are compared against observations and lower resolution models (12 and 50 km), which use parameterized convection. The simulations are integrated using the COSMO (Consortium for Small-Scale Modeling) model. The evaluation of snow water equivalent (SWE) in the simulation of present-day climate, driven by the ERA-Interim reanalysis, against an observational dataset, reveals that the high-resolution simulation clearly outperforms simulations with grid spacing of 12 and 50 km. The latter simulations underestimate the cumulative amount of SWE over Switzerland over the whole annual cycle by 33% (12 km simulation) and 56% (50 km simulation) while the high-resolution simulation shows a spatially and temporally averaged difference of less than 1%. Scenario simulations driven by GCM MPI-ESM-LR (2081–2090 RCP8.5 vs. 1991–2000) reveal a strong decrease of SWE over the Alps, consistent with previous studies. Previous studies had found that the relative decrease becomes gradually smaller with elevation, but this finding was limited to low and intermediate altitudes (as a 12 km simulation resolves the topography up to 2500 m). In the current study we find that the height gradient reverses sign, and relative reductions in snow cover increases above 3000 m asl, where important parts of the cryosphere are present. In addition, the simulations project a transition from permanent to seasonal snow cover at high altitudes, with potentially important impacts to Alpine permafrost. This transition and the more pronounced decline of SWE emphasize the value of the higher grid spacing. Overall, we show that high-resolution climate models offer a promising approach in improving the simulation of snow cover in Alpine terrain.

Recent studies indicate that the surface mass balance will dominate the Greenland Ice Sheet's (GrIS) contribution to 21st century sea level rise. Consequently, it is crucial to understand the liquid water balance (LWB) of the ice sheet and its response to increasing surface melt. We therefore analyse a firn simulation conducted with the SNOWPACK model for the GrIS and over the period 1960–2014 with a special focus on the LWB and refreezing. Evaluations of the simulated refreezing climate with GRACE and firn temperature observations indicate a good model–observation agreement. Results of the LWB analysis reveal a spatially uniform increase in surface melt (0.16 m w.e. a⁻¹) during 1990–2014. As a response, refreezing and run-off also indicate positive changes during this period (0.05 and 0.11 m w.e. a⁻¹, respectively), where refreezing increases at only half the rate of run-off, implying that the majority of the additional liquid input runs off the ice sheet. This pattern of refreeze and run-off is spatially variable. For instance, in the south-eastern part of the GrIS, most of the additional liquid input is buffered in the firn layer due to relatively high snowfall rates. Modelled increase in refreezing leads to a decrease in firn air content and to a substantial increase in near-surface firn temperature. On the western side of the ice sheet, modelled firn temperature increases are highest in the lower accumulation zone and are primarily caused by the exceptional melt season of 2012. On the eastern side, simulated firn temperature increases are more gradual and are associated with the migration of firn aquifers to higher elevations.

In this study, we focus on the model detection in the Antarctic Peninsula (AP) of so-called perennial firn aquifers (PFAs) that are widespread in Greenland and Svalbard and are formed when surface meltwater percolates into the firn pack in summer, which is then buried by snowfall and does not refreeze during the following winter. We use two snow models, the Institute for Marine and Atmospheric Research Utrecht Firn Densification Model (IMAU-FDM) and SNOWPACK, and force these (partly) with mass and energy fluxes from the Regional Atmospheric Climate MOdel (RACMO2.3p2) to construct a 1979–2016 climatology of AP firn density, temperature, and liquid water content. An evaluation using 75 snow temperature observations at 10 m depth and density profiles from 11 firn cores shows that output of both snow models is sufficiently realistic to warrant further analysis of firn characteristics. The models give comparable results: in 941 model grid points in either model, covering ∼28 000 km², PFAs existed for at least 1 year in the simulated period, most notably in the western AP. At these locations, surface meltwater production typically exceeds 200 $\mathrm{mm}\mathrm{w}.\mathrm{e}.{\mathrm{yr}}^{-\mathrm{1}}$ , with accumulation for most locations $>\mathrm{1000}\mathrm{mm}\mathrm{w}.\mathrm{e}.{\mathrm{yr}}^{-\mathrm{1}}$ . Most persistent and extensive are PFAs modelled on and around Wilkins Ice Shelf. Here, both meltwater production and accumulation rates are sufficiently high to sustain a PFA on 49 % of the ice shelf area in (up to) 100 % (depending on the model) of the years in the 1979–2016 period. Although this PFA presence is confirmed by recent observations, its extent in the models appears underestimated. Other notable PFA locations are Wordie Ice Shelf, an ice shelf that has almost completely disappeared in recent decades, and the relatively warm north-western side of mountain ranges in Palmer Land, where accumulation rates can be extremely high, and PFAs are formed frequently. PFAs are not necessarily more frequent in areas with the largest melt and accumulation rates, but they do grow larger and retain more meltwater, which could increase the likelihood of ice shelf hydrofracturing. We find that not only the magnitude of melt and accumulation is important but also the timing of precipitation events relative to melt events. Large accumulation events that occur in the months following an above-average summer melt event favour PFA formation in that year. Most PFAs are predicted near the grounding lines of the (former) Prince Gustav, Wilkins, and Wordie ice shelves. This highlights the need to further investigate how PFAs may impact ice shelf disintegration events through the process of hydrofracturing in a similar way as supraglacial lakes do.

Terrain parameters like topographic horizon and sky view factor (SVF) are used in numerous fields and applications. In atmospheric and climate modelling, such parameters are utilised to parameterise the effect of terrain geometry on radiation exchanges between the surface and the atmosphere. Ideally, these parameters are derived from a high-resolution digital elevation model (DEM) because inferring them from coarser elevation data induces a smoothing effect. Computing topographic horizon with conventional algorithms, however, is slow because large amounts of non-local terrain data have to be processed. We propose a new and more efficient method, which is based on a high-performance ray-tracing library. The new algorithm can speed up horizon calculation by 2 orders of magnitude relative to a conventional approach. By applying terrain simplification to remote topography, the ray-tracing-based algorithm can also be applied with very high-resolution (<5 m) DEM data, which would otherwise induce an excessive memory footprint. The topographic horizon algorithm is accompanied by an SVF algorithm, which was verified to work accurately for all terrain – even very steep and complex terrain. We compare the computational performance and accuracy of the new horizon algorithm with two reference methods from the literature and illustrate its benefits. Finally, we illustrate how sub-grid SVF values can be efficiently computed with the newly derived horizon algorithm for a wide range of target grid resolutions (1–25 km).

The Hengduan Mountains (HM) are located on the southeastern edge of the Tibetan Plateau and feature high mountain ridges (>6,000 m MSL) separated by deep valleys. The HM region also features an exceptionally high biodiversity, believed to have emerged from the topography interacting with the climate. To investigate the role of the HM topography on regional climate, we conduct simulations with the regional climate model COSMO at high horizontal resolutions (at ∼12 km and a convection‐permitting scale of ∼4.4 km) for the period 2001–2005. We conduct one control simulation with modern topography and two idealized experiments with modified topography, inspired by past geological processes that shaped the mountain range. In the first experiment, we reduce the HM's elevation by applying a spatially non‐uniform scaling to the topography. The results show that, following the uplift of the HM, the local rainy season precipitation increases by ∼25%. Precipitation in Indochina and the Bay of Bengal (BoB) also intensifies. Additionally, the cyclonic circulation in the BoB extends eastward, indicating an intensification of the East Asian summer monsoon. In the second experiment, we remove deep valleys by applying an envelope topography to quantify the effects of terrain undulation with high amplitude and frequency on climate. On the western flanks of the HM, precipitation slightly increases, while the remaining fraction of the mountain range experiences ∼20% less precipitation. Simulations suggest an overall positive feedback between precipitation, erosion, and valley deepening for this region, which could have influenced the diversification of local organisms.