by Rico Kronenberg and Lisa Marie Oehlschlägel
This manual describes how to run the BROOK90 as R implementation, based on the model by Federer, what files are necessary to download and what changes need to be done in the R-scripts.
Figure 1: BROOK90 Scheme.
- Overview of available files
- Download the data
- Run the programm
- The output
- Temporal model resolution (daily or subdaily)
These files are available for usage:
- html-Files
- to download and read the full documentation
- plot output
- example plot of calculated values
- Documentation including the written R scripts to run the programm
- html files are easier to read, Rmd-files for running the R-scripts
- precipitation data
- P_WBG_0060_Minuten_1999_point.txt
- meteorological data
- data_input_WB_1999_point.txt
- catchment parameters
- 6 txt-files where the characteristics of the catchment can be changed
Data files have to be formated as described HERE. The minimal amount of meteorological information needed to run the model is daily precipitation, mininmal and maximal temperatures. Other variables have to be set to zero when not available.
As the second step the files can be downloaded and unpacked.
To open the Rmd-files it is necessary to open R Studio first and to install the package for R-Markdown.
install.packages("rmarkdown")
The next step is to open the Rmd-script B90V4.Rmd. There the "projectpath" has to be changed in the first R-chunk. How this must be done is described in the file. Then MainProg.Rmd has to be opened. Important changes that have to be done in MainProg.Rmd are:
- change the "projectpath", it has to be the same as in B90V4.Rmd
- the input data (meteorological data, precipitatin data, catchment parameters) have to be in the same folder
- is already done if you unpack the zip-folder
- change catchment parameters (located in "Input_data" - folder) in these files, if necessary:
- canopy.txt
- fixed.txt
- flow_standart.txt
- initial.txt
- location.txt
- soil.txt
- if you are NOT using the example data, change the names of your input data files (meteoFile, precFile,...) and check the form of your data with consistency_meteoFile.Rmd.
The next step is to open again the script B90V4.Rmd. At the end of this script there is a code chunk to plot the output data. In this case precipitation, evaporation, calculated and measured flow are plotted. If other data should be shown, add or exchange the output data with "timeseries_" in its name. Finally B90V4.Rmd can also be started with "Run All" (click on the little arrow button next to "Run") and then the model is running and producing the output.
The programm calculates 24 daily time series of different water balance components in the considered catchment as output of BROOK90 at the current state. The table shows the explanation of the output data visible with the name "timeseries_..." and their corresponding shortcut:
| Shortcut | Explanation [d] |
|---|---|
| adef | available water deficit in root zone |
| awat | available soil water in root zone |
| evp | evapotranspiration |
| flow | total flow |
| irvp | evaporation rate of intercepted rain |
| isvp | evaporation rate of intercepted snow |
| mesld | measured flow |
| pint | average potential interception for day |
| prec | precipitation |
| ptran | average potential transpiration rate for day |
| rfald | rainfall rate |
| rintd | rain interception |
| rnet | rain reaching soil surface |
| rthrd | rain throughfall |
| sfald | snowfall |
| sintd | snow interception |
| slfld | input to soil surface |
| slvp | soil evaporation |
| smltd | snowmelt |
| snow | water equivalent of snow on the ground |
| snvp | evaporation from snowpack |
| sthrd | snow throughfall |
| swat | total soil water in all layers |
| trand | transpiration |
As mentioned in Run the programm, all of this time series can be plotted if you add or exchange time series at the end of the script B90V4.Rmd.
The default setting in this version is in daily resolution but the model can be driven by subdaily data like hourly data. This can be changed in the MainProg.Rmd (chunk2). For hourly values NPINT has to be changed from 1 to 24 and SUBDAYDATA has to be set TRUE.
| Parameter | Default setting | Description |
|---|---|---|
| NPINT | 1 | number of precipitation intervals per day, 1 = 1 day, 24 = 24 hours per day |
| SUBDAYDATA | FALSE | FALSE or TRUE |
The following lines show the part of MainProg.Rmd where the changes can be done:
NPINT <- 1 #change to 24
SUBDAYDATA <- FALSE #change to TRUE
RRD <- 0.55
UWD <- 3.0
The change of the temporal resolution varies the output. The results for daily and hourly resolution is shown in the next two figures.
Figure 2: Temporal model resolution daily
Figure 3: Temporal model resolution hourly