## 1. Introduction

The distribution of precipitation in complex terrain has usually been estimated using one of two methods: the interpolation of rain gauge station data and the so-called upslope method. In the former method, sparse station data are interpolated using inverse distance weighting or spline fitting, sometimes including an altitude or aspect correction (Daly et al. 1994; Hutchinson 1998). No influence of horizontal scale is included. In the latter method, the terrain slope and wind speed are used to estimate the condensation rate above the terrain (Collier 1975; Rhea 1978; Smith 1979; Neiman et al. 2002). Precipitation is computed from the assumption that condensed water falls immediately to the ground. Again, physical scale is not taken into account.

Several quasi-analytic models have included scale-dependent processes. Hobbs et al. (1973) and Bader and Roach (1977) have described raindrop/snow fall speeds and washout of cloud water. More recently, Alpert and Shafir (1989), Sinclair (1994), and Smith (2003a, hereafter S03) have included the advection of hydrometeors in their upslope models. The drift distances in these models, from 5 to 25 km, improve the agreement with station and radar data.

Simple upslope-advection models such as S03 have two limitations. First, they assume that the terrain-induced vertical air velocity penetrates upward through the moist layer without change. In fact, it is known from mountain wave theory that such motions may either decay or oscillate with altitude, depending on the horizontal scale and aspect ratio of the terrain and the static stability and wind speed in the environment. In either case, the upslope estimate may far exceed the actual condensation rate. Fraser et al. (1973), Durran and Klemp (1982b), and Kuligowski and Barros (1999) discuss the role of moist wave dynamics in a vertical plane.

A second problem relates to leeside descent. In upslope models, it is usually assumed that only upslope regions influence precipitation. With instantaneous conversion and fallout, the windward slope receives rain while the flat and downslope regions are dry. When time delays are included, the upslope condensate is distributed downstream, without regard for local terrain. This approach neglects the evaporation of cloud water and hydrometeors caused by descending air. The total precipitation is overestimated by this assumption. Models such as Alpert and Shafir's introduce an adjustable coefficient to reduce the prediction into a reasonable range.

The errors induced by these two problems are strongly scale dependent. For smooth hills with scales of 100 km, the upslope estimates are usually quite reasonable. On the other hand, when the terrain is rising and falling with scales of 20 km or less, both assumptions fail and the model can overestimate the total precipitation by a factor of 5 or greater. In the Italian Alps, Smith et al. (2003) showed that the upslope rainfall estimate exceeded the incoming moisture flux, an obvious violation of water conservation. The scale dependence of these errors means that the investigator must smooth the terrain by an “appropriate” amount before applying the model. The choice of smoothing scale introduces an arbitrary factor into the analysis.

Full numerical mesoscale models are now being used frequently to predict orographic precipitation. These models are replete with physical scales, in both the governing equations and the parameterizations. Numerical models are difficult to fully diagnose, however, because of their complexity and slow execution speed. Usually, less than 100 runs can be carried out to investigate their properties. Bruintjes et al. (1994), Miglietta and Buzzi (2001), Jiang (2003), and Smith et al. (2003) have been partially successful in diagnosing model behavior, but their range of parameter values was sharply restricted.

The goal of this paper is to develop a model of orographic precipitation with the following characteristics:

analytically tractable so that its properties can be easily understood;

applicable to actual complex terrain and arbitrary wind direction so that it can be tested against real data;

reduces to the classical upslope model so that it can be compared with earlier work;

includes the basic physical elements: airflow dynamics, condensed water conversion, advection and fallout, and downslope evaporation, leading to a theory of precipitation efficiency.

The first criterion suggests that our model should be mathematically linear. The challenge is to retain linearity while satisfying the other three criteria. Problems of cloud formation do not lend themselves to linear formulation because of the threshold nature of water vapor saturation (e.g., Barcilon et al. 1980; Durran and Klemp 1982b; Barcilon and Fitzjarrald 1985). Linear models must also ignore other nonlinearities such as moist airflow blocking (Jiang 2003) and cloud physics bifurcations (Jiang and Smith 2003). Nevertheless, we expect that our model will have a substantial domain of applicability.

The development of the model is described in four parts. In section 2, the equations for advection, conversion, and fallout of condensed water are introduced. In section 3, the role of airflow dynamics is discussed, following classical 3D mountain wave theory. In section 4, a treatment of evaporation and precipitation efficiency is formulated. In section 5, the combined effects of all these processes are discussed. The last two sections of the paper consider applications. In section 6, the influence of topographic scale is summarized. In section 7, we display the properties of the linear model with real terrain. The appendix describes our thermodynamic formulation.

## 2. Governing equations for advection, conversion, and fallout

In this section, we formulate a system of equations that describes the advection of condensed water by the mean wind.

### a. Formulation

*S*(

*x,*

*y*)] arising from forced ascent. While conventional “upslope” models assume instantaneous fallout of condensed water, we allow delays and downstream advection. We begin, following Smith (S03), by postulating a pair of steady-state advection equations describing the vertically integrated cloud water density [

*q*

_{c}(

*x,*

*y*)] and hydrometeor density [

*q*

_{s}(

*x,*

*y*)].

*τ*

_{c}is the time constant for conversion from cloud water to hydrometeors (i.e., rain or snow) and

*τ*

_{f}is the time constant for hydrometeor fallout. The physics of conversion and fallout in orographic clouds was discussed by Hobbs et al. (1973), Rauber (1992), and others. As our model is vertically integrated, we use average values of the time constants, representative of the whole column. Estimates of the time constants range from 200 to 2000 s (Jiang and Smith 2003; Smith et al. 2003). The advecting wind vector

**U**=

*U*

**i**+

*V*

**j**is assumed to be independent of time and space.

*S*(

*x,*

*y*) is the sum of a background rate of cloud water generation and local variations created by terrain-forced uplift. A simple proposal for

*S*(

*x,*

*y*) is the upslope model

*S*

*x,*

*y*

*S*

_{∞}

*C*

_{w}

**U**

*h*

*x,*

*y*

*C*

_{w}=

*ρ*

_{Sref}(Γ

_{m}/

*γ*)

*h*(

*x,*

*y*) is the terrain, and

*S*

_{∞}is the background nonorographic large-scale vertically integrated condensation rate. Equation (3) assumes saturated conditions and a terrain-forced vertical air velocity

*w*(

*x,*

*y*) =

**U**· ∇

*h*(

*x,*

*y*) that is independent of altitude. The final term in (1) is the conversion of cloud water to hydrometeors. It appears as a source term in (2). The final term in (2) is the loss of hydrometeors associated with precipitation

*P*

*x,*

*y*

*q*

_{s}

*x,*

*y*

*τ*

_{f}

To clarify the properties of (1)–(4), we consider three trivial solutions. First, in the absence of terrain *h*(*x,* *y*), the uniform source *S*_{∞} will create constant values of cloud water column density (*q*_{c} = *S*_{∞}*τ*_{c}), hydrometeor column density (*q*_{s} = *S*_{∞}*τ*_{f}), and precipitation (*P*_{∞} = *S*_{∞}). Because our equations are linear, the background precipitation can be added to any orographic component.

Second, in the special case of *S*_{∞} = 0 and *τ*_{c} = ∞ (i.e., no conversion), (1) and (3) can be integrated to give *q*_{c}(*x,* *y*) = *C*_{w}*h*(*x,* *y*) and *q*_{s} = 0, indicating the existence of cloud water over high terrain, but no precipitation. We call this the “lenticular” solution as steady clouds exist over each mountain peak with equal condensation and evaporation at their leading and trailing edges.

Third, if *S*_{∞} = 0 and *τ*_{f} = ∞ (i.e., no fallout), forward and reverse conversion can occur but the total condensed water behaves simply; that is, *q*_{c} + *q*_{s} = *C*_{w}*h*(*x,* *y*). Rainwater can be created but then it is stored in the parcel until, upon descent, negative *S* drives *q*_{c} negative, which in turn evaporates *q*_{s}. See section 4 for further discussion of evaporation.

*k,*

*l*) are the components of the horizontal wavenumber vector and

*σ*=

*Uk*+

*Vl*is the intrinsic frequency (see also Smith 2003b for examples of solving the advection equation with Fourier methods). The transform (5) can be inverted to obtain the precipitation distribution

*P*(

*x,*

*y*) using

*S*(

*x,*

*y*) is positive definite, so is

*P*(

*x,*

*y*). This property is proven in the next subsection. In section 4, we extend the definition of

*q*

_{c}to include negative values. A second property of (5), (6) is that, if

*τ*

_{c}=

*τ*

_{f}= 0 (i.e., fast conversion and fallout), then

*P̂*(

*k,*

*l*) =

*Ŝ*(

*k,*

*l*) and thus

*P*(

*x,*

*y*) =

*S*(

*x,*

*y*).

### b. Response to a point source

*S*(

*x,*

*y*) =

*Bδ*(

*x*)

*δ*(

*y*) in axis-parallel flow; for example,

*U*> 0,

*V*= 0. In applied mathematics, such a solution is called the Green's function. From (5), (6), using contour integration and the residue theorem,

*P*

*x,*

*y*

*Bδ*

*y*

*e*

^{−x/Uτ}

*Uτ*

*P*

*x,*

*y*

*Bδ*

*y*

*xe*

^{−x/Uτ}

*Uτ*

^{2}

*Uτ*downstream of the source. While (3), (5) represents a convenient model for distributing precipitation downstream, it has two serious disadvantages. First, if only positive values of

*S*are retained, every upslope source region causes a swath of precipitation with perfect efficiency. Smith et al. (2003) showed that, for complex terrain, this approach can seriously overestimate the precipitation. The same water molecules are lifted repeatedly and, according to (3), (5), every lifting event produces precipitation. With such a model, reasonable total precipitation estimates can be achieved only by smoothing the terrain to 30 km or more so that multiple lifting is minimized.

The second problem with (3), (5) is that it takes no account of airflow dynamics: including phenomena such as lateral influence of terrain, upwind tilt of the forced ascent, and the decay of forced ascent aloft.

## 3. Airflow dynamics and the condensed water source term

In this section we examine the effect of airflow dynamics on the vertically integrated condensation rate.

### a. Formulation

*C*

_{w}is the thermodynamics uplift sensitivity factor (see appendix) and

*H*

_{w}is the depth of the moist layer. To describe the dynamics of forced ascent, we use well-known results from linear Boussinesq mountain wave theory (e.g., Queney 1947; Sawyer 1962; Smith 1979, 2002; Wurtele et al. 1996). For each Fourier component (

*k,*

*l*), the vertical velocity takes a complex exponential form

*ŵ*

*k,*

*l,*

*z*

*ŵ*

*k,*

*l,*

*e*

^{imz}

*N*

^{2}

_{m}

*σ*

^{2}), the airflow is irrotational and

*m*=

*i*(

*k*

^{2}+

*l*

^{2})

^{1/2}. In the hydrostatic limit (

*N*

^{2}

_{m}

*σ*

^{2}), (12) reduces to

*l*= 0, and

*σ*=

*Uk*so that (13) reduces to

*ŵ*

_{0}(

*k,*

*l*) =

*ŵ*(

*k,*

*l,*0) =

*iσĥ*(

*k,*

*l*) from mountain wave theory, (15) yields a transfer function relating the transform of the terrain and the condensed water source function

In the denominator of (16), the square bracket contains the effect of vertical velocity variations up through the moist layer. If *σ*^{2} > *N*^{2}_{m}*m* is positive imaginary and the forced ascent decays with height. The denominator of (16) is real, positive, and greater than unity, reducing the amount of condensation. If the static stability (*N*) is large or the intrinsic frequency *σ* is small, *m* will be real. In this case, the vertical structure oscillates with height and partial cancellation occurs. The denominator of (16) is complex with a magnitude exceeding unity. Again, the total condensation is reduced. In both these cases, the depth of the moist layer (*H*_{w}) plays a key role. Only if |*mH*_{w}| ≪ 1, will the vertical motion penetrate through the moist layer without shift or decay.

The effect of static stability enters several ways in (16). For example, an increased static stability (i.e., a smaller magnitude of *γ*) increases the available water vapor (*C*_{w} and *H*_{w}), but decreases the depth of the lifting [i.e., an increased *m* from (A13) and (12)].

In the following subsections we examine five different examples of how airflow dynamics influences cloud water condensation.

### b. Sinusoidal terrain

*λ*), specified by

*h*

*x*

*Ae*

^{ikx}

*k*= 2

*π*/

*λ.*If the coefficient

*A*is real, (17) is

*h*(

*x*) =

*A*cos(

*kx*). Equation (17) can be used directly in (16), using generalized functions, to give

*m.*In the irrotational case (12) gives

*m*=

*i*|

*k*|, so the vertically integrated condensation is, from (18)

*k,*the condensation rate does not increase significantly as the increasing terrain slope is cancelled by the decreasing penetration.

*N*

_{m}=

*Uk,*the vertical wavenumber is

*m*= 0 [from (12)] and complete penetration occurs; that is, from (18)

*S*

*x*

*C*

_{w}

*kU*

*Aie*

^{ikx}

*A*is real,

*S*(

*x*) = −

*C*

_{w}

*kUA*sin(

*kx*).

*Ĥ*=

*N*

_{m}

*H*

_{w}/

*U*is the nondimensional moist layer depth. If the coefficient

*A*is real, (21) can be written

*Ĥ*can be interpreted as the ratio of the moist layer depth to the penetration depth of the forced ascent

*d*=

*U*/

*N*

_{m}. A typical value for

*Ĥ*

^{−1}

^{−1}

_{dyn}= 2

^{−1/2}= 71%. If

*Ĥ*≫ 1, the forced ascent is only felt by a small lower part of the moist layer, and the net condensation is very small (PE

_{dyn}≪ 1).

### c. Rise to a plateau

*h*(

*x*) = (

*A*/

*π*) tan

^{−1}(

*x*/

*a*), with total rise Δ

*h*=

*A.*The vertical velocity near the ground is

*Ĥ*is independent of

*k*in (26), the denominators can be brought outside the integrals. Evaluating (26) gives

*S*(

*x*). According to (27), see Fig. 2a, one effect of

*Ĥ*is to shift the condensation pattern upstream. The negative source region in (27), located at

*x*>

*a*/

*Ĥ,*is associated with a region of descending air in the mountain wave. Equation (27) satisfies (3) when

*Ĥ*= 0. The total source is

*Ĥ*= 1, (29) gives PE

_{dyn}= 50%. It is interesting that, even with permanent lifting to a plateau, the condensation is significantly reduced due to the lack of penetration of the vertical motion field aloft.

### d. Witch-of-Agnesi ridge

*πAa.*Its Fourier transform is

*S*(

*x*). Equation (33) satisfies (3) when

*Ĥ*= 0. Compared to the raw upslope function (i.e.,

*Ĥ*= 0), this formula shifts the condensation upstream and reduces its magnitude (Fig. 2b). It crosses zero and becomes negative, slightly upstream of the ridge peak. It also generates a small positive source region downwind, associated with ascent in a mountain wave aloft. This source function (33) has a zero net integral, but its upslope part can be isolated as

*S*

_{upslope}reduces to the expected raw upslope value. The dynamical precipitation efficiency is again given by (29).

### e. Isolated Gaussian hill

*h*

*x,*

*y*

*Ae*

^{−[(x/ax)2+(y/ay)2]}

*πAa*

_{x}

*a*

_{y}. The aspect ratio for southerly winds is the ratio of the east–west width (

*a*

_{x}) to the north–south width (

*a*

_{y}); that is,

*R*=

*a*

_{x}/

*a*

_{y}. In the hydrostatic case, the source distribution generated by (13), (16) is independent of scale, but dependent on

*R*and

*Ĥ.*To describe the influence of aspect ratio, we examine its influence on both the precipitation along the upwind centerline and over the full upslope region.

*y*axis) we define the upslope centerline total source as

*C*

_{w}

*VA.*Actual values of

*S*

_{CL}, determined from (36) using fast Fourier transform (FFT) evaluation of (13), (16), are shown in Table 1, normalized by the reference value (i.e.,

*S*

_{CL}/

*C*

_{w}

*VA*). From Table 1, we see that the circular hill (

*R*= 1) and the wind-aligned ridge (

*R*< 1) produce much less condensation on the centerline than a crosswind ridge (

*R*> 1) of the same height. This effect can be explained by recalling how lateral structure in a mountain wave field influences vertical scale. For a wind-aligned ridge, the dominant terrain wavenumber vectors have a large ratio of the lateral to the windward components. According to (12) or (13), the role of the lateral component is to increase

*m,*decreasing the vertical penetration of the forced ascent. Thus, the depth of the ascending air is only a fraction of the nominal depth (i.e.,

*U*/

*N*

_{m}). It follows that a crosswind ridge produces ascent that penetrates well aloft, while a wind-aligned ridge lifts only a shallow layer of air. If the moist layer is very thin, this rapid decay of vertical velocity (

*w*) has little influence. For a typical value of

*Ĥ*= 1, however, the reduction by 3D effect is quite significant. For a circular hill (

*R*= 1), the condensation over the upslope centerline is only 41% of the raw upslope value.

*π*

*C*

_{w}

*UaA*for a circular Gaussian hill (

*a*=

*a*

_{x}=

*a*

_{y}) with

*Ĥ*= 0. The normalized values of

*S*

_{upslope}are given in Table 2 as a function of nondimensional moist layer depth. With

*Ĥ*= 1, the total upslope condensation has dropped to 63% of the raw “no-dynamics” value. Tables 1 and 2 show that both centerline and total upstream condensation are reduced by dynamics. The total upslope condensation is reduced less due to lateral spreading of the forced ascent.

The effect of lateral spreading is also found when the wind blows through a mountain pass. While the local gap-forced ascent may be zero in such a case, the hills on either side produce ascent above the gap. If the gap is narrow and *Ĥ* is one or larger, the gap will experience nearly as much precipitation as the adjacent hill slopes.

### f. Influence of the Coriolis force

Until now we have neglected the influence of the Coriolis force. This neglect is consistent with our focus on relatively small horizontal scales from 100 m to 100 km. Other assumptions of our model break down at scales of 200 km and larger, where the Coriolis force begins to dominate the dynamics. For example, the assumption of uniform background flow or the neglect of surface evaporation would fail at synoptic scales. Furthermore, the focus of our study, the dynamic and cloud-delay controls on precipitation, are not very important at the larger scales. The vertical penetration of broad-scale uplift is substantial and the cloud processes are relatively quick for very broad scale terrain.

*d*=

*fa*/

*N*

_{m}. If

*f*= 10

^{−4}s

^{−1},

*a*= 500 km, and

*N*

_{m}= 0.005 s

^{−1}, the penetration depth is 10 km, well in excess of the typical moist layer depth.

## 4. Downslope evaporation

In this section, we describe how the linear model determines precipitation efficiency by treating downslope evaporation.

### a. Formulation

We now consider the role of negative source regions (i.e., evaporation) associated with descent. Negative source regions were encountered in all the examples in section 3, even when the terrain rises monotonically to a plateau (24). Negative source regions in orographic flow could only be avoided if the background condensation rate (*S*_{∞}) were larger than the largest negative orographic perturbation.

*S*

_{trun}

*x,*

*y*

*S,*

A possible evaporation scenario is illustrated in Fig. 1. In Fig. 1a, the sum of the two uplifts exceeds the net uplift (i.e., *S*_{1} + *S*_{2} > *S*_{3}). In cases with repeated uplift, the sum of sources can even exceed the incoming flux (i.e., *S*_{1} + *S*_{2} > *F*). If the total condensate falls immediately to the ground or drifts without evaporation, the model precipitation can exceed the influx (i.e., *P*_{1} + *P*_{2} > *F*)—an unphysical result. In reality, some fraction of *S*_{1} must evaporate over the next downslope region to allow condensation on the next uplift. In Fig. 1b, smooth ascent guarantees *S*_{3} < *F* and *P*_{3} < *F* without evaporation.

*P*

_{trun}

*x,*

*y*

*P,*

*P*includes the orographic and the background precipitation (

*P*

_{∞}). From (1), (2), and (4) we see that negative values for

*P*can arise where descent (

*S*< 0) has first driven

*q*

_{c}negative. Then, through the conversion term,

*q*

_{s}goes negative also while

*q*

_{c}is driven back toward zero. This sequence represents the evaporation of cloud water followed by the evaporation of hydrometeors. In cases when

*q*

_{c}can change sign, it is better to think of this quantity as the “degree of supersaturation.” When positive, it represents the column density of condensed cloud water. When negative, it represents the degree of subsaturation, that is, the amount of rainwater that could be evaporated to return

*q*

_{c}to zero. In (1) and (2), we use the same time scale

*τ*

_{c}to represent the conversion of cloud water to hydrometeors and the evaporation of hydrometeors by subsaturated air.

In cases with weak *P*_{∞}, rapid conversion and fallout and persistent descent and drying our linear “saturated” formulation becomes problematic. When *q*_{c} < 0 and *q*_{s} ≤ 0, the conversion term in (1) continues to remoisten the air even though there are no hydrometeors to evaporate. The truncation (40) assigns *P* = 0 to all such scenarios.

The airflow formulation in section 3 is also problematic in persistent dry descending flow, as we used the moist stability in the wave dynamics (14). While Jiang showed that dry leeside descent had little influence on windward ascent on single hills, strongly dried air encountering a second hill will not be treated well by our model.

### b. Sinusoidal terrain

*β*

_{c}=

*Ukτ*

_{c}and

*β*

_{f}=

*Ukτ*

_{f}. This function is shifted downstream compared to the source function (20) and it has reduced amplitude. The truncation condition (40) keeps only the positive parts of the sine wave in (41). A cloud-physics precipitation efficiency can be defined as the ratio of the positive areas of the precipitation and source sine waves. From (20) and (41), if

*P*

_{∞}= 0,

_{cloud}

*β*

^{2}

_{c}

*β*

^{2}

_{f}

^{−1/2}

*τ*

_{c}=

*τ*

_{f}= 1000 s,

*U*= 10 m s

^{−1},

*λ*= 20 km, and the terrain wavenumber

*k*= 2

*π*/

*λ*= 2

*π*/20 000 so that

*β*

_{c}=

*β*

_{f}= 3. The cloud efficiency (42) is PE

_{cloud}= 0.1. This reduction in efficiency is due to the effect of downslope evaporation.

### c. Triangle ridge

*G*is given by (7). For simplicity, we consider a triangle ridge of height

*A*centered at

*x*= 0 with total width 2

*d.*The cross-sectional area of the ridge is

*Ad.*The reference source functions on the windward and leeward slopes are

*S*

_{ref}

*C*

_{w}

*Udh*

*dx*

*c,*

*c*=

*C*

_{w}

*UA*/

*d.*For simplicity, we also set one of the time constants equal to zero so that the appropriate Green's function is (8). Using (8), (43), the sum of all the upwind sources gives an expression for the precipitation on the windward slope

*P*(

*x*) = 0. According to (45), the precipitation spills over the peak, reaching to the cutoff position

*x*

_{c}

*Uτ*

*e*

^{−d/Uτ}

*d.*

*x*= −

*d*to 0 and (45) from

*x*= 0 to

*x*=

*x*

_{c}gives the area total precipitation (

*P*

*d*/

*Uτ*≫ 1 and (47) becomes

*d*/

*Uτ*≪ 1, (47) gives

_{cloud}

*Uτ*

*d*

## 5. Combined airflow dynamics, advection, time delays, and downslope evaporation

### a. Mathematical properties and controlling parameters

*ĥ*(

*k,*

*l*) and the precipitation field

*P̂*(

*k,*

*l*). As before, the coefficient

*C*

_{w}is given by (A9). The first factor in the denominator of (49) describes airflow dynamics. The second and third factors describe cloud delays and advection. The reduction in precipitation efficiency due to downslope evaporation is not contained explicitly in (49), but is present when (40) is applied after the inverse Fourier transform.

The combined theory (49) contains a number of dimensional parameters and variables (see Table 3) and four nondimensional control parameters (Table 4). As the model is linear, the mountain height (*A*) appears as a coefficient rather than a control parameter.

An interesting property of (49) is that the dynamics and cloud-delay factors in the denominator have a similar form. The appearance of *i* = *m*(*k,* *l*) and *σ*(*k,* *l*). Note that the two cloud time scales, *τ*_{c} and *τ*_{f}, are mathematically analogous. When *m* = *τ*_{c} = *τ*_{f} = 0, (49) reduces to the standard upslope model with no airflow dynamics and no condensed water advection.

In the remainder of this section, we will apply (49) to two idealized hill shapes to illustrate its properties: that is, a triangle ridge and a 3D Gaussian hill.

### b. Triangle ridge

The triangle ridge is a useful example, as the raw upslope condensation rate is constant over the windward slope. Thus, it is easy to see modification caused by airflow dynamics. Some insight into the pure cloud physics effect for a triangle ridge was given by (47). The combined influence of full dynamics and cloud time delays [(13) and (49)] is shown in Fig. 3, with parameters: *T*_{0} = 280 K, *γ* = −5.8°C km^{−1}, *U* = 15 m s^{−1} so Γ_{m} = −6.5°C km^{−1}, *N*_{m} = 0.005 s^{−1}, *ρ*_{Sref}^{−3}, *H*_{w} = 2500 m. Also *A* = 500 m and *τ*_{c} = *τ*_{f} = 1000 s.

In both parts of Fig. 3, the effect of airflow dynamics is to reduce the total condensation and shift the maximum upwind, close to the “slope break” of the triangle ridge. The source term turns negative slightly upstream of the hill crest. The effect of cloud delay reduces the precipitation further and shifts the precipitation peak downstream. For the narrower ridge, the precipitation maximum (2.96 mm h^{−1}) is close to the hilltop. For the wider ridge, the maximum (2.33 mm h^{−1}) is on the windward slope, about two-thirds of the way toward the hill crest. For both ridges, there is downstream condensation in a lee wave, but leeside descent and evaporation prevent precipitation. Especially for the narrow ridge, the ratio of total precipitation to total raw upslope condensation is quite small.

### c. Isolated Gaussian hill

The isolated circular Gaussian hill (35 with *a*_{x} = *a*_{y}) is useful for showing lateral spreading, downstream drift, and the mountain wave contribution. Recall that the pure role of airflow dynamics over the isolated Gaussian hill was discussed in section 3e. Combined effects are shown in Fig. 4.

In Fig. 4a, we see that the raw upslope precipitation is confined to the windward slopes. With airflow dynamics (Fig. 4b), the condensation is much more widespread in the upstream region. Strong negative condensation values are seen in a butterfly pattern downstream. By far the strongest negative values are just downstream of the peak due to mountain-wave-induced descent there. Farther downstream, a significant positive region of condensation is present associated with a wave cloud aloft.

After the cloud advection and truncation have been applied (Fig. 4c), the precipitation peak shifts to the hilltop. Some spillover is seen. A dry region is maintained on the lee slopes. A modest region of precipitation appears well downstream due to ascent in the mountain wave above the lee slope. Most of this wave-induced precipitation has been eliminated by the leeside drying. In the real world, there is little evidence for lee wave precipitation, so we have cause to question this small remaining precipitation predicted by the linear model (see Bruintjes et al. 1994).

## 6. Scale dependence

*U,*

*N*

_{m},

*H*

_{w},

*τ*

_{c},

*τ*

_{f}, as the mountain width (

*a*) varies from 100 m to 100 km. The ratios are

*S*

_{ref}is an area integral computed from (3), (40);

*S*

_{dyn}and

*P*

*S*

_{dyn}we set

*τ*

_{c}=

*τ*

_{f}= 0. The fixed parameters are

*U*= 15 m s

^{−1},

*N*

_{m}= 0.005 s

^{−1},

*H*

_{w}= 3000 m,

*τ*

_{c}= 1000 s,

*τ*

_{f}= 1000 s. Unlike Table 1, the

*S*

_{dyn}in Fig. 5 includes the lee wave contribution. Note that, unlike sinusoidal terrain [e.g., (50)], the effects of pure airflow dynamics and pure cloud delay are not multiplicative. The ratio PE

_{cloud}in (51) describes the influence of cloud delays in the presence of airflow dynamics.

In Fig. 5, the ratio PE_{dyn} rises quickly as width increases. For small-scale terrain, when *σ*^{2} > *N*^{2}_{m}*a* = 5 km, a significant peak occurs associated with a near match between the mountain width and the buoyancy scale; that is, *a* = *U*/*N*_{m} = 3 km. This peak arises in part from good penetration and in part from the positive source contribution in a train of untrapped dispersive lee waves. These lee wave source regions do not cause precipitation as the condensed water from each cloud is advected quickly into a region of evaporation. Essentially they are “lee-wave lenticular clouds”. This effect disappears and hydrostatic results obtain when the mountain width *a* exceeds 10 km for the chosen parameters, as orographic forcing at the resonant wavelength is much reduced.

We now discuss the curve for PE_{cloud} in Fig. 5. For small-scale terrain, the precipitation efficiency from cloud processes is small. Condensed water from the upslope region is advected quickly to the lee side and evaporates before it has a chance to convert to hydrometeors and fall out (section 4). As the terrain scale increases, the time spent by a parcel in the upslope region increases. Progressively more of the condensed water has a chance to precipitate. For wide terrain, the value of PE_{cloud} approaches unity as all the condensate precipitates.

The combined effect of airflow dynamics and cloud delay on PE is shown in Fig. 5; PE remains low until PE_{cloud} begins to grow. For wide hills, PE_{cloud} ≈ 1 and PE ≈ PE_{dyn}, so PE can never exceed the value PE_{dyn} set by airflow dynamics. That is, if water does not condense, it cannot precipitate.

Changing the time-delay values will influence the precipitation efficiency, as shown in Table 5. As only the time constants are changed in the table, the quantity PE_{dyn} does not change. In general, as the time delays are increased, the PE_{cloud} decreases. More condensed water is advected onto the lee slope where it evaporates. In Table 5, note that the entry for *τ*_{c} = *τ*_{f} = 1000 s agrees with Fig. 5.

_{m}/

*γ*) = 2,

*A*= 500 m, and

*H*

_{w}= 3000 m, then DR = 20%. For order of magnitude calculations, we can roughly combine hydrostatic dynamics (29) and box model advection (48) with (52) to obtain

The drying ratio is a convenient measure of the nonlinearity of the orographic precipitation system. A key assumption in linear theory is that the air remains near saturation everywhere. If the drying ratio is reasonably small (perhaps less than 0.3) this assumption may be useful. If DR exceeds 0.5 however, this assumption is less appropriate.

## 7. An application of the linear model

To illustrate the properties of the linear model, we present one example of a predicted precipitation pattern over real terrain. We select the Olympic Range in Washington State as it is compact, complex, and relatively well studied. It is one of the rainiest spots in North America, but with a definite rain shadow on the northeast side. Our intention in this section is not to test the model, but only to exhibit its behavior.

For the example, we consider a southwest wind with speed 15 m s^{−1} and a moist stability of 0.005 s^{−1}. The surface temperature and specific humidity are 280 K and 6.2 g kg^{−1}. The moist layer depth is 2.5 km. The cloud time delays are each 1000 s. The calculation was done on a large 1024 by 1024 grid with one kilometer resolution. The surrounding mountains were reduced with a Gaussian weighting function centered on the Olympics. Only a small portion of the computational domain is shown in Fig. 6. The calculation took about 8 seconds on a small workstation. The same computation could be done on a 256 × 256 grid in about 1 s.

In Fig. 6, the terrain has been smoothed with a 3000-m spectral filter for clarity of presentation, even though the terrain was only smoothed to 800 m for the FFT model run. The 6-h accumulated precipitation is shown in millimeters with a maximum value of about 26 mm just upwind of the highest peak, Mount Olympus (2428 m). Several features can be noted. Four tongues of high precipitation are associated with four southwestward directed ridges. Light precipitation is found well upstream of the mountains, even over the sea. There is some spillover, but mostly the northeast lee slopes are dry. The model predicts that the high peaks in the northeast part of the massif collect no precipitation.

The sensitivity of these results to the model parameters can be estimated from the previous sections. The comparison with the raw upslope model (*m* = *τ*_{c} = *τ*_{f} = 0) is the most striking (not shown). It gives strong spikes of precipitation directly over each southwestern facing slope including those in the northeast corner of the range. Peak precipitation exceeds those in Fig. 6 by an order of magnitude.

## 8. Conclusions

We have developed a linear model of orographic precipitation including airflow dynamics, cloud time scales and advection, and downslope evaporation. The model is easy to apply to complex terrain. Only four steps are required: Fourier transform the terrain *h*(*x,* *y*) to obtain *ĥ*(*k,* *l*), multiply by the transfer function (49), perform an inverse Fourier transform, and apply the positive cutoff (40). The input parameters are *P*_{∞}, *U,* *V,* *T*_{0}, *N*_{m}, *τ*_{c}, *τ*_{f}, and a measure of vertical structure, *N*_{m} or *γ.* There is no need to smooth the terrain. The model includes crude representations of physical processes that provide the proper weights to the different scales, thus providing a theory of precipitation efficiency.

Several strong assumptions were made in the model formulation, for example, linear steady wave dynamics, near saturation, constant wind and moist stability with height and location, constant time delays, equal hydrometer growth and decay times, etc. The model only treats the vertically integrated condensed water. These simplifications could limit the accuracy and practical application of the model. Also, the model is unsuitable for unstable atmospheres.

Both the airflow and condensed water formulations become problematic in situations with fast conversion and fallout and persistent descent. The drying ratio (DR) can be a useful indicator of this breakdown in the theory.

We examined the role of scale in orographic precipitation using the linear model. Four natural atmospheric length scales appear: buoyancy scale, moist-layer thickness, and two drift distances. Both the pattern of precipitation and the total amount of precipitation are controlled by the ratios of mountain width to these inherent atmospheric scales.

The location of maximum precipitation is determined by a competition between the upstream shift caused by dynamics and the downstream shift caused by cloud delays. As mountain width decreases, the location of the maximum precipitation shifts from the windward slope to the hilltop. For narrow mountains and high wind speeds, spillover is possible.

The amount of precipitation is determined by the vertical penetration of the forced ascent relative to the depth of the moist layer and by the speed at which cloud droplets can convert and fall. Small-scale hills produce uplift that decays rapidly aloft, condensing little water. The ascent from broader hills penetrates much more using gravity wave dynamics. This penetration is also limited, however, and the condensation rate can be well below the raw upslope value.

The time delays from cloud processes also reduce precipitation. Here there is competition between the speed of conversion and fallout, and the time it takes for air parcels to reach the descending air in the lee. The cloud water or hydrometeors may drift to the lee side of a hill and evaporate before they can convert and fall to the ground. The combined effect of airflow dynamics and cloud delays may reduce precipitation by as much as an order of magnitude below the raw upslope model prediction.

We applied the linear FFT model to the Olympics under a southwest wind. The model ran quickly, with reasonable results. Quantitative testing of the linear model is under way and will be presented separately.

## Acknowledgments

Useful comments were given by Steve Skubis, Mathew Fearon, Qingfang Jiang, Michael Kunz, Gerard Roe, and three helpful reviewers. This research was partially supported by the National Science Foundation, Division of Atmospheric Sciences (ATM-0112354).

## REFERENCES

Alpert, P., and H. Shafir, 1989: Meso-gamma-scale distribution of orographic precipitation: Numerical study and comparison with precipitation derived from radar measurements.

,*J. Appl. Meteor***28****,**1105–1117.Bader, M. J., and W. T. Roach, 1977: Orographic rainfall in warm sectors of depressions.

,*Quart. J. Roy. Meteor. Soc***103****,**269–280.Barcilon, A., and D. Fitzjarrald, 1985: A nonlinear steady model for moist hydrostatic mountain waves.

,*J. Atmos. Sci***42****,**58–67.Barcilon, A., J. C. Jusem, and S. Blumsack, 1980: Pseudo-adiabatic flow over a two-dimensional ridge.

,*Geophys. Astrophys. Fluid Dyn***16****,**19–33.Bruintjes, R. T., T. Clark, and W. D. Hall, 1994: Interactions between topographic airflow and cloud/precipitation development during the passage of a winter storm in Arizona.

,*J. Atmos. Sci***51****,**48–67.Collier, C. G., 1975: A representation of the effects of topography on surface rainfall within moving baroclinic disturbances.

,*Quart. J. Roy. Meteor. Soc***101****,**407–422.Daly, C., R. P. Neilson, and D. L. Phillips, 1994: A statistical–topographic model for mapping climatological precipitation over mountainous terrain.

,*J. Appl. Meteor***33****,**140–158.Durran, D. R., and J. B. Klemp, 1982a: On the effects of moisture on the Brunt–Väisälä frequency.

,*J. Atmos. Sci***39****,**2152–2158.Durran, D. R., and J. B. Klemp, 1982b: The effects of moisture on trapped mountain lee waves.

,*J. Atmos. Sci***39****,**2490–2506.Fraser, A. B., R. C. Easter, and P. V. Hobbs, 1973: A theoretical study of the flow of air and fallout of solid precipitation over mountainous terrain: Part I. Airflow model.

,*J. Atmos. Sci***30****,**801–812.Hobbs, P. V., R. C. Easter, and A. B. Fraser, 1973: A theoretical study of the flow of air and fallout of solid precipitation over mountainous terrain: Part II. Microphysics.

,*J. Atmos. Sci***30****,**813–823.Hutchinson, M. F., 1998: Interpolation of rainfall data with thin plate smoothing splines: II. Analysis of topographic dependence.

,*J. Geogr. Inf. Decis. Anal***2****,**168–185.Jiang, Q., 2003: Moist dynamics and orographic precipitation.

,*Tellus***55A****,**301–326.Jiang, Q., and R. B. Smith, 2003: Cloud timescales and orographic precipitation.

,*J. Atmos. Sci***60****,**1543–1559.Kuligowski, R. J., and A. F. Barros, 1999: High-resolution short term quantitative precipitation forecasting in mountainous regions using a nested model.

,*J. Geophys. Res***94****,**31553–31564.Lalas, D. P., and F. Einaudi, 1973: On the stability of a moist atmosphere in the presence of a background wind.

,*J. Atmos. Sci***30****,**795–800.Miglietta, M. M., and A. Buzzi, 2001: A numerical study of moist stratified flows over isolated topography.

,*Tellus***53A****,**481–499.Neiman, P. J., F. M. Ralph, A. B. White, D. E. Kingsmill, and P. O. G. Persson, 2002: The statistical relationship between upslope flow and rainfall in California's coastal mountain: Observations during CALJET.

,*Mon. Wea. Rev***130****,**1468–1492.Queney, P., 1947: Theory of perturbations in stratified currents with applications to airflow over mountain barriers. Dept. of Meteorology, Univ. of Chicago, Miscellaneous Rep. 23, 81 pp.

Rauber, R. M., 1992: Microphysical structure and evolution of a central Sierra Nevada orographic cloud system.

,*J. Appl. Meteor***31****,**3–24.Rhea, J. O., 1978: Orographic precipitation model for hydrometeorological use. Ph.D. thesis, Colorado State University, Dept. of Atmospheric Science, Atmospheric Science Paper 287, 198 pp.

Sawyer, J. S., 1962: Gravity waves in the atmosphere as a three-dimensional problem.

,*Quart. J. Roy. Meteor. Soc***88****,**412–425.Sinclair, M. R., 1994: A diagnostic model for estimating orographic precipitation.

,*J. Appl. Meteor***33****,**1163–1175.Smith, R. B., 1979: The influence of mountains on the atmosphere.

*Advances in Geophysics,*Vol. 27, Academic Press, 87–230.Smith, R. B., 2002: Stratified airflow over topography.

*Environmental Stratified Flows,*R. Grimshaw, Ed., Vol. 3,*Topics in Environmental Fluid Mechanics,*Kluwer, 119–159.Smith, R. B., 2003a: A linear upslope-time-delay model of orographic precipitation, mountain hydrology and water resources.

,*J. Hydrol***282****,**2–9.Smith, R. B., 2003b: Advection, diffusion and deposition from distributed sources.

,*Bound.-Layer Meteor***107****,**273–287.Smith, R. B., Q. Jiang, M. Fearon, P. Tabary, M. Dorninger, J. Doyle, and R. Benoit, 2003: Orographic precipitation and airmass transformation: An Alpine example.

,*Quart. J. Roy. Meteor. Soc***129****,**433–454.Wurtele, M. G., R. D. Sharman, and A. Datta, 1996: Atmospheric lee waves.

,*Annu. Rev. Fluid Mech***28****,**429–476.

## APPENDIX

### Simple Thermodynamics

*e*

_{s}

*T*

*e*

_{Sref}

*e*

^{LΔT/RυT2ref}

*e*

_{Sref}

*T*

_{ref}are the saturation vapor pressure and the temperature at the ground, and Δ

*T*is the difference between the air temperature and reference temperature. In (A2), the latent heat is

*L*= 2.5 × 10

^{6}J kg

^{−1}and the gas constant for vapor is

*R*

_{υ}= 461 J kg

^{−1}K

^{−1}. If the environmental lapse rate is

*γ,*the value of Δ

*T*=

*γz.*The vertical distribution of water vapor is [approximately, using (A1) and (A2) with (

*T*≈

*T*

_{ref})]

*ρ*

_{S}

*z*

*ρ*

_{Sref}

*e*

^{−z/Hw}

*ρ*

_{Sref}

*e*

_{S}(

*T*

_{ref})/

*RT*

_{ref}and

*γ*< 0, the column-integrated water vapor is

*ρ*

_{Sref}

*H*

_{w}and the horizontal flux of water vapor advected by a constant wind is

**F**

*ρ*

_{Sref}

*H*

_{w}

**U**

*x*-component

*F*

_{x}=

*ρ*

_{Sref}

*H*

_{w}

*U.*Note that, if the ambient lapse rate approaches isothermal, the water vapor scale height grows and so does the column water vapor and flux, for a given surface value of vapor density. Positive values of lapse rate

*γ*are not allowed here, as the column density of water vapor is not bounded.

*L*/(

*R*

_{υ}

*T*

_{ref}) ≈ 20 for typical conditions: (A6) can be written as a function of height using (A3)

*w*(

*z*), parcels cool at the rate

*w*Γ

_{m}. The local condensation rate is

*w*(

*dρ*

_{S}/

*dT*)Γ

_{m}. The column-integrated condensation rate is

*C*

_{w}

*ρ*

_{Sref}

_{m}

*γ*

_{m}/

*γ*> 1 if the atmosphere is statically stable.

*w*=

*U*(

*dh*/

*dx*), is constant with altitude, (A8) gives

*γ*= Γ

_{m}). The total rate of condensation over the windward slope is, from (A10),

*A*is the mountain height and

*x*

_{top}is the location of the hilltop. The ratio of the condensation rate to incoming vapor flux is [using (A5) and (A11)]

*γ.*

_{m}is an average moist adiabatic lapse rate. If we average Γ

_{m}over the first kilometer, surface temperature of

*T*

_{0}= 253, 263, 273, 283, and 293 K give values of about Γ

_{m}= −9, −8, −7, −6, −5 (°C km

^{−1}), respectively. If, for example, the environmental lapse rate is

*γ*= −4°C km

^{−1}and

*T*

_{0}= 273 K, (A13) gives

*N*

_{m}= 0.0104 s

^{−1}. The ratio Γ

_{m}/

*γ*= 1.75. The water vapor scale height (A4) in this case is

*H*

_{w}= 3.4 km.

The influence of wave dynamics on vertically integrated condensation. The cloud water source function *S*(*x*) (nondimensional part) is shown for two values of nondimensional moist-layer depth; *Ĥ* = 0 (dashed) and *Ĥ* = 1 (solid). (a) Monotonic rise to a plateau (27). (b) Witch-of-Agnesi ridge (33). Wave dynamics reduces the source function, shifts it upstream, and modifies its shape

Citation: Journal of the Atmospheric Sciences 61, 12; 10.1175/1520-0469(2004)061<1377:ALTOOP>2.0.CO;2

The influence of wave dynamics on vertically integrated condensation. The cloud water source function *S*(*x*) (nondimensional part) is shown for two values of nondimensional moist-layer depth; *Ĥ* = 0 (dashed) and *Ĥ* = 1 (solid). (a) Monotonic rise to a plateau (27). (b) Witch-of-Agnesi ridge (33). Wave dynamics reduces the source function, shifts it upstream, and modifies its shape

Citation: Journal of the Atmospheric Sciences 61, 12; 10.1175/1520-0469(2004)061<1377:ALTOOP>2.0.CO;2

The influence of wave dynamics on vertically integrated condensation. The cloud water source function *S*(*x*) (nondimensional part) is shown for two values of nondimensional moist-layer depth; *Ĥ* = 0 (dashed) and *Ĥ* = 1 (solid). (a) Monotonic rise to a plateau (27). (b) Witch-of-Agnesi ridge (33). Wave dynamics reduces the source function, shifts it upstream, and modifies its shape

Citation: Journal of the Atmospheric Sciences 61, 12; 10.1175/1520-0469(2004)061<1377:ALTOOP>2.0.CO;2

Precipitation rates (mm h^{−1}) over a triangle ridge with three different assumptions: precipitation patterns with no dynamics or cloud delays (3) (solid); dynamics only (16) (dotted); dynamics and cloud delays (49) (dashed). Two mountain half-widths are shown: (a) *a* = 15 km, (b) *a* = 40 km. Parameters are *U* = 15 m s^{−1}, *N*_{m} = 0.005 s^{−1}, *T*_{0} = 280 K, *H*_{w} = 2500 m, *τ*_{c} = *τ*_{f} = 1000 s, *A* = 500 m, *S*_{∞} = 0. Note the great reduction in total precipitation caused by dynamics and cloud delays. The narrower hill has a precipitation maximum close to the hilltop and a lower overall precipitation efficiency

Precipitation rates (mm h^{−1}) over a triangle ridge with three different assumptions: precipitation patterns with no dynamics or cloud delays (3) (solid); dynamics only (16) (dotted); dynamics and cloud delays (49) (dashed). Two mountain half-widths are shown: (a) *a* = 15 km, (b) *a* = 40 km. Parameters are *U* = 15 m s^{−1}, *N*_{m} = 0.005 s^{−1}, *T*_{0} = 280 K, *H*_{w} = 2500 m, *τ*_{c} = *τ*_{f} = 1000 s, *A* = 500 m, *S*_{∞} = 0. Note the great reduction in total precipitation caused by dynamics and cloud delays. The narrower hill has a precipitation maximum close to the hilltop and a lower overall precipitation efficiency

Precipitation rates (mm h^{−1}) over a triangle ridge with three different assumptions: precipitation patterns with no dynamics or cloud delays (3) (solid); dynamics only (16) (dotted); dynamics and cloud delays (49) (dashed). Two mountain half-widths are shown: (a) *a* = 15 km, (b) *a* = 40 km. Parameters are *U* = 15 m s^{−1}, *N*_{m} = 0.005 s^{−1}, *T*_{0} = 280 K, *H*_{w} = 2500 m, *τ*_{c} = *τ*_{f} = 1000 s, *A* = 500 m, *S*_{∞} = 0. Note the great reduction in total precipitation caused by dynamics and cloud delays. The narrower hill has a precipitation maximum close to the hilltop and a lower overall precipitation efficiency

Planform patterns of condensation and/or precipitation for a circular Gaussian hill: (a) condensation source (*S*) with no dynamics or cloud delays (3); (b) condensation source with dynamics (16); (c) precipitation with dynamics and cloud time-delays (49). The parameters of the calculation are the same as in Fig. 3, but with *a* = 15 km, Δ*x* = Δ*y* = 750 m. The contour intervals differ for the different parts: (a) 2 mm h^{−1} from 0.025 to 12.025, (b) 2 mm h^{−1} from −10 to +6, and (c) 0.4 mm h^{−1} from 0.025 to 2.025. Note that the actual precipitation (c) is less intense but more widespread than the raw upslope pattern (a). The precipitation field has a maximum near the hilltop and a weak region of wave cloud precipitation is located downwind (*x* = 75 km)

Planform patterns of condensation and/or precipitation for a circular Gaussian hill: (a) condensation source (*S*) with no dynamics or cloud delays (3); (b) condensation source with dynamics (16); (c) precipitation with dynamics and cloud time-delays (49). The parameters of the calculation are the same as in Fig. 3, but with *a* = 15 km, Δ*x* = Δ*y* = 750 m. The contour intervals differ for the different parts: (a) 2 mm h^{−1} from 0.025 to 12.025, (b) 2 mm h^{−1} from −10 to +6, and (c) 0.4 mm h^{−1} from 0.025 to 2.025. Note that the actual precipitation (c) is less intense but more widespread than the raw upslope pattern (a). The precipitation field has a maximum near the hilltop and a weak region of wave cloud precipitation is located downwind (*x* = 75 km)

Planform patterns of condensation and/or precipitation for a circular Gaussian hill: (a) condensation source (*S*) with no dynamics or cloud delays (3); (b) condensation source with dynamics (16); (c) precipitation with dynamics and cloud time-delays (49). The parameters of the calculation are the same as in Fig. 3, but with *a* = 15 km, Δ*x* = Δ*y* = 750 m. The contour intervals differ for the different parts: (a) 2 mm h^{−1} from 0.025 to 12.025, (b) 2 mm h^{−1} from −10 to +6, and (c) 0.4 mm h^{−1} from 0.025 to 2.025. Note that the actual precipitation (c) is less intense but more widespread than the raw upslope pattern (a). The precipitation field has a maximum near the hilltop and a weak region of wave cloud precipitation is located downwind (*x* = 75 km)

Three efficiency ratios PE_{dyn} (dashed), PE_{cloud} (dotted), and PE (solid) as a function of mountain half-width for a Gaussian ridge according to the current linear model (12), (49). Parameters are the same as in Figs. 3 and 4

Three efficiency ratios PE_{dyn} (dashed), PE_{cloud} (dotted), and PE (solid) as a function of mountain half-width for a Gaussian ridge according to the current linear model (12), (49). Parameters are the same as in Figs. 3 and 4

Three efficiency ratios PE_{dyn} (dashed), PE_{cloud} (dotted), and PE (solid) as a function of mountain half-width for a Gaussian ridge according to the current linear model (12), (49). Parameters are the same as in Figs. 3 and 4

A real-world application of the linear FFT model; the Olympic Mountains under a 15 m s^{−1} southwesterly airflow. The 6-h accumulated precipitation is shown shaded with a 2.5-mm contour interval. The maximum precipitation is 25.96 mm. The map projection is Lambert conformal. The atmospheric parameters in the model calculation are the same as in Figs. 3 and 4. The terrain is shown (dotted) with a 200-m contour interval. The coastline is shown with a dark solid line

A real-world application of the linear FFT model; the Olympic Mountains under a 15 m s^{−1} southwesterly airflow. The 6-h accumulated precipitation is shown shaded with a 2.5-mm contour interval. The maximum precipitation is 25.96 mm. The map projection is Lambert conformal. The atmospheric parameters in the model calculation are the same as in Figs. 3 and 4. The terrain is shown (dotted) with a 200-m contour interval. The coastline is shown with a dark solid line

A real-world application of the linear FFT model; the Olympic Mountains under a 15 m s^{−1} southwesterly airflow. The 6-h accumulated precipitation is shown shaded with a 2.5-mm contour interval. The maximum precipitation is 25.96 mm. The map projection is Lambert conformal. The atmospheric parameters in the model calculation are the same as in Figs. 3 and 4. The terrain is shown (dotted) with a 200-m contour interval. The coastline is shown with a dark solid line

Normalized upslope centerline condensation (36) influenced by hydrostatic airflow dynamics over a Gaussian hill

Normalized total upslope condensation (37) for a circular Gaussian hill (*R* = 1)

Some symbols used in the model

Nondimensional control variables for orographic precipitation

Influence of time constants on precipitation efficiencies for the Gaussian ridge (case *U* = 15 m s^{−1} , *N* = 0.005 s^{−1} , *a* = 15 km, *H _{w}* = 3 km,

*S*

_{∞}= 0)