rating: Load estimation with log log rating curve

Description

Estimate the load using a log-log rating curve between values of concentration and river flow at the time of sampling (Quilbe' et al., 2006). The estimation is performed on the time period spanned by flow records. It is also possible to estimate the load monthly or annually specifying optional argument.

Usage

rating(db, ncomp, period)

Arguments

An input data frame with at least three column. A column with DateTime in standard format (yyyy-mm-dd HH:MM:SS) labeled "datetime", a column with flow records in cubic meter per second labeled "flow", and columns with scattered concentration data in milligram per litre. Alternatively, it can be used the output matrix of the function db.union.

ncomp

Number of compounds in the input data frame for which the load must be estimated.

period

Optional argument. It can be "month" for a monthly estimation of load, or it can be "year" for annual estimation of load. If it misses, the calculation is performed on the time period spanned by the streamflow data.

Value

A matrix with the estimated load with regression method for the different compounds. If the optional argument 'period' misses, the matrix will have only one row with the load estimation done throughout the time period spanned by streamflow data. Otherwise, if the period is specified, the load is estimated in the different months/year and the load values are returned in different rows, named with the month/year considered. The load is expressed in gram per period of estimation.

References

Quilbe', R., Rousseau, A. N., Duchemin, M., Poulin, A., Gangbazo, G., & Villeneuve, J. P. (2006). Selecting a calculation method to estimate sediment and nutrient loads in streams: Application to the Beaurivage River (Quebec, Canada). Journal of Hydrology, 326(1-4), 295-310. https://doi.org/10.1016/j.jhydrol.2005.11.008

Examples

Run this code

# NOT RUN {
data("flow.data1","conc.data1")
union<-db.union(flow.data1, conc.data1)
CQregression(union,1)
reg<-rating(union, 2)
reg.month<-rating(union, 2, "month")
reg.year<-rating(union, 2, "year")
# }