reservoir_strategic: Assemble of reservoir parameter file for WASA for import into...

In tpilz/lumpR: Landscape Unit Mapping Program for R

Assemble of reservoir parameter file for WASA for import into database with db_fill()

Description

Function generates the reservoir parameter file from a pre-processed reservoir vector map and optional supplemental parameter file.

Usage

reservoir_strategic(res_vect = NULL, subbasin = NULL,
  res_file = NULL, dir_out = "./", reservoir_file = "reservoir.txt",
  overwrite = F, silent = F)

Arguments

res_vect	Name of reservoir vector map in GRASS location. Should be point instead of polygon feature (i.e. reservoir outlet locations; consider function reservoir_outlet)! For mandatory columns in attribute table see Details. In any case, column subbas_id containing the subbasin ID derived from subbasin will be added to the attribute table.
subbasin	Subbasin raster map in GRASS location. Can be created with calc_subbas.
res_file	tab-delimited file containing reservoir properties (fields see details) and key res_id. If set to NULL, these attributes must be contained in the attribute table of res_vect. If they are missing, they will be set to deafult values (see Details).
dir_out	Character string specifying output directory (will be created if it does not yet exist). The result file can later be imported to the database with db_fill()
reservoir_file	Output: File of parameters for the strategic reservoirs assigned to subbasins. To be filled into a database using db_fill(). This file is NOT directly compatible with WASA-SED!
overwrite	logical. Shall output of previous calls of this function be deleted? If FALSE the function returns an error if output already exists. Default: FALSE.
silent	logical. Shall the function be silent (also suppressing warnings of internally used GRASS functions)? Default: FALSE.

Details

For each reservoir that should be modelled explicitly within WASA the following information need to be collected and written into the vector file's attribute table or res_file. Column order is not important. Additional columns can be given but will be ignored:

res_id
Unique numeric reservoir identifier (if res_file is given, it also needs to be defined in the vector file's attribute table!).

name
OPTIONAL: name of the reservoir. Will be filled with <NA> if not given.

minlevel
Initial minimum level in the reservoir [m]. Value varies because of sediment accumulation. Default: 0.

maxlevel
Maximum water level in the reservoir [m]. Default: Estimated using Molle's equation (alpha = 2.7, k = 1500).

vol0
Initial volume of the reservoir [10^3 m^3]. Value varies because of sediment accumulation. Set to '-999' if information is not available. Default: Estimated using Molle's equation (alpha = 2.7, k = 1500).

storecap
Initial storage capacity of the reservoir [10^3 m^3]. Value varies because of sediment accumulation. DEfault: vol0

damflow
Target outflow discharge of the reservoir (90 % reliability) [m^3/s]. Default: 9.99

damq_frac
Fraction of Q90 released from the reservoir in regular years [-]. Default: 1

withdrawal
Water withdrawal discharge from the reservoir to supply the water use sectors [m^3/s]. Outflow discharge through the dam is not considered. Default: 0

damyear
Year of construction of the dam (YYYY). Default: 1900

maxdamarea
Initial maximum area of the reservoir [ha]. Value varies because of sediment accumulation. Default from GIS

damdead
Initial dead volume of the reservoir [10^3 m^3]. Value varies because of sediment accumulation. Default: 0.01 * storecap

damalert
Initial alert volume of the reservoir [10^3 m^3]. Value varies because of sediment accumulation. Default: 0.01 * storecap

dama, damb
Parameters of the area-volume relationship in the reservoir: area = dama * Vol^damb [-]. Values of reservoir area and volume are expressed in m^2 and m^3, respectively. Attention: These are not the Molle parameters! The correct values can be derived from them as shown in the default. Default: dama = (1/k * (alpha * k)^(alpha/(alpha-1)))^((alpha-1)/alpha), damb = (alpha-1)/alpha

q_outlet
Maximum outflow discharge released through the bottom outlets of the reservoir [m^3/s]. Default: 0.999

fvol_botm
Fraction of storage capacity that indicates the minimum storage volume for sediment release through the bottom outlets of the reservoir [-]. Default: 0.01

fvol_over
Fraction of storage capacity that indicates the minimum storage volume for water release through the spillway of the reservoir [-]. Default: 1

damc, damd
Parameters of the spillway rating curve of the reservoir: Qout = damc * Hv^damd [-]. Values of water height over the spillway and overflow discharges are expressed in m and m^3/s, respectively. Default: 1.6 * 30, 1.5

elevbottom
Bottom outlet elevation of the reservoir [m]. Currently ignored in WASA. Default: 99

The output file reservoir_file contains the additional column pid which is the corresponding subbasin ID determined from input subbasin.

Note

If you applied reservoir_outlet using the outlet locations of strategic reservoirs as drainage points, it might be necessary to use locations of the function's output points_processed instead of the true reservoir outlet locations as otherwise the reservoirs might get assigned to the wrong subbasins!

Author(s)

Tobias Pilz [e-mail]

References

lumpR package introduction with literature study and sensitivity analysis:
Pilz, T.; Francke, T.; Bronstert, A. (2017): lumpR 2.0.0: an R package facilitating landscape discretisation for hillslope-based hydrological models. Geosci. Model Dev., 10, 3001-3023, doi: 10.5194/gmd-10-3001-2017

WASA model in general:
Guentner, A. (2002): Large-scale hydrological modelling in the semi-arid North-East of Brazil. PIK Report 77, Potsdam Institute for Climate Impact Research, Potsdam, Germany.

Reservoir module of the WASA model:
Mamede, G. L. (2008): Reservoir sedimentation in dryland catchments: Modeling and management. PhD Thesis, University of Potsdam, Germany.

tpilz/lumpR documentation built on Aug. 5, 2023, 1:31 a.m.