Incorporation of parameter uncertainty into spatial interpolation using Bayesian trans-Gaussian kriging

Quantitative precipitation estimation (QPE) plays an important role in meteorological and hydrological applications. Ground-based telemetered rain gauges are widely used to collect precipitation measurements. Spatial interpolation methods are commonly employed to estimate precipitation fields covering non-observed locations. Kriging is a simple and popular geostatistical interpolation method, but it has two known problems: uncertainty underestimation and violation of assumptions. This paper tackles these problems and seeks an optimal spatial interpolation for QPE in order to enhance spatial interpolation through appropriately assessing prediction uncertainty and fulfilling the required assumptions. To this end, several methods are tested: transformation, detrending, multiple spatial correlation functions, and Bayesian kriging. In particular, we focus on a short-term and time-specific rather than a long-term and event-specific analysis. This paper analyzes a stratiform rain event with an embedded convection linked to the passing monsoon front on the 23 August 2012. Data from a total of 100 automatic weather stations are used, and the rainfall intensities are calculated from the difference of 15 minute accumulated rainfall observed every 1 minute. The one-hour average rainfall intensity is then calculated to minimize the measurement random error. Cross-validation is carried out for evaluating the interpolation methods at regional and local levels. As a result, transformation is found to play an important role in improving spatial interpolation and uncertainty assessment, and Bayesian methods generally outperform traditional ones in terms of the criteria.