AbstractThe German Bight was exposed to record high riverine discharges in June 2013, as a result of flooding of the Elbe and Weser rivers. Several anomalous observations suggested that the hydrodynamical and biogeochemical states of the system were impacted by this event. In this study, we developed a biogeochemical model and coupled it with a previously introduced high-resolution hydrodynamical model of the southern North Sea in order to better characterize these impacts and gain insight into the underlying processes. Performance of the model was assessed using an extensive set of in situ measurements for the period 2011–2014. We first improved the realism of the hydrodynamic model with regard to the representation of cross-shore gradients, mainly through inclusion of flow-dependent horizontal mixing. Among other characteristic features of the system, the coupled model system can reproduce the low salinities, high nutrient concentrations and low oxygen concentrations in the bottom layers observed within the German Bight following the flood event. Through a scenario analysis, we examined the sensitivity of the patterns observed during July 2013 to the hydrological and meteorological forcing in isolation. Within the region of freshwater influence (ROFI) of the Elbe–Weser rivers, the flood event clearly dominated the changes in salinity and nutrient concentrations, as expected. However, our findings point to the relevance of the peculiarities in the meteorological conditions in 2013 as well: a combination of low wind speeds, warm air temperatures and cold bottom-water temperatures resulted in a strong thermal stratification in the outer regions and limited vertical nutrient transport to the surface layers. Within the central region, the thermal and haline dynamics interactively resulted in an intense density stratification. This intense stratification, in turn, led to enhanced primary production within the central region enriched by nutrients due to the flood but led to reduction within the nutrient-limited outer region, and it caused a widespread oxygen depletion in bottom waters. Our results further point to the enhancement of the current velocities at the surface as a result of haline stratification and to intensification of the thermohaline estuarine-like circulation in the Wadden Sea, both driven by the flood event.