fusesma.sim: Implementation of the framework for hydrological modelling...
In RHydro: Classes and methods for hydrological modelling and analysis

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

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

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)

`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 (-)

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):

Rainfall error (rferr)
- additive = 11
- multiplicative = 12
Architecture of the upper soil layer (arch1)
- Single state = 21
- Separate tension storage = 22
- Cascading buckets = 23
Architecture of the lower soil layer (arch2)
- Baseflow reservoir of fixed size = 31
- Tension reservoir plus two parallel tanks = 32
- Baseflow reservoir of unlimited size (frac rate) = 33
- Baseflow reservoir of unlimited size (power recession) = 34
Runoff (qsurf)
- Unsaturated zone Pareto = 41
- Unsaturated zone linear = 42
- Saturated zone topographic = 43
Percolation (qperc)
- Drainage above field capacity = 51
- Gravity drainage = 52
- Saturated zone control = 53
Evaporation (esoil)
- Sequential = 61
- Root weighting = 62
Interflows (qintf)
- Interflow denied = 71
- Interflow allowed = 72
Routing (q_tdh)
- Routing denied = 81
- Routing allowed using Gamma distribution = 82

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.

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.

Claudia Vitolo, Imperial College London

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

fuserouting.sim

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)