Journal cover Journal topic
Atmospheric Chemistry and Physics An interactive open-access journal of the European Geosciences Union
Journal topic
ACP | Articles | Volume 20, issue 24
Atmos. Chem. Phys., 20, 15585–15616, 2020
© Author(s) 2020. This work is distributed under
the Creative Commons Attribution 4.0 License.
Atmos. Chem. Phys., 20, 15585–15616, 2020
© Author(s) 2020. This work is distributed under
the Creative Commons Attribution 4.0 License.

Research article 15 Dec 2020

Research article | 15 Dec 2020

Reappraising the appropriate calculation of a common meteorological quantity: potential temperature

Reappraising the appropriate calculation of a common meteorological quantity: potential temperature
Manuel Baumgartner1,2, Ralf Weigel2, Allan H. Harvey3, Felix Plöger4,5, Ulrich Achatz6, and Peter Spichtinger2 Manuel Baumgartner et al.
  • 1Zentrum für Datenverarbeitung, Johannes Gutenberg University Mainz, Mainz, Germany
  • 2Institute for Atmospheric Physics, Johannes Gutenberg University Mainz, Mainz, Germany
  • 3Applied Chemicals and Materials Division, National Institute of Standards and Technology, Boulder, CO, USA
  • 4Forschungszentrum Jülich GmbH, Institute of Energy and Climate Research (IEK-7), Jülich, Germany
  • 5Institute for Atmospheric and Environmental Research, University of Wuppertal, Wuppertal, Germany
  • 6Institut für Atmosphäre und Umwelt, Goethe-Universität Frankfurt, Frankfurt am Main, Germany

Correspondence: Manuel Baumgartner (

Back to toptop

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 temperature-dependent 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.

1 Introduction
Back to toptop

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 Kutzbach2016). 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” (Thomson1862)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 eddy-driven 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 Helmholtz1888). 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 (p0) 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:

(1) θ = T p 0 p γ - 1 γ ,

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 p0 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 Bezold1888).

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 dry-adiabatic 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 one-to-one 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 dry-adiabatic 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

(2) N 2 = g θ θ z ,

where g=9.81ms-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 θz. In the research field of dynamical meteorology, the potential vorticity (PV) is often used (Ertel1942; 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 three-dimensional 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 large-scale 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 dry-air potential temperature have been developed (Hauf and Höller1987; Emanuel1994; Marquet2011; Marquet and Geleyn2015) 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. Emanuel1994). 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 so-called pseudo-adiabatic 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. Spichtinger2014). 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 in-cloud conditions. Therefore, above clouds and further aloft, e.g. within the stratosphere, the conventional dry-air 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 Klemp2008; 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 dry-air situations and cloud conditions.

In any case, the use of the potential temperature requires the following preconditions to be fulfilled:

  1. θ 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

  2. θ 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 ideal-gas assumption, a re-assessment of the fundamental atmospheric quantity θ is suggested, which is based on the state of knowledge of air's thermodynamic properties, and this re-assessed θ 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 p0=1000hPa, as applied to the Earth's atmosphere, is frequently used as a reference pressure for the atmosphere of other planets (Catling2015, 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 Klemp2008). However, the temperature dependency of the isobaric heat capacity cp is not generally negligible, especially when taking “deep atmospheres, such as on Venus” (Catling2015, 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üller-Wodarg et al.2014).

Moreover, the potential temperature is a frequently used quantity in oceanography (e.g. McDougall et al.2003; Feistel2008), 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 re-assessment 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 atmosphere-related purposes. These studies from other disciplines motivate the need for a re-assessment 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 cp is recalled in Sect. 2. In Sect. 3 the assumption of a constant cp is discussed together with a synopsis of various cp values as provided in the literature. The temperature dependency of cp 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 temperature-dependent specific heat capacity, while Sect. 6 focuses on the influence of real-gas effects on the resulting potential temperature. Section 7 presents some implications of the use of θref, and concluding remarks are given in Sect. 8.

2 Derivation of the potential temperature for an ideal gas
Back to toptop

The Gibbs equation (see e.g. Kondepudi and Prigogine1998) is a general thermodynamic relation to describe the state of a system with m components and reads as

(3) T d S = d H - V d p - k = 1 m μ k d M k ,

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 (m3), μk the chemical potential of component k in joules per kilogram (J kg−1), Mk 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 dMk=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 Ma of dry air; note that lowercase letters indicate specific variables – e.g. h=H/Ma) leads to

(4) T d s = d h - V M a d p d s = 1 T d h - V M a T d p .

Furthermore, approximating dry air as an ideal gas leads to the following simplifications.

  • The ideal-gas law

    (5) p V = M a R a T

    can be applied with the specific gas constant Ra of dry air, which is

    (6) R a = R M mol,a = 8.31446261815324 J mol - 1 K - 1 0.0289586 kg mol - 1 ± 0.0000002 kg mol - 1 287.11350 J kg - 1 K - 1 , 287.11748 J kg - 1 K - 1 ,

    with the molar gas constant R in Jmol-1K-1 (Tiesinga et al.2020; Newell et al.2018) and Mmol,a the molar mass of dry air (Lemmon et al.2000), composed of nitrogen N2, oxygen O2, and argon Ar.

  • The specific enthalpy is given by

    (7) d h = c p d T ,

    with cp 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

(8) d s = c p T d T - R a d p p .

For isentropic changes of state, i.e. ds=0, Eq. (8) reduces to

(9) c p T d T = R a d p p .

Note that the assumption of dry air being an ideal gas does not imply that in Eq. (9) the specific heat capacity cp is constant. While statistical mechanics excludes any pressure dependence in the ideal-gas heat capacity, the general derivation (cf. Appendix A) permits a temperature dependence of cp. However, usually the temperature dependence is neglected in atmospheric physics, and, instead, cp is assumed as constant (see e.g. Ambaum2010, p. 48/49, where vibrational modes of the air molecules are neglected). Immediately below and in Sect. 3, the treatment of cp as a temperature-independent constant is discussed. The introduction of the temperature dependence then follows in Sect. 4.

Treating cp as a constant, rearrangement of Eq. (9) leads to

(10) d T T = R a c p d p p .

Integration of Eq. (10) over the range from ground-level pressure and temperature (p0, T0) to the pressure and temperature at a specific height (p, T) yields

(11) ln T T 0 = T 0 T d T T = R a c p p 0 p d p p = R a c p ln p p 0 ,

and, after another straightforward conversion, one arrives at

(12) ln T 0 T = R a c p ln p 0 p .

With the definition θcp=T0, Eq. (12) is transformed into the commonly used expression for determining the potential temperature

(13) θ c p = T p 0 p R a c p ,

for which the ground-level pressure p0 is arbitrary but usually set to p0=1000hPa. This choice coincides with the definition of the World Meteorological Organisation (WMO1966) and the standard-state pressure (Tiesinga et al.2020) but should not be confused with the standard atmosphere 101 325 Pa (Tiesinga et al.2020). In the following, θcp denotes the potential temperature based on a constant cp, and, when a specific value of cp is applied, the subscript cp in the potential temperature's notation is replaced by the corresponding cp value.

3 Examining the assumption of constant cp for dry air
Back to toptop

The general theory of thermodynamics, assuming dry air as an ideal gas, gives the expression

(14) c p = 1 + f 2 R a

for the constant specific heat capacity, which is based on the results of statistical mechanics and the equipartition theorem (e.g. Huang1987). In Eq. (14), the parameter f=ftrans+frot+fvib 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 ftrans, rotation frot, and vibration fvib. Assuming further that dry air exclusively consists of the linear molecules N2 and O2 (implying ftrans=3 and frot=2, while the contribution of Ar remains disregarded) and additionally neglecting the vibrational degrees of freedom (fvib=0), the general relation Eq. (14) reduces to

(15) c p = 1 + 3 + 2 2 R a = 7 2 R a .

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 cp may be appropriate. According to the WMO (1966), the recommended value for cp of dry air is 1005Jkg-1K-1, and, furthermore (ibid.), it is defined that γ=cpcv=75=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 cp, a single consistently used value of cp was not found. Instead, the specified values of cp 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 cp that ranges from 994 to 1011Jkg-1K-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 1010Jkg-1K-1 is considered.

Wegener and Wegener (1935, converted from units other than SI)Roedel and Wagner (2011, p. 66)(WMO1966)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)(WMO1966)

Table 1Synopsis of temperature-independent constant values given mainly in textbooks for the specific heat capacity cp of dry air from various sources (non-exhaustive). Note that the WMO (1966) indicates a minimum and maximum “range of actual values” together with their recommended value cp=1005Jkg-1K-1.

Download Print Version | Download XLSX

Figure 1Vertical profiles of (a) atmospheric pressure and (b) temperature as functions of height, corresponding to the US Standard Atmosphere.


Figure 2Computed vertical course of the potential temperature θcp based on the two extremes of constant values for the specific heat capacity cp provided in the literature including the historical extreme values (a; cf. Table 1), and (b) the absolute differences Δθcp=θ994θ1011 and Δθcp=θ1000θ1010 between the two resulting curves of θcp. The absolute difference Δθcp=θ1003.5θ1006.5 is also shown (green curve), corresponding to a more realistic interval of cp values.


These different values of constant cp scatter within a small range (below ±1.1 %) around the WMO's recommendation 1005Jkg-1K-1, which may seem negligible if cp 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 θcp, cf. Eq. (13), the specific heat capacity cp does not contribute linearly but rather as the denominator in the exponent. Thus, the variety of different cp values, although scattering within a small range, impacts the resulting θcp significantly. To illustrate this impact, a computation of θcp 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 Atmosphere1976). From the list of the different cp in Table 1, the extreme values were selected in order to initially illustrate the sensitivity of the resulting θcp to variations in cp in the range of ∼1%, as seen in the literature. In Fig. 2a, the individual profiles of θcp are shown for the extremes of the historic cp values (Table 1), while Fig. 2b illustrates the absolute differences Δθcp=θ994θ1011 (red curve), Δθcp=θ1000θ1010 (blue curve), and Δθcp=θ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 cp values in the literature (Table 1). At an altitude of 8.5 km, the difference Δθcp already exceeds 1 K (blue curve). The values of Δθcp 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 Δθcp reaches 43 K. The green curve corresponds to the more realistic cp interval 1005Jkg-1K-1±1.5Jkg-1K-1 as recommended by the WMO; the difference reaches approximately 13 K at the stratopause.

Figure 2 illustrates the possible spread of θcp based on a range of cp values from different literature references; hence, if one uses a different value for cp from the literature than that defined by WMO (1966), the difference θ1005-θcp might be significant. Since the cp values provided by some literature references are close to the value cp=1005Jkg-1K-1 recommended by the WMO (1966), the subsequent comparisons will be made to θ1005. The θcp values based on cp values other than 1005Jkg-1K-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 cp values (cf. Table 1), they nevertheless illustrate the sensitive response of θcp to even small variations in cp, 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 small-scale 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 θcp. 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” (Ooyama2001); 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 (cp(T); cf. Sect. 4) may be easily adopted. However, the crucial drawback of the entropy-based 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 cp vary, but also the values of the gas constant Ra vary slightly due to different historical approximations for the molar gas constant4 R and for the composition of dry air. The variation of values for Ra is typically only on the order of 0.1Jkg-1K-1, whereas the variability in cp is on the order of a few Jkg-1K-1 (cf. Table 1). Therefore, within the exponent of the expression (13) for θcp, the variability of cp has by far a stronger impact on the resulting θcp value than the variability of Ra.

However, accepting for a moment the WMO's definition (15) of cp (WMO1966), the variability of air's cp should naturally be constrained to certain limits. With the specific gas constant Ra=287.05Jkg-1K-1 (WMO1966), the WMO's definition leads to cp=1004.675Jkg-1K-1. In contrast, taking into account the uncertainty introduced in Ra by the molar mass of dry air, cf. Eq. (6), the resulting range for air's specific heat capacity is 1004.897Jkg-1K-1cp1004.912Jkg-1K-1. It may be surmised that the rounded value cp=1005Jkg-1K-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.

4 Accounting for the temperature dependence of air's specific heat capacity
Back to toptop

Next, while retaining the ideal-gas assumption, we consider the dependence of air's cp on temperature, mainly over the atmospherically relevant range (180 to 300 K). The temperature dependence of cp 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 cp 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 cp values initially decline, then pass a minimum, and subsequently increase again at lower temperatures (T<170K). The description of refined experiments and ascertainable data of air's cp(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 cp(T) of dry air (“air containing neither carbon-dioxide nor steam”, Awano1936), 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 3Variety of suggested values for the specific heat capacity of air. Ranges of constant values for cp (including the historical) together with the recommended value by the WMO (1966) are displayed as given in Table 1 (dashed lines). The parameterisations of air's cp0(T), assuming dry air as an ideal gas, accounting for its temperature dependence by Lemmon et al. (2000, solid magenta curve) and by Dixon (2007, solid cyan curve) are displayed. Discrete measurement and literature data at about 1000 hPa (i.e. as often specified, at “about one atmosphere”) are indicated by dots. In addition, the studies by Awano (1936) and Vasserman et al. (1966) provide data at other atmospheric pressures, as indicated by squares, diamonds, and triangles.


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 cp(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 cp(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 cp(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 cp(T) diverges significantly. Possible reasons for this include the following:

  • the measurements of cp(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 cp(T) data, (ii) a constant cp (e.g. 1005Jkg-1K-1; WMO1966), and (iii) the parameterised cp(T) is found only for a temperature interval ranging from 270 to 300 K. For air temperatures below 270 K, the constant value cp=1005Jkg-1K-1 is only comparable with the values from Vasserman et al. (1966) but fails to coincide with other parameterised or experimentally determined values of cp(T).

4.1 The temperature dependence of the ideal-gas specific heat capacity

As already indicated by the data depicted in Fig. 3, the specific heat capacity cp depends on the gas temperature. With regard to measured values, the lack of constancy may be due to real-gas effects or to a dependence of the ideal-gas heat capacity on temperature. In this section, we focus on the latter effect, denoting the ideal-gas isobaric specific heat capacity by cp0(T), where the superscript 0 indicates the underlying ideal-gas assumption. For an individual gas, there is always a contribution from the three translational degrees of freedom, cp,trans0=52Ri, where Ri 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

(16) c p , rot 0 = R i , for linear (e.g. diatomic) molecules , 3 2 R i , for nonlinear molecules .

As mentioned previously, at finite temperatures molecules also have contributions to cp0(T) from intramolecular vibrations (and, at high temperatures, excited electronic states). To arrive at a temperature-dependent parameterisation for the ideal-gas 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 three-component mixture: the diatomic gases nitrogen (N2) and oxygen (O2) and the monatomic gas argon (Ar).

To determine the contribution of N2 to cp0(T), both Bücker et al. (2002) and Lemmon et al. (2000) use the ideal-gas 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 Δcp0 of less than 0.02 %.

For the contribution of O2, 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 Reuck1987), after refitting it to more recently obtained data, thereby achieving an overall uncertainty Δcp0 of less than ±0.015 % for O2 (Bücker et al.2002). However, the difference in the resulting specific heat capacity contribution by O2 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 cp0 for O2 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 cp0=52RAr.

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 Δcp0±0.015% in the temperature range of 200KT3300K. 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 cp0(T)

Besides a comprehensive survey of the available experimental data for the specific heat capacity of air, Lemmon et al. (2000) also provide state-of-the-art 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:

  1. an empirical model-based equation of state for standard (dry) air considered as a pseudo-pure fluid and

  2. assembly of a mixture model from equations of state for each pure fluid.

Each approach allows calculating the thermodynamic properties, e.g. cp, of gas mixtures such as dry air, and both are real-gas models with the ideal-gas behaviour as a boundary condition. The major difference between the models is that the first approach considers air as a pseudo-pure fluid while the second, more rigorous approach treats air as a mixture composed of N2, O2, 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 Atmosphere1976), 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 (CO2) 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 pseudo-pure fluid model and the mixture model are implemented in NIST's REFPROP database (cf., last access: 2 November 2020) for various physical properties of fluids over a wide range of temperatures and pressures.

Both the pseudo-pure fluid model and the mixture model of Lemmon et al. (2000) use the same expression for the ideal-gas heat capacity, which is rigorously given as a sum of the pure-component contributions:

(17) C p 0 ( T ) R = x N 2 C p 0 ( T ) R N 2 + x Ar C p 0 ( T ) R Ar + x O 2 C p 0 ( T ) R O 2 ,

where xi denotes the molar fraction of species i, and Cp0 and the molar gas constant R are given in units of Jmol-1K-1.

Like Bücker et al. (2002), Lemmon et al. (2000) use the expression of Span et al. (2000) for the contribution of N2 to the heat capacity and adopt Cp0=52R for Ar. Together with the contribution by O2 according to the formulation by Schmidt and Wagner (1985), the expression provided by Lemmon et al. (2000, Eq. 18 therein) for the ideal-gas heat capacity of dry air is

(18) C p 0 ( T ) R = N 1 + N 2 T + N 3 T 2 + N 4 T 3 + N 5 T - 3 2 + N 6 N 9 2 T 2 exp N 9 T exp N 9 T - 1 2 + N 7 N 10 2 T 2 exp N 10 T exp N 10 T - 1 2 + 2 N 8 3 N 11 2 T 2 exp - N 11 T 2 3 exp - N 11 T + 1 2 ,

with the scalar coefficients Ni for dry air (ibid.),

(19) N 1 = 3.490888032 , N 2 = 2.395525583 × 10 - 6 , N 3 = 7.172111248 × 10 - 9 , N 4 = - 3.115413101 × 10 - 13 , N 5 = 0.223806688 , N 6 = 0.791309509 , N 7 = 0.212236768 , N 8 = 0.197938904 , N 9 = 3364.011 , N 10 = 2242.45 , N 11 = 11580.4 ,

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, Cp0(T)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 cp0(T) from an engineer's perspective

The parameterisation from Dixon (2007)

(20) c p ( T ) = 1002.5 + 275 × 10 - 6 T - 200 2

for 200KT450K 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 ideal-gas limit cp0(T) (Lemmon et al.2000) and Dixon's parameterisation substantially increases. This is most likely due to the chosen type of polynomial approximation (Dixon2007), which increasingly departs from the reference cp0(T) 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 cp(T) is significantly larger (cf. Tsilingiris2008) than the uncertainty of dry air's cp(T) that is discussed in the present work. Furthermore, the consideration of water vapour as a component of air requires very individual and case-specific computations of cp(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 ideal-gas behaviour introduces excessive uncertainty.

5 The θcp(T) from the temperature-dependent specific heat capacity of air
Back to toptop

Previously introduced approaches for computing the specific heat capacity of dry air call for a brief discussion on how to use the obtained cp(T) to derive the potential temperature. In the following, θcp(T) 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 cp(T) value into the conventionally used equation and defining Eq. (13) for θcp (WMO1966) may appear tempting but definitely leads to results inconsistent with θcp(T), which is based on the reference parameterisation of dry air's cp(T). Therefore, the thermodynamically consistent use of cp(T) in the derivation of θ is described in the following.

5.1 Derivation of θcp(T) based on the temperature-dependent 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 re-assessment of Eq. (9) leads to

(21) c p ( T ) T d T = R a d p p .

Integration of Eq. (21) from the basic state (p0,θcp(T)) to any other state (p, T) yields

(22) R a ln p p 0 = R a p 0 p d p p = θ c p ( T ) T c p ( T ) T d T ,

where θcp(T) 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

(23) F ( x ) = x T c p ( T ) T d T - R a ln p p 0 .

To arrive at the desired potential temperature θcp(T) for any given temperature and pressure, the equation 0=F(x) must be solved for the variable x, which is the desired θcp(T). 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 ideal-gas reference potential temperature θref is introduced, based on the formulation of the ideal-gas limit of dry air's specific heat capacity cp0(T) 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 cp(T) is to be replaced by cp0(T); i.e. for given p, T the reference potential temperature θref solves the equation

(24) 0 = F ( θ ref ) = θ ref T c p 0 ( T ) T d T - R a ln p p 0 .

The parameterisation of cp0(T) 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 cp0(T) 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 cp0(T)T may be expected to give nevertheless a good approximation even if the accuracy of cp0(T) 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 cp(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 cp 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. Deuflhard2011) provides a standard approach to numerically approximate the zero of a nonlinear equation. Proceeding from an initial guess x0, Newton's method constructs a sequence {xk}k∈ℕ defined by the recursion

(25) x k + 1 = x k - F ( x k ) F ( x k ) = x k - F ( x k ) - c p ( x k ) x k = x k c p ( x k ) c p ( x k ) + F ( x k ) = x k c p ( x k ) c p ( x k ) - R a ln p p 0 + x k T c p ( T ) T d T .

The constructed sequence {xk}k∈ℕ converges to the equation's desired zero. For the computations described here, the iteration is stopped as soon as the absolute difference xk+1-xk of two consecutive iterations falls below 10−8 K.

For the reference of air's specific heat capacity, cp0(T), the integral (23) turns out not to be explicitly solvable. Therefore, with each iteration, the solution of the integral xkTcp0(T)TdT is approximated by subdividing the entire integration range, [xk, 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 x0 for the Newton iteration, the conventional definition of θcp based on a constant specific heat capacity (WMO1966) is inserted:

(26) x 0 = T p 0 p R a 1005 J kg - 1 K - 1 = θ 1005 .

In the course of Newton's method, the sequence {xk}k∈ℕ will converge to the unique zero for any initial guess x0 due to the monotonicity of F(x). However, the right choice of the initial guess x0 substantially decreases the error of the first iteration x1, speeding up convergence to the desired zero of the function F(x). Therefore, it seems wise to use the conventional definition of θcp as the first guess for the Newton iteration (25).

Figure 4(a) Reference potential temperature θref together with the potential temperatures θ994,θ1000,θ1010 and θ1011 relying on constant cp values (the dashed lines depict the historical extremes for cp; cf. Table 1). (b) Relative differences θcp-θref/θref for the same choices cp994Jkg-1K-1,1000Jkg-1K-1,1010Jkg-1K-1,1011Jkg-1K-1 as in panel (a) between the reference potential temperature and the potential temperatures relying on constant cp values. For comparison, the relative difference θ1005-θref/θref is displayed, for which cp=1005Jkg-1K-1 corresponds to the WMO recommendation. In addition, also comparisons with θ1003.5,θ1004,θ1006.5 are included. All profiles are based on the values for temperature and pressure according to the US Standard Atmosphere. Note the linear axis scaling inside and the logarithmic scaling outside of the grey-shaded area in panel (b).


Solving the previously described root-finding problem by Newton's method over the comprehensive range of iteration steps (until the set requirement, i.e. xk+1-xk<10-8K, is fulfilled) finally leads to the reference potential temperature θref. This θref is based on the ideal-gas limit of dry air's specific heat capacity cp0(T), 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 θcp are shown based on the two (historical) extremes cp=994Jkg-1K-1 and cp=1011Jkg-1K-1 (dashed curves) and based on the range limits of more recent values cp=1000Jkg-1K-1 and cp=1010Jkg-1K-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 cp, with respect to the reference, i.e. Δθ/θref=θcp-θref/θref, is depicted in Fig. 4b. The comparison demonstrates that the θcp 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 cp1000Jkg-1K-1,1010Jkg-1K-1, the relative error level of 0.1 % is exceeded at altitudes about 5 km. While θ1000 continues to increasingly deviate from the reference, θ1010 re-enters and crosses the 0.1 % relative error interval (grey-shaded area) at altitudes between ∼19 and 21 km, before it reaches similar errors to the other θcp profiles that are based on a constant cp. Although the extreme values cp1000Jkg-1K-1,1010Jkg-1K-1 appear in recent literature, these values may be considered unrealistic. For this reason, Fig. 4b also shows the relative deviations for the values cp1003.5Jkg-1K-1,1004Jkg-1K-1,1005Jkg-1K-1, 1006.5Jkg-1K-1, which include the recommended value of the WMO (1966) and a more realistic range; i.e. cp=1005Jkg-1K-1±1.5Jkg-1K-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 cp=1004Jkg-1K-1; cf. Table 1). Until 15 km altitude, both constant cp values lead to errors of calculated θcp 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

(27) 0 = - T θ c p ( T ) T d T - R a ln p p 0

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 T,pθ may be applied mutatis mutandis to calculate the transformation θ,pT.

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 re-assess 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 xTcp0(T)TdT 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(x)=b0+b1lnx-b2+b3x+b4x2 and which is recognisable as an approximated primitive of cp0(T)T; 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 cp.

As previously discussed (cf. Sect. 5.1), the formulation of a new expression for the potential temperature based on the temperature-dependent specific heat capacity cp(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

(28) F ^ ( x ) = f ( T ) - f ( x ) - R a ln p p 0 .

Consequently, the resulting approximated reference potential temperature, i.e. the respective zero of the function F^(x), is denoted as θrefapprox.

The difference between the approximation result and the reference, i.e.

(29) θ ref - θ ref approx ,

is then referred to as the basic error of the approximation. Note that the replacement of the function F by F^ only circumvents the integration in F; the root-finding problem 0=F^(x) for the approximated reference potential temperature θrefapprox remains analytically not solvable.

Therefore, the second move towards a practical approximation procedure is to construct approximations θ(k) to the zero of F^(x) 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 θrefapprox, and, in the limit k→∞, the approximation error

(30) θ ref approx - θ ( k )

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 well-chosen θ(k) could result in an approximation error θrefθ(k) that is smaller than the basic error.

Figure 5Absolute basic error Δθ=θref-θrefapprox, cf. Eq. (29), from approximating the reference potential temperature along the US Standard Atmosphere (a) and for the extended pressure range 1000 to 0.5 hPa and temperature range 180 to 300 K (b). For orientation, the white solid line indicates the pT profile from the US Standard Atmosphere. The relative basic error θ|∕θref is shown in panel (c) for the extended pressure and temperature range.


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 θref-θrefapprox/θ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.

Figure 6(a) Relative error Δθ/θref=θ(2)-θref/θref of the second iterate θ(2), obtained with Newton's method for the ranges of pressure and temperature from 1000 to 0.5 hPa and from 180 to 300 K, respectively. Panels (b) and (c) exhibit the difference Δθ=θ1005-θref and relative difference Δθθref, respectively, on a logarithmic scale between the reference potential temperature θref and the potential temperature θ1005 based on a constant specific heat capacity (cp=1005Jkg-1K-1). For orientation, the white solid line indicates the pT profile from the US Standard Atmosphere.


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

(31) θ ref - θ ( 2 ) = θ ref - θ ref approx + θ ref approx - θ ( 2 ) ,

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.

Figure 7Flowchart guiding through the process of computing the approximation θ(2) by using Newton's formulation (C5) until its second iteration, wherein T (in K) and p (in hPa) are the atmospheric air conditions in terms of temperature and pressure, respectively, and p0 is set to 1000 hPa (WMO1966). Table C1 collects values of θref and the approximation θ(2) together with intermediate results for selected pairs of temperature and pressure to verify a computation according to this instruction.


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 θcp (WMO1966) based on a constant specific heat capacity cp=1005Jkg-1K-1. Notably, over a wide altitude range within the troposphere (i.e. for atmospheric pressures greater than ∼100 hPa), the absolute error Δθ=θ1005-θref 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

(32) c p ( θ ref ) d θ ref θ ref = c p ( T ) d T T - R a d p p ,

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 θrefapprox 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 x0 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 x1 and x2. 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.

6 The potential temperature for air as a real gas
Back to toptop

To account for real-gas 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 state-of-the-art reference equations of state for pure nitrogen (Span et al.2000), oxygen (Schmidt and Wagner1985), and argon (Tegeler et al.1999). The mixing rule and binary interaction parameters are taken from the GERG-2008 model (Kunz and Wagner2012). From its definition in terms of an isentropic process, the potential temperature θreal(T, p) is defined implicitly by

(33) s ( θ real , p 0 ) = s ( T , p ) ,

where s is the specific entropy. Calculating θreal(T, p) is a two-step 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 p0 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 Wagner1985; 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 ideal-gas calculations; the correlations for cp0(T) for N2 and O2 are extrapolations beyond 2000 K. However, since the same ideal-gas values are used in the real-gas calculations, any inaccuracy in cp0(T) will cancel when evaluating the difference between ideal-gas and real-gas values of θ.

Figure 8Difference θrealθref reflecting the deviation of the potential temperature θreal, based on the properties of air behaving as a real gas under variable temperature and pressure, from the herein derived potential temperature expression θref for the ideal-gas limit of the air's specific heat capacity cp0(T). (a) Difference along the profile of the US Standard Atmosphere. (b) Relative difference in pT coordinates covering any combination of atmospherically relevant temperatures and pressures.


Figure 8 illustrates the comparison between the real-gas potential temperature θreal and the ideal-gas reference potential temperature θref. Figure 8a shows the difference θrealθref along the pT profile of the US Standard Atmosphere and Fig. 8b accounts for any pT 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 30×10-5=0.03% for the relative difference. As may be anticipated from the deviation of cp0 shown in Fig. 3 at low temperatures both from the experimentally determined values (which may be inaccurate) and from the REFPROP data, the real-gas effect on the specific heat capacity of dry air tends to increase towards the coldest gas temperatures. However, the difference between the real- and ideal-gas 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 real-gas and the ideal-gas specific heat capacity becomes insignificant. Within the intermediate (stratospheric) region, the low pressures (and thus the low air densities) cause the ideal-gas 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 pT; 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.

7 Implications on the use of the potential temperature
Back to toptop

As previously shown, the newly defined reference potential temperature θref deviates most from the WMO-defined 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. θrefz and θ1005z. 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 well-founded, 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 N2 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 Nfalse2. 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 Klemp1982; Marquet and Geleyn2013; Wallace and Hobbs2006; Ambaum2010). Retaining the assumption of hydrostatic conditions, the defining formula yields

(34) N 2 = g T T z + g c p ( T ) ,

where the temperature-dependent specific heat capacity cp(T) was implied, and which quantifies the balance between the actual temperature stratification Tz and the dry-adiabatic lapse rate -gcp(T) (e.g. Holton2004).

Figure 9Vertical profiles of (a) the temperature up to 85 km altitude as typical for mid-latitudes in June (red curve) and December (blue curve). (b) Resulting wrong Brunt–Väisälä frequency Nfalse2 (dashed lines) and the true Brunt–Väisälä frequency N2 (solid lines) for the two temperature profiles from panel (a).


To illustrate the deviation of Nfalse2 from N2, 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 Ortland2003) data and are assumed as typical at mid-latitudes 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 Nfalse2 and the true Brunt–Väisälä frequency N2 are shown in Fig. 9b. Evidently, the values of Nfalse2 (dashed lines) deviate significantly from N2 (solid lines) and increasingly so towards higher altitudes above 15 km. However, the absolute deviation |N2-N10052|, using N10052 as calculated with θ1005 in accordance with Eq. (2), does not exceed 1.6×10-6s-2 (not shown), indicating that N10052 is a good representation of N2 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 cp, to avoid aberrations and erroneous conclusions.

7.2 The potential vorticity

Ertel's potential vorticity (e.g. Ertel1942; Hoskins et al.1985; Schubert et al.2004; Holton2004) may be defined as the potential vorticity of the dry-air potential temperature by

(35) PV ( θ ) = 1 ρ 2 Ω + × u θ .

In this definition, 2Ω+×u is the absolute vorticity, Ω denotes Earth's angular velocity, u the three-dimensional wind vector, and ρ the air density (see e.g. Hoskins et al.1985; Cotton et al.2011; Marquet2014). Since Eq. (35) represents the defining equation for Ertel's potential vorticity, the two potential vorticities

(36) PV ref = PV ( θ ref ) , PV 1005 = PV ( θ 1005 )

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

(37) PV ( θ ) = 2 sin ( ϕ ) ρ 2 π t E θ z

for a position on Earth with geographical latitude ϕ and tE=24 h, the duration of one rotation of the Earth.

Figure 10(a) Vertical profiles of the potential vorticity PVref computed with θref (solid lines), and PV1005 computed with θ1005 (dashed lines), for an atmosphere at rest along the temperature profiles from Fig. 9a for June (red lines) and December (blue lines). Since the temperature profiles are representative for mid-latitudes on the Northern Hemisphere, the geographical latitude in Eq. (37) was set to 52 N. (b) Relative deviation |PVref-PV1005|/PVref of the potential vorticity profiles from panel (a).


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 northern-hemispheric mid-latitudes, 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 N2 (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, WMO1957) 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 high-altitude research aircraft such as the G-550 HALO (e.g. Wendisch et al.2016; Voigt et al.2017), the NASA WB-57 or ER-2 (e.g. Murphy et al.2007; Dessler2002), the M-55 Geophysica (Curtius et al.2005; Borrmann et al.2010; Frey2011), balloon-borne platforms (Lary et al.1995; Vernier et al.2018), or satellite-based vertical profiles (e.g. Davies et al.2006; Spang et al.2005) require consideration of the systematic error in θ if calculated as θcp in compliance with the definition by the WMO (1966). The possibly inconsistent use of a constant cp value of 1004 or 1005Jkg-1K-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 high-altitude 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 spacecraft-bound vertical soundings (e.g. from ASTRO-SPAS, SCIAMACHY, or ENVISAT), the error in the potential temperature determined by θcp 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 dqdt supplied to a given air parcel, e.g. by radiative heating, and are given in units of Jkg-1s-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,

(38) AHR ref d q d t = d T d t = 1 c p 0 ( T ) d q d t , AHR 1005 d q d t = d T d t = 1 1005 J kg - 1 K - 1 d q d t .

Again, the distinction was made between the temperature-dependent cp0(T) and the constant cp=1005Jkg-1K-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

(39) AHR 1005 ( x ) - AHR ref ( x ) AHR ref ( x ) = c p 0 ( T ) 1005 J kg - 1 K - 1 - 1 .

Apart from the absolute heating rates for the change in absolute temperature, the change in potential temperature due to a diabatic heating rate dqdt 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, Chipperfield2006; 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

(40) d q T = c p ( T ) T d T - R a d p p .

Comparing the right-hand 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

(41) d q T = c p ( θ ref ) d θ ref θ ref .

Consequently, the following two diabatic heating rates,

(42) d θ ref d t = θ ref c p 0 ( θ ref ) T d q d t = HR ref d q d t , d θ 1005 d t = θ 1005 ( 1005 J kg - 1 K - 1 ) T d q d t = HR 1005 d q d t ,

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

(43) HR 1005 ( x ) - HR ref ( x ) HR ref ( x ) = θ 1005 θ ref c p 0 ( θ ref ) 1005 J kg - 1 K - 1 - 1 .

Figure 11(a) Monthly averaged temperatures profiles for 52N. (b) The relative differences between the absolute heating rates, defined in Eq. (39). (c) The relative differences between the heating rates (Eq. 43) for the potential temperatures θ1005 and θref. (d) The resulting relative deviation of the potential temperatures θ1005,θref after 24 h of heating with constant diabatic heating dqdt and the resulting heating rates HRref, HR1005 at constant pressure.


In order to judge the magnitudes of the relative differences (Eqs. 39 and 43), the monthly averaged temperature profiles from ERA-Interim (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 dqdt 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 AHR1005 or AHRref are smaller than 0.044 K. However, the differences in the potential temperatures θ1005,θref, computed with the heating rates HRref, HR1005, 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 temperature-dependent instead of the constant cp. 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 upward-propagating seasonal signal in tropical stratospheric water vapour (the so-called 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 temperature-dependent cp is comparably small. However, in cases of smaller inter-model differences the additional cp-related uncertainty needs to be assessed.

Note that the determination of absolute temperatures T1005,Tref which correspond to the resulting potential temperatures θ1005,θref after 24 h differs by less than 0.014 K (not shown).

8 Summary and conclusions
Back to toptop

Under the assumption that dry air is an ideal gas, a re-assessment 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 state-of-the-art parameterisation of the ideal-gas 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 real-gas 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 θcp (with constant cp=1005Jkg-1K-1, as recommended by the World Meteorological Organisation, WMO1966) 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 ideal-gas properties of dry air requires the integration of Gibbs' equation and the subsequent solution of the resulting nonlinear equation. With a constant cp, 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 cp(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 θref-θrefapprox inherent in the integral's basic approximation and the error θrefapprox-θ(k) 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 higher-order polynomial.

One of the foremost implications of the re-assessed 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 re-assessed 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 1010Jkg-1K-1; see Table 1).

Significant errors and biases may arise if, for instance, the conventional derivation of θ (WMO1966) is used together with values for air's specific gas constant (Ra) or air's specific heat capacity (cp), 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 (p0) may cause an additional deviation of the resulting θ. Thus, the re-assessment 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 well-founded, 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 re-assessment could take into account the current state of knowledge regarding the accuracy of thermodynamic variables and substance-related 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 ideal-gas 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.

Appendix A: Derivation of the specific heat capacity from thermodynamics
Back to toptop

In the following, the derivation of the air's specific heat capacities CV, Cp (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 ideal-gas law

(A1) p V = N R T ,

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

(A2) d U = d Q - p d V ,

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

(A3) d Q - p d V = U T V , N d T + U V T , N d V

and subsequently

(A4) d Q = U T V , N d T + p + U V T , N d V .

If the system's volume is held constant, Eq. (A4) represents the definition of the constant-volume heat capacity CV in molar units, i.e.

(A5) d Q = U T V , N d T = C V ( p , T ) d T .

Alternatively, assuming the system's pressure as constant, its volume is variable with total derivative

(A6) d V = V T p , N d T + V p T , N d p = 0 = V T p , N d T

and therefore

(A7) d Q = U T V , N d T + p + U V T , N d V = U T V , N d T + p + U V T , N V T p , N d T = U T V , N + p + U V T , N V T p , N d T = C p ( p , T ) d T ,

defining the isobaric molar heat capacity Cp. 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. Fay1965).

Using this result, together with Eq. (A7) and the ideal-gas law (A1), it follows

(A8) C p = U T V , N + p + U V T , N V T p , N = U T V , N + p V T p , N = U T V , N + T p V p , N = U T V , N + T N R T p , N = U T V , N + N R .

In the previous computations, there is no restriction on the temperature dependence of the internal energy U(T). Therefore, even by assuming ideal-gas behaviour, the specific heat capacity Cp in Eq. (A8) is in general a function of temperature.

Appendix B: Sensitivity of the conventional definition of θ to perturbations of cp
Back to toptop

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 cp as a variable, the sensitivity of θcp (Eq. 13) to a small perturbation δ of cp is described by its Taylor expansion

(B1) θ c p + δ = θ c p + θ c p c p δ + O δ 2 = θ c p - θ c p R a c p 2 ln p 0 p δ + O δ 2 .

For any constant value of the specific heat capacity cp and for a minor perturbation δ, the second summand within the expansion (B1) remains small for small values of lnp0p. If the interval between the two pressure levels is very narrow, i.e. pp0, the expression lnp0p 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. lnp0p-, 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.

Appendix C: Approximate computation of the reference potential temperature
Back to toptop

This section summarises the detailed steps of approximating the function F(x), defined in Eq. (23), by F^(x), 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),

(C1) h ( x ) = T 1 x c p ( T ) T d T ,

with T1=180 K, the function F(x) may be rearranged as

(C2) F ( x ) = x T c p ( T ) T d T - R a ln p p 0 = h ( T ) - h ( x ) - R a ln p p 0 .

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

(C3) f ( x ) = b 0 + b 1 ln x - b 2 + b 3 x + b 4 x 2 .

Notably, if cp is constant, this function reduces to an exact primitive of the integrand cpT with b3=b4=0. Moreover, in this case, the resulting root-finding 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 cp(T) in Eq. (C1) is replaced by the ideal-gas limit of air's specific heat capacity cp0(T). The integration interval [T1, x] with T1x2000K 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 least-squares problem, the coefficients in Eq. (C3) for the approximation of h(x) by the function f(x) are estimated as

(C4) b 0 = - 4072.2121328563667 , b 1 = 797.09247926609601 , b 2 = 29.587047521428016 , b 3 = 0.41981158226925142 , b 4 = - 5.1008025097060311 × 10 - 5 .

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 ±1Jkg-1K-1. Exclusively at temperatures below 190 K, the approximation error rapidly rises above 1Jkg-1K-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).

Figure C1(a) Numerically evaluated function h(x) together with its approximation f(x); (b) the absolute approximation error h(x)−f(x).


Table C1Values of the new reference potential temperature θref, together with the first two iterates θ(1), θ(2) using Newton's method and the first iterate θHouseholder(1) using Householder's method for five pairs of temperature and pressure along the US Standard Atmosphere. The computed values are rounded to three digits.

Download Print Version | Download XLSX

C2 Finalised approximation of the reference potential temperature

As discussed in Sect. 5.1, the new formulation of the potential temperature based on the temperature-dependent specific heat capacity cp(T) requires solving the root-finding problem 0=F(x), where the function F(x) is defined in Eq. (23). However, since F(x) contains an integral that complicates the root-finding 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 F^(x) as defined in Eq. (28), and the zero of the equation 0=F^(x) is denoted as θrefapprox.

The equation 0=F^(x) is still not analytically solvable, so Newton's method is once more required. Using again x0=θ1005 as the initial guess, cf. Eq. (26), the iteration sequence for Newton's method is given by the recursion

(C5) x k + 1 = x k - F ^ ( x k ) F ^ ( x k ) = x k - f ( T ) - f ( x k ) - R a ln p p 0 - f ( x k ) = x k - R a ln p p 0 - f ( T ) + f ( x k ) f ( x k ) .

Instead of this standard formulation of Newton's method (C5), Householder's formulation

(C6) x k + 1 = x k - F ^ ( x k ) F ^ ( x k ) - F ^ ′′ ( x k ) 2 F ^ ( x k ) F ^ ( x k ) F ^ ( x k ) 2 = x k - R a ln p p 0 - f ( T ) + f ( x k ) f ( x k ) - f ′′ ( x k ) 2 f ( x k ) R a ln p p 0 - f ( T ) + f ( x k ) f ( x k ) 2

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

(C7) f ( x ) = b 1 x - b 2 + b 3 + 2 b 4 x , f ′′ ( x ) = 2 b 4 - b 1 x - b 2 2 .

The final step on the way to formulate a new expression for the potential temperature requires defining one of the iterates xk 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 θHouseholder(1) 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.

Appendix D: Approximation error for the reference potential temperature
Back to toptop

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 θref-θrefapprox and the approximation error that results from the approximation sequence θrefapprox-θ(k), 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 F^(x), and the approximation sequences θ(k) converge to θrefapprox for k→∞. Consequently, the approximation error θrefapprox-θ(k) tends to zero for k→∞.

Figure D1Total relative error along the US Standard Atmosphere arising from the iteration process by declaring (a) the first iterate θ(1) or (b) the second iterate θ(2) as the final approximation to the reference potential temperature θref. Red curves: iterates computed using Newton's method (C5); blue curves: iterates computed using Householder's method (C6). Note the different range of the abscissae.


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 θref-θ(1)/θref of the first iterate (Fig. D1a) and θref-θ(2)/θ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 θref-θ(2)/θref relative to the ultimate reference potential temperature (indicated by a smaller distance to the dashed zero-line above 45 km altitude). Nevertheless, the relative approximation error, θrefapprox-θ(2)/θ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.

Figure D2Relative error of the second iterates θ(2) with (a) Newton's method and (b) Householder's method for the ranges of pressure and temperature from 1000 to 0.5 hPa and from 180 to 300 K, respectively. (c) The absolute error arising from the first iterate θ(1) with Householder's method. The white solid line indicates the pT profile from the US Standard Atmosphere. Note the different ranges of the Δθ scales.


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 (∼50km). 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.

Appendix E: The derivative of the reference potential temperature
Back to toptop

As discussed in Sect. 5.1, the new reference potential temperature is defined as the zero of the function

(E1) F ( x , p , T ) = x T c p ( T ) T d T - R a ln p p 0

for given values of pressure p and temperature T; see Eq. (23). More precisely, for varying p, T, a function (p,T)θref(p,T) is implicitly defined by the equation

(E2) F ( θ ref ( p , T ) , p , T ) = 0 .

According to the implicit function theorem (e.g. Protter and Morrey1985, chap. 7), Eq. (E2) is uniquely solvable for θref(p, T); i.e. the function (p,T)θref(p,T) actually exists as a differentiable function of (p, T), if the condition Fθ0 holds. According to Eq. (E1), this partial derivative equals

(E3) F θ ( θ ref , p , T ) = - c p ( θ ref ) θ ref ,

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

(E4) θ ref p ( p , T ) , θ ref T ( p , T ) = - F θ ( θ ref , p , T ) - 1 F p ( θ ref , p , T ) , F T ( θ ref , p , T ) = θ ref c p ( θ ref ) - R a p , c p ( T ) T = - R a c p ( θ ref ) θ ref p , θ ref T c p ( T ) c p ( θ ref ) .

Note that these partial derivatives coincide with the partial derivatives of θcp in the case of a constant specific heat capacity. Using the partial derivatives (Eq. E4), the total differential of θref may be written as

(E5) d θ ref = θ ref p d p + θ ref T d T = - R a c p ( θ ref ) θ ref p d p + θ ref T c p ( T ) c p ( θ ref ) d T


(E6) c p ( θ ref ) d θ ref θ ref = c p ( T ) d T T - R a d p p .
Data availability
Back to toptop
Data availability. 

The data used in this study are available from the corresponding author upon request (

Author contributions
Back to toptop
Author contributions. 

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 real-gas potential temperatures. FP contributed the implications on diabatic heating. MB and RW wrote the manuscript with contributions and reviews from all authors.

Competing interests
Back to toptop
Competing interests. 

The authors declare that they have no conflict of interest.

Back to toptop

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 (, 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 (MS-GWaves) and through grants AC 71/12-2 and SP 1163/5-2.

Financial support
Back to toptop
Financial support. 

This research has been supported by the Deutsche Forschungsgemeinschaft (grant no. 257899354, AC 71/12-2, and SP 1163/5-2) and the Bundesministerium für Bildung, Wissenschaft, Forschung und Technologie (grant no. 01LG1205A).

Review statement
Back to toptop
Review statement. 

This paper was edited by Timothy Garrett and reviewed by Pascal Marquet and one anonymous referee.

Back to toptop

Ambaum, M. H. P.: Thermal Physics of the Atmosphere, John Wiley & Sons, Ltd., Chichester, UK,, 2010. a, b, c

Awano, S.: JS-Diagrams for Air, Report of Aeronautical Research Institute, Tokyo Imperial University, 11, available at: (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,<1046:TCOEPT>2.0.CO;2, 1980. a

Borchert, S., Zhou, G., Baldauf, M., Schmidt, H., Zängl, G., and Reinert, D.: The upper-atmosphere extension of the ICON general circulation model (version: ua-icon-1.0), Geosci. Model Dev., 12, 3541–3569,, 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: in-situ measurements of submicron particle abundance and volatility, Atmos. Chem. Phys., 10, 5573–5592,, 2010. a

Brasseur, G. P. and Solomon, S.: Aeronomy of the Middle Atmosphere, Springer,, 2005. a

Brusseau, M., Pepper, I. L., and Gerba, C. P.: Environmental and Pollution Science, 3rd Edn., Elsevier,, 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,, 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.,, 2015. a, b

Chang, X., Zhao, W., Zhang, Z., and Su, Y.: Sap flow and tree conductance of shelter-belt in arid region of China, Agr. Forest Meteorol., 138, 132–141,, 2006. a

Chipperfield, M. P.: New version of the TOMCAT/SLIMCAT off-line chemical transport model: Intercomparison of stratospheric tracer experiments, Q. J. Roy. Meteor. Soc., 132, 1179–1203,, 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,, 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 MIPAS-E satellite measurements in winter 2002/2003, Atmos. Chem. Phys., 6, 3149–3161,, 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., Monge-Sanz, B. M., Morcrette, J.-J., Park, B.-K., Peubey, C., de Rosnay, P., Tavolato, C., Thépaut, J.-N., and Vitart, F.: The ERA-Interim reanalysis: configuration and performance of the data assimilation system, Q. J. Roy. Meteor. Soc., 137, 553–597,, 2011. a

Dessler, A. E.: The effect of deep, tropical convection on the tropical tropopause layer, J. Geophys. Res.-Atmos., 107, ACH 6-1–ACH 6-5,, 2002. a

Deuflhard, P.: Newton Methods for Nonlinear Problems, Vol. 35 of Springer Series in Computational Mathematics, Springer-Verlag, Berlin Heidelberg,, 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 Brunt-Väisälä Frequency, J. Atmos. Sci., 39, 2152–2158,<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, Addison-Wesley, 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, Deep-Sea Res. Pt. I, 55, 1639–1671,, 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 High-Resolution Rovibronic Spectra and Definitive Ideal-Gas Thermochemistry of the 16O2 Molecule, J. Phys. Chem. Ref. Data, 48, 023101,, 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,, 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,<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.: Stratosphere-troposphere exchange, Rev. Geophys., 33, 403–439,, 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., Friedl-Vallon, 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,, 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 Physikalisch-Technischen 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 GERG-2008 Wide-Range Equation of State for Natural Gases and Other Mixtures: An Expansion of GERG-2004, J. Chem. Eng. Data, 57, 3032–3091,, 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.: Three-dimensional tracer initialization and general diagnostics using equivalent PV latitude–potential-temperature coordinates, Q. J. Roy. Meteor. Soc., 121, 187–210,, 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,, 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 Properties-REFPROP, Version 10.0, National Institute of Standards and Technology,, 2018. a, b, c

Li, C., Ingersoll, A. P., and Oyafuso, F.: Moist Adiabats with Multiple Condensing Species: A New Theory with Application to Giant-Planet Atmospheres, J. Atmos. Sci., 75, 1063–1072,, 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 FIRE-I data flights, Q. J. Roy. Meteor. Soc., 137, 768–791,, 2011. a

Marquet, P.: On the definition of a moist-air potential vorticity, Q. J. Roy. Meteor. Soc., 140, 917–929,, 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,, an english translation is available at arXiv: (last access: 2 November 2020), 2019. a

Marquet, P. and Geleyn, J.-F.: On a general definition of the squared Brunt-Väisälä frequency associated with the specific moist entropy potential temperature, Q. J. Roy. Meteor. Soc., 139, 85–100,, 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,, 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,<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üller-Wodarg, 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,, 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 NA for the revision of the SI, Metrologia, 55, L13–L16,, 2018. a, b

Ooyama, K. V.: A Thermodynamic Foundation for Modeling the Moist Atmosphere, J. Atmos. Sci., 47, 2580–2593,<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,<2073:ADATFF>2.0.CO;2, 2001. a, b

Poisson, S. D.: Traité de mécanique, Bachelier, Paris, 2nd Edn., available at: (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 long-lived trace species in the Chemical Lagrangian Model of the Stratosphere (CLaMS), Geosci. Model Dev., 7, 2895–2916,, 2014. a

Protter, M. H. and Morrey, C. B.: Intermediate Calculus, Undergraduate Texts in Mathematics, 2nd Edn., Springer, Berlin, Heidelberg,, 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

Reinke-Kunze, 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,, 2007. a

Roebuck, J. R.: The Joule-Thomson Effect in Air, P. Am. Acad. Arts Sci., 60, 537–596,1925. a

Roebuck, J. R.: The Joule-Thomson 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,, 2011. a

Scheel, K. and Heuse, W.: Die spezifische Wärme der Luft bei Zimmertemperatur und bei tiefen Temperaturen, Ann. Phys., 342, 79–95,, 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,, 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 twenty-one of Ertel's papers on geophysical fluid dynamics, Meteorol. Z., 13, 527–576,, 2004. a, b

Schumann, U.: Atmospheric Physics; Background – Methods – Trends, Research Topics in Aerospace, Springer-Verlag, Berlin Heidelberg,, 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 time-split nonhydrostatic atmospheric model for weather research and forecasting applications, J. Comput. Phys., 227, 3465–3485,, 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/TN-468+STR), University Corporation for Atmospheric Research,, 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,, 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,, 2005. a

Spichtinger, P.: Shallow cirrus convection – a source for ice supersaturation, Tellus A, 66, 19937,, 2014. a

Stamnes, K., Thomas, G. E., and Stamnes, J. J.: Radiative Transfer in the Atmosphere and Ocean, Cambridge University Press,, 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., Kromp-Kolb, 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.: Stratosphere-troposphere exchange: A review, and what we have learned from STACCATO, J. Geophys. Res.-Atmos., 108, 8516,, 2003. a

Swinbank, R. and Ortland, D. A.: Compilation of wind data for the (UARS) Reference Atmosphere Project, J. Geophys. Res., 108, 4615,, 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,, 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,, 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: (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 lce-Liquid Water Potential Temperature as a Thermodynamic Variable In Deep Atmospheric Models, Mon. Weather Rev., 109, 1094–1102,<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,, 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: (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,, 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.: Ml-Cirrus the Airborne Experiment on Natural Cirrus and Contrail Cirrus with the High-Altitude Long-Range Research Aircraft Halo, B. Am. Meteorol. Soc., 98, 271–288,, 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,, 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 fast-flying aircraft, Atmos. Meas. Tech., 9, 5135–5162,, 2016. a

Wendisch, M. and Brenguier, J.-L.: Airborne Measurements for Environmental Research: Methods and Instruments, Wiley‐VCH, Weinheim,, 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 ACRIDICON-CHUVA campaign: Studying tropical deep convective clouds and precipitation over Amazonia using the new German research aircraft HALO, B. Am. Meteorol. Soc., 97, 1885–1908,, 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,, 1989. a

Witkowski, A. W.: I. Thermodynamic properties of air, The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, 42, 1–37,, 1896. a, b

WMO: Meteorology – A three-dimensional science, WMO Bull., 6, 134–138, 1957. a

WMO: International Meteorological Tables, WMO-No.188.TP97, Secretariat of the World Meteorological Organization, Geneva, Switzerland, available at: (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 Non-hydrostatic) modelling framework of DWD and MPI-M: Description of the non-hydrostatic dynamical core, Q. J. Roy. Meteor. Soc., 141, 563–579,, 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 (Reinke-Kunze2013) and they married in the year 1913 (Hallam1975).


Cf. Bauer (1908), where, for the first time, the potential temperature and the entropy are set in a relationship.


These early applications of entropy in meteorology are also documented in Marquet (2019).


The value of R is now defined exactly, cf. Tiesinga et al. (2020) and Newell et al. (2018), and is used in Eq. (6).

Publications Copernicus
Short summary
The potential temperature is routinely used in atmospheric science. We review its derivation and suggest a new potential temperature, based on a temperature-dependent parameterization of the dry air's specific heat capacity. Moreover, we compare the new potential temperature to the common one and discuss the differences which become more important at higher altitudes. Finally, we indicate some consequences of using the new potential temperature in typical applications.
The potential temperature is routinely used in atmospheric science. We review its derivation and...
Final-revised paper