Numerical Analysis on the Effects of Binary Interaction between Typhoons Tembin and Bolaven in 2012

,e binary interaction is one of the most challenging factors to improve the forecast accuracy of multiple tropical cyclones (TCs) in close vicinity. ,e effect of binary interaction usually results in anomalous track and variable intensity of TCs. A typical interaction type, one-way influence mode, has been investigated by many studies which mainly focused on the anomalous track and recordbreaking precipitation, such as typhoons Morakot and Goni. In this paper, a typical case of this type, typhoons Tembin and Bolaven, occurred in the western North Pacific in August 2012, was selected to study how one typhoon impacts the track and intensity of the other one. ,e vortex of Tembin or Bolaven and the monsoon circulation were removed by a TC bogus scheme and a low-pass Lanczos filter, respectively, to carry out the numerical experiments. ,e results show that the presence of monsoon made the binary interaction more complex by affecting the tracks and the translation speeds of the TCs. ,e influence of Bolaven on the track of Tembin was more significant than the influence of Tembin on the track of Bolaven, and the looping track of Tembin was also affected by the external surrounding circulation associated with Bolaven. ,e absence of Tembin was not conducive to the development of Bolaven due to stronger vertical wind shear condition and the less kinetic energy being transported to upper troposphere. Note that the above analysis demonstrates the interacting processes between TCs and sheds some light on the prediction of binary TCs.


Introduction
e anomalous track and intensity of tropical cyclone (TC) caused by binary typhoon effect is an important topic in atmospheric sciences. In the early 1920s, Fujiwhara [1] found that two vortices would interact each other when they are close to each other. Brand [2], who analyzed the relationship between angular change and separation distance for 22 pairs of binary typhoons in the western North Pacific, suggested that the binary typhoon effect was related to the scale and intensity of typhoons and further indicated that two typhoons start to rotate when the separation distance is less than 750 miles. Lander and Holland [3] used a reference frame with unweighted approximation to centroid to make a detailed analysis of the behaviour of multiple-vortex interactions in the western North Pacific.
ey presented a common model of binary interaction that contains the major elements of approach, sudden capture, and then a period of steady cyclonic orbit, followed by a sudden escape or a merger. It further indicated that the vortices would be captured when the separation distance is less than 600 miles. Since many cases of binary typhoons whose tracks are quite different from the mutual advection (Fujiwhara type), Carr et al. [4] and Carr and Elsberry [5] categorized binary typhoon interaction processes and proposed four conceptual models: the direct interaction conceptual model, which includes one-way influence, mutual interaction, and merge modes; the semidirect interaction conceptual model; the indirect interaction conceptual model; and binary cyclone contribution to reverseoriented trough formation (a monsoon trough with a southwest-northeast orientation, which will occur when two or more TCs locating along a zonally oriented monsoonal cloud band move relative to one another [6]). However, due to the lack of observations, it is difficult to explain how the binary typhoon effect works between the TCs. For this reason, the numerical simulation is considered one of the available methods to investigate it.
Xu et al. [7] studied a binary typhoons case, Morakot and Goni, which occurred in 2009. ey pointed out that the existence of Goni played an important role in maintaining the intensity of Morakot, resulting in the heavy rainfall increasing over Taiwan, as compared to the experiment removing Goni's circulation. But, Wu et al. [8] drew an opposite conclusion for the same binary typhoons case. ey carried out a set of experiments with modifying the intensity of Morakot in its initial conditions to evaluate the factors affecting the interaction between the two TCs. e results showed that Morakot had a greater impact on the motion and was conducive to intensifying the strength of Goni, whereas Goni, located upstream of the southwesterly flow, intercepted part of the water vapor being transported to Morakot, making Morakot weaken. Xu et al. [9] applied the same experiment design to study another binary typhoons case, Bopha and Saomai, which occurred in 2006. ey indicated that Saomai's track and intensity responded to the change in the intensity of Bopha and Bopha had opposite effects at different development stages of Saomai. e interaction between typhoons Tembin and Bolaven resulted in strong winds and heavy rain in both Taiwan and Fujian provinces of China in August 2012. Reviewing the evolution of Tembin and Bolaven (Figure 1(a)) shows that the former was characterized by a looping track and slow translation speeds, while the latter was characterized by northwestward motion. It is interesting that the tracks of Tembin and Bolaven are quite similar to Morakot and Goni's, especially the looping tracks of Tembin and Goni. According to the studies of Carr et al., the two pairs of binary typhoons can be both considered as the one-way influence of the direct interaction mode in which the track of a smaller TC that is embedded in the circulation of a larger TC has an orbiting motion but causes no significant track alteration of the larger TC due to its much smaller circulation. In other words, the track of a smaller TC (i.e., Tembin and Goni) that is dominated by a larger TC (i.e., Bolaven and Morakot) has a cyclonic looping motion. Previous studies [10,11] on Tembin and Bolaven mainly focused on the complicated tracks, but these were insufficient for the investigation of the two typhoons' intensities. As mentioned above, for the oneway influence of interaction mode, such as Morakot and Goni, the answer to the question of how the binary typhoon effect affects the track and intensity of TCs remains unclear. erefore, it is worthy carrying out a study on this interaction mode to investigate how one typhoon impacts the track and intensity of another one.

Data.
e data used as the initial and boundary conditions in this study are the final (FNL) global analysis data (resolution: 1°× 1°) from the National Center for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) with the time resolution of 6 hours and the real-time global analysis sea surface temperature (SST) dataset (RTG_SST, 0.5°× 0.5°) from the National Center for Environmental Prediction/National Oceanic and Atmospheric Administration (NCEP/NOAA) with the time resolution of 24 hours. e best track of the TCs is obtained from the International Best Track Archive for Climate Stewardship (IBTrACS) [12] with the 6-hour time resolution, which contains the information about the central location of the TCs, the minimum pressure, the maximum  wind near the TCs' center, the maximum wind radius, etc. 2.2. TC Bogus Scheme. TC bogus is a TC initial processing scheme of Fredrick et al., and the method description partly reproduces their wording [13]. e scheme can either introduce a new vortex or remove an existing one. e vortex removal is done as follows. Firstly, use the best track location of the TC and 400 km as the search radius to look for the maximum relative vorticity in the pressure field nearest the surface. e grid point where the maximum relative vorticity is located serves as the center of the vortex to be removed. e vorticity and divergence outside the radius of 300 km of the vortex's center are set to zero by the scheme. en, use the successive overrelaxation method to solve the following equations for the stream function, velocity potential, nondivergent wind, and divergent wind: where ψ, ζ, v ψ , χ, δ, and v χ represent the stream function, relative vorticity, nondivergent wind, velocity potential, horizontal divergent, and divergent wind, respectively. Once the divergent wind is calculated, it is subtracted from the first-guess U and V wind fields. Next, the geopotential height anomaly, the geostrophic wind, and temperature anomaly are subtracted from the first-guess field. Finally, the surface and sea-level pressure perturbations are subtracted. So, only the background flow where the first-guess storm was originally located is left.

Numerical Simulation Settings.
With the purpose of studying the one-way influence of the interaction mode for typhoons Tembin and Bolaven, the advanced Weather Research and Forecasting Model (WRF, V3.6) is used for the numerical simulation. e model integration time is from 1200 UTC 23 August to 1200 UTC 27 August 2012. Euler mass coordinate, Runge-Kutta 3rd-order time integration option, and two-way nesting with double-nested grids are applied. e model is centered at 28.0°N, 125.0°E with horizontal resolutions of 27 km for D01 and 9 km for D02, 31 eta levels in the vertical direction, and a model top level at 50 hPa. e model configurations include a 241 × 203 grid in the outer nest (D01) and a 268 × 274 grid in the inner nest (D02). e physical parameterization schemes used in the numerical simulation include the WRF single-moment 6-class microphysics scheme [14], the rapid radiative transfer model longwave radiation scheme [15], the Dudhia shortwave radiation scheme [16], the uni ed Noah landsurface model land-surface scheme [17], the Yonsei University planetary boundary layer scheme [18], and the Kain-Fritsch convective scheme [19]. e experiment setting above is denoted as CTL experiment. To remove the Tembin and Bolaven vortices, TC bogus scheme provided by the WRF modeling system is applied for the sensitivity experiments (denoted as NoT experiment and NoB experiment, respectively) with the same experiment settings as the CTL experiment.

Synoptic Situation Analysis.
e synoptic evolution was shown by the 500 hPa geopotential height elds from the European Centre for Medium-Range Weather Forecasts (ECMWF) ERA-Interim reanalyses data from 0000 UTC 24 August to 0000 UTC 27 August ( Figure 2). ere were one trough and one ridge in midhigh latitudes regions locating in China and Japan, respectively. e circulation center of subtropical high (the contour of 576 dagpm) was located in the ocean and extended westward from 0000 UTC 24 August to 0000 UTC 25 August, with the westernmost contours of 576 dagpm locating at 110°E. During this time, the two TCs moved westward along the periphery of the subtropical high and the distance between typhoons Tembin and Bolaven gradually became smaller. After 0000 UTC 26 August, the trough started to deepen and move eastward and the intensity of the subtropical high weakened with its main body retreating eastward and stretching southward in a meridional pattern. During this period, Tembin moved southward and then turned northeastward, and Bolaven constantly moved northward. It should be noticed that when Tembin made a looping track, i.e., from 0000 UTC 25 August to 1200 UTC 27 August, its translation speed was slower than Bolaven's. To clearly describe the behaviour of interacting tropical vortices, the centroid-relative motion for the two typhoons is displayed in Figure 1(b). e two TCs started to approach the centroid at 0000 UTC 20 August, and a sudden onset of cyclonic orbit occurred after 1800 UTC 22 August, followed by a sudden capture. e distance between the two TCs tended to become slightly larger, and a period of stable cyclonic orbit began at 0600 UTC 26 August and ended at 0600 UTC 28 August. A sudden escape occurred when the separation distance gradually became larger after 0600 UTC 28 August.
At the point of interest, typhoons Tembin and Bolaven started to cyclonically orbit when the separation distance was about 800 miles (the dash circle in Figure 1(b) indicates Brand's typical cyclonic-orbit onset distance of 750 miles, about 13 degrees in latitude). In addition, the actual orbit track was not distinctly cyclonically curved. ese characteristics indicate that the tracks of typhoons Tembin and Bolaven were likely to be a ected by the monsoon over the western North Paci c. To verify that, an experiment (denoted as NoM experiment) applying a low-pass Lanczos lter with a 10-day cuto period to the wind elds to remove the low-frequency ows (original elds minus the result from the low-pass Lanczos lter) [20] is carried out to determine the binary interaction behaviour in the absence of the monsoon. e results will be discussed in Section 3.3.

Simulation Veri cation.
e best track data from the Japan Meteorological Agency (JMA) and Joint Typhoon Warning Center (JTWC) are used to verify the simulation results of the track and intensity of the two typhoons. As shown in Figure 3, the model succeeded well in simulating the tracks of Tembin and Bolaven in the CTL experiment, especially Tembin's looping track, but the translation speed of the two typhoons (in white circle) was faster than the best track. On the other hand, as shown in Figure 4, although the  Figure 1: Tracks of Typhoons Tembin and Bolaven (a) and a "Fujiwhara" diagram for the binary interaction between them (b) from 0000 UTC 20 August to 0000 UTC 1 September (interval: 6 hours; red dots: Tembin; blue dots: Bolaven). In (b), distances are in degrees latitude. Dots are centroid-relative positions. e major elements of "Approach", "Capture", "Orbit", and "Escape" are labelled. e radius of dash circle is 6.5 degrees. e origin of coordinate denotes the geographical center between the two vortices.
di erences in the minimum sea-level pressure (MSLP) and maximum surface wind speed (MSWS) from the JMA and JTWC were signi cant, even at the initial time, the evolutionary trends of the MSLP and MSWS were nearly the same. e intensity evolutionary trends of the two typhoons in the CTL experiment were successfully recreated by the numerical simulation. Note that because the model initial condition was modi ed in the two sensitivity experiments, the model was unstable at the rst 24 h integration time.
erefore, the veri cation of the two sensitivity experiments using the TC bogus scheme (NoT and NoB experiments) is shown by the multiple levels at 1200 UTC 24 August ( Figure 5). e vortex of Tembin was removed completely at the sea level and 500 hPa isobaric level. However, there was still weak circulation of Bolaven at the sea level after applying the TC bogus scheme. e reason for that may be due to the size of Bolaven which was much bigger, and the scheme could not remove this size of vortex completely. All in all, this left-over vortex of Bolaven (actually, it was the remnant of monsoon which will be discussed in Section 3.3) in the NoB experiment had little impact on the Tembin since its intensity was much weaker.

Binary Typhoon Track.
e tracks of Tembin and Bolaven in the two sensitivity experiments are shown in Figure 3. It is remarkable that Tembin made a steady westward movement in the NoB experiment, which was completely di erent from the observations. Compared with the observations, Bolaven shifted westward in the NoT experiment after 0000 UTC 26 August, which was likely to make a greater in uence on eastern China. On the other hand, Tembin made an evident southward motion in the NoM experiment without the monsoon obstructing it to move south as compared with observations and the CTL experiment ( Figure 6(a)), and the orbit path for the TCs was more cyclonically curved ( Figure 6(b)). Meanwhile, the translation speed of Bolaven was greatly reduced in the NoM experiment ( Figure 6(a)). e southwesterly low-level jet (wind speed over 12 m s −1 in Figure 6(f )) to the south and southeast of the TCs was removed in the NoM experiment so that the circulations of the two TCs both became smaller. Also, the strength of the anticyclone circulation in the south of the TCs was weaker than that of the CTL (Figure 6(c)) and NoT experiments ( Figure 5(f )), and the subtropical high in the north of the TCs was divided into two parts with one locating over the Eastern China and the other locating over Japan ( Figure 6(d)). Under this synoptic situation, Bolaven moved northward much more slowly as the steering ow was not distinct. ese results suggest that the monsoon had an impact on both the tracks and the translation speeds of the TCs, especially for the smaller vortex of Tembin which was also a ected by both the existence of Bolaven. It is worth mentioning that the left-over vortex in the NoB experiment was the remnant of monsoon which was broad and belonged to the low-frequency ows and the TC bogus scheme could not completely remove it. Many previous researches have demonstrated that the environmental steering ow has best correlation with TC motion [21][22][23][24][25][26]. To access the di erences in the trajectories of the two typhoons in the CTL and the two sensitivity experiments, the environmental steering ow is calculated as the deep-layer mean (from 850 hPa to 300 hPa) wind averaged over a 5°-7°annulus around the TC's center, that is, the steering ow at each level is V(p) 7 5 2π 0 Vr dr dθ 7 5 2π 0 r dr dθ , (2) and the deep-layer-mean steering ow is where V is the meridional component V x or the zonal component V y of the environmental steering ow. As seen in Figure   A PV tendency diagnostic approach proposed by Wu and Wang [27] was applied to evaluate the TC motion in this study. Previous studies [28,29] have shown that the PV tendency was useful to diagnose the track of TC, especially the anomaly track, indicating that TC moves to the region where the PV tendency reaches maximum. e PV and the vorticity equation in the pressure coordinate can be expressed, respectively, as where P, g, θ, ζ, and f are the potential vorticity, gravitational acceleration, potential temperature, relative vorticity, and the Coriolis parameter, x, y, and p are the zonal, meridional, Substituting equations (5) into (6), the PV tendency (PVT) equation can be written as zP zt where F 1 represents the horizontal PV advection, F 3 represents the diabatic heating, F 4 represents the friction, Equation (7) suggests that the PV tendency results from horizontal (HA) and vertical (VA) PV advection, diabatic heating (DH), and friction. HA represents that the external large-scale circulation contributes to the PVT, and DH represents that the PV is a ected by the internal convective activity of TC. Since the VA is smaller than any other terms in the PVT equation and TC mainly moves over ocean, the contributions from the VA and friction are ignored in this study. e vertical integrated HA, DH, and PVT in the air column are calculated and used to evaluate the di erent tracks of Tembin in the CTL and NoB experiments. In the CTL experiment, the distribution of PVT was northeastsouthwest direction at 0600 UTC 25 August (Figure 8(c)) and after that time, the distribution turned to south-north direction. During this period, Tembin moved from westward to southward (black dot represents the current TC's position, and purple dot represents the next 6 h TC's position), which indicates that Tembin tended to move to the regions with maximum positive PVT and the change in Tembin's track was correctly re ected by the change in the PVT distribution. Note that the pattern of HA was consistent with that of PVT but the distribution of DH was asymmetric and the order of magnitude was smaller, suggesting that the contribution from the DH to the PVT was smaller than that of the HA. At 1800 UTC 26 August, the distributions of PVT, HA, and DH were both turned to east-west direction and Tembin moved eastward (Figures 8(d)-8(f )). After that time, the distribution of DH became symmetric (Figure 8(h)), and the environmental steering ow and the translation speed started to increase gradually (Figures 7(c) and 7(d). In the NoB experiment, the distributions of PVT and HA were nearly the same with the results from the CTL experiment before 0600 UTC 26

Advances in Meteorology
August (Figures 9(a) and 9(c)). However, the distributions of them turned to east-west direction after that time (Figures 9(d) and 9(f ), and Tembin moved westward steadily (Figure 3), which was distinctly different to the results from the CTL experiment (Figures 8(d) and 8(f )). During this period, the distribution of DH was asymmetric (Figures 9(b) and 9(e)), which restrained the westward movement of Tembin, making its translation speed decrease from 1200 UTC 25 August to 1200 UTC 26 August. e influence of DH on the track of Tembin reduced after 1800 UTC 26 August with the magnitude of DH decreasing sharply (Figure 9(h)). e southerly greatly reduced in the south of Bolaven in the absence of the monsoon (Figure 6(d)). When Bolaven moved poleward with increasing intensity, the formation of the south side of the subtropical high began to alter and its ridge line moved southward (Figure 6(c)), which was evidently observed in the NoT experiment ( Figure 5(f )). However, the ridge line of the subtropical high was a constant east-west direction in the NoB experiment ( Figure 5(e)) and the corresponding environmental steering flow of Tembin (Figures 7(c) and 7(d)) was northwestward, making Tembin move westward steadily.

Binary Typhoon
Intensity. e area of the water vapor associated with the southwesterly flow over the South China Sea in the CTL run was larger than that in the NoT run. It suggests that Tembin had an ability to drag considerable amount of water vapor from the South China Sea. And, there was air flow between the two typhoons in the CTL experiment ( Figure 10). Obviously, the water vapor carried by the southwesterly flow was divided into two parts, where one was absorbed by Tembin and the other became one of Bolaven's water vapor sources. However, Tembin was located at upstream of the southwesterly, which interrupted part of the water vapor being transported to Bolaven with few water vapor inflow as compared with that in the NoT run. Simultaneously, the distribution of the water vapor flux near Bolaven's center in the NoT experiment was more symmetrical than that in the CTL run. ese results suggest that the existence of Tembin was not conducive to the supply of water vapor to Bolaven. e simulation results of the intensity of Bolaven in the CTL and NoT experiment (the diamond line in Figures 4(b) and 4(d) showed that Bolaven's intensity was higher after 1200 UTC 24 August but lower after 0000 UTC 27 August. ese differences in the intensity of Bolaven indicate that the existence of Tembin seemed to have an impact on the intensity of Bolaven, which is discussed below.
Bolaven intensified simultaneously with the magnitude of wind in low troposphere, which was directly reflected by the increased kinetic energy. e kinetic energy budget is used to investigate the reason for the variations of the intensity of Bolaven in the CTL and NoT experiments. e kinetic energy budget can be expressed as [30] zK zt where K � (1/2)(u 2 + v 2 ) is the horizontal kinetic energy per unit mass of air parcel, V h is the horizontal wind vector, C tc is the translation vector of Bolaven, and F is the friction. Simplifying equation (12) results in zK zt e left-hand side represents the local rate of change in kinetic energy per unit mass of air parcel. e terms on the right-hand side represent the ways in which kinetic energy can be changed and include the generation of kinetic energy (G k ), horizontal flux of kinetic energy (HF k ), vertical flux of kinetic energy (VF k ), and kinetic energy dissipation (R) due to friction. e temporal evolution of kinetic energy budget terms in the CTL and NoT experiments shows that the local rate of change in kinetic energy was basically positive, indicating that the intensity of Bolaven increased gradually. And, the major contribution to the kinetic energy increase in low troposphere came from G k and to the kinetic energy dissipation came from VF k (not shown). e differences in the kinetic energy budget terms between the CTL and NoT experiments are represented in Figure 11. e magnitude of positive values of G k was bigger in the NoT experiment during the investigating period (Figure 11(a)), suggesting that the source of kinetic energy, that is, the available potential energy converted into kinetic energy by pressure gradient force, was greater than that in the CTL experiment. HF k reflects the kinetic energy exchange between Bolaven and the surrounding atmosphere. e significant differences in HF k between the two experiments were observed in low troposphere (below 850-hPa) (Figure 11(b)), which was related to the variations in the low-level jet. A considerable amount of energy was transported out of the area of Bolaven by HF k whose values were negative below 900 hPa in both experiments after 1200 UTC 24 August, and this kinetic energy output was higher in the CTL experiment from 1200 UTC 24 August to 0000 UTC 26 August. Above 900 hPa, HF k was positive in the NoT experiment during the entire period, but its magnitude was smaller after 0000 UTC 27 August. ese results suggest that the existence of Tembin tended to hamper the kinetic energy being transported to Bolaven by the low-level jet, which corresponded to the analyses of water vapor flux ( Figure 10). VF k represents the vertical transport of kinetic energy from low levels to upper levels, and negative values found in both experiments in the low layer (950-700 hPa) dominate throughout the period (not shown). It indicates that the transport of kinetic energy aloft by the ascending motion was the source of kinetic energy to the upper layer (100-300 hPa) with the positive values of VF k . e frictional dissipation effect in both experiments was produced below 850 hPa; however, the positive values of R were found at the layer between 850 hPa and 300 hPa during the entire time.
e positive values indicate that a transfer of energy from subgrid to grid scales of motion is a source of kinetic energy, which was reported by previous studies [31][32][33]. At certain moments, such as active convection and mutual stage of TC, the conversion of kinetic energy in a subgrid scale to that in the grid scale could be a source of kinetic energy and the values of R could be positive. e positive values of R observed in middle troposphere reflect the contribution from active convection to the increase in the kinetic energy at the upper level. Note that, as shown in Figure 11, the main differences in the kinetic energy budget terms were both occurred from 0000 UTC 26 August to 0000 UTC 27 August in the low layer (950-850 hPa). e contributions from the kinetic energy budget terms in the NoT experiment seem to be more beneficial to the development of Bolaven, for example, more kinetic energy was generated by G k in low troposphere (Figure 11(a)), and more kinetic energy was transported to the upper troposphere by ascending motion (Figure 11(c)). But, actually, the intensity of Bolaven was weaker than that in the CTL experiment except after 0000 UTC 26 August.
A series of observations and numerical simulations have shown that vertical wind shear has an important effect on TC Advances in Meteorology intensity [34][35][36][37][38]. e area-averaged (1000 km × 1000 km) horizontal winds over Bolaven and the area-averaged (1000 km × 1000 km) magnitude of deep-layer shear vector between 200 hPa and 850 hPa, as well as the low-layer shear vector between 850 hPa and 1000 hPa, were calculated to investigate the reason why the intensity of Bolaven was weaker in the NoT experiment before 0000 UTC 26 August. e results (Figure 12) show that the magnitude of the  averaged horizontal wind at lower levels were evidently greater in the NoT experiment than in the CTL run (Figures 12(a) and 12(b)). Although the difference in the wind direction at all levels in the two experiments was minor, the low-layer shear was stronger than that in the CTL experiment due to an increase in wind speed at low levels in the NoT experiment. Note that the low-layer shear and the deep-layer shear in the NoT experiment were both greater as compared with the results from the CTL experiment from 1200 UTC 24 August to 0000 UTC 26 August (Figures 12(c) and 12(d)). e upper circulation would be blown away from the low-level center of TC due to a stronger vertical wind shear, which is highly not conducive to the development of ascending motion and the intensity of TC. As shown in Figure 11(c), the magnitude of negative values of VF k in the middle and upper layers (800-300 hPa) was bigger in the CTL experiment than in the NoT experiment before 0000 UTC 26 August, indicating that the ascending motion was so strong that more kinetic energy was transported to the upper level, especially at 1200 UTC 25 August. After that time, the vertical wind shear weakened in the NoT experiment (Figures 12(c) and 12(d)) and the strength of ascending motion became stronger so that more kinetic energy was transported to the upper levels ( Figure 11c) and the intensity of Bolaven became stronger (Figures 4(b) and 4(d)).

Conclusions
e binary typhoon effect is a challenge in typhoon research because of its difficulty to predict. Previous studies on how the one-way interaction mode works is still unclear, such as the case of typhoons Morakot and Goni in 2009. In this report, a typical case of one-way interaction type, typhoons Tembin and Bolaven, was analyzed to investigate the binary typhoon effect on the track and intensity of the two typhoons by using the WRF model. e results are as follows: (1) e tracks of Tembin and Bolaven were diagnosed by the environmental steering flow, showing that there was a good relationship between the environmental steering flow and the tracks of the two TCs. And, the impact of the existence of Tembin on the track of Bolaven was negligible due to relatively slight differences in the environmental steering flow and translation speed of Bolaven. However, the results from the potential vorticity tendency demonstrated that Tembin tended to move to the region with maximum positive PVT and the main contributions to its looping track were the horizontal PV advection. In other words, the track of Tembin was affected by the external surrounding circulation associated with Bolaven. In addition, the monsoon over the western North Pacific also affected the tracks and the translation speeds of the TCs, especially for the smaller vortex of Tembin, which makes the tracks of TCs more complex. (2) e analysis of kinetic energy budget suggests that the development of Bolaven owed to the generation of kinetic energy in low troposphere and the vertical transport of kinetic energy to upper troposphere. e existence of Tembin not only interrupted the water vapor being transported to Bolaven but also affected the strength of vertical wind shear. Without the vortex of Tembin, the relatively bigger magnitude of the lowlayer and deep-layer wind shear was not conducive to the vertical transport of kinetic energy by ascending motion, resulting in the much weaker intensity of Bolaven during the looping track of Tembin.
is study demonstrates clearly the mechanisms of the one-way interaction mode that affects the track and intensity of binary typhoons. ese conclusions can also shed some light on the improvement of the TC prediction when the binary typhoon effect occurs. More work related to other typical categories in binary typhoon effect research could be done in the future, such as merging cases to illustrate the physical mechanisms of binary typhoon systems in more detail.

Data Availability
Most of the data series used to support the findings of this study are obtained from the Research Data Archive, which is available on the web at https://rda.ucar.edu/. e European Centre for Medium-Range Weather Forecasts (ECMWF) ERA-Interim reanalyses data are also available on the web at http://apps.ecmwf.int/datasets/. And, the best tracks of tropical cyclones are collected from the International Best Track Archive for Climate Stewardship, which is available on the web at https://climatedataguide.ucar.edu/climate-data/ ibtracs-tropical-cyclone-best-track-data.

Advances in Meteorology
Conflicts of Interest e authors declare that they have no conflicts of interest.