## 1. Introduction

The track, structure, and intensity of a tropical cyclone during and after extratropical transition (ET) are highly dependent on the details of the interaction with the midlatitude flow. Before recurvature, the relative locations of a tropical cyclone and a preexisting midlatitude trough determine whether the tropical cyclone accelerates into the midlatitudes ahead of the trough, or becomes embedded in the equatorward flow behind the trough. If the relative locations are incorrectly forecast, 72-h track errors in excess of 1000 km may result (Beven 2000). Following recurvature, the interaction with the midlatitude flow can modify the flow upstream of the tropical cyclone (Bosart and Lackmann 1995), contribute to enhanced ridging downstream (Henderson et al. 1999; Riemer et al. 2008), and excite a Rossby wave train, resulting in downstream development (Harr and Dea 2009). Ensemble prediction systems often exhibit a large variability in the development downstream of an ET event (Harr et al. 2008), indicating a reduced predictability. This variability has been related to the interaction of the tropical cyclone and the upper-level jet stream during the ET event (Anwender et al. 2008).

In Scheck et al. (2011, hereafter Part I) we investigated the interaction of a tropical cyclone with a zonally oriented front in a barotropic nondivergent framework. The tropical cyclone excites a localized frontal wave that is characterized by strong ridging downstream of the cyclone. The circulation associated with the frontal wave has a strong influence on the vortex motion. The interaction is sensitive to the vortex structure, especially to the presence of anticyclonic relative vorticity around the cyclonic vortex core. For stronger jets the interaction is more marked, in the sense that the tropical cyclone approaches the jet more rapidly.

Resonant waves (i.e., waves propagating with a phase speed matching the zonal translation speed of the cyclone) are of particular importance. A continuous acceleration of the cyclone accompanied by a continuous amplification or damping of the frontal wave is possible only for wave modes close to resonance. In general, the spectrum of frontal waves excited by the cyclone is dominated by wave modes that are resonant in the initial phase of the interaction process. As a result, these modes have the strongest influence on the cyclone motion.

In contrast to the idealized straight jets studied in Part I, real jets can be strongly perturbed. Frontal waves characterized by a trough upstream and a ridge downstream of the cyclone are associated with a circulation that advects the cyclone northward, bringing it closer to the front. The inverse configuration with an upstream ridge and a downstream trough acts to advect the cyclone southward, away from the front. The largest effects on the cyclone track are to be expected for near-resonant waves that stay in phase with the cyclone for a long time and cause a continuous meridional acceleration of the latter. Near-resonant frontal waves should actually be a common case, as every stationary (orography) or slow-moving source of frontal deformation will excite standing or slow-moving waves.

Therefore we investigate in this second paper a more realistic scenario in which preexisting frontal waves modify the cyclone–front interaction. Furthermore, the impact of these waves on the predictability of cyclone tracks is discussed. The paper is organized as follows. A short description of the model setup is given in section 2. In section 3 we discuss the influence of preexisting frontal waves on the vortex–jet interaction and on the predictability. Finally, section 4 contains a summary of our results.

## 2. Model description

### a. Numerical methods and grid

In this study we use a nondivergent barotropic model on the *β* plane with second-order diffusion. We employ an explicit finite difference code to solve these equations numerically (see Part I for details).

In our calculations the numerical domain extends *L _{x}* = 27 000 km in the zonal direction and

*L*= 9000 km in the meridional direction. For the model runs of Part I typically 3000 × 1001 grid points were used (corresponding to a grid length of 9 km). To allow for the large number of simulations that are required to explore the parameter space, most of the model runs performed for this study are computed with a reduced resolution of 1500 × 501 zones. This results in some damping of the cyclones but does not alter the results discussed below qualitatively (as tested for several cases).

_{y}We use a channel geometry with periodic boundaries in the zonal direction and walls at *y* = ±*L _{y}*/2. The value of

*β*is 2 × 10

^{−11}m

^{−1}s

^{−1}, which corresponds to a latitude of approximately 30°.

### b. Front and cyclone representations

*η*=

*f*±

*η*

_{0}separated by an interface at

*y*=

*y*(

_{f}*x*). The Coriolis parameter is approximated as

*f*=

*f*

_{0}+

*βy*. For convenience, the constant

*f*

_{0}, which has no influence on the calculations, is set to zero. Thus we adopt an absolute vorticity distribution

*y*≡ 0 km), Eq. (1) gives rise to the zonal velocity profile

_{f}*u*

_{0}=

*u*

For the model runs in this study we set the jet parameters to *u*_{0} = 40 m s^{−1} and *η*_{0} = 4 × 10^{−5} s^{−1} and use the SUD vortex described in Part I, which reaches the maximum wind speed of 40 m s^{−1} at a radius of 100 km. The start position of the vortex is 1500 km south of the front. Thus the initial state is set up as in model run BS from Part I, except that frontal waves are already present at *t* = 0 days.

*k*= 2

_{m}*πm*/

*L*is the wavenumber,

_{x}*ϕ*the initial wave phase, and

_{m}^{i}*b*the initial wave amplitude of mode

_{m}^{i}*m*. In some cases several wave modes are excited.

*m*propagates with a phase speed

*λ*=

_{m}*L*/

_{x}*m*is the wavelength (see Part I).

*ν*=

_{m}*b*

_{m}^{i}η_{0}/2. Unless otherwise noted, we use a frontal displacement amplitude

*b*= 250 km. For wave mode

_{m}^{i}*m*= 5 this corresponds to an initial steering flow of 1.7 m s

^{−1}at the center of the cyclone.

As the initial zonal cyclone position is *X _{c}* = 0, a value of

*ϕ*= 0 corresponds to the configuration with a trough upstream and a ridge downstream of the cyclone, which causes a northward flow at the cyclone center. The values

_{m}^{i}*ϕ*=

_{m}^{i}*π*/2,

*π*, and 3

*π*/2 correspond to eastward, southward, and westward flows, respectively. Note that in the meridional direction the circulation is damped exponentially with an

*e*-folding scale of

*λ*/2

_{m}*π*.

The flow steering the cyclone contains contributions associated with the vorticity step [Eq. (2)], the preexisting frontal wave, the frontal waves excited by the cyclone in the course of the model run, and the anticyclones sheared from the vortex in the initial phase (see Part I).

## 3. Interaction with a perturbed front

The effects of preexisting frontal waves are demonstrated using several sets of model runs. Within a set of calculations the mode number *m* and the initial wave amplitude *b _{m}^{i}* are kept constant and the initial wave phase

*ϕ*is varied between 0 and 2

_{m}^{i}*π*. Each set consists of about 40–50 model runs. The naming convention for the model runs is BS

*m*-[

*ϕ*/

_{m}^{i}*π*] and a full set of calculations is denoted by BS

*m*-

*X*. As in Part I, the letters B and S indicate that the runs are performed on the

*β*plane and that the SUD vortex is used, respectively.

The mode *m* = 5 is closest to resonance at the initial time and thus of special interest for the setup considered here. Given its phase speed, *c*_{5} = 5.6 m s^{−1}, and the initial zonal translation velocity of the cyclone for a straight jet, *U _{c}*(

*t*= 0 days) = 2.5 m s

^{−1}, the time scale over which this mode can grow [

*τ*

_{5}≈ 10 days; see Part I, Eq. (23)] is rather long. Because of the preexisting frontal wave,

*U*(

_{c}*t*= 0 days) varies with longitude in the model runs of set BS5-

*X*, with

*ν*

_{5}≈ 1.7 m s

^{−1}.

### a. Sensitivity to the wave phase

As a first step we investigate calculation set BS5-*X*, in which the wavelength of the preexisting frontal wave, *λ*_{5} = *L _{x}*/5 = 5400 km, is close to the initial resonant wavelength. In Fig. 1 the evolution of the cyclone position is shown for the cases

*ϕ*

_{5}

^{i}= 0,

*π*, 1.43

*π*, and 1.44

*π*and for the run with an unperturbed front (run BS). Not unexpectedly, the cyclone in run BS5-0.0 is advected toward the front faster than in run BS and crosses the interface 4 days earlier (Fig. 1a). The cyclone in run BS5-1.0, which is initially advected southward by the wave circulation, reaches the front 7 days later than in run BS. However, BS5-1.0 is not the most extreme case in terms of the time needed to reach the front. In model run BS5-1.43 the cyclone does not cross the interface at all in the 20 days included in the calculation.

An unexpected sensitivity to the initial conditions is found when comparing runs BS5-1.43 and BS5-1.44, in which the preexisting frontal wave is shifted only *λ*_{5} × 0.01 = 54 km farther upstream. The start positions of the cyclones in these runs are close to the trough axis (*ϕ*_{5}^{i} = 1.5*π*). After about 5 days the meridional locations of the cyclones in the two runs diverge dramatically. While in run BS5-1.43 the cyclone moves southward, it is advected rapidly northward in run BS5-1.44 and crosses the front at about the same time as in the unperturbed case BS. The huge differences in the temporal evolution of the meridional cyclone coordinate lead also to vastly different zonal cyclone positions (Fig. 1b). The cyclones that approach the front are advected downstream with translation speeds of up to 30 m s^{−1}, whereas the cyclones that stay farther away from the front attain only rather low zonal translation speeds. The difference in the cyclone position of the runs BS5-1.43 and BS5-1.44 is initially very small but increases rapidly from *t* ≈ 7 days. For *t* > 9 days the distance grows by about 2000 km day^{−1}. Thus, a preexisting frontal wave can have a strong influence on the cyclone track and in some cases this influence can be very sensitive to the position of the cyclone with respect to the wave.

The sensitivity of the cyclone track to the initial conditions varies strongly with the initial wave phase *ϕ*_{5}^{i}. For phases that are associated with a northward advection of the cyclone, the sensitivity is rather low. In all these cases the cyclone crosses *y* = 0 km after 6–7 days (Fig. 2). In contrast, in cases where the cyclone is initially advected southward the sensitivity is much higher and cases can be found in which the cyclone stays south of the front for more than 20 days (Fig. 2). The most interesting feature, however, is that there are several critical phases at which the solution changes abruptly—most notably at *ϕ*_{5}^{i} = 1.44*π*, between the runs BS5-1.44 and BS5-1.43, but also around *ϕ*_{5}^{i} ≈ *π* Further model runs have shown that at these critical phases shifting the frontal wave by a distance smaller than the grid length (in this case Δ*x* = 18 km) is sufficient to switch between vastly different solutions, (e.g., BS5-1.43-like solutions in which the cyclone remains south of the front and solutions in which the cyclone reaches the front in less than 10 days, as in run BS5-1.44). This behavior can be interpreted as a bifurcation in the solution space.

*π*/2 <

*ϕ*

_{c,m}<

*π*/2 (around a trough–ridge transition) cyclones are advected northward by the wave circulation; otherwise they are advected southward. South of troughs (

*π*<

*ϕ*

_{c,m}< 2

*π*) cyclones are advected eastward relative to the basic-state flow, and south of ridges (0 <

*ϕ*

_{c,m}<

*π*) they are advected westward. The phase

*ϕ*

_{c,m}is indicated by the shading in Fig. 3, which displays the evolution of the zonal cyclone position for all runs of set BS5-

*X*.

The cyclones starting at phases *ϕ*_{5}^{i} < 0.5*π* or *ϕ*_{5}^{i} > 1.5*π* are initially advected northward, thereby increasing their zonal translation velocity sufficiently to stay in the northward advecting range (white background in Fig. 3) for several days and reach the front within 9 days. Cyclones starting at phases slightly smaller than *ϕ*_{5}^{i} = 1.5*π* [i.e., upstream of the trough center, where the preexisting frontal wave causes a southward directed flow (gray background) initially] are still able to overtake the trough and reach the northward advection zone. This is possible for two reasons. At *t* = 0 days the cyclones are already moving with a translation speed close to the phase speed of the wave (5.6 m s^{−1}) because the wave circulation is eastward near the trough center and thus *U*_{c}(*t* = 0 days) ≈ *u**Y _{c}*)+

*ν*

_{5}= 4.4 m s

^{−1}. The increase in zonal velocity required to overtake the ridge is provided by the anticyclones, which cause initially a northward advection (see Part I) and thus also a zonal acceleration of the cyclone.

However, at a critical phase *ϕ*_{crit} ≈ 1.44*π* the cyclone attains a zonal velocity matching the phase speed of the wave at a time when it is south of the center of the trough (i.e., on the trough axis). In this position the meridional velocity component at the cyclone center caused by the wave is zero. Thus, the cyclone remains at the same meridional position and its zonal translation speed remains constant and equal to the phase speed of the wave. In the frame of reference of the wave this is a stationary situation—cyclone and wave move with the same zonal velocity and their separation remains constant. In runs with *ϕ*_{5}^{i} = *ϕ*_{crit} the cyclone can stay close to this position for a rather long time. Very small changes in the initial wave phase determine whether the cyclone eventually moves toward the next trough–ridge transition downstream of the trough, where it is attracted by the front, or toward the ridge–trough transition upstream of the trough, where it is repelled from the front. The fact that very small changes in the initial state determine which of two vastly different solutions will occur gives rise to the bifurcation.

*V*=

_{c}*υ*(

_{f}*X*,

_{c}*Y*) = 0 and

_{c}*ϕ*= 3

_{c}*π*/2) where the background flow and the circulation of the preexisting frontal wave are parallel, and a point R south of the ridge center (

*ϕ*=

_{c}*π*/2), where they are antiparallel. For the set of calculations considered here Eq. (8) yields meridional coordinates

*Y*

_{cT}≈ −1400 km and

*Y*

_{cR}≈ −1100 km for points T and R, respectively. The bifurcation at

*ϕ*≈ 1.44

*π*occurs close to the trough center and is thus clearly related to point T. The bifurcation at

*ϕ*≈

*π*(see Fig. 2) is also located close to the trough center.

^{1}No unusual spreading of the tracks is visible for cyclones that attain positions south of the ridge center (dotted lines in Fig. 3), so there are no bifurcations that occur close to point R.

The fact that bifurcations occur at point T but not at point R is related to the different flows in the vicinity of the two points. Close to point T the solutions diverge. As sketched in Fig. 4, cyclones that are located slightly downstream of point T are attracted toward the front and accelerated in the zonal direction, whereas cyclones slightly upstream of point T are advected farther away from the front and slow down. Point T can be approached from the southeast because a cyclone located downstream and farther away from the front than point T moves more slowly than the point itself and experiences a northward directed flow component from the wave. Both effects will bring the cyclone closer to point T.

In contrast to point T, point R cannot be reached by cyclones from the south. To get closer to the front, the cyclone must be located upstream of point R. However, in this position the cyclone is moving more slowly than point R, so that by the time the cyclone has reached the meridional position of point R, the latter has already propagated downstream. If a cyclone is located downstream and south of point R, it will move southward and slow down until it is overtaken by point R. Subsequently this cyclone will move toward the front and also pass point R on the upstream side (see Fig. 4). Thus point R cannot be reached by a cyclone moving northward, but the flow around point R causes cyclones to revolve around it (in the frame of reference of the frontal wave).

In accordance with these considerations, cyclone tracks diverging near point T and cyclones tracks orbiting around point R are visible in Fig. 5 (top panel), which displays the cyclone tracks of simulation set BS5-*X* in the frame moving with the wave. These cyclone tracks are influenced not only by the preexisting frontal wave, but also by frontal waves excited by the cyclone and by the anticyclones. The locations of points R and T indicated in Fig. 5, however, are computed under the simplifying assumption that the influence of the anticyclones and of the additional frontal waves excited by the cyclone can be neglected. Thus it is assumed that the cyclone is advected like a tracer particle in a velocity field (*u _{p}*,

*υ*) that is a superposition of the stationary jet and the flow caused by the propagating preexisting frontal wave. The limits of this approximation are evident from Fig. 5 (bottom panel), which displays cyclone tracks together with contour lines of the streamfunction Ψ

_{p}*that gives rise to (*

_{p}*u*,

_{p}*υ*) in the frame of the wave. The points R and T are located at maxima and saddle points of Ψ

_{p}*, respectively. Overall the streamlines are similar to the cyclone tracks—they diverge near point T and orbit around point R, and the angle between the direction of the cyclone motion and the closest streamline is relatively small for most of the time and most of the tracks.*

_{p}However, clear differences between the tracks and the streamlines are visible. These differences are related to the anticyclones and the frontal waves excited by the cyclone. During the formation of the anticyclones in the first few days the cyclone is advected northward (see Part I). Because of this effect the tracks diverge at a point slightly upstream of the position marked with T, where the southward flow component caused by the wave compensates for the northward flow due to the anticyclones. At later times, and in particular when the cyclone comes close to the interface, the waves excited by the cyclone can no longer be neglected compared to the preexisting wave, and the tracks no longer follow the streamlines. Therefore, considering the cyclone as a tracer particle can explain some basic features of the tracks, but this is not useful to determine the cyclone tracks quantitatively. When the amplitude of the preexisting frontal wave is reduced, the cyclone-induced waves become important earlier and eventually the bifurcation should vanish. However, the bifurcation of tracks near point T is still clearly visible even for a preexisting wave amplitude of only 100 km.

Figure 5 illustrates also why a bifurcation occurs at *ϕ*_{5}^{i} ≈ *π* (see Fig. 2). The cyclone starting at *ϕ*_{5}^{i} = *π* circles around point R, approaches point T from the northwest, and is deflected toward the front. The cyclone starting at *ϕ*_{5}^{i} = 1.02*π* takes the same route initially but is then deflected to the south. For slightly higher values of *ϕ*_{5}^{i} the cyclone is advected northward before it can come closer to point T. Thus, a rather large area downstream of point R cannot be reached by any cyclone.

### b. Influence of perturbation wavelength

The wavelength of the *m* = 5 frontal wave mode used above is slightly smaller than the initial resonant wavelength and the phase speed is slightly faster than the initial zonal translation velocity of the cyclone. In the following, the influence of waves with larger or smaller wavelengths is discussed and cases are investigated in which several wave modes are excited. For this purpose, several additional sets of calculations are considered (see Table 1 for a list of all sets discussed in this section). In particular, we consider several calculation sets similar to BS5-*X*, but with frontal waves of the modes *m* = 4, 7, 10, and 15.

The wave amplitude *b _{m}^{i}* = 250 km is the same in all runs, which means that at a fixed distance from the front the modes with the largest wavelengths are associated with the highest wind speeds [see Eqs. (5) and (6)]. Therefore the

*m*= 4 wave has a strong impact on the cyclone track—the time required by the cyclone to approach the front varies between 8 and 14 days (Fig. 6). However, for all phases the solution changes smoothly with

*ϕ*and there are no bifurcations. This fundamental difference is related to the fact that

_{m}^{i}*m*= 5 waves propagate in the same direction as the fluid itself, whereas

*m*= 4 waves propagate in the upstream direction. The bifurcations for

*m*= 5 can arise because two kinds of solutions can develop from very similar initial conditions: one in which the cyclone overtakes the wave, and one in which the wave overtakes the cyclone. This situation can be encountered only if the cyclone can actually obtain a zonal translation speed matching the phase speed of the wave. As the zonal translation speed of the cyclone is close to the zonal speed of the background flow (unless the cyclone is close to the front), this requires that the phase speed is contained in the zonal speed profile of the jet. This is not the case for

*m*= 4(

*c*

_{4}< 0,but

*u*

For smaller wavelengths [i.e., higher values of *m* and faster (positive) phase speeds], *U _{c}* =

*c*is again possible. With decreasing wavelength the frontal wave circulation decays faster in the meridional direction [Eqs. (5) and (6)] and the meridional positions of points R and T move closer to the front. Thus the strongest effects of preexisting frontal waves with smaller wavelengths, including bifurcations, are expected to occur only after the cyclone has moved sufficiently close to the front. Indeed, for the

_{m}*m*≥ 7 cases the cyclone position in the first 7 days hardly varies for different phases

*ϕ*. During this period the cyclone motion is not yet influenced significantly by the preexisting frontal wave, but the cyclone circulation excites waves with the current resonant wavelength. This wave component causes a northward acceleration of the cyclone, making it considerably more difficult for the preexisting wave to advect the cyclone southward. Consequently, the range of phases for which the cyclone does not cross the front during the simulated 20 days becomes narrower. While for

_{m}^{i}*m*= 5 a considerable fraction

*ϕ*

_{5}

^{i}∈ [0.85

*π*, 1.43

*π*] causes the cyclone to reach the front only after 10 days or later (Fig. 2), this is the case for

*m*= 7 only when

*ϕ*

_{7}

^{i}∈ [0.6

*π*, 0.7

*π*]. For higher mode numbers there is more time for the cyclone to excite frontal waves before the bifurcation process takes place. Thus the influence of the preexisting wave becomes weaker compared to the other waves, and the impact of the bifurcation process on the cyclone track becomes less severe. At

*m*= 10 the bifurcation process is able to change the time required by the cyclone to reach the front by about 3 days (Fig. 7) but does not lead to such strongly diverging evolutions as in the case of runs BS5-1.43 and BS5-1.44 (Fig. 1).

For preexisting frontal waves with very small wavelengths and high phase speeds, bifurcations become impossible again because the wave is always faster than the cyclone. In our calculations with a jet speed *u*_{0} = 40 m s^{−1} the maximum zonal translation velocity obtained by the cyclone is about 30 m s^{−1}, which corresponds to *m* ≈ 15 for the *β*-plane cases discussed here (see Fig. 13 in Part I). For *m* = 10 bifurcations are still present whereas for *m* = 15 no bifurcations are encountered (Fig. 6). Thus, a rather narrow band of wavelengths smaller than the initial resonant wavelength has a strong influence on the sensitivity of the solution on the initial conditions. Assuming that the zonal translation velocity of the cyclone is in the range [0, *u*_{0}/2], bifurcations can only be caused by waves with a phase speed in the same range. The corresponding range of wavelengths can be computed using the relation *u*_{0} = *η*_{0}^{2}/2*β* (see Part I, section 3a) and the dispersion relation [Part I, Eq. (6)]. The maximum wavelength (i.e., the wavelength of the wave with phase speed zero) is thus *λ _{b}* =

*π*

^{−1}, respectively. The minimum wavelength is about

*λ*/2.

_{b}To examine whether bifurcations can still occur when preexisting waves containing several wave modes are present, several further sets of calculations are considered. In a first experiment, set BS5a-*X*, the preexisting frontal wave consists of an *m* = 4 and an *m* = 5 component with equal amplitudes (250 km). The phase of the *m* = 5 component is varied continuously, whereas for the *m* = 4 mode only some phases are investigated. We show only the case *ϕ*_{4}^{i} = 1.5*π*, which should be the least favorable for the occurrence of bifurcations, as the *m* = 4 wave causes the maximum northward advection for this phase. Yet bifurcations are present even in this case (Fig. 8). The bifurcation is suppressed only when the amplitude of the *m* = 5 mode is reduced to 100 km and *ϕ*_{4}^{i} = 1.5*π* (set BS5b-*X*, Fig. 8). If either the amplitude is increased or a more favorable phase is chosen, bifurcations are present.

Finally, in set BS5c-*X* four wave modes (*m* = 4, 5, 6, 8) are excited in the preexisting frontal wave, each with an amplitude of 250 km. In this case several bifurcations are present (Fig. 8). These results demonstrate that bifurcations due to waves that are in resonance with the cyclone at some time are not easily suppressed by the presence of other resonant or nonresonant waves.

### c. Influence of frontal wave–induced shear

In the previous sections we showed that the main features of the interaction between the cyclone and the frontal wave can be explained by treating the cyclone as a point vortex. A secondary effect was the modification of the cyclone structure by the background shear associated with the jet, leading to the formation of anticyclones. The tropopause evolution, however, is sensitive to the interaction with the upper-level anticyclone of a tropical cyclone that has a much larger horizontal scale than the cyclonic core. Thus it is important to consider the structural evolution of an anticyclone as a function of its position relative to the frontal wave. Because of its large horizontal scale, the anticyclone is more sensitive to the horizontal shear of the jet and the frontal wave than the cyclonic vortex considered previously.

We investigate this interaction with a set of calculations, BN5-*X*, for which the anticyclonic part of the ANTI vortex (see section 3b in Part I) is used. In this anticyclonic vortex (*υ*_{max} = −8 m s^{−1}, *r*_{max} = 500 km) the velocity gradient is of a size comparable to the background shear. Furthermore, for the *m* = 5 frontal wave with amplitude 250 km in BN5-*X* the shear of the frontal wave circulation is also of the same order as the background shear. Consequently, the vortex structure can be changed considerably because of the shear and the changes depend on the initial position of the vortex.

Two extreme cases are shown in Fig. 9. In run BN5-0.6 the anticyclone is initially located south of a ridge, in the vicinity of point R. In this region the shear of the frontal wave circulation partially compensates for the background shear. The shape of the anticyclone is basically preserved as it circles around point R (Fig. 9, top panel). In run BN5-1.5 the anticyclone starts south of a trough, close to point T. Because of the strong shear deformation caused by the jet profile and the wave circulation in this region, the vortex is deformed to an elongated shape and then disrupted into two filaments within several days. The two filaments are advected to the northeast and the southwest, respectively, and end up south of the two neighboring ridges (Fig. 9, bottom panel). In this model run the anticyclone cannot influence the frontal wave substantially, as its vorticity distribution soon extends over a zonal range comparable to the wavelength of the frontal wave. In contrast, the anticyclone in run BN5-0.6, which remains compact and stays in phase with the wave, can influence the wave more strongly and causes a phase shift of about 500 km in 10 days (Fig. 9, top panel).

## 4. Summary

We investigated the interaction of a tropical cyclone and a perturbed tropopause front using an idealized model. In a barotropic framework a vorticity step is used to model the front and gives rise to a jetlike velocity profile. To study the influence of preexisting troughs and ridges on the cyclone–front interaction, we performed model runs in which a sinusoidal perturbation was applied to the jet in the initial state. In these cases a high sensitivity of the cyclone track to the initial cyclone position is observed. A bifurcation of the cyclone tracks occurs close to a point on the trough axis where the zonal translation speed of the cyclone matches the phase speed of the wave and the meridional translation speed is zero. Arbitrarily small displacements from this location determine whether a cyclone is advected northward and accelerated in the zonal direction or is advected southward and slowed down. The difference in the cyclone position for these two cases increases by up to 2000 km day^{−1}. A second point on the ridge axis fulfilling the same conditions does not give rise to bifurcations, but the flow in its vicinity causes cyclones to revolve around it (in the frame moving with the wave). To a first approximation the cyclone behaves like a tracer particle advected with the flow due to the jet and the preexisting frontal wave. Two factors cause the cyclone motion to deviate from that of a tracer particle. First, the anticyclones sheared away from the outer region of the vortex modify the track in the first few days. Second, frontal waves excited by the cyclone influence the motion, particularly in later phases.

The strongest sensitivity to the initial cyclone position is found for perturbation wavelengths slightly smaller than the initial resonant wavelength. Waves with larger wavelengths still have a strong influence on the cyclone motion, but they cannot cause bifurcations. For smaller wavelengths the bifurcation occurs closer to the front and leads to less severe differences in the cyclone tracks. For very small wavelengths bifurcations become impossible again. Thus, bifurcations can occur only for a limited range of wavelengths of the preexisting wave (in our standard case between about 3000 and 6000 km). If several wave modes are excited, this leads in general to additional bifurcation points.

For weak, broad vortices with parameters typical for anticyclonic outflows of real tropical cyclones, a preexisting frontal wave can also have a strong influence on the vortex structure. The shear deformation of the background flow is stronger south of troughs and weaker south of ridges. Weak vortices that come close to the bifurcation point are disrupted under the influence of the strong shear deformation and their remnants are advected into the regions south of ridges.

Despite the strong idealizations, the model runs performed for this paper show a variety of processes that are important for the interaction of real cyclones with jets. The strongly idealized model has the advantage that the causes and effects of these processes are easier to analyze than in three-dimensional full-physics models. For instance, through the idealized approach we were able to identify bifurcations in the cyclone tracks and explain why and where they occur. These scenarios are interesting as test problems for sensitivity and predictability studies. In a future study we will investigate the behavior of singular vectors for these cases, in particular for cyclone trajectories that come close to a bifurcation point. Furthermore, the scenarios discussed in this paper will be used as test problems for adaptive grid refinement techniques, which are evaluated for the use in tropical cyclone–related problems.

## Acknowledgments

This study was carried out as part of the German Research Foundation (DFG) Priority Program SPP 1276 “MetStröm: Multiple Scales in Fluid Mechanics and Meteorology.” We acknowledge the support of Tobias Hahn, Björn Rocker, and Vincent Heuveline of the Institute of Applied and Numerical Mathematics, Karlsruhe Institute of Technology, for assistance in implementing the model on graphics hardware. This study benefitted from discussions with Huw Davies. We are grateful to Pat Harr, Chris Rozoff, and an anonymous reviewer for their helpful comments on a previous version of the manuscript.

## REFERENCES

Anwender, D., P. A. Harr, and S. C. Jones, 2008: Predictability associated with the downstream impacts of the extratropical transition of tropical cyclones: Case studies.

,*Mon. Wea. Rev.***136****,**3226–3247.Bell, G. I., 1990: Interaction between vortices and waves in a simple model of geophysical flow.

,*Phys. Fluids***2A****,**575–586. doi:10.1063/1.857757.Beven, J., 2000: Tropical cyclone report: Hurricane Alberto (2000). National Hurricane Center. [Available online at http://www.nhc.noaa.gov/2000alberto.html].

Bosart, L. F., and G. M. Lackmann, 1995: Postlandfall tropical cyclone reintensification in a weakly baroclinic environment: A case study of Hurricane David (September 1979).

,*Mon. Wea. Rev.***123****,**3268–3291.Harr, P. A., and J. M. Dea, 2009: Downstream development associated with the extratropical transition of tropical cyclones over the western North Pacific.

,*Mon. Wea. Rev.***137****,**1295–1319.Harr, P. A., D. Anwender, and S. C. Jones, 2008: Predictability associated with the downstream impacts of the extratropical transition of tropical cyclones: Methodology and a case study of Typhoon Nabi (2005).

,*Mon. Wea. Rev.***136****,**3205–3225.Henderson, J. M., G. M. Lackmann, and J. R. Gyakum, 1999: An analysis of Hurricane Opal’s forecast track errors using quasigeostrophic potential vorticity inversion.

,*Mon. Wea. Rev.***127****,**292–307.Riemer, M., S. C. Jones, and C. A. Davis, 2008: The impact of extratropical transition on the downstream flow: An idealized modelling study with a straight jet.

,*Quart. J. Roy. Meteor. Soc.***134****,**69–91. doi:10.1002/qj.189.Scheck, L., S. C. Jones, and M. Juckes, 2011: The resonant interaction of a tropical cyclone and a tropopause front in a barotropic model. Part I: Zonally oriented front.

,*J. Atmos. Sci.***68****,**405–419.Schwierz, C., S. Dirren, and H. C. Davies, 2004: Forced waves on a zonally aligned jet stream.

,*J. Atmos. Sci.***61****,**73–87.

Time evolution of the meridional cyclone position as a function of the initial wave phase for all model runs of calculation set BS5-*X*. The contour interval is 250 km; the thick line corresponds to *Y _{c}* = 0 km. Lines corresponding to negative values of

*Y*are dashed. The gray shaded area indicates the initial phase range for which the flow across the cyclone center is southward at

_{c}*t*= 0 days. The small vertical lines at the top of the panel indicate the phases for which model runs have been performed.

Citation: Journal of the Atmospheric Sciences 68, 3; 10.1175/2010JAS3483.1

Time evolution of the meridional cyclone position as a function of the initial wave phase for all model runs of calculation set BS5-*X*. The contour interval is 250 km; the thick line corresponds to *Y _{c}* = 0 km. Lines corresponding to negative values of

*Y*are dashed. The gray shaded area indicates the initial phase range for which the flow across the cyclone center is southward at

_{c}*t*= 0 days. The small vertical lines at the top of the panel indicate the phases for which model runs have been performed.

Citation: Journal of the Atmospheric Sciences 68, 3; 10.1175/2010JAS3483.1

Time evolution of the meridional cyclone position as a function of the initial wave phase for all model runs of calculation set BS5-*X*. The contour interval is 250 km; the thick line corresponds to *Y _{c}* = 0 km. Lines corresponding to negative values of

*Y*are dashed. The gray shaded area indicates the initial phase range for which the flow across the cyclone center is southward at

_{c}*t*= 0 days. The small vertical lines at the top of the panel indicate the phases for which model runs have been performed.

Citation: Journal of the Atmospheric Sciences 68, 3; 10.1175/2010JAS3483.1

(top) Frontal displacement for *t* = 0 days. (bottom) Time evolution of the zonal cyclone position *ϕ _{c}* as a function of the initial wave phase for all model runs of calculation set BS5-

*X*. Dashed lines mark trough axes; dotted lines indicate ridge axes. In gray shaded areas the flow across the cyclone center associated with the frontal wave is southward; in white areas it is northward. Thick lines are used for the cases

*ϕ*/

^{i}*π*= 1.0, 1.02, 1.43, and 1.44.

Citation: Journal of the Atmospheric Sciences 68, 3; 10.1175/2010JAS3483.1

(top) Frontal displacement for *t* = 0 days. (bottom) Time evolution of the zonal cyclone position *ϕ _{c}* as a function of the initial wave phase for all model runs of calculation set BS5-

*X*. Dashed lines mark trough axes; dotted lines indicate ridge axes. In gray shaded areas the flow across the cyclone center associated with the frontal wave is southward; in white areas it is northward. Thick lines are used for the cases

*ϕ*/

^{i}*π*= 1.0, 1.02, 1.43, and 1.44.

Citation: Journal of the Atmospheric Sciences 68, 3; 10.1175/2010JAS3483.1

(top) Frontal displacement for *t* = 0 days. (bottom) Time evolution of the zonal cyclone position *ϕ _{c}* as a function of the initial wave phase for all model runs of calculation set BS5-

*X*. Dashed lines mark trough axes; dotted lines indicate ridge axes. In gray shaded areas the flow across the cyclone center associated with the frontal wave is southward; in white areas it is northward. Thick lines are used for the cases

*ϕ*/

^{i}*π*= 1.0, 1.02, 1.43, and 1.44.

Citation: Journal of the Atmospheric Sciences 68, 3; 10.1175/2010JAS3483.1

Schematic view of the flow steering the cyclone in the presence of a frontal wave. The thin solid arrows at the bottom indicate the direction of the flow associated with the frontal wave (*u _{f}* ,

*υ*). Lighter shading is used for a northward flow due to the frontal wave and darker shading for southward flow. The dashed arrows on the left depict the background jet flow

_{f}*u*

*c*; the dashed arrow shows the flow velocity of the unperturbed jet; and the thin solid arrow shows the flow velocity due to the wave. At these points the sum of the background flow velocity and the wave-induced flow velocity is equal to

*c*. Thin dotted arrows indicate cyclone tracks in the frame moving with the phase speed of the wave.

Citation: Journal of the Atmospheric Sciences 68, 3; 10.1175/2010JAS3483.1

Schematic view of the flow steering the cyclone in the presence of a frontal wave. The thin solid arrows at the bottom indicate the direction of the flow associated with the frontal wave (*u _{f}* ,

*υ*). Lighter shading is used for a northward flow due to the frontal wave and darker shading for southward flow. The dashed arrows on the left depict the background jet flow

_{f}*u*

*c*; the dashed arrow shows the flow velocity of the unperturbed jet; and the thin solid arrow shows the flow velocity due to the wave. At these points the sum of the background flow velocity and the wave-induced flow velocity is equal to

*c*. Thin dotted arrows indicate cyclone tracks in the frame moving with the phase speed of the wave.

Citation: Journal of the Atmospheric Sciences 68, 3; 10.1175/2010JAS3483.1

Schematic view of the flow steering the cyclone in the presence of a frontal wave. The thin solid arrows at the bottom indicate the direction of the flow associated with the frontal wave (*u _{f}* ,

*υ*). Lighter shading is used for a northward flow due to the frontal wave and darker shading for southward flow. The dashed arrows on the left depict the background jet flow

_{f}*u*

*c*; the dashed arrow shows the flow velocity of the unperturbed jet; and the thin solid arrow shows the flow velocity due to the wave. At these points the sum of the background flow velocity and the wave-induced flow velocity is equal to

*c*. Thin dotted arrows indicate cyclone tracks in the frame moving with the phase speed of the wave.

Citation: Journal of the Atmospheric Sciences 68, 3; 10.1175/2010JAS3483.1

(top) Cyclone tracks (solid lines) in the frame moving with the phase speed of the wave. For the cases *ϕ*_{5}^{i}/*π* = 1.0, 1.02, 1.43, and 1.44 thick lines are used and the start of each day is marked with a dot. The initial meridional position of the cyclones is *Y _{c}* = −1500 km. The dashed line is the preexisting frontal wave. The location of points R and T (see Fig. 4) is indicated. (bottom) The thick lines are as at top; the thin lines are contour lines of the wave streamfunction in the frame moving with the wave.

Citation: Journal of the Atmospheric Sciences 68, 3; 10.1175/2010JAS3483.1

(top) Cyclone tracks (solid lines) in the frame moving with the phase speed of the wave. For the cases *ϕ*_{5}^{i}/*π* = 1.0, 1.02, 1.43, and 1.44 thick lines are used and the start of each day is marked with a dot. The initial meridional position of the cyclones is *Y _{c}* = −1500 km. The dashed line is the preexisting frontal wave. The location of points R and T (see Fig. 4) is indicated. (bottom) The thick lines are as at top; the thin lines are contour lines of the wave streamfunction in the frame moving with the wave.

Citation: Journal of the Atmospheric Sciences 68, 3; 10.1175/2010JAS3483.1

(top) Cyclone tracks (solid lines) in the frame moving with the phase speed of the wave. For the cases *ϕ*_{5}^{i}/*π* = 1.0, 1.02, 1.43, and 1.44 thick lines are used and the start of each day is marked with a dot. The initial meridional position of the cyclones is *Y _{c}* = −1500 km. The dashed line is the preexisting frontal wave. The location of points R and T (see Fig. 4) is indicated. (bottom) The thick lines are as at top; the thin lines are contour lines of the wave streamfunction in the frame moving with the wave.

Citation: Journal of the Atmospheric Sciences 68, 3; 10.1175/2010JAS3483.1

Time required by the cyclone to reach *Y*_{c}(*t*) = 0km as a function of the initial wave phase for calculation sets BS4-*X* (dashed), BS5-*X* (thick solid), BS7-*X* (dashed–dotted), BS10-*X* (dotted), and BS15-*X* (thin solid). Here *Y*_{c}(*t*) is the meridional cyclone coordinate, averaged over 24 h. Nearly vertical lines indicate bifurcations and are absent in the cases *m* = 4 and *m* = 15.

Citation: Journal of the Atmospheric Sciences 68, 3; 10.1175/2010JAS3483.1

Time required by the cyclone to reach *Y*_{c}(*t*) = 0km as a function of the initial wave phase for calculation sets BS4-*X* (dashed), BS5-*X* (thick solid), BS7-*X* (dashed–dotted), BS10-*X* (dotted), and BS15-*X* (thin solid). Here *Y*_{c}(*t*) is the meridional cyclone coordinate, averaged over 24 h. Nearly vertical lines indicate bifurcations and are absent in the cases *m* = 4 and *m* = 15.

Citation: Journal of the Atmospheric Sciences 68, 3; 10.1175/2010JAS3483.1

Time required by the cyclone to reach *Y*_{c}(*t*) = 0km as a function of the initial wave phase for calculation sets BS4-*X* (dashed), BS5-*X* (thick solid), BS7-*X* (dashed–dotted), BS10-*X* (dotted), and BS15-*X* (thin solid). Here *Y*_{c}(*t*) is the meridional cyclone coordinate, averaged over 24 h. Nearly vertical lines indicate bifurcations and are absent in the cases *m* = 4 and *m* = 15.

Citation: Journal of the Atmospheric Sciences 68, 3; 10.1175/2010JAS3483.1

Evolution of the meridional cyclone position for model runs BS (dashed), BS10-1.67, and BS10-1.68.

Citation: Journal of the Atmospheric Sciences 68, 3; 10.1175/2010JAS3483.1

Evolution of the meridional cyclone position for model runs BS (dashed), BS10-1.67, and BS10-1.68.

Citation: Journal of the Atmospheric Sciences 68, 3; 10.1175/2010JAS3483.1

Evolution of the meridional cyclone position for model runs BS (dashed), BS10-1.67, and BS10-1.68.

Citation: Journal of the Atmospheric Sciences 68, 3; 10.1175/2010JAS3483.1

As in Fig. 6, but for calculation sets BS5a-*X* (solid), BS5b-*X* (dashed), and BS5c-*X* (dashed–dotted).

Citation: Journal of the Atmospheric Sciences 68, 3; 10.1175/2010JAS3483.1

As in Fig. 6, but for calculation sets BS5a-*X* (solid), BS5b-*X* (dashed), and BS5c-*X* (dashed–dotted).

Citation: Journal of the Atmospheric Sciences 68, 3; 10.1175/2010JAS3483.1

As in Fig. 6, but for calculation sets BS5a-*X* (solid), BS5b-*X* (dashed), and BS5c-*X* (dashed–dotted).

Citation: Journal of the Atmospheric Sciences 68, 3; 10.1175/2010JAS3483.1

Evolution of the *η* = 0 s^{−1} contour line for model runs (top) BN5-0.5 and (bottom) BN5-1.6.

Citation: Journal of the Atmospheric Sciences 68, 3; 10.1175/2010JAS3483.1

Evolution of the *η* = 0 s^{−1} contour line for model runs (top) BN5-0.5 and (bottom) BN5-1.6.

Citation: Journal of the Atmospheric Sciences 68, 3; 10.1175/2010JAS3483.1

Evolution of the *η* = 0 s^{−1} contour line for model runs (top) BN5-0.5 and (bottom) BN5-1.6.

Citation: Journal of the Atmospheric Sciences 68, 3; 10.1175/2010JAS3483.1

Parameters of the sets of runs with preexisting frontal waves. In all runs SUD vortices and *β*-plane jets with *u*_{0} = 40 m s^{−1} are used.