Abstract. Information on the frequency and intensity of extreme precipitation is required by public authorities, civil security departments, and engineers for the design of buildings and the dimensioning of water management and drainage schemes. Especially for sub-daily resolutions, at which many extreme precipitation events occur, the observational data are sparse in space and time, distributed heterogeneously over Europe, and often not publicly available. We therefore consider it necessary to provide an impact-orientated data set of 10-year rainfall return levels over Europe based on climate model simulations and evaluate its quality. Hence, to standardize procedures and provide comparable results, we apply a high-resolution single-model large ensemble (SMILE) of the Canadian Regional Climate Model version 5 (CRCM5) with 50 members in order to assess the frequency of heavy-precipitation events over Europe between 1980 and 2009. The application of a SMILE enables a robust estimation of extreme-rainfall return levels with the 50 members of 30-year climate simulations providing 1500 years of rainfall data. As the 50 members only differ due to the internal variability in the climate system, the impact of internal variability on the return level values can be quantified. We present 10-year rainfall return levels of hourly to 24 h durations with a spatial resolution of 0.11∘ (12.5 km), which are compared to a large data set of observation-based rainfall return levels of 16 European countries. This observation-based data set was newly compiled and homogenized for this study from 32 different sources. The rainfall return levels of the CRCM5 are able to reproduce the general spatial pattern of extreme precipitation for all sub-daily durations with Spearman's rank correlation coefficients >0.76 for the area covered by observations. Also, the rainfall intensity of the observational data set is in the range of the climate-model-generated intensities in 60 % (77 %, 78 %, 83 %, 78 %) of the area for hourly (3, 6, 12, 24 h) durations. This results in biases between −16.3 % (hourly) to +8.2 % (24 h) averaged over the study area. The range, which is introduced by the application of 50 members, shows a spread of −15 % to +18 % around the median. We conclude that our data set shows good agreement with the observations for 3 to 24 h durations in large parts of the study area. However, for an hourly duration and topographically complex regions such as the Alps and Norway, we argue that higher-resolution climate model simulations are needed to improve the results. The 10-year return level data are publicly available (Poschlod, 2020; https://doi.org/10.5281/zenodo.3878887).

No abstract

Abstract. 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 areas. In addition, we present a unique assessment of the uncertainty when using only 50 years of data (what is typically available from observations). Training the impact function with short records leads to a general underestimation of the return levels as water level extremes are not well sampled. Also, the marginal distributions of the 50-year time series of the surge show high variability. Moreover, compounding effects tend to be underestimated when using 50-year slices to estimate the dependence pattern between predictors. Overall, the internal variability of the climate system is identified as a major source of uncertainty in the multivariate statistical model.

Abstract. Extreme daily rainfall is an important trigger for floods in Bavaria. The dimensioning of water management structures as well as building codes is based on observational rainfall return levels. In this study, three high-resolution regional climate models (RCMs) are employed to produce 10- and 100-year daily rainfall return levels and their performance is evaluated by comparison to observational return levels. The study area is governed by different types of precipitation (stratiform, orographic, convectional) and a complex terrain, with convective precipitation also contributing to daily rainfall levels. The Canadian Regional Climate Model version 5 (CRCM5) at a 12 km spatial resolution and the Weather and Forecasting Research (WRF) model at a 5 km resolution both driven by ERA-Interim reanalysis data use parametrization schemes to simulate convection. WRF at a 1.5 km resolution driven by ERA5 reanalysis data explicitly resolves convectional processes. Applying the generalized extreme value (GEV) distribution, the CRCM5 setup can reproduce the observational 10-year return levels with an areal average bias of +6.6 % and a spatial Spearman rank correlation of ρ=0.72. The higher-resolution 5 km WRF setup is found to improve the performance in terms of bias (+4.7 %) and spatial correlation (ρ=0.82). However, the finer topographic details of the WRF-ERA5 return levels cannot be evaluated with the observation data because their spatial resolution is too low. Hence, this comparison shows no further improvement in the spatial correlation (ρ=0.82) but a small improvement in the bias (2.7 %) compared to the 5 km resolution setup. Uncertainties due to extreme value theory are explored by employing three further approaches. Applied to the WRF-ERA5 data, the GEV distributions with a fixed shape parameter (bias is +2.5 %; ρ=0.79) and the generalized Pareto (GP) distributions (bias is +2.9 %; ρ=0.81) show almost equivalent results for the 10-year return period, whereas the metastatistical extreme value (MEV) distribution leads to a slight underestimation (bias is −7.8 %; ρ=0.84). For the 100-year return level, however, the MEV distribution (bias is +2.7 %; ρ=0.73) outperforms the GEV distribution (bias is +13.3 %; ρ=0.66), the GEV distribution with fixed shape parameter (bias is +12.9 %; ρ=0.70), and the GP distribution (bias is +11.9 %; ρ=0.63). Hence, for applications where the return period is extrapolated, the MEV framework is recommended. From these results, it follows that high-resolution regional climate models are suitable for generating spatially homogeneous rainfall return level products. In regions with a sparse rain gauge density or low spatial representativeness of the stations due to complex topography, RCMs can support the observational data. Further, RCMs driven by global climate models with emission scenarios can project climate-change-induced alterations in rainfall return levels at regional to local scales. This can allow adjustment of structural design and, therefore, adaption to future precipitation conditions.

This study assesses the change of the seasonal runoff characteristics in 98 catchments in central Europe between the reference period of 1981–2010, and in the near future (2011–2040), mid future (2041–2070) and far future (2071–2099). Therefore, a large ensemble of 50 hydrological simulations featuring the model WaSiM-ETH driven by a 50-member ensemble of the Canadian Regional Climate Model, version 5 (CRCM5) under the emission scenario Representative Concentration Pathway (RCP 8.5) is analyzed. A hierarchical cluster analysis is applied to group the runoff characteristics into six flow regime classes. In the study area, (glacio-)nival, nival (transition), nivo-pluvial and three different pluvial classes are identified. We find that the characteristics of all six regime groups are severely affected by climate change in terms of the amplitude and timing of the monthly peaks and sinks. According to our simulations, the monthly peak of nival regimes will occur earlier in the season and the relative importance of rainfall increases towards the future. Pluvial regimes will become less balanced with higher normalized monthly discharge during January to March and a strong decrease during May to October. In comparison to the reference period, 8% of catchments will shift to another regime class until 2011–2040, whereas until 2041–2070 and 2071–2099, 23% and 43% will shift to another class, respectively.