perfindices: Calculate 10 indices

In RHydro: Classes and methods for hydrological modelling and analysis

Description

This routine calculates 10 performance and stramflow indices

Usage

perfindices(Po,Qo,Qs,Qb)

Arguments

Po	observed precipitation [mm/timestep]
Qo	observed streamflow [mm/timestep]
Qs	simulated streamflow [mm/timestep]
Qb	simulated base flow [mm/timestep]

Details

perfindices() calculates 6 performance measures and 4 catchment descriptors (streamflow-based metrics). It uses simulated and observed streamflow, observed precipitation (to calculate the rainfall-runoff coefficient and streamflow elasticity) and simulated baseflow (to calculate the baseflow index, slope of duration curve, streamflow elastivity).

The 6 performance measures are:

• LAGTIME between observed and simulated discharge (from package tiger)

• MEAN ERROR between observed and simulated discharge

• MEAN ABSOLUTE ORDINAL ERROR between observed and simulated discharge (from package qualV)

• Complement to 1 of NASH-SUTCLIFFE EFFICIENCY between observed and simulated discharge (from package qualV)

• ROOT MEAN SQUARE ERROR between observed and simulated discharge (from package qualV)

• Difference in HIGH PULSE COUNT between observed and simulated discharge

All the above measures tend to zero as Qo approaches Qs.

The 4 catchment descriptors are:

• RAINFALL_RUNOFF coefficient (for Qs)

• BASEFLOW INDEX (for Qs)

• SLOPE of FLOW DURATION CURVE (for Qs)

• STREAMFLOW ELASTICITY (for Qs)

Value

The function returns a list of 10 indices.

Author(s)

Claudia Vitolo, Imperial College London; Till Francke, University of Potsdam

See Also

runoff_coeff, fdc_slope, sel,codebfi

Examples

data(modlist)
data(mopex)

# returns the instantaneous runoff
X <- fusesma.sim (fuse.DATA,mid=5,modlist,deltim=1,
                      states=FALSE,fluxes=TRUE,
                      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)
                        
# returns the routed runoff
Qrout <- fuserouting.sim(X$U, mid=5, modlist, fuse.parameters$timedelay, deltim=1)

Po <- fuse.DATA[,"P"]
Qo <- fuse.DATA[,"Q"]
Qs <- Qrout
Qb <- X$qbase_2
  
indices   <- perfindices(Po,Qo,Qs,Qb)

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