Research article 15 Dec 2020
Research article  15 Dec 2020
Reappraising the appropriate calculation of a common meteorological quantity: potential temperature
 ^{1}Zentrum für Datenverarbeitung, Johannes Gutenberg University Mainz, Mainz, Germany
 ^{2}Institute for Atmospheric Physics, Johannes Gutenberg University Mainz, Mainz, Germany
 ^{3}Applied Chemicals and Materials Division, National Institute of Standards and Technology, Boulder, CO, USA
 ^{4}Forschungszentrum Jülich GmbH, Institute of Energy and Climate Research (IEK7), Jülich, Germany
 ^{5}Institute for Atmospheric and Environmental Research, University of Wuppertal, Wuppertal, Germany
 ^{6}Institut für Atmosphäre und Umwelt, GoetheUniversität Frankfurt, Frankfurt am Main, Germany
 ^{1}Zentrum für Datenverarbeitung, Johannes Gutenberg University Mainz, Mainz, Germany
 ^{2}Institute for Atmospheric Physics, Johannes Gutenberg University Mainz, Mainz, Germany
 ^{3}Applied Chemicals and Materials Division, National Institute of Standards and Technology, Boulder, CO, USA
 ^{4}Forschungszentrum Jülich GmbH, Institute of Energy and Climate Research (IEK7), Jülich, Germany
 ^{5}Institute for Atmospheric and Environmental Research, University of Wuppertal, Wuppertal, Germany
 ^{6}Institut für Atmosphäre und Umwelt, GoetheUniversität Frankfurt, Frankfurt am Main, Germany
Correspondence: Manuel Baumgartner (manuel.baumgartner@unimainz.de)
Hide author detailsCorrespondence: Manuel Baumgartner (manuel.baumgartner@unimainz.de)
The potential temperature is a widely used quantity in atmospheric science since it is conserved for dry air's adiabatic changes of state. Its definition involves the specific heat capacity of dry air, which is traditionally assumed as constant. However, the literature provides different values of this allegedly constant parameter, which are reviewed and discussed in this study. Furthermore, we derive the potential temperature for a temperaturedependent parameterisation of the specific heat capacity of dry air, thus providing a new reference potential temperature with a more rigorous basis. This new reference shows different values and vertical gradients, in particular in the stratosphere and above, compared to the potential temperature that assumes constant heat capacity. The application of the new reference potential temperature is discussed for computations of the Brunt–Väisälä frequency, Ertel's potential vorticity, diabatic heating rates, and for the vertical sorting of observational data.
According to the book Thermodynamics of the Atmosphere by Alfred Wegener (1911), the first published use of the expression potential temperature in meteorology is credited to Wladimir Köppen (1888)^{1} and Wilhelm von Bezold (1888), both following the conclusions of Hermann von Helmholtz (1888) (see also Kutzbach, 2016). Even prior to the introduction of the entropy, Poisson (1833) and Thomson (1862) used the “adiabatic equation”, the basis of what is understood today as “potential temperature”^{2}, to describe adiabatic processes, e.g. the coincident variation of temperature and pressure on the movement of air, which is “independent of the effects produced by the radiation or conduction of heat” (Thomson, 1862)^{3}. Approximately 26 years later, von Helmholtz perceived that within the atmosphere the heat exchange between air masses of different temperatures, which are relatively moved, is insufficiently explained by heat transfer due only to radiation and convection. He argued that wind phenomena (e.g. the trade winds), storm events, and the atmospheric circulation were more intense, of larger extent, and more persistent than observed if the air's heat exchange within the discontinuity region (the friction surface of the different air masses) was not mainly due to eddydriven mixing. On his way to analytically describe the heat exchange of different air masses within the atmosphere, in May 1880, von Helmholtz introduced the air's immanent heat while its absolute temperature changes with changing pressure (von Helmholtz, 1888). In essence, von Helmholtz concluded that the temperature gained by a volume of dry air due to its adiabatic descent from a certain initial pressure level (p) to ground pressure (p_{0}) corresponds to the air's immanent heat. In November of the same year, in agreement with von Helmholtz and probably inspired by a presentation that was given in June by Köppen (1888), this property was renamed and reintroduced as the air's potential temperature (θ in the following) by von Bezold (1888) with the following definition for strictly adiabatic changes of state:
where T and p are the absolute temperature and pressure, respectively, of an air parcel at a certain initial (pressure) altitude level. The quantities θ and p_{0} are corresponding values of the same air parcel's absolute temperature and pressure if the air was exposed to conditions at ground level. The dimensionless coefficient γ, nowadays called the isentropic exponent, was specified as 1.41 (von Bezold, 1888).
Moreover, in the same publication, von Bezold concluded that for moist air's adiabatic changes of state, its potential temperature remains unchanged as long as the change of state occurs within dryadiabatic limits; and further, if there is condensation and precipitation, the potential temperature changes by a magnitude that is determined by the amount of water that falls out of the air parcel. From a modern perspective, it is clear that the air parcel is an isolated thermodynamic system, and adiabatic processes correspond to processes with conserved entropy (i.e. isentropic processes). The description of the immanent heat is then equivalent to the thermodynamic state function entropy, which corresponds to potential temperature of dry air in a onetoone relationship.
In general, the potential temperature has the benefit of providing a practicable vertical coordinate (equivalent to the pressure level or the altitude above, e.g. sea level) to visualise and analyse the vertical distribution and variability of (measured) data related to any type of atmospheric parameter. Admittedly, the use of the potential temperature as a vertical coordinate is initially less intuitive than applying altitude or pressure coordinates. Indeed, the potential temperature bears a certain abstractness to describe an air parcel's state at a certain altitude level by its imaginary dryadiabatic descent to ground conditions. However, one major advantage of using the potential temperature as a vertical coordinate is that the (measured) data are sortable with respect to the entropy state at which the atmospheric samples were taken. Thus, comparing repeated measurements of an atmospheric parameter on an isentropic surface or layer excludes any diabatic change in the probed air mass.
Apart from characterising the isentropes, the vertical profiles of the potential temperature (θ as a function of height z) are used as the reference for evaluating the atmosphere's actual vertical temperature gradient, which allows characterising its static stability. Notably, von Bezold (1888) already proposed the potential temperature as an atmospheric stability criterion. In its basic formulation, the potential temperature exclusively refers to the state of dry air, and thus the potential temperature characterises the atmosphere's static stability with respect to vertical displacements of a dry air parcel. In meteorology, the static stability parameter is expressed in terms of the (squared) Brunt–Väisälä frequency N, often written in the form
where $g=\mathrm{9.81}\phantom{\rule{0.125em}{0ex}}\mathrm{m}\phantom{\rule{0.125em}{0ex}}{\mathrm{s}}^{\mathrm{2}}$ is the gravitational acceleration. The potential temperature twice enters the formulation of the stability parameter, as the denominator (θ^{−1}) and as the vertical gradient $\frac{\partial \mathit{\theta}}{\partial z}$. In the research field of dynamical meteorology, the potential vorticity (PV) is often used (Ertel, 1942; Hoskins et al., 1985; Schubert et al., 2004). The PV is proportional to the scalar product of the atmosphere's vorticity (the air's local spinning motion) and its stratification (the air's tendency to spread in layers of diminished exchange). More concretely, the PV is the scalar product of the absolute vorticity vector and the threedimensional gradient of θ, i.e. not only the potential temperature's vertical gradient but also its partial derivatives on the horizontal plane add to the resulting PV, although, particularly at stratospheric altitudes, the vertical gradient constitutes the dominant contribution. For the analytical description of a fluid's motion within a rotational system, as is the atmosphere, the PV provides a quantity that varies exclusively due to diabatic processes. Frequently, the PV is used to define the tropopause height (usually at 2 PV units; see e.g. Gettelman et al., 2011) or the edge of a largescale cyclone such as the polar winter vortex on specific θ levels (cf. Curtius et al., 2005).
While for a dry atmosphere (i.e. with little or no water vapour) the potential temperature is the correct conserved quantity (corresponding to entropy) for reversible processes, for an atmosphere containing water in two or more phases (vapour, liquid, and/or solid phases) energy transfers due to phase changes play a major role. Thus, the formulation of the potential temperature has to be extended since entropy is still the right quantity for reversible processes, including phase changes. Starting from the equation for the moist specific entropy, as derived from the first law of thermodynamics and the Gibbs equation, further extensions of the dryair potential temperature have been developed (Hauf and Höller, 1987; Emanuel, 1994; Marquet, 2011; Marquet and Geleyn, 2015) to account for phase changes and deviations from thermodynamic equilibrium, e.g. by irreversible processes. By assuming only reversible processes (i.e. conserved entropy), approximate formulas can be derived (e.g. Emanuel, 1994). However, in the case of large hydrometeors, liquid or solid particles are removed due to gravitational acceleration, leading to an irreversible process; hence the formulas based on the assumption of a reversible process are no longer applicable. Sometimes for this situation a socalled pseudoadiabatic potential temperature is defined, assuming instantaneous removal of hydrometeors from the air parcel; usually, meaningful approximations to this quantity are given, since generally it cannot be derived from first principles. Equivalent potential temperature including phase changes for vapour and liquid water is often used for the determination of convective instabilities. The general formulation can be easily adapted for an ice equivalent potential temperature, i.e. for reversible processes in pure ice clouds (see e.g. Spichtinger, 2014). Although the latent heat of sublimation is larger than the latent heat of vaporisation, the absolute mass content of water vapour decreases exponentially with decreasing temperature, leading to only small corrections due to phase changes in pure ice clouds.
At altitudes above the clouds' top, within the upper troposphere and across the tropopause, the air is substantially dried out compared to tropospheric incloud conditions. Therefore, above clouds and further aloft, e.g. within the stratosphere, the conventional dryair potential temperature may suffice to provide a meaningful vertical coordinate. Moreover, the potential temperature or the virtual potential temperature, which includes water vapour, are commonly used as prognostic variables in numerical models for the formulations of the energy equation (e.g. Skamarock et al., 2005; Skamarock and Klemp, 2008; Zängl et al., 2015; Borchert et al., 2019). Thereby, very often both variants, the potential temperature as well as the equivalent potential temperature, are involved to account for dryair situations and cloud conditions.
In any case, the use of the potential temperature requires the following preconditions to be fulfilled:

θ should be based on a rigorous derivation to ensure its validity as a function of atmospheric altitude in order not to corrupt its character as a vertical coordinate that allows for appropriately comparing (measured) atmospheric parameters, and

θ should approximate to the greatest possible extent the true entropy state of a probed air mass and should preferably account for the implied dependencies on atmospheric variables, even under the assumption that air behaves as an ideal gas,
with the aim that the potential temperature behaves as a rational physical variable. Thus, still abiding by the idealgas assumption, a reassessment of the fundamental atmospheric quantity θ is suggested, which is based on the state of knowledge of air's thermodynamic properties, and this reassessed θ is comprehensively examined concerning its ability to hold also for atmospheric conditions above the troposphere.
In principle, the concept of the potential temperature is transferable to all systems of thermally stratified fluids such as a planetary gas atmosphere or an ocean, to investigate heat fluxes (advection or diffusion) or the static stability of the fluid. In astrophysics, the potential temperature is used almost identically as in atmospheric sciences to describe dynamic processes and thermodynamic properties (e.g. static stability or vorticity) in the atmosphere of planets other than the Earth. Here, the same value p_{0}=1000 hPa, as applied to the Earth's atmosphere, is frequently used as a reference pressure for the atmosphere of other planets (Catling, 2015, Table 4), whereby the formulations of the specific heat capacity require adaptations to account for the individual gas composition of the respective planetary atmosphere. In order to simulate the weather in the atmosphere of other planets, the Weather Research and Forecasting Model (WRF) was extended to planetWRF (Richardson et al., 2007), and the governing equations considered within the WRF model include a prognostic equation for the potential temperature (Skamarock et al., 2005; Skamarock and Klemp, 2008). However, the temperature dependency of the isobaric heat capacity c_{p} is not generally negligible, especially when taking “deep atmospheres, such as on Venus” (Catling, 2015, p. 436) into account or the temperature lapse rates on other planets (Li et al., 2018). The atmosphere of Saturn's moon Titan, the only known moon with a substantial atmosphere, was comprehensively studied with frequent application of the potential temperature based on profile measurement of temperature and pressure in Titan's atmosphere by the Huygens probe (MüllerWodarg et al., 2014).
Moreover, the potential temperature is a frequently used quantity in oceanography (e.g. McDougall et al., 2003; Feistel, 2008), while here the consideration of seawater's salinity and its impact on the specific heat capacity of seawater implies additional complexity. In particular, McDougall et al. (2003) suggests a reassessment of the potential temperature as applied in oceanography to approximate the adiabatic lapse rate; thus this study bears certain parallels to the present investigation aiming at the reappraisal of the potential temperature for atmosphererelated purposes. These studies from other disciplines motivate the need for a reassessment of the potential temperature for the atmospheric sciences. Thus, the approach provided herein proposes a modified calculation of the widely used quantity of the potential temperature by additionally accounting for the current state of knowledge concerning air's properties.
The study is organised as follows. The derivation of the potential temperature for an ideal gas with constant specific heat capacity c_{p} is recalled in Sect. 2. In Sect. 3 the assumption of a constant c_{p} is discussed together with a synopsis of various c_{p} values as provided in the literature. The temperature dependency of c_{p} is examined in Sect. 4, and a parameterisation is given. Section 5 is devoted to the definition and computation of a new reference potential temperature θ_{ref} based on the temperaturedependent specific heat capacity, while Sect. 6 focuses on the influence of realgas effects on the resulting potential temperature. Section 7 presents some implications of the use of θ_{ref}, and concluding remarks are given in Sect. 8.
The Gibbs equation (see e.g. Kondepudi and Prigogine, 1998) is a general thermodynamic relation to describe the state of a system with m components and reads as
where T denotes the absolute temperature in kelvin (K), S the entropy in joules per kelvin (J K^{−1}), H the enthalpy in joules (J), V the volume in cubic metres (m^{3}), μ_{k} the chemical potential of component k in joules per kilogram (J kg^{−1}), M_{k} the mass of component k in kilograms (kg), and p the static pressure in pascals (Pa). Assuming no phase conversion or chemical reaction within the system, the mass of each component does not change; hence dM_{k}=0 for each component k.
In the following, dry air is assumed to be the single component in the system. Expressing the Gibbs equation in its specific form (i.e. division by the total mass M_{a} of dry air; note that lowercase letters indicate specific variables – e.g. $h=H/{M}_{\mathrm{a}}$) leads to
Furthermore, approximating dry air as an ideal gas leads to the following simplifications.

The idealgas law
$$\begin{array}{}\text{(5)}& pV={M}_{\mathrm{a}}{R}_{\mathrm{a}}T\end{array}$$can be applied with the specific gas constant R_{a} of dry air, which is
$$\begin{array}{}\text{(6)}& \begin{array}{rl}{R}_{\mathrm{a}}& ={\displaystyle \frac{R}{{M}_{\text{mol,a}}}}\\ & ={\displaystyle \frac{\mathrm{8.31446261815324}\phantom{\rule{0.125em}{0ex}}\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{mol}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}}{\mathrm{0.0289586}\phantom{\rule{0.125em}{0ex}}\mathrm{kg}\phantom{\rule{0.125em}{0ex}}{\mathrm{mol}}^{\mathrm{1}}\pm \mathrm{0.0000002}\phantom{\rule{0.125em}{0ex}}\mathrm{kg}\phantom{\rule{0.125em}{0ex}}{\mathrm{mol}}^{\mathrm{1}}}}\\ & \in \left[\mathrm{287.11350}\phantom{\rule{0.125em}{0ex}}\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}},\phantom{\rule{0.125em}{0ex}}\mathrm{287.11748}\phantom{\rule{0.125em}{0ex}}\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}\right],\end{array}\end{array}$$with the molar gas constant R in $\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{mol}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}$ (Tiesinga et al., 2020; Newell et al., 2018) and M_{mol,a} the molar mass of dry air (Lemmon et al., 2000), composed of nitrogen N_{2}, oxygen O_{2}, and argon Ar.

The specific enthalpy is given by
$$\begin{array}{}\text{(7)}& \mathrm{d}h={c}_{p}\phantom{\rule{0.125em}{0ex}}\mathrm{d}T,\end{array}$$with c_{p} the specific heat capacity of dry air.
Based on these assumptions, the change in the specific entropy (within the fluid dry air) is given by
For isentropic changes of state, i.e. ds=0, Eq. (8) reduces to
Note that the assumption of dry air being an ideal gas does not imply that in Eq. (9) the specific heat capacity c_{p} is constant. While statistical mechanics excludes any pressure dependence in the idealgas heat capacity, the general derivation (cf. Appendix A) permits a temperature dependence of c_{p}. However, usually the temperature dependence is neglected in atmospheric physics, and, instead, c_{p} is assumed as constant (see e.g. Ambaum, 2010, p. 48/49, where vibrational modes of the air molecules are neglected). Immediately below and in Sect. 3, the treatment of c_{p} as a temperatureindependent constant is discussed. The introduction of the temperature dependence then follows in Sect. 4.
Treating c_{p} as a constant, rearrangement of Eq. (9) leads to
Integration of Eq. (10) over the range from groundlevel pressure and temperature (p_{0}, T_{0}) to the pressure and temperature at a specific height (p, T) yields
and, after another straightforward conversion, one arrives at
With the definition ${\mathit{\theta}}_{{c}_{p}}={T}_{\mathrm{0}}$, Eq. (12) is transformed into the commonly used expression for determining the potential temperature
for which the groundlevel pressure p_{0} is arbitrary but usually set to p_{0}=1000 hPa. This choice coincides with the definition of the World Meteorological Organisation (WMO, 1966) and the standardstate pressure (Tiesinga et al., 2020) but should not be confused with the standard atmosphere 101 325 Pa (Tiesinga et al., 2020). In the following, ${\mathit{\theta}}_{{c}_{p}}$ denotes the potential temperature based on a constant c_{p}, and, when a specific value of c_{p} is applied, the subscript c_{p} in the potential temperature's notation is replaced by the corresponding c_{p} value.
The general theory of thermodynamics, assuming dry air as an ideal gas, gives the expression
for the constant specific heat capacity, which is based on the results of statistical mechanics and the equipartition theorem (e.g. Huang, 1987). In Eq. (14), the parameter $f={f}_{\mathrm{trans}}+{f}_{\mathrm{rot}}+{f}_{\mathrm{vib}}$ is equal to the total number of degrees of freedom of the gas molecules of which dry air consists. The individual contributions to f comprise the degrees of freedom of translation f_{trans}, rotation f_{rot}, and vibration f_{vib}. Assuming further that dry air exclusively consists of the linear molecules N_{2} and O_{2} (implying f_{trans}=3 and f_{rot}=2, while the contribution of Ar remains disregarded) and additionally neglecting the vibrational degrees of freedom (f_{vib}=0), the general relation Eq. (14) reduces to
Although the neglect of vibrational excitation, particularly at very low temperatures, seems plausible and appropriate, errors are already introduced by this assumption for the temperature range relevant in the atmosphere.
In atmospheric sciences, for the majority of computations that require the specific heat capacity of dry air, a constant value of c_{p} may be appropriate. According to the WMO (1966), the recommended value for c_{p} of dry air is 1005 $\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}$, and, furthermore (ibid.), it is defined that $\mathit{\gamma}=\frac{{c}_{p}}{{c}_{v}}=\frac{\mathrm{7}}{\mathrm{5}}=\mathrm{1.4}$ (cf. Eq. 1). This definition is consistent with the general thermodynamic theory together with all aforementioned additional assumptions and results in Eq. (15) as well.
Even assuming a universally valid constant c_{p}, a single consistently used value of c_{p} was not found. Instead, the specified values of c_{p} vary among different textbooks and other sources. In Table 1, some of the available values of constant specific heat capacity for dry air are compiled, indicating a variability of c_{p} that ranges from 994 to 1011 $\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}$. However, the extremes in Table 1 are from old references of historical interest only; to reflect recently stated values the narrower range 1000 to 1010 $\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}$ is considered.
Wegener and Wegener (1935, converted from units other than SI)Roedel and Wagner (2011, p. 66)(WMO, 1966)Tripoli and Cotton (1981, the Appendix therein)Holton (2004, p. 491)Wallace and Hobbs (2006, p. 75)Schumann (2012)Wendisch and Brenguier (2013, p. 24)Liou (2002, Appendix F)Ambaum (2010, table “Useful Data”)Pruppacher and Klett (2010, converted from units other than SI; p. 489)Curry and Webster (1998, p. 62)WMO (1966)Bohren et al. (1998, p. 384)Houghton (2002, p. 275)Zdunkowski and Bott (2003, p. 705)Brasseur and Solomon (2005, p. 426)Seinfeld and Pandis (2006, p. 1178)Cotton et al. (2011, Table 2.1)Bolton (1980)Emanuel (1994, Appendix 2)Wendisch and Brenguier (2013, p. 69)Stamnes et al. (2017, p. 14)Chang et al. (2006)Tiwary and Williams (2019, beneath Eq. 8.8; possibly a typo, as indicated by inconsistencies on reproducing their conclusions based on this value)Brusseau et al. (2019, p. 59)(WMO, 1966)These different values of constant c_{p} scatter within a small range (below ±1.1 %) around the WMO's recommendation 1005 $\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}$, which may seem negligible if c_{p} contributes only as a linear coefficient within an equation (e.g. in the expression of a correction factor; cf. Weigel et al., 2016). However, in the formulation of the potential temperature ${\mathit{\theta}}_{{c}_{p}}$, cf. Eq. (13), the specific heat capacity c_{p} does not contribute linearly but rather as the denominator in the exponent. Thus, the variety of different c_{p} values, although scattering within a small range, impacts the resulting ${\mathit{\theta}}_{{c}_{p}}$ significantly. To illustrate this impact, a computation of ${\mathit{\theta}}_{{c}_{p}}$ by using Eq. (13) was based on the values of static pressure (p; cf. Fig. 1a) and absolute temperature (T; cf. Fig. 1b) corresponding to the US Standard Atmosphere (United States Committee on Extension to the Standard Atmosphere, 1976). From the list of the different c_{p} in Table 1, the extreme values were selected in order to initially illustrate the sensitivity of the resulting ${\mathit{\theta}}_{{c}_{p}}$ to variations in c_{p} in the range of ∼1 %, as seen in the literature. In Fig. 2a, the individual profiles of ${\mathit{\theta}}_{{c}_{p}}$ are shown for the extremes of the historic c_{p} values (Table 1), while Fig. 2b illustrates the absolute differences $\mathrm{\Delta}{\mathit{\theta}}_{{c}_{p}}={\mathit{\theta}}_{\mathrm{994}}$–θ_{1011} (red curve), $\mathrm{\Delta}{\mathit{\theta}}_{{c}_{p}}={\mathit{\theta}}_{\mathrm{1000}}$–θ_{1010} (blue curve), and $\mathrm{\Delta}{\mathit{\theta}}_{{c}_{p}}={\mathit{\theta}}_{\mathrm{1003.5}}$–θ_{1006.5} (green curve). The absolute error exhibited with the blue curve in Fig. 2b is based on the extremes of most recently referred c_{p} values in the literature (Table 1). At an altitude of 8.5 km, the difference $\mathrm{\Delta}{\mathit{\theta}}_{{c}_{p}}$ already exceeds 1 K (blue curve). The values of $\mathrm{\Delta}{\mathit{\theta}}_{{c}_{p}}$ reach approximately 1.2 K at 10 km and rise further, above 4 K, with increasing altitude up to 20 km. At 50 km, approximately where the stratopause is located, which is the chosen upper height limit for this investigation, the computed $\mathrm{\Delta}{\mathit{\theta}}_{{c}_{p}}$ reaches 43 K. The green curve corresponds to the more realistic c_{p} interval $\mathrm{1005}\phantom{\rule{0.125em}{0ex}}\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}\pm \mathrm{1.5}$ $\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}$ as recommended by the WMO; the difference reaches approximately 13 K at the stratopause.
Figure 2 illustrates the possible spread of ${\mathit{\theta}}_{{c}_{p}}$ based on a range of c_{p} values from different literature references; hence, if one uses a different value for c_{p} from the literature than that defined by WMO (1966), the difference ${\mathit{\theta}}_{\mathrm{1005}}{\mathit{\theta}}_{{c}_{p}}$ might be significant. Since the c_{p} values provided by some literature references are close to the value c_{p}=1005 $\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}$ recommended by the WMO (1966), the subsequent comparisons will be made to θ_{1005}. The ${\mathit{\theta}}_{{c}_{p}}$ values based on c_{p} values other than 1005 $\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}$ are only used to illustrate respective deviations. Although the curves in Fig. 2b depict extremes in the deviation of potential temperatures, as they are based on the extremes of c_{p} values (cf. Table 1), they nevertheless illustrate the sensitive response of ${\mathit{\theta}}_{{c}_{p}}$ to even small variations in c_{p}, on the order of 1 %. Further proof of this sensitivity from the mathematical perspective is provided in Appendix B. The impact of this sensitivity becomes important at altitudes of ∼10 km and above, thus where the use of the potential temperature becomes increasingly meaningful. Here, and in particular above the cloud tops, the smallscale and comparatively fast tropospheric dynamics (causing vertical transport and implying diabatic processes) become diminished, while further above, towards the stratosphere, an increasingly layered vertical structure of the atmosphere is taking over.
As indicated above, the reason for this sensitivity to small variations of air's specific heat capacity is that it affects the exponent of the equation for ${\mathit{\theta}}_{{c}_{p}}$. The studies of Ooyama (1990, 2001) document an interesting attempt to formulate for example the energy balance equations for the moist atmosphere, wherein entropy replaces the more common formulation using the potential temperature. This substitution avoids the use of the potential temperature, which “is merely an exponential transform of the entropy expressed in units of temperature” (Ooyama, 2001); thus, within this equation, air’s specific heat capacity is implied exclusively as a linear coefficient. Consequently, a parameterisation for the temperature dependence of the specific heat capacity (c_{p}(T); cf. Sect. 4) may be easily adopted. However, the crucial drawback of the entropybased equations is that to gain a numerical model for weather forecast purposes, the parameterisations of most of the physical processes within the atmosphere would require a reformulation.
It should be noted that not only do literature values of air's specific heat capacity c_{p} vary, but also the values of the gas constant R_{a} vary slightly due to different historical approximations for the molar gas constant^{4} R and for the composition of dry air. The variation of values for R_{a} is typically only on the order of 0.1 $\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}$, whereas the variability in c_{p} is on the order of a few $\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}$ (cf. Table 1). Therefore, within the exponent of the expression (13) for ${\mathit{\theta}}_{{c}_{p}}$, the variability of c_{p} has by far a stronger impact on the resulting ${\mathit{\theta}}_{{c}_{p}}$ value than the variability of R_{a}.
However, accepting for a moment the WMO's definition (15) of c_{p} (WMO, 1966), the variability of air's c_{p} should naturally be constrained to certain limits. With the specific gas constant R_{a}=287.05 $\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}$ (WMO, 1966), the WMO's definition leads to c_{p}=1004.675 $\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}$. In contrast, taking into account the uncertainty introduced in R_{a} by the molar mass of dry air, cf. Eq. (6), the resulting range for air's specific heat capacity is $\mathrm{1004.897}\phantom{\rule{0.125em}{0ex}}\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}\le {c}_{p}\le \mathrm{1004.912}$ $\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}$. It may be surmised that the rounded value c_{p}=1005 $\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}$ as recommended by the WMO (1966) had the main goal to simplify certain calculations, which at the time may have been mostly done by hand.
Next, while retaining the idealgas assumption, we consider the dependence of air's c_{p} on temperature, mainly over the atmospherically relevant range (180 to 300 K). The temperature dependence of c_{p} is, of course, not a new finding. Experimental approaches for determining the calorimetric properties of air and the temperature dependence of a fluid's specific heat capacity are described by Witkowski (1896), who investigated the change in the mean c_{p} as a function of temperature intervals between room temperature (as a fixed reference) and various warmer and colder temperatures, for atmospheric pressures and slightly beyond. Despite the potentially high uncertainty of the experimental results from these times, Witkowski (1896) already indicated that with decreasing temperature the experimentally determined c_{p} values initially decline, then pass a minimum, and subsequently increase again at lower temperatures (T<170 K). The description of refined experiments and ascertainable data of air's c_{p}(T) for temperatures below 293 K is summarised by Scheel and Heuse (1912), Jakob (1923), and Roebuck (1925, 1930), illustrating in comprehensive detail the experimental effort and providing the resulting data. The review by Awano (1936) compiled and compared the data of c_{p}(T) of dry air (“air containing neither carbondioxide nor steam”, Awano, 1936), and he attested – at that time – the previously mentioned studies to constitute “the most reliable experiments”. During the decades following these experiments, further insights were gained and landmarks were reached, which are summarised in the comprehensive survey by Lemmon et al. (2000) of the progress of modern formulations for the thermodynamic properties of air and about the experiments the previous formulations were based on.
Figure 3 illustrates the range of suggested constant values for the specific heat capacity as indicated in Table 1 (dashed curves) together with the measurements that were made to obtain air's behaviour as a function of temperature and pressure. Note that Fig. 3 includes data at other atmospheric pressures, indicated by squares, diamonds, and triangles. In the same figure, calculated values of c_{p}(T) of dry air are displayed, resulting from the equation of state which was derived from experimental p, V, and T data by Vasserman et al. (1966), who provided an extensive review of previous experimental and theoretical works and of the state of knowledge at that time. In addition, Fig. 3 exhibits two different parameterisations, by Lemmon et al. (2000) and by Dixon (2007, see p. 376 in his book – the accuracy is “within 0.1 % from 200 to 450 K”), which account for the temperature dependence of the specific heat capacity c_{p}(T). The parameterisation by Lemmon et al. (2000), to be discussed in detail in Sect. 4.2, is valid for dry air assumed as an ideal gas, whereas this distinction is not explicitly made in Dixon (2007). Moreover, Fig. 3 contains discrete values of dry air's c_{p}(T) extracted from the database REFPROP (Reference Fluid Thermodynamic and Transport Properties Database by NIST, the National Institute of Standards and Technology, Lemmon et al., 2018), which is based on parameterisations resulting from thermodynamic considerations discussed later.
The measurement data, as well as the parameterisations, clearly indicate a dependence of air's specific heat capacity on the temperature. At temperatures above 300 K, the data points by Jakob (1923) are surprisingly well captured by the parameterisations, while below 270 K the course of the parameterised and measured c_{p}(T) diverges significantly. Possible reasons for this include the following:

the measurements of c_{p}(T) have a precision likely no better than 1 % (in particular the historical measurements), and there could be systematic errors, especially at low temperatures;

the measured data reflect the true thermodynamic behaviour of the real gas rather than that of an ideal gas.
However, it is immediately obvious from Fig. 3 that a good agreement among (i) the experimentally determined c_{p}(T) data, (ii) a constant c_{p} (e.g. 1005 $\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}$; WMO, 1966), and (iii) the parameterised c_{p}(T) is found only for a temperature interval ranging from 270 to 300 K. For air temperatures below 270 K, the constant value c_{p}=1005 $\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}$ is only comparable with the values from Vasserman et al. (1966) but fails to coincide with other parameterised or experimentally determined values of c_{p}(T).
4.1 The temperature dependence of the idealgas specific heat capacity
As already indicated by the data depicted in Fig. 3, the specific heat capacity c_{p} depends on the gas temperature. With regard to measured values, the lack of constancy may be due to realgas effects or to a dependence of the idealgas heat capacity on temperature. In this section, we focus on the latter effect, denoting the idealgas isobaric specific heat capacity by ${c}_{p}^{\mathrm{0}}\left(T\right)$, where the superscript 0 indicates the underlying idealgas assumption. For an individual gas, there is always a contribution from the three translational degrees of freedom, ${c}_{p,\mathrm{trans}}^{\mathrm{0}}=\frac{\mathrm{5}}{\mathrm{2}}{R}_{i}$, where R_{i} is the specific gas constant of the gas. If the molecule is assumed to be a rigid rotor, there is also a rotational contribution given by
As mentioned previously, at finite temperatures molecules also have contributions to ${c}_{p}^{\mathrm{0}}\left(T\right)$ from intramolecular vibrations (and, at high temperatures, excited electronic states). To arrive at a temperaturedependent parameterisation for the idealgas specific heat capacity of dry air, the compounds' individual contributions, considering all degrees of freedom, need to be parameterised and then combined according to each compound's proportion in the mixture. For the following, dry air is considered a threecomponent mixture: the diatomic gases nitrogen (N_{2}) and oxygen (O_{2}) and the monatomic gas argon (Ar).
To determine the contribution of N_{2} to ${c}_{p}^{\mathrm{0}}\left(T\right)$, both Bücker et al. (2002) and Lemmon et al. (2000) use the idealgas heat capacity from the reference equation of state of Span et al. (2000) that compares well with the findings from other studies within an uncertainty $\mathrm{\Delta}{c}_{p}^{\mathrm{0}}$ of less than 0.02 %.
For the contribution of O_{2}, Lemmon et al. (2000) use the formulation given by Schmidt and Wagner (1985). Alternatively, Bücker et al. (2002) provide a slightly different formulation from the International Union of Pure and Applied Chemistry (IUPAC, Wagner and de Reuck, 1987), after refitting it to more recently obtained data, thereby achieving an overall uncertainty $\mathrm{\Delta}{c}_{p}^{\mathrm{0}}$ of less than ±0.015 % for O_{2} (Bücker et al., 2002). However, the difference in the resulting specific heat capacity contribution by O_{2} between the two approaches (Lemmon et al., 2000, or Bücker et al., 2002) is comparatively small. The recent work of Furtenbacher et al. (2019) leads to values of ${c}_{p}^{\mathrm{0}}$ for O_{2} with even smaller uncertainties, but the differences from the values used here are negligible in our context.
For a monatomic gas such as Ar, vibrational and rotational contributions to the heat capacity do not exist, and Bücker et al. (2002) consider that argon's excited electronic states are relevant only at temperatures above 3500 K. Hence, the contribution of Ar to the specific heat capacity of air reduces to ${c}_{p}^{\mathrm{0}}=\frac{\mathrm{5}}{\mathrm{2}}{R}_{\mathrm{Ar}}$.
The approach by Bücker et al. (2002) additionally considers the contribution of further constituents of air, such as water, carbon monoxide, carbon dioxide, and sulfur dioxide. These authors provide an analytical expression for specific heat capacity, accounting for this more complex but proportionally invariant air composition which is specified to deviate from the used reference by $\mathrm{\Delta}{c}_{p}^{\mathrm{0}}\le \pm \mathrm{0.015}\phantom{\rule{0.125em}{0ex}}\mathrm{\%}$ in the temperature range of $\mathrm{200}\phantom{\rule{0.125em}{0ex}}\mathrm{K}\le T\le \mathrm{3300}\phantom{\rule{0.125em}{0ex}}\mathrm{K}$. At atmospheric altitudes above the clouds' top, i.e. on average above ∼11 km, the air is assumed to have lost most of its water and is deemed as dry. Furthermore, for the following, trace gases that contribute to air's composition by molar fractions of less than that of Ar are neglected.
4.2 NIST's parameterisation of ${c}_{p}^{\mathrm{0}}\left(T\right)$
Besides a comprehensive survey of the available experimental data for the specific heat capacity of air, Lemmon et al. (2000) also provide stateoftheart knowledge for other thermodynamic properties (isochoric heat capacity, speed of sound, vapour–liquid equilibrium, etc.). Additionally, they give two approaches to derive air's thermodynamic properties, including the vapour–liquid equilibrium:

an empirical modelbased equation of state for standard (dry) air considered as a pseudopure fluid and

assembly of a mixture model from equations of state for each pure fluid.
Each approach allows calculating the thermodynamic properties, e.g. c_{p}, of gas mixtures such as dry air, and both are realgas models with the idealgas behaviour as a boundary condition. The major difference between the models is that the first approach considers air as a pseudopure fluid while the second, more rigorous approach treats air as a mixture composed of N_{2}, O_{2}, and Ar, in molar fractions of 0.7812, 0.2096, and 0.0092, respectively, following Lemmon et al. (2000, their Table 3). This fractional composition of dry air is assumed to be constant from ground level up to 80 km height (United States Committee on Extension to the Standard Atmosphere, 1976), and its fractional composition would have to be shifted significantly to cause a serious deviation of the resulting potential temperature. The contribution to the composition by carbon dioxide (CO_{2}) and of any other trace species was assumed to be negligible. The validity of both approaches is specified for various states of dry air, from its solidification point (59.75 K) up to temperatures of 1000 K, and for pressures up to 100 MPa and even much further beyond the pressure range that is relevant for atmospheric investigations. Both the pseudopure fluid model and the mixture model are implemented in NIST's REFPROP database (cf. https://www.nist.gov/srd/refprop, last access: 2 November 2020) for various physical properties of fluids over a wide range of temperatures and pressures.
Both the pseudopure fluid model and the mixture model of Lemmon et al. (2000) use the same expression for the idealgas heat capacity, which is rigorously given as a sum of the purecomponent contributions:
where x_{i} denotes the molar fraction of species i, and ${C}_{p}^{\mathrm{0}}$ and the molar gas constant R are given in units of $\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{mol}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}$.
Like Bücker et al. (2002), Lemmon et al. (2000) use the expression of Span et al. (2000) for the contribution of N_{2} to the heat capacity and adopt ${C}_{p}^{\mathrm{0}}=\frac{\mathrm{5}}{\mathrm{2}}R$ for Ar. Together with the contribution by O_{2} according to the formulation by Schmidt and Wagner (1985), the expression provided by Lemmon et al. (2000, Eq. 18 therein) for the idealgas heat capacity of dry air is
with the scalar coefficients N_{i} for dry air (ibid.),
which is specified as valid for temperatures from 60 to 2000 K. Because the underlying calculations are based on rigorous statistical mechanics and accurate spectroscopic data, $\frac{{C}_{p}^{\mathrm{0}}\left(T\right)}{R}$ should be accurate to within 0.01 % throughout this range, as discussed by Span et al. (2000).
The parameterisation (18) provides the isobaric specific heat capacity of dry air, considered as a mixture of ideal gases. This represents a more rigorous and accurate behaviour than assuming it to be a constant.
4.3 The parameterisation of ${c}_{p}^{\mathrm{0}}\left(T\right)$ from an engineer's perspective
The parameterisation from Dixon (2007)
for $\mathrm{200}\phantom{\rule{0.125em}{0ex}}\mathrm{K}\le T\le \mathrm{450}\phantom{\rule{0.125em}{0ex}}\mathrm{K}$ is not explicitly described to be based on particular assumptions or data sets. The author indicates his suggested parameterisation to hold within 0.1 % for temperatures between 200 and 450 K. For elevated air temperatures, the deviation between the idealgas limit ${c}_{p}^{\mathrm{0}}\left(T\right)$ (Lemmon et al., 2000) and Dixon's parameterisation substantially increases. This is most likely due to the chosen type of polynomial approximation (Dixon, 2007), which increasingly departs from the reference ${c}_{p}^{\mathrm{0}}\left(T\right)$ for gas temperatures exceeding 450 K.
Concerning the thermophysical properties of humid air, the study by Tsilingiris (2008) provides further insight. Its purpose was to evaluate the transport properties as a function of different levels of the relative humidity and as a function of temperature (from 273 to 373 K) for the gas mixture of air with water vapour at a constant pressure (1013 hPa). The atmospherically relevant pressure range below 1013 hPa and temperatures smaller than 273 K were not considered. Although this study focused on providing a comprehensive account of moisture within air, mainly for technical purposes and engineering calculations, the possible usefulness of these findings to atmospheric investigations is also apparent. However, the impact of water vapour on the resulting gas mixture's c_{p}(T) is significantly larger (cf. Tsilingiris, 2008) than the uncertainty of dry air's c_{p}(T) that is discussed in the present work. Furthermore, the consideration of water vapour as a component of air requires very individual and casespecific computations of c_{p}(T) of moist air, as water vapour is among the most variable constituents of the atmosphere.
The effort required to produce an analytical formulation for gas properties which best reflects the true gas behaviour may indicate that for engineering purposes (pneumatic shock absorbers, engines' combustion efficiency, improvements of turbofan/turboprop propulsion, aerodynamics, material sciences, etc.), especially where pressures exceed atmospheric, the assumption of idealgas behaviour introduces excessive uncertainty.
Previously introduced approaches for computing the specific heat capacity of dry air call for a brief discussion on how to use the obtained c_{p}(T) to derive the potential temperature. In the following, ${\mathit{\theta}}_{{c}_{p}\left(T\right)}$ denotes the derived potential temperature that accounts for the temperature dependence of dry air's specific heat capacity. Furthermore, it should be noted that simply substituting any c_{p}(T) value into the conventionally used equation and defining Eq. (13) for ${\mathit{\theta}}_{{c}_{p}}$ (WMO, 1966) may appear tempting but definitely leads to results inconsistent with ${\mathit{\theta}}_{{c}_{p}\left(T\right)}$, which is based on the reference parameterisation of dry air's c_{p}(T). Therefore, the thermodynamically consistent use of c_{p}(T) in the derivation of θ is described in the following.
5.1 Derivation of ${\mathit{\theta}}_{{c}_{p}\left(T\right)}$ based on the temperaturedependent specific heat capacity of dry air
In the derivation of the potential temperature (cf. Sect. 2), we note that, until reaching the expression for isentropic changes of state (9), no assumption was made about the specific heat capacity. As soon as the temperature dependence of the specific heat capacity comes into play, the reassessment of Eq. (9) leads to
Integration of Eq. (21) from the basic state $({p}_{\mathrm{0}},\phantom{\rule{0.125em}{0ex}}{\mathit{\theta}}_{{c}_{p}\left(T\right)})$ to any other state (p, T) yields
where ${\mathit{\theta}}_{{c}_{p}\left(T\right)}$ is the desired potential temperature.
The rearrangement of Eq. (22) makes evident that the desired potential temperature is a zero of the function F(x), given by
To arrive at the desired potential temperature ${\mathit{\theta}}_{{c}_{p}\left(T\right)}$ for any given temperature and pressure, the equation 0=F(x) must be solved for the variable x, which is the desired ${\mathit{\theta}}_{{c}_{p}\left(T\right)}$. Equation (23) has at most only one real zero, since its integrand is strictly positive, which means F(x) is strictly monotonic.
In the following, the idealgas reference potential temperature θ_{ref} is introduced, based on the formulation of the idealgas limit of dry air's specific heat capacity ${c}_{p}^{\mathrm{0}}\left(T\right)$ in accordance with Eq. (18) as formulated by Lemmon et al. (2000). This reference potential temperature θ_{ref} represents the zero of F(x) in Eq. (23), wherein c_{p}(T^{′}) is to be replaced by ${c}_{p}^{\mathrm{0}}\left({T}^{\prime}\right)$; i.e. for given p, T the reference potential temperature θ_{ref} solves the equation
The parameterisation of ${c}_{p}^{\mathrm{0}}\left({T}^{\prime}\right)$ is stated to give accurate values for temperatures from 60 to 2000 K (cf. Sect. 4.2); thus values of θ_{ref} should not exceed 2000 K, since otherwise ${c}_{p}^{\mathrm{0}}\left({T}^{\prime}\right)$ within the integrand in Eq. (24) is evaluated outside of its range of validity. However, due to the division by T^{′}, the value of the integrand $\frac{{c}_{p}^{\mathrm{0}}\left({T}^{\prime}\right)}{{T}^{\prime}}$ may be expected to give nevertheless a good approximation even if the accuracy of ${c}_{p}^{\mathrm{0}}\left({T}^{\prime}\right)$ is decreased; hence values θ_{ref}>2000 K should not be discarded.
It may be noted that further variants of a reference potential temperature are derivable by replacing c_{p}(T^{′}) in Eq. (23) by any other expression of the specific heat capacity of air which may appear sufficiently accurate. The steps to compute or approximate the zero of the function (23), described in this study, are independent of the chosen heat capacity formulation.
Unfortunately, for a straightforward solution of the integral (23), the suggested parameterisation of c_{p} is too complex and an analytically insolvable nonlinear equation 0=F(x) could result. Thus, an approximation of the equation's desired zero is required. Newton's method (cf. e.g. Deuflhard, 2011) provides a standard approach to numerically approximate the zero of a nonlinear equation. Proceeding from an initial guess x_{0}, Newton's method constructs a sequence {x_{k}}_{k∈ℕ} defined by the recursion
The constructed sequence {x_{k}}_{k∈ℕ} converges to the equation's desired zero. For the computations described here, the iteration is stopped as soon as the absolute difference $\left{x}_{k+\mathrm{1}}{x}_{k}\right$ of two consecutive iterations falls below 10^{−8} K.
For the reference of air's specific heat capacity, ${c}_{p}^{\mathrm{0}}\left(T\right)$, the integral (23) turns out not to be explicitly solvable. Therefore, with each iteration, the solution of the integral $\underset{{x}_{k}}{\overset{T}{\int}}\frac{{c}_{p}^{\mathrm{0}}\left({T}^{\prime}\right)}{{T}^{\prime}}\phantom{\rule{0.125em}{0ex}}\mathrm{d}{T}^{\prime}$ is approximated by subdividing the entire integration range, [x_{k}, T], into intermediate intervals with respective size of at most 0.1 K and by applying Simpson's rule on each subinterval.
As a first guess x_{0} for the Newton iteration, the conventional definition of ${\mathit{\theta}}_{{c}_{p}}$ based on a constant specific heat capacity (WMO, 1966) is inserted:
In the course of Newton's method, the sequence {x_{k}}_{k∈ℕ} will converge to the unique zero for any initial guess x_{0} due to the monotonicity of F(x). However, the right choice of the initial guess x_{0} substantially decreases the error of the first iteration x_{1}, speeding up convergence to the desired zero of the function F(x). Therefore, it seems wise to use the conventional definition of ${\mathit{\theta}}_{{c}_{p}}$ as the first guess for the Newton iteration (25).
Solving the previously described rootfinding problem by Newton's method over the comprehensive range of iteration steps (until the set requirement, i.e. $\left{x}_{k+\mathrm{1}}{x}_{k}\right<{\mathrm{10}}^{\mathrm{8}}$ K, is fulfilled) finally leads to the reference potential temperature θ_{ref}. This θ_{ref} is based on the idealgas limit of dry air's specific heat capacity ${c}_{p}^{\mathrm{0}}\left(T\right)$, which refers to the current thermodynamic state of knowledge, and, thus, we use θ_{ref} as our reference for the potential temperature in the following. For evaluating the results, the air temperature and pressure from the US Standard Atmosphere are used once more to set up the vertical profiles of the potential temperature. Figure 4a exhibits the resulting reference profile, i.e. θ_{ref} (red curve). Additionally, for comparison with the reference, further potential temperature profiles ${\mathit{\theta}}_{{c}_{p}}$ are shown based on the two (historical) extremes c_{p}=994 $\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}$ and c_{p}=1011 $\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}$ (dashed curves) and based on the range limits of more recent values c_{p}=1000 $\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}$ and c_{p}=1010 $\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}$ (solid green and magenta curves) of given constant values of air's specific heat capacity (cf. Table 1). Clearly, in particular at elevated altitudes, the courses of θ_{1000} and θ_{1010} significantly deviate from the reference. To quantitatively evaluate the match between the different profiles, the relative difference of the profiles based on a constant c_{p}, with respect to the reference, i.e. $\mathrm{\Delta}\mathit{\theta}/{\mathit{\theta}}_{\mathrm{ref}}=\left({\mathit{\theta}}_{{c}_{p}}{\mathit{\theta}}_{\mathrm{ref}}\right)/{\mathit{\theta}}_{\mathrm{ref}}$, is depicted in Fig. 4b. The comparison demonstrates that the ${\mathit{\theta}}_{{c}_{p}}$ profiles significantly depart from the reference by about 300 K at 50 km altitude, corresponding to a relative difference of about 16 %. With both extremes of the recent constant values ${c}_{p}\in \left\{\mathrm{1000}\phantom{\rule{0.125em}{0ex}}\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}},\mathrm{1010}\phantom{\rule{0.125em}{0ex}}\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}\right\}$, the relative error level of 0.1 % is exceeded at altitudes about 5 km. While θ_{1000} continues to increasingly deviate from the reference, θ_{1010} reenters and crosses the 0.1 % relative error interval (greyshaded area) at altitudes between ∼19 and 21 km, before it reaches similar errors to the other ${\mathit{\theta}}_{{c}_{p}}$ profiles that are based on a constant c_{p}. Although the extreme values ${c}_{p}\in \left\{\mathrm{1000}\phantom{\rule{0.125em}{0ex}}\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}},\mathrm{1010}\phantom{\rule{0.125em}{0ex}}\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}\right\}$ appear in recent literature, these values may be considered unrealistic. For this reason, Fig. 4b also shows the relative deviations for the values ${c}_{p}\in \left\{\mathrm{1003.5}\phantom{\rule{0.125em}{0ex}}\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}},\mathrm{1004}\phantom{\rule{0.125em}{0ex}}\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}},\mathrm{1005}\phantom{\rule{0.125em}{0ex}}\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}\right.$, $\left.\mathrm{1006.5}\phantom{\rule{0.125em}{0ex}}\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}\right\}$, which include the recommended value of the WMO (1966) and a more realistic range; i.e. ${c}_{p}=\mathrm{1005}\phantom{\rule{0.125em}{0ex}}\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}\pm \mathrm{1.5}$ $\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}$. Notably, up to an altitude of 15 km, the reference potential temperature is comparably well matched by both the recommended θ_{1005} and θ_{1004} (based on the frequently used alternative c_{p}=1004 $\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}$; cf. Table 1). Until 15 km altitude, both constant c_{p} values lead to errors of calculated ${\mathit{\theta}}_{{c}_{p}}$ which remain comparatively small within the 0.1 % relative error interval. However, above ∼17.5 km, both θ_{1004} and θ_{1005} exceed the 0.1 % relative error interval, and further aloft their relative error with respect to the reference θ_{ref} increases rapidly.
In the context of numerical models of the atmosphere, the energy balance equation is occasionally formulated based on the potential temperature θ; thus θ constitutes a prognostic model variable. In such a case, the temperature T needs to be calculated from a given pair of values of pressure p and potential temperature θ. Using once more the defining Eq. (22), for given θ a zero of the function
is to be computed. Since Eq. (27) corresponds to the function F defined in Eq. (23) with the exception of a negative sign, the identical approximation procedure as outlined above in this section for the calculation of $\left(T,\phantom{\rule{0.125em}{0ex}}p\right)\mapsto \mathit{\theta}$ may be applied mutatis mutandis to calculate the transformation $\left(\mathit{\theta},\phantom{\rule{0.125em}{0ex}}p\right)\mapsto T$.
In any case, a certain effort is required to implement the new formulation of the potential temperature in an atmospheric model, as this equation should be based on the implicit definition Eq. (22), and such a goal may be the subject of future endeavours.
5.2 Approximations of the reference potential temperature
Of course, the previously described procedure to compute the potential temperature may appear to be anything but practical. Indeed, due to the complications inherent with

the requirement to numerically solve the integral in the function F(x) and

the need to use Newton's method for an iteration sequence to approach the zero of F(x),
a convenient approach to reassess the conventional definition of the potential temperature is not provided at all. This motivates the development of a more practical approximation of the reference potential temperature. To arrive at a practicable approximation procedure, the two principal steps in the suggested procedure are briefly outlined in the following, whereas the comprehensive details and intermediate derivation steps are found in Appendix C.
Proceeding from the definition (23) of the function F(x), the computation of the integral $\underset{x}{\overset{T}{\int}}\frac{{c}_{p}^{\mathrm{0}}\left({T}^{\prime}\right)}{{T}^{\prime}}\phantom{\rule{0.125em}{0ex}}\mathrm{d}{T}^{\prime}$ becomes the first obstacle to a practical approximation. Therefore, a plausible initial step is to replace the integral by an expression that is easier to treat. This expression may be proposed as f(T)−f(x), where the function f is defined as $f\left(x\right)={b}_{\mathrm{0}}+{b}_{\mathrm{1}}\mathrm{ln}\left(x{b}_{\mathrm{2}}\right)+{b}_{\mathrm{3}}x+{b}_{\mathrm{4}}{x}^{\mathrm{2}}$ and which is recognisable as an approximated primitive of $\frac{{c}_{p}^{\mathrm{0}}\left({T}^{\prime}\right)}{{T}^{\prime}}$; see Appendix C1. The choice of the functional form of f is motivated by the exact primitive of the integral in the case of a constant c_{p}.
As previously discussed (cf. Sect. 5.1), the formulation of a new expression for the potential temperature based on the temperaturedependent specific heat capacity c_{p}(T) requires finding the zero of the equation 0=F(x), where the function F(x) is defined in Eq. (23). Replacing the exact integral in Eq. (23) by the difference f(T)−f(x) means that F(x) is substituted by the function
Consequently, the resulting approximated reference potential temperature, i.e. the respective zero of the function $\widehat{F}\left(x\right)$, is denoted as ${\mathit{\theta}}_{\mathrm{ref}}^{\mathrm{approx}}$.
The difference between the approximation result and the reference, i.e.
is then referred to as the basic error of the approximation. Note that the replacement of the function F by $\widehat{F}$ only circumvents the integration in F; the rootfinding problem $\mathrm{0}=\widehat{F}\left(x\right)$ for the approximated reference potential temperature ${\mathit{\theta}}_{\mathrm{ref}}^{\mathrm{approx}}$ remains analytically not solvable.
Therefore, the second move towards a practical approximation procedure is to construct approximations θ^{(k)} to the zero of $\widehat{F}\left(x\right)$ by using Newton's method; see Appendix C2. Newton's method is an iterative procedure; the notation θ^{(k)} refers to the kth computed iterate. Hence, θ^{(k)} constitutes an approximation to ${\mathit{\theta}}_{\mathrm{ref}}^{\mathrm{approx}}$, and, in the limit k→∞, the approximation error
vanishes. Two formulations of Newton's method are distinguished in Appendix C2, i.e. the principal application of Newton's method, and its further derivative, called Householder's method. Both formulations require the stipulation of one of the iterates θ^{(k)} as sufficient to obtain a result of appropriate accuracy. The higher the number of iterations, of course, the smaller is the error (30), whereas the basic error (29) remains unaffected by the number of iterations. Hence, in any case, the basic error (29) is to be accepted as at least implied in the final approximation, even though a wellchosen θ^{(k)} could result in an approximation error θ_{ref}−θ^{(k)} that is smaller than the basic error.
The various errors implied in the proposed approximation procedure combining for the approximation's total error, as well as accompanying details, are discussed in Appendix D. In brief, Fig. 5a illustrates the basic error (29) based on the pressure and temperature profiles of the US Standard Atmosphere, as these provide atmospherically meaningful averages of realistic temperature–pressure data pairs. Based on the parameters of the US Standard Atmosphere, the basic error inherent with the approximation remains below 1.25 K up to altitudes of 50 km. Thus, regarding the subsequent iteration process, a substantial improvement of the error compared to ∼1.5 K is not to be expected for the total error of approximating the reference potential temperature.
An error analysis exclusively based on the US Standard Atmosphere is constrained to specific combinations of the air's pressure and temperature, potentially suppressing latent errors that may emerge if certain fluctuations of the real atmosphere's temperature and pressure profiles are considered. Thus, the error analysis is extended to an atmospheric pressure (p) and temperature (T) range, from 1000 to 0.5 hPa and from 180 to 300 K, such that the conditions within the entire troposphere and stratosphere, including the stratopause, are covered. Figure 5b illustrates the absolute basic error (29) for the extended ranges of pressure and temperature while Fig. 5c illustrates the relative basic error $\left{\mathit{\theta}}_{\mathrm{ref}}{\mathit{\theta}}_{\mathrm{ref}}^{\mathrm{approx}}\right/{\mathit{\theta}}_{\mathrm{ref}}$. The contours in Fig. 5b and c mainly highlight two regions: at ∼100 hPa where Δθ never rises above 0.75 K, which corresponds to a maximum relative basic error of 0.15 %, and in a pressure range from ∼5 to 1 hPa where a Δθ of 1.25 K is never exceeded, corresponding to relative errors of at most 0.1 %. Note that the entire Δθ scale ranges up to 3 K, which may only be reached at pressures below 0.8 hPa combined with temperatures above 280 K.
As previously discussed, the basic error is unavoidable and is to be accepted when applying the suggested substitution for the integral in the definition of the function F(x) in Eq. (23). However, as outlined in Appendix C2, the second iterate θ^{(2)} of Newton's method (principal application) may thoroughly suffice for the final approximation to the reference potential temperature θ_{ref}, as this iteration level already features an approximation error (30) which is negligibly small. Figure 6a illustrates the total relative error of the suggested approximation θ^{(2)} with respect to the ultimate reference θ_{ref} for the extended ranges of pressure and temperature. Indeed, the contour pattern in Fig. 6a and the basic relative approximation error shown in Fig. 5c are remarkably similar. Thus, the iteration process itself imparts only a minor contribution to the total error compared to the basic approximation error.
The total approximation error, which is
is dominated by the unavoidable basic error (first bracket) and augmented by a negligible error inherent to the iteration (second bracket), also supporting the conclusion that the second iterate of Newton's method is an appropriate approximation procedure. Figure 7 presents stepwise instructions for the computation of the second iterate approximation to the reference potential temperature and may serve as a guide to follow the numerous equations and intermediate analytical steps described throughout the derivations in Appendix C.
For completeness, Fig. 6b and c exhibit a final comparison by means of the logarithmic difference and the logarithmic relative difference between the reference potential temperature θ_{ref} and the conventional definition ${\mathit{\theta}}_{{c}_{p}}$ (WMO, 1966) based on a constant specific heat capacity c_{p}=1005 $\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}$. Notably, over a wide altitude range within the troposphere (i.e. for atmospheric pressures greater than ∼100 hPa), the absolute error $\mathrm{\Delta}\mathit{\theta}=\left{\mathit{\theta}}_{\mathrm{1005}}{\mathit{\theta}}_{\mathrm{ref}}\right$ remains below 1 K, cf. Fig. 6b, corresponding to a relative error Δθ∕θ_{ref} of at most 0.1 %. However, in the pressure range below ∼100 hPa, deviations of the real atmospheric conditions from those of the US Standard Atmosphere could increase the absolute error Δθ from a few K to up to 10 K, corresponding to an increase in the relative error to 1 %. Further critical pressure levels are at ∼20 and ∼5 hPa, where the error's magnitude increases to several tens and several hundreds of K, respectively. At a pressure of 0.5 hPa, an absolute error Δθ of up to 500 K is reached, which corresponds to a relative error of 10 % or even more.
5.3 Implementation aspects
The use of the new reference potential temperature θ_{ref} in a numerical model requires additional computational effort to perform corresponding calculations. Hereafter, two aspects are briefly discussed: (i) the formulation of the model equations, which include θ_{ref}, and (ii) the calculation of θ_{ref}.
Although it is beyond the scope of the present study to provide a general derivation of an appropriate energy equation based on θ_{ref} for atmospheric models, a formulation of the total derivative of θ_{ref} is given by
where the details of its derivation are given in Appendix E. The total derivative of θ_{ref} may be useful, since the governing equations are commonly formulated as differential equations.
The calculation of both the reference potential temperature θ_{ref} and its approximation ${\mathit{\theta}}_{\mathrm{ref}}^{\mathrm{approx}}$ on the basis of given values of pressure p and temperature T requires an iterative procedure. The additional computational effort inherent with these calculations depends on the number of iterations. If, however, the second iteration θ^{(2)} already represents an appropriate approximation of θ_{ref} (cf. Sect. 5.2), then the flowchart in Fig. 7 immediately conveys the additional computational effort to be expected. The calculation of the starting value x_{0} is identical to computing θ_{1005}. An additional effort results from the evaluation of the functions f (three times) and f^{′} (two times), respectively, and the combination (two times) of obtained values to determine x_{1} and x_{2}. Since each of these evaluations causes additional numerical steps, the computational effort to obtain θ^{(2)} is in total about 7 times more than the calculation of the conventional θ_{1005}, while the algorithmic complexity is constant.
To account for realgas effects (that cause a behaviour other than that of an ideal gas; cf. Sect. 4) on the potential temperature, we use the model embedded in REFPROP (Lemmon et al., 2018), a standard reference database from NIST. This model treats air as a mixture and employs stateoftheart reference equations of state for pure nitrogen (Span et al., 2000), oxygen (Schmidt and Wagner, 1985), and argon (Tegeler et al., 1999). The mixing rule and binary interaction parameters are taken from the GERG2008 model (Kunz and Wagner, 2012). From its definition in terms of an isentropic process, the potential temperature θ_{real}(T, p) is defined implicitly by
where s is the specific entropy. Calculating θ_{real}(T, p) is a twostep process. First, the specific entropy s is computed at temperature T and pressure p. Then, the temperature θ_{real} that gives the same entropy s at the ground pressure p_{0} is found. This is an iterative calculation, but it is accomplished automatically within the REFPROP software (Lemmon et al., 2018).
One caveat should be mentioned regarding the computed potential temperatures. The range of validity of the equations of state for the air components (Span et al., 2000; Schmidt and Wagner, 1985; Tegeler et al., 1999) extends only up to 2000 K. At very high altitudes, computed values of θ_{real} exceed this limit. While all the equations extrapolate in a physically realistic way, their quantitative accuracy is less certain above 2000 K. This caveat also applies to the idealgas calculations; the correlations for ${c}_{p}^{\mathrm{0}}\left(T\right)$ for N_{2} and O_{2} are extrapolations beyond 2000 K. However, since the same idealgas values are used in the realgas calculations, any inaccuracy in ${c}_{p}^{\mathrm{0}}\left(T\right)$ will cancel when evaluating the difference between idealgas and realgas values of θ.
Figure 8 illustrates the comparison between the realgas potential temperature θ_{real} and the idealgas reference potential temperature θ_{ref}. Figure 8a shows the difference θ_{real}−θ_{ref} along the p–T profile of the US Standard Atmosphere and Fig. 8b accounts for any p–T combination of extended range but shows the relative difference instead. The difference between θ_{real} and θ_{ref} never exceeds 0.1 K for the absolute difference or $\mathrm{30}\times {\mathrm{10}}^{\mathrm{5}}=\mathrm{0.03}\phantom{\rule{0.125em}{0ex}}\mathrm{\%}$ for the relative difference. As may be anticipated from the deviation of ${c}_{p}^{\mathrm{0}}$ shown in Fig. 3 at low temperatures both from the experimentally determined values (which may be inaccurate) and from the REFPROP data, the realgas effect on the specific heat capacity of dry air tends to increase towards the coldest gas temperatures. However, the difference between the real and idealgas approaches results in essentially no substantial difference between the resulting θ's, neither at ground conditions (for any temperature at ∼1000 hPa) nor at very high altitudes (at pressures below ∼1 hPa). While the negligible difference between θ_{real} and θ_{ref} near ground levels is less surprising, the diminished difference at higher altitudes reflects that in this region the potential temperature reaches such high values that the difference between the realgas and the idealgas specific heat capacity becomes insignificant. Within the intermediate (stratospheric) region, the low pressures (and thus the low air densities) cause the idealgas assumption to be an accurate approximation even at low temperatures. In general, the degree to which a gas can be treated as ideal is primarily a function of the (molar) density. For an ideal gas, the density is proportional to the quotient $\frac{p}{T}$; this is almost true also for real air. Hence, declining pressures together with rising temperatures both make the air's behaviour increasingly close to ideal.
As previously shown, the newly defined reference potential temperature θ_{ref} deviates most from the WMOdefined potential temperature θ_{1005} at stratospheric altitudes and above (cf. Fig. 6). More particularly, not only do the values from both θ definitions differ, but also their vertical derivatives, i.e. $\frac{\partial {\mathit{\theta}}_{\mathrm{ref}}}{\partial z}$ and $\frac{\partial {\mathit{\theta}}_{\mathrm{1005}}}{\partial z}$. Whether such deviations have a significant effect on an application is very case dependent and requires detailed examination and specific appraisal. Below, four typical applications of the potential temperature were selected and are examined regarding the quantitative effect on the results due to deviations of the introduced reference potential temperature compared to the conventional and commonly used θ_{1005}. The purpose of this examination is to document the magnitude of errors to allow a wellfounded, individual decision for each application of the potential temperature as to whether it is worth applying the more rigorous calculation in the particular context.
7.1 The Brunt–Väisälä frequency
The formula for the (squared) Brunt–Väisälä frequency N^{2} is often given in the form of Eq. (2), i.e. a formula involving the potential temperature θ. The substitution of θ in Eq. (2) by the new reference potential temperature θ_{ref} may be tempting, but it is erroneous, and the resulting quantity is denoted as ${N}_{\mathrm{false}}^{\mathrm{2}}$. The Brunt–Väisälä frequency is not defined by Eq. (2), since this formula results from various simplifications in its derivation, e.g. by assuming hydrostatic conditions and a constant specific heat capacity. Consequently, the substitution of θ_{ref} in Eq. (2) leads to a wrong formula for the Brunt–Väisälä frequency that does not correctly consider the temperature dependence of dry air's specific heat capacity.
The Brunt–Väisälä frequency is the oscillation frequency of an air parcel due to a local density perturbation (see e.g. Durran and Klemp, 1982; Marquet and Geleyn, 2013; Wallace and Hobbs, 2006; Ambaum, 2010). Retaining the assumption of hydrostatic conditions, the defining formula yields
where the temperaturedependent specific heat capacity c_{p}(T) was implied, and which quantifies the balance between the actual temperature stratification $\frac{\partial T}{\partial z}$ and the dryadiabatic lapse rate $\frac{g}{{c}_{p}\left(T\right)}$ (e.g. Holton, 2004).
To illustrate the deviation of ${N}_{\mathrm{false}}^{\mathrm{2}}$ from N^{2}, vertical profiles of both variables were calculated based on the temperature profiles shown in Fig. 9a. The temperature data are taken from the Upper Atmosphere Research Satellite Reference Atmosphere Project (URAP; see Swinbank and Ortland, 2003) data and are assumed as typical at midlatitudes during June and December. The temperature profiles extend up to altitudes of 85 km and thus cover the entire stratosphere and most of the mesosphere. The hydrostatic assumption allowed for computing pressure profiles along the URAP values for the vertical temperature distribution. Subsequently, the reference potential temperature θ_{ref} and its vertical derivative were calculable. The resulting vertical profiles for ${N}_{\mathrm{false}}^{\mathrm{2}}$ and the true Brunt–Väisälä frequency N^{2} are shown in Fig. 9b. Evidently, the values of ${N}_{\mathrm{false}}^{\mathrm{2}}$ (dashed lines) deviate significantly from N^{2} (solid lines) and increasingly so towards higher altitudes above 15 km. However, the absolute deviation ${N}^{\mathrm{2}}{N}_{\mathrm{1005}}^{\mathrm{2}}$, using ${N}_{\mathrm{1005}}^{\mathrm{2}}$ as calculated with θ_{1005} in accordance with Eq. (2), does not exceed $\mathrm{1.6}\times {\mathrm{10}}^{\mathrm{6}}\phantom{\rule{0.125em}{0ex}}{\mathrm{s}}^{\mathrm{2}}$ (not shown), indicating that ${N}_{\mathrm{1005}}^{\mathrm{2}}$ is a good representation of N^{2} along these temperature profiles.
For equations involving the potential temperature, however, it should be emphasised that the substitution of θ by θ_{ref} rarely succeeds and that instead the entire derivation of the equations requires careful consideration of the assumptions, such as the constancy of c_{p}, to avoid aberrations and erroneous conclusions.
7.2 The potential vorticity
Ertel's potential vorticity (e.g. Ertel, 1942; Hoskins et al., 1985; Schubert et al., 2004; Holton, 2004) may be defined as the potential vorticity of the dryair potential temperature by
In this definition, $\mathrm{2}\mathbf{\Omega}+\mathrm{\nabla}\times \mathit{u}$ is the absolute vorticity, Ω denotes Earth's angular velocity, u the threedimensional wind vector, and ρ the air density (see e.g. Hoskins et al., 1985; Cotton et al., 2011; Marquet, 2014). Since Eq. (35) represents the defining equation for Ertel's potential vorticity, the two potential vorticities
based on the new reference potential temperature θ_{ref} and θ_{1005}, respectively, are considered. To provide a first comparison of these potential vorticities, u=0 is assumed, i.e. an atmosphere at rest. Additionally, the potential temperature is assumed as horizontally constant. Consequently, Eq. (35) reduces to
for a position on Earth with geographical latitude ϕ and t_{E}=24 h, the duration of one rotation of the Earth.
Using the temperature profiles from Fig. 9a together with the values of the potential temperatures θ_{ref} and θ_{1005}, the evaluation of the two potential vorticities (36) by (37) yields the potential vorticity profiles shown in Fig. 10a while their relative deviations are shown in Fig. 10b. Since the temperature profiles are representative for the northernhemispheric midlatitudes, the geographical latitude ϕ in Eq. (37) was set to 52^{∘} N. At tropospheric altitudes, the relative deviation between θ_{ref} and θ_{1005} is small and never exceeds ∼1 %, while it continuously increases towards higher altitudes. According to these profiles, the relative deviation exceeds 10 % at 30 km and reaches 100 % at the highest altitudes.
It is noteworthy, however, that the computations of both N^{2} (cf. Sect. 7.1 and Fig. 9b) and PV (Fig. 10b) are based on the specific temperature profiles from URAP (cf. Sect. 7.1 and Fig. 9a) and thus are not of general validity. The selection of these temperature profiles was entirely arbitrary and exclusively aimed at illustrating possible implications of the use of the developed reference potential temperature. The resulting and indicated deviations are ultimately subject to individual assessment on applying θ_{ref}.
7.3 Vertical sorting of data
For atmospheric investigations, e.g. in the region of the upper troposphere and lower stratosphere (UT/LS), it is common practice to set vertical profiles of atmospheric parameters in relation to the potential temperature as a vertical coordinate. This way, the increasingly isentropic stratification of the atmosphere above the UT is taken into account. The transport of an air mass along isentropic surfaces, i.e. surfaces of constant potential temperature and entropy, is to be regarded as adiabatic. Hence, the air's composition and properties within the same isentrope interval, regardless of the observation location, are better comparable than they would be if based on other isopleths (i.e. height or pressure coordinates). Investigations of air mass compositions over time and from different regions at the same θ level largely exclude that, during its transport history, the air had experienced vertical displacement and/or diabatic processes (radiative heating, condensation/evaporation) which would result in energy conversion. The tropopause height is often used as a reference height in the θ coordinate system in connection with the vertical sorting of observational data, whereby the assignment of tropospheric and stratospheric processes is made or exchanges across the tropopause are investigated (Holton et al., 1995; Stohl et al., 2003). Consequently, the tropopause height is also determined by the potential vorticity (e.g. Gettelman et al., 2011, and cf. Sect. 7.2) if the conventional tropopause definitions (cold point or lapse rate, WMO, 1957) do not allow for clearly determining the tropopause height, e.g. in the Asian monsoon anticyclone (cf. Höpfner et al., 2019) or in the polar winter vortex (Wilson et al., 1989; Weigel et al., 2014). The conventional definition of θ implies a systematic error in the vertical sorting of observational data in the θ coordinate system, independent of the measurement platform. Investigations with highaltitude research aircraft such as the G550 HALO (e.g. Wendisch et al., 2016; Voigt et al., 2017), the NASA WB57 or ER2 (e.g. Murphy et al., 2007; Dessler, 2002), the M55 Geophysica (Curtius et al., 2005; Borrmann et al., 2010; Frey, 2011), balloonborne platforms (Lary et al., 1995; Vernier et al., 2018), or satellitebased vertical profiles (e.g. Davies et al., 2006; Spang et al., 2005) require consideration of the systematic error in θ if calculated as ${\mathit{\theta}}_{{c}_{p}}$ in compliance with the definition by the WMO (1966). The possibly inconsistent use of a constant c_{p} value of 1004 or 1005 $\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}$ (or any other) in different and compared data sets, which could be due to different literature references for this value (cf. Table 1), will not be explored here. At altitudes between 15 and 20 km (ceiling of highaltitude research aircraft), an overestimation by about 0.1 %–0.5 % is to be expected for the potential temperature according to the conventional definition (cf. Fig. 4b). At altitudes of 30–35 km, an overestimation by up to 2 %–5 % results. Whether this error is significant or small compared to the uncertainty of ambient temperature and pressure measurement aboard the respective aircraft is left to individual judgement in the course of data processing. In the case of spacecraftbound vertical soundings (e.g. from ASTROSPAS, SCIAMACHY, or ENVISAT), the error in the potential temperature determined by ${\mathit{\theta}}_{{c}_{p}}$ exceeds 10 % at altitudes above 40 km, as shown in Fig. 4b. Finally, we note that the specified errors apply exclusively along the vertical profile of the US standard atmosphere and that deviations of the actual temperature profile from the US standard atmosphere, e.g. warmer temperatures, could lead to larger errors (cf. Fig. 6).
7.4 Diabatic heating rates
Diabatic heating rates refer to the rate of energy $\frac{\mathrm{d}q}{\mathrm{d}t}$ supplied to a given air parcel, e.g. by radiative heating, and are given in units of $\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{s}}^{\mathrm{1}}$. This energy supply causes a temperature change in an air parcel at a rate which hereafter is referred to as the absolute heating rate,
Again, the distinction was made between the temperaturedependent ${c}_{p}^{\mathrm{0}}\left(T\right)$ and the constant c_{p}=1005 $\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}$ specific heat capacity. From the defining Eq. (38), the relative difference between these absolute heating rates, where x designates an arbitrary diabatic heating rate, is
Apart from the absolute heating rates for the change in absolute temperature, the change in potential temperature due to a diabatic heating rate $\frac{\mathrm{d}q}{\mathrm{d}t}$ is of interest. For example, it is the change in potential temperature that modifies the altitude of modelled trajectories in Lagrangian chemical transport models based on isentropic coordinates rather than the change in absolute temperature (e.g. the SLIMCAT model, Chipperfield, 2006; or CLaMS model, Pommrich et al., 2014).
Taking the relation T ds=dq for the specific entropy into account, Gibbs' Eq. (8) may be rewritten as
Comparing the righthand side of this equation to the total derivative of the new reference potential temperature θ_{ref} (see Appendix E for the detailed computation and Eq. E6 for the result) Eq. (40) amounts to
Consequently, the following two diabatic heating rates,
for the potential temperatures θ_{ref} and θ_{1005} may be defined. Denoting again by x an arbitrary diabatic heating rate, the relative difference between the heating rates (Eq. 42) is
In order to judge the magnitudes of the relative differences (Eqs. 39 and 43), the monthly averaged temperature profiles from ERAInterim (Dee et al., 2011) data for 52^{∘} N geographical latitude are used; see Fig. 11a. The relative differences of the absolute heating rates (Eq. 39) are shown in Fig. 11b, and the difference appears to be small. However, the relative differences of the heating rates (Eq. 43) in Fig. 11c are much larger, as relative deviations exceeding 50 % are reached in the upper stratosphere and lower mesosphere (at pressures below 1 hPa). Additionally, the temperatures that resulted after 24 h of heating with a constant heating rate $\frac{\mathrm{d}q}{\mathrm{d}t}$ as given in the (averaged) data set were computed, where a constant pressure is assumed for simplicity. As may be anticipated from the small deviations in Fig. 11b, the differences in the final absolute temperatures by using the absolute heating rates AHR_{1005} or AHR_{ref} are smaller than 0.044 K. However, the differences in the potential temperatures ${\mathit{\theta}}_{\mathrm{1005}}^{\ast},\phantom{\rule{0.125em}{0ex}}{\mathit{\theta}}_{\mathrm{ref}}^{\ast}$, computed with the heating rates HR_{ref}, HR_{1005}, are much larger (Fig. 11d) and amount to about 3 % at 10 hPa and about 15 % at 1 hPa. For transport calculations done in isentropic coordinates, these differences are of the same order of magnitude as the deviations resulting from the use of the temperaturedependent instead of the constant c_{p}. It remains to be decided on individual application whether this additional effect in the calculation is significant.
A standard diagnostic for the speed of the stratospheric circulation is the time lag of the upwardpropagating seasonal signal in tropical stratospheric water vapour (the socalled tape recorder, Mote et al., 1996). Here, differences between calculations (done in isentropic coordinates) based on different current meteorological reanalysis data sets amount to about 10 %–30 % for the signal's upward propagation below about 10 hPa (Tao et al., 2019), such that the additional deviation from using the temperaturedependent c_{p} is comparably small. However, in cases of smaller intermodel differences the additional c_{p}related uncertainty needs to be assessed.
Note that the determination of absolute temperatures ${T}_{\mathrm{1005}}^{\ast},\phantom{\rule{0.125em}{0ex}}{T}_{\mathrm{ref}}^{\ast}$ which correspond to the resulting potential temperatures ${\mathit{\theta}}_{\mathrm{1005}}^{\ast},\phantom{\rule{0.125em}{0ex}}{\mathit{\theta}}_{\mathrm{ref}}^{\ast}$ after 24 h differs by less than 0.014 K (not shown).
Under the assumption that dry air is an ideal gas, a reassessment of computing the potential temperature was introduced that accounts for the hitherto unconsidered temperature dependence of air's specific heat capacity. The new reference potential temperature θ_{ref} was introduced, which is thermodynamically consistent and based on a stateoftheart parameterisation of the idealgas specific heat capacity of dry air from the National Institute of Standards and Technology (NIST). This reference potential temperature was compared to a potential temperature θ_{real} wherein the realgas behaviour of dry air is considered. In the range of temperatures from 180 to 300 K and the range of pressures from 1000 to 0.5 hPa, covering the atmospheric conditions of roughly the entire troposphere and stratosphere, the relative differences between θ_{ref} and θ_{real} are smaller than 0.03 % and may be considered negligible. Consequently, θ_{ref} even provides a reasonable approximation to the potential temperature of the real gas.
The difference between the newly derived reference potential temperature θ_{ref} and the conventionally determined potential temperature ${\mathit{\theta}}_{{c}_{p}}$ (with constant c_{p}=1005 $\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}$, as recommended by the World Meteorological Organisation, WMO, 1966) increases with altitude, e.g. Δθ≥1 K at pressures p≤60 hPa.
Derivation of a potential temperature that is consistent with thermodynamics and that accounts for the idealgas properties of dry air requires the integration of Gibbs' equation and the subsequent solution of the resulting nonlinear equation. With a constant c_{p}, both analytical steps are straightforward, resulting in the conventional expression (13) as suggested by WMO (1966). However, if instead the temperature dependence of air's specific heat capacity c_{p}(T) is considered, the integrals as well as the equations are not analytically solvable, and, thus, the solution must be approximated. Both approximations were performed and described in detail. The integral was treated with the basic approximation and the solution of the nonlinear equation was approximated by the second iterate of Newton's method. As an alternative to Newton's classical method, a modified formulation of Householder's iteration method is provided, featuring accelerated convergence properties.
The suggested approximation steps to obtain a reference potential temperature have two main sources of error: the error ${\mathit{\theta}}_{\mathrm{ref}}{\mathit{\theta}}_{\mathrm{ref}}^{\mathrm{approx}}$ inherent in the integral's basic approximation and the error ${\mathit{\theta}}_{\mathrm{ref}}^{\mathrm{approx}}{\mathit{\theta}}^{\left(k\right)}$ of the kth Newton iterate. The latter error approaches zero as k→∞, whereas the error resulting from the basic approximation remains well below 0.1 % (along the US Standard Atmosphere) for values of θ_{ref} of up to ∼2000 K, hence up to stratopause altitudes. To keep this low error level also for θ_{ref}>2000 K, the approximation may require an extension by means of a higherorder polynomial.
One of the foremost implications of the reassessed potential temperature's definition concerns the use of θ as a vertical coordinate for the sorting, grouping, and comparison of (measured) data, e.g. along or across isentropes. Thereby, the reassessed potential temperature constitutes a more accurate consideration of the air's actual properties. This particularly concerns for example the specific heat capacity which is conventionally assumed as constant and for which various values are given depending on the textbook consulted (offering a range from 1000 to 1010 $\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}$; see Table 1).
Significant errors and biases may arise if, for instance, the conventional derivation of θ (WMO, 1966) is used together with values for air's specific gas constant (R_{a}) or air's specific heat capacity (c_{p}), which better comply with the most recent state of knowledge. Moreover, the use of the standard pressure 1013.25 instead of 1000 hPa as defined by WMO (1966) and consistently used herein as ground level pressure (p_{0}) may cause an additional deviation of the resulting θ. Thus, the reassessment of θ's definition could largely diminish such errors and biases and improve the comparability of data.
In addition to the vertical sorting of data, implications of the new reference potential temperature were discussed for several other applications in which the potential temperature is used. On the one hand, results may appear mostly unaffected by using θ_{ref} instead of the conventional θ_{1005}, such as the values of the Brunt–Väisälä frequency or the temperature change in air parcels due to diabatic heating. On the other hand, it was illustrated that any formula which involves the potential temperature needs to be carefully reviewed to see if its derivation relies on the assumed constancy of the specific heat capacity. If this is the case, substituting θ_{ref} for all occurrences of θ within the particular formula may lead to a wrong computation.
In contrast, examples were shown where the computation of Ertel's potential vorticity and the rate of change of potential temperature in response to diabatic heating yields different results by the use of θ_{ref} instead of θ_{1005}. The differences increased with altitude; hence they become more important for applications within the stratosphere and above.
It should be emphasised that all these examples were based on assuming particular profiles of temperature and pressure together with other assumptions. Moreover, only a limited number of examples could be investigated, while the applications of potential temperature are numerous. Consequently, a wellfounded, individual decision is required for each application of the potential temperature as to whether it is worth applying the more rigorous calculation in the particular context.
On the one hand, such a reassessment could take into account the current state of knowledge regarding the accuracy of thermodynamic variables and substancerelated properties. On the other hand, this way, the conceptional abstractness already inherent in θ is not further complicated by a misleading selection of parameters or reputed constants. There is no doubt that the conventional method is suitable for the description of most processes occurring within the troposphere. However, at stratospheric or even mesospheric altitudes, the neglect of the temperature dependence of the idealgas heat capacity in the conventional definition increasingly distorts the resulting absolute values as well as the vertical course of the potential temperature. Ultimately, it seems obvious to profit from the computing capacities available today and from the known higher accuracy of physical variables and atmospheric parameters to carry out a reappraisal of the potential temperature, a useful (but not always consistently used) meteorological quantity.
In the following, the derivation of the air's specific heat capacities C_{V}, C_{p} (capital letters indicate molar units) at constant volume and pressure, respectively, is summarised, mainly following the textbook exposition by Kondepudi and Prigogine (1998). We start with the idealgas law
with p the pressure, V the volume of the system, N the amount of gas within the volume, T the temperature, and R the universal gas constant. Additionally, the first law of thermodynamics is
with U the internal energy of the system and dQ specifying the change in heat. Insertion of the total derivative of the internal energy U in Eq. (A2), and assuming the system as thermodynamically closed, i.e. the molar amount N remains conserved (dN=0), leads to
and subsequently
If the system's volume is held constant, Eq. (A4) represents the definition of the constantvolume heat capacity C_{V} in molar units, i.e.
Alternatively, assuming the system's pressure as constant, its volume is variable with total derivative
and therefore
defining the isobaric molar heat capacity C_{p}. In general, this quantity depends on pressure as well as on temperature. However, if the gas is assumed as ideal, an important conclusion from the statistical description of an ideal gas is the fact that the internal energy U must be independent of the pressure (see e.g. Fay, 1965).
Using this result, together with Eq. (A7) and the idealgas law (A1), it follows
In the previous computations, there is no restriction on the temperature dependence of the internal energy U(T). Therefore, even by assuming idealgas behaviour, the specific heat capacity C_{p} in Eq. (A8) is in general a function of temperature.
This section explores, from a mathematical perspective, the sensitivity of the potential temperature formulation (13) based on a constant specific heat capacity. Considering the specific heat capacity c_{p} as a variable, the sensitivity of ${\mathit{\theta}}_{{c}_{p}}$ (Eq. 13) to a small perturbation δ of c_{p} is described by its Taylor expansion
For any constant value of the specific heat capacity c_{p} and for a minor perturbation δ, the second summand within the expansion (B1) remains small for small values of $\mathrm{ln}\left(\frac{{p}_{\mathrm{0}}}{p}\right)$. If the interval between the two pressure levels is very narrow, i.e. p≈p_{0}, the expression $\mathrm{ln}\left(\frac{{p}_{\mathrm{0}}}{p}\right)$ approximately equals ln (1)=0. Contrarily, if the pressure approaches very low values, i.e. p→0 Pa, the logarithmic expression diverges to negative infinity, i.e. $\mathrm{ln}\left(\frac{{p}_{\mathrm{0}}}{p}\right)\to \mathrm{\infty}$, implying that the impact of the second summand intensifies with decreasing pressure, i.e. for increasing altitudes. Moreover, this may explain why the deviation between θ_{1000} and θ_{1010}, as illustrated in Fig. 2b, remains comparatively small within the troposphere and systematically increases with rising altitude, i.e. decreasing pressure levels.
This section summarises the detailed steps of approximating the function F(x), defined in Eq. (23), by $\widehat{F}\left(x\right)$, defined in Eq. (28) (Sect. C1), as well as the approximations of the solutions of the resulting nonlinear equations by Newton's method (Sect. C2).
C1 Reformulating the function F(x)
Proceeding from the definition of a function h(x),
with T_{1}=180 K, the function F(x) may be rearranged as
The advantage of this reformulation of F(x) is the inclusion of h(x), consisting of an integral with fixed lower bound and a sole variable upper bound. This way, the function h(x) is numerically solvable, and subsequently h(x) can be substituted by an approximation f(x) that is defined as
Notably, if c_{p} is constant, this function reduces to an exact primitive of the integrand $\frac{{c}_{p}}{{T}^{\prime}}$ with ${b}_{\mathrm{3}}={b}_{\mathrm{4}}=\mathrm{0}$. Moreover, in this case, the resulting rootfinding problem 0=F(x) is exactly solvable and finally leads to the known conventional definition (13) of the potential temperature.
As a further step, the function h(x) is numerically approximated, while c_{p}(T) in Eq. (C1) is replaced by the idealgas limit of air's specific heat capacity ${c}_{p}^{\mathrm{0}}\left(T\right)$. The integration interval [T_{1}, x] with ${T}_{\mathrm{1}}\le x\le \mathrm{2000}\phantom{\rule{0.125em}{0ex}}\mathrm{K}$ is traversed in steps of at most 0.001 K, while each step of the integration process is carefully approximated by using Simpson's rule.
By solving a leastsquares problem, the coefficients in Eq. (C3) for the approximation of h(x) by the function f(x) are estimated as
In Fig. C1a the function h(x) is graphed together with the approximation f(x), as well as the respective deviations h(x)−f(x) in Fig. C1b. Evidently, the absolute error inherent to the approximations is comparatively small as, over the entire temperature range above 190 K, the approximation error never exceeds ±1 $\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}$. Exclusively at temperatures below 190 K, the approximation error rapidly rises above 1 $\mathrm{J}\phantom{\rule{0.125em}{0ex}}{\mathrm{kg}}^{\mathrm{1}}{\mathrm{K}}^{\mathrm{1}}$, bearing in mind that such absolute temperatures are only occasionally found in the atmosphere within a relatively narrow altitude interval at the cold point tropopause. Moreover, the difference between f(x) and h(x) appears negligible as the profiles almost ideally coincide (cf. Fig. C1a).
C2 Finalised approximation of the reference potential temperature
As discussed in Sect. 5.1, the new formulation of the potential temperature based on the temperaturedependent specific heat capacity c_{p}(T) requires solving the rootfinding problem 0=F(x), where the function F(x) is defined in Eq. (23). However, since F(x) contains an integral that complicates the rootfinding process, this integral is substituted by the difference f(T)−f(x), where f is given in Sect. C1. Therefore, F(x) is replaced by the function $\widehat{F}\left(x\right)$ as defined in Eq. (28), and the zero of the equation $\mathrm{0}=\widehat{F}\left(x\right)$ is denoted as ${\mathit{\theta}}_{\mathrm{ref}}^{\mathrm{approx}}$.
The equation $\mathrm{0}=\widehat{F}\left(x\right)$ is still not analytically solvable, so Newton's method is once more required. Using again x_{0}=θ_{1005} as the initial guess, cf. Eq. (26), the iteration sequence for Newton's method is given by the recursion
Instead of this standard formulation of Newton's method (C5), Householder's formulation
may be used, which allows for reducing the computation time due to its accelerated convergence speed. For completeness, the required derivatives f^{′} and f^{′′} in the recursion formulas (C5) and (C6) are
The final step on the way to formulate a new expression for the potential temperature requires defining one of the iterates x_{k} as appropriate enough for the approximations that result from applying the different methods:

the standard of Newton's method (C5), simply referred to as Newton's method in the sequel, or

Householder's method (C6).
While the mathematical expressions in Eqs. (C5) and (C6) are of increasing complexity, the convergence rate of the approximating sequence increases with rising mathematical complication. The preferred method is determined by the accuracy required; i.e. better accuracy is necessarily associated with elevated computational effort for the approximation method. A discussion of the approximation errors is found in Appendix D.
Table C1 collects values of the new reference potential temperature θ_{ref}, together with the first two iterates θ^{(1)}, θ^{(2)} using Newton's method (C5) and the first iterate ${\mathit{\theta}}_{\mathrm{Householder}}^{\left(\mathrm{1}\right)}$ using Householder's method (C6) for five pairs of temperature and pressure along the US Standard Atmosphere, cf. Fig. 1, which allows verification of computations. The first height is chosen midway along the linearly decreasing temperature profile within the troposphere, while the other heights correspond to the kinks of the temperature profile.
The following aims at a comprehensive investigation of the errors inherent with approximating the ultimate reference potential temperature θ_{ref}. As discussed in Sect. 5.2, the total error is a combination of the basic error ${\mathit{\theta}}_{\mathrm{ref}}{\mathit{\theta}}_{\mathrm{ref}}^{\mathrm{approx}}$ and the approximation error that results from the approximation sequence ${\mathit{\theta}}_{\mathrm{ref}}^{\mathrm{approx}}{\mathit{\theta}}^{\left(k\right)}$, where θ^{(k)} denotes the kth iterate of the approximation sequence which is computed in accordance with either Newton's or Householder's method. The formulations of Newton's (Eq. C5) and Householder's (Eq. C6) method require replacing the function F(x) by $\widehat{F}\left(x\right)$, and the approximation sequences θ^{(k)} converge to ${\mathit{\theta}}_{\mathrm{ref}}^{\mathrm{approx}}$ for k→∞. Consequently, the approximation error ${\mathit{\theta}}_{\mathrm{ref}}^{\mathrm{approx}}{\mathit{\theta}}^{\left(k\right)}$ tends to zero for k→∞.
The analysis of the approximation error is initially based on the pressure and temperature profiles of the US Standard Atmosphere. Figure D1 shows the total relative errors $\left({\mathit{\theta}}_{\mathrm{ref}}{\mathit{\theta}}^{\left(\mathrm{1}\right)}\right)/{\mathit{\theta}}_{\mathrm{ref}}$ of the first iterate (Fig. D1a) and $\left({\mathit{\theta}}_{\mathrm{ref}}{\mathit{\theta}}^{\left(\mathrm{2}\right)}\right)/{\mathit{\theta}}_{\mathrm{ref}}$ of the second iterate (Fig. D1b), computed with Newton's or Householder's method. The first iterate still causes the approximation to have significant errors, especially at altitudes above 35 km. However, the second iterate with either Newton's or Householder's method yields results with negligible approximation error. Hence, the total error of the approximation procedure is dominated by the unavoidable basic error and may be deduced from the provided figures whenever the total error profile nearly congruently follows the profile of the basic error (cf. Figs. D1b and 5a).
It may be noted that Householder's method achieves a significantly lower error level than Newton's method due to its accelerated rate of convergence. Compared to the first iterate approximations, computation up to the second iterate (cf. Fig. D1b) achieves, in general, a considerable improvement for both methods, and both second iterate approximations approach the basic error quite closely (cf. Fig. D1b). As is also evident from Fig. D1b, compared to Householder's method, the second iterate with Newton's method results in a smaller total relative error $\left({\mathit{\theta}}_{\mathrm{ref}}{\mathit{\theta}}^{\left(\mathrm{2}\right)}\right)/{\mathit{\theta}}_{\mathrm{ref}}$ relative to the ultimate reference potential temperature (indicated by a smaller distance to the dashed zeroline above 45 km altitude). Nevertheless, the relative approximation error, $\left({\mathit{\theta}}_{\mathrm{ref}}^{\mathrm{approx}}{\mathit{\theta}}^{\left(\mathrm{2}\right)}\right)/{\mathit{\theta}}_{\mathrm{ref}}$, is larger compared to the second iterate with Householder's method. So, luckily, the second iterate with Newton's method provides a better approach to the reference potential temperature than that with Householder's method.
As with the discussion of the basic error in Sect. 5.2, the analysis of the total error should include all possible combinations of pressure and temperature in order to take into account fluctuations in the real atmosphere that deviate from the profile of the US Standard Atmosphere. Therefore, the extended analysis of the approximation error is summarised in Fig. D2. The upper panels illustrate the total relative error of the second iterate for Newton's (Fig. D2a) and Householder's method (Fig. D2b). As previously shown, further iteration with either method does not improve the approximation quality. The contour patterns in these panels show a remarkable similarity to the contours for the relative error of the basic approximation in Fig. 5c. Also here (upper panels of Fig. D2), two regions are highlighted by the contours, i.e. at ∼100 hPa and in a pressure range from ∼5 to 1 hPa, featuring the same impact on Δθ∕θ_{ref} of identical strength as the basic error. This result may not be surprising, since the second iteration step with both methods, Newton's and Householder's, was already proven to approach the approximation comparatively well, without worsening the total error level (cf. Fig. D1b).
Consequently, concerning the required number of iterations and the method to use, the second iteration of Newton's method can be recommended to deliver appropriate results, with a relative error of less than 0.3 %, up to the stratopause level (∼50 km). Householder's method features an accelerated convergence rate, and its use up to its first iterate θ^{(1)} may be already appropriate for certain applications. According to the total error of Householder's method up to its first iterate θ^{(1)} (Fig. D2c), the resulting relative error remains below 7 % to a pressure level of ∼50 hPa and Δθ stays below 0.3 % to pressures of ∼2 hPa. Thus, Fig. D2 may serve as guidance to decide how many iterations with one or the other method best meet the individual accuracy requirements.
As discussed in Sect. 5.1, the new reference potential temperature is defined as the zero of the function
for given values of pressure p and temperature T; see Eq. (23). More precisely, for varying p, T, a function $(p,\phantom{\rule{0.125em}{0ex}}T)\mapsto {\mathit{\theta}}_{\mathrm{ref}}(p,\phantom{\rule{0.125em}{0ex}}T)$ is implicitly defined by the equation
According to the implicit function theorem (e.g. Protter and Morrey, 1985, chap. 7), Eq. (E2) is uniquely solvable for θ_{ref}(p, T); i.e. the function $(p,\phantom{\rule{0.125em}{0ex}}T)\mapsto {\mathit{\theta}}_{\mathrm{ref}}(p,\phantom{\rule{0.125em}{0ex}}T)$ actually exists as a differentiable function of (p, T), if the condition $\frac{\partial F}{\partial \mathit{\theta}}\ne \mathrm{0}$ holds. According to Eq. (E1), this partial derivative equals
being strictly negative, since the specific heat capacity is always positive. Moreover, the implicit function theorem states that the derivatives of the implicit function θ_{ref}(p, T) are given by
Note that these partial derivatives coincide with the partial derivatives of ${\mathit{\theta}}_{{c}_{p}}$ in the case of a constant specific heat capacity. Using the partial derivatives (Eq. E4), the total differential of θ_{ref} may be written as
or
The data used in this study are available from the corresponding author upon request (manuel.baumgartner@unimainz.de).
MB, RW, and PS conceived, designed, and carried out the main part of the research. UA contributed the URAP data and implications on breaking heights of gravity waves. AHH gave advice about the heat capacity and performed the calculations of realgas potential temperatures. FP contributed the implications on diabatic heating. MB and RW wrote the manuscript with contributions and reviews from all authors.
The authors declare that they have no conflict of interest.
We thank Eric W. Lemmon for his advice on the equation of state of dry air, Vera Bense for fruitful discussions on gravity wave breaking, Gergely Bölöni for providing us the URAP vertical profiles in Sect. 7, Miklós Szakáll for the hint concerning Titan's atmosphere, and Heini Wernli for his interest in an early version of the draft. We thank Pascal Marquet and an anonymous reviewer for their very valuable comments and hints which led to a significantly improved manuscript. Additionally, we gratefully acknowledge the contributions of Timothy Garret as editor in the review process of this article. Manuel Baumgartner and Peter Spichtinger acknowledge support by the Deutsche Forschungsgemeinschaft (DFG) within the Transregional Collaborative Research Centre TRR165 Waves to Weather (https://www.wavestoweather.de, last access: 2 November 2020), projects B7 and Z2. Ralf Weigel received financial support by the Bundesministerium für Bildung und Forschung (BMBF) under the joint ROMIC project SPITFIRE (01LG1205A). Ulrich Achatz and Peter Spichtinger acknowledge partial support by the DFG through the research unit Multiscale Dynamics of Gravity Waves (MSGWaves) and through grants AC 71/122 and SP 1163/52.
This research has been supported by the Deutsche Forschungsgemeinschaft (grant no. 257899354, AC 71/122, and SP 1163/52) and the Bundesministerium für Bildung, Wissenschaft, Forschung und Technologie (grant no. 01LG1205A).
This paper was edited by Timothy Garrett and reviewed by Pascal Marquet and one anonymous referee.
Ambaum, M. H. P.: Thermal Physics of the Atmosphere, John Wiley & Sons, Ltd., Chichester, UK, https://doi.org/10.1002/9780470710364, 2010. a, b, c
Awano, S.: JSDiagrams for Air, Report of Aeronautical Research Institute, Tokyo Imperial University, 11, available at: https://jaxa.repo.nii.ac.jp/?action=repository_uri&item_id=35290&file_id=31&file_no=1 (last access: 11 December 2020), 1936. a, b, c
Bauer, L. A.: The relation between “potential temperature” and “entropy”, Phys. Rev., 26, 177–183, 1908. a
Bohren, C., Albrecht, B., and Albrecht, P.: Atmospheric Thermodynamics, Oxford University Press, New York, USA, Oxford, UK, 1998. a
Bolton, D.: The Computation of Equivalent Potential Temperature, Mon. Weather Rev., 108, 1046–1053, https://doi.org/10.1175/15200493(1980)108<1046:TCOEPT>2.0.CO;2, 1980. a
Borchert, S., Zhou, G., Baldauf, M., Schmidt, H., Zängl, G., and Reinert, D.: The upperatmosphere extension of the ICON general circulation model (version: uaicon1.0), Geosci. Model Dev., 12, 3541–3569, https://doi.org/10.5194/gmd1235412019, 2019. a
Borrmann, S., Kunkel, D., Weigel, R., Minikin, A., Deshler, T., Wilson, J. C., Curtius, J., Volk, C. M., Homan, C. D., Ulanovsky, A., Ravegnani, F., Viciani, S., Shur, G. N., Belyaev, G. V., Law, K. S., and Cairo, F.: Aerosols in the tropical and subtropical UT/LS: insitu measurements of submicron particle abundance and volatility, Atmos. Chem. Phys., 10, 5573–5592, https://doi.org/10.5194/acp1055732010, 2010. a
Brasseur, G. P. and Solomon, S.: Aeronomy of the Middle Atmosphere, Springer, https://doi.org/10.1007/1402038240, 2005. a
Brusseau, M., Pepper, I. L., and Gerba, C. P.: Environmental and Pollution Science, 3rd Edn., Elsevier, https://doi.org/10.1016/C20170004809, 2019. a
Bücker, D., Span, R., and Wagner, W.: Thermodynamic Property Models for Moist Air and Combustion Gases, J. Eng. Gas Turb. Power, 125, 374–384, https://doi.org/10.1115/1.1520154, 2002. a, b, c, d, e, f, g
Catling, D. C.: 10.13 – Planetary Atmospheres, in: Treatise on Geophysics, edited by: Schubert, G., 429–472, Elsevier, Oxford, 2nd Edn., https://doi.org/10.1016/B9780444538024.001858, 2015. a, b
Chang, X., Zhao, W., Zhang, Z., and Su, Y.: Sap flow and tree conductance of shelterbelt in arid region of China, Agr. Forest Meteorol., 138, 132–141, https://doi.org/10.1016/j.agrformet.2006.04.003, 2006. a
Chipperfield, M. P.: New version of the TOMCAT/SLIMCAT offline chemical transport model: Intercomparison of stratospheric tracer experiments, Q. J. Roy. Meteor. Soc., 132, 1179–1203, https://doi.org/10.1256/qj.05.51, 2006. a
Cotton, W. R., Bryan, G. H., and van den Heever, S. C.: Storm and Cloud Dynamics, 2nd Edn., Academic Press, San Diego, CA, USA, 2011. a, b
Curry, J. and Webster, P.: Thermodynamics of Atmospheres and Oceans, Vol. 65 of International Geophysics, Elsevier, London, UK, 1998. a
Curtius, J., Weigel, R., Vössing, H.J., Wernli, H., Werner, A., Volk, C.M., Konopka, P., Krebsbach, M., Schiller, C., Roiger, A., Schlager, H., Dreiling, V., and Borrmann, S.: Observations of meteoric material and implications for aerosol nucleation in the winter Arctic lower stratosphere derived from in situ particle measurements, Atmos. Chem. Phys., 5, 3053–3069, https://doi.org/10.5194/acp530532005, 2005. a, b
Davies, S., Mann, G. W., Carslaw, K. S., Chipperfield, M. P., Remedios, J. J., Allen, G., Waterfall, A. M., Spang, R., and Toon, G. C.: Testing our understanding of Arctic denitrification using MIPASE satellite measurements in winter 2002/2003, Atmos. Chem. Phys., 6, 3149–3161, https://doi.org/10.5194/acp631492006, 2006. a
Dee, D. P., Uppala, S. M., Simmons, A. J., Berrisford, P., Poli, P., Kobayashi, S., Andrae, U., Balmaseda, M. A., Balsamo, G., Bauer, P., Bechtold, P., Beljaars, A. C. M., van de Berg, L., Bidlot, J., Bormann, N., Delsol, C., Dragani, R., Fuentes, M., Geer, A. J., Haimberger, L., Healy, S. B., Hersbach, H., Hólm, E. V., Isaksen, L., Kållberg, P., Köhler, M., Matricardi, M., McNally, A. P., MongeSanz, B. M., Morcrette, J.J., Park, B.K., Peubey, C., de Rosnay, P., Tavolato, C., Thépaut, J.N., and Vitart, F.: The ERAInterim reanalysis: configuration and performance of the data assimilation system, Q. J. Roy. Meteor. Soc., 137, 553–597, https://doi.org/10.1002/qj.828, 2011. a
Dessler, A. E.: The effect of deep, tropical convection on the tropical tropopause layer, J. Geophys. Res.Atmos., 107, ACH 61–ACH 65, https://doi.org/10.1029/2001JD000511, 2002. a
Deuflhard, P.: Newton Methods for Nonlinear Problems, Vol. 35 of Springer Series in Computational Mathematics, SpringerVerlag, Berlin Heidelberg, https://doi.org/10.1007/9783642238994, 2011. a
Dixon, J. C.: The Shock Absorber Handbook, 2nd Edn., John Wiley & Sons Ltd, Chichester, 2007. a, b, c, d, e
Durran, D. R. and Klemp, J. B.: On the Effects of Moisture on the BruntVäisälä Frequency, J. Atmos. Sci., 39, 2152–2158, https://doi.org/10.1175/15200469(1982)039<2152:OTEOMO>2.0.CO;2, 1982. a
Emanuel, K. A.: Atmospheric Convection, Oxford University Press, New York, USA, Oxford, UK, 1994. a, b, c
Ertel, H.: Ein neuer hydrodynamischer Wirbelsatz, Meteorol. Z., 59, 277–281, 1942. a, b
Fay, J. A.: Molecular Thermodynamics, Engineering Sciences, AddisonWesley, Reading, Massachusetts, 1965. a
Feistel, R.: A Gibbs function for seawater thermodynamics for −6 to 80 ^{∘}C and salinity up to 120 g kg^{−1}, DeepSea Res. Pt. I, 55, 1639–1671, https://doi.org/10.1016/j.dsr.2008.07.004, 2008. a
Frey, W.: Airborne in situ measurements of ice particles in the tropical tropopause layer, PhD thesis, Universität Mainz, Mainz, 2011. a
Furtenbacher, T., Horváth, M., Koller, D., Sólyom, P., Balogh, A., Balogh, I., and Császár, A. G.: MARVEL Analysis of the Measured HighResolution Rovibronic Spectra and Definitive IdealGas Thermochemistry of the ^{16}O_{2} Molecule, J. Phys. Chem. Ref. Data, 48, 023101, https://doi.org/10.1063/1.5083135, 2019. a
Gettelman, A., Hoor, P., Pan, L. L., Randel, W. J., Hegglin, M. I., and Birner, T.: The Extratropical Upper Troposphere and Lower Stratosphere, Rev. Geophys., 49, RG3003, https://doi.org/10.1029/2011RG000355, 2011. a, b
Hallam, A.: Alfred Wegener and the Hypothesis of Continental Drift, Sci. Am., 232, 88–97, 1975. a
Hauf, T. and Höller, H.: Entropy and Potential Temperature, J. Atmos. Sci., 44, 2887–2901, https://doi.org/10.1175/15200469(1987)044<2887:EAPT>2.0.CO;2, 1987. a
Holton, J. R.: An Introduction to Dynamic Meteorology, Vol. 88 of International Geophysics Series, Elsevier Academic Press, 4th Edn., 2004. a, b, c
Holton, J. R., Haynes, P. H., McIntyre, M. E., Douglass, A. R., Rood, R. B., and Pfister, L.: Stratospheretroposphere exchange, Rev. Geophys., 33, 403–439, https://doi.org/10.1029/95RG02097, 1995. a
Höpfner, M., Ungermann, J., Borrmann, S., Wagner, R., Spang, R., Riese, M., Stiller, G., Appel, O., Batenburg, A. M., Bucci, S., Cairo, F., Dragoneas, A., FriedlVallon, F., Hünig, A., Johansson, S., Krasauskas, L., Legras, B., Leisner, T., Mahnke, C., Möhler, O., Molleker, S., Müller, R., Neubert, T., Orphal, J., Preusse, P., Rex, M., Saathoff, H., Stroh, F., Weigel, R., and Wohltmann, I.: Ammonium nitrate particles formed in upper troposphere from ground ammonia sources during Asian monsoons, Nat. Geosci., 12, 608–612, https://doi.org/10.1038/s4156101903858, 2019. a
Hoskins, B. J., McIntyre, M. E., and Robertson, A. W.: On the use and significance of isentropic potential vorticity maps, Q. J. Roy. Meteor. Soc., 111, 877–946, 1985. a, b, c
Houghton, J.: The Physics of Atmospheres, Cambridge University Press, Cambridge, UK, 2002. a
Huang, K.: Statistical Mechanics, 2nd Edn., John Wiley & Sons, New York, 1987. a
Jakob, M.: Die spezifische Wärme der Luft im Bereich von 0 bis 200 at und von −80 bis 250^{∘}, Mitteilungen aus der PhysikalischTechnischen Reichsanstalt, Zeitschr. f. techn. Physik, 12, 460–468, 1923. a, b
Kondepudi, D. and Prigogine, I.: Modern Thermodynamics, John Wiley & Son, Chichester, 1998. a, b
Köppen, W. P.: Über Luftmischung und potentielle Temperatur, in: Anlehnung an die neueste Abhandlung von Herrn v. Helmholtz, oral presentation at Hamburg's local agency of the German Meteorological Society referenced in Wegener (1911), 1888. a, b
Kunz, O. and Wagner, W.: The GERG2008 WideRange Equation of State for Natural Gases and Other Mixtures: An Expansion of GERG2004, J. Chem. Eng. Data, 57, 3032–3091, https://doi.org/10.1021/je300655b, 2012. a
Kutzbach, G.: The Thermal Theory of Cyclones: A History of Meteorological Thought in the Nineteenth Century, Meteorological Monographs, American Meteorological Society, Boston, USA, 2016. a
Lary, D. J., Chipperfield, M. P., Pyle, J. A., Norton, W. A., and Riishøjgaard, L. P.: Threedimensional tracer initialization and general diagnostics using equivalent PV latitude–potentialtemperature coordinates, Q. J. Roy. Meteor. Soc., 121, 187–210, https://doi.org/10.1002/qj.49712152109, 1995. a
Lemmon, E. W., Jacobsen, R. T., Penoncello, S. G., and Friend, D. G.: Thermodynamic Properties of Air and Mixtures of Nitrogen, Argon, and Oxygen From 60 to 2000 K at Pressures to 2000 MPa, J. Phys. Chem. Ref. Data, 29, 331–385, https://doi.org/10.1063/1.1285884, 2000. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o
Lemmon, E. W., Bell, I. H., Huber, M. L., and McLinden, M. O.: NIST Standard Reference Database 23: Reference Fluid Thermodynamic and Transport PropertiesREFPROP, Version 10.0, National Institute of Standards and Technology, https://doi.org/10.18434/T4/1502528, 2018. a, b, c
Li, C., Ingersoll, A. P., and Oyafuso, F.: Moist Adiabats with Multiple Condensing Species: A New Theory with Application to GiantPlanet Atmospheres, J. Atmos. Sci., 75, 1063–1072, https://doi.org/10.1175/JASD170257.1, 2018. a
Liou, K. N.: An Introduction to Atmospheric Radiation, Vol. 84 of International Geophysics, 2nd Edn., Elsevier, San Diego, USA, 2002. a
Marquet, P.: Definition of a moist entropy potential temperature: application to FIREI data flights, Q. J. Roy. Meteor. Soc., 137, 768–791, https://doi.org/10.1002/qj.787, 2011. a
Marquet, P.: On the definition of a moistair potential vorticity, Q. J. Roy. Meteor. Soc., 140, 917–929, https://doi.org/10.1002/qj.2182, 2014. a
Marquet, P.: Le troisième principe de la thermodynamique ou une définition absolue de l'entropie – Partie 2: Définitions et applications en météorologie et en climat, https://doi.org/10.4267/2042/70559, an english translation is available at arXiv: https://arxiv.org/abs/1904.11699 (last access: 2 November 2020), 2019. a
Marquet, P. and Geleyn, J.F.: On a general definition of the squared BruntVäisälä frequency associated with the specific moist entropy potential temperature, Q. J. Roy. Meteor. Soc., 139, 85–100, https://doi.org/10.1002/qj.1957, 2013. a
Marquet, P. and Geleyn, J.F.: Formulations of moist thermodynamics for atmospheric modelling, chap. 22, in: Parameterization of Atmospheric Convection, edited by: Plant, R. S. and Yano, J.I., Imperial College Press, London, UK, 221–274, https://doi.org/10.1142/9781783266913_0026, 2015. a
McDougall, T. J., Jackett, D. R., Wright, D. G., and Feistel, R.: Accurate and Computationally Efficient Algorithms for Potential Temperature and Density of Seawater, J. Atmos. Ocean. Techn., 20, 730–741, https://doi.org/10.1175/15200426(2003)20<730:AACEAF>2.0.CO;2, 2003. a, b
Mote, P. W., Rosenlof, K. H., McIntyre, M. E., Carr, E. S., Gille, J. G., Holton, J. R., Kinnersley, J. S., Pumphrey, H. C., Russell III, J. M., and Waters, J. W.: An atmospheric tape recorder: The imprint of tropical tropopause temperatures on stratospheric water vapor, J. Geophys. Res., 101, 3989–4006, 1996. a
MüllerWodarg, I., Griffith, C. A., Lellouch, E., and Cravens, T. E.: Titan: Interior, Surface, Atmosphere, and Space Environment, Vol. 14 of Cambridge Planetary Science, Cambridge University Press, New York, USA, 2014. a
Murphy, D. M., Cziczo, D. J., Hudson, P. K., and Thomson, D. S.: Carbonaceous material in aerosol particles in the lower stratosphere and tropopause region, J. Geophys. Res.Atmos., 112, D04203, https://doi.org/10.1029/2006JD007297, 2007. a
Newell, D. B., Cabiati, F., Fischer, J., Fujii, K., Karshenboim, S. G., Margolis, H. S., de Mirandés, E., Mohr, P. J., Nez, F., Pachucki, K., Quinn, T. J., Taylor, B. N., Wang, M., Wood, B. M., and Zhang, Z.: The CODATA 2017 values of h, e, k, and N_{A} for the revision of the SI, Metrologia, 55, L13–L16, https://doi.org/10.1088/16817575/aa950a, 2018. a, b
Ooyama, K. V.: A Thermodynamic Foundation for Modeling the Moist Atmosphere, J. Atmos. Sci., 47, 2580–2593, https://doi.org/10.1175/15200469(1990)047<2580:ATFFMT>2.0.CO;2, 1990. a
Ooyama, K. V.: A Dynamic and Thermodynamic Foundation for Modeling the Moist Atmosphere with Parameterized Microphysics, J. Atmos. Sci., 58, 2073–2102, https://doi.org/10.1175/15200469(2001)058<2073:ADATFF>2.0.CO;2, 2001. a, b
Poisson, S. D.: Traité de mécanique, Bachelier, Paris, 2nd Edn., available at: https://gallica.bnf.fr/ark:/12148/bpt6k9605452x.texteImage (last access: 2 November 2020), 1833. a
Pommrich, R., Müller, R., Grooß, J.U., Konopka, P., Ploeger, F., Vogel, B., Tao, M., Hoppe, C. M., Günther, G., Spelten, N., Hoffmann, L., Pumphrey, H.C., Viciani, S., D'Amato, F., Volk, C. M., Hoor, P., Schlager, H., and Riese, M.: Tropical troposphere to stratosphere transport of carbon monoxide and longlived trace species in the Chemical Lagrangian Model of the Stratosphere (CLaMS), Geosci. Model Dev., 7, 2895–2916, https://doi.org/10.5194/gmd728952014, 2014. a
Protter, M. H. and Morrey, C. B.: Intermediate Calculus, Undergraduate Texts in Mathematics, 2nd Edn., Springer, Berlin, Heidelberg, https://doi.org/10.1007/9781461210863, 1985. a
Pruppacher, H. R. and Klett, J. D.: Microphysics of Clouds and Precipitation, Vol. 18 of Atmospheric and Oceanographic Sciences Library, Kluwer Academic Publishers, Dordrecht, 2010. a
ReinkeKunze, C.: Alfred Wegener: Polarforscher und Entdecker der Kontinentaldrift, Lebensgeschichten aus der Wissenschaft, Birkhäuser Basel, 2013. a
Richardson, M. I., Toigo, A. D., and Newman, C. E.: PlanetWRF: A general purpose, local to global numerical model for planetary atmospheric and climate dynamics, J. Geophys. Res.Planet., 112, E09001, https://doi.org/10.1029/2006JE002825, 2007. a
Roebuck, J. R.: The JouleThomson Effect in Air, P. Am. Acad. Arts Sci., 60, 537–596,1925. a
Roebuck, J. R.: The JouleThomson Effect in Air. Second Paper, P. Am. Acad. Arts Sci., 64, 287–334, 1930. a
Roedel, W. and Wagner, T.: Physik unserer Umwelt: Die Atmosphäre, Springer Berlin Heidelberg, https://doi.org/10.1007/9783642157295, 2011. a
Scheel, K. and Heuse, W.: Die spezifische Wärme der Luft bei Zimmertemperatur und bei tiefen Temperaturen, Ann. Phys., 342, 79–95, https://doi.org/10.1002/andp.19113420106, 1912. a
Schmidt, R. and Wagner, W.: A new form of the equation of state for pure substances and its application to oxygen, Fluid Phase Equilibr., 19, 175–200, https://doi.org/10.1016/03783812(85)870163, 1985. a, b, c, d
Schubert, W., Ruprecht, E., Hertenstein, R., Ferreira, R. N., Taft, R., Rozoff, C., Ciesielski, P., and Kuo, H.C.: English translations of twentyone of Ertel's papers on geophysical fluid dynamics, Meteorol. Z., 13, 527–576, https://doi.org/10.1127/09412948/2004/00130527, 2004. a, b
Schumann, U.: Atmospheric Physics; Background – Methods – Trends, Research Topics in Aerospace, SpringerVerlag, Berlin Heidelberg, https://doi.org/10.1007/9783642301834, 2012. a
Seinfeld, J. H. and Pandis, S. N.: Atmospheric Chemistry and Physics: From Air Pollution to Climate Change, Wiley, Hoboken, New Jersey, USA, 2006. a
Skamarock, W. C. and Klemp, J. B.: A timesplit nonhydrostatic atmospheric model for weather research and forecasting applications, J. Comput. Phys., 227, 3465–3485, https://doi.org/10.1016/j.jcp.2007.01.037, 2008. a, b
Skamarock, W. C., Klemp, J. B., Dudhia, J., Gill, D. O., Barker, D. M., Wang, W., and Powers, J. G.: A Description of the Advanced Research WRF Version 2 (No. NCAR/TN468+STR), University Corporation for Atmospheric Research, https://doi.org/10.5065/D6DZ069T, 2005. a, b
Span, R., Lemmon, E. W., Jacobsen, R. T., Wagner, W., and Yokozeki, A.: A Reference Equation of State for the Thermodynamic Properties of Nitrogen for Temperatures from 63.151 to 1000 K and Pressures to 2200 MPa, J. Phys. Chem. Ref. Data, 29, 1361–1433, https://doi.org/10.1063/1.1349047, 2000. a, b, c, d, e
Spang, R., Remedios, J. J., Kramer, L. J., Poole, L. R., Fromm, M. D., Müller, M., Baumgarten, G., and Konopka, P.: Polar stratospheric cloud observations by MIPAS on ENVISAT: detection method, validation and analysis of the northern hemisphere winter 2002/2003, Atmos. Chem. Phys., 5, 679–692, https://doi.org/10.5194/acp56792005, 2005. a
Spichtinger, P.: Shallow cirrus convection – a source for ice supersaturation, Tellus A, 66, 19937, https://doi.org/10.3402/tellusa.v66.19937, 2014. a
Stamnes, K., Thomas, G. E., and Stamnes, J. J.: Radiative Transfer in the Atmosphere and Ocean, Cambridge University Press, https://doi.org/10.1017/9781316148549, 2017. a
Stohl, A., Bonasoni, P., Cristofanelli, P., Collins, W., Feichter, J., Frank, A., Forster, C., Gerasopoulos, E., Gäggeler, H., James, P., Kentarchos, T., KrompKolb, H., Krüger, B., Land, C., Meloen, J., Papayannis, A., Priller, A., Seibert, P., Sprenger, M., Roelofs, G. J., Scheel, H. E., Schnabel, C., Siegmund, P., Tobler, L., Trickl, T., Wernli, H., Wirth, V., Zanis, P., and Zerefos, C.: Stratospheretroposphere exchange: A review, and what we have learned from STACCATO, J. Geophys. Res.Atmos., 108, 8516, https://doi.org/10.1029/2002JD002490, 2003. a
Swinbank, R. and Ortland, D. A.: Compilation of wind data for the (UARS) Reference Atmosphere Project, J. Geophys. Res., 108, 4615, https://doi.org/10.1029/2002JD003135, 2003. a
Tao, M., Konopka, P., Ploeger, F., Yan, X., Wright, J. S., Diallo, M., Fueglistaler, S., and Riese, M.: Multitimescale variations in modeled stratospheric water vapor derived from three modern reanalysis products, Atmos. Chem. Phys., 19, 6509–6534, https://doi.org/10.5194/acp1965092019, 2019. a
Tegeler, C., Span, R., and Wagner, W.: A New Equation of State for Argon Covering the Fluid Region for Temperatures from the Melting Line to 700 K at Pressures up to 1000 MPa, J. Phys. Chem. Ref. Data, 28, 779–850, https://doi.org/10.1063/1.556037, 1999. a, b
Thomson, W.: On the convective equilibrium of temperature in the atmosphere, Manch. Lit. Philos. Soc., 2, 170–176, 1862. a, b
Tiesinga, E., Mohr, P. J., Newell, D. B., and Taylor, B. N.: The 2018 CODATA Recommended Values of the Fundamental Physical Constants, Web Version 8.1, Database developed by: Baker, J., Douma, M., and Kotochigova, S., available at: http://physics.nist.gov/constants (last access: 2 November 2020), 2020. a, b, c, d
Tiwary, A. and Williams, I.: Air Pollution: Measurement, Modelling and Mitigation, 4th Edn., CRC Press, Boca Raton, FL, USA, 2019. a
Tripoli, G. J. and Cotton, W. R.: The Use of lceLiquid Water Potential Temperature as a Thermodynamic Variable In Deep Atmospheric Models, Mon. Weather Rev., 109, 1094–1102, https://doi.org/10.1175/15200493(1981)109<1094:TUOLLW>2.0.CO;2, 1981. a
Tsilingiris, P. T.: Thermophysical and transport properties of humid air at temperature range between 0 and 100 ^{∘}C, Energ. Convers. Manage., 49, 1098–1110, https://doi.org/10.1016/j.enconman.2007.09.015, 2008. a, b
United States Committee on Extension to the Standard Atmosphere: U.S. standard atmosphere, 1976, National Oceanic and Atmospheric Administration, U.S. Government Printing Office, Washington D. C., available at: https://ntrs.nasa.gov/archive/nasa/casi.ntrs.nasa.gov/19770009539.pdf (last access: 2 November 2020), 1976. a, b
Vasserman, A. A., Kazavchinskii, Y. Z., and Rabinovich, V. A.: Thermophysical Properties of Air and Air Components, Academy of Sciences of the USSR, Israel Program for Scientific Translation Ltd., Jerusalem 1971, 1966. a, b, c
Vernier, J.P., Fairlie, T. D., Deshler, T., Ratnam, M. V., Gadhavi, H., Kumar, B. S., Natarajan, M., Pandit, A. K., Raj, S. T. A., Kumar, A. H., Jayaraman, A., Singh, A. K., Rastogi, N., Sinha, P. R., Kumar, S., Tiwari, S., Wegner, T., Baker, N., Vignelles, D., Stenchikov, G., Shevchenko, I., Smith, J., Bedka, K., Kesarkar, A., Singh, V., Bhate, J., Ravikiran, V., Rao, M. D., Ravindrababu, S., Patel, A., Vernier, H., Wienhold, F. G., Liu, H., Knepp, T. N., Thomason, L., Crawford, J., Ziemba, L., Moore, J., Crumeyrolle, S., Williamson, M., Berthet, G., Jégou, F., and Renard, J.B.: BATAL: The Balloon Measurement Campaigns of the Asian Tropopause Aerosol Layer, B. Am. Meteorol. Soc., 99, 955–973, https://doi.org/10.1175/bamsd170014.1, 2018. a
Voigt, C., Schumann, U., Minikin, A., Abdelmonem, A., Afchine, A., Borrmann, S., Boettcher, M., Bucuchholz, B., Bugliaro, L., Costa, A., Curtius, J., Dollner, M., Dörnbrack, A., Dreiling, V., Ebert, V., Ehrlich, A., Fix, A., Forster, L., Frank, F., Futterer, D., Giez, A., Graf, K., Grooss, J. U., Gross, S., Heimerl, K., Heinold, B., Huneke, T., Järvinen, E., Jurkat, T., Kaufmann, S., Kenntner, M., Klingebiel, M., Klimach, T., Kohl, R., Krämer, M., Krisna, T. C., Luebke, A., Mayer, B., Mertes, S., Molleker, S., Petzold, A., Pfeilsticker, K., Port, M., Rapp, M., Reutter, P., Rolf, C., Rose, D., Sauer, D., Schafer, A., Schlage, R., Schnaiter, M., Schneider, J., Spelten, N., Spichtinger, P., Stock, P., Walser, A., Weigel, R., Weinzierl, B., Wendisch, M., Werner, F., Wernli, H., Wirth, M., Zahn, A., Ziereis, H., and Zoger, M.: MlCirrus the Airborne Experiment on Natural Cirrus and Contrail Cirrus with the HighAltitude LongRange Research Aircraft Halo, B. Am. Meteorol. Soc., 98, 271–288, https://doi.org/10.1175/BamsD1500213.1, 2017. a
von Bezold, W.: Zur Thermodynamik der Atmosphaere, 1189–1206, Verlag der königlichen Akademie der Wissnschaften, Berlin, 1888. a, b, c, d
von Helmholtz, H.: Über atmosphaerische Bewegungen, 647–663, Verlag der königlichen Akademie der Wissnschaften, Berlin, 1888. a, b
Wagner, W. and de Reuck, K. M.: Oxygen, International Thermodynamic Tables of the Fluid State, Vol. 9 of IUPAC Thermodynamic Tables Project, Blackwell Science, Oxford, UK, 1987. a
Wallace, J. M. and Hobbs, P. V.: Atmospheric Science, Vol. 92 of International Geophysics Series, 2nd Edn., Academic Press, San Diego, USA, 2006. a, b
Wegener, A.: Thermodynamik der Atmosphäre, J. A. Barth, Leipzig, Germany, 1911. a
Wegener, A. and Wegener, K.: Vorlesung über Physik der Atmosphäre, J. A. Barth, Leipzig, 1935. a
Weigel, R., Volk, C. M., Kandler, K., Hösen, E., Günther, G., Vogel, B., Grooß, J.U., Khaykin, S., Belyaev, G. V., and Borrmann, S.: Enhancements of the refractory submicron aerosol fraction in the Arctic polar vortex: feature or exception?, Atmos. Chem. Phys., 14, 12319–12342, https://doi.org/10.5194/acp14123192014, 2014. a
Weigel, R., Spichtinger, P., Mahnke, C., Klingebiel, M., Afchine, A., Petzold, A., Krämer, M., Costa, A., Molleker, S., Reutter, P., Szakáll, M., Port, M., Grulich, L., Jurkat, T., Minikin, A., and Borrmann, S.: Thermodynamic correction of particle concentrations measured by underwing probes on fastflying aircraft, Atmos. Meas. Tech., 9, 5135–5162, https://doi.org/10.5194/amt951352016, 2016. a
Wendisch, M. and Brenguier, J.L.: Airborne Measurements for Environmental Research: Methods and Instruments, Wiley‐VCH, Weinheim, https://doi.org/10.1002/9783527653218, 2013. a, b
Wendisch, M., Pöschl, U., Andreae, M. O., Machado, L. A. T., Albrecht, R., Schlager, H., Rosenfeld, D., Martin, S. T., Abdelmonem, A., Afchine, A., Araùjo, A., Artaxo, P., Aufmhoff, H., Barbosa, H. M. J., Borrmann, S., Braga, R., Buchholz, B., Cecchini, M. A., Costa, A., Curtius, J., Dollner, M., Dorf, M., Dreiling, V., Ebert, V., Ehrlich, A., Ewald, F., Fisch, G., Fix, A., Frank, F., Fütterer, D., Heckl, C., Heidelberg, F., Hüneke, T., Jäkel, E., Järvinen, E., Jurkat, T., Kanter, S., Kästner, U., Kenntner, M., Kesselmeier, J., Klimach, T., Knecht, M., Kohl, R., Kölling, T., Krämer, M., Krüger, M., Krisna, T. C., Lavric, J. V., Longo, K., Mahnke, C., Manzi, A. O., Mayer, B., Mertes, S., Minikin, A., Molleker, S., Münch, S., Nillius, B., Pfeilsticker, K., Pöhlker, C., Roiger, A., Rose, D., Rosenow, D., Sauer, D., Schnaiter, M., Schneider, J., Schulz, C., de Souza, R. A. F., Spanu, A., Stock, P., Vila, D., Voigt, C., Walser, A., Walter, D., Weigel, R., Weinzierl, B., Werner, F., Yamasoe, M. A., Ziereis, H., Zinner, T., and Zöger, M.: The ACRIDICONCHUVA campaign: Studying tropical deep convective clouds and precipitation over Amazonia using the new German research aircraft HALO, B. Am. Meteorol. Soc., 97, 1885–1908, https://doi.org/10.1175/BAMSD1400255.1, 2016. a
Wilson, J. C., Loewenstein, M., Fahey, D. W., Gary, B., Smith, S. D., Kelly, K. K., Ferry, G. V., and Chan, K. R.: Observations of condensation nuclei in the Airborne Antarctic Ozone Experiment: Implications for new particle formation and polar stratospheric cloud formation, J. Geophys. Res.Atmos., 94, 16437–16448, https://doi.org/10.1029/JD094iD14p16437, 1989. a
Witkowski, A. W.: I. Thermodynamic properties of air, The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, 42, 1–37, https://doi.org/10.1080/14786449608620887, 1896. a, b
WMO: Meteorology – A threedimensional science, WMO Bull., 6, 134–138, 1957. a
WMO: International Meteorological Tables, WMONo.188.TP97, Secretariat of the World Meteorological Organization, Geneva, Switzerland, available at: https://library.wmo.int/doc_num.php?explnum_id=7997 (last access: 11 December 2020), 1966. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o, p, q, r, s, t, u, v, w
Zängl, G., Reinert, D., Rípodas, P., and Baldauf, M.: The ICON (ICOsahedral Nonhydrostatic) modelling framework of DWD and MPIM: Description of the nonhydrostatic dynamical core, Q. J. Roy. Meteor. Soc., 141, 563–579, https://doi.org/10.1002/qj.2378, 2015. a
Zdunkowski, W. and Bott, A.: Dynamics of the Atmosphere: A Course in Theoretical Meteorology, Cambridge University Press, Cambridge, UK, 2003. a
Wegener mentioned a talk given by Köppen in a footnote on p. 111. In the publication year (1911) of Wegener's book, Köppen's daughter Else got engaged to Alfred Wegener (ReinkeKunze, 2013) and they married in the year 1913 (Hallam, 1975).
Cf. Bauer (1908), where, for the first time, the potential temperature and the entropy are set in a relationship.
The value of R is now defined exactly, cf. Tiesinga et al. (2020) and Newell et al. (2018), and is used in Eq. (6).
 Abstract
 Introduction
 Derivation of the potential temperature for an ideal gas
 Examining the assumption of constant c_{p} for dry air
 Accounting for the temperature dependence of air's specific heat capacity
 The ${\mathit{\theta}}_{{c}_{p}\left(T\right)}$ from the temperaturedependent specific heat capacity of air
 The potential temperature for air as a real gas
 Implications on the use of the potential temperature
 Summary and conclusions
 Appendix A: Derivation of the specific heat capacity from thermodynamics
 Appendix B: Sensitivity of the conventional definition of θ to perturbations of c_{p}
 Appendix C: Approximate computation of the reference potential temperature
 Appendix D: Approximation error for the reference potential temperature
 Appendix E: The derivative of the reference potential temperature
 Data availability
 Author contributions
 Competing interests
 Acknowledgements
 Financial support
 Review statement
 References
 Abstract
 Introduction
 Derivation of the potential temperature for an ideal gas
 Examining the assumption of constant c_{p} for dry air
 Accounting for the temperature dependence of air's specific heat capacity
 The ${\mathit{\theta}}_{{c}_{p}\left(T\right)}$ from the temperaturedependent specific heat capacity of air
 The potential temperature for air as a real gas
 Implications on the use of the potential temperature
 Summary and conclusions
 Appendix A: Derivation of the specific heat capacity from thermodynamics
 Appendix B: Sensitivity of the conventional definition of θ to perturbations of c_{p}
 Appendix C: Approximate computation of the reference potential temperature
 Appendix D: Approximation error for the reference potential temperature
 Appendix E: The derivative of the reference potential temperature
 Data availability
 Author contributions
 Competing interests
 Acknowledgements
 Financial support
 Review statement
 References