A Numerical Study on Impact of Taiwan Island Surface Heat Flux on Super Typhoon Haitang (2005)

Three to four tropical cyclones (TCs) by average usually impact Taiwan every year. This study, using the Developmental Tested Center (DTC) version of the Hurricane WRF (HWRF) model, examines the effects of Taiwan’s island surface heat fluxes on typhoon structure, intensity, track, and its rainfall over the island. The numerical simulation successfully reproduced the structure and intensity of super Typhoon Haitang. The model, especially, reproduced the looped path and landfall at nearly the right position. Sensitive experiments indicated that Taiwan’s surface heat fluxes have significant influence on the super TyphoonHaitang. Compared to sensible heat (SH) fluxes, latent heat (LH) is the dominant factor affecting the intensity and rainfall, but they showed opposite effects on intensity and rainfall. LH (SH) flux of Taiwan Island intensified (weakened) Typhoon Haitang’s intensity and structure by transferringmore energy from (to) surface. However, only LH played amajor role in the looped path before the landfall of the Typhoon Haitang.


Introduction
TCs are one of the most dangerous natural hazards to human society and the environment [1][2][3][4].Severe TCs produce destructive winds, high surges, torrential rains, and severe floods, usually resulting in serious property damage and loss of life [3,[5][6][7][8].Thus, it is of great interest to accurately predict TC's track and intensity.However, although the skill of position (track) forecasts has increased significantly, the accuracy of intensity forecasts has still to be improved [9][10][11][12][13][14]. Apparently, tropical cyclone intensity change interweaves multiscale nonlinear interactions among different physical processes, which are often poorly resolved or parameterized by global and even regional models [15][16][17].Such processes include surface fluxes of moisture and heat, moist convection, radiative transfer, and cloud microphysics [13,18].
Surface fluxes of moisture and heat are critical for mesoscale meteorological models [13,19].Modeling studies of sensitivity of TC development and intensity to sea/land surface conditions [20][21][22] were performed through alterations of the exchange coefficients of heat, moisture, and momentum at the surface.These studies revealed the relation between tropical cyclones development and maintenance of the model storm.The air-sea exchange of SH appears to be far less important.The studies of Tuleya and Kurihara [23] indicate that suppression of evaporation is the most important factor in the decay of a storm upon landfall.A reasonable change in surface temperature, moisture, and distribution of surface roughness can make a considerable difference in behavior [24].
Three to four tropical cyclones (TCs) by average usually impact Taiwan every year.The aforementioned studies emphasize role of surface conditions of land and ocean in affecting TCs, and the effect of surface heat fluxes of Taiwan Island is the issue that still needs to be addressed.In this study, we aimed to extend these previous studies by numerical simulations of Typhoon Haitang to explore the influence of surface heating of Taiwan Island on typhoon behavior.Super Typhoon Haitang was selected because of its cyclonic looping motion around Taiwan and strong intensity.
The track deflection of super typhoon Haitang was simulated in Jian and Wu [25] utilizing the advanced research version of the Weather Research and Forecasting model (WRF) with a 4 km fine mesh.It is found that the northerly advection flow results in a sharp southward turn of the westward-moving storm.The advection wind vectors rotate 2 Advances in Meteorology cyclonically in time and well match the motion of the simulated vortex during the looping process.Lowering the Taiwan terrain elevations to 70% or 40% of those in the control experiment reduces the southward track deflection and loop amplitude.Reducing elevation to 10% does not show a looping track and thus demonstrates the key role of the terrain-induced channeling effect.
Despite the work in Jian and Wu [25], the role of surface heat fluxes over Taiwan Island to Haitang remains to be investigated.In this study, some numerical experiments are designed to examine the track, intensity, and rainfall over Taiwan associated with surface heat fluxes over the island in two particular aspects.One is the role of heat fluxes of Taiwan Island in the intensity and structure of Haitang, and the other is its impact on the cyclonic looping motion of Super Typhoon Haitang before its landfall in Taiwan, which poses a rather unique scientific and forecasting problem [25].The rest of the paper is organized as follows.Section 2 describes the overview of Typhoon Haitang.Sections 3 and 4 show the model description and experiment design.The results of simulations and the detailed effects of heat fluxes on intensity, rainfall, and track were presented in Section 5. Section 6 gives summary.

Overview of Typhoon Haitang
Super typhoon Haitang developed from a poorly organized depression at UTC 10 July 2005.By 0000 UTC 12, it had reached tropical storm (17 m s −1 ) and upgraded to typhoon (32.7-41.4m s −1 ) at 1800 UTC 13.It turned westsouthwestward after 0600 UTC 12 and west-northwestward 3 days later heading toward Taiwan Island, largely steered by the easterlies associated with the subtropical high to its north [26].By 06 UTC July 16, it strengthened into a super typhoon (⩾51.0 m s −1 ) and became a threat to Taiwan Islands.Before making landfall in Taiwan at 0000 UTC 18, Haitang experienced sharp southward in its track and a loop to the east of Taiwan.

Model Description
The Hurricane Weather Research and Forecast System (HWRF) Model [27,28], an advanced hurricane prediction system based on Weather Research and Forecast software framework [29], was developed at NCEP's Environmental Modeling Center (EMC) for understanding the hurricane forecast issues and improvement of hurricane prediction.This HWRF modeling system, based on the Nonhydrostatic Mesoscale Model dynamic core [30,31], is a highresolution coupled air-sea-land prediction model with a movable nested grid capable of following the moving TC [27].The configuration of the HWRF system consists of outermost (parent) domain and two movable nested grids with resolutions of 27, 9, and 3 km, respectively, and 42 vertical levels.This HWRF modeling system has been used for real-time operational forecasts at NCEP since 2007 and was made accessible to research community in April 2010.The rotated latitude-longitude projection is utilized to maximize efficiency.The discretization grid is based on the Arakawa-E grid.The parameterization schemes are briefly described in references [27,28,32].HWRF uses a TC initialization and relocation algorithm to initialize the TC based on the observed position and intensity.The Gridpoint Statistical Interpolation (GSI), a 3-Dimensional Variational (3DVAR) Data Assimilation (DA) system, is implemented for initializing the environmental circulation outside the TC vortex area.The vortex initialization, relocation, and GSI DA produce the initial conditions to better match the actual TC observations.Therefore, a composite synthetic vortex, which is generated from bogus vortex or cycled from the vortex of the previous cycle, is assimilated in the observed TC location.The initialization of the HWRF Advances in Meteorology

Experimental Design
A series of 54 h integration experiments are carried out to explore the impact of Taiwan Island heat flux on Haitang.NCEP/NCAR FNL reanalysis [33]

Result
Figure 1 shows comparison of the best track and the simulated track (CTRL) of Typhoon Haitang.Haitang initially moves northwestward to Taiwan and then turns suddenly southward with a cyclonic looping path before landfall.It makes landfall around 0600 UTC 18 July and leaves at about 1800 UTC 18 July.The CTRL experiment reproduces reasonably the cyclonic looping path and the time of landing in and leaving Taiwan.Figure 2 compares the TC intensity between the control experiment and JTWC analysis.Such a comparison indicates that the simulated intensity agrees well with that from JTWC during most of the integration period, despite of an underestimated intensification rate for the 6∼18 h of simulated time.In general, control experiment shows good track simulations for Haitang.The entire simulation in CTRL shows the fairly accurate track and intensity evolutions.Figure 3 shows the horizontal distributions of reflectivity from the radar observation and the 3 km CTRL simulation.The model simulation reproduced the observed typhoon structure reasonably well, including eye, eyewall, and rainbands.The azimuthal phase of the simulated rainbands, especially, matches well that of the observation.Since the CTRL reproduced reasonably well Haitang's intensity, structure, and movement, the CTRL simulation is studied to understand the impact of Taiwan Island surface heat flux on the intensity and structure of and the looping motion of Haitang.
Taking an average over the period of model integration (Figure 4), the LH fluxes are positive over the Taiwan Island.In contrast, SH fluxes are negative over most of Taiwan, indicating an opposite direction of the surface energy transport between the typhoon and the boundary.The highest values of the LH flux are found in the east and west sides, while the lowest are in the central area.Similarly, distribution of the SH flux resembles that of the LH flux.To evaluate the impact of Taiwan surface heat fluxes on the typhoon, it is also meaningful to compare radar reflectivity from sensitivity simulations.Figure 6 reveals that all the simulations, including the CTRL, produce some rainbands from south of the storm to its north.It is important to note that the highest differences among the experiments are seen north of Taiwan.As radar reflectivity is a variable that is sensitive to the hydrometeor variables, whose structure and intensity changes are closely related to tropical cyclone intensity [34], the difference of radar reflectivity may be   7) is consistent with results above, that is, the LH (SH) trend to weaken (strengthen) wind at typhoon center (Figures 7(b) and 7(c)).Figure 8 shows the West-east cross sections of equivalent potential temperature at 1200 UTC 18 July 2005.The result from the CTRL experiment shows that an area of high   appeared from surface to 100 hPa over the typhoon center.Low   areas located at the west and east of the typhoon center.The gradient of the   is large between two low areas and the typhoon center, indicating strong upward movement in the area of the typhoon eyewall [14].Such distribution of   suggests the warm and wet structure of typhoon [5].Comparatively, No hflux shows larger gradient, as a result of the cutoff of SH fluxes from typhoon to the island surface, which is favorable for sustaining of large gradient   structure of typhoon.By contrast, without LH in No qflux (Figure 8(c)), the gradient is lower than that in CTRL.No flux (Figure 8(d)) also shows the same trend as No qflux.

Impact of Heat Fluxes on the Typhoon Intensity
By comparing the 54 h accumulated rainfall over Taiwan produced from the sensitivity experiments and the CTRL run (Figure 9), it is found that the transfer of both latent and sensible heat fluxes over Taiwan Island has an impact on the intensity of typhoon rainfall (Figures 9(b), 9(c), and  9(d)).Without the latent heat effect (Figure 9(c)), typhoon precipitation would be reduced significantly over north of the Taiwan Island but increased to the south.The area averaged precipitation decreased from 344 mm in CTRL to 306 mm in No qflux and increased to 347 mm in No hflux over Taiwan Island.SH show little impact on the island (Figure 10).From the sensitivity experiments, it is suggested that the LH flux plays a dominant role in intensity of typhoon rainfall.

Summary
In this study, a mesoscale numerical model particularly designed for the purpose of hurricane research and forecast, HWRF, is used to quantify the influence of Taiwan Island's heat flux on the track and structure of Typhoon Haitang that made landfall at Taiwan.The HWRF system consists of two movable nested grids with resolutions of 27, 9, and 3 km.CTRL experiment is performed with NCEP/NCAR FNL reanalysis.The experiment provided a reasonable replication of the track of typhoon Haitang during the 54 h integration, especially reproducing the looping motion east of Taiwan and the landfall location.The overall evolution of intensity also compared favorably with observations, although it underestimates intensification rate for the 6∼18 h of simulated time.Furthermore, radar reflectivity was very similar to what was observed.The model simulation reproduced the observed typhoon structure reasonably well, including eye, eyewall, and rainbands.The future work is intended to study several other loopedpaths around Taiwan and to summarize their common characteristics.Furthermore, we hope to reveal the roles of island surface heat flux played in the structure of TC.The goal of this research is to improve the understanding of and ability to predict anomalous track of typhoon around Taiwan Island.
Figure 5 shows MSLP of the sensitivity experiments and the control run.It shows that both the LH flux (Exp.No qflux) and the SH flux (Exp.No hflux) have an effect on typhoon intensity.The result in Exp.No qflux illustrates that LH flux from the Taiwan Island is favorable for typhoon maintenance.Without considering the effect on the LH flux of Taiwan, Typhoon Haitang would be weakened.This might be related to increased TC energy transported from the island surface and enhanced gain of energy when the typhoon made landfall over Taiwan.In contrast, Typhoon Haitang lost energy and attenuated during landfall.In short, Haitang is attenuated in No qflux but deepens in No hflux.Taking integration results until 54 h, the MSLP in Exp.No hflux (without SH) is decreased by 7.3 hPa while it increased by 12.8 and 10.8 hPa in Exp.No qflux and Exp.No flux, respectively.These results indicate that heating fluxes from Taiwan affect the development of super typhoon Haitang with opposite impact.
Figure 11   shows MSLP of the modeled TCs.Typhoon Haitang showed an asymmetric structure, when it is crossing Taiwan.High

Figure 9 :
Figure 9: (a) Horizontal maps of the 54 h accumulated rainfall (units: mm) over Taiwan from CTRL experiment during 0000 UTC 17 July to 0600 UTC 19 July 2005; (b) differences of the 54 h accumulated rainfall (units: mm) between No hflux and CTRL experiment; (c) and (d) are the same as (b) but for No qflux and No flux.