Liang Li, Zhao Linna, Qi Dan, et al. The experiment of hydrologic probabilistic forecast based on the precipitation forecast calibrated by bayesian model averaging. J Appl Meteor Sci, 2013, 24(4): 416-424.
Citation: Liang Li, Zhao Linna, Qi Dan, et al. The experiment of hydrologic probabilistic forecast based on the precipitation forecast calibrated by bayesian model averaging. J Appl Meteor Sci, 2013, 24(4): 416-424.

The Experiment of Hydrologic Probabilistic Forecast Based on the Precipitation Forecast Calibrated by Bayesian Model Averaging

  • Received Date: 2012-09-04
  • Rev Recd Date: 2013-04-09
  • Publish Date: 2013-08-31
  • Based on 24-h accumulated precipitation data of the Huaihe Basin from 1 June to 31 August in 2008 and the corresponding ensemble forecast of 24 h, 48 h, 72 h from T213 model, the method of Bayesian Model Averaging (BMA) is used to calibrate quantitative precipitation forecasts of 15 members from the ensemble forecast based on the training data of 30 days. The calibrated results are verified by continuous ranked probability score (CRPS) and mean absolute error (MAE). Second, the Dapoling—Wangjiaba Catchment in the upper stream of the Huaihe River, which is subdivided into Dapoling—Xixian Catchment and Xixian-Wangjiaba Catchment, is investigated with the hydrological simulation experiment. The Xixian and Wangjiaba hydrologic stations in the upper stream of the Huaihe River are selected as representative stations. The rain process occurring from 23 June to 3 August in 2008 is investigated to simulate the runoff tendency. Then, the 25th and the 75th percentiles of 24-h, 48-h and 72-h precipitation ensemble forecast which are calibrated by BMA are used to force the Variable Infiltration Capacity (VIC) hydrological model respectively to obtain the corresponding runoff, and finally the simulate results are analyzed comparing with daily runoff observations.Results show that the precision of 24-h, 48-h, 72-h precipitation forecast of BMA model is improved after the calibration. The raw ensemble forecast of 24 h is calibrated well by BMA model. As the leading hours increase, the calibration of 48 h and 72 h is as good as that of 24 h. Although BMA calibrates on the raw ensemble forecast, the improvement of calibrated forecast depends on the accuracy of raw ensemble forecast. The valid interval given by BMA model, namely the interval from the 25th percentile to the 75th one of ensemble forecast, is more likely to contain the true value of observations according to the verification analysis of hydrological probabilistic forecast. From this aspect, the performance of BMA forecast outperform deterministic forecast. It can improve the accuracy of forecast and reduce the error by BMA, describing forecast uncertainty in the form of a probability distribution. Known from the analysis of verification index of hydrologic probabilistic forecast, it shows that the hydrological simulation forced by the calibrated precipitation is almost consistent with the runoff tendency of observations. It is effective to grasp the trend of runoff change. It indicates that the precipitation forecast calibrated by BMA can be established coupling with the VIC hydrological model as well as increasing the forecast accuracy significantly. It can meet the more and more objective, quantitative needs of decision-making service, and improve the benefit of weather forecast greatly.
  • Fig. 1  The location and coverage of Huaihe Basin (a) and the coverage of Dapoling—Wangjiaba Catchment with the illustration of interpolation (b)

    (Ⅰ represents the coverage of Dapoling—Xixian Catchment, Ⅱ represents the coverage of Xixian—Wangjiaba Catchment)

    Fig. 2  The schematic chart of hydrologic probabilistic forecast based on precipitation forecast calibrated by Bayesian Model Averaging

    Fig. 3  The comparison of Box-Whisker plot and observation plot between probabilistic forecast calibrated by Bayesian Model Averaging (a) and raw ensemble forecast (b) at Xixian Station from 21 to 30 in July 2008

    Fig. 4  The comparison of Box-Whisker plot and observation plot between probabilistic forecast calibrated by Bayesian Model Averaging (a) and raw ensemble forecast (b) at Wangjiaba Station from 21 to 30 in July 2008

    Fig. 5  The comparison between runoff simulations forcing by probabilistic forecast calibrated by Bayesian Model Averaging (BMA) and runoff observations at Xixian Station of 24 h (a), 48 h (b) and 72 h (c)

    Fig. 6  The comparison between runoff simulations forcing by probabilistic forecast calibrated by Bayesian Model Averaging (BMA) and runoff observations at Wangjiaba Station of 24 h (a), 48 h (b) and 72 h (c)

    Table  1  The comparison of CRPS and MAE scores between raw ensemble forecast (REF) and probabilistic forecast calibrated by Bayesian Model Averaging (BMA) for 24 h, 48 h and 72 h

    预报日期 评分 24 h 48 h 72 h
    REF BMA REF BMA REF BMA
    07-21 CRPS 2.05 0.62 0.71 1.21 0.86 1.52
    MAE 2.52 0.43 0.88 0.70 1.02 1.16
    07-22 CRPS 9.36 9.13 9.82 8.73 8.99 8.73
    MAE 10.84 10.02 11.02 10.96 10.74 10.71
    07-23 CRPS 31.45 45.59 37.79 47.12 40.26 48.42
    MAE 38.62 56.36 45.28 56.64 49.28 57.47
    07-24 CRPS 12.01 11.52 12.32 13.31 14.72 12.07
    MAE 14.80 14.55 14.97 14.84 18.34 17.05
    07-25 CRPS 5.11 3.16 4.25 3.31 3.52 2.96
    MAE 6.81 3.41 5.31 3.71 4.58 3.39
    07-26 CRPS 8.63 4.23 4.38 3.56 4.50 3.86
    MAE 10.78 4.96 5.38 4.28 5.74 4.70
    07-27 CRPS 6.97 4.74 5.75 4.79 5.33 4.82
    MAE 8.91 5.96 7.25 6.32 6.65 6.42
    07-28 CRPS 7.14 4.14 6.36 3.60 3.52 3.04
    MAE 9.53 4.66 8.63 3.93 4.89 3.31
    07-29 CRPS 2.82 1.18 1.56 1.08 3.20 1.63
    MAE 3.64 0.86 2.17 0.49 4.62 0.85
    07-30 CRPS 8.91 6.39 8.15 6.97 7.85 6.59
    MAE 11.22 8.85 8.78 8.64 8.94 8.35
    DownLoad: Download CSV

    Table  2  The statistics of verification index of hydrologic probabilistic forecast

    站点 预报时效/h 确定性系数 洪峰相对误差 峰现时间差/d
    第25百分位数 第75百分位数 第25百分位数 第75百分位数
    息县 24 0.56 0.70 -0.41 -0.36 0
    48 -1.13 -1.85 -0.60 -0.50 1
    72 -0.56 -2.65 -0.81 -0.70 2
    王家坝 24 0.58 0.60 0.16 0.26 -1
    48 0.71 0.82 0.04 0.17 -1
    72 0.34 0.55 0.19 0.03 0
      注:峰现时间差为正,表示预测提前;峰现时间差为负,表示预测滞后。
    DownLoad: Download CSV
  • [1]
    张亚萍, 程明虎, 徐慧, 等.雷达定量测量在佛子岭流域径流模拟中的应用.应用气象学报, 2007, 18(3):295-305. http://qikan.camscma.cn/jams/ch/reader/view_abstract.aspx?file_no=20070351&flag=1
    [2]
    王莉莉, 陈德辉, 赵琳娜.GRAPES气象-水文模式在一次洪水预报中的应用.应用气象学报, 2012, 23(3):274-284. http://qikan.camscma.cn/jams/ch/reader/view_abstract.aspx?file_no=20120303&flag=1
    [3]
    陈丽娟, 张培群, 赵振国.松嫩辽流域夏季面雨量预测因子探讨.应用气象学报, 2005, 16(5):663-669. doi:  10.11898/1001-7313.20050513
    [4]
    赵琳娜, 包红军, 田付友, 等.水文气象研究进展.气象, 2012, 38(2):147-154. doi:  10.11898/1001-7313.20120203
    [5]
    李泽椿, 毕宝贵, 朱彤, 等.近30年中国天气预报业务进展.气象, 2004, 30(12):4-10. doi:  10.3969/j.issn.1000-0526.2004.12.002
    [6]
    马培迎.应用贝叶斯原理修正降水概率预报.气象科技, 1999, 1:45-48. doi:  10.3969/j.issn.1007-9033.1999.02.015
    [7]
    Raftery A E, Gneiting T, Balabdaoui F, et al.Using Bayesian model averaging to calibrate forecast ensembles.Mon Wea Rev, 2005, 133:1155-1174. doi:  10.1175/MWR2906.1
    [8]
    陈朝平, 冯汉中, 陈静.基于贝叶斯方法的四川暴雨集合概率预报产品释用.气象, 2012, 36(5):32-39. http://www.cnki.com.cn/Article/CJFDTOTAL-QXXX201005006.htm
    [9]
    陈法敬, 矫梅燕, 陈静.一种温度集合预报产品释用方法的初步研究.气象, 2011, 37(1):14-20. doi:  10.7519/j.issn.1000-0526.2011.01.002
    [10]
    Mclean J S, Adrian E R, Rilmann G, et al.Probabilistic quantitative precipitation forecasting using Bayesian Model Averaging.MonWea Rev, 2007, 135:3209-3220.
    [11]
    李新, 程国栋, 卢玲.空间内插方法比较.地球科学进展, 2000, 15(3):260-265. http://www.cnki.com.cn/Article/CJFDTOTAL-QHYH200402003.htm
    [12]
    高歌, 龚乐冰, 赵珊珊, 等.日降水量空间插值方法研究.应用气象学报, 2007, 18(5):732-736. http://qikan.camscma.cn/jams/ch/reader/view_abstract.aspx?file_no=200705111&flag=1
    [13]
    秦涛, 付宗堂.ArcGIS中几种空间内插方法的比较.物探化探计算技术, 2007, 29(1):72-75. http://www.cnki.com.cn/Article/CJFDTOTAL-WTHT200701016.htm
    [14]
    赵琳娜, 吴昊, 田付友, 等.基于TIGGE资料的流域概率性降水预报评估.气象, 2010, 36(7):133-142. doi:  10.7519/j.issn.1000-0526.2010.07.020
    [15]
    Hamill T M, Whitaker J S, Wei X.Ensemble re-forecasting improving medium-range forecast skill using retrospective forecasts.MonWea Rev, 2004, 132:1434-1447.
    [16]
    吴洪宝, 王盘兴, 林开平.广西6、7月份若干日内最大日降水量的概率分布.热带气象学报, 2004, 20(5):586-592. http://www.cnki.com.cn/Article/CJFDTOTAL-RDQX200405014.htm
    [17]
    梁莉, 赵琳娜, 巩远发, 等.淮河流域汛期20 d内最大日降水量概率分布.应用气象学报, 2011, 22(4):421-428. doi:  10.11898/1001-7313.20110404
    [18]
    Liang Li, Zhao Linna, GongYuanfa, et al.Probability distribution of summer daily precipitation in the Huaihe Basin of China based on Gamma distribution.Acta Meteor Sinica, 2012, 26(1):72-84. doi:  10.1007/s13351-012-0107-2
    [19]
    Wu C F J.On the convergence properties of the EM algorithm.Ann Stat, 1983, 11:95-103. doi:  10.1214/aos/1176346060
    [20]
    Wood E F, Lettenmaier D P, Zartarian V G.A land-surface drology parameterization with subgrid variability for general circulation models.J Geophys Res, 1992:97.
    [21]
    Liang Xu, Xie Zhenghui.Important factors in land-atmosphere interactions:Surface runoff generations and interactions between surface and groundwater.Global Planetary Change, 2003, 38:101-114. doi:  10.1016/S0921-8181(03)00012-2
    [22]
    Zhao Linna, Qi Dan, Tian Fuyou, et al.Probabilistic flood prediction in the upper Huaihe Catchment using TIGGE data.Acta Meteor Sinica, 2012, 26(1):62-71. doi:  10.1007/s13351-012-0106-3
    [23]
    林建, 谢正辉, 陈锋, 等.2006年汛期VIC水文模型模拟结果分析.气象, 2008, 34(3):69-77. doi:  10.7519/j.issn.1000-0526.2008.03.011
    [24]
    徐虹, 朱爱华, 张宏.降水概率预报的评分和经济效益评估.陕西气象, 1999, 1:1-3. http://www.cnki.com.cn/Article/CJFDTOTAL-SXQI901.000.htm
  • 加载中
  • -->

Catalog

    Figures(6)  / Tables(2)

    Article views (5045) PDF downloads(2113) Cited by()
    • Received : 2012-09-04
    • Accepted : 2013-04-09
    • Published : 2013-08-31

    /

    DownLoad:  Full-Size Img  PowerPoint