2. The Semicompressible Equations
As discussed in [
1], the Type I equations are
where density,
, is evaluated using Conservative Temperature, Θ, Absolute Salinity,
and the hydrostatic pressure,
, determined from the reference density profile,
(at most a function of
z) according to
The quantity
is viscosity and
.
De Szoeke and Samelson [
3] have shown that the hydrostatic Boussinesq equations are equivalent to the compressible hydrostatic Navier-Stokes equations written in pressure coordinates. Equation (1a) improves upon the classic hydrostatic balance in the hydrostatic limit through the presence of the sound speed. The density approximation it represents is a more accurate approximation of true density than
, density evaluated at the hydrostatic pressure. Equation (1) is also more flexible in that they can be used in non-hydrostatic settings.
The notation
in Equation (1) denotes the non-advective contributions to
X. Absolute Salinity, being a purely conserved quantity (ignoring chemistry, see [
4]), the form is clear
where
is the diffusive flux of Absolute Salinity. We purposely leave the form for
unspecified, as part of the present exercise is to determine it. Other notation is standard.
The equations in Equation (1) satisfy the energy equation
with potential energy played by the quantity
The quantity
is the viscous dissipation of kinetic energy,
is the Type I buoyancy, the quantity
Rdenotes diffusive terms and Equation (1b) demonstrates the presence of compressibility in these equations.
The Type II equations are
where
is either a chosen reference density profile yielding the sound speed profile
or a profile constructed from a chosen sound speed profile,
. In either case, density,
, is a function of depth only and the hydrostatic pressure is computed from it
Compressibility appears in these equations in Equation (7b) where the time independence of
has been recognized. The Type II equations satisfy the energy constraint
where the role of potential energy is now played by
The quantity
is the Type II buoyancy. The advantages and disadvantages of these equations are discussed in [
1]. Extracting Equations (1) and (7) from the full Navier–Stokes equations is akin to substituting the Boussinesq density for the full density in most places, and retaining a reference density profile where needed for potential energy consistency.
3. Semicompressible Thermodynamics
We now consider how to augment (1) and (7) with a consistent thermodynamic behavior. The First Law of Thermodynamics is
where
T is temperature,
is specific enthalpy,
η is specific entropy, specific volume is denoted by
,
P is pressure and
μ is the relative chemical potential of salt in seawater. The quantity
ϵ is the heat of viscous energy dissipation, and
is a generalized heat flux. Enthalpy is defined by
where
e is specific internal energy. Enthalpy from Equation (12) is seen to depend naturally on the variables
η,
, and
P, and in this form plays the role of a thermodynamic potential from which all thermodynamic variables can be obtained by differentiation. For example, derivatives of enthalpy with respect to its natural variables are
T,
μ, and
ν, respectively. The derivation of Equation (12) is discussed extensively in the TEOS-10 manual ([
4]; see also [
5,
6]).
Equation (12) is exact. The semicompressible equations consider modifications to the thermodynamic variables caused mainly by the difference between static and full pressure, and we analyze Equation (12) from that perspective. In what follows, we will adopt entropy, salinity, and pressure as state variables, and will often cast pressure dependency in terms of static pressure. For example,
where
is static pressure and
p the dynamic pressure, given by
. The superscript (
) will denote quantities evaluated at the static pressure. Specific volume, density,
etc. will be written in a manner similar to Equation (14). Note, according to Equation (12), Equation (14) is equivalent to
3.1. Type I Thermodynamics
We first develop what is a consistent statement of Type I thermodynamics. We begin by considering the form of Equation (12) to
;
i.e.,
where
and similarly for
. Evaluating the derivative of the static pressure, Equation (16) becomes
While Equation (18) is literally correct to
, it will turn out that a less accurate approximation, to
, is required on the right hand side of Equation (12)
The reason for this will become clear later, here we stress only that it assures thermodynamic consistency of the Type I equations.
3.2. Type I Mechanical Energy Equation
The mechanical energy equation derived from Equation (1a) is
where
is the kinetic energy density. Introducing the gravitational potential
and using Equation (18)
Using Equation (12)
so
where
The quantity
is recognized as the potential enthalpy defined by [
6], and the last integral in Equation (23) is the so-called dynamic enthalpy of [
7]. The quantity
is a constant reference surface pressure. The middle integral in Equation (23) reduces upon inspection to
.
Thus, Equation (21) becomes
Lastly, we exploit Equation (19) to obtain
which is in conservative form. This implies the Type I equations have a well-formed energy principle, provided they are augmented by Equations (18) and (19) as their thermodynamic equations. Note that the entropy in Equation (19) can be rewritten as
which is entropy “non-conservation” in its familiar form. Hence, the Type I equations also have a consistent entropy variable.
At this point, a practical difficulty is that entropy is not a typical ocean modelling variable. Instead, McDougall’s Conservative Temperature
McDougall ([
6]) is preferable and so it is necessary to work out a Conservative Temperature equation to complement (27). From Equation (18), one finds
By definition
so
where
θ is potential temperature and
. Potential enthalpy is thus related to entropy via
From Equation (23)
so using Equation (32)
The underlined terms in Equation (34) are quite small, according to [
8]. Neglecting them leads to the statement that Conservative Temperature is very accurately portrayed as a conserved quantity; however, all of those terms must be included in order to ensure full energy conservation.
It is seen, therefore, that the Type I equations yield energy conservation, but the cost is the insertion of a large quantity (
) into the energy equation. A similar result was labelled in [
7] as a “crippling” disadvantage of the Boussinesq-Bernoulli equation. It might thus be more useful to work with a partial equation when considering energy;
i.e., one that doesn’t explicitly include Conservative Temperature, but in any case, the full potential enthalpy Equation (34) is required. We have also found the Type I Semicompressible equations have a proper entropy variable and conserve total energy, and so are fully thermodynamically consistent. These equations are very similar to the pseudo-incompressible equation set for general equations of state in [
2].
3.3. Type II Equations: Potential, Kinetic, and Thermal Energies
The Type II equations are quite similar to the Type I equations except for the appearance of the reference density profile, rather than the Boussinesq density, in the mass conservation equation and in front of the momentum acceleration. Most of the analysis carries over in a manner similar to the Type I analysis, with the one important distinction that the velocity divergence diagnosed from mass conservation has no diabatic contributions
Therefore, the mechanical energy equation for the system analogous to Equation (21) is
with the underlined terms cancelling. Taking the material derivative of Equation (11) and substituting in Equation (36) yields
To put Equation (37) in conservative form requires the potential enthalpy equation
which is equivalent to Equation (33). This leads to
The Type II Conservative Temperature equation becomes
The most significant difference relative to the Type I set is in Equation (35), where the diabatic contributions to expansion are absent. Like for Type I, it is still necessary to retain all terms for Conservative Temperature, rather than the approximate form recommended in [
6]. We conclude that the Type II equations also have a consistent mechanical and thermal energy structure. Our Type II equations are very similar to the thermodynamically consistent anelastic set derived in [
9].
3.4. A Comparison to the Boussinesq Set
To ensure full energetic consistency for the semicompressible equations, it has been necessary to carry along some additional non-conservative effects in the calculations. We now analyze the full energetics behavior of the Boussinesq set to gauge any additional computational burden of the semicompressible set relative to them. If we adopt Conservative Temperature, salinity, and pressure as the thermodynamic coordinates and follow the procedure outlined in [
7], we arrive at the Boussinesq version of the Bernoulli equation
(see Equation (48) in [
7]) where
is the Boussinesq buoyancy and
is the Boussinesq reference density. Dynamic enthalpy as defined in [
7] is
Returning to a form of the potential enthalpy consistent with the Boussinesq approximation
where
is potential enthalphy, and taking its material derivative yields
Using Equation (44) to eliminate the pressure integrals in Equation (41) turns it into the conservative Boussinesq energy equation
The Boussinesq Conservative Temperature equation is given by
The point of this section is that to obtain a conservative energy equation out of the Boussinesq set, it is necessary to carry along the same additional diffusive terms as for the semi-compressible equations. A similar result has recently been derived in [
10,
11]. This is not done in any Boussinesq ocean circulation model of which we are aware.
4. Summary
The Type I and II semicompressible equations discussed in [
1], which are analogues of the pseudo-incompressible and anelastic equations known in meteorology, have undergone an examination of their thermodynamic structure. To summarize, the Types I and II equations are fully consistent from an energetic and thermodynamic perspective. To obtain this, it is necessary to include a more complete equation for Conservative Temperature as a part of the equation sets than is normally done. Interestingly, it turns out the Boussinesq equations must also carry analogous quantities to be in purely conservative energetic form. We thus see the semicompressible equations remain computationally competitive with the Boussinesq equations while admitting some degree of compressibility. Other results here include the associated definitions of the dynamic enthalpy variables belonging to the sets.
From a simple computer throughput perspective, this analysis recommends Type II relative to Type I because they are slightly less computationally intensive. However, the Type II equations represent the oceanic density structure by a single vertical profile. Any such profile cannot perform well globally, so the “cost” of the simplicity is a less accurate density representation. This promises to be the most troubling when computing sea level changes via the continuity equation. The additional accuracy provided by the Boussinesq density in the Type I set might well compensate for the (slight) additional complexity.
Any further exploration with these equation sets should bear in mind the above points.