CN112596127A - Novel method for calculating typhoon potential generation index - Google Patents

Novel method for calculating typhoon potential generation index Download PDF

Info

Publication number
CN112596127A
CN112596127A CN202011336283.9A CN202011336283A CN112596127A CN 112596127 A CN112596127 A CN 112596127A CN 202011336283 A CN202011336283 A CN 202011336283A CN 112596127 A CN112596127 A CN 112596127A
Authority
CN
China
Prior art keywords
typhoon
wind field
wind
field
atmospheric
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202011336283.9A
Other languages
Chinese (zh)
Inventor
哈瑶
钟中
朱益民
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
National University of Defense Technology
Original Assignee
National University of Defense Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN202011336283.9A priority Critical patent/CN112596127A/en
Publication of CN112596127A publication Critical patent/CN112596127A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01WMETEOROLOGY
    • G01W1/00Meteorology
    • G01W1/10Devices for predicting weather conditions

Abstract

The invention relates to a new method for calculating a typhoon potential generation index, which is characterized in that a typhoon wind field is quantitatively separated from atmospheric reanalysis data, and 850hPa absolute vorticity of a troposphere low layer and a vertical wind shear of 200hPa minus 850hPa of a troposphere high and low layer are calculated; calculating 600hPa relative humidity and typhoon potential intensity; and inputting the calculation result into the original GPI index so as to finally obtain the improved GPI index. The novel method provided by the invention optimizes the traditional potential generation index GPI algorithm on the basis of physical significance and a calculation scheme, and calculates the GPI index by using the element variables without typhoon weather system data such as a typhoon wind field and the like, thereby obtaining the index of potential contribution of a purer atmospheric marine environment background field to typhoon generation.

Description

Novel method for calculating typhoon potential generation index
Technical Field
The invention relates to the technical field of typhoon climate science, in particular to a novel method for calculating a typhoon potential generation index.
Background
Typhoon (also known as Tropical Cyclone (TC)) is often used to estimate the combined impact of various environmental field elements in the atmospheric ocean on the generation of marine-generated weather scale atmospheric systems, especially characterizing the strong convective and dangerous weather system generation contributions like typhoons.
Typhoons are mainly generated and developed on tropical ocean surfaces in summer and autumn, and for typhoons affecting China, most of the typhoons move to the west, the northwest or the north after being generated in the northwest pacific and the south sea of the tropical sea, and a few of the typhoons move to the east or the south. The typhoon is likely to be generated and continuously developed and enhanced under certain atmospheric marine environmental field conditions, the classical atmospheric science theory considers that key factors influencing the generation and development of the typhoon comprise a wide warm ocean surface, small vertical shear of the wind speed of a troposphere, a rotating parameter larger than a certain numerical value and certain small disturbance to a low layer of the troposphere, and after the conditions are met, the small disturbance on the tropical ocean surface can be used as an embryo of the typhoon development and can develop and strengthen in the next few days, namely the typhoon generation. Therefore, according to the physical conditions of the atmospheric environmental field, the factors of the meteorological environment that affect the generation and development of the typhoon mainly include the sea surface temperature, the relative vorticity of the troposphere low layer, the relative humidity of the troposphere, the vertical movement of the atmosphere, the vertical wind shear, and the like.
Based on the basic theory, the concept of tropical cyclone potential generation index (GPI) is proposed by two meteorologists, Emanuel and Nolan, and the index can better describe the climatic state space-time characteristics generated by typhoon, and the GPI expression of the index is shown in the following section. This index relates to several large scale environmental factors closely related to typhoon generation, including the absolute vorticity (η) of 850hPa in the troposphere, the relative humidity (RH-relative humidity) of 600hPa in the troposphere, the cyclone potential strength (V)pot) Wind shear (V) of 200hPa and 850hPashear-vertical shear). The parameter of potential strength is a diagnostic variable which depends on the difference between the temperature at which the sea gas is thermally unstable and the adiabatically rising boundary layer gas block reaches neutral buoyancy and the sea surface temperature, which relates to the sea surface temperature, sea level air pressure, vertical atmospheric temperature and the mixing ratio, and which later takes into account the effect of dissipative heating. The GPI index is applied to typhoon generation and mobile diagnosis analysis and forecast prediction, and is also applied to output data information of various weather modes for researching the weather change of tropical cyclone. A large number of results confirm that the pattern results and the observation results show a certain degree of coincidence in the seasonal variation and the generation position of typhoon, and thus prove to have a certain degree of influence on typhoon activitiesCapabilities are described.
The most similar technical scheme of the invention is the calculation mode of the above mentioned classic GPI, and the expression and calculation method of the GPI comprises the following steps:
Figure BDA0002797232000000021
where eta is 850hPa absolute vorticity, H is 600hPa relative humidity, VpotIs the typhoon latent strength defined by Emanuel, which is determined by Sea Surface Temperature (SST), air pressure in the vertical direction, air temperature and specific humidity and can be accurately calculated, VshearAre the magnitudes of 200hPa and 850hPa vertical wind shear. The magnitude of the GPI calculation result respectively represents the influence of the environmental conditions of the atmospheric ocean elements on the typhoon generating activity. The large numerical value represents strong contribution, typhoon is easy to generate and develop, and frequently occurs to a scale system of the flowing weather; the small numerical value indicates weak contribution, activities such as typhoon generation and the like are restrained by environmental conditions, and strong convection weather such as typhoon and the like is not easy to develop. All data (such as wind field, humidity field, temperature field of each level of the atmosphere, sea surface temperature and the like) for calculating the GPI index can use conventional reanalysis data in the field of atmospheric science research or output data of various atmospheric modes. For example, the national center for atmospheric environmental research (NCEP/NCAR) re-analyzes data sets including high altitude 17-layer wind fields, relative humidity, vertical velocity, and temperature fields, as well as reconstructed sea surface temperature data, global sea surface temperature grid point data, provided by the National Oceanic and Atmospheric Administration (NOAA).
The traditional technical scheme is simpler and more intuitive. The data set used in the field of atmospheric science contains all parameter variables in GPI formula, and the data of all kinds of elements are directly substituted into the GPI formula, so that the result of classical GPI index space-time distribution can be calculated.
In the existing method for calculating the GPI index, the physical significance and the calculation mode of part of calculation item parameters are not accurate enough, so that the calculated GPI cannot completely and truly reflect the contribution of atmospheric marine environment field conditions generated by typhoon. The main points are as follows: the GPI index contains 2 important parameters, namely the tropospheric low-rise 850hPa absolute vorticity and the 200hPa minus 850hPa vertical wind shear (200hPa minus 850hPa latitudinal wind). Because the calculation of the contribution of each meteorological ocean element in the typhoon generation environment diagnosis is usually directed at the periods of excessive typhoon generation and active activity in summer, autumn and the like, at this time, when the absolute vorticity item of 850hPa of troposphere low-level and the shear of 850hPa vertical wind subtracted from 200hPa are calculated, the wind field of the troposphere high-level and low-level analysis data already contains a large amount of wind field information of typhoon systems, and particularly in the areas and periods of strong typhoon and ultra-strong typhoon activity, the influence of the wind field on the calculation result is more obvious, so that the GPI index calculated under the condition is actually the possibility of 'diagnosing' generation of the typhoon systems by using the change characteristic of the typhoon systems, the diagnosis result has larger deviation, and even has wrong logic contrary to the 'potential generation index' design original purpose. The index analysis is utilized to research the related problems of typhoon activities in the atmospheric science field, so that the contribution of an environmental field to typhoon generation is exaggerated in the analysis result, and meanwhile, the judgment of the potential typhoon generation possibility is deviated, so that the accuracy of the contribution diagnosis of the environmental field to typhoon generation is greatly reduced.
Disclosure of Invention
The invention aims to provide a novel method for calculating a typhoon potential generation index, which is characterized in that when a GPI index is calculated, the influence from a typhoon system is removed from the terms of 850hPa absolute vorticity of a convection layer and 200hPa minus 850hPa vertical wind shear of a high and low layer of the convection layer, and the GPI index is calculated by using element variables which do not contain typhoon weather system data such as a typhoon field and the like, so that an index of potential contribution of a purer atmospheric ocean environment background field to typhoon generation is obtained, and the problems in the background technology are solved.
In order to achieve the purpose, the technical scheme of the invention is as follows:
a new method of calculating a typhoon potential generation index, the new method comprising the steps of:
the first step is as follows: obtaining a whole-layer wind field structure of an atmospheric circulation background wind field from atmospheric reanalysis data, wherein the whole-layer wind field structure comprises a typhoon wind field component and an atmospheric circulation background wind field component, quantitatively separating the typhoon wind field from the atmospheric reanalysis data, the typhoon wind field comprises a non-rotational motion wind field and a non-dispersive motion wind field, and removing the typhoon wind field from the whole-layer wind field structure;
the second step is that: computing absolute vorticity eta of troposphere lower layer 850hPa by utilizing atmospheric circulation background wind fieldNOTCAnd convection layer high and low layer 200hPa minus 850hPa vertical wind shear Vshear-NOTC
Flow layer low layer 850hPa absolute vorticity etaNOTCThe calculation formula is as follows:
Figure BDA0002797232000000041
wherein
Figure BDA0002797232000000042
And
Figure BDA0002797232000000043
respectively showing the latitude wind and the longitude wind of the atmospheric circulation background after eliminating the typhoon wind field on the 850hPa level;
convection layer high and low layer 200hPa minus 850hPa vertical wind shear Vshear-NOTCThe calculation formula of (2) is as follows:
Figure BDA0002797232000000044
or
Figure BDA0002797232000000045
Wherein
Figure BDA0002797232000000046
And
Figure BDA0002797232000000047
respectively represents the atmospheric ring after eliminating the typhoon field on the 200hPa levelFlow background weft and warp, wherein
Figure BDA0002797232000000048
And
Figure BDA0002797232000000049
respectively showing the latitude wind and the longitude wind of the atmospheric circulation background after eliminating the typhoon wind field on the 850hPa level;
the third step: calculation of 600hPa relative humidity H and typhoon potential intensity V using conventional datapot
The fourth step: using the parameter η obtained in the second stepNOTCAnd Vshear-NOTCAnd the parameters H and V obtained in the third steppotCalculating GPI to finally obtain an improved GPI index;
the GPI calculation formula obtained by the new method is as follows:
Figure BDA0002797232000000051
in the above new process, in the first step: the atmospheric circulation background wind field is quantitatively separated from the atmospheric re-analysis data by a typhoon wind field separation method and is applied to each standard level of the atmospheric re-analysis data.
In the above new process, in the first step: the method for separating the typhoon wind field comprises the steps of determining the center of the position of the typhoon in the analysis data, and then performing wind field decomposition on an area which takes the point as the center of a circle and takes a certain distance as the radius.
In the new method, the radius range of the wind field decomposition is set to 800km, and reanalysis data with a resolution of 1 DEG x 1 DEG or more is used.
In the new method, the whole-layer wind field structure comprises a typhoon wind field component and an atmospheric circulation background wind field component, and the typhoon wind field comprises a non-rotational motion wind field and a non-scattered motion wind field.
In the new method, the principle of the wind field decomposition is that the plane of incompressible fluid has no rotational motion and has velocity potential and flow function, the horizontal wind field in the reanalysis data can be regarded as the superposition of no rotational motion, no scattered motion and residual error terms, the velocity fields without rotational motion and no scattered motion are superposed to obtain the typhoon velocity field by calculating the velocity potential and flow function of the reanalysis data wind field, and the residual error terms of the reanalysis velocity field and the typhoon velocity field are the atmospheric circulation background wind field.
Compared with the prior art, the invention has the beneficial effects that: when the GPI index is calculated, the influence from the typhoon system is removed from the terms of the convection layer 850hPa absolute vorticity and the convection layer high and low layer 200hPa minus 850hPa vertical wind shear, the GPI index is calculated by using the element variables without typhoon weather system data such as a typhoon wind field, and the index of potential contribution of the atmospheric marine environment background field to the typhoon generation is further obtained.
Specifically, the method removes the wind field data information of typhoon in the atmospheric circulation, calculates two parameters of absolute vorticity and vertical wind shear by only using the atmospheric background wind field, and substitutes the parameters into the GPI index. The improved GPI index can be considered as potential contribution of a pure environmental background field to generation and activities of a typhoon weather system due to the fact that the improved GPI index does not contain typhoon information interference, and the indication capability of the improved GPI index to the generation relation of the atmospheric marine environmental conditions, typhoons and other marine strong convection weather scale disturbance systems is greatly improved, so that the purpose of objectively and quantitatively and accurately evaluating the influence of the atmospheric marine background field on the typhoon activities is achieved.
Drawings
FIG. 1 is a flow chart of an embodiment of the present invention;
FIG. 2 is a schematic diagram of a computational process implemented by the present invention;
FIG. 3 is a schematic view of a positioning analysis typhoon center implemented in the present invention;
FIG. 4 is a graph of an 850hPa wind field and a vorticity field of the Western Pacific ocean at 27, 12 months 8, 1997 in which (a) the reanalyzed data fields, (b) the typhoon field, (c) the ambient air circulation background wind field has a contour interval of 10-5s-1(ii) a The average 850hPa wind field and vorticity field in 6-10 months of the northwest Pacific in 1997, (d) reanalysis of the data field, (e) typhoon field, (f) atmospheric circulation background wind field, and the interval between the contour lines is 10-6s-1(ii) a The solid and dashed lines represent positive and negative vorticity, respectively.
Detailed Description
The technical solution of the present invention is further described in detail with reference to the accompanying drawings and examples.
As shown in fig. 1 and 2, a new method of calculating a typhoon potential generation index includes the steps of:
the first step is as follows: and obtaining the whole layer wind field structure of the atmospheric circulation background wind field from the atmospheric re-analysis data, and quantitatively separating the typhoon wind field from the atmospheric re-analysis data.
The whole-layer wind field structure comprises a typhoon wind field component and an atmospheric circulation background wind field component, the typhoon wind field comprises a non-rotary motion wind field and a non-dispersion motion wind field, and finally the typhoon wind field is removed from the whole-layer wind field structure.
The atmosphere reanalysis data serve for numerical weather forecast or weather pattern analysis and the like, the numerical weather forecast is more common, and most of weather forecast seen by people at ordinary times is estimation and judgment of a weather numerical pattern simulation result. The core of the numerical model is a set of atmospheric dynamics equations, and in order to obtain a better result, the solution of the equations needs accurate boundary conditions and initial conditions, which are provided by input data, namely, atmospheric re-analysis data.
In order to eliminate the space-time distribution characteristics of a typhoon wind field and the variation rule of the strength of the space-time distribution characteristics, and obtain the homogeneous distribution of the atmospheric environment under different climatic backgrounds, so that the accurate contribution of factors such as the vorticity of the atmospheric circulation background wind field, the vertical wind shear and the like in the environmental field to the generation and the activity of the typhoon can be more accurately analyzed, the first step of the invention is to remove the typhoon vortex wind field from the atmospheric reanalysis data.
The principle of removing the typhoon vortex wind field is that the incompressible fluid plane can be decomposed without rotating motion, the relevant wind field information of the typhoon can be quantitatively separated from the analysis data after the decomposition, and the method is applied to each standard level of the atmosphere re-analysis data, so that the whole layer wind field structure of the atmosphere circumfluence background wind field can be obtained. The method is a conventional means in typhoon numerical simulation and forecast, and achieves the purpose of simplifying the nonlinear action between atmospheric circulation systems with different scales by decomposing the global wind field data on different atmospheric pressure surfaces into typhoon components and background components.
The specific method for removing the typhoon vortex wind field comprises the steps of firstly determining the center of the position of the typhoon in the analysis data, and then carrying out wind field decomposition on an area which takes the point as the center of a circle and takes a certain distance as the radius. The process of determining the center of the analytic typhoon is shown in fig. 3, and the grid point of the reanalysis data closest to the positioning position of the optimal path set of the typhoon in the analytic data is used as the real typhoon center, the grid point with the maximum positive vorticity on the reanalysis data is searched by taking 600km as the radius (r _ search), the searched grid point is used as the center of the analytic typhoon, and the atmosphere reanalysis data is decomposed into typhoon wind field components and background wind field components within the range of 800km radius (r _ vor).
The gray triangles in FIG. 3 indicate the position of the best path set data positioning typhoon, and since this position is the actual position of the typhoon activity, its probability is not at the grid point of the analysis data. The gray large round point represents the center of the real typhoon, namely the grid point closest to the position of the typhoon in the optimal path set, the black large round point represents the center of the searched analysis typhoon, the dotted line ring area represents the maximum vorticity search area, the solid line ring area represents the analysis typhoon area, the black small round point represents the grid point position where the analysis data is decomposed, and the decomposed analysis field is located on the black small round point in the large circle. In the figure, the reanalysis data with 2.5 DEG x 2.5 DEG spatial resolution is taken as an example, and the method is similar if reanalysis data with 1 DEG x 1 DEG spatial resolution is used.
After the wind field range is determined, the wind field decomposition can be carried out on the target area, and a background wind field and a typhoon wind field are separated.
The principle of decomposing and analyzing the data wind field is that the plane of incompressible fluid has no rotational motion and has velocity potential and flow function, the horizontal wind field in the reanalysis data can be regarded as the superposition of the non-rotational motion wind field, the non-scattered motion wind field and the residual error term, the velocity fields of the non-rotational motion wind field and the non-scattered motion wind field can be superposed to obtain the typhoon wind field by calculating the velocity potential and the flow function of the reanalysis data wind field, and the residual error term of the reanalysis velocity field and the typhoon wind field is the atmospheric circulation background wind field.
The specific calculation scheme is that in the range of 800km radius with the center of the typhoon to be analyzed as the center of a circle, a formula is utilized
Figure BDA0002797232000000081
And
Figure BDA0002797232000000082
calculating a non-dispersion motion wind field, wherein psi represents a flow function of non-dispersion wind, and zeta represents relative vorticity;
Figure BDA0002797232000000083
is a translational motion velocity field. At the same time, using the formula
Figure BDA0002797232000000084
And
Figure BDA0002797232000000085
calculating a non-rotational motion wind field, wherein x represents the velocity potential, delta represents the divergence,
Figure BDA00027972320000000812
is a rotational motion velocity field. Typhoon wind field VTC=vψ+vχEliminating environmental background wind field V after typhoon wind fieldE=V-VTC
The concrete solving steps are as follows:
s1, atmospheric wind field of known meteorological data reanalysis data
Figure BDA0002797232000000086
(radial wind component v and latitudinal wind component u), calculating relative vorticity ζ and divergence δ:
Figure BDA0002797232000000087
s2, according to
Figure BDA0002797232000000088
And
Figure BDA0002797232000000089
solving psi and chi;
s3, according to
Figure BDA00027972320000000810
And
Figure BDA00027972320000000811
v is obtainedψAnd vχ(ii) a (k is a unit vector in the z direction)
S4,VTC=vψ+vχFurther, V is obtainedE=V-VTC
The above calculation method is prior art and is used herein to verify the feasibility of the separation typhoon technique.
The specific flow of the typhoon wind field part and the atmospheric circulation background wind field part is divided into meteorological data by utilizing a typhoon wind field separating technology. The analysis data of each standard layer in the atmosphere vertical direction are processed by the method, and the typhoon wind field component of the whole layer of each time-order troposphere and the atmospheric circulation background wind field component can be obtained. The present inventor utilizes Fortran program source code to achieve the purpose of decomposing the atmospheric re-analysis data into typhoon vortex wind field component and atmospheric circulation background wind field component.
It should be noted that, since reanalysis data used in atmospheric science research at present have different spatial resolutions (e.g., reanalysis data of resolutions such as 2.5 ° × 2.5 °, 1 ° × 1 °), when the reanalysis data has low spatial resolution, the spatial resolution of the decomposed typhoon field component is also low, so that the representativeness of the information such as typhoon intensity, influence range and the like is insufficient, in order to solve the problem of weak climatological characteristics of the typhoon field component possibly caused by the method to a certain extent, the invention sets the radius range of the decomposition analysis data field to be a larger threshold (for example 800km), and simultaneously, reanalysis data with higher resolution (such as 1 degree multiplied by 1 degree and above) is proposed to be used, so that more typhoon wind field information in the grid point data is reflected in the separated typhoon components, and a typhoon vortex wind field and an atmospheric environment field background wind field can be more accurately represented. For the selection of the analysis area threshold, the spatial range of the decomposed atmospheric wind field can be determined by adopting a mode of dynamically positioning and analyzing the typhoon scale according to the satellite observation data, so that the purpose of eliminating the typhoon system in the environmental field more accurately is achieved.
In order to verify the effectiveness of the method for eliminating the typhoon field, the inventor provides experimental implementation data. Fig. 4 shows the state before and after decomposing the 850hPa wind field of the troposphere low layer in the northwest pacific and south sea areas and the vorticity field thereof into typhoon vortex and atmospheric background circulation at 8, 27, 12 in 1997.
As can be seen from fig. 4(a), three typhoon systems exist simultaneously on the positive vorticity zone distributed in the latitudinal direction around 20 ° N at this moment, which is also the position of the east asian monsoon wind channel, which is a typical potential field favorable for the generation and development of typhoons, so that the typhoon area is the maximum area of the positive vorticity in the monsoon wind channel. Fig. 4(b) shows the 850hPa wind field and the vorticity field with three typhoons separated, and a comparison with fig. 4(a) shows that the positive vorticity is mainly caused by typhoons. Fig. 4(c) shows the background field portion after the typhoon vortex is removed, and it can be seen that the positive vorticity in the monsoon groove is greatly reduced, but the relevant features of the monsoon groove still remain in the background data.
To illustrate the effect of decomposition on long-time scale climate characterization, FIG. 4(d) shows the mean distribution of 850hPa wind field and vorticity field in the season of typhoon (6-10 months) in 1997. FIG. 4(e) is a seasonal average of isolated typhoon circulation, which is mainly manifested by a significant positive vorticity distribution in the southeast region of the North-West Pacific, which is most active in typhoon activity in the summer and autumn, corresponding to the strongest early nino in 1997
Figure BDA0002797232000000101
The modulation effect on the activity of the northwest Pacific typhoon can be accurately reversedReflecting the climatic properties of the activity in the typhoon season (note: El)
Figure BDA0002797232000000102
During the year, western pacific typhoons frequently move and strong typhoons frequently occur). Fig. 4(f) shows the seasonal average of the background field with the typhoon field removed, and compared to fig. 4(e), the wind speed and positive vorticity in the southeast region of the pacific northwest are reduced by the seasonal average of the typhoon field, while the main large-scale circulating system is still evident. Therefore, the re-analysis data is processed by the processing method, so that the purposes of separating the typhoon related information from the analysis field and effectively retaining the large-scale circulation characteristics in the background field are achieved.
The second step is that: computing absolute vorticity eta of troposphere lower layer 850hPa by utilizing atmospheric circulation background wind fieldNOTCAnd convection layer high and low layer 200hPa minus 850hPa vertical wind shear Vshear-NOTC
The first step is utilized to obtain the data of the large-scale environmental wind field of each layer of the troposphere after eliminating the typhoon wind field, and the absolute vorticity eta of 850hPa can be calculatedNOTCAnd convection layer high and low layer 200hPa minus 850hPa vertical wind shear Vshear-NOTCThese two parameters. Wherein, the absolute vorticity eta of the flow layer lower layer 850hPaNOTCThe calculation formula is as follows:
Figure BDA0002797232000000103
wherein therein
Figure BDA0002797232000000104
And
Figure BDA0002797232000000105
respectively showing the latitude wind and the longitude wind of the atmospheric circulation background after eliminating the typhoon wind field on the 850hPa level;
convection layer high and low layer 200hPa minus 850hPa vertical wind shear Vshear-NOTCThe calculation formula of (2) is as follows:
Figure BDA0002797232000000106
or
Figure BDA0002797232000000111
Wherein
Figure BDA0002797232000000112
And
Figure BDA0002797232000000113
respectively representing the atmospheric circulation background latitude wind and longitude wind after eliminating the typhoon field on the 200hPa level, wherein
Figure BDA0002797232000000114
And
Figure BDA0002797232000000115
respectively showing the atmospheric circulation background latitude wind and longitude wind after eliminating the typhoon field on the 850hPa level.
The third step: calculation of 600hPa relative humidity H and typhoon potential intensity V using conventional datapot
The relative humidity H of 600hPa does not need to be calculated, and can be directly acquired from meteorological data in a centralized manner. Potential intensity of typhoon VpotThe calculation method can be obtained by referring to the conventional method for calculating the GPI index formula, and all parameters and variables can be directly obtained from meteorological data in a centralized manner.
The fourth step: using the parameter η obtained in the second stepNOTCAnd Vshear-NOTCAnd the parameters H and V obtained in the third steppotCalculating GPI to finally obtain an improved GPI index;
the GPI calculation formula obtained by the new method is as follows:
Figure BDA0002797232000000116
thus, an improved GPI index can be obtained.
The invention aims to eliminate the influence from a typhoon system in terms of 850hPa absolute vorticity of a troposphere and 200hPa minus 850hPa vertical wind shear of a troposphere when calculating a GPI index, and calculate the GPI index by using element variables without typhoon weather system data such as a typhoon wind field and the like, thereby obtaining an index of potential contribution of a more pure atmospheric marine environment background field to typhoon generation.
Specifically, the method removes the wind field data information of typhoon in the atmospheric circulation, calculates two parameters of absolute vorticity and vertical wind shear by only using the atmospheric background wind field, and substitutes the parameters into the GPI index. The improved GPI index can be considered as potential contribution of a pure environmental background field to generation and activities of a typhoon weather system due to the fact that the improved GPI index does not contain typhoon information interference, and the indication capability of the improved GPI index to the generation relation of the atmospheric marine environmental conditions, typhoons and other marine strong convection weather scale disturbance systems is greatly improved, so that the purpose of objectively and quantitatively and accurately evaluating the influence of the atmospheric marine background field on the typhoon activities is achieved.
In practical application, the improved GPI index can be calculated by removing the typhoon wind field and other information from the analysis data and calculating the GPI using the data from which the typhoon wind field is removed, thereby achieving the similar purpose and expected calculation result to those of the present invention. At present, in addition to the method mentioned in the present invention, it is possible to remove the typhoon wind field part in the analysis data by using a method combining smoothing and filtering techniques, which is as follows:
decomposing an initial wind field V into a basic field VBAnd a disturbance field VD:V=VB+VD
And then decomposing the disturbance field into a background disturbance field and a typhoon disturbance field, wherein the typhoon disturbance field is a typhoon field component to be eliminated. Disturbing typhoon in disturbing field VTyphoonRemoving to obtain environment wind field V with typhoon removedE:VE=VB+(VD-VTC)=V-VTC
Here, V denotes a vector wind field, and includes two components, i.e., a latitudinal wind u and a latitudinal wind V. The wind field in the technical route is directly represented by the latitudinal and longitudinal winds (u and v in lower case), so it is felt as if it did not correspond.
Formula VE=V-VTCIndicating a corresponding subtraction of vector components (warp and weft), i.e.
Weft wind: u. ofE=u-uTC
The meridian wind: v. ofE=v-vTC
The two can be combined together by formula VE=V-VTCDenoting, upper case V denotes the wind field vector.
To determine the location of the analytic typhoon (typhoon center), the influence of the distance can be taken into account, adding a distance weight. Calculating the disturbance wind speed V by using the following formula with the grid point closest to the observation typhoon position as the centerDijA centroid position within a range.
Figure BDA0002797232000000131
Figure BDA0002797232000000132
Figure BDA0002797232000000133
Wherein D is 200km, D is 400km,
Figure BDA0002797232000000134
respectively, the longitude and latitude, Delta S, of the center of the typhoonijIs the area corresponding to the grid or grid,
Figure BDA0002797232000000135
is the wind speed of the typhoon field, r is the distance of the analysis point from the center position of the typhoon.
The method for removing the typhoon wind field part in the analysis data by using the method of combining the smoothing and filtering technologies is an existing calculation formula, and is not repeated here.
The above-mentioned embodiments, objects, technical solutions and advantages of the present invention are further described in detail, it should be understood that the above-mentioned embodiments are only specific embodiments of the present invention, and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.

Claims (6)

1. A new method for calculating a typhoon potential generation index, the new method comprising the steps of:
the first step is as follows: obtaining the whole layer wind field structure of the atmospheric circulation background wind field from the atmospheric re-analysis data, and quantitatively separating the typhoon wind field from the atmospheric re-analysis data;
the second step is that: computing absolute vorticity eta of troposphere lower layer 850hPa by utilizing atmospheric circulation background wind fieldNOTCAnd convection layer high and low layer 200hPa minus 850hPa vertical wind shear Vshear-NOTC
Flow layer low layer 850hPa absolute vorticity etaNOTCThe calculation formula is as follows:
Figure FDA0002797231990000011
wherein
Figure FDA0002797231990000012
And
Figure FDA0002797231990000013
respectively showing the latitude wind and the longitude wind of the atmospheric circulation background after eliminating the typhoon wind field on the 850hPa level;
convection layer high and low layer 200hPa minus 850hPa vertical wind shear Vshear-NOTCThe calculation formula of (2) is as follows:
Figure FDA0002797231990000014
or
Figure FDA0002797231990000015
Wherein
Figure FDA0002797231990000016
And
Figure FDA0002797231990000017
respectively representing the atmospheric circulation background latitude wind and longitude wind after eliminating the typhoon field on the 200hPa level, wherein
Figure FDA0002797231990000018
And
Figure FDA0002797231990000019
respectively showing the latitude wind and the longitude wind of the atmospheric circulation background after eliminating the typhoon wind field on the 850hPa level;
the third step: calculation of 600hPa relative humidity H and typhoon potential intensity V using conventional datapot
The fourth step: using the parameter η obtained in the second stepNOTCAnd Vshear-NOTCAnd the parameters H and V obtained in the third steppotCalculating GPI to finally obtain an improved GPI index;
the GPI calculation formula obtained by the new method is as follows:
Figure FDA00027972319900000110
2. a new method of calculating a typhoon potential generation index according to claim 1, characterized in that in the first step: the atmospheric circulation background wind field is quantitatively separated from the atmospheric re-analysis data by a typhoon wind field separation method and is applied to each standard level of the atmospheric re-analysis data.
3. A new method of calculating a typhoon potential generation index according to claim 2, characterized in that in the first step: the method for separating the typhoon wind field comprises the steps of determining the center of the position of the typhoon in the analysis data, and then performing wind field decomposition on an area which takes the point as the center of a circle and takes a certain distance as the radius.
4. A new method of calculating a typhoon potential generation index according to claim 3, characterized in that: the radius range of the wind field decomposition is set to 800km, and reanalysis data with a resolution of 1 DEG x 1 DEG or more is used.
5. A new method of calculating a typhoon potential generation index according to claim 3, characterized in that: the principle of the wind field decomposition is that the plane non-rotational motion of the incompressible fluid has a velocity potential and a flow function, the horizontal wind field in the reanalysis data can be regarded as the superposition of the non-rotational motion, the non-scattered motion and the residual error term, the velocity fields of the non-rotational motion and the non-scattered motion are superposed to obtain a typhoon velocity field by calculating the velocity potential and the flow function of the reanalysis data wind field, and the residual error term of the reanalysis velocity field and the typhoon velocity field is the atmospheric circulation background wind field.
6. The method of claim 1, wherein the method further comprises the steps of: the whole-layer wind field structure comprises a typhoon wind field component and an atmospheric circulation background wind field component, and the typhoon wind field comprises a non-rotary motion wind field and a non-dispersion motion wind field.
CN202011336283.9A 2020-11-25 2020-11-25 Novel method for calculating typhoon potential generation index Pending CN112596127A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011336283.9A CN112596127A (en) 2020-11-25 2020-11-25 Novel method for calculating typhoon potential generation index

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011336283.9A CN112596127A (en) 2020-11-25 2020-11-25 Novel method for calculating typhoon potential generation index

Publications (1)

Publication Number Publication Date
CN112596127A true CN112596127A (en) 2021-04-02

Family

ID=75183741

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011336283.9A Pending CN112596127A (en) 2020-11-25 2020-11-25 Novel method for calculating typhoon potential generation index

Country Status (1)

Country Link
CN (1) CN112596127A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113868581A (en) * 2021-10-15 2021-12-31 江苏海洋大学 Method for judging marine phytoplankton outbreak caused by tropical cyclone

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002286863A (en) * 2001-03-28 2002-10-03 Songai Hokenryoritsu Santeikai System of evaluating wind damage due to typhoon
CN103824280A (en) * 2013-12-30 2014-05-28 杭州师范大学 Typhoon center extraction method
CN104036462A (en) * 2014-06-06 2014-09-10 南京大学 Method to estimate intensity of tropical cyclone across the Northwest Pacific Ocean by using geostationary meteorological satellite image
CN105095655A (en) * 2015-07-17 2015-11-25 南京信息工程大学 Perturbation method based on self-adaptive sensitive factors
CN109583653A (en) * 2018-11-30 2019-04-05 南京信息工程大学 The extended peroid forecasting procedure of NORTHWESTERN PACIFIC TYPHOON based on statistical model
CN109960882A (en) * 2019-03-28 2019-07-02 苏州浪潮智能科技有限公司 Forecast that tropical cyclone generates method, apparatus, equipment and the storage medium of frequency
CN110286424A (en) * 2019-04-15 2019-09-27 南京大学 Objective weather classifying method based on numerical statistic
US20200170200A1 (en) * 2018-10-02 2020-06-04 John Barton Huber System and method for modifying intensity or path of a tropical cyclone
CN111338005A (en) * 2020-03-02 2020-06-26 中国人民解放军国防科技大学 Prediction method for generation frequency of tropical cyclone in northwest Pacific ocean on monthly scale

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002286863A (en) * 2001-03-28 2002-10-03 Songai Hokenryoritsu Santeikai System of evaluating wind damage due to typhoon
CN103824280A (en) * 2013-12-30 2014-05-28 杭州师范大学 Typhoon center extraction method
CN104036462A (en) * 2014-06-06 2014-09-10 南京大学 Method to estimate intensity of tropical cyclone across the Northwest Pacific Ocean by using geostationary meteorological satellite image
CN105095655A (en) * 2015-07-17 2015-11-25 南京信息工程大学 Perturbation method based on self-adaptive sensitive factors
US20200170200A1 (en) * 2018-10-02 2020-06-04 John Barton Huber System and method for modifying intensity or path of a tropical cyclone
CN109583653A (en) * 2018-11-30 2019-04-05 南京信息工程大学 The extended peroid forecasting procedure of NORTHWESTERN PACIFIC TYPHOON based on statistical model
CN109960882A (en) * 2019-03-28 2019-07-02 苏州浪潮智能科技有限公司 Forecast that tropical cyclone generates method, apparatus, equipment and the storage medium of frequency
CN110286424A (en) * 2019-04-15 2019-09-27 南京大学 Objective weather classifying method based on numerical statistic
CN111338005A (en) * 2020-03-02 2020-06-26 中国人民解放军国防科技大学 Prediction method for generation frequency of tropical cyclone in northwest Pacific ocean on monthly scale

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
赵华睿: "西北太平洋台风环流分离方法比较", 《中国优秀硕士学位论文全文数据库 (基础科学辑)》 *
邓诗茹: "种再分析资料中热带气旋潜在生成指数分析", 《气象科学》 *
钱伊恬等: "2015年西北太平洋台风季提早展开:2015/2016超级厄尔尼诺的影响", 《大气科学学报》 *
黄小刚: "削去分析台风方法的对比研究", 《应用气象学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113868581A (en) * 2021-10-15 2021-12-31 江苏海洋大学 Method for judging marine phytoplankton outbreak caused by tropical cyclone
CN113868581B (en) * 2021-10-15 2023-08-04 江苏海洋大学 Method for judging sea surface phytoplankton outbreak caused by tropical cyclone

Similar Documents

Publication Publication Date Title
Tilinina et al. Comparing cyclone life cycle characteristics and their interannual variability in different reanalyses
Bosilovich et al. Water vapor tracers as diagnostics of the regional hydrologic cycle
Blockley et al. Recent development of the Met Office operational ocean forecasting system: an overview and assessment of the new Global FOAM forecasts
Langland et al. Initial condition sensitivity and error growth in forecasts of the 25 January 2000 East Coast snowstorm
Schenkel et al. An examination of tropical cyclone position, intensity, and intensity life cycle within atmospheric reanalysis datasets
Carrier et al. Impact of assimilating ocean velocity observations inferred from Lagrangian drifter data using the NCOM-4DVAR
CN102682335B (en) Neural network method for precisely determining tropospheric delay in region
JP2007304080A (en) Device, method and program for predicting gas condition, and diffusion state prediction system
De Mey Data Assimilation at the Oceanic Mesoscale: A Review (gtSpecial IssueltData Assimilation in Meteology and Oceanography: Theory and Practice)
Wu et al. Potential vorticity diagnosis of the factors affecting the track of Typhoon Sinlaku (2008) and the impact from dropwindsonde data during T-PARC
CN108764527B (en) Screening method for soil organic carbon library time-space dynamic prediction optimal environment variables
Beadling et al. Assessing the quality of Southern Ocean circulation in CMIP5 AOGCM and Earth system model simulations
CN114910980A (en) Tropical cyclone gale wind circle forecasting method based on subjective path strength forecasting and parameterized wind field model
CN116522600A (en) Complex wind field wind speed simulation method, device, equipment and medium
CN112596127A (en) Novel method for calculating typhoon potential generation index
Mohanty Tropical cyclones in the Bay of Bengal and deterministic methods for prediction of their trajectories
CN113326624A (en) Method and system for predicting height of atmospheric boundary layer in desert area
Huang et al. Track deflection of Typhoon Maria (2018) during a westbound passage offshore of northern Taiwan: Topographic influence
Chamberlain et al. Acoustic float tracking with the Kalman smoother
CN116609859A (en) Weather disaster high-resolution regional mode forecasting system and method
Halperin et al. Diagnosing conditions associated with large intensity forecast errors in the Hurricane Weather Research and Forecasting (HWRF) Model
Stauffer et al. A field-coherence technique for meteorological field-program design for air quality studies. Part I: Description and interpretation
Zhang et al. Spatial interpolation of air temperature with ANUSPLIN in Three Gorges Reservoir Area
Kang et al. Wind resource assessment and potential development of wind farms along the entire coast of South Korea using public data from the Korea meteorological administration
Bothwell et al. A kinematics-based GIS methodology to represent and analyze spatiotemporal patterns of precipitation change in IPCC A2 scenario

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20210402

RJ01 Rejection of invention patent application after publication