The co-occurrence of (not necessarily extreme) precipitation and surge can lead to extreme inland water levels in coastal areas. In a previous work the positive dependence between the two meteorological drivers was demonstrated in a managed water system in the Netherlands by empirically investigating an 800-year time series of water levels, which were simulated via a physical-based hydrological model driven by a regional climate model large ensemble. In this study, we present an impact-focused multivariate statistical framework to model the dependence between these flooding drivers and the resulting return periods of inland water levels. This framework is applied to the same managed water system using the aforementioned large ensemble. Composite analysis is used to guide the selection of suitable predictors and to obtain an impact function that optimally describes the relationship between high inland water levels (the impact) and the explanatory predictors. This is complex due to the high degree of human management affecting the dynamics of the water level. Training the impact function with subsets of data uniformly distributed along the range of water levels plays a major role in obtaining an unbiased performance. The dependence structure between the defined predictors is modelled using two- and three-dimensional copulas. These are used to generate paired synthetic precipitation and surge events, transformed into inland water levels via the impact function. The compounding effects of surge and precipitation and the return water level estimates fairly well reproduce the earlier results from the empirical analysis of the same regional climate model ensemble. Regarding the return levels, this is quantified by a root-mean-square deviation of 0.02 m. The proposed framework is able to produce robust estimates of compound extreme water levels for a highly managed hydrological system. Even though the framework has only been applied and validated in one study area, it shows great potential to be transferred to other ...