skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Simulated Greenland Surface Mass Balance in the GISS ModelE2 GCM: Role of the Ice Sheet Surface
Abstract The 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.  more » « less
Award ID(s):
1713072
PAR ID:
10375147
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Earth Surface
Volume:
124
Issue:
3
ISSN:
2169-9003
Page Range / eLocation ID:
p. 750-765
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract. Surface mass loss from the Greenland ice sheet (GrIS) hasaccelerated over the past decades, mainly due to enhanced surface meltingand liquid water runoff in response to atmospheric warming. A large portionof runoff from the GrIS originates from exposure of the darker bare ice inthe ablation zone when the overlying snow melts, where surface albedo playsa critical role in modulating the energy available for melting. In thisregard, it is imperative to understand the processes governing albedovariability to accurately project future mass loss from the GrIS. Bare-icealbedo is spatially and temporally variable and contingent on non-linearfeedbacks and the presence of light-absorbing constituents. An assessment ofmodels aiming at simulating albedo variability and associated impacts onmeltwater production is crucial for improving our understanding of theprocesses governing these feedbacks and, in turn, surface mass loss fromGreenland. Here, we report the results of a comparison of the bare-iceextent and albedo simulated by the regional climate model ModèleAtmosphérique Régional (MAR) with satellite imagery from theModerate Resolution Imaging Spectroradiometer (MODIS) for the GrIS below70∘ N. Our findings suggest that MAR overestimates bare-ice albedoby 22.8 % on average in this area during the 2000–2021 period with respectto the estimates obtained from MODIS. Using an energy balance model toparameterize meltwater production, we find this bare-ice albedo bias canlead to an underestimation of total meltwater production from the bare-icezone below 70∘ N of 42.8 % during the summers of 2000–2021. 
    more » « less
  2. Modèle Atmosphérique Régional (MAR) is a regional climate model that is fully coupled to a one-dimensional surface-atmosphere energy and mass transfer scheme, SISVAT (Soil Ice Snow Vegetation Atmosphere Transfer) (Fettweis et al., 2005, 2020; Lefebre et al., 2005). SISVAT employs a multilayered snowpack model, CROCUS, that simulates meltwater production, percolation, and refreeze (Brun et al., 1989), while also accounting for changes in albedo due to snow metamorphism (Brun et al., 1992). MAR has been extensively verified over the Greenland Ice Sheet and is therefore particularly well suited for analyses of Greenland ice sheet surface mass balance (Fettweis et al., 2011; Fettweis et al., 2020; Lefebre et al. 2005; Mattingly et al. 2020). Brun, E., Martin, E., Simon, V., Gendre, C., and Coléou, C. (1989). An energy and mass model of snow cover suitable for operational avalanche forecasting. Journal of Glaciology, 35, 333. https://doi.org/10.1017/S0022143000009254 Brun, E., David, P., Sudul, M., and Brunot, G. (1992). A numerical model to simulate snow-cover stratigraphy for operational avalanche forecasting. Journal of Glaciology, 38(128), 13–22. https://doi.org/10.3189/S0022143000009552 Fettweis, X., Gallée, H., Lefebre, F., and van Ypersele, J.-P. (2005). Greenland surface mass balance simulated by a regional climate model and comparison with satellite-derived data in 1990–1991. Climate Dynamics, 24(6), 623–640. https://doi.org/10.1007/s00382-005-0010-y Fettweis, X., Tedesco, M., van den Broeke, M., and Ettema, J. (2011). Melting trends over the Greenland ice sheet (1958–2009) from spaceborne microwave data and regional climate models. The Cryosphere, 5(2), 359–375. https://doi.org/10.5194/tc-5-359-2011 Fettweis, X., Hofer, S., Krebs-Kanzow, U., Amory, C., Aoki, T., Berends, C. J., et al. (2020). GrSMBMIP: intercomparison of the modelled 1980–2012 surface mass balance over the Greenland Ice Sheet. The Cryosphere, 14(11), 3935–3958. https://doi.org/10.5194/tc-14-3935-2020 Lefebre, F., Fettweis, X., Gallée, H., Van Ypersele, J.-P., Marbaix, P., Greuell, W., and Calanca, P. (2005). Evaluation of a high-resolution regional climate simulation over Greenland. Climate Dynamics, 25(1), 99–116. https://doi.org/10.1007/s00382-005-0005-8 Mattingly, K. S., Mote, T. L., Fettweis, X., van As, D., Van Tricht, K., Lhermitte, S., et al. (2020). Strong summer atmospheric rivers trigger Greenland ice sheet melt through spatially varying surface energy balance and cloud regimes. Journal of Climate, 33(16), 6809–6832. https://doi.org/10.1175/JCLI-D-19-0835.1 
    more » « less
  3. The simulation of ice sheet-climate interaction such as surface massbalance fluxes are sensitive to model grid resolution. Here we simulatethe multicentury evolution of the Greenland Ice Sheet (GrIS) and itsinteraction with the climate using the Community Earth System Modelversion 2.2 (CESM2.2) including an interactive GrIS component (theCommunity Ice Sheet Model v2.1 [CISM2.1]) under an idealized warmingscenario (atmospheric CO2 increases by 1% yr−1 until quadrupling thepre-industrial level and then is held fixed). A variable-resolution (VR)grid with 1/4◦ regional refinement over broader Arctic and 1◦ resolutionelsewhere is applied to the atmosphere and land components, and theresults are compared to conventional 1◦ lat-lon grid simulations toinvestigate the impact of grid refinement. An acceleration of GrIS massloss is found at around year 110, caused by rapidly increasing surfacemelt as the ablation area expands with associated albedo feedback andincreased turbulent fluxes. Compared to the 1◦ runs, the VR run featuresslower melt increase, especially over Western and Northern Greenland,which slope gently towards the peripheries. This difference patternoriginates primarily from the weaker albedo feedback in the VR run,complemented by its smaller cloud longwave radiation. The steeper VRGreenland surface topography favors slower ablation zone expansion, thusleading to its weaker albedo feedback. The sea level rise contributionfrom the GrIS in the VR run is 53 mm by year 150 and 831 mm by year 350,approximately 40% and 20% smaller than the 1◦ runs, respectively. 
    more » « less
  4. Abstract. Snowfall is the major source of mass for the Greenland ice sheet (GrIS) but the spatial and temporalvariability of snowfall and the connections between snowfall and mass balance have so far been inadequatelyquantified. By characterizing local atmospheric circulation and utilizing CloudSat spaceborne radarobservations of snowfall, we provide a detailed spatial analysis of snowfall variability and its relationshipto Greenland mass balance, presenting first-of-their-kind maps of daily spatial variability in snowfallfrom observations across Greenland. For identified regional atmospheric circulation patterns, we show that thespatial distribution and net mass input of snowfall vary significantly with the position and strength ofsurface cyclones. Cyclones west of Greenland driving southerly flow contribute significantly more snowfall thanany other circulation regime, with each daily occurrence of the most extreme southerly circulation patterncontributing an average of 1.66 Gt of snow to the Greenland ice sheet. While cyclones east of Greenland,patterns with the least snowfall, contribute as little as 0.58 Gt each day. Above 2 km on the ice sheet wheresnowfall is inconsistent, extreme southerly patterns are the most significant mass contributors, with up to1.20 Gt of snowfall above this elevation. This analysis demonstrates that snowfall over the interior ofGreenland varies by up to a factor of 5 depending on regional circulation conditions. Using independentobservations of mass changes made by the Gravity Recovery and Climate Experiment (GRACE), we verify that thelargest mass increases are tied to the southerly regime with cyclones west of Greenland. For occurrences of thestrongest southerly pattern, GRACE indicates a net mass increase of 1.29 Gt in the ice sheet accumulation zone(above 2 km elevation) compared to the 1.20 Gt of snowfall observed by CloudSat. This overall agreementsuggests that the analytical approach presented here can be used to directly quantify snowfall masscontributions and their most significant drivers spatially across the GrIS. While previous research hasimplicated this same southerly regime in ablation processes during summer, this paper shows that ablation massloss in this circulation regime is nearly an order of magnitude larger than the mass gain from associatedsnowfall. For daily occurrences of the southerly circulation regime, a mass loss of approximately 11 Gt isobserved across the ice sheet despite snowfall mass input exceeding 1 Gt. By analyzing the spatialvariability of snowfall and mass changes, this research provides new insight into connections between regionalatmospheric circulation and GrIS mass balance. 
    more » « less
  5. Abstract. Quantifying the total liquid water amounts (LWAs) in the Greenland ice sheet (GrIS) is critical for understanding GrIS firn processes, mass balance, and global sea level rise. Although satellite microwave observations are very sensitive to ice sheet melt and thus can provide a way of monitoring the ice sheet melt globally, estimating total LWA, especially the subsurface LWA, remains a challenge. Here, we present a microwave retrieval of LWA over Greenland using enhanced-resolution L-band brightness temperature (TB) data products from the National Aeronautics and Space Administration (NASA) Soil Moisture Active Passive (SMAP) satellite for the 2015–2023 period. L-band signals receive emission contributions deep in the ice sheet and are sensitive to the liquid water content (LWC) in the firn column. Therefore, they can estimate the surface-to-subsurface LWA, unlike higher-frequency signals (e.g., 18 and 37 GHz bands), which are limited to the top few centimeters of the surface snow during the melt. We used vertically polarized TB (V-pol TB) with empirically derived thresholds to detect liquid water and identify distinct ice sheet zones. A forward model based on radiative transfer (RT) in the ice sheet was used to simulate TB. The simulated TB was then used in an inversion algorithm to estimate LWA. Finally, the retrievals were compared with the LWA obtained from two sources. The first source was a locally calibrated ice sheet energy and mass balance (EMB) model, and the second source was the Glacier Energy and Mass Balance (GEMB) model within NASA's Ice-sheet and Sea-level System Model (ISSM). Both models were forced by in situ measurements from six automatic weather stations (AWSs) of the Programme for Monitoring of the Greenland Ice Sheet (PROMICE) and the Greenland Climate Network (GC-Net) located in the percolation zone of the GrIS. The retrievals show generally good agreement with both the references, demonstrating the potential for advancing our understanding of ice sheet physical processes to better project Greenland's contribution to the global sea level rise in response to the warming climate. 
    more » « less