Sea ice growth processes

Sea ice is a complex composite composed primarily of pure ice in various states of crystallization along with air bubbles and included pockets of brine. Understanding its growth processes is important both for climate scientists for use in simulations as well remote sensing specialists since the composition and microstructural properties of the ice ultimately affect how it interacts with electromagnetic radiation.

Sea ice growth models for predicting the ice distribution and extent are also valuable for shipping concerns. An ice growth model can be combined with remote sensing measurements in an assimilation model as a means of generating more accurate ice charts.

Overview

Several formation mechanisms of sea ice have been identified. At its earliest stages, sea ice
consists of elongated, randomly oriented crystals. This is called frazil and mixed with water in
the unconsolidated state is known as grease ice. If wave and wind conditions are calm, these
crystals will consolidate at the surface and by selective pressure, begin to grow preferentially
in the downward direction, forming nilas. In more turbulent conditions, the frazil will consolidate
by mechanical action to form pancake ice, which has a more random structure

Another common formation mechanism,
especially in the Antarctic where precipitation over sea ice is high, is from snow deposition: on
thin ice, the snow will weigh down the ice enough to cause flooding. Subsequent freezing will
form ice with a much more granular structure.

One of the more interesting processes to occur within consolidated ice packs is changes in
the saline content. As the ice freezes, most of the salt content gets rejected and forms highly
saline brine inclusions between the crystals. With decreasing temperatures in the ice sheet, the
size of the brine pockets decreases while the salt content goes up. Since ice is less dense than
water, increasing pressure causes some of the brine to be ejected from both the top and bottom,
producing the characteristic ‘C’-shaped salinity profile of first year ice.


Brine will also
drain through vertical channels, particularly in the melt season. Thus multi-year ice will tend
to have both lower salinity and lower density than first-year ice

Vertical growth

The downward growth of consolidated ice in calm conditions is determined by the rate of heat transfer, ''Q^*'', at the ice-water interface. If we assume that the ice is in thermal equilibrium both with itself and its surroundings and that the weather conditions are known, then we can determine ''Q^*'' by solving the following equation:

:$h Q^* = k (T_s - T_w)$

for ''T_s'', the surface temperature. The water temperature, ''T_w'', is assumed to be at or near freezing, while the ice thickness, ''h'', is assumed to be known and we can approximate the thermal conductivity, ''k'', as an average over the layers (which have different salinities) or simply use the value for pure ice. The net heat flux is a total over four components:

:
Q^* = Q_E \left (T_s) \right+ Q_H(T_s) + Q_(T_s^4) + Q_

which are latent, sensible, longwave and shortwave fluxes, respectively. For a description of the approximate parameterizations, see determining surface flux under sea ice thickness. The equation can be solved using a numerical root-finding algorithm such as bisection: the functional dependencies on surface temperature are given, with ''e'' being the equilibrium vapor pressure.

While Cox and Weeks


assume thermal equilibrium, Tonboe


uses a more complex thermodynamic model based on numerical solution of the heat equation. This would be appropriate when the ice is thick or the weather conditions are changing rapidly.

The rate of ice growth can be calculated from heat flux by the following equation:

:$g=\frac$

where ''L'' is the latent heat of fusion for water and $\rho$ is the density of ice. The growth rate in turn determines the saline content of the newly frozen ice. Empirical equations for determining the initial brine entrapment in sea ice have been derived by Cox and Weeks and Nakawo and Sinha


and take the form:

:$S=S_0 f(g)$

where ''S'' is ice salinity, ''S₀'' is the salinity of the parent water and ''f'' is an empirical function of ice growth rate, e.g.:

:$f(g) = \frac$

where ''g'' is in cm/s.

Salt content

Brine entrapped in sea ice will always be at or near freezing since any departure will either cause some of the water in the brine to freeze, or melt some of the surrounding ice. Thus, brine salinity is variable and can be determined based strictly on temperature—see freezing point depression. References and

contain empirical formulas relating sea ice temperature to brine salinity.

The relative brine volume, ''V_b'', is defined as the fraction of brine relative to the total volume. It too is highly variable, however its value is more difficult to determine since changes in temperature may cause some of the brine to be ejected or move within the layers, particularly in new ice. Writing equations relating the salt content of the brine, the total salt content, the brine volume, the density of the brine and the density of the ice and solving for brine volume produces the following relation:

:$V_b=\frac$

where ''S'' is sea ice salinity, ''S_b'' is brine salinity, $\rho_i$ is the density of the ice and $\rho_b$ is brine density. Compare with this empirical formula from Ulaby et al.:

:$V_b = 10^S\left(-\frac+0.532\right)$
where ''T'' is ice temperature in degrees Celsius and ''S'' is ice salinity in parts per thousand.

In new ice, the amount of brine ejected as the ice cools can be determined by assuming that the total volume stays constant and subtracting the volume increase from the brine volume. Note that this is only applicable to newly formed ice: any warming will tend to generate air pockets as the brine volume will increase more slowly than the ice volume decreases, again due to the density difference.
Cox and Weeks provide the following formula determining the ratio of total ice salinity between temperatures, ''T₁'' and ''T₂'' where ''T₂'' < ''T₁'':

:
\frac=\frac\frac
\exp \left \lbrace \frac{\rho_i \left _b(T_1) - S_b(T_2)\right \right \rbrace

where ''c''=0.8 kg m⁻³ is a constant.
As the ice goes through constant warming and cooling cycles, it becomes progressively more porous through ejection of the brine and drainage through the resulting channels.

The figure above shows a scatter plot of salinity versus ice thickness for ice cores taken from the Weddell Sea, Antarctica, with an exponential fit of the form, $S = \exp(a h + b)$, overlaid, where ''h'' is ice thickness and ''a'' and ''b'' are constants.

Horizontal motion

The horizontal motion of sea ice is quite difficult to model because ice is a non-Newtonian fluid.
Sea ice will deform primarily at fracture points which in turn will form at the points of greatest stress and lowest strength, or where the ratio between the two is a maximum. Ice thickness, salinity and porosity will all affect the strength of the ice. The motion of the ice is driven primarily by ocean currents, though to a lesser extent by wind. Note that stresses will not be in the direction of the winds or currents, but rather will be shifted by coriolis effects—see, for instance, Ekman spiral.

See also

* Sea ice
* Sea ice thickness
* Sea ice concentration
* Sea ice emissivity modelling

References

Sea ice
Climatology