In this paper, we describe the construction of an efficient probabilistic parameterization that could be used in a coarse-resolution numerical model in which the variation of moisture is not properly resolved. An Eulerian model using a coarse-grained field on a grid cannot properly resolve regions of saturation—in which condensation occurs—that are smaller than the grid boxes. Thus, in the absence of a parameterization scheme, either the grid box must become saturated or condensation will be underestimated. On the other hand, in a stochastic Lagrangian model of moisture transport, trajectories of parcels tagged with humidity variables are tracked, and small-scale moisture variability can be retained; however, explicitly implementing such a scheme in a global model would be computationally prohibitive. One way to introduce subgrid-scale saturation into an Eulerian model is to assume the humidity within a grid box has a probability distribution. To close the problem, this distribution is conventionally determined by relating the required subgrid-scale properties of the flow to the grid-scale properties using a turbulence closure. Here, instead, we determine an assumed probability distribution by using the statistical moments from a stochastic Lagrangian version of the system. The stochastic system is governed by a Fokker–Planck equation, and we use that, rather than explicitly following the moisture parcels, to determine the parameters of the assumed distribution. We are thus able to parameterize subgrid-scale condensation in an Eulerian model in a computationally efficient and theoretically well-founded way. In two idealized advection–condensation problems, we show that a coarse Eulerian model with the subgrid parameterization is well able to mimic its Lagrangian counterpart.
Water vapor is carried around as a tracer, normally condensing when the vapor pressure reaches the saturation value given by the Clausius–Clapeyron relation. (Strictly, the Clausius–Clapeyron relation determines the vapor pressure for water, or other condensible, in thermodynamic equilibrium with liquid water or ice, and if neither are present, then the water vapor can become supersaturated, but for the purposes of this paper, we will assume condensation occurs upon saturation.) Furthermore, condensation normally occurs quickly upon saturation, and it is common in the atmospheric sciences to refer to the “fast condensation limit.” The liquid water produced by condensation may then, in simple models, be assumed to fall to the ground, or more realistically, it may be advected by the flow before the drops coalesce and fall as rain or form ice, as determined in models by more or less complicated microphysical parameterization schemes.
Putting aside the complications of microphysical effects, the simplest advection–condensation model consistent with this picture postulates that as an air parcel is advected by the large-scale wind, and away from evaporation sources, the moisture content of the parcel remains constant except when it exceeds the local saturation limit, at which point the excessive water vapor is removed by condensation. Previous work (e.g., Salathé and Hartmann 1997; Pierrehumbert and Roca 1998; Galewsky et al. 2005; Dessler and Minschwaner 2007) has applied this idea to reconstruct large-scale features of the atmospheric moisture distribution. In these studies, the trajectory of a parcel is traced backward in time to the location where the parcel is last saturated, for example, the point at which it last encounters convection or the lower boundary layer. The specific humidity of the parcel at its present location is then given by the minimum saturation specific humidity encountered along the trajectory. The success of these studies highlights the Lagrangian nature of the large-scale transport and condensation of atmospheric moisture. Indeed, Pierrehumbert et al. (2007) suggested that the proper approach to represent moisture transport in climate models is to take the stochastic Lagrangian viewpoint whereby the fluctuations in the trajectories of moist parcels are parameterized by random processes.
Over the past few decades, stochastic Lagrangian models, which describe the trajectories of air parcels using a model of random velocity, have been developed to study turbulent transport in the atmosphere (Wilson and Sawford 1996; Thomson and Wilson 2013). Under a Markov assumption, the parcel position and velocity are random variables satisfying some stochastic differential equations. Equivalently, the model can also be specified by a Fokker–Planck equation that governs the joint probability density function (PDF) of position and velocity. Applying this approach to moisture transport, various studies (O’Gorman and Schneider 2006; Pierrehumbert et al. 2007; Sukhatme and Young 2011; Beucler 2016; Tsang and Vanneste 2017) have investigated theoretically the advection–condensation of water vapor by evolving an ensemble of particles, each carrying its own set of dynamical and thermodynamical variables obeying stochastic model equations. Whereas the stochastic Lagrangian description does have the advantage of, in principle, retaining local fluctuations at small scales, it also comes with a high computational cost—it is simply impractical to carry around a very large number of Lagrangian particles representing moist air parcels. One possible way to address this problem is to use a hybrid parcel-in-cell method (Dritschel et al. 2018), but below, we will describe a qualitatively different approach, in which the resulting equations are Eulerian (and so can be efficiently solved) but the underlying parameterization is explicitly based on a Lagrangian description.
The conventional practice is to represent atmospheric water vapor as a coarse-grained field on a numerical grid, writing the equations of motion in the Eulerian form as a partial differential equation (PDE), for example,
In this equation, u is velocity, q is specific humidity, S is a moisture source, C represents the effects of condensation, and D(x, t) is a diffusivity. The condensation term is zero until saturation occurs. In reality, D would be the molecular diffusivity and is very small indeed, so that the specific humidity of an unsaturated parcel is essentially conserved. However, in a coarse-resolution model—such as a climate model with a horizontal resolution measured in kilometers—D is often a parameterized diffusivity much larger than the molecular one. It also cannot be small for numerical reasons. (A semi-implicit, semi-Lagrangian scheme may be stable at low resolution without a high explicit diffusivity, but these methods can also be diffusive or inaccurate.) Furthermore, if condensation is only allowed to occur at saturation, then the effects of diffusion are in many circumstances such as to make large regions saturated, producing moisture fields noticeably different from a Lagrangian model (Pierrehumbert et al. 2007; Vallis 2017, chapter 18). In a climate model with a moisture equation similar to (1), condensation and rainfall will only occur when a grid box is entirely saturated. This has long been recognized to be in many circumstances quite unrealistic (e.g., Sommeria and Deardorff 1977), and because of the strong dependence of the absorption of outgoing longwave radiation on water content, such a misrepresentation can be especially significant in the modeling of Earth’s radiation budget.
The problem with the Eulerian approach, as noted by Pierrehumbert et al. (2007), is that the coarse graining that is in practice required does not commute with the highly nonlinear condensation process. Our first goal is in fact to demonstrate theoretically how this causes an Eulerian model without subgrid-scale condensation to produce large regions of saturation compared to its Lagrangian counterpart. A possible solution to such problems, sometimes used in cloud modeling (Tompkins 2002; Jakob and Miller 2002), is to suppose that the specific humidity (and possibly other thermodynamical variables) inside a given grid box is not single valued but has a probability distribution, thus introducing local fluctuations into the system. Then part of the box may be saturated even though the average specific humidity over the box is less than the saturation limit, and some fraction of the water vapor content, as determined by the probability function, may then be removed by condensation. A probabilistic parameterization of subgrid-scale condensation along these lines was proposed by Sommeria and Deardorff (1977) and Mellor (1977) to model moist convection in the boundary layer. They assumed the total mixing ratio and the liquid potential temperature have a joint Gaussian PDF and determined the cloud fraction within a grid cell from such a PDF. Bougeault (1981) later used this “assumed PDF” method with several different PDF shapes to model the trade wind cumulus layer. Since then, numerous variations have been developed and employed in atmospheric numerical models. For example, one of the cloud schemes in the Met Office Unified Model is the Smith (1990) scheme, which uses a triangular PDF (Wilson et al. 2008); a somewhat more complicated scheme has been used at ECMWF (Tiedtke 1993); and various other, sometimes still more complicated (and computationally intensive) schemes have been proposed (e.g., Lappen and Randall 2001; Tompkins 2002; Golaz et al. 2002; Kuwano-Yoshida et al. 2010; Bogenschutz and Krueger 2013). A crucial step in these probabilistic schemes is to determine the parameters of the prescribed PDF, such as width and skewness. This is often done by linking the PDF parameters to various eddy fluxes or correlation functions, and turbulence closure models are then used to predict these correlations from the resolved scales. The difficulty with these approaches is, of course, that the parameterization is only as good as the turbulence closure it is based upon.
Evidently, then, both Lagrangian and Eulerian approaches have advantages and shortcomings—the former is accurate but impractical, and the latter is practical but less accurate, with the contrast stemming from the fundamental differences in the representation of particle motion and condensation of the two formulations. In this paper, we seek to combine these two approaches. The idea is to use information extracted from a corresponding stochastic Lagrangian model to derive a parameterized Eulerian model that can produce similar results to the stochastic Lagrangian model but at a fraction of the computational cost. We aim to achieve two goals. The first is to provide a sound theoretical basis to the heuristic probabilistic schemes that are in common use. The second is to describe a systematic way whereby a probabilistic parameterization for the condensation in an Eulerian model may be derived, for example, to provide a C in the Eulerian equation, (1).
The premise of our method is that the small-scale velocity V′ of a moist parcel can be modeled as a random process. We may then represent the moist dynamics by a stochastic Lagrangian model in which an ensemble of moist parcels is advected by the velocity V + V′, where V is the deterministic large-scale parcel velocity. Since V′ is random, at each location and time, the stochastic system produces a PDF of the humidity, . Now, it is expensive to obtain by performing a Monte Carlo simulation of the stochastic differential equations or by solving the high-dimensional Fokker–Planck equation governing . Instead, we use an assumed PDF as surrogate for and require the moments of to match those of derived from the Fokker–Planck equation. The so determined will then be used in a probabilistic parameterization of C for an Eulerian model such as (1). Thus, our scheme involves two steps. First, an appropriate stochastic Lagrangian model must be constructed, and second, the Fokker–Planck equation—as an alternative to turbulence closures—must be used to derive parameters for an assumed PDF. We carry out this procedure in two idealized advection–condensation problems and show that a coarse Eulerian model with the subgrid parameterization is, in fact, well able to mimic its Lagrangian counterpart. Because of the idealized nature of these problems, we are able to solve the Lagrangian model directly, by Monte Carlo simulations of moist particles advected by a large-scale field and a random component, and so provide a true test of the methodology.
The paper is organized as follows. In section 2, we present the basic Lagrangian and Eulerian methodologies using a model of moisture transport in an overturning cell and show that an Eulerian model tends to produce saturated air. Section 3 gives the details of probabilistic parameterization of condensation. We then compare results from Eulerian models with and without parameterization to those of Lagrangian models for a steady flow in section 4 and for an unsteady flow in section 5. In section 6, we discuss the use of an underlying stochastic Lagrangian model to parameterize condensation in coarse-grained atmospheric models, and we conclude the paper in section 7.
2. Lagrangian particles versus Eulerian fields
a. Advection–condensation in an overturning cell
We consider the advection of moist air in a square domain on the x–y plane. Condensation occurs as water vapor is transported by a prescribed velocity through a saturation specific humidity field . We assume the velocity has an incompressible large-scale component and a turbulent component at the small scales. In this section, as a crude model with some similarities to the Hadley cell, we take as a steady overturning flow with streamfunction
Figure 1 shows the streamlines of u in a schematic of the system. We assume varies only with the altitude y and is independent of time. Specifically, we assume a linear temperature profile in y:
together with gives
and we define
Here, we set and . Hence, and . We assume there is an evaporation source S located at the bottom boundary to maintain the specific humidity along at .
The interplay between large-scale coherent flow, small-scale turbulence, and condensation in this model roughly reproduces several interesting features of the atmosphere (Tsang and Vanneste 2017): a humid bottom boundary layer that resembles the planetary boundary layer, a narrow region of intense condensation along reminiscent of the tropics, and a relative humidity minimum at the center of the cell.
b. Deterministic coarse-grained field formulation
For a deterministic Eulerian formulation of the advection–condensation problem described above, the specific humidity is represented by a coarse-grained field whose time evolution is governed by the PDE:
Above, u is the large-scale velocity, and unresolved small-scale turbulence is represented by the diffusion term with constant eddy diffusivity . In accord with the advection–condensation paradigm, molecular diffusion is assumed to be negligible. The condensation C may be written as
where is the condensation time scale and is the Heaviside step function. For most of this paper, we employ the rapid-condensation limit of and implement C as a rule to prevent supersaturation:
The source S is implemented as a boundary condition:
At the other boundaries, we have the no-flux conditions:
We solve (7) for the field using the split-step approach. Given at time , we obtain an intermediate (supersaturated) field by time stepping forward the advection–diffusion equation
to . (Here and elsewhere in the paper, we write the equations with partial derivatives with respect to time, such as , with the understanding that the procedure takes us from to .) We use a semi-Lagrangian scheme for the advection and the alternative direction implicit method for the diffusion. We then carry out the condensation using (9) to produce the moisture field at time :
Figure 2c shows the steady-state relative humidity field,
at a large time t from a 5132 simulation with .
c. Stochastic particle formulation
We now turn to a stochastic Lagrangian model of the system. Inside the square domain, the moist air is now represented by an ensemble of air parcels. The domain boundaries are reflective. The parcels are initially uniformly distributed over the domain and will remain so for an incompressible advecting flow.
Let be the position of a parcel and Q be its specific humidity. Consider as random variables, the advection–condensation of each moist parcel is described by the following set of stochastic differential equations:
The resolved large-scale velocity u is once again given by (2); and are Wiener processes, and is the associated Brownian diffusivity. Thus, the small-scale turbulent velocity is modeled as white noise, denoted by . In other words, the turbulent velocity of each parcel is a Gaussian random variable at any instance of time t and is uncorrelated in time with correlation function:
Exchange of moisture between parcels, which may be important in some situations (Haynes and Anglade 1997), is not included in this model, and each parcel evolves independently. For finite condensation rate , C is given by
and in the limit of , we have
The action of the source S at is that it resets the specific humidity of air parcels to the local saturation value upon hitting the bottom boundary. For a detailed analysis of this stochastic system, we refer the readers to Tsang and Vanneste (2017).
We perform Monte Carlo simulation of (15) using the Euler–Murayama method (Higham 2001). At , 106 saturated parcels are uniformly distributed over the domain. Figure 2a shows a snapshot of a subset of the parcels after the system has reached a statistically steady state. The color indicates the relative humidity of each parcel:
To visualize the spatial distribution of moisture over the domain, we construct a bin-averaged field from the Monte Carlo data by dividing the domain into a uniform gird of square bins. We then average over all parcels inside the bin centered at to obtain . As the velocity u is steady and we are interested in the statistically steady distribution, we further average over time to obtain . Figure 2b shows corresponding to the simulation in Fig. 2a; 5132 bins have been used. An interpretation of this averaging procedure is that many parcels with different contribute to a single observation of taken over a small area about .
d. Noncommutation between condensation and coarse graining
Let us now compare the relative humidity field calculated from the two formulations. As shown clearly in Figs. 2b and 2c, the Lagrangian and the Eulerian models produce starkly different results. The Eulerian model has the unrealistic feature that a large part of the domain is fully saturated with r = 1. Figure 2e plots the difference in the relative humidity field from the two models. Unsurprisingly, the largest discrepancy occurs in the rising half of the cellular flow where most of the condensation happens. Generally, the saturated region in the Eulerian model will shrink as decreases (e.g., cf. Figs. 2c and 7a). However, regardless of the value of , the boundary at x = 0 will remain saturated. This is fundamentally different from the results of the Lagrangian model. Pierrehumbert et al. (2007) had observed similar behavior in simple one-dimensional models and attributed it to the loss of local fluctuation in a coarse-grained field representation, and Vallis (2017) qualitatively described similar behavior in a two-dimensional model. Here, we investigate this effect quantitatively in the two-dimensional case.
Figure 2a clearly shows parcels with a broad range of specific humidity coexist within a small area. When the moisture distribution is represented by a coarse-grained field, such subgrid-scale fluctuation is averaged out, leading the system to bias toward saturation. Mathematically, this is because the condensation process and the coarse-graining process do not commute. To elucidate, we examine from a theoretical viewpoint how one goes from Fig. 2a to Fig. 2b. To this end, it is more convenient to momentarily revert to a small but nonzero condensation time . For an ensemble of parcels described by the random variables obeying the stochastic differential equations, (15), the joint PDF of specific humidity and position satisfies the Fokker–Planck equation (Pavliotis 2014),
supplemented by appropriate boundary conditions in the domain . Above, we have used the incompressibility condition and
The mean specific humidity at a given position is the conditional expectation value,
where is the conditional probability density for a parcel to have specific humidity given it is located at . So the bin-averaged field in Fig. 2b is a numerical approximation to .
We now derive the evolution equation of . By the definition of conditional PDF, in (23) is related to P by
is the marginal PDF that gives the probability that a parcel is located at regardless of its specific humidity. Integrating (21) over yields
Note that the boundary term involving C from the integration vanishes. This is because by (22) and also as , no parcel can have inside the domain: for all . It then follows from (21) and (26) that satisfies
Since the parcels are uniformly distributed at t = 0 in our setup, (26) implies for all t; hence, the parcels remain uniformly distributed, and the term involving in (27) vanishes. This simply means we are concerned with a constant air density.
Multiplying (27) by and integrating over , we finally get the equation for :
On the other hand, the governing equation, (7), of in the Eulerian formulation for nonzero reads
Comparing (28) with (29), we see that the differences in and q stem from the condensation term. In (28), condensation for each individual parcel is considered before their contributions to are added up. Thus, local fluctuations are accounted for. In (29), only the coarse-grained value q is available, and condensation only happens when , causing the system to retain more moisture as seen in Fig. 2c. Figure 3 illustrates this noncommutation between condensation and coarse graining pictorially with an example.
3. Probabilistic parameterization of condensation
a. An effective condensation
We have seen in previous sections that modeling water vapor distribution using a coarse-grained field is prone to producing saturation. On the other hand, the Lagrangian approach is able to produce more realistic results, albeit with higher computational cost, by accounting for the effects of subgrid-scale moisture fluctuation on condensation. Here, we ask the question, If we regard the Lagrangian model as “truth,” how do we construct an Eulerian PDE-based model that might be used in its place to give similar results? Comparing (28) and (29) suggests naturally the answer is to replace the condensation term in the Eulerian equation, (29), by an effective condensation:
where is an approximation to the “true” conditional PDF in the Lagrangian model. Equation (30) resembles the formula for liquid water content in a conventional probabilistic subgrid-scale cloud scheme (Sommeria and Deardorff 1977). To specify , we take the assumed PDF approach by assuming a functional form for that contains a small number of parameters. These parameters are then determined by matching the moments of to those of up to a certain order. Note that is not governed by an evolution equation and there is the freedom to assume different functional forms at different times. We shall explain the detailed procedure through examples in sections 4 and 5.
b. Numerical implementation
To adopt the above representation of condensation into our numerical framework, first recall from the discussion around (12) that we employ the split-step algorithm and first solve the advection–diffusion step to obtain the intermediate field . This is then followed by solving the condensation step:
In the limit , we can consistently set in the effective condensation, (30), where is the time step of the simulation. Assume at the end of the advection–diffusion step is known, and denote it by . Discretizing (31) in time with C given by (30) leads to the condensation formula:
which gives the value of the specific humidity at the end of one full time step.
Before we proceed further, we give a physical interpretation to (32) and also set the stage for specifying in the next sections. The idea is to interpret the value of the specific humidity at a given grid point after the advection–diffusion step as the mean from an ensemble of parcels with specific humidity distribution ; that is,
Note that some of these imagined parcels can have specific humidity higher than even if . This is illustrated in the top panel of Fig. 4. Next, we carry out rapid condensation on this ensemble to reduce the specific humidity of all supersaturated parcels to . The distribution after condensation is
where α is fixed by the normalization condition . Figure 4 shows a schematic of this parameterized condensation. Finally, the specific humidity field at time is given by
from which (32) follows.
4. A steady overturning flow
In our first example of applying the condensation parameterization, we use the system introduced in section 2 where moist air in a square cell is advected by the steady overturning flow u in (2). The coarse-grained specific humidity field in the Eulerian formulation is governed by the PDE (7). The fairly large diffusivity of magnifies the susceptibility to saturation in the coarse-grained model and puts the condensation parameterization to stringent test.
Our first task is to choose an ansatz for . For this particular setup, subsidence of dry parcels from the upper part of the domain significantly affect the distribution of specific humidity. The driest parcels are created at the top boundary. They roam through the domain and maintain their dryness of until they hit the localized moisture source at the bottom boundary. As a consequence, we expect to be composed of a dry spike of amplitude β (Sukhatme and Young 2011; Tsang and Vanneste 2017) and a continuous part :
In part for simplicity and in part because we expect the distribution of specific humidity to be smooth over a range of values (as we show later), we assume a top-hat shape for at all times. Referring to the top panel of Fig. 4, is centered at a with width 2σ. Normalization condition gives
Thus, is fixed by the three parameters that generally vary with both position and time. For comparison, Fig. 5 shows the true (time averaged) PDF from the Monte Carlo simulation of Fig. 2a. To obtain the specific humidity field after condensation, we substitute (36) into the condensation formula given in (32), or equivalently (35). Depending on the proportion of supersaturated parcels in the distribution, in other words, the location of relative to , we have three cases:
We discuss how to determine in the next sections with further technical details concerning some exceptional cases given in appendix A.
a. Amplitude of the dry spike β
Because the dry parcels with simply move around the domain without undergoing condensation, it is particularly easy to calculate the amplitude of the dry spike β. Recalling from (24) that , we substitute into (21). Noting that , we find β satisfies
Because of rapid condensation, parcels at the top boundary always have . At the bottom boundary where the source is located, there is zero probability that . Hence, the boundary conditions at the top and bottom are
b. Center of , a
The value of is determined by the intermediate field obtained at the end of the advection–diffusion stage in the split-step algorithm described around (12). As mentioned before, is interpreted as the mean of the distribution . Substitution of (36) into (33) yields
Note that no extra prognostic equation is introduced here.
c. Width of , σ
The width of describes the subgrid-scale fluctuation of the specific humidity about its mean value before the action of condensation in each time step. To determine , we use the second moment,
from the stochastic model. Consider the advection–diffusion of the ensemble of parcels without condensation from time to and assume the initial condition is known. During this time, μ evolves to an intermediate value according to
which follows from (21). The boundary conditions are
and vanishing normal derivative at all other boundaries. Knowing , we set the value of σ in by requiring
After rapid condensation, the conditional PDF of specific humidity of the imagined ensemble becomes given by (34) and depicted in the bottom panel of Fig. 4. Therefore, the initial condition for the next iteration is
Let us now summarize the full procedure. Given , , and at time (with spatial arguments momentarily suppressed for clarity), we time step forward the three advection–diffusion equations: (12) for the moisture q itself, (39) for the amplitude of the dry spike β, and (43) for the second moment μ of the “true” distribution . This gives , , and , which in turn allows us to calculate a and σ from (41) and (46), respectively, and hence fully specify . Finally, the action of parameterized condensation depicted in Fig. 4 gives in (38) and in (47). Note that the full Fokker–Planck equation is not solved (nor could it be). Rather, there are only as many evolution equations as there are parameters in the assumed PDF.
We run the parameterized system until it reaches the steady state. We first examine the spatial structure of the PDF parameters in order to gain further insights into the parameterization process. Figure 6a plots the steady-state dry spike amplitude β. As expected from the boundary condition and the circulating flow pattern, along the top and east edges, while along the bottom and west boundaries. A more surprising feature is that for much of the area away from the boundaries. This means that in the central area, roughly half of the parcels in the imagined ensemble have the minimum specific humidity . This evinces the importance of subsidence of dry parcels by the random velocity. Figure 6b shows how close the center of is to the local saturation limit at different positions . We find that for all with the first inequality follows directly from (41). Figure 6c plots the measure of supersaturation of the ensemble. Recalling the schematic in Fig. 4, we see that inside the red supersaturated region where , some of the imagined parcels are about 1%–3% over the local saturation limit. This is the region where the condensation parameterization is in action.
We now assess the effectiveness of the condensation parameterization. Figure 2d shows the relative humidity field of the parameterized system at a late time. In contrast to from the unparameterized model with the same simulation parameters shown in Fig. 2c, does not have large areas of complete saturation and approximately resembles the bin-averaged field from the Lagrangian model in Fig. 2b. The most visibly noticeable discrepancy appears inside the boundary layer near x = 0. Figure 2f plots and shows that the biggest difference is located between such a boundary layer and the central dry region. Comparing Fig. 2f to Fig. 2e and noting the difference in the color scales, we can see quantitatively the improvement due to the parameterization. Figure 7 plots r, , and for the case of small eddy diffusivity and shows the condensation parameterization is similarly effective.
For further comparison, we plot the variation of the relative humidity along y at a fixed for κ =10−1 and 10−2 in Figs. 8a and 8b, respectively. For both values of κ, for all y. Near the top and bottom boundaries, and virtually have the same values, whereas elsewhere. We also examine the total moisture content in the system by calculating the mean specific humidity. For the Lagrangian formulation, we have
where N is the total number of parcels and indicates averaging over many snapshots in the statistically steady state. In the Eulerian formulation, using the steady solution at some large time , we compute
Figure 8c plots the mean specific humidity versus κ for the three models studied here. The moisture content of the unparameterized Eulerian system is the highest, and it also increases the fastest with κ. With the condensation parameterization implemented, the mean specific humidity in the Eulerian model and the rate at which it increases with κ are both reduced to nearly the same as those in the Lagrangian model.
An important quantity in atmospheric moisture transport is the vertical moisture flux F. As the advecting velocity in our present system is steady, we focus on the steady-state flux. Hence, for the Eulerian formulation, we compute
where is some large time in the simulations. In the Lagrangian formulation, we estimate by monitoring over a long period of time in the statistically steady state the specific humidity Q of those parcels crossing a given altitude y. We relegate the detailed implementation of this diagnostic to appendix B. Figure 9a plots the horizontal profile of the vertical moisture flux across for different models at . Figure 9b shows the same for the case of . Generally, there is a large positive flux associated with the rising arm of the overturning cell for and a small negative flux in the descending arm for . For both values of κ, the unparameterized system has the largest flux in magnitude because of its high moisture content. When the condensation is parameterized in the Eulerian model, the magnitude of the flux is reduced, and the profile becomes close to that of the Lagrangian model. Figure 9c plots the total vertical moisture flux across ,
for different κ. The total flux generally increases with κ. Not surprisingly, the unparameterized Eulerian model produces the largest at all κ. With the condensation parameterized, in the Eulerian model is reduced by about 50%. The smallest is observed in the Lagrangian model.
5. An unsteady channel flow
For our second example, we apply the condensation parameterization to a time-dependent flow. We introduce a configuration that roughly mimics the transport of moisture by baroclinic eddies along moist isentropic surfaces in midlatitudes (Vallis 2017). Consider a channel of width π in the y direction and periodic in the x direction. The streamfunction of the unsteady velocity in the channel is taken to be
The values of the parameters are , , , , , , and . We choose to ensure all streamlines are open and wrap around the periodic x direction. Figure 10c shows several streamlines of (52) at one instance of time. The waviness of the streamlines, controlled by , varies with time as the whole pattern propagates eastward. We use the saturation profile given in (5) with y interpreted as the meridional direction. Hence, we take and , which gives and . We once again assume an evaporation source that saturates air parcels is located along . At , we have . The domain is initially saturated. Advected by the time-periodic velocity, the moisture field eventually reaches a time-periodic state that varies at the same frequency as the velocity.
We first consider the Lagrangian formulation of the problem. As in the previous example, we perform Monte Carlo simulation of the stochastic system (15) and calculate the bin-averaged relative humidity field from the data. Figure 10c shows for at a late time after the transient, specifically, . We see that the areas along the top and bottom edges are close to saturation. Large regions of low relative humidity are formed in the middle of the channel. These dry patches are separated by tongues of humid air erupting periodically from the top and bottom boundary layers as the general large-scale pattern propagates eastward. Interestingly, the jets of humid air emerging from the bottom boundary are filamentous, creating sharp gradients in humidity. Similar features have been reported in more complex models of moisture decay on isentropic surfaces (Yang and Pierrehumbert 1994).
Turning to the Eulerian formulation, we recall that the system is now governed by the PDE (7). Figure 10a shows the relative humidity field obtained from a solution of (7) with unparameterized rapid condensation and Fig. 10b shows for the case when condensation is parameterized. We once again use the ansatz (36) in our parameterization with the three parameters determined by the same procedure described in the previous section. In both figures, and , that is, the same diffusivity and time instance as in Fig. 10c. All three relative humidity fields in Fig. 10 display the same general structure of high and low values. However, large areas of complete saturation can be seen in the unparameterized Eulerian model. Furthermore, the minimum relative humidity inside the dry patches is about 20%–30% higher than those in the other two models. Figure 11 shows the meridional relative humidity profile obtained by averaging over the zonal direction x and time t after the initial transient. The difference in the magnitude of the relative humidity minimum in the profile from the three models is obvious. The high moisture content in the unparameterized system is also evident in Fig. 12a, which plots the time evolution of the mean specific humidity. Interestingly, Fig. 12a also shows that the mean specific humidity oscillates with a much larger amplitude in the unparameterized system. Figure 12b shows how the specific humidity changes with time at one particular location over a single period of variation. We see that the evolution from the parameterized Eulerian model approximately follows the one from the Lagrangian model. This demonstrates that these two systems are not only close to each other in the average sense but actually have similar spatiotemporal behavior.
6. Parameterization in atmospheric models
Probabilistic (or statistical) schemes are often used in atmospheric general circulation models (GCMs), and sometimes cloud-resolving models, to parameterize subgrid-scale moisture variability. As discussed in the introduction, these schemes often employ turbulence closures to obtain the moments required to fix the assumed PDF. In view of our results that an Eulerian model with probabilistic condensation, namely, (30), can successfully mimic a Lagrangian model, we suggest a strategy that makes use of a stochastic Lagrangian model instead of turbulence closures. Of course, at a fundamental level, the two methodologies are not so different, for there is a close relationship between stochastic Lagrangian models and turbulence models, in particular, second-moment closures (Pope 1994b). However, our method avoids the “intermediate” step of constructing a closure. Stochastic Lagrangian models are also often used as models of turbulent diffusion (Rodean 1996) and the dispersion of passive, nonreactive scalars in the atmosphere (Wilson and Sawford 1996). It may also be noted (e.g., Pope 1994a) that the Lagrangian framework is especially fit for modeling reactive flows, and condensation can be considered mathematically as a form of reaction.
To construct a parameterization that might be used in an atmospheric GCM, one would first construct a stochastic Lagrangian model of water vapor transport, such as (15), for the atmospheric flow under consideration. Imagine an ensemble of moist parcels advected by the flow. Each parcel carries a set of thermodynamical variables (e.g., specific humidity and potential temperature) that evolves because of moist processes such as condensation and evaporation. The parcel moves with velocity V + V′. The large-scale velocity V is interpolated from the Eulerian velocity field u provided by the atmospheric model. To complete the model, we assume the salient properties of the small-scale parcel velocity V′ can be captured by a suitably chosen random process—so that our stochastic system is a good representation of the moisture dynamics. This is a nontrivial issue, and the details will depend on the particular application. For example, in some cases, the Markovian (i.e., memoryless), Brownian process is a sufficiently good model, while in other cases, it may be necessary to consider time-correlated or non-Markovian random processes.
If computational cost were not a constraint, we could perform Monte Carlo simulation (as in Fig. 2b) or solve for the governing PDF P using the Fokker–Planck equation of the stochastic system. The mean humidity field can then be computed. However, the large number of parcels required to obtain good statistics and the high dimension of the Fokker–Planck equation render these propositions impractical (and, in any case, were computational costs not a consideration, one could perform extremely high-resolution Eulerian simulations without the need to parameterize subgrid-scale motion). Instead, we apply the effective condensation , given in (30), to the evolution PDE for the Eulerian humidity field q in the atmospheric model; for example,
At each time step, the parameters of the assumed PDF embedded in are determined, as described in section 4, by matching a certain number of moments of to those of P. Investigation in sections 4 and 5 suggests that (54) will produce similar results to the stochastic Lagrangian model. Therefore, we can use (54) in place of the Lagrangian model to parameterize the actual moisture transport.
In the procedure presented here, the stochastic Lagrangian model forms the foundation of an integrated parameterization scheme. It provides the theoretical basis for the effective condensation, that is (30) (as discussed in section 3), as well as fixing the parameters in . It is in the explicit use of an underlying Lagrangian model that our method differs from other parameterization schemes in which moisture variability is not incorporated so directly. A possible advantage of the approach is the flexibility to incorporate different Lagrangian dynamics into the parameterization through the stochastic model (Wilson and Sawford 1996; Sawford 2001). Results from atmospheric tracer experiments (Stohl 1998) or novel theoretical transport models such as anomalous fractional diffusion (Goulart et al. 2017) may also be adopted into the scheme.
The potential disadvantage of the scheme is that it requires extra prognostic equations in addition to the one for the humidity q. Generally, the total number of equations equals the number of undetermined parameters in the assumed PDF; thus, in our examples where has three parameters, two additional equations are introduced, namely, (39) and (43). However, these prognostic equations are solved at the same resolution as the other variables in the atmospheric model. Modern GCMs often have a very large number of prognostic equations, especially if the model has an aerosol scheme, so that the additional expense of our scheme would be relatively small.
7. Summary and conclusions
The representation of subgrid-scale condensation of moisture in climate or weather models is a matter of both theoretical interest and considerable practical concern. Using the simple advection–condensation model, (7), we have shown that, without any condensation parameterization, a coarse-grained PDE model tends to retain excessive moisture and develop large regions of high humidity. Fundamentally, this is because the nonlinear condensation process and the coarse-graining operation do not commute, and local fluctuations are therefore lost when moisture is represented by a coarse-grained field, as illustrated in Fig. 3. On the other hand, the comparison in Fig. 2 shows that a Lagrangian formulation, where air parcels tagged with a humidity variable are tracked, is able to account for small-scale fluctuations, as found in nature.
It is, however, possible for an Eulerian model to produce results similar to the corresponding Lagrangian model if subgrid-scale moisture variability is properly introduced. Section 3 presents a way of achieving this using a probabilistic condensation parameterization given in (30). This mimics the Lagrangian representation of condensation in (28) by using an assumed PDF of humidity, with the parameters of the PDF being given through the use of the Fokker–Planck equation that governs the stochastic Lagrangian model. In both the simple single-celled circulation patterns shown in Fig. 7 and in the unsteady channel flow shown in Fig. 10, we see that this methodology reduces the excessive saturation in the Eulerian model, allowing it to produce moisture distribution close to that of the original Lagrangian model, obtained by a Monte Carlo simulation that explicitly follows the moist parcels. That the Eulerian model with a probabilistic parameterization can mimic the explicit Lagrangian model is a quite stringent test for the method.
The use of such a parameterization of condensation in a GCM trying to model real atmospheric flows would be rather more complex than our examples but would follow the same methodology. That is, presuming that trajectories in the atmosphere can be modeled by a stochastic Lagrangian model, a coarse-resolution Eulerian model that incorporates subgrid-scale Lagrangian information via the probabilistic condensation, (30), could be used to parameterize water vapor transport efficiently, as discussed in section 6. The first step is to construct a stochastic model analogous to (15) but with a more complete thermodynamics and, potentially, non-Markovian dynamics. The second step is to choose an ansatz for the PDF of the thermodynamic variables, a PDF that is determined by a small number of parameters, and to use that ansatz in the Fokker–Planck equation to determine those parameters. In the examples we computed, we chose the PDF to be a dry spike plus a continuous (top hat) component with a finite width, but other choices, with more free parameters, are possible. Since explicitly computing the Lagrangian model with a Monte Carlo simulation will not generally be possible in such cases, the efficacy of the choices will ultimately be determined by comparison with observation.
Testing this method in a range of models of varying complexity is the next step, starting from fairly idealized settings such as nonprecipitating moist Rayleigh–Bénard convection (Pauluis and Schumacher 2013) or a minimal precipitating convection model (Hernandez-Duenas et al. 2013). Note too that the general idea behind the method is not limited to the condensation “reaction”—its applicability to the parameterization of mix-down time in atmospheric chemical transport (Thuburn and Tan 1997) could also be explored.
YKT is indebted to Steven Böing for many helpful discussions. This work is supported by a Feasibility Grant from the EPSRC network Research on Changes of Variability and Environmental Risk (ReCoVER). YKT also thanks the EPSRC-funded Past Earth Network (Grant EP/M008363/1) for hosting a writing retreat during which part of this paper is written. GKV also acknowledges support from the Royal Society (Wolfson Foundation), the Leverhulme Trust, and NERC.
Exceptional Cases in the Determination of
As discussed in sections 4 and 5, we take the Lagrangian model, (15), as a good parameterization of the moist dynamics in the two systems considered there. We then determine the parameters of the assumed-shaped PDF in the effective condensation, (30), using information from the Lagrangian model. Although it rarely occurs in practice, two issues could, in principle, arise because generally does not satisfy the Fokker–Planck equation of the Lagrangian model. First, we would have in (46) if
If this occurs, we set . Second, when , and hence a, gets close to or , it is possible for a portion of (the continuous top-hat component of ) to lie outside the range . When this happens, we reduce σ so that either or . In the highly unlikely case that , we set and adjust β to make .
Estimation of Vertical Moisture Flux in Monte Carlo Simulations
Consider a Monte Carlo simulation using N parcels in a domain. Following Tsang and Vanneste (2017), we estimate the vertical moisture flux as follows. Assume between time t and , there are parcels crossing a given height y in either direction and whose x-position lies in . Let be the sign of ; then
The statistically steady is then obtained by averaging over t.