Abstract. Basal melt at the bottom of Antarctic ice shelves is a major control on glacier dynamics, as it modulates the amount of buttressing that floating ice shelves exert onto the ice streams feeding them. Three-dimensional ocean circulation numerical models provide reliable estimates of basal melt rates but remain too computationally expensive for century scale projections. Ice sheet modelers therefore routinely rely on simplified parameterizations either based on ice shelf depth or on more sophisticated box models. However, existing parameterizations do not accurately resolve the complex spatial patterns of sub-shelf melt rates that have been observed over Antarctica's ice shelves, especially in the vicinity of the grounding line, where basal melt is one of the primary drivers of grounding line migration. In this study, we couple the Potsdam Ice-shelf Cavity mOdel (PICO) to a buoyant Plume melt rate parameterization to create PICOP, a novel basal melt rate parameterization that is easy to implement in transient ice sheet numerical models and produces a melt rate field that is in excellent agreement with the spatial distribution and magnitude of observations for a wide variety of ocean basins. We test PICOP on the Amundsen Sea sector of West Antarctica, Totten and Moscow University ice shelves in Eastern Antarctica, and the Ronne-Filchner ice shelf and compare the results to PICO. We find that PICOP is able to reproduce the high melt rates near the grounding lines of Pine Island, Thwaites, and Totten glaciers (on the order of 100 m/yr) and removes the “banding” pattern observed in melt rates produced by PICO over the Ronne-Filchner ice shelf. PICOP resolves many of the issues contemporary basal melt rate parameterizations face and is therefore a valuable tool for those looking to make future projections of Antarctic glaciers.
Brief communication: PICOP, a new ocean melt parameterization under ice shelves combining PICO and a plume model