CN109655941A - 一种基于目标系统扰动的热带气旋集合预报方法 - Google Patents
一种基于目标系统扰动的热带气旋集合预报方法 Download PDFInfo
- Publication number
- CN109655941A CN109655941A CN201811499437.9A CN201811499437A CN109655941A CN 109655941 A CN109655941 A CN 109655941A CN 201811499437 A CN201811499437 A CN 201811499437A CN 109655941 A CN109655941 A CN 109655941A
- Authority
- CN
- China
- Prior art keywords
- field
- initial
- scale
- tropical cyclone
- disturbance
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01W—METEOROLOGY
- G01W1/00—Meteorology
- G01W1/10—Devices for predicting weather conditions
Landscapes
- Environmental & Geological Engineering (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Atmospheric Sciences (AREA)
- Biodiversity & Conservation Biology (AREA)
- Ecology (AREA)
- Environmental Sciences (AREA)
- Cyclones (AREA)
Abstract
本发明提供一种基于目标系统扰动的热带气旋集合预报方法,包括下列步骤:1)确定模式初始背景场;2)对初始背景场进行尺度分离,得到大尺度背景场和中小尺度扰动场;3)对分离出的小尺度场进行扰动;4)将分离出的大尺度场与扰动后的中小尺度场进行叠加,得到基于目标系统扰动的热带气旋集合预报初始成员;5)将步骤4)得到的初始集合成员作为初始场,开展热带气旋集合预报。本发明通过将背景场分解为缓变的大尺度场和快变的中小尺度场,并对中小尺度场进行扰动得到基于目标系统扰动的热带气旋集合预报初始成员,解决了目前热带气旋集合预报中多采用全球模式中长期预报扰动技术的问题,使初始集合成员更有针对性。
Description
技术领域
本发明属于数值天气预报技术领域,尤其涉及一种基于目标系统扰动的热带气旋集合预报方法。
背景技术
对热带气旋数值预报而言,目前形成的共识是由于初值误差和模式误差的存在,单一初值单一模式的确定性预报难以做出稳定准确的预报结果。为了减小各种误差对预报结果的影响,集合预报的方法被应用于热带气旋的数值预报研究和业务预报中来。目前研究和业务中常见的热带气旋集合预报方法可分为三类,即扰动初始环境场、扰动热带气旋初始位置和强度、扰动模式。目前热带气旋的集合数值预报尚处在发展阶段,还存在一些亟待解决的问题,包括:(1)目前使用较多的初值扰动方法没有针对性,如BGM、SV方法等,主要用于中高纬度地区的中期天气预报,而低纬地区的天气系统,特别是热带气旋等,在动力和热力机制上与中高纬地区天气系统有很大差别,导致被应用到热带气旋的集合数值预报中效果大打折扣。(2)热带气旋路径数值预报对积云参数化方案敏感性不强,不同参数化方案之间的差异,不足以对热带台风内部对流结构造成决定性的影响,这导致基于物理过程扰动的集合预报发散度不够。(3)多模式多初值的超级集合预报方法在实际工作中还不便开展,对大多数业务机构来说,能够同时获取的初始场和数值模式是有限的。
对热带气旋的路径而言可以分为两类,一类是稳定的正常路径,其余是突然转向的异常路径,目前数值模式的发展水平对正常路径的热带气旋已经具有很好的预报能力,而应重点针对异常路径热带气旋的产生原因开展专门的设计,进而设计有针对性的热带气旋集合预报技术。从宏观角度来说,制约热带气旋运动过程的科学问题是两个方面:一是热带气旋系统自身的结构,包括大小、强度等;另一方面是环境场的引导作用。对环境场的扰动,如前所述主要是采用MCF、LAF、BGM、SV等方法,这些方法适合全球模式,对热带气旋的预报来说缺乏针对性。台风由热带地区生成往中纬度运动过程中,主要受到中纬度环流系统的影响,因此相较于对整个初始场进行扰动,只对中纬度相关环流系统进行有针对性的扰动,可能会得到更合理的集合预报初始成员。
中纬度的环流系统可以分解为稳定的大尺度背景场和快变的中小尺度扰动场,而热带气旋路径的突变往往和热带气旋周边中小尺度扰动场的变化有关。因此可将热带气旋周边的中小尺度扰动场看作对热带气旋路径突变有重要影响的目标系统,通过对其进行分离合扰动,从而得到有针对性的基于目标系统扰动的热带气旋集合预报方法。
发明内容
本发明的目的是提供一种基于目标系统扰动的热带气旋集合预报方法,用以解决目前还缺乏专门针对热带气旋数值预报的集合成员构造方法的问题。
为实现上述目的,本发明采用以下技术方案:
一种基于目标系统扰动的热带气旋集合预报方法,包括下列步骤:
1)根据需要确定数值初始背景场,该背景场在水平方向为经纬网格,其水平分辨率为Δφ;
2)对初始背景场进行尺度分离,得到大尺度背景场和中小尺度扰动场。假设某要素场F=FB+FD,其中FB为大尺度背景场,FD为中小尺度扰动场,为了对初始背景场进行尺度分离,分别对初始场F在纬向和经向进行滤波,其滤波公式为:
其中,下标λ、分别为经向和纬向的某一格点,n=1,2,3,......,M,m=2,3,4,2,5,6,2,8,9,2。
3)利用步骤2)得到中小尺度扰动场FD=F-FB,将FD看作目标系统并对其进行扰动,其扰动方法为将目标系统乘以一系数s,该系数取值范围设定为-1≤s≤1。
4)将分离出的大尺度场与扰动后的中小尺度场进行叠加,得到基于目标系统扰动的热带气旋集合预报初始成员。其中某一集合成员计算式为FE=F+s×FD。
5)假设通过步骤4)得到N个初始集合成员,则将这N个初始集合成员分别作为模式初始场,通过台风初始化后即可开展热带气旋路径和强度集合预报。
本发明具有如下有益效果:本发明提供了一种客观有效的方法,通过背景场尺度分离得到对热带气旋异常路径影响显著的中小尺度系统,即目标系统。并通过对该目标系统的合理扰动得到初始集合成员,解决了目前在初值扰动方法中多采用适用于全球模式中长期预报的集合成员构造方法,而没有专门针对热带气旋初值扰动的初始集合成员构造方法的问题。
应当理解,前述构思以及在下面更加详细地描述的额外构思的所有组合只要在这样的构思不相互矛盾的情况下都可以被视为本公开的发明主题的一部分。另外,所要求保护的主题的所有组合都被视为本公开的发明主题的一部分。
结合附图从下面的描述中可以更加全面地理解本发明教导的前述和其他方面、实施例和特征。本发明的其他附加方面例如示例性实施方式的特征和/或有益效果将在下面的描述中显见,或通过根据本发明教导的具体实施方式的实践中得知。
附图说明
附图不意在按比例绘制。在附图中,在各个图中示出的每个相同或近似相同的组成部分可以用相同的标号表示。为了清晰起见,在每个图中,并非每个组成部分均被标记。现在,将通过例子并参考附图来描述本发明的各个方面的实施例,其中:
图1是基于目标系统扰动的热带气旋集合预报方法流程图;
图2-1~图2-3是背景场尺度分离实例,该要素场为海平面气压场;其中:
图2-1为原始场背景;图2-2为分离出的大尺度环境场;图2-3为中小尺度扰动场;
图3是对中小尺度扰动场进行扰动后,500hPa高度场上5880和5820位势米等高线分布图;
图4是以0920号台风为例利用本发明提供的基于目标系统扰动的集合预报方法开展热带气旋路径集合预报的结果。
具体实施方式
为了更了解本发明的技术内容,特举具体实施例并配合所附图式说明如下。
在本公开中参照附图来描述本发明的各方面,附图中示出了许多说明的实施例。本公开的实施例不必定意在包括本发明的所有方面。应当理解,上面介绍的多种构思和实施例,以及下面更加详细地描述的那些构思和实施方式可以以很多方式中任意一种来实施,这是因为本发明所公开的构思和实施例并不限于任何实施方式。另外,本发明公开的一些方面可以单独使用,或者与本发明公开的其他方面的任何适当组合来使用。
本发明提供的基于目标系统扰动的热带气旋集合预报方法主要包括下列步骤:
1)根据需要确定数值初始背景场,该背景场在水平方向为经纬网格,背景场的水平分辨率为Δφ;
2)对采用滤波方法对初始背景场进行尺度分离,得到大尺度背景场和中小尺度扰动场;
方法为:假设某要素场F=FB+FD,其中FB为大尺度背景场,FD为中小尺度扰动场,为了对初始背景场进行尺度分离,分别对初始场F在纬向和经向进行滤波,其滤波方法为:
其中,下标λ、分别为经向和纬向的某一格点,n=1,2,3,......,M,m=2,3,4,2,5,6,2,8,9,2;
3)利用步骤2)得到中小尺度扰动场FD=F-FB,将FD看作目标系统并对其进行扰动,其扰动方法为将目标系统乘以扰动系数s,该系数取值范围设定为-1≤s≤1;
4)将分离出的大尺度场与扰动后的中小尺度场进行叠加,得到基于目标系统扰动的热带气旋集合预报初始成员,某一集合预报初始成员FE=F+s×FD;
5)设通过步骤4)得到N个初始集合成员,则将这N个初始集合成员分别作为模式初始场,通过台风初始化后即可开展热带气旋路径和强度集合预报。
本例中,以0920号台风路径集合预报为例:
步骤1)以0920号台风路径集合预报为例,其预报背景场采用美国环境预报中心NCEP再分析资料,其水平分辨率为1°。
步骤2)对初始场采用步骤2)提供的公式开展尺度分离,图2-1~2-3为2009年10月20日00时(世界时)海平面气压场尺度分离实例,图2-1为原始场背景;图2-2为分离出的大尺度环境场;图2-3为中小尺度扰动场。其中Δφ=1、M=1。
步骤3)将分离出的中小尺度场看作是目标系统,并对其进行扰动。其中扰动系数s取值为-1.0、-0.9、-0.8、-0.7、-0.6、-0.5、-0.4、-0.3、-0.2、-0.1、0.0、0.1、0.2、0.3、0.4、0.5、0.6、0.7、0.8、0.9、1.0。图3是对目标系统进行扰动后,2009年10月20日00时(世界时)500hPa高度场上5880和5820位势米等高线分布图;
步骤4)将扰动场与大尺度背景场叠加,得到初始集合成员。
步骤5)将步骤4)得到的初始集合成员作为初始场开展热带气旋路径和强度预报,图4为采用基于目标系统扰动的热带气旋集合预报方法得到的0920号台风路径集合预报结果。
虽然本发明已以较佳实施例揭露如上,然其并非用以限定本发明。本发明所属技术领域中具有通常知识者,在不脱离本发明的精神和范围内,当可作各种的更动与润饰。因此,本发明的保护范围当视权利要求书所界定者为准。
Claims (1)
1.一种基于目标系统扰动的热带气旋集合预报方法,其特征在于,包括下列步骤:
1)根据需要确定数值初始背景场,该背景场在水平方向为经纬网格,背景场的水平分辨率为Δφ;
2)对采用滤波方法对初始背景场进行尺度分离,得到大尺度背景场和中小尺度扰动场;
方法为:假设某要素场F=FB+FD,其中FB为大尺度背景场,FD为中小尺度扰动场,为了对初始背景场进行尺度分离,分别对初始场F在纬向和经向进行滤波,其滤波方法为:
其中,下标λ、分别为经向和纬向的某一格点,
3)利用步骤2)得到中小尺度扰动场FD=F-FB,将FD看作目标系统并对其进行扰动,其扰动方法为将目标系统乘以扰动系数s,该系数取值范围设定为-1≤s≤1;
4)将分离出的大尺度场与扰动后的中小尺度场进行叠加,得到基于目标系统扰动的热带气旋集合预报初始成员,某一集合预报初始成员FE=F+s×FD;
5)设通过步骤4)得到N个初始集合成员,则将这N个初始集合成员分别作为模式初始场,通过台风初始化后即可开展热带气旋路径和强度集合预报。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811499437.9A CN109655941A (zh) | 2018-12-09 | 2018-12-09 | 一种基于目标系统扰动的热带气旋集合预报方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811499437.9A CN109655941A (zh) | 2018-12-09 | 2018-12-09 | 一种基于目标系统扰动的热带气旋集合预报方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109655941A true CN109655941A (zh) | 2019-04-19 |
Family
ID=66113037
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811499437.9A Pending CN109655941A (zh) | 2018-12-09 | 2018-12-09 | 一种基于目标系统扰动的热带气旋集合预报方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109655941A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110443993A (zh) * | 2019-06-12 | 2019-11-12 | 中国科学院大气物理研究所 | 一种适用于模式预报enso的方法 |
CN111353217A (zh) * | 2020-02-19 | 2020-06-30 | 中山大学 | 台风路径预测的方法、设备、计算机设备和存储介质 |
CN113012285A (zh) * | 2021-02-01 | 2021-06-22 | 中国气象局广州热带海洋气象研究所(广东省气象科学研究所) | 一种针对南海台风分尺度混合的涡旋初始场构造方法 |
CN113485995B (zh) * | 2021-07-23 | 2024-03-22 | 中国科学院大气物理研究所 | 一种集合预报系统初值的扰动方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20160060793A (ko) * | 2014-10-21 | 2016-05-31 | 군산대학교산학협력단 | 능동형 정보수집 스크립트 및 수치모형을 이용한 태풍해일 자동 예측 방법 |
CN107238873A (zh) * | 2017-05-05 | 2017-10-10 | 池德龙 | 一种关于气旋状态的统计分析预报系统及调节方法 |
CN108647806A (zh) * | 2018-04-04 | 2018-10-12 | 南京信息工程大学 | 一种考虑大尺度环流背景场的延伸期预报方法 |
-
2018
- 2018-12-09 CN CN201811499437.9A patent/CN109655941A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20160060793A (ko) * | 2014-10-21 | 2016-05-31 | 군산대학교산학협력단 | 능동형 정보수집 스크립트 및 수치모형을 이용한 태풍해일 자동 예측 방법 |
CN107238873A (zh) * | 2017-05-05 | 2017-10-10 | 池德龙 | 一种关于气旋状态的统计分析预报系统及调节方法 |
CN108647806A (zh) * | 2018-04-04 | 2018-10-12 | 南京信息工程大学 | 一种考虑大尺度环流背景场的延伸期预报方法 |
Non-Patent Citations (2)
Title |
---|
熊春晖 等: "《集合-变分数据同化方法的发展和应用》", 《地球科学进展》 * |
谭燕 等: "《多种扰动组合的热带气旋路径集合预报技术研究》", 《高原气象》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110443993A (zh) * | 2019-06-12 | 2019-11-12 | 中国科学院大气物理研究所 | 一种适用于模式预报enso的方法 |
CN111353217A (zh) * | 2020-02-19 | 2020-06-30 | 中山大学 | 台风路径预测的方法、设备、计算机设备和存储介质 |
CN113012285A (zh) * | 2021-02-01 | 2021-06-22 | 中国气象局广州热带海洋气象研究所(广东省气象科学研究所) | 一种针对南海台风分尺度混合的涡旋初始场构造方法 |
CN113012285B (zh) * | 2021-02-01 | 2023-12-01 | 中国气象局广州热带海洋气象研究所(广东省气象科学研究所) | 一种针对南海台风分尺度混合的涡旋初始场构造方法 |
CN113485995B (zh) * | 2021-07-23 | 2024-03-22 | 中国科学院大气物理研究所 | 一种集合预报系统初值的扰动方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109655941A (zh) | 一种基于目标系统扰动的热带气旋集合预报方法 | |
Jaiser et al. | Impact of sea ice cover changes on the Northern Hemisphere atmospheric winter circulation | |
Alamaniotis et al. | Evolutionary multiobjective optimization of kernel-based very-short-term load forecasting | |
Martínez‐Alvarado et al. | The dichotomous structure of the warm conveyor belt | |
Jia et al. | Evaluation of machine learning models for predicting daily global and diffuse solar radiation under different weather/pollution conditions | |
Jouini et al. | The Sicily Channel surface circulation revisited using a neural clustering analysis of a high‐resolution simulation | |
Meidani et al. | Long-lead streamflow forecasting in the southwest of Iran by sea surface temperature of the Mediterranean Sea | |
Gangopadhyay et al. | Feature-oriented regional modeling and simulations in the Gulf of Maine and Georges Bank | |
Dai et al. | The outflow–rainband relationship induced by environmental flow around tropical cyclones | |
CN113779760A (zh) | 基于可预测气候模态的动力-统计结合季节气候预测方法 | |
Yoon et al. | Combinational optimization of the WRF physical parameterization schemes to improve numerical sea breeze prediction using micro-genetic algorithm | |
Ling et al. | Different effects of tropical cyclones generated in the South China Sea and the northwest Pacific on the summer South China Sea circulation | |
Durán-Grados et al. | The influence of emissions from maritime transport on air quality in the strait of gibraltar (Spain) | |
Zhang et al. | Seasonal velocity variations over the entire Kuroshio path part I: Data analysis and numerical experiments | |
Wu et al. | Do environmental regulation and foreign direct investment drive regional air pollution in China? | |
Leps et al. | Sensitivity of limited area atmospheric simulations to lateral boundary conditions in idealized experiments | |
Yin et al. | Kick risk forecasting and evaluating during drilling based on autoregressive integrated moving average model | |
Sun et al. | The role of planetary waves in the downward influence of stratospheric final warming events | |
Long et al. | Scenario changes in the climatology of winter midlatitude cyclone activity over eastern North America and the Northwest Atlantic | |
CN114021445B (zh) | 一种基于随机森林模型的海洋涡旋混合非局地性预测方法 | |
Yang et al. | Non-linear interactions between tides and storm surges during extreme weather events over the eastern Canadian shelf | |
Kriebel et al. | Coastal upwelling prediction with a mixture of neural networks | |
Tsegas et al. | A metamodelling implementation of a two-way coupled mesoscale–microscale flow model for urban area simulations | |
Yamashiki et al. | Sensitivity Analysis on Lake Biwa under the A1B SRES climate change scenario using Biwa-3D Integrated Assessment Model: part I–projection of lake temperature | |
Horton et al. | Automated input variable selection for analog methods using genetic algorithms |
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: 20190419 |
|
RJ01 | Rejection of invention patent application after publication |