Holistic flood risk management (FRM) aims at minimizing the flood risk (FR) for an entire river basin. To this end, the probability of flooding and its consequences must also be taken into account in the planning of flood risk management measures. If the FR is to be limited by the flooding probability value (e.g. by changing the structural reliability), knowledge of the variability of the limit state is required. As an example, this study quantifies for the first time the variability of the probability of a landside slope failure on homogeneous river levees, using the Monte Carlo simulation. In the slip circle method developed by Krey that is used for this purpose, four geometric and six material input parameters are included as being uncertain, of which four emerge in the result as being relevant influencing variables of the failure probability. These analyses extend our understanding of landward slope failure by the variability of the limit state; the results allow a targeted adjustment of the local failure probability to limit the flood risk. In the future, the methodology presented is to be supplemented by the integration of parameter dependencies and, once in use, expanded to include structured cross-sections or additional failure mechanisms.