CAMELS_bootjack: Jackknife after bootstrap for all CAMELS sites

In gumboot: Bootstrap Analyses of Sampling Uncertainty in Goodness-of-Fit Statistics

Jackknife after bootstrap for all CAMELS sites

Description

Hydrologic model simulations can be produced using input-response data from the 671 catchments in the CAMELS dataset (Catchment Attributes and MEteorology for Large-sample Studies). Newman et al. (2015) and Addor et al. (2017) provide details on the hydrometeorological and physiographical characteristics of the CAMELS catchments. The CAMELS catchments are those with minimal human disturbance (i.e., minimal land use changes or disturbances, minimal water withdrawals), and are hence almost exclusively smaller, headwater-type catchments (median basin size of 336 km^2^). The CAMELS data used for the large-domain model simulations are publicly available at the National Center for Atmospheric Research at https://ral.ucar.edu/solutions/products/camels.

Usage

CAMELS_bootjack(
  CAMELS_sites = NULL,
  NetCDF_file = NULL,
  sim_var = "kge",
  GOF_stat = c("NSE", "KGE"),
  nSample = 1000,
  waterYearMonth = 10,
  startYear = NULL,
  endYear = NULL,
  minDays = 100,
  minYears = 10,
  seed = NULL,
  bootYearFile = NULL,
  quiet = FALSE
)

Arguments

CAMELS_sites	Required. Data frame of CAMELS sites. Must contain a field called hcdn_site. The data frame hcdn_conus_sites will work. You can subset this data frame if you want to use fewer sites.
NetCDF_file	Required. NetCDF file containing CAMELS modelled and gauged flows.
sim_var	Required. Name of variable containing simulated flows in NetCDF.
GOF_stat	Required. Name(s) of simulation goodness of fit statistic(s) to be calculated. Currently both NSE and KGE are supported.
nSample	Required. Number of samples for bootstrapping.
waterYearMonth	Required. Month of beginning of water year. Default is 10 (October). If the calendar year is required, set waterYearMonth = 13.
startYear	Optional. First year of data to be used. If NULL then not used.
endYear	Optional. Last year of data to be used. If NULL then not used.
minDays	Required. Minimum number of days per year with valid (i.e. greater than 0) flows. Default is 100.
minYears	Required. Minimum number years to be used. Default is 10.
seed	Optional. If NULL (the default) then no seed is specified for the random number generator used for the bootstrapping. If a value is specified then the bootstrapping will always use the same set of pseudo-random numbers.
bootYearFile	Optional. If NULL (the default) the years used for the bootstrapping are neither output nor input. If a file is specified, and it it does not already exist, then the bootstrap years will be written to a .csv file as a table with the dimensions of years x nSample. If a file is specified, and it _does_ exist, then the years will be read in, and used for the bootstrapping.
quiet	Optional. If FALSE (the default) a progress bar is displayed. If TRUE, it is not.

Value

Returns a data frame containing the following variables:

CAMELS_site

CAMELS site number

lat

CAMELS site latitude

lon

CAMELS site longitude

GOF_stat

Goodness of fit statistics (i.e. NSE or KGE)

seJack

standard error of jacknife

seBoot

standard error of bootstrap

p05, p50, p95

the 5th, 50th and 95th percentiles of the estimates

score

the jackknife score

biasJack

the bias of the jackknife

biasBoot

the bias of the bootstrap

seJab

the standard error of the jackknife after bootstrap

Author(s)

Martyn Clark and Kevin Shook

References

N. Addor, A. Newman, M. Mizukami, and M. P. Clark, 2017. Catchment attributes for large-sample studies. Boulder, CO: UCAR/NCAR. \Sexpr[results=rd]{tools:::Rd_expr_doi("10.5065/D6G73C3Q")}

Addor, N., Newman, A. J., Mizukami, N. and Clark, M. P.: The CAMELS data set: catchment attributes and meteorology for large-sample studies, Hydrol. Earth Syst. Sci., 21, 5293–5313, \Sexpr[results=rd]{tools:::Rd_expr_doi("10.5194/hess-21-5293-2017")}, 2017.

See Also

read_CAMELS

Examples

## Not run: 
camels <- CAMELS_bootjack(CAMELS_sites = sites, NetCDF_file = "CAMELS_flow.nc")

## End(Not run)

gumboot documentation built on Oct. 18, 2023, 5:08 p.m.