Jiang Wenjing, Liang Xudong. Application of PFI-4DVar data assimilation technique to nowcasting of numerical model. J Appl Meteor Sci, 2020, 31(5): 543-555. DOI:  10.11898/1001-7313.20200503.
Citation: Jiang Wenjing, Liang Xudong. Application of PFI-4DVar data assimilation technique to nowcasting of numerical model. J Appl Meteor Sci, 2020, 31(5): 543-555. DOI:  10.11898/1001-7313.20200503.

Application of PFI-4DVar Data Assimilation Technique to Nowcasting of Numerical Model

DOI: 10.11898/1001-7313.20200503
  • Received Date: 2020-06-16
  • Rev Recd Date: 2020-07-30
  • Publish Date: 2020-09-30
  • Nowcasting is mainly based on radar echo or satellite image extrapolation method. However, the prediction ability of extrapolation method decreases with time, because this method cannot describe the physical mechanism during the occurrence, development and extinction of severe convective weather systems. Considering that the prediction ability of numerical model improves with time, the nowcasting system should be based on numerical forecast model. And appropriate data assimilation technology can be used to produce a more accurate initial field, making the integral forecast results closer to the reality. The PFI-4DVar assimilation method (four-dimensional variational technology under physical filter initialization) can filter in the process of assimilation rather than model integration, thus shortening the model spin-up time and getting a more dynamic and physically coordinated analysis field. Therefore, PFI-4DVar assimilation method not only improves model prediction results, but also makes initial field closer to observations, which is very suitable for nowcasting.Using WRF model and WRFDA assimilation system, effects of PFI-4DVar on prediction ability of numerical nowcasting are explored. Through the precipitation case in North China on 11 August 2018, prediction results in control and assimilation tests are discussed. According to ETS scores, the precipitation prediction of assimilation test is closer to the observation compared with control test. The water vapor in assimilated ground and sounding data, the dynamic field in assimilated radar radial wind data and the appropriate cumulus parameterization scheme make the amplitude of divergence in high-level and convergence in low-level in analysis field of assimilation test much stronger than those in background field, thus creating vertical motion. Moreover, the precipitation of assimilation test is mainly caused by process of cumulus.A batch test is carried out on 17 precipitation cases of North China in August 2018. It shows that PFI-4DVar can significantly improve the prediction ability for short precipitation (especially large order precipitation) and timely predict the fall area of heavy rain or rainstorm. After assimilation, ETS scores of 6-hour accumulated precipitation (greater than 25.0 mm) in batch test increase from 0.125 to 0.190, and ETS scores of 6-hour accumulated precipitation (greater than 60.0 mm) increase from 0.016 to 0.081. PFI-4DVar significantly improves the precipitation nowcasting. Calculations are reduced by selecting 12-minute assimilation time window, which greatly saves computational resources. And the time of assimilation test is shortened, ensuring the time efficiency of 6-hour forecast. Therefore, PFI-4DVar can improve and enhance the prediction ability of precipitation nowcasting.
  • Fig. 1  Observed six-hour accumulated precipitation during 1200-1800 UTC on 11 Aug 2018

    Fig. 2  Six-hour accumulated precipitation in observation, CTL experiment(without data assimilation), and PFI experiment(with data assimilation) on 11 Aug 2018

    Fig. 3  One-hour accumulated precipitation in observation, CTL(control experiment without data assimilation) and PFI(data assimilation experiment) during 1200-1300 UTC, 1300-1400 UTC and 1400-1500 UTC on 11 Aug 2018

    Fig. 4  Section of relative humidity in background and analysis fields on 11 Aug 2018

    Fig. 5  Winds in the background, analysis fields and analysis increments at 850 hPa, 700 hPa and 500 hPa on 11 Aug 2018

    Fig. 6  The same as in Fig. 5, but for divergence

    Fig. 7  Averaged convective parameterization(RAINC) and averaged grid-resolvable(RAINNC) precipitation accumulated within one-hour from the 1st to the 3rd hour in observation, CTL(without data assimilation), and PFI experiment(with data assimilation) of seventeen cases in Aug 2018

    Table  1  Seventeen selected precipitation cases in North China in Aug 2018

    日期 起始时刻 6 h最大降水量/mm 有降水产生的站点总数 6 h累积降水量不大于60 mm的站点数
    2018-08-03 00:00 116.6 3536 67
    2018-08-03 12:00 131.3 3837 37
    2018-08-07 00:00 66.4 3642 2
    2018-08-08 00:00 137.0 3825 65
    2018-08-09 00:00 61.5 2479 1
    2018-08-11 12:00 137.3 3522 71
    2018-08-12 12:00 166.8 6616 112
    2018-08-13 00:00 115.1 7608 22
    2018-08-13 12:00 191.5 6624 268
    2018-08-14 00:00 172.0 4614 194
    2018-08-14 12:00 255.7 2572 67
    2018-08-16 00:00 97.0 5095 12
    2018-08-17 12:00 231.3 7707 219
    2018-08-18 00:00 258.0 7528 315
    2018-08-19 00:00 283.8 6946 174
    2018-08-19 12:00 161.0 6097 177
    2018-08-30 00:00 108.5 4473 14
    DownLoad: Download CSV

    Table  2  ETS scores of one-hour accumulated precipitation forecasts of seventeen cases in Aug 2018

    预报时间 试验 ETS评分
    0.1 mm 1.5 mm 7.0 mm 13.0 mm 40.0 mm
    第1小时 CTL 0.204 0.087 0.008 0 0
    3DV 0.203 0.081 0.007 0 0
    PFI 0.183 0.136 0.073 0.027 0
    第2小时 CTL 0.240 0.215 0.085 0.007 0
    3DV 0.242 0.201 0.082 0.009 0
    PFI 0.210 0.211 0.133 0.076 0.010
    第3小时 CTL 0.216 0.195 0.056 0.004 0
    3DV 0.221 0.188 0.069 0.009 0
    PFI 0.206 0.215 0.131 0.039 0.010
    DownLoad: Download CSV

    Table  3  ETS scores of six-hour accumulated precipitation forecasts of seventeen cases in Aug 2018

    试验 ETS评分
    0.1 mm 4.0 mm 13.0 mm 25.0 mm 60.0 mm
    CTL 0.202 0.245 0.187 0.125 0.016
    3DV 0.202 0.222 0.163 0.108 0.017
    PFI 0.181 0.241 0.210 0.190 0.081
    DownLoad: Download CSV

    Table  4  ETS scores of one-hour accumulated precipitation forecasts from the 1st to the 3rd hour on 11 Aug 2018

    预报时间 ETS评分
    0.1 mm 1.5 mm 7.0 mm
    CTL PFI CTL PFI CTL PFI
    第1小时 -0.025 0.049 -0.004 0.008 0.0 0.0
    第2小时 -0.024 0.024 -0.038 0.015 -0.015 0.014
    第3小时 0.002 0.089 -0.024 0.053 -0.014 0.042
    DownLoad: Download CSV
  • [1]
    Ligda M G.Horizontal Motion of Small Precipitation Areas as Observed by 443 Radar.Tech Rep 21, Department of Meteorology, MIT, 1953:60.
    [2]
    Kessler E.Computer program for calculating average lengths of weather radar echoes and pattern bandedness.J Atmos Sci, 1966, 23:569-574. doi:  10.1175/1520-0469(1966)023<0569:CPFCAL>2.0.CO;2
    [3]
    Barclay P A, Wilk K E.Severe thunderstorm radar echo motion and related weather events hazardous to aviation operations.Essa Tech Mem, 1970, 46:63. http://cn.bing.com/academic/profile?id=2f05b99ce7d381c9cf690160563e76ae&encoded=0&v=paper_preview&mkt=zh-cn
    [4]
    陈明轩, 俞小鼎, 谭晓光, 等.对流天气临近预报技术的发展与研究进展.应用气象学报, 2004, 15(6):754-766. http://qikan.camscma.cn/article/id/20040693
    [5]
    Rinehart R E, Garvey E T.Three-dimensional storm motion detection by conventional weather radar.Nature, 1978, 273:287-289. doi:  10.1038/273287a0
    [6]
    Rinehart R E.A pattern recognition technique for use with conventional weather radar to determine internal storm motions.J Atmos Technol, 1981, 13:119-134.
    [7]
    王改利, 刘黎平, 阮征.多普勒雷达资料在暴雨临近预报中的应用.应用气象学报, 2007, 18(3):388-395;417. http://qikan.camscma.cn/article/id/20070363
    [8]
    陈明轩, 王迎春, 俞小鼎.交叉相关外推算法的改进及其在对流临近预报中的应用.应用气象学报, 2007, 18(5):690-701. http://qikan.camscma.cn/article/id/200705105
    [9]
    Browning K A, Collier C G.An integrated radar-satellite nowcasting system in the UK//Nowcasting.New York: Academic Press, 1982: 47-61.
    [10]
    王丽荣, 卞韬, 苏运涛, 等.晴空回波在强对流天气临近预报中的应用.应用气象学报, 2010, 21(5):606-613. http://qikan.camscma.cn/article/id/20100510
    [11]
    Dixon M, Wiener G.TITTAN:Thunderstorm identification, tracking, analysis and nowcasting-A radar-based methodology.J Atmos Oceanic Technol, 1993, 10:785-797. doi:  10.1175/1520-0426(1993)010<0785:TTITAA>2.0.CO;2
    [12]
    Wilson J W, Crook N A, MuellerC K, et al.Nowcasting thunderstorms:A status report.Bull Amer Meteor Soc, 1998, 79:2079-2100. doi:  10.1175/1520-0477(1998)079<2079:NTASR>2.0.CO;2
    [13]
    Wolfson M M, Forman G E, HallowellR G, et al.Consolidated Storm Prediction for Aviation (CoSPA)//13th Conference on Aviation, Range and Aerospace Meteorology.New Orleans, LA, Amer Meteor Soc, 2008: 58-62.
    [14]
    Benjamin S G, Hu M, Weygandt S, et al.Rapid Updating NWP: Integrated Assimilation of Radar/Sat/METAR Cloud Data for Initial Hydrometer/Divergence to Improve Hourly Updated Short-range Forecasts from RUC/RR/HRRR.WMO Symposium on Nowcasting, Whistler, Canada, 2009.
    [15]
    Benjamin S G, Weygandt S S, Brown J M.A North American hourly assimilation and model forecast cycle-The rapid refresh.Mon Wea Rev, 2016, 144:1669-1694. doi:  10.1175/MWR-D-15-0242.1
    [16]
    俞小鼎, 王秀明, 李万莉, 等.雷暴与强对流临近预报.北京:气象出版社, 2020.
    [17]
    冯佳宁, 端义宏, 徐晶, 等.雷达资料同化对2015年台风彩虹数值模拟改进.应用气象学报, 2017, 28(4):399-413. doi:  10.11898/1001-7313.20170402
    [18]
    余贞寿, 冀春晓, 杨程, 等.同化风廓线雷达资料对浙江降水预报改进评估.应用气象学报, 2018, 29(1):97-110. doi:  10.11898/1001-7313.20180109
    [19]
    Abhilash S, Sahai A K, Mohankumar K, et al.Assimilation of doppler weather radar radial velocity and reflectivity observations in WRF-3DVAR system for short-range forecasting of convective storms.Pure & Applied Geophysics, 2012, 169(11):2047-2070. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=1036ae8d0a7accd7bd52091eac2cf33e
    [20]
    Takano I, Segami A.Assimilation and initialization of a mesoscale model for improved spin-up of precipitation.Analytical Chemistry, 2011, 83(3):5145-5152. http://cn.bing.com/academic/profile?id=4fbd2a8b99a7d2f2b3a0cb66c94e3e1c&encoded=0&v=paper_preview&mkt=zh-cn
    [21]
    Lynch P, Huang X Y.Initialization of the HIRLAM model using a digital filter.Mon Wea Rev, 1992, 120:1019-1034. doi:  10.1175/1520-0493(1992)120<1019:IOTHMU>2.0.CO;2
    [22]
    Polavarapu S, Tanguay M, Fillion L.Four-dimensional variational data assimilation with digital filter initialization.Mon Wea Rev, 2000, 128:2491-2510. doi:  10.1175/1520-0493(2000)128<2491:FDVDAW>2.0.CO;2
    [23]
    Wee T K, Kuo Y H.Impact of a digital filter as a weak constraint in MM54DVAR:An observing system simulation experiment.Mon Wea Rev, 2004, 132:543-559. doi:  10.1175/1520-0493(2004)132<0543:IOADFA>2.0.CO;2
    [24]
    Sun J, Crook N A.Real-time low-level wind and temperature analysis using single WSR-88D data.Wea Forecasting, 2001, 16:117-132. doi:  10.1175/1520-0434(2001)016<0117:RTLLWA>2.0.CO;2
    [25]
    Kleist D T, Parrish D F, Derber J C, et al.Improving incremental balance in the GSI 3DVAR analysis system.Mon Wea Rev, 2009, 137:1046-1060. doi:  10.1175/2008MWR2623.1
    [26]
    Vendrasco E P, Sun J, Herdies D L, et al.Constraining a 3DVAR radar data assimilation system with large-scale analysis to improve short-range precipitation forecasts.J Climate Appl Meteor, 2016, 55:673-690. doi:  10.1175/JAMC-D-15-0010.1
    [27]
    Liang X, Wang B, Chan J C L, et al.Tropical cyclone forecasting with model-constrained 3D-Var Ⅰ:Description.Quart J Roy Meteor Soc, 2007, 133:147-153. doi:  10.1002/qj.9
    [28]
    黄伟, 梁旭东.台风涡旋循环初始化方法及其在GRAPES-TCM中的应用.气象学报, 2010, 68(3):365-375. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=qxxb201003008
    [29]
    Peng W, Liang X D, Zhang X, et al.Application of physical filter initialization in 4DVar.Mon Wea Rev, 2017, 145:2201-2216. doi:  10.1175/MWR-D-16-0274.1
    [30]
    Sun J, Wang H, Tong W, et al.Comparison of the impacts of momentum control variables on high-resolution variational data assimilation and precipitation forecasting.Mon Wea Rev, 2016, 144:149-169. doi:  10.1175/MWR-D-14-00205.1
    [31]
    Parrish D F, Derber J C.The National Meteorological Center's spectral statistical-interpolation analysis system.Mon Wea Rev, 1992, 120:1747-1763. doi:  10.1175/1520-0493(1992)120<1747:TNMCSS>2.0.CO;2
    [32]
    王曼, 李华宏, 段旭, 等.WRF模式三维变分中背景误差协方差估计.应用气象学报, 2011, 22(4):482-492. http://qikan.camscma.cn/article/id/20110411
    [33]
    Liang X D, Xie Y X, Yin J F, et al.An IVAP based dealiasing method for radar velocity data quality control.J Atmos Oceanic Technol, 2019, 36:2069-2085. doi:  10.1175/JTECH-D-18-0216.1
    [34]
    Donner L J.An initialization for cumulus convection in numerical weather prediction models.Mon Wea Rev, 1988, 116(2):377-385. doi:  10.1175/1520-0493(1988)116<0377:AIFCCI>2.0.CO;2
  • 加载中
  • -->

Catalog

    Figures(7)  / Tables(4)

    Article views (2130) PDF downloads(122) Cited by()
    • Received : 2020-06-16
    • Accepted : 2020-07-30
    • Published : 2020-09-30

    /

    DownLoad:  Full-Size Img  PowerPoint