comp_GSIMdurfreq: Compute intermittency statistics for GSIM gauging stations

In NaiaraLopezRojo/globalIRmap: What the Package Does (One Line, Title Case)

Usage

comp_GSIMdurfreq(
  path_mo,
  path_sea,
  in_gaugep,
  maxgap,
  mdurthresh = 1,
  windowsize = 100,
  fullwindow = FALSE,
  monthsel = NULL,
  verbose = FALSE
)

Arguments

path_mo	(character) file path to a GSIM-formatted monthly streamflow time series table (e.g., ""C:/globalIRmap/data/GSIM/GSIM_indices/TIMESERIES/monthly/MA_0000010.mon")
path_sea	(character) file path to a GSIM-formatted seasonal streamflow time series table (e.g., ""C:/globalIRmap/data/GSIM/GSIM_indices/TIMESERIES/seasonal/MA_0000010.seas")
in_gaugep	(table; data.frame or data.table) of gauges' hydro-environmental attributes (including mean monthly temperature data)
maxgap	(integer) maximum number of days with missing data beyond which a year is not used in the computation of statistics
mdurthresh	(numeric) threshold of mean annual number of zero-flow days beyond which to classify gauge as intermittent.
windowsize	(integer) window size to check for zero-flow days.
fullwindow	(logical) whether years for which the window is truncated (e.g., beginning and end of time series) are taken in account in moving window analysis.
monthsel	(integer vector) selected months to compute the statistics over
verbose	whether to print input path

Value

One row data.table with 110 columns:

• gsim_no - (char) unique identifier for the gauge

• firstYear - (num) first year on full record

• lastYear - (num) last year on full record

• totalYears - (int) total number of years on full record
  For three subsets of the time series post-1800, post-1961, and post-1971 (e.g., suffix "mDur_o1800"):

  • firstYear_kept - (num) first year on record with < maxgap missing days

  • lastYear_kept - (num) first year on record with < maxgap missing days

  • totalYears_kept - (int) total number of years with < maxgap missing days

  • totaldays - (num) total number of days with discharge data

  • sumDur - (int) total number of days with discharge = 0

  • mDur - (num) mean number of days/year with discharge = 0

  • intermittent - (factor): binary flow intermittence class. 1: non-perennial (if mDur >= 1, i.e., if gauging station recorded zero-flow for at least one day per year on average); 0: perennial. '

  • movinginter - (logical): whether there is at least one zero-flow day in every windowsize-year (e.g., 20-year) moving window across the record. (with at least one day of flow between periods)

  • monthly statistics - (numeric) long-term mean monthly number of zero-flow days (column name example: 'Jan_mdur_o1800')

  • winteronlyir - Indicates whether gauge is only non-perennial during winter months. if intermittent == 1 AND the average annual number of zero-flow days during warm months < 1. Warm months are those with mean monthly catchment air temperature >= 10 (WorldClim v2; Fick and Hijmans 2017). 0: either perennial, or non-perennial outside of winter months. Only including years with <= maxgap missing daily discharge values.

Determine general characteristics of the whole time series and of the subset of years that have less than a given threshold of missing data as well as intermittency statistics. The intermittency statistics can be computed for a subset of months of the year (e.g. only winter months)

NaiaraLopezRojo/globalIRmap documentation built on Dec. 17, 2021, 5:19 a.m.