Greenland Ice Sheet Surface Mass Research Articles

The Effect of Physically Based Ice Radiative Processes on Greenland Ice Sheet Albedo and Surface Mass Balance in E3SM

AbstractA significant portion of surface melt on the Greenland Ice Sheet (GrIS) is due to dark ice regions in the ablation zone, where solar absorption is influenced by the physical properties of the ice, light absorbing constituents (LACs), and the overlying crustal surface or melt ponds. Earth system models (ESMs) typically prescribe the albedo of ice surfaces as a constant value in the visible and near‐infrared spectral regions. This work advances ESM ice radiative transfer modeling by (a) incorporating a physically based radiative transfer model (SNow, ICe and Aerosol Radiation model Adding‐Doubling Version 4; SNICAR‐ADv4) into the Energy Exascale Earth System Model (E3SM), (b) determining spatially and temporally varying bare ice physical properties over the GrIS ablation zone from satellite observations to inform SNICAR‐ADv4, and (c) assessing the impacts on simulated GrIS albedo and surface mass balance associated with modeling of more realistic bare ice albedo. GrIS‐wide bare ice albedo in E3SMv2 is overestimated by ∼4% in the visible and ∼7% in the near‐infrared wavelengths compared to the Moderate Resolution Imaging Spectroradiometer. Our bare ice physical property retrieval method found that LACs, ice crustal surfaces, and melt ponds reduce visible albedo by 30% in the bare ice region of the GrIS ablation zone. The realistic bare ice albedo reduces surface mass balance by ∼145 Gt, or 0.4 mm of sea‐level equivalent between 2000 and 2021 compared to the default E3SM. This work highlights the importance of simulating bare ice albedo accurately and realistically to improve our ability to quantify changes in the GrIS surface mass and radiative energy budgets.

A stochastic parameterization of ice sheet surface mass balance for the Stochastic Ice-Sheet and Sea-Level System Model (StISSM v1.0)

Abstract. Many scientific and societal questions that draw on ice sheet modeling necessitate sampling a wide range of potential climatic changes and realizations of internal climate variability. For example, coastal planning literature demonstrates a demand for probabilistic sea level projections with quantified uncertainty. Further, robust attribution of past and future ice sheet change to specific processes or forcings requires a full understanding of the space of possible ice sheet behaviors. The wide sampling required to address such questions is computationally infeasible with sophisticated numerical climate models at the resolution required to accurately force ice sheet models. Stochastic generation of climate forcing of ice sheets offers a complementary alternative. Here, we describe a method to construct a stochastic generator for ice sheet surface mass balance varying in time and space. We demonstrate the method with an application to Greenland Ice Sheet surface mass balance for 1980–2012. We account for spatial correlations among glacier catchments using sparse covariance techniques, and we apply an elevation-dependent downscaling to recover gridded surface mass balance fields suitable for forcing an ice sheet model while including feedback from changing ice sheet surface elevation. The efficiency gained in the stochastic method supports large-ensemble simulations of ice sheet change in a new stochastic ice sheet model. We provide open source Python workflows to support use of our stochastic approach for a broad range of applications.

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.

High‐Latitude Stratospheric Aerosol Injection to Preserve the Arctic

AbstractStratospheric aerosol injection (SAI) has been shown in climate models to reduce some impacts of global warming in the Arctic, including the loss of sea ice, permafrost thaw, and reduction of Greenland Ice Sheet (GrIS) mass; SAI at high latitudes could preferentially target these impacts. In this study, we use the Community Earth System Model to simulate two Arctic‐focused SAI strategies, which inject at 60°N latitude each spring with injection rates adjusted to either maintain September Arctic sea ice at 2030 levels (“Arctic Low”) or restore it to 2010 levels (“Arctic High”). Both simulations maintain or restore September sea ice to within 10% of their respective targets, reduce permafrost thaw, and increase GrIS surface mass balance by reducing runoff. Arctic High reduces these impacts more effectively than a globally focused SAI strategy that injects similar quantities of SO2 at lower latitudes. However, Arctic‐focused SAI is not merely a “reset button” for the Arctic climate, but brings about a novel climate state, including changes to the seasonal cycles of Northern Hemisphere temperature and sea ice and less high‐latitude carbon uptake relative to SSP2‐4.5. Additionally, while Arctic‐focused SAI produces the most cooling near the pole, its effects are not confined to the Arctic, including detectable cooling throughout most of the northern hemisphere for both simulations, increased mid‐latitude sulfur deposition, and a southward shift of the location of the Intertropical Convergence Zone. For these reasons, it would be incorrect to consider Arctic‐focused SAI as “local” geoengineering, even when compared to a globally focused strategy.

Impact of Grids and Dynamical Cores in CESM2.2 on the Surface Mass Balance of the Greenland Ice Sheet

AbstractSix different configurations, a mixture of grids and atmospheric dynamical cores available in the Community Earth System Model, version 2.2 (CESM2.2), are evaluated for their skill in representing the climate of the Arctic and the surface mass balance of the Greenland Ice Sheet (GrIS). The finite‐volume dynamical core uses structured, latitude‐longitude grids, whereas the spectral‐element dynamical core is built on unstructured meshes, permitting grid flexibility such as quasi‐uniform grid spacing globally. The 1°–2° latitude‐longitude and quasi‐uniform unstructured grids systematically overestimate both accumulation and ablation over the GrIS. Of these 1°–2° grids, the latitude‐longitude grids outperform the quasi‐uniform unstructured grids because they have more degrees of freedom to represent the GrIS. Two Arctic‐refined meshes, with 1/4° and 1/8° refinement over Greenland, were developed for the spectral‐element dynamical core and are documented here as newly supported configurations in CESM2.2. The Arctic meshes substantially improve the simulated clouds and precipitation rates in the Arctic. Over Greenland, these meshes skillfully represent accumulation and ablation processes, leading to a more realistic GrIS surface mass balance. As CESM is in the process of transitioning away from conventional latitude‐longitude grids, these new Arctic‐refined meshes improve the representation of polar processes in CESM by recovering resolution lost in the transition to quasi‐uniform grids, albeit at increased computational cost.

Reanalysis Surface Mass Balance of the Greenland Ice Sheet Along K‐Transect (2000–2014)

AbstractAccurate estimates of surface mass balance over the Greenland ice sheet (GrIS) would contribute to understanding the cause of recent changes and would help to better estimate the future contribution of the GrIS to sea‐level rise. Given the limitations of in‐situ measurements, modeling, and remote sensing, it is critical to explore the opportunity to merge the available data to better characterize the spatial and temporal variation of the GrIS surface mass balance (SMB). This work utilizes a particle batch smoother data assimilation technique that yields SMB estimates that benefit from the snow model Crocus and a 16‐day albedo product derived from satellite remote sensing data. Comparison of the results against in‐situ SMB measurements shows that the assimilation of the albedo product reduces the root mean square error of the posterior estimates of SMB by 51% and reduces bias by 95%.

Influence of Arctic sea-ice loss on the Greenland ice sheet climate

The Arctic is the region on Earth that is warming the fastest. At the same time, Arctic sea ice is reducing while the Greenland ice sheet (GrIS) is losing mass at an accelerated pace. Here, we study the seasonal impact of reduced Arctic sea ice on GrIS surface mass balance (SMB), using the Community Earth System Model version 2.1 (CESM2), which features an advanced, interactive calculation of SMB. Addressing the impact of sea-ice reductions on the GrIS SMB from observations is difficult due to the short observational records. Also, signals detected using transient climate simulations may be aliases of other forcings. Here, we analyze dedicated simulations from the Polar Amplification Model Intercomparison Project with reduced Arctic sea ice and compare them with preindustrial sea ice simulations while keeping all other forcings constant. In response to reduced sea ice, the GrIS SMB increases in winter due to increased precipitation, driven by the more humid atmosphere and increasing cyclones. In summer, surface melt increases due to a warmer, more humid atmosphere providing increased energy transfer to the surface through the sensible and latent heat fluxes, which triggers the melt-albedo feedback. Further, warming occurs throughout the entire troposphere over Baffin Bay. This deep warming results in regional enhancement of the 500 hPa geopotential heights over the Baffin Bay and Greenland, increasing blocking and heat advection over the GrIS’ surface. This anomalous circulation pattern has been linked to recent increases in the surface melt of the GrIS.

Brief communication: Reduction in the future Greenland ice sheet surface melt with the help of solar geoengineering

Abstract. The Greenland Ice Sheet (GrIS) will be losing mass at an accelerating pace throughout the 21st century, with a direct link between anthropogenic greenhouse gas emissions and the magnitude of Greenland mass loss. Currently, approximately 60 % of the mass loss contribution comes from surface melt and subsequent meltwater runoff, while 40 % are due to ice calving. In the ablation zone covered by bare ice in summer, most of the surface melt energy is provided by absorbed shortwave fluxes, which could be reduced by solar geoengineering measures. However, so far very little is known about the potential impacts of an artificial reduction in the incoming solar radiation on the GrIS surface energy budget and the subsequent change in meltwater production. By forcing the regional climate model MAR with the latest CMIP6 shared socioeconomic pathways (SSP) future emission scenarios (SSP245, SSP585) and associated G6solar experiment from the CNRM-ESM2-1 Earth system model, we estimate the local impact of a reduced solar constant on the projected GrIS surface mass balance (SMB) decrease. Overall, our results show that even in the case of a low-mitigation greenhouse gas emissions scenario (SSP585), the Greenland surface mass loss can be brought in line with the medium-mitigation emissions scenario (SSP245) by reducing the solar downward flux at the top of the atmosphere by ∼ 40 W/m2 or ∼ 1.5 % (using the G6solar experiment). In addition to reducing global warming in line with SSP245, G6solar also decreases the efficiency of surface meltwater production over the Greenland ice sheet by damping the well-known positive melt–albedo feedback. With respect to a MAR simulation where the solar constant remains unchanged, decreasing the solar constant according to G6solar in the MAR radiative scheme mitigates the projected Greenland ice sheet surface melt increase by 6 %. However, only more constraining geoengineering experiments than G6solar would allow us to maintain a positive SMB until the end of this century without any reduction in our greenhouse gas emissions.

Surface melting over the Greenland ice sheet derived from enhanced resolution passive microwave brightness temperatures (1979–2019)

Abstract. Surface melting is a major component of the Greenland ice sheet surface mass balance, and it affects sea level rise through direct runoff and the modulation of ice dynamics and hydrological processes, supraglacially, englacially and subglacially. Passive microwave (PMW) brightness temperature observations are of paramount importance in studying the spatial and temporal evolution of surface melting due to their long temporal coverage (1979–present) and high temporal resolution (daily). However, a major limitation of PMW datasets has been the relatively coarse spatial resolution, which has historically been of the order of tens of kilometers. Here, we use a newly released PMW dataset (37 GHz, horizontal polarization) made available through a NASA “Making Earth System Data Records for Use in Research Environments” (MeASUREs) program to study the spatiotemporal evolution of surface melting over the Greenland ice sheet at an enhanced spatial resolution of 3.125 km. We assess the outputs of different detection algorithms using data collected by automatic weather stations (AWSs) and the outputs of the Modèle Atmosphérique Régional (MAR) regional climate model. We found that sporadic melting is well captured using a dynamic algorithm based on the outputs of the Microwave Emission Model of Layered Snowpack (MEMLS), whereas a fixed threshold of 245 K is capable of detecting persistent melt. Our results indicate that, during the reference period from 1979 to 2019 (from 1988 to 2019), surface melting over the ice sheet increased in terms of both duration, up to 4.5 (2.9) d per decade, and extension, up to 6.9 % (3.6 %) of the entire ice sheet surface extent per decade, according to the MEMLS algorithm. Furthermore, the melting season started up to 4.0 (2.5) d earlier and ended 7.0 (3.9) d later per decade. We also explored the information content of the enhanced-resolution dataset with respect to the one at 25 km and MAR outputs using a semi-variogram approach. We found that the enhanced product is more sensitive to local-scale processes, thereby confirming the potential of this new enhanced product for monitoring surface melting over Greenland at a higher spatial resolution than the historical products and for monitoring its impact on sea level rise. This offers the opportunity to improve our understanding of the processes driving melting, to validate modeled melt extent at high resolution and, potentially, to assimilate these data in climate models.

GrSMBMIP: intercomparison of the modelled 1980–2012 surface mass balance over the Greenland Ice Sheet

Abstract. Observations and models agree that the Greenland Ice Sheet (GrIS) surface mass balance (SMB) has decreased since the end of the 1990s due to an increase in meltwater runoff and that this trend will accelerate in the future. However, large uncertainties remain, partly due to different approaches for modelling the GrIS SMB, which have to weigh physical complexity or low computing time, different spatial and temporal resolutions, different forcing fields, and different ice sheet topographies and extents, which collectively make an inter-comparison difficult. Our GrIS SMB model intercomparison project (GrSMBMIP) aims to refine these uncertainties by intercomparing 13 models of four types which were forced with the same ERA-Interim reanalysis forcing fields, except for two global models. We interpolate all modelled SMB fields onto a common ice sheet mask at 1 km horizontal resolution for the period 1980–2012 and score the outputs against (1) SMB estimates from a combination of gravimetric remote sensing data from GRACE and measured ice discharge; (2) ice cores, snow pits and in situ SMB observations; and (3) remotely sensed bare ice extent from MODerate-resolution Imaging Spectroradiometer (MODIS). Spatially, the largest spread among models can be found around the margins of the ice sheet, highlighting model deficiencies in an accurate representation of the GrIS ablation zone extent and processes related to surface melt and runoff. Overall, polar regional climate models (RCMs) perform the best compared to observations, in particular for simulating precipitation patterns. However, other simpler and faster models have biases of the same order as RCMs compared with observations and therefore remain useful tools for long-term simulations or coupling with ice sheet models. Finally, it is interesting to note that the ensemble mean of the 13 models produces the best estimate of the present-day SMB relative to observations, suggesting that biases are not systematic among models and that this ensemble estimate can be used as a reference for current climate when carrying out future model developments. However, a higher density of in situ SMB observations is required, especially in the south-east accumulation zone, where the model spread can reach 2 m w.e. yr−1 due to large discrepancies in modelled snowfall accumulation.

Global Warming Threshold and Mechanisms for Accelerated Greenland Ice Sheet Surface Mass Loss.

The Community Earth System Model version 2.1 (CESM2.1) is used to investigate the evolution of the Greenland ice sheet (GrIS) surface mass balance (SMB) under an idealized CO forcing scenario of 1% increase until stabilization at 4 pre‐industrial at model year 140. In this simulation, the SMB calculation is coupled with the atmospheric model, using a physically based surface energy balance scheme for melt, explicit calculation of snow albedo, and a realistic treatment of polar snow and firn compaction. By the end of the simulation (years 131–150), the SMB decreases with 994 Gt yr with respect to the pre‐industrial SMB, which represents a sea‐level rise contribution of 2.8 mm yr. For a threshold of 2.7‐K global temperature increase with respect to pre‐industrial, the rate of expansion of the ablation area increases, the mass loss accelerates due to loss of refreezing capacity and accelerated melt, and the SMB becomes negative 6 years later. Before acceleration, longwave radiation is the most important contributor to increasing energy for melt. After acceleration, the large expansion of the ablation area strongly reduces surface albedo. This and much increased turbulent heat fluxes as the GrIS‐integrated summer surface temperature approaches melt point become the major sources of energy for melt.

The Spatiotemporal Variability of Cloud Radiative Effects on the Greenland Ice Sheet Surface Mass Balance

AbstractTo better understand and quantify the impact of clouds on the Greenland Ice Sheet surface mass balance (SMB), we study the spatiotemporal variability of the cloud radiative effect (CRE). The total CRE is separated in short‐term and long‐term impacts by performing multiple simulations with the SNOWPACK model for 2001‐+2010. The annual total CRE is 16.8 ± 4.5 W m−2, reducing the SMB with −157 ± 3.8 Gt yr−1. Summer cloud radiative cooling is −6.4 ± 5.7 W m−2 in the ablation area, increasing the SMB with 121 ± 2.2 Gt yr−1. The annual integrated impact is cloud‐reduced SMB of −36 Gt yr−1. The short‐term effect dominates the opposing long‐term effects through the albedo‐melt feedback. A long‐term warming effect decreases the albedo and so preconditions the surface for enhanced (summer) melt. The impact of the CRE, determined by spatial, temporal and initial conditions, explains existing conflicted views on the role of cloud radiation and emphasizes the need for accurate cloud and albedo representations in future studies.

Remapping of Greenland ice sheet surface mass balance anomalies for large ensemble sea-level change projections

Abstract. Future sea-level change projections with process-based stand-alone ice sheet models are typically driven with surface mass balance (SMB) forcing derived from climate models. In this work we address the problems arising from a mismatch of the modelled ice sheet geometry with the geometry used by the climate model. We present a method for applying SMB forcing from climate models to a wide range of Greenland ice sheet models with varying and temporally evolving geometries. In order to achieve that, we translate a given SMB anomaly field as a function of absolute location to a function of surface elevation for 25 regional drainage basins, which can then be applied to different modelled ice sheet geometries. The key feature of the approach is the non-locality of this remapping process. The method reproduces the original forcing data closely when remapped to the original geometry. When remapped to different modelled geometries it produces a physically meaningful forcing with smooth and continuous SMB anomalies across basin divides. The method considerably reduces non-physical biases that would arise by applying the SMB anomaly derived for the climate model geometry directly to a large range of modelled ice sheet model geometries.

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.

Modeling a modern-like <i>p</i>CO<sub>2</sub> warm period (Marine Isotope Stage KM5c) with two versions of an Institut Pierre Simon Laplace atmosphere–ocean coupled general circulation model

Abstract. The mid-Piacenzian warm period (3.264 to 3.025 Ma) is the most recent geological period with present-like atmospheric pCO2 and is thus expected to have exhibited a warm climate similar to or warmer than the present day. On the basis of understanding that has been gathered on the climate variability of this interval, a specific interglacial (Marine Isotope Stage KM5c, MIS KM5c; 3.205 Ma) has been selected for the Pliocene Model Intercomparison Project phase 2 (PlioMIP 2). We carried out a series of experiments according to the design of PlioMIP2 with two versions of the Institut Pierre Simon Laplace (IPSL) atmosphere–ocean coupled general circulation model (AOGCM): IPSL-CM5A and IPSL-CM5A2. Compared to the PlioMIP 1 experiment, run with IPSL-CM5A, our results show that the simulated MIS KM5c climate presents enhanced warming in mid- to high latitudes, especially over oceanic regions. This warming can be largely attributed to the enhanced Atlantic Meridional Overturning Circulation caused by the high-latitude seaway changes. The sensitivity experiments, conducted with IPSL-CM5A2, show that besides the increased pCO2, both modified orography and reduced ice sheets contribute substantially to mid- to high latitude warming in MIS KM5c. When considering the pCO2 uncertainties (+/-50 ppmv) during the Pliocene, the response of the modeled mean annual surface air temperature to changes to pCO2 (+/-50 ppmv) is not symmetric, which is likely due to the nonlinear response of the cryosphere (snow cover and sea ice extent). By analyzing the Greenland Ice Sheet surface mass balance, we also demonstrate its vulnerability under both MIS KM5c and modern warm climate.

Recent precipitation decrease across the western Greenland ice sheet percolation zone

Abstract. The mass balance of the Greenland Ice Sheet (GrIS) in a warming climate is of critical interest in the context of future sea level rise. Increased melting in the GrIS percolation zone due to atmospheric warming over the past several decades has led to increased mass loss at lower elevations. Previous studies have hypothesized that this warming is accompanied by a precipitation increase, as would be expected from the Clausius–Clapeyron relationship, compensating for some of the melt-induced mass loss throughout the western GrIS. This study tests that hypothesis by calculating snow accumulation rates and trends across the western GrIS percolation zone, providing new accumulation rate estimates in regions with sparse in situ data or data that do not span the recent accelerating surface melt. We present accumulation records from sixteen 22–32 m long firn cores and 4436 km of ground-penetrating radar, covering the past 20–60 years of accumulation, collected across the western GrIS percolation zone as part of the Greenland Traverse for Accumulation and Climate Studies (GreenTrACS) project. Trends from both radar and firn cores, as well as commonly used regional climate models, show decreasing accumulation rates of 2.4±1.5 % a−1 over the 1996–2016 period, which we attribute to shifting storm tracks related to stronger atmospheric summer blocking over Greenland. Changes in atmospheric circulation over the past 20 years, specifically anomalously strong summertime blocking, have reduced GrIS surface mass balance through both an increase in surface melting and a decrease in accumulation rates.

Regional grid refinement in an Earth system model: impacts on the simulated Greenland surface mass balance

Abstract. In this study, the resolution dependence of the simulated Greenland ice sheet surface mass balance (GrIS SMB) in the variable-resolution Community Earth System Model (VR-CESM) is investigated. Coupled atmosphere–land simulations are performed on two regionally refined grids over Greenland at 0.5∘ (∼55 km) and 0.25∘ (∼28 km), maintaining a quasi-uniform resolution of 1∘ (∼111 km) over the rest of the globe. On the refined grids, the SMB in the accumulation zone is significantly improved compared to airborne radar and in situ observations, with a general wetting (more snowfall) at the margins and a drying (less snowfall) in the interior GrIS. Total GrIS precipitation decreases with resolution, which is in line with best-available regional climate model results. In the ablation zone, CESM starts developing a positive SMB bias with increased resolution in some basins, notably in the east and the north. The mismatch in ablation is linked to changes in cloud cover in VR-CESM, and a reduced effectiveness of the elevation classes subgrid parametrization in CESM. Overall, our pilot study introduces VR-CESM as a new tool in the cryospheric sciences, which could be used to dynamically downscale SMB in scenario simulations and to force dynamical ice sheet models through the CESM coupling framework.

Simulated Greenland Surface Mass Balance in the GISS ModelE2 GCM: Role of the Ice Sheet Surface

AbstractThe rate of growth or retreat of the Greenland and Antarctic ice sheets remains a highly uncertain component of future sea level change. Here we examine the simulation of Greenland ice sheet surface mass balance (GrIS SMB) in a development branch of the ModelE2 version of the NASA Goddard Institute for Space Studies (GISS) general circulation model (GCM). GCMs are often limited in their ability to represent SMB compared with polar region regional climate models. We compare ModelE2‐simulated GrIS SMB for present‐day (1996–2005) simulations with fixed ocean conditions, at a spatial resolution of 2° latitude by 2.5° longitude (~200 km), with SMB simulated by the Modèle Atmosphérique Régionale (MAR) regional climate model (1996–2005 at a 25‐km resolution). ModelE2 SMB agrees well with MAR SMB on the whole, but there are distinct spatial patterns of differences and large differences in some SMB components. The impacts of changes to the ModelE2 surface are tested, including a subgrid‐scale representation of SMB with surface elevation classes. This has a minimal effect on ice sheet‐wide SMB but corrects local biases. Replacing fixed surface albedo with satellite‐derived values and an age‐dependent scheme has a larger impact, increasing simulated melt by 60%–100%. We also find that lower surface albedo can enhance the effects of elevation classes. Reducing ModelE2 surface roughness length to values closer to MAR reduces sublimation by ~50%. Further work is required to account for meltwater refreezing in ModelE2 and to understand how differences in atmospheric processes and model resolution influence simulated SMB.

Brief communication: Impact of the recent atmospheric circulation change in summer on the future surface mass balance of the Greenland Ice Sheet

Abstract. Since the 2000s, a change in the atmospheric circulation over the North Atlantic resulting in more frequent blocking events has favoured warmer and sunnier weather conditions over the Greenland Ice Sheet (GrIS) in summer, enhancing the melt increase. This circulation change is not represented by general circulation models (GCMs) of the Coupled Model Intercomparison Project Phase 5 (CMIP5), which do not predict any circulation change for the next century over the North Atlantic. The goal of this study is to evaluate the impact of an atmospheric circulation change (as currently observed) on projections of the future GrIS surface mass balance (SMB). We compare GrIS SMB estimates simulated by the regional climate model MAR forced by perturbed reanalysis (ERA-Interim with a temperature correction of +1, +1.5, and +2 ∘C at the MAR lateral boundaries) over 1980–2016 to projections of the future GrIS SMB from MAR simulations forced by three GCMs over selected periods for which a similar temperature increase of +1, +1.5, and +2 ∘C is projected by the GCMs in comparison to 1980–1999. Mean SMB anomalies produced with perturbed reanalysis over the climatologically stable period 1980–1999 are similar to those produced with MAR forced by GCMs over future periods characterised by a similar warming over Greenland. However, over the 2 last decades (2000–2016) when an increase in the frequency of blocking events has been observed in summer, MAR forced by perturbed reanalysis suggests that the SMB decrease could be amplified by a factor of 2 if such atmospheric conditions persist compared to projections forced by GCMs for the same temperature increase but without any circulation change.