## 1. Introduction

Within the earth’s atmosphere, water can be found in gaseous, liquid, and solid states. Water molecules undergo multiple phase transitions as they are transported by atmospheric flows. These phase transitions are associated with significant release or absorption of latent heat, and play an important role in the atmospheric circulation. For example, evaporation of water accounts for approximately 80% of the energy exchange between the atmosphere and oceans (Peixoto and Oort 1992). Phase transitions are at the core of the complex interplay between thermodynamics and fluid dynamics that is an intrinsic aspect of many weather phenomena ranging from convection to hurricanes and midlatitude storms.

Moist air can be treated as an ideal mixture of dry air, water vapor, and condensed water. The state of moist air can be uniquely described by the combination of entropy *S*, total pressure *p*, specific humidity *q _{υ}*, and liquid water content

*q*. However, not all possible combinations of these four state variables can be realized within the atmosphere. Indeed, air parcels within the lower 100 km of the atmosphere are considered to be in local thermodynamic equilibrium. This assumption implies that, when water vapor and condensed water are both present, the partial pressure of the water vapor in the parcel is equal to the saturation vapor pressure. Consistently with Gibbs’s phase rule, the requirement that the different water phases be in thermodynamic equilibrium reduces by one the number of variables necessary to describe the state of moist air. Here, we introduce the total water content

_{l}*q*=

_{T}*q*+

_{υ}*q*as the mass of water in all phases per unit mass of moist air, so that the state of a moist air parcel is uniquely determined by its entropy

_{l}*S*, pressure

*p*, and total water content

*q*. In particular, the partitioning of the total water content

_{T}*q*between the water vapor

_{T}*q*and liquid water

_{υ}*q*can be obtained once the three state variables,

_{l}*S*,

*p*, and

*q*, are known. For example, at constant pressure and total water content, higher entropy corresponds to higher temperature and water vapor content but lower condensed water content. At sufficiently large entropy, low pressure, or low total water content, all of the water will be in the vapor phase, with no condensed water present. Air parcels with only water vapor and no condensed phase are referred to as unsaturated parcels, while parcels with both phases present are referred to as saturated parcels.

_{T}*T*/∂

*z*). For a typical tropical sounding, it is roughly 0.01 K m

^{−1}near the surface but decreases abruptly to a value on the order of 0.004 K m

^{−1}above the cloud base. This abrupt change in the lapse rate is a reflection of the difference in behavior between unsaturated and saturated air parcels. The temperature profile in regions of active convection in the tropics follows closely that of a parcel raised adiabatically from the surface. For such an adiabatic ascent, the temperature

*T*is a function of the entropy

*S*, total water content

*q*, and pressure

_{T}*p:*Below the cloud base, the lapse rate corresponds to the adiabatic expansion of an unsaturated parcel withwith

*ρ*the mass of air per unit volume,

*g*the gravitational acceleration, and

*C*the heat capacity at constant pressure. The second equality arises from the hydrostatic balance: (∂

_{p}*p*/∂

*z*) = −

*ρg*. Above the cloud base, the lapse rate corresponds to the adiabatic expansion of the same parcel that is now saturated. Emanuel (1994) and others have computed the lapse rate for a saturated ascent:where

*L*is the latent heat of vaporization;

_{υ}*R*and

_{d}*R*are the ideal gas constants of dry air and water vapor, respectively; and

_{υ}*ϵ*= (

*R*

_{d}/

*R*

_{υ}) is their ratio. The approximation lies in neglecting the contribution of the condensed water and replacing the mixing ratio by the specific humidity. The saturated lapse rate is always lower than the unsaturated lapse rate because the cooling due to the parcel’s expansion is partially compensated for by the latent heat release due to condensation. Observations in the tropics reflect the fact that, for a parcel raised adiabatically, the lapse rate is given by the unsaturated value Γ

*below the cloud base, and abruptly drops to its saturated value Γ*

_{u}*above it.*

_{m}*F*(

*S*,

*p*,

*q*,

_{υ}*q*) with

_{l}*q*the specific humidity and

_{υ}*q*the liquid water content. The state property

_{l}*F*is a smooth function and can be determined for any combination of

*S*,

*p*,

*q*, and

_{υ}*q*. However, as discussed above, not all of these combinations can be realized due to the requirement that moist air be in local thermodynamic equilibrium. For an unsaturated air parcel, only water vapor is present,

_{l}*q*=

_{υ}*q*and

_{T}*q*= 0. Partial derivatives must be taken assuming that

_{l}*dq*= 0 and

_{l}*dq*=

_{υ}*dq*; for example,In contrast, for a saturated parcel with both water vapor and liquid water present, the thermodynamic equilibrium between the two phases is expressed by requiring the chemical potential of both water vapor

_{T}*μ*and liquid water

_{υ}*μ*to be equal:

_{l}*μ*=

_{υ}*μ*. In this case, partial derivatives are taken assuming both

_{l}*dq*+

_{υ}*dq*=

_{l}*dq*and

_{T}*dμ*=

_{υ}*dμ*; for example,Local thermodynamic equilibrium leads to different constraints on saturated and unsaturated parcels. Partial derivatives are taken along different directions in four-dimensional space (

_{l}*S*,

*q*,

_{υ}*q*,

_{l}*p*) depending on whether the air is saturated or unsaturated. As a result, the partial derivatives of the state equation are discontinuous at the saturation point.

The Boussinesq approximation (Oberbeck 1879; Boussinesq 1903; Spiegel and Veronis 1960; Mihaljan 1962) offers a practical simplification of the Navier–Stokes equations for a compressible fluid by replacing the continuity equation with a nondivergence constraint on the mass flux. This not only simplifies the analytic treatment of the equations, but also filters out fast-propagating sound waves. The use of the Boussinesq approximation to study atmospheric motion is limited by the fact that it requires the vertical extent of the domain to be small in comparison to the density-scale height. This problem is addressed in the anelastic approximation by allowing for a height-dependent reference density profile (Batchelor 1953; Ogura and Phillips 1962; Dutton and Fitchl 1969; Lipps and Hemler 1982; Durran 1989; Vallis 2006). The original derivations of the Boussinesq and anelastic approximations require the atmosphere to be close to a dry isentropic reference state in hydrostatic balance, and an equation of state based on the ideal gas law, with density proportional to the potential temperature.

The primary motivation of this work is to ensure that the use of the anelastic approximation is fully consistent even for a fluid exhibiting a highly nonlinear thermodynamic state law such as moist air. Section 2 discusses an implementation anelastic approximation for a moist atmosphere. It assumes that the atmospheric density distribution is a small perturbation around a hydrostatic reference profile with constant entropy and total water content. The impacts of density and pressure on the momentum equation are linearized, so that density fluctuations enter the momentum equations only as a vertical buoyancy acceleration. A central feature of this version of the anelastic approximation is that it keeps the full nonlinear dependency of the buoyancy on the entropy and total water content.

The central result presented in section 3 can be seen as an extension of Ingersoll’s (2005) work on nonlinear equations of state. It is shown that the anelastic approximation conserves the sum of the kinetic and thermodynamic energy. The thermodynamic energy here is equal to the difference between the enthalpy of the parcel estimated at the reference pressure and the enthalpy of the reference state at the same level. The temperature and chemical potential are given by the temperature and chemical potential of a moist air parcel with the same entropy and total water content but they are estimated at the reference pressure. Furthermore, the thermodynamic relations between these various quantities under the anelastic approximation are the same as these for moist air except that the pressure has been replaced by the reference state pressure. As the anelastic approximation discussed here uses entropy as a prognostic variable, this results implies that this anelastic model simultaneously obey the first and second laws of thermodynamics.

Sections 4 and 5 address some implications of this result. For example, analysis of the entropy and energy budgets such as that in Pauluis and Held (2002) can provide an estimate of the work produced by the atmosphere. In section 4, it is shown that, under the anelastic approximation, the same considerations yield a constraint on the vertically integrated buoyancy flux. This is illustrated here with a simple example that compares a direct calculation of the buoyancy flux with the predictions from the thermodynamic constraints. Section 5 discusses the possibility of using enthalpy or temperature as a prognostic variable instead of entropy, with the evolution equation obtained from the thermodynamic energy equation. The entropy in this model would be given by the entropy of the parcel evaluated at the reference state pressure. Such a model, while explicitly obeying the first law of thermodynamics, would also implicitly obey the second law.

## 2. The anelastic approximation

We are interested here in deriving the anelastic approximation for a fluid in which the specific volume *α*(*S*, *q _{T}*,

*p*) is an arbitrary function of three state variables: here, the entropy

*S*, total water content

*q*and pressure

_{T}*p*. As discussed in the introduction, moist air can be treated as an ideal mixture of dry air, water vapor, and condensed water, with the water phases in thermodynamic equilibrium.

^{1}This later requirement implies that the partitioning of the total water between its gas and condensed phases can be uniquely determined from the three state variables

*S*,

*p*, and

*q*, consistently with Gibb’s phase rule. The state law for moist air exhibits a strong nonlinear behavior in that its partial derivatives are discontinuous at the saturation point. In this section, we derive a version of the anelastic approximation that maintains the nonlinear dependency of the specific volume on both the entropy and total water content.

_{T}*S*

_{0}and total water content

*q*

_{T}_{0}. The entropy here includes contributions from the water vapor and condensed water. The reference state pressure

*p*

_{0}(

*z*) is in hydrostatic balance:with

*α*

_{0}(

*z*) =

*α*[

*S*

_{0},

*q*

_{T}_{0},

*p*

_{0}(

*z*)] the specific volume in the reference state. The total pressure is then given by

*p*(

*x*,

*y*,

*z*,

*t*) =

*p*

_{0}(

*z*) +

*δp*(

*x*,

*y*,

*z*,

*t*), and the specific volume is expanded asThe contributions of the pressure and specific volume terms in the momentum equations are linearized:Note that in (5), the dependency of the specific volume on pressure has been linearized but the full nonlinear dependency of the specific volume on the two other states variables,

*S*and

*q*, has been kept. This has some specific advantages over versions based on a fully linearized equation of state. First, while the pressure and specific volume perturbations must be small for the expansion of the momentum equation to remain accurate, a nonlinear equation of state allows for large fluctuations of the entropy or the total water content as long as their combined impacts on the specific volume is small. Second, this makes it possible to handle discontinuous partial derivatives and other complex aspects of the equation of state. Finally, we will establish in section 3 that preserving the nonlinearity dependency of the specific volume on entropy and the total water content also preserves some key thermodynamics relationships.

_{T}*B*is given byThe linearization of the momentum equation is based on the Taylor expansion of the specific volume

*α*and pressure terms

*p*, and is only accurate for small fluctuations in these quantities.

*U*/

*c*

_{s}), with

*U*a velocity scale for the flow and

*c*the speed of the sound. Hence, only flows with low Mach numbers (Mach ≪ 1) can be accurately represented under the anelastic approximation. This requirement is separate from the previous requirement that the density fluctuations be small.

_{s}*ρ*

_{0}

**u**:The equation differs from the kinetic energy for a compressible fluid in that it uses the density of the reference state rather than the actual density of the fluid, and in that the net production of kinetic energy is given by the buoyancy flux

*ρ*

_{0}

*Bw*.

## 3. Thermodynamic energy in the anelastic approximation

**∇**· (

**u**

*δp*) is the divergence of the flux, it integrates to zero over a closed domain and does not affect the energy conservation. The objective here is to find a function of entropy and total water content (and possibly height and pressure) whose tendency cancels out the buoyancy term when added to the kinetic energy equation. Following Ingersoll (2005) and Vallis (2006, chapter 2.5), we define the thermodynamic energy asHere, the function Φ

_{0}can be any arbitrary function of the entropy and humidity.

*is obtained by differentiating (14):where*E

*ρ*

_{0}and adding the kinetic energy equations, yields the following equation for the conservation of the perturbation energy:Equation (14) defines the thermodynamic energy

*up to an additive function of entropy and water content, and the same holds for its derivatives*E

*μ*

_{υ}and

*μ*

_{d}are the chemical potentials for water vapor and dry air. The Maxwell relationships are obtained by comparing the second-order derivatives of enthalpy

*H*(

*S*,

*q*,

_{T}*p*). These relationships make it possible to rewrite the partial derivative of the buoyancy with respect to entropy:where Γ = −(∂

*T*/∂

*z*)

_{S,qT}(

*S*,

*q*

_{T},

*z*) is the adiabatic lapse rate, which is a function of the entropy

*S*, the total water content

*q*, and height

_{T}*z*. Similarly, the partial derivative the of buoyancy, respectively, to the total water mixing ratio isThe vertical dependencies of

_{0}to be given by,then its partial derivatives areUsing these expressions in (23) and (24) results in the functions

*T*and the chemical potential

*μ*and

_{υ}*μ*are evaluated at the reference state pressure

_{d}*p*

_{0}(

*z*). This can be compared with the thermodynamic relationship between enthalpy, entropy, and water content for moist air:The thermodynamic energy in the anelastic approximation

*H*at the reference state pressure

*p*

_{0}(

*z*). By construction, the thermodynamic energy of the reference is

*S*

_{0},

*q*,

_{T}*z*) = 0. Hence, the thermodynamic energy is equal to the enthalpy perturbation defined as the enthalpy difference between the parcel and the reference state:

*H*+

*gz*of the reference state is thus constant, and the thermodynamic energy (31) can then be written aswhere

*C*= −

*H*[

*S*

_{0},

*q*

_{T}_{0},

*p*

_{0}(

*z*

_{0})] −

*gz*

_{0}can be treated as an arbitrary integration constant. The thermodynamic energy in the anelastic approximation can thus be interpreted either as the perturbation enthalpy or as the moist static energy (enthalpy plus geopotential energy).

*δT*is the difference between these two temperatures, the relative temperature error scales aswhere

*U*is the velocity scale,

*c*is the speed of sound, and Mach = (

_{s}*u*/

*c*

_{s}) is the Mach number. Hence, the temperature in the anelastic approximation is close to the actual temperature as long as the Mach number remains low. A similar scaling also applies for the error in the chemical potential and other thermodynamic variables. As a low Mach number is a necessary condition for the validity of the anelastic approximation, the thermodynamic consistency is guaranteed in the general conditions of use of the anelastic approximation.

The equivalence between perturbation energy _{0}(*S*, *q _{T}*) =

*H*[

*S*,

*q*,

_{T}*p*

_{0}(

*z*)] in (14). Different choices for Φ

_{0}can be perfectly justified and would yield alternative forms of the perturbation energy that would differ from the enthalpy by a function of entropy and total water. These multiple forms of the perturbation energy can be thought of as capturing the thermodynamic behavior of other “fluids” for which the buoyancy would have the same dependency to the state variables

*S*and

*q*within the atmospheric domain, but with different state relationships outside the range of states considered in the atmospheric domain. The nonuniqueness of the perturbation energy is thus a reflection of the wide range of hypothetical fluids that would have the same behavior under the anelastic approximation.

_{T}## 4. Thermodynamic constraints on the buoyancy flux

In the previous section, it has been established that total energy is conserved under the anelastic approximation, with the thermodynamic energy given by the moist static energy evaluated at the reference state pressure *p*_{0}(*z*). In this section, we show that the thermodynamics constrain the buoyancy flux in a manner similar to the way the entropy budget constrains the work performed in the atmosphere.

*W*given by the domain integral of the buoyancy flux:For an atmosphere in statistical equilibrium, the work

*W*can be determined if the sources and sinks of entropy and water are known. Multiplying (15) by

*ρ*

_{0}and integrating over the entire domain yieldsafter taking advantage of the fact that the integral of

*ρ*

_{0}(

*D*

*DT*) over the atmospheric domain vanishes in statistical equilibrium. Equation (37) makes it possible to directly estimate the production of kinetic energy from the knowledge of the entropy and water sources. A similar expression can be obtained for a compressible atmosphere, and differs from (37) in that the actual density replaces the reference density and that the temperature and chemical potential are estimated at the actual pressure. Hence, understanding how the production of kinetic energy is constrained by the thermodynamics under the anelastic approximation should translate directly to a fully compressible case.

*Q*

_{in}at the lower boundary at temperature

*T*

_{in}and an energy sink

*Q*

_{out}at the upper boundary at temperature

*T*

_{out}. There is no source or sink of water here. The main difference from the traditional Rayleigh–Bénard problem is that a portion of the atmospheric domain is saturated. To keep the discussion as simple as possible, it is also assumed that the work

*W*is not dissipated but rather extracted from the system, which is equivalent to omitting the frictional heating and the corresponding irreversible entropy production. In this case, the only entropy sources and sinks in the system are the external energy sources and sinks. For the system to be in statistical equilibrium, the total entropy remains constant over time, which requires that (

*Q*

_{in}/

*T*

_{in}) + (

*Q*

_{out}/

*T*

_{out}) = 0. The work performed by the atmosphere is obtained by using (37):This corresponds to the work performed by a Carnot cycle acting between the temperatures

*T*

_{in}and

*T*

_{out}. The same expression applies for the work performed by a compressible atmosphere.

In deriving (41), the adiabatic lapse rate is approximated by that of the reference state. Given the fact that the lapse rate is discontinuous when a parcel goes through saturation, this is only accurate for levels at which the atmosphere is either entirely unsaturated or entirely saturated. This would be the case through most of the atmosphere for a stratus-topped boundary layer, but (41) cannot be applied to a problem such as cumulus convection where saturated and unsaturated air parcels coexist within a thick atmospheric layer. Nevertheless, the relation between the vertically integrated buoyancy flux and thermodynamic forcing (37) holds even in the case where the fluctuations of the entropy and total water content are large.

Equation (41) also indicates that the buoyancy flux in moist convection can be discontinuous. Indeed, below the cloud base where all air parcels are unsaturated, the buoyancy flux is equal to the entropy flux multiplied by the dry-adiabatic lapse rate; while above the cloud base, it is given by the moist-adiabatic lapse rate. The buoyancy flux exhibits a jump at the cloud base even when the flux of the entropy and water vapor are smooth. This is a consequence of the discontinuity in the partial derivatives of the state equation for moist air—in this case, the adiabatic lapse rate—as a parcel goes through saturation.

*q*is the specific humidity,

_{υ}*R*is the ideal constant for the water vapor, and

_{υ}*R*is the ideal gas constant for dry air. Such an expression for the buoyancy flux can be obtained directly from the definition of the buoyancy (9) and the ideal gas law. Similarly, the entropy flux can be approximated by

_{d}*C*the heat capacity of the moist air at constant pressure. Comparing this expression with (42) shows that the buoyancy flux here is equal to the entropy flux multiplied by the dry-adiabatic lapse rate Γ

_{p}_{d}= (

*g*/

*C*

_{p}).

*L*the latent heat of the vaporization and

_{υ}*q*

_{s}_{0}the saturation specific humidity in the reference state. For a saturated atmospheric layer, the entropy flux becomeswhile the buoyancy flux isThe ratio between these two quantities is indeed the moist-adiabatic lapse rate (2).

## 5. Thermodynamic energy as a prognostic variable

The anelastic framework discussed in the previous sections uses entropy and total water content as prognostic variables to describe the thermodynamic properties of moist air. Other thermodynamic state variables such as temperature, enthalpy, chemical potential, or specific humidity can be obtained using state relationships for moist air. Relying on moist entropy as a prognostic variable is not a new proposition. Ogura and Phillips (1962) mention their unsuccessful attempt to base their derivation of the anelastic approximation based on moist entropy. Ingersoll (2005) also uses the entropy of salty water as a prognostic variable for the anelastic approximation. Nevertheless, many authors (Ogura and Phillips 1962; Lipps and Hemler 1982; Durran 1989; Bannon 1996) have derived the anelastic approximations for different combinations of thermodynamic variables, most commonly temperature or potential temperature. Here, it is shown that the thermodynamic energy Eq. (15) can be used to obtain an alternative prognostic equation for these other prognostic variables.

*T*and

*μ*−

_{υ}*μ*, the last two terms on the right-hand side of (15) can be combined into

_{d}*TṠ*+ (

*μ*

_{υ}−

*μ*

_{d})

*q̇*

_{T}=

*Q*. Here,

*Q*is the external heating rate, which includes contributions from radiative heating, frictional heating, and molecular diffusion but not latent heat release. Equation (15) thus yields a prognostic equation for the thermodynamic energy:Given the definition of the thermodynamic energy (33), its rate of change is given bywhere

*C*is the rate at which the water vapor condenses (evaporation corresponding to a negative value for

*C*). The thermodynamic energy equation then becomesafter using

*ρ*

_{0}

*w*(

*g*+

*B*) =

*ρ*

_{0}

*g*(

*α*/

*α*

_{0})

*w*= −

*ρ*

_{0}

*α*(

*Dp*

_{0}/

*Dt*). Equation (50) corresponds to the first law of thermodynamics in which the pressure has been replaced by the reference state pressure

*p*

_{0}(

*z*).

Equation (50) has been derived from the tendency equations for entropy and total water. Any two of these three equations are sufficient to determine the thermodynamic evolution of the fluid. On the one hand, we can use entropy and total water as prognostic variables as in section 2 and 3, with the anelastic system conserving energy and the first law of thermodynamic is being implicitly obeyed. On the other hand, we can also use temperature and total water content as prognostic variables, with the temperature tendency obtained from Eq. (50). In this case, energy conservation is explicitly enforced. In addition, the system would have an implicit entropy equation that is consistent with the second law of thermodynamics. The entropy in the anelastic approximation would be defined as the entropy of parcel of moist air at the same temperature and water content, but estimated at the reference state pressure: *S* = *S*[*T*, *q _{T}*,

*p*

_{0}(

*z*)].

*θ*=

*π*

^{−1}

*T*, where

*π*= (

*p*/

*p*

_{ref})

^{(Rd/Cpd)}is the Exner function

*using the specific heat C*

_{pd}

*and gas constant R*. The quantity

_{d}for dry air*p*

_{ref}is an arbitrary pressure. While Eq. (51) is perfectly appropriate for a dry atmosphere, it fails to account for the effect of water vapor on both the heat capacity and gas constant, and as such is inconsistent with the first law of thermodynamics.

*T*= (

*p*

_{0}/

*p*

_{ref})

^{(Rd/Cpd)}

*θ*=

*π*

_{0}

*θ*. In the case of Lipps and Hemler (1982), the potential temperature tendency isUsing this expression in the potential temperature, Eq. (51) yieldsThis equation differs from the thermodynamic energy Eq. (48) in two respects. First, it uses the heat capacity of dry air

*C*

_{pd}rather than that of moist air

*C*. Second, the term multiplying the rate of change of pressure (

_{p}*R*

_{d}

*T*/

*p*

_{0}) =

*α*(

*p*

_{0},

*T*,

*q*

_{T}= 0) is equal to the specific volume of a dry air parcel. This omits the impacts of water on the density of moist air. Lipps and Hemler (1982) indeed observe that their formulation of the anelastic approximation does not conserve energy. This fact however should be attributed to their use of an incorrect form of the potential temperature equation for a moist atmosphere, and is not indicative of an intrinsic limitation of the anelastic approximation. Indeed, an anelastic model based on Eq. (50) and the full nonlinear equation of state for buoyancy not only conserves energy but also obeys the second law of thermodynamics.

## 6. Conclusions

This paper has investigated the use of a complex equation of state in conjunction with the anelastic approximation. The version of the anelastic approximation discussed here differs slightly from the previous work of Ogura and Phillips (1962), Dutton and Fitchl (1969), Lipps and Hemler (1982), and Durran (1989) in that it uses moist entropy and total water content as prognostic variables and has a reference state with constant entropy and total water content. Buoyancy is computed with the full nonlinear equation of state for the specific volume. This anelastic approximation remains accurate even for large fluctuations of the entropy and total water as long as the perturbations in the specific volume and total pressure remain small. While this paper has mainly focused on the behavior of a moist atmosphere, the results reported here are valid for a wide range of fluids with complex equations of state. In particular, the use of the anelastic approximation is justified for studying deep convection in a moist atmosphere in which the equation of state is highly nonlinear due to the discontinuity of its partial derivatives at saturation, or in the ocean, in which case the nonlinearity results from thermobaric effects.

The key finding here is that the anelastic approximation conserves energy, with the thermodynamic energy given by the difference between the enthalpy of the parcel at the reference pressure and the enthalpy of the reference state at the same level. Furthermore, the thermodynamic relationships between the thermodynamic energy, entropy, and total water content for an anelastic fluid are the same as for a fully compressible fluid, with the exception that the reference state pressure must be used instead of the total pressure. The thermodynamic consistency of the anelastic approximation holds as long as the Mach number of the flow is small. This result indicates that, even for complex fluids, anelastic models can be constructed so that they simultaneously obey the first and second laws of thermodynamics.

The thermodynamic energy equation can also serve as an alternative prognostic equation without affecting the thermodynamic consistency. In an anelastic model, temperature or enthalpy can replace moist entropy as a prognostic variable, with the evolution equation derived from the thermodynamic energy equation. Such a model would explicitly follow the first law of thermodynamics, while implicitly obeying the second law of thermodynamics, with the entropy given by the moist entropy of the parcel at the reference state pressure.

These findings ensure that anelastic models offer an appropriate framework within which to study the thermodynamic behavior of the atmosphere. In a compressible atmosphere, the work produced by the circulation can be determined from the energy and entropy budgets (Pauluis and Held 2002). Here, the same approach has been applied to an anelastic model to determine the buoyancy flux from the thermodynamic constraints imposed by the first and second laws of thermodynamics. In numerical models, ensuring mass and energy conservation is usually straightforward under the anelastic approximation, but presents several difficult challenges under the compressible Navier–Stokes equations. Hence, while there has been a strong push toward the development and use of fully compressible atmospheric codes, anelastic models still offer a very attractive platform, particularly for studying the maintenance of the atmospheric circulation over long time scales.

## Acknowledgments

I thank Bjorn Stevens, Dale Durran, Geoff Vallis, and two anonymous reviewers for their insightful comments on this manuscript. This work was supported by NSF Grant ATM-0545047.

## REFERENCES

Bannon, P., 1996: On the anelastic approximation for a compressible atmosphere.

,*J. Atmos. Sci.***53****,**3618–3628.Batchelor, G. K., 1953: The condition for dynamical similarity of motions in a frictionless perfect-gas atmosphere.

,*Quart. J. Roy. Meteor. Soc.***79****,**224–235.Boussinesq, J., 1903:

*Theorie Analytique de la Chaleur*. Vol. 2. Gauthier-Villars, 170–172.Durran, D. R., 1989: Improving the anelastic approximation.

,*J. Atmos. Sci.***46****,**1453–1461.Dutton, J., , and G. Fitchl, 1969: Approximate equations of motion for gases and liquids.

,*J. Atmos. Sci.***26****,**241–254.Emanuel, K. A., 1994:

*Atmospheric Convection*. Oxford University Press, 580 pp.Ingersoll, A., 2005: Boussinesq and anelastic approximations revisited: Potential energy release during thermobaric instability.

,*J. Phys. Oceanogr.***35****,**1359–1369.Lipps, F., , and R. Hemler, 1982: A scale analysis of deep moist convection and some related numerical calculations.

,*J. Atmos. Sci.***39****,**2192–2210.Mihaljan, J. M., 1962: A rigorous exposition of the Boussinesq approximations applicable to a layer of thin fluid.

,*J. Astrophys.***136****,**1126–1133.Oberbeck, A., 1879: Üeber die wärmleitung der flüssigkeiten bei berücksichtigung der strömungen infolge von temperaturdifferenzen. (On the thermal conduction of liquid taking into account flows due to temperature differences).

,*Ann. Phys. Chem. Neue Folge***7****,**271–292.Ogura, Y., , and N. A. Phillips, 1962: Scale analysis of deep and shallow convection in the atmosphere.

,*J. Atmos. Sci.***19****,**173–179.Pauluis, O., , and I. M. Held, 2002: Entropy budget of an atmosphere in radiative–convective equilibrium. Part I: Maximum work and frictional dissipation.

,*J. Atmos. Sci.***59****,**125–139.Peixoto, J. P., , and A. H. Oort, 1992:

*Physics of Climate*. AIP Press, 520 pp.Spiegel, E. A., , and G. Veronis, 1960: On the Boussinesq approximation for a compressible fluid.

,*Astrophys. J.***131****,**442–447.Vallis, G. K., 2006:

*Atmospheric and Oceanic Fluid Dynamics. Fundamentals and Large-Scale Circulation*. Cambridge University Press, 745 pp.

^{1}

Precipitation is made of large water droplets and ice crystals that are not in thermodynamic equilibrium with the surrounding air. Condensed water in the anelastic model discussed here is limited to small cloud droplets and crystals for which the thermodynamic equilibrium assumption is appropriate. A very rough treatment of rain and snow within this framework is to altogether remove the condensed water when it starts falling. A more comprehensive approach requires the addition of at least one state variable for precipitating water, but does not introduce any particular difficulty as far as the anelastic approximation is concerned.