Geostatistical Enhancement of European Hydrological Predictions (GEEHP)

This tutorial will guide the user through the main steps for the implementation of the procedure developed in the GEEHP experiment. For more details see the experiment protocol.

Retrieve input data and scripts from GitHub

All R-scripts and other resources can be retrieved from GitHub. Download the project folder into your working directory. A quick method for Linux users is:

git clone https://github.com/alessio-pugliese/geehp

Create a working directory (wd) in your local machine, or make the downloaded folder from git-hub your wd. Then, download streamflow data and catchment boundaries of Tyrol region from GitHub:

  1. EHYPE catchment shape files;

  2. EHYPE streamflows;

  3. Observed catchment shape files;

  4. Observed streamflows (IT);

  5. Observed streamflows (AU);

The code will look at all the shape file into a sub-folder of your wd, thus create a sub-folder called “shp” where all the GIS files must be stored.

Run the code

Cross-validation of the regional model for predicting FDCs.

Install and load all the required packages:

library(rtop)
library(maptools)
library(rgdal)
library(pracma)
library(nortest)
library(rgeos)

then, run the leave-one-out cross-validation algorithm with the command below. This commands will provide some graphical outputs, so make sure your data fit the code correctly.

source(file = "GEEHP_FDCs_LOOCV.R")

Predict FDCs on EHYPE sites and enhance simulations

The script GEEHP_selected_catchments.r computes FDCs at EHYPE catchments assuming these sites are ungauged, then predicted FDCs are used to enhance EHYPE simulations, running the procedure. Each EHYPE catchment selected for testing the procedure is paired with a gauged catchment and each pair assumes that either the EHYPE or the observed catchment represent the same basin. The table below report selected catchment pairs.

Country EHYPE OBS A_OBS A_EHYPE Diff_A Dist Quality
IT 9001070 3675 207.4 215.0 3.6% 1.6 NA
IT 9740616 5715 419.7 428.3 2.0% 2.1 NA
IT 9783013 6735 1929.4 2148.0 11.3% 9.9 NA
IT 9001106 6979 2899.1 3152.9 8.8% 12.5 NA
IT 9782936 7989 4199.3 3770.0 10.2% 14.4 NA
IT 9783018 8555 7000.6 7033.4 0.5% 6.8 NA
AT 9608296 201236 394.2 410.3 4.1% 3.9 1
AT 9608974 201277 270.7 274.0 1.2% 2.5 1
AT 9608705 201392 512.0 478.0 6.6% 5.1 0
AT 9616258 201434 785.2 892.2 13.6% 8.6 0
AT 9602501 201624 853.6 854.2 0.1% 1.5 0

In this phase the regional FDCs model must be validated and checked for reliability beforehand. The command below will provide a lot of graphical outputs, so make sure your data fit the code correctly. The script provide graphical comparison between EHYPE simulations with corrected EHYPE (cEHYPE). Plus, a focus on the best improvement among all pairs is also provided.

source(file = "GEEHP_selected_catchments.r")