This paper explores the use of the constant grid flux form forward-in-time upstream-biased advection schemes for the advection of temperature and salinity in ocean modeling. The constant grid flux form schemes are shown to be an improvement over the traditional central differencing commonly used in ocean models. In addition, nonoscillatory versions of the scheme, which employ flux limiters, are explored. The limiters are based on total variation diminishing concepts and are applied to higher-order (in space) versions of the constant grid flux form scheme. The constant grid flux form schemes are Crowley-type upstream-biased Eulerian advection schemes. They are mass conserving and possess small amplitude and phase errors. The flux limiters prevent the under- and overshooting associated with the numerical dispersion of the unlimited schemes. The limited schemes are easy to implement, efficient, and nonoscillatory. Of these schemes the third-order and fifth-order versions employing the PDM limiter ar... Abstract This paper explores the use of the constant grid flux form forward-in-time upstream-biased advection schemes for the advection of temperature and salinity in ocean modeling. The constant grid flux form schemes are shown to be an improvement over the traditional central differencing commonly used in ocean models. In addition, nonoscillatory versions of the scheme, which employ flux limiters, are explored. The limiters are based on total variation diminishing concepts and are applied to higher-order (in space) versions of the constant grid flux form scheme. The constant grid flux form schemes are Crowley-type upstream-biased Eulerian advection schemes. They are mass conserving and possess small amplitude and phase errors. The flux limiters prevent the under- and overshooting associated with the numerical dispersion of the unlimited schemes. The limited schemes are easy to implement, efficient, and nonoscillatory. Of these schemes the third-order and fifth-order versions employing the PDM limiter ar...