Greenland Ice Sheet Surface Mass Balance Research Articles

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.

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.

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.

Modelling the climate and surface mass balance of polar ice sheets using RACMO2 – Part 1: Greenland (1958–2016)

Abstract. 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.

Reconstructions of the 1900–2015 Greenland ice sheet surface mass balance using the regional climate MAR model

Abstract. With the aim of studying the recent Greenland ice sheet (GrIS) surface mass balance (SMB) decrease relative to the last century, we have forced the regional climate MAR (Modèle Atmosphérique Régional; version 3.5.2) model with the ERA-Interim (ECMWF Interim Re-Analysis; 1979–2015), ERA-40 (1958–2001), NCEP–NCARv1 (National Centers for Environmental Prediction–National Center for Atmospheric Research Reanalysis version 1; 1948–2015), NCEP–NCARv2 (1979–2015), JRA-55 (Japanese 55-year Reanalysis; 1958–2014), 20CRv2(c) (Twentieth Century Reanalysis version 2; 1900–2014) and ERA-20C (1900–2010) reanalyses. While all these forcing products are reanalyses that are assumed to represent the same climate, they produce significant differences in the MAR-simulated SMB over their common period. A temperature adjustment of +1 °C (respectively −1 °C) was, for example, needed at the MAR boundaries with ERA-20C (20CRv2) reanalysis, given that ERA-20C (20CRv2) is ∼ 1 °C colder (warmer) than ERA-Interim over Greenland during the period 1980–2010. Comparisons with daily PROMICE (Programme for Monitoring of the Greenland Ice Sheet) near-surface observations support these adjustments. Comparisons with SMB measurements, ice cores and satellite-derived melt extent reveal the most accurate forcing datasets for the simulation of the GrIS SMB to be ERA-Interim and NCEP–NCARv1. However, some biases remain in MAR, suggesting that some improvements are still needed in its cloudiness and radiative schemes as well as in the representation of the bare ice albedo. Results from all MAR simulations indicate that (i) the period 1961–1990, commonly chosen as a stable reference period for Greenland SMB and ice dynamics, is actually a period of anomalously positive SMB (∼ +40 Gt yr−1) compared to 1900–2010; (ii) SMB has decreased significantly after this reference period due to increasing and unprecedented melt reaching the highest rates in the 120-year common period; (iii) before 1960, both ERA-20C and 20CRv2-forced MAR simulations suggest a significant precipitation increase over 1900–1950, but this increase could be the result of an artefact in the reanalyses that are not well-enough constrained by observations during this period and (iv) since the 1980s, snowfall is quite stable after having reached a maximum in the 1970s. These MAR-based SMB and accumulation reconstructions are, however, quite similar to those from Box (2013) after 1930 and confirm that SMB was quite stable from the 1940s to the 1990s. Finally, only the ERA-20C-forced simulation suggests that SMB during the 1920–1930 warm period over Greenland was comparable to the SMB of the 2000s, due to both higher melt and lower precipitation than normal.

Regional Greenland accumulation variability from Operation IceBridge airborne accumulation radar

Abstract. The mass balance of the Greenland Ice Sheet (GrIS) in a warming climate is of critical interest to scientists and the general public in the context of future sea-level rise. An improved understanding of temporal and spatial variability of snow accumulation will reduce uncertainties in GrIS mass balance models and improve projections of Greenland's contribution to sea-level rise, currently estimated at 0.089 ± 0.03 m by 2100. Here we analyze 25 NASA Operation IceBridge accumulation radar flights totaling > 17 700 km from 2013 to 2014 to determine snow accumulation in the GrIS dry snow and percolation zones over the past 100–300 years. IceBridge accumulation rates are calculated and used to validate accumulation rates from three regional climate models. Averaged over all 25 flights, the RMS difference between the models and IceBridge accumulation is between 0.023 ± 0.019 and 0.043 ± 0.029 m w.e. a−1, although each model shows significantly larger differences from IceBridge accumulation on a regional basis. In the southeast region, for example, the Modèle Atmosphérique Régional (MARv3.5.2) overestimates by an average of 20.89 ± 6.75 % across the drainage basin. Our results indicate that these regional differences between model and IceBridge accumulation are large enough to significantly alter GrIS surface mass balance estimates. Empirical orthogonal function analysis suggests that the first two principal components account for 33 and 19 % of the variance, and correlate with the Atlantic Multidecadal Oscillation (AMO) and wintertime North Atlantic Oscillation (NAO), respectively. Regions that disagree strongest with climate models are those in which we have the fewest IceBridge data points, requiring additional in situ measurements to verify model uncertainties.

Evaluation of the updated regional climate model RACMO2.3: summer snowfall impact on the Greenland Ice Sheet

Abstract. We discuss Greenland Ice Sheet (GrIS) surface mass balance (SMB) differences between the updated polar version of the RACMO climate model (RACMO2.3) and the previous version (RACMO2.1). Among other revisions, the updated model includes an adjusted rainfall-to-snowfall conversion that produces exclusively snowfall under freezing conditions; this especially favours snowfall in summer. Summer snowfall in the ablation zone of the GrIS has a pronounced effect on melt rates, affecting modelled GrIS SMB in two ways. By covering relatively dark ice with highly reflective fresh snow, these summer snowfalls have the potential to locally reduce melt rates in the ablation zone of the GrIS through the snow-albedo-melt feedback. At larger scales, SMB changes are driven by differences in orographic precipitation following a shift in large-scale circulation, in combination with enhanced moisture to precipitation conversion for warm to moderately cold conditions. A detailed comparison of model output with observations from automatic weather stations, ice cores and ablation stakes shows that the model update generally improves the simulated SMB-elevation gradient as well as the representation of the surface energy balance, although significant biases remain.

Inferring the variation of climatic and glaciological contributions to West Greenland iceberg discharge in the twentieth century

Iceberg discharge is a major component of the mass balance of the Greenland Ice Sheet (GrIS). While bulk estimates of discharge variation over time exist, inferred remotely from measurements of grounding line ice velocities or surface mass balance calculations, few detailed measurements of discharge itself from individual marine-terminating glaciers existed until recent years. Recently, it has been shown, through a combination of ocean–iceberg modelling and non-linear system identification, that the century-long record of iceberg numbers crossing 48oN in the West Atlantic is a good first-order proxy for discharge from at least south and west Greenland. Here, we explore the varying relative importance of ice sheet, oceanic and climatic forcing of iceberg discharge from these areas over the twentieth century, by carrying out sensitivity studies of a non-linear auto-regressive mathematical model of the 48oN time series. We find that the relationships are mainly non-linear, with the contribution of the GrIS surface mass balance to iceberg discharge likely to be dominant in the first half of the century. This period is followed by several decades where oceanic temperature effects are most important in determining the model variation in iceberg discharge. In recent decades, all physical processes play a non-negligible part in explaining the iceberg discharge and the model suggests that the glacial response time to environmental changes may have decreased.

Coupled simulations of Greenland Ice Sheet and climate change up to A.D. 2300

AbstractRecent observations indicate a high sensitivity of the Greenland Ice Sheet (GrIS) to climate change. We examine the coupling between the GrIS surface mass balance, elevation, and dynamical flow with one of the few coupled GrIS and atmosphere‐ocean general circulation models. Bidirectional coupling from the early Holocene reveals a growing present‐day GrIS in the absence of anthropogenic forcing. We identify atmospheric sources of biases in the simulated present‐day GrIS and assess the GrIS sensitivity to future greenhouse gas forcing through three Representative Concentration Pathways and their extensions and to climate variability. The elevation‐surface mass balance feedback contributes to future GrIS mass loss with 8–11% (by 2100), depending on the forcing scenario, and 24–31% (by 2300). Climate variability causes a 2.5 times spread in the magnitude of the simulated present‐day GrIS mass trends in a three‐member ensemble. Our results represent a first step toward more advanced higher resolution coupled modeling of GrIS and climate evolution.

Sensitivity of Greenland Ice Sheet surface mass balance to perturbations in sea surface temperature and sea ice cover: a study with the regional climate model MAR

Abstract. During recent summers (2007–2012), several surface melt records were broken over the Greenland Ice Sheet (GrIS). The extreme summer melt resulted in part from a persistent negative phase of the North Atlantic Oscillation (NAO), favoring warmer atmospheric conditions than normal over the GrIS. Simultaneously, large anomalies in sea ice cover (SIC) and sea surface temperature (SST) were observed in the North Atlantic, suggesting a possible connection. To assess the direct impact of 2007–2012 SIC and SST anomalies on GrIS surface mass balance (SMB), a set of sensitivity experiments was carried out with the regional climate model MAR forced by ERA-Interim. These simulations suggest that perturbations in SST and SIC in the seas surrounding Greenland do not considerably impact GrIS SMB, as a result of the katabatic wind blocking effect. These offshore-directed winds prevent oceanic near-surface air, influenced by SIC and SST anomalies, from penetrating far inland. Therefore, the ice sheet SMB response is restricted to coastal regions, where katabatic winds cease. A topic for further investigation is how anomalies in SIC and SST might have indirectly affected the surface melt by changing the general circulation in the North Atlantic region, hence favoring more frequent warm air advection towards the GrIS.

Evaluation of the Surface Representation of the Greenland Ice Sheet in a General Circulation Model

Simulated surface conditions of the Goddard Earth Observing System model, version 5 (GEOS-5), atmospheric general circulation model (AGCM) are examined for the contemporary Greenland Ice Sheet (GrIS). A surface parameterization that explicitly models surface processes including snow compaction, meltwater percolation and refreezing, and surface albedo is found to remedy an erroneous deficit in the annual net surface energy flux and provide an adequate representation of surface mass balance (SMB) in an evaluation using simulations at two spatial resolutions. The simulated 1980–2008 GrIS SMB average is 24.7 ± 4.5 cm yr−1 water-equivalent (w.e.) at ½° model grid spacing, and 18.2 ± 3.3 cm yr−1 w.e. for 2° grid spacing. The spatial variability and seasonal cycle of the ½° simulation compare favorably to recent studies using regional climate models, while results from 2° integrations reproduce the primary features of the SMB field. In comparison to historical glaciological observations, the coarser-resolution model overestimates accumulation in the southern areas of the GrIS, while the overall SMB is underestimated. These changes relate to the sensitivity of accumulation and melt to the resolution of topography. The GEOS-5 SMB fields contrast with available corresponding atmospheric models simulations from phase 5 of the Coupled Model Intercomparison Project (CMIP5). It is found that only a few of the CMIP5 AGCMs examined provide significant summertime runoff, a dominant feature of the GrIS seasonal cycle. This is a condition that will need to be remedied if potential contributions to future eustatic change from polar ice sheets are to be examined with GCMs.

Future climate warming increases Greenland ice sheet surface mass balance variability

AbstractThe integrated surface mass balance (SMB) of the Greenland ice sheet (GrIS) has large interannual variability. Long‐term future changes to this variability will affect GrIS dynamics, freshwater fluxes, regional oceanography, and detection of changes in ice volume trends. Here we analyze a simulated 1850–2100 GrIS SMB time series from the Community Earth System Model, currently the only global climate model that realistically simulates GrIS SMB. We find a significant increase in interannual integrated SMB variability over time, which we attribute primarily to a shift to a high‐variability melt‐dominated SMB regime due to GrIS ablation area growth. We find temporal increases to characteristic ablation and accumulation area‐specific SMB variabilities to be of secondary importance. Since ablation area SMB variability is driven largely by variability in summer surface melt, variability in the climate processes regulating the energy fluxes that control melting will likely increasingly determine future GrIS SMB variability.

Estimating the Greenland ice sheet surface mass balance contribution to future sea level rise using the regional atmospheric climate model MAR

Abstract. To estimate the sea level rise (SLR) originating from changes in surface mass balance (SMB) of the Greenland ice sheet (GrIS), we present 21st century climate projections obtained with the regional climate model MAR (Modèle Atmosphérique Régional), forced by output of three CMIP5 (Coupled Model Intercomparison Project Phase 5) general circulation models (GCMs). Our results indicate that in a warmer climate, mass gain from increased winter snowfall over the GrIS does not compensate mass loss through increased meltwater run-off in summer. Despite the large spread in the projected near-surface warming, all the MAR projections show similar non-linear increase of GrIS surface melt volume because no change is projected in the general atmospheric circulation over Greenland. By coarsely estimating the GrIS SMB changes from GCM output, we show that the uncertainty from the GCM-based forcing represents about half of the projected SMB changes. In 2100, the CMIP5 ensemble mean projects a GrIS SMB decrease equivalent to a mean SLR of +4 ± 2 cm and +9 ± 4 cm for the RCP (Representative Concentration Pathways) 4.5 and RCP 8.5 scenarios respectively. These estimates do not consider the positive melt–elevation feedback, although sensitivity experiments using perturbed ice sheet topographies consistent with the projected SMB changes demonstrate that this is a significant feedback, and highlight the importance of coupling regional climate models to an ice sheet model. Such a coupling will allow the assessment of future response of both surface processes and ice-dynamic changes to rising temperatures, as well as their mutual feedbacks.