## 1. Introduction

Gravity waves in the atmosphere and ocean have been studied extensively, and much is known about how they propagate, evolve, and affect the ambient flow [e.g., see texts by Lighthill (1978), Gill (1982), Lindzen (1990), and Nappo (2002) for a review of basic gravity wave dynamics]. Particular attention has been devoted to understanding how gravity waves interact with the mean flow in the stratosphere where vertically propagating waves may impose a drag due to convergence of momentum flux [see review by Fritts and Alexander (2003) and references therein]. A problem that is complimentary to the wave–mean flow interaction problem is the wave–forcing interaction problem. The latter has received relatively less scrutiny in the literature (exceptions noted below) but may be of significance to mesoscale dynamics in the troposphere. The wave–forcing interaction problem is of particular relevance to convectively generated waves that, unlike topographically forced waves, originate from sources which are themselves dynamic. Given that space and time scales characterizing waves and convection are similar, source and spectrum may be mutually dependent.

For convenience, the relationship between gravity waves and their sources is usually analyzed as though it proceeds in one direction. Either one analyzes the waves generated by convection, which are assumed to have no subsequent influence on the convection itself (e.g., Bretherton and Smolarkiewicz 1989; Nicholls et al. 1991; Chagnon and Bannon 2005a,b), or one examines the role of waves in the initiation of convection but not in the subsequent evolution of convection (e.g., Marsham and Parker 2006). There have been a few exceptions in which the mutual interaction has been considered. The wave-CISK hypothesis of Lindzen (1974) could be considered the most extreme type of interaction in which no dynamical distinction between waves and convection is made. Cohen and Craig (2004) hypothesized that gravity waves play an important role in establishing equilibrium spatial and temporal scales in fields of convection simulated using a cloud-resolving model but stopped short of identifying the actual mechanisms involved. Lin (1987) demonstrated that the linear wave response to a diabatic heating applied in a sheared flow could produce a maximum vertical velocity that is either in or out of phase with the heating, depending on the Richardson number of the background flow. Most recently, Robinson et al. (2008) found that the amplitude and depth of the response to patches of surface heating can be optimized for specific horizontal scales of the heating that excite a resonant response. A question that follows from these studies is to what extent an evolving internal forcing such as convection might be influenced by interfering tendencies between dynamical response and forcing.

This paper considers the relationship between high frequency waves and the forcing from which they originate. The effect of waves on forcing (and vice versa) may be thought of as a problem complimentary to the effect of the environment (e.g., wind shear) on waves. As discussed above, the latter has received much attention for its significance to stratospheric dynamics, while the former is relevant to the dynamics of the troposphere. The problem examined in this paper considers a wave forcing mechanism such as that due to the release of a local instability (e.g., convection). The nature of this instability and the origin of the forcing mechanism are not addressed in this paper. If the environment is statically stable, then the forcing will generate a dynamical response involving waves. If the response induces dynamical perturbations in the vicinity of the forcing, then the response may subsequently do work on the forcing itself, thereby acting as either a resistance or conductance to the forcing. Under statically stable conditions, the effect is usually a resistance. This dynamical resistance depends on the properties of the forcing. A forcing that generates a wave response that minimizes the dynamical resistance imposed by the resulting wave field may be characterized as efficient. The analysis contained in this paper explores the effect of the dynamical response on the energetic efficiency of the forcing across a wide parameter space describing the forcing and environment. Although the scope of this paper is limited to the introduction and basic demonstration of this concept, the analysis generalizes the specific examples of enhanced forcing efficiency provided by Lin (1987) and Robinson et al. (2008).

In the next section a linear model is introduced that will allow for a computation of the energetic efficiency of a forcing that excites acoustic and gravity waves in a compressible, inviscid, nonhydrostatic atmosphere linearized about a basic state with arbitrary background flow and static stability. We will not assume to know anything about the actual process that generates the waves. Instead, we will examine the dependence of the forcing efficiency and dynamical resistance on the properties of the forcing as well as the properties of the environment into which the forcing is applied. Section 3 will demonstrate several examples of how energetic efficiency varies with forcing geometry against a background flow. Specifically, experiments are designed to test the effect of forcing duration, forcing-relative flow speed, wind shear, and multiple versus isolated forcing elements. Additionally, total energy calculations are performed along gravity wave particle ray paths in order to clarify some of the results obtained in the linear model simulations. Section 4 discusses the results and suggests routes for further study.

## 2. Forcing efficiency in a linear model

A linear model will be used to examine gravity waves, dynamical resistance, and forcing efficiency throughout this investigation. The eigensolutions of the linear system comprise acoustic and gravity waves. Waves are excited by arbitrary external forcing terms in the momentum and heat equations. These forcing terms may be associated with “fast” and/or diabatic processes taking place within a cumulus cloud. It will be shown that the total energy excited by the time-dependent forcing is dependent on the work performed by the waves upon the forcing.

*f*plane, assuming horizontal homogeneity in the

*y*direction, linearized about a hydrostatic basic state (denoted by a subscript

*s*) with arbitrary velocity

*u*areandwhere the operator ∂/∂

_{s}*t*is defined as

_{a}*H*=

_{s}*RT*/

_{s}*g*is the scale height,

*R*= 287 J kg

^{−1}K

^{−1}is the dry ideal gas constant,

*g*= 9.816 m s

^{−2}is the vertical component of acceleration due to gravity,

*F*,

_{u}*F*,

_{υ}*F*, and

_{w}*F*. Although these forcing terms are likely to be interdependent in the real atmosphere (e.g., Song et al. 2003), the present investigation imposes no such constraints (see section 4 for a discussion of this issue).

_{θ}*F*, but the procedure may be generalized to any other type of forcing. Suppose the forcing is turned on at time

_{w}*t*= 0 and shut off at time

*t*=

*τ*. Let brackets denote a spatial integral over the full domain. The total energy generation isIntegrating the vertical momentum equation (3), the vertical velocity perturbations may generally be written aswhere the function

*D*represents the effects of the internal dynamics on the evolution of the velocity perturbations. Substituting (15) into (14) yieldswhereandFor the impulsive forcing (Dirac delta function in time), the total energy generation

_{w}*E*

_{gen}can be computed directly from

*E*

_{fcg}. For the finite duration forcing (

*τ*> 0)

*E*

_{gen}cannot be computed without first considering the dynamical response

*E*

_{dyn}, which may either amplify or impede the total energy generation relative to

*E*

_{fcg}. If the dynamical response

*D*is directed opposite to the forcing

_{w}*F*, then the dynamical resistance

_{w}*E*

_{dyn}will be positive, thus reducing the total energy generated by the forcing. We may also define a forcing efficiency:Note that an efficiency

*ϵ*= 1 corresponds to the case in which there is no dynamical response during the application of the forcing. An impulsive forcing has

*ϵ*= 1. Under statically stable conditions the forcing efficiency is typically less than one. Nonetheless, configurations of the basic state and forcing exist that minimize the resistance and maximize the efficiency.

The following section demonstrates the forcing efficiency *ϵ* as a function of properties of the forcing and the basic state. Although a solution is known for (1)–(5) in the special case that the initial state is isothermal and at rest (Chagnon and Bannon 2005a,b), a solution for the general case is unknown. We therefore compute *E*_{dyn} by solving the system (1)–(5) numerically in the following way.

- A discrete Fourier transform is applied to (1)–(5), resulting in system of equations that depend on horizontal wavenumber
*k*as well as height*z*and time*t*. - Vertical derivatives are represented as centered finite differences accurate to fifth order (i.e., using two points on either side of the central level) except near the boundaries where forward and backward representations of second-order accuracy are used.
- The resulting time-dependent equations are solved using a third-order Adams–Bashforth scheme, except at the first time step when a forward step is used.
- An inverse fast Fourier transform is finally applied to the horizontal wavenumber solution to obtain the full solution in physical space.

*z*= 0 and 25 km. Over the range of forcing parameters used in this investigation, the upper boundary is sufficiently far from the tropospheric wave source to ensure that no waves reflected from the upper boundary can interfere with the forcing at a later time. For example, waves generated at

*z*= 5 km elevation would have to travel at a vertical speed exceeding 30 m s

^{−1}in order to traverse the necessary 40-km distance to arrive back at the forcing elevation within 20 min. Such a speed is larger than the fastest vertical group speeds (approximately 25 m s

^{−1}) associated with the deepest nonhydrostatic modes (i.e., vertical wavelength 20 km) propagating against a basic state with

*N*= 0.02 s

_{s}^{−1}. Note also that acoustic waves contribute insignificantly to solutions when the forcing duration exceeds approximately 3 min (Chagnon and Bannon 2005a). The numerical solution procedure conserves the total perturbation energy to within one-tenth of a percent for simulations of one hour duration in which the shear is zero.

An advantage of using total energy (and the energetic efficiency *ϵ*) as a basis for analysis is that it provides one single number to characterize an entire experiment evolving in space and time. It is therefore a convenient diagnostic for exploring a vast parameter space, as is done in the following section. Furthermore, total energy is conserved following termination of forcing, except for the case in which there is a background shear. If the environment contains shear, then the total energy of wave perturbations is not conserved due to the work done by/against the mean flow. We proceed with this caveat in mind, noting that the convenience of using the energetic efficiency to investigate dynamical resistance justifies its application here.

## 3. Forcing efficiency examples

*x*,

_{i}*z*of the formwhere

_{i}*S*(

*x*,

*z*) is a function describing the geometry of the forcing,

*T*(

*t*) is a function describing the time dependence of the forcing, and

*δF*is the forcing amplitude. The forcing geometry is Gaussian:where the parameters

*a*and

*d*describe the width and depth of the forcing, respectively. The time dependence is a sine-squared function:where

*τ*is the forcing duration and

*H*is the Heaviside step function. Note that the Gaussian forcing (21) used in the experiments below is consistent with both upper and lower boundary conditions. The heating is centered in the domain such that the value of

*S*(

*x*,

*z*) is the same on the left and right lateral boundaries, thus avoiding problems with the discrete Fourier transform. The vertical velocity and heating are fixed at zero on the lower boundary to avoid spurious wave generation at the boundary; in all experiments, the heating elevation

*z*is chosen to be at least twice the heating depth

_{i}*d*so as to avoid large heating gradients near the lower boundary. A Gaussian heating geometry has been chosen rather than heating of finite extent in order to avoid the Gibbs phenomenon and the spurious generation of grid-scale waves.

Three sets of experiments are performed in this section. The first is a control experiment in which the forcing is stationary and there is no background wind. The second set of experiments considers the efficiency of a forcing that is moving relative to the background flow, including cases with and without environmental shear. The third experiment, presented at the end of this section, examines an ensemble of forcing elements (as opposed to a single isolated forcing). In the case of multiple forcing elements, the dynamical response from nearby forcing events may interfere with neighboring forcing elements. All experiments presented below correspond to the case where *N _{s}* = 0.02 s

^{−1}is constant. Experiments in which

*N*varied with elevation were performed (e.g., “tropopause-like” doubling of

_{s}*N*; reduced stability layer above forcing) but did not differ in significant qualitative respects from the results shown below. More exotic experiments such as examination of wave ducts in the vicinity of the forcing would constitute an interesting basis for investigation but are deferred to future studies in favor of presenting as simple and basic an analysis as possible for the present.

_{s}### a. Control experiment: Isolated forcing, no background wind

The demonstration of forcing efficiency begins with the simple experiment in which waves are generated by an isolated heating applied to a resting atmosphere. The control experiment demonstrates the dependence of total energy generation on the duration of forcing. Figure 1 presents the energy efficiency as a function of *τ* during two stages of adjustment—one stage lasting only a few minutes during which acoustic waves are generated, followed by a second stage lasting tens of minutes characterized by the emergence of high-frequency gravity waves. A longer duration forcing generates less total energy than a forcing applied rapidly. Similar results have been demonstrated analytically (see Sotack and Bannon 1999; Chagnon and Bannon 2005a). The loss of energy at longer duration is due to the destructive interference of waves whose period is shorter than the time scale of the forcing. Such behavior is typical in statically stable systems for which the dynamical response is normally directed against an imposed forcing. To explain this behavior, Vadas and Fritts (2001) made an analogy to a system comprising a pendulum fixed to a moveable base. If the base is accelerated rapidly from a resting position (i.e., on a time scale shorter than the pendulum period), then an oscillation may be excited. If the forcing time scale is longer than the pendulum period, then higher frequency oscillations will interfere destructively with the forcing, resulting in an oscillation of lesser amplitude. In the following subsection it is demonstrated that these destructive influences can be overcome if the forcing and environment are configured such that the change in wave phase is minimized relative to the forcing.

### b. Forcing against a background flow

This subsection considers the effect of a forcing-relative background flow *u _{s}* on

*ϵ*of an isolated forcing. We first examine an unsheared constant

*u*=

_{s}*U*and then examine a flow having constant shear Λ. The unsheared experiment corresponds to the case where the forcing is moving through a static fluid at a speed

*U*. Figure 2 presents the forcing efficiency for varying background flow speed as a function of

*k*. To produce the field contoured in Fig. 2, the total energy following a prescribed heating was decomposed into its spectral contributions via Fourier transform along

*k*. The efficiency (19) was then calculated as the total energy in

*k*normalized by the Fourier transformed integral (18) evaluated at

*k*. The parameter

*d*was varied between a shallow case (500 m, Fig. 2a) and a moderate-depth case (2 km, Fig. 2b). The heating width, elevation, and duration were fixed at

*a*= 5 km,

*z*= 4 km, and

_{i}*τ*= 20 min, respectively.

Figure 2 demonstrates that a moving forcing is potentially more efficient than a stationary forcing. The enhancement occurs at larger forcing-relative flow speeds *U* when *k* is decreased. In the shallow heating case (Fig. 2a) the envelope of enhanced efficiency covers a wide range of speeds. In the moderate depth case (Fig. 2b) the envelope of enhanced efficiency is narrower and confined to higher flow speeds. Because the solution at each wavenumber shown in Fig. 2 contains contributions from both positive and negative vertical wavenumbers and frequencies (and hence wave particles propagating in all quadrants relative to the source), the efficiency contoured in Fig. 2 is symmetric about the *k* and *U* axes. Consequently, only one quadrant is shown in Fig. 2.

_{ray}. If the particle phase varies slowly along the ray path, then the wave particle performs net work in the direction of the forcing. This net work results in a larger quantity of total energy generated along the ray path. In this scenario, TE

_{ray}is positive. If the particle undergoes one or more full cycles along the path, then the particle performs zero net work in the direction of the forcing and TE

_{ray}is zero. If the particle performs net work directed opposite to the forcing, then TE

_{ray}is negative. Consider gravity waves in an incompressible atmosphere satisfying the dispersion relation:where

*ω*is the flow-relative wave frequency, and

*m*is the vertical wavenumber. For a constant background

*U*, the total change in phase along a ray path (assuming constant

*N*) iswhere

_{s}*d*

_{fcg}is the vertical path length through the forcing region,

*c*is the vertical component of group velocity, and

_{gz}*c*is the horizontal component of the phase velocity. The total work performed by the wave particle along the ray path in the direction of the forcing, found by evaluating (A5), iswhere

_{x}*A*is a constant proportional to the forcing magnitude and initial wave amplitude. Figure 3 contours TE

_{ray}given by (25) as a function of horizontal wavenumber

*k*and

*U*for a range of initial particle vertical wavelengths

*λ*= 2

_{z}*π*/

*m*and forcing regions depths

*d*

_{fcg}. All curves have been normalized by the coefficient in (25) such that a value of one corresponds to stationary phase along the ray path. Only that quadrant corresponding to wave particles moving into the background flow is shown in Fig. 3. Features common to all of the panels in Fig. 3 include (i) a primary envelope of enhanced efficiency corresponding to minimal total phase change whose width and mean value along the

*U*axis increase for decreasing

*k*and (ii) additional bands of maxima and minima flanking the primary envelope, which correspond to particles that have undergone multiple cycles (e.g., phase given by

*nπ*). Comparing the panels in the columns from right to left reveals that the width of the envelope is larger for a shallower forcing (i.e., shallower particle path). For a shallow forcing, the particle can exit the forcing region at an earlier time in its phase cycle and is therefore less likely to interfere destructively with the forcing. Comparing the panels in the rows from top to bottom reveals that the maximum efficiency occurs at higher

*U*for particles having deeper

*λ*. Waves with deep vertical wavelength have a larger horizontal phase speed and therefore achieve stationary phase at higher forcing relative flow speeds. These features of the analytic wave particle calculation shown in Fig. 3 describe the general characteristics of the efficiency calculated in the linear model simulations (Fig. 2).

_{z}*z*represents the elevation of the heating maximum. Figure 4 presents the forcing efficiency for varying Λ as a function of

_{i}*k*for a shallow heating (

*d*= 500 m, Fig. 4a) and a moderate depth heating (

*d*= 2 km, Fig. 4b). In the shallow heating case (Fig. 4a) the shear has very little impact on the forcing efficiency after 20 minutes. The forcing efficiency is a decreasing function of wavenumber. Note that higher wavenumbers correspond to wave particles having higher frequencies, which undergo a complete cycle in a shorter period of time and may therefore interfere destructively with the forcing at an earlier time than particles of lower wavenumber. In the moderate depth case (Fig. 4b) bands of both maximum and minimum efficiency are evident that slope upward toward larger shear magnitude at lower horizontal wavenumber. As in Fig. 2, the efficiency contoured in Fig. 4 is symmetric about

*k*and

*U*; hence, only one quadrant is shown.

*U*is the forcing-relative flow speed at the wave particle’s launch elevation and erf the error function. Figure 5 contours the total energy given by (27) as a function of Λ and horizontal wavenumber for a range of

_{o}*d*

_{fcg}and

*λ*. Both positive and negative shear values are shown in Fig. 5 because a significant response is evident in both regions of the parameter space (note again that the spectra from the linear model simulations in Fig. 4 contain contributions from both halves). The efficiency along ray paths in the shear case (Fig. 5) shares many properties in common with the constant background flow case (Fig. 3). For example, a deeper heating (larger

_{z}*d*

_{fcg}) is more destructive than a shallow heating; hence, values of enhanced efficiency occur across a narrower envelope of shear magnitude. Furthermore, a particle having deeper

*λ*experiences enhanced efficiency at larger values of shear. The differences between the effects of shear and a constant background flow are essentially due to the differences in ray path geometries. In the shear case the ray paths are curved. Particles may therefore spend unequal proportions of their lifetime at particular elevations within the forcing. Although the particles are assumed to initially be in phase with the forcing, it is possible for some particles to spend enough time having phase directed against the forcing that the quantity TE

_{z}_{ray}may become negative. Note that this does not imply that the total energy is negative; TE

_{ray}is a measure of the work performed by the wave particle in the direction of the forcing and is thus related to −

*E*

_{dyn}.

### c. Ensemble of multiple forcing elements

*a*.

_{c}Figure 6 plots the efficiency of an ensemble of elements normalized by the efficiency of an isolated forcing for the three experiments described above. The first experiment (solid curves, Fig. 6a) demonstrates that an ensemble of warming elements is more efficient than an isolated warming. Each warming element generates a wave response that is dominated by a horizontally propagating signal of positive potential temperature anomaly (e.g., see Bretherton and Smolarkiewicz 1989). These leading fronts interfere constructively with neighboring warming elements upon arrival. Closely spaced, the elements are conducive to constructive interference and efficiency enhancement. The second experiment (dashed curves, Fig. 6a) demonstrates that the efficiency of an ensemble of alternating cooling and warming elements is reduced relative to the isolated forcing. In this case the leading wave fronts are out of phase and therefore interfere destructively with their nearest neighboring elements upon arrival. The third experiment (Fig. 6b) demonstrates the efficiency of an ensemble of hybrid elements. The response contains both positive and negative anomalies propagating through forcings of both positive and negative sign. The forcing efficiency of the hybrid ensemble can be either enhanced or reduced relative to the isolated forcing, depending on the configuration of the ensemble and the time into the simulation. Figure 6 demonstrates generally that an ensemble of multiple forcing elements can be either more or less efficient than an isolated forcing depending on the phase of the wave response arriving at neighboring elements. The interaction is most pronounced for ensembles characterized by a large degree of symmetry.

## 4. Discussion

The purpose of this paper has been to demonstrate the concepts of dynamical resistance and forcing efficiency as a means of analyzing the relationship between gravity waves and their sources. If a forcing is applied on a time scale similar to that of the waves it excites, then the dynamical perturbations associated with the waves may interfere constructively or destructively with the forcing itself. It is hypothesized that this effect, identified as the dynamical resistance, may explain some properties of the gravity wave energy spectrum corresponding to a given forcing. If there is little or no resistance to the forcing, then the forcing can generate more energy. Depending on the characteristics of the environment into which the forcing is applied, waves may remain nearly in phase with the forcing. For example, if a forcing is moving through a fluid at a constant speed, then gravity waves having phase speed similar to the forcing-relative flow speed will exert minimal resistance upon the forcing. The significance of this phase locking behavior to gravity waves and convection in nature might be viewed from either one of two possible perspectives. These perspectives, described below, represent two extremes concerning the extent to which waves could interact with their sources.

In the first perspective, we suppose that the evolving characteristics of convection are entirely unaffected by the dynamical resistance imposed by the convectively generated waves. In this scenario, peaks in the wave energy spectrum would correspond to the scales at which gravity waves are optimally phase locked relative to the forcing. This result would challenge the typical assumption that the initial spectrum of waves generated by convection is determined exclusively by the spatial characteristics of the forcing alone. Similarly, Holton et al. (2002) argued that the dominant vertical scale of waves generated by convection depends not only on the depth of the cloud but also on its duration. Chagnon and Gray (2008) have recently identified the existence of multiple discrete peaks in the gravity wave spectrum above convection both in observations from a wind-profiling Doppler radar and in high-resolution convection-permitting numerical simulations. These peaks could not be explained as having originated due to a parametric instability. As part of a future investigation, we are examining the possibility that these spectral peaks correspond to scales at which gravity wave particles experience minimal phase change relative to the background forcing-relative flow near the tropopause.

In the second perspective, we suppose that an evolving field of convection responds to the dynamical resistance by adjusting toward a more efficient configuration. This is the perspective assumed by Robinson et al. (2008) and by Cohen and Craig (2004) to explain the dominant scales present in ensembles of convection. Figure 7 schematically illustrates one example of how a dynamical resistance might influence the structure of an isolated forcing. Suppose a cloud has formed due to the release of some available energy associated with an instability, denoted *E*_{ins} (which might be related to, say, the convective available potential energy in a cloud or the production of subcloud layer entropy). Suppose that the instability induces maximum acceleration along a direction *α*_{ins}. As this instability produces motion it will excite a dynamical response in the surrounding environment. The dynamical response may in turn do work upon the region containing the instability, denoted *E*_{dyn}. An interesting question to consider is whether there exists a forcing angle *α* that allows the forcing to draw energy from the instability while minimizing the resistance from the dynamical response. Could this angle deviate from *α*_{ins}?

The scope of this investigation has been limited to the introduction and basic demonstration of dynamical resistance, and the context has been highly idealized. Nature does not obey the linear, 2D, inviscid dynamics that have been explored in this paper. A more general formulation of the concepts of dynamical resistance and forcing efficiency are needed before we can assess their relevance to convection and gravity waves in the real world. An attractive candidate is the Eulerian available energy formalism described in Bannon (2005) and applied to the analysis of atmospheric adjustment in Edson and Bannon (2008).

A major shortcoming of the analysis contained in this paper is that no constraints have been placed on the evolution and character of the forcing. In fact, the analysis has made no reference to the actual origin of the forcing at all. An interesting next step might involve the optimization of energetic efficiency subject to a few such constraints. For example, Song et al. (2003) examined the gravity wave forcing in fully nonlinear cloud-resolving numerical simulations. In addition to diabatic heating, the wave forcing was shown to comprise significant contributions from nonlinear heat and momentum fluxes. These nonlinear momentum and heat flux terms were shown to be out of phase with the diabatic heating, thus reducing the total wave forcing. Chun et al. (2008) have proposed that the reduction in total wave forcing could be incorporated into gravity wave drag schemes via a “nonlinearity factor” multiplying the heat and momentum flux terms; a large nonlinearity factor was shown to yield a reduction in wave amplitude above the forcing. In light of these studies, a prudent application of the dynamical resistance calculation would be to explore the efficiency of a forcing comprised of both diabatic and nonlinear heat and momentum flux terms. Furthermore, the nonlinear heat and momentum flux terms may incorporate the effects of waves themselves. Future formulations of the dynamical resistance and forcing efficiency should accommodate the direct effect of waves on the forcing.

## Acknowledgments

This work was supported by the UK National Centre for Atmospheric Science under the Weather Directorate. The author is grateful to Michael Reeder and Peter Bannon for helpful discussions. The author wishes to thank the two anonymous reviewers whose feedback led to a significant improvement on the original draft.

## REFERENCES

Bannon, P. R., 2005: Eulerian available energetics in moist atmospheres.

,*J. Atmos. Sci.***62****,**4238–4252.Bretherton, C. S., , and P. K. Smolarkiewicz, 1989: Gravity waves, compensating subsidence, and detrainment around cumulus clouds.

,*J. Atmos. Sci.***46****,**740–759.Chagnon, J. M., , and P. R. Bannon, 2005a: Wave response during hydrostatic and geostrophic adjustment. Part I: Transient dynamics.

,*J. Atmos. Sci.***62****,**1311–1329.Chagnon, J. M., , and P. R. Bannon, 2005b: Wave response during hydrostatic and geostrophic adjustment. Part II: Potential vorticity conservation and energy partitioning.

,*J. Atmos. Sci.***62****,**1330–1345.Chagnon, J. M., , and S. L. Gray, 2008: Analysis of convectively-generated gravity waves in mesoscale model simulations and wind-profiler observations.

,*Quart. J. Roy. Meteor. Soc.***134****,**663–676. doi:10.1002/qj.239.Chun, H-Y., , H-J. Choi, , and I-S. Song, 2008: Effects of nonlinearity on convectively forced internal gravity waves: Application to a gravity wave drag parameterization.

,*J. Atmos. Sci.***65****,**557–575.Cohen, B. G., , and G. C. Craig, 2004: The response time of a convective cloud ensemble to a change in forcing.

,*Quart. J. Roy. Meteor. Soc.***130****,**933–944.Edson, A. R., , and P. R. Bannon, 2008: Nonlinear atmospheric adjustment to momentum forcing.

,*J. Atmos. Sci.***65****,**953–969.Fritts, D. C., , and M. J. Alexander, 2003: Gravity wave dynamics and effects in the middle atmosphere.

,*Rev. Geophys.***41****,**1003. doi:10.1029/2001RG000106.Gill, A. E., 1982:

*Atmosphere–Ocean Dynamics*. Academic Press, 662 pp.Holton, J. R., , J. H. Beres, , and X. Zhou, 2002: On the vertical scale of gravity waves excited by localized thermal forcing.

,*J. Atmos. Sci.***59****,**2019–2023.Lighthill, J., 1978:

*Waves in Fluids*. Cambridge University Press, 516 pp.Lin, Y-L., 1987: Two-dimensional response of a stably stratified shear flow to diabatic heating.

,*J. Atmos. Sci.***44****,**1375–1393.Lindzen, R. S., 1974: Wave-CISK in the tropics.

,*J. Atmos. Sci.***31****,**156–179.Lindzen, R. S., 1990:

*Dynamics in Atmospheric Physics*. Cambridge University Press, 310 pp.Marsham, J. H., , and D. J. Parker, 2006: Secondary initiation of multiple bands of cumulonimbus over southern Britain. II: Dynamics of secondary initiation.

,*Quart. J. Roy. Meteor. Soc.***132****,**1053–1072. doi:10.1256/qj.05.152.Nappo, C. J., 2002:

*An Introduction to Atmospheric Gravity Waves*. International Geophysics Series, Vol. 85, Academic Press, 276 pp.Nicholls, M., , R. Pielke, , and W. Cotton, 1991: Thermally forced gravity waves in an atmosphere at rest.

,*J. Atmos. Sci.***48****,**1869–1884.Robinson, F. J., , S. C. Sherwood, , and Y. Li, 2008: Resonant response of deep convection to surface hot spots.

,*J. Atmos. Sci.***65****,**276–286.Song, I-S., , H-Y. Chun, , and T. P. Lane, 2003: Generation mechanisms of convectively forced internal gravity waves and their propagation to the stratosphere.

,*J. Atmos. Sci.***60****,**1960–1980.Sotack, T., , and P. Bannon, 1999: Lamb’s hydrostatic adjustment for heating of finite duration.

,*J. Atmos. Sci.***56****,**71–81.Vadas, S. L., , and D. C. Fritts, 2001: Gravity wave radiation and mean responses to local body forces in the atmosphere.

,*J. Atmos. Sci.***58****,**2249–2279.

## APPENDIX

### Wave Particle Perspective

*x*

_{0},

*z*

_{0}at an initial time. The particle follows a trajectory (i.e., a ray path), defined by the parametric curve

*x*(

_{p}*t*),

*z*(

_{p}*t*), traveling through the broader region of forcing

*S*(

*x*,

*z*). Relative to the background flow, the particle has a horizontal phase speed

*c*. Let us consider the effect of this wave particle on the energetic efficiency of the forcing

_{x}*F*. The field perturbations produced by the wave particle along the ray are characterized by a wave functionwhere the phase is given byin which

*k*and

*m*are horizontal and vertical wavenumbers,

*ω*is the wave frequency,

*δx*(

_{p}*t*) =

*x*(

_{p}*t*) −

*x*

_{0}, and

*δz*(

_{p}*t*) =

*z*(

_{p}*t*) −

*z*

_{0}. We assume that the medium through which the particle propagates only varies in the

*z*direction. Consequently,

*k*and

*ω*may be treated as constants that are fixed for the particle, but

*m*may vary in

*z*. Because the background is flowing at speed

*u*(

_{s}*z*), the particle wave frequency may be expressed aswhere

*ω*

_{intr}is the frequency relative to the moving flow (i.e., the intrinsic frequency). The work performed by the wave particle in the direction of the forcing along the ray path is given byThe integral (A4) may be thought of as the total work done by the wave particle in the direction of the forcing

*F*as it travels along the ray path through the forcing region. The quantity TE

_{ray}is similar to −

*E*

_{dyn}, defined in (17). For convenience, let us assume that the wave perturbation is in phase with the forcing at the initial time and that the spatial structure of

*S*(

*x*,

*z*) is given by a constant value

*Ŝ*along the ray path. The work performed by the wave particle along the ray path (A4) may then be written aswhere

*Ŝψ̂*is a positive real number and the change in phase is given byGiven a dispersion relation defining

*ω*

_{intr}and a background flow

*u*, the integral TE

_{s}_{ray}may be computed by analytical or numerical means. If the background flow is constant, then the integral (A5) is maximized for wave particle having an intrinsic phase velocity

*c*equal in magnitude and opposite the background flow,

_{x}*c*= −

_{x}*u*, for which the phase is stationary along the ray path.

_{s}