Intercomparison of Empirical Formulations of Maximum Wind Radius in Parametric Tropical Storm Modeling over Zhoushan Archipelago

: The Holland (2010) parametric wind model has been extensively utilized in tropical cyclone and storm surge-related coastal hazard mitigation and management studies. The only remaining input parameter, the radius of maximum wind speed ( R m ), is usually generated by previously proposed empirical relations which are, however, sensitivity to study areas in producing better performed numerical results. In order to acquire optimal R m formulations over the region of Zhoushan Archipelago, East China Sea, 16 empirical relations were compiled into the Holland (2010) model to produce time series of the pressure, wind speed, and wind direction in comparison to observational records taken at three stations during the tropical cyclone events of Ampil and Rumbai. Their respective agreements were evaluated by error metrices including the root mean square error, correlation coefﬁcient, mean bias error, and scatter index, whilst the overall performances of the 16 formulations were ranked according to a proposed comprehensive error. In the following order, the R m formulations of Lu (2012), Zhou (2005), Kato (2018), and Jiang (2008) ranked the best for both events in terms of their minimum comprehensive errors; however, recommendations on the application of speciﬁc empirical formulations for the region of Zhoushan Archipelago are also provided herein from the perspective of conservation and accuracy. Y.B.


Introduction
The East China Sea coast is frequently devastated by typhoon-induced storm surges and disastrous waves during tropical cyclone season each year due to its gentle continental shelf and low-lying coastal landform. Direct economic loss from major historical events steadily increases year on year according to the Bulletin of China Marine Disaster, regardless of the construction of modern levee and revetment along the meandering coastlines. Moreover, the potential catastrophic impact on the coastal community in the future is currently far from mitigated since a northwestward tendency of upcoming typhoon tracks has been plausibly predicted over the Western North Pacific according to representative concentration pathway scenarios [1,2]. Therefore, the accurate and efficient description of wind fields for historical or hypothetical typhoon events plays a fundamental role in the assessment of coastal damage and the planning of a mitigation strategy.
Numerical weather prediction (NWP) models, such as the Fifth-Generation Penn State/NCAR Mesoscale Model (MM5), the Weather Research and Forecasting (WRF) Model, and the Rapid Refresh Model (RAP) have been extensively utilized to predict tropical cyclone tracks and generate high-resolution wind fields at a regional scale [3][4][5].
Under specific boundary conditions, these models solve equations of atmospheric dynamics and thermodynamics while considering a large number of physical processes, including the radiation, surface layer, planet boundary condition, cumulus, microphysics, turbulence, and diffusion to yield the temperature, wind, pressure, precipitation, and other atmospheric variables in support of atmospheric research and operational forecasting [6][7][8]. However, the NWP models are computationally intensive, especially for multi-layer scenarios, and have low bias in modeling storm intensity [9,10]. Their performance is extremely sensitive to not only physical parameters but also the initialization and boundary conditions in the simulation of mesoscale weather events [11,12]. As an immediate result of these model characteristics, the NWP models are especially inefficient for the reconstruction of a large number of historical tropical cyclone events with extratropical transition and landing in probabilistic ensemble hindcasting to fulfill risk assessment. Moreover, without specific boundary conditions in hypothetical scenarios, they are not capable of generating the pressure and wind field of typhoons or hurricanes with a designed track and enforced intensity for engineering applications.
With features such as their low computational cost and compact physical balance, parametric wind models based on empirical formulations and a limited set of parameters serve as an economical solution. This type of model-such as the Rankine vortex model [13], Sea, Lake, and Overland Surges from Hurricanes model [14], Emanuel's model [15,16], and Holland's gradient wind model [17,18]-assumes simplified physical balance which is generalized from tropical cyclone dynamics, and they can reconstruct symmetric or asymmetric surface wind fields of a typhoon or hurricane from the evolution of an axisymmetric stationary vortex with a forward motion. On the premise of ensuring consistency between the radial wind profile and the regional observation data, modifications to the two key parameters, namely the maximum sustained wind speed and the radius of maximum wind, give rise to diverse parametric wind models [19][20][21][22]. These models are broadly applied in studies of hurricane/typhoon-induced storm surges and hazardous waves [23][24][25], natural hazard assessment [26][27][28], hazard mitigation and emergency management [29][30][31], etc.
Among the aforementioned parametric wind models, Holland's gradient wind model [18] (hereafter referred to as H80) is commonly utilized to produce radial axisymmetric wind and pressure distributions through a rectangular hyperbola relationship in many applications involving various ocean-atmospheric processes with tropical cyclone wind forcing [32][33][34]. The H80 model is only determined by two key parameters, namely the radius of maximum winds and the Holland parameter B. The subsequent simplicity of usage and computational efficiency enable this model to reconstruct historical tropical cyclones in various typhoon/hurricane-prone regions with reasonable accuracy by using the information of track, central pressure, and maximum wind speed [3,25,35,36], and then readily generate a large number of statistical scenarios for evaluating the hazardous impacts from hypothetical or ensemble events [37][38][39]. Maintaining the rectangular hyperbolic form of the H80 model, Holland et al. [17,40] (hereafter referred to as H10) introduced an additional exponential parameter x in the wind profile to address the relationship between the surface and gradient layers and improved the agreement between the external wind velocity and the observed data. Therefore, more and more researchers have begun to adopt the H10 model for studies of tropical cyclone-induced natural hazards due to its better performance in comparison to the H80 model [41][42][43][44][45].
In the H10 model, the original scaling parameter B is determined by the pressure drop, central pressure rate change, latitude of cyclone center, and cyclone translation speed. The newly introduced shape parameter x is associated with the pressure drop. Given the track information, central pressure, and the maximum wind speed-usually provided by today's meteorological agencies for historical tropical cyclones-the above variables can be readily obtained, and the only unknown remaining in the computation of the wind profile is the radius of maximum wind speed R m . There are many existing parameterized formulations of R m proposed from statistical analysis or practical experience in different tropical cyclone events. However, it is indicated that these formulations show strong region-dependent characteristics [46], and one equation typically utilized in one region should not be cross-adopted into other regions without careful validation. Consequently, the selection of proper empirical relations of R m becomes crucial in accurately describing the wind profiles of typhoons across East China Sea.
In this study, we implement the H10 model and evaluate the performance of 16 empirical relations of the radius of maximum wind speed in modeling the wind fields of tropical storm Ampil (No.1810) and Rumbia (No.1817) approaching Zhoushan Archipelago, Zhejiang Province, China in 2018. Optimal relations are selected from the rank of comprehensive error and recommendations for the application of specific relations are also advised. In the following, the Holland parametric wind models, the 16 empirical relations of the radius of maximum wind speed, the two selected typhoon events, and the error matric for error analysis are introduced in Section 2. A comparison of pressure, wind speed, and wind direction; error analysis leading to a rank among all 16 relations; and discussion on the implementation of the optimal relations are elaborated in Section 3.

Holland Parametric Wind Model for Tropical Cyclones
The H10 model is built upon the H80 model to depict the radial axisymmetric wind and pressure profiles at the nominal height of 10m above sea surface, with better performance in the description of external wind speed. In this section, we briefly introduce both models to demonstrate their differences in formulation and parameters.
Based on the index pressure distribution model from Schloemer [47], and together with introduction of the Holland parameter B, the pressure profile (in hPa) at the sea surface between the minimum central pressure P c and the peripheral or ambient environment pressure P n in the H80 model is written as where the symbol r denotes the radial distance measured from the typhoon center to the target location, and R m the radius of the maximum wind speed. The associated gradient wind velocity profile is formulated as in which ρ a is air density, and the Coriolis parameter f is expressed in terms of latitude φ as f = 2Ωsin(φ), with the earth rotation rate as Ω. If the maximum wind speed V m is given from best track database or real-time observation, the gradient wind velocity profile can be simplified into The non-dimensional Holland parameter B, usually varies from 0.5 to 2.5, is critical in relating the pressure field to the wind field. Its value can be obtained from parametric relations involving pressure drop, radius of maximum wind speed, or maximum wind velocity [48][49][50][51] and, here, we select the one from Harper and Holland [52], who expressed B in terms of central pressure P c as to complete the formulation of the H80 model. Maintaining the pressure profile as in Equation (1), Holland et al. [17] proposed a redesigned wind profile formulation to better account for horizontal and vertical wind profile of tropical cyclones. The formulation can be written in terms of pressure drop with a newly introduced shape parameter x as or in terms of maximum wind speed V m as where the expression of V m is defined similarly to the one in H80 as ρ ms e (P ns − P cs ) 0.5 (7) in which the symbol e is the base of the natural logarithm, ρ ms is the surface air density ρ s at the radius of the maximum wind speed R m , the subscript s denotes surface values at the nominal height of 10 m, and the shape parameter x will be defined below. The Equations (5) and (6) are equivalent when the exponent 0.5 in the Equation (7) is regarded as the shape parameter x, and the Equation (7) is invoked only when the maximum wind speed is not known beforehand. The surface scaling parameter b s in Equations (5)- (7) functions in a similar way to the Holland parameter B, but considers more impact factors. It is suggested that if the maximum surface wind speed V ms and the surface central pressure P cs can be observed independently, then the surface scaling parameter b s is determined as b s = V 2 ms ρ ms e 100(P ns − P cs ) (8) Otherwise, if remote assessment such as Dvorak analysis is involved in the measurement, the following empirical relation driving by temporal evolution of the surface central pressure P cs is recommended to be b s = −4.4 × 10 −5 (P ns − P cs ) 2 + 0.01(P ns − P cs ) + 0.03 ∂P cs ∂t − 0.014 ϕ + 0.15 v n t + 1.0 (9) in which the symbol ϕ denotes the absolute value of latitude in degree, v t represents the tropical cyclone translation speed in m/s, its exponent n equals 0.6 − 0.0028(P ns − P cs ), and ∂P cs ∂t varies in hPa per hour. If invoking Equations (5), (7) and (8) is necessary in computation, the air density at any location on the surface wind plane can be approximated as where the surface pressure profile P s (r) is similarly defined as the Equation (1) for the H80 model, R = 286.9 J kg −1 K −1 denotes the gas constant for dry air, and the virtual surface temperature T vs measured in K depends on surface temperature T s in • C and surface moisture q s (g/kg) in the following relation Sustainability 2021, 13, 11673

of 23
The surface temperature is estimated by T s = SST − 1, with the SST denoting sea surface temperature; the surface moisture q s is formulated in terms of T s and P s (r) as where the surface relative humidity RH s is usually assumed to be 0.9 if no direct observation is available. With the input of SST and surface pressure profile, the variable density is proved to reduce modeling errors for intense cyclones, but is not critical if a constant value is selected. The shape parameter x used in Equations (5) and (6) is newly introduced in the H10 model, together with the scaling parameter b s , to improve the overall wind profile agreement with observations. The original definition of x forms a piecewise linear function as where an observational wind record measured at the radius of r n is utilized to best fit the computed wind profile in the radial direction with an adjusted exponent x n . The proposed shape parameter x maintains the same value as in the H80 model for the cyclone core section and linearly varies in the external section under calibration with a distant observation. This original definition may become inconvenient when no external wind records are available.
In an attempt to yield a parametric hurricane/typhoon model suitable for forecast and ensemble evaluation, in this paper, we assume x equals 0.6 [1 − (P ns − P cs )/215], which is implemented by Kalourazi et al. [43] and Murty et al. [44] in analyzing tropical cyclones at the Indian coast and Gulf of Mexico and produces reasonable results. Adjusted by inflow angle, which is typically defined by the National Weather Service (NWS), the pressure gradient wind field described by either the Equations (2) and (3) in the H80 model or the Equations (5) and (6) in the H10 model forms the axisymmetric cyclostrophic wind velocity. Additional translation velocity is necessarily superposed to synthesize the wind field of a moving tropical cyclone that is asymmetric in reality. Here, we adopt the translation velocity V t (r) proposed by Jelesnianski [14] in the form of to complete the synthesis.

Empirical Formulations of Maximum Wind Radius
The radius of maximum wind speed is the radial distance from the cyclone center to the location, where the wind speed reaches the maximum of the wind field. It is usually utilized to measure the horizontal scale of the cyclone. It is shown from the governing equations of the H10 model that the radius of maximum wind speed plays a critical role in determining the pressure distribution and wind profile. Since observational records of the radius of the maximum wind speed are difficult to obtain and scarce in number, many empirical relations based on regional tropical cyclone databases are established between R m and other parametersto yield a proper approximation in numerical modeling. In this study, we select 16 formulas proposed from the year of 1959 to 2017 to compare their performance in modeling typhoon near Zhejiang coast. Each formulation is briefly introduced below.

Graham's Empirical Formulation (1959)
Graham and Nunn [53] investigated historical tropical cyclone events across the east coast of the United States and the Gulf of Mexico, compared the influence of central  (15) where the symbol φ stands for geographical latitude, V t the translation speed, and P c the centralpressure of the tropical cyclone. The unit for R m is kilometers, and the one for P c is hPa.

Zhang's Empirical Formulation (2015)
In modeling typhoon events over the South China Sea, Zhang [54] explored the influence of four typhoon parameters on the wind profile across four gradient wind models and confirmed the effectiveness of these models. In particular, among several candidate formulations, he suggested an empirical relationship for Holland models: where the P n denotes the peripheral environment pressure and P c the pressure of cyclone center.

Vickery's Empirical Formulation I (2000a, Vickery-L in Short)
While introducing a new technique for assessing hurricane risk in the United States, Vickery et al. [50] adopted a storm reconstruction method to model the entire path of each hurricane when it crossed the ocean and landed. The method considers the hurricane central pressure as a function of sea surface temperature, and its effectiveness is verified by site measurements. After comparison with the collected data in two regions, the one written in the form as was selected because of its better performance in yielding minimum variance and standard deviation in comparison to the two other proposed formulations.

Vickery's Empirical Formulation II (2000b, Vickery-HRD in Short)
Vickery et al. [55] utilized a dynamic model equipped with a planetary boundary layer (PBL) to calculate hurricane wind fields, which are calibrated by aircraft detection and HRD observation data. Then, they performed a statistical analysis on the characteristic parameters of hurricanes and solved a steady-state wind field on a set of nested rectangular grids using the finite difference scheme. A total of 1560 hurricanes are modeled the with Holland B parameter ranging from 0.5 to 2.5, a pressure drop between 1 and 150 millibars, a radius of the maximum wind speed between 8 and 150 km, and a translational velocity between 0 to 40 m/s. The Fourier fitting method is implemented in the final step to give the proposed R m formulation as where variables follow the prior definitions.

Willoughby's Empirical Formulation (2004)
Based on 493 records of the U.S. NOAA Hurricane Research Division (HRD)'s flightlevel archives from 606 instances of flight detection of tropical cyclones over the Atlantic and eastern Pacific Ocean from 1977 to 2000, Willoughby and Ryan [56] conducted a least square analysis to fit the Holland80 model with these data and examined wind profile in relation to key parameters' distribution. It was found through a frequency distribution histogram of the maximum wind speed radius that In R m is a more appropriate dependent variable in depicting the wind profile. Meanwhile, a multivariate normal distribution is formed to obtain an exponential relationship between the radius of maximum wind speed and the maximum wind speed: where the flight-level maximum wind speed V f max can be adjusted to surface level according to Powell and Black [57] for study in this paper.

Kato's Empirical Formulation (2005)
Kato [58] evaluated the mean and standard deviation of the radius of maximum wind speed in comparison to the previously measured typhoon wind field data and assessed the associated storm surge impact along Japan coast. It was observed that when the central pressure is lower than 950 hpa, the radius of maximum wind speed equals 0.769 P c − 650.55, and when the central pressure is higher than 950 hpa, the R m equals 1.633 P c − 1471.35. The relation can be simplified as which is applied in the following study.

Powell's Empirical Formulation (2005)
Taking into account NOAA HRD's surface wind real-time analysis data archive from 1995 to 2002, HDR's aircraft observation archive from 1984 to 1987 for the west coast of the U.S. and Mexico, the National Hurricane Center's best track archive from 1988 to 1999, and storm archive data collected from 1900 to 1983, Powell et al. [48] conducted a multi-variable analysis and the observed logarithm of the radius of maximum wind speed follows normal distribution. An empirical formulation is proposed to relate R m with pressure drop and latitude in terms of which has been implemented in a parametric wind model developed for the state of Florida to make probabilistic risk assessments on insured residential properties damaged by hurricanes.

Knaff's Empirical Formulation (2007)
Based on the measured radii of maximum wind speed in different tropical cyclone events archived in the U.S. National Hurricane Center (NHC) and Joint Typhoon Warning Center (JTWC), Knaff et al. [59] assumed free parameters in the modified ranking vortex equation as functions of parameters available in the best track data, namely, latitude, cyclone translation speed, maximum wind speed, and the established CLIPER model for continuously computing R m at different latitudes. The proposed formulation contains three constants which vary with the specific latitude of input location, and the formulation suitable for the study area here is where the three constants m 0 , m 1 , and m 2 are 38, −0.1167 and −0.004, respectively.

Li's Empirical Formulation (2007)
Restrained by limited and discontinuous observational data recorded in non-uniform format along China coast, statistical analysis is of great importance for the determination of typhoon parameters for the specific study region because of its targeted and restrictive characteristics. Li [60] applied nonparametric tests, including the Kolmogorov-Smirnov test (K-S test) and Chi-square test, on the probabilistic distribution of critical typhoon parameters, and proposed the radius of maximum wind speed for typhoon-prone area along the south and east coasts of China as (23) in relation to pressure drop through least square regression.

Jiang's Empirical Formulation (2008)
Atkinson and Holiday (1977) analyzed 28 years of maximum wind records measured at coastal and island stations over the western North Pacific and proposed an exponential relationship between central pressure and maximum wind speed. Referring to structure of the Atkinson and Holiday's formulation [61], Jiang [62] investigated tropical cyclone events from 1949 to 2002 archived in the China National Typhoon Yearbook and proposed a similar relationship with the modified parameters in a power exponential form as which is validated by tropical cyclone No.198710 and 200601 in wind field and wave field modeling.

Lu's Empirical Formulation (2012)
Landing at Zhejiang province, China, on 10 August 2006, the super typhoon Saomai (international designation: 0608) devastated the east China coast with 60 m/s maximum wind speed and caused hundreds of deaths and about 20 billion RMB in direct economic loss. Lu [63] used the WRF-ARW model with 1km spatial resolution and modeled the evolution of the pressure distribution and wind field for a duration of 37 h. Statistically validated by onsite measurements, the continuous model output allows the establishment of a relationship between the radius of maximum wind speed and pressure drop in the form as which facilitate development of the Holland B parameter in Lu's study.

Lin's Empirical Formulation (2013)
Based on the observational records of 26 historical typhoon events monitored by ground-based meteorological stations at Hainan province, China, Lin and Fang [64] determined an optimal empirical relationship for computing the B parameter in the H80 model and then analyzed the characteristics of its spatial distribution over the Northwest Pacific to facilitate the application of the H80 model in China. To provide the necessary input of the maximum wind speed radius for computing the B parameter, they compiled 6426 records of tropical cyclones events between the years of 2001 and 2009 archived by Joint Typhoon Warning Center (JTWC) and found a relation between radius of maximum wind speed and pressure drop in the form of R m = −18.04 ln(P n − P c ) + 110.22 (26) through regression analysis.

Kan's Empirical Formulation (2014)
To evaluate storm surge impact on Putuo District, Zhoushan, China, Kan [65] selected No. 1211 and No. 9806 typhoons as study cases and compared the wind field accuracy between Fujita, Jelesnianski and Meng's models. The radius of maximum wind speed required by those three models are provided by two methods. If the radius of moderate gale (Beaufort wind force scale 7) is available in recorded typhoon information, then one Sustainability 2021, 13, 11673 9 of 23 tenth of the radius is regarded as the maximum wind speed radius. If not, the R m is expressed as a second-degree polynomial equation in terms of central pressure as (27) where the empirical constant R k is usually between 30 to 60, depending on research area.
In the paper, we utilize the Equation (27) with R k = 40 as one of sixteen empirical relations for computing the maximum wind speed radius.
2.2.14. FEMA's Empirical Formulation (2005) Based on the landfalling storm information provided by the U.S. National Weather Service, a statistical model proposed by Vickery et al. [55] describes the relationship between the radius of maximum wind speed and both the central pressure and latitude of the storm center. The U.S. Federal Emergency Management Agency (FEMA) made appropriate improvements to the statistical model, and the revised model provides a series of updates to the maximum wind radius information for selected storms occurring in the Atlantic basin, including Hurricane Mitch (1998), Brett (1999), Floyd (1999), and Gilbert (1988) [66]. The empirical relation derived from the revised model is expressed as which is generally suitable for tropical cyclones with maximum wind speed radius of 8~150 km.

Fang's Empirical Formulation (2017)
Based on the best track records of tropical cyclones in the northwest Pacific Ocean from 2001 to 2015 released by the Joint Typhoon Warning Center (JTWC), Fang et al. [67] removed tropical depression and land-based cyclone data to yield a dataset with 5438 records and investigated the relationship between the maximum wind speed radius and both the geographical latitude and central low pressure by referring to Lin [64]. Additionally, a large number of wind, wave, and satellite altimeter observations in five typhoon events over the East China Sea are compared with numerical model results to statistically yield an empirical relation with 95% confidence as which is determined by pressure drop and latitude.
which is implemented in its proposed integrated model to produce a reasonable storm surge level under different tidal conditions.

Two Typhoon Events over Zhoushan Archipelago in 2018
As a better reflection of current global climate status, only the latest tropical cyclones passing the Zhoushan archipelago are considered in this study in order to train a better performing H10 model. Among a few candidate events which occurred in recent years, we chose only representative events in which time span coincided with our available observational records. As a result, taking place in 2018 and passing through Zhoushan archipelago, two tropical storm events, Ampil and Rumbia, were selected in this study to implement the parametric wind model and compare the measurements. Their tracks passed right between station 1 and 2, making them ideal to capture typhoon characteristics on both sides of their forward direction (Figure 1).

Computational Setup and Metric of Error Analysis
The Holland10 model is utilized to reconstruct the evolution of Ampil from 0 a.m., 18 July to 11 p.m., 24  To give comprehensive measurements on the performance of each formulation of R m in determining wind and pressure fields during the two historical typhoon events, four common metrics in error analysis, including mean error or bias, root mean square error (RMSE), correlation coefficient (COR), mean bias error (Bias), and scatter index (SI), are implemented here:

Mean bias error: Bias
Root mean square error: Correlation coefficient: Scatter index: where x i denotes the ith measured value in an observational record with n elements, y i denotes the ith computed value in a model dataset with n elements, and x and y are their arithmetic means, respectively. The computed values of bias, RMSE, and SI from comparison of pressure and wind speed are considered in identifying the optimal R m formulations, since these three indexes, different from COR, proportionally reflect the accumulation of error and the wind direction is not sensitive to applications of various R m values in the parametric typhoon model.

Results and Discussion
Sixteen empirical formulas of maximum wind speed radius embedded in the H10 model give a spectrum of pressure and wind profiles which are evaluated by observational records at three stations through error analysis to yield the optimal formula suitable for typhoon modeling over Zhoushan archipelago. In this section, time series of atmospheric pressure, wind speed, and wind direction are compared between numerical outputs and measurements together with error metrics. Then, optimal R m formulations are determined through comprehensive error analysis followed by a discussion on the result and overall procedure.

Comparison of Pressure at Observational Stations
Tropical cyclone features clear sky, warm temperature, and low atmospheric pressure in its central region or the eye. The larger the gradient between the central pressure and the ambient pressure, the stronger the wind speed. When the typhoon center transitions into a new location, the R m value of each empirical formulation is computed based on model input and then fed into Equation (1) to generate a pressure profile over the computational domain and pressure values at three stations. Figure 2 shows time series comparisons of pressure between the model results and the measurements at each station for tropical storm Ampil and Rumbia, respectively. The red solid line indicates a continuous fluctuation of the measured pressure over a constant baseline of ambient pressure. The center of tropical storm Ampil passed through Zhoushan archipelago in the middle of 22 July and produced a pressure trough at each station around the same time. The lowest pressure values measured at station 1, 2, and 3 are 986 hPa, 993 hPa, and 990 hPa, respectively. On the way to reaching or leaving the pressure trough, all sixteen R m formulations show a sharp drop or rebound, in contrast to a moderate decrease and a mild recovery in the measurements. Both the 3-or 6-h interval of the center location and the strong interaction between the typhoon and on-land atmospheric system prior to landing are likely to contribute the evident difference. In terms of overall fit to the measurements computed as error metrics in Tables 1-3, Lu's R m formulation performs the best, with the lowest RMSE = 3.02, Bias = 0.9233, and SI = 0.003 and highest COR = 0.85 averaged across three stations, and Zhou's (RMSE = 3.6567, Bias = 2.4333, SI = 0.003, COR = 0.8267) and Kato's (RMSE = 3.7567, Bias = 2.5433, SI = 0.003, COR = 0.8233) formulations rank the second and third, respectively. However, in terms of the extreme value that the R m formulation can reach, Lu tends to overestimate the minimum pressure the most, although this drawback taking place for a short time range does not deteriorate the error metrics much.
The center of tropical cyclone Rumbia passed through Zhoushan archipelago around 17 August. Since the intensity of Rumbia was relatively weak, it produced a relatively broad pressure trough at station 1, with a minimum value of 992 hPa, a less distinct trough at station 2 with a minimum value of 993 hPa, and a normal pressure drop at station 3 with a minimum value of 989 hPa. There was a low-pressure system following typhoon Rumbia causing mild drops in measurements at all three stations after August 20. Among the 16 empirical formulations, Lu's R m again yielded the overall best performance with averaged RMSE = 2.86, COR = 0.93, Bias = 1.5, and negligible difference in SI across all three stations. Zhou's (RMSE = 3.25, COR = 0.91, Bias = 1.89) and Kato's (RMSE = 3.44, COR = 0.91, Bias = 2.2) formulations rank second and third best, respectively, in reproducing atmospheric pressure in the event of Rumbia.
formulation performs the best, with the lowest RMSE = 3.02, Bias = 0.9233, and SI = 0.003 and highest COR = 0.85 averaged across three stations, and Zhou's (RMSE = 3.6567, Bias = 2.4333, SI = 0.003, COR = 0.8267) and Kato's (RMSE = 3.7567, Bias = 2.5433, SI = 0.003, COR = 0.8233) formulations rank the second and third, respectively. However, in terms of the extreme value that the formulation can reach, Lu tends to overestimate the minimum pressure the most, although this drawback taking place for a short time range does not deteriorate the error metrics much. The center of tropical cyclone Rumbia passed through Zhoushan archipelago around 17 August. Since the intensity of Rumbia was relatively weak, it produced a relatively broad pressure trough at station 1, with a minimum value of 992 hPa, a less distinct trough at station 2 with a minimum value of 993 hPa, and a normal pressure drop at station 3 with a minimum value of 989 hPa. There was a low-pressure system following typhoon Rumbia causing mild drops in measurements at all three stations after August 20. Among the 16 empirical formulations, Lu's again yielded the overall best performance with averaged RMSE = 2.86, COR = 0.93, Bias = 1.5, and negligible difference in SI across all three stations. Zhou's (RMSE = 3.25, COR = 0.91, Bias = 1.89) and Kato's (RMSE = 3.44, COR = 0.91, Bias = 2.2) formulations rank second and third best, respectively, in reproducing atmospheric pressure in the event of Rumbia.   . Lower central pressure and steeper pressure gradient can usually yield higher wind speed. Tropical cyclones Ampil and Rumbia belong to the scale of the tropical storms, and the wind magnitude of Ampil is slightly higher than the one of Rumbia, owing to its lower minimum central pressure. In the Holland10 model, the radius of maximum wind speed R m not only influences the wind profile from the pressure field through Equation (1), but also modifies the wind profile through Equations (5) and (6) directly.  Figure 3 shows time series of wind speed computed from the H10 model with 16 R m formulations input, respectively, in comparison to measurements at station 1, 2, and 3. In the event of Ampil, the maximum values of wind speed measured at station 1, 2, and 3 were 29 m/s, 21 m/s, and 21 m/s, respectively. The peaks in the time series of the measurements arrived several hours early at all three stations compared to those in the model. This time lag can also be observed from the time series of pressure comparison in Figure 2(a1-a3). The systematic error is possibly associated with less accurate track location and a deformation of wind profile from circular to elliptical shape before landfall. After reaching the peak, the measured wind speed does not recover to the environmental level at the same rate of the computed results because of the involvement of other background atmospheric processes. The computed error indices listed in Tables 4-6   In the event of Rumbia, the measured wind speed slowly reduced to about 10 m/s after the peak arrival at around 0 o clock on August 17 and then was sustained due to the involvement of another atmospheric process; this trend is consistent with the measured pressure shown in Figure 2(b1-b3). Considering the process of the tropical cyclone alone, the modelled wind speed falls quicker than the measurements. A systematic sudden drop in the modelled wind speed is probably related to a 120 • swerve of the center location close to station 3 (Figure 1d

Comparison of Wind Direction at Observational Stations
Wind direction, defined as the arctangent of the east/west component vs. the north/south component, is a sensitive parameter with many complicated influencing factors. The tracks of tropical storms Ampil and Rumbia passed right between station 1 and station 2 (or 3) and produced distinct characteristics of wind direction variation on measurements from each group of stations (Figure 1). In the event of Ampil, station 1 experienced south wind which transitioned into east wind counterclockwise around 22 July and then quickly swung back to SSW clockwise with undulations around 24 July. Stations 2 and 3 also experienced southern wind in the beginning, but it turned to north wind counterclockwise on about 22 July and then continued the rotation to SSW with direction variation on 24 July. Under the same environment setup, all 16 R m formulations contribute to the east/west and north/south components equally and generate basically the same variation in wind direction. The evident direction differences between the modeled and the measured prior to the tropical storm arrival are due to the dominant background wind field, which is not considered in the H10 model. In the event of Rumbia, station 1 experienced eastern wind which transitioned into south or SSW wind clockwise. Meanwhile, station 2 and 3 experienced northern wind on August 16 and then the wind direction turned to south counterclockwise. All 16 R m formulations yield the same numerical results, which basically follow the general trend of the measured wind direction with visible error in different time spans. Since all R m formulations perform the same in the time series comparison in Figure 4, we do not consider wind direction in the following error analysis. measurements arrived several hours early at all three stations compared to those in the model. This time lag can also be observed from the time series of pressure comparison in Figure 2(a1-a3). The systematic error is possibly associated with less accurate track location and a deformation of wind profile from circular to elliptical shape before landfall. After reaching the peak, the measured wind speed does not recover to the environmental level at the same rate of the computed results because of the involvement of other background atmospheric processes. The computed error indices listed in Tables 4-6

Error Analysis
The 16 R m formulations generate a spectrum of time series of pressure and wind speed at three stations through the running of the H10 model for tropical cyclone events Ampil and Rumbia to be compared with the observational records. Their performance is evaluated by four indices: RMSE, COR, Bias, and SI (Tables 1-6). The COR measures the similarity in trend in the time series form, and low correlation does not necessarily mean high deviation in individual magnitude. Moreover, the time series generated by the 16 R m formulations shown in Figures 2 and 3 mainly differ in magnitude, not in trend, and the COR cannot efficiently reflect their performance. Therefore, we do not consider the index COR in the following error analysis. Among the three indices remaining, the value of SI is too small and is not at the same magnitude as RMSE and Bias. To produce the three indices at the same amount of weight in evaluating, we multiply 1000 for SI in pressure comparison and 10 in wind speed comparison. Therefore, for each event (Ampil or Rumbia), and each physical variable (pressure or wind speed), the error listed in Table 7 for all three stations is calculated from Tables 1-6 as where x = 1000 when SI is measuring pressure and x = 10 when SI is measuring wind speed. If we consider both events and physical variables as having the same amount of weight, the average of the two error values provide a comprehensive error for each event listed in Table 7.
In the Ampil event, Lu's formulation performs the best in reproducing the pressure time series, but falls short of reaching good agreement with the wind speed. The first place in pressure with dominant advantage and seventh place in wind speed give it the top rank in the list of comprehensive error among the 16 R m formulations. Zhou's formulation ranks first in wind speed comparison and second in pressure, leading to a second place in the comprehensive error list. Kato's formulation provides the same best performance as Zhou's formulation in wind speed comparison and ranks third in the pressure comparison, resulting in a third place among all 16 formulations. Jiang's formulation follows immediately after Kato's and ranks fourth place in the comprehensive error list. In the Rumbia event, Lu, Zhou, Kato, and Jiang's formulation rank the first, second, third, and fourth place, respectively, in both the wind speed and pressure comparison, leading to the same positions as in the comprehensive error list. The Ampil and Rumbia events were independent tropical storm events in 2018 and had completely different characteristics in terms of distribution of pressure and the wind speed recorded at three stations. Nevertheless, both events give rise to the same rank in the comprehensive error, which strongly supports the reasonability of the proposed methodology and credibility of the optimal formulation list.

Discussion
The investigation of the evidence in support of the lowest comprehensive error being produced by Lu's R m formulation reveals that the plausible first rank is simply an occasional incident. Since the tropical cyclones Ampil and Rumbia are relatively weak in intensity, the measured pressure and wind speed are heavily saturated by background atmospheric processes, especially on the sections 24 h before and after the peak/trough region, and the characteristics of the tropical storms were distorted to a certain extent. Meanwhile, it so happened that Lu's R m formulation usually yields the best performance on saturated sections. This coincidence results in the best agreement on pressure primarily and leads to the top rank in terms of comprehensive error for Lu's formulation. It can also be implied that the radius of the maximum wind speed from Lu's prediction for Ampil and Rumbia must be significantly larger than the realistic one. From a conservative perspective, Lu's R m formulation is still recommended to outline the worst case in hazard mitigation. Meanwhile, from the perspective of accuracy in modeling physical variables, Zhou, Kato, and Jiang, with realistic predictions on R m and steadily good performance, are recommended to be optimal formulations for the Holland10 modeling of tropical cyclones over Zhoushan archipelago in the East China Sea. results, which basically follow the general trend of the measured wind direction with visible error in different time spans. Since all formulations perform the same in the time series comparison in Figure 4, we do not consider wind direction in the following error analysis. formulations and onsite measurements at three stations for typhoon Ampil and Rumbia: (e1-e3) pressure comparison at stations 1-3 for tropical storm Ampil; (f1-f3) pressure comparison at stations 1-3 for tropical storm Rumbia.

Error Analysis
The 16 formulations generate a spectrum of time series of pressure and wind speed at three stations through the running of the H10 model for tropical cyclone events Ampil and Rumbia to be compared with the observational records. Their performance is evaluated by four indices: RMSE, COR, Bias, and SI (Tables 1-6). The COR measures the similarity in trend in the time series form, and low correlation does not necessarily mean high deviation in individual magnitude. Moreover, the time series generated by the 16 formulations shown in Figures 2 and 3 mainly differ in magnitude, not in trend, and

Conclusions
In comparison to numerical weather prediction models, parametric wind models are more efficient in reconstructing pressure and wind fields of tropical cyclone events and have been extensively adopted in storm surge-related coastal hazards and emergency management. As one of the representative realizations, the Holland10 model, inheriting essentials of its predecessor H80 model, has been implemented more frequently. As the only undetermined input parameter sensitive to interest regions, the radius of maximum wind speed can effectively influence the model performance, but optimal empirical formulations in generating R m are not yet known over the Zhoushan archipelago, in the East China Sea. In this paper, 16 R m formulations are compiled into the H10 model to produce a time series of pressure, wind speed, and wind direction at three observational stations in the historical tropical cyclones events of Ampil and Rumbia, both of which took place in 2018. Based on the recorded data at three stations, the performance of numerical results is evaluated through error metrices, including RMSE, Bias, COR and SI, and optimal R m formulations are selected and recommended in terms of the opposite perspective.
The time series comparison in pressure and wind speed among all 16 proposed R m empirical formulations indicates that the radius of maximum wind speed mainly alters the magnitude, instead of the variation pattern. Additionally, the R m computed from different formulations yields the same wind direction, owing to the same contribution on each velocity component. Based on the comprehensive error synthesized from RMSE, Bias, and SI, Lu's, Zhao's, and Kato's formulations yield the top three minimum error rankings in comparison to onsite measurements and Jiang's formulation ranks the fourth. However, the relatively weak nature of the tropical cyclone events Ampil and Rumbia induces saturation into the background environment in measurements, and the R m performance evaluated under this baseline should be perceived cautiously. In Lu's formulation, the radius of maximum wind speed is expressed in terms of pressure drop in an exponential form, and the corresponding time series at each station usually constitutes the lower bound in pressure and upper bound in wind speed. The corresponding feature indicates that Lu's formulation provides the most extreme prediction on the radius of maximum wind speed. It is recommended that Lu's formulation can be applied to produce the most conservative results, and Zhao's, Kato's, and Jiang's formulations are suitable to depict the tropical cyclone more accurately. Data Availability Statement: Data used in the study cannot be publicized due to copyright restrictions.