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.

The calving of icebergs accounts for a significant fraction of the mass loss from both the Antarctic and Greenland ice sheets. Iceberg melting affects the water properties impacting sea ice formation, local circulation and biological activity. Laboratory experiments have investigated the effects of the Earth’s rotation on iceberg melting and the possible formation of Taylor columns (TCs) underneath icebergs. It is found that at high Rossby number, $Ro$, when rotation is weak compared to advection, iceberg melting is unaffected by the background rotation. As $Ro$ decreases, the melt rate shows an increasing trend, which is expected to reverse for very low $Ro$. This behaviour is explained by considering the integrated horizontal velocity at the base of the iceberg. For moderate $Ro$, a partial TC is formed and its effective blocking accelerates the flow under the remainder of the iceberg, which increases the melt rate since the melting is proportional to the flow velocity. It is expected that for very low $Ro$ the melt rate decreases because, with the expansion of the TC, the region of flow acceleration occurs away from the base of the iceberg. For low free stream velocity the freshwater produced by the ice melting introduces another dynamical effect. It is observed that there is a threshold free stream velocity below which the melt rate is constant. This is explained with the formation of a gravity current at the base of the iceberg that insulates it from the free flow and controls its melting.

The propagation of full-depth lock-exchange bottom gravity currents past a submerged array of circular cylinders is investigated using laboratory experiments and large eddy simulations. Firstly, to investigate the front velocity of gravity currents across the whole range of array density $\unicode[STIX]{x1D719}$ (i.e. the volume fraction of solids), the array is densified from a flat bed ($\unicode[STIX]{x1D719}=0$) towards a solid slab ($\unicode[STIX]{x1D719}=1$) under a particular submergence ratio $H/h$, where $H$ is the flow depth and $h$ is the array height. The time-averaged front velocity in the slumping phase of the gravity current is found to first decrease and then increase with increasing $\unicode[STIX]{x1D719}$. Next, a new geometrical framework consisting of a streamwise array density $\unicode[STIX]{x1D707}_{x}=d/s_{x}$ and a spanwise array density $\unicode[STIX]{x1D707}_{y}=d/s_{y}$ is proposed to account for organized but non-equidistant arrays ($\unicode[STIX]{x1D707}_{x}\neq \unicode[STIX]{x1D707}_{y}$), where $s_{x}$ and $s_{y}$ are the streamwise and spanwise cylinder spacings, respectively, and $d$ is the cylinder diameter. It is argued that this two-dimensional parameter space can provide a more quantitative and unambiguous description of the current–array interaction compared with the array density given by $\unicode[STIX]{x1D719}=(\unicode[STIX]{x03C0}/4)\unicode[STIX]{x1D707}_{x}\unicode[STIX]{x1D707}_{y}$. Both in-line and staggered arrays are investigated. Four dynamically different flow regimes are identified: (i) through-flow propagating in the array interior subject to individual cylinder wakes ($\unicode[STIX]{x1D707}_{x}$: small for in-line array and arbitrary for staggered array; $\unicode[STIX]{x1D707}_{y}$: small); (ii) over-flow propagating on the top of the array subject to vertical convective instability ($\unicode[STIX]{x1D707}_{x}$: large; $\unicode[STIX]{x1D707}_{y}$: large); (iii) plunging-flow climbing sparse close-to-impermeable rows of cylinders with minor streamwise intrusion ($\unicode[STIX]{x1D707}_{x}$: small; $\unicode[STIX]{x1D707}_{y}$: large); and (iv) skimming-flow channelized by an in-line array into several subcurrents with strong wake sheltering ($\unicode[STIX]{x1D707}_{x}$: large; $\unicode[STIX]{x1D707}_{y}$: small). The most remarkable difference between in-line and staggered arrays is the non-existence of skimming-flow in the latter due to the flow interruption by the offset rows. Our analysis reveals that as $\unicode[STIX]{x1D719}$ increases, the change of flow regime from through-flow towards over- or skimming-flow is responsible for increasing the gravity current front velocity.

The advance of the front of a dense gravity current propagating in a rectangular channel and V-shaped valley both horizontally and up a shallow slope is examined through theory, full-depth lock–release laboratory experiments and hydrostatic numerical simulations. Consistent with theory, experiments and simulations show that the front speed is relatively faster in the valley than in the channel. The front speed measured shortly after release from the lock is 5–22 % smaller than theory, with greater discrepancy found in upsloping V-shaped valleys. By contrast, the simulated speed is approximately 6 % larger than theory, showing no dependence on slope for rise angles up to ${\it\theta}=8^{\circ }$. Unlike gravity currents in a channel, the current head is observed in experiments to be more turbulent when propagating in a V-shaped valley. The turbulence is presumably enhanced due to the lateral flows down the sloping sides of the valley. As a consequence, lateral momentum transport contributes to the observed lower initial speeds. A Wentzel–Kramers–Brillouin like theory predicting the deceleration of the current as it runs upslope agrees remarkably well with simulations and with most experiments, within errors.

We discuss laboratory experiments investigating mixing in a density-driven current flowing down a sloping bottom, in a rotating homogenous fluid. A systematic study spanning a wide range of Froude, 0.8 < Fr < 10, and Reynolds, 10 < Re < 1400, numbers was conducted by varying three parameters: the bottom slope; the flow rate; and the density of the dense fluid. Different flow regimes were observed, i.e. waves (non-breaking and breaking) and turbulent regimes, while changing the above parameters. Mixing in the density-driven current has been quantified within the observed regimes, and at different locations on the slope. The dependence of mixing on the relevant non-dimensional numbers, i.e. slope, Fr and Re, is discussed. The entrainment parameter, E, was found to be dependent not only on Fr, as assumed in previous studies, but also on Re. In particular, mixing increased with increasing Fr and Re. For low Fr and Re, the magnitude of the mixing was comparable to mixing in the ocean. For large Fr and Re, mixing was comparable to that observed in previous laboratory experiments that exhibited the classic turbulent entrainment behaviour.