|
|
|
|
|
see [code](https://github.com/Ibasam/IBASAM/blob/master/IBASAM/R/river_climate_model.R)
|
|
|
|
|
|
The climate in the simulated river was described by daily water temperature and flow. As explained in the [input](howtouse) section, the time series of temperature and flow of a real river can be used for parameterization purpose or specific simulations.
|
|
|
|
|
|
However, for further simulations and to obtain realistic replicates of time series of these two environmental factors, a random autoregressive sinusoidal function was used for each of these factors:
|
... | ... | @@ -21,8 +18,10 @@ function (m, a, cs, alp, bet, d, nbd) |
|
|
where *m* and *a* were the mean and amplitude of the environmental parameter (temperature or flow) over a year, *d* was the Julian day, *cs* was a shifting parameter to adjust the sinusoidal signal and *error(dd)* was an autoregressive random error with *alp* as a non-zero parameter and *esp* as a normal random number centered on 0 and of variance *bet*.
|
|
|
|
|
|
The Scorff River data series (over 15 yrs) allowed a parameterization of these two functions:
|
|
|
(1) for the river water temperature (T(t));
|
|
|
(2) for the logarithm of the flow (flow(t) in m~3/s) divided by the module of the river (5.03 m~3/s).
|
|
|
1. for the river water temperature (T(t));
|
|
|
2. for the logarithm of the flow (flow(t) in m^3/s) divided by the module of the river (5.03 m^3/s).
|
|
|
|
|
|
Link to the [code](https://github.com/Ibasam/IBASAM/blob/master/IBASAM/R/river_climate_model.R)
|
|
|
|
|
|
``` R
|
|
|
river_climate_model <- function (nYears, CC_Temp, CC_Amp)
|
... | ... | |