Interactive comment on “ Typhoon effect on Kuroshio and Green Island wake : a modelling study ”

General Comments The paper is a study of the effect of one real typhoon and four simulated typhoons on the ocean currents near Green Island. The study region lies in the path of the Kuroshio with depths extending down to 5000m. The study makes use of a finite element model based on the shallow water equations. It reads as a reasonable paper but the description of the model is poor and the paper make a basic error in using a barotropic model for a region of ocean where baroclinic effects are likely to be important. The paper is not acceptable in it present form. More importantly the use of a barotropic model means that even if other parts of the paper were considerably improved the C1713


Introduction
According to Zdravkovich [1], during flow past a bluff obstacle with the Reynolds number Re = u ∞ L/v t in the range of 50 and 800, a well-organized downstream wake occurs, where u ∞ is the characteristic flow velocity, L the characteristic length, and v t the eddy viscosity of fluid.Flow becomes periodic with the detachment of the free shear layers and consequent alternate shedding of vortices.The periodic phenomenon is referred to as the vortex shedding, where the anti-symmetric clockwise and counterclockwise wakes pattern is called the von Karman vortex streets [1].
The Kuroshio, a western boundary current of the sub-tropic North Pacific Ocean, originates from the North Equatorial Current and flows northward to the eastern coast of Taiwan.The passage of the Kuroshio mainstream is parallel to the eastern shoreline of Taiwan.Green Island is located at (121 • 28 E, 22 • 35 N) and is 40 km off the southeastern coast of Taiwan.The climatology of the Kuroshio velocity as revealed from in-situ measurements shows that Green Island is approximately located in the mainstream of the energetic Kuroshio passage and acts as an obstacle to the stream.Hence, downstream island wakes are prone to occur.Characteristics of the vortex streets in downstream Green Island can be found from the moderate resolution imaging spectroradiometer (MODIS) satellite images and in-situ measurements [17].The seasonal variation of mainstream current patterns, and the transport of the Kuroshio in the east of Taiwan, has been numerically investigated [18].Spatial and temporal scales of downstream Green Island wakes due to passing of the Kuroshio have been numerically studied [19].The seasonal wind-, heat flux-, and buoyancy-driven circulation and thermohaline structure of the Caspian Sea is numerically investigated [20].Salient features of the mesoscale dynamics, mixing and transport of the Caspian Sea are well reproduced, and compared well with the observation.
Typhoons, also known as tropical low-pressure cyclones, occur in the tropical and subtropical seas, and are a frequent cause of serious disasters in coastal countries and regions.Taiwan is located in the Northwest Pacific-one of the most typhoon-active areas of the world.An annual average of four to five typhoons hit the island according to the records of the Central Weather Bureau (CWB) of Taiwan (http://www.cwb.gov.tw/) over the past century .Typhoons typically occur in the period from late April to December, but with most storms concentrated in July through September.In recent years, the number of typhoons and the maximum wind speed have increased and the central pressure of typhoons have decreased due to the climate change.
Typhoons have a significant impact on the physical and chemical characteristics of the water of the Kuroshio.The storms disturb the Kuroshio's path and reduce its flow rate, while strengthening upwelling and lowering the surface temperature.The impact of typhoons on the upper water layers is quite dramatic in the Kuroshio.Suda [21] suggested that typhoons are the main cause for the Kuroshio meandering near Japan.Sun et al. [22] conducted a numerical analysis for the impact of typhoons on the Kuroshio's path south of Japan, and found that violent disturbances on the sea surface caused by a typhoon produces a strong upwelling, further strengthening the counterclockwise eddy and causing the Kuroshio to deviate southward from its original path by about 2 • of latitude.Chen et al. [23] suggested that, in the sea of eastern Taiwan, typhoons with a diameter exceeding 200 km (a large counterclockwise eddy) can reduce significantly the Kuroshio's flow volume and cause a 180 • change of the Kuroshio flow direction.
Taiwan has an excellent marine current energy resource as it is an island country and experiences abundant marine current flows.However, this resource has yet to be explored.The Kuroshio is known for its strong and steady flow.It could be a potential source of renewable energy as it flows steadily all year round, and Green Island is the potential plant site for Kuroshio current energy in Taiwan.Chen [24] proposed a conceptual design for the Kuroshio power plant.The assessment of the potential Kuroshio power test site has been performed using the in-situ measurements and modeling by Hsu et al. [25].Safety, maintenance and operations of the Kuroshio power plant may be subjected to impacts from earthquakes, typhoons, climate change, and other natural factors.To prevent the potential damages caused by typhoon-driven waves, meandering vortex shedding, and the cavitation occurring on the surface of turbine blades, turbines and platforms should be submerged several tens of meters below the water surface.The effect of monsoon on the Green Island wakes has been previously studied [25,26].The effect of typhoon on the Kuroshio and Green Island wakes is presented using a high resolution (250-2250 m) shallow water equations (SWEs) model in this study.

Theory and Numerical Method
Shallow-water equations (SWEs) for describing shallow water flows and the space-time least-squares finite-element method (STLSFEM) used to solve SWEs are briefly introduced in this section.

Shallow-Water Equations
The two-dimensional vertical-averaged SWEs are a simple form for describing the horizontal structure of the ocean dynamics [27][28][29][30].SWEs can be expressed in the conservative or non-conservative form.Comparison of the conservative and non-conservative form of SWEs can be found in Toro [31,32].The two-dimensional viscous SWEs in a non-conservative form read where η is the water surface elevation, t is time, H is the still water depth, h = η + H is the total water depth, u is the x-component velocity, v is the y-component velocity, Z b is the bottom height, and v t is the eddy viscosity which is an important parameter for modeling the Reynolds stresses for the fluid [33].τ b x and τ b y representing the xand y-component bottom frictional shear stresses are expressed as where C Db is the coefficient of bottom friction ranging between 0.0025 to 0.0040 for ocean modeling [34].
Holland's wind field model [35] is employed to compute the typhoon wind field.The equation of the gradient of wind speed is where B is a scaling number, P n is the ambient pressure, P c is the pressure of typhoon center, R mw is the radius of maximum wind speed, r is the radius to typhoon center, and ρ air = 1.2 kg/m 3 is the air density.
The wind shear stresses on the water surface are usually expressed in terms of the wind speed at 10 m above the water surface W 10 .Powell [36] suggested that the relationship between the gradient of wind speed W g and W 10 can be expressed as W 10 = 0.8W g .Therefore, τ s x and τ s y representing the xand y-component wind shear stresses are given as where W 10x and W 10y are the xand y-component of W 10 , C Ds is the coefficient of surface stresses [37].

Space-Time Least-Squares Finite-Element Method
We use the two-dimensional SWEs to illustrate the space-time least-squares finite-element method (STLSFEM).The two-dimension SWEs read where the source terms of the momentum equations are The Newton method is applied to linearize the nonlinear terms in Equation (7), and the resulting equations are where the symbol "~" denotes the value from the previous iteration or time step.In the STLSFEM, the unknowns u = {η, u, v} T are approximated by the polynomial interpolations where M(x) and N(t) are the space and time interpolation functions, respectively.Substituting the approximations Equation ( 9) into Equation (10), the residuals can be written as Note that a piecewise continuous linear interpolation function for time is used, where the subscripts 1 and 2 denote the number of temporal local nodes, and superscripts n and n + 1 denote the values of the space-time element at t = t n and t = t n+1 , respectively.
Upon applying the least square method, we thus have min where are the residuals, Ω xt is the element of the space-time.Equation ( 11) can be rewritten in the form Details of the STLSFEM can be found in Gunzburger [38], Hughes and Hulbert [39], Jiang [40], Laible and Pinder [41], Liang and Hsu [42].

Computed Results and Discussion
In this section, the effect of typhoon Soulik and wind field of Holland's model [35] on the Kuroshio and Green Island wakes is presented.

Typhoon Soulik
A high resolution (250-2250 m) SWEs model is used to investigate the hydrodynamic characteristics of the Kuroshio and downstream Green Island wakes.Green Island is located at (121 • 28 E, 22 • 350 N), 40 km off the southeastern coast of Taiwan.Figure 1a depicts the study area, the location of Lanyu Island and Green Island, and the bathymetry of the adjacent sea waters.The bathymetry of the study area varies from 400 to 5000 m, with an average value of about 2000 m.The bathymetry around the Green Island is within the range from 500 to 4000 m, with an average value of 2000 m.There are sea-mountains in the northwest of the Green Island where water depth varies from 200 to 400 m.

Computed Results and Discussion
In this section, the effect of typhoon Soulik and wind field of Holland's model [35] on the Kuroshio and Green Island wakes is presented.

Typhoon Soulik
A high resolution (250-2250 m) SWEs model is used to investigate the hydrodynamic characteristics of the Kuroshio and downstream Green Island wakes.Green Island is located at (121°28′ E, 22°350′ N), 40 km off the southeastern coast of Taiwan.Figure 1a depicts the study area, the location of Lanyu Island and Green Island, and the bathymetry of the adjacent sea waters.The bathymetry of the study area varies from 400 to 5000 m, with an average value of about 2000 m.The bathymetry around the Green Island is within the range from 500 to 4000 m, with an average value of 2000 m.There are sea-mountains in the northwest of the Green Island where water depth varies from 200 to 400 m.However, Kuroshio is a sub-surface flow where flow mainly occurs at the top 400 m layer of water with an average current speed around 1.0 m/s; water below 800 m is essentially motionless [17].Therefore, we limit the bathymetry ranging from 10 to 360 m as the redefined bathymetry from the ETOPO1 (1 Arc-Minute Global Relief Model) in the simulations.Figure 1b  The flow field data of HYCOM (hybrid coordinate ocean model) on 7 November 2014 after interpolation to the model grids, shown in Figure 1b, is utilized to initiate the flows.The choice of this particular flow field as the initial condition is based on the hypothesis proposed by Zheng and Zheng [43] in which the downstream Green Island wakes are prone to occur when Kuroshio mainstream heads on the island.Figure 2 shows the boundaries of the study domain.AB is the land boundary where no flux and free-slip boundary condition is applied, BC is the open boundary where the radiation boundary conditions are specified, and CD, DO, and OA are the velocity boundaries using the flow field data from HYCOM on 7 November 2014.
Zheng [43] in which the downstream Green Island wakes are prone to occur when Kuroshio mainstream heads on the island.Figure 2 shows the boundaries of the study domain.AB is the land boundary where no flux and free-slip boundary condition is applied, BC is the open boundary where the radiation boundary conditions are specified, and CD, DO , and OA are the velocity boundaries using the flow field data from HYCOM on 7 November 2014. Figure 3 shows a composite ERS-1 SAR (European Remote Sensing satellite1, Synthetic Aperture Radar) image of downstream island vortex trains of Lanyu Island and Green Island taken with one week interval sequentially between 17:01 UT 25 September (lower) and 19:02 UT 2 October (top), 1996.Two downstream island vortex trains are clearly seen.One a small distance from the coast, called the Lanyu Island (also called the Orchid Island) vortex train, occurs downstream of Lanyu Island (a little beyond the southern boundary of the image), and consists of three vortexes, two cyclonic appearing as bright shading and one anticyclonic appearing as dark shading vortexes.Another near the eastern coast of Taiwan, called the Green Island vortex train, occurs downstream of Green Island, and consists of three pairs of cyclonic vortexes (bright center) and anticyclonic vortexes (dark center).Spatial and temporal scales of downstream Green Island wakes due to passing of the Kuroshio has been numerically studied [19].The main difference between Liang et al. [19] and the present study lies in a constant inflow used in Liang et al. is replaced by a spatial varying inflow from HYCOM in the present study.Two downstream island vortex trains are clearly seen.One a small distance from the coast, called the Lanyu Island (also called the Orchid Island) vortex train, occurs downstream of Lanyu Island (a little beyond the southern boundary of the image), and consists of three vortexes, two cyclonic appearing as bright shading and one anticyclonic appearing as dark shading vortexes.Another near the eastern coast of Taiwan, called the Green Island vortex train, occurs downstream of Green Island, and consists of three pairs of cyclonic vortexes (bright center) and anticyclonic vortexes (dark center).Spatial and temporal scales of downstream Green Island wakes due to passing of the Kuroshio has been numerically studied [19].The main difference between Liang et al. [19] and the present study lies in a constant inflow used in Liang et al. is replaced by a spatial varying inflow from HYCOM in the present study.
The SWEs model is applied to simulation the typhoon Soulik event.Typhoon Soulik hit the Taiwan area during 7 to 14 July 2013.Figure 4 depicts the track of typhoon Soulik.Typhoon Soulik, developed from a tropical depression on 7 July 2013 through a tropical storm, moderate typhoon, severe typhoon, and then decreased its strength and made landfall on Taiwan on 7 December 2013.The typhoon Soulik data with temporal intervals of every 6 h are from the CWB of Taiwan, which are computed by the NFS (non-hydrostatic forecast system).NFS uses the structured orthogonal meshes; however, the SWEs model uses the unstructured triangular meshes.Therefore, wind data needs to be interpolated into the computational nodes of the SWEs model.The IDW (inverse distance weighted interpolation [44]) is employed to interpolate the data from NFS into nodal points of the SWEs model.In this study, the 5-km resolution typhoon Soulik data is used.The maximum wind speed of typhoon Soulik is 51 m/s, categorized as a severe typhoon.The closest distance of typhoon to Green Island is about 250 km.The maximum wind speed of the model area is 17 m/s and the average value is smaller than 10 m/s.The SWEs model is applied to simulation the typhoon Soulik event.Typhoon Soulik hit the Taiwan area during 7 to 14 July 2013.Figure 4 depicts the track of typhoon Soulik.Typhoon Soulik, developed from a tropical depression on 7 July 2013 through a tropical storm, moderate typhoon, severe typhoon, and then decreased its strength and made landfall on Taiwan on 7 December 2013.The typhoon Soulik data with temporal intervals of every 6 h are from the CWB of Taiwan, which are computed by the NFS (non-hydrostatic forecast system).NFS uses the structured orthogonal meshes; however, the SWEs model uses the unstructured triangular meshes.Therefore, wind data needs to be interpolated into the computational nodes of the SWEs model.The IDW (inverse distance weighted interpolation [44]) is employed to interpolate the data from NFS into nodal points of the SWEs model.In this study, the 5-km resolution typhoon Soulik data is used.The maximum wind speed of typhoon Soulik is 51 m/s, categorized as a severe typhoon.The closest distance of typhoon to Green Island is about 250 km.The maximum wind speed of the model area is 17 m/s and the average value is smaller than 10 m/s.
Figure 5a depicts the global view of the study area, where the enclosed region near the Green Island indicates the sub-domain that will be presented in the succeeding plots of the downstream recirculation and island wakes.Figure 5a depicts the global view of the study area, where the enclosed region near the Green Island indicates the sub-domain that will be presented in the succeeding plots of the downstream recirculation and island wakes.Figure 5b depicts the local view of the sub-domain as well as the location of the cross section Q 0 Q 1 and monitor point P at downstream Green Island, where time series of flow quantities are presented and analyzed.Length of Q 0 Q 1 is 20 km.It is located about 22 km downstream the Green Island.
Figure 6 depicts the typhoon wind field and flow streamlines as well as vectors of wind (red) and flow (blue) along cross section Q 0 Q 1 at S 0 -S 4 five time instances with a 12 h of increment, shown in Figure 5, where S 0 represents 00:00 of 12 July, S 1 12:00 of 12 July, S 2 00:00 of 13 July, S 3 12:00 of 13 July, and S 4 00:00 of 14 July 2013, respectively.As can be seen, the downstream island wakes are well reproduced.The flow field is affected by the typhoon, especially in the shallow water area and the lee of the islands.Interactions of Kuroshio and downstream island wakes with typhoon are obvious.When typhoon approaches the study area, because of the cyclonic rotation and southward winding of typhoon, Kuroshio current decreases.At S 3 , we can see the moment when the center of the typhoon is closest to the Green Island.The wind starts to change its direction and to blow northward.Comparing flow vectors along cross section Q 0 Q 1 at S 2 with that at S 1 and S 3 of Figure 6, it is noticed that downstream island wakes are significantly affected.Kuroshio currents increase after S 2 due to the northward winding and resumes its normal flow condition after typhoon passes far away.5, where S0 represents 00:00 of 12 July, S1 12:00 of 12 July, S2 00:00 of 13 July, S3 12:00 of 13 July, and S4 00:00 of 14 July, 2013, respectively.As can be seen, the downstream island wakes are well reproduced.The flow field is affected by the typhoon, especially in the shallow water area and the lee of the islands.Interactions of Kuroshio and downstream island wakes with typhoon are obvious.5, where S0 represents 00:00 of 12 July, S1 12:00 of 12 July, S2 00:00 of 13 July, S3 12:00 of 13 July, and S4 00:00 of 14 July, 2013, respectively.As can be seen, the downstream island wakes are well reproduced.The flow field is affected by the typhoon, especially in the shallow water area and the lee of the islands.Interactions of Kuroshio and downstream island wakes with typhoon are obvious.5, where S0 represents 00:00 of 12 July, S1 12:00 of 12 July, S2 00:00 of 13 July, S3 12:00 of 13 July, and S4 00:00 of 14 July, 2013, respectively.As can be seen, the downstream island wakes are well reproduced.The flow field is affected by the typhoon, especially in the shallow water area and the lee of the islands.Interactions of Kuroshio and downstream island wakes with typhoon are obvious.
a severe typhoon, 51 m/s < Vmax, and S 0 -S 4 are the five time instances with a 12 h of increment to show the flow fields, see Figure 6, for example, S 0 represents 00:00 of 12 July, S 1 12:00 of 12 July, S 2 00:00 of 13 July, S 3 12:00 of 13 July and S 4 00:00 of 14 July 2013, respectively.(Modified from CWB of Taiwan).5, where S0 represents 00:00 of 12 July, S1 12:00 of 12 July, S2 00:00 of 13 July, S3 12:00 of 13 July, and S4 00:00 of 14 July, 2013, respectively.As can be seen, the downstream island wakes are well reproduced.The flow field is affected by the typhoon, especially in the shallow water area and the of typhoon, Kuroshio current decreases.At S3, we can see the moment when the center of the typhoon is closest to the Green Island.The wind starts to change its direction and to blow northward.Comparing flow vectors along cross section 1 0 Q Q at S2 with that at S1 and S3 of Figure 6, it is noticed that downstream island wakes are significantly affected.Kuroshio currents increase after S2 due to the northward winding and resumes its normal flow condition after typhoon passes far away.Since TOROS data are the surface flow quantities, they are larger than the depth-averaged model predictions, in general.Moreover, the left one-third of the cross section (0-20 km) is right in the lee of Green Island.The meandering downstream island wakes are observed in the modelling, but not present in the TOROS data, because of its coarse spatial resolution (~10 km).However, computed results show Kuroshio currents increase when flow is in the same direction as the counterclockwise rotation of typhoon, and vice versa.This finding is consistent with the TOROS measured datasets [45].
(a)  Since TOROS data are the surface flow quantities, they are larger than the depth-averaged model predictions, in general.Moreover, the left one-third of the cross section (0-20 km) is right in the lee of Green Island.The meandering downstream island wakes are observed in the modelling, but not present in the TOROS data, because of its coarse spatial resolution (~10 km).However, computed results show Kuroshio currents increase when flow is in the same direction as the counterclockwise rotation of typhoon, and vice versa.This finding is consistent with the TOROS measured datasets [45].Since TOROS data are the surface flow quantities, they are larger than the depth-averaged model predictions, in general.Moreover, the left one-third of the cross section (0-20 km) is right in the lee of Green Island.The meandering downstream island wakes are observed in the modelling, but not present in the TOROS data, because of its coarse spatial resolution (~10 km).However, computed results show Kuroshio currents increase when flow is in the same direction as the counterclockwise rotation of typhoon, and vice versa.This finding is consistent with the TOROS measured datasets [45]. (a)

Effect of Model Typhoon on Green Island Wakes
Figure 8 depicts flow streamlines around the Green Island of no wind forcing in a period of vortex shedding.There are small recirculations in the lee of Green Island.Its size is about the size of the island (~8 km).A well-organized coherent and alternating asymmetric meandering eddies in the lee of the Green Island is reproduced, which has been observed in field measurement [17] and satellite images [19,43].In this study, simulation using the boundary conditions from HYCOM produces a more realistic Kuroshio flow and downstream Green Island wakes.

Effect of Model Typhoon on Green Island Wakes
Figure 8 depicts flow streamlines around the Green Island of no wind forcing in a period of vortex shedding.There are small recirculations in the lee of Green Island.Its size is about the size of the island (~8 km).A well-organized coherent and alternating asymmetric meandering eddies in the lee of the Green Island is reproduced, which has been observed in field measurement [17] and satellite images [19,43].In this study, simulation using the boundary conditions from HYCOM produces a more realistic Kuroshio flow and downstream Green Island wakes.According to the statistical analysis of typhoons from 1911 to 2010 by the CWB of Taiwan, 83.6% of typhoons pass the eastern Taiwan coastal waters and 13.5% pass the western Taiwan coastal waters.Therefore, two typical tracks of typhoons, namely the SN typhoon (17% of total) moving from southwest to northeast and the EW typhoon (67% of total) moving from east to west, shown in  According to the statistical analysis of typhoons from 1911 to 2010 by the CWB of Taiwan, 83.6% of typhoons pass the eastern Taiwan coastal waters and 13.5% pass the western Taiwan coastal waters.Therefore, two typical of typhoons, namely the SN typhoon (17% of total) moving from southwest to northeast and the EW typhoon (67% of total) moving from east to west, shown in Figure 9, are chosen to investigate the effect of typhoon on the Kuroshio and Green Island wakes.Two typical values of the typhoon moving speed, 2.5 m/s (9 km/h) and 5.0 m/s (18 km/h), are considered.Value of the radius of maximum wind speed R mw = 50 km and pressure of typhoon center P c = 950 hPa are used in the Holland's wind field model.The resulting maximum of W 10 is about 47 m/s.In order to better quantify the net influence of typhoon on Green Island wakes, we subtract the flow field of SN typhoon case and EW typhoon case with no wind forcing case.The net influence of typhoon on the flow field is defined by Figures 10 and 11 show the "net" u-and v-contours of SN typhoon at T0-T2 three time instances as indicated in Figure 16.It is noted that the effect of typhoon on the y-component velocity are more pronounced than on the x-component velocity, because the direction of moving SN typhoon and the mainstream of Kuroshio is the same, in general.Due to the counterclockwise rotation of typhoon, flow accelerates (decelerates) when it flows in the favorable (adverse) direction of wind.For example, Figures 14 and 15 show that y-component velocity increases in the right of SN typhoon's track and decreases in the left of SN typhoon's track, the so called rightward bias phenomenon.The rightward bias phenomenon of u contours is not as evident as v contours, since the northeastward Kuroshio and x-component wind forcing are not in the same or opposite direction.
Comparing Figure 10 with Figure 11, we notice that the impact of the slow-moving typhoon (WT = 2.5 m/s) is more significant than the fast-moving typhoon (WT = 5.0 m/s), especially in the ycomponent velocity, when the typhoon moves in favor of the Kuroshio mainstream direction.When the typhoon moves slowly, it takes a longer time to travel the same distance.Therefore, more momentum and energy are input to water, and water takes more time to respond to the wind forcing.In order to better quantify the net influence of typhoon on Green Island wakes, we subtract the flow field of SN typhoon case and EW typhoon case with no wind forcing case.The net influence of typhoon on the flow field is defined by Figures 10 and 11 show the "net" uand v-contours of SN typhoon at T 0 -T 2 three time instances as indicated in Figure 16.It is noted that the effect of typhoon on the y-component velocity are more pronounced than on the x-component velocity, because the direction of moving SN typhoon and the mainstream of Kuroshio is the same, in general.Due to the counterclockwise rotation of typhoon, flow accelerates (decelerates) when it flows in the favorable (adverse) direction of wind.For example, Figures 14 and 15 show that y-component velocity increases in the right of SN typhoon's track and decreases in the left of SN typhoon's track, the so called rightward bias phenomenon.The rightward bias phenomenon of u contours is not as evident as v contours, since the northeastward Kuroshio and x-component wind forcing are not in the same or opposite direction.
Comparing Figure 10 with Figure 11, we notice that the impact of the slow-moving typhoon (W T = 2.5 m/s) is more significant than the fast-moving typhoon (W T = 5.0 m/s), especially in the y-component velocity, when the typhoon moves in favor of the Kuroshio mainstream direction.When the typhoon moves slowly, it takes a longer time to travel the same distance.Therefore, more momentum and energy are input to water, and water takes more time to respond to the wind forcing.In order to better quantify the net influence of typhoon on Green Island wakes, we subtract the flow field of SN typhoon case and EW typhoon case with no wind forcing case.The net influence of typhoon on the flow field is defined by Figures 10 and 11 show the "net" u-and v-contours of SN typhoon at T0-T2 three time instances as indicated in Figure 16.It is noted that the effect of typhoon on the y-component velocity are more pronounced than on the x-component velocity, because the direction of moving SN typhoon and the mainstream of Kuroshio is the same, in general.Due to the counterclockwise rotation of typhoon, flow accelerates (decelerates) when it flows in the favorable (adverse) direction of wind.For example, Figures 14 and 15 show that y-component velocity increases in the right of SN typhoon's track and decreases in the left of SN typhoon's track, the so called rightward bias phenomenon.The rightward bias phenomenon of u contours is not as evident as v contours, since the northeastward Kuroshio and x-component wind forcing are not in the same or opposite direction.
Comparing Figure 10 with Figure 11, we notice that the impact of the slow-moving typhoon (WT = 2.5 m/s) is more significant than the fast-moving typhoon (WT = 5.0 m/s), especially in the ycomponent velocity, when the typhoon moves in favor of the Kuroshio mainstream direction.When the typhoon moves slowly, it takes a longer time to travel the same distance.Therefore, more momentum and energy are input to water, and water takes more time to respond to the wind forcing.the moment when the typhoon passes the Green Island.The third time instance (T 2 ), represented by the blue vertical dashed line, indicates the time instance when the typhoon is located 50 km northeast of Green Island.This shows that the typhoon has little influence on the flow field near P when it approaches P and is far away from P. However, the influence of the typhoon increases as the typhoon approaches P, and remains significant in a short period after the typhoon has passed far away from P.
Atmosphere 2018, 9, x FOR PEER REVIEW 19 of 22 Green Island.The second time instance (T1), represented by the green vertical dashed line, indicates the moment when the typhoon passes the Green Island.The third time instance (T2), represented by the blue vertical dashed line, indicates the time instance when the typhoon is located 50 km northeast of Green Island.This shows that the typhoon has little influence on the flow field near P when it approaches P and is far away from P. However, the influence of the typhoon increases as the typhoon approaches P, and remains significant in a short period after the typhoon has passed far away from P.

Conclusions
Safety, maintenance, and operations of Kuroshio power plant may be subjected to impacts from earthquakes, typhoons, climate change, and other natural factors.Typhoons have a significant impact on the physical and chemical characteristics of the water of the Kuroshio.A shallow water model based on the shallow water equations and the space-time least-squares finite-element method has been developed.The model has been applied to study the hydrodynamics of Kuroshio and Green Island wakes with a small computational domain (72 km × 156 km) and coarse mesh resolution [19,25] (500-3500 m).The spatial and temporal scales of meandering downstream Green Island wakes have been quantified and compared with the in-situ measurements and satellite images.Special efforts of the present study are made to refine the model, including to integrate Holland's wind field model [35], to use a larger computational domain (124 km× 220 km) with a finer mesh resolution (250-2250 m), and to use flow field from HYCOM as the initial and boundary conditions to drive the flows.

Conclusions
Safety, maintenance, and operations of Kuroshio power plant may be subjected to impacts from earthquakes, typhoons, climate change, and other natural factors.Typhoons have a significant impact on the physical and chemical characteristics of the water of the Kuroshio.A shallow water model based on the shallow water equations and the space-time least-squares finite-element method has been developed.The model has been applied to study the hydrodynamics of Kuroshio and Green Island wakes with a small computational domain (72 km × 156 km) and coarse mesh resolution [19,25] (500-3500 m).The spatial and temporal scales of meandering downstream Green Island wakes have been quantified and compared with the in-situ measurements and satellite images.Special efforts of the present study are made to refine the model, including to integrate Holland's wind field model [35], to use a larger computational domain (124 km× 220 km) with a finer mesh resolution (250-2250 m), and to use flow field from HYCOM as the initial and boundary conditions to drive the flows.
A high resolution (250-2250 m) shallow-water model is used to investigate effect of typhoon on the hydrodynamics of the Kuroshio and Green Island wakes.In the first case, we simulate the typhoon Soulik event to validate the applicability of the SWEs model.Computed results indicate that salient characteristics of Kuroshio and meandering downstream island wakes seem less affected by the typhoon because typhoon Soulik is 250 km away Green Island and wind speed is smaller than 10 m/s near the region of Green Island.However, Kuroshio currents increase when flow is in the same direction as the counterclockwise rotation of typhoon, and vice versa.This finding is consistent with the TOROS observed dataset [45].
In the second case, the SWEs model, forced by the Kuroshio and Holland's wind field model, successfully reproduces the downstream recirculation and meandering vortex streets.Kuroshio and the downstream Green Island wakes are found to be significantly affected by the moving typhoon that passes the Green Island directly, especially for the shallow waters and the lee of the islands.Computed results clearly reveal the rightward bias phenomenon-due to the counterclockwise rotation of the typhoon, Kuroshio currents increase in the right of the moving typhoon's track and decrease in the left of the moving typhoon's track.Computed results also show that the slower the moving typhoon, the more significant is the impact of the typhoon on the Kuroshio and downstream Green Island wakes.

Figure 1 .
Figure 1.(a) Schematic illustration of the study area and the bathymetry derived from the ETOPO1 and (b) computational meshes with fine meshes used in the coastal waters around Lanyu Island and Green Island.

Figure 1 .
Figure 1.(a) Schematic illustration of the study area and the bathymetry derived from the ETOPO1 and (b) computational meshes with fine meshes used in the coastal waters around Lanyu Island and Green Island.
depicts the computational meshes.It contains 50,842 nodes and 101,125 triangles.Spatial resolution varies from 250 to 2250 m.Fine meshes are employed near Lanyu Island and Green Island since flow field is expected to change dramatically in these areas.The increment of time step ∆t = 45 s is used in the computations.The eight-day simulation of typhoon Soulik takes about 40 h on a personal computer equipped with an Intel(R) Core TM i7-4700HQ CPU@2.40GHz.

Figure 2 .
Figure 2. Study domain and boundaries, as well as streamlines from the hybrid coordinate ocean model) (HYCOM).

Figure 2 .
Figure 2. Study domain and boundaries, as well as streamlines from the hybrid coordinate ocean model) (HYCOM).

Figure 3
Figure3shows a composite ERS-1 SAR (European Remote Sensing satellite1, Synthetic Aperture Radar) image of downstream island vortex trains of Lanyu Island and Green Island taken with one week interval sequentially between 17:01 UT 25 September (lower) and 19:02 UT 2 October (top), 1996.Two downstream island vortex trains are clearly seen.One a small distance from the coast, called the Lanyu Island (also called the Orchid Island) vortex train, occurs downstream of Lanyu Island (a little beyond the southern boundary of the image), and consists of three vortexes, two cyclonic appearing as bright shading and one anticyclonic appearing as dark shading vortexes.Another near the eastern coast of Taiwan, called the Green Island vortex train, occurs downstream of Green Island, and consists of three pairs of cyclonic vortexes (bright center) and anticyclonic vortexes (dark center).Spatial and temporal scales of downstream Green Island wakes due to passing of the Kuroshio has been numerically studied[19].The main difference between Liang et al.[19] and the present study lies in a constant inflow used in Liang et al. is replaced by a spatial varying inflow from HYCOM in the present study.The SWEs model is applied to simulation the typhoon Soulik event.Typhoon Soulik hit the Taiwan area during 7 to 14 July 2013.Figure4depicts the track of typhoon Soulik.Typhoon Soulik, developed from a tropical depression on 7 July 2013 through a tropical storm, moderate typhoon, severe typhoon, and then decreased its strength and made landfall on Taiwan on 7 December 2013.The typhoon Soulik data with temporal intervals of every 6 h are from the CWB of Taiwan, which are computed by the NFS (non-hydrostatic forecast system).NFS uses the structured orthogonal meshes; however, the SWEs model uses the unstructured triangular meshes.Therefore, wind data needs to be interpolated into the computational nodes of the SWEs model.The IDW (inverse distance weighted interpolation[44]) is employed to interpolate the data from NFS into nodal points of the SWEs model.In this study, the 5-km resolution typhoon Soulik data is used.The maximum wind speed of typhoon Soulik is 51 m/s, categorized as a severe typhoon.The closest distance of typhoon to Green Island is about 250 km.The maximum wind speed of the model area is 17 m/s and the average value is smaller than 10 m/s.

Figure 3 .
Figure 3. (a) A composite ERS-1 SAR image of island-induced ocean vortex trains of Lanyu Island and Green Island taken at 17:01 UT 25 September (lower) and 19:02 UT 2 October (top), 1996.(b) Zoomed Green Island and head of vortex train (modified from Figure 2 of Zheng and Zheng [43]).Blue circle represents for cyclonic (bright shading) vortexes and red circle represents for anticyclonic (dark shading) vortexes.
Figure 5b depicts the local view of the sub-domain as well as the location of the cross section 1 0 Q Q and monitor point P at downstream Green Island, where time series of flow quantities are presented and analyzed.Length of 1 0 Q Q is 20 km.It is located about 22 km downstream the Green Island.

Figure 3 .
Figure 3. (a) A composite ERS-1 SAR image of island-induced ocean vortex trains of Lanyu Island and Green Island taken at 17:01 UT 25 September (lower) and 19:02 UT 2 October (top), 1996.(b) Zoomed Green Island and head of vortex train (modified from Figure 2 of Zheng and Zheng [43]).Blue circle represents for cyclonic (bright shading) vortexes and red circle represents for anticyclonic (dark shading) vortexes.

Figure 5 .
Figure 5. (a) Global view of the study area, where the enclosed region near the Green Island indicates the sub-domain that will be presented in the succeeding plots of the downstream recirculation and island wakes, and (b) local view of the sub-domain as well as location of the cross section

Figure 6
Figure 6 depicts the typhoon wind field and flow streamlines as well as vectors of wind (red) and flow (blue) along cross section 1 0 Q Q at S0-S4 five time instances with a 12 h of increment, shown in Figure5, where S0 represents 00:00 of 12 July, S1 12:00 of 12 July, S2 00:00 of 13 July, S3 12:00 of 13 July, and S4 00:00 of 14 July, 2013, respectively.As can be seen, the downstream island wakes are well reproduced.The flow field is affected by the typhoon, especially in the shallow water area and the lee of the islands.Interactions of Kuroshio and downstream island wakes with typhoon are obvious.
denotes the tropical depression with maximum wind speed Vmax < 17.2 m/s; Atmosphere 2018, 9, x FOR PEER REVIEW 8 of 22

Figure 5 .
Figure 5. (a) Global view of the study area, where the enclosed region near the Green Island indicates the sub-domain that will be presented in the succeeding plots of the downstream recirculation and island wakes, and (b) local view of the sub-domain as well as location of the cross section

Figure 6
Figure 6 depicts the typhoon wind field and flow streamlines as well as vectors of wind (red) and flow (blue) along cross section 1 0 Q Q at S0-S4 five time instances with a 12 h of increment, shown in Figure5, where S0 represents 00:00 of 12 July, S1 12:00 of 12 July, S2 00:00 of 13 July, S3 12:00 of 13 July, and S4 00:00 of 14 July, 2013, respectively.As can be seen, the downstream island wakes are well reproduced.The flow field is affected by the typhoon, especially in the shallow water area and the lee of the islands.Interactions of Kuroshio and downstream island wakes with typhoon are obvious.

Figure 6
Figure 6 depicts the typhoon wind f and flow (blue) along cross section 1 0 Q Q in Figure 5, where S0 represents 00:00 of 1 July, and S4 00:00 of 14 July, 2013, respecti reproduced.The flow field is affected by lee of the islands.Interactions of Kuroshio

Figure 5 .
Figure 5. (a) Global view of the study area, where the enclosed region near the Green Island indicates the sub-domain that will be presented in the succeeding plots of the downstream recirculation and island wakes, and (b) local view of the sub-domain as well as location of the cross section

Figure 6
Figure 6 depicts the typhoon wind field and flow streamlines as well as vectors of wind (red) and flow (blue) along cross section 1 0 Q Q at S0-S4 five time instances with a 12 h of increment, shown in Figure5, where S0 represents 00:00 of 12 July, S1 12:00 of 12 July, S2 00:00 of 13 July, S3 12:00 of 13 July, and S4 00:00 of 14 July, 2013, respectively.As can be seen, the downstream island wakes are well reproduced.The flow field is affected by the typhoon, especially in the shallow water area and the lee of the islands.Interactions of Kuroshio and downstream island wakes with typhoon are obvious.

Figure 5 .
Figure 5. (a) Global view of the study area, where the enclosed region near the Green Island indicates the sub-domain that will be presented in the succeeding plots of the downstream recirculation and island wakes, and (b) local view of the sub-domain as well as location of the cross section

Figure 6
Figure 6 depicts the typhoon wind field and flow streamlines as well as vectors of wind (red) and flow (blue) along cross section 1 0 Q Q at S0-S4 five time instances with a 12 h of increment, shown in Figure5, where S0 represents 00:00 of 12 July, S1 12:00 of 12 July, S2 00:00 of 13 July, S3 12:00 of 13 July, and S4 00:00 of 14 July, 2013, respectively.As can be seen, the downstream island wakes are well reproduced.The flow field is affected by the typhoon, especially in the shallow water area and the

Figure 5 .
Figure 5. (a) Global view of the study area, where the enclosed region near the Green Island indicates the sub-domain that will be presented in the succeeding plots of the downstream recirculation and island wakes, and (b) local view of the sub-domain as well as location of the cross section Q 0 Q 1 and monitor point P at downstream Green Island.Length of Q 0 Q 1 is 20 km.It locates about 22 km downstream the Green Island.

Figure
Figure 7a illustrates the location of 22.84° N (y = 150 km) cross section, where black dots denote the position of TOROS data, and Figure7b,c depicts the comparison of the TOROS-observed surface currents with the computed flow vectors along the 22.84° N cross section at S0-S4: five time instances.Since TOROS data are the surface flow quantities, they are larger than the depth-averaged model predictions, in general.Moreover, the left one-third of the cross section (0-20 km) is right in the lee of Green Island.The meandering downstream island wakes are observed in the modelling, but not present in the TOROS data, because of its coarse spatial resolution (~10 km).However, computed results show Kuroshio currents increase when flow is in the same direction as the counterclockwise rotation of typhoon, and vice versa.This finding is consistent with the TOROS measured datasets[45].

Figure 6 .
Figure 6.(a) W 10 vectors and contours (unit: m/s), (b) flow streamlines, and (c) wind vectors (red) and flow vectors (blue) of cross section Q 0 Q 1 at S 0 -S 4 five time instances of the typhoon Soulik, see Figure 5.

Figure
Figure 7a illustrates the location of 22.84 • N (y = 150 km) cross section, where black dots denote the position of TOROS data, and Figure 7b,c depicts the comparison of the TOROS-observed surface currents with the computed flow vectors along the 22.84 • N cross section at S 0 -S 4 : five time instances.Since TOROS data are the surface flow quantities, they are larger than the depth-averaged model predictions, in general.Moreover, the left one-third of the cross section (0-20 km) is right in the lee of Green Island.The meandering downstream island wakes are observed in the modelling, but not present in the TOROS data, because of its coarse spatial resolution (~10 km).However, computed results show Kuroshio currents increase when flow is in the same direction as the counterclockwise rotation of typhoon, and vice versa.This finding is consistent with the TOROS measured datasets[45].

Figure
Figure 7a illustrates the location of 22.84° N (y = 150 km) cross section, where black dots denote the position of TOROS data, and Figure7b,c depicts the comparison of the TOROS-observed surface currents with the computed flow vectors along the 22.84° N cross section at S0-S4: five time instances.Since TOROS data are the surface flow quantities, they are larger than the depth-averaged model predictions, in general.Moreover, the left one-third of the cross section (0-20 km) is right in the lee of Green Island.The meandering downstream island wakes are observed in the modelling, but not present in the TOROS data, because of its coarse spatial resolution (~10 km).However, computed results show Kuroshio currents increase when flow is in the same direction as the counterclockwise rotation of typhoon, and vice versa.This finding is consistent with the TOROS measured datasets[45].

Figure 7 .
Figure 7. (a) Location of the 22.84 • N (y = 150 km) cross section, where black dots denote the position of TOROS data.(b) and (c) depict comparison of flow vectors of TOROS data and model predictions along 22.84 • N cross section at S 0 -S 4 : five time instances, see Figure 5.

Figure 8 .
Figure 8. Flow streamlines of the Green Island wakes in the no wind forcing case.

Figure 9 ,
are chosen to investigate the effect of typhoon on the Kuroshio and Green Island wakes.Two typical values of the typhoon moving speed, 2.5 m/s (9 km/h) and 5.0 m/s (18 km/h), are considered.Value of the radius of maximum wind speed Rmw = 50 km and pressure of typhoon center Pc = 950 hPa are used in the Holland's wind field model.The resulting maximum of W10 is about 47 m/s.

Figure 8 .
Figure 8. Flow streamlines of the Green Island wakes in the no wind forcing case.

Figure 9 .
Figure 9. Schematic diagram of the study domain and track of SN and EW typhoon.

Figure 9 .
Figure 9. Schematic diagram of the study domain and track of SN and EW typhoon.

Atmosphere 2018, 9 , 22 Figure 9 .
Figure 9. Schematic diagram of the study domain and track of SN and EW typhoon.

Figure 10 .
Figure 10.Contours of (a) W10, "net" (b) u-, and (c) v-contours of (unit: m/s) SN typhoon with moving speed WT = 2.5 m/s at T0-T2 three time instances indicated by the three vertical lines of Figure 16a.Figure 10.Contours of (a) W 10 , "net" (b) u-, and (c) v-contours of (unit: m/s) SN typhoon with moving speed W T = 2.5 m/s at T 0 -T 2 three time instances indicated by the three vertical lines of Figure 16a.

Figure 10 .Figure 10 .
Figure 10.Contours of (a) W10, "net" (b) u-, and (c) v-contours of (unit: m/s) SN typhoon with moving speed WT = 2.5 m/s at T0-T2 three time instances indicated by the three vertical lines of Figure 16a.Figure 10.Contours of (a) W 10 , "net" (b) u-, and (c) v-contours of (unit: m/s) SN typhoon with moving speed W T = 2.5 m/s at T 0 -T 2 three time instances indicated by the three vertical lines of Figure 16a.

Figure 11 .
Figure 11.Contours of (a) W10, "net" (b) u-, and (c) v-contours of (unit: m/s) SN typhoon with moving speed WT = 5.0 m/s at T0-T2 three time instances indicated by the three vertical lines of Figure 16b.

Figure 12
Figure 12 plots the wind field as well as wind vectors (red) and flow (blue) vectors along

Figure 11 .
Figure 11.Contours of (a) W 10 , "net" (b) u-, and (c) v-contours of (unit: m/s) SN typhoon with moving speed W T = 5.0 m/s at T 0 -T 2 three time instances indicated by the three vertical lines of Figure 16b.

Figure 12 Figure 11 .
Figure12plots the wind field as well as wind vectors (red) and flow (blue) vectors along Q 0 Q 1 at T 0 -T 2 three time instances.The flow field and downstream island wakes are significantly affected by the northeastward typhoon.Kuroshio currents increase as the typhoon approaches and decrease as the typhoon passes away.

Figure 12
Figure 12 plots the wind field as well as wind vectors (red) and flow (blue) vectors along

Figures 13
Figures 13 and 14 plot the "net" u-and of EW typhoon at T0-T2 three time instances indicated in Figure 16.The rightward bias feature due to the counterclockwise rotation of typhoon is evident, especially for the slow moving typhoon and y-component velocities.Comparing Figures 10 and 11 with Figures 13 and 14, the impact of SN typhoon on the Kuroshio and downstream Green Island wakes is apparently more significant than EW typhoon.The effect of the typhoon strengthens when the typhoon moves in favor of the Kuroshio and weakens when the typhoon moves against the Kuroshio.

Figures 13 Figure 12 .
Figures 13 and 14 plot the "net" uand v-contours of EW typhoon at T 0 -T 2 three time instances indicated in Figure16.The rightward bias feature due to the counterclockwise rotation of typhoon is evident, especially for the slow moving typhoon and y-component velocities.Comparing Figures10 and 11with Figures13 and 14, the impact of SN typhoon on the Kuroshio and downstream Green Island wakes is apparently more significant than EW typhoon.The effect of the typhoon strengthens when the typhoon moves in favor of the Kuroshio and weakens when the typhoon moves against the Kuroshio.

Figures 13
Figures 13 and 14 plot the "net" u-and v-contours of EW typhoon at T0-T2 three time instances indicated in Figure 16.The rightward bias feature due to the counterclockwise rotation of typhoon is evident, especially for the slow moving typhoon and y-component velocities.Comparing Figures 10 and 11 with Figures 13 and 14, the impact of SN typhoon on the Kuroshio and downstream Green Island wakes is apparently more significant than EW typhoon.The effect of the typhoon strengthens when the typhoon moves in favor of the Kuroshio and weakens when the typhoon moves against the Kuroshio.

Figure 13 .
Figure 13.Contours of (a) W10, "net" (b) u-, and (c) v-contours (unit: m/s) of EW typhoon with moving speed WT = 2.5 m/s at three time instances T0-T2 indicated by three vertical lines of Figure 16a.

Figure 13 .Figure 13 .
Figure 13.Contours of (a) W 10 , "net" (b) u-, and (c) v-contours (unit: m/s) of EW typhoon with moving speed W T = 2.5 m/s at three time instances T 0 -T 2 indicated by the three vertical lines of Figure 16a.

Figure 14 .
Figure 14.Contours of (a) W 10 , "net"(b) u-, and (c) v-contours of (unit: m/s) EW typhoon with moving speed W T = 5.0 m/s at T 0 -T 2 three time instances indicated by the three vertical lines of Figure 16b.

Figure 15 22 Figure 14 .
Figure15plots the wind field as well as wind vectors (red) and flow vectors (blue) along Q 0 Q 1 at T 0 -T 2 three time instances of the EW typhoon case.The flow field and downstream island wakes are also significantly affected by the EW typhoon.Kuroshio currents decrease as the typhoon approaches and increase in a short period as the typhoon passes away.

Figure 15 Figure 15 .
Figure 15 plots the wind field as well as wind vectors (red) and flow vectors (blue) along

Figure 16
Figure 16 shows a comparison of time series of u and v of the monitor point P of no wind forcing case and SN typhoon case with (a) WT = 2.5 m/s and (b) WT = 5.0 m/s, respectively.Three vertical lines are drawn and represent the three time instances.The first time instance (T0), represented by the red vertical dashed line, indicates the time instance when the typhoon is located 50 km southwest of

Figure 15 .
Figure 15.(a) W 10 -contours and (b) wind vectors (red) and flow vectors (blue) of EW typhoon along cross section Q 0 Q 1 with moving speed W T = 2.5 m/s at T 0 -T 2 three time instances indicated by the three vertical lines of Figure 16a.

Figure 16
Figure 16 shows a comparison of time series of u and v of the monitor point P of no wind forcing case and SN typhoon case with (a) W T = 2.5 m/s and (b) W T = 5.0 m/s, respectively.Three vertical lines are drawn and represent the three time instances.The first time instance (T 0 ), represented by the red vertical dashed line, indicates the time instance when the typhoon is located 50 km southwest of Green Island.The second time instance (T 1 ), represented by the green vertical dashed line, indicates

Figure 16 .
Figure 16.Comparison of u and v of monitor point P of no wind, SN, and EW typhoon with the moving speed WT = (a) 2.5 m/s (upper row) and (b) 5.0 m/s (lower row), respectively.

Figure 16 .
Figure 16.Comparison of u and v of monitor point P of no wind, SN, and EW typhoon with the moving speed W T = (a) 2.5 m/s (upper row) and (b) 5.0 m/s (lower row), respectively.