LIM2 is a two-level thermodynamic-dynamic sea ice model. Sensible heat storage and vertical heat conduction within snow and ice are determined by a three-layer model (one layer for snow and two layers for ice). The effect of the subgrid-scale snow and ice thickness distributions is implicitly accounted for through an effective thermal conductivity. The storage of latent heat inside the ice resulting from the trapping of shortwave radiation by brine pockets is taken into account. Vertical and lateral growth/decay rates of the ice are obtained from prognostic energy budgets at both the bottom and surface boundaries of the snow-ice cover and in leads. The surface albedo depends on the state of the surface, the snow and ice thicknesses, and the cloudiness. When the load of snow is large enough to depress the snow-ice interface under the water level, seawater is supposed to infiltrate the entirety of the submerged snow and to freeze there, forming snow ice. Ice dynamics are simulated by assuming that sea ice behaves as a two-dimensional viscous-plastic or elastic-viscous-plastic continuum in dynamical interaction with atmosphere and ocean. The ice momentum and transport equations are written in curvilinear, orthogonal coordinates and are numerically solved on a B- or C-grid. For the transport equation, a numerical method that conserves the second-order moments of the spatial distribution of the advected quantity is used. The code can be downloaded from the NEMO website



Please quote the following two references in all scientific publications using LIM2 :

Fichefet, T., and M.A. Morales Maqueda, 1997 : Sensitivity of a global sea ice model to the treatment of ice thermodynamics and dynamics. Journal of Geophysical Research, 102, 12,609-12,646, doi:10.1029/97JC00480.

Bouillon, S., M.A. Morales Maqueda, V. Legat, and T. Fichefet, 2009 : An elastic-viscous-plastic sea ice model formulated on Arakawa B and C grids. Ocean Modelling, 27, 174-184, doi : 10.1016/j.ocemod.2009.01.004.


Typical namelists :