fusesma.sim: Implementation of the framework for hydrological modelling...

Description Usage Arguments Details Value Author(s) References See Also Examples

Description

Loss module derived from the Fortran version of FUSE by Martyn Clark (2011).

Usage

1
2
3
4
5
fusesma.sim(DATA,mid,modlist,deltim, states, fluxes, 
    fracstate0, rferr_add, rferr_mlt, 
    frchzne,fracten,maxwatr_1,percfrac,fprimqb,qbrate_2a,qbrate_2b,
    qb_prms,maxwatr_2,baserte,rtfrac1,percrte,percexp,sacpmlt,
    sacpexp,iflwrte,axv_bexp,sareamax,loglamb,tishape,qb_powr)

Arguments

DATA

data.frame containing observations. It consists of 3 columns: Rainfall (P), Potential Evapo-Transpiration (E) and Streamflow (Q)

mid

model id number in Model List 2011(see below for details)

modlist

list of availabe models

deltim

observation time step (days)

states

list of internal state variables

fluxes

list of fluxes

fracstate0

initial store fraction (initialization)

rferr_add

additive rainfall error (mm)

rferr_mlt

multiplicative rainfall error (-)

frchzne

fraction tension storage in recharge zone (-)

fracten

fraction total storage in tension storage (-)

maxwatr_1

depth of the upper soil layer (mm)

percfrac

fraction of percolation to tension storage (-)

fprimqb

fraction storage in 1st baseflow reservoir (-)

qbrate_2a

baseflow depletion rate 1st reservoir (day-1)

qbrate_2b

baseflow depletion rate 2nd reservoir (day-1)

qb_prms

baseflow depletion rate (day-1)

maxwatr_2

depth of the lower soil layer (mm)

baserte

baseflow rate (mm day-1)

rtfrac1

fraction of roots in the upper layer (-)

percrte

percolation rate (mm day-1)

percexp

percolation exponent (-)

sacpmlt

SAC model percltn mult for dry soil layer (-)

sacpexp

SAC model percltn exp for dry soil layer (-)

iflwrte

interflow rate (mm day-1)

axv_bexp

ARNO/VIC "b" exponent (-)

sareamax

maximum saturated area (-)

loglamb

mean value of the topographic index (m)

tishape

shape param for the topo index Gamma dist (-)

qb_powr

baseflow exponent (-)

Details

fusesma.sim() is a function to generate an ensemble of SOIL MOISTURE ACCOUNTING models. It is compatible with the HYDROMAD framework (see hydromad package: http://hydromad.catchment.org/).

fusesma.sim() can simulate several model structures. The default list is a modlist dataframe contained in the data folder of this package. In modlist each row identifies a model structure, and each column identifies a different element of the model structure:

first column mid = model id number
second column rferr = type of rainfall error (optional)
third column arch1 = architecture of the upper soil layer
fourth column arch2 = architecture of the lower soil layer
fifth column qsurf = surface runoff
sixth column qperc = vertical drainage
seventh column esoil = evapotranspiration
eighth column qintf = interflow
ninth column q_tdh = routing

For each element of the model structure, several model decisions can be made (see Clark et al. 2011 for details):

For instance, model 5 is identified by the following string (from the table above reading first the row number then column number for each model decision): 5 11 22 33 41 51 62 71 82

The parameter set varies depending on the selected model structure. Ranges of parameter values are in fusesma.ranges

Flow can then be routed using the function fuserouting.sim (which is based on the Gamma function) or any other routing function.

Value

The function returns an array of simulated "instantaneous" discharges. If necessary, fuserouting.sim can be run to obtain routed discharges using a two parameter Gamma distribution.

Author(s)

Claudia Vitolo, Imperial College London

References

Clark M. P., SlaterA. G., Rupp D. E., Woods R. A., Vrugt J. A., Gupta H. V., Wagener T. and Hay L. E. (2008), Framework for Understanding Structural Errors (FUSE): A modular framework to diagnose differences between hydrological models, Water Resour. Res. 44 p. 91-94

Clark M. P., McMillan H. K., Collins D. B. G., Kavetski D. and Woods R. A. (2011), Hydrological field data from a modeller's perspective: Part 2: process-based evaluation of model hypotheses. Hydrological Processes, 25: 523-543. doi: 10.1002/hyp.7902

See Also

fuserouting.sim

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
data(modlist)
data(mopex)

# returns the instantaneous runoff
Qinst <- fusesma.sim (fuse.DATA,mid=5,modlist,deltim=1,
                      states=FALSE,fluxes=FALSE,
                      fracstate0=0.25,
                      fuse.parameters$rferr_add,fuse.parameters$rferr_mlt,
                      fuse.parameters$frchzne,fuse.parameters$fracten,
                      fuse.parameters$maxwatr_1,fuse.parameters$percfrac,
                      fuse.parameters$fprimqb,fuse.parameters$qbrate_2a,
                      fuse.parameters$qbrate_2b,fuse.parameters$qb_prms,
                      fuse.parameters$maxwatr_2,fuse.parameters$baserte,
                      fuse.parameters$rtfrac1,fuse.parameters$percrte,
                      fuse.parameters$percexp,fuse.parameters$sacpmlt,
                      fuse.parameters$sacpexp,fuse.parameters$iflwrte,
                      fuse.parameters$axv_bexp,fuse.parameters$sareamax,
                      fuse.parameters$loglamb,fuse.parameters$tishape,
                      fuse.parameters$qb_powr)

RHydro documentation built on May 2, 2019, 6:24 p.m.