A second-order accurate in time IMplicit-EXplicit (IMEX) integration scheme for sea ice dynamics

Current sea ice models use numerical schemes based on a splitting in time between the momentum and continuity equations. Because the ice strength is explicit when solving the momentum equation, this can create unrealistic ice stress gradients when using a large time step. As a consequence, noise develops in the numerical solution and these models can even become numerically unstable at high resolution. To resolve this issue, we have implemented an iterated IMplicit-EXplicit (IMEX) time integration method. This IMEX method was developed in the framework of an already implemented Jacobian-free Newton-Krylov solver. The basic idea of this IMEX approach is to move the explicit calculation of the sea ice thickness and concentration inside the Newton loop such that these tracers evolve during the implicit integration. To obtain second-order accuracy in time, we have also modified the explicit time integration to a second-order Runge-Kutta approach and by introducing a second-order backward difference method for the implicit integration of the momentum equation. These modifications to the code are minor and straightforward. By comparing results with a reference solution obtained with a very small time step, it is shown that the approximate solution is second-order accurate in time. The new method permits to obtain the same accuracy as the splitting in time but by using a time step that is 10 times larger. Results show that the second-order scheme is more than five times more computationally efficient than the splitting in time approach for an equivalent level of error.

[1]  Ron Kwok,et al.  Elastic‐decohesive constitutive model for sea ice , 2006 .

[2]  David M. Holland,et al.  A comparison of the Jacobian-free Newton-Krylov method and the EVP model for solving the sea ice momentum equation with a viscous-plastic formulation: A serial algorithm study , 2012, J. Comput. Phys..

[3]  R. Reynolds,et al.  The NCEP/NCAR 40-Year Reanalysis Project , 1996, Renewable Energy.

[4]  Peter Lemke,et al.  Results of the Sea Ice Model Intercomparison Project: Evaluation of sea ice rheology schemes for use in climate simulations , 2000 .

[5]  Homer F. Walker,et al.  Choosing the Forcing Terms in an Inexact Newton Method , 1996, SIAM J. Sci. Comput..

[6]  William H. Lipscomb,et al.  Ridging, strength, and stability in high-resolution sea ice models , 2007 .

[7]  Jean-François Lemieux,et al.  Numerical convergence of viscous‐plastic sea ice models , 2009 .

[8]  Z. Bai,et al.  A globally convergent Newton-GMRES method for large sparse systems of nonlinear equations , 2007 .

[9]  Jean-Marc Molines,et al.  Evaluation of high-resolution sea ice models on the basis of statistical and scaling properties of Arctic sea ice drift and deformation , 2009 .

[10]  Sylvain Bouillon,et al.  A new modeling framework for sea-ice mechanics based on elasto-brittle rheology , 2011, Annals of Glaciology.

[11]  Elizabeth C. Hunke,et al.  Viscous–Plastic Sea Ice Dynamics with the EVP Model: Linearization Issues , 2001 .

[12]  Keguang Wang,et al.  Modeling linear kinematic features in pack ice , 2009 .

[13]  D. Menemenlis,et al.  On the formulation of sea-ice models. Part 1: Effects of different solver implementations and parameterizations , 2010 .

[14]  Dale R. Durran,et al.  Implicit–Explicit Multistep Methods for Fast-Wave–Slow-Wave Problems , 2012 .

[15]  William D. Hibler,et al.  Modeling Pack Ice as a Cavitating Fluid , 1992 .

[16]  Tamara G. Kolda,et al.  An overview of the Trilinos project , 2005, TOMS.

[17]  John K. Dukowicz,et al.  Comments on ''Stability of the Viscous-Plastic Sea Ice Rheology'' , 1997 .

[18]  Sylvain Bouillon,et al.  The elastic–viscous–plastic method revisited , 2013 .

[19]  Pierre Rampal,et al.  Scaling properties of sea ice deformation from buoy dispersion analysis , 2008 .

[20]  Yousef Saad,et al.  A Flexible Inner-Outer Preconditioned GMRES Algorithm , 1993, SIAM J. Sci. Comput..

[21]  Hrvoje Jasak,et al.  A strength implicit correction scheme for the viscous-plastic sea ice model , 2004 .

[22]  Rick M. Rauenzahn,et al.  A second order self-consistent IMEX method for radiation hydrodynamics , 2010, J. Comput. Phys..

[23]  D. Keyes,et al.  Jacobian-free Newton-Krylov methods: a survey of approaches and applications , 2004 .

[24]  Paul F. Tupper,et al.  Improving the numerical convergence of viscous-plastic sea ice models with the Jacobian-free Newton-Krylov method , 2010, J. Comput. Phys..

[25]  William D. Hibler,et al.  On an efficient numerical method for modeling sea ice dynamics , 1997 .

[26]  W. Hibler.  A Dynamic Thermodynamic Sea Ice Model , 1979 .

[27]  Jean-François Lemieux,et al.  A parallel Jacobian-free Newton-Krylov solver for a coupled sea ice-ocean model , 2014, J. Comput. Phys..

[28]  Dana A. Knoll,et al.  A fully second order implicit/explicit time integration technique for hydrodynamics plus nonlinear heat conduction problems , 2010, J. Comput. Phys..

[29]  L. Mysak,et al.  Modeling Sea Ice as a Granular Material, Including the Dilatancy Effect , 1997 .

[30]  Stephen J. Thomas,et al.  Using the preconditioned Generalized Minimum RESidual (GMRES) method to solve the sea‐ice momentum equation , 2008 .

[31]  William H. Lipscomb,et al.  Modeling Sea Ice Transport Using Incremental Remapping , 2004 .