hydromt_wflow.WflowSedimentModel.setup_reservoirs#

WflowSedimentModel.setup_reservoirs(reservoirs_fn: str | Path | GeoDataFrame, timeseries_fn: str | Path | DataFrame = None, min_area: float = 1.0, **kwargs)[source]#

Generate maps of lake areas and outlets.

Also generates well as parameters with average reservoir area, demand, min and max target storage capacities and discharge capacity values.

The data is generated from features with min_area [km2] (default is 1 km2) from a database with reservoir geometry, IDs and metadata.

Data requirements for direct use (ie wflow parameters are data already present in reservoirs_fn) are reservoir ID ‘waterbody_id’, area ‘ResSimpleArea’ [m2], maximum volume ‘ResMaxVolume’ [m3], the targeted minimum and maximum fraction of water volume in the reservoir ‘ResTargetMinFrac’ and ‘ResTargetMaxFrac’ [-], the average water demand ResDemand [m3/s] and the maximum release of the reservoir before spilling ‘ResMaxRelease’ [m3/s].

In case the wflow parameters are not directly available they can be computed by HydroMT based on time series of reservoir surface water area. These time series can be retrieved from either the hydroengine or the gwwapi, based on the Hylak_id the reservoir, found in the GrandD database.

The required variables for computation of the parameters with time series data are reservoir ID ‘waterbody_id’, reservoir ID in the HydroLAKES database ‘Hylak_id’, average volume ‘Vol_avg’ [m3], average depth ‘Depth_avg’ [m], average discharge ‘Dis_avg’ [m3/s] and dam height ‘Dam_height’ [m]. To compute parameters without using time series data, the required variables in reservoirs_fn are reservoir ID ‘waterbody_id’, average area ‘Area_avg’ [m2], average volume ‘Vol_avg’ [m3], average depth ‘Depth_avg’ [m], average discharge ‘Dis_avg’ [m3/s] and dam height ‘Dam_height’ [m] and minimum / normal / maximum storage capacity of the dam ‘Capacity_min’, ‘Capacity_norm’, ‘Capacity_max’ [m3].

Adds model layers:

  • wflow_reservoirareas map: reservoir IDs [-]

  • wflow_reservoirlocs map: reservoir IDs at outlet locations [-]

  • ResSimpleArea map: reservoir area [m2]

  • ResMaxVolume map: reservoir max volume [m3]

  • ResTargetMinFrac map: reservoir target min frac [m3/m3]

  • ResTargetFullFrac map: reservoir target full frac [m3/m3]

  • ResDemand map: reservoir demand flow [m3/s]

  • ResMaxRelease map: reservoir max release flow [m3/s]

  • reservoirs geom: polygon with reservoirs and wflow reservoir parameters

Parameters:

  • reservoirs_fn (str) –

    Name of data source for reservoir parameters, see data/data_sources.yml.

    • Required variables for direct use: [‘waterbody_id’, ‘ResSimpleArea’, ‘ResMaxVolume’, ‘ResTargetMinFrac’, ‘ResTargetFullFrac’, ‘ResDemand’, ‘ResMaxRelease’]

    • Required variables for computation with timeseries_fn: [‘waterbody_id’, ‘Hylak_id’, ‘Vol_avg’, ‘Depth_avg’, ‘Dis_avg’, ‘Dam_height’]

    • Required variables for computation without timeseries_fn: [‘waterbody_id’, ‘Area_avg’, ‘Vol_avg’, ‘Depth_avg’, ‘Dis_avg’, ‘Capacity_max’, ‘Capacity_norm’, ‘Capacity_min’, ‘Dam_height’]

  • timeseries_fn (str {'gww', 'hydroengine', 'none'}, optional) – Download and use time series of reservoir surface water area to calculate and overwrite the reservoir volume/areas of the data source. Timeseries are either downloaded from Global Water Watch ‘gww’ (using gwwapi package) or JRC ‘jrc’ (using hydroengine package). By default None.

  • min_area (float, optional) – Minimum reservoir area threshold [km2], by default 1.0 km2.