1. Introduction
Skillful prediction of typhoon (TY) intensity and sudden track change remains a challenge for numerical models (Carr and Elsberry 1995; Harr et al. 1996; Rappaport et al. 2009; DeMaria et al. 2014; Emanuel and Zhang 2016; Ge et al. 2018; Yan et al. 2019). The challenge is partly related to our limited understanding of the complicated interactions between a typhoon and nearby systems, such as easterly waves, monsoon throughs, and upper-level cold lows (UTCLs) or troughs (Molinari and Vollaro 1989; Lander 1994; Wang and Wu 2004; Davidson et al. 2008; Hendricks et al. 2011; Leroux et al. 2016; Peirano et al. 2016; Fischer et al. 2019). Among them, UTCLs are characterized by upper-level cyclonic potential vorticity (PV) anomalies and occur frequently over the tropical and subtropical western North Pacific (WNP) (Molinari and Vollaro 1989; Patla et al. 2009; Li et al. 2012) during the typhoon season (June–October) (Wei et al. 2016). As indicated by the name, a UTCL has a cold-core structure in the upper troposphere. The cold-core anomalies peak at 300–400 hPa, with the maximum wind speed around 250 hPa (Kelly and Mock 1982; Chen and Chou 1994; Chen et al. 2001; Campetella and Possia 2007). Although most pronounced in the upper troposphere, a UTCL may extend down to the lower troposphere or even to the surface. The mean lifespan of a UTCL is 6.3 days (Chen and Chou 1994). Some UTCLs develop in a tropical upper-tropospheric trough (TUTT; Sadler 1975) and are also known as TUTT cells, while some others are cutoff lows originated from the midlatitudes westerlies (Molinari and Vollaro 1989).
Patla et al. (2009) proposed a conceptual model on how the track and intensity of typhoons may be affected by a UTCL. The authors identified several factors that affect the typhoon–UTCL interaction, including the intensity, horizontal scale and depth of the UTCL, the typhoon intensity, and the separation distance between them. Using reanalysis data, Li et al. (2012) analyzed the sudden track change of Typhoon Meranti (2010). The authors showed that a UTCL modified the environmental vertical wind shear (VWS) and horizontal vorticity advection around the typhoon, leading to its sudden northward turning. Wen et al. (2019) studied the interaction between Typhoon Meranti and a UTCL using the European Centre for Medium-Range Weather Forecasts (ECMWF) ensemble forecasts. By comparing the success and failure groups, they found that the group with better track predications had a better representation of the UTCL structure. Additionally, past studies have also indicated that the typhoon’s intensity could be affected by upper-level troughs (Molinari et al. 1995; Bosart et al. 2000; Leroux et al. 2013; Chen et al. 2015). Komaromi and Doyle (2018) simulated the interaction between a typhoon and an upper-level trough in an idealized framework. Their results indicated that the equatorward outflow is stronger than the poleward outflow in the simulated typhoon in the absence of an upper-level trough. In contrast, the presence of an upper-level trough north of the typhoon modifies the inertial stability such that the outflow channel shifts poleward. The reduction of inertial stability may lead to enhanced upper-level divergence that promotes typhoon’s intensification. However, the interaction is sensitive to the relative distance between the typhoon and the trough, and the typhoon weakens due to the enhanced vertical wind shear when the trough is too close.
While previous studies have provided valuable insights on the typhoon–UTCL interaction, the impacts of UTCL on the intensity and track evolution of typhoons are mostly qualitative in the studies based on observational data and real-storm simulations. The objective of this study is to quantitatively assess the impacts of a UTCL on the track and intensity evolution of a typhoon using semi-idealized simulations. A notable difference between a trough and UTCL is that the propagation of a trough is largely determined by the midlatitude westerlies, while the motion of a cutoff low may be strongly affected by its interaction with a typhoon, which introduces an additional degree of freedom. We will focus on Typhoon Jongdari in 2018. The typhoon has an unusual counterclockwise track and rapid intensity changes likely tied to the interaction with a UTCL. We will employ the piecewise potential vorticity inversion (PPVI) method to separate the flow component associated with the UTCL (Davis and Emanuel 1991). The effects of UTCL on the track and intensity change of Typhoon Jongdari are quantified through the analysis of a series of numerical simulations.
The remaining sections of this paper are organized as follows. An overview of Typhoon Jongdari is presented in section 2. The numerical model, experiment designs, and PPVI method are described in section 3. The effects of a UTCL on the track and intensity of Typhoon Jongdari are analyzed in sections 4 and 5, respectively. The summary and discussion are given in section 6.
2. Overview of Typhoon Jongdari
Typhoon Jongdari is a strong, long-lived typhoon with an unusual track in late July and early August 2018. Typhoon Jongdari formed over the Northwest Pacific at 0500 UTC 25 July. From 26 to 30 July, Typhoon Jongdari experienced a rare counterclockwise turning southeast of Japan. Jongdari made landfall in Japan on 29 July, produced heavy precipitation, and then moved toward eastern China after a small looping track (Fig. 1a). The typhoon had a relatively large meridional translation speed before reaching Japan, with its center moved from 22° to 35°N in 3 days. It intensified rapidly on 26 July and reached its peak intensity at 1200 UTC 27 July, with the minimum central pressure of 950 hPa (Fig. 1b) and the maximum surface wind speed of 45 m s−1 (Fig. 1c) according to the JTWC best track data. The storm then weakened rapidly.


(a)The observed tracks of Typhoon Jongdari (black) and the UTCL (red) from 0000 UTC 26 Jul to 0000 UTC 1 Aug (Typhoon Jongdari and the UTCL are collocated at 1200 UTC 31 Jul). (b) Time series of the minimum central sea level pressure (CMSLP; units: hPa) of Typhoon Jongdari. (c) Time series of the maximum surface wind speed (Vmax; units: m s−1) of Typhoon Jongdari. (d) A longitude–pressure cross section of meridional wind (contours; units: m s−1) and temperature anomalies (shading; units: K) of the UTCL along 33°N at 0000 UTC 26 Jul. The temperature anomalies in (d) are defined with respect to the average over 120°–150°E along 33°N.
Citation: Monthly Weather Review 149, 5; 10.1175/MWR-D-20-0271.1

(a)The observed tracks of Typhoon Jongdari (black) and the UTCL (red) from 0000 UTC 26 Jul to 0000 UTC 1 Aug (Typhoon Jongdari and the UTCL are collocated at 1200 UTC 31 Jul). (b) Time series of the minimum central sea level pressure (CMSLP; units: hPa) of Typhoon Jongdari. (c) Time series of the maximum surface wind speed (Vmax; units: m s−1) of Typhoon Jongdari. (d) A longitude–pressure cross section of meridional wind (contours; units: m s−1) and temperature anomalies (shading; units: K) of the UTCL along 33°N at 0000 UTC 26 Jul. The temperature anomalies in (d) are defined with respect to the average over 120°–150°E along 33°N.
Citation: Monthly Weather Review 149, 5; 10.1175/MWR-D-20-0271.1
(a)The observed tracks of Typhoon Jongdari (black) and the UTCL (red) from 0000 UTC 26 Jul to 0000 UTC 1 Aug (Typhoon Jongdari and the UTCL are collocated at 1200 UTC 31 Jul). (b) Time series of the minimum central sea level pressure (CMSLP; units: hPa) of Typhoon Jongdari. (c) Time series of the maximum surface wind speed (Vmax; units: m s−1) of Typhoon Jongdari. (d) A longitude–pressure cross section of meridional wind (contours; units: m s−1) and temperature anomalies (shading; units: K) of the UTCL along 33°N at 0000 UTC 26 Jul. The temperature anomalies in (d) are defined with respect to the average over 120°–150°E along 33°N.
Citation: Monthly Weather Review 149, 5; 10.1175/MWR-D-20-0271.1
Due to the unusual counterclockwise track and rapid intensity changes, Typhoon Jongdari was a forecasting challenge for operational numerical models. Lei et al. (2020) showed that both the ECMWF and NCEP ensemble forecasts had low skill for the track and intensity forecast of Jongdari (2018) when the forecasts were initiated before the typhoon and UTCL started interacting, but the forecast skill initialized after the typhoon–UTCL interaction was improved. The contrast suggests that the interaction between the storm and UTCL is an important source of uncertainty. Close examinations of the National Centers for Environmental Prediction (NCEP) Global Forecast System (GFS) Final (FNL) analysis (1° × 1°) show that the UTCL moved southwestward and approached Typhoon Jongdari, and that these two systems were nearly collocated on 1200 UTC 31 July (Fig. 1a; the UTCL center is defined as the position of the maximum relative vorticity at 250 hPa). The UTCL had a maximum wind speed of about 20 m s−1 at 250 hPa, and the diameter of the outer closed circulation was about 2000 km at 0000 UTC 26 July (Fig. 1d). The maximum cold core occurred around 300 hPa, and the cold anomalies extended downward to the surface. These characteristics are consistent with those documented in previous studies for a UTCL (Sadler 1975; Chen and Chou 1994; Patla et al. 2009; Li et al. 2012). The question here is how much the UTCL affected Typhoon Jongdari’s track and intensity. A set of numerical experiments are conducted to investigate this question.
3. Methodology and experiment designs
a. Experiment designs
The Advanced Research version of the Weather Research and Forecasting Model (WRF-ARW), version 3.9.1 (Skamarock et al. 2008), is used in this study. Three two-way interactive domains are used, with horizontal resolutions of 27, 9, and 3 km, respectively. The outer domain covers 0°–57°N, 95°–174°E. The model has 55 sigma levels with more vertical levels in the boundary layer and near the tropopause, and the model top is set at 50 hPa. The initial and boundary conditions are obtained from 6-hourly NCEP GFS FNL analysis with 1° × 1° resolution. The Kain–Fritsch convective scheme (Kain and Fritsch 1993) is applied to the outer two domains, and turned off for the inner domain of 3 km, in which microphysics is used. In addition, the Rapid Radiation Transfer Model (RRTM) longwave radiation (Mlawer et al. 1997), Dudhia shortwave radiation (Dudhia 1989) schemes, YSU planetary boundary scheme (Hong et al. 2006), and the unified Noah LSM surface scheme (Tewari et al. 2004) are used.
A total of seven numerical experiments are conducted to examine the impacts of UTCL on the track and intensity evolution of Typhoon Jongdari (Table 1). In the control experiment (CTL), the UTCL and the TY are present in the GFS analysis. In the sensitivity experiment denoted as RCL, the UTCL circulation is removed in the initial and boundary fields using the PPVI method (described in section 3b). As such, the comparison between CTL and RCL runs helps reveal the influences of UTCL on the intensity and track of Typhoon Jongdari. The CTL and RCL runs both employ the WRF single-moment 6-class (WSM-6) microphysics scheme (Hong and Lim 2006). To test the sensitivity of the prediction to microphysics schemes, additional experiments were conducted with the Purdue Lin microphysics scheme (LIN) and the ETA Grid-scale Cloud and Precipitation microphysics scheme (Ferrier) for both the CTL and RCL runs (Lin et al. 1983; Ferrier et al. 2002). Both LIN and WSM6 schemes are 6-class single-moment schemes, containing the following species: water vapor, cloud water, cloud ice, snow, rain, and graupel. The WSM6 can be regarded as an updated LIN scheme with the reduced sensitivity to the model time step and is more suitable for high-resolution simulations. In contrast, the Ferrier scheme is a 4-class double-moment scheme and includes water vapor, cloud water, rain, and ice (cloud ice/snow/graupel). Additionally, another experiment, “Rland,” is carried out to examine the effects of terrain. In this experiment, the Japanese Archipelago is replaced with ocean of SST 300.7 K, which is the averaged SST over 10-km radius of the typhoon center at 60 h (right before landfall) and at 66 h (after the typhoon reemerges over the ocean). All the simulations were integrated from 0000 UTC 26 July to 0000 UTC 30 July.
Description of different experiments.



b. Piecewise potential vorticity inversion
PV thinking is a useful approach to study vortex evolution (Hoskins et al. 1985). PV is defined based on absolute vorticity and static stability (potential temperature gradient), and contains the information on both thermodynamic and dynamic states of the atmosphere. Based on the PV invertibility principle, the balanced wind and mass fields can be determined from the PV distribution given proper boundary conditions. In particular, piecewise potential vorticity inversion is a useful diagnostic technique (Davis 1992; Davis and Emanuel 1991). Using PPVI, the total PV field is divided into discrete portions corresponding to different dynamical features of interest, and the PV field for each portion is then inverted separately to obtain the balanced wind and mass fields associated with individual dynamical features. PPVI has been widely used to isolate the balanced circulations associated with individual entities in order to quantitatively examine their characteristics and interactions among them (Wu and Emanuel 1995a,b; Shapiro 1996, 1999; Wu et al. 2003, 2009; Yang et al. 2008). Here we will apply PPVI to isolate the circulation components associated with the UTCL and Typhoon Jongdari and assess how the UTCL affects the intensity and track evolution of Typhoon Jongdari.
The perturbation fields are defined as
To examine how well this method removes the balanced UTCL circulation, the total flow fields before and after the removal of UTCL at 0000 UTC 26 July are compared in Fig. 2. The oblique cross section of PV going through the centers of the typhoon and the UTCL (red line in Fig. 2d) shows that Typhoon Jongdari is characterized by cyclonic PV extending throughout the troposphere, with a maximum around 500 hPa (Fig. 2a). In contrast, the UTCL is mainly confined in the middle to upper troposphere, and its intensity increases nearly monotonically from 500 to 200 hPa. Corresponding to the upper-level PV feature, there are low geopotential anomalies, cold temperature anomalies and a cyclonic circulation at 400 hPa associated with the UTCL, while the circulation component of the typhoon is much weaker at the same level (Figs. 2b–d). In the initial conditions for the RCL experiment (Figs. 2e–h), the PV feature associated with the UTCL is largely gone, and the signals in the geopotential height, temperature and wind fields at 400 hPa are also substantially reduced. Overall, the PPVI technique does a reasonable job removing the UTCL from the initial conditions of the RCL run.


(a)–(d) Total fields and (e)–(h) fields after the removal of the UTCL using the PPVI method. (a),(e) The oblique cross sections of PV (units: 0.01 PVU) in two experiments that pass through both the typhoon and UTCL centers in the CTL run [as shown by the red lines in (d) and (h)]. The other columns show (b),(f) geopotential height (shaded; units: gpm), (c),(g) air temperature (shaded; units: K), and (d),(h) wind vectors at 400 hPa (units: m s−1).
Citation: Monthly Weather Review 149, 5; 10.1175/MWR-D-20-0271.1

(a)–(d) Total fields and (e)–(h) fields after the removal of the UTCL using the PPVI method. (a),(e) The oblique cross sections of PV (units: 0.01 PVU) in two experiments that pass through both the typhoon and UTCL centers in the CTL run [as shown by the red lines in (d) and (h)]. The other columns show (b),(f) geopotential height (shaded; units: gpm), (c),(g) air temperature (shaded; units: K), and (d),(h) wind vectors at 400 hPa (units: m s−1).
Citation: Monthly Weather Review 149, 5; 10.1175/MWR-D-20-0271.1
(a)–(d) Total fields and (e)–(h) fields after the removal of the UTCL using the PPVI method. (a),(e) The oblique cross sections of PV (units: 0.01 PVU) in two experiments that pass through both the typhoon and UTCL centers in the CTL run [as shown by the red lines in (d) and (h)]. The other columns show (b),(f) geopotential height (shaded; units: gpm), (c),(g) air temperature (shaded; units: K), and (d),(h) wind vectors at 400 hPa (units: m s−1).
Citation: Monthly Weather Review 149, 5; 10.1175/MWR-D-20-0271.1
4. Simulated results and track analysis
a. Overview of the simulated track and intensity evolution
The track and intensity evolution of Jongdari in the CTL and RCL runs are shown in Fig. 3. The simulated track of Typhoon Jongdari in CTL is nearly in line with the JTWC best track (Fig. 1a), with the storm taking a counterclockwise path and then making landfall in Japan on 29 July. The track error is less than 70 km before the landfall in Japan. The storm intensifies in the first 42 h and then weakens, which is also largely consistent with the best track (Fig. 1b) but with some quantitative differences. In contrast, when the UTCL is removed, the storm takes a nearly zonal path and tracks westward between 20° and 25°N at a slow speed in RCL. The simulated Typhoon Jongdari in RCL remains over the open ocean throughout the 3-day period, and intensifies steadily up to 72 h. The large differences between the CTL and RCL runs suggest that the typhoon–UTCL interaction is responsible for the unusual movement of Jongdari, as will be analyzed in the subsequent sections.


(a) TY tracks and (b) the time series of TY intensity in terms of CMSLP in the different runs from 0000 UTC 26 Jul to 0000 UTC 30 Jul. Shading in (a) shows the time-independent skin temperature field used in all simulations. The black box in (a) depicts the innermost domain, and the plotting domain roughly represents the intermediate domain. Different line styles and colors represent different simulations as shown in the legends.
Citation: Monthly Weather Review 149, 5; 10.1175/MWR-D-20-0271.1

(a) TY tracks and (b) the time series of TY intensity in terms of CMSLP in the different runs from 0000 UTC 26 Jul to 0000 UTC 30 Jul. Shading in (a) shows the time-independent skin temperature field used in all simulations. The black box in (a) depicts the innermost domain, and the plotting domain roughly represents the intermediate domain. Different line styles and colors represent different simulations as shown in the legends.
Citation: Monthly Weather Review 149, 5; 10.1175/MWR-D-20-0271.1
(a) TY tracks and (b) the time series of TY intensity in terms of CMSLP in the different runs from 0000 UTC 26 Jul to 0000 UTC 30 Jul. Shading in (a) shows the time-independent skin temperature field used in all simulations. The black box in (a) depicts the innermost domain, and the plotting domain roughly represents the intermediate domain. Different line styles and colors represent different simulations as shown in the legends.
Citation: Monthly Weather Review 149, 5; 10.1175/MWR-D-20-0271.1
Also shown in Fig. 3 are the track and intensity evolution of the experiments using different microphysics schemes (Table 1). The Lin scheme produces the strongest storm in both the control group and RCL group of experiments. This is consistent with Maw and Min (2017) finding that Lin scheme tends to produce a stronger tropical cyclone than the Ferrier or WSM6 scheme due to more heating between 300 and 100 hPa. Despite the sensitivity to microphysics, the storm intensity differences between the groups of simulations with and without a UTCL are robust. Interestingly, the track evolution is not sensitive to the choice of the microphysics scheme in this case. Previous studies (e.g., Fovell et al. 2009; Choudhury and Das 2017) suggested that different microphysics schemes may produce different outer wind structures for a typhoon and can potentially affect its track by modifying the beta gyres. However, the sensitivity may vary from case to case. In the case of Jongdari, the steering flow seems to play a dominant role in the storm motion and make the track insensitive to the microphysics options. For brevity, our subsequent analysis will focus on the first set of experiments using the WSM6 microphysics (i.e., CTL versus RCL run). We will examine the track differences between the CTL and RCL runs in section 4 and the intensity differences in section 5.
b. Potential vorticity budget analysis
In Eq. (5), the left-hand-side term is the local PV tendency (PVT); the first and second right-hand-side terms represent the horizontal PV advection (HAD), which includes the advection by the large-scale steering flow, the beta effect, and the nonlinear self-advection (Wu and Emanuel 1995a,b); the third term is the vertical PV advection (VAD); the fourth term is the diabatic heating term (DH); and the residual term F* represents friction and subgrid processes (including diffusion; not shown). We will focus on the wavenumber-1 structure because previous studies suggested that the wavenumber-1 PVT well reflects the typhoon’s motion (Holland 1983; Wu and Wang 2000; Chan and Williams 1987; Bi et al. 2015). When applying Eq. (5) to the wavenumber-1 pattern, the residual term also includes the interaction of the wavenumber-1 pattern with the other wavenumber components.
Figure 4 shows the PVT and different terms for the wavenumber-1 component averaged over 850–600 hPa in both experiments (CTL and RCL runs) at 24 and 48 h. The two snapshots are chosen to represent the typhoon–UTCL interaction at different distances and when the typhoon has different intensities. As shown in Fig. 4, the HAD term has a pattern similar to the PVT and makes the major contribution to the PVT in both simulations. Compared to the HAD term, the VAD and DH terms are more concentrated in the inner-core region of the typhoon. Specifically, HAD term in the CTL run shows a positive tendency northeast of the storm center at 24 h, causing the northeastward motion of the storm at the early stage. This is consistent with the advective impacts of the UTCL, which is located northeast of the typhoon. The DH term has a northwest–southeast-oriented dipole pattern with a negative center in the northwest. This DH pattern is largely consistent with the vertical velocity or convection pattern. Due to the impacts of the vertical wind shear, enhanced convection occurs on downshear left side, which is to the southeast of the storm center at 24 h as reflected by the reflectivity pattern (Fig. 5a). The VAD term is determined by both the vertical velocity distribution and the vertical PV gradient. DH partially cancel HAD and VAD within the radius of 150 km in the CTL run. In the RCL run, the PVT term shows strong positive values to the northwest of the storm center, consistent with the northwestward track of the storm. DH is consistent with the enhanced convection south of the storm center under vertical wind shear (Fig. 5c). The VAD and DH terms in the RCL run are weaker than those in the CTL run, consistent with the weaker storm intensity in RCL at the time.


Columns show (left to right) the total PV tendency (PVT), the horizontal advection (HAD), the vertical advection (VAD), and the diabatic heating (DH) terms (units: 0.01 PVU s−1). Rows show the (a) CTL and (b) RCL runs at t = 24 h and the (a) CTL and (b) RCL runs at t = 48 h.
Citation: Monthly Weather Review 149, 5; 10.1175/MWR-D-20-0271.1

Columns show (left to right) the total PV tendency (PVT), the horizontal advection (HAD), the vertical advection (VAD), and the diabatic heating (DH) terms (units: 0.01 PVU s−1). Rows show the (a) CTL and (b) RCL runs at t = 24 h and the (a) CTL and (b) RCL runs at t = 48 h.
Citation: Monthly Weather Review 149, 5; 10.1175/MWR-D-20-0271.1
Columns show (left to right) the total PV tendency (PVT), the horizontal advection (HAD), the vertical advection (VAD), and the diabatic heating (DH) terms (units: 0.01 PVU s−1). Rows show the (a) CTL and (b) RCL runs at t = 24 h and the (a) CTL and (b) RCL runs at t = 48 h.
Citation: Monthly Weather Review 149, 5; 10.1175/MWR-D-20-0271.1


Wavenumber-1 reflectivity (dBZ; shaded) in (a),(b) CTL and (c),(d) (b) RCL run averaged over 850–600 hPa. The black vector indicates the direction of environmental vertical wind shear averaged over 500 km of the typhoon center.
Citation: Monthly Weather Review 149, 5; 10.1175/MWR-D-20-0271.1

Wavenumber-1 reflectivity (dBZ; shaded) in (a),(b) CTL and (c),(d) (b) RCL run averaged over 850–600 hPa. The black vector indicates the direction of environmental vertical wind shear averaged over 500 km of the typhoon center.
Citation: Monthly Weather Review 149, 5; 10.1175/MWR-D-20-0271.1
Wavenumber-1 reflectivity (dBZ; shaded) in (a),(b) CTL and (c),(d) (b) RCL run averaged over 850–600 hPa. The black vector indicates the direction of environmental vertical wind shear averaged over 500 km of the typhoon center.
Citation: Monthly Weather Review 149, 5; 10.1175/MWR-D-20-0271.1
At 48 h, the PVT terms in the CTL and RCL indicate northwestward and westward propagation, respectively, consistent with their storm tracks. Compared to 24 h, the PVT terms at 48 h have become strong in both the CTL and RCL runs, which is mainly due to the stronger HAD terms. In addition, the DH and VAD terms in the CTL run are still stronger than those in the RCL run, consistent with the stronger typhoon intensity in the CTL run. The DH terms have a similar spatial pattern as the convection distribution under the impacts of vertical wind shear in both experiments (Fig. 5). The VAD term tends to cancel the DH term in the RCL run while the VAD and DH terms have a nearly 90° phase shift in the CTL run. Overall, the HAD terms have a similar spatial pattern and magnitude to the PVT terms in both runs, indicating the important role of steering flow in determining the typhoon track, which will be further examined in section 4c.
c. Steering flow


(top) The time–height profiles of the steering flow averaged within a 7° × 7° box following the typhoon center, which is defined as the location of minimum SLP. (bottom) The deep-layer steering flow components averaged over 800–200 hPa: “T” denotes the TY translation speed, “Vs” is the total steering flow, “Vcl” is the steering flow component induced by the UTCL, and “Vm” is the steering flow after removing the UTCL. (a) The CTL run and (b) the RCL run. The vectors represent the zonal and meridional flow components, with the vector scale shown in the lower right in the units of m s−1.
Citation: Monthly Weather Review 149, 5; 10.1175/MWR-D-20-0271.1

(top) The time–height profiles of the steering flow averaged within a 7° × 7° box following the typhoon center, which is defined as the location of minimum SLP. (bottom) The deep-layer steering flow components averaged over 800–200 hPa: “T” denotes the TY translation speed, “Vs” is the total steering flow, “Vcl” is the steering flow component induced by the UTCL, and “Vm” is the steering flow after removing the UTCL. (a) The CTL run and (b) the RCL run. The vectors represent the zonal and meridional flow components, with the vector scale shown in the lower right in the units of m s−1.
Citation: Monthly Weather Review 149, 5; 10.1175/MWR-D-20-0271.1
(top) The time–height profiles of the steering flow averaged within a 7° × 7° box following the typhoon center, which is defined as the location of minimum SLP. (bottom) The deep-layer steering flow components averaged over 800–200 hPa: “T” denotes the TY translation speed, “Vs” is the total steering flow, “Vcl” is the steering flow component induced by the UTCL, and “Vm” is the steering flow after removing the UTCL. (a) The CTL run and (b) the RCL run. The vectors represent the zonal and meridional flow components, with the vector scale shown in the lower right in the units of m s−1.
Citation: Monthly Weather Review 149, 5; 10.1175/MWR-D-20-0271.1
The deep-layer steering flow is also calculated for the RCL simulation. It also shows a good agreement with the typhoon motion. In particular, the steering flow does not show a strong poleward component after the removal of the UTCL. Although the environmental flow field from the CTL run with the UTCL winds removed (Vm) is close to the environmental flow field from the RCL run (Vs), the storm tracks in the two experiments are different and contribute to the differences between Vm from the CTL and Vs from the RCL, which are calculated along the storm tracks. The results suggest that the track differences between the two simulations are mainly due to the direct impacts of the UTCL on the steering flow, and the differences can be further amplified when a storm moves to a different location and is subject to different environmental flow.
d. Fujiwhara effect
The impacts of the UTCL on Jongdari’s track can be explained by the Fujiwhara effect (Fujiwhara 1921; 1923). The Fujiwhara effect has been widely used to interpret the interaction between two tropical cyclones. Additionally, Bi et al. (2015) suggested that the interaction between a typhoon and a monsoon gyre can be explained by the Fujiwhara effect as well. The interaction between the two involved cyclonic systems is sensitive to their sizes, depth, intensities, and the environmental conditions such as the vertical wind shear, moisture, and sea surface temperature (Chang 1983; Ritchie and Holland 1993; Wang and Holland 1995; Wu et al. 2003, 2009, 2010; Yang et al. 2008; Peng and Reynolds 2005, 2006). The interaction is also sensitive to the distance between two vortices. When the distance is less than a critical distance (i.e., ~1400 km), the two vortices rotate around the vorticity centroid (Lander and Holland 1993; Liou et al. 2016); when the relative distance is shorter (i.e., ~750 km), the two vortices attract each other (Brand 1970).
Figures 7a–e shows the time evolution of the UTCL-related horizontal flow averaged over 800–200 hPa and the typhoon and UTCL centers in the CTL run. The typhoon (UTCL) center is defined as the location of the maximum relative vorticity at 850 (250) hPa. The typhoon is southwest of the UTCL at 0000 UTC 26 July. Jongdari and the UTCL rotate counterclockwise around each other over the following few days, and the distance between the two systems decreases with time. Figure 7f shows the orbiting tracks of the two cyclonic vortices with respect to the midpoint between the two vorticity centers. Although the typhoon and UTCL have very different vertical structures, the orbiting motions of the two vortices suggest that the Fujiwhara effect is still largely valid: the two systems rotate around each other cyclonically while moving closer together from 26 to 30 July. The relative distance decreases from more than 1400 km on 26 July to about 600 km on 30 July. As they move closer to each other, the intensity of UTCL (represented by the 250-hPa PV) and the typhoon both decrease with time (Fig. 3b). We will examine the intensity evolution of Typhoon Jongdari in the next section.


Interaction between the UTCL and Typhoon Jongdari in CTL. (a)–(e) The 250-hPa PV (shading; units: 0.1 PVU) and averaged 800–200-hPa wind field (vectors; units: m s−1) associated with the UTCL. The blue star and the red TY symbol represent the UTCL and typhoon centers, respectively. (f) The relative positions of the UTCL and typhoon with respect to their vorticity centroid at different times with 6-h intervals.
Citation: Monthly Weather Review 149, 5; 10.1175/MWR-D-20-0271.1

Interaction between the UTCL and Typhoon Jongdari in CTL. (a)–(e) The 250-hPa PV (shading; units: 0.1 PVU) and averaged 800–200-hPa wind field (vectors; units: m s−1) associated with the UTCL. The blue star and the red TY symbol represent the UTCL and typhoon centers, respectively. (f) The relative positions of the UTCL and typhoon with respect to their vorticity centroid at different times with 6-h intervals.
Citation: Monthly Weather Review 149, 5; 10.1175/MWR-D-20-0271.1
Interaction between the UTCL and Typhoon Jongdari in CTL. (a)–(e) The 250-hPa PV (shading; units: 0.1 PVU) and averaged 800–200-hPa wind field (vectors; units: m s−1) associated with the UTCL. The blue star and the red TY symbol represent the UTCL and typhoon centers, respectively. (f) The relative positions of the UTCL and typhoon with respect to their vorticity centroid at different times with 6-h intervals.
Citation: Monthly Weather Review 149, 5; 10.1175/MWR-D-20-0271.1
5. Storm intensity and the impacts of the upper-level cold low
As described in section 4a, the simulated intensity of Typhoon Jongdari in the CTL run increases initially and then decreases, similar to that observed. The impacts of the ULCL on the typhoon intensity evolution are examined in this section.
a. Upper-level outflow
Figure 8 shows the 200-hPa horizontal flow in CTL and RCL averaged during 12–24 h, a time period when the storm intensity clearly diverges between the two simulations (Fig. 3b). A prominent feature in the CTL run is a strong southwesterly jet north of the typhoon center. The jet is located east of the UTCL. Although the strong wind does not extend above the typhoon center, it is connected to the outflow of the typhoon. Figure 8a also shows a region of reduced inertial stability north of the typhoon center and near the UTCL in the CTL run. In other words, the UTCL modifies the upper-level absolute vorticity and reduces the inertial stability north of the typhoon, which contributes to stronger poleward outflow (Fig. 8a). It helps enhance the secondary circulation of the typhoon and favors its development (Rodgers et al. 1991; Shi et al. 1997; Fischer et al. 2017). In contrast, Fig. 8b shows that the outflow in RCL is much weaker and is directed equatorward toward the region of lower inertial stability (Rappin et al. 2011; Barrett et al. 2016; Komaromi and Doyle 2018).


The 200-hPa flow (vectors; units: m s−1), the inertial stability (shading; units:10−5 s−1), and the jet (contours; wind speed > 20 m s−1 with an interval of 5 m s−1) averaged over 12–24 h. The center of the domain represents the typhoon center, and the tick marks show the distance from the typhoon center in km.
Citation: Monthly Weather Review 149, 5; 10.1175/MWR-D-20-0271.1

The 200-hPa flow (vectors; units: m s−1), the inertial stability (shading; units:10−5 s−1), and the jet (contours; wind speed > 20 m s−1 with an interval of 5 m s−1) averaged over 12–24 h. The center of the domain represents the typhoon center, and the tick marks show the distance from the typhoon center in km.
Citation: Monthly Weather Review 149, 5; 10.1175/MWR-D-20-0271.1
The 200-hPa flow (vectors; units: m s−1), the inertial stability (shading; units:10−5 s−1), and the jet (contours; wind speed > 20 m s−1 with an interval of 5 m s−1) averaged over 12–24 h. The center of the domain represents the typhoon center, and the tick marks show the distance from the typhoon center in km.
Citation: Monthly Weather Review 149, 5; 10.1175/MWR-D-20-0271.1
b. Eddy fluxes convergence
The left column in Fig. 9 shows the radius–time plots of EFC in the CTL and RCL runs. The EFC is vertically averaged between 100 and 300 hPa. As expected, the EFC shows marked differences between CTL and RCL runs (Figs. 9a,d). In the CTL run, the EFC maximum occurs around the radius of 1200 km during the first few hours of the model integration. It then propagates inward and moves to ~300 km by 40 h (Fig. 9a), accompanied by an increase in the low-level tangential wind. The time–radius plot of the outflow averaged over 100–200 hPa is shown in Fig. 9b, and the vertical motion averaged over 500–200 hPa is shown in Fig. 9c. Since the upper-level outflow not only influences convection but also responds to convection, one should not expect a perfect agreement between EFC and υr due to the stochastic nature of convection. Nevertheless, it is worth noting that the inward propagation of the EFC is accompanied by the inward propagation of upper-level outflow (Fig. 9b). In addition, strong upward motion in the inner-core region occurs from 18 h onward, about 10 h after the inward propagation of EFC (Fig. 9c). The strong upward motion is concomitant with the intensification of Typhoon Jongdari in the CTL run prior to 40 h, then the storm weakens, and w declines accordingly (Figs. 9c and 3b; also see the tangential wind in Fig. 9a). In contrast, with the removal of the UTCL, the EFC and the outflow at the large radii in RCL are much weaker than those in the CTL prior to 48 h, and there is no apparent inward propagation of EFC or upper-level outflow. The storm undergoes a steady but slow intensification before 48 h, and then intensifies more rapidly and has strong upward motions after 48 h (Figs. 9f and 3b; also see the tangential wind in Fig. 9d). The results suggest that the UTCL can enhance the storm intensity via eddy angular momentum convergence, consistent with the balanced vortex theory (e.g., Molinari and Vollaro 1989).


Time–radius plots of azimuthally averaged (a),(d) 300–100-hPa layer-mean eddy flux convergence (EFC; shaded; units: m s−1 day−1) and 800-hPa tangential wind (Vt; contours; units: m s−1); (b),(e) 200–100-hPa layer-mean radial wind (Vr; shaded; units: m s−1); and (c),(f) 500–200-hPa layer-mean vertical wind (w; shaded; units: m s−1) in the (a)–(c) CTL and (d)–(f) RCL runs.
Citation: Monthly Weather Review 149, 5; 10.1175/MWR-D-20-0271.1

Time–radius plots of azimuthally averaged (a),(d) 300–100-hPa layer-mean eddy flux convergence (EFC; shaded; units: m s−1 day−1) and 800-hPa tangential wind (Vt; contours; units: m s−1); (b),(e) 200–100-hPa layer-mean radial wind (Vr; shaded; units: m s−1); and (c),(f) 500–200-hPa layer-mean vertical wind (w; shaded; units: m s−1) in the (a)–(c) CTL and (d)–(f) RCL runs.
Citation: Monthly Weather Review 149, 5; 10.1175/MWR-D-20-0271.1
Time–radius plots of azimuthally averaged (a),(d) 300–100-hPa layer-mean eddy flux convergence (EFC; shaded; units: m s−1 day−1) and 800-hPa tangential wind (Vt; contours; units: m s−1); (b),(e) 200–100-hPa layer-mean radial wind (Vr; shaded; units: m s−1); and (c),(f) 500–200-hPa layer-mean vertical wind (w; shaded; units: m s−1) in the (a)–(c) CTL and (d)–(f) RCL runs.
Citation: Monthly Weather Review 149, 5; 10.1175/MWR-D-20-0271.1
c. Impacts of vertical wind shear


Time series of (a) the ventilation index, (b) entropy deficit (Xm), (c) potential intensity (MPI; units: m s−1), and (d) vertical wind shear (VWS; units: m s−1) in the CTL (black solid) and RCL (red dashed) simulations averaged within 500 km from the typhoon center.
Citation: Monthly Weather Review 149, 5; 10.1175/MWR-D-20-0271.1

Time series of (a) the ventilation index, (b) entropy deficit (Xm), (c) potential intensity (MPI; units: m s−1), and (d) vertical wind shear (VWS; units: m s−1) in the CTL (black solid) and RCL (red dashed) simulations averaged within 500 km from the typhoon center.
Citation: Monthly Weather Review 149, 5; 10.1175/MWR-D-20-0271.1
Time series of (a) the ventilation index, (b) entropy deficit (Xm), (c) potential intensity (MPI; units: m s−1), and (d) vertical wind shear (VWS; units: m s−1) in the CTL (black solid) and RCL (red dashed) simulations averaged within 500 km from the typhoon center.
Citation: Monthly Weather Review 149, 5; 10.1175/MWR-D-20-0271.1
Overall the UTCL has two direct, competing influences on the typhoon intensity: strong vertical wind shear that hinders typhoon intensification (Gray 1968; DeMaria and Kaplan 1994; Riemer et al. 2010), and enhanced upper-level outflow that strengthens the secondary circulation and facilitates storm intensification (e.g., Holland and Merrill 1984; Komaromi and Doyle 2018). The impacts of the UTCL on the typhoon intensity in the CTL run are consistent with previous studies: a UTCL may strengthen storm intensity at the early stage via enhanced upper-level outflow, but the negative impacts of strong vertical wind shear become dominant and result in storm stagnant or weakening when the UTCL moves too close to the storm (e.g., DeMaria and Kaplan 1994; Komaromi and Doyle 2018). In addition, the UTCL steers the typhoon to a different track and affects the typhoon intensity indirectly via the different SST and other thermal conditions along the altered track.
d. Terrain effects
Terrain may affect typhoon track via the separation of low-level and upper-level vortex centers and cause track deflection or discontinuity (Bender et al. 1985; Yeh and Elsberry 1993a,b; Wu 2001). Additionally, typhoons usually weaken quickly after landfall (Yu et al. 2017). In the CTL run, Typhoon Jongdari moves over Japan shortly after 60 h, and its intensity decreases rapidly (Fig. 3b). This is also the time period when the typhoon and UTCL are very close to each other, and the UTCL may weaken the typhoon by inducing strong vertical wind shear as suggested by previous studies (DeMaria et al. 1993; Komaromi and Doyle 2018). To quantify the effect of UTCL on the typhoon evolution after 60 h, the “Rland” experiment was carried out by replacing the Japanese Archipelago with ocean. As shown in Fig. 11, the simulated storm track in Rland is similar to that in the CTL run. Small track differences occur only after the storm has passed the location of the Japanese Archipelago, and the track in the CTL run turns slightly leftward compared to the track in the Rland run, which is consistent with typhoon track deflections reported in previous studies (Lin et al. 1999; Hsu et al. 2018). The storm intensity in the Rland run is close to that in the CTL run before 60 h, but diverges sharply afterward. In the Rland run, the typhoon still weakens after 60 h (albeit at a slower rate) even without the presence of Japanese Archipelago. This suggests that the weakening of Typhoon Jongdari at the later stage in the CTL run cannot be completely attributed to landfall, and another possible factor is the moderate vertical wind shear induced by the UTCL when it gets close to the typhoon (Fig. 10d).


(a) Typhoon tracks and (b) time series of CMSLP in the CTL and Rland simulations.
Citation: Monthly Weather Review 149, 5; 10.1175/MWR-D-20-0271.1

(a) Typhoon tracks and (b) time series of CMSLP in the CTL and Rland simulations.
Citation: Monthly Weather Review 149, 5; 10.1175/MWR-D-20-0271.1
(a) Typhoon tracks and (b) time series of CMSLP in the CTL and Rland simulations.
Citation: Monthly Weather Review 149, 5; 10.1175/MWR-D-20-0271.1
6. Summary and discussion
The impacts of an upper-level cold low (UTCL) on the track and intensity evolution of Typhoon Jongdari (2018) is examined in this study using numerical simulations. Typhoon Jongdari had an unusual counterclockwise looping track before making landfall over Japan. The CTL run that is initialized using the GFS analysis captures the track and intensity change reasonably well. In the RCL run, the UTCL is removed from the initial and boundary conditions using the piecewise potential vorticity inversion (PPVI) method, and the track and intensity evolution of the storm is substantially different from that in the CTL run. The PV budget analysis suggests that the motion of the storm is mainly controlled by horizontal PV advection, and that the storm motion is also affected by the vertical motion and diabatic heating fields, which are modified by the interaction between the UTCL and typhoon as well. The analysis of the steering flow in combination with PPVI suggests that the UTCL strongly modulates the steering flow and contributes to the looping track of the storm. Typhoon Jongdari and the UTCL rotate counterclockwise around each other as they move close together, a typical evolution under the Fujiwhara effect.
Further analysis showed that the UTCL could enhance the upper-level outflow around Typhoon Jongdari through eddy flux convergence of angular momentum and by reducing the inertial stability, and thus contributes to the intensification of Typhoon Jongdari. Typhoon Jongdari weakens at the later stage of the interaction in the CTL run due to strong ventilation. The strong ventilation can be attributed to increased vertical wind shear and reduced SST, which are related to the close proximity of the UTCL to the typhoon as well as a more northern typhoon track due to the impacts of the UTCL on the environmental steering flow. In addition, an experiment in which the Japanese Archipelago is replaced with ocean shows that the Typhoon weakens at a slow rate when it reaches Japan area and suggests that the environmental conditions, including the UTCL, and the land/terrain effects both contribute to the weakening of Typhoon Jongdari at its later stage.
This study provides quantitative evidence for the strong impacts of a UTCL on the track and intensity evolution of a nearby typhoon. The analysis of this particular case suggests that realistic representation of the upper-level low (or trough) and its interaction with a typhoon is important for skillful prediction of typhoon track and intensity. Our conclusion is based on one case study. Idealized simulations with UTCLs of different size and intensity and at different relative locations from the typhoon would provide more insight into the complex nature of the typhoon–UTCL interaction, which merits further study.
Acknowledgments
This work was jointly sponsored by the Science and Technology Innovation Project of Ningbo (Grant 2019B10025), the National Key R&D Program of China (2017YFC1502000), the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD), and Graduate research and innovation projects of Jiangsu Province (KYCX20_0912).
REFERENCES
Alaka, M. A., 1961: The occurrence of anomalous winds and their significance. Mon. Wea. Rev., 89, 482–494, https://doi.org/10.1175/1520-0493(1961)089<0482:TOOAWA>2.0.CO;2.
Barrett, B. S., E. R. Sanabia, S. C. Reynolds, J. K. Stapleton, and A. L. Borrego, 2016: Evolution of upper tropospheric outflow in Hurricanes Iselle and Julio (2014) in the Navy Global Environmental Model (NAVGEM) analyses and in satellite and dropsonde observations. J. Geophys. Res. Atmos., 121, 13 273–13 286, https://doi.org/10.1002/2016JD025656.
Bender, M. A., R. E. Tuleya, and Y. Kurihara, 1985: A numerical study of the effect of a mountain range on a landfalling tropical cyclone. Mon. Wea. Rev., 113, 567–583, https://doi.org/10.1175/1520-0493(1985)113<0567:ANSOTE>2.0.CO;2.
Bi, M. T., M. Li, and X. Shen, 2015: Interactions between Typhoon Megi (2010) and a low-frequency monsoon gyre. J. Atmos. Sci., 72, 2682–2702, https://doi.org/10.1175/JAS-D-14-0269.1.
Bister, M., and K. Emanuel, 2002: Low frequency variability of tropical cyclone potential intensity. 1. Interannual to interdecadal variability. J. Geophys. Res., 107, 4801, https://doi.org/10.1029/2001JD000776.
Bosart, L. F., W. E. Bracken, J. Molinari, C. S. Velden, and P. G. Black, 2000: Environmental influences on the rapid intensification of Hurricane Opal (1995) over the Gulf of Mexico. Mon. Wea. Rev., 128, 322–352, https://doi.org/10.1175/1520-0493(2000)128<0322:EIOTRI>2.0.CO;2.
Brand, S., 1970: Interaction of binary tropical cyclones of the western North Pacific Ocean. J. Appl. Meteor., 9, 433–441, https://doi.org/10.1175/1520-0450(1970)009<0433:IOBTCO>2.0.CO;2.
Bryan, G., 2008: On the computation of pseudoadiabatic entropy and equivalent potential temperature. Mon. Wea. Rev., 136, 5239–5245, https://doi.org/10.1175/2008MWR2593.1.
Campetella, C. M., and N. E. Possia, 2007: Upper-level cut-off lows in southern South America. Meteor. Atmos. Phys., 96, 181–191, https://doi.org/10.1007/s00703-006-0227-2.
Carr, L. E., and R. L. Elsberry, 1995: Monsoonal interactions leading to sudden tropical cyclone track changes. Mon. Wea. Rev., 123, 265–290, https://doi.org/10.1175/1520-0493(1995)123<0265:MILTST>2.0.CO;2.
Chan, J. C., and W. M. Gray, 1982: Tropical cyclone movement and surrounding flow relationships. Mon. Wea. Rev., 110, 1354–1374, https://doi.org/10.1175/1520-0493(1982)110<1354:TCMASF>2.0.CO;2.
Chan, J. C., and R. T. Williams, 1987: Analytical and numerical studies of the beta-effect in tropical cyclone motion. Part I: Zero mean flow. J. Atmos. Sci., 44, 1257–1265, https://doi.org/10.1175/1520-0469(1987)044<1257:AANSOT>2.0.CO;2.
Chan, J. C., F. M. Ko, and Y. M. Lei, 2002: Relationship between potential vorticity tendency and tropical cyclone motion. J. Atmos. Sci., 59, 1317–1336, https://doi.org/10.1175/1520-0469(2002)059<1317:RBPVTA>2.0.CO;2.
Chang, S. W. J., 1983: A numerical study of the interaction between two tropical cyclones. Mon. Wea. Rev., 111, 1806–1817, https://doi.org/10.1175/1520-0493(1983)111<1806:ANSOTI>2.0.CO;2.
Chen, G., and L. F. Chou, 1994: An investigation of cold vortices in the upper troposphere over the western North Pacific during the warm season. Mon. Wea. Rev., 122, 1436–1448, https://doi.org/10.1175/1520-0493(1994)122<1436:AIOCVI>2.0.CO;2.
Chen, T. C., and Coauthors, 2001: Summer upper-level vortex over the North Pacific. Bull. Amer. Meteor. Soc., 82, 1991–2006, https://doi.org/10.1175/1520-0477-82.9.1991.
Chen, X., Y. Wang, and K. Zhao, 2015: Synoptic flow patterns and large-scale characteristics associated with rapidly intensifying tropical cyclones in the South China Sea. Mon. Wea. Rev., 143, 64–87, https://doi.org/10.1175/MWR-D-13-00338.1.
Choudhury, D., and S. Das, 2017: The sensitivity to the microphysical schemes on the skill of forecasting the track and intensity of tropical cyclones using WRF-ARW model. J. Earth Syst, 126, 57, https://doi.org/10.1007/s12040-017-0830-2.
Davidson, N. E., C. M. Nguyen, and M. J. Reeder, 2008: Downstream development during the rapid intensification of Hurricanes Opal and Katrina: The distant trough interaction problem. 28th Conf. on Hurricanes and Tropical Meteorology, Orlando, FL, Amer. Meteor. Soc., 9B.4, https://ams.confex.com/ams/28Hurricanes/techprogram/paper_138060.htm.
Davis, C. A., 1992: Piecewise potential vorticity inversion. J. Atmos. Sci., 49, 1397–1411, https://doi.org/10.1175/1520-0469(1992)049<1397:PPVI>2.0.CO;2.
Davis, C. A., and K. A. Emanuel, 1991: Potential vorticity diagnostics of cyclogenesis. Mon. Wea. Rev., 119, 1929–1953, https://doi.org/10.1175/1520-0493(1991)119<1929:PVDOC>2.0.CO;2.
DeMaria, M., 1996: The effect of vertical shear on tropical cyclone intensity change. J. Atmos. Sci., 53, 2076–2088, https://doi.org/10.1175/1520-0469(1996)053<2076:TEOVSO>2.0.CO;2.
DeMaria, M., and J. Kaplan, 1994: A Statistical Hurricane Intensity Prediction Scheme (SHIPS) for the Atlantic basin. Wea. Forecasting, 9, 209–220, https://doi.org/10.1175/1520-0434(1994)009<0209:ASHIPS>2.0.CO;2.
DeMaria, M., J.-J. Baik, and J. Kaplan, 1993: Upper-level eddy angular momentum fluxes and tropical cyclone intensity change. J. Atmos. Sci., 50, 1133–1147, https://doi.org/10.1175/1520-0469(1993)050<1133:ULEAMF>2.0.CO;2.
DeMaria, M., C. R. Sampson, J. A. Knaff, and K. D. Musgrave, 2014: Is tropical cyclone intensity guidance improving? Bull. Amer. Meteor. Soc., 95, 387–398, https://doi.org/10.1175/BAMS-D-12-00240.1.
Dudhia, J., 1989: Numerical study of convection observed during the winter monsoon experiment using a mesoscale two-dimensional model. J. Atmos. Sci., 46, 3077–3107, https://doi.org/10.1175/1520-0469(1989)046<3077:NSOCOD>2.0.CO;2.
Emanuel, K., and F. Zhang, 2016: On the predictability and error sources of tropical cyclone intensity forecasts. J. Atmos. Sci., 73, 3739–3747, https://doi.org/10.1175/JAS-D-16-0100.1.
Ferrier, B. S., and Coauthors, 2002: Implementation of a new grid-scale cloud and precipitation scheme in the NCEP Eta Model. 19th Conf. on Weather Analysis and Forecasting/15th Conf. on Numerical Weather, Seattle, WA, Amer. Meteor. Soc., 280–283.
Fischer, M. S., B. H. Tang, and K. L. Corbosiero, 2017: Assessing the influence of upper-tropospheric troughs on tropical cyclone intensification rates after genesis. Mon. Wea. Rev., 145, 1295–1313, https://doi.org/10.1175/MWR-D-16-0275.1.
Fischer, M. S., B. H. Tang, and K. L. Corbosiero, 2019: A climatological analysis of tropical cyclone rapid intensification in environments of upper-tropospheric troughs. Mon. Wea. Rev., 147, 3693–3719, https://doi.org/10.1175/MWR-D-19-0013.1.
Fovell, R. G., K. L. Corbosiero, and H. Kuo, 2009: Cloud microphysics impact on hurricane track as revealed in idealized experiments. J. Atmos. Sci., 66, 1764–1778, https://doi.org/10.1175/2008JAS2874.1.
Frank, W. M., and E. A. Ritchie, 2001: Effects of vertical wind shear on the intensity and structure of numerically simulated hurricanes. Mon. Wea. Rev., 129, 2249–2269, https://doi.org/10.1175/1520-0493(2001)129<2249:EOVWSO>2.0.CO;2.
Fujiwhara, S., 1921: The mutual tendency toward symmetry of motion and its application as a principle in meteorology. Quart. J. Roy. Meteor. Soc., 47, 287–292, https://doi.org/10.1002/qj.49704720010.
Fujiwhara, S., 1923: On the growth and decay of vortical systems. Quart. J. Roy. Meteor. Soc., 49, 75–104, https://doi.org/10.1002/qj.49704920602.
Ge, X. Y., Z. Yan, M. Peng, M. Bi, and T. Li, 2018: Sensitivity of tropical cyclone track to the vertical structure of a nearby monsoon gyre. J. Atmos. Sci., 75, 2017–2028, https://doi.org/10.1175/JAS-D-17-0201.1.
George, J. E., and W. M. Gray, 1976: Tropical cyclone motion and surrounding parameter relationships. J. Appl. Meteor., 15, 1252–1264, https://doi.org/10.1175/1520-0450(1976)015<1252:TCMASP>2.0.CO;2.
Gray, W. M., 1968: Global view of the origin of tropical disturbances and storms. Mon. Wea. Rev., 96, 669–700, https://doi.org/10.1175/1520-0493(1968)096<0669:GVOTOO>2.0.CO;2.
Hanley, D. E., J. Molinari, and D. Keyser, 2001: A composite study of the interactions between tropical cyclones and upper-tropospheric troughs. Mon. Wea. Rev., 129, 2570–2584, https://doi.org/10.1175/1520-0493(2001)129<2570:ACSOTI>2.0.CO;2.
Harr, P. A., R. L. Elsberry, and J. C. Chan, 1996: Transformation of a large monsoon depression to a tropical storm during TCM-93. Mon. Wea. Rev., 124, 2625–2643, https://doi.org/10.1175/1520-0493(1996)124<2625:TOALMD>2.0.CO;2.
Hendricks, E. A., M. S. Peng, X. Ge, and T. Li, 2011: Performance of a dynamic initialization scheme in the Coupled Ocean–Atmosphere Mesoscale Prediction System for Tropical Cyclones (COAMPS-TC). Wea. Forecasting, 26, 650–663, https://doi.org/10.1175/WAF-D-10-05051.1.
Holland, G. J., 1983: Tropical cyclone motion: Environmental interaction plus a beta effect. J. Atmos. Sci., 40, 328–342, https://doi.org/10.1175/1520-0469(1983)040<0328:TCMEIP>2.0.CO;2.
Holland, G. J., and R. T. Merrill, 1984: On the dynamics of tropical cyclone structural changes. Quart. J. Roy. Meteor. Soc., 110, 723–745, https://doi.org/10.1002/qj.49711046510.
Hong, S.-Y., and J.-O. J. Lim, 2006: The WRF single-moment 6-class microphysics scheme (WSM6). J. Korean Meteor. Soc., 42, 129–151.
Hong, S.-Y., Y. Noh, and J. Dudhia, 2006: A new vertical diffusion package with an explicit treatment of entrainment processes. Mon. Wea. Rev., 134, 2318–2341, https://doi.org/10.1175/MWR3199.1.
Hoskins, B. J., M. E. Mclntyre, and A. W. Robertson, 1985: On the use and significance of isentropic potential-vorticity maps. Quart. J. Roy. Meteor. Soc., 111, 877–946, https://doi.org/10.1002/qj.49711147002.
Hsu, L. H., S. H. Su, and G. Robert, 2018: On typhoon track deflections near the east coast of Taiwan. Mon. Wea. Rev., 146, 1495–1510, https://doi.org/10.1175/MWR-D-17-0208.1.
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.
Kelly, W. E., and D. Mock, 1982: A diagnostic study of upper tropospheric cold lows over the western North Pacific. Mon. Wea. Rev., 110, 471–480, https://doi.org/10.1175/1520-0493(1982)110<0471:ADSOUT>2.0.CO;2.
Komaromi, W. A., and J. D. Doyle, 2018: On the dynamics of tropical cyclone and trough interactions. J. Atmos. Sci., 75, 2687–2709, https://doi.org/10.1175/JAS-D-17-0272.1.
Lander, M., 1994: Description of a monsoon gyre and its effects on the tropical cyclones in the western North Pacific during August 1991. Wea. Forecasting, 9, 640–654, https://doi.org/10.1175/1520-0434(1994)009<0640:DOAMGA>2.0.CO;2.
Lander, M., and G. J. Holland, 1993: On the interaction of tropical-cyclone-scale vortices: Observations. Quart. J. Roy. Meteor. Soc., 119, 1347–1361, https://doi.org/10.1002/qj.49711951406.
Lei, L., Y. Ge, Z. Tan, and X. Bao, 2020: An evaluation and improvement of tropical cyclone prediction in the western North Pacific basin from global ensemble forecasts. Sci. China Earth Sci., 63, 12–26, https://doi.org/10.1007/s11430-019-9480-8.
Leroux, M. D., M. Plu, D. Barbary, F. Roux, and P. Arbogast, 2013: Dynamical and physical processes leading to tropical cyclone intensification under upper-level trough forcing. J. Atmos. Sci., 70, 2547–2565, https://doi.org/10.1175/JAS-D-12-0293.1.
Leroux, M. D., M. Plu, and F. Roux, 2016: On the sensitivity of tropical cyclone intensification under upper-level trough forcing. Mon. Wea. Rev., 144, 1179–1202, https://doi.org/10.1175/MWR-D-15-0224.1.
Li, Y., L. Guo, Y. Ying, and S. Hu, 2012: Impacts of upper-level cold vortex on the rapid change of intensity and motion of Typhoon Meranti (2010). J. Trop. Meteor., 18, 207–219.
Lin, Y. L., R. D. Rarley, and H. D. Orville, 1983: Bulk parameterization of the snow field in a cloud model. J. Climate Appl. Meteor., 22, 1065–1092, https://doi.org/10.1175/1520-0450(1983)022<1065:BPOTSF>2.0.CO;2.
Lin, Y. L., J. Han, D. W. Hamilton, and C.-Y. Huang, 1999: Orographic influence on a drifting cyclone. J. Atmos. Sci., 56, 534–562, https://doi.org/10.1175/1520-0469(1999)056<0534:OIOADC>2.0.CO;2.
Liou, Y. A., J.-C. Liu, M.-X. Wu, Y.-J. Lee, C.-H. Cheng, C.-P. Kuei, and R.-M. Hong, 2016: Generalized empirical formulas of threshold distance to characterize cyclone-cyclone interactions. IEEE Trans. Geosci. Remote Sens., 54, 3502–3512, https://doi.org/10.1109/TGRS.2016.2519538.
Maw, K. W., and J. Min, 2017: Impacts of microphysics schemes and topography on the prediction of the heavy rainfall in western Myanmar associated with tropical cyclone ROANU (2016). Adv. Meteor., 2017, 3252503, https://doi.org/10.1155/2017/3252503.
Mlawer, E. J., S. J. Taubman, P. D. Brown, M. J. Iacono, S. A. Clough, 1997: Radiative transfer for inhomogeneous atmosphere: RRTM, a validated correlated-k model for the longwave. J. Geophys. Res., 102, 16 663–16 682, https://doi.org/10.1029/97JD00237.
Molinari, J., and D. Vollaro, 1989: External influences on hurricane intensity. Part I: Outflow layer eddy momentum fluxes. J. Atmos. Sci., 46, 1093–1105, https://doi.org/10.1175/1520-0469(1989)046<1093:EIOHIP>2.0.CO;2.
Molinari, J., and D. Vollaro, 2014: Symmetric instability in the outflow layer of a major hurricane. J. Atmos. Sci., 71, 3739–3746, https://doi.org/10.1175/JAS-D-14-0117.1.
Molinari, J., S. Skubis, and D. Vollaro, 1995: External influences on hurricane intensity. Part III: Potential vorticity structure. J. Atmos. Sci., 52, 3593–3606, https://doi.org/10.1175/1520-0469(1995)052<3593:EIOHIP>2.0.CO;2.
Nguyen, L. T., and J. Molinari, 2015: Simulation of the downshear reformation of a tropical cyclone. J. Atmos. Sci., 72, 4529–4551, https://doi.org/10.1175/JAS-D-15-0036.1.
Patla, J. E., D. Stevens, and G. M. Barnes, 2009: A conceptual model for the influence of TUTT cells on tropical cyclone motion in the northwest Pacific Ocean. Wea. Forecasting, 24, 1215–1235, https://doi.org/10.1175/2009WAF2222181.1.
Peirano, C. M., K. L. Corbosiero, and B. H. Tang, 2016: Revisiting trough interactions and tropical cyclone intensity change. Geophys. Res. Lett., 43, 5509–5515, https://doi.org/10.1002/2016GL069040.
Peng, M. S., and C. A. Reynolds, 2005: Double trouble for typhoon forecasters. Geophys. Res. Lett., 32, L02810, https://doi.org/10.1029/2004GL021680.
Peng, M. S., and C. A. Reynolds, 2006: Sensitivity of tropical cyclones forecasts as revealed by singular vectors. J. Atmos. Sci., 63, 2508–2528, https://doi.org/10.1175/JAS3777.1.
Rappaport, E. N., and Coauthors, 2009: Advances and challenges at the National Hurricane Center. Wea. Forecasting, 24, 395–419, https://doi.org/10.1175/2008WAF2222128.1.
Rappin, E. D., M. C. Morgan, and G. J. Tripoli, 2011: The impact of outflow environment on tropical cyclone intensification and structure. J. Atmos. Sci., 68, 177–194, https://doi.org/10.1175/2009JAS2970.1.
Riemer, M., M. T. Montgomery, and M. E. Nicholls, 2010: A new paradigm for intensity modification of tropical cyclones: Thermodynamic impact of vertical wind shear on the inflow layer. Atmos. Chem. Phys., 10, 3163–3188, https://doi.10.5194/ACP-10-3163-2010.
Ritchie, E. A., and G. J. Holland, 1993: On the interaction of tropical cyclone-scale vortices. II: Discrete vortex patches. Quart. J. Roy. Meteor. Soc., 119, 1363–1379, https://doi.org/10.1002/QJ.49711951407.
Rodgers, E. B., S. W. Chang, J. Stout, J. Steranka, and J.-J. Shi, 1991: Satellite observations of variations in tropical cyclone convection caused by upper-tropospheric troughs. J. Appl. Meteor., 30, 1163–1184, https://doi.org/10.1175/1520-0450(1991)030<1163:SOOVIT>2.0.CO;2.
Sadler, J. C., 1975: The monsoon circulation and cloudiness over the GATE area. Mon. Wea. Rev., 103, 369–387, https://doi.org/10.1175/1520-0493(1975)103<0369:TMCACO>2.0.CO;2.
Schubert, W. H., and J. J. Hack, 1982: Inertial stability and tropical cyclone development. J. Atmos. Sci., 39, 1687–1697, https://doi.org/10.1175/1520-0469(1982)039<1687:ISATCD>2.0.CO;2.
Shapiro, L. J., 1996: The motion of Hurricane Gloria: A potential vorticity diagnosis. Mon. Wea. Rev., 124, 2497–2508, https://doi.org/10.1175/1520-0493(1996)124<2497:TMOHGA>2.0.CO;2.
Shapiro, L. J., 1999: Potential vorticity asymmetries and tropical cyclone motion. Mon. Wea. Rev., 127, 124–131, https://doi.org/10.1175/1520-0493(1999)127<0124:PVAATC>2.0.CO;2.
Shi, J. J., S. W. Chang, and S. Raman, 1997: Interaction between Hurricane Florence (1988) and an upper-tropospheric westerly trough. J. Atmos. Sci., 54, 1231–1247, https://doi.org/10.1175/1520-0469(1997)054<1231:IBHFAA>2.0.CO;2.
Skamarock, W. C., and Coauthors, 2008: A description of the Advanced Research WRF version 3. NCAR Tech. Note NCAR/TN-475+STR, 113 pp., https://doi.org/10.5065/D68S4MVH.
Tang, B., and K. Emanuel, 2012: A ventilation index for tropical cyclones. Bull. Amer. Meteor. Soc., 93, 1901–1912, https://doi.org/10.1175/BAMS-D-11-00165.1.
Tewari, M. F., and Coauthors, 2004:Implementation and verification of the unified Noah land surface model in the WRF model. 20th Conf on Weather Analysis and Forecasting/16th Conf on Numerical Weather Prediction, Seattle, WA, Amer. Meteor. Soc., 14.2a, https://ams.confex.com/ams/84Annual/techprogram/paper_69061.htm.
Velden, C. S., and L. M. Leslie, 1991: The basic relationship between tropical cyclone intensity and the depth of the environmental steering layer in the Australian region. Wea. Forecasting, 6, 244–253, https://doi.org/10.1175/1520-0434(1991)006<0244:TBRBTC>2.0.CO;2.
Wang, Y., and G. J. Holland, 1995: On the interaction of tropical cyclone-scale vortices. IV: Baroclinic vortices. Quart. J. Roy. Meteor. Soc., 121, 95–126, https://doi.org/10.1002/qj.49712152106.
Wang, Y., and C. C. Wu, 2004: Current understanding of tropical cyclone structure and intensity changes—A review. Meteor. Atmos. Phys., 87, 257–278, https://doi.org/10.1007/s00703-003-0055-6.
Wei, N., Y. Li, D.-L. Zhang, Z. Mai, and S.-Q. Yang, 2016: A statistical analysis of the relationship between upper-tropospheric cold low and tropical cyclone track and intensity change over the western North Pacific. Mon. Wea. Rev., 144, 1805–1822, https://doi.org/10.1175/MWR-D-15-0370.1.
Wen, D., and Coauthors, 2019: An ensemble analysis on abrupt north turning of Typhoon Meranti (1010) under the influence of an upper tropospheric cold low (in Chinese). Chin. J. Atmos. Sci., 43, 730–740.
Wu, C. C., 2001: Numerical simulation of Typhoon Gladys (1994) and its interaction with Taiwan terrain using the GFDL hurricane model. Mon. Wea. Rev., 129, 1533–1549, https://doi.org/10.1175/1520-0493(2001)129<1533:NSOTGA>2.0.CO;2.
Wu, C. C., and K. Emanuel, 1995a: Potential vorticity diagnostics of hurricane movement. Part I: A case study of Hurricane Bob (1991). Mon. Wea. Rev., 123, 69–92, https://doi.org/10.1175/1520-0493(1995)123<0069:PVDOHM>2.0.CO;2.
Wu, C. C., and K. Emanuel, 1995b: Potential vorticity diagnostics of hurricane movement. Part II: Tropical Storm Ana (1991) and Hurricane Andrew (1992). Mon. Wea. Rev., 123, 93–109, https://doi.org/10.1175/1520-0493(1995)123<0093:PVDOHM>2.0.CO;2.
Wu, C. C., and Y. Kurihara, 1996: A numerical study of the feedback mechanisms of hurricane–environment interaction on hurricane movement from the potential vorticity perspective. J. Atmos. Sci., 53, 2264–2282, https://doi.org/10.1175/1520-0469(1996)053<2264:ANSOTF>2.0.CO;2.
Wu, C. C., T. S. Huang, W. P. Huang, and K.-H. Chou, 2003: A new look at the binary interaction: Potential vorticity diagnosis of the unusual southward movement of Tropical Storm Bopha (2000) and its interaction with Super Typhoon Saomai (2000). Mon. Wea. Rev., 131, 1289–1300, https://doi.org/10.1175/1520-0493(2003)131<1289:ANLATB>2.0.CO;2.
Wu, C. C., K. W. Kevin, and Y.-Y. Lo, 2009: Numerical study of the rainfall event due to the interaction of Typhoon Babs (1998) and the northeasterly monsoon. Mon. Wea. Rev., 137, 2049–2064, https://doi.org/10.1175/2009MWR2757.1.
Wu, C. C., K. K. W. Cheung, J.-H. Chen, and C.-C . Chang, 2010: The impact of Tropical Storm Paul (1999) on the motion and rainfall associated with Tropical Storm Rachel (1999) near Taiwan. Mon. Wea. Rev., 138, 1635–1650, https://doi.org/10.1175/2009MWR3021.1.
Wu, L. G., and B. Wang, 2000: A potential vorticity tendency diagnostic approach for tropical cyclone motion. Mon. Wea. Rev., 128, 1899–1911, https://doi.org/10.1175/1520-0493(2000)128<1899:APVTDA>2.0.CO;2.
Yan, Z. Y., X. Ge, M. Peng, and T . Lim, 2019: Does monsoon gyre always favour tropical cyclone rapid intensification? Quart. J. Roy. Meteor. Soc., 145, 2685–2697, https://doi.org/10.1002/qj.3586.
Yang, C. C., C.-C. Wu, K.-H. Chou, and C.-Y . Lee, 2008: Binary interaction between Typhoons Fengshen (2002) and Fungwong (2002) based on the potential vorticity diagnosis. Mon. Wea. Rev., 136, 4593–4611, https://doi.org/10.1175/2008MWR2496.1.
Yeh, T. C., and R. L. Elsberry, 1993a: Interaction of typhoons with the Taiwan orography. Part I: Upstream track deflections. Mon. Wea. Rev., 121, 3193–3212, https://doi.org/10.1175/1520-0493(1993)121<3193:IOTWTT>2.0.CO;2.
Yeh, T. C., and R. L. Elsberry, 1993b: Interaction of typhoons with the Taiwan orography. Part II: Continuous and discontinuous tracks across the island. Mon. Wea. Rev., 121, 3213–3233, https://doi.org/10.1175/1520-0493(1993)121<3213:IOTWTT>2.0.CO;2.
Yu, Z. F., Y. Wang, H. Xu, N. Davidson, Y. Chen, Y. Chen, and H. Yu, 2017: On the relationship between intensity and rainfall distribution in tropical cyclones making landfall over China. J. Appl. Meteor. Climatol., 56, 2883–2901, https://doi.org/10.1175/JAMC-D-16-0334.1.
