Highly resolved and accurate daily precipitation data are required for impact models to perform adequately and correctly measure the impacts of high-risk events. In order to produce such data, bias correction is often needed. Most of those statistical methods correct the probability distributions of daily precipitation by modeling them with either empirical or parametric distributions. A recent semi-parametric model based on a penalized Berk–Jones (BJ) statistical test, which allows for automatic and personalized splicing of parametric and non-parametric distributions, has been developed. This method, called the Stitch-BJ model, was found to be able to model daily precipitation correctly and showed interesting potential in a bias correction setting. In the present study, we will consolidate these results by taking into account the seasonal properties of daily precipitation in an out-of-sample context and by considering dry days probabilities in our methodology. We evaluate the performance of the Stitch-BJ method in this seasonal bias correction setting against more classical models such as the Gamma, Exponentiated Weibull (ExpW), Extended Generalized Pareto (EGP) or empirical distributions. Results show that a seasonal separation of data is necessary in order to account for intra-annual non-stationarity. Moreover, the Stitch-BJ distribution was able to consistently perform as well as or better than all the other considered models over the validation set, including the empirical distribution, which is often used due to its robustness. Finally, while methods for correcting dry day probabilities can be easily applied, their relevance can be discussed as temporal and spatial correlations are often neglected.

We are interested in creating statistical methods to provide informative summaries of random fields through the geometry of their excursion sets.To this end, we introduce an estimator for the length of the perimeter of excursion sets of random fields on $\mathbb R^2$ observed over regular square tilings. The proposed estimator acts on the empirically accessible binary digital images of the excursion regions and computes the length of a piecewise linear approximation of the excursion boundary. The estimator is shown to be consistent as the pixel size decreases, without the need of any normalization constant, and with neither assumption of Gaussianity nor isotropy imposed on the underlying random field. In this general framework, even when the domain grows to cover $\mathbb R^2$, the estimation error is shown to be of smaller order than the side length of the domain. For affine, strongly mixing random fields, this translates to a multivariate Central Limit Theorem for our estimator when multiple levels are considered simultaneously. Finally, we conduct several numerical studies to investigate statistical properties of the proposed estimator in the finite-sample data setting.