hydromt_wflow.WflowSedimentModel.setup_reservoirs#

WflowSedimentModel.setup_reservoirs(reservoirs_fn: str | Path | GeoDataFrame, overwrite_existing: bool = False, duplicate_id: str = 'error', min_area: float = 1.0, trapping_default: float = 1.0, output_names: dict = {'reservoir_area__count': 'reservoir_area_id', 'reservoir_location__count': 'reservoir_outlet_id', 'reservoir_surface__area': 'reservoir_area', 'reservoir_water_sediment~bedload__trapping_efficiency': 'reservoir_trapping_efficiency'}, geom_name: str = 'meta_reservoirs', **kwargs)[source]#

Generate maps of reservoir areas and outlets.

Also generates parameters with average reservoir area, and trapping efficiency for large particles.

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

Adds model layers:

  • reservoir_area_id map: reservoir IDs [-]

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

  • reservoir_area map: reservoir area [m2]

  • reservoir_trapping_efficiency map: reservoir bedload trapping efficiency coefficient [-] (0 for natural lakes, 0-1 depending on the type of dam)

  • meta_reservoirs geom: polygon with reservoirs and parameters

  • reservoirs geom: polygon with all reservoirs as in the model

Parameters:

  • reservoirs_fn (str) –

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

    • Required variables: ['waterbody_id', 'Area_avg']

    • Optional variables: ['reservoir_trapping_efficiency']

  • overwrite_existing (bool, optional) – If True, overwrite existing reservoirs in the model grid. Default is False.

  • duplicate_id ({"error", "skip"}, optional) – Action to take if duplicate reservoir IDs are found when merging with existing reservoirs. Options are "error" (default) to raise an error, or "skip" to skip adding new reservoirs.

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

  • trapping_default (float, optional) – Default trapping efficiency coefficient for large particles [0-1]. Default is 1.0, meaning 100% of large particles (sand to gravel) are trapped (e.g. in a gravity dam). For other dam types, the natural deposition rates from Camp are used.

  • output_names (dict, optional) – Dictionary with output names to be used in the model NetCDF input files. Keys should be Wflow.jl variable names, values the names in the NetCDF file.

  • geom_name (str, optional) – Name of the reservoirs geometry in the staticgeoms folder. Default is "meta_reservoirs" (for meta_reservoirs.geojson).

  • kwargs (dict, optional) – Additional keyword arguments passed to hydromt.DataCatalog.get_rasterdataset().