## 1. Introduction

This work has two objectives: 1) to propose an energy budget adapted from Lorenz’s (1955, hereafter L55) analysis in order to precisely quantify the growth of an African easterly wave (AEW) in a limited domain and 2) to apply this energy budget to the AEW that spawned Hurricane Helene (2006).

AEWs are wavelike disturbances observed in West Africa during the summer months. They have wavelengths of 3000–5000 km and periods of 3–5 days. Their amplitude is maximum at the level of the African easterly jet (AEJ) near 3000-m altitude (e.g., Erickson 1963; Burpee 1972). Norquist et al. (1977), following the energetic analysis in L55, showed that the growth of AEWs over West Africa is likely to be the consequence of a both barotropic and baroclinic instability of the AEJ. AEWs usually reach their maximal amplitude near the Guinean coast and start to decrease as they reach the tropical Atlantic Ocean, where baroclinic growth is reduced (e.g., Norquist et al. 1977; Thompson et al. 1979).

Several decades of observations have confirmed that a large proportion of Atlantic hurricanes evolve from AEWs (e.g., Avila and Clark 1989). However, few AEWs have a cyclogenetic evolution and the physical processes leading to such developments are not yet fully understood, although several hypotheses have been made. Landsea and Gray (1992) found a positive correlation between the Atlantic cyclonic activity and West Sahel rainfall variability during the summer months. This suggests that developing cases are the consequence of an interaction between AEWs and mesoscale convective systems (MCSs) that occur over West Africa during the monsoon period (we refer to the work of Laing et al. 2008 for an up-to-date description of West African MCSs during the monsoon season). Using statistics of AEWs in a 20-yr European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis dataset, Thorncroft and Hodges (2001) proposed that the cyclogenetic evolution of AEWs depends on their low-level amplitude when they leave the West African coast. Berry and Thorncroft (2005) studied the AEW associated with the pre-Alberto (2000) disturbance over West Africa. They found that the cyclonic vorticity structure of this AEW merged with cyclonic vortices of convective origin over the Guinea Highlands, thus resulting in its intensification. They finally hypothesized that this interaction process between the AEW and MCSs led to the cyclogenesis of Alberto (2000), one day after the disturbance left the West African coast.

The interaction between AEWs and convection has also been investigated quantitatively with the L55 energetic analysis in several idealized case studies (e.g., Thorncroft 1995) and realistic simulations (e.g., Hsieh and Cook 2007). These authors found that the latent heat released by convection increases the baroclinic growth of AEWs, but they did not separate the developing cases from the nondeveloping ones. Using composite statistics on a 60-yr National Centers for Environmental Prediction (NCEP) reanalysis dataset, Aviles (2004) applied the L55 energetic analysis to differentiate developing AEWs from nondeveloping ones. Surprisingly, she found that “the sources of eddy kinetic energy of developing AEWs are not as clearly barotropic and baroclinic as for nondeveloping ones” (136–137). Indeed, the cyclogenetic evolution of AEWs does not depend as much on their energetic growth over West Africa as on their more or less favorable internal structure when they leave the West African coast (“anticyclonic flow in the upper levels and surrounding the positive vorticity center, deep rising motion, weak and small cold core, a moist environment;” Aviles 2004, p. 140) and, to lesser extent, on a favorable environment in the Atlantic basin (e.g., Gray 1968).

Based on 15-yr ECMWF reanalyses, Kiladis et al. (2006) evaluated the adiabatic forcing of vertical motion with the divergence of **Q** vectors (Hoskins et al. 1978). They found “strong correspondence between the dynamically forced vertical motion at 850 hPa and deep convection for the day −2 through day +1 stages of the AEW evolution along 10°N” (2226–2227). In other words, vertical motions associated with AEWs in this region are mostly forced adiabatically, resulting in a simple relationship between convection and dynamics (convection occurs preferentially near the wave’s trough; see, e.g., Payne and McGarry 1977). However, the divergence of **Q** vectors “does a poor job in accounting for the more complicated vertical motion field for the 15°N, 17.5°W base point. This is in line with the likely important role of diabatic processes in the energetics of African easterly waves over some regions” (Kiladis et al. 2006, p. 2227). Kiladis et al. (2006) concluded that the relationship between AEWs and convection could be precisely analyzed with vorticity and energy budgets.

Here, we study the energetics of the West African disturbance that spawned Hurricane Helene about 1000 km east of Cape Verde Islands on 16 September 2006, during NASA’s African Monsoon Multidisciplinary Analysis (NAMMA) and AMMA’s third special observing period (SOP-3). Our goal is to quantify the respective contributions of large-scale and internal processes involved in this developing case study. This case study has been numerically simulated with the French nonhydrostatic, anelastic model Méso-NH (Lafore et al. 1998), whose main characteristics are given in section 2. The results of this simulation are analyzed in section 3 and compared to *Meteosat-9* images in the water vapor channel at 7.3 *μ*m. In section 4, we show the limits of L55’s analysis when applied to the energetics of a disturbance in a finite domain, and we propose an alternate formulation for hydrostatic, compressible, and anelastic frameworks in a limited domain. In section 5, a scale analysis applied to the simulated pre-Helene disturbance allows us to separate the most important terms in the proposed energy equations, and an energy budget including potential, zonal, and eddy kinetic energy is derived. This energy budget obtained for the pre-Helene simulated disturbance is discussed in section 6. Additional ideas for improving this energy budget are discussed in section 7. Section 8 gives conclusions and some perspectives on this work.

## 2. Simulation with Méso-NH

The pre-Helene disturbance was associated with a well-defined AEW and a growing MCS on 9 September 2006 northwest of Burkina Faso (Fig. 1). Successive convective developments associated with the AEW occurred during the two following days. On 12 September, the offshore convective redevelopment south of Cape Verde Islands spawned a tropical depression, which became Hurricane Helene four days later. Following the hypothesis of Berry and Thorncroft (2005) for the case study of Alberto (2000), the cyclogenesis of Helene (2006) could have resulted from the interaction between an AEW and MCSs.

To further investigate the processes involved in this cyclogenesis, a 4-day Méso-NH simulation, starting at 0000 UTC 9 September 2006, was conducted. During these four days, the pre-Helene disturbance propagated over about 3000 km in longitude. Because of computer limitations, it was not possible to run a simulation on this domain and for this duration at a resolved convection scale. Hence, we used a single model configuration with a horizontal resolution of 24 km and parameterized shallow and deep convection with the Kain–Fritsch–Bechtold scheme (Kain and Fritsch 1993; Bechtold et al. 2001). The simulated domain is 4800 km × 4800 km horizontally (Fig. 1). The vertical grid has 66 levels from the surface to 28 km with a grid spacing of 60 m near the surface up to 600 m at the tropopause level, and the orography is taken into account with the Gal-Chen and Sommerville (1975) vertical coordinate. This simulation is coupled with ECMWF operational analyses every 6 h at the boundaries of the simulated domain. Microphysics is parameterized with a one-moment mixed scheme and six classes of liquid and ice hydrometeors (Pinty and Jabouille 1998). Turbulence is parameterized with the 1D scheme with a 1.5-order closure proposed by Bougeault and Lacarrère (1989). Subgrid condensation is represented with the scheme proposed by Chaboureau and Bechtold (2005), and radiative processes with the radiation scheme used at ECMWF (Gregory et al. 2000).

The model outputs are saved every hour. The diagnosed data contain the dynamic and thermodynamic variables as well as the model-derived brightness temperature in the water vapor channel 6 (7.3 *μ*m) of *Meteosat-9* (Saunders et al. 2005; Chaboureau et al. 2000, 2002). The three-dimensional variables are interpolated from the Méso-NH original grid to constant altitude levels to compute the budgets of potential and kinetic energy.

## 3. Qualitative analysis

The Hovmöller diagram computed with *Meteosat-9* brightness temperature in the water vapor channel 6 (7.3 *μ*m) (Fig. 2a) shows that the pre-Helene disturbance was associated with three convective events: on 10 September over West Africa, on 11 September near the Guinean coast, and on 12 September off shore. The *Meteosat-9* images in Figs. 3a–g allow us to describe this convective activity more precisely. An MCS grew northwest of Burkina Faso on the afternoon of 9 September, propagated westward during the following night and morning, and then dissipated at 1200 UTC 10 September. A squall line developed in front of this system during the afternoon of 10 September, but it quickly dissipated over the Guinean Highlands. An MCS formed behind the previous squall line on 11 September morning, north of the Guinean Highlands, and dissipated over the nearby Atlantic during the afternoon. Finally, an oceanic MCS developed on the morning of 12 September south of Cape Verde Islands, which really started the cyclogenesis of Helene (2006).

The Hovmöller diagram with model-derived brightness temperature (Fig. 2b) shows comparable convective redevelopments, although not exactly at the same times and locations and with slightly warmer brightness temperatures. The model-derived brightness temperature images in Figs. 3h–n display this simulated convective activity more precisely. The MCS starting northwest of Burkina Faso on 9 September is simulated approximately at the right time and location. It spawned another MCS on 10 September afternoon, but this one did not have the shape of a squall line and it dissipated north of the Guinean Highlands during the night of 10 September. Another MCS redeveloped over Dakar (Senegal) on the afternoon of 11 September and intensified gradually. On the afternoon of 12 September, this simulated MCS showed evidence of tropical cyclogenesis (organized convection, cyclonically rotating cloud strips), but it was located about 200 km northeast of the actual MCS at this time. Hence, Méso-NH simulated a comparable but slightly different evolution compared to the actual pre-Helene disturbance. Despite these differences, we consider that a quantitative analysis of the energetic growth of this simulated developing system can provide valuable information on the processes occurring in the real atmosphere.

The Hovmöller diagram of the ECMWF analyzed relative vorticity at 700 hPa (Fig. 2c) shows the zonally moving AEW trough, associated with the previously described convective activity. We therefore suggest the synoptic AEW could have interacted with convection. The Hovmöller diagram of the model-derived relative vorticity at 3000 m (Fig. 2d) shows smaller-scale structures with high values embedded in the cyclonic vorticity zone of the AEW, which might be a consequence of such interactions. Horizontal cross sections of the simulated relative vorticity at 3000 m (Figs. 3o–u) reveal small-scale cyclonic vortices, probably of convective origin, that merged with the synoptic cyclonic vorticity structure of the AEW, thus intensifying and organizing it, especially near the West African coast. This is consistent with hypothesis of Berry and Thorncroft (2005) that interactions between AEWs and MCSs over the Guinean Highlands play an important role in tropical cyclogenesis off West Africa.

## 4. Energy budget

Burpee (1972) was the first to apply the L55 energetic analysis to study the growth of AEWs over West Africa. His work was based on 8-yr radiosonde data over the African continent. By evaluating the correlations [*υ***T**] and [*υ***u**] (see the list of symbols in the appendix), he found 1) a sensible heat transport directed down the temperature gradient at each level above the surface, which corresponds to a production of eddy available potential energy by zonal available potential energy, and 2) a horizontal momentum transport with largest magnitude at 700 hPa directed equatorward during the summer season, although the meridional shear of zonal momentum could not be evaluated and the sign of the barotropic conversion between zonal and eddy kinetic energy remained undetermined. Moreover, the correlation [*ω***T**] could not be evaluated and the baroclinic production of eddy kinetic energy was not quantified in this analysis.

A more complete energy budget was realized by Norquist et al. (1977). Their work was based on the data from the third Global Atmospheric Research Program (GARP) Atlantic Tropical Experiment (GATE III), which was held over West Africa and the eastern tropical Atlantic during the summer of 1974. Their energy budget was done over a limited domain but they neglected the boundary terms, assuming the processes responsible for the growth and/or decay of AEWs occur within their limited domain. They found that the growth of AEWs is both barotropic and baroclinic over West Africa but mostly barotropic over the ocean. This was confirmed by Thompson et al. (1979). Moreover, Norquist et al. (1977) supposed that latent heat released by convective systems over West Africa plays an important role in the baroclinic growth, without quantifying it. In an idealized case study of global circulation, Thorncroft (1995) also suggested that diabatic convective processes participate significantly in the baroclinic growth of AEWs. There was no boundary term in his budget because it was done at global scale. On the contrary, from an idealized case study in a limited domain with moist processes, Paradis et al. (1995) deduced that the baroclinic conversion is mainly compensated by a boundary term representing the redistribution of eddy kinetic energy by the pressure perturbation and that, consequently, it is not associated with production of eddy kinetic energy. Recently, Hsieh and Cook (2007) applied the L55 energetic analysis to a realistic simulation of AEWs during the summer season. They used a wavelet analysis to better understand their evolution. They found that the energy of AEWs over West Africa is mainly maintained by baroclinic conversion and to a smaller extent by barotropic conversion, and that baroclinic growth is sustained mainly by convective heating. Finally, the boundary terms in their budget have a minor contribution and the sink of AEWs’ energy is principally due to the frictional dissipation.

These analyses (except Paradis et al. 1995) that are based on the L55 energetic framework confirmed more or less the importance of convective processes in the growth of AEWs. The main problem with the L55 analysis is that it was originally proposed to study the maintenance of the general circulation and as such cannot strictly be applied to a limited domain. Thorncroft (1995), for example, avoids this problem by studying the energetics of a global circulation in an idealized case study. However, this cannot be done for a real case study. Consequently, there is a need to adapt the L55 analysis to a limited domain.

### a. Lorenz energetic analysis

*A*of the general circulation is proportional to the variance of potential temperature

*θ*[Eq. (8) in L55]:

*A*〉 is obtained [Eq. (16) in L55]:

*C*〉 is the baroclinic conversion term and 〈

*G*〉 is the diabatic source term. The horizontal kinetic energy 〈

*K*〉 of the general circulation is defined as [Eq. (11) in L55]

*K*〉:

*D*

_{friction}〉 is the energy dissipation by friction.

*ω*) in pressure coordinates is equal to zero when averaged for the general circulation. This leads to another prognostic equation for 〈

*K*〉 [Eq. (20) in L55]:

*C*〉 between 〈

*A*〉 and 〈

*K*〉 in the general circulation is proportional to the correlation 〈

*ω̃*

*T̃*〉. This means that at global scale the ascent of warm air (or descent of cold air) is associated with a decrease of 〈

*A*〉 and, conversely, an increase of 〈

*K*〉.

### b. Limits of the Lorenz energetic analysis

*E*〉 defined as

*E*is proportional to the enthalpy (

*H*=

*ρ*). In pressure coordinates, the thermodynamic equation is

_{a}c_{pa}T*E*〉:

*E*〉 and 〈

*K*〉 is proportional to the correlation 〈

*ω̃*

*T̃*〉, so it is equivalent to conversion between 〈

*A*〉 and 〈

*K*〉. This equivalence should also be valid at the local scale if

*A*is properly defined, which is not an easy task. We do not follow Marquet (1993, 1995, 2003a) and Shepherd (1993), who tried to define available potential energy as a local quantity. Here, we use the total potential energy

*E*or the enthalpy

*H*in relation with the kinetic energy, and the conversion term between

*E*and

*K*is referred to as the “horizontal baroclinic conversion.”

*ω*in pressure coordinates,

*ρ*in the definition of Φ in hydrostatic balance,

_{a}*E*〉:

*E*〉 and 〈

*K*〉 is proportional to the horizontal pressure work. It is interesting to note that the horizontal pressure work can only result from the ageostrophic component of the wind. When the circulation is in geostrophic balance and the wind is parallel to the isobars, there is no horizontal baroclinic conversion between

*E*and

*K*.

Marquet (2003b) called the horizontal pressure work the “ageostrophic conversion term” and compared it to the baroclinic conversion term of L55 for several case studies in a limited domain. He found that local changes in kinetic energy were directly related to the ageostrophic conversion term, whereas the baroclinic conversion term of L55 was “out of phase” and compensated by an additional term: the total (i.e., horizontal plus vertical) work of pressure forces.

*E*

*K*

*C*

*C*

*E*

*K*

### c. Energetic analysis including the vertical pressure work

The vertical pressure work CVP in (14) is also a source term in the equation of vertical kinetic energy, although it is not possible to write a budget of vertical kinetic energy in hydrostatic balance, as supposed in L55. Thus, we developed a more general energy budget taking into account the vertical kinetic energy, based on the nonhydrostatic, compressible equations given by Bannon (2002).

*ρ*instead of the dry

_{m}*ρ*, even though it is not a conservative quantity. This is indeed very important for the energetic analysis, given that we have seen in section 4b that the pressure work (linked to pressure gradients) is a measure of the rate of change of potential energy and that both dry air and water vapor actually contribute to the total pressure. From (18) we obtain the following prognostic equation for

_{a}*ρ*:

_{m}*ρ*is equal to the sum of its flux divergence plus a term quantifying the variation of mass of moist air due to changes in the water vapor content.

_{m}*K*= ½

^{C}*ρ*(

_{m}*u*

^{2}+

*υ*

^{2}) and the vertical kinetic energy KV

*= ½*

^{C}*ρ*

_{m}w^{2}:

*K*is equal to the sum of flux divergence BK

^{C}*and source terms: the horizontal pressure work CP*

^{C}*; the partial work of the Coriolis force CC*

^{C}*, which corresponds to a conversion between*

^{C}*K*and KV

^{C}*; the horizontal work of friction*

^{C}*K*due to variations of the water vapor content

^{C}*is equal to the sum of the flux divergence BKV*

^{C}*and the source terms: the vertical pressure work CVP*

^{C}*, the partial work of the Coriolis force CC*

^{C}*, the work of the gravitational force*

^{C}*due to variations of water vapor content*

^{C}*H*=

_{m}^{C}*ρ*), the thermodynamic equation of Bannon (2002) becomes

_{a}c_{pm}T*H*is equal to the sum of flux divergence BH

_{m}^{C}*and source terms: the pressure tendency*

^{C}*G*; the total pressure work (CP

_{p}^{C}*+ CPV*

^{C}*); the term associated with changes of water phase*

^{C}*H*due to variations of water vapor content

_{m}^{C}Equations (23), (24), and (25) show that *H _{m}^{C}* can be converted into (produced by) horizontal and/or vertical kinetic energy through the horizontal CP

*and vertical CVP*

^{C}*baroclinic conversions, the sum of these terms being equal to the total pressure work. This is a generalization of the result found in section 4b. When integrated over a spatial (finite or global) domain, (23), (24), and (25) constitute an energy budget including moist enthalpy and horizontal and vertical kinetic energies, schematically displayed in Fig. 4.*

^{C}### d. Energetic analysis including the anelastic approximation

*ρ*

_{aref}, is

*r*

_{hydro}=

*r*+

_{l}*r*) because they have been written with the moist air density, which includes hydrometeors. At first order, development of

_{i}*r*

_{hydro}allows us to separate the contribution of hydrometeors in the pressure terms. The terms (d4), (e4), and (f4) are the components of the Coriolis force; (d5), (e5), and (f5) the components of the frictional force; (d6), (e6), and (f6) the variation of the anelastic version of moist momentum due to variations of water vapor content; and (f7) is the gravitational force.

*K*= ½(1 +

^{A}*r*)

_{v}*ρ*

_{aref}(

*u*

^{2}+

*υ*

^{2}) and vertical kinetic energy KV

*= ½(1 +*

^{A}*r*)

_{υ}*ρ*

_{aref}

*w*

^{2}:

*π*

_{ref}is the Exner function of the reference state and

*c*

_{hydro}=

*r*, where

_{l}c_{l}+ r_{i}c_{i}*r*(

_{l }*r*) and

_{i}*c*(

_{l}*c*) are the mixing ratio and specific heat of liquid water (ice). In (31), term (a) is proportional to the diabatic heating rate, and term (b) corresponds to the moist correction in absence of any phase change. Using (25) and (31), the equation for the moist enthalpy can be written as

_{i}*p*

_{ref}is the pressure of the reference state. Although they are written differently, the same terms appear in (25) and (32). The additional term

*p*

_{ref}is used in the moist correction term. The fact that CP

*and CVP*

^{A}*appear in (29), (30), and (32) implies that the horizontal (vertical) baroclinic conversion between*

^{A}*H*and

_{m}^{A}*K*(KV

^{A}*) is correctly described with this set of anelastic equations. Hence, when integrated over the spatial (global or finite) domain, (29), (30), and (32) constitute an energetic budget encompassing anelastic versions of moist enthalpy and horizontal and vertical kinetic energies. For simplicity, the superscript*

^{A}*A*’s are omitted hereafter.

## 5. Application and scale analysis

The energy budget presented in section 4d is applied to the results of the Méso-NH simulation of the pre-Helene disturbance within domain D defined as 9°–18°N, 26°W–5°E horizontally (see Fig. 3) and 1000–18 000 m above mean sea level. Then, the terms BK, BKV, and BH of (30), (31), and (32) are boundary terms that represent the contribution of energy exchanges between atmospheric flows inside and outside domain D.

### a. Budget of moist enthalpy

The budget of *H _{m}* [Eq. (32); Figs. 5 and 6] is nearly balanced by seven “first-order” terms referred to as “primary processes” (<±0.1 W m

^{−3}in Fig. 6). The evolution of

*H*is characterized by a diurnal cycle with positive tendencies during the day and negative ones during the night, between the surface and 12 000 m. This diurnal cycle is mostly forced by radiative heating

_{m}*G*

_{diab}, with some influence of the semidiurnal variation of

*G*, in relation to the atmospheric tide (e.g., Dai and Wang 1999) and the coupling with the ECMWF analyses every 6 h. In domain D, the mean vertical velocity is always positive owing to convective and radiative processes with larger values during evening, as can be deduced from Fig. 6c showing positive values of CPV, proportional to (−

_{P}*w*∂

*p*/∂

*z*), with

*p*decreasing with height. The boundary term BH is dominated by the vertical flux divergence [−∂(

*Hw*)/∂

*z*], which is positive because

*w*is positive and

*H*decreases with height. This means that in this modeled case study, the production of moist enthalpy by diabatic processes, including radiation and latent heat release, is mostly dissipated by a direct vertical baroclinic conversion (CPV > 0).

The term *G*_{wv}, related to humidity changes, is mainly negative, which is coherent with water vapor globally condensing within the domain (*G*_{phase} > 0). The diurnal positive values of *G*_{wv} at low levels correspond to surface evaporation at the hottest time of the day. They are associated with negative values of *G*_{phase} and enhanced positive values of *G*_{diab} (see Figs. 6d,e,g).

The budget of *H _{m}* also involves two “second-order” terms referred to as “secondary processes” (<±10

^{−3}W m

^{−3}in Fig. 7). One of them is the generation of moist enthalpy by hydrometeors G

_{hydro}, which mostly results from the upward transport of nonprecipitating liquid and ice particles that carry sensible heat from lower levels and transmit it to the colder air above. It can be noted that the maximum of

*G*

_{hydro}during the night of 9–10 September is correlated with the maximum of

*G*

_{phase}, associated with strong convective activity. The other secondary process is the horizontal baroclinic conversion term CP with two periods with positive values (−CP < 0 in Fig. 7c) on 9 and late 12 September, during which energy is transferred from the moist enthalpy to the kinetic energy through a direct horizontal baroclinic conversion, and a period of negative values (−CP > 0 in Fig. 7c) during 10–12 September, during which energy is transferred from the kinetic energy to the moist enthalpy through an inverse horizontal baroclinic conversion. An analysis of the associated physical processes is given below.

The additional term *G*_{residual} (not shown) in the budget of *H _{m}* is negligible compared to the other ones presented above. Its amplitude is about 50 times smaller (<±2 10

^{−5}W m

^{−3}) than the secondary terms displayed in Fig. 7. This result gives us confidence in the accuracy of the anelastic approximation of Méso-NH in terms of energy conservation.

### b. Budget of horizontal kinetic energy

The budget of *K* [Eq. (29)] is approximately balanced at first order by three primary processes (<±10^{−3} W m^{−3} in Fig. 8). The horizontal kinetic energy increases on 9 September and partly on 10 September, decreases on the 11th and 12th, and then increases again during the evening of the 12th. This evolution is apparently controlled by the horizontal baroclinic conversion CP, with some influence from the negative friction term *D*_{friction} and from the positive advection term BK in the low to mid levels.

The first-order residual of the kinetic energy budget is almost balanced by the second-order conversion term *CC* between horizontal and vertical kinetic energy due to the work of the Coriolis force (<±0^{−5} W m^{−3} in Fig. 9). This term is proportional to −*uw* and, owing to the mean upward motions, it is negative (positive) in the low (mid to upper) levels because of the presence of southwesterly (easterly) winds. It can be noted that the higher extension of negative values on 10 and 11 September could be related to the upward transport of southwesterly momentum, in conjunction with enhanced convective activity diagnosed by *G*_{phase} in Fig. 6. The remaining terms—the dissipation of *K* by condensation/deposition of water vapor *D*_{wv} and the momentum forcing associated with hydrometeors *D*_{hydro}—are negligible.

### c. Budget of vertical kinetic energy

The budget of KV [Eq. (30)] is characterized by a very weak tendency with a near balance between two first-order terms (<±0.1 W m^{−3} in Fig. 10), referred to as primary processes—the vertical baroclinic conversion term CVP and the work of the gravitational force *D*_{gravity}. It is a well-known result that at first order, changes in potential energy are balanced by the work associated with the vertical displacement of the atmospheric column.

The residual of these first-order terms is mostly balanced by two secondary processes (<±10^{−3} W m^{−3} in Fig. 11): the vertical frictional dissipation DV_{friction} and the dissipation by vertical momentum forcing of moist air associated with hydrometeors DV_{hydro}. They have comparable contributions, with the second showing a maximum in relation with convective activity during the night of 9–10 September.

The budget of KV also involves “third-order” terms referred to as “tertiary processes” (<±10^{−5} W m^{−3} in Fig. 12): the conversion CC between *K* and KV due to the work of the Coriolis force, the boundary term BKV, and the tendency of KV. The vertically integrated contribution of BKV is negligible as a result of opposing values in the lower and upper troposphere, probably related to a net inflow of ascending air in the low levels and a net outflow of subsiding air in the upper levels. The tendency of KV is small, but not negligible locally. However, it varies rapidly from positive to negative values, and vice versa, so that KV remains close to zero. The vertical dissipation DV_{wv} associated with water vapor condensation and/or deposition of water vapor (not shown) is about two orders of magnitude smaller than the third-order terms in the budget of KV, so it can be neglected.

### d. Synthesis

Several processes at different orders of magnitude have been described in the energy budgets (see Table 1 for a classification of these processes with respect to their magnitude). In this study, we are concerned with the horizontal kinetic energy of the simulated pre-Helene disturbance, so we only need to consider the primary processes in the budget of *K*. At this first order of magnitude, the conversion term between *K* and KV can be neglected, so the budget of KV is disregarded (as well as term CC in the budget of *K*). Moreover, the vertical baroclinic conversion term CPV can be considered as a “dissipation” term in the budget of *H _{m}*, in relation with the vertical reorganization of the atmosphere through the work of the gravitational force DV

_{gravity}. However, CPV cannot be approximated by −DV

_{gravity}because CP and the sum CPV + DV

_{gravity}have comparable magnitudes. The introduction of water vapor in the boundary terms BK and BH resulted in additional dissipation terms

*D*

_{wv}and

*G*

_{wv}because of nonconservation of the water vapor mass in the atmosphere. At the considered order of magnitude,

*D*

_{wv}is negligible, so the water vapor contribution will not be considered in the expressions of

*K*, BK, CP, and

*D*. However,

*G*

_{wv}is a primary process in the budget of enthalpy

*H*and it cannot be neglected for the abovementioned reason. The contributions of water vapor in the other primary processes—tendency of

_{m}*H*, CPV,

_{m}*G*

_{phase},

*G*

_{diab}, and

*G*—have to be kept as well.

_{p}### e. Budget of eddy kinetic energy

The energetic cycle proposed in L55 was aimed at studying the maintenance of the global circulation, more particularly the transfers between kinetic energy of the zonal mean flow, kinetic energy of the eddy circulations, and available potential energy of the global atmosphere. Consequently, L55 separated the zonal and eddy components of *A* and *K*. The eddy components characterize the perturbations with respect to the zonal mean in the global atmosphere. In the L55 analysis, the conversion between zonal kinetic energy *K _{Z}* and eddy kinetic energy

*K*is referred to as “barotropic growth” and the conversion between eddy available potential energy

_{E}*A*and

_{E}*K*is referred to as “baroclinic growth.” Note that

_{E}*A*can be produced through a conversion of zonal available potential energy (

_{E}*A*) and/or through diabatic processes. The L55 energetic cycle is an attractive concept for studying wave disturbances but, as explained above, it cannot be directly applied to a limited domain.

_{Z}Applying the theoretical results of section 4, *A* is replaced by *H _{m}* and is not separated in zonal and eddy components, and the energy budget is done in a limited domain. Separating the eddy and zonal components of

*H*in the analysis would be meaningful on condition that an unambiguous relation exists between the eddy potential energy and the eddy kinetic energy, which is not straightforward. The conversion between

_{m}*K*(

_{E}*K*) and

_{Z}*H*is referred to as the “eddy (zonal) baroclinic conversion.”

_{m}*K*and

_{Z}*K*as

_{E}*H*used for the present energy budget is

_{m}*K*and

_{Z}*K*can be derived from (29) where the second-order terms and the water vapor contribution in the first-order terms are neglected, using the property that

_{E}*ρ*

_{aref}depends only on the vertical coordinate:

*G*

_{residual}and adding terms

*G*

_{phase},

*G*

_{diab},

*G*

_{hydro},

*G*, and

_{p}*G*

_{wv}into

*G*, (32) becomes

*, and BK*

_{Z}*are the boundary terms of the budgets of*

_{E}*H*,

_{m}*K*and

_{Z}*K*, respectively, in the limited domain. All the other terms characterize sinks/sources or conversions of energy within this domain. Herein,

_{E}*G*is the production/destruction of

*H*by diabatic processes; CPV is the vertical baroclinic conversion term between H

_{m}_{m}and the vertical kinetic energy KV, and is considered as a sink in the moist enthalpy budget; CPZ and CPE quantify the zonal and eddy horizontal baroclinic production/destruction of

*K*and

_{Z}*K*, respectively; CPE (CPZ) is the work of eddy (zonal) pressure force by the eddy (zonal) circulation; and DZ and DE are the frictional dissipation of

_{E}*K*and

_{Z}*K*. When the zonal (eddy) circulation is in geostrophic balance, there is no zonal (eddy) baroclinic conversion. Note that CK corresponds to the barotropic conversion between

_{E}*K*and

_{Z}*K*; it is the sum of four terms although it is dominated by the first one: term (1) in (46), related to the eddy meridional transport of eddy zonal momentum. Thus, CK is considered as the conversion between the energy of the zonal mean flow and

_{E}*K*.

_{E}## 6. Results

### a. Eddy vorticity

As seen in Fig. 3, the midlevel trough and ridge of an AEW (characterized by large-scale cyclonic and anticyclonic circulations, respectively) are found in domain D, although the pre-Helene disturbance is associated with the midlevel trough only. Hence, it is necessary to objectively determine the respective contributions of the cyclonic and anticyclonic circulations to the eddy perturbation with respect to the zonal mean in domain D.

*ξ*, is defined as the vertical vorticity of the eddy horizontal wind,

_{E}### b. Eddy geostrophic vorticity

*ξ*, which is equal to opposite of the horizontal divergence of the eddy pressure force normalized by the Coriolis parameter:

_{G}*ξ*and

_{E}*ξ*are equal when the atmospheric flow is geostrophic. The comparison between the mean profiles of

_{G}*ξ*and

_{E}*ξ*in Fig. 14 shows that the eddy circulation was approximately geostrophic, which was not certain for such a large perturbation at low latitudes (9°–18°N). On the afternoon of 12 September, however,

_{G}*ξ*was slightly weaker than

_{E}*ξ*between 1000 and 8000 m. This suggests that during this cyclogenesis event pressure decreased first, and then geostrophic adjustment of the wind field to the mass field led to an increase of

_{G}*ξ*.

_{E}### c. Eddy kinetic energy

The evolution of the vertical profile of *K _{E}* (Fig. 15) shows a maximum around 4000 m, which corresponds to the energetic signature of the simulated AEW. On the morning of 9 September,

*K*increased between 3000 and 5000 m and then remained approximately constant until the end of the day. On the morning of 10 September a strong increase of

_{E}*K*in the low levels was followed by a decrease of similar intensity during the afternoon. However, as seen in Fig. 14, at that time the main contribution to the eddy circulation came from the AEW ridge, which played no role in cyclogenesis. From the afternoon of 11 September,

_{E}*K*increased constantly between 1000 and 5000 m, especially at low levels, in conjunction with the increasing cyclonic vorticity (Fig. 14).

_{E}### d. Energy budget

Evolution of the vertically integrated terms involved in the budget of *K _{E}* is shown in Fig. 16; their vertical profiles are shown in Fig. 17. On 9 September, the first stage of

*K*growth was mostly due to eddy baroclinic conversion (CPE > 0), with a positive contribution of barotropic conversion (CK) below 4000 m altitude and a negative one from net advection BKE, except in the lowest levels and between 5000 and 6000 m. This situation is similar to that analyzed by Norquist et al. (1977), who used the inadequate L55 formulation for limited domains but also found that AEW growth over the West African continent is both barotropic and baroclinic. This also meets the idealized case study by Thorncroft (1995) and the space–time spectral analysis by Hsieh and Cook (2007), although these authors did not consider specifically the relationship between energetic conversions and cyclogenetic evolution off shore. Aviles (2004) found smaller baroclinic and barotropic conversions over West Africa for developing cases of AEWs compared to nondeveloping ones. However, the eddy baroclinic conversion in the proposed energy budget does not quantify the same processes as in L55, so it is difficult to compare the present result with previous ones. A new interpretation of this eddy baroclinic conversion term is proposed here.

_{E}In particular, it can be noticed that pressure and geostrophic vorticity increased while convective activity was relatively intense during the night of 9–10 September (cf. *G*_{phase} in Fig. 6d and *ξ _{G}* in Fig. 14b). Hence, the eddy circulation was mostly down the eddy pressure gradient in the low to midtroposphere (CPE > 0). This suggests that the wind field was adjusting to the pressure field through horizontal displacement of air. The eddy kinetic energy produced by this horizontal displacement of air through positive eddy pressure work (CPE > 0) was partly exported to the outer flow, as can be deduced from the negative net advection BKE.

The positive barotropic conversion (CK) in the low to midtroposphere from 10 to 12 September was a consequence of the increasingly cyclonic curvature of the easterly flow, in response to the pressure decrease associated with convective activity (Fig. 6d). The negative eddy baroclinic conversion (CPE) indicates that the work of the eddy pressure forces was then negative (i.e., the eddy air circulation was mostly oriented up the eddy pressure gradient). This resulted from the presence of the abovementioned high pressure zone (centered on 29°N, 27°W at 1200 UTC 11 September), which induced a strong northeasterly divergent flow in the western part of domain D (Figs. 3r,s).

The fact that baroclinic conversion became positive again on late 12 September could be the signature of strong convective activity and associated pressure decrease in an environment where the AEJ was weak and the northwestern anticyclonic circulation had moved farther west (at 1200 UTC 12 September, it was centered on 28°N, 32°W). Actually, barotropic conversion decreased on late 12 September because the flow was no longer zonal, but a closed cyclonic and convergent circulation formed instead (Fig. 3u). Thus, some of the kinetic energy created by the positive baroclinic conversion on late 12 September was restored to the outer flow by weak net advection. At that time the system became tropical depression Helene and was able to take enough energy from the sea to be energetically self-sustained (Emanuel 1986).

The observed pressure decrease in this developing system was probably of a convective origin. The next step to investigate this hypothesis would be to more precisely analyze the relation among convection, pressure changes, energy conversions, and the adjustments between mass and wind fields.

## 7. Perspective for the energetic analysis

*H*in the anelastic case (

_{m}^{A}*H*=

_{m}^{A}*ρ*

_{aref}

*c*). In Fig. 18, we show horizontal cross sections at 3000 m of the normalized quantities

_{pm}T*R*/

_{a}*C*(Figs. 18a–d) and

_{pa}H_{m}^{A}*R*/

_{a}*C*(Figs. 18e–h) so they can be compared with the total pressure (Figs. 18i–l). Pressure decreased in association with the intensification of pre-Helene disturbance, with minimum values of 711 hPa at 0000 UTC 11 September, 709 hPa at 0000 UTC 12 September, and 702 hPa at 0000 UTC 13 September. This is coherent with the increase of

_{pa}H_{m}^{C}*ξ*observed during this period (Fig. 14b).

_{G}The evolution of *H _{m}^{C}* (Figs. 18e–h) is comparable with the evolution of pressure because the moisture correction is small, and the budget of

*H*can be used to analyze the pressure decrease in the developing AEW. However, the evolution of

_{m}^{C}*H*(Figs. 18a–d), which is the adequate thermodynamic variable in Méso-NH, is related to the temperature only, since the density

_{m}^{A}*ρ*

_{aref}is constant at each level. This implies that the processes responsible for the variations of pressure cannot be thoroughly analyzed with the moist enthalpy equation in the anelastic approximation.

Assuming that the processes involved in the budget of *H _{m}^{C}* should have magnitudes comparable to those in the budget of

*H*, we suggest that in a numerical simulation using compressible equations, the same primary processes would be directly related to the variations of pressure. The residual term of these primary processes would still be linked to the eddy baroclinic conversion. Consequently, the budget of

_{m}^{A}*H*would allow us to quantify the interactions between the vertical circulation forced by convection, pressure decrease, energetic conversions, and adjustment between mass and wind fields.

_{m}^{C}## 8. Conclusions and perspectives

The pre-Helene disturbance was associated with an AEW and several MCSs over West Africa during the period before cyclogenesis over the tropical eastern Atlantic near the Cape Verde Islands. A 4-day numerical simulation with Méso-NH using convective parameterization does not exactly reproduce the observed series of these MCSs, but the simulated disturbance shows a relatively similar evolution. The development of this simulated AEW is analyzed with an energy budget adapted from the L55 analysis. Because the energetic analysis of L55 cannot be applied directly in a limited domain, we developed a concept of energy transformations in the atmosphere that can be used locally.

First, considering that there is no robust and versatile definition of available potential energy *A* for a limited domain, we use the total potential energy *E* or the enthalpy *H* to quantify the thermal atmospheric reservoir. We verified that the conversions between *E* (or *H*) and *K*, and between *A* and *K*, are equivalent at global scale—the original framework of L55.

Second, we found that in the same hydrostatic framework as L55, the so-called horizontal baroclinic conversion between enthalpy and horizontal kinetic energy can be expressed as the horizontal pressure work resulting from the ageostrophic component of the wind. This means that the horizontal pressure force is a measure of the amount of potential energy that can be converted into horizontal kinetic energy. This result was partly given by Marquet (2003b), although his “ageostrophic conversion term” did not appear explicitly in his budget of available potential energy. According to Marquet (2003b), there is a “potential-energy component which plays the role of vertical distribution of energy by transforming baroclinic into ageostrophic conversions, via unknown processes to be discovered” (p. 2492). We found that this ageostrophic conversion is actually the horizontal baroclinic conversion between enthalpy and horizontal kinetic energy, so there is no unknown source or sink in this energetic conversion. This definition of horizontal baroclinic conversion is different from that proposed by L55, but it is valid at global and local scales.

Third, we showed that in the most general case (e.g., using compressible equations) the potential energy of the atmosphere can be characterized by the moist enthalpy and that its conversion in total kinetic energy is still proportional to the total pressure work. A positive horizontal (vertical) pressure work corresponds to a growth of horizontal (vertical) kinetic energy. These energy conversions are referred to as horizontal and vertical baroclinic conversions, respectively.

Fourth, the energetic analysis has been adapted to the anelastic equations of Méso-NH. An application to the simulated pre-Helene disturbance confirms that the redefined baroclinic conversions are correctly described with this set of equations. At first order, the tendency of moist enthalpy results mostly from diabatic heating, net advection, and vertical baroclinic conversion. These so-called primary processes show a diurnal cycle similar to that of radiative heating. Convective heating also follows a diurnal cycle with a maximum during the evening. However, convection is also active during the day and the mean vertical velocity in the domain is always positive, although the maximum values occur during the evening. Therefore, the vertical baroclinic conversion is always positive with maximum values during the evening. Variations of pressure due to atmospheric tides and variations of water vapor content have a small but nonnegligible contribution compared to the previous primary processes, so they are considered as primary processes as well. The residual of these primary processes is balanced by the much smaller horizontal baroclinic conversion and total heating rate of hydrometeors. Since the primary processes in the budget of moist enthalpy are partly caused by the vertical circulation of convective origin, convection is associated with pressure decrease, and horizontal baroclinic conversion is related to the adjustment of the wind field to the pressure change, we suggest that there is a physical link between these processes. Other processes theoretically identified in the budget of *H _{m}* are negligible.

The tendency of horizontal kinetic energy *K* results mainly from three primary processes: horizontal baroclinic conversion, net advection, and frictional dissipation. The production/destruction of *K* is explained by the horizontal baroclinic conversion, with a negative frictional term and weak net advection. Other processes are negligible. In particular, the conversion term between *K* and the vertical kinetic energy KV, related to the work of the Coriolis force, is a secondary process, so the budget of KV is disregarded and the vertical baroclinic conversion term between *H _{m}* and KV is considered as a sink term in the budget of

*H*.

_{m}Fifth, this simplified energetic analysis is used to derive an energy budget similar to L55. Here, *K* is separated in the zonal mean *K _{Z}* and the eddy contribution

*K*. Likewise, the horizontal baroclinic conversion CP is separated in zonal CPZ and eddy CPE components to quantify the baroclinic growths of

_{E}*K*and

_{Z}*K*, respectively; CPE (CPZ) is proportional to the work of eddy (zonal) pressure forces by the eddy (zonal) circulation. When the zonal (eddy) circulation is in geostrophic balance there is no zonal (eddy) baroclinic conversion. The conversion between

_{E}*K*and

_{Z}*K*is classically referred to as the barotropic conversion. Unlike in the analysis in L55, however, the enthalpy

_{E}*H*is not separated into zonal and eddy components, so the energy budget proposed here involves only three energies:

_{m}*H*,

_{m}*K*, and

_{Z}*K*. This energy budget also includes boundary terms, diabatic production of

_{E}*H*, dissipation of

_{m}*H*by vertical baroclinic conversion, and frictional dissipations of

_{m}*K*and

_{Z}*K*. It is therefore adapted to quantify the respective contribution of large-scale and internal processes to the evolution of eddy kinetic energy of a disturbance such as an AEW in a limited domain.

_{E}Our main finding with the energy budget of *K _{E}* is that the intensification of the simulated pre-Helene disturbance during its continental stage resulted from both barotropic and baroclinic processes at low and mid levels. This positive baroclinic conversion was associated with a horizontal displacement of air down the eddy pressure gradient in response to a pressure decrease in the trough. The resulting positive eddy pressure work was partly balanced by a negative net advection of eddy kinetic energy. From the afternoon of 11 September over the near-tropical Atlantic, the intensification of the simulated pre-Helene disturbance was associated with a positive barotropic conversion between 1000 and 5000 m and a negative eddy baroclinic conversion between 2000 and 5000 m. The positive barotropic conversion resulted from the decreasing intensity of the AEW entering the considered domain through its eastern limit and the intensifying large-scale motion in the trough region. The negative baroclinic conversion was associated with a horizontal displacement of air up the eddy pressure gradient forced by an external process (i.e., the anticyclonic circulation to the northwest). Finally, on late 12 September, the midlevel tropospheric flow was no longer zonal so barotropic conversion weakened, positive baroclinic conversion started, and some eddy kinetic energy was restored to the outer flow. By then the observed developing AEW was already a tropical depression and was able to take enough energy from the sea to be energetically self-sustained (Emanuel 1986).

Sixth, we showed that the eddy circulation of the simulated pre-Helene disturbance was approximately geostrophic. We propose that the energetic conversions responsible for the observed development might have resulted from an adjustment of the wind field to the pressure decrease induced by convection. In the tropics, because the Coriolis force is not strong enough to lead to a rapid evolution to geostrophic equilibrium, the eddy circulation is therefore convergent and tends to fill up the depression. This is what probably occurred during the continental stage of simulated pre-Helene disturbance. Over the ocean, however, the pressure decrease induced by convection was accompanied by a positive barotropic conversion as the cyclonic curvature of the easterly flow intensified in the low pressure system. The eddy circulation was globally divergent with regard to the westward propagation of an anticyclonic circulation, centered northwest of the trough. In consequence, the wind field adjusted geostrophically to the pressure perturbation induced by convection, a relatively rare event in the tropics. This geostrophic adjustment of the wind following a pressure decrease of convective origin could be the main process responsible for this simulated cyclogenetic evolution off shore. However, in the anelastic case the budget of moist enthalpy cannot be explicitly related to pressure variations, so the link between energetic conversions and this geostrophic adjustment is not explicit, and this hypothesis cannot be further investigated.

At a smaller scale, Rogers and Fritsch (2001) among others have shown that successive convective developments in a mesoscale vortex can reduce the local Rossby radius, which facilitates geostrophic adjustment of the wind field to the mass field. It would be interesting to analyze the adjustment occurring in the pre-Helene disturbance in terms of interactions between the vertical circulation forced by convection, the pressure decrease, and other processes described by the energetic conversion terms. We suggest that a pressure budget deduced from the moist enthalpy Eq. (25) would allow us to quantify this. A numerical simulation resolving the compressible equations of Bannon (2002) can be used for that purpose. In the anelastic case, however, the link between variations of pressure and convection could be investigated with a diagnostic equation of pressure (an example of such a pressure diagnostic in the anelastic case is provided by Smith and Bannon 2008 for a two-dimensional idealized case).

These results will be compared with energy budgets for other developing and nondeveloping cases of AEWs over West Africa and the near-tropical eastern Atlantic. We will also simulate some stages of the pre-Helene and other West African disturbances with Méso-NH at higher resolution with explicit convection. The analysis of the vorticity production by MCSs and their interaction with the developing AEW, at the convective scale and mesoscale, will certainly complement the results obtained from energy budgets at synoptic scale.

## Acknowledgments

The present work is part of the first author’s thesis at Université Paul Sabatier Toulouse III, France. Based on a French initiative, AMMA was built by an international scientific group and is currently funded by a large number of agencies, especially from France, the United Kingdom, the United States, and Africa. It has been the beneficiary of a major financial contribution from the European Community’s Sixth Framework Research Programme. Detailed information on scientific coordination and funding is available on the AMMA International web site (see http://www.amma-international.org). Numerical simulations were conducted on the CNRS/IDRIS computer under Grants 070591 and 080591. We thank Didier Gazen and Juan Escobar for the technical support in Méso-NH, Dr. Jean-Pierre Chaboureau for the scientific support in Méso-NH, and Dr. Jean-Pierre Pinty for valuable discussions on cloud physics.

## REFERENCES

Avila, L. A., and G. B. Clark, 1989: Atlantic tropical systems of 1988.

,*Mon. Wea. Rev.***117****,**2260–2265.Aviles, L. B., 2004: African easterly waves: Evolution and relationship to Atlantic tropical cyclones. Ph.D. thesis, Dept. of Atmospheric Sciences, University of Illinois at Urbana–Champaign, 212 pp.

Bannon, P. R., 2002: Theoretical foundations for models of moist convection.

,*J. Atmos. Sci.***59****,**1967–1982.Bechtold, P., E. Bazile, F. Guichard, P. Mascart, and E. Richard, 2001: A mass-flux convection scheme for regional and global models.

,*Quart. J. Roy. Meteor. Soc.***127****,**869–886.Berry, G. J., and C. Thorncroft, 2005: Case study of an intense African easterly wave.

,*Mon. Wea. Rev.***133****,**752–766.Bougeault, P., and P. Lacarrère, 1989: Parameterization of orography-induced turbulence in a mesobeta-scale model.

,*Mon. Wea. Rev.***117****,**1872–1890.Burpee, R. W., 1972: The origin and structure of easterly waves in the lower troposphere of North Africa.

,*J. Atmos. Sci.***29****,**77–90.Chaboureau, J-P., and P. Bechtold, 2005: Statistical representation of clouds in a regional model and the impact on the diurnal cycle of convection during Tropical Convection, Cirrus and Nitrogen Oxides (TROCCINOX).

,*J. Geophys. Res.***110****,**D17103. doi:10.1029/2004JD005645.Chaboureau, J-P., J-P. Cammas, P. Mascart, J-P. Pinty, C. Claud, R. Roca, and J-J. Morcrette, 2000: Evaluation of a cloud system life-cycle simulated by Meso-NH during FASTEX using METEOSAT radiances and TOVS-3I cloud retrievals.

,*Quart. J. Roy. Meteor. Soc.***126****,**1735–1750.Chaboureau, J-P., J-P. Cammas, P. Mascart, J-P. Pinty, and J-P. Lafore, 2002: Mesoscale model cloud scheme assessment using satellite observations.

,*J. Geophys. Res.***107****,**4301. doi:10.1029/2001JD000714.Dai, A., and J. Wang, 1999: Diurnal and semidiurnal tides in global surface pressure fields.

,*J. Atmos. Sci.***56****,**3874–3891.Durran, D. R., 1989: Improving the anelastic approximation.

,*J. Atmos. Sci.***46****,**1453–1461.Emanuel, K. A., 1986: An air–sea interaction theory for tropical cyclones. Part I: Steady-state maintenance.

,*J. Atmos. Sci.***43****,**585–605.Erickson, C. O., 1963: An incipient hurricane near the West African coast.

,*Mon. Wea. Rev.***91****,**61–68.Gal-Chen, T., and R. C. J. Sommerville, 1975: On the use of a coordinate transformation for the solution of the Navier–Stokes equations.

,*J. Comput. Phys.***17****,**209–228.Gray, W. M., 1968: Global view of the origin of tropical disturbances and storms.

,*Mon. Wea. Rev.***96****,**669–700.Gregory, D., J-J. Morcrette, C. Jakob, A. M. Beljaars, and T. Stockdale, 2000: Revision of convection, radiation and cloud schemes in the ECMWF Integrated Forecasting System.

,*Quart. J. Roy. Meteor. Soc.***126****,**1685–1710.Hoskins, B. J., I. Draghici, and H. C. Davies, 1978: A new look at the

*ω*-equation.,*Quart. J. Roy. Meteor. Soc.***104****,**31–38.Hsieh, J-S., and K. H. Cook, 2007: A study of the energetics of African easterly waves using a regional climate model.

,*J. Atmos. Sci.***64****,**421–440.Kain, J. S., and J. M. Fritsch, 1993: Convective parameterization for mesoscale models: The Kain–Fritsch scheme.

*The Representation of Cumulus Convection in Numerical Models, Meteor. Monogr.,*No. 46, Amer. Meteor. Soc., 165–170.Kiladis, G. N., C. D. Thorncroft, and N. M. J. Hall, 2006: Three-dimensional structure and dynamics of African easterly waves. Part I: Observations.

,*J. Atmos. Sci.***63****,**2212–2230.Lafore, J. P., and Coauthors, 1998: The Méso-NH Atmospheric Simulation System. Part I: Adiabatic formulation and control simulations.

,*Ann. Geophys.***16****,**90–109.Laing, A. G., R. Carbone, V. Levizzani, and J. Tuttle, 2008: The propagation and diurnal cycles of deep convection in northern tropical Africa.

,*Quart. J. Roy. Meteor. Soc.***134****,**93–109.Landsea, C. W., and W. M. Gray, 1992: The strong association between western Sahelian monsoon rainfall and intense Atlantic hurricanes.

,*J. Climate***5****,**435–453.Lorenz, E. N., 1955: Available potential energy and the maintenance of the general circulation.

,*Tellus***7****,**157–167.Marquet, P., 1993: Exergy in meteorology: Definition and properties of moist available enthalpy.

,*Quart. J. Roy. Meteor. Soc.***119****,**567–590.Marquet, P., 1995: On the concept of pseudo-energy of T. G. Shepherd.

,*Quart. J. Roy. Meteor. Soc.***121****,**455–459.Marquet, P., 2003a: The available-enthalpy cycle. I: Introduction and basic equations.

,*Quart. J. Roy. Meteor. Soc.***129****,**2445–2466.Marquet, P., 2003b: The available-enthalpy cycle. II: Applications to idealized baroclinic waves.

,*Quart. J. Roy. Meteor. Soc.***129****,**2467–2494.Norquist, D. C., E. E. Recker, and R. J. Reed, 1977: The energetics of African wave disturbances as observed during phase III of GATE.

,*Mon. Wea. Rev.***105****,**334–342.Paradis, D., J-P. Lafore, and J-L. Redelsberger, 1995: African easterly waves and convection. Part I: Linear simulations.

,*J. Atmos. Sci.***52****,**1657–1679.Payne, S. W., and M. M. McGarry, 1977: The relationship of satellite inferred convective activity to easterly waves over West Africa and the adjacent ocean during phase III of GATE.

,*Mon. Wea. Rev.***105****,**413–420.Pinty, J. P., and P. Jabouille, 1998: A mixed-phase cloud parameterization for use in mesoscale non-hydrostatic model: Simulations of a squall line and of orographic precipitations. Preprints,

*Conf. on Cloud Physics,*Everett, WA, Amer. Meteor. Soc., 217–220.Rogers, R. F., and J. M. Fritsch, 2001: Surface cyclogenesis from convectively driven amplification of midlevel mesoscale convective vortices.

,*Mon. Wea. Rev.***129****,**605–637.Saunders, R., M. Matricardi, P. Brunel, S. English, P. Bauer, U. O’Keeffe, P. Francis, and P. Rayer, 2005: RTTOV-8 Science and validation report. Met Office Tech Rep. NWPSAF-MO-TV-007, 46 pp.

Shepherd, T. G., 1993: A unified theory of available potential energy.

,*Atmos.–Ocean***31****,**1–26.Smith, J. W., and P. R. Bannon, 2008: A comparison of compressible and anelastic models of deep dry convection.

,*Mon. Wea. Rev.***136****,**4555–4571.Thompson, R. M., S. W. Payne, E. E. Recker, and R. J. Reed, 1979: Structure and properties of the synoptic-scale wave disturbances in the intertropical convergence zone of the eastern Atlantic.

,*J. Atmos. Sci.***36****,**53–72.Thorncroft, C. D., 1995: An idealized study of African easterly waves. III: More realistic basic states.

,*Quart. J. Roy. Meteor. Soc.***121****,**1589–1614.Thorncroft, C. D., and K. Hodges, 2001: African easterly wave variability and its relationship to Atlantic tropical cyclone activity.

,*J. Climate***14****,**1166–1179.

## APPENDIX

### Mathematical Operators

### List of Symbols

Superscripts *C* and *A* refer to the compressible and anelastic cases. Subscripts *υ*, *Z*, *E* refer to the vertical component, zonal mean, and eddy component, respectively. Vectors are in bold.

Classification of the processes involved in the energy budget of section 4d (see Fig. 4) with respect to their magnitude. The middle column indicates the processes’ names and the right column gives their order of magnitude in W m^{−3}.