CAMELS_bootjack function

Jackknife after bootstrap for all CAMELS sites

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.

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.

Returns

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

Examples

## Not run:

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

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.

See Also

read_CAMELS

Author(s)

Martyn Clark and Kevin Shook