keff in the Sea Ice Zone

Processes leading turbulence production in the sea ice zone.

Processes leading turbulence production in the sea ice zone.


Version 2 has been released . This version fixes an inconsistency in the way that we represented the gas transfer velocity from fetch-limited waves. Field data from 2015 and 2016 have given new insights into the TKE dissipation caused by these short period waves and the new formulation for keff is based on an entirely consistent forcing from shear, convection and fetch-limited waves. This eliminates a pre-existing bias in open water and at high wind speeds.

Download the new version now!!

The Matlab routines contained in this distribution calculate the effective gas transfer velocity (keff) for sea-ice-covered waters where turbulence can arise from (1) ice-water drag, (2) air-water drag and (3) buoyant convection. The routines are based on the parameter model described in:

Reference: Loose, B., W. R. McGillis, D. Perovich, C. J. Zappa, and P. Schlosser (2014), A parameter model of gas exchange for the seasonal sea ice zone. Ocean Science, 10(4), doi:10.5194/os-10-1-2014.


The program Saltywatr1.0.0 is collection of Matlab routines that implement a numerical solution of the shallow-water flow equations for rivers, estuaries, and the coastal sea. Imagine a two-dimensional planar flow, as in a shallow estuary, with spatial variations in velocity, water height, and salt. The model calculates the above flow variables using the conservation equations for water, salt, and momentum. To simulate flow driven by differ- ences in salt concentration, the momentum equation is coupled to the conservation of salt through the hydrostatic pressure term. Saltywatr is capable of modeling domains of irreg- ular geometry, such as estuaries and rivers.

Download: Saltywatr1.0.0.tar
Manual: Saltywatr_manual