Ocean mode

PALM allows for studying the ocean mixed layer (OML) by using the ocean mode where the sea surface is defined at the top of the model, so that negative values of z indicate the depth. Hereafter, we keep the terminology and use the word surface and index 0 for variables at the sea surface and top of the ocean model. For a list of ocean specific parameters, see the following Table 3.

List of ocean model parameters

The ocean mode differs from the atmospheric mode by a few modifications, which are handled in the code by case distinctions, so that both versions share the same basic code. In particular, seawater buoyancy and static stability depend not only on θ, but also on the salinity Sa. In order to account for the effect of salinity on density, a prognostic equation is added for Sa (in PSU, practical salinity unit):

\begin{align*}
 \frac{\partial {\text{Sa}}}{\partial t} = -  \frac{\partial u_j
    {\text{Sa}}}{\partial x_j} - \frac{\partial}{\partial
    x_j}\left(\overline{u_j^{\prime\prime}{\text{Sa}}^{\prime\prime}}\right) +
  \Psi_{\text{Sa}}\;,
\end{align*}

where ΨSa represents sources and sinks of salinity. Furthermore, θv is replaced by potential density ρθ in the buoyancy term of the momentum equation (see Sect. governing equations)

\begin{align*}
 + g \frac{\theta_\mathrm{v} -
    \langle\theta_{\mathrm{v}}\rangle}{\langle\theta_{\mathrm{v}}\rangle}
  \delta_{i3}\,\rightarrow\,- g \frac{\rho_{\theta} -
    \langle\rho_{\theta}\rangle}{\langle\rho_{\theta}\rangle}
  \delta_{i3}\;,
\end{align*}

in the stability-related term of the SGS-TKE equation (see Sect. turbulence closure)

\begin{align*}
+ \frac{g}{\theta_{\mathrm{v},0}}\overline{u_3^{\prime\prime}
    {\theta_{\mathrm{v}}}^{\prime\prime}}\,\rightarrow\,+
  \frac{g}{\rho_{\theta,0}}\overline{u_3^{\prime\prime}
    {\rho_{\theta}}^{\prime\prime}}
\end{align*}

as well as in the calculation of the mixing length (see Sect. turbulence closure)

\begin{align*}
  & \left(\frac{g}{\theta_{\mathrm{v},0}}
    \frac{\partial{\theta_{\mathrm{v}}}}{\partial
      z}\right)^{-\frac{1}{2}}\,\rightarrow\,\left(\frac{g}{\rho_{\theta,0}}
    \frac{\partial{\rho_{\theta}}}{\partial
      z}\right)^{-\frac{1}{2}}\;.
\end{align*}

ρθ is calculated from the equation of state of seawater after each time step using the algorithm proposed by Jackett et al. (2006). The algorithm is based on polynomials depending on Sa, θ, and p (see Jackett et al., 2006, Table A2). At the moment, only the initial values of p enter this equation.
The ocean is driven by prescribed fluxes of momentum, heat and salinity at the top. The boundary conditions at the bottom of the model can be chosen as for atmospheric runs, including the possibility to use topography at the sea bottom.
Starting from revision r...., the ocean mode can also account for the effect of surface waves (Langmuir circulation and wave-breaking) as described in Noh et al. (2004). The ocean mode in its current state was recently used for simulations of the ocean mixed layer by Esau (2014), who investigated indirect air-sea interactions by means of the atmosphere-ocean coupling scheme that are described in Sect. coupled atmosphere-ocean simulations. Note that most previous PALM studies of the OML used the atmospheric mode, subsequent inversion of the z-axis and appropriate normalization of the results, instead of using the relatively new ocean mode (e.g., Noh et al., 2004, 2009).

The ocean mode is switched on by providing the FORTRAN namelist &ocean_parameters in the input parameter file.

References

  • Jackett DR, Mc Dougall TJ, Feistel R, Wright DG, Griffies SM. 2006. Algortihms for density, potential temperature, conservative temperature, and the freezing temperature of seawater. J. Atmos. Ocean. Tech. 23: 1709–1728.
  • Noh Y, Min HS, Raasch S. 2004. Large eddy simulation of the ocean mixed layer: the effects of wave breaking and Langmuir circulation. J. Phys. Oceanogr. 34: 720–735.
  • Esau, I. 2014. Indirect air-sea interactions simulated with a coupled turbulence-resolving model. Ocean Dynam. 64: 689–705. doi
  • Noh Y, Goh G, Raasch S, Gryschka M. 2009. Formation of a diurnal thermocline in the ocean mixed layer simulated by LES. J. Phys. Oceanogr. 39: 1244–1257.
Last modified 5 years ago Last modified on Oct 2, 2018 8:53:56 AM

Attachments (1)

Download all attachments as: .zip