Search for: All records

Abstract Ice shelves regulate ice sheet dynamics, with their stability influenced by horizontal flow and vertical flexure. MacAyeal and others (2021) developed the theoretical foundation for a coupled flow-flexure model (the “M21 model”), combining the Shallow Shelf Approximation with thin-beam flexure, providing a computationally efficient tool for studying phenomena like ice shelf rumpling and lake drainage. However, the M21 model relies on proprietary software, is unstable under compressive flow conditions, and does not incorporate fracture processes critical for capturing ice-shelf damage evolution. We present an open-source version of the M21 model addressing these limitations. Using the free Python librariesFiredrakeandicepack, we introduce a plastic failure mechanism, effectively limiting bending stresses and thereby stabilizing the model. This enhancement expands the viscous M21 model into a viscoplastic flow-flexure-fracture (3F) framework. We validate the 3F model through test cases replicating key ice shelf phenomena, including marginal rumpling and periodic surface meltwater drainage. By offering this tool as open-source software, we aim to enable broader adoption, with the ultimate aim of representing surface meltwater induced flow-flexure-fracture processes in large-scale ice sheet models. 

Abstract Global Navigation Satellite System (GNSS) observations and ground-based timelapse photography obtained over the record-high 2019/2020 melt season are combined to characterise the flexure and fracture behaviour of a previously formed doline on George VI Ice Shelf, Antarctica. The GNSS timeseries shows a downward vertical displacement of the doline centre with respect to the doline rim of ~60 cm in response to loading from a central meltwater lake. The GNSS data also show a tens-of-days episode of rapid-onset, exponentially decaying horizontal displacement, where the horizontal distance between the doline rim and its centre increases by ~70 cm. We interpret this event as the initiation and/or widening of a fracture, aided by stress perturbations associated with meltwater loading in the doline basin. Viscous flexure modelling indicates that the meltwater loading generates tensile surface stresses exceeding 75 kPa. This, together with our timelapse photos of circular fractures around the doline, suggests the first such documentation of meltwater-loading-induced ‘ring fracture’ formation on an ice shelf, equivalent to the fracture type proposed as part of the chain-reaction lake drainage process involved in the 2002 breakup of the Larsen B Ice Shelf.