CAMELS_bootjack: 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
)

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

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.

Author

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. 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, tools:::Rd_expr_doi("10.5194/hess-21-5293-2017"), 2017.

Examples

Run this code

if (FALSE) {
camels <- CAMELS_bootjack(CAMELS_sites = sites, NetCDF_file = "CAMELS_flow.nc")
}

Run the code above in your browser using DataLab