Please note that this product is experimental and is not operationally supported, meaning that there may be some gaps in the availability on the EFAS web interface due to outages.

The radar-based river flash flood impact layer forecasts 1-arcminute (~1.4 km) river grid cells where there is a risk of impacts within the next 5 days due to river flash flooding, which refers to the rapid rise of small rivers (with a drainage area <=1000 km2) due to intense rainfall. The forecasts are generated every hour using a blending of ensemble precipitation forecast inputs derived from the latest radar observations and NWP forecasts from the ECMWF ensemble (the blending procedure is described in the methodology section of Radar-based Accumulated Precipitation 80th Percentile). The precipitation forecasts are used to predict grid cells where there is a probability of flash flooding exceeding the 2 year return period threshold. The flash flood probability is overlaid with a static dataset which summarises the exposure of population and critical infrastructure. The result of this overlay is used to categorise flash flood impacts on an impact matrix as being either low (yellow), medium (orange), high (red) or severe (purple), river grid cells affected by flash flooding will be highlighted in one of these four colours (Figure 1).

The forecasts have a 1-hour time step for the first 6-hours of the forecast, and a 6-hour time step thereafter up to a maximum lead time of 120-hours (5-days).

Visualising on the EFAS Web Interface

On the EFAS-IS web interface, the layer can be selected from the Flash Flood product tab. Once loaded, the user can use a drop down menu (Figure 1, bottom left corner) to select forecasts generated during the past 5 days. An animation slider (Figure 1, bottom left corner) allows each forecast time step to be visualised.  Please note, this layer is only available for the current and previous 5 days.

Figure 1. The Radar-based river flash flood impact layer shown on the EFAS web interface for the Dublin area on 4th September 2022 at 04:00 UTC (left), and the impact matrix (right) which is used to assign the colour of each river grid cell affected by flash flooding.

Suggested Workflow for Analysing the Radar-based River Flash Flood Impact Layer

When viewing the animated Radar-based river flash flood impact layer, it can be difficult to identify specific locations and time steps which may be affected by flash flooding. This is partly because when zoomed out it is very difficult to see the shaded 1-arcminute grid cells. To solve this, the following workflow is recommended:

  1.  Load in the Radar-based river flash flood impact catchment summary layer for the forecast lead time horizon window of interest e.g. load in the Radar-based River Flash Flood Impact Catchment summary 0-6h if you are interested in the first 6 hours of the forecast leadtime horizon 
  2. Use this layer to identify the specific sub-catchments which are highlighted as being at risk of flash flooding during this forecast horizon
  3. Zoom in closer to the highlighted subcatchments of interest
  4. Load in the Radar-based river flash flood impact layer, use the animation slider to identify the specific 1-arcminute grid cells and time steps where flash flood impacts are possible

The visibility of the Radar-based river flash flood impact layer can be improved by changing the map background to show country borders only, this can be toggled from the top right of the webviewer, it is the icon to the right of the WMS icon

Methodology

There are two main steps in the generation of the river flash flood impact layer: 1) river flash flood hazard forecasting using blended radar and NWP precipitation, 2) the prediction of river flash flood impact (Figure 2).

Figure 2. The impact matrix upon which flash flood hazard probability (y-axis) is combined with exposure (x-axis) to decide the impact category as shown by the colour of each cell. The definition of each impact category is shown on the right hand side.

Step 1: River Flash Flood Hazard Forecasts

A probabilistic forecast of river flash flood hazard is generated firstly by blending hourly nowcasts of precipitation derived from radar observations with NWP precipitation forecast from the ECMWF ensemble. More details about how the two datasets are blended can be found in the methodology section of: Radar-based Accumulated Precipitation 80th Percentile. Radar data are not available in all areas covered in the EFAS spatial domain, a map of radar coverage can be found in the Radar Coverage layer. The resulting blended 51-member ensemble precipitation forecasts have a spatial resolution of 2 km, and 1-hour time steps for the first 12 hours lead time which increases to 6-hours thereafter for a maximum lead time of 120-hours (5 days).

Next, at each time step of the forecast, at each grid cell with an upstream drainage area <=1000 km2, the rainfall is summarised as the total amount of rain which was forecasted to fall across the area upstream of the grid cell over the preceding time period equal to the river's time of concentration, Tc. The time of concentration refers to the time between rain falling and the river reaching its peak, it was computed at each grid cell using the method of Giandotti, 1934 which was derived for small rural catchments:

Tc = (4√A + 1.5L) / 0.8 √H

where A = upstream drainage area (km2), L = the slope length (km) and H = difference between the mean elevation in the area upstream of the grid cell and the elevation of the grid cell (m) 

The summarised rainfall is compared against a climatological threshold, which is the summarised rainfall associated with the 2-year return period. This climatological threshold is derived from the following datasets:

  • 8-year gauge-adjusted OPERA radar rainfall data
  • A dataset of 20-year reforecasts obtained with ECMWF Integrated Forecasting System (IFS)

At each location on the river network and each forecast time step, the total number of ensemble members which exceed the summarised rainfall associated with the 2-year return period is computed and then divided by the total number of ensemble members (51) to give a probability of exceedance (%).

Step 2: River Flash Flood Impact Prediction

The next step is to compute the exposure, at each river location where the exceedance probability of the 2-year return period of accumulated rainfall (step 1).

The exposure data accounts for population and critical infrastructure in the form of health, education, transport, and energy generation facilities. Population data were obtained from the GHSL (Global Human Settlement Layer), the critical infrastructure data were from HARCI-EU within EU member states and OpenStreetMap for non-EU member states. Data for each of these five categories was harmonised and combined with equal weighting to create a combined exposure layer whose values ranged from 1.0-2.0.

The river flash flood impact level is computed at each time step by intersecting the river flash flood hazard forecast with the combined exposure dataset. The low, medium, and high values for flash flood hazard likelihood on the y-axis of the impact matrix (Figure 2) indicate where there is a 5%-50%, 50%-80%, and >80% probability of exceeding the 2-year return period threshold.

The categories for low, medium, and high exposure, on the x-axis of the impact matrix (Figure 2), relate to combined exposure values of 1.0-1.3, 1.3-1.6 and 1.6-2.0 respectively. These categories were chosen based on the statistical distribution of exposure values across the EFAS domain, and consequently mean 81.0% the exposure values are classified as low exposure, 8.0% as medium exposure and 1.2% as high exposure. This reflects the reality that most grid points in Europe have a low population density and with few exposed critical infrastructures

This combination of flash flood hazard level with exposure gives impact levels for each grid cell on the river network at each time step of the 120 hour forecast period. 

  • No labels