Combining observations and a high-resolution climate model for the generation of spatial rainfall return levels

Rainfall return levels are guiding hazard protection, insurance models, infrastructure design, construction, and planning of cities. When deriving information about the frequency and intensity of extremes by fitting extreme value models to pointwise observations, the regionalization of these models is challenging. Rain gauges are distributed unevenly, where some regions suffer from data scarcity in space and time. To address this, topographical and/or climatological covariates are often used for the spatial interpolation. On the other hand, high-resolution climate simulations are available to provide spatial information on rainfall extremes. However, these simulations are still governed by model biases, where the bias adjustment of extremes at ungauged locations is also inducing uncertainty.   

In this study, we propose a combination of observations and a high-resolution convection-permitting climate model simulation in the framework of smooth spatial Generalized Extreme Value (GEV) models in order to estimate spatial rainfall return levels. We choose a study area over southern Germany with complex terrain, which is densely monitored with 1132 rain gauges providing more than 30-year daily rainfall observations. There, a 30-year simulation of the Weather and Forecasting Research (WRF) model is available at 1.5 km resolution driven by ERA5 reanalysis data. We combine observations and covariates from the WRF simulation in the spatial GEV and refer to this approach as sGEV-WRF.

We want to answer three research questions to assess the added value of the proposed framework:

• Is it worth the effort? Does the sGEV-WRF improve the generation of rainfall return levels compared to the WRF alone?
• Does the WRF simulation as covariate add value? Can the sGEV-WRF outperform a topography-only spatial GEV?
• Does the dynamical downscaling at high resolution add value? Can sGEV-WRF outperform a spatial GEV based on observations and covariates from the coarser resolution ERA5?

By evaluating the percentage bias, mean absolute error, and root-mean-square error, we show that the combination of observations and WRF can improve the representation of 10-year and 100-year return levels of daily rainfall.

In addition, we aim to assess the performance of this framework under data-scarce conditions. Therefore, we devise an extensive cross-validation study. We select 5%, 10%, 20%, 50%, 80%, 90%, and 95% of all 1132 rain gauges to re-build the spatial GEV models with 1000 random folds each. We show that the performance is robust under these conditions, highlighting the potential for the application in data-scarce regions. Furthermore, in a non-stationary setup with climate model future projections, it can serve as a reliable tool to assess climate change effects on heavy rainfall.