Bayesian typhoon precipitation prediction with a mixture of ensemble forecast-based and historical event-based prediction functions

Abstract

Statistical post-processing has been widely adopted in probabilistic forecasting, especially when a predictive distribution of weather variables is of interest. Among the available post-processing methods, a common tool is the ensemble Bayesian model averaging. Challenges arise when it is applied to heavy rain prediction, including the precipitation produced by tropical cyclones. The issues are related to the rare extreme event and limited available data for model training. To address this, we propose in this research a Bayesian mixture model, Bmix, based on information from forecast members and historical tropical cyclones, rather than solely on immediately observed precipitation data. This approach provides a grid-specific Bayesian predictive distribution and an easy-to-derive probabilistic forecast with posterior samples generated from the Markov chain Monte Carlo algorithm. This post-processing procedure can be executed on each grid simultaneously, reducing the computational time significantly. In addition, a categorized Bmix (C.Bmix) is implemented to accommodate different scenarios such as heavy or torrential rain in applications. The proposed approaches are demonstrated on the 24-hour accumulated precipitation forecast of a category 4-equivalent super typhoon Dujuan landing in Taiwan in 2015; while two others, Matmo and Soudelor, are adopted to provide historical information. Over the Taiwan main island, a total of 1282 grids are considered, each with a resolution of 5 km and 20 ensemble member forecasts from the Weather Research and Forecasting Ensemble Prediction System (WEPS) at a 6-hour interval during the typhoon period. The analyses show that the proposed mixture approaches present larger percentages of positive CRPSS values than traditional BMA. For instance, at the time when Dujuan made landfall, the percentages are 92.2 % for Bmix, 95.2 % for C.Bmix, and 15.4 % for BMA when the true precipitation amount exceeds 200 mm/24 h. The proposed methods outperform in the probability prediction and the overall pattern of precipitation.