We present a theoretical and experimental study of viscous flows injected into a porous medium that is confined vertically by horizontal impermeable boundaries and filled with an ambient fluid of different density and viscosity. General three-dimensional equations describing such flows are developed, showing that the dynamics can be affected by two separate contributions: spreading due to gradients in hydrostatic pressure, and that due to the pressure drop introduced by the injection. In the illustrative case of a two-dimensional injection of fluid at a constant volumetric rate, the injected fluid initially forms a viscous gravity current insensitive both to the depth of the medium and to the viscosity of the ambient fluid. Beyond a characteristic time scale, the dynamics transition to being dominated by the injection pressure, and the injected fluid eventually intersects the second boundary to form a second moving contact line. Three different late-time asymptotic regimes can emerge, depending on whether the viscosity of the injected fluid is less than, equal to or greater than that of the ambient fluid. With a less viscous injection, the flow undergoes a slow decay towards a similarity solution in which the two contact lines extend linearly in time with differing prefactors. Perturbations from this long-term state are shown to decay algebraically with time. Equal viscosities result in both contact lines approaching the same leading-order asymptotic position but with a first-order correction to the distance between them that expands as t 1/2 due to gravitational spreading. For a more viscous injection, the distance between the contact lines approaches a constant value, with perturbations decaying exponentially. Data from a new series of laboratory experiments confirm these theoretical predictions.

I present a theoretical and experimental study of floating viscous fluid films introduced into a channel of finite length, motivated by the flow of glacial ice shelves. The dynamics are characterized by a mixture of viscous extensional stresses, transverse shear stresses and a driving buoyancy force. A theory based on a width-integrated model is developed and investigated using analytical, asymptotic and numerical methods. With fluid introduced at a constant rate, the flow is found to approach a steady state with two possible asymptotic forms depending on the length of the channel. For channel lengths less than half the width, the flow is similar to a purely extensional one-dimensional flow, characterized by concave surface profiles and being insensitive to the position of the channel exit (or calving front). Greater lengths result in a more complex asymptotic structure in which the flow adjusts over a short distance towards a prevailing flow of universal dimensionless form. In complete contrast to the extensional regime, the prevailing flow is controlled by the position of the channel exit. Data from a new laboratory experiment involving particle velocimetry of a floating fluid film compares well with the predicted along-channel velocity. Motivated by glaciological application, the analysis is generalized to power-law rheologies and the results used to classify the flow regimes of a selection of ice shelves. The prediction for the frontal speed is in good agreement with geophysical data, indicating that the universal profile predicted by the theory is common in nature.

We present a theoretical and experimental study of viscous gravity currents introduced at the surface of a denser inviscid fluid layer of finite depth inside a vertical HeleShaw cell. Initially, the viscous fluid floats on the inviscid fluid, forming a self-similar, buoyancy-driven current resisted predominantly by the viscous stresses due to shear across the width of the cell. Once the viscous current contacts the base of the cell, the flow can be considered in two regions: a grounded region in which the current lies in full contact with the base; and a floating region. The subsequent advance of the grounding line separating these regions is shown to be controlled by the thickening of the current associated with balancing the local shear stresses. An understanding of the flow transitions is developed using asymptotic and numerical analysis of a model based on lubrication theory.

Abstract

A long-standing open question in glaciology concerns the propensity for ice sheets that lie predominantly submerged in the ocean (marine ice sheets) to destabilise under buoyancy. This paper addresses the processes by which a buoyancy-driven mechanism for the retreat and ultimate collapse of such ice sheets – the marine ice sheet instability – is suppressed by lateral stresses acting on its floating component (the ice shelf). The key results are to demonstrate the transition between a mode of stable (easily reversible) retreat along a stable steady-state branch created by ice-shelf buttressing to tipped (almost irreversible) retreat across a critical parametric threshold. The conditions for triggering tipped retreat can be controlled by the calving position and other properties of the ice-shelf profile and can be largely independent of basal stress, in contrast to principles established from studies of unbuttressed grounding-line dynamics. The stability and recovery conditions introduced by lateral stresses are analysed by developing a method of constructing grounding-line stability (bifurcation) diagrams, which provide a rapid assessment of the steady-state positions, their natures and the conditions for secondary grounding, giving clear visualisations of global stabilisation conditions. A further result is to reveal the possibility of a third structural component of a marine ice sheet that lies intermediate to the fully grounded and floating components. The region forms an extended grounding area in which the ice sheet lies very close to flotation, and there is no clearly distinguished grounding line. The formation of this region generates an upsurge in buttressing that provides the most feasible mechanism for reversal of a tipped grounding line. The results of this paper provide conceptual insight into the phenomena controlling the stability of the West Antarctic Ice Sheet, the collapse of which has the potential to dominate future contributions to global sea-level rise.