thresh: Rainfall threshold for event-based landslide occurrences

View source: R/thresh.R

threshR Documentation

Rainfall threshold for event-based landslide occurrences

Description

This function uses statistical approaches to calculate landslide rainfall thresholds and to quantify their uncertainties.

Usage

thresh(
  Re,
  D,
  method = c("LS", "QR", "NLS"),
  prob.thresh = 0.05,
  trans.log10 = TRUE,
  bootstrapping = TRUE,
  R = 1000,
  use.bootMedian = FALSE,
  use.normal = FALSE,
  nls.pw = 10,
  nls.tw = 10,
  nls.method = c("tw", "pw"),
  nls.bounds = list(lower = c(t = 0, alpha = 0, gamma = 0), upper = c(t = 500, alpha =
    100, gamma = 10)),
  seed = 123,
  cores = 1
)

Arguments

Re

vector containing the rain event variable, e.g. cumulated event rainfall (in mm) or intensity (mm/h)

D

vector containing the duration of the rainfall events

method

method to compute threshold. Either 'LS' for least square, 'QR' for quantile regression, or 'NLS' for non-linear least squares Method. Default: 'LS'

prob.thresh

exceedance probability level. Default: 0.05 (5 percent)

trans.log10

log-transformation of input vectors Re and D. Default: TRUE

bootstrapping

If TRUE bootstrapping is performed. Default: TRUE

R

the number of bootstrap replicates, see boot::boot() for more information. Default: 1000

use.bootMedian

Threshold is delineated from median of bootstrapping result (i.e., regression is overwritten). Default: FALSE

use.normal

Use normal approximation of residuals. This method is described in literature. Default: FALSE

nls.pw

For NLS method, size of point-wise ED-pair filtering. Default: 10

nls.tw

For NLS method, size of time-wise ED-pair filtering. Default: 10

nls.method

Define method for NLS ED-pair filtering. Default: 'tw'

nls.bounds

bounds for finding NLS starting parameter. Default: list(lower = c(t = 0, alpha = 0, gamma = 0), upper = c(t = 500, alpha = 100, gamma = 10))

seed

replicable bootstrapping. Default: 123

cores

If cores > 1 than parallelisation for bootstrapping is initialized via future back-end. Default: 1

Value

list object containing threshold metrics and settings (see description).

Note

The use of method = "NLS" is not recommended.

  • Brunetti, M. T., Peruccacci, S., Rossi, M., Luciani, S., Valigi, D., & Guzzetti, F. (2010). Rainfall thresholds for the possible occurrence of landslides in Italy. Natural Hazards and Earth System Sciences, 10(3), 447.

  • Peruccacci, S., Brunetti, M. T., Luciani, S., Vennari, C., & Guzzetti, F. (2012). Lithological and seasonal control on rainfall thresholds for the possible initiation of landslides in central Italy. Geomorphology, 139, 79-90.

  • Rossi, M., Luciani, S., Valigi, D., Kirschbaum, D., Brunetti, M. T., Peruccacci, S., & Guzzetti, F. (2017). Statistical approaches for the definition of landslide rainfall thresholds and their uncertainty using rain gauge and satellite data. Geomorphology, 285, 16-27.

  • Guzzetti, F., Peruccacci, S., Rossi, M., & Stark, C. P. (2007). Rainfall thresholds for the initiation of landslides in central and southern Europe. Meteorology and atmospheric physics, 98(3-4), 239-267.


raff-k/RainSlide documentation built on July 19, 2023, 6:26 a.m.