Abstract. Rainfall thresholds, namely rainfall intensity-duration conditions beyond which the probability of debris flow occurrence is considered significant, can be used as a forecasting tool in debris-flow early warning system. Many uncertainties may affect the thresholds calibration and, in turn, the reliability and effectiveness of this tool. The purpose of this study is to assess the uncertainty in the determination of the rainfall threshold for stony debris flow based on the Back Dynamical Approach (BDA) (Rosatti et al., 2019), an innovative method to estimate the rainfall duration and averaged intensity strictly related to measured debris flow. The uncertainty analysis has been computed performing two Monte Carlo cascade simulations: (i) to assess the variability in the estimate of rainfall conditions due to the uncertainty of some of the BDA parameters and (ii) to quantify the impact of this variability on the threshold parameters, obtained by using the frequentist method. Then, the deviation between these analysis outcomes and the values obtained in Rosatti et al. (2019) has been examined. The results highlight that the variability in the rainfall condition estimate is strongly related to the debris flow characteristics and the hyetograph shape. Depending on these features, the spreading of the obtained distributions can take both low and high values. Instead, the threshold parameters are characterised by a low statistical spreading. Finally, the consistency between the outcome of this study and the results obtained in Rosatti et al. (2019) has been proved and the critical issues related to the rainfall condition estimation have been discussed.