## 1. Introduction

In a recent paper, Mrowiec et al. (2016, hereafter MPZ) investigated the thermodynamic properties of a three-dimensional hurricane simulation. The MPZ paper focuses on isentropic analysis based on the conditional averaging of the mass transport with respect to the equivalent potential temperature *θ*_{e}.

The methodology employed in MPZ originated in the concept of a thermal or Carnot heat engine applied to hurricanes by Riehl (1950), Kleinschmidt (1951), Emanuel (1986, 1988, 1991, hereafter E86, E88, E91), and Gray (1994) and to a study of steam cycle and entropy budgets described in Pauluis (2011, hereafter P11). In these papers, it is assumed that the equivalent potential temperature *θ*_{e} or the saturation value *θ*_{es} are logarithmic measurements of the moist-air entropy, according to *s* ≈ *S*_{0m} + *C*_{pm} ln(*θ*_{e}) or *s** ≈ *S*_{0m} + *C*_{pm} ln(*θ*_{es}), which are valid for unsaturated or saturated moist air, respectively. It is thus assumed that isentropic surfaces are represented by constant values of *θ*_{e} or *θ*_{es} despite the terms *S*_{0m} and *C*_{pm}, which may depend on the water content and may thus vary with space or time, making the link between *s* and *θ*_{e} or *θ*_{es} unclear, at the very least.

The equivalent potential temperatures are sometimes replaced by the moist static energy counterparts *C*_{pm}*T* + *L*_{υ}*q*_{υ} + *gz* or *C*_{pm}*T* + *L*_{υ} *q*_{sw} + *gz*, where *C*_{pm} may depend on water vapor or condensed water contents, *L*_{υ} is the latent heat of vaporization, and *q*_{υ} or *q*_{sw} is the unsaturated or saturated specific content of water vapor, respectively.

The use of static energies or equivalent potential temperatures as a proxy for the moist-air entropy has been described and discussed at some length in other studies over the years: Malkus (1958), Miller (1958), Green et al. (1966), Darkow (1968), Madden and Robitaille (1970, 1972), Levine (1972), Betts and Dugan (1973), Betts (1973, 1974, 1975), Rotunno and Emanuel (1987), Emanuel et al. (1987), E88, Emanuel (1989), Ooyama (1990), Peixoto et al. (1991), E94, Emanuel (1995), Rennó and Ingersoll (1996), Emanuel and Bister (1996), Bister and Emanuel (1998), Pauluis et al. (2000), Rennó (2001), Pauluis and Held (2002a,b), Goody (2003), Emanuel (2003, 2004), Bannon (2005), Romps (2008), Smith et al. (2008), Pauluis et al. (2008, 2010), Romps and Kuang (2010), Romps (2010), Emanuel and Rotunno (2011), Raymond (2013), and Pauluis and Mrowiec (2013).

In the present paper, it is shown that the way in which the moist entropy and the equivalent potential temperatures are defined may lead to opposite results in studies of isentropic processes in hurricanes. It is shown that, the more the total water varies with space, the more the isentropes differ, leading to large discrepancies in computations and plots of 1) the isentropes themselves; 2) the changes with space of the equivalent potential temperatures and the moist-air entropies for a thermodynamic heat cycle; 3) the isentropic streamfunction Ψ defined in Pauluis and Mrowiec (2013) and studied in MPZ; 4) the heat input and work function studied in E88, E91, Emanuel (1995), Rennó and Ingersoll (1996), Pauluis and Held (2002a), and P11; and 5) the efficiency of a (moist) steam cycle defined in P11.

The definitions and the values of the heat input and work functions, the streamfunction, the efficiency, and the moist-air entropy itself cannot be subject to arbitrariness. This paper thus advocates using the absolute value of the moist-air entropy to study the thermodynamic properties of atmospheric processes or to perform isentropic analyses.

It is explained in appendix A that the third law of thermodynamics leads to such a physically founded definition of the absolute entropy, which can be computed for all substances. The third-law values are determined by using experimental calorimetric values of the specific heat and with the hypothesis that the same (zero) value of entropy is reached for the most stable crystalline state of all substances when approaching absolute zero temperature. Appendix A also explains that the third-law values are very close to values of entropies determined from theoretical methods based on statistical and quantum physics, for both mono- and polyatomic molecules. Accordingly, the term “third law” entropy will herein denote either the theoretical or the experimental absolute value of the entropy.

Appendix B recalls that several applications of the third law of thermodynamics to atmospheric studies already exist: Hauf and Höller (1987), Bannon (2005), and Marquet (2011, hereafter M11), a process recently improved by Marquet and Geleyn (2015) and Marquet (2015b, 2016b,c). It is thus possible to study the absolute moist-air entropy defined by *s*(*θ*_{s}) = *s*_{ref} + *c*_{pd} ln(*θ*_{s}), where both *s*_{ref} and *c*_{pd} are constant and where the third-law potential temperature *θ*_{s} is synonymous with the absolute moist-air entropy.

The paper is organized as follows. The definitions of most existing equivalent potential temperatures *θ*_{e} or *θ*_{es} are recalled in section 2a, together with the third-law value *θ*_{s}. Several of the existing moist-air entropies are listed in section 2b, including the third-law value *s*(*θ*_{s}). A dataset derived from a simulation of Hurricane Dumile is presented in section 2c and isentropic surfaces corresponding to *θ*_{e} and *θ*_{s} are computed, plotted, and compared in sections 3a and 3b for a cross section of this hurricane. Differences in the associated moist-air entropies are described in section 3c and it is shown that the two isentropic streamfunctions computed in section 3d for both *θ*_{e} and *θ*_{s} exhibit large differences. The heat input and work function are computed in section 3e for a steam cycle in the so-called temperature–entropy diagram. The efficiencies of such steam cycles are computed and compared in section 3f for both *θ*_{e} and *θ*_{s}. Finally, conclusions are presented in section 4.

## 2. Data and method

### a. Moist-air potential temperatures

*θ*

_{e}defined by Eq. (4.5.11) in Emanuel (1994, hereafter E94) can be written aswhere

*T*is the temperature,

*p*

_{0}is the standard pressure,

*p*

_{d}is the dry-air pressure, and

*R*

_{d}and

*R*

_{υ}are the gas constant of water vapor and dry air, respectively. The specific heat

*c*

_{pd}+

*r*

_{t}

*c*

_{l}depends on the values for dry air (

*c*

_{pd}) and liquid water (

*c*

_{l}). The mixing ratios

*r*

_{υ}and

*r*

_{t}represent water vapor and total water, respectively. The quantity

*L*

_{υ}is the latent heat of vaporization. The relative humidity with respect to liquid water

*H*

_{l}=

*e*/

*e*

_{sw}is the ratio of the water vapor pressure (

*e*) to the saturated value (

*e*

_{sw}).

*p*is the total pressure and

*r*

_{sw}(

*T*,

*p*) is the saturation mixing ratio at temperature

*T*and pressure

*p*.

*θ*=

*T*(

*p*

_{0}/

*p*)

^{κ}with

*κ*=

*R*

_{d}/

*c*

_{pd}≈ 0.286.

*θ*

_{e}on page 1860 of MPZ can be written asThe difference between

*θ*

_{e/MPZ}and

*θ*

_{e/E94}is that

*p*

_{d}is replaced by

*p*in the Exner function. The differences between

*θ*

_{e/MPZ}and

*θ*

_{e/E86}are that

*r*

_{sw}is replaced by

*r*

_{υ}and the additional term depending on

*H*

_{l}is included. The differences between

*θ*

_{e/MPZ}and

*θ*

_{e/B73}are that

*c*

_{pd}is replaced by

*q*

_{υ}is replaced by

*r*

_{υ}, and the additional term depending on

*H*

_{l}is included.

Even though the potential temperatures truly considered in MPZ are *θ*_{e/E94}, it is interesting to consider *θ*_{e/MPZ} given by Eq. (5) in order to demonstrate that the a priori small change of *p*_{d} into *p* may have a large impact on the definition of isentropic processes.

*λ*=

*c*

_{pυ}/

*c*

_{pd}− 1 ≈ 0.837,

*η*=

*R*

_{υ}/

*R*

_{d}≈ 1.608,

*ε*=

*R*

_{d}/

*R*

_{υ}≈ 0.622,

*δ*=

*η*− 1 ≈ 0.608, and

*γ*=

*R*

_{υ}/

*c*

_{pd}≈ 0.46. The term Λ

_{r}= [(

*s*

_{υ})

_{r}− (

*s*

_{d})

_{r}]/

*c*

_{pd}≈ 5.87 depends on the third-law values of the reference entropies for water vapor and dry air, and

*L*

_{s}is the latent heat of sublimation. The water vapor, liquid, ice, and total specific contents

*q*

_{υ},

*q*

_{l},

*q*

_{i}, and

*q*

_{t}=

*q*

_{υ}+

*q*

_{l}+

*q*

_{i}replace the mixing ratios involved in most of the previous formulations.

The water vapor and dry-air reference entropies (*s*_{υ})_{r} = *s*_{υ}(*T*_{r}, *e*_{r}) and (*s*_{d})_{r} = *s*_{d} (*T*_{r}, *p*_{r} − *e*_{r}) depend on the reference temperature *T*_{r} = 273.15 K, the total reference pressure *p*_{r} = 1000 hPa, and the water vapor partial pressure *e*_{r}. However, it is shown in M11 that *θ*_{s/M11} defined in Eq. (6) is independent of *T*_{r} and *p*_{r} if the reference mixing ratio *r*_{r} = *εe*_{r}/(*p*_{r} − *e*_{r}) ≈ 3.82 g kg^{−1} corresponds to the saturating pressure value at *T*_{r}: *e*_{r} = *e*_{sw}(*T*_{r}) = 6.11 hPa.

The first two lines of Eq. (6) are derived in (40) in M11. The four terms in the last line of Eq. (6) derived in Marquet (2016c) are improvements with respect to M11. They take possible nonequilibrium processes into account, such as under- or supersaturation with respect to liquid water (*H*_{l} ≠ 1) or ice (*H*_{i} ≠ 1) and possible temperatures of rain *T*_{l}, or snow, *T*_{i}, which may differ from the *T* for dry air and water vapor.

The advantage of the term *q*_{l} replaces *r*_{υ} in the exponent, making no impact in clear air, or under- or supersaturated moist regions (where *H*_{l} ≠ 1 and *r*_{υ} may be large, but where *q*_{l} = 0) and having less impact in cloud in under or supersaturated regions (where *H*_{l} ≠ 1 but where, typically, *q*_{l} < *r*_{υ}).

*θ*

_{s}are derived in Marquet (2015b, 2016b,c) and lead towhere

*r*

_{*}≈ 0.0124 kg kg

^{−1}. Both (

*θ*

_{s})

_{1}and (

*θ*

_{s})

_{2}must be multiplied by the last line of Eq. (6) if nonequilibrium processes are to be described.

### b. The moist-air entropies

*s*

_{ref}≈ 1139 J K

^{−1}kg

^{−1}and

*c*

_{pd}≈ 1004 J K

^{−1}kg

^{−1}are constant, making

*θ*

_{s}a true equivalent of the specific moist-air entropy. The formulation (11) is expressed “per unit of dry air,” in order to be better compared with the entropies computed in other studies such as E94 or MPZ.

*s*

_{ref}or

*c*

_{pd}(often both of them) replaced by quantities that depend on the total-water mixing ratio

*r*

_{t}. This is true in Eq. (4.5.10) in E94, which can be written aswhere

*p*

_{0}≈ 1000 hPa is a constant standard value. The division of

*s*by

*q*

_{d}means that the entropy in Eq. (12) is expressed “per unit mass of dry air.” The reference values of entropies disagree with the third law in E94, the consequence being that several terms are missing or are set to zero in Eq. (12). These missing terms may impact the specific entropy if

*q*

_{t}varies in space or time, since these missing terms must be multiplied by

*q*

_{d}= 1 −

*q*

_{t}to compute

*s*from

*s*/

*q*

_{d}given by Eq. (12). Moreover, since

*c*

_{pd}+

*r*

_{t}

*c*

_{l}depends on

*r*

_{t}, changes in

*s*and

*s*/

*q*

_{d}cannot be represented by

*θ*

_{e/E94}in Eq. (12) for varying values of

*r*

_{t}. This prevents

*θ*

_{e/E94}from being a true equivalent of the specific moist-air entropy,

*s*, for hurricanes where properties of saturated regions (large values of

*r*

_{t}) are to be compared with nonsaturated ones (small values of

*r*

_{t}).

*T*

_{0}= 273.15 K is a constant standard value. Again, it is an entropy expressed per unit mass of dry air and the specific heat

*r*

_{t}, twice preventing

*θ*

_{e/MPZ}from being a true equivalent of the specific moist entropy.

*r*

_{t}on the term

*θ*

_{es/E86}given by Eq. (2), yielding

### c. The dataset for Hurricane Dumile

The aim of the following sections is to show that varying the values for *r*_{υ}, has significant impact on the definition and the plot of the isentropic surfaces. To do this, all the moist-air equivalent potential temperatures and the entropies given by Eqs. (1)–(15) are computed for a series of 15 points selected arbitrarily and plotted in the west–east cross sections indicated on Figs. 1–4 for Hurricane Dumile.

At 0000 UTC 3 January 2013, Hurricane Dumile was located northwest of Réunion Island and east of Madagascar, near latitude 18.5°S and longitude 54.25°E. The pressure is used as a vertical coordinate and the black regions close to 1000 hPa represent the east coast of Madagascar on the left and the center of Dumile on the right.

The following figures are plotted for a 12-h forecast employing the French model Aire Limitée Adaptation Dynamique Développement International (ALADIN), with a resolution of about 8 km. The ALADIN–Réunion model is as described in Montroty et al. (2008) but with two main changes: (i) the diagnostic turbulence scheme (Cuxart et al. 2000) is based on a 1.5-order prognostic TKE equation and (ii) the shallow convection is based on Bechtold et al. (2001).

The ALADIN model is not perfect, owing to the uncertainties concerning the description of the dynamics (nonconservative schemes) and the thermodynamics (more or less accurate parameterizations of microphysics, turbulence, and clouds). Nonetheless, it is assumed that the quality of the ALADIN model is sufficient to describe, in broad terms, relevant spatial variations of the thermodynamic variables.

The two cross sections marked in Figs. 1 and 2 are plotted with the pseudoadiabatic potential temperature *H*_{l}. The use of *θ*_{e} recalled in section 2a, which will be compared in later sections.

The 15 points describe a moist-air steam cycle similar to the one described in P11 and inspired by the Carnot heat engine studied by E86, E88, E91, and Emanuel (2004). The basic thermodynamic conditions (*p*, *T*, *r*_{υ}, *q*_{l} = *q*_{i} = 0) of these points are listed in Table 1. The dry descent follows a path of almost constant relative humidity between 60% and 70% (points 11–15).

Pressure (hPa), temperature (K), water vapor mixing ratios (g kg^{−1}), relative humidity (relative to liquid water, %), and pseudoadiabatic potential temperature (K) for the 15 points indicated in Figs. 1–3.

The impacts of condensed water and of nonequilibrium terms will not be tested with these 15 points, chosen such that *T*_{l} = *T*_{i} = *T* everywhere and *H*_{l} = *H*_{i} = 1 and *q*_{l} = *q*_{i} = 0 for the just-saturated regions (points 5–7). However, these impacts are expected to be small in comparison with those induced by the large changes in *q*_{t} = *q*_{υ}, which vary from 3 to 17 g kg^{−1}.

The eyewall and the core of Hurricane Dumile are similar to the Figs. 16 and 12 (top) plotted in Hawkins and Imbembo (1976) for Hurricane Inez, where *θ*_{e} was likely computed as an equivalent for *θ*_{e/B73} are similar to those for *θ*_{e/B73} almost overlap.

## 3. Results

### a. Impact on plotting isentropic surfaces

The cross section indicated in Fig. 4 for *θ*_{s} exhibits great differences in comparison with Figs. 1 and 3 valid for *θ*_{e/B73} and with Fig. 2a for *θ*_{e/MPZ} in MPZ. To facilitate comparisons, the isolines of *θ*_{s} and *θ*_{e/B73} are plotted in Fig. 5 for the complete west–east cross section.

Generally speaking, the isentropes traced with *θ*_{s/M11} are smoother than the isolines plotted with *θ*_{e/B73}. This can be interpreted as a moderate impact of the total-water content *q*_{t} on the moist-air entropy variable *θ*_{s/M11} ≈ *θ*_{l} exp(6*q*_{t}), where the factor 6 is about two-thirds of that included in *θ*_{e/E94} ≈ *θ*_{l} exp(9*q*_{t}). The general aspect of the isolines plotted in Fig. 5 is pretty similar to those that can be observed on polar–equator sections depicted in zonal averages (not shown).

The vertical changes in potential temperatures in the core region are characterized by minimum values of *θ*_{s} close to the surface with increasing values above, up to 100 hPa. In contrast, maximum values are observed for *θ*_{e/B73} at 1000 hPa, with values decreasing up to 600 hPa and increasing values above.

At some distance from the center, *θ*_{s} exhibits midtropospheric minimum values in the layers from 900 to 750 hPa while the minimum values of *θ*_{e/B73}, and *θ*_{e/MPZ} are located higher, close to 700 hPa. The mixing in *θ*_{s} is greater in the boundary layer below 950 hPa, where the vertical gradient in *θ*_{s} is smaller than those for *θ*_{e/B73}.

The isentropes plotted with *θ*_{s} in Fig. 5 almost coincide with the isolines of *θ*_{e/B73} in the dry regions where *q*_{t} is small, namely in the high troposphere above the 200-hPa level and west of the core region (longitudes < 50°). The specific humidity *q*_{t} is larger east of the core region because the hurricane is not symmetric (not shown). The isentropes are thus different from the isolines of *θ*_{e/B73} at all levels below 300 hPa and for longitudes between 56° and 60°. Consequently, the vertical tilts of the isolines of *θ*_{s} and *θ*_{e/B73} can be very different locally, and especially within the eyewalls of the hurricane (close to or above the points 4–9, where largest values of *q*_{t} are simulated).

The isentropes computed with *θ*_{s} are therefore not compatible with the isolines of *θ*_{e/B73} and the larger the specific content *q*_{t} and the relative humidity *H*_{l} are in Fig. 2, the more the isentropes differ in Fig. 5. It can be seen in Fig. 5 that the ascending branches of the eyewalls of Hurricane Dumile cannot follow both the isentropes plotted with the third-law value *θ*_{s} and the isolines of *θ*_{e/B73}. This demonstrates that the different ways of computing the entropy, with either one of the formulations of *θ*_{e} or with *θ*_{s}, cannot be true simultaneously. Only one of them can be applied to plot moist-air isentropes and, thus, to achieve isentropic analyses for the atmosphere.

### b. Impacts on computing potential temperatures

All moist-air potential temperatures given by Eqs. (1)–(8) are plotted in Fig. 6 and listed in Table 2 (except *θ*_{s1/M11}).

Three groups of loops are shown in Fig. 6: (i) the three definitions for *θ*_{s}, (*θ*_{s})_{1}, and (*θ*_{s})_{2} are plotted in dark on the left; (ii) the three equivalent potential temperatures *θ*_{e/MPZ}, *θ*_{e/E94}, and *θ*_{e/B73} are plotted in light gray in the center; and (iii) the two saturated equivalent potential temperatures *θ*_{es/E86} and *θ*_{es/E94} are plotted in dark on the right.

Clearly, (*θ*_{s})_{1}, (*θ*_{s})_{2}, and *θ*_{s} remain close to each other with an accuracy of ±0.8 K for (*θ*_{s})_{1}, and with (*θ*_{s})_{2} almost overlapping *θ*_{s} (errors are smaller than 0.08 K). This is new proof that (*θ*_{s})_{1} and (*θ*_{s})_{2} are accurate increasing-order approximations for *θ*_{s}.

The equivalent formulations *θ*_{e/MPZ}, *θ*_{e/E94}, and *θ*_{e/B73} exhibit larger discrepancies, especially in the warm and moist ascent where differences of ±2.5 K are observed between 900 and 700 hPa. Moreover, the dry descent for the saturated versions *θ*_{es/E86} and *θ*_{es/E94} are 10 K warmer than that of other definitions of *θ*_{e}, making the loops for *θ*_{es/E86} and *θ*_{es/E94} (where *r*_{υ} is replaced by the larger saturating value *r*_{sw}) much narrower than the others.

The saturated values *θ*_{es/E86} and *θ*_{es/E94} remain close to one another. The differences of about 1 K in the dry descents can be understood by computing the factor *H*_{l} ≈ 60%, *R*_{υ} /*R*_{υ} /*c*_{pd} ≈ 0.46 and *r*_{sw} ≈ 15 g kg^{−1}, leading to (0.60)^{−0.46×0.015} ≈ 1.0035 and to 340 × 1.0035 ≈ 341.2 K.

Moderate differences of about ±3.5 K are observed between the five potential temperatures at elevated levels, at 450 hPa where ≈ 3 g kg^{−1} is small. The impact is much larger at low level where *r*_{υ} > 15 g kg^{−1}, with the ascent values of *θ*_{s} about 15 K colder than those for *θ*_{e} and *θ*_{es} at 950 hPa.

“Isentropic” surfaces or regions cannot, therefore, be the same when diagnosed using values of either *θ*_{s}, *θ*_{e}, or *θ*_{es}. The comparisons described in MPZ, which were based on analyses of altitude–*θ*_{e/MPZ} diagrams (MPZ’s Figs. 4 and 6–9), might be invalid since changes in the vertical direction may be of sign opposite the one for *θ*_{s}.

The solid disks plotted in Fig. 6 and the corresponding values in Table 1 show that ascents between the levels 800 and 450 hPa correspond to a small increase in *θ*_{s} of about +1 K, while they correspond to a large decrease in *θ*_{e/MPZ} of about −6 K. These differences in signs and in magnitudes between the properties of updrafts and downdrafts are not compatible with the differences in equivalent potential temperatures of about 5 K described in MPZ (p. 1867, right).

Another way to analyze the difference between *θ*_{s} and *θ*_{e/MPZ} or *θ*_{e/E94} is to consider the gap between descending and ascending values at 800 hPa: Δ*θ*_{s} ≈ 17 K and Δ*θ*_{e/MPZ} ≈ *θ*_{e/E94} ≈ 21 K. The difference is even larger for Δ*θ*_{e/B73} ≈ 24 K and is much smaller for Δ*θ*_{es/E86} ≈ 12 K.

The important consequence of these findings is that changes in moist-air entropy represented by either *θ*_{s} or any of the definitions of *θ*_{e} cannot be simultaneously positive or negative, because it would then be impossible to decide whether or not turbulence, convection, or radiation processes would increase or decrease the moist-air entropy in the atmosphere. Moreover, since isentropic processes or changes in entropy must be observable facts, it is impossible to consider that all definitions given by Eqs. (1)–(8) are equivalent; at most only one them can correspond to real atmospheric processes.

### c. Impacts on computing moist-air entropies

The issue of computing the relevant moist-air entropy is even harder than choosing one of the equivalent potential temperatures studied in the previous section, namely either *θ*_{s} or one of the versions of *θ*_{e} or *θ*_{es}. Since the aim of E86, E88, E91, Emanuel (2004), Pauluis et al. (2010), P11, or MPZ is to analyze meteorological properties in moist-air isentropic coordinates, a comparison of values of the moist-air entropy itself is needed. Let us therefore plot, in the temperature–entropy diagram of Fig. 7, the values of the six moist-air entropies considered in section 2b, with three of them listed in Table 3.

The moist-air entropies (J K^{−1} kg^{−1}) are *s*(*θ*_{s/M11}) computed by Eq. (10) with an offset of −6850 J K^{−1} kg^{−1}, (*s*/*q*_{d})(*θ*_{es/E86}) computed by Eq. (15) with an offset of −5650 J K^{−1} kg^{−1}, and (*s*/*q*_{d})(*θ*_{e/MPZ}) computed by Eq. (13) with an offset of +550 J K^{−1} kg^{−1}. The absolute temperature (K) is in the second column in order to follow the loops plotted in Fig. 7 with *T* as ordinates.

The differences between the loops (Carnot or steam cycles) are large. Some of the loops are very narrow, whereas others are wide and have a flared shape (a large gap between ascending and descending regions). Some of the loops are almost vertical (with a small change of less than ±20 J K^{−1} kg^{−1} in entropy between the surface and the upper air), whereas others exhibit a pronounced tilt (with a large decrease of 400 J K^{−1} kg^{−1} between the warm, moist regions and the cold, dry ones). The third-law entropy *s*(*θ*_{s}) increases between points 2 and 10, whereas all the other formulations lead to a decrease between these two points.

Therefore, since the moist-air entropy is a state function, and since isentropic processes or changes in entropy must be observable facts, it is impossible to consider that all definitions given by Eqs. (10)–(15) are equivalent: at most only one of them can correspond to real atmospheric processes. Otherwise, it would be possible to modify the second law by imposing arbitrarily an increase or a decrease in entropy for an ascending parcel of moist air.

### d. Impact on isentropic streamfunctions

*θ*

_{e/E94};

*z*) is studied in Pauluis and Mrowiec (2013) and MPZ and the alternative third-law value Ψ(

*θ*

_{s/M11};

*z*) is tested in this section. The isentropic streamfunction defined in MPZ can be written aswhere

*θ*

_{e/s}stands either for

*θ*

_{e/E94}or

*θ*

_{s/M11}, respectively. The isentropic mass flux 〈

*ρw*〉(

*r*,

*θ*,

*z*) is computed according toThe first rectangular function

*D*(

*θ*) is different from zero and equal to 1/Δ

*θ*= 1 K

^{−1}for |

*θ*−

*θ*′(

*x*,

*y*,

*z*)| < Δ

*θ*/2 and for

*θ*=

*θ*

_{e/E94}or

*θ*

_{s/M11}varying from

*θ*

_{min}= 320 K to

*θ*

_{max}= 360 K in steps of Δ

*θ*(41 bin values). The second rectangular function

*D*(

*r*) is different from zero and equal to 1/Δ

*r*= 1/15 km

^{−1}for |

*r*−

*r*′(

*x*,

*y*,

*z*)| < Δ

*r*/2 and for

*r*varying from 7.5 to 487.5 km in steps of Δ

*r*(33 bin values). The radial distances

*r*and

*r*′(

*x*,

*y*,

*z*) are computed with respect to the center of Hurricane Dumile and the average value

*z*) is computed for each vertical level for

*r*< 495 km.

The diagram plotted in Fig. 8 for Ψ(*θ*_{e/E94}; *z*) and for Hurricane Dumile is similar to the one plotted in Fig. 1b in Pauluis and Mrowiec (2013) and in Fig. 7a in MPZ for an idealized hurricane. The locations of the minimum values of the streamlines correspond to zero mean vertical mass flux and they are roughly aligned with the mean values of *z*) (solid line).

The arrows give some hints on the average anticlockwise circulations. Values of *θ*_{e/E94} are conserved only for the descending branch between 4 and 2 km and for values close to 338 K. Values of *θ*_{e/E94} decrease in the ascending branch below the freezing level at 5 km and they increase above this level up to the top of the troposphere.

The diagram plotted in Fig. 9 for the third-law streamfunction Ψ(*θ*_{s/M11}; *z*) exhibits large differences with Fig. 8, in particular below the freezing level and for the larger values of *q*_{t}.

Values of *θ*_{s/M11} are conserved for the ascending branch between 1.5 and 4 km and for values close to 335 K. As in Fig. 8, the locations of the minimum values of the streamlines are roughly aligned with the mean values of *z*) (solid line). Moreover, a dashed straight line has been added on Fig. 9 to show an intriguing linear organization from 2 to 12 km, which follows the location of all minimum values of Ψ(*θ*_{s/M11}; *z*). The upward and downward motions in Hurricane Dumile are located below and above this straight line and they correspond to monotonic increases and decreases in *θ*_{s/M11}, respectively. This simple linear organization is not observed for Ψ(*θ*_{e/E94}; *z*) in Fig. 8.

The underlying hypothesis justifying the study of the isentropic streamfunction is the fact that adiabatic and reversible processes lead to isentropic states, namely that any change in the moist-air entropy must correspond to diabatic sources or sinks. Accordingly, the diabatic processes of mixing and entrainment are considered in MPZ to explain 1) the decrease in *θ*_{e/E94} by following the streamlines of Ψ(*θ*_{e/E94}; *z*) from the surface up to the melting level and 2) the increase in *θ*_{e/E94} above the melting level. Clearly, if a similar increase of *θ*_{s/M11} is observed in Fig. 9 in dry regions above the freezing level, the decrease of *θ*_{s/M11} with *z* below the freezing level does not exist, except for a small decrease close to the ground (below 0.5 km) which is likely due to the impact of surface heating.

Comparing Figs. 8 and 9 reveals that *θ*_{e/E94} ≈ *θ*_{l} exp(9*q*_{t}) overestimates the impact of humidity in comparison with the third-law value *θ*_{s/M11} ≈ *θ*_{l} exp(6*q*_{t}). The ratio 2/3 for the coefficients in factor of *q*_{t} in the exponential functions reveals new patterns and it is thus important to use the absolute entropy *s*(*θ*_{s/M11}) when computing the isentropic streamfunction given by Eqs. (16) and (17).

The differences between the potential temperatures *θ*_{e/E94} and *θ*_{s/M11} described by studying the outputs of the ALADIN model may depend on the various approximations that are made in all numerical models to describe moist-air processes and thermodynamic equations.

However, it is indicated in M11 that similarly large differences between *θ*_{e/E94} and *θ*_{s/M11} are observed for several vertical profiles of potential temperatures computed from the radial flights in stratocumulus for the First International Satellite Cloud Climatology Project (ISCCP) Regional Experiment (FIRE-I), Atlantic Stratocumulus Transition Experiment (ASTEX), Eastern Pacific Investigation of Climate Processes in the Coupled Ocean–Atmosphere System (EPIC), Second Dynamics and Chemistry of the Marine Stratocumulus field study (DYCOMS II) campaigns.

Moreover, it is possible to validate the previous results derived from numerical outputs of ALADIN for Hurricane Dumile by comparing these results with several observed soundings into Hurricane Earl (Wang et al. 2015). Figures 10a and 10b show that values of *θ*_{s/M11} near the ground are colder by about 11 K than those at altitude at 400 hPa or 10-km height, whereas values of *θ*_{e/E94} are warmer by about 10 K. In addition, the height of minimum value and the dispersion of the profiles are both smaller with *θ*_{s/M11} than with *θ*_{e/E94.}

The shape of individual and average profiles in Figs. 10a and 10b are thus similar to the descending parts of the streamfunctions in Figs. 8 and 9 and of the loops of potential temperatures in Fig. 6. This confirms that the numerical outputs from the ALADIN model is a relevant laboratory to investigate isentropic processes and that large differences do exist between the potential temperatures *θ*_{e/E94} and *θ*_{s/M11}.

### e. Impacts on heat input and work functions

Many alternative definitions for the heat input or work functions have been considered in E86, E91, E94, Rennó and Ingersoll (1996), Pauluis and Held (2002a,b), and Goody (2003) and have been applied to different Carnot or steam cycles. The aim of this section is to use the outputs of Hurricane Dumile to compare the numerical values of the heat input and work functions computed for several of these definitions.

*e*

_{i}is the specific internal energy,

*h*the specific enthalpy,

*ρ*the density, and

*μ*

_{k}=

*h*

_{k}−

*Ts*

_{k}the specific Gibbs’ function for

*k*= 0, 1 2, 3 denoting dry-air, water vapor, liquid water, and ice, respectively.

*W*= −

*ρ*

^{−1}

*dp*can then be computed by integrating Eq. (18) along a closed circuit. Since the integral of the total differential of enthalpy

*W*can be written asThe last integral is derived from the last integral in Eq. (18) for the steam cycle considered in Figs. 1 to 4, where

*q*

_{l}=

*q*

_{i}= 0,

*q*

_{υ}=

*q*

_{t}, and

*q*

_{d}+

*q*

_{t}= 1, leading to

*dq*

_{l}=

*dq*

_{i}= 0,

*dq*

_{d}= −

*dq*

_{υ}= −

*dq*

_{t}, and thus to

*μ*

_{d}

*dq*

_{d}+

*μ*

_{υ}

*dq*

_{υ}= (

*μ*

_{υ}−

*μ*

_{d})

*dq*

_{t}. The same equation for

*W*holds true if liquid or solid water exist, provided that changes of phases are reversible, which implies (

*μ*

_{υ}−

*μ*

_{l})

*dq*

_{l}= (

*μ*

_{υ}−

*μ*

_{i})

*dq*

_{i}= 0.

*μ*

_{k}=

*h*

_{k}−

*Ts*

_{k}, the Gibbs Eq. (18) can be rewritten asThis relation clearly shows that, since

*dh*

_{k}=

*c*

_{pk}

*dT*(for all species),

*ds*

_{k}=

*c*

_{pk}

*dT*/

*T*−

*R*

_{k}

*dp*

_{k}/

*p*

_{k}(for gases) and

*ds*

_{k}=

*c*

_{pk}

*dT*/

*T*(for liquids or solids), the reference values of enthalpies and entropies have no impact and may remain undetermined if the aim is to compute −

*dp*/ρ.

Accordingly, the reference enthalpies *h*_{v0} and *h*_{d0} can be discarded in the computations of the enthalpies of water vapor and dry air contained in the last integral of Eq. (19), which depends on (*μ*_{υ} − *μ*_{d}) = (*h*_{υ} − *h*_{d}) − *T*(*s*_{υ} − *s*_{d}) . This property holds true because *h*_{v0} − *h*_{d0}) *dq*_{t} = (*h*_{v0} − *h*_{d0}) *dq*_{t} = 0 since *dq*_{t} ≡ 0 for a closed loop.

In contrast, the reference entropies must be taken into account in the computations of *T*(*s*_{d0} − *s*_{v0})*dq*_{t} = (*s*_{d0} − *s*_{v0})*Tdq*_{t}, because *T*(*q*_{t})*dq*_{t} ≠ 0 if the area of the loop is different from zero in a *T*–*q*_{t} diagram. However, the impact of the reference entropies *s*_{d0} and *s*_{v0} on the last integral in Eq. (19) must cancel out with the impact on the integral *T*(*s*) *ds*. This means that a change in the reference entropies must modify the values of the two integrals on the right-hand side of Eq. (19), although their sum, which is equal to *W* = −*dp*/*ρ*, must not depend on *s*_{d0} or *s*_{v0}.

*h*and

*s*replaced by

*h*/

*q*

_{d}and

*s*/

*q*

_{d}. Furthermore, the reference entropies for dry air and liquid water are assumed to vanish at

*T*

_{0}= 0°C in P11, leading toThe work functions

*W*or

*W*

_{(d)}are approximated in the older studies on Carnot cycles by considering only the first integral on the right-hand sides of Eqs. (19) or (21)—see E91 (Eqs. 4 and 10) and Rennó and Ingersoll [1996, Eq. (3)]—yieldingIn Eqs. (23) and (24)

*W*

_{H}and

*W*

_{H(d)}are commonly called “heat input.” They are equal to the areas enclosed by the loops in the temperature–entropy diagram of Fig. 7 and they are independent of the global offset chosen for each loop. Values of the heat input

*W*

_{H}and

*W*

_{H(d)}listed in Table 4 are calculated using the simple but robust first-order Riemann method.

The heat inputs *W*_{H} [or *W*_{H(d)}] (J kg^{−1}) computed for the six moist-air entropies considered in section 2b. Here *W*_{H} [or *W*_{H(d)}] represents the areas of the six loops shown in Fig. 7 (Riemann integrals). The wind scale *V* = (2*W*_{H})^{1/2} (m s^{−1}).

Large factors, of 0.5 and 4, are observed between *W*_{H} ≈ 870 J kg^{−1}, computed with the third-law value *s*(*θ*_{s}), and *W*_{H(d)} ≈ 442 J kg^{−1} or *W*_{H(d)} ≈ 3405 J kg^{−1}, computed with *s*(*θ*_{es/E86})/*q*_{d} or *s*(*θ*_{e/E94})/*q*_{d}, respectively. Furthermore, the impact of *c*_{pd}, leading to a factor of more than 6 for *W*_{H(d)} computed with the same potential temperature *θ*_{es/E86} (442 vs 2758 J kg^{−1}).

The impact on *W*_{H} of the definition per unit of dry air versus the definition per unit of moist air (specific value) can be evaluated by comparing *s* and *s*/*q*_{d} for the same third-law value *θ*_{s/M11}. The impact 1528 − 870 = 658 J kg^{−1} is large, leading to an increase of more than 75% for the definition per unit of dry air.

The impact of the choice of the different formulations for the moist-air entropy can be evaluated differently, by computing the wind scale *W*_{H} = *V* ^{2}/2, with *V* becoming a crude proxy for the surface wind that a perfect Carnot engine might produce. The last line in Table 4 shows that *V* would vary from about 30 to more than 80 m s^{−1}.

Since temperatures can vary by 30 K along the Carnot cycle, and for an accuracy of about 1 J K^{−1} kg^{−1} for entropies, errors in computations of *W*_{H} and *W*_{H(d)} are on the order of 30 J kg^{−1}. This means that the observed differences, of the order of hundreds or thousands of joules per kilogram, are significant and, since the moist-air entropy is a state variable, the third-law formulation computed with *θ*_{s/M11} must be used to compute both *W*_{H} and *W*_{H(d)}.

Independent computation of the integrals on the right-hand sides of Eqs. (19) and (21) is thus impossible without knowledge of the third-law reference entropies, because only the sum of these integrals is independent of these reference values. It is thus necessary to compute the full work function and to switch from a Carnot cycle to a steam cycle, for which changes in water vapor have a large impact via the last integrals of Eqs. (19) and (21).

It is assumed in E88 that part of the energy available from the steam cycle is actually used to lift water and is not available for the generation of kinetic energy. However, insofar as Eqs. (19) and (21) are considered without including these possible effects, the values of all integrals appearing for *W* and *W*_{(d)} are listed in Table 5. The left- and right-hand sides of either *W* or *W*_{(d)} are equal within an accuracy of about 40 J kg^{−1}, which is better than 10% for *W* and 6% for *W*_{(d)}. This accuracy can be improved by using a Simpson second-order method, leading to 2.5% for *W* and 0.5% for *W*_{(d)} (not shown). This confirms the theoretical methods for establishing both Eqs. (19) and (21). In particular, it proves that the reference enthalpies and entropies of dry air and water vapor do not have an impact on *W* and *W*_{(d)}.

It seems somewhat surprising that *W*_{(d)} can be (703 − 455)/455 = 55% higher than *W*. However, this large difference is not a numerical artifact; it merely comes from the replacement of *ρ* by *ρ*_{d} in the computation of *W*_{(d)}. To demonstrate this result, values of *ρ* and *ρ*_{d} are shown in Fig. 11 for the 15 points of the steam cycle of Hurricane Dumile. To facilitate visual comparisons, the terms *F*(*p*) = 0.7(*p* − 950)/450 and *F*(*p*) − 0.014 have been added to 1/*ρ* and 1/*ρ*_{d}, respectively. These terms are the same for a given pressure and they do not affect the computation of the differences between upward and downward values.

Although the replacement of 1/*ρ* by 1/*ρ*_{d} leads to a relative change of less than 2% if *q*_{υ} < 20 g kg^{−1}, the differences between upward and downward values are much larger: (1/*ρ*_{d})_{up} − (1/*ρ*_{d})_{down} is more than 40% larger than (1/*ρ*_{d})_{up} − (1/*ρ*_{d})_{down} for all levels between 900 and 600 hPa. Therefore, the area of the loop *W*_{(d)} is logically about 45% larger than *W*, and the mixing-ratio view expressed per unit mass of dry air leads to a value *W*_{(d)} that largely overestimates the work function *W* defined in thermodynamics with specific values. The consequence is the need to consider *W* for studying steam cycles, since the use of *W*_{(d)} would lead to large errors and to erroneous physical meanings if *W*_{(d)} was to be converted into potential or kinetic energy.

The aim of the present paper is not to provide the most general formulation for the work functions *W* given by Eq. (19) or (21). Indeed, *W* may be subject to some uncertainty with, for example, some additional terms that are added to Eq. (21) in E88. These additional terms are intended to take into account possible irreversible entropy sources ranging from the fallout of condensed phase water to irreversible mixing across water concentration gradients to freezing of supercooled water. Moreover, it is assumed in E88 that part of the energy available from the steam cycle is actually used to lift water and is not available for the generation of kinetic energy.

Differently, it is shown in the present paper that once a formulation for *W* is chosen, with or without the additional terms of E88, and provided that it is indeed derived from the Gibbs equation, the value of the work function may depend on the reference entropies of dry air and water vapor. And in this case, it is essential to use the values derived from the third law of thermodynamics.

### f. Impacts on the efficiency of steam cycles

It was shown in the previous section that it is impossible to compute the integrals on the right-hand sides of Eqs. (19) and (21) independently. Nevertheless, calculations of this kind are often used in atmospheric science to evaluate the efficiency factor of a thermal engine.

The efficiency *E* of a tropical cyclone considered as a thermal engine was defined in Palmén and Riehl (1957) as “the ratio of the mechanical energy produced to the heat released” and *E* was evaluated to be about 2%–3%.

*W*=

*Q*

_{1}=

*Q*

_{2}and

*W*

_{(d)}=

*Q*

_{1(d)}=

*Q*

_{2(d)}, whereandThe “energy input” is defined in P11 as

*Q*

_{in}=

*δQ*

^{+}, where

*δQ*

^{+}= max(0,

*δQ*) is different from zero only for positive values of

*δQ*. The mechanical efficiency is then defined and computed in P11 by the ratiowhere all quantities are expressed per unit mass of dry air. From Eq. (22) in P11, this ratio can be approximated by

*χ*

_{S}≈ (

*T*

_{in}−

*T*

_{out})/1540 K, where

*T*

_{in}is the highest and

*T*

_{out}the lowest temperature of the cycle and where (

*R*

_{d}

*L*

_{υ})/(

*R*

_{υ}

*c*

_{pd}) ≈ 1540 K. The steam cycle considered for Hurricane Dumile corresponds to

*T*

_{in}≈ 296 K and

*T*

_{out}≈ 265 K, leading to

*χ*

_{S}≈ 0.020.

Table 6 shows that *W* = *Q*_{1} ≈ *Q*_{2} and *W*_{(d)} = *Q*_{1(d)} ≈ *Q*_{2(d)} with the same accuracy of about 40 J kg^{−1} previously noted for the Riemann method. The differences *Q*_{1} ≠ *Q*_{1(d)} and *Q*_{2} ≠ *Q*_{2(d)} are due to the use of mixing ratios and of quantities expressed per unit mass of dry air in P11. This means that the way the work function is defined and computed may largely impact the diagnosis of the external heating required with any transformation of moist air and may modify the numerator of the efficiency factor given by Eq. (29).

External heating (*Q*), energy inputs (*Q*_{in}), and mechanical efficiency (*χ*). All integrals are computed by using the first-order Riemann method.

Similarly, the way *δQ*^{+} and *Q*_{in} are computed may impact the denominator of Eq. (29). The computations of *δQ*^{+} for *Q*_{1}, *Q*_{2}, *Q*_{1(d)}, or *Q*_{2(d)} are especially uneasy because they rely on accurate definitions of all the enthalpies (*h* of *h*/*q*_{d}), the entropies (*s* or *s*/*q*_{d}), and the Gibbs functions (*μ*_{υ}, *μ*_{d}, or *g*_{υ}). Moreover, whereas *dh* = 0 for a closed cycle, the integral of *dh* is different from zero if it is computed only for those points where *δQ*_{1} > 0.

*h*/

*q*

_{d}), entropy (

*s*/

*q*

_{d}), and Gibbs functions (

*g*

_{υ}) appearing in Eqs. (27) and (28) are those defined in P11. The specific entropy involved in Eqs. (26) is the third-law value defined in M11 and recalled in Eq. (10). The specific enthalpy in Eq. (25) is the one defined in Marquet (2015a,c):where

*L*

_{h}(

*T*

_{r}) ≈ 2603 kJ kg

^{−1}and

*h*

_{ref}≈ 256 kJ kg

^{−1}. The Gibbs functions

*μ*

_{υ}and

*μ*

_{d}are computed by

*h*

_{υ}−

*Ts*

_{υ}and

*h*

_{d}−

*Ts*

_{d}, where

*h*

_{d}(

*T*

_{r}) ≈ 530 kJ kg

^{−1}and

*h*

_{υ}(

*T*

_{r}) ≈ 3133 kJ kg

^{−1}.

The latent heat *L*_{h}(*T*) defined in Eq. (31) depends on the values of *L*_{h}(*T*_{r}), in the same way as *L*_{υ} depends on *L*_{υ}(*T*_{r}) ≈ 2501 kJ kg^{−1}, and in the same way as the entropy potential temperature *θ*_{s/M11} defined in Eq. (6) depends on the third-law values Λ_{r}, (*s*_{υ})_{r}, and (*s*_{d})_{r}.

However, *L*_{h}(*T*) is independent of the reference temperature *T*_{r} because *L*_{h}(*T*_{r}) + (*c*_{pυ} − *c*_{pd})*T*_{r} is a constant, in the same way as *L*_{υ} and *L*_{s} are independent of the value *T*_{r} used to compute them, and in the same way as *θ*_{s/M11} does not depend on *T*_{r}, *p*_{r}, and *r*_{r} = *r _{sw}*(

*T*,

_{r}*p*) = ε

_{r}*e*/(

_{r}*p*−

_{r}*e*) provided that

_{r}*e*=

_{r}*e*(

_{sw}*T*).

_{r}The efficiency factors shown in the last column of Table 6 are almost the same for *Q*_{1} and *Q*_{2} on the one hand and for *Q*_{1(d)} and *Q*_{2(d)} on the other hand. However, the values, respectively close to 2.2% and 3.4%, are not compatible. Although these two efficiency values are close to those derived in Palmén and Riehl (1957), only the former value is close to the approximate value *χ*_{S} ≈ 0.020 derived in P11. The explanation for this important difference seems to come from the energy input terms *Q*_{in/1} and *Q*_{in/2}, which are only 1% smaller than *Q*_{in/1(d)} and *Q*_{in/2(d)}, whereas the full work functions *Q*_{1(d)} and *Q*_{2(d)} are 50% larger than *Q*_{1} and *Q*_{2}.

The definitions of Palmén and Riehl (1957) or P11 for the efficiency factor may not be the most relevant ones. It is thus probably necessary either to change these definitions in order to be independent in the choice of the reference values of entropies and enthalpies or to use the third-law potential temperature, with all specific quantities expressed per unit of moist air, in order to comply with the recommendations of open-system thermodynamics.

## 4. Discussion and conclusions

Isentropic analysis is probably a powerful tool for investigating moist-air energetics by plotting moist-air isentropes or by computing isentropic mass fluxes. However, the quality and the realism of such an analysis rely on a clear definition of the moist-air entropy.

It has been shown, here, that the way the potential temperatures *θ*_{s}, *θ*_{e}, or *θ*_{es} are defined as “equivalents” of the moist-air entropy significantly impacts the computations and plots of isentropic surfaces, making “isentropic” analyses similar to those published in MPZ or P11 uncertain.

It has been pointed out that the streamfunction, the heat input, and the efficiency factor may be largely modified by the way in which the entropy and the enthalpy are defined: *s* or *s*/*q*_{d}, *h* or *h*/*q*_{d}; modified reference values for entropies and enthalpies; the choice of *θ*_{s}, *θ*_{e}, or *θ*_{es}, the choice of *c*_{pd} or

The issue associated with the “per unit mass of dry air” view can be understood as follows. If isentropic processes are defined as in E94, MPZ, and P11, with constant values of *s*/*q*_{d}, the definitions of the geopotential, the wind components, or the kinetic energy one should modify accordingly by plotting, for instance, *gz*/*q*_{d}, *u*/*q*_{d}, *υ*/*q*_{d}, or (*u*^{2} + *υ*^{2})/(2*q*_{d}). These definitions are unusual in the thermodynamics of open systems. Moreover, if *s*/*q*_{d} could be defined within a global constant *C*, the specific value *s* would depend on *q*_{d}*C*, which varies with *q*_{d} and renders the integral *C* is an unknown term. This is not realistic.

The third law of thermodynamics, which can be based on calorimetric, quantum, or statistical physics, can now be considered as well established. It provides a fortunate opportunity to compute, describe, analyze, and possibly understand several features of atmospheric energetics in a new way.

It is important to clearly differentiate the wish to define conserved variables like *θ*_{l}, *θ*_{e}, or *θ*_{es}, on the one hand, from the need to define the specific moist-air entropy by computing the third-law value *θ*_{s}, on the other hand.

These two aspects are complementary because a given process which might correspond to the law of conservation of *θ*_{e}, for instance, may correspond to a change of entropy that can be precisely computed by using the third-law value *s*(*θ*_{s}). Conversely, a nontrivial isentropic process where *θ*_{s} is conserved, but where both *θ*_{l} and *q*_{t} vary with time and space, might not correspond to a conservation of any of *θ*_{e}, or *θ*_{es}.

Moist-air isentropic surfaces are not subject to uncertainty in the natural world and it is likely that the third-law definitions *θ*_{s} and *s*(*θ*_{s}) given by Eqs. (6) and (10) are the most relevant, since they are based on general thermodynamic principles and use specific values expressed per unit mass of moist air, as with all other variables in fluid dynamics.

I thank the two reviewers for their comments, which helped to improve the manuscript. The data for Hurricane Earl were provided by NCAR/EOL under the sponsorship of the National Science Foundation (https://data.eol.ucar.edu/).

# APPENDIX A

## The History of the Third Law of Thermodynamics

It is surprising that so many potential temperatures may serve to compute the moist-air entropy at the same time. Since the moist-air entropy is a state function, the status of the isentropic or diabatic evolution of a given parcel of moist air is an observable feature and, in the words of Richardson (1922, p. 157), “approximations are not here permissible.” The only way to solve this problem is to rely on general thermodynamics where the aim of the third law is precisely to provide an absolute and nonambiguous definition of the entropy state function. According to reviews of the old treatises and papers of thermodynamics (Lewis and Randall 1923; Wilks 1961; Barkan 1999; Coffey 2006; Klimenko 2012), there are three ways to express the third law.

*H*and Δ

*S*are the enthalpy and entropy of reaction and

*T*is the temperature. Negative values of Δ

*G*would make the reaction proceed spontaneously because, according to the second law of thermodynamics, the total entropy of the system and its surroundings, (Δ

*S*)

_{tot}, would increase.

Gibbs (1878a, 150–152) explicitly introduced two constants of integration appearing in the computation of the specific energy, entropy, free energy, and free enthalpy for an ideal gas. He then arrived at the definitions of the entropy, free energy, and free enthalpy of a mixture of ideal gases [Eqs. (278), (279) and (293); p. 156 and 163]. These definitions depend on the constants of integration, which, a priori, are different for each individual gas and are multiplied by the individual variable concentrations. Moreover, since Δ*S* is multiplied by the variable temperature in the computation of Δ*G* given by Eq. (A1), the free enthalpy computed in Eq. (293) of Gibbs (1878a) must depend on the constants of integration for the specific entropies.

*H*, Δ

*S*, and Δ

*G*was next considered by Le Chatelier (1888), Lewis (1899), Richards (1902), van’t Hoff (1904), and Haber (1905) with an alternative version of Eq. (A1). They all focused on the free energy

*A*= Δ

*G*on both sides of Eq. (A1) by using the Helmholtz equationswhere

*A*represents the maximum available work (“maximaler Arbeit”). Comparison of Eqs. (A1) and (A2) shows that the derivative at constant pressure

*dA*/

*dT*is equal to −Δ

*S*(the opposite of the entropy of reaction).

*H*(

*T*) for different absolute temperatures from laboratory experiments, and then of finding the value of

*A*(

*T*)/

*T*by an integration of Eq. (A2) between

*T*

_{0}and

*T*, leading to

The knowledge of *A*(*T*) from Eq. (A3) should allow the spontaneity of any chemical reaction to be predicted, for a given temperature *T*; it depends on the sign of *A*(*T*) = Δ*G*. However, the constant of integration *A*(*T*_{0})/*T*_{0} is multiplied by *T* and influences the value of *A*(*T*) and, thus, possibly the sign of Δ*G*(*T*). This cannot be true because it would bring some arbitrariness into the spontaneity of the chemical reaction, making it possible to modify at will the behavior of chemical reactions like H_{2} + O_{2} = H_{2}O in the atmosphere by choosing any arbitrary value for the constant of integration *A*(*T*_{0})/*T*_{0} for the temperature *T*_{0}.

*A*, leading to

*A*(

*T*

_{0})/

*T*

_{0}→ 0 as

*T*

_{0}→ 0 K in Eq. (A3). This is the third law of thermodynamics in its original formulation suggested by Nernst.

Nernst tested the consequences of his theorem by computing the theoretical values of the affinities or the reaction rates of several chemical reactions. The good agreement he showed with experimental results supported his theorem. The hope of Nernst was then to demonstrate his theorem starting from the observed anomalous behavior of specific heats, which tend to zero at low temperature and with *C*(*T*)/*T* remaining finite at absolute zero, which implied that the entropy itself defined by *ds* = [*C*(*T*)/*T*]*dT* remained finite at 0 K.

In response to Einstein’s objections to the heat theorem, expressed during the First Solvay Congress in Brussels in the fall of 1911, Nernst (1912) suggested the alternative principle of unattainability of *T* = 0 (“Unerreichbarkeit des absoluten Nullpunktes”), namely that absolute zero temperature cannot be reached in a finite time interval and in a finite number of steps. A proof for this second way to express the third law of thermodynamics is given in the recent paper Masanes and Oppenheim (2017).

Planck made the heat theorem more general in his third (1911) and fifth (1917) German editions of his *Treatise on Thermodynamics* (Planck 1917). The third formulation of the third law expressed by Planck was clearly based on the entropy and was written in two parts:

- “The gist of the (heat) theorem is contained in the statement that, as the temperature diminishes indefinitely the entropy of a chemical homogeneous body of finite density approaches indefinitely near to a definite value, which is independent of the pressure, the state of aggregation and of the special chemical modification” (p. 273).
- “…without loss of generality, we may write [this definite value] lim(
*T*→ 0) entropy = 0. We can now in this sense speak of an absolute value of the entropy” (p. 274).

Planck wrote in the preface of the fifth edition in 1917, “The theorem in its extended form has in the interval (since 1911) received abundant confirmation and may now be regarded as well established.” However, it was discovered after 1917 that a residual entropy may exist for anomalous species like Ice-lh at 0 K owing to quantum effects explained by Pauling (1935) and Nagle (1966) (proton disorder and remaining randomness of hydrogen bonds at 0 K). For these reasons, Planck’s formulation must now be amended and applied to the “most stable form of pure crystalline solid substances.”

*T*

_{0}and pressure

*p*

_{0}can be written asThe integral of the piecewise function

*c*

_{p}(

*T*,

*p*

_{0}) is computed from 0 K to

*T*

_{0}for all solid(s), liquid, and gaseous states. The sum over

*k*is extended to all transitions of phases occurring at

*T*

_{k}with the latent heats

*L*(

*T*

_{k}).

The validity of the third law is now considered as established by the relevance of the computations of Δ*G*, which provide accurate predictions of the constants of chemical reactions and of the spontaneity of these reactions with respect to temperature. This is why thermochemical tables include absolute values for entropy *S*, whereas only relative standard enthalpies of reactions

According to the results already available in Kelley (1932) and in more recent tables, the absolute values of entropies determined from the calorimetric method (A5) and with the third law are in good agreement with the other available methods: 1) by using the statistical-physics method—that is, by computing the translational partition function and the corresponding Sackur–Tetrode equation valid for monoatomic gases only, and thus for Argon in the atmosphere (Grimus 2013); 2) by computing the partition functions valid for more complex molecules (like O_{2}, N_{2}, H_{2}O, and CO_{2}), including the impact of possible rotational, vibrational, or anharmonic vibration corrections or electronic states of molecules (Gordon and Barnes 1932; Gordon 1934, 1935, and subsequent works); 3) by analyzing residual-ray frequencies and infrared absorption spectra of crystals or by using calculations from spectroscopic data for gases; and 4) by measuring both Δ*G* and Δ*H*, leading from Eq. (A1) to Δ*S* = (Δ*H* − Δ*G*)/*T*.

It is important to note that the third law corresponds to a similar hypothesis, which is implicitly assumed in Boltzmann’s entropy written as *S* = *S*_{0} + *k* ln(*W*), namely that *S*_{0} is a constant, which is usually set to zero. It is clearly explained in Schrödinger (1989, lectures delivered in 1944 at Dublin) that the second part of the third-law assumption suggested by Planck (1917), to write *S*_{0} = 0 without loss of generality, must not be regarded as the essential thing. This would create confusion and would draw attention away from the point really at issue, namely that *S*_{0} is a universal constant that has no physical meaning only if it is independent of any internal physical parameters of any species, no matter what the value of *S*_{0} might be.

The entropies of the five most frequent atmospheric species are listed in Table A1 for the values computed with either the third-law quantum statistical-physics (theoretical) method or with the third-law calorimetric (experimental) method. Good agreement can be observed between theoretical and experimental values. The larger discrepancy observed for H_{2}O and GR96 (45.132 vs 44.31) is not observed for MG15 (45.2). The explanation for the larger value observed for O_{2} and MG15 (49.7 vs 49.031) is given in appendix B.

The first three columns list the standard entropies for five atmospheric gases at 298.15 K and *p*_{0} = 1000 hPa published in Chase (1998, C98) and Gokcen and Reddy (1996, GR96). The statistical-physics method is based on *S* = *k*ln(*W*) + *S*_{0} with *S*_{0} = 0. The calorimetric method is based on Eq. (A5) with the third law *S*(0 K, *p*_{0}) = 0. The fourth column corresponds to values computed with *S*(0 K, *p*_{0}) = 0 and the calorimetric datasets described in Marquet (2015c). Units are cal K^{−1} mol^{−1} throughout.

# APPENDIX B

## Applications of the Third Law to Atmospheric Studies

The apparent paradox is that, while there is no chemical reaction, the third law discovered in thermochemistry must be applied if the moist-air entropy is to be computed in the atmosphere.

*R*

_{d}and

*R*

_{υ}and specific heats

*c*

_{pd}and

*c*

_{pυ}are assumed to be constant in the range of atmospheric temperatures (from 170 to 340 K), the entropies for dry air, water vapor, and moist-air mixing can be written aswhere

*p*=

*p*

_{d}+

*e*,

*q*

_{d}= 1 −

*q*

_{υ}, and

*c*

_{p}=

*q*

_{d}

*c*

_{pd}+

*q*

_{υ}

*c*

_{pυ}and with the use of the partial pressures

*p*,

*e*,

*p*

_{d0}, and

*e*

_{0}, which automatically takes the entropy of mixing of the two gases into account.

The impact of the last bracketed terms in Eq. (B4) is similar to the impact of the last bracketed terms in Eq. (A3). These terms are the product of the constants of integration *A*(*T*_{0})/*T*_{0}, *s*_{d0}, and *s*_{v0} by the variable terms *T*(*x*, *y*, *z*), *q*_{d}(*x*, *y*, *z*), and *q*_{υ}(*x*, *y*, *z*), respectively.

Therefore, the same consequences as described for thermochemistry hold true for the atmosphere: in order to determine if a process is diabatic or isentropic, or for making isentropic analyses, it is necessary to compute local values of moist-air entropy with Eq. (B4) and, in particular, *q*_{d}(*x*, *y*, *z*)*s*_{d0} and *q*_{υ}(*x*, *y*, *z*)*s*_{v0}. Consequently, the two reference values *s*_{d0} and *s*_{v0} must be known and the use of the third-law values would be in full agreement with general thermodynamics recommendations.

Bjerknes (2009) may well have been the first to imagine in 1904 to using entropy in weather prediction, as a prognostic variable in one of his seven basic equations. It is recalled in Marquet (2016a) that the consequences of the third law on atmospheric thermodynamics was described soon after in Richardson (1922, 158–160), who was already aware of the problem of the constant of integration in the entropy. Richardson suggested that “the most natural way of reckoning the entropy of the water substance would be to take it as zero at the absolute zero of temperature.” This corresponds to the use of the third-law values for *s*_{d0} and *s*_{v0}.

However, Richardson was not able to continue accurate computations of the moist-air entropy in 1922 because values of *c*_{p}(*T*) were not available at that time for an absolute temperature varying from near zero to 350 K. These measurements were made later on for all atmospheric species by using the magnetic refrigeration method to attain extremely low temperatures far below 1 K (Giauque 1964 in 1949; Tiselius 1949). It is now possible to find the absolute values of entropies for all atmospheric species: N_{2}, O_{2}, Ar, H_{2}O, CO_{2}, etc. in thermodynamic tables. All these third-law values were already available in Kelley (1932).

The use made of entropy by Rossby et al. (1937) for a moist-air isentropic analysis was based on the use of the dry-air value *θ*. This seems unrealistic.

The first application of the third law for computing the moist-air entropies in atmospheric science was probably that of Hauf and Höller (1987), who wrote that “the entropy reference value (*q*_{d}*s*_{d0} + *q*_{υ}*s*_{v0}) is not a constant,” and that “the values of the zero entropies have to be determined experimentally or by quantum statistical considerations,” namely from the third law or spectroscopic methods. They used the numerical values *s*_{v0} ≈ 10 320 J K^{−1} kg^{−1} and *s*_{d0} ≈ 6675 J K^{−1} kg^{−1} for *T*_{0} = 273.15 K and *p*_{d0} = *e*_{0} = 1000 hPa.

Another application of the third law to the atmosphere was that of Bannon (2005), who used the reference entropies given by Chase (1998) with the same value of *s*_{v0} but with a smaller value *s*_{d0} ≈ 6612 J K^{−1} kg^{−1} for dry air.

The term Λ_{r} appearing in Eq. (6) depends on the standard values at 273.15 K and 1000 hPa considered in Hauf and Höller (1987) including a correction made in M11 to account for the impact of changes in partial pressures: *e _{r}* ≈ 6.11 hPa for water vapor and

*p*

_{dr}≈ 1000 − 6.11 ≈ 994 hPa for dry air. The reference entropies are

*s*

_{υr}≈ 12 673 J K

^{−1}kg

^{−1}and

*s*

_{dr}≈ 6777 J K

^{−1}kg

^{−1}, leading to Λ

_{r}= [(

*s*

_{υ})

_{r}− (

*s*

_{d})

_{r}]/

*c*

_{pd}≈ 5.87.

The third-law reference entropies given by A5 are explicitly computed in Marquet (2015c) and Marquet and Geleyn (2015) for N_{2}, O_{2}, and H_{2}O. A larger value *s*_{d0} ≈ 6848 J K^{−1} kg^{−1} for dry air is derived by considering the change in *c*_{p}(*T*, *p*_{0}) for O_{2} at the second-order solid *α*–*β* transition occurring at 23.85 K, forming a kind of Dirac pulse with no latent heat (Fagerstroem and Hollis Hallett 1969). The corresponding term Λ_{r} ≈ 5.80 is 1.1% smaller than the one considered in M11. The reference entropies considered in Chase (1998) lead to a 2.7% larger value, Λ_{r} ≈ 6.03, where the second-order transition for O_{2} at 23.85 K is not taken into account in Joint Army–Navy–Air Force (JANAF) tables.

It is shown in Marquet (2015c, section 2.4) and Marquet and Geleyn (2015, section 5.3) that the linear combination (1 − *a*)*s*(*θ*_{e}) + *as*(*θ*_{l}) described in appendix C of Pauluis et al. (2010) can lead to the third-law value of entropy *s*(*θ*_{s/M11}) if the weighting factor is set to the value *a* = [(*s*_{d})_{r} − (*s*_{l})_{r}]*T*_{r}/*L*_{υ}(*T*_{r}) ≈ 0.356. Another value for *a* would lead to a definition of the specific moist-air entropy which would disagree with the third law. The case *a* = 1 corresponds to *S*(*θ*_{l}) if it is assumed that both (*s*_{d})_{r} and (*s*_{υ})_{r} vanish at 0°C (*T* = 272.15 K). The case *a* = 0 corresponds to *S*(*θ*_{e}) if it is assumed that both (*s*_{d})_{r} and (*s*_{l})_{r} vanish at 0°C.

The definitions of entropy of moist air selected in the IAPWS releases and guidelines (http://www.iapws.org/) are mostly based on reference-state conditions of vanishing entropies of both dry air and liquid water at the standard ocean state 273.15 K and 1013.25 hPa (IAPWS 2010; Feistel et al. 2010). However, an absolute (third law) definition for entropy of ice-lh is envisaged in Feistel and Wagner (2006). It is based on Gordon (1934) and a value of 10 320.7 J K^{−1} kg^{−1}, which is close to that taken in Hauf and Höller (1987) and M11 for water vapor at 273.15 K and 1000 hPa.

The advantage of Eq. (10) is that both *s*_{ref} and *c*_{pd} are constant, making *θ*_{s} a true equivalent of *s*. Other definitions of *θ*_{l} or *θ*_{e}, and even *θ*_{S} defined in Hauf and Höller (1987), have the disadvantage of using values of *s*_{ref} and *c*_{p} as factors of the logarithm, which both depend on *q*_{t} or *r*_{t}, thus preventing these potential temperatures from being an equivalent of the moist-air entropy if *q*_{t} is not a constant.

## REFERENCES

Bannon, P. R., 2005: Eulerian available energetics in moist atmospheres.

,*J. Atmos. Sci.***62**, 4238–4252, doi:10.1175/JAS3516.1.Barkan, D. K., 1999:

*Walther Nernst and the Transition to Modern Physical Science.*Cambridge University Press, 288 pp.Bechtold, P., E. Bazile, F. Guichard, P. Mascart, and E. Richard, 2001: A mass-flux convection scheme for regional and global models.

,*Quart. J. Roy. Meteor. Soc.***127**, 869–886, doi:10.1002/qj.49712757309.Betts, A. K., 1973: Non-precipitating cumulus convection and its parameterization.

,*Quart. J. Roy. Meteor. Soc.***99**, 178–196, doi:10.1002/qj.49709941915.Betts, A. K., 1974: Further comments on “A comparison of the equivalent potential temperature and the static energy.”

,*J. Atmos. Sci.***31**, 1713–1715, doi:10.1175/1520-0469(1974)031<1713:FCOCOT>2.0.CO;2.Betts, A. K., 1975: Parametric interpretation of trade-wind cumulus budget studies.

,*J. Atmos. Sci.***32**, 1934–1945, doi:10.1175/1520-0469(1975)032<1934:PIOTWC>2.0.CO;2.Betts, A. K., and F. J. Dugan, 1973: Empirical formula for saturation pseudoadiabats and saturation equivalent potential temperature.

,*J. Appl. Meteor.***12**, 731–732, doi:10.1175/1520-0450(1973)012<0731:EFFSPA>2.0.CO;2.Bister, M., and K. A. Emanuel, 1998: Dissipative heating and hurricane intensity.

,*Meteor. Atmos. Phys.***65**, 233–240, doi:10.1007/BF01030791.Bjerknes, V., 2009: Das problem der Wettervorhersage, betrachtet vom Standpunkte der Mechanik und der Physik (The problem of weather prediction, considered from the viewpoints of mechanics and physics).

,*Meteor. Z.***18**, 663–667, doi:10.1127/0941-2948/2009/416.Chase, J. M. W., 1998:

*JANAF Thermochemical Tables.*4th ed. American Chemical Society, 1951 pp.Coffey, P., 2006: Chemical free energies and the third law of thermodynamics.

,*Hist. Stud. Phys. Biol. Sci.***36**, 365–396, doi:10.1525/hsps.2006.36.2.365.Cuxart, J., P. Bougeault, and J.-L. Redelsperger, 2000: A turbulence scheme allowing for mesoscale and large-eddy simulations.

,*Quart. J. Roy. Meteor. Soc.***126**, 1–30, doi:10.1002/qj.49712656202.Darkow, G. L., 1968: The total energy environment of severe storms.

,*J. Appl. Meteor.***7**, 199–205, doi:10.1175/1520-0450(1968)007<0199:TTEEOS>2.0.CO;2.de Groot, S., and P. Mazur, 1984:

*Non-Equilibrium Thermodynamics.*Dover Publications, 510 pp.Emanuel, K. A., 1986: An air-sea interaction theory for tropical cyclones. Part I: Steady-state maintenance.

,*J. Atmos. Sci.***43**, 585–605, doi:10.1175/1520-0469(1986)043<0585:AASITF>2.0.CO;2.Emanuel, K. A., 1988: The maximum intensity of hurricanes.

,*J. Atmos. Sci.***45**, 1143–1155, doi:10.1175/1520-0469(1988)045<1143:TMIOH>2.0.CO;2.Emanuel, K. A., 1989: The finite-amplitude nature of tropical cyclogenesis.

,*J. Atmos. Sci.***46**, 3431–3456, doi:10.1175/1520-0469(1989)046<3431:TFANOT>2.0.CO;2.Emanuel, K. A., 1991: The theory of hurricanes.

,*Annu. Rev. Fluid Mech.***23**, 179–196, doi:10.1146/annurev.fl.23.010191.001143.Emanuel, K. A., 1994:

*Atmospheric Convection*. Oxford University Press, 580 pp.Emanuel, K. A., 1995: The behavior of a simple hurricane model using a convective scheme based on subcloud-layer entropy equilibrium.

,*J. Atmos. Sci.***52**, 3960–3968, doi:10.1175/1520-0469(1995)052<3960:TBOASH>2.0.CO;2.Emanuel, K. A., 2003: Tropical cyclones.

,*Annu. Rev. Earth Planet. Sci.***31**, 75–104, doi:10.1146/annurev.earth.31.100901.141259.Emanuel, K. A., 2004: Tropical cyclone energetics and structure.

*Atmospheric Turbulence and Mesoscale Meteorology*, E. Fedorovich, R. Rotunno, and B. Stevens, Eds., Cambridge University Press, 165–192, doi:10.1017/CBO9780511735035.010.Emanuel, K. A., and M. Bister, 1996: Moist convective velocity and buoyancy scales.

,*J. Atmos. Sci.***53**, 3276–3285, doi:10.1175/1520-0469(1996)053<3276:MCVABS>2.0.CO;2.Emanuel, K. A., and R. Rotunno, 2011: Self-stratification of tropical cyclone outflow. Part I: Implications for storm structure.

,*J. Atmos. Sci.***68**, 2236–2249, doi:10.1175/JAS-D-10-05024.1.Emanuel, K. A., M. Fantini, and A. J. Thorpe, 1987: Baroclinic instability in an environment of small stability to slantwise moist convection. Part I: Two-dimensional models.

,*J. Atmos. Sci.***44**, 1559–1573, doi:10.1175/1520-0469(1987)044<1559:BIIAEO>2.0.CO;2.Fagerstroem, C.-H., and A. C. Hollis Hallett, 1969: The specific heat of solid oxygen.

,*J. Low Temp. Phys.***1**, 3–12, doi:10.1007/BF00628329.Feistel, R., and W. Wagner, 2006: A new equation of state for H2O ice Ih.

,*J. Phys. Chem. Ref. Data***35**, 1021–1047, doi:10.1063/1.2183324.Feistel, R., D. G. Wright, H.-J. Kretzschmar, E. Hagen, S. Herrmann, and R. Span, 2010: Thermodynamic properties of sea air.

,*Ocean Sci.***6**, 91–141, doi:10.5194/os-6-91-2010.Giauque, W. F., 1964: Some consequences of low temperature research in chemical thermodynamics.

*Nobel lectures: Chemistry 1942-1962*, Elsevier, 227–250, http://www.nobelprize.org/nobel_prizes/chemistry/laureates/1949/giauque-lecture.pdf.Gibbs, J. W., 1878a: On the equilibrium of heterogeneous substances.

*Transactions of the Connecticut Academy of Arts and Sciences*, Vol. III, Connecticut Academy of Arts and Sciences, 108–248, 343–524, http://radarmet.atmos.colostate.edu/AT620/papers/Gibbs1875-1878-Equilibrium_of_Heterogeneous_Substances.pdf.Gibbs, J. W., 1878b: On the equilibrium of hererogeneous substances (continued).

,*Amer. J. Sci. Arts***16**, 441–458, http://brbl-dl.library.yale.edu/vufind/Record/3447347.Gokcen, N. A., and R. G. Reddy, 1996: The third law of thermodynamics.

*Thermodynamics*, Springer, 109–118, doi:10.1007/978-1-4899-1373-9_7.Goody, R., 2003: On the mechanical efficiency of deep, tropical convection.

,*J. Atmos. Sci.***60**, 2827–2832, doi:10.1175/1520-0469(2003)060<2827:OTMEOD>2.0.CO;2.Gordon, A. R., 1934: The calculation of thermodynamic quantities from spectroscopic data for polyatomic molecules; the free energy, entropy and heat capacity of steam.

,*J. Chem. Phys.***2**, 65–72, doi:10.1063/1.1749422.Gordon, A. R., 1935: The calculation of the free energy of polyatomic molecules from spectroscopic data. II.

,*J. Chem. Phys.***3**, 259–265, doi:10.1063/1.1749651.Gordon, A. R., and C. Barnes, 1932: The entropy of steam, and the water-gas reaction.

,*J. Phys. Chem.***36**, 1143–1151, doi:10.1021/j150334a006.Gray, S. L., 1994: Theory of mature tropical cyclones: A comparison between Kleinschmidt (1951) and Emanuel (1986). Joint Centre for Mesoscale Meteorology Tech. Rep. 40, 50 pp.

Green, J. S. A., F. H. Ludlam, and J. F. R. McIlveen, 1966: Isentropic relative-flow analysis and the parcel theory.

,*Quart. J. Roy. Meteor. Soc.***92**, 210–219, doi:10.1002/qj.49709239204.Grimus, W., 2013: 100th anniversary of the Sackur–Tetrode equation.

,*Ann. Phys.***525**, A32–A35, doi:10.1002/andp.201300720.Haber, F., 1905:

*Thermodynamik Technischer Gasreaktionen: Sieben Vorlesungen*. R. Oldenburg, 321 pp., https://archive.org/details/thermodynamikte00habegoog.Hauf, T., and H. Höller, 1987: Entropy and potential temperature.

,*J. Atmos. Sci.***44**, 2887–2901, doi:10.1175/1520-0469(1987)044<2887:EAPT>2.0.CO;2.Hawkins, H. F., and S. M. Imbembo, 1976: The structure of a small, intense hurricane—Inez 1966.

,*Mon. Wea. Rev.***104**, 418–442, doi:10.1175/1520-0493(1976)104<0418:TSOASI>2.0.CO;2.IAPWS, 2010: Guideline on an equation of state for humid air in contact with seawater and ice, consistent with the IAPWS formulation 2008 for the thermodynamic properties of seawater. International Association for the Properties of Water and Steam Publ., 21 pp., http://www.iapws.org/relguide/SeaAir.pdf.

Kelley, K. K., 1932: Contributions to the data on theoretical metallurgy: I. The entropies of inorganic substances. Bureau of Mines Bulletin 350, 63 pp., http://digicoll.manoa.hawaii.edu/techreports/PDF/USBM-350.pdf.

Kleinschmidt, E., 1951: Grundlagen einer Theorie der tropischen Zyklonen.

,*Arch. Meteor. Geophys. Bioklimatol.***4A**, 53–72, doi:10.1007/BF02246793.Klimenko, A. Y., 2012: Teaching the third law of thermodynamics.

,*Open Thermodyn. J.***6**, 1–14, doi:10.2174/1874396X01206010001.Le Chatelier, M. H., 1888:

*Recherche Expérimentales et Téoriques sur les Équilibres Chimiques*. Dunod, 229 pp., http://iris.univ-lille1.fr/handle/1908/3467.Levine, J., 1972: Comments on “A comparison of the equivalent potential temperature and the static energy.”

,*J. Atmos. Sci.***29**, 201–202, doi:10.1175/1520-0469(1972)029<0201:COCOTE>2.0.CO;2.Lewis, G. N., 1899: The development and application of a general equation for free energy and physico-chemical equilibrium.

,*Proc. Amer. Acad. Arts Sci.***35**, 3–38, doi:10.2307/25129891.Lewis, G. N., and M. Randall, 1923:

*Thermodynamics and the Free Energy of Chemical Substances*. McGraw-Hill, 653 pp.Madden, R. A., and F. E. Robitaille, 1970: A comparison of the equivalent potential temperature and the static energy.

,*J. Atmos. Sci.***27**, 327–329, doi:10.1175/1520-0469(1970)027<0327:ACOTEP>2.0.CO;2.Madden, R. A., and F. E. Robitaille, 1972: Reply.

,*J. Atmos. Sci.***29**, 202–203, doi:10.1175/1520-0469(1972)029<0202:R>2.0.CO;2.Malkus, J. S., 1958: On the structure and maintenance of the mature hurricane eye.

,*J. Meteor.***15**, 337–349, doi:10.1175/1520-0469(1958)015<0337:OTSAMO>2.0.CO;2.Marquet, P., 2011: Definition of a moist entropy potential temperature: Application to FIRE-I data flights.

,*Quart. J. Roy. Meteor. Soc.***137**, 768–791, doi:10.1002/qj.787.Marquet, P., 2015a: Definition of total energy budget equation in terms of moist-air enthalpy surface flux. Research activities in atmospheric and oceanic modelling, WCRP Rep. 12/2015, 16–17, http://www.wcrp-climate.org/WGNE/BlueBook/2015/chapters/BB_15_S4.pdf.

Marquet, P., 2015b: An improved approximation for the moist-air entropy potential temperature . Research activities in atmospheric and oceanic modelling, WCRP Rep. 12/2015, 14–15, http://www.wcrp-climate.org/WGNE/BlueBook/2015/chapters/BB_15_S4.pdf.

Marquet, P., 2015c: On the computation of moist-air specific thermal enthalpy.

,*Quart. J. Roy. Meteor. Soc.***141**, 67–84, doi:10.1002/qj.2335.Marquet, P., 2016a: Comments on “MSE minus CAPE is the true conserved variable for an adiabatically lifted parcel.”

,*J. Atmos. Sci.***73**, 2565–2575, doi:10.1175/JAS-D-15-0299.1.Marquet, P., 2016b: Etude de l’énergétique de l’air humide et des paramétrisations de l’atmosphère: Propriété de l’exergie, de l’enthalpie utilisable, de l’entropie et de l’enthalpie. Institut National Polytechnique de Toulouse, 307 pp., https://tel.archives-ouvertes.fr/tel-01504276/file/Marquet_HDR_2017_short.pdf.

Marquet, P., 2016c: The mixed-phase version of moist-entropy. Research activities in atmospheric and oceanic modelling, WCRP Rep. 15/2016, 7–8, http://www.wcrp-climate.org/WGNE/BlueBook/2016/documents/Sections/BB_16_S4.pdf.

Marquet, P., and J.-F. Geleyn, 2015: Formulations of moist thermodynamics for atmospheric modelling.

*Current Issues and New Theories*, R. S. Plant and J.-I. Yano, Eds., Vol. II,*Parameterization of Atmospheric Convection*, World Scientific, 221–274, doi:10.1142/9781783266913_0026.Masanes, L., and J. Oppenheim, 2017: A general derivation and quantification of the third law of thermodynamics.

,*Nat. Commun.***8**, 14538, doi:10.1038/ncomms14538.Miller, B. I., 1958: On the maximum intensity of hurricane.

,*J. Meteor.***15**, 184–195, doi:10.1175/1520-0469(1958)015<0184:OTMIOH>2.0.CO;2.Montroty, R., F. Rabier, S. Westrelin, G. Faure, and N. Viltard, 2008: Impact of wind bogus and cloud- and rain-affected SSM/I data on tropical cyclone analyses and forecasts.

,*Quart. J. Roy. Meteor. Soc.***134**, 1673–1699, doi:10.1002/qj.308.Mrowiec, A. A., O. M. Pauluis, and F. Zhang, 2016: Isentropic analysis of a simulated hurricane.

,*J. Atmos. Sci.***73**, 1857–1870, doi:10.1175/JAS-D-15-0063.1.Nagle, J. F., 1966: Lattice statistics of hydrogen bonded crystals. I. The residual entropy of ice.

,*J. Math. Phys.***7**, 1484–1491, doi:10.1063/1.1705058.Nernst, W., 1907: Integration of the equation of the reaction isochore, preliminary discussion of the undetermined constant of integration, and of the relation between the total and the free energies at very low temperatures.

*Experimental and Theoretical Applications of Thermodynamics to Chemistry*, W. Nernst, Ed., Lectures at the Yale University, Charles Scribner’s Sons, 39–52, https://archive.org/details/experimentaland04nerngoog.Nernst, W., 1912: Thermodynamik und spezifische Wärme.

,*Sitzungsberichte der Königlich Preussischen Akademie der Wissenschaften zu Berlin***6**, 134–140.Ooyama, K. V., 1990: A thermodynamic foundation for modeling the moist atmosphere.

,*J. Atmos. Sci.***47**, 2580–2593, doi:10.1175/1520-0469(1990)047<2580:ATFFMT>2.0.CO;2.Palmén, E., and H. Riehl, 1957: Budget of angular momentum and energy in tropical cyclones.

,*J. Meteor.***14**, 150–159, doi:10.1175/1520-0469(1957)014<0150:BOAMAE>2.0.CO;2.Pauling, L., 1935: The structure and entropy of ice and of other crystals with some randomness of atomic arrangement.

,*J. Amer. Chem. Soc.***57**, 2680–2684, doi:10.1021/ja01315a102.Pauluis, O., 2011: Water vapor and mechanical work: A comparison of Carnot and steam cycles.

,*J. Atmos. Sci.***68**, 91–102, doi:10.1175/2010JAS3530.1.Pauluis, O., and I. M. Held, 2002a: Entropy budget of an atmosphere in radiative–convective equilibrium. Part I: Maximum work and frictional dissipation.

,*J. Atmos. Sci.***59**, 125–139, doi:10.1175/1520-0469(2002)059<0125:EBOAAI>2.0.CO;2.Pauluis, O., and I. M. Held, 2002b: Entropy budget of an atmosphere in radiative–convective equilibrium. Part II: Latent heat transport and moist processes.

,*J. Atmos. Sci.***59**, 140–149, doi:10.1175/1520-0469(2002)059<0140:EBOAAI>2.0.CO;2.Pauluis, O., and A. A. Mrowiec, 2013: Isentropic analysis of convective motions.

,*J. Atmos. Sci.***70**, 3673–3688, doi:10.1175/JAS-D-12-0205.1.Pauluis, O., V. Balaji, and I. M. Held, 2000: Frictional dissipation in a precipitating atmosphere.

,*J. Atmos. Sci.***57**, 989–994, doi:10.1175/1520-0469(2000)057<0989:FDIAPA>2.0.CO;2.Pauluis, O., A. Czaja, and R. Korty, 2008: The global atmospheric circulation on moist isentropes.

,*Science***321**, 1075–1078, doi:10.1126/science.1159649.Pauluis, O., A. Czaja, and R. Korty, 2010: The global atmospheric circulation in moist isentropic coordinates.

,*J. Climate***23**, 3077–3093, doi:10.1175/2009JCLI2789.1.Peixoto, J. P., A. H. Oort, M. De Almeida, and A. Tomé, 1991: Entropy budget of the atmosphere.

,*J. Geophys. Res.***96**, 10 981–10 988, doi:10.1029/91JD00721.Planck, M., 1917:

*Treatise on Thermodynamics*(in German). Dover, 297 pp.Raymond, D. J., 2013: Sources and sinks of entropy in the atmosphere.

,*J. Adv. Model. Earth Syst.***5**, 755–763, doi:10.1002/jame.20050.Rennó, N. O., 2001: Comments on “Frictional dissipation in a precipitating atmosphere.”

,*J. Atmos. Sci.***58**, 1173–1177, doi:10.1175/1520-0469(2001)058<1173:COFDIA>2.0.CO;2.Rennó, N. O., and A. P. Ingersoll, 1996: Natural convection as a heat engine: A theory for CAPE.

,*J. Atmos. Sci.***53**, 572–585, doi:10.1175/1520-0469(1996)053<0572:NCAAHE>2.0.CO;2.Richards, T. W., 1902: The significance of changing atomic volume. III. The relation of changing heat capacity to change of free energy, heat of reaction, change of volume, and chemical affinity.

,*Proc. Amer. Acad. Arts Sci.***38**, 293–317, doi:10.2307/20021755.Richardson, L. F., 1922:

*Weather Prediction by Numerical Process*. Cambridge University Press, 236 pp.Riehl, H., 1950: A model of hurricane formation.

,*J. Appl. Phys.***21**, 917–925, doi:10.1063/1.1699784.Romps, D. M., 2008: The dry-entropy budget of a moist atmosphere.

,*J. Atmos. Sci.***65**, 3779–3799, doi:10.1175/2008JAS2679.1.Romps, D. M., 2010: A direct measure of entrainment.

,*J. Atmos. Sci.***67**, 1908–1927, doi:10.1175/2010JAS3371.1.Romps, D. M., and Z. Kuang, 2010: Do undiluted convective plumes exist in the upper tropical troposphere?

,*J. Atmos. Sci.***67**, 468–484, doi:10.1175/2009JAS3184.1.Rossby, C.-G., D. P. Keily, J. W. W. Osmin, and J. Namias, 1937: Isentropic analysis.

,*Bull. Amer. Meteor. Soc.***18**, 201–209.Rotunno, R., and K. A. Emanuel, 1987: An air–sea interaction theory for tropical cyclones. Part II: Evolutionary study using a nonhydrostatic axisymmetric numerical model.

,*J. Atmos. Sci.***44**, 542–561, doi:10.1175/1520-0469(1987)044<0542:AAITFT>2.0.CO;2.Schrödinger, E., 1989:

*Statistical Thermodynamics*. Dover, 95 pp.Smith, R. K., M. T. Montgomery, and S. Vogl, 2008: A critique of Emanuel’s hurricane model and potential intensity theory.

,*Quart. J. Roy. Meteor. Soc.***134**, 551–561, doi:10.1002/qj.241.Tiselius, A., 1949: The Nobel Prize in Chemistry 1949: Award ceremony speech. Noble Media AB, http://www.nobelprize.org/nobel_prizes/chemistry/laureates/1949/press.html.

van’t Hoff, J. H., 1904: Einfluss der Änderung der spezifischen Wärme auf die Umwandlungsarbeit (The influence of the change in specific heat on the work of conversion). Festschrift Ludwig Boltzmann gewidmet zum sechzigsten

*Geburtstage 20 Februar 1904*(*Commemorative of the Sixtieth Birthday of Ludwig Boltzmann on 20 February 1904*), J. A. Barth, 233–241, https://archive.org/details/festschriftludw00meyegoog.Wang, J. J., and Coauthors, 2015: A long-term, high-quality, high-vertical-resolution GPS dropsonde dataset for hurricane and other studies.

,*Bull. Amer. Meteor. Soc.***96**, 961–973, doi:10.1175/BAMS-D-13-00203.1.Wilks, J., 1961:

*The Third Law of Thermodynamics*. Oxford University Press, 142 pp.