The Effect Of Rheology On Seasonal Sea-Ice Simulations

Abstract

On the seasonal time scales relevant to numerical investigations of climate, the rheology used in large-scale sea ice models significantly affects the ice thickness build-up and ice velocity fields. Plastic rheologies with a normal flow rule have been used to-date in seasonal dynamic thermodynamic simulations. These rheologies have proved useful in simulating discontinuous slip near the coast while still supplying relatively robust velocity fields in the central Arctic Basin. However, as indicated by limited numerical sensitivity studies with different types of elliptical yield curves, the amount of shear strength significantly affects the ice build-up and can possibly cause a stoppage of the ice outflow through Fram Strait. In addition to the shear strength issue, there is also the possibility that non-normal flow rule rheologies, such as the Mohr Coulomb failure criterion used in soil mechanics, may cause somewhat different types of flow patterns, especially in the Fram Strait region. However, to date no seasonal simulations with such non-normal flow rule rheologies have been carried out. In order to investigate the role of different rheologies on the large-scale flow patterns in the Arctic Basin, a more general numerical scheme than that of Hibler (1979) is developed, which allows the simulation of the dynamic thermodynamic behavior of sea ice with a wide variety of different non linear rheologies. Using this numerical scheme, comparative simulations are carried out to seasonal equilibrium with several variations of the Mohr Coulomb rheology and compared to the more standard Elliptical yield curve results. In particular, the main control Mohr Coulomb case is a capped rheology in which the shear strength is taken to be proportional to the compressive stress. In this capped case only shear flow is allowed until a maximum allowable compressive stress is reached. This cap strength is parameterized to be a function of the ice thickness and compactness. For comparison, a simulation with a very large cap strength is also carried out, and an experiment with a similar compressive cap but much lower shear strength. Overall the results are analyzed to determine the sensitivity of the ice build-up to flow rule and shear strength magnitude. In addition special attention is given to the character of the flow and stoppage (if any) through Fram Strait.