Landfast sea ice material properties derived from ice bridge simulations using the Maxwell elasto-brittle rheology

• Published in: The Cryosphere
• Main Authors: M. Plante, B. Tremblay, M. Losch, J.-F. Lemieux
• Format: Article in Journal/Newspaper
• Language: English
• Published: Copernicus Publications 2020
• Subjects: geo; envir; Arctic; Sea ice; The Cryosphere
• Online Access: https://doi.org/10.5194/tc-14-2137-2020; https://tc.copernicus.org/articles/14/2137/2020/tc-14-2137-2020.pdf; https://doaj.org/article/0f4fd09506534a16b10752c494bad83f

The Maxwell elasto-brittle (MEB) rheology is implemented in the Eulerian finite-difference (FD) modeling framework commonly used in classical viscous-plastic (VP) models. The role of the damage parameterization, the cornerstone of the MEB rheology, in the formation and collapse of ice arches and ice bridges in a narrow channel is investigated. Ice bridge simulations are compared with observations to derive constraints on the mechanical properties of landfast sea ice. Results show that the overall dynamical behavior documented in previous MEB models is reproduced in the FD implementation, such as the localization of the damage in space and time and the propagation of ice fractures in space at very short timescales. In the simulations, an ice arch is easily formed downstream of the channel, sustaining an ice bridge upstream. The ice bridge collapses under a critical surface forcing that depends on the material cohesion. Typical ice arch conditions observed in the Arctic are best simulated using a material cohesion in the range of 5–10 kN m−2. Upstream of the channel, fracture lines along which convergence (ridging) takes place are oriented at an angle that depends on the angle of internal friction. Their orientation, however, deviates from the Mohr–Coulomb theory. The damage parameterization is found to cause instabilities at large compressive stresses, which prevents the production of longer-term simulations required for the formation of stable ice arches upstream of the channel between these lines of fracture. Based on these results, we propose that the stress correction scheme used in the damage parameterization be modified to remove numerical instabilities.