Andrew Wells

Abstract:

Thermodynamic ablation of ice in contact with the ocean is an essential element of ice sheet and ocean interactions but is challenging to model and quantify. Building on earlier observations of sea ice ablation, a variety of recent theoretical, experimental and observational studies have considered ice ablation in contrasting geometries, from vertical to near-horizontal ice faces, and reveal different scaling behaviour for predicted ablation rates in different dynamical regimes. However, uncertainties remain about when the contrasting results should be applied, as existing model parameterisations do not capture all relevant regimes of ice–ocean ablation. To progress towards improved models of ice–ocean​ interaction, we synthesise current understanding into a classification of ablation types. We examine the effect of the classification on the parameterisation of turbulent fluxes from the ocean towards the ice, and identify the dominant processes next to ice interfaces of different orientation. Four ablation types are defined: melting and dissolving based on ocean temperatures, and shear-controlled and buoyancy-controlled regimes based on the dynamics of the near-ice molecular sublayer. We describe existing observational and modelling studies of sea ice, ice shelves, and glacier termini, as well as laboratory studies, to show how they fit into this classification. Two sets of observations from the Ross and Ronne Ice Shelf cavities suggest that both the buoyancy-controlled and shear-controlled regimes may be relevant under different oceanographic conditions. Overall, buoyancy-controlled dynamics are more likely when the molecular sublayer has lower Reynolds number, and shear for higher Reynolds number, although the observations suggest some variability about this trend.

Abstract:

Marine-terminating glaciers, such as those along the coastline of Greenland, often release meltwater into the ocean in the form of subglacial discharge plumes. Though these plumes can dramatically alter the mass loss along the front of a glacier, the conditions surrounding their genesis remain poorly constrained. In particular, little is known about the geometry of subglacial outlets and the extent to which seawater may intrude into them. Here, the latter is addressed by exploring the dynamics of an arrested salt wedge – a steady-state, two-layer flow system where salty water partially intrudes a channel carrying fresh water. Building on existing theory, we formulate a model that predicts the length of a non-entraining salt wedge as a function of the Froude number, the slope of the channel and coefficients for interfacial and wall drag. In conjunction, a series of laboratory experiments were conducted to observe a salt wedge within a rectangular channel. For experiments conducted with laminar flow (Reynolds number Re < 800), good agreement with theoretical predictions are obtained when the drag coefficients are modelled as being inversely proportional to Re. However, for fully turbulent flows on geophysical scales, these drag coefficients are expected to asymptote toward finite values. Adopting reasonable drag coefficient estimates for this flow regime, our theoretical model suggests that typical subglacial channels may permit seawater intrusions of the order of several kilometres. While crude, these results indicate that the ocean has a strong tendency to penetrate subglacial channels and potentially undercut the face of marine-terminating glaciers.