A New Vortex Initialization Scheme Coupled with WRF-ARW

The ability of numerical simulations to predict typhoons has been improved in recent decades. Although the track prediction is satisfactory, the intensity prediction is still far from adequate. Vortex initialization is an efficient method to improve the estimations of the initial conditions for typhoon forecasting. In this paper, a new vortex initialization scheme is developed and evaluated. The scheme requires only observational data of the radius of maximumwind and the max wind speed in addition to the global analysis data. This scheme can also satisfy the vortex boundary conditions, which means that the vortex is continuously merged into the background environment. The scheme has a low computational cost and has the flexibility to adjust the vortex structure. It was evaluated with 3 metrics: track, center sea-level pressure (CSLP), and maximum surface wind speed (MWSP). Simulations were conducted using the WRF-ARW numerical weather prediction model. Super and severe typhoon cases with insufficiently strong initial MWSP were simulated without and with the vortex initialization scheme. The simulation results were compared with the 6-hourly observational data from Hong Kong Observatory (HKO). The vortex initialization scheme improved the intensity (CSLP and MWSP) prediction results. The scheme was also compared with other initialization methods and schemes.


Introduction
To reduce the damage caused by tropical cyclones (TC), it is very important to predict their track and intensity correctly.In recent decades, track prediction has been steadily improved [1].However, there has been little improvement in predicting the intensity of typhoons [1].During 1989-2012, TC intensity errors from the best available models have decreased at 1%-2% yr −1 at 24, 48, and 72 h [2].
As mentioned by Marks et al. [3], the primary reason for the slow progress in improvement of intensity prediction is that the intensity depends on the inner-core dynamics, whereas the track prediction depends more on the large-scale environments.The inner-core dynamics are mainly associated with vortex structure, intensity, and vortex size.It has been suggested that the accurate representation of the inner-core structure of tropical cyclones in the initial conditions is as important as the representation of the largescale environment [4].To improve the representation of the initial vortex, a bogus vortex is often applied because the initial vortices provided by the large-scale analysis are often too weak or misplaced [5].For example, for the super typhoon Genevieve (2014) at 12:00 (UTC) on 9 Aug 2014, the center sea-level pressure (CSLP) derived from FNL (NCEP final analyses) data is around 993 hPa, whereas the observational CSLP from HKO is 940 hPa.
Vortex initialization is an efficient method to provide better estimates of the initial conditions for numerical simulation of tropical cyclones.The advantages of the use of bogus vortices are the ability to customize the vortex structure and low computational cost, compared with dynamical initialization and variational data assimilation methods [6].
In Kwon and Cheong [7], an initialization method with an idealized 3D bogus vortex is developed for track and intensity predictions.In their scheme, a 3D axisymmetric bogus vortex is designed empirically in terms of analytic functions.The bogus vortex is used as a disturbance field and merged into the environmental field.The surface pressure profile of the bogus vortex is given by an exponent decay function.The temperature is derived from the 3D geopotential by using the 2 Advances in Meteorology hydrostatic balance equation.The relative humidity is given by an empirical function, which is not dynamically consistent with the other variables.In our scheme, this problem is overcome by solving the temperature and the humidity simultaneously (Equation (50)).
Reed and Jablonowski [8] used analytic functions to derive the initial conditions.The initial axisymmetric vortex is in hydrostatic and gradient wind balance.The first step is to provide the analytic background specific humidity profile: where   is set as 15 km, which approximates the tropopause height;  0 is the specific humidity at the surface, which is set as 21 g kg −1 ; and   is the specific humidity in the upper atmosphere, which is set as 10 −8 g kg −1 . 1 = 3000 m and  2 = 8000 m are two constants to control the specific humidity profile.The virtual temperature profile is given by where  V0 is the virtual potential temperature at the surface, and Γ is the virtual temperature lapse rate, which is set as 0.007 K m −1 . V =  V0 − Γ  is the virtual temperature in the upper atmosphere.The background temperature profile can then be derived as The background pressure profile () is computed using the hydrostatic equation and the ideal gas law.To obtain the perturbation profile of temperature, an empirical exponent decay function is used to modify the pressure perturbation.Then, by the hydrostatic equation and the ideal gas law, the temperature is derived.Finally, the tangential velocity V  (, ) is derived by the gradient wind balance, which is given by where   = 2Ω sin(  ) is the Coriolis parameter at the constant latitude   and Ω = 7.292 × 10 −5 s −1 is the rotational speed of Earth.  is set to be 10 ∘ . is distance to the TC center.  = 287.04J kg −1 K −1 is the ideal gas constant for dry air. V and  are virtual temperature and pressure.The analytic formulation of this scheme enables it to be applied in various weather and climate models.However, the specific humidity is constructed as a horizontally homogeneous variable, which may not be realistic for the tropical cyclones.Also, the specific humidity is not adjustable, which means that the same specific humidity field must be used for all the tropical cyclones.Therefore, the scheme may be lacking flexibility to simulate different tropical cyclones.In Rappin et al. [6], a highly configurable vortex initialization scheme is developed, which can be easily manipulated to generate different initial vortex structures.The vortex construction starts with parameterizations of the radial wind structure.The tangential wind () of the MR (modified Rankine) vortex is given by where  is the decaying factor, RMW is the radius of maximum wind,  max is the tangential velocity at RMW, and  is the distance from the TC center.
The vertical structure of tangential wind speed follows Stern and Nolan [9]: where  is set as 2 and  is set as 3175 m below  max and 4762.5 m above  max [9]. max is the height of the boundary layer top.Finally, the hydrostatic balance and gradient wind balance are used to determine temperature and pressure.One advantage of this scheme is that it makes use of ( 6), which considers the vertical structure of tangential wind speed.By using this equation, the 3D wind field can be derived from a 2D wind structure (e.g., the MR profile) directly, which makes the scheme highly configurable.The equations for constructing horizontal and vertical profiles introduced in Rappin et al. [6] may also be helpful in designing new schemes.In our new scheme, we further improve the wind construction by making the decaying factor and RMW as functions of .
In this paper, a vortex initialization scheme is developed for the weather research and forecasting (WRF) model with the advanced research WRF (ARW) core.Further, in our scheme, the variables satisfy geophysical formulas, which ensures that the variables are dynamically consistent.The analytic formulations of this scheme enable it to be conducted efficiently.The problems in other schemes mentioned above can also be solved by our scheme.Notice that these problems are mainly associated with two aspects: (1) the variables should be dynamically consistent with each other; (2) the TCs constructed by a scheme should be adjustable and realistic.
The vortex initialization problem can be viewed as follows: given the data sets of (, V, , , , etc.) from global analysis data (e.g., FNL and GFS (global forecast system) data), how can we correct them to new data ( * , V * ,  * ,  * ,  * , etc.)?Here, , V, , ,  represent EW-wind, NS-wind, pressure, temperature, and specific humidity, respectively; " * " is used to represent corrected data throughout this paper.In most situations, we need to provide some extra data besides the global analysis data to build the new vortex.For example, we may need the max wind speed and RMW.Although providing more data may help us to achieve better initial fields, it is difficult to obtain observational data on the ocean.In our scheme, we aim to construct a 3D vortex that can represent basic features of typhoon's structure.The extra inputs are vortex location, max wind speed, and the radius of the max wind.
The rest of the paper is organized as follows: in Section 2, the vortex initialization scheme is derived; in Section 3, numerical simulations for super and severe typhoons in the northwestern Pacific are conducted to evaluate the performance of the new scheme; the simulation results are discussed in Section 4; and in Section 5, conclusions are given.The wind construction can be disaggregated as size correction, intensity correction, and computation of the intensity correction factor.

Size and Intensity
Correction.The first step is "size correction" [13].If the max wind speed and the RMW are computed by the FNL data directly, the max wind speed is generally smaller than the observational data, and the RMW is generally too large.The typhoon size can be corrected by compressing the model grids. and  * are the distances from a grid point to the typhoon center, before and after size correction.The size correction is conducted by  * 1 ( * , , ) =  (, , ) , (if  * ≤  *  ) ; where  * 1 and V * 1 are the wind components after size correction. is the azimuthal angle and  is the height from the ground.This means that the wind speed with distance  is moved to the position with distance  * .In particular, the wind speeds with distance   (RMW) and   (vortex size) are moved to new positions with distance  *  and  *  , respectively.The relationship between  and  * is given by the assumption of linear compression as follows: It can deduced that In practice, the vortex boundary is defined as the contour line of a sea-level pressure value (e.g., 1013 hPa).After size correction, the new vortex is called  * 1 -vortex.The second step is "intensity correction" [13].As mentioned by Lord [14], the vortex initialization is usually carried out by first implanting a synthetic vortex into the large-scale analysis.Similarly, the intensity correction in this paper is conducted by superposing an axisymmetric vortex onto  * 1 -vortex.The axisymmetric vortex is called  2vortex.When adding  2 -vortex to  * 1 -vortex, the intensity is adjusted through the intensity correction factor.In the next subsection, the issues related to the computation of the intensity correction factor are addressed.
In the axisymmetric  2 -vortex, the tangential velocity is (, ), and the radial velocity is zero.(, ) is defined as follows: denote   () as the max tangential wind speed at height  in the axisymmetric vortex, which is the same as the max wind magnitude in the original vortex and  * 1 -vortex.(, ) is the tangential velocity on the circle of radius  (the center of the circle is the TC center).(, ) is defined in the form of an MR vortex: where () is the profile decay rate, which is determined later.For variable , define The physical meaning of ⟨⟩ is the azimuthal mean of variable .By using this symbol, define Then, the decay rate can be deduced: Azimuthal means are used to compute () because this is more representative than the use of single-point estimations.Then,  2 (, , ) = −(, ) sin(), V 2 (, , ) = (, ) cos().Here,  is the angle between -axis and the radial directional Advances in Meteorology vector;  2 and V 2 are the EW-and NS-wind components of  2 -vortex.In the final step, the wind is corrected by  * (, , ) =  * 1 (, , ) +  (, )  2 (, , ) ; where  * and V * are the functions of corrected wind and  is the intensity correction factor.

Computation of the Intensity Correction Factor.
In this subsection, the function of  (the intensity correction factor) is determined.To simplify notations, let  1 () = ( *  (), ),  2 () = ( *  (), ), and The following equation shows that the max wind speed at RMW is intensified to be the same as the observational max wind: where  *  () is the corresponding angle at RMW with height .The vertical profile of  obs () is specified later.Then the following equation can be deduced: The next step is to deduce  2 (), which is the intensity correction factor at  *  ().The reasoning is based on the fact that the typhoon vortex should be merged with the environment continuously.In other words, the vortex boundary values should be close to those of the large-scale background.Mathematically, on the vortex boundary of a typhoon, all of the variables should be continuous.Now, consider a radial direction with the corresponding angle   .By the boundary continuity, lim Note that lim → * −   * (,   , ) =  * ( *  ,   , ), which is the EW-wind component at the vortex boundary after correction.On the other hand, lim is the EW-wind component at the vortex boundary before correction.The same arguments can be applied to V. Therefore, at the boundary of the vortex, the wind magnitude after the change should be equal to the wind magnitude before the change.Choose point  0 on the vortex boundary before size correction.After size correction,  0 becomes .Consider the wind speed at : This can be expanded as Therefore, Because  and  2 are known, the above equation for  2 can be solved numerically.Here, an analytic approximation for  2 is introduced, which can make the computations more efficient.
where   is the tangential wind at , which is counterclockwise.  is the radial wind at , which is outward.In the above approximations, two facts are applied: (1) the wind's magnitude changes a little as  changes; and (2) the radial wind is small, compared with the tangential wind.Similarly, V 2 ( *  ,   , ) = V(  ,   , ).Therefore, In practice, azimuthal means are used to estimate  2 , and the approximation formula is as follows: The above estimation is deduced under an assumption of wind continuity.Further conclusions, considering the continuity of other variables, are demonstrated in later sections.Finally, the formula of (, ) is given.Because of sparse data, let  be a piecewise linear function: By following all the steps above, the 2D wind fields  * and V * can be constructed if the observational max wind speed is provided.In later sections, the vertical profile of  obs is demonstrated and the 3D wind field is deduced.

3D Wind Correction.
The inputs to correct the vertical layer at height  are  *  (),  *  (), and  obs ().Once these three profiles are provided, the 3D wind fields can be corrected.In practice,  *  () =   (), which can be deduced from the continuous vortex boundary conditions in later sections.In this section, the focus is on solving  *  () and  obs ().In Rappin et al. [6] and Stern and Nolan [9], the vertical wind structure is constructed as where () is the tangential wind speed at  and  max , (, ) is the tangential wind speed at  and ,  max is the height of the maximum wind,  is the vertical depth scale, and  is the parameter that controls the decay rate.This profile is an important result in Stern and Nolan [15].Change (25) to The reasons for this are as follows: (1) the radial wind is smaller than the tangential wind, so the magnitude of the horizontal wind is dominated by the tangential wind; (2) in a typhoon, the radial wind above 3 km is very small (typically less than 1 m s −1 ).Notice that new parameters introduced here are  obs ( max ), , and  (a height related to the decay rate). is set as 2 in Stern and Nolan [9].By estimating the observational data, Stern and Nolan [9] set parameter  as 3175 m below  max and 4762.5 m above  max .Stern and Nolan [15] also confirmed that the absolute angular momentum at the RMW for different heights was almost a constant.Now use  obs () to represent the tangential wind in reverse: Therefore, The 3D wind can be corrected from ( 26) and (28).

3D Pressure Correction.
The CSLP is an important metric, which is always used to evaluate the quality of tropical cyclone simulations.In this section, the gradient wind balance is used to deduce the pressure from the wind.This method is especially efficient when the wind fields are constructed analytically.By the gradient wind balance of the governing equation for the axisymmetric component, the following approximation can be deduced: where  is constructed with the same max wind speed and vortex boundary value as the original vortex.Therefore, its profile is used as the axisymmetric part of  * 1 -vortex.Because  * =  * 1 + 2 , (1+) is used as the axisymmetric part of the corrected vortex. corresponds to  * 1 , and  corresponds to  2 .These two profiles have the same directions on the same grids, so they can be superimposed.The sea-level pressure correction step is based on the gradient wind stream function (, ), which satisfies where  0 is the Coriolis parameter at the TC center.Then, The new gradient wind stream function is Let   be the environmental sea-level pressure.The profiles of the environmental pressure are given by the following equations: where Θ  is the perturbed virtual temperature, which can be obtained from the outputs of WRF-ARW.Let Δ(, , ) = (, , ) −   () and Δ * (, , ) =  * (, , ) −   () be the pressure perturbations before and after correction, respectively.Assuming that  and  0 are constants, then, This shows that the pressure perturbation is proportional to the gradient wind stream function.The new pressure perturbation is The explicit formulations are (, ) = 0,  >  *  () , (36c) Finally,  * (, , ) = [ * (, )/(, )][(, , )−  ()]+   () is the desired pressure profile.

3D Temperature and Humidity Correction.
In this subsection, a method to construct a 3D virtual temperature field ( * V ) is first introduced [13].The key technique is to divide the total virtual temperature field into an environmental field and a vortex field.The environmental field provides the base state, and the vortex field provides the perturbations.From the hydrostatic equation and the ideal gas law, the following equation is obtained: where   = 287.04J kg −1 K −1 .Let  be the height of the top of the typhoon vortex. 0 is a fixed height.Integrate the above equation along the segment (, ) × [ 0 , ].In the environment, where  0 ,   are the pressures at  0 and the top of the computational domain, respectively. V is the virtual temperature of the environment.In the total field, where "Δ" signifies perturbation in the vortex, and the subscript "0" signifies the value at  0 .Δ V is the enhancement made to the original environment by the vortex part.Then the equation derived in the total field is modified as Use this equation and the equation for the environment, and an approximation can be obtained for the perturbation of the virtual temperature: After correction, Define Γ(, ) =  * (, )/(, ).Recall that Δ * 0 = Γ(,  0 )Δ 0 .Therefore, For a fixed  0 , Δ * V () = Γ(,  0 )Δ V () for any  ∈ [ 0 , ] is one simple solution.In particular, Δ * V () = Γ(, )Δ V ().Now,  * V is known.Between virtual temperature and temperature, they are connected by where  = 0.608.During the corrections,  * / *  = /  is assumed, where  and   () are vapor pressure and saturation vapor pressure.By definition,  = (0.622)/( − 0.378), which is the specific humidity.Then,  * = ( The continuity of virtual temperature can also be deduced from the continuity of temperature and humidity: Equation ( 51) is equivalent to lim Because lim → * −   * (, ) = 0, lim → * −  (, ) = 0, lim Therefore,  2 = 0 is a necessary and sufficient condition to make lim From the continuity of wind, it can be deduced that *  () =   () can guarantee the above equations.The continuity of other variables can be deduced from the continuity of pressure and temperature.In conclusion, the vortex boundary conditions for all of the variables require that  *  () =   ().

A Real Case Example: Super Typhoon Phanfone (2014).
The super typhoon Phanfone ( 2014) is used as an example to illustrate the effects of the vortex initialization scheme on the initial conditions.The initial time is 00:00:00 (UTC) on 1 Oct, 2014.In Figures 1-4, the surface wind speed WSP  , surface pressure   , surface temperature   , and surface water vapor mixing ratio   before and after initialization are compared in the parent domain (D01), the middle domain (D02), and the inner domain (D03).The parent domain (D01) has a horizontal resolution of 27 km, with 199 × 172 grids; the middle nest domain (D02) has a resolution of 9 km, with 301 × 223 grids; and the inner domain (D03) has a resolution of 3 km, with 469 × 370 grids.The domains are shown in Figure 5.The variables in the inner domain are identically interpolated back to the middle domain and the parent domain.The modified RMW  *  ( max ) is set as 90 km (the sensitivity of this parameter is discussed later); the observational max wind speed  obs ( max ) is set as 60 kt.After initialization, strong wind around the typhoon center can be observed in Figure 1(f).Before initialization, the eye wall of the typhoon is not formed (Figure 1(e)).At the initial time, the observational CSLP is 975 hPa, whereas the CSLP before initialization is higher than 990 hPa (Figure 2(e)), which needs to be intensified.As shown in Figures 3 and 4, the temperature and moisture fields at the initial time are both intensified simultaneously.2014), because the initial MWSP derived from FNL is not strong enough for these cases.The cases are listed in chronological order, and the simulation periods are also listed in Table 1.The simulations start half a day before the typhoon stage.The duration for most simulations is 96 hours, except for Genevieve (for which the observational data starts from the 180th meridian) and Wutip (which ends with dissipation).The simulation results are compared with the observational data from the Hong Kong Observatory (HKO).

Typhoon Cases and Evaluation
The MAE is affected equally by all of the time points during the simulation period and can be used as a variable to evaluate a simulation's overall performance.The improvements of the initialization scheme are evaluated by where "no" signifies no initialization and "init" signifies initialization.In the figures and tables, the units for track, CSLP, and MWSP are km, hPa, and m s −1 , respectively.
For the other cases, the domains are horizontally shifted by changing the reference latitude and reference longitude in the WRF preprocessing system (WPS).For each case, D03 is designed to make the initial vortex located at its right-bottom corner, which helps the inner domain to cover the tracks as far as possible.All of the domains have 39 vertical layers.The Lambert projection method is applied.The time steps are 120 s, 40 s, and 40/3 s for D01, D02, and D03, respectively.To maintain stability, the time step in seconds must be 6 × Δ or less, where Δ is the grid size in km.The current settings for time steps can satisfy the stability constraint and control the computational cost at the same time.The pressure at the top of the domains is set as 50 hPa, which is enough to reach the tropopause height for tropical cyclone simulations.Twoway nesting is used for all of the simulations.Time varying sea surface temperature, sea-ice, vegetation fraction, and albedo are used.No data assimilations are used.For physical parameterizations, the YSU planetary boundary layer scheme is used.The microphysics scheme is the WSM6 scheme [16].
The long-wave and short-wave radiation schemes are the RRTMG schemes [17].The Grell-Freitas cumulus scheme [18] is applied in D01 and D02.The initial conditions and boundary conditions are generated by the FNL data (6-hourly, 1-degree resolution).The initial conditions can be further modified by the vortex initialization scheme, if the initial MWSP is not strong enough.The configurations of all the six cases simulated with the vortex initialization scheme are summarized in Table 2.
Although the observational RMW is currently unavailable, the input parameter RMW in the initialization scheme has a low sensitivity.As shown in Figure 6, the super typhoon Phanfone (2014) was simulated with three RMW values (80, 90, and 100 km).From the time series of track, CSLP, and MWSP, it can be observed that the simulation results are close to each other for different RMW values.

Results and Analysis
For Wutip, Fitow, Krosa, Genevieve, Phanfone, and Nuri, comparisons between without and with initialization are made for track, CSLP, and MWSP.The main metric to evaluate the simulations is the MAE, which is defined in (57).If the FNL data for one case can provide sufficiently strong initial MWSP but cannot provide a low enough CSLP, the vortex initialization scheme will still not be used because our scheme uses the gradient wind balance to deduce the pressure field from the wind field.If the wind field from FNL is stronger than the observation, the intensity correction factor will be negative, which further increases the center pressure.In the following subsections, we analyze the casespecific results and summarize the overall performance.

Analysis of Case-Specific Results.
In Table 3, the improvements of track, CSLP, and MWSP made by the vortex initialization scheme are illustrated for each case.The vortex initialization scheme has little impact on the track prediction results.This is because the track is largely dependent on the large-scale analysis, while vortex initialization can only modify the vortex part.For Fitow and Genevieve, the scheme increases the MAE of track by 4.6% and 3.4%, respectively.For Wutip and Nuri, the scheme results in larger improvements in the MAE of track.In general, the vortex initialization scheme improves the CSLP and MWSP prediction results.The scheme reduces the MAE of CSLP by 7.6%-30.2%and the MAE of MWSP by 7.9%-42.3%.

Analysis of Overall Performance.
In this section, the averaged terms are analyzed, with the averages taken for all the cases.Through the vortex initialization scheme, the averaged MAEs of track, CSLP, and MWSP are improved by 4.0%, 17.6%, and 14.5%, respectively (Table 4).In Table 4, averaged AEs and MAEs are reported.With the exception of 24 h track, 48 h track, and 72 h CSLP, the scheme improves all of the variables (Table 4).Table 5 presents comparisons between our scheme and some other schemes.
In the column "our scheme," a value in bold is the minimum of its row (which means that our scheme provides the best performance of the corresponding metric among all the schemes).
Advances in Meteorology 13 In Xiao et al. [10], the averaged track error of the bogus data assimilation (BDA) scheme is around 110 km, whereas the averaged MAE of track simulated by our scheme is 114.2 km; the averaged CSLP error is around 17 hPa (our scheme: 11.6 hPa, as shown in Table 5); and the averaged MWSP is around 9 m/s (our scheme: 5.7 m/s).The scheme in Kwon and Cheong [7] decreases the 48 h track error to 173.9 km on average (our scheme: 103.4 km); the 48 h error for CSLP is improved from 21.2 hPa to 10.0 hPa on average (our scheme: 15.0 hPa); and the 24 h CSLP error is improved from 37.3 hPa to 11.2 hPa on average (our scheme: 9.2 hPa).In Hsiao et al. [11], the track errors are 208, 95, and 312 km for the with-relocation runs at 24, 48, and 72 h, respectively (our scheme: 69.7, 103.4,and 163.4 km at 24, 48, and 72 h).In Cha and Wang [12], for the MAEs averaged for all forecasts, the dynamic initialization scheme reduces the absolute track error to 103.3 km (our scheme: 114.2 km) and reduces the absolute MWSP error to 7.3 m/s (our scheme: 5.7 m/s).

Summary and Conclusions
In this paper, a new vortex initialization scheme is developed and tested.Through this scheme, a new framework of initial vortex construction is introduced, which starts from 3D wind construction.The other variables, such as pressure, temperature, and humidity, are also constructed on the basis of dynamical and thermodynamical equations.
The scheme requires only observational data for the center location and max wind speed in addition to global analysis data.One advantage of this scheme is that it can estimate a 3D  TC vortex with only little observational data.The derivation of the scheme is mainly based on the dynamical balance constraints.Therefore, the scheme is highly analytic, which makes it satisfy the vortex boundary conditions and with low computational cost.Some highlights of this work include (1) determination of the intensity correction factor at the vortex boundary  2 (23); (2) using the absolute angular momentum at the RMW to derive the RMW function  *  () (28); (3) explicit analytic formulations of the gradient wind stream functions ((36a), (36b), (36c), (37a), (37b), and (37c)), which are useful for improving the computational speed of pressure construction; (4) the construction of temperature and   humidity is unified in (50); and (5) discussion of the vortex boundary continuity (a necessary and sufficient condition of vortex boundary continuity under the framework of this scheme is given in Section 2.5).
Some future work includes improving the parameters, for example,  and  introduced in (6), which is the estimation for vertical profile of tangential wind speed.One potential weakness of the scheme is about the RMW, which cannot be obtained from observations currently.
The scheme was evaluated via simulations for severe and super typhoons in the northwestern Pacific.The scheme can provide a better prediction of typhoon intensity: the averaged MAEs of CSLP and MWSP are improved by 17.6% and 14.5%, respectively.Compared with other methods and initialization schemes, the scheme can provide the best performance for the averaged MAEs of CSLP and MWSP.

Table 1 :
Simulation periods.The variables to be compared are track, CSLP, and maximum surface wind speed (MWSP).The observational data were available every 6 hours.Therefore, the simulated data are compared with observations every 6 hours.Consider the time series of observational data   and simulated data   , 1 ≤  ≤ . is the number of time points to be compared, that is, (length of simulation)/6 + 1.The definition of absolute error (AE) in this paper for -th time point is |  −   |, which is slightly different from the traditional definition (  −   ).The mean absolute error (MAE) is defined as

Table 2 :
Configurations in the vortex initialization scheme.

Table 3 :
Improvements of track, CSLP, and MWSP by the vortex initialization scheme.

Table 4 :
Averaged terms and improvements of track, CSLP, and MWSP by the vortex initialization scheme.

Table 5 :
Comparisons of our scheme with other schemes.