22 MONTHLY WEATHER REVIEW VOLUME 138
A Diagnostic Study of the Intensity of Three Tropical Cyclones in the Australian Region. Part II: An Analytic Method for Determining the Time Variation of the Intensity of a Tropical Cyclone*
FRANCE LAJOIE AND KEVIN WALSH School of Earth Sciences, University of Melbourne, Parkville, Australia
(Manuscript received 20 November 2008, in final form 14 May 2009)
ABSTRACT
The observed features discussed in Part I of this paper, regarding the intensification and dissipation of Tropical Cyclone Kathy, have been integrated in a simple mathematical model that can produce a reliable 15– 30-h forecast of (i) the central surface pressure of a tropical cyclone, (ii) the sustained maximum surface wind and gust around the cyclone, (iii) the radial distribution of the sustained mean surface wind along different directions, and (iv) the time variation of the three intensity parameters previously mentioned. For three tropical cyclones in the Australian region that have some reliable ground truth data, the computed central surface pressure, the predicted maximum mean surface wind, and maximum gust were, respectively, within 63 hPa and 62ms21 of the observations. Since the model is only based on the circulation in the boundary layer and on the variation of the cloud structure in and around the cyclone, its accuracy strongly suggests that (i) the maximum wind is partly dependent on the large-scale environmental circulation within the boundary layer and partly on the size of the radius of maximum wind and (ii) that all factors that contribute one way or another to the intensity of a tropical cyclone act together to control the size of the eye radius and the central surface pressure.
1. Introduction intensity and radial wind profile. The model uses input parameters that are readily available and simple to use. The problem of forecasting tropical cyclone intensity The principles upon which the new model is based are and wind field remains an important issue in tropical further described here and then expanded in the en- meteorology. A review of recent work is contained in suing sections. Lajoie and Walsh (2010, hereafter Part I). Based on the Because different atmospheric processes are involved sequence of observed characteristic features during in different parts of a tropical cyclone, and also because intensification and dissipation of Tropical Cyclone the radial wind profiles in different regions have distinct Kathy in Part I, we have developed a mathematical physical features (Shea and Gray 1973), analysis of the model for forecasting all aspects of tropical cyclone circulation is considered in four distinct regions, as shown intensity. As discussed in Part I, it was suggested that in Fig. 1. Region 1 is the circular area of radius r , the eye the inflow of a band of moist, near-equatorial air into e radius. Region 2 is the annular area between r and r . the cyclone center was responsible for bursts of in- e m Region 3 is the area bounded between radii r and r , tensification of the storm. In addition, it was proposed m c the radius of the inner core of the cyclone or of the that the inflow angle of the band was related to the rate central dense overcast (CDO) region, which on average of intensification. The concepts are used here to de- is 18 latitude, as it was for Kathy. The two last regions velop an accurate numerical model of tropical cyclone contain the eyewall cloud, the maximum wind, most of the destructive winds, and very heavy precipitation. * Supplemental information related to this paper is available at the Region 4 is between rc and a radius of 48 latitude. Journals Online Web site: http://dx.doi.org/10.1175/2010MWR2876.s1. To start the analysis, an initial radial wind profile in region 4 a few hours before the start of intensification, or Corresponding author address: Kevin Walsh, School of Earth at t 5 0, is obtained from a streamline and isotach anal- Sciences, University of Melbourne, Parkville VIC, 3010 Australia. ysis. The initial radial surface pressure profile is then cal- E-mail: [email protected] culated to balance the wind profile. Because (i) the feeder
DOI: 10.1175/2009MWR2876.1
Ó 2010 American Meteorological Society Unauthenticated | Downloaded 09/29/21 04:37 AM UTC JANUARY 2010 L A J O I E A N D W A L S H 23
In region 2, the wind is assumed to decrease linearly
from Vrm to zero at the cyclone center and the cyclo- strophic wind equation is used to obtain pre, the surface pressure at re. Since surface pressure does not vary ap- preciably with radius in the eye of a cyclone (see rele-
vant barograms on 285–286 in Riehl (1954)), pre 5 pc, the central surface pressure.
Parameters rm and re are determined from the Trop- ical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI; Simpson et al. 1988; Lonfat et al. 2004), TRMM Precipitation Radar (PR) data (Iguchi et al. 2000), or digitized high-resolution IR satellite cloud data (see Lajoie 2007; Lajoie and Walsh 2008). Finally, by using a reduction factor to convert gradient-level mean wind to mean surface winds, and the mean asymmetric FIG. 1. The inner and outer regions of a tropical cyclone. The wind distribution around a tropical cyclone, the hori- former is subdivided into the following: region 1, the eye where 0 # zontal distribution of the sustained mean surface wind as r # re; region 2 where re , r , rm (rm being the radius of maximum well as the maximum gust around the cyclone can be wind); region 3 where rm # r # rc (rc being the external radius of the estimated. CDO); and region 4, the outer region where r , r # 48 latitude. c Section 2 of this paper details the mathematical for- mulation of the prediction model. Section 3 gives results band originates at about 440 km from the cyclone center, of hindcast predictions of intensity changes for Tropical (ii) only moist near-equatorial air spirals through the Cyclone Kathy using the model. Section 4 shows calcu- feeder band into the cyclone core and no cool dry trade air lations of the estimated radial wind profile in Kathy. participates in the cyclonic circulation during intensifica- Section 5 gives hindcast intensity predictions for two tion, and (iii) the convective cloud suppression around the other cyclones, section 6 discusses issues associated with cyclone during the period of intense convection in the the intensification rate of cyclones, and section 7 pro- CDO is observed to extend to about that distance, it is vides concluding remarks. assumed that all physical processes that influence the cy- clone intensity take place within 440 km of the cyclone 2. Model formulation center and that only the cyclonic circulation starting from 440 km to the north and northeast of the cyclone center a. Region 4 (rc # r # r4) needs to be taken into consideration to determine the 1) CIRCULATION ANALYSIS IN REGION 4 change in its intensity or in its intensification rate. An energy equation is then applied to the circulation Streamlines around stationary axisymmetric tropical in region 4. It is shown that because of the sudden cyclones are logarithmic spirals (see Abdullah 1966; change of a4, the angle of inflow at a radius of 48 latitude, Lahiri 1981; Wong et al. 2007). According to Senn and to the north and east of the cyclone at t 5 0, gradient- Hiser (1957, 1959), the trajectories of air parcels are level winds subsequently increase at all radii in region 4 equiangular logarithmic spirals far from the cyclone (48 latitude . r $ 18 latitude). On assuming that air center, while close to the center the angle of inflow de- trajectories are logarithmic spirals, the distances trav- creases with decreasing radius. A similar assumption is eled by air parcels between selected radii can be de- used here. In region 4, the angle of inflow a (see Fig. 2), termined, and from the radial wind profile, Vrc the wind in the northerly airstream to the north and northeast of at r 5 rc or 18 latitude and the time at which air parcels the cyclone, which is assumed to be constant between leaving r 5 48 latitude reaches rc can also be determined. radii 48 and 28 latitude and varies as shown in Fig. 2 for In region 3, a theoretically determined gradient-level smaller radii. The maximum angle of inflow is 208 before wind profile, particularly developed for that region by intensification, 42.58 and later 338 during intensification, Riehl (1954, 1963) and empirically documented by Shea and 208 during dissipation as observed in Fig. 5 of Part I. and Gray (1973), is then used to determine the gradient- The angle of inflow a is assumed to decrease linearly to level winds at different selected radii and the maximum 208 for 28 latitude $ r $ 18 latitude and to decrease to 108 wind Vrm at rm. The gradient wind equation is then used for 18 latitude . r $ 08 latitude. The last assumption is to obtain the surface pressure gradient and the surface made to simplify the computations of pc without greatly pressure at selected radii and at rm. affecting the results.
Unauthenticated | Downloaded 09/29/21 04:37 AM UTC 24 MONTHLY WEATHER REVIEW VOLUME 138
The solid curves ending with an arrow represent the streamlines. AB is an arc of a circle of radius r. The quantity P is the location of an air parcel moving along the s direction with a speed V. Vector n is normal to the s direction, u is the angle POX, and O is the cyclone center. Following Malkus (1958, 1962) the equations of mo- tion of an air parcel at P in natural coordinates, assuming large-scale two-dimensional flow, can be written as dV ›V ›V w›V 1 ›p 5 1 V 1 5 À 1 F, dt ›t ›s ›z r ›s (1) V2 À1 ›p 1 ›p 1 fV5 5 cosa , (2) R r ›n r ›r r
FIG. 2. Radial variation of the angle of inflow a at different times. where V is the horizontal wind at the gradient level along the s direction, t is time, F is the frictional accel- It is further assumed, following Malkus et al. (1961), eration, R is the radius of curvature of the air trajecto- Malkus (1962), Gentry (1984), and also Lucas et al. (1994), ries, and f is the Coriolis parameter. The following that 1%–5% of air parcels in region 4 will rise in convective assumptions are used in Eq. (1): updrafts, and the rest will move on a quasi-horizontal tra- jectory to reach rc, the outside edge of the CDO. 1) In the expansion of dV/dt in Eq. (1), the term ›V/›t is The natural and polar coordinate systems used in the the local rate of change of V. For ›V/›t to change, the analysis of the circulation in region 4 are shown in Fig. 3. surface pressure gradient in region 4 must change
FIG. 3. Natural and polar coordinate system used to analyze the flow of near-equatorial air toward the cyclone. The curves are logarithmic spirals representing air trajectories. Air parcels at point P move along the s direction. The n is at right angle to s, O is the cyclone center, u is the angle that OP makes with a fixed arbitrary axis, and a is the angle of inflow (angle between s and the tangent at P to the arc of a circle of radius OP).
Unauthenticated | Downloaded 09/29/21 04:37 AM UTC JANUARY 2010 L A J O I E A N D W A L S H 25
FIG. 4. The two logarithmic spirals T1 and T2 represent trajectories of near-equatorial air with different angles of inflow at r 5 48 latitude distance. The DS is the distance between selected radii along the trajectory. Here Vr and pr are the wind speed and surface pressure, respectively, at r. Subscript 1 or 2 refers to T1 or T2.
with time, but we have shown in Fig. 14 of Part I that The reason for this assumption is because (i) only the surface pressure at Centre Island (assumed to be moist near-equatorial air originating at 440 km to the the same at any radius r in region 4) was constant north and northeast of the cyclone spirals through for about 9–12 h before a sudden surface pressure fall feeder bands into the cyclone core (no cool dry trade that lasted about 3 h. Thus, ›V/›t 5 0 for the 9 h, while air participates in the cyclonic circulation during in- during the following 3 h it will be shown that the tensification) and (ii) the convective cloud suppres- change of V is brought about by the term V(›V/›s). sion around the cyclone during the period of intense 2) Since we shall be using Eq. (1) only in region 4 and convection in the CDO is observed to extend to only at the gradient level (i.e., at the height of max- about that distance. imum wind, where ›V/›z vanishes), neglecting of the 6) The surface pressure field and hence the surface pre- term w›V/›z is entirely justified. ssure gradient outside the CDO remain constant while 3) The influence of the ocean–atmosphere surface drag ‘‘near-equatorial air’’ moves from r 5 48 latitude to
is negligible at the gradient level. r 5 rc or 18 latitude and changes when it moves from rc 4) The wind direction is almost constant with height near to the eyewall. the gradient level and the frictional force that is mainly Equation (1) therefore reduces simply to due to the vertical eddy transport of momentum by turbulence and convective elements is mostly along dV ›V 1 ›p 5 V 5 À 1 F. (1a) the s direction and negligible along the n direction. dt ›s r ›s 5) The surface- and gradient-level wind and the surface pressure at r 5 48 latitude distance from the cyclone The symbols used in the analysis are illustrated in Fig. 4,
center remain constant during the evolution of the which shows two trajectories: T1 with a4 of 208 and T2 cyclone. with a4 of 42.58. The former is the trajectory of air
Unauthenticated | Downloaded 09/29/21 04:37 AM UTC 26 MONTHLY WEATHER REVIEW VOLUME 138 parcels when Kathy was in a steady nonintensifying state Friction at the gradient level is due to turbulent mixing and the latter when it had started to intensify. Suppose and thermal updrafts at that level and can be regarded as 2 that at time t 5 0 a parcel of air is at point Po1 in Fig. 4 on a drag force, which is proportional to (V 2 ub) , where ub radius ro, where the surface pressure is po1, and where is the horizontal component of the wind inside and at the the air parcel is moving along trajectory T1 with an initial base of the updrafts (Malkus 1952). Since on the one velocity Vo1. Suffix ‘‘o’’ denotes the initial position or hand, air that is entrained from below the gradient level initial values of p or V. Suppose also that in time dt,theair tends to conserve its horizontal momentum (Austin and parcel has moved along T1 through a distance Ds1 and has Houze 1973), and on the other, surface winds are lighter reached the point Pf1 on radius rf (5r 2 dr) where the and proportional to V (actually 0.8V, see section 2b), 2 2 surface pressure is pf1 and where its speed is now Vf1. then F is proportional to V and can be written F 5 kV , Suffix f denotes the final position or the final values of p where k is a constant. If we assume that in a small finite and V after dt. Suffix 1 denotes that the air parcel has been distance Ds1, V increases linearly, then V 5 Vo 1 ls and moving on trajectory T1. If the air parcel had been ð ð Ds1 Ds1 moving on trajectory T its final position after dt would be 2 2 Fds5 k(Vo 1 ls) ds, 0 0 Pf2 (see Fig. 4) and the surface pressure and wind speed would be pf2 and Vf2. Now if the air density in region 4 is where l 5 (Vf 2 Vo)/Ds1. The friction term then be- assumed to remain constant and Eq. (1a) is integrated comes (for T1) with respect to s between s 5 0ands 5Ds1,weobtain ð ð Ds Ds 1 kDs 1 1 1 1 2 2 2 2 Fds5 (Vo1 1 Vf1 1 Vf1Vo1), (Vf1 À Vo1) 5 (po1 À pf1) À Fds. (3) 0 3 2 r 0
Equation (3) is an energy equation similar to that derived and Eq. (3) becomes (for T1) by Haltiner and Martin (1957) except for the vertical motion term that is neglected here for reasons discussed 2 2 2 2k Vf1 À Vo1 5 (po1 À pf1) À Ds1 above. Physically, Eq. (3) states that the increase in r 3 kinetic energy of the air parcel is simply the difference 2 2 3 (Vo1 1 Vf1 1 Vo1Vf1). (4) between the work done by the pressure gradient force and the energy dissipated by friction. In other words, Similarly, air parcels moving from point Po along tra- when the system is in a steady-state, nonintensifying jectory T2 instead of along trajectory T1, would arrive at condition, that is when a4 5 208, the work done by the point Pf2 with a velocity Vf2 given by the pressure force is then just enough to balance the energy dissipated by friction and to maintain the ex- 2 2 2 2k Vf2 À Vo2 5 (po2 À pf2) À Ds2 isting wind field. It is worth noting that Eq. (3) is used in r 3 region 4 only. In regions 2 and 3, because of the many 2 2 3 (Vo2 1 Vf2 1 Vo2Vf2), (5) complex atmospheric processes that take place there and that affect the wind speed we will use the radial where Ds2 is the distance between points Po2 and Pf2 mean wind profiles that have been determined by other along T2. By substituting for k in Eqs. (4) and (5), we investigators. obtain the following: