The Netherlands, Germany, and Belgium were hit by heavy and prolonged precipitation in July 2021. As time passed, weather warnings escalated, leading to evacuations due to predicted floods, including in the Rur catchment. It was difficult to forecast the flooding of the Rur, raising the question of which elements are crucial in a flood predictionmodel for the Rur river. This question is addressed by addressing both a hindcast of the 2021 flood event and creating forecasts based on the weather forecast of July 13, 2021.

The Rur river basin is characterised by topographic and geological variations, with the steep Eifel responding differently than the flat lowlands, and human intervention in the form of reservoirs and lignite mines. A hydrological Wflow_SBM model has been derived for the Rur river basin, encompassing these characteristics, along with a hydrodynamic ProMaIDes model for the downstream reach of the Rur. These models were compared to investigate various aspects: river routing, floodplain flow, tirbutary interactions, the influence of reservoirs, and the impact of reduced groundwater levels.

The results of the 2021 floods indicate that modelling flows in floodplains is crucial to shaping the flood wave, both in tributaries and the Rur itself. Additionally, the reservoir played a significant role in attenuating the flood wave, with the increase in the outflow of the reservoir primarily affecting the tail of the wave. The reduced groundwater level was simulated by adding a leakage termto the saturated subsurface zone, whose indirect effect is significantly greater than the leakage termitself. Moreover, the tributaries Worm and Inde, particularly, are influential in the Rur’s discharge. These characteristics are also evident in the simulated forecasts, although the spatial and temporal resolution is significantly lower for these meteorological predictions.

Finally, the characteristic response of the Rur demonstrates that not everymodel type is equally practical for flood forecasting. The dominant flow from the reservoirs is highly regulated and is unlikely to induce inundations downstream. Complex flow patterns in floodplains only become relevant in the Dutch Rur, which makes two-dimensional modelling particularly valuable here. Therefore, it is recommended to use a one-dimensional discharge model, incorporating delay effects from winter bed flows. When predicted discharges at the Stah station are exceeded, two-dimensional simulations may provide a solution, the model area reduced to the Dutch Rur, focussing on predictions where a critical value related to floodplain capacity (Qlimit = 300 m^3/s) is exceeded.