event_separation: Detect Rainfall Events with Specified Separation Time (Dry...

In mevr: Fitting the Metastatistical Extreme Value Distribution MEVD

Detect Rainfall Events with Specified Separation Time (Dry Time) and Minimum Rain Threshold

Description

This function identifies rainfall events in time-series data based on a minimum rainfall threshold and a separation time between events. It uses morphological operations to detect and separate events, and it can ignore short events that do not meet a minimum duration criterion. Based on Marra et al. 2024: "Predicting Extreme Sub-Hourly Precipitation Intensification Based on Temperature Shifts“. Hydrology and Earth System Sciences 28, Nr. 2: 375–89. \Sexpr[results=rd]{tools:::Rd_expr_doi("10.5194/hess-28-375-2024")}.

Usage

event_separation(
  data,
  separation_period = 360,
  time_resolution = 10,
  ignore_event_duration = 30,
  min_rain = 0.1
)

Arguments

data	A data frame with two columns, observation date and values. 'groupvar': A time variable (e.g., POSIXct or Date) indicating the timestamp of each observation. 'val': A numeric of the observed rainfall values.
separation_period	Numeric. The minimum dry period (in minutes) that separates two distinct rainfall events. Default is 360 minutes (6 hours). Values less than min_rain count as zero. Note that this should depend on local climatological conditions.
time_resolution	Numeric. The time step resolution of the data (in minutes). For example, if observations are recorded every 10 minutes, 'time_resolution' should be 10. Default is 10 minutes.
ignore_event_duration	Numeric. The minimum duration (in minutes) of rainfall required for an event to be considered valid. Events shorter than this will be ignored. Default is 30 minutes.
min_rain	Numeric. The minimum rainfall value (e.g., in mm) required for a data point to be considered part of an event. Default is 0.1.

Details

The function works by first dilating the rainfall data using a structuring element defined by the separation time, followed by erosion to remove noise. It then detects events based on the separation time and the minimum rainfall threshold, returning the time intervals where events occurred. Events that span more than two years or are shorter than the 'ignore_event_duration' are ignored.

Value

A list with the following elements:

• 'data': A modified version of the input data with an additional column 'is_event', indicating whether each time step is part of a detected event (1 for event, 0 for no event). This is useful in the next step ordinary_events.

• 'fromto': A tibble with two columns, 'from' and 'to', representing the indices of the start ('from') and end ('to') of each detected rainfall event.

• 'ts_res': The time resolution used in the function (in minutes).

Examples

## Not run: 
# Example data
mock_data <- data.frame(
  groupvar = seq.POSIXt(Sys.time(),
                        by = "10 min",
                        length.out = 100),
  val = c(rep(0, 10), runif(20, 0, 1), rep(0, 70))
)

# Detect rainfall events with a separation of 6 hours
# and a minimum rain threshold of 0.1
result <- event_separation(mock_data,
                           separation_period = 360,
                           time_resolution = 10,
                           min_rain = 0.1)
print(result$fromto)  # View detected events

## End(Not run)

mevr documentation built on Oct. 13, 2026, 5:07 p.m.