Simulating the Greenland Ice Sheet Surface Mass Balance using an Earth System Model

• Main Author: van Kampenhout, Leonardus
• Other Authors: Broeke, M.R. van den, Lenaerts, J.T.M.
• Format: Doctoral or Postdoctoral Thesis
• Language: English
• Published: Utrecht University 2020
• Subjects: Groenland; klimaatmodel; zeespiegelstijging; ijskap; modelontwikkeling; Greenland; Ice Sheet
• Online Access: https://dspace.library.uu.nl/handle/1874/400695

Observational evidence shows that the Greenland Ice Sheet (GrIS) is currently losing mass, and that this mass loss is accelerating. In this dissertation, we look at the representation of GrIS climate in the Community Earth System Model (CESM), in which ice sheet interactions with other parts of the climate system can be consistently resolved. The central metric is the surface mass balance (SMB), as a realistic representation of SMB is a prerequisite for gaining confidence in ice sheet mass loss projections. Our first study investigates the representation of polar snow and firn in CESM. We introduce new parametrizations for fresh snow density and snow compaction and increase the number of snow layers in the model. These model improvements largely resolve issues identified by previous work, including an insufficient refreezing capacity and spurious subsurface melt rates. Our second study concerns spatial resolution. Previous work on regional models has indicated that ~20 km resolution is needed to resolve steep SMB gradients near the GrIS margins, whereas CESM has a nominal resolution of 1° (~100 km). We apply, for the first time, regional grid refinement (55 and 28 km) to the Greenland region using Variable Resolution CESM. We find that although precipitation is greatly improved with enhanced resolution, surface runoff is not. Our third study evaluates the present-day GrIS climate in CESM version 2 (CESM2). We find that CESM2 overestimates the net radiative input at the surface in summer, primarily driven by a low biased albedo, which is compensating for a lack in sensible turbulent heat flux. GrIS ablation zones do not always overlap with that in RACMO2, but the difference in total ablation area is small (-11%). Integrated amounts of SMB components melt, refreezing and runoff are all bracketed by RACMO2 equivalents at the resolutions 11 km (native) and 1 km (statistically down- scaled), indicating a good performance. Our last study explores the future timing of GrIS SMB turning from positive to negative. If such ...