• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Which Features of the SST Forcing Error Most Likely Disturb the Simulated Intensity of Tropical Cyclones?

    2021-04-07 10:20:56JiaweiYAOWansuoDUANandXiaohaoQIN
    Advances in Atmospheric Sciences 2021年4期

    Jiawei YAO,Wansuo DUAN*,and Xiaohao QIN

    1State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics (LASG),Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029, China

    2University of Chinese Academy of Sciences, Beijing 100049, China

    ABSTRACT Among all of the sources of tropical cyclone (TC) intensity forecast errors,the uncertainty of sea surface temperature(SST) has been shown to play a significant role.In the present study,we determine the SST forcing error that causes the largest simulation error of TC intensity during the entire simulation period by using the WRF model with time-dependent SST forcing.The SST forcing error is represented through the application of a nonlinear forcing singular vector (NFSV)structure.For the selected 12 TC cases,the NFSV-type SST forcing errors have a nearly coherent structure with positive (or negative) SST anomalies located along the track of TCs but are especially concentrated in a particular region.This particular region tends to occur during the specific period of the TCs life cycle when the TCs present relatively strong intensity,but are still intensifying just prior to the mature phase,especially within a TC state exhibiting a strong secondary circulation and very high inertial stability.The SST forcing errors located along the TC track during this time period are verified to have the strongest disturbing effect on TC intensity simulation.Physically,the strong inertial stability of TCs during this time period induces a strong response of the secondary circulation from diabatic heating errors induced by the SST forcing error.Consequently,this significantly influences the subsidence within the warm core in the eye region,which,in turn,leads to significant errors in TC intensity.This physical mechanism explains the formation of NSFV-type SST forcing errors.According to the sensitivity of the NFSV-type SST forcing errors,if one increases the density of SST observations along the TC track and assimilates them to the SST forcing field,the skill of TC intensity simulation generated by the WRF model could be greatly improved.However,this adjustment is most advantageous in improving simulation skill during the time period when TCs become strong but are still intensifying just prior to reaching full maturity.In light of this,the region along the TC track but in the time period of TC movement when the NFSV-type SST forcing errors occur may represent the sensitive area for targeting observation for SST forcing field associated with TC intensity simulation.

    Key words:tropical cyclone,sea surface temperature,sensitivity,error growth,inertial stability

    1.Introduction

    A tropical cyclone (TC) is one of the most destructive synoptic scale systems on earth.It frequently brings great economic loss to coastal areas and even interior locations of countries (Peduzzi et al.,2012).It is therefore important to improve the skill of TC forecasts.With recent improvements of model simulation capability,to include data assimilation methods,great progress has been made resulting in greater skill regarding TC track forecasting.However,forecasting TC intensity still presents a huge challenge.The main reason for this is that TC intensity is not only dependent upon large-scale environmental factors such as vertical wind shear and sea surface temperature (SST) but also upon mesoscale/microscale processes represented by parameterizations such as moist convection,microphysics,etc.(Emanuel et al.,2004;Wang and Wu,2004;Zhang et al.,2011;Hakim,2013;Zhang and Tao,2013;Tao and Zhang,2014;Torn,2016).The uncertainties of these factors operating at different spatial scales and their interactions with one another,complicate and therefore limit the predictability of TC intensity.Emanuel and Zhang (2016) used the Coupled Hurricane Intensity Prediction System (CHIPS;Emanuel et al.,2004) to explore the sources of TC intensity forecasting uncertainties.They found that the TC intensity errors grow during the early stages of development and is dominated by initial intensity errors,while the errors of TC track and environmental shear become more pronounced in affecting TC intensity during later developmental periods.Torn (2016)compared the role of atmospheric uncertainties and oceanic uncertainties in TC intensity simulation error and found that although the atmospheric uncertainties dominate during the initial period,oceanic uncertainties become equally as important during later periods.One may conclude that the inclusion of oceanic uncertainties can help to improve the predictability of TC intensity.

    Previous studies reported the influence of SST on TC intensity either in theories or through observations.On one hand,the ocean is thought to be the energy source for TC genesis and intensification (Riehl,1950;Emanuel,1986,1988;Holland,1997).Specifically,the theory of maximum potential intensity (MPI),proposed by Emanuel (1986),treated TCs as heat engine.The MPI of TCs is a function of SST,outflow temperature and other relevant parameters,where the outflow temperature over tropical and subtropical oceans is strongly controlled by SST (Reid and Gage,1981).In particular,the outflow temperature was shown to be linearly correlated to the SST,especially when the SST is higher than about 24 C but smaller than 29 C (DeMaria and Kaplan,1994;Schade,2000) [see Fig.1 in Schade (2000)].Therefore,the MPI of TCs is almost exclusively determined by the SST.Moreover,the rate of TC intensification is also strongly affected by SST (?rnivec et al.,2016;Xu et al.,2016).Specifically,it is known that TCs absorb heat energy from the ocean and intensify due to the strong TC-ocean interaction which is potentially enhanced when TCs encounter warm oceanic eddies or rings (Lloyd and Vecchi,2011;Yablonsky and Ginis,2012;Kilic and Raible,2013;Ma et al.,2017).The rapid intensification of both Hurricane Opal and Katrina,occurred when they moved through warm oceanic eddies.In the case of Katrina,it was the warm Loop Current of the Gulf of Mexico.(Hong et al.,2000;Shay et al.,2000;Scharroo et al.,2005).More specifically,Hong et al.(2000) demonstrated that a 1K SST increase will induce a drop in minimum central pressure of a TC by 10 hPa through a warm eddy sensitivity experiment.On the other hand,due to entrainment/mixing and upwelling processes,SST cooling always occurs along the right side of the track and can be as large as 4 K (Price,1981;Schade and Emanuel,1999;Schade,2000;Davis et al.,2008).The cooler sea surface will inhibit the upward entropy flux and eventually reduce the intensity of TCs.Therefore,in order to improve the forecasting skill regarding TC intensity,a well-simulated SST forcing field is necessary.However,quite a lot of numerical models use a fixed SST forcing field and ignore the“SST cooling”,leading to an overestimation of the TC intensity (Winterbottom et al.,2012;Sun et al.,2014).Even if coupled models were used,the model errors of the atmospheric and oceanic components and its coupling frequency will reflect the error of a simulated SST(Davis et al.,2008;Strazzo et al.,2016;Scoccimarro et al.,2017).We may conclude that the SST error is inevitable.Therefore,it is necessary to estimate the potential effect of SST uncertainties upon TC intensity and then try to minimize this effect,ultimately improving the TC intensity forecast skill.

    Most of operational TC forecasting models are composed of numerical weather forecasting models with a fixed SST forcing field,rather than applying an ocean-atmospheric coupled model,although the latter is of great expectation.Therefore,we analyze these weather forecasting models to explore the effect of the uncertainties of SST forcing on TC intensity simulation errors.To make the SST forcing more realistic,we adopt,in the present study,the observed time-dependent SST,rather than a fixed SST,as an external forcing of TC system and consider the effect of the errors superimposed on these SSTs upon the simulation of TC intensity,in attempt to explore which error has the largest effect on TC intensity.Here,the error of the SST forcing may describe the uncertainties occurring in SST due to inaccurate SST observations or imperfect TC-ocean interactions as indicated by a coupled model.Based on this analysis,we naturally ask how to improve the accuracy of the SST forcing or the TC-ocean interaction factors associated with TC intensity.It is certainly true that an increase in SST observations will improve the SST forcing field;while improvements regarding the TC-ocean interaction factors require improvements within a coupled model,which is also dependent upon having sufficient observations.Obviously,both of these aspects rely upon increasing observations.The question then is,in order to improve the result of TC intensity simulation efficiently,in which regions should the density of SST observations be increased?

    The above question is related to a new observational strategy called“targeted observation”(Snyder,1996;Mu,2013).The tool of“targeted observation”was developed in the 1990s and originally proposed for an initial value-problem.Its general idea is as follows:to better predict an event at a future time

    t

    (verification time) in a focused area (verification area),additional observations are deployed at a future time

    t

    (targeted time,

    t

    t

    ) in some special areas (sensitive areas) where additional observations are expected to contribute most profoundly to reducing the prediction errors in the verification area [Mu (2013)].These additional observations can then be introduced into a model by a data assimilation system to form a more reliable initial state,which results in a more accurate prediction or simulation [see Fig.1 in Majumdar (2016) for a schematic example for targeted observation].

    Fig.1.The patterns of the NFSV-type SST forcing errors (K) of the selected 12 TC cases.The blue,green,yellow,red and purple dots indicate the TC intensity of the unperturbed run within [980,1000],[970,980],[960,970],[950,960]and [900,950](hPa).

    Generally,targeted observation is used to decrease the initial error (Peterson et al.,2006;Buizza et al.,2007;Wu et al.,2007;Qin and Mu,2012;Duan and Hu,2015;Zou et al.,2016;Zhang et al.,2017).To deal with uncertainties of external forcing upon simulation skill,Wen and Duan(2019) extended the idea of targeted observation to treat the forcing error and considered which observations are more helpful for reducing forcing error and improving simulation skill.In the present study,TC simulation is investigated from the perspective of the SST forcing influencing TC intensity.Therefore,the targeted observation associated with reducing the forcing error proposed by Wen and Duan (2019)can be reasonably adopted to deal with the external forcing of SST observations for TC intensity simulation.

    The key of targeted observation is to determine the sensitive area,defined here as the area where the simulation uncertainties are most sensitive to the forcing errors.In the present study,we would first identify the most sensitive SST forcing error and then subject the sensitive area to targeted observation associated with a TC intensity simulation.In order to identify the most sensitive error of the SST forcing,the approach of Nonlinear Forcing Singular Vector(NFSV) proposed by Duan and Zhou (2013) is useful.The NFSV represents the forcing error leading to the largest forecast/simulation error.The NFSV approach has been applied to the predictability studies of ENSO and Kuroshio Current effectively and succeeded in obtaining the most sensitive forcing error (Duan and Zhao,2015;Wen and Duan,2019).Through observation system simulation experiments(OSSEs),Wen and Duan (2019) showed that the region of large values in the NFSV-type errors represents the sensitive area for targeted observation associated with external forcing errors.In light of these successes,the present study uses the NFSV approach to determine the sensitive area for targeted observation associated with the TC intensity simulation.Thus,we seek answers to the following questions.1) What kind of SST forcing error can lead to the largest TC intensity errors? 2) What kind of properties does the sensitivity of TC intensity on SST error have? 3) Does the structure of NFSV-type errors indicate the sensitivity of TC intensity on SST errors? 4) Which region,in time and space,represents the sensitive area for target observation for TC intensity simulation?

    The arrangement of this paper is as follows:the settings of model,approach and algorithm is briefly introduced in section two.The NFSV-type SST forcing errors for 12 TC cases are calculated and corresponding sensitivity are shown in section three.In section four,the TC states responsible for the occurrence of the NFSV-type SST forcing errors are revealed and associated physical mechanisms are explored.In section five,the mechanism of the NFSV-type SST forcing error inhibiting TC intensity is discussed.Finally,a summary and a discussion are made in section six.

    2.Model,approach and algorithm

    2.1.Model

    The regional model used in the present study is the Advanced Research Weather Research and Forecasting(WRF) model (ARW) in its version 3.8.1 (Skamarock et al.,2008).This particular version of the WRF model is fully compressible and based upon the non-hydrostatic Euler equations and has often been used in studies of TCs.The model adopts the microphysics scheme of Lin et al.(1983) and the Kain-Fritsch scheme for cumulus parameterization (Kain,2004).The model considers the longwave and shortwave radiation by using the Rapid Radiative Transfer Model (RRTM;Mlawer et al.,1997) and the Dudhia scheme (Dudhia,1989).The boundary layer is parameterized by the Yonsei University scheme (Hong et al.,2006).The TC simulation conducted by the WRF is subject to SST forcing that is updated every 6 hours;that is to say,TC is forced by a time-dependent SST field.The WRF here adopts the horizontal resolution of 30 km×30 km without nesting and 24 levels in the vertical direction.The model top is set as 5000 Pa and the time step of simulation is 90 s.All TC cases in the present study are simulated for 120 hours.

    The atmospheric data (including variables associated with wind,pressure,cloud,soil,precipitation,etc.) adopted here is the FNL reanalysis from National Centers for Environmental Prediction (NCEP),whose resolution is 1×1 degree.The SST observation data is from NCEP and has resolution of 0.083×0.083 degree whose interpolated data is updated every six hours by the ungrib process,which assimilates in the SST data that is used to force the WRF model.

    2.2.Approach:Nonlinear Forcing Singular Vector

    The Nonlinear Forcing Singular Vector (NFSV) was proposed by Duan and Zhou (2013),which is an extension of (linear) forcing singular vector (FSV;Barkmeijer et al.,2003)in nonlinear regime.For convenience,the NFSV is briefly described as follows.

    Assume that the Eq.(1) describes a state equation

    where

    G

    (

    U

    ) and

    F

    (

    x

    ,

    t

    ) is model equation tendency and external forcing,and

    f

    (

    x

    ) is a forcing error of the forcing term.According to the definition of NFSV,it is the tendency perturbation that generates the largest deviation from the reference state in a nonlinear model at a given time based on a physical constraint condition.For a given forcing error,the NFSV can be understood as the forcing error that has the largest effect upon the simulation or prediction error at the given future time.This can be formulated into the following optimization Eq.(2):

    where

    J

    is cost function and the norms ||·||and ||·||measure the amplitude of the forcing error

    f

    and its resultant simulation error against the reference state,respectively.The forcing error,

    f

    ,is subject to the constraint radius δ;

    M

    (

    f

    ) and

    M

    (0) are the propagators of a nonlinear model with and without forcing error

    f

    (

    x

    ) from time 0 to

    t

    ,respectively;and

    U

    is the initial value of the reference state.By solving Eq.(2),the NFSV,denoted by

    f

    *in Eq.(2) can be obtained.

    2.3.Algorithm:Particle swarm optimization

    The particle swarm optimal (PSO) algorithm was initially proposed by Kennedy and Eberhart (1995) to imitate the process of bird foraging,but soon it was widely used to solve optimization problems and,in doing so,achieved great successes (Banks et al.,2008;Zheng et al.,2017).We will use this algorithm to calculate the optimization problem associated with NFSV.Next,we briefly describe the algorithm.

    To make the cost function as in Eq.(2) to have the largest value in the given constraint condition,a series of particles characterized by positions (denoted by

    X

    ;the forcing perturbations here) and velocities (denoted by

    V

    ;the iteration velocity) are randomly generated.Then the cost function is calculated with these particles.These particles will be updated by iterations according to the values of the cost function.Specifically,the iterations are realized by calculating the Eqs.(3) and (4).

    For a high-dimensional dynamic system,it is impossible to generate particles with the same number as the dimensions of the model dynamical system to calculate the NSFV by using the PSO algorithm.Thus,we apply an empirical orthogonal function (EOF) analysis to reduce the dimensions and obtain the representative particles required by the PSO.For the SST forcing errors associated with TC intensity,we adopt the following strategy to generate the particles of the PSO algorithm.

    (i) We take the SST field in the north-west Pacific region (i.e.10o?35oN,100o?150oE) every five days from 1 July to 30 September or from 1 June to 31 August during 2006?16.Which period is selected depends upon the season when TCs happen most frequently during the year.Then,we subtract the SST field on 1 July from that on 5 July (or on 1 June from that on 5 June).This method is also applied to the SST fields on 5 July from that on 10 July (or on 5 June from that on 10 June),on 10 July from that on 15 July (or on 10 June from that on 15 June),and so on.Then 198 SST forcing perturbations are obtained.

    (ii) An EOF analysis is applied to the 198 SST forcing perturbations and the leading 30 dominant modes are experimentally selected to yield the NFSV.The 30 dominant modes explain 80% of total variance of SST forcing perturbations which is believed to be enough for searching the NFSV using the PSO algorithm.

    (iii) We scale the leading 30 modes to have the same amplitude in terms of the adopted norm (see next section)and assign them as the initial positions of particles,while the initial velocities of particles are,at first,guessed as being equal to their initial positions.

    3.NFSV structure and its sensitivity

    In the present study,12 TC cases are selected according to the best-track data from the Japan Meteorology Agency (JMA) for exploring the sensitivity of TC intensities on the SST forcing errors.The 12 TC cases are chosen according to the following criteria:(1) the TCs which greatly influenced China,(2) the lifetimes of TCs were longer than 5 days and stayed over the ocean during the overwhelming majority of the simulation period,and (3) the TC whose tracks were minimally altered when the SST forcing was modified,with the largest deviation of this subset being less than 60 km when the SST forcing is modified,which almost excludes the impacts of the bias of TC tracks on TC intensities.The brief information concerning the selected 12 TCs is listed in Table 1.The geometric centers of the TCs are identified as the locations with the minimum sea level pressures (MSLPs),which are utilized to represent the TC intensities hereafter.The NFSV represents a special type of SST forcing error,with which the simulated TC intensities depart the most from,as compared to without it,as Eq (5)describes.

    Table 1.Twelve TC cases of investigation.

    The NFSV-type SST forcing errors are calculated for the predetermined 12 TC cases and plotted in Fig.1.It is shown that,although the tracks of these 12 TC cases differ a lot from each other,the NFSV-type SST forcing errors are always along the TC tracks.This indicates that the SST errors in the areas along the TC tracks,compared with those in the areas away from the TCs,are likely to exert a greater influence the TC intensities.Although the NFSV-type SST forcing errors are located along the TC tracks,they display the largest anomalies in different time periods (e.g.from 24 h to 48 h for Soulik,from 60 h to 96 h for Rammasun,and so on;see Fig.1) of the different TCs.Recalling the definition of the aforementioned NFSV scheme,the NFSV-type SST forcing errors represent the forcing errors that result in the largest TC intensity simulation errors over the 120 hrs.Then the distribution of the NFSV-type SST forcing error along the TC track may indicate that the TC intensity simulations that are significantly sensitive to the SST forcing errors occurring in a particular time period of TC movement.

    The regionally-averaged NFSV-type SST errors of the selected 12 TCs,within a radius of 300 km,centered at the central location of the TC,and at

    t

    are respectively calculated,which is plotted in Fig.3 as a function of the spread of the 22 perturbed runs in TC intensity for each TC during the interval,[

    t

    ,

    t

    +6].With the change of

    t

    ,the linear regression between the regionally-averaged NFSV-type SST forcing errors and the ensemble spread is calculated (see the line in Fig.3).It follows that the regionally-averaged NFSV-type SST forcing errors are significantly correlated with the spread of the intensity in the perturbed runs of the TCs,which demonstrates significance at the 0.01 level when subjected to a

    t

    -test.That is to say,the larger the ensemble spread during one time period of TC,the larger the corresponding NFSV-type SST forcing errors.It is therefore obvious that the NFSV-type SST forcing errors can identify the time period when the TC intensity simulations are highly sensitive to the SST forcing errors.

    Fig.2.The spread of TC intensity simulations perturbed by 22 SST forcing perturbations during [t0,t0+6]with t0being 0 h,6 h,12 h,...,114 h and the TC intensity of the unperturbed run.Red lines denote the TC intensity [indicated by minimum sea level pressure (MSLP);units:hPa]of the unperturbed run.The blue bars represent the spread of the 22 perturbed runs of the TCs during [t0,t0+6](units:hPa).

    Fig.3.The regionally-averaged NFSV-type SST forcing errors at t0as a function of the spreads of TC intensity simulations perturbed by the 22 SST forcing perturbations during [t0,t0+6]with t0being 0 h,6 h,12 h,...,114 h.The red line represents the linear regression line between NFSV-type SST forcing errors and spreads.

    The 22 randomly-selected SST forcing perturbations described above are further-superimposed on the unperturbed SST fields to force the TCs for the whole simulation period of 120 hrs.Then the total error of the TC intensity in each perturbed run,during the 120 hrs,is calculated for each TC.After this is done,the correlation coefficients are calculated at each grid point,for each TC between the total errors of the TC intensities in the perturbed runs and the corresponding member of the 22 SST forcing perturbations(see Fig.4).Figure 1 shows that the distributions of the correlation coefficients are very similar to those of the corresponding NFSV-type SST errors.More specifically,the spatial correlation coefficients between them can be as large as 0.64,on average,for the 12 TC cases (the details can be seen in Table 2),indicating that the larger the ensemble spread particular to a TC location,the larger the NFSV-type errors there.It can also be seen that the correlation coefficients are much larger along the TC track,which implies that the total error of the TC intensity is especially sensitive to the SST forcing errors along the TC track.The similarity between the distributions of the correlation coefficients and the corresponding NFSV-type SST errors indicates that the NFSV-type SST errors link the sensitivity of the TC intensity simulation errors to the SST forcing errors in space.Furthermore,this relationship confirms that the SST errors in the areas along the TC tracks,especially those during the time period when the NFSV-type SST forcing errors attain large values,may significantly influence the TC intensities.

    4.TC states responsible for the occurrence of the NFSV-type SST forcing errors

    It is clear now that the NFSV-type SST forcing errors can identify the particular time period of the TC movement when the TC intensities exhibit high sensitivity to the SST forcing errors;but such time periods are dependent on the individual TC cases.The issue then becomes,whether or not these different time periods for different TCs correspond to common physical and environmental states of TCs.That is to say,what physical and environmental factors determine the sensitivity displayed by the NFSV-type SST forcing errors? Since the NFSV-type SST forcing errors cause the largest simulation errors of the TC intensity,we choose to explore the contributing factors of the NFSV-type SST errors by addressing which states of TCs are favorable for the SST forcing errors that yield large TC intensity errors.We go on to explain the physical processes responsible for the TC state that is consistent with the formation of the NFSV-type SST forcing errors.

    4.1.Which states of the TCs are favorable for the SST forcing error causing large intensity error?

    We classify the lifetimes of TC movement into two categories of time periods according to the sensitivities of the simulated TC intensities in perturbed runs to the SST forcing perturbations.Specifically,for the selected 22 SST forcing perturbations and the time periods [

    t

    ,

    t

    +6]as in section 3,if a time period possesses a spread (among the 22 perturbed runs in TC intensities) greater than 1.5 hPa (i.e.the mean value of the ensemble spreads during [

    t

    ,

    t

    +6]with the changing

    t

    ),this time period is categorized as a relatively high sensitivity period (referred to as“H-Sen”hereafter);conversely,the other time periods are categorized as relatively low sensitivity periods (referred to as“L-Sen”hereafter).For the environmental factors,the SST,relative humidity (RH),vertical wind shear (VWS),and translation speed are considered.During the period [

    t

    ,

    t

    +6](noting that

    t

    changes),all of the above environmental factors are calculated at

    t

    ,

    t

    +3 hrs,and

    t

    +6 hrs for the unperturbed run,respectively.At each of these timings,the SST is regionally-averaged in a circular domain with a radius of 300 km centered at the simulated TC center;it is then further averaged consistent with the three timings.This averaged SST represents the environmental SST of the TCs in the period [

    t

    ,

    t

    +6].The VWS is similarly calculated,but denotes the difference of the regionally-averaged horizontal wind between 200 hPa and 850 hPa;and the RH is vertically averaged in the layers between 1 km and 6 km in the aforementioned circular domain.The translation speed,is calculated by taking the difference of the centers of the TC at both

    t

    and

    t

    +6 hrs and then dividing this distance by the 6 hrs time interval.All of the above calculations are classified according to H-Sen and L-Sen and the results are shown in Table 3.It is found that the TCs in the H-Sen periods,compared with those in the L-Sen periods,have much more humid innercores,move profoundly slower over warmer SSTs,and are accompanied by stronger VWS.In previous studies,all of these environmental factors were shown to benefit the intensification of TCs (Mei et al.,2012;Zhang and Tao,2013;Tao and Zhang,2014;Walker et al.,2014;Torn,2016;Zhao and Chan,2017).Nevertheless,when we calculate the correlation coefficients between the spread of the 22 perturbed runs in TC intensity during the period [

    t

    ,

    t

    +6],with the changing

    t

    ,and the environmental factors of SST,RH,VWS,and translation speed in the unperturbed run,they are shown to be very low and thus weakly correlated [see Fig.5 (a?d)].This suggests that the environmental factors of TCs cannot be responsible for the relatively high sensitivity in the HSen period of TCs.

    Fig.4.The spatial correlation coefficients (shaded) between the 22 SST forcing perturbations and the sum of the absolute values of intensity errors during the whole simulation period of 120 hrs for the selected 12 TC cases.The blue,green,yellow,red and purple dots represent the simulated TC intensity within [980,1000],[970,980],[960,970],[950,960]and[900,950](hPa).

    The characteristics of the TCs themselves in the H-Sen and L-Sen periods are also examined and associated inflow,vertical velocity,and inertial stability influencing TC intensity are calculated,respectively.Specifically,the inflow is calculated by taking the regionally-averaged radial component of the horizontal wind in a circular area centered at the simulated TC center,with the radii between 50 km and 300 km and vertically-averaged from 0 km to 1 km;the vertical velocity is estimated by a similar scheme but estimates the vertical wind by taking the vertically-and regionally-averaged component from 0 km to 15 km within a circular area centered at the simulated TC center with the radii between 50 km and 150 km;and the inertial stability is calculated by the formula I=(f+ξ)(f+2v/r) (where

    f

    is the Coriolis parameter,v represents the tangential wind velocity,ξ denotes the relative vorticity,and r is the radius),which is then vertically and regionally-averaged from 0 km to 15 km and in a circular area,centered at the simulated TC center,with a radius of 100 km.The above results are also listed in Table 3 according to H-Sen and L-Sen.It is shown that the TC intensity in the H-Sen period is about 957.9 hPa on average,which is stronger than the average of 965.9 hPa in the L-Sen,and is found to be significant at the 0.0016 level using a t-test.This indicates that the TC intensities show higher sensitivity to the SST forcing errors when they are much stronger.We also note similar results from Table 3 with respect to the inflow,vertical velocity,and inertial stability;that is,these variables are much larger in H-Sen period,which is consistent with the presence of relatively strong TCs in the H-Sen periods.Furthermore,we calculate the correlation coefficients between the spread of the 22 perturbed runs in TC intensity during the period [

    t

    ,

    t

    +6]with the changing

    t

    and corresponding inflow,vertical velocity,inertial stability,and MSLP of the unperturbed run.We find that they,compared to those for the environmental factors,are more significantly correlated.This indicates that the TC intensities are much more sensitive to the SST forcing perturbations when the TCs exhibit strong intensity,larger inflow,large vertical velocity,and strong inertial stability (see Fig.5).Moreover,from Fig.2,it is shown that the TC cases tend to have the largest forecast spread when the TCs exhibit relatively strong intensity but are still intensifying,just prior to being fully mature,with the notable exception of TC Noul.To facilitate the discussion,such spatio-temporal period is hereafter called the Critical Time and Phase (CTP).

    Table 2.The spatial correlation coefficient between intensity errors and the 198 SST perturbations for 12 TC cases.

    Table 3.TC states and their environmental factors during H-Sen and L-Sen.

    When we further examine the corresponding secondary circulation and inertial stability,it is found that the CTP period shows strong secondary circulation but the highest inertial stability for TCs.As an example,Fig.6 plots the evolution of inertial stability,inflow,and vertical velocity of the TC Soulik (201307).It is shown that the inertial stability and the vertical velocity reach the largest values during the time period from 24 h to 48 h,while the inflow is still increasing in this period.Here,the time period from 24 h to 48 h fits the CTP of TC Soulik (also see Fig.2).Therefore,we conclude that the strongest sensitivity of TC intensities to SST forcing perturbations occurs during the CTP of the TCs,especially with a TC state of the highest inertial stability and the strong secondary circulation.Such time periods can also be seen in Fig.1,i.e.the period with the deepest color in the shaded area,during which the SST errors associated with the NFSV are of the largest assigned values.Particularly,the time period with the deepest color in the shaded area for the TC Soulik is about from 24 h to 48 h,which coincides with the CTP of the TC Soulik Therefore,the NFSV-type SST forcing errors locating at the TC track not only describe the sensitivity of the TC intensities to the SST forcing errors (see section 3) but also capture the time period of the TC movement when the TC intensity presents the strongest sensitivity to the SST forcing errors.Based upon these results,we advance the hypothesis that the TC states of high inertial stability and the strong secondary circulation are responsible for the NFSV-type SST forcing errors being dominated by the errors occurring during the CTP.

    4.2.Why do the SST forcing errors in the H-Sen period influence the TC intensity more significantly?

    Fig.5.The scatter plots between spreads of TC intensity simulation yielded by the 22 SST forcing perturbations and corresponding (a) SST forcing (units:K),(b) RH (units:%),(c) translation speed (units:km h?1),(d) VWS (units:m s?1),(e) vertical velocity (units:m s?1),(f) inflow (units:m s?1),(g)inertial stability (units:10?3 s?2) and (h) MSLP (units:hPa).The red lines represent the linear regression.The red dots are for the H-Sen period and the black dots are for the L-Sen periods.The SST,RH,translation speed,VWS,vertical velocity,inflow,inertial stability and MSLP are calculated as in Table 3 (see section 4).

    The eleven TC cases with the strongest sensitivity in their respective CTP are used to address the nature of the physical mechanisms which explain the enhanced sensitivity.To facilitate the discussion,we take TC Soulik as an example.It has been previously mentioned that the distribution of the NFSV-type SST forcing errors reflects the sensitivity of TC intensities to the SST forcing errors;and the NFSV-type SST forcing errors tend to have the largest anomalies during the time period from 24 h to 48 h of TC movement when the TC is during the CTP (see the NFSV-type SST forcing error pattern for the TC Soulik in Fig.1).This time period is coincident with the strongest sensitivity of TC intensity on the SST forcing errors and is particularly associated with a strong secondary circulation and very highest inertial stability of the unperturbed run.The issue of concern is whether or not the along-track SST forcing errors that occurred during this unique and relative stage and timing in TC evolution,significantly differ from the response of TCs that exist at another timings in the TC growth cycle.

    The errors of the TC intensities in the experiments of UN-all,UN-Sen,and UN-non-Sen with respect to the unperturbed run are shown in Fig.7a.It can be seen that the error of the TC intensity in the UN-Sen experiment is significantly larger than that in the UN-non-Sen experiments and accounts for a forecast error of almost 80% compared to that in the UN-all experiment (Fig.7b).This indicates that the total error of TC intensity in the UN-all experiment is mainly caused by the SST forcing errors during the H-Sen period.In addition,we can notice from Fig.7a that the growth rate of the TC intensity error in the UN-Sen experiment is much larger than that in the UN-non-Sen experiment.Specifically,during the time periods from 0 h to 24 h and from 96 h to 120 h of the L-Sen periods in the UN-non-Sen experiment,the intensity errors are constrained to be than 6 hPa;while in the UN-Sen experiment,the intensity errors increase by about 10 hPa during the time period from 24 h to 48 h of the H-Sen period,which is the equivalent time interval of 24 hrs as the former two periods(from 0 h to 24 h and from 96 h to 120 h,in the L-Sen period).Furthermore,the time period from 24 h to 48 h of the H-Sen period coincides with the interval when the NFSV-type SST forcing error of the TC Soulik achieves the largest values and the TC intensities are most sensitive to the SST forcing perturbations.The rapid increase of the intensity error during the time period from 24 h to 48 h can also be seen in the UN-all experiment if one uses the slope of the error evolutionary curve to measure the error growth (see Fig.7a).In particular,the time period from 24 h to 48 h shows the fastest growth of intensity error for the TC Soulik in the UN-all experiment,therefore,contributing the most to the total error of the TC intensity during the entire simulation period.This may explain why the NFSV-type SST forcing errors occur in the H-Sen period,particularly during the time period from 24 h to 48 h in the case of TC Soulik.

    Fig.6.The evolution of (a) inertial stability (units:10?3 s?2),(b) inflow (units:m s?1) and (c) vertical velocity (units:m s?1)for the TC Soulik.The inertial stability is regionally-averaged in a round area centered at the simulated TC center within 50 km;the inflow is an azimuthal mean of radial velocity,which is vertically-averaged in the layers between 0 km and 1.5 km;and the vertical velocity is vertically-averaged in the layers between 3 km to 7.5 km.The black lines in (b) and (c) denote the RMW at 2.0 km.

    Fig.7.(a) The evolutionary curves of TC intensity of unperturbed run (black line) for the TC Soulik and its error evolutionary curve of TC intensity perturbed by the NFSV-type SST forcing error (orange);and the error evolutionary curves of TC intensity in the UN-all (blue),UN-Sen (purple),and UN-non-Sen (red) experiments;(b)the total error of the intensity [see Eq.(5)]for the TC simulation perturbed by the NFSV-type SST forcing error(orange),and the TC simulation in the UN-all (blue),UN-Sen (purple),and UN-non-Sen (red) experiments.Note that the maximum intensity error measured by MSLP forced by NFSV-type SST forcing errors can reach up to about 35 hPa during 72?96 h,which is comparable with the 23 hPa of the annual mean of root mean square errors (RMSEs)for TC center pressure forecasts with lead time 72 hrs.

    Next,we will explore how the SST forcing errors in the H-and L-Sen periods influences TC intensity by perturbing the processes influencing TC intensity.To facilitate the calculation,we select only the 18 h,114 h,and 42 h of the unperturbed run as representative of two L-Sen and one H-Sen periods to calculate the TC processes influencing TC intensity and associated simulation errors,where these three timings are all relative to the initial time of their respective time periods for 18 hrs.

    Fig.8.TC Soulik:azimuthal mean of surface latent heat flux errors (black lines,units:W m?2),water vapor errors (shaded,units:g kg?1),and diabatic heating errors (red contours,units:K h?1,contour interval:0.5 K h?1) are plotted at (a) 18 h,(b)42 h,and (c) 114 h;azimuthal mean of potential temperature errors (shaded,units:K),vertical velocity errors (red contours,units:m s?1,contour interval:0.1 m s?1),and inflow errors (blue contours,units:m s?1,contour interval:1 m s?1) are shown at (d) 18 h,(e) 42 h,(f) 114 h.

    4.3.The effect of inertial stability

    The results in Table 3 and Fig.5 show that the greatest sensitivities of TC intensity on SST forcing errors are associated with conditions of high inertial stability.We then pose the following two questions.How does the inertial stability affect TC intensities and through what mechanism does the high inertial stability in H-Sen period cause the SST forcing errors to yield such large intensity errors? Here,we show that the high inertial stability in the H-Sen period is favorable for the large growth of errors in TC intensity that is caused by the SST forcing errors.Furthermore,we suggest that this is caused by the response of the secondary circulation to the heat forcing.In the developments to follow,we address this issue by analyzing the Sawyer-Eliassen (SE)equation (Montgomery et al.,2006;Chen et al.,2018)through sensitivity experiments of the unperturbed run associated with TC Soulik with respect to the TC state in the UNSen experiment.

    Fig.9.The azimuthal mean of (a) vertical velocity (units:0.06 m s?1) and (c) inflow (units:m s?1) in the SE equation for the unperturbed run TC Soulik at time 42 h;the differences of the SE solution in (b) vertical velocity (units:0.006 m s?1) between the run in (a) and that with the coefficient C in the UN-Sen experiment and in (d) the inflow(0.1 m s?1) between the run in (a) and that with the coefficient C in the UN-Sen experiment.

    5.The mechanism of NFSV-type SST errors affecting TC intensity

    In section 4,we showed that the NFSV-type errors occur during the time period (i.e.the H-Sen period above)when the TCs present strong intensities,strong secondary circulations and high inertial stability.Moreover,the NFSVtype SST forcing errors tend to have the largest anomalies when the TCs are in the CTP,which identifies the time period when the intensity of TCs exhibit the strongest sensitivity to the SST forcing errors.In sections 4.2 and 4.3,we analyze TC Soulik as an example to explain why the aforementioned TC states are favorable for the growth of the intensity errors which are caused by the NFSV-type SST forcing error and we furthermore advance a physical mechanism which logically explains the the response of the secondary circulation to the NFSV-type SST forcing error.All evidence leads to the conclusion that the strongest sensitivity of TC intensity to the SST forcing errors along the TC track during the period when the TCs are of the CTP state.Therefore,if we manage to implement additional along-track SST observations during this particular time period,it would help to obtain a much more accurate SST forcing field for the WRF model,ultimately,reducing the simulation uncertainties of the TC intensity.If one uses a coupled model to simulate the TC,the along-track SST in this time period should be better simulated so as to greatly improve the TC intensity simulation skill.

    In the present section,we continue to use the TC Soulik to explore how the NFSV-type SST forcing errors perturb the TC intensities by influencing the processes associated with TC intensification.From Fig.7a,it is shown that,from about 24 h to 48 h,the intensity error caused by the NFSV-type SST forcing error produces larger growth rates(measured by the slope of the error evolutionary curve) than that in the UN-all experiment,which then causes larger intensity errors during the mature phase (i.e.from 48 h to 96 h) of TC Soulik despite the SST forcing errors in the UN-all experiment -having the same amplitude as in the NFSV-type SST forcing errors.We have known that the time period from 24 h to 48 h corresponds to the one when the NFSV-type SST forcing error possesses the largest anomalies and that the TC intensities are most sensitive to the SST forcing errors,which,as revealed in section 4,can explain why the NFSVtype SST forcing errors cause much larger intensity errors of a TC.We now pose the question,how does the NFSVtype SST forcing error influence the processes associated with the TC intensity and finally perturb the TC intensity?

    Since we use the minimum MSLP to measure the intensity of TC,it is required to figure out the mechanism of the NFSV-type SST forcing error resulting in the change of the MSLP of the TC.As mentioned above,the MSLP of the TC is strongly correlated with the potential temperature in the upper layers of the TC.Therefore,we derive the potential temperature (PT) error tendency equation by examining the difference between the PT tendency equation component of the WRF model associated with unperturbed run and perturbed run mentioned above.The equation is given in Eq.(12):

    Here,the overbar denotes the azimuthal mean,the star signifies the deviation from the azimuthal mean,the prime indicates the error caused by SST forcing error,and

    R

    is the residual term associated with unresolved processes as in section 4.However,for the residual term

    R

    ,we cannot exactly separate the role of each of its inclusive processes.That is to say,the total effect of

    R

    is not of clear physics.For simplicity,we only considered the role of the terms with clear physics and do not analyze

    R

    here.The meanings of the other terms are listed in Table 4.

    With the NFSV-type SST forcing error disturbed,we calculate the terms of Eq.(12) (confined in the eye region of TC within a round area centered at the simulated TC center with a radius of 50 km) for the TC Soulik.We find that the terms II,V and VII are much larger,which indicates that the diabatic heating error,eddy process error,and vertical advection of the PT by the vertical velocity error play an important role in generating PT error.These three terms are plotted in Fig.10.It is found that the term VII (i.e.the vertical advection of the azimuthal mean of PT by the azimuthal mean of vertical velocity error) is negative especially in the upper layers (i.e.about 9 km above the surface for TC Soulik).Out of the three terms,Term VII is the largest when measured according to absolute value.This indicates that the term VII plays a dominate role in affecting PT,further noting that the vertical advection denoted by the term VII contributes to suppressing PT growth,especially when the TC is still intensifying just prior to the mature phase [i.e.from 30 to 48 h,which is within the window of the most sensitive period from 24 h to 48 h identified by the NFSV-type SST error].Since the positive direction of the vertical PT gradient vector is upward,the negative contribution of the term VII to the PT growth mainly results from the positive error of vertical velocity (see Fig.10f).As we know,in the eye region of TC,the vertical velocity is generally downward.Thus,the positive vertical velocity error reduces downward advection,and limits the amount of high entropy air parcels that enter the warm core which results in a decrease of TC intensity,ultimately yielding a negative error of TC intensity.Compared to term VII,term V,the error of eddy component of radial PT advection is much smaller,which indicates the SST forcing error has a negligible effect upon the asymmetric structure of the vortex associated with TC.Term 2,the diabatic heating error,is also relatively small,which may be due to the latent heat release occurring in the convection region instead of the eye region.In fact,as shown in section 4,the important effect of diabatic heating error,in the convection region,upon TC intensity simulation is uncertain.Therefore,the change of diabatic heating induced by the NSFV-type SST forcing error affects the warm core of TC in an indirect way.

    We plot in Fig.11a the evolution of the absolute value of negative latent heat flux errors caused by the NFSV-type SST forcing error for TC Soulik.It is obvious that the latent heat flux error increases from 0 h to 48 h,especially within the interval of 24 h to 48 h due to the increase in surface wind consistent with the strong intensity of the TC during this period.Correspondingly,the water vapor error shows similar evolutionary behavior,with the largest reduction during the same sensitive period,from 24 h to 48 h.Then the water vapor advected into the inner core and upper layers is reduced as a consequence of the error,which causes large negative diabatic heating errors in the upper layers especially during the mature phase (i.e.from 48 h to 96 h) of TC Soulik,where the strong secondary circulation of the unperturbed run,as clarified in section 4,enhances the diabatic heating error in the upper layers induced by the water vapor error.Fig.11b shows negative diabatic heating errors and vertical velocity errors forced by diabatic heating.It is shown that both the vertical velocity and the secondary circulation weaken,due to the effect of negative diabatic heating error.Despite the fact that the maximum of the negative diabatic heating errors do not occur during the most sensitive period of TC Soulik (see Fig.11b),the strongest inertial stability of the unperturbed run during this period (see Fig.6a) triggers a response of the secondary circulation to the diabatic heating and significantly increases the errors occurring in the secondary circulation,ultimately decreasing the vertical velocity,to the largest extent,during the most sensitive period from 24 h to 48 h.Particularly,we can see from Fig.6a that the inertial stability is highest during the period from 24 h to 48 h and correspondingly,the vertical velocity error is the largest during this period.In response,the subsidence in the inner core is increased by a positive error (see Fig.11c),especially from 24 h to 48 h.This results in reducing potential temperature of the upper layers of the inner core which leads to negative errors in potential temperature which accumulate rapidly from 24 h to 48 h(see Fig.11c).The potential temperature then exhibits oscillatory behavior during the mature phase (i.e.from 48 h to 96 h) before dropping quickly from 96 h to 120 h,which coincides with the evolutionary behavior of the TC intensity error caused by the NFSV-type SST forcing error shown in Fig.7a.Thus,the above mechanism interprets how the NFSV-type SST forcing errors perturb the intensity error.

    Table 4.The meaning of the terms in the PT error tendency equation.

    Fig.10.The time-height cross section of terms in potential temperature (PT) error tendency equation.PT error is caused by NFSV-type SST errors.All the terms are calculated by regional average within a round area centered at the simulated TC center with a radius of 50 km.(a) the term I;(b) the sum of rhs of Eq.(12) except the residual term;(c) the term II;(d) the term V;(e) the term VII;(f) the vertical velocity error (units:0.1 m s?1) caused by NFSV-type SST errors.

    Fig.11.Evolution of (a) surface latent heating errors(black line;units:W m?2) and water vapor errors (shaded;units:g kg?1) regionally-averaged in a round centered at the simulated TC center within a radius of 300 km,(b)vertical velocity errors (red contours;units:m s?1,contour interval:0.1 m s?1) and diabatic heating errors (shaded;units:10 K h?1) regionally-averaged in a ring area centered at the simulated TC center with the radii between 50 km and 150 km,and (c) subsidence errors (shaded;units:0.1 m s?1) and potential temperature errors (blue contours;units:K,contour interval:2 K) regionally-averaged in a round centered at the simulated TC center within a radius of 50 km,which are all yielded by NFSV-type SST forcing errors.

    Apart from the TC Soulik,we also explore the other ten TC cases with the strongest sensitivity occuring in the CTP of TCs.They have mechanisms similar to TC Soulik in the NFSV-type SST forcing errors affecting the TC intensity,except that some cases show signs opposite to those of the TC Soulik in the NFSV-type SST forcing errors and their resultant TC intensity uncertainties.Therefore,we can summarize the physical mechanisms as follows.When a NFSV-type SST forcing error with negative anomalies (or positive anomalies,depending on TC cases) occurs,it causes the surface latent heat flux and water vapor in low layers to decrease(increase).This further leads to a reduction (an increase) in diabatic heating which forces a weaker (stronger) secondary circulation and causes the downward vertical velocity in the eye region to decrease (increase).Eventually,the warm core of the TC becomes weaker (stronger) and the TC intensity tends to be under-(over-) estimated.In this process,the unperturbed run of the TC tends to present relatively strong intensity with the highest inertial stability and strong secondary circulation during the time period when the NFSV-type SST forcing error occurs.This property of the unperturbed run greatly enhances the latent heat flux errors in the low layers due to the effect of SST forcing error which goes on to increase the diabatic heating error,thereby,promoting a significant response from the secondary circulation before finally increasing the TC intensity error,most notably during the CTP.

    6.Summary and discussion

    The present study explores which TC features are affected the most by the SST forcing error that leads to TC intensity simulation uncertainties,with the intent of providing useful insights and ideas concerning target observation for SST forcing associated with TC intensity simulation.The approach of applyinga nonlinear forcing singular vector (NFSV) upon 12 TC cases for their 120-hour simulations are used to address this concern.The results show that the SST forcing errors of the NFSV structure often cause the largest simulation error of TC intensity.The NFSV-type SST forcing errors tend to be distributed along the TC track but are mainly concentrated during the time period when the TC is of strong intensity with strong secondary circulation and inertial stability,and exhibit positive or negative SST errors (which are dependent on the particular TC cases).Specifically,the NFSV-type SST forcing error of the TC tends to be dominated by the errors during the time period when the TC is in the CTP.Analysis of both the spatial and temporal distribution of those domains exhibiting model sensitivity (high correlation between SST forcing errors and TC intensity simulation errors) reveals that the critical time period of the TC response is consistent with the CTP of TC.This is identified by the NFSV-type SST forcing errors,which,by design,recognizes the time when the TC intensity is most sensitive to the SST forcing error.Therefore,if one manages to deploy additional SST observations during that particular time period of TC movement and assimilate them to the model SST forcing field,the TC intensity simulation level will be greatly improved.Even if a coupled model is used,the particular time period for targeting SST observations identified here will enable us to know when the SST should be preferentially well-simulated for improving the TC intensity simulation.

    By sensitivity experiments and analysis of the SE equation,we show that the high inertial stability of the TC in the CTP determines the degree of the response of the secondary circulation to the SST forcing errors,which finally resolves to what extent the SST forcing errors influence the TC intensity.This indicates that if the SST forcing errors happen to occur during the CTP,the TC intensity will respond through the amplification or muting of the secondary circulation,thereby exerting a control upon subsidence in the core.The high inertial stability present during the CTP of the lifecycle of the TC greatly enhances this process and increases uncertainty of the TC intensity simulation at this time.Clearly,this argument explains the formation of the NFSVtype SST forcing errors for the TC intensity simulation.

    By tracing the evolution of the TC intensity error,we suggest a mechanism that the NFSV-type SST forcing errors influence TC intensity.In particular,when the NFSV-type SST forcing error occurs along track and during the aforementioned CTP,the following events are likely to occur if the SST forcing error is negative (positive):it will cause both the low-layer atmospheric temperature and water vapor to be lower (higher) and consequently,the diabatic heating to decrease (increase).At this time,the secondary circulation,dominantly forced by the diabatic heating,becomes weaker(stronger) which causes the downward vertical velocity in the eye region of TC to become smaller (larger),finally causing the warm core of the TC to weaken (strengthen).Since the minimum MSLP is a proxy for TC intensity,coupled with the fact that this is strongly tied to the strength of the warm core,it follows that the TC intensity is under-(over-)estimated which presents large uncertainty.Additionally,the very large inertial stability of the TC the sensitive time period promotes a vigorous response of the secondary circulation to the SST forcing errors and causes a rapid accumulation of the TC intensity errors during this period,ultimately contributing the most to the total error of the TC intensity.This mechanism further illustrates that the NFSV-type SST forcing errors determine the time period and the region(where and when) targeted additional observations of SST forcing are needed to promote a more accurate TC intensity simulation.

    Due to the destructive effects of TCs and associated danger,their direct observations,especially those of the ocean component,are often difficult to obtain and therefore very valuable.Even if some of the observations are retained,they are subject to large uncertainties.We have therefore researched,which physical variable,and within which region and time period,would be optimally suited for increasing the spatial density of observations with the intent of improving TC simulations and subsequently deepening the understanding of the TC system.We conclude that SST observations should be enhanced along the TC track,while the TC is strong,but still intensifying,just prior to the mature phase.This particular time period and region close to the eye of the TC provides a great logistical challenge to instrumentation personnel.The issue now becomes one of engineering and deployment.Can the observing equipment and instruments reach the sensitive area during stormy conditions and how will the instruments be deployed? The intent of this paper is restricted to isolating the spatio-temporal domain where the observations are needed and does not address the logistical issues related to instrument deployment,which is,at first glance,difficult and beyond the scope of the present study.

    The present study used a version of the WRF model that has a horizontal resolution of 30 km×30 km.This coarse resolution,in of itself,can cause the simulation of TC intensity and structure to exhibit large biases.For example,the uncertainties of the RMW (radius of maximum wind) and the radius of the gale force wind cannot be adequately simulated after the SST forcing errors are superimposed.It then becomes apparent that the NFSV and its resultant sensitive area/ time period for target observations are based on biased TC simulations.Therefore,the evidence presented here discloses which time period and in which region the observations of the SST forcing are particularly important for improving the accuracy of the TC intensity simulation that contains biases.Despite this limitation,we put forward that the sensitive period for TC intensity simulation demonstrated here is still instructive,despite the model resolution that leaves much to be desired.Moreover,the primary results from our study are evidenced from 11 TC cases,all of which have the aforementioned sensitive time period/area in common.The reason why the TC Noul does not show coherent result is still unclear,which needs to be explored in future.

    It is known that versions of the WRF model which use finer resolutions require more computational costs for the simulation of TC intensity.Furthermore,the computation of the NFSV is also expensive.When these two factors are combined,the TC simulations based upon a high-resolution WRF model with the NFSV computations,provide a challenge to computational resources.However,considering that the NFSV approach is useful for constructing theories for improving forecasting skill regarding TC intensity,we deem the computational investment worthy.We remain hopeful that the WRF model with finer resolution,together with a much more efficient algorithm for the NFSV approach will be used to accurately identify the sensitive period for target observation for not only SST forcing,but also for other relevant atmospheric variables,such as,initial wind field,initial temperature field,initial pressure field,and initial moisture field associated with TC intensity simulation.Such efforts can then fine tune the strategy for implementing target observations for TC intensity simulations.These considerations represent our subject areas for future investigations.

    Acknowledgements

    .The authors appreciate the anonymous reviewers very much for their very useful comments and suggestions.This work was jointly sponsored by the National Nature Scientific Foundations of China (Grant Nos.41930971) and the National Key Research and Development Program of China (Grant No.2018YFC1506402) and the National Nature Scientific Foundations of China (Grant No.41575061).

    一级毛片 在线播放| 国产熟女午夜一区二区三区 | 国产女主播在线喷水免费视频网站| 欧美精品人与动牲交sv欧美| 一本—道久久a久久精品蜜桃钙片| 噜噜噜噜噜久久久久久91| 日日爽夜夜爽网站| 亚洲精品第二区| 久久久久久久久久久丰满| a级一级毛片免费在线观看| 五月开心婷婷网| 亚洲精品乱码久久久久久按摩| 亚洲欧美成人综合另类久久久| 国产国拍精品亚洲av在线观看| 一级片'在线观看视频| 午夜影院在线不卡| 中文字幕精品免费在线观看视频 | 丝瓜视频免费看黄片| 性色av一级| 亚洲av综合色区一区| 国产69精品久久久久777片| 国产精品国产三级国产专区5o| 黄色视频在线播放观看不卡| 人妻少妇偷人精品九色| 人妻少妇偷人精品九色| 狂野欧美白嫩少妇大欣赏| 色婷婷久久久亚洲欧美| 国产老妇伦熟女老妇高清| 亚洲欧美日韩卡通动漫| 国产乱来视频区| 涩涩av久久男人的天堂| 99国产精品免费福利视频| 成人毛片60女人毛片免费| 国产极品粉嫩免费观看在线 | 国产精品一区二区在线不卡| 欧美日韩av久久| 亚洲av不卡在线观看| 两个人免费观看高清视频 | 曰老女人黄片| 99热这里只有是精品50| 九九在线视频观看精品| 亚洲综合色惰| 亚洲国产最新在线播放| 十八禁高潮呻吟视频 | av不卡在线播放| 欧美成人精品欧美一级黄| 久久精品熟女亚洲av麻豆精品| 老司机影院毛片| 国产熟女午夜一区二区三区 | 一本一本综合久久| 乱人伦中国视频| 少妇的逼水好多| 黑丝袜美女国产一区| 精华霜和精华液先用哪个| 亚洲精品国产成人久久av| 欧美 亚洲 国产 日韩一| 黑人巨大精品欧美一区二区蜜桃 | 亚洲人成网站在线观看播放| 边亲边吃奶的免费视频| 国产成人免费观看mmmm| 最后的刺客免费高清国语| 欧美成人精品欧美一级黄| 久久99一区二区三区| 少妇被粗大猛烈的视频| 噜噜噜噜噜久久久久久91| 国产精品一区www在线观看| 中文在线观看免费www的网站| 国产黄片美女视频| 国产一区二区在线观看av| 日日摸夜夜添夜夜爱| 青春草国产在线视频| 国产淫语在线视频| 亚洲性久久影院| 成人毛片60女人毛片免费| 十分钟在线观看高清视频www | 国产免费又黄又爽又色| 永久网站在线| 又大又黄又爽视频免费| 黄色毛片三级朝国网站 | 日韩不卡一区二区三区视频在线| 亚洲欧美日韩卡通动漫| 中文字幕亚洲精品专区| 久久精品国产亚洲av涩爱| 男女啪啪激烈高潮av片| 一边亲一边摸免费视频| 国产精品久久久久久精品电影小说| 日韩一区二区三区影片| 熟女av电影| 日本黄大片高清| 欧美精品一区二区大全| av一本久久久久| 丝袜脚勾引网站| 久久韩国三级中文字幕| 91成人精品电影| av有码第一页| 极品人妻少妇av视频| 免费观看无遮挡的男女| 中文在线观看免费www的网站| 人人妻人人澡人人爽人人夜夜| 日本色播在线视频| 91精品国产国语对白视频| 丝袜在线中文字幕| 看免费成人av毛片| 精品国产乱码久久久久久小说| 视频中文字幕在线观看| 国产精品久久久久久久久免| 三级经典国产精品| 亚洲一级一片aⅴ在线观看| 午夜91福利影院| 久久精品国产a三级三级三级| 欧美激情极品国产一区二区三区 | 亚洲av欧美aⅴ国产| 九色成人免费人妻av| 女性被躁到高潮视频| 免费少妇av软件| 午夜福利网站1000一区二区三区| 国产国拍精品亚洲av在线观看| a 毛片基地| 秋霞伦理黄片| 日本与韩国留学比较| 国国产精品蜜臀av免费| 久久精品夜色国产| 人人妻人人爽人人添夜夜欢视频 | 国产精品女同一区二区软件| 另类精品久久| 国产黄片美女视频| 日韩大片免费观看网站| 日韩亚洲欧美综合| 80岁老熟妇乱子伦牲交| 国产永久视频网站| 大码成人一级视频| 老司机影院成人| 亚洲三级黄色毛片| 久久99精品国语久久久| 国产精品女同一区二区软件| 晚上一个人看的免费电影| 看免费成人av毛片| 五月玫瑰六月丁香| 少妇 在线观看| 亚洲精品一二三| 亚洲欧美精品专区久久| 亚洲成人手机| 99久久精品一区二区三区| 日韩人妻高清精品专区| 97在线视频观看| 亚洲欧美日韩另类电影网站| 久久鲁丝午夜福利片| 26uuu在线亚洲综合色| 高清不卡的av网站| 观看免费一级毛片| 日日啪夜夜撸| 人人妻人人添人人爽欧美一区卜| 嘟嘟电影网在线观看| 国产欧美亚洲国产| 国产男人的电影天堂91| 人妻少妇偷人精品九色| 久久亚洲国产成人精品v| 亚洲天堂av无毛| 日日摸夜夜添夜夜爱| 涩涩av久久男人的天堂| 亚洲精品久久久久久婷婷小说| 国产精品一区二区在线不卡| 高清在线视频一区二区三区| 国产永久视频网站| 青青草视频在线视频观看| 国产精品久久久久久久久免| 三级国产精品欧美在线观看| 天堂8中文在线网| 三上悠亚av全集在线观看 | 丝袜脚勾引网站| 熟女av电影| 欧美bdsm另类| 日韩亚洲欧美综合| 日本黄色日本黄色录像| 日韩制服骚丝袜av| 久久人人爽人人片av| 亚洲欧美日韩东京热| 精品一区二区免费观看| 99热6这里只有精品| 尾随美女入室| 99热国产这里只有精品6| 六月丁香七月| 一级毛片aaaaaa免费看小| 国产黄色免费在线视频| 国产黄片视频在线免费观看| 另类精品久久| 寂寞人妻少妇视频99o| 有码 亚洲区| 久久人人爽人人片av| 欧美激情极品国产一区二区三区 | 亚洲av成人精品一二三区| 亚洲综合精品二区| 人人妻人人看人人澡| 欧美日韩在线观看h| 毛片一级片免费看久久久久| 青春草视频在线免费观看| 亚洲激情五月婷婷啪啪| 亚洲高清免费不卡视频| 两个人免费观看高清视频 | 如何舔出高潮| 欧美日韩一区二区视频在线观看视频在线| 成人影院久久| 欧美3d第一页| 欧美最新免费一区二区三区| 免费久久久久久久精品成人欧美视频 | av天堂久久9| 亚洲精品一区蜜桃| 妹子高潮喷水视频| 少妇精品久久久久久久| 伊人久久精品亚洲午夜| 日韩av不卡免费在线播放| 美女国产视频在线观看| 久久精品久久精品一区二区三区| 亚洲精品成人av观看孕妇| 男人添女人高潮全过程视频| 寂寞人妻少妇视频99o| 久久久国产一区二区| 午夜激情福利司机影院| 亚洲丝袜综合中文字幕| 欧美区成人在线视频| 亚洲美女视频黄频| 丰满迷人的少妇在线观看| 一级黄片播放器| 日韩一区二区视频免费看| 在线看a的网站| 国产成人免费观看mmmm| 精品少妇内射三级| 国产深夜福利视频在线观看| 老司机影院毛片| www.色视频.com| 精品少妇久久久久久888优播| 一级av片app| 99视频精品全部免费 在线| 伦理电影免费视频| 久久久a久久爽久久v久久| 亚洲一区二区三区欧美精品| 免费观看的影片在线观看| 亚洲丝袜综合中文字幕| 久久精品国产自在天天线| 春色校园在线视频观看| 久久久久久久亚洲中文字幕| 亚洲电影在线观看av| 亚洲欧洲精品一区二区精品久久久 | 街头女战士在线观看网站| 国产成人a∨麻豆精品| 亚洲三级黄色毛片| 最新中文字幕久久久久| 51国产日韩欧美| 黄色一级大片看看| 丰满乱子伦码专区| videos熟女内射| 亚洲精品亚洲一区二区| 精品一区在线观看国产| 大片免费播放器 马上看| 丰满少妇做爰视频| 欧美成人精品欧美一级黄| 黄色一级大片看看| 美女主播在线视频| 91久久精品国产一区二区三区| 久久精品夜色国产| 熟女av电影| 国产真实伦视频高清在线观看| 久久国产乱子免费精品| 又粗又硬又长又爽又黄的视频| 国产精品一区二区在线不卡| 最近手机中文字幕大全| 乱码一卡2卡4卡精品| 97超碰精品成人国产| 赤兔流量卡办理| 国产在线视频一区二区| 亚洲熟女精品中文字幕| 精品一区在线观看国产| 少妇熟女欧美另类| 尾随美女入室| 看非洲黑人一级黄片| 卡戴珊不雅视频在线播放| 超碰97精品在线观看| 久久久欧美国产精品| 国内少妇人妻偷人精品xxx网站| 又粗又硬又长又爽又黄的视频| 久久久久网色| 免费人成在线观看视频色| 人妻夜夜爽99麻豆av| 日韩大片免费观看网站| 97在线视频观看| 国产精品99久久99久久久不卡 | 日本与韩国留学比较| 国产亚洲精品久久久com| 亚洲欧美日韩东京热| 国产亚洲午夜精品一区二区久久| 自线自在国产av| a级毛片在线看网站| 91精品伊人久久大香线蕉| 欧美变态另类bdsm刘玥| 亚洲国产欧美在线一区| 丁香六月天网| 最黄视频免费看| 亚洲国产色片| 伦理电影大哥的女人| 亚洲精品国产色婷婷电影| 免费黄频网站在线观看国产| 深夜a级毛片| 成年人午夜在线观看视频| 九草在线视频观看| 免费大片18禁| 在线播放无遮挡| 精品午夜福利在线看| 久久久久久久大尺度免费视频| av网站免费在线观看视频| 能在线免费看毛片的网站| 国产毛片在线视频| 男女边吃奶边做爰视频| 极品少妇高潮喷水抽搐| 久久精品国产自在天天线| 在线免费观看不下载黄p国产| 99精国产麻豆久久婷婷| av卡一久久| 91精品一卡2卡3卡4卡| 一本色道久久久久久精品综合| 亚洲精品中文字幕在线视频 | 亚洲高清免费不卡视频| 黄色一级大片看看| 亚洲精品久久午夜乱码| 国产精品国产av在线观看| 国产一区亚洲一区在线观看| 亚洲国产欧美日韩在线播放 | 免费在线观看成人毛片| 精品人妻偷拍中文字幕| 亚洲va在线va天堂va国产| 99国产精品免费福利视频| 97在线视频观看| 久久久a久久爽久久v久久| 看非洲黑人一级黄片| 国产乱来视频区| 亚洲欧美一区二区三区黑人 | 午夜福利在线观看免费完整高清在| 91久久精品国产一区二区成人| 中文资源天堂在线| 黑人猛操日本美女一级片| 老女人水多毛片| 免费观看a级毛片全部| 久久精品夜色国产| 久久久久久久久久久久大奶| 成年人午夜在线观看视频| 精品国产一区二区三区久久久樱花| 久久精品久久精品一区二区三区| 欧美激情国产日韩精品一区| 欧美性感艳星| 亚洲欧美成人精品一区二区| 一级a做视频免费观看| 97在线人人人人妻| 免费播放大片免费观看视频在线观看| 欧美日韩亚洲高清精品| 国产精品福利在线免费观看| 青春草国产在线视频| 交换朋友夫妻互换小说| 国产淫片久久久久久久久| 精品国产露脸久久av麻豆| 青青草视频在线视频观看| 久久毛片免费看一区二区三区| 国产亚洲91精品色在线| 亚洲国产毛片av蜜桃av| 免费不卡的大黄色大毛片视频在线观看| 亚洲国产成人一精品久久久| 欧美日韩精品成人综合77777| 国产精品一二三区在线看| 午夜福利网站1000一区二区三区| 我要看黄色一级片免费的| 欧美最新免费一区二区三区| 王馨瑶露胸无遮挡在线观看| 夫妻性生交免费视频一级片| 久久国产亚洲av麻豆专区| 亚洲无线观看免费| 久久久久人妻精品一区果冻| 欧美老熟妇乱子伦牲交| 91久久精品国产一区二区三区| 国产淫语在线视频| 九色成人免费人妻av| 国产成人精品一,二区| 美女国产视频在线观看| videos熟女内射| 国产一区二区在线观看av| 国产爽快片一区二区三区| 免费观看a级毛片全部| 日本黄色片子视频| 国产精品成人在线| 欧美3d第一页| 一区在线观看完整版| 国产69精品久久久久777片| 这个男人来自地球电影免费观看 | 国产又色又爽无遮挡免| av天堂久久9| 另类亚洲欧美激情| 久久久久人妻精品一区果冻| 亚洲无线观看免费| 观看美女的网站| kizo精华| 欧美性感艳星| 十八禁网站网址无遮挡 | 久久国产精品大桥未久av | 在线观看av片永久免费下载| 中文乱码字字幕精品一区二区三区| 久久人妻熟女aⅴ| 午夜老司机福利剧场| 女性生殖器流出的白浆| 91精品国产国语对白视频| 九草在线视频观看| 热re99久久国产66热| 成年人免费黄色播放视频 | 80岁老熟妇乱子伦牲交| 亚洲精品国产色婷婷电影| 欧美变态另类bdsm刘玥| 我要看日韩黄色一级片| 亚洲欧美成人精品一区二区| 久久午夜综合久久蜜桃| 久久久久久久久久成人| 日日爽夜夜爽网站| 亚洲成色77777| 99久久综合免费| 国产精品一区二区在线观看99| 91精品伊人久久大香线蕉| 国产黄片美女视频| 能在线免费看毛片的网站| 99视频精品全部免费 在线| 一区二区三区乱码不卡18| 在线天堂最新版资源| 婷婷色综合www| 夫妻性生交免费视频一级片| 最近手机中文字幕大全| 久久99一区二区三区| 麻豆成人午夜福利视频| 纵有疾风起免费观看全集完整版| 国产精品99久久99久久久不卡 | 亚洲av电影在线观看一区二区三区| 美女中出高潮动态图| 国产精品免费大片| 精品一品国产午夜福利视频| 国产一区亚洲一区在线观看| h视频一区二区三区| 亚州av有码| 欧美精品高潮呻吟av久久| 一本色道久久久久久精品综合| 国产色爽女视频免费观看| 亚州av有码| 成人综合一区亚洲| 激情五月婷婷亚洲| 91精品一卡2卡3卡4卡| 国产精品女同一区二区软件| 九草在线视频观看| videossex国产| 肉色欧美久久久久久久蜜桃| 曰老女人黄片| 欧美亚洲 丝袜 人妻 在线| 特大巨黑吊av在线直播| 少妇猛男粗大的猛烈进出视频| 国产在线男女| 日韩伦理黄色片| 亚洲欧美日韩另类电影网站| 大话2 男鬼变身卡| 国产成人免费无遮挡视频| 日韩欧美一区视频在线观看 | 日本欧美国产在线视频| h日本视频在线播放| 亚洲精品国产成人久久av| 久久精品熟女亚洲av麻豆精品| 亚洲av男天堂| 国产亚洲一区二区精品| 黑丝袜美女国产一区| 天美传媒精品一区二区| 成人综合一区亚洲| 日本wwww免费看| 超碰97精品在线观看| 成人毛片a级毛片在线播放| 好男人视频免费观看在线| 久久久久久久久久成人| 欧美国产精品一级二级三级 | 色婷婷av一区二区三区视频| 久久久久久久大尺度免费视频| 这个男人来自地球电影免费观看 | 久久女婷五月综合色啪小说| 国产精品一区二区在线不卡| 男女无遮挡免费网站观看| 久久久久久久久大av| 国产69精品久久久久777片| 国产免费一级a男人的天堂| 观看美女的网站| 国产欧美亚洲国产| 国产在视频线精品| 伊人亚洲综合成人网| 国产视频内射| 日韩不卡一区二区三区视频在线| 欧美精品高潮呻吟av久久| 国产亚洲午夜精品一区二区久久| 国产成人一区二区在线| 国产黄片视频在线免费观看| 国产男女超爽视频在线观看| 国产淫片久久久久久久久| 蜜臀久久99精品久久宅男| 国产精品久久久久久精品古装| 少妇人妻 视频| 观看免费一级毛片| 欧美另类一区| 国产毛片在线视频| 18禁在线无遮挡免费观看视频| 日本黄大片高清| 狂野欧美白嫩少妇大欣赏| 国产真实伦视频高清在线观看| 亚洲美女搞黄在线观看| 五月玫瑰六月丁香| 熟女av电影| 卡戴珊不雅视频在线播放| 免费观看av网站的网址| 免费黄色在线免费观看| 精华霜和精华液先用哪个| 精品亚洲成a人片在线观看| 午夜激情久久久久久久| 麻豆成人av视频| 啦啦啦视频在线资源免费观看| 全区人妻精品视频| 精品人妻熟女av久视频| 日本午夜av视频| 久久精品国产鲁丝片午夜精品| 在线观看三级黄色| 尾随美女入室| 国模一区二区三区四区视频| 午夜福利网站1000一区二区三区| 国产无遮挡羞羞视频在线观看| 菩萨蛮人人尽说江南好唐韦庄| 夜夜爽夜夜爽视频| 搡老乐熟女国产| 看十八女毛片水多多多| 免费观看的影片在线观看| 午夜日本视频在线| 成年av动漫网址| 亚洲精品乱久久久久久| 少妇被粗大的猛进出69影院 | 草草在线视频免费看| 青春草视频在线免费观看| 免费大片18禁| 九九在线视频观看精品| 亚洲精品第二区| 久久精品久久久久久噜噜老黄| 日本欧美视频一区| 一级毛片久久久久久久久女| 黄色日韩在线| 国产片特级美女逼逼视频| 国产精品无大码| 精品亚洲成国产av| 狂野欧美激情性xxxx在线观看| 国产一区二区三区综合在线观看 | av卡一久久| 国产女主播在线喷水免费视频网站| 免费看av在线观看网站| 最近最新中文字幕免费大全7| 欧美一级a爱片免费观看看| 美女xxoo啪啪120秒动态图| 国产高清不卡午夜福利| 男人和女人高潮做爰伦理| 日日啪夜夜撸| 国产乱人偷精品视频| 精品久久久久久久久av| 精品人妻熟女毛片av久久网站| 高清欧美精品videossex| 这个男人来自地球电影免费观看 | 日本免费在线观看一区| 久久精品国产亚洲网站| 日本vs欧美在线观看视频 | 毛片一级片免费看久久久久| 欧美精品人与动牲交sv欧美| 久久久久久久国产电影| 精品久久久久久久久av| 国产在线免费精品| 国产精品熟女久久久久浪| 少妇人妻一区二区三区视频| 亚洲国产毛片av蜜桃av| 国产综合精华液| 全区人妻精品视频| 日本黄色片子视频| 国产精品国产三级国产av玫瑰| 国产伦理片在线播放av一区| 日韩免费高清中文字幕av| 国产精品偷伦视频观看了| 91成人精品电影| 一级二级三级毛片免费看| 国产有黄有色有爽视频| 亚洲精品第二区| 亚洲国产精品专区欧美| 日韩在线高清观看一区二区三区| 日本av手机在线免费观看| 中文字幕人妻丝袜制服| 韩国高清视频一区二区三区| 国产精品国产三级国产av玫瑰| 人妻一区二区av| 久久国内精品自在自线图片| 少妇人妻一区二区三区视频| 久久国产精品男人的天堂亚洲 | 中国国产av一级| 国产av码专区亚洲av| 韩国高清视频一区二区三区| 日韩三级伦理在线观看| 十八禁网站网址无遮挡 | 久久国产乱子免费精品| 日本欧美视频一区| 久久97久久精品| 美女大奶头黄色视频| 免费黄网站久久成人精品| 一级毛片电影观看| 国产高清不卡午夜福利| a级毛色黄片| 国产欧美另类精品又又久久亚洲欧美| 久久精品国产亚洲av涩爱| 亚洲国产色片| 搡女人真爽免费视频火全软件| 18+在线观看网站| 亚洲欧美清纯卡通|