The moist-air entropy can be used to analyze and better understand the general circulation of the atmosphere or convective motions. Isentropic analyses are commonly based on studies of different equivalent potential temperatures, all of which are assumed to fully represent the entropy of moist air. It is, however, possible to rely either on statistical physics or the third law of thermodynamics when defining and computing the absolute entropy of moist air and to study the corresponding third-law potential temperature, which is different from the previous ones. The third law assumes that the entropy for the most stable crystalline state of all substances is zero when approaching absolute zero temperature.
This paper shows that the way all these moist-air potential temperatures are defined has a large impact on (i) the plotting of the isentropes for a simulation of Hurricane Dumile, (ii) the changes in moist-air entropy computed for a steam cycle defined for this hurricane, (iii) the analyses of isentropic streamfunctions computed for this hurricane, and (iv) the computations of the heat input, the work function, and the efficiency defined for this steam cycle.
The moist-air entropy is a state function and the isentropic analyses must be completely determined by the local moist-air conditions. The large differences observed between the different formulations of moist-air entropy are interpreted as proof that the isentropic analyses of moist-air atmospheric motions must be achieved by using the third-law potential temperature defined from general thermodynamics.
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 ≈ S0m + Cpm ln(θe) or s* ≈ S0m + Cpm 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 S0m and Cpm, 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 CpmT + Lυqυ + gz or CpmT + Lυ qsw + gz, where Cpm may depend on water vapor or condensed water contents, Lυ is the latent heat of vaporization, and qυ or qsw 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) = sref + cpd ln(θs), where both sref and cpd 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
where , T is the temperature, p0 is the standard pressure, pd is the dry-air pressure, and Rd and Rυ are the gas constant of water vapor and dry air, respectively. The specific heat = cpd + rtcl depends on the values for dry air (cpd) and liquid water (cl). The mixing ratios rυ and rt represent water vapor and total water, respectively. The quantity Lυ is the latent heat of vaporization. The relative humidity with respect to liquid water Hl = e/esw is the ratio of the water vapor pressure (e) to the saturated value (esw).
where p is the total pressure and rsw(T, p) is the saturation mixing ratio at temperature T and pressure p.
where the dry-air potential temperature is θ = T(p0/p)κ with κ = Rd/cpd ≈ 0.286.
The equivalent potential temperature θe on page 1860 of MPZ can be written as
The difference between θe/MPZ and θe/E94 is that pd is replaced by p in the Exner function. The differences between θe/MPZ and θe/E86 are that rsw is replaced by rυ and the additional term depending on Hl is included. The differences between θe/MPZ and θe/B73 are that cpd is replaced by , qυ is replaced by rυ, and the additional term depending on Hl 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 pd into p may have a large impact on the definition of isentropic processes.
where λ = cpυ/cpd − 1 ≈ 0.837, η = Rυ/Rd ≈ 1.608, ε = Rd/Rυ ≈ 0.622, δ = η − 1 ≈ 0.608, and γ = Rυ/cpd ≈ 0.46. The term Λr = [(sυ)r − (sd)r]/cpd ≈ 5.87 depends on the third-law values of the reference entropies for water vapor and dry air, and Ls is the latent heat of sublimation. The water vapor, liquid, ice, and total specific contents qυ, ql, qi, and qt = qυ + ql + qi replace the mixing ratios involved in most of the previous formulations.
The water vapor and dry-air reference entropies (sυ)r = sυ(Tr, er) and (sd)r = sd (Tr, pr − er) depend on the reference temperature Tr = 273.15 K, the total reference pressure pr = 1000 hPa, and the water vapor partial pressure er. However, it is shown in M11 that θs/M11 defined in Eq. (6) is independent of Tr and pr if the reference mixing ratio rr = εer/(pr − er) ≈ 3.82 g kg−1 corresponds to the saturating pressure value at Tr: er = esw(Tr) = 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 (Hl ≠ 1) or ice (Hi ≠ 1) and possible temperatures of rain Tl, or snow, Ti, which may differ from the T for dry air and water vapor.
The advantage of the term in Eq. (6) compared with in Eq. (1) or (5) is that ql replaces rυ in the exponent, making no impact in clear air, or under- or supersaturated moist regions (where Hl ≠ 1 and rυ may be large, but where ql = 0) and having less impact in cloud in under or supersaturated regions (where Hl ≠ 1 but where, typically, ql < rυ).
where 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
The moist-air entropy is computed in M11 from the third law of thermodynamics. It can be written as
where both sref ≈ 1139 J K−1 kg−1 and cpd ≈ 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.
Other definitions of “moist-air entropy” are derived with either sref or cpd (often both of them) replaced by quantities that depend on the total-water mixing ratio rt. This is true in Eq. (4.5.10) in E94, which can be written as
where p0 ≈ 1000 hPa is a constant standard value. The division of s by qd 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 qt varies in space or time, since these missing terms must be multiplied by qd = 1 − qt to compute s from s/qd given by Eq. (12). Moreover, since = cpd + rtcl depends on rt, changes in s and s/qd cannot be represented by θe/E94 in Eq. (12) for varying values of rt. 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 rt) are to be compared with nonsaturated ones (small values of rt).
Similarly, the moist-air entropy defined in section 3 in MPZ can be written as
where T0 = 273.15 K is a constant standard value. Again, it is an entropy expressed per unit mass of dry air and the specific heat depends on varying values of rt, twice preventing θe/MPZ from being a true equivalent of the specific moist entropy.
To better analyze the impact of rt on the term , and thus on the definition of the moist-air entropy, two kinds of “saturated equivalent entropy” are defined. They are based on the definition of θ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 and the relative humidity Hl. The use of gives a clear, unambiguous definition of thermal properties, differing from the uncertain and multiple definitions 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υ, ql = qi = 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).
The impacts of condensed water and of nonequilibrium terms will not be tested with these 15 points, chosen such that Tl = Ti = T everywhere and Hl = Hi = 1 and ql = qi = 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 qt = 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 . This is a kind of validation of the modeling of Hurricane Dumile by ALADIN. The same high values of observed in the eye of Inez are simulated in Fig. 1 for Dumile at the lower and upper levels. The same area of minimum and lower relative humidity is observed in the core region in the 500–700-hPa layer. Figure 3 shows that the general patterns for θe/B73 are similar to those for shown in Fig. 1, but with change of labels for the potential temperature units. This confirms that the isolines defined by constant values of either or θe/B73 almost overlap.
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 and θ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 qt on the moist-air entropy variable θs/M11 ≈ θl exp(6qt), where the factor 6 is about two-thirds of that included in θe/E94 ≈ θl exp(9qt). 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 and θ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 and θe/B73.
The isentropes plotted with θs in Fig. 5 almost coincide with the isolines of θe/B73 in the dry regions where qt is small, namely in the high troposphere above the 200-hPa level and west of the core region (longitudes < 50°). The specific humidity qt 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 qt are simulated).
The isentropes computed with θs are therefore not compatible with the isolines of or θe/B73 and the larger the specific content qt and the relative humidity Hl 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
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 rsw) 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 with Hl ≈ 60%, Rυ /≈ Rυ /cpd ≈ 0.46 and rsw ≈ 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 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
An isentropic streamfunction Ψ(θ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 as
where θe/s stands either for θe/E94 or θs/M11, respectively. The isentropic mass flux 〈ρw〉(r, θ, z) is computed according to
The 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.
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(9qt) overestimates the impact of humidity in comparison with the third-law value θs/M11 ≈ θl exp(6qt). The ratio 2/3 for the coefficients in factor of qt 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.
According to de Groot and Mazur (1984) and for a state of local equilibrium, the Gibbs differential equation for open systems can be written as
where ei is the specific internal energy, h the specific enthalpy, ρ the density, and μk = hk −Tsk the specific Gibbs’ function for k = 0, 1 2, 3 denoting dry-air, water vapor, liquid water, and ice, respectively.
The work produced per unit mass of moist air W = −ρ−1dp can then be computed by integrating Eq. (18) along a closed circuit. Since the integral of the total differential of enthalpy vanishes along such a loop, W can be written as
The last integral is derived from the last integral in Eq. (18) for the steam cycle considered in Figs. 1 to 4, where ql = qi = 0, qυ = qt, and qd + qt = 1, leading to dql = dqi = 0, dqd = −dqυ = −dqt, and thus to μddqd + μυdqυ = (μυ − μd)dqt. The same equation for W holds true if liquid or solid water exist, provided that changes of phases are reversible, which implies (μυ − μl)dql = (μυ − μi)dqi = 0.
With the help of , , and μk = hk − Tsk, the Gibbs Eq. (18) can be rewritten as
This relation clearly shows that, since dhk = cpkdT (for all species), dsk = cpkdT/T − Rkdpk/pk (for gases) and dsk = cpkdT/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 hv0 and hd0 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υ − hd) − T(sυ − sd) . This property holds true because (hv0 − hd0) dqt = (hv0 − hd0) dqt = 0 since dqt ≡ 0 for a closed loop.
In contrast, the reference entropies must be taken into account in the computations of T(sd0 − sv0)dqt = (sd0 − sv0)Tdqt, because T(qt)dqt ≠ 0 if the area of the loop is different from zero in a T–qt diagram. However, the impact of the reference entropies sd0 and sv0 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 sd0 or sv0.
The integrals in Eq. (19) are computed in E94, Goody (2003), and P11 with all extensive quantities expressed per unit of dry air, and thus with the specific values of h and s replaced by h/qd and s/qd. Furthermore, the reference entropies for dry air and liquid water are assumed to vanish at T0 = 0°C in P11, leading to
The 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)]—yielding
In Eqs. (23) and (24) WH and WH(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 WH and WH(d) listed in Table 4 are calculated using the simple but robust first-order Riemann method.
Large factors, of 0.5 and 4, are observed between WH ≈ 870 J kg−1, computed with the third-law value s(θs), and WH(d) ≈ 442 J kg−1 or WH(d) ≈ 3405 J kg−1, computed with s(θes/E86)/qd or s(θe/E94)/qd, respectively. Furthermore, the impact of is very large compared to that of cpd, leading to a factor of more than 6 for WH(d) computed with the same potential temperature θes/E86 (442 vs 2758 J kg−1).
The impact on WH 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/qd 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 WH = 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 WH and WH(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 WH and WH(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%.
The efficiency of a hurricane considered as a steam cycle is defined in P11 by taking account of the last integrals in Eqs. (19) and (21), which are intended to model the effects due to the open cycle in which mass is added and removed. These Eqs. (19) and (21) can be rewritten as W = Q1 = Q2 and W(d) = Q1(d) = Q2(d), where
The “energy input” is defined in P11 as Qin = δ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 ratio
where all quantities are expressed per unit mass of dry air. From Eq. (22) in P11, this ratio can be approximated by χS ≈ (Tin − Tout)/1540 K, where Tin is the highest and Tout the lowest temperature of the cycle and where (RdLυ)/(Rυcpd) ≈ 1540 K. The steam cycle considered for Hurricane Dumile corresponds to Tin ≈ 296 K and Tout ≈ 265 K, leading to χS ≈ 0.020.
Table 6 shows that W = Q1 ≈ Q2 and W(d) = Q1(d) ≈ Q2(d) with the same accuracy of about 40 J kg−1 previously noted for the Riemann method. The differences Q1 ≠ Q1(d) and Q2 ≠ Q2(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).
Similarly, the way δQ+ and Qin are computed may impact the denominator of Eq. (29). The computations of δQ+ for Q1, Q2, Q1(d), or Q2(d) are especially uneasy because they rely on accurate definitions of all the enthalpies (h of h/qd), the entropies (s or s/qd), 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 δQ1 > 0.
The enthalpy (h/qd), entropy (s/qd), 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 Lh(Tr) ≈ 2603 kJ kg−1 and href ≈ 256 kJ kg−1. The Gibbs functions μυ and μd are computed by hυ −Tsυ and hd −Tsd, where hd(Tr) ≈ 530 kJ kg−1 and hυ(Tr) ≈ 3133 kJ kg−1.
The latent heat Lh(T) defined in Eq. (31) depends on the values of Lh(Tr), in the same way as Lυ depends on Lυ(Tr) ≈ 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 (sd)r.
However, Lh(T) is independent of the reference temperature Tr because Lh(Tr) + (cpυ − cpd)Tr is a constant, in the same way as Lυ and Ls are independent of the value Tr used to compute them, and in the same way as θs/M11 does not depend on Tr, pr, and rr = rsw(Tr, pr) = εer/(pr − er) provided that er = esw(Tr).
The efficiency factors shown in the last column of Table 6 are almost the same for Q1 and Q2 on the one hand and for Q1(d) and Q2(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 Qin/1 and Qin/2, which are only 1% smaller than Qin/1(d) and Qin/2(d), whereas the full work functions Q1(d) and Q2(d) are 50% larger than Q1 and Q2.
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/qd, h or h/qd; modified reference values for entropies and enthalpies; the choice of θs, θe, or θes, the choice of cpd or as a factor of the logarithm; etc.
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/qd, the definitions of the geopotential, the wind components, or the kinetic energy one should modify accordingly by plotting, for instance, gz/qd, u/qd, υ/qd, or (u2 + υ2)/(2qd). These definitions are unusual in the thermodynamics of open systems. Moreover, if s/qd could be defined within a global constant C, the specific value s would depend on qdC, which varies with qd and renders the integral indeterminate, because it would depend on , where 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 qt 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/).
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.
In the mid- and late nineteenth century, an important problem in chemistry was how to improve the understanding of chemical reactions and to predict their spontaneity. The answer given by Gibbs (1878a,b) was to compute the sign of the Gibbs function (or free enthalpy) of reaction
where Δ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.
The problem of finding the numerical values of Δ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 equations
where 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).
The advantage of Eq. (A2) is that it provides the actual possibility of measuring the heat 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 T0 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(T0)/T0 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 H2 + O2 = H2O in the atmosphere by choosing any arbitrary value for the constant of integration A(T0)/T0 for the temperature T0.
This problem was solved by Nernst (1907) in 1906 by assuming (page 49) a “new hypothesis,” which is nowadays called the “heat theorem”:
The heat theorem was not expressed in terms of the entropy but, instead, was based on the Helmholtz function A, leading to A(T0)/T0 → 0 as T0 → 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).
The “finite density” hypothesis clearly excludes applications to perfect gases at 0 K and the third law must only be applied to the solid state.
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.”
Accordingly, the standard third-law values of entropies at temperature T0 and pressure p0 can be written as
The integral of the piecewise function cp(T, p0) is computed from 0 K to T0 for all solid(s), liquid, and gaseous states. The sum over k is extended to all transitions of phases occurring at Tk with the latent heats L(Tk).
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 are provided.
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 O2, N2, H2O, and CO2), 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 = S0 + k ln(W), namely that S0 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 S0 = 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 S0 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 S0 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 H2O and GR96 (45.132 vs 44.31) is not observed for MG15 (45.2). The explanation for the larger value observed for O2 and MG15 (49.7 vs 49.031) is given in 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.
The close link between the two subjects can be understood by analyzing the simple case of mixing dry air and water vapor, both of them considered as perfect gases. If the gas constants Rd and Rυ and specific heats cpd and cpυ 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 as
where p = pd + e, qd = 1 − qυ, and cp = qdcpd + qυcpυ and with the use of the partial pressures p, e, pd0, and e0, 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(T0)/T0, sd0, and sv0 by the variable terms T(x, y, z), qd(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, qd(x, y, z)sd0 and qυ(x, y, z)sv0. Consequently, the two reference values sd0 and sv0 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 sd0 and sv0.
However, Richardson was not able to continue accurate computations of the moist-air entropy in 1922 because values of cp(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: N2, O2, Ar, H2O, CO2, 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 (qdsd0 + qυsv0) 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 sv0 ≈ 10 320 J K−1 kg−1 and sd0 ≈ 6675 J K−1 kg−1 for T0 = 273.15 K and pd0 = e0 = 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 sv0 but with a smaller value sd0 ≈ 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: er ≈ 6.11 hPa for water vapor and pdr ≈ 1000 − 6.11 ≈ 994 hPa for dry air. The reference entropies are sυr ≈ 12 673 J K−1 kg−1 and sdr ≈ 6777 J K−1 kg−1, leading to Λr = [(sυ)r − (sd)r]/cpd ≈ 5.87.
The third-law reference entropies given by A5 are explicitly computed in Marquet (2015c) and Marquet and Geleyn (2015) for N2, O2, and H2O. A larger value sd0 ≈ 6848 J K−1 kg−1 for dry air is derived by considering the change in cp(T, p0) for O2 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 O2 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 = [(sd)r − (sl)r]Tr/Lυ(Tr) ≈ 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 (sd)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 (sd)r and (sl)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 sref and cpd 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 sref and cp as factors of the logarithm, which both depend on qt or rt, thus preventing these potential temperatures from being an equivalent of the moist-air entropy if qt is not a constant.
A comment/reply has been published regarding this article and can be found at http://journals.ametsoc.org/doi/abs/10.1175/JAS-D-18-0059.1 and http://journals.ametsoc.org/doi/abs/10.1175/JAS-D-18-0126.1