Ice base slope effects on the turbulent ice shelf-ocean boundary current

Abstract Efforts to parameterize ice shelf basal melting within climate models are limited by an incomplete understanding of the influence of ice base slope on the turbulent ice shelf–ocean boundary current (ISOBC). Here, we examine the relationship between ice base slope, boundary current dynamics, and melt rate using 3D, turbulence-permitting large-eddy simulations (LESs) of an idealized ice shelf–ocean boundary current forced solely by melt-induced buoyancy. The range of simulated slopes (3%–10%) is appropriate to the grounding zone of small Antarctic ice shelves and to the flanks of relatively wide ice base channels, and the initial conditions are representative of warm-cavity ocean conditions. In line with previous studies, the simulations feature the development of an Ekman boundary layer adjacent to the ice, overlaying a broad pycnocline. The time-averaged flow within the pycnocline is in thermal wind balance, with a mean shear that is only weakly dependent on the ice base slope angle α , resulting in a mean gradient Richardson number 〈Ri g 〉 that decreases approximately linearly with sin α . Combining this inverse relationship with a linear approximation to the density profile, we derive formulations for the friction velocity, thermal forcing, and melt rate in terms of slope angle and total buoyancy input. This theory predicts that melt rate varies like the square root of slope, which is consistent with the LES results and differs from a previously proposed linear trend. The derived scalings provide a potential framework for incorporating slope dependence into parameterizations of mixing and melting at the base of ice shelves. Significance Statement The majority of Antarctica’s contribution to sea level rise can be attributed to changes in ocean-driven melting at the base of ice shelves (the floating extensions of the Antarctic ice sheet). Turbulent ocean currents and melting are strongest where the ice base is steeply sloped, but few studies have systematically examined this effect. We use an idealized ice shelf–ocean model to examine how variations in ice base slope influence ocean mixing and ice melting. We derive a formula predicting that melting varies like the square root of the ice base slope, and this scaling is supported by the simulations. These results provide a potential framework for improving the representation of ice shelf melting in climate models.