loading page

Simulating Temperature Transients in First Year Sea Ice
  • Nitay Ben Shachar,
  • Mark McGuinness,
  • Joe Trodahl
Nitay Ben Shachar
Victoria University of Wellington, School of Mathematics and Statistics

Corresponding Author:[email protected]

Author Profile
Mark McGuinness
Victoria University of Wellington, School of Mathematics and Statistics
Author Profile
Joe Trodahl
Victoria University of Wellington, School of Chemical and Physical Sciences
Author Profile


A dynamic one-dimensional model of the freezing of first-year sea ice with a snow cover is presented along with Matlab code. The model is based on historical models in the literature, and is primarily intended to design a computational tool that accurately simulates transient temperature changes within growing sea ice. This may complement temperature data at various depths and times in growing sea ice. The model is validated by comparing the simulated ice temperature profile and growth rate to thermistor temperature profiles measured in 1997 in McMurdo Sound, Antarctica. This comparison shows evidence of platelet ice growth. Despite this, an accuracy of ±0.4 o C is achieved throughout the entire ice temperature profile, when comparing simulated transient temperatures to data. The model is also used, with the inclusion of a 15mm thick snow layer, to correct a previously observed discrepancy between ice growth rates measured via banding structures in ice cores and the predictions of a steady state ice growth model.