Fracture-induced softening for large-scale ice dynamics

Abstract

Abstract. Floating ice shelves can exert a retentive and hence stabilizing force onto the inland ice sheet of Antarctica. However, this effect has been observed to diminish by fracture-coupled dynamic processes within the protective ice shelves leading to accelerated ice flow and hence to a sea-level contribution. In order to better understand the role of fractures in ice dynamics we apply a large-scale continuum representation of fractures and related fracture growth into the prognostic Parallel Ice Sheet Model (PISM). To this end we introduce a higher-order accuracy advection scheme for the transport of the two-dimensional fracture density across the regular computational grid. Dynamic coupling of fractures and ice flow is attained by a reduction of effective ice viscosity proportional to the inferred fracture density. This formulation implies the possibility of a non-linear threshold behavior due to self-amplified fracturing in shear regions triggered by small variations in damage threshold. As a result of prognostic flow simulations, flow patterns with realistically large across-flow velocity gradients in fracture-weakened regions as seen in observations are reproduced. This model framework is expandable to grounded ice streams and accounts for climate-induced effects on fracturing and hence on the ice-flow dynamics. It further allows for an enhanced fracture-based calving parameterization.