Conservative temperature explained

Conservative temperature

(\Theta) is a thermodynamic property of seawater. It is derived from the potential enthalpy and is recommended under the TEOS-10 standard (Thermodynamic Equation of Seawater - 2010) as a replacement for potential temperature as it more accurately represents the heat content in the ocean.^{[1] [2]}

Motivation

Conservative temperature was initially proposed by Trevor McDougall in 2003. The motivation was to find an oceanic variable representing the heat content that is conserved during both pressure changes and turbulent mixing. In-situ temperature

T is not sufficient for this purpose, as the compression of a water parcel with depth causes an increase of the temperature despite the absence of any external heating. Potential temperature \theta can be used to combat this issue, as it is referenced to a specific pressure and so ignores these compressive effects. In fact, potential temperature is a conservative variable in the atmosphere for air parcels in dry adiabatic conditions, and has been used in ocean models for many years.^[3] However, turbulent mixing processes in the ocean destroy potential temperature, sometimes leading to large errors when it is assumed to be conservative.^[4]

By contrast, the enthalpy of the parcel is conserved during turbulent mixing. However, it suffers from a similar problem to the in-situ temperature in that it also has a strong pressure dependence. Instead, potential enthalpy is proposed to remove this pressure dependence. Conservative temperature is then proportional to the potential enthalpy.^[5]

Derivation

Potential enthalpy

The fundamental thermodynamic relation is given by:^[6]

$$dh\; -\; \backslash fracdp\; =\; T\; \backslash ,\; d\backslash sigma\; +\; \backslash mu\; \backslash ,\; dS$$

where

h is the specific enthalpy, p is the pressure, \rho is the density, T is the temperature, \sigma is the specific entropy, S is the salinity and \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:

$$\backslash left(\backslash right)\_\; =\; \backslash frac$$

By integrating this equation, the potential enthalpy

h ^{0} is defined as the enthalpy at a reference pressure p _{r} :

$$h^0(S,\; \backslash ,\; \backslash theta,\; \backslash ,\; p\_r)\; =\; h(S,\; \backslash ,\; \backslash theta,\; \backslash ,\; p)\; -\; \backslash int^p\_\; \backslash frac\; dp\text{'}$$

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

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

$$\backslash Theta\; =\; \backslash frac$$

where

C ^{0} _{p} = 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.^[7]

Conservative properties of potential enthalpy

Conservation form

The first law of thermodynamics can be written in the form:^[8]

$$\backslash rho\; \backslash left(-\; (p\_0\; +\; p)\backslash frac\; \backslash right)\; =\; -\; \backslash nabla\; \backslash cdot\; \backslash mathbf\; +\; \backslash rho\; \backslash epsilon\_M$$

or equivalently:

$$\backslash rho\; \backslash left(-\; \backslash frac\; \backslash right)\; =\; -\; \backslash nabla\; \backslash cdot\; \backslash mathbf\; +\; \backslash rho\; \backslash epsilon\_M$$

where

\epsilon denotes the internal energy, F _{Q} represents the flux of heat and \rho \epsilon _{M} is the rate of dissipation, which is small compared to the other terms and can therefore be neglected. The operator {D \over Dt} = {\partial \over \partial t} + u ⋅ \nabla is the material derivative with respect to the fluid flow u , and \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:

$$=\; -\; \backslash frac\; -\; \backslash int^p\_\; \backslash frac\; dp\text{'}\; +\; \backslash int^p\_\; \backslash frac\; dp\text{'}$$

where

\tilde{\alpha} = - (1)/(\rho) \left( { \partial \rho \over \partial \theta } \right) _{S,} and \tilde{\beta} = (1)/(\rho) \left( { \partial \rho \over \partial S } \right) _{\theta,} . 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:

$$=\; -\; \backslash frac$$

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

$$\backslash rho\; =\; -\; \backslash nabla\; \backslash cdot\; \backslash mathbf$$

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

S _{A} as the primary salinity variable, ^[9]

• Introduction of preformed salinity as a conservative variable under biogeochemical processes,^[10]
• Defining all oceanic variables with respect to the Gibbs function.^[11]

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.

Notes and References

1. Book: IOC. SCOR. IAPSO. amp. The international thermodynamic equation of seawater – 2010: Calculation and use of thermodynamic properties . 2010 . Intergovernmental Oceanographic Commission, UNESCO (English) . 196pp .
2. McDougall . Trevor J. . 2003. Potential Enthalpy: A Conservative Oceanic Variable for Evaluating Heat Content and Heat Fluxes . Journal of Physical Oceanography. 33 . 5 . 945–963 . 10.1175/1520-0485(2003)033<0945:PEACOV>2.0.CO;2. 2003JPO....33..945M . free.
3. 2015-04-01. Observational and energetics constraints on the non-conservation of potential/Conservative Temperature and implications for ocean modelling. Ocean Modelling. en. 88. 26–37. 10.1016/j.ocemod.2015.02.001. 1463-5003. Tailleux . Rémi . 2015OcMod..88...26T . 26958618 .
4. Graham. Felicity S.. McDougall. Trevor J.. 2013-05-01. Quantifying the Nonconservative Production of Conservative Temperature, Potential Temperature, and Entropy. Journal of Physical Oceanography. en. 43. 5. 838–862. 10.1175/jpo-d-11-0188.1. 2013JPO....43..838G . 0022-3670. free.
5. Book: Paulowicz, Rich . April 2013 . TEOS-10 Primer: What every oceanographer needs to know . 10pp . Department of Earth, Ocean and Atmospheric Sciences, University of British Columbia . 16 February 2023 .
6. Warren. Bruce A.. August 2006 . The First Law of Thermodynamics in a salty ocean. Progress in Oceanography. en . 70 . 2–4 . 149–167 . 10.1016/j.pocean.2006.01.001. 2006PrOce..70..149W . 1912/1289 . 0079-6611. free.
7. 2003-07-01. A new extended Gibbs thermodynamic potential of seawater . Progress in Oceanography . en . 58. 1. 43–114. 10.1016/S0079-6611(03)00088-0 . 0079-6611. Feistel . Rainer .
8. Davis. Russ E.. 1994-04-01. Diapycnal Mixing in the Ocean: Equations for Large-Scale Budgets. Journal of Physical Oceanography . EN. 24. 4. 777–800. 10.1175/1520-0485(1994)024<0777:DMITOE>2.0.CO;2 . 0022-3670. free. 1994JPO....24..777D .
9. Wright. D. G.. Pawlowicz. R.. McDougall. T. J.. Feistel. R.. Marion. G. M.. 2011-01-06. Absolute Salinity, Density Salinity and the Reference-Composition Salinity Scale: present and future use in the seawater standard TEOS-10. Ocean Science. English. 7. 1. 1–26. 10.5194/os-7-1-2011. 2011OcSci...7....1W . 1812-0784. free.
10. Pawlowicz. R.. Wright. D. G.. Millero. F. J.. 2011-06-01. The effects of biogeochemical processes on oceanic conductivity/salinity/density relationships and the characterization of real seawater. Ocean Science. English. 7. 3. 363–387. 10.5194/os-7-363-2011. 2011OcSci...7..363P . 1812-0784. free.
11. 2008-12-01. A Gibbs function for seawater thermodynamics for −6 to 80 °C and salinity up to 120 g kg–1. Deep Sea Research Part I: Oceanographic Research Papers. en. 55. 12. 1639–1671. 10.1016/j.dsr.2008.07.004. 0967-0637. Feistel . Rainer . 2008DSRI...55.1639F .