HTML
--> --> --> -->3.1. Numerical model
The numerical simulations in this study were conducted and performed with the WRF model (version 3.7; Skamarock and Klemp, 2008). Figure 2a shows the three nested domains with 9-, 3- and 1-km horizontal resolutions. All domains had 35 vertical (η) levels, with the finest resolution near the boundary layer and the model top at 50 hPa. The terrain datasets and land-use information were obtained from the U. S. Geological Survey with a 30-s temporal resolution. The WRF single-moment six-class microphysics (Hong et al., 2004), rapid radiative transfer model for longwave radiation (Mlawer and Clough, 1997), Dudhia shortwave radiation scheme (Dudhia, 1989), and Yonsei University planetary boundary layer scheme (Noh et al., 2003; Hong et al., 2006) were adopted in all the domains. In addition, the modified version of the Kain-Fritsch convective parameterization scheme (Kain, 2004) was applied in the outermost domain, but no cumulus scheme was used in the inner domains.2
3.2. Initial conditions and experimental design
There was independent integration in each domain because of the computational limitations: one-way interpolation was used between the domains of 9- (D01), 3- (D02) and 1-km grids (D03) (Fig. 2a). The bogus scheme was not used in each domain. Two experiments were designed as follows:Control experiment (CTRL): (i) in the 9-km D01, the initial and boundary conditions were provided by the 0.5° operational analyses of NCEP/GFS at 6-h intervals. Firstly, D01 ran freely from 0000 to 2000 UTC 6 August 2009, and then the horizontal wind retrieved from QuikSCAT data based on Ma and Tan (2010) were assimilated to obtain a new initial model field for the subsequent integration from 2000 UTC 6 to 1200 UTC 9 August. (ii) In the 3-km D02, the initial and boundary conditions were provided by the outputs of D01, and was run for 64 h from 2000 UTC 6 to 1200 UTC 9 August. (iii) The 1-km D03 was run for 42 h from 0600 UTC 7 to 0000 UTC 9 August, with the initial and boundary conditions from the outputs of D02.
Sensitivity experiment (FLAT): the experimental design was the same as CTRL, except that the overall terrain height on Taiwan Island was set to 1 m. The experiment adopted the methodology of (Yang et al., 2008) and (Tang et al., 2012) to deal with the vertical volume that was originally occupied by the terrain.
Note that all analyses and discussions in the following sections are based on the output of the 1-km grid (D03), except for the simulated tracks (Fig. 2b) and the calculation of steering flow from D01, and circulation centers and parcel trajectories from D02.
-->
4.1. Verification of CTRL
4.1.1. Track and intensityThe track of the simulated minimum sea level pressure (MSLP) center in CTRL is shown in Fig. 2b, as well as the best-track data issued by the CMA. In general, Morakot's track was reproduced well, especially regarding the two landfall sites——namely, on the eastern coast of Hualian County for the first landfall and on the northeastern coast of Fujian Province in China for the second landfall. A difference appeared over northwestern Taiwan in which the simulated storm presented a looping motion (Fig. 2b). This was responsible for the determination of the MSLP center, which may have been easily influenced by the downslope adiabatic warming and depression effect over the northwestern CMR, where SLs always appear. As for the typhoon intensity, the weakening process after landfall in Taiwan was also captured by the numerical simulation in the 1-km D03 (Fig. 2c).
4.1.2. Accumulative rainfall
As identified by (Wang et al., 2012), Morakot's key feature was its asymmetric structure, with an intense and long-lasting rainband over the southern Taiwan Strait, leading to the heavy rainfall in southern Taiwan. Therefore, good reproduction of this rainband is an important evaluation standard in terms of a successful numerical simulation. Figures 3a and b compare the 20-h accumulative rainfall simulated by CTRL and observations from Taiwan Central Weather Bureau (CWB) during Morakot's crossing over Taiwan from 1600 UTC 7 to 1200 UTC 8 August. In general, the model reproduced the heavy rainfall over southern Taiwan well, particularly with respect to the rainfall maxima observed at the border of Kaohsiung and Pingtung (near 23.0°N), and three poor rainfall areas in Taipei, Yilan and Changhua (Figs. 3a and b). However, the rainfall amount was a little overpredicted, e.g., the simulated rainfall maximum of 1338 mm versus the observed maximum of 1064 mm. Note that the simulated rainfall over the southern tip was greatly overestimated, as were the simulations reported by (Van Nguyen and Chen, 2011) and (Wang et al., 2012). Although this overestimated rainfall cannot be verified, due to the few rain gauges in these remote mountainous areas, the rainfall estimations from the Tropical Rainfall Measuring Mission satellite (Wang et al., 2012, Fig. 4) and operational radar in Taiwan (Figs. 4f-h) partly support that there should have been heavy rainfall over the southern tip. In general, there were two principal rain bands contributing to the heavy rainfall over southern Taiwan: (i) a high-reflectivity rainband, as displayed by Line1 in Fig. 4a, corresponding to a southwest-northeast-oriented rainfall area during 1200 UTC 7 to 0000 UTC 8 August (Fig. 3c); and (ii) Line2 in Figs. 4b and c, associated with a quasi-east-west-oriented rainfall area near 23.0°N during 0000-1200 UTC 8 August (Fig. 3d).
Figure3. (a, b) 20-h accumulative rainfall (units: mm) from 1600 UTC 7 to 1200 UTC 8 August (a) simulated by CTRL and (b) observed by Taiwan CWB; (c) 8-h rainfall from 1200 UTC 7 to 0000 UTC 8 August; and (d) 12-h rainfall from 0000 to 1200 UTC 8 August. The locations of the rain gauges at Yuyao shan (23.0°N, 120.7°E), Majia (22.7°N, 120.65°E) and Fangliao (22.4°N, 120.6°E) are marked in (a).
Figure4. (a-d) 850-hPa wind vectors (units: m s-1) and radar reflectivity (color shading; units: dBZ) simulated by CTRL and (e-h) corresponding radar reflectivity observations at 1200 UTC 7, 1800 UTC 7, 0000 UTC 8 and 1200 UTC 8 August. "Line1" and "Line2" denote two intense rain bands.
2
4.2. Evolution of Morakot's structure
Figure 5 shows the streamline field at 925, 700 and 500 hPa during Morakot's crossing of Taiwan. Prior to landfall, the storm centers were separated in the vertical direction. When the low-level center (below 4 km, lower than the highest altitude of the CMR) approached the central eastern coast of Taiwan, the mid-level center (above 4 km, denoted at 500 hPa herein) was located more than 100 km toward the southwestern coast (Figs. 5a and b). The dislocation between the low-level center and the mid-level center was associated with the strong southwesterly wind (Hong et al., 2010; Van Nguyen and Chen, 2011). The steering flows, defined by the area-averaged wind inside an annulus between 300 and 700 km from the simulated MSLP centers at 850, 700 and 500 hPa, showed that the low-level wind (850 and 700 hPa) changed from southeasterly to southwesterly, while a quasi-easterly and gradually weakening wind prevailed at 500 hPa during 2000 UTC 6 to 0600 UTC 7 August (Fig. 6a), which closely agrees with in the findings of Chien and Kuo (2011, see their Fig. 11). As a result, the initially overlapped storm centers had completely separated in the vertical direction by 0600 UTC 7 August (Fig. 6b). Similar to the methods used in (Jian et al., 2006) and (Van Sang et al., 2008), the circulation centers of the simulated storm in Fig. 6b were defined by the weighted-average VCs within a radius of 50 km from the streamline centers at 850 hPa and 500 hPa. Moreover, a strong northeasterly wind shear, derived from the steering flow at 500 hPa and 850 hPa (Fig. 6a), made the mid-level center of Morakot above the CMR tilt to the southwest relative to the low-level center (Fig. 5a and Fig. 6b), similar to that in previous studies (Wang and Holland, 1996; Reasor et al., 2004, Wu et al., 2011). Consequently, when the low-level center made landfall on the central eastern coast of Taiwan, the mid-level center was far away the southeastern coast (Fig. 5a). Then, the low-level center was filled over the eastern coast of Taiwan (Figs. 5a-c and 6b), whereas the low-level northerly wind of the storm flowed around the northern tip of the CMR and met the southwesterly monsoon to result in a confluent flow over the southern Taiwan Strait (denoted by CZ1 in Figs. 5a and b), corresponding to Line1 in Figs. 4a and b. The backward trajectories of air parcels released within CZ1 at 500 m above ground level (AGL), calculated by the Flexible Particle Dispersion Model with WRF (Stohl et al., 2005), demonstrate that this confluent flow was caused by Morakot's northerly wind and the southwesterly monsoon (Fig. 7a). By about 1400 UTC 7 August, a VC had emerged from the northern edge of CZ1, when the confluent flow was blocked by the southern CMR (Fig. 5b). After that, the VC moved northward along the western CMR and approached central western Taiwan (Fig. 6b), before finally coupling with the primary mid-level center of the storm passing over the CMR to reestablish a new and vertically stacked typhoon at 0000 UTC 8 August (Figs. 5c-e).Figure5. Time series of simulated streamline at 925, 700 and 500 hPa in CTRL at (a) 1000 UTC 7, (b) 1400 UTC 7, (c) 1600 UTC 7, (d) 2000 UTC 7, (e) 0000 UTC 8, and (f) 0600 UTC 8 August. The black stars denote the primary centers; the solid (dashed) line indicates the vertical extent of the typhoon (secondary) center; CZ1 denotes a convergence zone (or confluent flow) between the northerly wind of the storm and the southwesterly wind; and VC indicates the location of the vorticity center.
Figure6. (a) Environmental steering flow (units: m s-1) simulated in D01 by CTRL, defined by the area-averaged wind inside an annulus between 300 and 700 km from the circulation centers at three levels (850, 700 and 500 hPa) from 2000 UTC 6 to 0600 UTC 7 August. (b) Tracks of hourly circulation centers at 850 hPa (line with circles) and 500 hPa (line with squares) in CTRL from 2000 UTC 6 to 1200 UTC 8 August, with triangles denoting the tracks of the SCs at 850 hPa. The circulation centers are defined by the weighted-average vertical VC within a 50-km radius away from the streamline center, and the terrain height is indicated by shading with an interval of 500 m.
Figure7. The 8-h (a) backward and (b) forward trajectories (time interval: 1 h) for 10 air parcels inside D02 in CTRL at 1400 UTC 7 August. The positions of parcels were chosen at 500 m AGL within CZ1 (displayed in Fig. 5b).
2
4.3. Formation and development of the SC
4.3.1. FormationThe evolution of the simulated MSLP and 10-m wind vectors are shown in Fig. 8. Figures 8a-c show that the typhoon center was gradually filled by the land after landfall at Hualien, and two thermodynamic SLs meanwhile emerged on the leeside——namely, SL1 in the northwest, relative to the easterly wind passing over the northern CMR; and SL2 at the southeastern coast, relative to the westerly wind passing over the southern CMR (Figs. 8a-c)——consistent with the analysis of NCEP/GFS (Fig. 1a), which were caused by the downslope adiabatic warming and depression effect. In addition, a dynamic VC was produced over southwestern Taiwan, characterized by a low-pressure center and closed isobars (Figs. 8a-c). The VC was produced due to the interaction between the confluent flow and the CMR. Because the VC was initiated on the windslope relative to the westerly wind, there was an absence of a warm core inside it until 1500 UTC 7 August (Fig. 9a). Therefore, the VC was not a well-defined SC in this formation stage.
Figure8. Simulated MSLP (contour interval: 2 hPa) and 10-m wind vectors (units: m s-1) in CTRL at (a) 1400 UTC 7, (b) 1500 UTC 7, (c) 1600 UTC 7, (d) 2000 UTC 7, and (e) 0000 UTC 8 August. "TC", "SLn (n=1,2,3)," and "VC" denote the primary typhoon center, secondary lows, and vorticity center, respectively. Lines SN and EW (S'N' and E'W') cross the VC at 1500 UTC 7 August (0000 UTC 8 August).
Figure9. Simulated horizontal wind vectors (units: m s-1) and temperature (color shading; units: °C) at the lowest model level, superposed with the terrain height (contour interval: 500 m) at (a) 1500 UTC and (b) 2100 UTC 7 August. The black and green lines denote the hourly backward trajectories released at 500 m AGL inside the SC.
To highlight the formation process, the hourly horizontal wind vectors and absolute vorticity field at the lowest model level are given in Fig. 10. At 1200 UTC 7 August, a confluent flow appeared between the northerly wind of the storm and the southwesterly wind, resulting in a quasi-westerly wind near 22.8°N (Fig. 10a). When this westerly wind was blocked by the southern CMR, a part of the airflow deflected northward to become southerly wind along the western slope (Figs. 10b and c). This southerly wind then deflected westward again when it encountered the wider and higher orographic barriers near 23.2°N, before finally merging with the northerly wind to produce a closed circulation (Fig. 10d). The forward trajectories of air parcels released within the confluent flow confirm this closed circulation (Fig. 7b). Figure 11 shows vertical cross sections of temperature and horizontal wind crossing the SC (along lines EW and SN in Fig. 8b) at 1500 UTC 7 August. The results indicate that a strong westerly flow was located near 22.7°N (to the south of the VC). As a result, a positive shear vorticity -? u/? y>0 would have been produced at the northern edge of the strong westerly flow (Fig. 11a). Moreover, when the westerly wind was blocked by the southern CMR to result in southerly flow along the western slope of the CMR, another component of positive shear vorticity ? v/? x>0 could be produced (Fig. 11b). Consequently, there was considerable positive vertical vorticity ? v/? x-? u/? y>0 over southwestern Taiwan, even though it looked partially asymmetric at this time. The vertical cross sections of zonal wind and meridional wind at line AA' (displayed in Fig. 10d) show that there was reversed flow (easterly wind) to the north of the VC (Fig. 12), which also reconfirmed that closed circulation appeared over southwestern Taiwan at 1500 UTC 7 August.
Figure10. Simulated horizontal wind vectors (full barb is 10 m s-1) and absolute vorticity (shading; units: 10-5 s-1) at the lowest model level from 1200 UTC to 15 UTC 7 August. The blue contours denote the terrain height (contour interval: 500 m), and green line AA’ is located at 23.2°N.
Figure11. Vertical cross sections of temperature (shading; units: °C) and (a) zonal wind (contour interval: 4 m s-1) along line SN, and (b) meridional wind along line EW (displayed in Fig. 8b) at 1500 UTC 7 August; (c) zonal wind along line S'N', and (d) meridional wind along line E'W' (displayed in Fig. 8e) at 0000 UTC 8 August. The black dots below the horizontal coordinate denote the location of the SC.
Figure12. Vertical cross sections of zonal wind vectors (shading; <0 m s-1) and meridional wind (contour interval: 4 m s-1) along line AA' (displayed in Fig. 10d) in CTRL at (a) 1300 UTC 7 and (b) 1500 UTC 7 August.
4.3.2. Development
Although the VC featured closed circulation with high vertical vorticity at 1500 UTC 7 August, the absence of a warm core and symmetric structure meant that it was still not a well-defined SC at this time. As the VC drifted northward along the western slope of the CMR, due to the steering of the southerly flow and complete departure from the confluent flow, its asymmetry gradually disappeared, ultimately becoming rather symmetric in structure over central western Taiwan at 0000 UTC 8 August (Figs. 11c and d). Similarly, the temperature inside the VC gradually increased as it moved northward during 1400 UTC 7 to 0000 UTC 8 August (Fig. 13a). Consequently, the VC has a warm core (Fig. 13b). The backward trajectories show that the warm air, produced by the adiabatic warming on the downslope of the northern CMR during the passage of the easterly wind across northern Taiwan, was continuously advected into the VC (Fig. 9b), contributing to the increase in temperature inside the VC (Fig. 13). So, by 0000 UTC 8 August, the VC had developed to become a well-defined SC with a near-symmetric, closed-circulation and warm-core structure.
Figure13. (a) Time-height (AGL) section of the area-averaged temperature (shading; units: °C) within 50 km from the circulation center of the SC (shown in Fig. 6b) and (b) the evolution of surface temperature difference between being averaged within 50 km and inside an annulus between 50 and 100 km from the circulation center of the SC in CTRL.
2
4.4. Orographic effects
A sensitivity experiment with the CMR removed was performed to examine the effect of orography on the structural evolution of Morakot during its crossing of the CMR. As shown by Fig. 2b, prior to landfall, the simulated track in FLAT looked similar to that in CTRL, apart from a more northward track closely before landfall because of the absence of topographic channel impact, as proposed by Yeh and Elsberry (1993a, 1993b). After landfall, FLAT also simulated a more northward track than that observed and in CTRL. This is because the lower mountain height led to a more northward steering flow and track, in agreement with previous studies (Yang et al., 2008; Liang et al., 2011; Wang et al., 2012).Figure 14 shows the streamline field of the simulated storm in FLAT, every 6 h from 1200 UTC 7 to 1200 UTC 8 August. Similar to the simulation in CTRL, the storm centers also inclined southward with altitude prior to landfall, and the low-level northerly wind of the storm also coupled with the southwesterly wind to result in a confluent flow over the southern Taiwan Strait, but no SC formed over southwestern Taiwan (Fig. 14a), because of the absence of closed circulation due to the orographic blocking.
Figure14. As in Fig. 5 but for FLAT every 6 h from 1200 UTC 7 to 12 UTC 8 August.
2
4.5. Discussion on the evolution of the SC
Based on the above analysis, the evolution of the SC inside Morakot can be divided into two stages: formation and development. As the primary center moved toward the eastern coast of Taiwan, its cyclonic wind had flowed around the northern tip of Taiwan and met the strong southwesterly wind to result in a confluent flow (CZ1 in Fig. 5b). The westerly wind within the CZ1 would deflect northward to induce a southerly wind along the western slope of the CMR, due to being blocked by the CMR. As a result, a VC with high vertical vorticity and closed circulation was produced over southwestern Taiwan (Fig. 10d and Figs. 11a and b). Because the VC initially emerged in the relatively cold area on the windslope relative to the westerly wind, it had no warm core in the initial formation stage (Fig. 9a). As the VC departed from its formation region and drifted northward along the western coast of Taiwan, similar to the findings of (Lee et al., 2008) and Lin et al. (1999, 2006), warm air produced within the wake flow of an easterly wind passing over the northern CMR was continuously transported into the VC to produce a warm core (Fig. 9b). Consequently, the VC became a well-defined SC with a warm-core and near-symmetric structure. In summary, although the SC inside Morakot was initially produced by a dynamic process, its development was closely associated with the thermodynamic process documented by previous studies.But what caused the uniqueness of the formation of the SC inside Morakot? Comparing the wind fields of typhoons Morakot (2009), Mindulle (2004) and Dot (1990), the largest difference is that for Morakot there was a quasi-east-west-oriented, strong and long-lasting confluent flow sustained over the southern Taiwan Strait (Fig. 10; Jian et al., 2006, Fig. 5; Lee et al., 2008, Fig. 12), which was likely a very important factor in producing the VC inside. Moreover, as identified by (Chien and Kuo, 2011), the maintenance of the confluent flow was because of two unique features of Morakot: the strong southwesterly flow and the slow translation speed of the TC. Firstly, the southwesterly wind was so strong that it encountered the northerly wind of the storm to result in a confluent flow over the southern Taiwan Strait prior to landfall. Secondly, the translation of Morakot crossing Taiwan was sufficiently slow that there was abundant time for interaction between the confluent flow and the southern CMR to produce a VC.