Introduction

The Growth Curve for your Site of Interest is estimated using the pooled or Enhanced Single Site (ESS) site method, as presented in detail by Kjeldsen et al., 2008.

If you have at-site AM data, the ESS procedure can be applied to give additional weight in the pooling procedure to the moments of the at-site growth curve. WINFAP defaults to a recommendation of using the ESS procedure if at-site data meet the pooling criteria (except that pertaining to the URBEXT2000). If the data do not meet the criteria, you can still select to ‘use at-site data’ for the Growth Curve thereby entering your station into the pooling group with an additional weighting.

The 'Pooling & QMED Analysis' dashboard provides all the tools required to develop and assess the pooling group, as shown below.

image.png

Defining a pooling group

The estimated statistical properties of your Site of Interest's growth curve are estimated in the pooling procedure by taking a weighted average of the L-moments of the distributions of Annual Maxima (AM) data from a pooling group of hydrologically similar gauging stations.

When starting the Pooled & QMED Analysis method, a number of options are presented to you for use in developing the pooling group.

image.png

You should consider the following settings:

  • Pooling group size (years of AM data). This defaults to 500 years.
  • Only include stations suitable for pooling. This is selected by default.
  • Use 'small catchment' recommendations as default. This is advised for catchments below 25km2 and is available for catchments below or equal to 40km2. If selected the small catchments distance measure (SDM) will be used to create the pooling group and only 1 station will selected by default for use in the Donor adjustment of QMED.
  • Include stations with URBEXT2000 < 'X'. This defaults to 0.03 which is categorised as 'essentially rural'.
  • The pooling group is, by default, selected from stations in the default NRFA Peak Flow Dataset, identified in the Peak Flow Registry. Additional stations can be added through the identification of a supplementary folder.

These settings can be reviewed from the main dashboard via the information icon next to Pooling Group Details.

The pooling group consists of multiple gauging stations across the UK which are judged as hydrologically similar to the target catchment. Similarity is identified using the FEH similarity distance measure (SDM). By default, for stations with a catchment area above 25km2 hydrological similarity is based on work completed by Kjeldsen et al., 2008 and is assessed with regards to the catchment descriptors: AREA, SAAR, FARL and FPEXT. For stations with catchment areas below 40km2 the option is provided to use the equations developed as part of the Small Catchments project (Stewart et. al., 2019), which is based on AREA and SAAR only with a lower weighting associated with AREA. This option is selected as default for catchments less than 25km2.

The standard FEH similarity distance measure (SDM) is estimated using the following equation as described in Kjeldsen et al., 2008:

SDM = \sqrt{3.2 \bigg[ {{Ln(Area_i) -Ln(Area_j)} \over {1.28}} \bigg]^2 + 0.5 \bigg[ {{Ln(SAAR_i) -Ln(SAAR_j)} \over {0.37}} \bigg]^2 + 0.1 \bigg[ {{FARL_i -FARL_j} \over {0.05}} \bigg]^2 + 0.2 \bigg[ {{FPEXT_i -FPEXT_j} \over {0.04}} \bigg]^2 }

The following equation can be used to estimate the SDM within small catchments and is described in Stewart et. al., 2019:

SDM = \sqrt{\bigg[ {{Ln(Area_i) -Ln(Area_j)} \over {1.264}} \bigg]^2 + \bigg[ {{Ln(SAAR_i) -Ln(SAAR_j)} \over {0.349}} \bigg]^2 }

The pooling group members are presented in an expandable list in the sidebar.

Pooling Group Criteria

In addition to options presented to the user during initiation of pooled analysis, stations must satisfy the following criteria to be included in the pooling group:

  • There are at least 5 years of AM data. It is preferable to have at least 8 years of AM data.
  • Catchment descriptors (AREA, SAAR, FARL, FPEXT) are known.
  • Catchment area ≥ 0.5 km².

You are encouraged to review the pooling group to ensure that the best stations are selected to predict the Growth Curve for your Site of Interest. Consideration should be given to the quality of the peak flow data for stations in the pooling group. The dashboard provides comprehensive options for reviewing and refining the pooling group:

Data for a site can be edited via the Station Manager or Single Site Analysis Dashboard. To ensure consistency between the Site of Interest and pooling group members, it is recommended that any edits to the Site of Interest data are completed prior to selecting 'Pooled & QMED Analysis'. Whilst any site edits implemented (e.g. annual maxima edits) are used in the calculations for existing pooling groups, the membership of any existing pooling groups is not updated. It is therefore also recommended that any existing pooling groups are regenerated following any edits to pooling group member site catchment descriptors (e.g. SAAR) that are used in selecting stations for the pooling group. Following these recommendations will prevent any potential inconsistent results from changes to the default pooling membership.

At-site data

Details of the at-site data are provided. For a gauged site this will provide information on the URBEXT2000, the suitability criteria and the record length available at the site.

Where the site has an URBEXT2000 greater than the threshold set for stations to be included in the pooling group a warning will be presented. If you are using the ESS method, where at-site data is used within the pooling method, you should make an assessment of whether it is appropriate to include the station within the pooling group. Given the improvement in estimates associated with the use of local observed data it will often be the case that it will be appropriate to retain the station within the pooling group.

image.png

Urbanisation

For the Site of Interest and those in the pooling group, you have the option to de-urbanise the L-moments for use in the Growth Curve development procedure. See the section on Urbanisation for more details. The default settings for this can be over-ridden but warnings are provided where selections are inconsistent.

Pooling Group Table

You can inspect various parameters of your pooling group stations in the Pooling Group Table. The data are tabulated in two tabs: 'AM Data' and 'Catchment Descriptors'.

AM Data tab

  • Distance (SDM): this is the similarity distance measure between your Site of Interest and the pooling-group station calculated using the respective descriptors. Smaller values are associated with stations that are more hydrologically similar to your Site of Interest.
  • Years of data: the number of years of AM data in the record.
  • QMED AM: the observed QMED calculated from the AM data.
  • L-CV Observed, L-SKEW Observed: sample L-moments of the AM series.
  • L-CV Deurbanised, L-SKEW Deurbanised: deurbanised sample L-moments of the AM series.
  • Discordancy: a measure of how different the station's L-moment ratios are from the pooled L-moment ratios. A higher value indicates a ‘more different’ dataset.

A colour and symbol coding system has been adopted to assist in the identification of:

  • Stations with short records (default GREY).
  • Stations deemed to be not suitable for use in a pooling group in the NFRA Peak Flow Data (default YELLOW).
  • Stations deemed to be not suitable for use in a pooling group or for the estimation of QMED in the NFRA Peak Flow Data (default MAGENTA).
  • Stations with a significant discordancy measure (default RED).
  • Stations that have been added through the supplementary folder are identified by a +.
  • Stations which have been edited in the project are identified by a *.

You can select to use the observed or deurbanised L-moments. The values used in the analysis are shown in bold. See urbanisation procedures for more details on default and recommended settings.

Stations that have been added through the supplementary folder (+) or edited in the project (*) are also identified using the same symbols in the sidebar lists. Rejected stations are listed below the main table in pale grey.

Catchment Descriptors tab

  • Distance (SDM): this is the similarity distance measure between your Site of Interest and the pooling-group station calculated using the respective descriptors in the similarity equation. Smaller values are associated with stations that are more hydrologically similar to your Site of Interest.
  • AREA: catchment area (km²).
  • SAAR: standard annual average rainfall (mm), a measure of catchment wetness.
  • FPEXT: an index of floodplain coverage within the catchment.
  • FARL: an index for assessing the relative important of flood attenuation due to reservoirs and lakes.
  • URBEXT2000: an index of urbanisation within the catchment based on landuse in 2000.
  • BFIHOST19: an index of the permeability of the catchment. BFIHOST19 is not used in either similarity measure but can be useful in identifying catchments that may have very different flood mechanisms to those within the Site of Interest's catchment.
  • A dynamic link to the NRFA station's page.

Rejected stations are listed below the main table in pale grey.

Modify Pooling Group

'Add Station' and 'Reject Station' allows you to add and reject stations from the pooling group. Rejected stations are listed below the main table in pale grey. Note that where at-site data has been rejected this will not be presented within the list of rejected stations.

Add Station

On selection of 'Add Station' a list of the appropriate stations in the station manager is displayed.

image.png

You can change which stations are displayed using the Station Selection Criteria, and order by the 'similarity distance' or 'Station number'.

All stations previously rejected are displayed in grey. Stations that have been modified or from the supplementary data folder are also highlighted.

When a station is added that has a higher SDM than other stations available, all stations with an SDM lower are considered to have been rejected.

Reject Station

To remove a station from the pooling group, highlight the station in the pooling table and select the 'Reject Station' button.

Review Pooling Group

The 'Review Pooling Group' section provides the following graphical options:

  • Station Record Parameters: a UK map of the gauging stations; the distribution of the members of the pooling group with respect to flood behaviour; individual flood growth curves for each member of the pooling group.
  • 3D L-Moment Graph: a three-dimensional plot of the L-CV, L-skewness and L-kurtosis for each of the pooling group members and your Site of Interest (if the Site of Interest has observed data).
  • Catchment Descriptor Graphs: the ‘location’ of your Site of Interest in the data space of a variety of catchment descriptor variables. Your Site of Interest is shown as a cross (X) and the members of the pooling group as dots on a histogram.
  • All Analysis Graphs: a combination of all graphs from the Station Record Parameters and Catchment Descriptor Graphs.

Exploratory Data Analysis

The two Exploratory Data Analysis methods, ‘Goodness of Fit’ and ‘Heterogeneity’, are intended to assist in your choices about pooling group membership. They should not be used in isolation to refine pooling groups; expert judgement is essential. Please note that both these methods use random sampling hence, for a given site and settings, values will vary.

The background to the diagnostic tools is presented in Flood Estimation Handbook Volume 3, Chapter 16.

Distribution Selection

The growth curve can be estimated using a number of different distributions.

Following publication of Kjeldsen et. al., 2017, the Kappa 3 distribution, a 4 parameter distribution, has been added to the list of distributions which can be used in the growth curve estimation process. The GLO, GEV, and GPA distributions are specific applications of the Kappa distribution where ξ is the location parameter, α the scale parameter and κ the shape parameter. The 4th order parameter, 'h', is set to -1, 0, and 1 respectively for each of these distributions. Based on the research, the Kappa 3 distribution uses a fixed 'h' of -0.4; the optimum value found for UK catchments.

By default the Generalised Logistic, Generalised Extreme Value and Kappa 3 distributions are selected. Kjeldsen et al., 2010 showed that the Generalised Logistic distribution is, in the main, the preferred distribution in the UK, however other distributions should also be considered.

The Goodness of Fit tools can aid in the selection of the appropriate distribution.