We use a sea ice model to reproduce ice growth observations from two buoys deployed on coastal sea ice and analyze the improvements brought by new physics that represent the presence of saline liquid water in the ice interior. We find that the new physics with default parameters degrade the model performance, with overly rapid ice growth and overly early snow flooding on top of the ice. The performance is largely improved by simple modifications to the ice growth and snow-flooding algorithms.

We developed a standard viscous–plastic sea-ice model based on the numerical framework called smoothed particle hydrodynamics. The model conforms to the theory within an error of 1 % in an idealized ridging experiment, and it is able to simulate stable ice arches. However, the method creates a dispersive plastic wave speed. The framework is efficient to simulate fractures and can take full advantage of parallelization, making it a good candidate to investigate sea-ice material properties.

In some shallow seas, grounded ice ridges contribute to stabilizing and maintaining a landfast ice cover. A scheme has already proposed where the keel thickness varies linearly with the mean thickness. Here, we extend the approach by taking into account the ice thickness and bathymetry distributions. The probabilistic approach shows a reasonably good agreement with observations and previous grounding scheme while potentially offering more physical insights into the formation of landfast ice.

Sea ice motion is a versatile parameter for monitoring the Arctic climate system. In this contribution, we use data from drifting buoys, winds, and ice thickness to parameterize the motion of sea ice in a free drift regime – i.e., flowing freely in response to the forcing from the winds and ocean currents. We show that including a dependence on sea ice thickness and taking into account a climatology of the surface ocean circulation significantly improves the accuracy of sea ice motion estimates.

We propose a generalized form for the damage parameterization such that super-critical stresses can return to the yield with different final sub-critical stress states. In uniaxial compression simulations, the generalization improves the orientation of sea ice fractures and reduces the growth of numerical errors. Shear and convergence deformations however remain predominant along the fractures, contrary to observations, and this calls for modification of the post-fracture viscosity formulation.

Deformations in the Arctic sea ice cover take the shape of narrow lines. High-resolution sea ice models recreate these deformation lines. Recent studies have shown that the most widely used sea ice model creates fracture lines with intersection angles larger than those observed and cannot create smaller angles. In our work, we change the way sea ice deforms post-fracture. This change allows us to understand the link between the sea ice model and intersection angles and create more acute angles.

Canada's coastlines include diverse ocean environments. In response to the strong need to support marine activities and security, we present the first pan-Canadian operational regional ocean analysis system. A novel online tidal harmonic analysis method is introduced that uses a sliding-window approach. Innovations are compared to those from the Canadian global analysis system. Particular improvements are found near the Gulf Stream due to the higher model grid resolution.

Sea ice plays a crucial role in global energy and water budgets. To get a better simulation of sea ice, we coupled a sea ice model with an atmospheric and ocean model to form a fully coupled system. The sea ice simulation results of this coupled system demonstrated that a two-way coupled model has better performance in terms of sea ice, especially in summer. This indicates that sea-ice–ocean–atmosphere interaction plays a crucial role in controlling Arctic summertime sea ice distribution.