Objective Center-Finding Algorithm for Tropical Cyclones in Numerical Models

: Precise center-detection of tropical cyclones (TCs) is critical for dynamic analysis in high resolution model data. The existence of both smaller scale perturbations and larger scale circulations could reduce the accuracy of center positioning. In this study, an objective center-ﬁnding algorithm is developed based on a two-dimensional Fourier ﬁlter and a vorticity centroid algorithm. This proposed algorithm is able to automatically adjust its parameters according to the scale of the target vortex instead of using artiﬁcially prescribed parameters in previous research. What’s more, this new algorithm has been optimized and validated by a hundred idealized vortexes with different sizes and small-scale perturbations. A high-resolution simulation of Typhoon Soudelor (2015) was used to evaluate the performance of the new algorithm, and the proposed objective center-ﬁnding algorithm was found able to detect a precise and reliable center.


Introduction
Detecting and tracking the rotational center of tropical cyclones (TC) are critical for objective TC-intensity estimation [1,2], the detection of rapid intensification [3], dynamic diagnosis [4][5][6] and TC climatological descriptions [7][8][9][10][11]. Studies have shown that the overall conclusions could be heavily sensitive to the accuracy of the center [12,13], because positioning the center is the initial step in most of the applications. Thus, the need has arisen for the improvement and sophistication of the center-finding algorithm.
At present, most algorithms used for detecting TC center are based on Doppler radar [14], visible and infrared imageries [13,[15][16][17], microwave radiometers and synthetic aperture radar (SAR) [18][19][20]. For instance, Wood [21] estimated the rotational center by fitting the Doppler velocity field to an idealized Rankine vortex. Similarly, by fitting multiple-Doppler radar or model data to a simple analytical model, Potvin et al. [22,23] presented an objective algorithm to detect mesoscale convective vortices. Lee et al. proposed the velocity track display (VTD) technique to extract real-time tropical cyclone circulations using a single Doppler radar and retrieved the kinematic structure of Typhoon Alex (1987) with an improved VTD-simplex algorithm [24,25]. Detecting methods applied to satellite images were designed based on the structure and organization of cloud features. Wimmers and Velden constructed an automated rotational center hurricane eye retrieval (ARCHER) method using a combination of the brightness temperature gradients in the TC banding patterns and the ring-shaped edge of a possible eye [18]. Recently, they improved the ARCHER algorithm by considering a larger set of satellite based imageries [13]. Zheng et al. proposed two newly developed algorithms by using different image operators to extract eyes from SAR images and found that the TC centers from best track (BT) data were closer to the locations extracted from SAR images than those from infrared (IR) images [26]. However, there are potential uncertainties in the center position calculated by algorithms of the VTD family when asymmetric convections embedded in the primary TC circulation. A larger Fourier filter, truncated wave size and searching radius of the centroid algorithm. This algorithm is designed to automatically adjust these parameters according to the scale of a target vortex, rather than using prescribed settings.
The rest of this paper is organized as follows. Section 2 introduces the components of the new algorithm. The construction of idealized vortexes regarded as benchmarks for optimizing and validating the algorithm is described in Section 3. The new algorithm is optimized and validated using constructed benchmarks in Section 4. The new algorithm is then applied to a high-resolution model simulation of Typhoon Soudelor (2015) in Section 5. Conclusions are given at last in Section 6.

Description of the Objective Center-Finding Algorithm
The new algorithm aims to find the center of a symmetric vortex in complicated situations such as asymmetric convections at different small scales, large scale circulations and other cyclonic vortex in a comparable scale. According to the low-order model of Potvin et al., the wind fields are decomposed into a modified combined Rankine vortex and a background flow [22,23]. The latter is considered as the sum of a spatially constant flow, linearly sheared flow, and linearly divergent flow in the low-order model. Considering that vorticity induced by a spatially constant flow equals to zeros, and that by a linearly sheared flow is constant, the low-order model could be simplified by a low-order vorticity model. Therefore, this new algorithm prefers to vorticity field instead of wind fields to detect center.
The vorticity field containing at least a cyclonic vortex could be decomposed into three parts, a modified Rankine vortex (symmetric vortex), asymmetric perturbations and vorticity from large scale circulations. Vorticity from large scale circulation is supposed constant within a 'small' area. In other words, it could be eliminated by removing the domain mean. Meanwhile, the asymmetric vorticity perturbations in a TC vortex are supposed to be at smaller scale than the corresponding symmetric part. Therefore, retrieving the major vorticity field of the symmetric vortex by eliminating asymmetric perturbations is the goal of the objective center-finding algorithm. Besides, this new algorithm is also designed to be able to detect multiple TCs even when they are close to each other. The procedure is as follows.

•
Step 1. Detection potential cyclonic vortex Detecting potential vortex is the fundamental step towards finding the location of centers in the existence of multiple cyclonic vortexes. Similar to the criteria used in previous research, the grid point whose smoothed pressure or geopotential height is lower than its eight surroundings is determined as a candidate for a potential cyclonic vortex center [41,42]. For simplicity, the field is smoothed by a nine-point averaged operator. Meanwhile, the radius of maximum azimuthal mean tangential wind (RMMTW) and maximum azimuthal tangential wind (MMTW) are calculated for each candidate. A candidate will be disqualified once if its MMTW is found lower than a prescribed value (for instance, category 1 vortex defined by the Saffir-Simpson hurricane wind scale was used in this manuscript). The precise center position of potential vortex will be detected by the following two steps, individually.

•
Step 2. Removing small scale perturbations Vorticity in a square area centered at each candidate is selected, and it is on this area where a two-dimensional Fourier transformation is carried on. The domain size varies among different candidates as it is supposed to be relative to the estimated RMW of the target vortex. The Fourier coefficients are computed as follows: where F is the Fourier coefficients. M and N are domain sizes in direction x and y, respectively. As the domain is chosen as a square area, M equals to N. The wavenumber in x and y direction is represented by u and v. Besides, the space resolution is supposed to be equal in both directions.
To remove smaller-scale perturbations, Fourier coefficients are set to be 0 where its wave length (λ u,v ) is shorter than a prescribed scale (λ 0 ), where λ u,v is defined as ( M·dx u ) 2 + ( N·dy v ) 2 ). The filtered coefficients are defined in Equation (2).
The filtered vorticity field is obtained by applying an inverting Fourier transformation with the modified Fourier coefficients. The inverted Fourier transformation is described in Equation (3). • Step 3. Determining the vortex center A centroid algorithm is applied to the filtered vorticity filed to determine the final center of the vortex. The centroid algorithm is briefly illustrated in Equation (4), where x and y represent the location of center in x and y direction, respectively. R represents the searching radius, and ξ with a subscript i represents vorticity at each grid within the searching radius. There are only two parameters (an initial guess center and a searching radius R) involved in a centroid algorithm. The guess center is chosen as the position of the candidate calculated in step 1. Meanwhile, the searching radius is objectively defined according to the scale of target vortex which will be illustrate in the Section 4.
This objective center-finding algorithm applies pressure field, wind fields and vorticity field together to determine centers of cyclonic vortexes. By involving the first step, the algorithm is able to detect multiple cyclones even when they are close. A two-dimensional Fourier filter based on automatically adjusted wave-length is applied to eliminate small-scale perturbations, and retrieve the primary symmetric circulation of each target vortex in step 2. Finally, the center is detected by applying a vorticity centroid algorithm in step 3. It should be noted that parameters involved in the new algorithm are automatically adjusted according to the scale of a target vortex rather than artificially defined as constant values.

Construction of Idealized Vortexes
The center of a tropical cyclone is usually defined by the location of minimum wind or minimum pressure. Meanwhile, the local maxima of relative vorticity could also be used to identify the center. As there is no standard definition of a TC's center, and centers defined by a variety of center-finding algorithms are not identical, defining a center under idealized vortexes offer an option to evaluate the performance of center detecting algorithm. Idealized cyclonic circulations were supposed to be constructed by four individual parts in a cylindrical coordinate-the symmetric tangential wind, the symmetric radial wind, environment flow, and the isolated nondivergent asymmetric perturbations. The center of a cyclonic vortex is defined as the center of its symmetric vortex.
The environment flow is supposed to be a summary of a linear sheared flow and a uniform flow which introduces a constant vorticity at each point. It is also supposed that the symmetric radial wind is represented by a lognormal function profile, which has no influence on vorticity. Besides, according to previous research, adding the symmetric radial wind does not affect the computation of RMMTW and MMTW [34]. Therefore, neither the symmetric radial wind or the environment flow is initialized when an idealized vortex is constructed, as they have no influence on the performance of a vorticity based center-finding algorithm. The initialization of the idealized vortex is illustrated as follows.
Firstly, the symmetric vortexes are constructed by vorticity-based profiles. The profile that was used to generate a Rankine-like vortex in previous studies smoothly connects radial intervals of constant vorticity with relatively low values at the eye, and the highest value at the eyewall annulus [43,44]. The formula is described in Equation (5), where S(x) = 1 − 3x 2 + 2x 3 is the cubic Hermitian polynomial. Values for each of the ζ i , r i , d i , r i and d i can be changed to represent a hurricane-like vortex with a variety of strengths and various distributions. It is the radial distribution of vorticity rather than its amplitude that needs to be paid attention in validating the center finding algorithm. Therefore, ζ 1 , ζ 2 and ζ 3 are set to 5 × 10 −4 s −1 , 2.5 × 10 −3 s −1 and 1.0 × 10 −4 s −1 , respectively. The symmetric tangential winds were initialized on a high-resolution cylindrical grid, which was 0.5 km in radius and π/360 in radians. Once initialized, the fields on the cylindrical grid were interpolated onto a Cartesian grid using a bi-cubic interpolation. The Cartesian grid was set as 1 km in space to emulate the finer structure of tropical cyclones. The tangential winds were also transformed into Cartesian winds. It was on this Cartesian grid that the center was detected by center-finding algorithms.
Secondly, to simulate the asymmetric structure of a cyclonic vortex, idealized vortexes are constructed by adding isolated asymmetries to symmetric vortexes. Asymmetries are defined as Gaussian-type perturbations in the streamfunction given as follows: where A is the wave amplitude and b is the radial location of asymmetries relative to the vortex center. c, n, λ, and ϕ represent the radial width, azimuthal wavenumber, azimuth angle, and wave phase, respectively.
To initialize the nondivergent asymmetries, the symmetric vortex was first converted into the streamfunction by inverting the relationship, ζ(r) = ∇ψ(r). The distribution of ψ in the prescribed high-resolution cylindrical grid was numerically computed by supposing the streamfunction at the outer boundary to be 0. Then, the full wind and vorticity were calculated from the total streamfunction filed, including both ψ and ψ , by solving a Laplace equation. Once initialized, the fields on the cylindrical grid were interpolated onto the Cartesian grid.
To generate a variety of different vortexes, parameters controlling the radial location and the range and strength of Gaussian-type asymmetries were 'randomly' chosen under certain restrictions. So does the number of asymmetries. Meanwhile, parameter r 2 of the symmetric vortex is also 'randomly' chosen to generate symmetric vortexes at different scales. The restrictions for each of these parameters are listed in Table 1.  A hundred different vortexes were constructed as benchmarks to optimize and validate the performance of the objective center-finding algorithm. The RMWs of these vortexes varied between 30 km to 45 km. Meanwhile, large differences were found in the amplitude and radial structure of asymmetric perturbations. The first 50 vortexes (hereafter, training group) were used for optimization, and the last vortexes (hereafter, testing group) were used for validation in Section 4.
For instance, the wind structures of one in the hundred vortexes were displayed in Figure 1. A polygon vorticity field was found around the real center. The maximum vorticity was 2.64 × 10 −3 s −1 , and the maximum wind speed was 33.4 m s −1 . Negative perturbations in vorticity were located just outside the RMW. This vortex could represent a category 1 tropical cyclone with small scale convections at RMW. It should be noted that a pressure or geo-potential height field was needed to detect the candidates of potential cyclonic vortexes. However, neither of them was provided by the initialization procedure. Therefore, the streamfunction was used instead of a pressure field.

Algorithm Optimization and Validation
There are three parameters involved in the proposed algorithm: (1) the domain size where two-dimensional Fourier transformation is carried on, (2) the truncated scale used to filter perturbations, and (3) the searching radius for a vorticity centroid algorithm. The object of this section is to find the optimum values of the three parameters with fifty realistic vortexes in the training group. Then, the optimized algorithm is validated and compared with a centroid algorithm with the testing group.

Method
Centers are calculated by a variety of new algorithms with different settings of the three parameters. Table 2 shows the range of each parameter, and a total number of 252 testing algorithms are generated accordingly. It should be noted that the ranges of each parameter could be set larger and finer, the main conclusions remain unchanged.
The center position of each vortex in the training group is detected by the 252 testing algorithms, simultaneously. Considering the size of the vortexes in the training group varies, the three parameters are normalized by the RMW of the target vortex to eliminate difference in the scale of the target vortex.
Statistics of the position error are carried on according to the three normalized parameters rather than their physical values. The medium, 25% and 75% percentiles of the position error are used to evaluate the accuracy, in addition to the group mean and standard deviation. Table 2. Ranges of parameters for generating testing algorithms. Units are in kilometers. The interval of each parameter is displayed in blanks.

The Domain Size
The performance of the proposed algorithm is evaluated by the domain size at first. Figure 2a is the 'boxplot' of position error where the medium value in each group is marked in a red line [45]. The edges of the box are the 25th percentile and 75th percentile, the whiskers extend to the most extreme position error excluding the outliers. Meanwhile, the outliers are defined as follows: where Q1 and Q3 are the 25th and 75th percentiles, respectively. According to Figure 2a, the medium position error decreases from 0.76 km to 0.3 km with an increasing of normalized domain size, and it varies a little when normalized domain size is larger than 8 RMW. The medium position decreases 20% when the domain size increases from 7 RMW to 8 RMW. Meanwhile, the decline rate of medium error is less than 10% once the domain size is lager than 8 RMW. The lowest position error seems to be not sensitive to the domain size as it is close to zero at each group. So does the 25th percentile of position error when the domain size is larger than 8 RMW. The 75th percentile of position error decreases a little when the normalized domain size is larger than 10 RMW. Both the group mean and standard deviation in Figure 2b share a similar distribution with the medium position error. The group mean of position error decreases from 5.5 km to 0.41 km when the domain size increases from 6 RMW to 13 RMW. Besides, the group mean and standard deviation are much larger when domain size is smaller than 7 RMW. Therefore, it can be concluded that the performance of the proposed algorithm is sensitive to choice of the domain size which is also relative to the scale of the symmetric vortex. The larger the domain size is, the better the algorithm performs.
In practice, as the 'true' RMW of the target vortex is unknown, estimated RMW of the target vortex is considered as an alternative. The mean difference between the estimated RMW and the 'true' RMW is found quite small (less than 0.11 RMW), which is consistent with conclusions of Ryglicki et al. [34]. Considering that the large scale circulation is supposed to be uniformed in a 'small' area, it is reasonable that the domain size is set as ten times of the estimated RMW for both the computation efficiency and accuracy.

The Truncated Scale
Position error of the proposed algorithm is examined with the training vortexes based on the truncated scale. As shown in Figure 3a, the medium position error increases from 0.37 km to 0.75 km when the normalized truncated scale increases from 1.5 RMW to 3.5 RMW. Meanwhile, the minimum position error does not vary significantly. The 25th percentile of position error is close to 0.2 km when the normalized truncated scale is less than 2 RMW, and it increases with truncated scale after then. The 75th percentile of position error is close to 0.7 km when the normalized truncated scale is between 1.0 RMW and 2.0 RMW. After that it increases with the growing of the truncated scale. Both the group mean and standard deviation grow with the increasing of the truncated scale once the truncated scale is larger than 1.5 RMW. Therefore, it can be concluded that the performance of the proposed algorithm is sensitive to the choice of truncated scale as well. Perturbations whose wave length are smaller than 1.5 RMW should be eliminated to retrieve the primary symmetric vortex and improve the accuracy of the new algorithm. Therefore, the optimized truncated scale is chosen to be 1.5 times of the estimated RMW of the target vortex in the proposed algorithm.

The Searching Radius
The statistics of the searching radius in the proposed algorithm are shown in Figure 4. As shown in Figure 4a, the medium position error decreases from 1.32 km to nearly 0.3 km when the normalized searching radius increases from 1.25 RMW to 3.75 RMW. The medium position error varies a little when the searching radius is larger than 2.25 RMW. The 25th percentile of position error is less than 0.05 km in all groups. Meanwhile, difference in the 75th percentile of position error is much significant. The 75th percentiles of position error decreases from 3.7 km to 1.2 km when the searching radius increases from 1.25 RMW to 2.75 RMW. Then, it increases with the growing of the searching radius. Besides, the group mean position error decreases with the growing of searching radius until 2.25 RMW in Figure 4b. So does the standard deviation. It should also be noticed that the standard deviation is much larger at 2.75 RMW than that at 2.25 RMW. The reason for this significant increasing in standard deviation is that there are more outliers in the group of 2.75 RMW. Besides, more outliers indicate that the performance of the new algorithm is less reliable when the searching radius is set close to 2.75 RMW. Therefore, it is suggested that the optimized searching radius is chosen as 2.25 times of the estimated RMW in the proposed algorithm.  It should be realized that statistical accuracy of the new algorithm is found to be sensitive to the settings of the three parameters which are found relative to the scale of the target vortex. Therefore, the three parameters are not optimized in their physical values, but in the relationship with the scale of target vortex instead. In this way, the new algorithm could objectively and automatically adjust the physical values of each parameter according to the statistical optimized relationships.

Algorithm Validation
The accuracy of the proposed optimized algorithm (hereafter FVC) is examined with the fifty vortexes in the testing group. Based on the findings in previous sections, the domain size, truncated scale and searching radius in the optimized algorithm are objectively set to 10, 1.5 and 2.25 times of the estimated RMW of the target vortex, respectively. In other words, the physical values of the three parameters are objectively and automatically defined based on the size of the target vortex. A vorticity centroid algorithm (hereafter VC) is chosen as a compared algorithm. In previous research, the domain size of a vorticity centroid algorithm may vary from 60 km to 300 km [27,38]. As the asymmetric perturbations are initialized within 3 times of the RMW of each vortex which is about 150 km, it is reasonable to set 150 km × 150 km as the domain size of the compared algorithm. Figure 5 displays the center locations of fifty testing vortexes detected by the FVC method and VC method. There are 33 centers detected by the FVC whose position error is less than 0.5 km. Meanwhile, only 26 centers detected by VC are found with position error less than 0.5 km. Centers detected by VC method could have larger position error, as the maximum position error in VC is about two times of that in FVC. Table 3 shows statistical results of the two algorithm. The mean position error of FVC is 0.4 km, which is much smaller than that of VC. The standard deviation of the VC is as large as 1.06 km. Meanwhile, the performance of FVC seems to be more reliable as its standard deviation is only 0.354 km. The medium position errors of the two algorithm are close to each other, so does the minimum position errors. Therefore, it can be concluded that the proposed objective center-finding algorithm is better than a vorticity centroid algorithm in both accuracy and reliability.

Application to High-Resolution Model Data
The objective center-finding algorithm is then applied to a high-resolution (2-km horizontal grid spacing) numerical simulation of Typhoon Soudelor (2015). Soudelor formed as a tropical depression near Pohnpei on 29 July 2015. Then it strengthened slowly before starting rapid intensification on 2 August. The storm made its landfall in eastern Taiwan on 8 August. Then, it moved inland and degraded to a tropical depression by 9 August. The proposed algorithm is applied to the model simulations and its performance is compared with other commonly used algorithms which are shown in Table 4. Table 4. Center-finding algorithms used in comparison.

Model Setup
the Advanced Research Weather Research and Forecasting Model (ARW-WRF version 3.3.8) was used to simulate the evolution of Typhoon Soudelor. The model domain was two-way interactive and triply nested with horizontal grid spacings of 18, 6, and 2 km, respectively. The three meshes were set with sizes of 311 × 251, 271 × 271 and 211 × 211. There were 40 σ levels in the vertical, hydrostatic pressure as the vertical coordinate extending from the surface to the model top at 50 hPa. The two nested inner meshes automatically moved following the TC during the model integration, while the outermost 18-km mesh was fixed. Other configurations of the ARW-WRF is displayed in Table 5 similar to those by Wang et al. [46,47]. The model was initialized at 0000 UTC 3 August 2015 and integrated for 120 h up to 0000 UTC 8 August. The Final Analysis (FNL) data with horizontal resolution of 1 • × 1 • was used for model initial and lateral boundary conditions.

Verification of the Simulation
The track and intensity of the simulated Soudelor are compared with those from the JTWC (Joint Typhoon Warning Center) best track data in Figure 6.
It seems that the model simulated the northwest ward movement quite well. Whereas, the ability to simulate the evolution of intensity is not good enough. Figure 7 shows the evolution of maximum wind speed and minimum surface level pressure.
The Soudelor intensified much more rapidly from 0000 UTC to 1200 UTC 3 August in our simulation. After the rapid intensification (RI), the modeled Soudelor maintained as a category 4 hurricane for nearly 12 h (from 1200 UTC 3 to 0000 UTC 4 August), and then decreased to a category 2 hurricane until 0600 UTC 7 August. It re-intensified before making its landfall in Taiwan at 1800 UTC 7 August. The simulated Soudelor moved faster and landed 6 h earlier than the real storm. The simulation was able to capture the RI, decay, re-intensification and landfall of a tropical cyclone. It is not our intention to clarify every details of the simulation, and attentions are paid on the performance of the proposed algorithm.   Figure 8 shows the tracks detected by different center-finding algorithms displayed in Table 4. The centers detected by the six different center-finding algorithms are very close until 0000 UTC 6 August before the re-intensification occured. In other words, the dispersion of center locations is relatively small during the RI, the mature and the decaying period. Whereas, significant differences among the six algorithms exist when the simulated Soudelor re-intensified and made its landfall in Taiwan. Therefore, it seems that it is not the strength of a tropical cyclone that determines the dispersion of the center locations. Besides, the track determined by FVC is found smoother that others.

Performance of the Proposed Algorithm
Six snap shots are selected from the 120 h simulation to further evaluate the performance of the proposed algorithm (FVC). The snap shots are selected to represent the structure of simulated Soudelor at different stages. The surface level pressure, vorticity and winds at 10 m above the sea level are displayed in Figure 9 with the center locations detected by the six algorithms as well.
The detected centers are close to each other as the maximum difference of the center locations is less than 4.2 km in Figure 9a when the Soudelor intensifies rapidly. When Soudelor becomes mature, the dispersion of center locations is also relatively small as the maximum difference is nearly 3.59 km in Figure 9b. The center detected by the FVC algorithm is also found close to others in Figure 9c when the Soudelor decreases to a category 2 tropical cyclone. It is difficult to tell which algorithm is more accuracy from Figure 9a-c as the detected centers are close to each other and they all locate in small wind speed areas. However, the dispersion of center locations is much larger in Figure 9d-f. The vorticity field in Figure 9d,e has a triangle distribution rather than circular distributions in Figure 9a-c. The gradient of SLP is much weaker at inner core as well. The centers detected by Pmin, PC100 and PC150 are found located in areas with wind speed higher than 6 m/s in Figure 9d. Besides, the performance of VC100 is worse as the its center locates in a much higher wind speed area. However, centers detected by VC150 and FVC seems to be more reasonable as they all locate at small wind speed area. Similar to Figure 9d, the centers detected by VC150 and FVC seems to be more reasonable than those detected by other algorithms in Figure 9e. Therefore, the VC150 and FVC are more appropriate when applied to an asymmetric storm with a polygon vorticity field. It is notable that the maximum difference of center locations could be as large as 70 km in Figure 9f when the simulated Soudelor makes its landfall in Taiwan. Due to the influence of topography by Central Mountain of Taiwan, the fields of Soudelor become more complicated. Centers detected by the Pmin, PC100 and PC150 are found north to others. Meanwhile, the VC100 and VC150 detect the centers in the high wind speed areas to the west of the Central Mountain. The center detected by the FVC seems to be more reasonable as it located in a small wind speed area and it appeares to be the center of the large scale circulation. Therefore, it is thought that the FVC is able to detect a reasonable center for a tropical cyclone even when the storm is weak and asymmetric. Figure 9. The horizontal structures of simulated Soudelor including the surface level pressure (black lines), vorticity (shaded) and winds (wind-bars) at 10 m above the sea level. The location of detected centers are represented by different markers displayed at left. (a-f) display the fields at time 0900 UTC 3, 2300 UTC 3, 1700 UTC 5, 1900 UTC 6, 0000 UTC 7, and 2100 UTC 7 August, respectively.

Conclusions
An objective center-finding algorithm was proposed to detect centers of cyclonic vortexes in numerical model datas. The main advantages of the proposed algorithm are summarized as follows: • Pressure, winds and vorticity fields are involved in the new algorithm rather than using a single field in previous algorithms. The participating of multiple fields provides more physical factors of a TC. • A two-dimensional Fourier filter is involved in the new algorithm to eliminate perturbations whose wave length are smaller than a certain threshold.

•
Most importantly, the parameters involved in the new algorithm could automatically change their physical values according to the scale of target vortex.
As there is no standard definition of the center in a real TC, the objective center-finding algorithm was optimized and validated under idealized vortexes with known centers. A hundred idealized vortexes were constructed by combining symmetric vortex with different sizes and asymmetric perturbations with different structures as benchmarks. The first fifty vortexes (training group) were used for optimization, and the rest vortexes (testing group) were used for validation.
Statistics on position error of fifty idealized vortexes in the training group were used to examine the influence of the settings for the three parameters on the accuracy of algorithm. It is found that the performance of the new algorithm is sensitive to the choice of parameters. Meanwhile, the scale of symmetric vortex has an influence on the performance. Using evaluating indicators such as group mean, standard deviation, medium, 25th and 75th percentiles, the performance of the new algorithm was found improved when the domain size, truncated scale and searching radius were set close to 10, 1.5 and 2.75 times of the RMW of the target vortex, respectively. A link between the parameters of a center-finding algorithm and the scale of the target TC is built objectively through statistical analysis. In practice, the estimated RMW of a target vortex calculated in step 1 was used instead of the unknown 'real' RMW. Therefore, the physical values of the three parameters in this new algorithm were objectively and automatically defined according to the scale of the target vortex.
The optimized algorithm is validated using the fifty idealized vortexes in the testing group. The averaged position error was 0.4 km which is much smaller than that of a vorticity centroid algorithm with a 150 km × 150 km domain size. The new algorithm was also proved more reliable as it is smaller in standard deviation and maximum position error.
The performance of the proposed objective center-finding algorithm is examined and compared with other commonly used algorithms in the WRF simulation of Typhoon Soudelor (2015). The center locations detected by different algorithms are close to each other when the storm is symmetric. Whereas, large dispersion of algorithms exists when the storm has a polygon vorticity field or is influenced by topography. The FVC is found to be able to detect a reasonable center for a tropical cyclone at different stages including RI, decaying, re-intensification and making a landfall.
The objective center-finding algorithm is optimized and validated with idealized vortexes where their symmetric parts are constructed stronger than a category 1 vortex. Future study should focus on the performance of the new algorithm when it is applied to real TCs, especially when extra-tropical transformation occurs. As the new algorithm is developed able to detect multiple vortexes, it is necessary to examine its performance by considering the merging of two TCs.

Conflicts of Interest:
The authors declare no conflict of interest.