Skip to contents

Analyzes seasonal response changes using cumulative precipitation-runoff (R-Q) differences, identifying breakpoints between dry and wet periods.

Usage

rq_slopes(
  Q,
  R,
  hdays,
  start = 15,
  end = 183,
  bp = mean(c(start, end)),
  intercept = TRUE
)

Arguments

Q

Numeric vector of streamflow values L³/T

R

Numeric vector of precipitation values L/T

hdays

Numeric vector of hydrological day indices (1-365/366)

start

Integer first day to analyze (default: 15)

end

Integer last day to analyze (default: 183)

bp

Initial guess for breakpoint day (default: midpoint between start/end)

intercept

Logical indicating whether to estimate intercepts (TRUE) or fix at 0 (FALSE)

Value

Named numeric vector containing:

  • bp: Breakpoint day

  • bp_strength: Magnitude of slope change

  • slp_dry: Dry period slope

  • b_dry: Dry period intercept

  • slp_wet: Wet period slope

  • b_wet: Wet period intercept