Modelling rainfall with a Bartlett–Lewis process: new developments
Abstract. The use of Poisson cluster processes to model rainfall time series at a range of scales now has a history of more than 30 years. Among them, the randomised (also called modified) Bartlett–Lewis model (RBL1) is particularly popular, while a refinement of this model was proposed recently (RBL2; Kaczmarska et al., 2014). Fitting such models essentially relies upon minimising the difference between theoretical statistics of the rainfall signal and their observed estimates. The first statistics are obtained using closed form analytical expressions for statistics of the orders 1 to 3 of the rainfall depths, as well as useful approximations of the wet–dry structure properties. The second are standard estimates of these statistics for each month of the data. This paper discusses two issues that are important for the optimal model fitting of RBL1 and RBL2. The first issue is that, when revisiting the derivation of the analytical expressions for the rainfall depth moments, it appears that the space of possible parameters is wider than has been assumed in past papers. The second issue is that care must be exerted in the way monthly statistics are estimated from the data. The impact of these two issues upon both models, in particular upon the estimation of extreme rainfall depths at hourly and sub-hourly timescales, is examined using 69 years of 5 min and 105 years of 10 min rainfall data from Bochum (Germany) and Uccle (Belgium), respectively.