Repository logo

Infoscience

  • English
  • French
Log In
Logo EPFL, École polytechnique fédérale de Lausanne

Infoscience

  • English
  • French
Log In
  1. Home
  2. Academic and Research Output
  3. Journal articles
  4. Comparison of different numerical approaches to the 1D sea-ice thermodynamics problem
 
Loading...
Thumbnail Image
research article

Comparison of different numerical approaches to the 1D sea-ice thermodynamics problem

Dupont, Frederic
•
Vancoppenolle, Martin
•
Tremblay, Louis-Bruno
Show more
2015
Ocean Modelling

The vertical one-dimensional sea-ice thermodynamic problem using the principle of conservation of enthalpy is revisited here using (1) the Bitz and Lipscomb (1999) finite-difference approach (FD), (2) a reformulation of the sigma-level transformation of Huwald et al. (2005b) (FV) and (3) a Finite Element approach also in sigma coordinates (FE). These three formulations are compared in terms of physics, numerics, and performance, in order to identify the best choice for large-scale climate models. The BL99 formulation sequentially treats the diffusion of heat and the changes in the vertical position of the ice-snow layers. In contrast, the FV sigma-level transformation elegantly treats both simultaneously. The original FV formulation suffers however from slow convergence. The convergence can nonetheless be improved significantly with a few simple modifications to the original code. The three formulations are compared following the experimental protocol of the Sea Ice Model Intercomparison Project for ice thermodynamics (SIMIP2). It is found that all formulations converge to the same solution. The FD approach, however, suffers from the added cost of the remapping step at large number of ice layers we include in the appendix an optimized version of the FD code–written by one of the reviewer–that resolves this issue. Finally the FE formulation results in a sub-surface temperature over-estimation at low resolution, a problem which disappears at high resolution. Hence, only FD and FV are found suitable for climate models.

  • Files
  • Details
  • Metrics
Loading...
Thumbnail Image
Name

OCEMOD942.pdf

Type

Publisher's Version

Access type

openaccess

Size

1.02 MB

Format

Adobe PDF

Checksum (MD5)

1921838a9a1a2fdfe401b3d883f5b776

Loading...
Thumbnail Image
Name

mmc1.doc

Access type

openaccess

Size

14.5 KB

Format

Microsoft Word

Checksum (MD5)

af2b8cffce5d8edd91ee0b71a1a97372

Logo EPFL, École polytechnique fédérale de Lausanne
  • Contact
  • infoscience@epfl.ch

  • Follow us on Facebook
  • Follow us on Instagram
  • Follow us on LinkedIn
  • Follow us on X
  • Follow us on Youtube
AccessibilityLegal noticePrivacy policyCookie settingsEnd User AgreementGet helpFeedback

Infoscience is a service managed and provided by the Library and IT Services of EPFL. © EPFL, tous droits réservés