calcREIC: Robust Event Identification Criteria (REIC)

In hydroEvents: Extract Event Statistics in Hydrologic Time Series

Robust Event Identification Criteria (REIC)

Description

Calculates the REIC value for a set of rainfall–runoff events given rainfall, quickflow, and event data.

Usage

calcREIC(rainfall, quickflow, events, n_rainfall, area)

Arguments

rainfall	Numeric vector. Rainfall data (mm) for the study period.
quickflow	Numeric vector. Quickflow (direct runoff) data (ML) corresponding to the same period as rainfall.
events	Data frame. Identified paired rainfall–runoff events, including rainfall and runoff events' start and end.
n_rainfall	Integer. The number of rainfall events detected.
area	Numeric. Catchment area in square kilometres.

Details

The Robust Event Identification Criteria (REIC) calculates the REIC value for a set of paired rainfall–runoff events based on two criteria outlined in Mohammadpour Khoie et al. (2025). It helps identify rainfall–runoff events that are both consistent and reliable (i.e., with fewer runoff coefficients greater than one).

The REIC metric is easily transferable and can be applied to various sets of rainfall–runoff events identified using different parameter values, helping users select the best-performing configurations.

Value

A numeric vector of REIC values.

References

Mohammadpour Khoie, M.M., Guo, D., & Wasko, C. (2025). Improving the consistency of hydrologic event identification. Environmental Modelling & Software, 106521.

See Also

pairEvents, plotPairs

Examples

## Not run: 
# Example usage
rain_events <- eventPOT(dataLoch, threshold = 1, min.diff = 1)
flow_events <- eventMaxima(dataBassRiver, delta.y = 200, delta.x = 1, 
threshold = 0)
paired_events <- pairEvents(rain_events, flow_events, lag = 5, type = 1)
quickflow <- baseflowA(dataBassRiver, alpha = 0.925, passes = 3)$bf

REIC_values <- calcREIC(dataLoch, quickflow, paired_events, 
n_rainfall = nrow(rain_events), area = 3012)

## End(Not run)

hydroEvents documentation built on Dec. 11, 2025, 5:06 p.m.