- Articles & preprints
- Submission
- Policies
- Peer review
- Editorial board
- About
- EGU publications
- Manuscript tracking

Journal cover
Journal topic
**The Cryosphere**
An interactive open-access journal of the European Geosciences Union

Journal topic

- Articles & preprints
- Submission
- Policies
- Peer review
- Editorial board
- About
- EGU publications
- Manuscript tracking

- Articles & preprints
- Submission
- Policies
- Peer review
- Editorial board
- About
- EGU publications
- Manuscript tracking

**Research article**
08 Jan 2018

**Research article** | 08 Jan 2018

Frazil-ice growth rate and dynamics in mixed layers and sub-ice-shelf plumes

^{1}Atmospheric, Oceanic and Planetary Physics, Department of Physics, University of Oxford, Clarendon Laboratory, Parks Road, Oxford, OX1 3PU, UK^{2}Department of Earth Sciences, University of Oxford, South Parks Road, Oxford, OX1 3AN, UK

^{1}Atmospheric, Oceanic and Planetary Physics, Department of Physics, University of Oxford, Clarendon Laboratory, Parks Road, Oxford, OX1 3PU, UK^{2}Department of Earth Sciences, University of Oxford, South Parks Road, Oxford, OX1 3AN, UK

**Correspondence**: David W. Rees Jones (david.reesjones@earth.ox.ac.uk)

**Correspondence**: David W. Rees Jones (david.reesjones@earth.ox.ac.uk)

Abstract

Back to toptop
The growth of frazil or granular ice is an important mode of ice formation in the cryosphere. Recent advances have improved our understanding of the microphysical processes that control the rate of ice-crystal growth when water is cooled beneath its freezing temperature. These advances suggest that crystals grow much faster than previously thought. In this paper, we consider models of a population of ice crystals with different sizes to provide insight into the treatment of frazil ice in large-scale models. We consider the role of crystal growth alongside the other physical processes that determine the dynamics of frazil ice. We apply our model to a simple mixed layer (such as at the surface of the ocean) and to a buoyant plume under a floating ice shelf. We provide numerical calculations and scaling arguments to predict the occurrence of frazil-ice explosions, which we show are controlled by crystal growth, nucleation, and gravitational removal. Faster crystal growth, higher secondary nucleation, and slower gravitational removal make frazil-ice explosions more likely. We identify steady-state crystal size distributions, which are largely insensitive to crystal growth rate but are affected by the relative importance of secondary nucleation to gravitational removal. Finally, we show that the fate of plumes underneath ice shelves is dramatically affected by frazil-ice dynamics. Differences in the parameterization of crystal growth and nucleation give rise to radically different predictions of basal accretion and plume dynamics, and can even impact whether a plume reaches the end of the ice shelf or intrudes at depth.

Download & links

How to cite

Back to top
top
How to cite.

Rees Jones, D. W. and Wells, A. J.: Frazil-ice growth rate and dynamics in mixed layers and sub-ice-shelf plumes, The Cryosphere, 12, 25–38, https://doi.org/10.5194/tc-12-25-2018, 2018.

1 Introduction

Back to toptop
Frazil-ice formation is an extremely rapid mode of ice growth occurring as the initial phase of ice growth in turbulent waters. Frazil ice forms as a suspension of crystals in oceans, lakes, rivers, and sub-glacial ice streams from liquid water supercooled beneath its freezing temperature (Martin and Kauffman, 1981; Lawson et al., 1998). Supercooled water at the surface of the ocean occurs when it is cooled efficiently by the atmosphere. Such conditions can occur in gaps in the ice pack (called leads) and in extensive areas of open water (called polynyas), as observed by Skogseth et al. (2009). In some Antarctic regions, frazil-ice growth in supercooled water also contributes to the accretion of platelet ice on the underside of sea ice (e.g. Gough et al., 2012; Langhorne et al., 2015).

Frazil ice can also form underneath floating ice shelves at the margins of
the Antarctic continent. Plumes of relatively fresh, buoyant “ice shelf
water” (ISW) flow along the underside of the ice shelves. These rise over a
depth range of about a kilometre, a range associated with significant
variation of the pressure-dependent freezing temperature of seawater, which
varies by −0.76^{∘}C km^{−1} with depth (Millero and Leung, 1976).
Consequently, the temperature of a rising plume can fall beneath the
in situ freezing temperature (Lewis and Perkin, 1986), triggering the
formation of frazil ice. Some of the ice precipitates onto the base of the
ice shelf, where it forms so-called marine ice, which has a granular texture.
The presence of marine ice was inferred and subsequently observed by drilling
boreholes through the ice shelf (Engelhardt and Determann, 1987; Oerter et al., 1992). Frazil-ice
formation can affect the dynamics of these plumes by changing their buoyancy
directly (because ice is less dense than water) and by changing their
temperature and salinity.

It is just becoming possible to assess the role of frazil-ice formation on sea ice and ocean conditions through large-scale models (e.g. Galton-Fenzi et al., 2012; Wilchinsky et al., 2015; Smedsrud and Martin, 2015). Such models rely on previous theoretical work concerning frazil-ice dynamics, which was pioneered by Daly (1984). Models of frazil-ice dynamics have been applied to the study of frazil in the upper ocean (Svensson and Omstedt, 1994, 1998; Heorton et al., 2017) and also to the study of frazil ice beneath ice shelves (Jenkins and Bombosch, 1995; Khazendar and Jenkins, 2003; Smedsrud and Jenkins, 2004; Holland and Feltham, 2005; Jordan et al., 2014, 2015). The theory of frazil-ice dynamics involves parameterizations of a number of physical processes that affect the evolution of a population of ice crystals. In this paper, we revisit the theory of frazil-ice dynamics, taking into account new understanding of the microphysics of crystal growth (Rees Jones and Wells, 2015), before suggesting likely implications for these large-scale models.

In a recent paper, Rees Jones and Wells (2015) presented numerical evidence that the growth rate of ice crystals has been significantly underestimated in some previous studies as detailed below. In this section, we briefly review this finding and explain the underlying physical ideas.

Frazil ice is observed to consist of roughly disk-shaped crystals that
typically have a much greater radius *R* than thickness *H*
(McFarlane et al., 2014). Crystal growth is predominantly radial, with attachment
kinetics limiting growth in the basal plane and maintaining the disk-shaped
geometry for crystals of modest size (Fujioka and Sekerka, 1974). The radial growth
rate *G* of a frazil crystal depends on the rate at which the latent heat
released by crystal growth is transported away from the crystal. In general,
the radial growth rate can be written in the form

$$\begin{array}{}\text{(1)}& {\mathit{\rho}}_{\mathrm{i}}LG=(\mathit{Nu}\phantom{\rule{0.125em}{0ex}}{k}_{\mathrm{l}}\mathrm{\Delta}T/H)f,\end{array}$$

where *ρ*_{i} is the density of ice; *L* is the latent heat of
solidification; *Nu* is the crystal Nusselt number, which equals 1
for purely diffusive growth and can be enhanced by flow; *k*_{l} is the thermal
conductivity of the liquid phase; Δ*T* is the amount of supercooling
below the in situ freezing temperature; and *f* is a dimensionless geometric
factor. A helpful way to interpret Eq. (1) is to rearrange
it into an expression for the rate of crystal-mass growth, with the
ice-crystal mass *M*=*ρ*_{i}*π**R*^{2}*H*. We find

$$\begin{array}{}\text{(2)}& L{\displaystyle \frac{\mathrm{d}M}{\mathrm{d}t}}=\mathit{Nu}\phantom{\rule{0.125em}{0ex}}\phantom{\rule{0.125em}{0ex}}{k}_{\mathrm{l}}\mathrm{\Delta}T\mathrm{2}\mathit{\pi}Rf\propto \mathit{Nu}\phantom{\rule{0.125em}{0ex}}\phantom{\rule{0.125em}{0ex}}{k}_{\mathrm{l}}\mathrm{\Delta}T{\displaystyle \frac{A}{{\mathit{\delta}}_{\mathrm{T}}}}f.\end{array}$$

The right-hand side is the product of the area for heat transfer *A* and the
heat flux scale *k*_{l}Δ*T*∕*δ*_{T}, where *δ*_{T} is a
thermal boundary layer thickness. Numerical calculations of the temperature
distribution around an ice crystal (an example is shown in
Fig. 1) show that *δ*_{T}∝*H* near the crystal edges, which have an area *A*∝*R**H*. However,
*δ*_{T}∝*R* near the crystal faces, which have an area
*A*∝*R*^{2}. In either case, the ratio $A/{\mathit{\delta}}_{\mathrm{T}}\propto R$. Thus the scaling argument suggests *f*∝1 (cf.
Eq. 2). It is interesting to note that the mass growth
rate of spherical crystals is also proportional to crystal radius *R*, so the
rate of latent heat release seems to depend on crystal size *R* but not on
the details of the geometry.

We now consider three possible parameterizations of crystal growth, which we
denote ${f}_{\mathrm{1},\mathrm{2},\mathrm{3}}$. Numerical calculations of the heat transfer by diffusion
from a disk-shaped crystal (Rees Jones and Wells, 2015) show that the growth rate
depends logarithmically on aspect ratio ${f}_{\mathrm{1}}(h=H/\mathrm{2}R)=\mathrm{1}/[\mathrm{0.9008}-\mathrm{0.2634}\mathrm{log}(h\left)\right]$, which is similar to 1. Some previous studies are
also consistent with the scaling *f* ∼ 1. For example,
Svensson and Omstedt (1994) and Jenkins and Bombosch (1995) take *f*_{2}=1. By contrast, some
later studies are inconsistent with the scaling argument. For example,
Smedsrud and Jenkins (2004), Holland et al. (2007), and Galton-Fenzi et al. (2012) take
*A*∝*R**H* and *δ*_{T}∝*R*, which gives a growth
rate proportional to ${f}_{\mathrm{3}}\equiv H/R\ll \mathrm{1}$, i.e. a very much smaller growth
rate. A further complication arises in that it is sometimes additionally
assumed that the crystal aspect ratio $h=H/\mathrm{2}R$ is constant, rather than the
crystal thickness *H* being constant. In this case, *f*_{3}≡2*h*, which is
a constant, like *f*_{2}, but very much smaller (e.g. Smedsrud and Jenkins, 2004, take
*h*=0.02). These papers are illustrative of a wider range of studies
(e.g.
Svensson and Omstedt, 1998; Khazendar and Jenkins, 2003; Holland and Feltham, 2005; Jordan et al., 2014, 2015; Wilchinsky et al., 2015; Smedsrud and Martin, 2015);
recently it appears that growth law *f*_{3} has been used most commonly, if not
exclusively. In summary, numerical calculations show that there is only weak
dependence on aspect ratio: *f*_{1} is typically close to *f*_{2}; however, *f*_{1}
is some 10–100 times greater than *f*_{3}, as illustrated in Fig. 2.

The presence of salt in seawater reduces the crystal growth rate because the
supercooling is reduced and salt rejected by the growing crystal needs to
diffuse away. Numerical calculations performed to investigate these effects
(Rees Jones and Wells, 2015) support the scaling argument used to account for the
effect of salt by Galton-Fenzi et al. (2012), which in turn was based on
Holland and Jenkins (1999). For practical modelling purposes, the supercooling needs
to be adjusted for the salt content of seawater, and the Nusselt number
should be reduced to account for salt diffusion. A good approximation based
on our numerical calculations is $\mathit{Nu}={\left[\mathrm{1}+\mathrm{1.4}\phantom{\rule{0.125em}{0ex}}\times \phantom{\rule{0.125em}{0ex}}(-aS{k}_{\mathrm{l}})/\left({D}_{\mathrm{S}}{\mathit{\rho}}_{\mathrm{i}}L\right)\right]}^{-\mathrm{1}},$ where *a*<0 is the rate of change of
freezing temperature with salinity *S*, and *D*_{S} is the diffusivity of salt
in water.

2 Frazil-ice dynamics

Back to toptop
How does the growth rate of an ice crystal affect the overall ice production rate of a system? To address this question we need to investigate “frazil-ice dynamics”. We follow the comprehensive framework of the influential reviews of Daly (1984, 1994), which accounts for the evolution of a crystal size distribution in time, in space, and in crystal size space. The evolution occurs through crystal nucleation, growth, flocculation, break-up, and transport by fluid motion. There is a high degree of uncertainty in the rate of each of these processes, which in turn drives uncertainty in predictions of crucial, environmentally relevant quantities, such as the total ice production rate.

In this section we set out continuum equations that describe the evolution of
frazil ice in a general framework that can be applied to a wide range of
specific situations, before later focussing on examples of ice growth in a
mixed layer and under ice shelves. Suppose that the size of a crystal can be
characterized by a single length scale *R*, the radius of a disk-shaped
crystal. We introduce the crystal number density *n*, which is defined as the
number of crystals per unit volume of mixture per unit length in crystal size
space. Other quantities can be derived from *n*. For example, the crystal
concentration density *c*=*n**V*, where *V*=*π**R*^{2}*H* is the volume of a
disk-shaped crystal of thickness *H*, and the total crystal concentration
$C={\int}_{\mathrm{0}}^{\mathrm{\infty}}c\phantom{\rule{0.33em}{0ex}}\mathrm{d}R$. Note that *C* is the volume occupied by ice
crystals per unit volume of mixture. The total number density
$N={\int}_{\mathrm{0}}^{\mathrm{\infty}}n\phantom{\rule{0.33em}{0ex}}\mathrm{d}R$. The density *n* is a function of time *t*, position
** x**, and crystal size

$$\begin{array}{ll}\text{(3)}& {\displaystyle \frac{\partial n}{\partial t}}+& {\displaystyle}\mathrm{\nabla}\cdot \left(\mathit{u}n\right)-\mathrm{\nabla}\cdot \left({D}_{\mathrm{c}}\mathrm{\nabla}n\right)={\displaystyle}-& {\displaystyle \frac{\partial}{\partial R}}\left(Gn\right)-W{\displaystyle \frac{\partial n}{\partial z}}-{\displaystyle \frac{\mathrm{1}}{V}}{\displaystyle \frac{\partial}{\partial R}}\left(BVn\right)+\dot{N}\mathit{\delta}\left(R\right),\end{array}$$

where ** u** is the fluid velocity and

The first term represents crystal growth, where *G* is the radial crystal
growth rate discussed in Sect. 1.2. For compactness,
we rewrite Eq. (1) as

$$\begin{array}{}\text{(4)}& G={G}_{\mathrm{0}}f,\end{array}$$

where ${G}_{\mathrm{0}}=\mathit{Nu}\phantom{\rule{0.125em}{0ex}}{k}_{\mathrm{l}}\mathrm{\Delta}T/\left({\mathit{\rho}}_{\mathrm{i}}LH\right)$, and *f* is given by one
of the three growth laws. The effect of this term is to shift the crystal
size distribution to larger radii *R*, without increasing the total number of
crystals. Thus growth increases crystal concentration but not the number of
crystals.

The second term represents removal due to buoyant crystal rise, where *W* is
an effective crystal rise speed. It is well established that larger crystals
rise faster. Recent experimental observations and the parameterization of
crystal rise speed are discussed in McFarlane et al. (2014). For the simplest
treatment, we use a linear relationship,

$$\begin{array}{}\text{(5)}& W={W}_{\mathrm{0}}R,\end{array}$$

with *W*_{0}=16 s^{−1}, because more complicated
parameterizations do not fit the data much better than this simple fit.
Indeed, such a relationship is consistent with the crystal rise being a
Stokes settling velocity under the assumption that crystal thickness is
constant. The drag is proportional to *μ**W**R*, and the buoyancy is
proportional to Δ*ρ**g**R*^{2}*H*, where Δ*ρ* is the density
difference between ice and water. Thus balancing drag and buoyancy yields *W*∝*R*.

The third term represents the net effect of the processes of flocculation and
break-up, where *B* is the rate. Positive *B* corresponds to flocculation
greater than break-up. Note that this term is constructed to conserve crystal
volume, which is physically appropriate. To see this, multiply Eq. (3) by *V* and integrate from *R*=0 to *R*=∞. The total
volume of ice is unaffected by the flocculation term. To our knowledge, this
term has received relatively little attention within the frazil-ice
literature. One exception, Svensson and Omstedt (1994), includes it and takes

$$\begin{array}{}\text{(6)}& B={B}_{\mathrm{0}}{R}^{\mathrm{2}}.\end{array}$$

As a technical aside, we note that Svensson and Omstedt (1994) describe their
flocculation law as linear. However, this linearity applies only to the
particular discrete set of equations they present, which use logarithmically
spaced size classes. At the continuum level, the quadratic Eq. (6) applies. There is no direct evidence for the form of this
relationship, although Svensson and Omstedt (1994) found it helpful in fitting some
experimental data. Their choice matches the intuition that larger crystals
might flocculate more readily since they are more likely to come into near
contact with other crystals. However, it does not account for the fact that
flocculation should increase with frazil concentration. A fuller treatment
would take *B* as an integral of an interaction kernel *K* multiplied by
number density over crystal radius, $B={\int}_{\mathrm{0}}^{\mathrm{\infty}}Kn\phantom{\rule{0.125em}{0ex}}\mathrm{d}R$. This kind of
approach has proved fruitful in the theory of sea-ice thickness and floe-size
distributions (Thorndike, 2000; Godlovitch et al., 2011; Horvat and Tziperman, 2015; Toppaladoddi and Wettlaufer, 2015). In
view of the considerable uncertainties in parameterizing flocculation, we
neglect this process in all of our calculations (*B*=0). Indeed, even the
sign of *B* is uncertain, as it not clear whether flocculation or break-up
dominates (and the balance of these processes may well depend on the fluid
dynamical conditions). If break-up dominates (perhaps in more turbulent
environments), setting *B*=0 might overestimate the number of large crystals.
Conversely, if flocculation dominates, setting *B*=0 might underestimate the
number of large crystals.

The fourth term represents crystal nucleation, where $\dot{N}$ is nucleation
rate. We use the mathematical construct of a delta function *δ*(*R*) in
Eq. (3) to express the fact that nucleated crystals are
extremely small. By integrating Eq. (3) from
*R*=0 to $R=\stackrel{\mathrm{\u0303}}{\mathit{\u03f5}}>\mathrm{0}$, it can be shown that the
nucleation flux balances the growth of small crystals, and

$$\begin{array}{}\text{(7)}& \underset{\stackrel{\mathrm{\u0303}}{\mathit{\u03f5}}\to {\mathrm{0}}^{+}}{lim}{\left.Gn\right|}_{R=\stackrel{\mathrm{\u0303}}{\mathit{\u03f5}}}=\dot{N},\end{array}$$

since the other terms give rise to contributions that are proportional to
$\stackrel{\mathrm{\u0303}}{\mathit{\u03f5}}$ and vanish in the limit $\stackrel{\mathrm{\u0303}}{\mathit{\u03f5}}\to {\mathrm{0}}^{+}$.
After some primary nucleation event, nucleation is assumed to be dominated by
secondary nucleation, sometimes called collisional breeding. Indeed,
Daly (1984) argues that homogenous and heterogenous nucleation are
extremely unlikely to occur in natural systems because the levels of
supercooling achieved are less than 1 ^{∘}C. We follow, for
instance, Svensson and Omstedt (1994) and suppose that collisions between crystals cause
microscopic pieces of ice to break off, which in turn become new crystals with
very small radius. The total nucleation rate depends on the the rate at which
a volume is swept out by a crystal and the crystal number density. We write

$$\begin{array}{}\text{(8)}& \dot{N}=\stackrel{\mathrm{\u0303}}{n}\underset{\mathrm{0}}{\overset{\mathrm{\infty}}{\int}}\mathit{\pi}{R}^{\mathrm{2}}{U}_{r}n\left(R\right)\phantom{\rule{0.125em}{0ex}}\mathrm{d}R,\end{array}$$

where

$$\begin{array}{}\text{(9)}& {U}_{r}=\sqrt{\mathrm{4}\mathit{\u03f5}{R}^{\mathrm{2}}/\mathrm{15}\mathit{\nu}+({W}_{\mathrm{0}}R{)}^{\mathrm{2}}}\equiv {U}_{\mathrm{0}}R\end{array}$$

is an effective collisional velocity scale taken to be the geometric mean of
a velocity scale based on turbulent motions (*ϵ* is the turbulent
dissipation rate, *ν* is the kinematic viscosity) and one based on buoyant
crystal rise. We define ${U}_{\mathrm{0}}=\sqrt{\mathrm{4}\mathit{\u03f5}/\mathrm{15}\mathit{\nu}+{W}_{\mathrm{0}}^{\mathrm{2}}}$ and use a
value $\mathit{\nu}=\mathrm{2}\times {\mathrm{10}}^{-\mathrm{6}}$ m^{2} s^{−1}. The nucleation efficiency
scale $\stackrel{\mathrm{\u0303}}{n}=min(N,{\stackrel{\mathrm{\u0303}}{n}}_{\mathrm{max}})$, where
${\stackrel{\mathrm{\u0303}}{n}}_{\mathrm{max}}$ is a calibration parameter that limits the
efficiency of secondary nucleation. Smedsrud (2002) points out that some
of the nucleated crystals will be below the so-called “critical size” for
crystals to grow, so it is plausible that $\stackrel{\mathrm{\u0303}}{n}<N$, but it must be
conceded that this parameterization is rather ad hoc. We use this
formulation primarily for consistency with previous studies, to allow us to
isolate the effect of crystal growth rate. It is simply a continuous version
of that used by, for example, Svensson and Omstedt (1994), Smedsrud (2002), Smedsrud and Jenkins (2004), and Holland and Feltham (2005).
Before the efficiency cap is reached, secondary nucleation is a quadratic in
the number of crystals, leading to very rapid growth in crystal number.

Equation (3) can be discretized in radial space to
facilitate numerical solution, following, for instance, Svensson and Omstedt (1994). The spatial
problem is a standard advection–diffusion problem, so we do not discuss here
how to discretize the left-hand side of Eq. (3) and
focus on the crystal interaction terms on the right-hand side. Let *R*_{i} be a
discrete set of points in radial space, where $\mathrm{1}\phantom{\rule{0.125em}{0ex}}\le \phantom{\rule{0.125em}{0ex}}i\phantom{\rule{0.125em}{0ex}}\le \phantom{\rule{0.125em}{0ex}}M$. We introduce
the notation *W*_{i}=*W*(*R*_{i}), *G*_{i}=*G*(*R*_{i}), and ${V}_{i}=V\left({R}_{i}\right)=\mathit{\pi}{R}_{i}^{\mathrm{2}}H$. We work
in terms of the total number of particles in size class *i*, denoted *m*_{i},
which evolves according to

$$\begin{array}{ll}\text{(10)}& {\displaystyle \frac{\partial {m}_{i}}{\partial t}}& {\displaystyle}=-{\mathrm{\Gamma}}_{i}{m}_{i}+{\mathrm{\Gamma}}_{i-\mathrm{1}}{m}_{i-\mathrm{1}}-{W}_{i}{\displaystyle \frac{\partial {m}_{i}}{\partial z}}-{\mathit{\alpha}}_{i}{m}_{i}\phantom{\rule{0.125em}{0ex}}\phantom{\rule{0.125em}{0ex}}\phantom{\rule{0.125em}{0ex}}(i\ge \mathrm{2}),{\displaystyle \frac{\partial {m}_{i}}{\partial t}}& {\displaystyle}=-{\mathrm{\Gamma}}_{i}{m}_{i}-{W}_{i}{\displaystyle \frac{\partial {m}_{i}}{\partial z}}+\sum _{j=\mathrm{2}}^{j=M}\stackrel{\mathrm{\u0303}}{n}\mathit{\pi}{R}_{j}^{\mathrm{2}}{U}_{r}\left({R}_{j}\right)\phantom{\rule{0.125em}{0ex}}\phantom{\rule{0.125em}{0ex}}\phantom{\rule{0.125em}{0ex}}(i=\mathrm{1}),\end{array}$$

where

$$\begin{array}{}\text{(11)}& {\displaystyle}{\mathrm{\Gamma}}_{i}& {\displaystyle}={\displaystyle \frac{{G}_{i}\mathrm{2}\mathit{\pi}{R}_{i}H}{{V}_{i+\mathrm{1}}-{V}_{i}}},\text{(12)}& {\displaystyle}{\mathit{\alpha}}_{i}& {\displaystyle}=\stackrel{\mathrm{\u0303}}{n}\mathit{\pi}{R}_{i}^{\mathrm{2}}{U}_{r}\left({R}_{i}\right){\displaystyle \frac{{V}_{\mathrm{1}}}{{V}_{i}}},\phantom{\rule{2em}{0ex}}(i\ge \mathrm{2}).\end{array}$$

The discrete distribution *n*(*R*_{i}) can be recovered: ${n}_{i}={m}_{i}/\mathrm{\Delta}{R}_{i}$,
where $\mathrm{\Delta}{R}_{i}={R}_{i+\mathrm{1}}-{R}_{i}$. We note that Eq. (10) is only a first-order discretization in radial
space, so an alternative approach could be to use a second-order
discretization. This numerical representation is conservative, and we use a
formulation of secondary nucleation (in terms of *α*_{i}) that conserves
crystal volume even when *V*_{1} is non-zero. Note that in the limit *R*_{1}→0 and Δ*R*_{i}→0 we recover the
continuum equations discussed above.

Equation (10) is equivalent to Eq. (1) in Svensson and Omstedt (1994). They demonstrate that this model is capable of reproducing the main features of the laboratory experiments of Michel (1963) and Carstens (1966), so we do not include any experimental comparison here. However, we discuss (Sect. 3.3) how such consistency is insufficient to fully validate the model. For practical purposes, we find it advantageous to use a logarithmically spaced set of crystal sizes and test the accuracy of our discretization by increasing the number of size classes to 1024. We find that good accuracy can be achieved with 128 classes, but accuracy noticeably degrades beneath this (cf. Holland and Feltham, 2005). Software code to reproduce the calculations in the paper is available (Rees Jones, 2017).

3 Frazil ice in a mixed layer

Back to toptop
The upper layer of a lake or ocean can sometimes be approximated as a
well-mixed layer, an approximation that can also be applied to the laboratory
experiments of Michel (1963) and Carstens (1966). We assume that
background turbulent stirring is sufficient to keep the layer well mixed such
that all physical quantities (temperature and crystal size distribution) are
uniform over the layer. Such turbulence might be driven, for example in the
oceans, by wind, waves, and buoyancy-driven convection. A turbulent flow is
mechanically driven in laboratory experiments. Thus we only need to solve
evolution equations for average physical quantities across the layer. This
approximation also significantly simplifies Eq. (3)
while still retaining the key frazil-ice dynamics. Averaging Eq. (3) over the mixed layer of depth *D* yields

$$\begin{array}{}\text{(13)}& {\displaystyle \frac{\partial n}{\partial t}}=-{\displaystyle \frac{\partial}{\partial R}}\left(Gn\right)-\mathit{\gamma}n+\dot{N}\mathit{\delta}\left(R\right),\end{array}$$

where $\mathit{\gamma}=W/D$ is an effective gravitational removal term. In reality,
crystal concentration would tend to decrease with depth (Svensson and Omstedt, 1998)
because of crystal buoyancy. Nevertheless, $\mathit{\gamma}=W/D$ is an appropriate
scaling relationship because removal increases with crystal rise velocity *W*
and decreases with mixed-layer depth *D* because turbulent eddies act to mix
crystals down to that depth range. This type of depth-integrated
representation of the process of gravitational removal has been used
successfully in previous studies of turbulent, particle-laden gravity
currents (Bonnecaze et al., 1993).

The temperature of the mixed layer or tank evolves due to heat extraction to
the atmosphere per unit volume *Q* and release of the latent heat of
solidification:

$$\begin{array}{}\text{(14)}& {\mathit{\rho}}_{\mathrm{l}}{c}_{\mathrm{l}}{\displaystyle \frac{\mathrm{d}T}{\mathrm{d}t}}=-Q+\mathrm{2}\mathit{\pi}\mathit{Nu}{k}_{\mathrm{l}}\mathrm{\Delta}T\underset{\mathrm{0}}{\overset{\mathrm{\infty}}{\int}}fnR\phantom{\rule{0.33em}{0ex}}\mathrm{d}R.\end{array}$$

*ρ*_{l} is the density of the liquid phase, and *c*_{l} is
the specific heat capacity of the liquid phase. There is an implicit assumption that the ice removed through gravitational
settling does not inhibit heat loss to the atmosphere (by ice accumulation at
the surface); otherwise *Q* would decrease over time. Note that, in this
section, we neglect the depth dependence of the freezing temperature, which
affects the supercooling Δ*T*. This is a good approximation provided
the mixed layer is relatively shallow, but it would not be appropriate for mixed
layers deeper than *O*(100 m).

In a typical experiment, a relatively small number of crystals are seeded into supercooled water, for example by running a saw blade over a block of ice. Over time, the number of crystals undergoes a period of rapid growth, producing an optically dense suspension (Hanley and Tsang, 1984). Svensson and Omstedt (1994) include a figure from S. F. Daly (1992, personal communication) showing a period of rapid growth in the number of crystals: the total number of crystals increased by 4 orders of magnitude over around 250 s. The frazil-ice explosion was observed to reduce the supercooling in the mixed layer to a small residual amount.

Our goal in this section is to ascertain the conditions under which such a
frazil-ice explosion can occur and hence to determine conditions for their
occurrence in geophysical settings as well as laboratory experiments. To
motivate our approach, we consider the time evolution of a mixed layer seeded
with some small initial concentration and cooled beneath freezing by a
constant flux *Q*. The initial size distribution of crystals is taken to be
uniform on [0,2*R*_{0}], and we vary the total number of crystals to vary the
initial concentration. Throughout this section, we fix the crystal growth law
*f*_{2}=1. We present an example of such a calculation in Fig. 3. In one calculation, with slightly less ice initially
present (blue curve in Fig. 3), all of the ice is
removed (by gravitational rise), and supercooling continues to build.
Eventually we would expect heterogenous and later homogenous nucleation to
occur (Daly, 1984), but we do not model these processes. In the other
calculation, with slightly more ice initially present (red curve in Fig. 3), the ice concentration increases rapidly before
attaining a steady state in which supercooling is almost exhausted (see
Sect. 3.4). We consider this an example of a “frazil-ice
explosion” of the kind observed in experiments. A greater initial seeding
concentration of ice always makes an explosion more likely, so we investigate
the minimum initial concentration (or equivalently number of crystals, if the
initial size is fixed) required to trigger an explosion as a function of the
other parameters of the system.

We summarize the results of our investigations in Fig. 4.
Increasing the turbulent intensity *ϵ* (Fig. 4a)
increases the rate of secondary nucleation, since crystals are more likely to
collide, which promotes frazil explosions. Increasing the mixed-layer depth
*D* (Fig. 4b) reduces the rate at which crystals are removed
gravitationally, which again promotes frazil explosions. A slightly weaker
effect (note the different scale on the axis) is that increasing the cooling
rate *Q* (Fig. 4c) promotes frazil explosions. The direct
mechanism is that higher cooling promotes ice growth, increasing the frazil
concentration. However, there is also an important indirect mechanism: ice
growth shifts the crystal size distribution to larger crystal sizes, which
are more likely to collide, leading to greater secondary nucleation. This
effect is somewhat offset by the fact that larger crystals are also more
effectively removed by gravitational rise.

These mechanisms can be understood more quantitatively by scaling analysis. First, we integrate Eq. (13) across crystal sizes to obtain an evolution equation for the total number density of crystals (recalling the growth shifts the size distribution but does not change the total number of crystals):

$$\begin{array}{}\text{(15)}& {\displaystyle \frac{\mathrm{d}N}{\mathrm{d}t}}=\dot{N}-\underset{\mathrm{0}}{\overset{\mathrm{\infty}}{\int}}n\mathit{\gamma}\phantom{\rule{0.125em}{0ex}}\mathrm{d}R.\end{array}$$

If gravitation removal were to act alone, we would find that

$$\begin{array}{}\text{(16)}& {\displaystyle \frac{\mathrm{d}N}{\mathrm{d}t}}=-{\displaystyle \frac{{W}_{\mathrm{0}}}{D}}\stackrel{\mathrm{\u203e}}{R}N,\end{array}$$

where $\stackrel{\mathrm{\u203e}}{R}$ is the mean crystal size. Thus crystals are removed
exponentially on a settling timescale
$\mathit{\tau}=D/{W}_{\mathrm{0}}\stackrel{\mathrm{\u203e}}{R}\phantom{\rule{0.125em}{0ex}}\approx \phantom{\rule{0.125em}{0ex}}\mathrm{300}$ s (based on *D*=1 m and $\stackrel{\mathrm{\u203e}}{R}=\mathrm{0.2}$ mm, the
initial average crystal radius), which is commensurate with the evolution
timescale observed in Fig. 3.

Second, we consider a balance between secondary nucleation and gravitational removal. We expect a frazil explosion when the secondary nucleation (Eq. 8) is much greater than gravitational removal:

$$\begin{array}{}\text{(17)}& {\displaystyle}{N}^{\mathrm{2}}{U}_{\mathrm{0}}\mathit{\pi}{\stackrel{\mathrm{\u203e}}{R}}^{\mathrm{3}}& {\displaystyle}\gg {\displaystyle \frac{N{W}_{\mathrm{0}}\stackrel{\mathrm{\u203e}}{R}}{D}},\text{(18)}& {\displaystyle}\Rightarrow N& {\displaystyle}\gg {N}_{\mathrm{crit}.}\sim {\displaystyle \frac{{W}_{\mathrm{0}}}{{U}_{\mathrm{0}}{\stackrel{\mathrm{\u203e}}{R}}^{\mathrm{2}}D}}.\end{array}$$

If $\stackrel{\mathrm{\u203e}}{R}$ were given by the initial average crystal radius, then in
terms of the external parameters of the system shown in Fig. 4 we would naively expect *N*_{crit.} to decrease
with turbulent intensity (inversely proportional to *U*_{0}) and mixed-layer depth
(inversely proportional to *D*), and be independent of *Q*. The first
prediction (Fig. 4a) is supported by the numerical results.
However, the second prediction (Fig. 4b) and third
prediction (Fig. 4c) are not (the dashed curves do not agree
with the numerical results). The resolution of these discrepancies lies in
recognizing that the average crystal size $\stackrel{\mathrm{\u203e}}{R}$ is not a constant
external parameter (i.e. set by the initial condition as a consequence of the
seeding strategy) but rather depends on crystal growth.

We now suppose that the average crystal size is determined by the amount a
crystal can grow over a crystal removal timescale *τ*, i.e.

$$\begin{array}{}\text{(19)}& \stackrel{\mathrm{\u203e}}{R}\sim G\mathit{\tau}.\end{array}$$

This is a good approximation provided *G**τ* is much larger than the
initial crystal size. The growth rate *G* is proportional to the
supercooling; in particular $G=\mathit{Nu}f\phantom{\rule{0.125em}{0ex}}{k}_{\mathrm{l}}\mathrm{\Delta}T/\left({\mathit{\rho}}_{\mathrm{i}}LH\right)$.
We can estimate the supercooling from the heat balance Eq. (14), in which the crystal growth term is negligible until
the frazil explosion occurs. We find

$$\begin{array}{ll}{\displaystyle}{\mathit{\rho}}_{\mathrm{l}}{c}_{\mathrm{l}}\mathrm{\Delta}T& {\displaystyle}\sim Q\mathit{\tau},\\ \text{(20)}& {\displaystyle}\Rightarrow G& {\displaystyle}\sim {\displaystyle \frac{\mathit{Nu}f\phantom{\rule{0.125em}{0ex}}{k}_{\mathrm{l}}Q\mathit{\tau}}{{\mathit{\rho}}_{\mathrm{l}}{c}_{\mathrm{l}}{\mathit{\rho}}_{\mathrm{i}}LH}}.\end{array}$$

We substitute Eq. (20) into Eq. (19) and recall that $\mathit{\tau}=D/{W}_{\mathrm{0}}\stackrel{\mathrm{\u203e}}{R}$. Rearranging for $\stackrel{\mathrm{\u203e}}{R}$, we find

$$\begin{array}{}\text{(21)}& {\stackrel{\mathrm{\u203e}}{R}}^{\mathrm{3}}\sim {\displaystyle \frac{\mathit{Nu}f\phantom{\rule{0.125em}{0ex}}{k}_{\mathrm{l}}Q}{{\mathit{\rho}}_{\mathrm{l}}{c}_{\mathrm{l}}{\mathit{\rho}}_{\mathrm{i}}LH}}{\left({\displaystyle \frac{D}{{W}_{\mathrm{0}}}}\right)}^{\mathrm{2}}.\end{array}$$

We then substitute this estimate for $\stackrel{\mathrm{\u203e}}{R}$ into Eq. (18) and obtain

$$\begin{array}{}\text{(22)}& {N}_{\mathrm{crit}.}\sim {\displaystyle \frac{\mathrm{1}}{{U}_{\mathrm{0}}}}{\left({\displaystyle \frac{{W}_{\mathrm{0}}}{D}}\right)}^{\mathrm{7}/\mathrm{3}}{\left({\displaystyle \frac{\mathit{Nu}f\phantom{\rule{0.125em}{0ex}}{k}_{\mathrm{l}}Q}{{\mathit{\rho}}_{\mathrm{l}}{c}_{\mathrm{l}}{\mathit{\rho}}_{\mathrm{i}}LH}}\right)}^{-\mathrm{2}/\mathrm{3}}.\end{array}$$

Equation (22) is very appealing because it can explain nearly
all the results presented in Fig. 4a, b, and c (the solid curves agree with the numerical results much
better than the dashed curves). The heat flux result is perhaps slightly
affected by the initial crystal size distribution at small *Q*, but overall
the agreement is very good.

In terms of our crystal growth rate, our scaling argument in Eq. (22) suggests that the faster growth laws would necessitate a smaller initial concentration of ice to trigger a frazil explosion, something that we observe in numerical experiments (cf. Fig. 6).

In conclusion, we find that the explosive growth of frazil ice requires a sufficiently large number of seed crystals. Seed crystals might be supplied from the atmosphere as sea spray freezes, from broken-off pieces of an ice shelf above a plume, or (perhaps unlikely) by sediment acting as nuclei for crystal growth. Such growth is promoted by high turbulent intensity, a deeper mixed layer, and strong cooling rate (or larger seed crystals).

Figure 5 shows an example of how the crystal size distribution (CSD) evolves when a frazil explosion occurs. Initially, the larger seed crystals are removed gravitationally, while crystals are nucleated at the smallest size due to collisional breeding. These crystals grow. Note the “travelling-wave” type solutions evident at 100 and 200 s with the radius of crystals increasing over time. Indeed, there are travelling-wave solutions to Eq. (13) if crystal growth is the only process that affects the CSD evolution. Finally, a steady-state distribution is achieved, which we discuss in more detail in Sect. 3.4.

We next consider the impact of different parameterizations of crystal growth
*f*_{1−3}. One main experimental measurement is mixed-layer temperature as a
function of time. We find that this observable is sensitive to the crystal
growth rate, as shown in Fig. 6. Faster crystal growth
means a faster increase in crystal concentration, with the peak growth rate
occurring several hundred seconds earlier. This in turn means that the peak
supercooling is lower, because of the latent heat liberated by crystal
growth. These differences are experimentally detectable.

Our new parameterization produces broadly similar transient evolution curves to the older model of Svensson and Omstedt (1994). It is therefore encouraging to note that Svensson and Omstedt (1994) were able to use their model to explain the experimental observations of degree of supercooling. However, demonstration of consistency with experiments does not conclusively show that a parameterization of crystal growth is correct, because other factors also affect the predicted supercooling, such as the size distribution of the initial seed crystals (which was not controlled in the experiments of Michel (1963) and Carstens (1966) that Svensson and Omstedt (1994) used to test their model) as shown in Fig. 7. Larger seed crystals grow more slowly and achieve greater maximum supercooling, which produces similar predictions to using a slower growth-rate law. This suggests that it is worthwhile for experimentalists to try to measure crystal sizes, as well as supercooling, in order to discriminate between models.

In conclusion, we have shown that crystal growth rate significantly affects
the transient evolution of the crystal size distribution. Further experimental
observations are needed to discriminate between models. Geophysically, we
note that the differences between models occur on timescales of a few hundred
seconds. This timescale is proportional to mixed-layer depth, so a deeper
mixed layer would be associated with even longer transient frazil-ice
dynamics. The transient differences are therefore likely to be most
significant to systems where the frazil ice is subject to processes that act
on similar or shorter timescales to the transient relief of supercooling.
(For processes that act on longer timescales, the frazil-ice dynamics would
have equilibrated to the steady states discussed in the next section.) For
example, a lateral current of 0.1 m s^{−1} would take 100 s to move
material across a lead that is 10 m wide. These numbers offer some
indication that these transient model differences may well be geophysically
significant. Indeed, we show an example in the context of Ice Shelf Water
plumes in Sect. 4.

We observed that the crystal size distribution evolves to a steady state. In this section we study these steady states by numerically integrating our transient model to reach a steady state for each of the three growth laws, and by finding analytical steady-state solutions of the governing equations for two of the growth laws. We present an example of numerically obtained steady states in Fig. 8. Changing the growth law subtly shifts the crystal size distribution.

In order to understand better the physical processes involved in maintaining this steady state, we analyse the steady-state solutions of Eq. (13), namely

$$\begin{array}{}\text{(23)}& {\displaystyle \frac{\partial}{\partial R}}\left({G}_{\mathrm{0}}fn\right)+{\mathit{\gamma}}_{\mathrm{0}}Rn=\dot{N}\mathit{\delta}\left(R\right),\end{array}$$

where ${\mathit{\gamma}}_{\mathrm{0}}={W}_{\mathrm{0}}/D$. We start with the growth law *f*=*f*_{2} (a constant).
First, we integrate Eq. (23) when *R*>0 to
obtain

$$\begin{array}{}\text{(24)}& n={n}_{\mathrm{0}}\mathrm{exp}\left(-{\displaystyle \frac{{\mathit{\gamma}}_{\mathrm{0}}}{\mathrm{2}{G}_{\mathrm{0}}{f}_{\mathrm{2}}}}{R}^{\mathrm{2}}\right).\end{array}$$

Second, at $R={\mathrm{0}}^{+}$, Eq. (7) implies that

$$\begin{array}{}\text{(25)}& {G}_{\mathrm{0}}{f}_{\mathrm{2}}n(R={\mathrm{0}}^{+})=\mathit{\pi}{U}_{\mathrm{0}}{\stackrel{\mathrm{\u0303}}{n}}_{\mathrm{max}}\underset{\mathrm{0}}{\overset{\mathrm{\infty}}{\int}}n\left(R\right){R}^{\mathrm{3}}\phantom{\rule{0.125em}{0ex}}\mathrm{d}R,\end{array}$$

where we assume that the total number of crystals exceeds ${\stackrel{\mathrm{\u0303}}{n}}_{\mathrm{max}}$ so $\stackrel{\mathrm{\u0303}}{n}={\stackrel{\mathrm{\u0303}}{n}}_{\mathrm{max}}$. This is reasonable because there is a very large number of crystals after a frazil-ice explosion has occurred. Equation (25) can be manipulated by substituting in Eq. (24) and integrating to show that ${G}_{\mathrm{0}}{f}_{\mathrm{2}}={\mathit{\gamma}}_{\mathrm{0}}^{\mathrm{2}}/\mathrm{2}\mathit{\pi}{U}_{\mathrm{0}}{\stackrel{\mathrm{\u0303}}{n}}_{\mathrm{max}}.$ This expression allows the steady-state supercooling to be calculated since ${G}_{\mathrm{0}}=\mathit{Nu}{k}_{\mathrm{l}}\mathrm{\Delta}T/{\mathit{\rho}}_{\mathrm{i}}LH$. Third, we use the overall heat balance from Eq. (14) in steady state,

$$\begin{array}{}\text{(26)}& Q=\mathrm{2}\mathit{\pi}\mathit{Nu}{k}_{\mathrm{l}}\mathrm{\Delta}T\underset{\mathrm{0}}{\overset{\mathrm{\infty}}{\int}}fnR\phantom{\rule{0.33em}{0ex}}\mathrm{d}R,\end{array}$$

to determine the unknown prefactor *n*_{0}. Finally, we calculate the average
crystal size (mean) $\stackrel{\mathrm{\u203e}}{R}$, the total number of crystals *N*, and the
total crystal concentration *C*.

We then repeat the analysis for the growth law $f={f}_{\mathrm{3}}\equiv H/R$. We report the results in Table 1.

We conclude from this analysis that the average crystal radius is insensitive
to the crystal growth rate. This initially surprising result can be
understood by considering that the balance between growth and precipitation
at large crystal sizes gives ${G}_{\mathrm{0}}f\sim {\mathit{\gamma}}_{\mathrm{0}}{\stackrel{\mathrm{\u203e}}{R}}^{\mathrm{2}}$, while the
balance between growth and nucleation of the smallest crystals gives ${G}_{\mathrm{0}}f\sim {U}_{\mathrm{0}}{\stackrel{\mathrm{\u0303}}{n}}_{\mathrm{max}}{\stackrel{\mathrm{\u203e}}{R}}^{\mathrm{4}}$. The growth-rate-dependent
term *G*_{0}*f* can be eliminated between these equations, and

$$\begin{array}{}\text{(27)}& \stackrel{\mathrm{\u203e}}{R}\sim {\left({\displaystyle \frac{{\mathit{\gamma}}_{\mathrm{0}}}{{U}_{\mathrm{0}}{\stackrel{\mathrm{\u0303}}{n}}_{\mathrm{max}}}}\right)}^{\mathrm{1}/\mathrm{2}},\end{array}$$

in agreement with the expressions in Table 1. Therefore
average crystal size depends on (1) secondary nucleation (affected by
turbulent intensity through *U*_{0} and efficiency of secondary nucleation
through ${\stackrel{\mathrm{\u0303}}{n}}_{\mathrm{max}}$, where more secondary nucleation means
smaller crystals) and on (2) gravitational removal (a larger gravitational
removal rate prefactor *γ*_{0} means larger crystals). The first effect is
readily understood: secondary nucleation creates tiny crystals. The second is
more subtle because gravitational removal tends to remove larger crystals.
However, secondary nucleation increases more rapidly as a function of crystal
radius than gravitational removal. Thus enhanced gravitational settling
enhances the removal of large crystals and mutes their efficiency in driving
secondary nucleation, leading to the scaling given in Eq. (27). In geophysical settings and laboratory experiments,
the crystal rise velocity, mixed-layer depth, and turbulent intensity can be
measured much more easily than the efficiency of secondary nucleation. We
therefore suggest choosing this parameter to match with observations of
average crystal size. For example, choosing the reduced value
${\stackrel{\mathrm{\u0303}}{n}}_{\mathrm{max}}=\mathrm{4}\phantom{\rule{0.125em}{0ex}}\times \phantom{\rule{0.125em}{0ex}}{\mathrm{10}}^{\mathrm{5}}$ m^{−3} would give an
average crystal size of about 0.5 mm.

The total crystal concentration *C* is also insensitive to the crystal growth
rate. We can show this by continuing our scaling analysis as follows. From
Eq. (26), we estimate

$$\begin{array}{ll}{\displaystyle}Q& {\displaystyle}\sim \mathit{Nu}{k}_{\mathrm{l}}\mathrm{\Delta}Tf\stackrel{\mathrm{\u203e}}{R}N,\\ {\displaystyle}& {\displaystyle}\sim {\mathit{\rho}}_{\mathrm{i}}LH{G}_{\mathrm{0}}f\stackrel{\mathrm{\u203e}}{R}N,\\ {\displaystyle}& {\displaystyle}\sim {\mathit{\rho}}_{\mathrm{i}}LH{\mathit{\gamma}}_{\mathrm{0}}{\stackrel{\mathrm{\u203e}}{R}}^{\mathrm{3}}N,\\ \text{(28)}& {\displaystyle}& {\displaystyle}\sim {\mathit{\rho}}_{\mathrm{i}}L{\mathit{\gamma}}_{\mathrm{0}}\stackrel{\mathrm{\u203e}}{R}C,\end{array}$$

where we have used ${G}_{\mathrm{0}}f\sim {\mathit{\gamma}}_{\mathrm{0}}{\stackrel{\mathrm{\u203e}}{R}}^{\mathrm{2}}$ from the growth-versus-settling balance. If we define a surface heat flux scale ${Q}_{\mathrm{surf}.}=QD$ and recall ${\mathit{\gamma}}_{\mathrm{0}}\stackrel{\mathrm{\u203e}}{R}={W}_{\mathrm{0}}\stackrel{\mathrm{\u203e}}{R}/D$, we find

$$\begin{array}{}\text{(29)}& C\sim {\displaystyle \frac{{Q}_{\mathrm{surf}.}}{{\mathit{\rho}}_{\mathrm{i}}L{W}_{\mathrm{0}}\stackrel{\mathrm{\u203e}}{R}}}.\end{array}$$

Thus at steady state, the total amount of frazil ice is determined by a balance between the surface heat flux and the rate of export of latent heat in the form of frazil ice that is removed gravitationally. This steady-state balance is unaffected by crystal growth rate (at least in the absence of advective processes).

4 Frazil-laden plume underneath an ice shelf

Back to toptop
Frazil ice also forms in plumes of ISW beneath floating ice shelves. A plume is fed by the discharge of subglacial meltwater at the start of the shelf and by melting from the shelf itself. These meltwaters are relatively fresh, so the plume rises buoyantly. The plume entrains ocean waters, resulting in an intermediate temperature and salinity called ISW. A full examination of the dynamics of these plumes is beyond the scope of this paper, and we refer the reader to previous studies by Jenkins (1991), Jenkins and Bombosch (1995), and Smedsrud and Jenkins (2004). Instead we focus more narrowly by considering a simple case study that illustrates the possible impact of different treatments of frazil-ice processes on the dynamics of an ISW plume. A linear ice shelf rises from a depth of 1400 m below sea level to a depth of 285 m below sea level over a horizontal distance of 600 km. The ambient seawater is treated as an approximation to High Salinity Shelf Water (HSSW) with a linear stratification. Jenkins and Bombosch (1995) conceived this setting as a simple configuration that is representative of a large Antarctic ice shelf. The plume becomes supercooled as it ascends the ice shelf base because of the fall in pressure and consequent change in the freezing temperature. This supercooling leads to a combination of frazil-ice formation and direct basal freezing. Frazil ice increases the plume buoyancy and so accelerates the plume. Thus we might naively expect that faster crystal growth would lead to higher frazil concentrations and faster-flowing plumes. In this section, we show that this expectation is confounded by complex feedbacks between plume dynamics and frazil-ice processes.

The plume model accounts for the evolution of plume depth *D*, the
depth-averaged plume velocity *U*, temperature *T*, and salinity *S* as a
function of distance *s* along the ice shelf. Note that we also average the
freezing temperature over the depth of the plume. The frazil-ice dynamics
part of the model is essentially the same as that described in
Eq. (3) but integrated over the depth of the plume. The
depth-averaged frazil crystal size distribution evolves according to

$$\begin{array}{}\text{(30)}& {\displaystyle \frac{\partial \left(DUn\right)}{\partial s}}=-D{\displaystyle \frac{\partial}{\partial R}}\left(Gn\right)-p\left(R\right)n+D\dot{N}\mathit{\delta}\left(R\right),\end{array}$$

where *p*(*R*) is the rate at which frazil precipitates onto the base of the
ice shelf.

We retain the approach of Smedsrud and Jenkins (2004) as far as possible. The full set
of governing equations is described in that paper. Software code to reproduce
the calculations in the paper is available (Rees Jones, 2017). Note that
our thermal calculation includes an estimate of the conductive heat flux into
the ice shelf (based on the thermal boundary layer parameterization of
Holland and Jenkins (1999), using a core ice-shelf temperature of −15^{∘}C; A. Jenkins, personal communication, 2014).
We use a large number of crystal size
classes in the discrete calculation (1000), to ensure the crystal size
distribution is well resolved. By contrast Smedsrud and Jenkins (2004) use only 10
classes. This affects the quantitative results but not the qualitative
behaviour of the system. One important difference compared to the mixed-layer
models (Sect. 3) is that the precipitation rate depends
both on crystal rise velocity and on the plume velocity, because
precipitation from a turbulent plume occurs when crystal buoyancy exceeds the
turbulent shear stress acting to keep it in suspension. Thus precipitation
occurs when the plume velocity *U* is less than some critical velocity *U*_{c}
that can be expressed in terms of a critical Shields number.

The dynamics of an ISW plume can be very sensitive to frazil-ice processes.
Our numerical investigation found two basic types of behaviour. Sometimes
frazil ice precipitates out over a relatively short distance of *O*(10 km) and the plume itself is barely affected by the frazil. At other
times the frazil ice is sustained over *O*(100 km) and the plume is
rendered more buoyant. We illustrate this range of behaviour and explain the
underlying physical mechanisms by varying the rates of secondary nucleation
and crystal growth (Fig. 9). In terms of secondary
nucleation, we consider no nucleation ${\stackrel{\mathrm{\u0303}}{n}}_{\mathrm{max}}=\mathrm{0}$ m^{−3}, intermediate nucleation
${\stackrel{\mathrm{\u0303}}{n}}_{\mathrm{max}}=\mathrm{500}$ m^{−3} comparable to Smedsrud and Jenkins (2004), and high nucleation
${\stackrel{\mathrm{\u0303}}{n}}_{\mathrm{max}}=\mathrm{4}\phantom{\rule{0.125em}{0ex}}\times \phantom{\rule{0.125em}{0ex}}{\mathrm{10}}^{\mathrm{6}}$ m^{−3} comparable to
Svensson and Omstedt (1994). In terms of crystal growth, we contrast a slow growth law
and a fast growth law. For a slow growth law, we use Smedsrud and Jenkins (2004), one
of the class of growth laws we labelled *f*_{3} previously. For a fast growth
law, we use Rees Jones and Wells (2015), labelled *f*_{1} previously. Calculations
with the growth law *f*_{2} introduced in Sect. 1.2 are
very similar to the results with *f*_{1}.

Our sensitivity experiments (Fig. 9) show that secondary nucleation is needed to sustain the frazil-ice population. We would also expect a continuous source of small seed crystals to have a similar effect, were the source sufficiently large. In calculations without nucleation, the crystals precipitate out, and the total concentration remains small, insufficient to affect the plume dynamics. The faster-growing crystals precipitate more over a shorter distance (dashed blue curve, panel e), because larger crystals rise faster and are more difficult to keep in suspension. After the frazil ice precipitates out of the plume, supercooling increases (blue curves, panel d), leading to a high rate of direct basal freezing (blue curves, panel g).

By contrast, a high nucleation rate triggers rapid growth of frazil ice,
which relieves the supercooling in the plume (red curves, panels c,
d, f). This behaviour is analogous to the frazil-ice
explosion we observed previously (Sect. 3.2), and it occurs
when secondary nucleation exceeds crystal removal by precipitation. The
increased frazil concentration leads to a more buoyant plume, causing it to
accelerate and have a slightly smaller depth *D* (red curves, panels
a, b). Precipitation is relatively unimportant (red curves,
panel e) as a result of a positive feedback: a faster-flowing plume
keeps crystals suspended more easily. Furthermore, nucleation produces small
crystals, which again are kept in suspension more easily. A faster crystal
growth rate is associated with a faster increase in crystal concentration
along the slope, although similar quasi-steady states are reached after the
supercooling is almost exhausted. As we found previously (Eq. 29), the quasi-steady ice concentration reflects the
overall energy balance of the system, rather than the growth dynamics.

The case of intermediate nucleation rate illustrates the surprising interplay between nucleation, growth, and precipitation of crystals. The calculation with a faster growth rate initially leads to a greater concentration of frazil ice, but the ice concentration is eventually overtaken by the slower growth rate calculation (green curves, panel c). Faster growth leads to larger crystals, which in turn are more readily precipitated (dashed green curve, panel e). This means that the crystal concentration eventually decreases, reducing the plume buoyancy and causing it to decelerate (dashed green curve, panel a). In this case, the plume thickness starts to increase rapidly as the plume begins to intrude at depth (dashed green curve, panel b). By contrast, the case of slower growth rate eventually reaches a crystal concentration comparable to the calculations with a higher crystal nucleation rate.

In terms of the large-scale dynamics, different parameterizations of crystal growth rate and nucleation can be the difference between a plume that is reinvigorated by frazil ice and reaches the end of the shelf and a plume that decelerates and intrudes at depth. This behaviour is likely to affect the ocean circulation and water mass transformation in the shelf seas around Antarctica. The differences between models could in principle be observed by considering the amount of frazil precipitation relative to basal freezing. The total amount of frazil formation also differs between the models (panels c, f). These differences are surprising: faster growth can lead to less frazil-ice formation in total if it is removed from suspension before it can multiply. This suggests that small-scale frazil-ice processes, which are hard to constrain in models, can have major implications for our understanding of the dynamics of plumes of ISW beneath Antarctica's floating ice shelves.

5 Conclusions

Back to toptop
The theory of frazil-ice dynamics pioneered by Daly (1984) encompasses the nucleation, growth, and removal of frazil ice. It describes the evolution of the size distribution of a population of crystals. We have applied this theory to understand ice formation in a supercooled ocean mixed layer and in a plume of ISW underneath a floating ice shelf. Understanding frazil-ice processes is significant to our understanding of ice–ocean interaction in the earliest, most explosive phase of ice growth. We have identified critical conditions for a self-sustained frazil-ice explosion, which occurs when secondary nucleation exceeds crystal removal. Crystal growth rate affects such explosions by changing the crystal size distribution and also alters the transient evolution of frazil ice, promoting faster increases in frazil concentration. We determined steady-state crystal size distributions and found that these were relatively insensitive to crystal growth rate but were sensitive to secondary nucleation and crystal removal. Thus measurement of crystal sizes could be used to estimate the nucleation rate indirectly. Finally, we showed that the parameterization of crystal growth rate and nucleation can dramatically affect the fate of plumes of supercooled ice shelf water, with implications for ice accretion on ice shelves and ocean circulation. Although our understanding of crystal growth rate has advanced recently, our understanding of crystal nucleation remains limited. Our calculations suggest that this is potentially a significant uncertainty, and it is a topic ripe for future research.

Code availability

Back to toptop
Code availability.

Please see https://github.com/davidreesjones/frazil-dynamics for software code to reproduce calculations and figures in the paper (Rees Jones, 2017).

Competing interests

Back to toptop
Competing interests.

The authors declare that they have no conflict of interest.

Acknowledgements

Back to toptop
Acknowledgements.

We thank Grae Worster for comments on an earlier version of this
paper and Adrian Jenkins for discussing previous models of
frazil-laden plumes. This publication arises from research funded by the John Fell Oxford University Press (OUP) Research Fund,
and Andrew J. Wells also acknowledges financial support through the research programme of the European Union FP7 award PCIG13-GA-2013-618610 SEA-ICE-CFD.
David W. Rees Jones acknowledges research funding through the NERC Consortium Grant NE/M000427/1 and NERC Standard Grant NE/I026995/1.
We would like to thank the Isaac Newton Institute for Mathematical Sciences for its hospitality during the programme Melt in
the Mantle, which was supported by EPSRC grant number EP/K032208/1.

Edited by: Christian Haas

Reviewed by: two anonymous referees

References

Back to toptop
Bonnecaze, R. T., Huppert, H. E., and Lister, J. R.: Particle-driven gravity currents, J. Fluid Mech., 250, 339–369, https://doi.org/10.1017/S002211209300148X, 1993. a

Carstens, T.: Experiments with supercooling and ice formation in flowing water, Geofys. Publ. Norway, 26, 3–18, 1966. a, b, c

Daly, S. F.: Frazil ice dynamics, CRREL Monograph, 84, 46 pp., 1984. a, b, c, d, e, f

Daly, S. F.: Report on frazil ice, Tech. Rep. 94-23, USA Cold Regions Research and Engineering Laboratory, CRREL Special Report, Hanover, New Hampshire, USA, 1994. a

Engelhardt, H. and Determann, J.: Borehole evidence for a thick layer of basal ice in the central Ronne Ice Shelf, Nature, 327, 318–319, https://doi.org/10.1038/327318a0, 1987. a

Fujioka, T. and Sekerka, R. F.: Morphological stability of disc crystals, J. Cryst. Growth, 24–25, 84–93, https://doi.org/10.1016/0022-0248(74)90284-X, 1974. a

Galton-Fenzi, B. K., Hunter, J. R., Coleman, R., Marsland, S. J., and Warner, R. C.: Modeling the basal melting and marine ice accretion of the Amery Ice Shelf, J. Geophys. Res., 117, C09031, https://doi.org/10.1029/2012JC008214, 2012. a, b, c

Godlovitch, D., Illner, R., and Monahan, A.: Smoluchowski coagulation models of sea ice thickness distribution dynamics, J. Geophys. Res., 116, C12005, https://doi.org/10.1029/2011JC007125, 2011. a

Gough, A., Mahoney, A., Langhorne, P., Williams, M., Robinson, N., and Haskell, T.: Signatures of supercooling: McMurdo Sound platelet ice, J. Glaciol., 58, 38–50, https://doi.org/10.3189/2012JoG10J218, 2012. a

Hanley, T. O. and Tsang, G.: Formation and properties of frazil in saline water, Cold Reg. Sci. Technol., 8, 209–221, https://doi.org/10.1016/0165-232X(84)90052-1, 1984. a

Heorton, H. D. B. S., Radia, N., and Feltham, D. L.: A model of sea ice formation in leads and polynyas, J. Phys. Oceanogr., 47, 1701–1718, https://doi.org/10.1175/JPO-D-16-0224.1, 2017. a

Holland, D. M. and Jenkins, A.: Modeling thermodynamic ice–ocean interactions at the base of an ice shelf, J. Phys. Oceanogr., 29, 1787–1800, https://doi.org/10.1175/1520-0485(1999)029<1787:MTIOIA>2.0.CO;2, 1999. a, b

Holland, P. R. and Feltham, D. L.: Frazil dynamics and precipitation in a water column with depth-dependent supercooling, J. Fluid Mech., 530, 101–124, https://doi.org/10.1017/S002211200400285X, 2005. a, b, c, d

Holland, P. R., Feltham, D. L., and Daly, S. F.: On the Nusselt number for frazil ice growth – a correction to “Frazil evolution in channels” by Lars Hammar and Hung-Tao Shen, J. Hydraul. Res., 45, 421–424, https://doi.org/10.1080/00221686.2007.9521775, 2007. a

Horvat, C. and Tziperman, E.: A prognostic model of the sea-ice floe size and thickness distribution, The Cryosphere, 9, 2119–2134, https://doi.org/10.5194/tc-9-2119-2015, 2015. a

Jenkins, A.: A one-dimensional model of ice shelf-ocean interaction, J. Geophys. Res., 96, 20671–20677, https://doi.org/10.1029/91JC01842, 1991. a

Jenkins, A. and Bombosch, A.: Modeling the effects of frazil ice crystals on the dynamics and thermodynamics of Ice Shelf Water plumes, J. Geophys. Res., 100, 6967–6981, https://doi.org/10.1029/94JC03227, 1995. a, b, c, d

Jordan, J. R., Holland, P. R., Jenkins, A., Piggott, M. D., and Kimura, S.: Modeling ice-ocean interaction in ice-shelf crevasses, J. Geophys. Res., 119, 995–1008, https://doi.org/10.1002/2013JC009208, 2014. a, b

Jordan, J. R., Kimura, S., Holland, P. R., Jenkins, A., and Piggott, M. D.: On the Conditional Frazil Ice Instability in Seawater, J. Phys. Oceanogr., 45, 1121–1138, https://doi.org/10.1175/JPO-D-14-0159.1, 2015. a, b

Khazendar, A. and Jenkins, A.: A model of marine ice formation within Antarctic ice shelf rifts, J. Geophys. Res., 108, 3235, https://doi.org/10.1029/2002JC001673, 2003. a, b

Langhorne, P. J., Hughes, K. G., Gough, A. J., Smith, I. J., Williams, M. J. M., Robinson, N. J., Stevens, C. L., Rack, W., Price, D., Leonard, G. H., Mahoney, A. R., Haas, C., and Haskell, T. G.: Observed platelet ice distributions in Antarctic sea ice: An index for ocean-ice shelf heat flux, Geophys. Res. Lett., 42, 5442–5451, https://doi.org/10.1002/2015GL064508, 2015. a

Lawson, D. E., Strabser, J. C., Evenson, E. B., Alley, R. B., Larson, G. J., and Arcone, S. A.: Glaciohydraulic supercooling: a freeze-on mechanism to create stratified, debris-rich basal ice: I. Field evidence, J. Glaciol., 44, 547–563, https://doi.org/10.3189/S0022143000002069, 1998. a

Lewis, E. L. and Perkin, R. G.: Ice pumps and their rates, J. Geophys. Res., 91, 11756–11762, https://doi.org/10.1029/JC091iC10p11756, 1986. a

Martin, S. and Kauffman, P.: A field and laboratory study of wave damping by grease ice, J. Glaciol., 27, 283–313, 1981. a

McFarlane, V., Loewen, M., and Hicks, F.: Laboratory measurements of the rise velocity of frazil ice particles, Cold Reg. Sci. Technol., 106–107, 120–130, https://doi.org/10.1016/j.coldregions.2014.06.009, 2014. a, b

Michel, B.: Theory of formation and deposit of frazil ice, Eastern Snow Conference, Proc. Annual Meeting, Quebec, 1963. a, b, c

Millero, F. J. and Leung, W. H.: The thermodynamics of seawater at one atmosphere, Am. J. Sci., 276, 1035–1077, https://doi.org/10.2475/ajs.276.9.1035, 1976. a

Oerter, H., Kipfstuhl, J., Determann, J., Miller, H., Wagenbach, D., Minikin, A., and Graft, W.: Evidence for basal marine ice in the Filchner-Ronne ice shelf, Nature, 358, 399–401, https://doi.org/10.1038/358399a0, 1992. a

Rees Jones, D. W.: Frazil-dynamics code, https://doi.org/10.5281/zenodo.1063501, https://github.com/davidreesjones/frazil-dynamics, 2017. a, b, c

Rees Jones, D. W. and Wells, A. J.: Solidification of a disk-shaped crystal from a weakly supercooled binary melt, Phys. Rev. E, 92, 022406, https://doi.org/10.1103/PhysRevE.92.022406, 2015. a, b, c, d, e, f, g, h

Skogseth, R., Nilsen, F., and Smedsrud, L. H.: Supercooled water in an Arctic polynya: Observations and modeling, J. Glaciol., 55, 43–52, https://doi.org/10.3189/002214309788608840, 2009. a

Smedsrud, L. H.: A model for entrainment of sediment into sea ice by aggregation between frazil-ice crystals and sediment grains, J. Glaciol., 48, 51–61, https://doi.org/10.3189/172756502781831520, 2002. a, b

Smedsrud, L. H. and Jenkins, A.: Frazil ice formation in an ice shelf water plume, J. Geophys. Res., 109, C03025, https://doi.org/10.1029/2003JC001851, 2004. a, b, c, d, e, f, g, h, i, j

Smedsrud, L. H. and Martin, T.: Grease ice in basin-scale sea-ice ocean models, Ann. Glaciol., 56, 295–306, https://doi.org/10.3189/2015AoG69A765, 2015. a, b

Svensson, U. and Omstedt, A.: Simulation of supercooling and size distribution in frazil ice dynamics, Cold Reg. Sci. Technol., 22, 221–233, https://doi.org/10.1016/0165-232X(94)90001-9, 1994. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o, p

Svensson, U. and Omstedt, A.: Numerical simulations of frazil ice dynamics in the upper layers of the ocean, Cold Reg. Sci. Technol., 28, 29–44, https://doi.org/10.1016/S0165-232X(98)00011-1, 1998. a, b, c

Thorndike, A.: Sea ice thickness as a stochastic process, J. Geophys. Res., 105, 1311–1313, https://doi.org/10.1029/1999JC900271, 2000. a

Toppaladoddi, S. and Wettlaufer, J. S.: Theory of the sea ice thickness distribution, Phys. Rev. Lett., 115, 148–501, https://doi.org/10.1103/PhysRevLett.115.148501, 2015. a

Wilchinsky, A. V., Heorton, H. D. B. S., Feltham, D. L., and Holland, P. R.: Study of the impact of ice formation in leads upon the sea ice pack mass balance using a new frazil and grease ice parameterization, J. Phys. Oceanogr., 45, 2025–2047, https://doi.org/10.1175/JPO-D-14-0184.1, 2015. a, b

Short summary

Frazil or granular ice grows rapidly from turbulent water cooled beneath its freezing temperature. We analyse numerical models of a population of ice crystals to provide insight into the treatment of frazil ice in large-scale models and hence in the environment. We determine critical conditions for explosively rapid frazil growth. We show that frazil-ice processes impact whether a plume of ice shelf water beneath an Antarctic ice shelf intrudes at depth or reaches the end of the shelf.

Frazil or granular ice grows rapidly from turbulent water cooled beneath its freezing...

The Cryosphere

An interactive open-access journal of the European Geosciences Union