reservoir_lumped: Generation of WASA parameter files for simulation of lumped...

In tpilz/lumpR: Landscape Unit Mapping Program for R

Generation of WASA parameter files for simulation of lumped reservoirs ("small reservoirs")

Description

Function generates the WASA parameter files lake.dat and lake_number.dat from a pre-processed reservoir vector and a subbasin raster file stored in a GRASS location.

Usage

reservoir_lumped(res_vect = NULL, subbas = NULL,
  res_vect_classified = NULL, dir_out = "./",
  reservoirs_small_classes_file = "reservoirs_small_classes.dat",
  r_subbas_contains_reservoirs_small_file = "r_subbas_contains_reservoirs_small.dat",
  res_param = data.frame(class = 1:5, f_vol_init = 0.2, class_change = 0,
  alpha_Molle = 2.7, damk_Molle = 1500, damc_hrr = c(7, 14, 21, 28, 35),
  damd_hrr = 1.5), keep_temp = F, overwrite = F, silent = F)

Arguments

res_vect	Name of reservoir vector map in GRASS location. Can be point or polygon feature (e.g. reservoir centroids). If it has the columns 'volume' with information on volume in [m^3] or column 'area' with information on lake area in [m^2], these will be used for computation. If none is encountered, 'area' can be generated from the area of the polygons, if present.
subbas	Name of subbasin raster map in GRASS location. Can be created with calc_subbas.
res_vect_classified	Output: Name for the vector reservoir map to be created in GRASS location. This is a point map based on res_vect with information of area or volume (whichever is missing), ID of the subbasin in which the reservoir is located, and the classified size class appended to the attribute table. If NULL (default) it will not be created.
dir_out	Character string specifying output directory (will be created if not available and files will be overwritten if overwrite = TRUE.
reservoirs_small_classes_file	Parameters for the reservoir size classes. See db_fill.
r_subbas_contains_reservoirs_small_file	Distribution of small reservoirs in subbasins. See db_fill.
res_param	A data.frame object containing parameters for the reservoir size classes. The default parameter set is adjusted to semi-arid Brazil. See Details.
keep_temp	logical. Set to TRUE if temporary files shall be kept in the GRASS location, e.g. for debugging or further analyses. Default: FALSE.
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

This function creates WASA input files which are needed to run the model with option doacudes=.T. (simulate small unlocated reservoirs).

The default values for res_param set were estimated by Molle (1989) for the semi-arid NE of Brazil and need to be adjusted if applied to some other region!

lake_file / res_param
Specification of parameters for the reservoir size classes. Note that lake_file and res_param have the same order, but different header names! If information on 'maxlake0' / 'vol_max' is not available, you can specify 'area_max' in res_param, i.e. the maximum area of reservoir size classes in m^2. This is internally converted to volume by relationship of Molle (1989) using parameters 'alpha_Molle' and 'damk_Molle'. If neither is given (default), the values will be estimated from the 20 percentiles of the value distribution of area or volume in res_vect by interpolation between both values using a logarithmic relationship.

Reservoir_class-ID / class
ID of reservoir size class.

maxlake0 / vol_max
Upper limit of reservoir size class in terms of volume in m^3.

lake_vol0_factor / f_vol_init
Fraction of storage capacity that indicates the initial water volume in the reservoir size classes (dimensionless).

lake_change / class_change
Factor that indicates yearly variation in the number of reservoirs of the size classes (dimensionless).

alpha_Molle, damk_Molle
Parameters of the area-volume relationship in the reservoir size classes: Area = alpha_Molle * damk_Molle * (Volume / damk_Molle)^( (alpha_Molle - 1) / alpha_Molle). Unit of Area: m^2, unit of Volume: m^3.

damc_hrr, damd_hrr
Parameters of the spillway rating curve in the reservoir size classes: Overflow = damc_hrr * Height^(damd_hrr). Unit of Overflow: m^3/s, unit of Height (over spillway): m.

Value

The two output files specified in reservoirs_small_classes_file and r_subbas_contains_reservoirs_small_file; and the vector map res_vect_class in GRASS

Note

Prepare GRASS location and necessary spatial objects in advance and start GRASS session in R using initGRASS.

Use reservoir_strategic to prepare the input for large / strategic reservoirs.

Points in res_vect not overlapping with any subbas will be silently removed during processing!

Polygons overlapping over various subbasins will be classified to the subbasin containing the polygon's centroid. If that occurs frequently, you should consider running reservoir_outlet in advance where reservoir outlet locations are estimated instead of using the centroid.

Author(s)

Tobias Pilz [e-mail], Till Francke [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.

Reservoir parameter set herein given as standard values:
Molle, F. (1989): Evaporation and infiltration losses in small reservoirs. Serie Hydrologia, 25, SUDENE / ORSTOM, Recife, Brazil, in Portuguese.

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