Conservative temperature

Derivation

Potential enthalpy

The fundamental thermodynamic relation is given by:

$${\displaystyle dh-\frac{1}{\rho }dp=Td\sigma +\mu dS}$$

where ${\displaystyle h}$ is the specific enthalpy, ${\displaystyle p}$ is the pressure, ${\displaystyle \rho }$ is the density, ${\displaystyle T}$ is the temperature, ${\displaystyle \sigma }$ is the specific entropy, ${\displaystyle S}$ is the salinity and ${\displaystyle \mu }$ is the relative chemical potential of salt in seawater.

During a process that does not lead to the exchange of heat or salt, entropy and salinity can be assumed constant. Therefore, taking the partial derivative of this relation with respect to pressure yields:

$${\displaystyle {\left(\frac{\mathrm{\partial }h}{\mathrm{\partial }p}\right)}_{S,\sigma }=\frac{1}{\rho }}$$

By integrating this equation, the potential enthalpy ${\displaystyle h^0}$ is defined as the enthalpy at a reference pressure ${\displaystyle p_r}$:

$${\displaystyle h^0(S,\theta ,p_r)=h(S,\theta ,p)-{\int }_{p_r}^p\frac{1}{\rho (S,\theta ,p^{\prime })}d{p}^{\prime }}$$

Here the enthalpy and density are defined in terms of the three state variables: salinity, potential temperature and pressure.

Conversion to conservative temperature

Conservative temperature ${\displaystyle \mathrm{\Theta }}$ is defined to be directly proportional to potential enthalpy. It is rescaled to have the same units (Kelvin) as the in-situ temperature:

$${\displaystyle \mathrm{\Theta }=\frac{h^0}{C_p^0}}$$

where ${\displaystyle C_p^0}$ = 3989.24495292815 J kg⁻¹K⁻¹ is a reference value of the specific heat capacity, chosen to be as close as possible to the spatial average of the heat capacity over the entire ocean surface.

Conservative properties of potential enthalpy

Conservation form

The first law of thermodynamics can be written in the form:

$${\displaystyle \rho \left(\frac{Dϵ}{Dt}-(p_0+p)\frac{1}{{\rho }^2}\frac{D\rho }{Dt}\right)=-\mathrm{\nabla }\cdot {\mathbf{F}}_{\mathbf{Q}}+\rho {ϵ}_M}$$

or equivalently:

$${\displaystyle \rho \left(\frac{Dh}{Dt}-\frac{1}{\rho }\frac{Dp}{Dt}\right)=-\mathrm{\nabla }\cdot {\mathbf{F}}_{\mathbf{Q}}+\rho {ϵ}_M}$$

where ${\displaystyle ϵ}$ denotes the internal energy, ${\displaystyle {\mathbf{F}}_{\mathbf{Q}}}$ represents the flux of heat and ${\displaystyle \rho {ϵ}_M}$ is the rate of dissipation, which is small compared to the other terms and can therefore be neglected. The operator ${\displaystyle \frac{D}{Dt}=\frac{\mathrm{\partial }}{\mathrm{\partial }t}+\mathbf{u}\cdot \mathrm{\nabla }}$ is the material derivative with respect to the fluid flow ${\displaystyle \mathbf{u}}$, and ${\displaystyle \mathrm{\nabla }}$ is the nabla operator.

In order to show that potential enthalpy is conservative in the ocean, it must be shown that the first law of thermodynamics can be rewritten in conservation form. Taking the material derivative of the equation of potential enthalpy yields:

$${\displaystyle \frac{D{h}^0}{Dt}=\frac{Dh}{Dt}-\frac{1}{\rho }\frac{Dp}{Dt}-\frac{D\theta }{Dt}{\int }_{p_r}^p\frac{\tilde{\alpha }(S,\theta ,p^{\prime })}{\rho (S,\theta ,p^{\prime })}d{p}^{\prime }+\frac{DS}{Dt}{\int }_{p_r}^p\frac{\tilde{\beta }(S,\theta ,p^{\prime })}{\rho (S,\theta ,p^{\prime })}d{p}^{\prime }}$$

where ${\displaystyle \tilde{\alpha }=-\frac{1}{\rho }{\left(\frac{\mathrm{\partial }\rho }{\mathrm{\partial }\theta }\right)}_{S,p}}$ and ${\displaystyle \tilde{\beta }=\frac{1}{\rho }{\left(\frac{\mathrm{\partial }\rho }{\mathrm{\partial }S}\right)}_{\theta ,p}}$. It can be shown that the final two terms on the right-hand side of this equation are as small or even smaller than the dissipation rate discarded earlier and the equation can therefore be approximated as:

$${\displaystyle \frac{D{h}^0}{Dt}=\frac{Dh}{Dt}-\frac{1}{\rho }\frac{Dp}{Dt}}$$

Combining this with the first law of thermodynamics yields the equation:

$${\displaystyle \rho \frac{D{h}^0}{Dt}=-\mathrm{\nabla }\cdot {\mathbf{F}}_{\mathbf{Q}}}$$

which is in the desired conservation form.

Comparison to potential temperature

Given that conservative temperature was initially introduced to correct errors in the oceanic heat content, it is important to compare the relative errors made by assuming that conservative temperature is conserved to those originally made by assuming that potential temperature is conserved. These errors occur from non-conservation effects that are due to entirely different processes; for conservative temperature heat is lost due to work done by compression, whereas for potential temperature this is due to surface fluxes of heat and freshwater. It can be shown that these errors are approximately 120 times smaller for conservative temperature than for potential temperature, making it far more accurate as a representation of the conservation of heat in the ocean.

Usage

TEOS-10 framework

Conservative temperature is recommended under the TEOS-10 framework as the replacement for potential temperature in ocean models. Other developments in TEOS-10 include:

• Replacement of practical salinity with the absolute salinity ${\displaystyle S_A}$ as the primary salinity variable,
• Introduction of preformed salinity as a conservative variable under biogeochemical processes,
• Defining all oceanic variables with respect to the Gibbs function.

Models

Conservative temperature has been implemented in several ocean general circulation models such as those involved in the Coupled Model Intercomparison Project Phase 6 (CMIP6). However, as these models have predominantly used potential temperature in previous generations, not all models have decided to switch to conservative temperature.