A Systematic Study on Berthing Capacity Assessment of Sanya Yazhou Fishing Port by Typhoon Prediction Model

: This paper sheds light on the effect of combination modes on the evaluation of berthing capacity for Sanya Yazhou Fishing Port (SYFP) under hypothetical typhoon conditions. By statistically analysing the maximum probability of moving speeds and directions of historical typhoons passing through the ﬁshing port, the representative typhoon path was determined with the nonparametric regression method. The designed typhoon wind ﬁelds of levels 12–17 were generated based on Holland’s parametric wind model. Then, the MIKE 21 BW model was used to obtain the high-precision wave distribution in the ﬁshing port. The boundary conditions (signiﬁcant wave height and peak period) of the MIKE 21 BW model were calculated by combining the MIKE 21 SW model with the designed typhoon wind ﬁelds. In SYFP, ships usually adopt the modes of multi-ship side-by-side and single anchor mooring during typhoons. In fair weather, approximately 158 vessels can be berthed if they are all large ones, while approximately 735 vessels can be moored if they are all small ones. However, with an increase in typhoon levels, the anchoring area for small vessels decreases. From the perspective of wave distribution in the ﬁshing port, the number of large vessels moored was hardly affected by typhoons. This can be attributed to the breakwater, which signiﬁcantly decreases the large wave height in the ﬁshing port. Finally, a study on the framework of a method for hazard assessment of berthing capacity in the coming typhoon-driven storm waves was set up.


Introduction
Typhoons are very severe natural disasters that can pose a serious threat to harbours [1,2]. A steady sea condition is a decisive factor for the management of berthed vessels in fishing ports [3][4][5]. Therefore, to ensure the safety of fishing ports and reduce the risk of typhoon disasters to the greatest extent, it is necessary to evaluate the berthing capacity of fishing ports [6][7][8]. There are several methods for studying the berthing capacity in fishing ports under typhoon fields [1,[9][10][11]. This study focuses on Sanya Yazhou Fishing Port (SYFP) as an example to establish a comprehensive evaluation system for the berthing capacity of fishing ports. Two primary components are used for this evaluation: the typhoon wind model and the wave model.
Typhoon wind fields, in particular, are important input parameters of the wave model, which can be simulated with different typhoon models [12]. Aside from this, the quality of the wind fields and historical climate are two important factors that affect the accuracy of the wave model [3,13]. In the past several decades, numerous wind models have been developed to hindcast and forecast typhoon surface wind fields [14]. Numerical simulations of hydrodynamics caused by typhoons and storms have been applied to the study of extreme hydrometeorology [6]. Holland developed a popular scheme, which has been used extensively for reconstructing typhoon wind fields in a wide range of applications [15]. Furthermore, an updated version of this approach was presented, which has less sensitivity to data errors and reproduces typhoon wind profiles with high temporal and spatial resolution [16].
In addition, a suitable wave model is a vital part of the evaluation of berthing capacity, which can in turn shape government policies, corporate budgets, and spending decisions [17]. Both the MIKE 21 SW model and tide-surge-wave coupled model have been applied to simulate the wave fields in many areas [18,19]. Hsiao demonstrated that the simulations of the maximum significant wave height have almost no discrepancies in the spatial distribution or the magnitude between the fully coupled and decoupled models [20]. The MIKE 21 SW model based on wave action balance equations is popular for simulating the wave distribution in large areas. In contrast, when considering a harbour shielded by a breakwater or artificial islands, the BW model of the enhanced Boussinesq-type equations is more suitable, as it can simulate the propagation of directional wave trains travelling from deep to shallow water [21,22]. Wang et al. used the MIKE 21 BW model to calculate the wave height distribution behind the breakwater in the Xiaoguoju fishing port project [23]. Further, Gou et al. studied the wave height distribution and berthing stability of different breakwater panning schemes to optimise the layout of the breakwater with the BW model [24].
By using the statistical results of the historical typhoon paths and the nearest record points from the historical typhoon paths to SYFP, this study used Holland's model to generate hypothetical designed typhoon wind fields of different levels. In light of this, the wave distribution of the model can be affected by the model itself. To study the evolution process of a given wave under hypothetical typhoons in SYFP, a high-precision BW model was established. This was used to calculate the generation, evolution, reflection, refraction, and diffraction of waves in SYFP combined with the MIKE 21 SW mode. This research is a systematic study on typhoon prevention assessment of SYFP and can be of great significance for formulating reasonable wind shelter measures.

Materials and Methods
To study the effect of typhoons on wave distribution in the fishing port, numerical models of different scales were built. In the large-scale model, the wind field generated by the Holland model drives the MIKE 21 SW model (a spectrum wave model) to calculate the wind waves. As boundary conditions for the MIKE 21 BW model (a Boussinesq wave model) are small scale, the calculated wave characteristics were used to obtain the wave distribution in the fishing port. A flowchart of the proposed method is shown in Figure 1.

Wind Field Model
According to Holland's parametric approach, the surface wind speed with a height of 10 m can be calculated using the following formula [16]: x bb

Wind Field Model
According to Holland's parametric approach, the surface wind speed with a height of 10 m can be calculated using the following formula [16]: where v ms is the maximum wind, R m is the radius of the maximum winds, and r is the distance between the cyclone centre and the calculating point (km). b s is a scaling parameter related to the proportion of the pressure gradient near the maximum wind radius. Alternatively, the maximum wind speed for a given pressure drop can be estimated as follows: where ∆p s is the pressure drop from a defined external pressure to the cyclone centre, and ∂p cs /∂t is the intensity change per hour. v f is the moving speed of the cyclone, ρ ms is the density of air, and ϕ is the absolute value of latitude in degrees. Typhoon parameters ∆p s and v f can be obtained from the China Meteorological Administration tropical cyclone database (CMA dataset) [25]. R m can be calculated using the following formula [26]: The asymmetric parametric wind model is calculated using the following formula [27]: where V(r, θ) denotes the asymmetric tangential wind field. ε represents the degree of asymmetry and α represents the azimuth that controls the location of the maximum wind speed. In this study, ε = V f /v gm , where v gm is the maximum wind speed of the gradient wind.

MIKE 21 SW Model
In horizontal Cartesian coordinates, the conservation equation for wave action can be written as: where N is the action density. c x , c y , c σ , and c θ are the propagation velocities of action in geographical (x, y) space, relative frequency, and wave propagation direction space dimensions, respectively. S is the source term of the energy density.
Here, S represents the superposition of source functions describing various physical phenomena: where S in represents the generation of energy by wind, S n1 is the wave energy transfer due to nonlinear wave-wave interaction, S ds is the dissipation of wave energy due to white capping, S bot is the dissipation due to bottom friction, and S sur f is the dissipation of wave energy due to depth-induced breaking.

MIKE 21 BW Model
The Boussinesq wave modules of MIKE 21 BW solve the enhanced Boussinesq equations [22]. These are expressed in one or two horizontal dimensions in terms of the free surface elevation (ξ), and the depth-integrated velocity components P and Q [22].
The Boussinesq equations read: Continuity equation: n ∂ξ ∂t Momentum equations in the horizontal direction: where ξ is the free surface elevation; p and q are flux densities in the directions x and y, respectively; h is the total water depth; D is the averaged water depth; and c is the Chezy resistance number, c = Mh 1/6 . In this study, M is the Manning coefficient, which is taken as 32 m 1/3 /s. Finally, E is the turbulence coefficient, and g is the gravitational acceleration.

Assessment of the Berthing Quantity of the Fishing Port
To evaluate the berthing capacity of SYFP under the influence of the designed typhoons, single anchor mooring and multiple vessels side-by-side with a single anchor mooring are used in SYFP for evaluation according to the "Master Design Code of Fishery Port (SC/T 9010-2000, 2000)" [28].

Single Anchor Mooring
This method was applied to the anchoring of fishing vessels in wind sheltering. The mooring area is a circular area with radius R, which can be calculated as follows: where R is the mooring radius, h 3 is the water depth of the anchorage at an extremely high water level, and L c is the full length of the representing designed vessel.

Multiple Vessels Abreast with a Single Anchor
This method applies to anchoring when the wind direction and water level in the harbour are stable. The anchorage area can be calculated as follows: where F 1 is the mooring area of each group with multiple vessels side-by-side and single anchor mooring; m 2 is the number of corresponding fishing vessels in each group, which can be taken as 2-6, where a small value is taken for the large vessels, while the large value for the small ones; lastly, B c represents the full width of the designed vessel type.
The number of berthed vessels in the harbour basin were estimated in combination with the type of representative vessels in accordance with Equations (12) and (13).

Wind Field Model
To obtain the representative designed typhoons, the typhoon paths with the maximum probability were determined based on the historical data of typhoons that passed through SYFP.
Firstly, the most likely starting point of the designed typhoon was determined. As shown in Figure Therefore, as seen in Figure 2, the red point represents the starting point, whose latitude and longitude coordinates are 114.95° E, 16.17° N.
The motion of the typhoon is described in Figure 3. Consistent with the record of the CMA, the position and central pressure of the typhoons were recorded every 6 h. Therefore, the next location of the typhoon centre record point could be determined based on the distance and direction of the movement for 6 h. The method used in this study fit the 6-h distance and angle according to the historical data of the typhoons within 2° of the position. The maximum probability density was then calculated as the distance and direction of the point moving downward. Typhoons were terminated when the number of data points within 2° of the point was less than 10. The nonparametric regression method was used to calculate the position with the highest density of the starting points. Consider X 1 , X 2 , . . . , X n as the sample from the population, and the kernel density estimate of the population density function as f (x). At any point, x is defined as [29]:f where K x−X i h is the kernel function and h is the window width. When the population obeys N(0, σ 2 ), the distribution, and the function K x−X i h is a Gaussian kernel function, the optimal window width is:ĥ Therefore, as seen in Figure 2, the red point represents the starting point, whose latitude and longitude coordinates are 114.95 • E, 16.17 The motion of the typhoon is described in Figure 3. Consistent with the record of the CMA, the position and central pressure of the typhoons were recorded every 6 h. Therefore, the next location of the typhoon centre record point could be determined based on the distance and direction of the movement for 6 h. The method used in this study fit the 6-h distance and angle according to the historical data of the typhoons within 2 • of the position. The maximum probability density was then calculated as the distance and direction of the point moving downward. Typhoons were terminated when the number of data points within 2 • of the point was less than 10.
At the nearest record point in this port, different intensities of central air pressure were assigned, with the maximum wind speeds corresponding to level 12-17 typhoons, respectively. Both the central air pressures at the starting and ending points were set to 1010 hPa. It should be noted that the typhoon had the greatest impact on the study area when it was close to the fishing port, while the rest of the time it had less impact. Therefore, the nearest point from the port was set as the strongest point of the typhoon. The increase and decrease in typhoon intensity were then simplified to a linear change. The time-history curves of the central pressure of typhoons of different levels are shown in Figure 4. At the nearest record point in this port, different intensities of central air pressure were assigned, with the maximum wind speeds corresponding to level 12-17 typhoons, respectively. Both the central air pressures at the starting and ending points were set to 1010 hPa. It should be noted that the typhoon had the greatest impact on the study area when it was close to the fishing port, while the rest of the time it had less impact. Therefore, the nearest point from the port was set as the strongest point of the typhoon. The increase and decrease in typhoon intensity were then simplified to a linear change. The time-history curves of the central pressure of typhoons of different levels are shown in Figure 4. Historical typhoon data in the study area suggests that typhoons below level 12 had little impact on the SYFP. Therefore, typhoons between levels 12 and 17 specified by the grade of tropical cyclones were taken as the designed typhoons for assessment, as shown in Figure 5a-f [30]. In the figure, the arrows indicate the wind direction, and the colours represent the wind speed.
An increase in typhoon levels leads to a gradual increase in the maximum wind speed. The maximum wind speeds of the level 12 and level 17 typhoons were 40 and approximately 60 m/s, respectively.  At the nearest record point in this port, different intensities of central air pressure were assigned, with the maximum wind speeds corresponding to level 12-17 typhoons, respectively. Both the central air pressures at the starting and ending points were set to 1010 hPa. It should be noted that the typhoon had the greatest impact on the study area when it was close to the fishing port, while the rest of the time it had less impact. Therefore, the nearest point from the port was set as the strongest point of the typhoon. The increase and decrease in typhoon intensity were then simplified to a linear change. The time-history curves of the central pressure of typhoons of different levels are shown in Figure 4. Historical typhoon data in the study area suggests that typhoons below level 12 had little impact on the SYFP. Therefore, typhoons between levels 12 and 17 specified by the grade of tropical cyclones were taken as the designed typhoons for assessment, as shown in Figure 5a-f [30]. In the figure, the arrows indicate the wind direction, and the colours represent the wind speed.
An increase in typhoon levels leads to a gradual increase in the maximum wind speed. The maximum wind speeds of the level 12 and level 17 typhoons were 40 and approximately 60 m/s, respectively. Historical typhoon data in the study area suggests that typhoons below level 12 had little impact on the SYFP. Therefore, typhoons between levels 12 and 17 specified by the grade of tropical cyclones were taken as the designed typhoons for assessment, as shown in Figure 5a-f [30]. In the figure, the arrows indicate the wind direction, and the colours represent the wind speed.
An increase in typhoon levels leads to a gradual increase in the maximum wind speed. The maximum wind speeds of the level 12 and level 17 typhoons were 40 and approximately 60 m/s, respectively.

Verification
The MIKE 21 SW model combined with the designed typhoon wind field model was used to analyse the wave distribution near the SYFP. As shown in Figure 6, the computational domain must be sufficiently large to simulate waves generated by typhoons travelling long distances to SYFP [31].
Two historical typhoons, Typhoon Doksuri June in 2012 and Typhoon Kai-Tak in August in 2012, were selected for verification of the numerical models. Figure 7 shows the validation of the simulated and measured significant wave heights in the study area during the onset of typhoons Doksuri and Kai-Tak. The location that recorded the measured data of the two typhoons are (21.

Verification.
The MIKE 21 SW model combined with the designed typhoon wind field model was used to analyse the wave distribution near the SYFP. As shown in Figure 6, the computational domain must be sufficiently large to simulate waves generated by typhoons travelling long distances to SYFP [31].  Figure 7 shows the validation of the simulated and measured significant wave heights in the study area dur- August in 2012, were selected for verification of the numerical models. Figure 7 shows the validation of the simulated and measured significant wave heights in the study area during the onset of typhoons Doksuri and Kai-Tak. The location that recorded the measured data of the two typhoons are (21.116 N, 112.617 E) and (22.284 N, 115.592 E). The simulated significant wave heights were in good agreement with the measured ones, indicating that the models were consistent with the actual situation. The SW model, therefore, provides a reasonable wave boundary for the next simulation calculation.

Setting of the Model
The range of the MIKE 21 BW model is shown in Figure 8. The water depth in the SYFP was obtained using a multi-beam sounding system in August 2017. The resolution of the measured water depth was approximately 10 m. For the purpose of this model, a

Setting of the Model
The range of the MIKE 21 BW model is shown in Figure 8. The water depth in the SYFP was obtained using a multi-beam sounding system in August 2017. The resolution of the measured water depth was approximately 10 m. For the purpose of this model, a structured rectangular grid was adopted. To obtain better simulation results, including the wave generation and evolution process in SYFP, the grid step was 10 m, as shown in Figure 9. There were 600 horizontal grids, 560 vertical grids, and 181,266 water grids.
J. Mar. Sci. Eng. 2021, 9, x FOR PEER REVIEW 10 of 16 structured rectangular grid was adopted. To obtain better simulation results, including the wave generation and evolution process in SYFP, the grid step was 10 m, as shown in Figure 9. There were 600 horizontal grids, 560 vertical grids, and 181,266 water grids.

Maximum Significant Wave Heights and Peak Wave Periods
The wave fields in large-scale areas were simulated using the SW model. The simulation results (including the maximum significant wave heights and peak wave periods) were obtained near the study area under the six designed typhoons. Then, the wave boundary conditions for the BW model of the six designed typhoons under the designed high water level were obtained, as shown in Table 1. Table 1. Wave boundary conditions for the BW model.

Level of Typhoon
Significant Wave Height (m) T (s) 12 2.21 10.08

Maximum Significant Wave Heights and Peak Wave Periods
The wave fields in large-scale areas were simulated using the SW model. The simulation results (including the maximum significant wave heights and peak wave periods) were obtained near the study area under the six designed typhoons. Then, the wave boundary conditions for the BW model of the six designed typhoons under the designed high water level were obtained, as shown in Table 1.

Wave Distribution in SYFP
The topography of the SYFP is complicated by a large proportion of shallow water areas. When the wave direction is in the WSW direction, the entrance of this port is almost parallel to the WSW direction. The positive wave incidence was adopted in this study to ensure that the wave propagates directly into the harbour basin [32]. In contrast, the port had a better shielding effect under other wave direction conditions. The calculation results of the wave distribution of the MIKE 21 BW model under the designed high water level are shown in Figure 10a-f.
With the propagation of directional wave trains travelling from deep to shallow water, wave refraction, diffraction, reflection, and breaking would occur, causing the wave energy to redistribute and create a new wave field. As shown in Figure 10, when the designed typhoon level was small, it had little influence on SYFP because of the good shielding of the harbour's gate. However, when the designed typhoon level increased, the values of the significant wave height were close to 1 m in local areas and 2 m at the entrance of the basin. When the waves were incident from the WSW direction, they propagated directly into the harbour basin because the entrance of SYFP was almost directly parallel to the direction of WSW. As such, if the wave travelled into the interior of the SYFP, the significant wave height in the inner most part of the harbour would be the largest. Similarly, the significant wave height in the centre of the study area was higher than that on the two sides.
The topography of the SYFP is complicated by a large proportion of shallow water areas. When the wave direction is in the WSW direction, the entrance of this port is almost parallel to the WSW direction. The positive wave incidence was adopted in this study to ensure that the wave propagates directly into the harbour basin [32]. In contrast, the port had a better shielding effect under other wave direction conditions. The calculation results of the wave distribution of the MIKE 21 BW model under the designed high water level are shown in Figure 10a The BW model can accurately simulate the wave distribution in a small range based on the numerical solution of Boussinesq-type equations in the time domain and the nonlinearity and dispersion of waves. This means that it is more suitable for the berthing quantity evaluation of SYFP as indicated in this study.

Berthing Quantity of SYFP
The average depth in the basin is 6 m, and the berthing area of SYFP is approximately 460,000 m 2 in fine weather. Under the influence of typhoons of different levels, the waves in the study area presented different distributions. According to the "Master Design Code of Fishery Port (SC/T 9010-2000)", small-and medium-sized fishing vessels can withstand the significant wave heights within 0.5 m, medium and large fishing vessels can bear the significant wave heights within 0.75 m, and large fishing vessels can withstand the significant wave heights within 1 m. However, when the significant wave height exceeds 1 m, the fishing vessels can no longer be berthed.
In this section, the MIKE 21 BW model is used to calculate the wave distributions and to estimate the number of fishing vessels that can be accommodated in SYFP under typhoons of different levels. Berthing areas of different typhoon levels are shown in Figure 11a-f. The wave distribution at different levels is divided into four areas according to wave height: 0-0.5 m (green), 0.5-0.75 m (light blue), 0.75-1 m (yellow), and >1 m (red).

Berthing Quantity of SYFP
The average depth in the basin is 6 m, and the berthing area of SYFP is approximately 460,000 m 2 in fine weather. Under the influence of typhoons of different levels, the waves in the study area presented different distributions. According to the "Master Design Code of Fishery Port (SC/T 9010-2000)", small-and medium-sized fishing vessels can withstand the significant wave heights within 0.5 m, medium and large fishing vessels can bear the significant wave heights within 0.75 m, and large fishing vessels can withstand the significant wave heights within 1 m. However, when the significant wave height exceeds 1 m, the fishing vessels can no longer be berthed. Under the influence of typhoons at level 12 or above, the significant wave height at the entrance of the basin exceeds 1 m (red area), so the fishing vessels cannot be berthed. Large vessels can only be docked anywhere except in the red area, medium and large fishing vessels can only be berthed in the light blue and green area, and small fishing vessels must try to avoid the northernmost end of the pier. At the entrance, there is an area of more than 10,000 m 2 where the significant wave height exceeds 1 m.
When a typhoon strikes, vessels usually adopt the mode of multi-vessel side-by-side and single anchor mooring. At this time, the number of vessels that can be berthed is approximately 158 if the vessels are all large, while it is approximately 735 if all vessels are small. The total number of vessels that can be berthed in the basin using different anchoring ways during fair weather is shown in Table 2. Mooring Area of Sin-Number of Ves- Under the influence of typhoons at level 12 or above, the significant wave height at the entrance of the basin exceeds 1 m (red area), so the fishing vessels cannot be berthed. Large vessels can only be docked anywhere except in the red area, medium and large fishing vessels can only be berthed in the light blue and green area, and small fishing vessels must try to avoid the northernmost end of the pier. At the entrance, there is an area of more than 10,000 m 2 where the significant wave height exceeds 1 m.
When a typhoon strikes, vessels usually adopt the mode of multi-vessel side-by-side and single anchor mooring. At this time, the number of vessels that can be berthed is approximately 158 if the vessels are all large, while it is approximately 735 if all vessels are small. The total number of vessels that can be berthed in the basin using different anchoring ways during fair weather is shown in Table 2. Based on the wave distribution in typhoons of level 12 or above, the significant wave height at the entrance is higher, and therefore, it is not suitable for berthing. The anchoring area of this fishing port is about 460,000 m 2 . Under the effect of different levels of typhoon, the waves in the harbour basin show different distributions. Table 3 shows the number of fishing vessels that can be berthed in SYFP under the influence of different levels of typhoons. It can be estimated that under typhoons of level 12 or above, the anchoring area of small vessels is decreasing. Therefore, the total number of small fishing vessels berthed is approximately 500-600. Due to the good shielding effect at the entrance, the anchorage area of the large vessel is within 450,000 m 2 and the total number of large vessels is approximately 156. The evaluation methods for the typhoon-resistant capability of fishing ports have been summarised in this paper. For the specific vessel type parameters, the typhoon-resistant capacity of the anchorage is improved when multiple vessels berth side-by-side with a single anchor when compared with single anchor mooring. However, if the representing parameters of the vessel type change, the corresponding results will also change.
From the perspective of the layout of the fishing port, the artificial islands and breakwater outside provide a better shielding effect for the harbour. This means that the wave effects weaken gradually as the waves propagate into the fishing port, which is caused by wave reflection and diffraction because of these buildings. Therefore, the arrangement of these structures stabilises the waves in most areas of the SYFP during typhoons [33]. At the same time, due to the influence of the port position, a typhoon moving to the northeast might have a great impact on the port when evaluating the berthing capacity of SYFP during the striking of typhoons. Aside from this, the most adverse wind direction for this port should also be considered.
It should be noted that the berthing capacity of a fishing port is not only related to the characteristics of fishing vessels, but also to the sediment conditions in the fishing port, the typhoon paths, and many other factors [34].

Conclusions
The berthing capability of Sanya Yazhou central fishing port was obtained under level 12-17 typhoons by assuming the WSW direction of the most unfavourable wave direction. Under the predicted typhoons of level 12 to 17, the anchoring area of small vessels at the entrance are: 385,000, 375,000, 3,700,000, 365,000, 350,000, and 315,000 m 2 , respectively. The anchorage area of the large vessel is within 450,000 m 2 .
This paper presents a method for evaluating the berthing capacity of fishing ports under hypothetical typhoons. The typhoon path with the maximum probability was determined based on the data of hypothetical typhoons passing through the port. In terms of the prediction of the starting point, the large sample information meant higher complexity, leading to an unclear existing trend. As such, applying the nonparametric regression method is of great significance in determining such a complex regression relationship. Typhoon wind fields were generated using the Holland model. The MIKE 21 SW model can provide wave boundary conditions for the MIKE 21 BW model by combining the typically designed typhoon field models. Finally, the assessment of the berthing capacity was concluded from the wave field distribution simulated by the MIKE 21 BW model. A study on the framework of a method for hazard assessment was set up for the berthing capacity in the coming typhoon-driven storm waves. In addition, a multi-factor coupling model can be established in the future including the sediment conditions in the fishing port, the typhoon paths, and other factors.