Ice shelf fracture parameterization in an ice sheet model
Abstract. Floating ice shelves exert a stabilizing force onto the inland ice sheet. However, this buttressing effect is diminished by the fracture process, which on large scales effectively softens the ice, accelerating its flow, increasing calving, and potentially leading to ice shelf breakup. Here, we explore how the application of a continuum damage model (CDM) to the prognostic ice sheet model BISICLES can account for the effects of fracture processes on viscous ice dynamics. Damage is created by the local stress field and advects downstream. This continuum damage model is coupled to the dynamical ice flow model by decreasing the effective viscosity proportional to the damage field. To evaluate the physical role of the fracture process on large-scale ice sheet dynamics and also discern the relative importance of the parameters used in the damage model, we carry out a suite of numerical experiments based on the MISMIP+ (Marine Ice Sheet Model Intercomparison Project) marine ice sheet geometry. We find that behavior of the simulated marine ice sheet is sensitive to fracture processes on the ice shelf. In the case of a geometry that produces strong lateral stress, the stiffness of ice around the grounding line is essential to ice sheet evolution, with softer or more damaged ice leading to thinning and grounding line retreat.