Search

Geologic archives constraining the variability of the Greenland ice sheet (GrIS) during the Holocene provide targets for ice sheet models to test sensitivities to variations in past climate and model formulation. Even as data–model comparisons are becoming more common, many models simulating the behavior of the GrIS during the past rely on meshes with coarse horizontal resolutions (≥10 km). In this study, we explore the impact of model resolution on the simulated nature of retreat across southwestern Greenland during the Holocene. Four simulations are performed using the Ice Sheet System Model (ISSM): three that use a uniform mesh and horizontal mesh resolutions of 20, 10, and 5 km, and one that uses a nonuniform mesh with a resolution ranging from 2 to 15 km. We find that the simulated retreat can vary significantly between models with different horizontal resolutions based on how well the bed topography is resolved. In areas of low topographic relief, the horizontal resolution plays a negligible role in simulated differences in retreat, with each model instead responding similarly to retreat driven by surface mass balance (SMB). Conversely, in areas where the bed topography is complex and high in relief, such as fjords, the lower-resolution models (10 and 20 km) simulate unrealistic retreat that occurs as ice surface lowering intersects bumps in the bed topography that would otherwise be resolved as troughs using the higher-resolution grids. Our results highlight the important role that high-resolution grids play in simulating retreat in areas of complex bed topography, but also suggest that models using nonuniform grids can save computational resources through coarsening the mesh in areas of noncomplex bed topography where the SMB predominantly drives retreat. Additionally, these results emphasize that care must be taken with ice sheet models when tuning model parameters to match reconstructed margins, particularly for lower-resolution models in regions where complex bed topography is poorly resolved. [ABSTRACT FROM AUTHOR]

Quantifying the Greenland Ice Sheet's future contribution to sea level rise is a challenging task that requires accurate estimates of ice sheet sensitivity to climate change. Forward ice sheet models are promising tools for estimating future ice sheet behavior, yet confidence is low because evaluation of historical simulations is challenging due to the scarcity of continental-wide data for model evaluation. Recent advancements in processing of Gravity Recovery and Climate Experiment (GRACE) data using Bayesianconstrained mass concentration ("mascon") functions have led to improvements in spatial resolution and noise reduction of monthly global gravity fields. Specifically, the Jet Propulsion Laboratory's JPL RL05M GRACE mascon solution (GRACE_JPL) offers an opportunity for the assessment of model-based estimates of ice sheet mass balance (MB) at ~300 km spatial scales. Here, we quantify the differences between Greenland monthly observed MB (GRACE_JPL) and that estimated by state-of-the-art, high-resolution models, with respect to GRACE_JPL and model uncertainties. To simulate the years 2003-2012, we force the Ice Sheet System Model (ISSM) with anomalies from three different surface mass balance (SMB) products derived from regional climate models. Resulting MB is compared against GRACE_JPL within individual mascons. Overall, we find agreement in the northeast and southwest where MB is assumed to be primarily controlled by SMB. In the interior, we find a discrepancy in trend, which we presume to be related to millennial-scale dynamic thickening not considered by our model. In the northwest, seasonal amplitudes agree, but modeled mass trends are muted relative to GRACE_JPL. Here, discrepancies are likely controlled by temporal variability in ice discharge and other related processes not represented by our model simulations, i.e., hydrological processes and ice-ocean interaction. In the southeast, GRACE_JPL exhibits larger seasonal amplitude than predicted by the models while simultaneously having more pronounced trends; thus, discrepancies are likely controlled by a combination of missing processes and errors in both the SMB products and ISSM. At the margins, we find evidence of consistent intra-annual variations in regional MB that deviate distinctively from the SMB annual cycle. Ultimately, these monthly-scale variations, likely associated with hydrology or ice-ocean interaction, contribute to steeper negative mass trends observed by GRACE_JPL. Thus, models should consider such processes at relatively high (monthly-to-seasonal) temporal resolutions to achieve accurate estimates of Greenland MB. [ABSTRACT FROM AUTHOR]

Calving is a major mechanism of ice discharge of the Antarctic and Greenland ice sheets, and a change in calving front position affects the entire stress regime of marine terminating glaciers. The representation of calving front dynamics in a 2-D or 3-D ice sheet model remains non-trivial. Here, we present the theoretical and technical framework for a level-set method, an implicit boundary tracking scheme, which we implement into the Ice Sheet System Model (ISSM). This scheme allows us to study the dynamic response of a drainage basin to user-defined calving rates.We apply the method to Jakobshavn Isbræ, a major marine terminating outlet glacier of the West Greenland Ice Sheet. The model robustly reproduces the high sensitivity of the glacier to calving, and we find that enhanced calving triggers significant acceleration of the ice stream. Upstream acceleration is sustained through a combination of mechanisms. However, both lateral stress and ice influx stabilize the ice stream. This study provides new insights into the ongoing changes occurring at Jakobshavn Isbræ and emphasizes that the incorporation of moving boundaries and dynamic lateral effects, not captured in flow-line models, is key for realistic model projections of sea level rise on centennial timescales. [ABSTRACT FROM AUTHOR]