Pool {UKFE} | R Documentation |
Create pooling group
Description
Function to develop a pooling group based on catchment descriptors
Usage
Pool(
CDs = NULL,
AREA,
SAAR,
FARL,
FPEXT,
N = 500,
exclude = NULL,
iug = FALSE,
UrbMax = 0.03,
DeUrb = FALSE
)
Arguments
CDs |
catchment descriptors derived from either GetCDs or CDsXML |
AREA |
catchment area in km2 |
SAAR |
catchment standard average annual rainfall (1961-1990) in mm |
FARL |
catchment flood attenuation from reservoirs & lakes |
FPEXT |
catchment floodplain extent. The proportion of the catchment that is estimated to be inundated by a 100-year flood |
N |
minimum Number of total gauged record years for the pooling group |
exclude |
sites to exclude from the pooling group. Either a single site reference or a vector of site references (numeric) |
iug |
iug stands for 'include urban gauge' - which refers to a gauged subject site if it's > UrbMax. It's a logical argument with default of FALSE. TRUE will over-ride the default and add the closest site in catchment descriptor space (should be the gauge of interest) to the pooling group if it has URBEXT2000 >= UrbMax |
UrbMax |
Maximum URBEXT2000 level with a default of 0.03. Any catchment with URBEXT2000 above this level will be excluded from the pooling group |
DeUrb |
logical argument with a default of FALSE. If true, the Lcv and LSkew of any site in the pooling group with URBEXT2000 > 0.03 will be de-urbanised |
Details
A pooling group is created from a CDs object, derived from GetCDs or CDsXML, or specifically with the catchment descriptors (see arguments). To change the default pooling group, one or more sites can be excluded using the 'exclude' option, which requires either a site reference or multiple site references in a vector. If this is done, the site with the next lowest similarity distance measure is added to the group (until the total number of years is at least N). Sites with URBEXT2000 (urban extent) > 0.03 are excluded by default and this can be adjusted with UrbMax. If a gauged assessment is required and the site of interest is > UrbMax it can be included by setting iug = TRUE. De-urbanise the Lcv and Lskew (L-moment ratios) for sites with URBEXT2000 > UrbMax by setting DeUrb = TRUE. If the user has more data available for a particular site within the pooling group, the Lcv and Lskew for the site can be updated after the group has been finalised. An example of doing so is provided below. The pooling method is outlined in Science Report: SC050050 - Improving the FEH statistical procedures for flood frequency estimation.
Value
A data.frame of the pooling group with site reference row names and 24 columns, each providing catchment & gauge details for the sites in the pooling group.
Author(s)
Anthony Hammond
Examples
#Get some catchment descriptors
CDs.73005 <- GetCDs(73005)
#Set up a pooling group object called Pool.73005 excluding sites 79005 & 71011.
#Then print the group to the console
Pool.73005 <- Pool(CDs.73005, exclude = c(79005, 71011))
Pool.73005
#Form a pooling group, called PoolGroup, with the catchment descriptors specifically
PoolGroup <- Pool(AREA = 1000, SAAR = 800, FARL = 1, FPEXT = 0.01)
#Form a pooling group using an urban catchment which is intended for enhanced
#single site estimation - by including it in the group.
CDs.39001 <- GetCDs(39001)
Pool.39001 <- Pool(CDs.39001, iug = TRUE, DeUrb = TRUE)
#Change the Lcv and LSkew of the top site in the pooling group to 0.19 & 0.18,
#respectively.
PoolUpdate <- LRatioChange(Pool.39001, SiteID = 39001, 0.19, 0.18)