Sea-ice strength and internal stresses from in-situ measurements

A knowledge of sea-ice dynamics is key to achieving proper Arctic climate simulations. The goal of this thesis is to lay the foundation for a three-dimensional ice dynamics model that considers both thermal and dynamic stresses. This thesis presents in-situ observations from buoys deployed in the Canadian Arctic Archipelago (CAA). In the CAA, the sea ice is landfast for approximately six months in winter. In a first step, we have analyzed in-situ data from a sea-ice stress buoy deployed in the Viscount Melville Sound. Results demonstrate that thermal stress is the dominant source of internal stress in the CAA with few short-lived important dynamic stress events caused by ice floe collisions in the free drift season prior to the landfast season. This is in contrast with similar internal sea-ice stress measurements made in the Arctic Ocean, in which both dynamic and thermal stresses are of similar magnitude. Prior to landfast ice onset, the thermal stresses are isotropic, as hypothesized in prior analyses of ice internal stress data measured in the Arctic Ocean. After landfast ice onset, however, the ther- mal stresses become anisotropic. Results from the buoy data, together with results from a 1.5D thermal stress model (forced with simulated internal sea-ice temperatures), demonstrate that the anisotropy in thermal stress arises from land confinement induced by the coastline in the direction of the short-axis of the channel. Results from the model are in good agreement with the observed stress in the direction of both principal stresses. They suggest that anisotropy in thermal stress could impact the mode of failure of sea ice in the CAA. The results also suggest that viscous creep stress relaxation is important and acts on time scale of several days, which is longer than the time scale (several hours) suggested from the previous measurements. In a second step, we derive estimates of the sea-ice compressive strength parameter (P∗) based on a simple force balance and known external forcing (surface air-ice and ice-ocean stresses) and whether sea ice drifts under the action of these external loads or not. Results from a proof of concept experiment using internally consistent data from a fully coupled ice-ocean model (the Regional Ice Operation Prediction System, RIOPS) demonstrate that it is indeed possible to estimate P∗, which is a known quantity in the model, from the simple force balance presented. When the same method is applied to in-situ observations and reanalysis data, the method only produces meaningful bounds of P∗ ( = 94.4 ± 4.4 kN/m2) when the pack ice is mostly composed of first-year ice with little multi-year ice present. This P∗ estimate is approximately three times the value currently used in the modelling community. This highlights the fact that the ice drift measured at a point may not be representative of the ice behaviour on average in a region. For instance, results suggest that larger tidal ocean currents in the region are enough to prevent a landfast ice cover to develop locally while the pack ice is mostly landfast on larger scales. The failures of the method when it is applied to some buoy data also suggest that there may be error in the surface forcing from the CGRF (Cana- dian Meteorological Center’s Global Deterministic Prediction System Reforecast). The analysis of these errors is left for future work. Results discussed in this thesis highlight the importance of thermal stresses in sea-ice models. The only forcing required by a thermal stress model is the internal temperature profile. Therefore, a thermal stress model could be implemented in current sea-ice thermodynamic models with a little effort. It would then be coupled with the dynamical part of the model by developing yield criteria for ice failure that are a function of the total (dynamical and thermal) stresses at a point rather than the depth average internal stress as currently done in the community.Cette thèse présente l’analyse d’observations enregistrées par des bouées déployées dans l’Archipel Arctique Canadien (AAC), où la glace de mer consiste en un couvert de glace de rive pendant approximativement six mois par année. En premier lieu, les contraintes internes enregistrées sur place par une bouée déployée dans le détroit du Vicomte de Melville sont analysées. Cette étude démontre que les contraintes thermiques dominent les forces internes de la banquise dans l’AAC, et ce malgré la présence de quelques évènements où de fortes contraintes dynamiques de courte durée ont été enregistrées. Ces résultats diffèrent des observations obtenues par des bouées similaires déployées dans l’océan Arctique, qui indiquent plutôt une équivalence en magnitude entre les contraintes de source thermique et de source dynamique. Pendant la période de dérive, les contraintes thermiques mesurées dans l’AAC présentent une isotropie qui concorde avec les analyses pour l’océan Arctique. Les contraintes thermiques deviennent cependant anisotropes après la transition vers le régime de glace de rive. En combinant ces résultats avec des simulations provenant d’un modèle de glace à 1.5 dimensions, il est démontré que cette anisotropie découle du confinement de la glace de rive par les berges du chenal. Les contraintes simulées par le modèle, dans chacun des axes principaux, concordent bien avec celles observées par les bouées. Les résultats de cette étude suggèrent que le mode de fracture de la glace de rive dans l’AAC peut être affecté par l’anisotropie des contraintes thermiques. De plus, il est observé que la relaxation visqueuse des contraintes due au fluage de la glace ne devient importante que sur une échelle temporelle de l’ordre de plusieurs jours, ce qui est plus long que l’ordre suggéré par les mesures en laboratoire. En second lieu, nous avons obtenu des valeurs approximatives pour le paramètre de résistance à la compression de la glace de mer (P*) à partir de l’équilibre des forces, en se basant sur l’état observé de la dérive de la glace (dérive ou glace de rive) et en utilisant un forçage connu. Afin de démontrer qu’il est possible de déterminer P* à partir du simple équilibre des forces, la méthode proposée est appliquée sur les données d’un modèle de glace entièrement couplé à l’océan (RIOPS), dans lequel P* est une variable connue. En appliquant cette méthode sur les données d’observations et de réanalyses, des limites physiques sont établies pour restreinte la plage de valeur de P* (= 94.4 ± 4.4 kN/m2) lorsque la glace est composée principalement de glace annuelle. Cette plage de valeur correspond approximativement au triple de la valeur couramment utilisée dans la communauté de modélisation. Les résultats cette étude suggèrent que les mesures de dérive à un point précis peuvent ne pas correspondre au comportement moyen de la glace dans une région. Par exemple, il est montré que dans une région principalement recouverte de glace de rive, la présence de courants de marée plus forts peut être suffisante pour prévenir le développement de la glace de rive localement. Dans le cas de certaines bouées, l’échec de la méthode proposée suggère la présence d’erreurs dans les forçage de surface du CGRF, et l’analyse de ces erreurs est laissé comme projet futur. Les résultats discutés dans cette thèse soulignent l’importance des contraintes thermiques pour la modélisation de la glace de mer. Le profil de température interne de la glace étant le seul forçage nécessaire dans un modèle de contraintes thermiques, il serait possible de l’implémenter dans un modèle thermodynamique de glace de mer. Ce modèle thermodynamique pourrait ensuite être couplé à une composante dynamique par le développement de limites de contraintes en fonctions des contraintes totales (dynamiques + thermiques) en un point, plutôt qu’en fonction des contraintes internes moyennées sur l’épaisseur, comme il est couramment utilisé dans la communauté