Sea ice and the ocean mixed layer over the Antarctic shelf seas

An ocean mixed-layer model has been incorporated into the Los Alamos sea ice model CICE to investigate regional variations in the surface-driven formation of Antarctic shelf waters. This model captures well the expected sea ice thickness distribution, and produces deep (> 500 m) mixed layers in the Weddell and Ross shelf seas each winter. This results in the complete destratification of the water column in deep southern coastal regions leading to high-salinity shelf water (HSSW) formation, and also in some shallower regions (no HSSW formation) of these seas. Shallower mixed layers are produced in the Amundsen and Bellingshausen seas. By deconstructing the surface processes driving the mixed-layer depth evolution, we show that the net salt flux from sea ice growth/melt dominates the evolution of the mixed layer in all regions, with a smaller contribution from the surface heat flux and a negligible input from wind stress. The Weddell and Ross shelf seas receive an annual surplus of mixing energy at the surface; the Amundsen shelf sea energy input in autumn/winter is balanced by energy extraction in spring/summer; and the Bellingshausen shelf sea experiences an annual surface energy deficit, through both a low energy input in autumn/winter and the highest energy loss in spring/summer. An analysis of the sea ice mass balance demonstrates the contrasting mean ice growth, melt and export in each region. The Weddell and Ross shelf seas have the highest annual ice growth, with a large fraction exported northwards each year, whereas the Bellingshausen shelf sea experiences the highest annual ice melt, driven by the advection of ice from the northeast. A linear regression analysis is performed to determine the link between the autumn/winter mixed-layer deepening and several atmospheric variables. The Weddell and Ross shelf seas show stronger spatial correlations (temporal mean – intra-regional variability) between the autumn/winter mixed-layer deepening and several atmospheric variables compared to the Amundsen and Bellingshausen. In contrast, the Amundsen and Bellingshausen shelf seas show stronger temporal correlations (shelf sea mean – interannual variability) between the autumn/winter mixed-layer deepening and several atmospheric variables

Continuum sea ice rheology determined from subcontinuum mechanics

[1] A method is presented to calculate the continuum-scale sea ice stress as an imposed, continuum-scale strain-rate is varied. The continuum-scale stress is calculated as the area-average of the stresses within the floes and leads in a region (the continuum element). The continuum-scale stress depends upon: the imposed strain rate; the subcontinuum scale, material rheology of sea ice; the chosen configuration of sea ice floes and leads; and a prescribed rule for determining the motion of the floes in response to the continuum-scale strain-rate. We calculated plastic yield curves and flow rules associated with subcontinuum scale, material sea ice rheologies with elliptic, linear and modified Coulombic elliptic plastic yield curves, and with square, diamond and irregular, convex polygon-shaped floes. For the case of a tiling of square floes, only for particular orientations of the leads have the principal axes of strain rate and calculated continuum-scale sea ice stress aligned, and these have been investigated analytically. The ensemble average of calculated sea ice stress for square floes with uniform orientation with respect to the principal axes of strain rate yielded alignment of average stress and strain-rate principal axes and an isotropic, continuum-scale sea ice rheology. We present a lemon-shaped yield curve with normal flow rule, derived from ensemble averages of sea ice stress, suitable for direct inclusion into the current generation of sea ice models. This continuum-scale sea ice rheology directly relates the size (strength) of the continuum-scale yield curve to the material compressive strength

Continuum sea ice rheology determined from subcontinuum mechanics

[1] A method is presented to calculate the continuum-scale sea ice stress as an imposed, continuum-scale strain-rate is varied. The continuum-scale stress is calculated as the area-average of the stresses within the floes and leads in a region (the continuum element). The continuum-scale stress depends upon: the imposed strain rate; the subcontinuum scale, material rheology of sea ice; the chosen configuration of sea ice floes and leads; and a prescribed rule for determining the motion of the floes in response to the continuum-scale strain-rate. We calculated plastic yield curves and flow rules associated with subcontinuum scale, material sea ice rheologies with elliptic, linear and modified Coulombic elliptic plastic yield curves, and with square, diamond and irregular, convex polygon-shaped floes. For the case of a tiling of square floes, only for particular orientations of the leads have the principal axes of strain rate and calculated continuum-scale sea ice stress aligned, and these have been investigated analytically. The ensemble average of calculated sea ice stress for square floes with uniform orientation with respect to the principal axes of strain rate yielded alignment of average stress and strain-rate principal axes and an isotropic, continuum-scale sea ice rheology. We present a lemon-shaped yield curve with normal flow rule, derived from ensemble averages of sea ice stress, suitable for direct inclusion into the current generation of sea ice models. This continuum-scale sea ice rheology directly relates the size (strength) of the continuum-scale yield curve to the material compressive strength

Sea ice–ocean feedbacks in the Antarctic shelf seas

Observed changes in Antarctic sea ice are poorly understood, in part due to the complexity of its interactions with the atmosphere and ocean. A highly simplified, coupled sea ice–ocean mixed layer model has been developed to investigate the importance of sea ice–ocean feedbacks on the evolution of sea ice and the ocean mixed layer in two contrasting regions of the Antarctic continental shelf ocean: the Amundsen Sea, which has warm shelf waters, and the Weddell Sea, which has cold and saline shelf waters. Modeling studies where we deny the feedback response to surface air temperature perturbations show the importance of feedbacks on the mixed layer and ice cover in the Weddell Sea to be smaller than the sensitivity to surface atmospheric conditions. In the Amundsen Sea the effect of surface air temperature perturbations on the sea ice are opposed by changes in the entrainment of warm deep waters into the mixed layer. The net impact depends on the relative balance between changes in sea ice growth driven by surface perturbations and basal-driven melting. The changes in the entrainment of warm water in the Amundsen Sea were found to have a much larger impact on the ice volume than perturbations in the surface energy budget. This creates a net negative ice albedo feedback in the Amundsen Sea, reversing the sign of this typically positive feedback mechanism

The effect of a new drag-law parameterization on ice shelf water plume dynamics

A drag law accounting for Ekman rotation adjacent to a flat, horizontal boundary is proposed for use in a plume model that is written in terms of the depth-mean velocity. The drag law contains a variable turning angle between the mean velocity and the drag imposed by the turbulent boundary layer. The effect of the variable turning angle in the drag law is studied for a plume of ice shelf water (ISW) ascending and turning beneath an Antarctic ice shelf with draft decreasing away from the grounding line. As the ISW plume ascends the sloping ice shelf-ocean boundary, it can melt the ice shelf, which alters the buoyancy forcing driving the plume motion. Under these conditions, the typical turning angle is of order -10° over most of the plume area for a range of drag coefficients (the minus sign arises for the Southern Hemisphere). The rotation of the drag with respect to the mean velocity is found to be significant if the drag coefficient exceeds 0.003; in this case the plume body propagates farther along and across the base of the ice shelf than a plume with the standard quadratic drag law with no turning angle. © 2007 American Meteorological Society

The rate- and state-dependence of sea ice friction

Ice loading on offshore structures is often governed by the behaviour of granular ice flows around the structures. This behaviour is in turn governed by frictional interactions between ice floes. Observed stick-slip behaviour in ice friction suggests a dependence on slip history, which is not accounted for in current ice friction models. We propose, by analogy with results from rock friction, a rate- and state-model for ice friction. In this paper we present the results of a series of metre-scale ice basin experiments, in which we determine this rate- and state- dependence. We then propose a simple parameterisation which gives a general rate- and state- model for sea ice friction. The model suggests that when the slip rate increases, the ice-ice friction reaches a peak before decaying to some steady-state value. To assess peak frictional loading, therefore, the slip history must be considered as well as the current slip velocity

A multithickness sea ice model accounting for sliding friction

A multithickness sea ice model explicitly accounting for the ridging and sliding friction contributions to sea ice stress is developed. Both ridging and sliding contributions depend on the deformation type through functions adopted from the Ukita and Moritz kinematic model of floe interaction. In contrast to most previous work, the ice strength of a uniform ice sheet of constant ice thickness is taken to be proportional to the ice thickness raised to the 3/2 power, as is revealed in discrete element simulations by Hopkins. The new multithickness sea ice model for sea ice stress has been implemented into the Los Alamos "CICE" sea ice model code and is shown to improve agreement between model predictions and observed spatial distribution of sea ice thickness in the Arctic. © 2006 American Meteorological Society

The positive effect of mirror visual feedback on arm control in children with Spastic Hemiparetic Cerebral Palsy is dependent on which arm is viewed

Mirror visual feedback has previously been found to reduce disproportionate interlimb variability and neuromuscular activity in the arm muscles in children with Spastic Hemiparetic Cerebral Palsy (SHCP). The aim of the current study was to determine whether these positive effects are generated by the mirror per se (i.e. the illusory perception of two symmetrically moving limbs, irrespective of which arm generates the mirror visual feedback) or by the visual illusion that the impaired arm has been substituted and appears to move with less jerk and in synchrony with the less-impaired arm (i.e. by mirror visual feedback of the less-impaired arm only). Therefore, we compared the effect of mirror visual feedback from the impaired and the less-impaired upper limb on the bimanual coupling and neuromuscular activity during a bimanual coordination task. Children with SHCP were asked to perform a bimanual symmetrical circular movement in three different visual feedback conditions (i.e. viewing the two arms, viewing only one arm, and viewing one arm and its mirror image), combined with two head orientation conditions (i.e. looking from the impaired and looking from the less-impaired body side). It was found that mirror visual feedback resulted in a reduction in the eccentric activity of the Biceps Brachii Brevis in the impaired limb compared to the condition with actual visual feedback from the two arms. More specifically, this effect was exclusive to mirror visual feedback from the less-impaired arm and absent when mirror visual feedback from the impaired arm was provided. Across conditions, the less-impaired arm was the leading limb, and the nature of this coupling was independent from visual condition or head orientation. Also, mirror visual feedback did not affect the intensity of the mean neuromuscular activity or the muscle activity of the Triceps Brachii Longus. It was concluded that the positive effects of mirror visual feedback in children with SHCP are not just the result of the perception of two symmetrically moving limbs. Instead, in order to induce a decrease in eccentric neuromuscular activity in the impaired limb, mirror visual feedback from the ‘unaffected’ less-impaired limb is required

A model for the consolidation of rafted sea ice

Rafting is one of the important deformation mechanisms of sea ice. This process is widespread in the north Caspian Sea, where multiple rafting produces thick sea ice features, which are a hazard to offshore operations. Here we present a one-dimensional, thermal consolidation model for rafted sea ice. We consider the consolidation between the layers of both a two-layer and a three-layer section of rafted sea ice. The rafted ice is assumed to be composed of layers of sea ice of equal thickness, separated by thin layers of ocean water. Results show that the thickness of the liquid layer reduced asymptotically with time, such that there always remained a thin saline liquid layer. We propose that when the liquid layer is equal to the surface roughness the adjacent layers can be considered consolidated. Using parameters representative of the north Caspian, the Arctic, and the Antarctic, our results show that for a choice of standard parameters it took under 15 h for two layers of rafted sea ice to consolidate. Sensitivity studies showed that the consolidation model is highly sensitive to the initial thickness of the liquid layer, the fraction of salt release during freezing, and the height of the surface asperities. We believe that further investigation of these parameters is needed before any concrete conclusions can be drawn about rate of consolidation of rafted sea ice features