On the Response of Zhejiang Coastal Waters to 12 Typhoons from 2011 to 2015

: Based on the hybrid coordinate ocean model (HYCOM) reanalysis data, the dynamical and thermal response of Zhejiang coastal waters to 12 typhoons from 2011 to 2015 was studied, and the relationship between the oceanic response and typhoon characteristics was analyzed. The HYCOM reanalysis data were validated by satellite-observed sea surface temperature data. Results show that all the 12 typhoons caused near-inertial waves, sea surface cooling and seaﬂoor warming, but signiﬁcant differences existed among them. Based on statistics, it was found that the domain-averaged near-inertial kinetic energy and seaﬂoor warming are signiﬁcantly correlated with the typhoon’s maximum wind speed, both of which show a positive relationship. The domain-averaged seaﬂoor warming also decreases with the typhoon’s translation speed. Moreover, there exists a positive correlation between the domain-averaged sea surface cooling and seaﬂoor warming. The different mechanisms that caused the seaﬂoor warming for the 12 typhoons were also investigated in this study.


Introduction
The Northwest Pacific Ocean is the region with the highest concentration of tropical cyclones (typhoons) around the world [1]. China is located to the west of the Northwest Pacific Ocean and has a long coastline. There are many typhoons landing on the southeastern coast of China per year, which cause severe damages. Given that there are many coastal facilities and marine ranches in the coastal waters of China, understanding the oceanic response to typhoons is of great importance.
Typhoons are a powerful weather system that has strong dynamical and thermal effects on the ocean as they pass over. Typhoons import a large amount of energy into the ocean in a short period through momentum exchange at the air-sea interface, causing strong near-inertial waves (NIWs) and turbulent mixing. This process can be divided into the forcing stage and relaxation stage. In the forcing stage, the wind stress of typhoons is strong, and the action time ranges from half a day to a day, which will induce NIWs [2][3][4][5]. In the relaxation stage, the NIWs caused by typhoons propagate to the deep ocean [6][7][8][9]. In terms of thermal response, a significant effect of typhoons on the ocean is the cooling of sea surface temperature (SST) and the warming of subsurface temperature (it is manifested as seafloor warming in coastal waters) [10]. There are three main processes through which typhoons cause SST cooling: turbulent mixing, advection (mainly upwelling), and air-sea heat exchange [11]. Both observations and numerical simulations suggest that nearly 75-90% of typhoon-induced SST cooling is attributed to the turbulent mixing at the bottom of the mixed layer, which brings the relatively colder water to the sea surface [12].
As part of the East China Sea, the Zhejiang coastal waters are frequently affected by typhoons. However, due to the lack of in situ observations, our knowledge about the response of Zhejiang coastal waters to typhoons is limited. Several previous studies [13][14][15] have reported that the hybrid coordinate ocean model (HYCOM) reanalysis data can reasonably reproduce the oceanic response to typhoons. Therefore, the HYCOM reanalysis results were adopted in this study to investigate the dynamical and thermal response of Zhejiang coastal waters to the 12 typhoons from 2011 to 2015, and the relationship between oceanic response and typhoon characteristics was revealed based on statistics, aiming to deepen our understanding of regional ocean dynamics. The paper is organized as follows. The HYCOM reanalysis data and corresponding data analysis methods are introduced in Section 2. The main results are presented in Section 3. Finally, a summary completes the paper in Section 4.
can reasonably reproduce the oceanic response to typhoons. Therefore, the HYCOM reanalysis results were adopted in this study to investigate the dynamical and thermal response of Zhejiang coastal waters to the 12 typhoons from 2011 to 2015, and the relationship between oceanic response and typhoon characteristics was revealed based on statistics, aiming to deepen our understanding of regional ocean dynamics. The paper is organized as follows. The HYCOM reanalysis data and corresponding data analysis methods are introduced in Section 2. The main results are presented in Section 3. Finally, a summary completes the paper in Section 4.

2.1.Typhoons
We selected 12 typhoons that passed over Zhejiang coastal waters from 2011 to 2015 for research. They are Meari and Muifa in 2011; Khunan, Haikui, Tembin and Bolaven in 2012; Leepi and Kongrey in 2013; Neoguri, Nakri and Fungwong in 2014; Chanhom in 2015. The time that they passed over Zhejiang coastal waters was mainly from June to September. Among them, Haikui, Fungwong and Chanhom landed on the Zhejiang coast. During the period in Zhejiang coastal waters (27-33° N, 120-130° E), Chanhom and Meari had the largest and smallest averaged maximum wind speeds, which were 43.33 and 12.82 m/s, respectively. The paths of the 12 typhoons are shown in Figure 1. Detailed information about the 12 typhoons is listed in Table 1. These data are available from the Tropical Cyclone Data Center of the China Meteorological Administration [16,17] (http://tcdata.typhoon.org.cn, accessed on 1 September 2021).

Satellite-Observed SST Data
The satellite-observed SST data during the passage of the 12 typhoons in the Zhejiang coastal waters were downloaded from the Remote Sensing Systems (http://www.remss.com, accessed on 1 September 2021) and compared with the HYCOM reanalysis data. These data have a spatial resolution of 1/4 • and a temporal interval of 1 day. Appendix A shows the comparison between the satellite-observed and HYCOM SST, which validates the reasonability of the HYCOM reanalysis data.

Methodology
The power spectrum analysis of horizontal velocity during each typhoon's transit period showed that there is a significant peak near the local inertial frequency. Although the local inertial frequency is close to those of diurnal tides in Zhejiang coastal waters, the peak near the local inertial frequency corresponds to typhoon-induced NIWs, because the HYCOM reanalysis results (version: GLBu0.08/expt_53.X) do not contain tidal forcing. According to the result of power spectrum analysis, the fourth-order Butterworth filter with cutoff frequencies of (0.77,1.03) cpd (0.79-1.05 times the local inertial frequency) was adopted to extract NIWs from the HYCOM reanalysis data.
To estimate the strength of typhoon-induced NIWs, the near-inertial kinetic energy (NIKE) was calculated as: where ρ 0 = 1024 kg/m 3 is the reference density, and u and v are the zonal and meridional near-inertial velocities, respectively. The depth-integrated NIKE was also calculated in this study.
Typhoons usually cause strong turbulent mixing in the ocean. To evaluate the intensity of turbulent mixing, the gradient Richardson number (Ri) was calculated: where N is the Brunt-Väisälä frequency and is calculated as: where g = 9.8 m/s 2 is the acceleration due to gravity, and ρ and z are the density and depth of seawater. Generally, Ri < 1/4 is regarded as the criterion for shear instability [18,19]. In this study, we first interpolated the velocities and densities onto uniform 5-m levels and then calculated the Ri at each level. The proportion of shear instability (Ri < 1/4) in the upper 50 m was used to assess the intensity of turbulent mixing. Since the HYCOM does not provide vertical velocity, it is diagnosed by calculating the divergence of horizontal velocities according to the continuity equation of the incompressible fluid [20]: ∂w ∂z = − ∂u ∂x + ∂v ∂y (4) Assuming that w at the seafloor is zero, we can estimate the direction of vertical velocity in the water column according to the divergence of horizontal velocities. If the divergence of horizontal velocities is always positive (negative) in the water column, w is negative (positive), which indicates the occurrence of downwelling (upwelling).

Dynamical and Thermal Response of Zhejiang Coastal Waters to Typhoon Muifa in 2011
The dynamical and thermal response of Zhejiang coastal waters to typhoon Muifa in 2011 was shown as an example.   Figure 3 shows the depth-integrated NIKE at 00:00 from 5 to 19 August with an interval of two days. On 5 August, Muifa-induced NIKE had appeared in this domain before Muifa entered. Muifa left the domain on 7 August, but the NIKE continuously increased until 11 August and reached the maximum of 6 kJ/m 2 . Then, the NIKE gradually dissipated. On 19 August, the NIKE was almost on the same level as that on 5 August. As shown in Figure 2b,d, Muifa induced NIWs at (29.2 • N, 128.7 • E) on 6 August, which gradually increased and reached the maximum (~0.2 m/s) around 11 August. The phase of Muifa-induced NIWs had an obvious pattern of upward propagation, which means that the vertical group speed was downward. In other words, the NIKE induced by Muifa propagated downward to the deep ocean. As shown in Figure 2, Muifa-induced NIWs could reach~800 m depth several days after its passage. Figure 3 shows the depth-integrated NIKE at 00:00 from 5 to 19 August with an interval of two days. On 5 August, Muifa-induced NIKE had appeared in this domain before Muifa entered. Muifa left the domain on 7 August, but the NIKE continuously increased until 11 August and reached the maximum of 6 kJ/m 2 . Then, the NIKE gradually dissipated. On 19 August, the NIKE was almost on the same level as that on 5 August. Moreover, it was found that the NIKE was significant where the water depth is greater than 300 m.  Figure 3 shows the depth-integrated NIKE at 00:00 from 5 to 19 August with an interval of two days. On 5 August, Muifa-induced NIKE had appeared in this domain before Muifa entered. Muifa left the domain on 7 August, but the NIKE continuously increased until 11 August and reached the maximum of 6 kJ/m 2 . Then, the NIKE gradually dissipated. On 19 August, the NIKE was almost on the same level as that on 5 August. Moreover, it was found that the NIKE was significant where the water depth is greater than 300 m. Zhejiang coastal waters not only had a dynamical response to Muifa but also exhibited a thermal response. Figure 4 shows the SST at 00:00 from 3 to 14 August with an interval of one day. Before Muifa entered the area on 3 August, the SST of the whole area Zhejiang coastal waters not only had a dynamical response to Muifa but also exhibited a thermal response. Figure 4 shows the SST at 00:00 from 3 to 14 August with an interval of one day. Before Muifa entered the area on 3 August, the SST of the whole area was relatively high (26-32 • C). On 5 and 6 August, slight SST cooling appeared. On 7 August when Muifa left the area, significant SST cooling appeared to the north of 31 • N with the minimum SST lower than 22 • C. With time going on, the SST cooling was enhanced. Until 14 August, significant SST cooling was still easily detected in the domain. To the north of 31 • N, the maximum SST cooling was approximately −9 • C, whereas to the south of 31 • N, the maximum SST cooling was approximately −4 • C. Figure 5 shows the seafloor temperature at 00:00 from 3 to 14 August, with an interval of one day. Before Muifa entered the area on 3 August, the seafloor temperature varied within a wide range from 0 to 30 • C. Due to the deep bathymetry, the seafloor temperature in the sea of 27-31 • N and 127-130 • E was about 0-10 • C. As the water depth becomes shallower, the seafloor became warmer. When Muifa entered the area on 6 August, the seafloor temperature of Zhejiang coastal waters became slightly warmer. Muifa left the area on 7 August, and the seafloor temperature in the area of 27-33 • N and 122-124 • E showed a significant increase. On 7 August, the seafloor warming reached a maximum of about 10 • C. Similar phenomenon has been reported in the area to the north of the Yangtze River [21]. Figure 6 shows the NIKE of the 12 typhoons in the Zhejiang coastal waters. Note that because typhoons Tembin and Bolaven were too close in time, we did not separate them. Strong NIKE was induced during each typhoon's transit. For each typhoon (or double typhoons Tembin and Bolaven), the strongest NIKE appeared to the right of the typhoon's track and the NIKE was significant where the water depth is greater than 300 m. Among the 12 typhoons, Neoguri caused the strongest NIKE (Figure 6h), followed by Tembin and Bolaven (Figure 6e) and Chanhom (Figure 6k), of which the maximum intensity was 45 kJ/m 2 , 25 kJ/m 2 and 25 kJ/m 2 , respectively. In contrast, the NIKE caused by the other typhoons was relatively small, mostly below 10 kJ/m 2 .   Figure 5 shows the seafloor temperature at 00:00 from 3 to 14 August, with an interval of one day. Before Muifa entered the area on 3 August, the seafloor temperature varied within a wide range from 0 to 30 °C. Due to the deep bathymetry, the seafloor temperature in the sea of 27-31° N and 127-130° E was about 0-10 °C. As the water depth becomes shallower, the seafloor became warmer. When Muifa entered the area on 6 August, the seafloor temperature of Zhejiang coastal waters became slightly warmer. Muifa left the area on 7 August, and the seafloor temperature in the area of 27-33° N and 122-124° E showed a significant increase. On 7 August, the seafloor warming reached a maximum of about 10 °C. Similar phenomenon has been reported in the area to the north of the Yangtze River [21].  Figure 6 shows the NIKE of the 12 typhoons in the Zhejiang coastal waters. Note that because typhoons Tembin and Bolaven were too close in time, we did not separate them. Strong NIKE was induced during each typhoon's transit. For each typhoon (or double typhoons Tembin and Bolaven), the strongest NIKE appeared to the right of the   Figure 7 shows the SST cooling of Zhejiang coastal waters caused by the 12 typhoons, which was calculated as the difference between the minimum SST within seven days after the passage of the typhoon and the SST before the typhoon. Among the 12 typhoons, Muifa caused the maximum SST cooling (Figure 7b), followed by Tembin and Bolaven (Figure 7e), of which the maximum SST cooling was −9 °C and −8 °C, respectively. In contrast, the SST cooling caused by the other typhoons was relatively small, mostly above −4 °C. For most typhoons (or double typhoons Tembin and Bolaven), the maximum SST cooling appeared to the right of the typhoon's track, but for Khunan and Kongrey, which have small values of the maximum wind speed (Table 1), the maximum SST cooling appeared to the left of the typhoon's track  Figure 7 shows the SST cooling of Zhejiang coastal waters caused by the 12 typhoons, which was calculated as the difference between the minimum SST within seven days after the passage of the typhoon and the SST before the typhoon. Among the 12 typhoons, Muifa caused the maximum SST cooling (Figure 7b), followed by Tembin and Bolaven (Figure 7e), of which the maximum SST cooling was −9 • C and −8 • C, respectively. In contrast, the SST cooling caused by the other typhoons was relatively small, mostly above −4 • C. For most typhoons (or double typhoons Tembin and Bolaven), the maximum SST cooling appeared to the right of the typhoon's track, but for Khunan and Kongrey, which have small values of the maximum wind speed (Table 1), the maximum SST cooling appeared to the left of the typhoon's track. Figure 8 shows the seafloor warming caused by the 12 typhoons, which was calculated as the difference between the maximum seafloor temperature within seven days after the passage of the typhoon and the seafloor temperature before the typhoon. The seafloor warming caused by the 12 typhoons was within 10 • C. Muifa caused the maximum seafloor warming (Figure 8b), which was followed by Haikui (Figure 8d). Although the characteristics of the 12 typhoons (Table 1) and the SST cooling caused by the 12 typhoons ( Figure 8) are different, the seafloor warming patterns induced by the 12 typhoons show a common feature: the seafloor warming was remarkable in the shallow waters near the Zhejiang coast. Figure 9 shows the correlation between the aforementioned oceanic response of the Zhejiang coastal waters and typhoon characteristics, including the maximum wind speed and translation speed ( Table 1). The SST cooling, seafloor warming and NIKE were averaged in the whole domain. From Figure 9a, it is easy to detect an obvious positive correlation between the SST cooling and the maximum wind speed, the same as [11]. In contrast, there is no significant change in the SST cooling with the increase of the typhoon's translation speed (Figure 9b), which suggests that there is no obvious correlation between them. As shown in Figure 9c, the seafloor warming shows an increasing trend with the maximum wind speed, which suggests a positive correlation between them. In contrast, the seafloor warming decreases with the increase of typhoon's translation speed (Figure 9d), indicating a negative correlation. As for the NIKE, it is clearly shown in Figure 9e that the NIKE becomes stronger with the increase of the maximum wind speed, whereas its relationship with the typhoon's translation speed is not obvious (Figure 9f). Moreover, the seafloor warming is also found to be positively correlated with the SST cooling ( Figure 10).  Figure 8 shows the seafloor warming caused by the 12 typhoons, which was calculated as the difference between the maximum seafloor temperature within seven days after the passage of the typhoon and the seafloor temperature before the typhoon. The seafloor warming caused by the 12 typhoons was within 10 °C. Muifa caused the maximum seafloor warming (Figure 8b), which was followed by Haikui (Figure 8d). Although the characteristics of the 12 typhoons (Table 1) and the SST cooling caused by the 12 typhoons ( Figure 8) are different, the seafloor warming patterns induced by the 12 typhoons show a common feature: the seafloor warming was remarkable in the shallow waters near the Zhejiang coast.   Figure 9 shows the correlation between the aforementioned oceanic response of the Zhejiang coastal waters and typhoon characteristics, including the maximum wind speed and translation speed ( Table 1). The SST cooling, seafloor warming and NIKE were averaged in the whole domain. From Figure 9a, it is easy to detect an obvious positive correlation between the SST cooling and the maximum wind speed, the same as [11]. In contrast, there is no significant change in the SST cooling with the increase of the typhoon's translation speed (Figure 9b), which suggests that there is no obvious correlation between them. As shown in Figure 9c, the seafloor warming shows an increasing trend   (Table 1) to sep rate their individual influence on the SST cooling, seafloor warming, and NIKE. The red line each subfigure denotes the linear fitting result. Note that because the NIKE caused by typho Neoguri is extremely high, it is not shown in (e,f) and is not considered in the fitting or correlati calculation. To study the significance of correlations between the oceanic response and typho  (Table 1) to separate their individual influence on the SST cooling, seafloor warming, and NIKE. The red line in each subfigure denotes the linear fitting result. Note that because the NIKE caused by typhoon Neoguri is extremely high, it is not shown in (e,f) and is not considered in the fitting or correlation calculation. To study the significance of correlations between the oceanic response and typhoon characteristics, the corresponding correlation coefficients and p-values were calculated and are listed in Table 2. It is generally considered that a small value of p (usually smaller than 0.05) means a significant correlation. As shown in Table 2, the correlation coeffi- To study the significance of correlations between the oceanic response and typhoon characteristics, the corresponding correlation coefficients and p-values were calculated and are listed in Table 2. It is generally considered that a small value of p (usually smaller than 0.05) means a significant correlation. As shown in Table 2, the correlation coefficient between the seafloor warming (NIKE) and the maximum wind speed is 0.70 (0.83) and the corresponding p-value is 0.02 (0.01), which means that the seafloor warming (NIKE) is significantly positively correlated with the maximum wind speed. The p-value between the seafloor warming and typhoon's translation speed (SST cooling) is 0.10 (0.08), which is only a little larger than 0.05, suggesting that the correlation between the seafloor warming and typhoon's translation speed (SST cooling) is generally acceptable, combined with the results shown in Figures 9 and 10. The correlation coefficient between the seafloor warming and typhoon's translation speed (SST cooling) is −0.54 (0.57), indicating a negative (positive) correlation between them. As for other correlations, the p-values are much greater than 0.05, indicating that these correlations are insignificant. These results are generally consistent with the results shown in Figures 9 and 10. Table 2. Correlation coefficients and p-values (in the bracket) between the oceanic response (SST cooling, seafloor warming and NIKE) and typhoon characteristics (the maximum wind speed and translation speed). It is generally recognized that both typhoon intensity and translation speed have impacts on the oceanic response, e.g., [22,23]. However, neither the NIKE nor the SST cooling in Zhejiang coastal waters shows a significant correlation with the typhoon's translation speed (Figure 9 and Table 2). We speculated that a possible cause is that the influence of the typhoon's translation speed is not isolated from other typhoon characteristics (e.g., typhoon intensity) and the influence of typhoon intensity on the oceanic response may be greater than that of the typhoon's translation speed in Zhejiang coastal waters.

Mechanisms of the Seafloor Warming in the Shallow Waters near the Zhejiang Coast
The oceanic turbulent mixing, advection and air-sea heat exchange control the SST cooling [11,12], of which the first two processes can also contribute to the seafloor warming. Attention was first paid to the influence of turbulent mixing on the seafloor warming. Generally, Ri < 1/4 is regarded as the criterion for shear instability and occurrence of turbulent mixing [18,19]. Therefore, the proportion of shear instability (Ri < 1/4) in the upper 50 m was used to assess the intensity of turbulent mixing in this study. Figure 11 shows the proportions of shear instability in the upper 50 m during the passage of typhoon Muifa as an example. It is clearly shown that when Muifa entered the study domain at 00:00 on 6 August (Figure 11a), strong turbulent mixing was found as the proportion of shear instability was nearly 100% at most of the area. Note that the proportion of shear instability in the typhoon eye was very small. This phenomenon lasted to 06:00 on 7 August when Muifa's center reached the northern boundary of the study domain (Figure 11f). Thereafter, the turbulent mixing decreased gradually. After 00:00 on 8 August, the proportion of shear instability in the upper 50 m at most areas was smaller than 40%. These results confirm the occurrence of enhanced turbulent mixing induced by typhoon Muifa. Figure 12 shows the proportions of shear instability in the upper 50 m for the 12 typhoons. It is clearly shown that for most typhoons (Muifa, Tembin, Bolaven, Neoguri, Nakri, Fungwong and Chanhom), strong turbulent mixing with the proportion of shear instability close to 100% appeared in the shallow waters near Zhejiang coast. For the other typhoons, the proportion of shear instability in the shallow waters near Zhejiang coast was much smaller. Specially, for typhoons Khunan and Leepi (Figure 12c,f), the proportion of shear instability in the shallow waters near Zhejiang coast was almost 0, which suggests that no enhanced turbulent mixing was caused by the two typhoons near Zhejiang coast. However, typhoons Khunan and Leepi indeed caused seafloor warming (Figure 8c,f), although it was almost the weakest among the 12 typhoons. This result suggests that in addition to the turbulent mixing, the advection (mainly downwelling) may play a role in the seafloor warming. Actually, according to the tracks of Khunan and Leepi, there existed north winds near the Zhejiang coast during the passage of the two typhoons. Given that the Zhejiang coast is nearly south-north, there should exist downwelling near the Zhejiang coast during the passage of Khunan and Leepi, according to the traditional Ekman theory. Figure 13 shows the divergence of horizontal velocities at (28 • N, 121.76 • E) during the passage of Khunan and Leepi. The water depth at the point is 20 m. It is clearly shown that after the passage of Khunan and Leepi, there exists a period (several days) in which the divergence of horizontal velocities in the whole water column was positive. According to Equation (4), it can be deduced that the vertical velocity in the water column in this period was negative, given that the vertical velocity was zero at the seafloor. In other words, downwelling appeared in this period, which accounts for the seafloor warming (Figure 13c,d). Actually, downwelling was also found for the other typhoons near the Zhejiang coast, which are not shown. Based on the above analysis, we can conclude that, for most typhoons in the Zhejiang coast waters, both typhoon-induced turbulent mixing and downwelling contribute to the seafloor warming; whereas for Khunan in 2012 and Leepi in 2013, the typhoon-induced turbulent mixing was not significant, and the seafloor warming was mostly related to downwelling.  Figure 12 shows the proportions of shear instability in the upper 50 m for the 12 typhoons. It is clearly shown that for most typhoons (Muifa, Tembin, Bolaven, Neoguri,Nakri, Fungwong and Chanhom), strong turbulent mixing with the proportion of shear instability close to 100% appeared in the shallow waters near Zhejiang coast. For the other typhoons, the proportion of shear instability in the shallow waters near Zhejiang coast was much smaller. Specially, for typhoons Khunan and Leepi (Figure 12c,f), the proportion of shear instability in the shallow waters near Zhejiang coast was almost 0, which suggests that no enhanced turbulent mixing was caused by the two typhoons near Zhejiang coast. However, typhoons Khunan and Leepi indeed caused seafloor warming (Figure 8c,f), although it was almost the weakest among the 12 typhoons. This result suggests that in addition to the turbulent mixing, the advection (mainly

Summary
Based on the HYCOM reanalysis data, the dynamical and thermal responses of the Zhejiang coastal waters to the 12 typhoons (Meari and Muifa in 2011; Khunan, Haikui, Tembin and Bolaven in 2012; Leepi and Kongrey in 2013; Neoguri, Nakri and Fungwong in 2014; Chanhom in 2015) were investigated in this study. Results indicate that all the 12 typhoons caused SST cooling, seafloor warming and NIKE. However, due to the different characteristics, the oceanic response of the Zhejiang coastal waters to the 12 typhoons was different. The SST cooling caused by the typhoons ranged from −9 to −1 • C, and the seafloor warming from 3 to 10 • C. For most typhoons, significant SST cooling appeared to the right of the typhoon's track, whereas the maximum SST cooling caused by Khunan and Kongrey was to the left of their tracks. In contrast, the seafloor warming caused by the 12 typhoons had a common feature, i.e., the seafloor warming was remarkable in the shallow waters near the Zhejiang coast. As for the typhoon-induced NIKE, large values were found to the right of the typhoon's track and in the region where the water depth is greater than 300 m.
To further understand the oceanic response of the Zhejiang coastal waters to typhoons, the correlation coefficients and corresponding p-values between the oceanic response (domain-averaged SST cooling, seafloor warming and NIKE) and typhoon characteristics (typhoon's maximum wind speed and translation speed) were calculated. Results indicate that the domain-averaged NIKE and seafloor warming are significantly positively correlated with the typhoon's maximum wind speed. The domain-averaged seafloor warming is also correlated with the typhoon's translation speed, but show a decreasing trend. Moreover, the domain-averaged seafloor warming is found to be correlated with the domain-averaged SST cooling, although their spatial patterns showed large differences.
In this study, the mechanisms to cause seafloor warming were also discussed. It is found that for most typhoons, the seafloor warming was caused by the combined effects of typhoon-induced turbulent mixing and downwelling. However, typhoons Khunan in 2012 and Leepi in 2013 did not cause strong turbulent mixing near the Zhejiang coast and downwelling was found to account for the seafloor warming caused by the two typhoons.
Because of the lack of in situ observations, we believe that the results of this study can deepen our understanding of the oceanic response of Zhejiang coastal waters to typhoons to some extent. Moreover, given that there are many marine ranches near the Zhejiang coast and some of the farmed fish and shellfish are sensitive to water temperature, the results of this study may be useful for the site selection and protection of the marine ranches.

Conflicts of Interest:
The authors declare no conflict of interest. Figure A1 shows the Muifa-induced SST cooling maps from HYCOM reanalysis data and satellite observations. On the whole, the HYCOM SST cooling is generally consistent with that from satellite observations; the maximum SST cooling appeared to the right of Muifa's track and the maximum SST cooling was approximately −9 • C. Figure A2 displays the comparison of domain-averaged SST cooling between HYCOM reanalysis results and satellite observations for the 10 typhoons except for Tembin and Bolaven. For most typhoons, the difference between the HYCOM and satellite-observed SST cooling was small and within ±0.5 • C. Only for the two typhoons (Muifa and Fungwong), the difference was larger than ±1 • C. This discrepancy is probably attributed to a combined effect of the lack of tidal forcing and inaccuracy of wind filed used in HYCOM. Note that the temporal intervals of HYCOM reanalysis data (3 h) and satellite observations (1 day) are different, which may be a cause for the discrepancy. Based on the aforementioned results, we think that the HYCOM reanalysis data can reasonably reproduce the thermal response of Zhejiang coastal waters to most typhoons and hence can be used in this study. most typhoons, the difference between the HYCOM and satellite-observed SST cooling was small and within ±0.5 °C. Only for the two typhoons (Muifa and Fungwong), the difference was larger than ±1 °C. This discrepancy is probably attributed to a combined effect of the lack of tidal forcing and inaccuracy of wind filed used in HYCOM. Note that the temporal intervals of HYCOM reanalysis data (3 h) and satellite observations (1 day) are different, which may be a cause for the discrepancy. Based on the aforementioned results, we think that the HYCOM reanalysis data can reasonably reproduce the thermal response of Zhejiang coastal waters to most typhoons and hence can be used in this study.      (Table 1) to separate their individual influence on SST cooling.