This work aims to increase our understanding of the nature of large scale features of sea ice from a dynamics point of view.Sea ice plays an important part in the exchange of heat and humidity between sea and air and thus is an important component of the climate system. Its physical presence also directly impacts the various forms of life such as diatoms, polar bears and humans alike.The dynamics of sea ice affect both weather and climate, through the large scale drift in the Arctic from the Siberian coast towards Fram Strait, through creation of cracks in the ice called leads or polynyas, and through ridging and other mechanical deformations of ice floes.In this work, we have focused on modelling of ridged ice for a number of reasons. Direct observations of the internal ice state is very difficult to perform and in general, observations of sea ice are either sparse or of limited information density. Ridged ice can be seen as the memory of high ice stress events, giving us a view on these highly dynamic events. Ridging is of major importance for the ice thickness distribution, as the thickest ice can only be formed through mechanical processes. Further, ridged ice is of direct interest for anyone conducting shipping through seasonal or perennial ice covered seas as it can form impenetrable barriers or in extreme even cases crush a ship caught within the ice pack. To this end, a multi-category sea ice model, the HELsinki Multi category Ice model (HELMI), was implemented into the Rossby Centre Ocean model (RCO). HELMI has explicit formulations for ridged and rafted ice, as well as sub-grid scale ice thickness distribution (a feature shared with other multi category models) and an ice strength based on energetics. These features give RCO better representation of sub-grid scale physics and gives us the possibility to study the deformed ice in detail. In paper I we look at the change in behaviour in the Arctic as the ice becomes more mobile, leading to a slight increase in modelled ridged ice volume in the central Arctic, despite a general trend of a decreasing ice cover.Paper II takes us to the Baltic Sea and the possibilities of modelling ridge ice concentration with a statistical model.In Paper III we investigate how the diminishing ice cover in future scenarios affects the biological activity in the Baltic Sea.Finally Paper IV investigates how the ice stress and the internal ice force can be interpreted in terms of ice compression on the ship scale. / <p>At the time of the doctoral defence the following paper was unpublished and had a status as follows: Paper 4: Manuscript</p>

The aim of this research was to study the relationship between rock slope stability and glacial processes. An in-depth analysis of our current understanding of how glaciated rock slopes develop instability and movement during deglaciation is presented; this shows that understanding is incomplete without an appreciation of the variable mechanical behaviour of glacier ice. In this thesis, I argue that: (1) The ductile behaviour of ice at low strain rates allows movement of rock slopes buttressed by ice. Field evidence and simple force models are used to explore rate of movement of ice-contact slopes and the conditions under which they evolve. The results indicate that large rockslides can move and deform glacial ice at rates of 10-2 to 102 m-yr. This implies that ice-contact slope movement may be important for slope evolution and the erosion and entrainment processes of glaciers; and (2) the elastic strength of glacier ice at the high strain rates associated with seismic shaking enables ice to modify the response of the surrounding rock to seismic shaking. To explore this, numerical analyses of the interaction between glacial erosion, glacier mass, topography, and earthquake shaking intensity are undertaken. Shaking of mountains of variable shape and with different levels of ice inundation is simulated using FLAC 6.0. The results suggest that complete inundation by ice can significantly reduce shaking intensity. This, in combination with glacial steepening of slopes, may make recently deglaciated slopes more prone to coseismic failure. In the final chapter of the thesis, I present a conceptual model of the evolution of slope stability during stages of glaciation and deglaciation. The model incorporates the ideas presented in the thesis. I then offer recommendations for how our understanding of these processes can be further advanced.

De récentes analyses statistiques de données satellitales et de bouées dérivantes ont révélé le caractère hautement hétérogène et intermittent de la déformation de la banquise Arctique, démontrant de ce fait que le schéma rhéologique visco-plastique utilisé traditionnellement en modélisation climatique et opérationnelle ne simule pas adéquatement le comportement dynamique des glaces ainsi que les efforts mécaniques en leur sein.Un cadre rhéologique alternatif, baptisé "Maxwell-Élasto-Fragile" (Maxwell-EB) est donc développé dans le but de reproduire correctement la dérive et la déformation des glaces dans les modèles continus de la banquise à l'échelle régionale et globale. Le modèle se base en partie sur un cade de modélisation élasto-fragile utilisé pour les roches et la glace. Un terme de relaxation visqueuse est ajouté à la relation constitutive d'élasticité linéaire ainsi qu'une viscosité effective, ou "apparente", laquelle évolue en fonction du niveau d'endommagement local du matériel simulé, comme son module d'élasticité. Ce cadre rhéologique permet la dissipation partielle des contraintes internes par le biais de déformations permanentes, possiblement grandes, le long de failles (ou "leads") lorsque le matériel est fortement endommagé ainsi que la conservation de la mémoire des contraintes associées aux déformations élastiques dans les zones où le matériel reste relativement peu endommagé.The schéma numérique du modèle Maxwell-EB est basé sur des méthodes de calcul variationnel et par éléments finis. Une représentation Eulérienne des équations du mouvement est utilisée et des méthodes dites Galerkin discontinues sont implémentées pour le traitement des processus d'advection.Une première série de simulations idéalisées et sans advection est présentée, lesquelles démontrent que la rhéologie Maxwell-Élasto-Fragile reproduit les caractéristiques principales du comportement mécanique de la banquise, c'est-à-dire la localisation spatiale, l'anisotropie et l'intermittence de la déformation ainsi que les lois d'échelle qui en découlent. La représentation adéquate de ces propriétés de la déformation se traduit par la présence de très forts gradients au sein des champs de contrainte, de déformation et du niveau d'endommagement simulés par le modèle. Des tests visant à évaluer la diffusion numérique découlant de l'advection de ces gradients extrêmes ainsi qu'à identifier certaines contraintes numériques du modèle sont ensuite présentés. De premières simulations en grandes déformations, incluant les processus d'advection, sont réalisées, lesquelles permettent une comparaison aux résultats d'une expérience de Couette annulaire sur de la glace fabriquée en laboratoire. Le modèle reproduit en partie le comportement mécanique observé. Par ailleurs, les différences entre les résultats des simulations et ceux obtenus en laboratoire permettent d'identifier certaines limitations, numériques et physiques, du modèle en grandes déformations. Finalement, le modèle rhéologique est utilisé pour modéliser la dérive et la déformation des glaces à l'échelle de la banquise Arctique. Des simulations idéalisées de l'écoulement de glace dans un chenal étroit sont présentées. Le modèle simule une propagation localisée de l'endommagement, définissant des failles en forme d'arche, et la formation de ponts de glace stables. / In recent years, analyses of available ice buoy and satellite data have revealed the strong heterogeneity and intermittency of the deformation of sea ice and have demonstrated that the viscous-plastic rheology widely used in current climate models and operational modelling platforms does not simulate adequately the drift, deformation and mechanical stresses within the ice pack.A new alternative rheological framework named ''Maxwell-Elasto-Brittle” (Maxwell-EB) is therefore developed in the view of reproducing more accurately the drift and deformation of the ice cover in continuum sea ice models at regional to global scales. The model builds on an elasto-brittle framework used for ice and rocks. A viscous-like relaxation term is added to a linear-elastic constitutive relationship together with an effective viscosity that evolves with the local level of damage of the material, like its elastic modulus. This framework allows for part of the internal stress to dissipate in large, permanent deformations along the faults/leads once the material is highly damaged while retaining the memory of small, elastic deformations over undamaged areas. A healing mechanism is also introduced, counterbalancing the effects of damaging over large time scales.The numerical scheme for the Maxwell-EB model is based on finite elements and variational methods. The equations of motion are cast in the Eulerian frame and discontinuous Galerkin methods are implemented to handle advective processes.Idealized simulations without advection are first presented. These demonstrate that the Maxwell-EB rheological framework reproduces the main characteristics of sea ice mechanics and deformation : the strain localization, the anisotropy and intermittency of deformation and the associated scaling laws. The successful representation of these properties translates into very large gradients within all simulated fields. Idealized numerical experiments are conducted to evaluate the amount of numerical diffusion associated with the advection of these extreme gradients in the model and investigate other limitations of the numerical scheme. First large-deformation simulations are carried in the context of a Couette flow experiment, which allow a comparison with the result of a similar laboratory experiment performed on fresh-water ice. The model reproduces part of the mechanical behaviour observed in the laboratory. Comparison of the numerical and experimental results allow identifying some numerical and physical limitations of the model in the context of large-deformation and laboratory-scale simulations. Finally, the Maxwell-EB framework is implemented in the context of modelling the drift and deformation of sea ice on geophysical scales. Idealized simulations of the flow of sea ice through a narrow channel are presented. The model simulates the propagation of damage along arch-like features and successfully reproduces the formation of stable ice bridges.