## 1. Introduction

The geostrophic adjustment problem for rotating fluids is an initial-value problem where an unbalanced flow evolves toward a balanced state defined by the equilibrium between Coriolis and pressure gradient forces. The relationship between the initial and the final balanced states was first studied by Rossby (1938) and a review on the geostrophic adjustment problem can be found in the paper by Blumen (1972). Perhaps the simplest example of geostrophic adjustment is the one on an infinity channel, where a rotating shallow water fluid is initially at rest with a step profile in temperature. The flow evolves toward a final state where the velocity field balances the temperature gradient (Gill 1982). In this case, the final state is achieved after inertia–gravity waves propagate away from the initial disturbance, carrying energy away from the system. Therefore, the end state energy is less than the initial state by the amount of energy carried away by inertia–gravity waves. Importantly, the final state can be obtained analytically from the initial condition, and so can the energy partition between geostrophic and ageostrophic modes.

The atmosphere can be seen as a mixed fluid composed of dry air and water, in which fluctuations in velocity and temperature are associated with fluctuations in precipitation through latent heating due to water phase transitions. The atmospheric geostrophic adjustment to stationary and transient heat sources has been previously investigated (Gill 1980; Silva Dias et al. 1983) and it has been shown that the energy contribution to the transient adjustment by gravity waves depends on the spatial and temporal distribution of the forcing. Although these models provide relevant simulations, they have the limitation that they do not include the feedback between circulation and precipitation. To address this limitation, Emanuel et al. (1994), based on the quasi-equilibrium (QE) assumption (Arakawa and Schubert 1974), argue that the effect of convection on the large-scale circulation is to reduce the effective static stability of the atmosphere. Building on this approach, Frierson et al. (2004, hereafter FMP) developed an idealized framework to study the feedback between water vapor and large-scale circulation that allows for interactions between precipitating and nonprecipitating regions. The model by FMP has been previously used to study the propagation of precipitating regions (Stechmann and Majda 2006; Pauluis et al. 2008) and to investigate the propagation of convectively coupled equatorial waves (Dias and Pauluis 2009).

In the present work, we use the model by FMP to study the impacts of convective lifetime on the geostrophic adjustment process. Observational studies such as Betts (1986) suggest a convective adjustment time between 2 and 12 h, while the Coriolis parameter varies from 0 s^{−1} at the equator to 10^{−4} s^{−1} in the midlatitudes. Thus, it is important to understand the effects of the separation between these two time scales, both for the transient and for the balanced flow. To investigate the mechanisms of the adjustment we integrate the model by FMP on an *f* plane and in one space dimension, analyzing the adjustment for different ratios of precipitation and rotation time scales in comparison to the dry geostrophic adjustment.

Physically, the ratio between convective and geostrophic adjustment time can be seen in two ways. On the one hand, when the convective adjustment time is fixed, the comparison between time scales can be thought as the effects of moist convection on the geostrophic adjustment closer or farther from the equator. On the other hand, when the rotation frequency is fixed, the impact on the geostrophic adjustment of a shorter or longer convective adjustment time can be compared at a fixed distance from the equator.

The effects of convective lifetime on balanced flows is presented in the context of a geostrophic adjustment process in an infinite channel that is initially precipitating on one-half of it and dry on the other half. This paper is organized as follows. Section 2 reviews the governing equations and equations for the moist and dry potential vorticity (PV) are derived. Unlike the traditional baroclinic adjustment problem, because dry PV is not conserved, the final balanced flow cannot be determined from the initial PV field. To illustrate that, the initial-value problem is presented in section 3. In particular, the transient behavior and the final state are described for different ratios of convective and geostrophic adjustment time. In section 4, the mechanisms of the adjustment are analyzed and we derive approximations for the final state in two limiting cases, fast and slow convective adjustment time compared to the earth’s rotation time scale. In the last section we summarize the main results.

## 2. Modeling framework

*U*=

*U*

*x*,

*y*,

*t*) +

*u*(

*x*,

*y*,

*t*)Ψ

_{u}(

*z*), the meridional wind is given by

*V*=

*V*

*x*,

*y*,

*t*) +

*υ*(

*x*,

*y*,

*t*)Ψ

_{υ}(

*z*), the temperature

*T*=

*T*

*z*) +

*T*(

*x*,

*y*,

*t*)Ψ

_{T}(

*z*), the vertical velocity

*W*=

*w*(

*x*,

*y*,

*t*)Ψ

*(*

_{w}*z*), and the moisture

*T*

*z*) and

*Q*

*z*) are the horizontally uniform temperature and moisture profiles, and

*q*(

*x*,

*y*,

*t*) is the vertically averaged moisture. The vertical velocity is assumed to be zero at the surface and tropopause. The first baroclinic vertical structure is then given by Ψ

*(*

_{u}*z*) = Ψ

*(*

_{υ}*z*) = cos(

*z*) and Ψ

*(*

_{T}*z*) = Ψ

*(*

_{w}*z*) = sin(

*z*), where the height of the tropopause

*H*was nondimensionalized to

_{T}*π*. The vertical structure functions were chosen so that the vertical derivative of Ψ

*and Ψ*

_{T}*are equal to Ψ*

_{w}*and Ψ*

_{u}*. In particular, the continuity equations can be written as*

_{υ}**∇**

*= (∂*

_{H}*, ∂*

_{x}*) and (*

_{y}*x*,

*y*,

*z*) represent the zonal, meridional, and upward distances.

Although keeping the barotropic wind *U*

### a. Governing equations

*y*and the Coriolis parameter is constant. The flow remains independent of

*y*for all times; hence, we present only the one-dimensional model. In this case, in the absence of the barotropic flow and on the

*f*plane, the nondimensional flow equations from FMP are

*t̂*= 1/

*f*and the length scale is

*L̂*=

*c*/

_{d}*f*, with the characteristic speed given by

*c*

_{d}^{2}=

*π*

^{−1}

*H*∂

_{T}C_{p}_{z}

*T*

*gH*, where

_{T}*C*is the heat capacity at constant pressure. The diabatic forcing in the temperature Eq. (1c) represents the precipitation rate

_{p}*P*. To close the system, in the FMP formulation, an equation for the vertically integrated moisture is derived:

*C*

_{p}T_{0}/

*L*, where

_{υ}*T*

_{0}is the temperature scale and

*L*is the latent heat of vaporization. A key assumption in this derivation is that the gross moisture stratification

_{υ}*Q̃*depends only on the reference moisture profile

*Q*

*z*), which is obtained from a mean sounding. Hence, the moisture gradient depends only on the background moisture, and the term

*Q̃*∂

*= −*

_{x}u*Q̃w*represents advection of moisture. The parameter

*Q̃*is positive since typical moisture soundings fall exponentially with height, and it is less than one for a typical mean sounding. Here, we fix

*Q̃*= 0.9, corresponding to a background moisture gradient of 1.75 g kg

^{−1}km

^{−1}.

In comparison to the QTCM, FMP further simplifies the model by assuming that the vertical average of moisture is only advected by the barotropic flow. The main model parameters and scales are summarized on Table 1. Note that since we consider that *U*

*q̂*over a given time scale. In regions where moisture is below saturation, convection is inactive and precipitation vanishes. Mathematically,

*τ*is the convective adjustment time. For a given saturation value

_{c}*q̂*a new moisture variable can be defined as

*q*′ =

*q*−

*q̂*, and with this change of variables an equivalent set of equations is obtained where

*q̂*= 0. Hence, hereafter we consider only the case

*q̂*= 0; that is, negative values of moisture should be interpreted as the vertically integrated moisture content being below saturation, in which case precipitation is inactive. Conversely, in a region where

*q*> 0, precipitation is active. The problem remains nonlinear because of the transition between precipitating

*P*> 0 and the nonprecipitating regions where

*P*= 0, which represents the interface between moist and dry regions. Importantly, in FMP energy equations are derived including a quantity that can be interpreted as a moist available potential energy, and it is shown that the entire system is dissipative, independent of the convective adjustment time.

### b. Moist and dry potential vorticity

*ζ*= ∂

*and the dry and moist potential vorticity are defined by*

_{x}υ*x*and combining with Eqs. (1c) and (2) yields

*P*≥ 0, the amplitude of the dry PV increases when precipitation is active, whereas the moist PV is conserved.

### c. Geostrophic balanced state

*u*= 0,

_{F}*υ*= −∂

_{F}*/*

_{x}T_{F}*f*,

*P*= 0, and

_{F}*F*stands for variables at the final steady state,

*c*= 1 is the nondimensional gravity wave speed, and

_{d}*c*=

_{m}*P*= 0, the final moisture distribution

_{F}*q*has to be below saturation. In addition, time integration of Eqs. (5) and (6) gives a relation between the initial and final potential vorticity

_{F}*I*stands for variables at initial time and

*P*〉, whereas using the moist potential vorticity it requires knowledge of the final distribution of moisture

*q*. In fact, the two formulations are equivalent since time integration of (1c) and (2) yields

_{F}Notice that if there is no precipitation—that is, if the initial value of *q* is sufficiently negative—the final temperature can be determined from the initial data. Thus, in contrast to the dry geostrophic adjustment problem, Eq. (8a) suggests that we cannot determine the final state given only the initial state of the flow. To illustrate the effects of moist convection in the geostrophic adjustment process, in the next section we numerically integrate the governing Eqs. (1), given a particular choice of initial conditions, and in section 4 we discuss the physical mechanisms of the adjustment process.

## 3. Numerical simulations

Numerical solutions are obtained utilizing the nonoscillatory balanced scheme introduced by Khouider and Majda (2005a) and Khouider and Majda (2005b). To avoid reflection of gravity waves at the boundaries, we work on a very large domain with the adjustment region (∼10*L _{d}*) located at its center. Moreover, we stop integrations before any reflected wave could enter the region of interest. The domain corresponds to a latitude line of size 2

*L*(i.e.,

_{c}*x*∈ [−

*L*,

_{c}*L*] where

_{c}*L*is large compared to the Rossby radius

_{c}*L*=

_{d}*c*/

_{d}*f*).

Since the adjustment depends on both the initial state and the convective adjustment time *τ _{c}*, here we compare an initial value problem for

*τ*= 0.1/

_{c}*f*,

*τ*= 1/

_{c}*f*, and

*τ*= 10/

_{c}*f*. In this problem we initialized the model with a motionless atmosphere, with a uniform distribution of temperature and a step profile in the initial moisture content. In particular, to the left of the origin the domain is initially supersaturated and to the right it is unsaturated. The portion of the domain that is initially unsaturated remains unsaturated throughout the adjustment process (see Fig. 1a).

The propagation of gravity waves away from the initial discontinuity is illustrated in Fig. 2. It is noticeable that the propagation speed is slower to the left than to the right of the discontinuity when *τ _{c}* = 0.1/

*f*(left panels) than when

*τ*= 10/

_{c}*f*(right panels). The initial discontinuity in the moisture content triggers precipitation to the left of the origin and the initial precipitation is stronger when

*τ*is shorter (bottom panels). The total amount of precipitation shows a similar behavior in both cases, except that when

_{c}*τ*is shorter it precipitates more, as can be seen in Fig. 1a, which is associated with a final moisture distribution below saturation even in the initially moist portion of the channel, as shown in Fig. 1b.

_{c}The final balanced state is shown in Fig. 3, where it is compared to the dry geostrophic balance solution (thin dotted line). The two moist cases have in common a stronger jet due to the stronger temperature gradient and the loss of symmetry due to the asymmetric precipitation distribution illustrated in Fig. 1a. Precipitation acts as a heat source; thus, the final temperature is higher than in the dry case everywhere in the channel. However, since precipitation is never active on the right portion of the domain, temperature is smooth and comparable to the dry analytical solution for all values of *τ _{c}*.

Although Fig. 3 indicates that the latent heat initially released at the discontinuity sharpens the final temperature front independently of the relation between *τ _{c}* and

*f*, the initial burst of precipitation has a distinct effect depending on the ratio between the two time scales. When

*τ*is short, the final adjusted solution shows a bump in temperature near the maximum precipitation location and the amplitude of the bump increases with decreasing

_{c}*τ*. When the convective adjustment is slow compared to rotation, the total precipitation and the final solution are smoother. In fact, the final temperature is similar to the dry solution for a shorter length scale to the left of the origin. In the next section we turn to the analysis of the two limiting cases: short and long

_{c}*τ*.

_{c}## 4. Adjustment in the slow and fast convective limit

Based on observational studies that have suggested a convective adjustment time between 2 and 12 h (Betts 1986), the parameter *τ _{c}f* ranges from 0.1 near the equator to 4.0 in the midlatitudes, as can be seen in Fig. 4. Hence, in theory, the limiting behaviors

*τ*≫ 1 (slow convective limit) and

_{c}f*τ*≪ 1 (fast convective limit) described in the previous section can be observed. In this section, we first explore the similarities between the slow and fast convective limits. Next, in sections 4b and 4c, we separately analyze the distinct characteristics of the two limiting cases. The focus of the discussion is on both the transient flow (

_{c}f*t*≪ 1/

*f*) and the final balanced state. Remarkably, we find that in both limits the final state flow can be estimated from the initial flow.

### a. Similarities of the moist geostrophic adjusted flow

The results from the previous section indicate that the adjustment to the unbalanced initial moisture distribution triggers precipitation to the left of the edge of the initially supersaturated region. Latent heat is released in this region and temperature rises; hence, the final balanced flow exhibits a sharper temperature front, and a sharper jet, than the equivalent dry system would have. In this section, in order to clarify the physical mechanisms of the moist geostrophic adjustment process, we carry on a careful investigation of the adjustment shortly after the beginning of the integration time (*t* ≪ 1/*f* ).

*P*is uniform to the left of the origin and vanishes to the right of it (

*P*= 0 for

*x*> 0). Using linear approximation in space and time, we obtain a scaling for the amplitude of the temperature to the left of the origin:

*t*≪ 1/

*f*and

*q*(0,

*x*) = −

*q*

_{0}sign(

*x*).

*c*, the temperature jump from Eq. (10) is spread over a region of length Δ

_{d}*x*=

*c*Δ

_{d}*t*. This implies that the wind speed at the origin can be approximated by

*Q*

*z*) is the background moisture profile, cos(

*z*) is the structure function for the horizontal velocity

*u*, and the gross moisture stratification

*Q̃*is defined as the integral term in the second line. The atmospheric transport of water vapor is in the same direction as the low-level flow and is given by the term

*Q̃u*. This accounts for the fact that the low-level flow from the dry region to the moist regions contains more water vapor than the return flow in the upper troposphere. That is, the initial precipitation sets up a direct circulation in which the low-level air flows from the dry region into the moist region (

*u*< 0), which is associated with transport of water vapor from the dry region to the moist regions.

*c*Δ

_{d}*t*, 0):

*t*<

*τ*), this inflow of water vapor increases the water content in the moist regions. However, over longer time scales, precipitation has to prevail: the water vapor that was transported into the moist regions will be precipitated out. Equations (14) and (15) further indicate that the adjustment process brings moisture into the warm region (

_{c}*x*< 0), which leads to enhanced precipitation. This initial circulation can be better understood with the schematic picture shown in Fig. 5. The arrows indicate the convective cell where the fluid flows toward the initially precipitating region at the surface, leading to more precipitation. The flow ascends to the left of the origin (

*w*= −∂

*), turning back to the initially dry region. It then subsides, drying out the initially dry region.*

_{x}uAnother interesting feature of the moist adjustment is that in the region where *q*(*x*, 0) < 0 the amplitude of oscillation of all variables are bounded throughout the integration time. As a result, when moisture starts sufficiently below saturation, it never reaches saturation (as shown in Fig. 1). Therefore, the final balanced solution in the dry region corresponds to the dry analytical solution, except that the temperature at the interface (*x* = 0) depends on the moist convective adjustment in the precipitating region.

Although these estimates hold for any *τ _{c}*, the results from the previous section and the dependence on

*τ*in Eq. (15) suggest that the initially unbalanced flow evolves toward distinct final states in the slow or fast convective limits. This distinct behavior is discussed in sections 4b and 4c.

_{c}### b. The slow convective limit

Because the precipitation dissipation rate is related to *τ _{c}*, precipitation is active for a long period of time when the convective adjustment time is long. In this section, we show that, as a result, in the slow convective limit, the distribution of the total amount of precipitating water 〈

*P*〉 in the system can be estimated, as well as the final balanced flow. The argument is built on the fact that in this case the coupling between temperature and moisture is weak and after

*t*≫ 1/

*f*time, most gravity wave activity must have propagated away from the region of interest, which implies that the convergence term ∂

*is small. Consequently, once all the water in the system has precipitated, variations in moisture are also small (∂*

_{x}u*≈*

_{t}q*Q̃*∂

*) and the final equilibrated moisture distribution has to be close to saturation. That is,*

_{x}u*q*→ 0 when

_{F}*τ*→ ∞.

_{c}f*q*= 0, we can compute a solution in which the moist PV is conserved to the left of the origin and the dry PV is conserved to the right of it:

_{F}*L*=

_{d}*c*/

_{d}*f*is the dry Rossby radius,

*L*=

_{m}*c*/

_{m}*f*is the moist Rossby radius, and

*α*and

*β*are determined by matching the two solutions at the interface. For illustration, when

*T*+

_{I}*q*=

_{I}*q*

_{0}if

*x*< 0 and

*T*= 0 if

_{I}*x*≥ 0, we have

*α*= −[

*c*/(1 +

_{m}*c*)]

_{m}*q*

_{0}and

*β*= [1/(1 +

*c*)]

_{m}*q*

_{0}. This analytical solution will be referred as the moist-dry analytical solution. In particular, the moist-dry analytical solution shows that the temperature at the center of the channel is higher than it would be in the absence of precipitation,

*τ*.

_{c}### c. The fast convective limit

*τ*. The fast convective adjustment time in the absence of rotation has been discussed by FMP and Pauluis et al. (2008) and corresponds to the strict quasi-equilibrium (SQE) theory presented there. In SQE, moisture is instantaneously relaxed toward the saturation profile (in this case, to zero). That is, whenever moisture is above saturation,

_{c}*= (1 −*

_{t}T*Q̃*)∂

*when*

_{x}u*P*> 0, and ∂

*= ∂*

_{t}T*otherwise. The coupling between temperature and moisture tendencies is stronger and more likely to affect the geostrophic balanced solution since a precipitating region influences the dynamics of gravity waves. In particular, Pauluis et al. (2008) shows that gravity waves moving across precipitating and nonprecipitating regions can be treated as moving between regions of different refractive indices; therefore, waves are partially reflected and partially transmitted when they encounter the boundary between a moist and dry region. As previously mentioned, there is a clear distinction between the final states in the two limits; namely, in Fig. 3 a bump in temperature can be seen to the left of the initial discontinuity, which appears to be a result of the initial burst of precipitation in that region.*

_{x}u*P*

_{0}. Note that to satisfy the moist PV equation in the initially moist portion of the channel, the final temperature has to depend on both dry and moist Rossby radius, whereas in the dry area it depends only on the dry Rossby radius. That is,

*α*= −[

*L*

_{m}L_{d}^{2}/(

*L*

_{d}^{2}−

*L*

_{m}^{2})]

*P*

_{0}. The parameters

*λ*and

*β*are obtained assuming that the balanced temperature is continuous and differentiable at

*x*= 0 and are functions of

*P*

_{0}. Note that the temperature at the interface must be warmer than in the dry adjustment:

*P*

_{0}obtained in the previous section (when

*τ*→ ∞) results in

_{c}*λ*= 0. Therefore, this formulation is consistent with the final temperature obtained for large

*τ*. Interestingly, if

_{c}*P*

_{0}is large enough, the matching temperature at the interface will be above the temperature away from the origin (in this case,

*q*

_{0}) and the final temperature will exhibit the bump observed in the numerical experiment.

*P*

_{0}that yields

*τ*, the precipitating region moves at speed

_{c}*c*, as shown in Fig. 7 and illustrated in Fig. 8. Using the estimated value of

_{m}*P*

_{0}, we compute the final temperature distribution and total amount of precipitation, shown in Fig. 9. Notice the local maximum in temperature to the left of the origin. In the appendix it is shown that the extra amount of precipitation (in comparison to the slow convective adjustment) comes from the fact that since the precipitation region propagates westward, precipitation is not only generated by drying out the initially dry region, but also results from the contribution from the region between the moving precipitating front and the origin.

Notice that this analytical estimate does not account for the precipitation generated by gravity waves and thus it cannot explain the complete behavior of the balanced flow. For instance, the analytical solution exhibits a bump in temperature to the left of the origin; however, its amplitude is smaller and slightly shifted to the right in comparison to the numerical solution. Moreover, the total analytical precipitation underestimates the total numerical precipitation to the left of the origin and overestimates at the origin. Nevertheless, we showed that the basic features of the balanced flow are well represented, depending primarily on the initial adjustment process.

## 5. Conclusions

We presented a theoretical study of a geostrophic adjustment process in an idealized moist atmosphere where convection is represented by a Betts and Miller (1986) type of relaxation scheme. In particular, the case of a motionless flow evolving toward a balanced state was analyzed when only half of the domain is precipitating. A difference between the geostrophic adjustment in the dry and moist cases lies in that the moist adjustment process is associated with a water vapor transport from the dry regions to the moist regions that enhances the precipitation, sharpens the temperature gradient, and intensifies the jet. In this formulation, the gross moist stratification *Q̃* determines the moist phase speed of the waves and, therefore, the length scale of the adjustment *L _{m}* =

*f*. Thus, the temperature front sharpens and the jet intensifies with decreasing values of

*Q̃*.

A more fundamental difference rises from how gravity waves and the rotational flow interact in the two adjustment processes. On the one hand, in the dry case gravity waves do not affect the potential vorticity. This makes it possible to compute the final balanced state based solely on the initial potential vorticity distribution. On the other hand, in the moist case gravity waves can generate precipitation, which in turn affects the potential vorticity distribution. This means that in the moist case the fast dynamics associated with the propagation of gravity waves can impact the slow dynamics tied to evolution of potential vorticity. As this coupling occurs through precipitation, one expects that details of the representation of convective processes should affect geostrophic adjustment in a moist atmosphere. This is evident in our numerical simulations from section 3, which show that adjusted flow strongly depends on the convective relaxation time.

To further understand the final flow dependence on the convective adjustment time, in section 4 analytical approximations are derived in two limits: the slow convective limit with *τ _{c}f* ≫ 1 and the fast convective limit with

*τ*≪ 1. In the slow convective limit, the coupling through precipitation is weak, gravity waves propagate away from the origin, and the final solution converges to a solution similar to the analytical dry solution, but with a reduced length scale

_{c}f*L*in the initially supersaturated portion of the domain. In both cases precipitation is enhanced by the initial westward inflow, which dries out the initially dry region and moistens the initially moist region. In the slow convective limit, the right edge of the precipitating region is stationary and it is dissipated as gravity waves propagate away from the origin. In the fast case, the precipitating region propagates westward and more moisture is transported toward the precipitating region when compared to the slow case. The warming effect of this extra amount of precipitation produces a local maximum to the left of the origin in the final temperature distribution.

_{m}Our study emphasizes the fact that precipitation acts as a coupling mechanism between gravity waves and potential vorticity. In particular, the long-term evolution of the flow was shown to depend critically on the behavior of convection, which in our study could be characterized by a single nondimensional parameter *τ _{c}f*. While our analysis is focused on a highly idealized problem and omits several other processes, such as radiation, evaporation, and vertical wind shear, our theoretical analysis provides a relatively simple framework to study the physical mechanisms by which precipitation affects geostrophic flows.

## Acknowledgments

The comments of two reviewers are gratefully acknowledged, as they have considerably improved the clarity of this paper. Juliana Dias and Olivier Pauluis are supported by the NSF under Grant AGS-0545047.

## REFERENCES

Arakawa, A., and W. H. Schubert, 1974: Interaction of a cumulus cloud ensemble with the large-scale environment, Part I.

,*J. Atmos. Sci.***31****,**674–701.Betts, A. K., 1986: A new convective adjustment scheme. Part I: Observational and theoretical basis.

,*Quart. J. Roy. Meteor. Soc.***112****,**677–691.Betts, A. K., and M. J. Miller, 1986: A new convective adjustment scheme. Part II: Single column tests using GATE wave, BOMEX, ATEX and arctic air-mass data sets.

,*Quart. J. Roy. Meteor. Soc.***112****,**693–709.Blumen, W., 1972: Geostrophic adjustment.

,*Rev. Geophys. Space Phys.***10****,**485–528.Dias, J., and O. Pauluis, 2009: Convectively coupled waves propagating along an equatorial ITCZ.

,*J. Atmos. Sci.***66****,**2237–2255.Emanuel, K. A., J. D. Neelin, and C. S. Bretherton, 1994: On large-scale circulations in convecting atmospheres.

,*Quart. J. Roy. Meteor. Soc.***120****,**1111–1143.Frierson, D. M. W., A. Majda, and O. Pauluis, 2004: Large scale dynamics of precipitation fronts in the tropical atmosphere: A novel relaxation limit.

,*Commun. Math. Sci.***2****,**591–626.Gill, A. E., 1980: Some simple solutions for heat-induced tropical circulation.

,*Quart. J. Roy. Meteor. Soc.***106****,**447–462.Gill, A. E., 1982:

*Atmosphere–Ocean Dynamics*. Academic Press, 662 pp.Khouider, B., and A. Majda, 2005a: A non-oscillatory balanced scheme for an idealized tropical climate model. Part I: Algorithm and validation.

,*Theor. Comput. Fluid Dyn.***19****,**331–354.Khouider, B., and A. Majda, 2005b: A non-oscillatory balanced scheme for an idealized tropical climate model. Part II: Nonlinear coupling and moisture effects.

,*Theor. Comput. Fluid Dyn.***19****,**355–375.Neelin, J. D., and N. Zeng, 2000: A quasi-equilibrium tropical circulation model—Formulation.

,*J. Atmos. Sci.***57****,**1741–1766.Pauluis, O., A. Majda, and D. M. Frierson, 2008: Propagation, reflection and transmission of precipitation fronts in the tropical atmosphere.

,*Quart. J. Roy. Meteor. Soc.***124****,**913–930.Rossby, C-G., 1938: On the mutual adjustment of pressure and velocity distributions in certain simple current systems.

,*J. Mar. Res.***1****,**239–263.Silva Dias, P. L., W. H. Schubert, and M. DeMaria, 1983: Large-scale response of the tropical atmosphere to transient convection.

,*J. Atmos. Sci.***40****,**2689–2707.Stechmann, S., and A. Majda, 2006: The structure of precipitation fronts for finite relaxation time.

,*Theor. Comput. Fluid Dyn.***20****,**377–404.

## APPENDIX

### The Total Amount of Precipitation in the Fast Convective Limit

*τ*→ 0. The first assumption we make stands on the empirical evidence that the total amount of precipitation 〈

_{c}*P*〉(

*x*) decays similarly in the fast and slow convective adjustment process (see Fig. 1a). Since in the case of the slow convective adjustment it was shown that 〈

*P*〉(

*x*) ∝

*e*

^{−|x|/Lm}, we assume that in the fast case 〈

*P*〉(

*x*) = (

*P*

_{0}/

*L*)

_{m}*e*

^{x/Lm}if

*x*< 0, when

*τ*is small, and we define the total amount of precipitation generated during the geostrophic adjustment process by

_{c}*τ*infinitely small, moisture is relaxed toward saturation instantaneously; therefore, precipitation is generated and dissipated very rapidly. Since the initial zonal flow moves toward the moist portion of the domain, it must propagate at the moist gravity wave speed

_{c}*c*, converging at

_{m}*x*= −

*c*where it generates a localized burst of precipitation. Thus, the precipitating region must be inside a region defined by −

_{m}t*c*≤

_{m}t*x*≤ −

*c*+

_{m}t*δ*, where

*δ*is a small positive distance (it corresponds to the gray area in Fig. 8b). Because precipitation is completely dissipated before the flow is balanced, it must satisfy

*t*

_{0}will be determined later. Recall that in SQE,

*P*= −

*Q̃*∂

*; thus,*

_{x}u*U*

_{0}(

*t*) =

*u*(−

*c*−

_{m}t*δ*,

*t*), that is, the zonal flow at the right edge of the precipitating region (as the scheme in Fig. 8 indicates). Notice that to determine

*P*

_{0}, it is necessary to estimate both

*t*

_{0}and

*U*

_{0}. To determine

*U*

_{0}(

*t*), we first derive the momentum budget equations

*x*between the characteristics

*t*= −

*x*/

*c*and

_{m}*t*=

*x*/

*c*(as illustrated in Fig. 8b). Computing the zonal average of

_{d}*u*(

*x*,

*t*), we obtain

*U*

*c*+

_{d}*c*)

_{m}*tU*

_{0}(

*t*) and substitution in Eq. (A6) yields

*u*(

*x*,

*t*) → 0 as

*t*→ 0 is proportional to sin(

*ft*)/(

*c*+

_{d}*c*)

_{m}*t*. To determine

*U*

_{0}(0), we use Eq. (A4)

*t*→ 0, which gives

*U*

_{0}(0) = −Δ

*T*/(

*c*+

_{m}*c*). As a result

_{d}*U*

_{0}(

*t*) switches sign at

*t*=

_{P}*π*/

*f*; since

*P*> 0, we obtain

*T*=

*q*

_{0}. That is, using the estimate in (10) when Δ

*t*=

*τ*implies that the initial discontinuity in temperature with moisture initially at saturation and the initial discontinuity in moisture with initially uniform temperature are equivalent in SQE.

_{c}Initial adjustment for *q*(0, *x*) = −*q*_{0} sign(*x*) (*q*_{0} = 0.5): (top) temperature, (middle) zonal velocity, and (bottom) precipitation: (left) *τ _{c}* = 0.1/

*f*, (right)

*τ*= 10/

_{c}*f*. The thin dotted line corresponds to variables at

*t*= 0, the thick line to

*t*= 2/

*f*, and the dashed line to

*t*= 4/

*f*.

Citation: Journal of the Atmospheric Sciences 67, 9; 10.1175/2010JAS3405.1

Initial adjustment for *q*(0, *x*) = −*q*_{0} sign(*x*) (*q*_{0} = 0.5): (top) temperature, (middle) zonal velocity, and (bottom) precipitation: (left) *τ _{c}* = 0.1/

*f*, (right)

*τ*= 10/

_{c}*f*. The thin dotted line corresponds to variables at

*t*= 0, the thick line to

*t*= 2/

*f*, and the dashed line to

*t*= 4/

*f*.

Citation: Journal of the Atmospheric Sciences 67, 9; 10.1175/2010JAS3405.1

Initial adjustment for *q*(0, *x*) = −*q*_{0} sign(*x*) (*q*_{0} = 0.5): (top) temperature, (middle) zonal velocity, and (bottom) precipitation: (left) *τ _{c}* = 0.1/

*f*, (right)

*τ*= 10/

_{c}*f*. The thin dotted line corresponds to variables at

*t*= 0, the thick line to

*t*= 2/

*f*, and the dashed line to

*t*= 4/

*f*.

Citation: Journal of the Atmospheric Sciences 67, 9; 10.1175/2010JAS3405.1

Comparison between final balanced state when *q*(0, *x*) = −*q*_{0} sign(*x*) for *τ _{c}f* = 0.1, 1, and 10, and the dry case (thin dotted line). All variables where normalized such that Δ

*q*(

*x*, 0) = 1.

Citation: Journal of the Atmospheric Sciences 67, 9; 10.1175/2010JAS3405.1

Comparison between final balanced state when *q*(0, *x*) = −*q*_{0} sign(*x*) for *τ _{c}f* = 0.1, 1, and 10, and the dry case (thin dotted line). All variables where normalized such that Δ

*q*(

*x*, 0) = 1.

Citation: Journal of the Atmospheric Sciences 67, 9; 10.1175/2010JAS3405.1

Comparison between final balanced state when *q*(0, *x*) = −*q*_{0} sign(*x*) for *τ _{c}f* = 0.1, 1, and 10, and the dry case (thin dotted line). All variables where normalized such that Δ

*q*(

*x*, 0) = 1.

Citation: Journal of the Atmospheric Sciences 67, 9; 10.1175/2010JAS3405.1

Scale *τ _{c}f* as a function of latitude for two fixed values of

*τ*. The top line corresponds to

_{c}*τ*= 12 h and the bottom line to

_{c}*τ*= 2 h.

_{c}Citation: Journal of the Atmospheric Sciences 67, 9; 10.1175/2010JAS3405.1

Scale *τ _{c}f* as a function of latitude for two fixed values of

*τ*. The top line corresponds to

_{c}*τ*= 12 h and the bottom line to

_{c}*τ*= 2 h.

_{c}Citation: Journal of the Atmospheric Sciences 67, 9; 10.1175/2010JAS3405.1

Scale *τ _{c}f* as a function of latitude for two fixed values of

*τ*. The top line corresponds to

_{c}*τ*= 12 h and the bottom line to

_{c}*τ*= 2 h.

_{c}Citation: Journal of the Atmospheric Sciences 67, 9; 10.1175/2010JAS3405.1

Schematic representation of the initial adjustment process. The light gray region represents the initially supersaturated region and the dark gray area represents the precipitating region due to the initial inflow toward the moist region. The arrows represent the convective cell generated by the initial unbalanced moisture distribution. On the far right side the vertical structure function for *u* and *w* are displayed.

Citation: Journal of the Atmospheric Sciences 67, 9; 10.1175/2010JAS3405.1

Schematic representation of the initial adjustment process. The light gray region represents the initially supersaturated region and the dark gray area represents the precipitating region due to the initial inflow toward the moist region. The arrows represent the convective cell generated by the initial unbalanced moisture distribution. On the far right side the vertical structure function for *u* and *w* are displayed.

Citation: Journal of the Atmospheric Sciences 67, 9; 10.1175/2010JAS3405.1

Schematic representation of the initial adjustment process. The light gray region represents the initially supersaturated region and the dark gray area represents the precipitating region due to the initial inflow toward the moist region. The arrows represent the convective cell generated by the initial unbalanced moisture distribution. On the far right side the vertical structure function for *u* and *w* are displayed.

Citation: Journal of the Atmospheric Sciences 67, 9; 10.1175/2010JAS3405.1

Comparison between mixed analytical solution (solid line) and numerical solution for large *τ _{c}* (dashed and dotted lines).

Citation: Journal of the Atmospheric Sciences 67, 9; 10.1175/2010JAS3405.1

Comparison between mixed analytical solution (solid line) and numerical solution for large *τ _{c}* (dashed and dotted lines).

Citation: Journal of the Atmospheric Sciences 67, 9; 10.1175/2010JAS3405.1

Comparison between mixed analytical solution (solid line) and numerical solution for large *τ _{c}* (dashed and dotted lines).

Citation: Journal of the Atmospheric Sciences 67, 9; 10.1175/2010JAS3405.1

(top) Time evolution of precipitation for short *τ _{c}*. The dashed lines correspond to the dry (

*x*=

*c*) and moist (

_{d}t*x*=

*c*) characteristics, and the horizontal line corresponds to

_{m}t*t*=

*π*/

*f*. (bottom) Three snapshots (at

*t*= 0.5/

*f*, 1/

*f*, and

*π*/

*f*) of the precipitation amplitude.

Citation: Journal of the Atmospheric Sciences 67, 9; 10.1175/2010JAS3405.1

(top) Time evolution of precipitation for short *τ _{c}*. The dashed lines correspond to the dry (

*x*=

*c*) and moist (

_{d}t*x*=

*c*) characteristics, and the horizontal line corresponds to

_{m}t*t*=

*π*/

*f*. (bottom) Three snapshots (at

*t*= 0.5/

*f*, 1/

*f*, and

*π*/

*f*) of the precipitation amplitude.

Citation: Journal of the Atmospheric Sciences 67, 9; 10.1175/2010JAS3405.1

(top) Time evolution of precipitation for short *τ _{c}*. The dashed lines correspond to the dry (

*x*=

*c*) and moist (

_{d}t*x*=

*c*) characteristics, and the horizontal line corresponds to

_{m}t*t*=

*π*/

*f*. (bottom) Three snapshots (at

*t*= 0.5/

*f*, 1/

*f*, and

*π*/

*f*) of the precipitation amplitude.

Citation: Journal of the Atmospheric Sciences 67, 9; 10.1175/2010JAS3405.1

Schematic representation of the precipitating region when *τ _{c}* ∼ 0, where

*U*

_{0}(

*t*) is the zonal flow at its edge at time

*t*. (a) Precipitating region at time

*t*in the (

*x*,

*z*) plane. (b) Precipitating region localized at

*x*= −

*c*in the

_{m}t*x*–

*t*plane (schematic description of top panel in Fig. 7).

Citation: Journal of the Atmospheric Sciences 67, 9; 10.1175/2010JAS3405.1

Schematic representation of the precipitating region when *τ _{c}* ∼ 0, where

*U*

_{0}(

*t*) is the zonal flow at its edge at time

*t*. (a) Precipitating region at time

*t*in the (

*x*,

*z*) plane. (b) Precipitating region localized at

*x*= −

*c*in the

_{m}t*x*–

*t*plane (schematic description of top panel in Fig. 7).

Citation: Journal of the Atmospheric Sciences 67, 9; 10.1175/2010JAS3405.1

Schematic representation of the precipitating region when *τ _{c}* ∼ 0, where

*U*

_{0}(

*t*) is the zonal flow at its edge at time

*t*. (a) Precipitating region at time

*t*in the (

*x*,

*z*) plane. (b) Precipitating region localized at

*x*= −

*c*in the

_{m}t*x*–

*t*plane (schematic description of top panel in Fig. 7).

Citation: Journal of the Atmospheric Sciences 67, 9; 10.1175/2010JAS3405.1

(top) Comparison between the numerical (solid line) and the estimated (dashed line) total precipitation for small *τ _{c}* final temperature. (bottom) As in (top), but for the final temperature.

Citation: Journal of the Atmospheric Sciences 67, 9; 10.1175/2010JAS3405.1

(top) Comparison between the numerical (solid line) and the estimated (dashed line) total precipitation for small *τ _{c}* final temperature. (bottom) As in (top), but for the final temperature.

Citation: Journal of the Atmospheric Sciences 67, 9; 10.1175/2010JAS3405.1

(top) Comparison between the numerical (solid line) and the estimated (dashed line) total precipitation for small *τ _{c}* final temperature. (bottom) As in (top), but for the final temperature.

Citation: Journal of the Atmospheric Sciences 67, 9; 10.1175/2010JAS3405.1

Model parameters and scales.