Hydrological models - data pre-processing

Overview

We provide a collection of scripts that pre-processes environmental data for use in several hydrological models:

PCR-GLOBWB

PCR-GLOBWB (PCRaster Global Water Balance) is a large-scale hydrological model intended for global to regional studies and developed at the Department of Physical Geography, Utrecht University (Netherlands). The recipe pre-processes ERA-Interim reanalyses data for use in the PCR-GLOBWB.

MARRMoT

MARRMoT (Modular Assessment of Rainfall-Runoff Models Toolbox) is a rainfall-runoff model comparison framework that allows objective comparison between different conceptual hydrological model structures https://github.com/wknoben/MARRMoT. The recipe pre-processes ERA-Interim and ERA5 reanalyses data for use in the MARRMoT.

MARRMoT requires potential evapotranspiration (evspsblpot). The variable evspsblpot is not available in ERA-Interim. Thus, we use the debruin function (De Bruin et al. 2016) to obtain evspsblpot using both ERA-Interim and ERA5. This function needs the variables tas, psl, rsds, and rsdt as input.

wflow_sbm and wflow_topoflex

Forcing data for the wflow_sbm and wflow_topoflex hydrological models can be prepared using recipe_wflow.yml.

HYPE

The hydrological catchment model HYPE simulates water flow and substances on their way from precipitation through soil, river and lakes to the river outlet. HYPE is developed at the Swedish Meteorological and Hydrological Institute. The recipe pre-processes ERA-Interim and ERA5 data for use in HYPE.

Available recipes and diagnostics

Recipes are stored in esmvaltool/recipes/hydrology

  • recipe_pcrglobwb.yml
  • recipe_marrmot.yml
  • recipe_wflow.yml
  • recipe_hype.yml

Diagnostics are stored in esmvaltool/diag_scripts/hydrology

  • pcrglobwb.py
  • marrmot.py
  • wflow.py
  • hype.py

User settings in recipe

  1. recipe_pcrglobwb.yml

    Required preprocessor settings:

    • start_year: 1979
    • end_year: 1979
  2. recipe_marrmot.yml

    There are two diagnostics, one for daily and one for hourly data.

    Required preprocessor settings:

    The settings below should not be changed.

    extract_shape:

    • shapefile: meuse_hydrosheds.shp (MARRMoT is a hydrological Lumped model that needs catchment-aggregated forcing data. The catchment is provided as a shapefile, the path can be relative to auxiliary_data_dir as defined in config-user.yml.).
    • method: contains
    • crop: true

    daily_statistics:

    • operator: mean (MARRMoT needs daily forcing data. Hourly forcing data are converted to daily values by mean operator).

    Required diagnostic script settings:

    • basin: Name of the catchment
  3. recipe_wflow.yml

    Required preprocessor settings:

    • extract_region: the region specified here should match the catchment
    • daily_statistics: if the frequency of the input data is not daily, it should be converted to daily using the preprocessor function daily_statistics with operator: mean.

    Required diagnostic script settings:

    • basin: name of the catchment
    • dem_file: netcdf file containing a digital elevation model with elevation in meters and coordinates latitude and longitude.
  4. recipe_hype.yml

    Required preprocessor settings:

    • start_year: 1979
    • end_year: 1979
    • shapefile: Meuse_HYPE.shp (expects shapefile with subcatchments)

    These settings should not be changed

    • method: contains
    • decomposed: true

Variables

  1. recipe_pcrglobwb.yml

    • tas (atmos, daily, longitude, latitude, time)
    • pr (atmos, daily, longitude, latitude, time)
  2. recipe_marrmot.yml

    • pr (atmos, daily or hourly mean, longitude, latitude, time)
    • psl (atmos, daily or hourly mean, longitude, latitude, time)
    • rsds (atmos, daily or hourly mean, longitude, latitude, time)
    • rsdt (atmos, daily or hourly mean, longitude, latitude, time)
    • tas (atmos, daily or hourly mean, longitude, latitude, time)
  3. recipe_wflow.yml

    • orog (fx, longitude, latitude)
    • pr (atmos, daily or hourly mean, longitude, latitude, time)
    • tas (atmos, daily or hourly mean, longitude, latitude, time)

    Either potential evapotranspiration can be provided:

    • evspsblpot(atmos, daily or hourly mean, longitude, latitude, time)

    or it can be derived from tas, psl, rsds, and rsdt using the De Bruin formula, in that case the following variables need to be provided:

    • psl (atmos, daily or hourly mean, longitude, latitude, time)
    • rsds (atmos, daily or hourly mean, longitude, latitude, time)
    • rsdt (atmos, daily or hourly mean, longitude, latitude, time)
  4. recipe_hype.yml

    • tas (atmos, daily or hourly, longitude, latitude, time)
    • tasmin (atmos, daily or hourly, longitude, latitude, time)
    • tasmax (atmos, daily or hourly, longitude, latitude, time)
    • pr (atmos, daily or hourly, longitude, latitude, time)

Observations and reformat scripts

Note: see headers of cmorization scripts (in esmvaltool/cmorizers/obs) for download instructions.

  • ERA-Interim (esmvaltool/cmorizers/obs/cmorize_obs_era_interim.py)
  • ERA5 (esmvaltool/cmorizers/obs/cmorize_obs_era5.py)

Output

  1. recipe_pcrglobwb.yml

  2. recipe_marrmot.yml

    The forcing data, the start and end times of the forcing data, the latitude and longitude of the catchment are saved in a .mat file as a data structure readable by MATLAB or Octave.

  3. recipe_wflow.yml

    The forcing data, stored in a single NetCDF file.

References

  • Sutanudjaja, E. H., van Beek, R., Wanders, N., Wada, Y., Bosmans, J. H. C., Drost, N., van der Ent, R. J., de Graaf, I. E. M., Hoch, J. M., de Jong, K., Karssenberg, D., López López, P., Peßenteiner, S., Schmitz, O., Straatsma, M. W., Vannametee, E., Wisser, D., and Bierkens, M. F. P.: PCR-GLOBWB 2: a 5 arcmin global hydrological and water resources model, Geosci. Model Dev., 11, 2429-2453, https://doi.org/10.5194/gmd-11-2429-2018, 2018.
  • De Bruin, H. A. R., Trigo, I. F., Bosveld, F. C., Meirink, J. F.: A Thermodynamically Based Model for Actual Evapotranspiration of an Extensive Grass Field Close to FAO Reference, Suitable for Remote Sensing Application, American Meteorological Society, 17, 1373-1382, DOI: 10.1175/JHM-D-15-0006.1, 2016.
  • Arheimer, B., Lindström, G., Pers, C., Rosberg, J. och J. Strömqvist, 2008. Development and test of a new Swedish water quality model for small-scale and large-scale applications. XXV Nordic Hydrological Conference, Reykjavik, August 11-13, 2008. NHP Report No. 50, pp. 483-492.
  • Lindström, G., Pers, C.P., Rosberg, R., Strömqvist, J., Arheimer, B. 2010. Development and test of the HYPE (Hydrological Predictions for the Environment) model – A water quality model for different spatial scales. Hydrology Research 41.3-4:295-319.