--- title: "FEH Method (FEH Rainfall Runoff Method) | Technical Reference" slug: "feh-method-feh-rainfall-runoff-method" description: "Find the data, theory and guidance for the FEH method (FEH rainfall runoff method)." tags: ["1D river modelling", "rainfall-runoff", "error messages and diagnostics"] updated: 2024-10-24T20:19:03Z published: 2024-10-24T20:19:03Z canonical: "help.floodmodeller.com/feh-method-feh-rainfall-runoff-method" stale: true --- > ## Documentation Index > Fetch the complete documentation index at: https://help.floodmodeller.com/llms.txt > Use this file to discover all available pages before exploring further. # FEH Method (FEH Rainfall Runoff Method) The FEH Rainfall Runoff Method Boundary (FEHBDY) derives an inflow hydrograph from a catchment or sub-catchment. The hydrograph then becomes a boundary condition equivalent to a Flow Time Boundary. Alternatively, the rainfall component of the FEH Boundary can be used as a direct rainfall boundary (equivalent to a REBDY), when used in conjunction with a Lateral Inflow unit, by selecting the 'Hyetograph' Boundary Type Option. ## Theory and Guidance ### Introduction The FEH Boundary is a rainfall-runoff model based on procedures described in the [Flood Estimation Handbook (1999)](/v1/docs/feh-method-feh-rainfall-runoff-method#references). It is essential that users of the FEH Boundary are fully conversant with the appropriate parts of the Flood Estimation Handbook - neither the software nor this documentation fully reproduce the guidance or definition of the methods presented in the Handbook. In particular, Volumes 1 and 4 of the Handbook include important guidance on the choice of methods for flood frequency estimates and on restrictions in the applicability of the methods. For example, design flow hydrographs produced by the FEH rainfall-runoff method can be highly inaccurate if parameter estimation is based solely on catchment descriptors. For up-to-date information on the Flood Estimation Handbook please visit the [Centre for Ecology and Hydrology website](https://fehweb.ceh.ac.uk/). The FEH Boundary is designed to be flexible so that the components of the rainfall-runoff model (such as baseflow, percentage runoff, unit hydrograph) can either be estimated from the catchment descriptors or can be directly provided based on event data (usually the preferred option). In addition, the generated flow hydrographs can be scaled to match the flood peak obtained by the FEH statistical methods. The rainfall-runoff model is based on unit hydrograph (UH) theory. For the required location, effective rainfall is transformed into a runoff hydrograph by convolution with the unit hydrograph. Please refer to Volume 4 of the FEH or standard hydrology texts for general details of unit hydrograph theory. In the absence of flood event data, the parameters of the FEH rainfall-runoff model can be estimated from catchment descriptors. These catchment descriptors and the parameters of the rainfall frequency model can be obtained from the [Centre for Ecology and Hydrology website](https://fehweb.ceh.ac.uk/) and/or the FEH CD-ROM (if appropriate). A facility exists within Flood Modeller to import all relevant catchment descriptors and rainfall parameters from the standard FEH CD-ROM CSV (comma separated variable) output file. The hydrograph produced by the FEH Boundary can be checked prior to use by either clicking on the 'Calculated hydrograph' tab in the FEHBDY properties dialog in the Flood Modeller interface, or by running a 'Boundary mode' simulation. A summary of the input data and rainfall-runoff component values is written to an ASCII file with the extension '.zzb'. The hydrograph and mass balance is tabulated in an ASCII file with the extension 'zzh'. When a boundary mode run has been undertaken, the hydrographs at all hydrological boundaries can be viewed (and compared) using the Flood Modeller Interface time series feature. You can generate a datafile with: - only one FEHBDY node - many FEHBDY nodes each representing versions of the same catchment (eg different return periods or model component options) - many FEHBDY nodes each representing different sub-catchments, optionally linked to additional Flood Modeller units such as Muskingum-Cunge routing units, hydrodynamic river sections, reservoirs/ponds, structures, junctions etc. The FEHBDY unit contains a number of options to control the method of calculation. Within the Flood Modeller ASCII data file, these options are referred to in the form xxFLAG (e.g. PRFLAG, UHFLAG etc). This format is also used below to describe the options (the Flood Modeller Interface does not require knowledge of these 'flags'). ### Equations The FEH Boundary calculates total runoff as a function of the following components: P = precipitation; either calculated for a specific return period (ERFLAG = 'FEHER') or 'observed' (ERFLAG = 'OBSER'). This rainfall depth will use a rainfall profile either based on winter data (RPFLAG = 'WINRP'), summer data (RPFLAG = 'SUMRP'), or observed/derived values (RPFLAG = 'OBSRP'). UH = unit hydrograph; either calculated (UHFLAG = 'FSRUH') or observed/derived (UHFLAG = 'OBSUH') PR = percentage runoff; either calculated (PRFLAG ='FEHPR') or observed/derived (PRFLAG = 'OBSPR') BF = baseflow, either calculated (BFFLAG= 'F16BF') or observed/derived (BFFLAG='OBSBF') Where the above components are 'calculated' then the following equations are used (new terms are either defined after introduction or are listed in the subsequent Data section). #### Rainfall If using the standard FEH rainfall profile (ERFLAG=FEHER): The design rainfall depth is calculated from the depth-duration-frequency (DDF) relationships defined in [FEH Volume 2 Section 2.3](/v1/docs/feh-method-feh-rainfall-runoff-method#references). The rainfall depth is multiplied by the areal reduction factor to produce areal rainfall. #### Unit Hydrograph parameters The unit hydrograph can either be derived from local data and specified as a time series of UH ordinates (usually preferred), or the standard FEH triangular UH can be used. Where the standard FEH UH is used (UHFLAG=FSRUH) then Tp can either be specified (TPFLAG=OBSTP) or calculated (TPFLAG = FEHTP) from:       $T_{p}(0) = 4.27DPSBAR - 0.35PROPWET^{-0.8}DPLBAR^{0.54}(1+URBEXT) - 5.77$ $T_{p}(t) = T_{p}(0) + t/2$ $Q_{p} = \frac{2.2 CAREA}{T_{p}(t)}$ $T_{B} = 2.52 T_{p}(t)$ where $T_{p}(0)$= time to peak of the instantaneous unit hydrograph (hrs) $T_{p}(t)$= time to peak of t-hour unit hydrograph (hrs) (generally referred to as Tp) $Q_{p}$ = peak flow of t-hour unit hydrograph (m3/s) $t$ = data interval (hrs) $T_{B}$ = time base of UH (hrs) and DPSBAR, DPLBAR, PROPWET, URBEXT and CAREA are defined in the [Data section](/v1/docs/feh-method-feh-rainfall-runoff-method#data) and are usually obtained from the [Centre for Ecology and Hydrology website](https://fehweb.ceh.ac.uk/) and/or the FEH CD-ROM (if appropriate). #### Total Percentage Runoff If using the FEH method and SPR to calculate percent runoff (PRFLAG = FEHPR):     $PR = PR_{RURAL}(1-0.615 URBEXT) + 70(0.615 URBEXT)$ , if $PR_{RURAL} \gt 70$ $PR = PR_{RURAL}$ , otherwise. where $PR_{RURAL} = SPR + DPRCWI + DPRRAIN$ $DPRCWI = 0.25 (CWI - 125)$ $DPRRAIN = 0.45 (P - 40)^{0.7}$ , if $P \gt 40$mm $DPRRAIN = 0$ , otherwise $P$ = precipitation (mm) $DPR$ = Dynamic Percentage Runoff and SPR and CWI are as defined in the [Data section](/v1/docs/feh-method-feh-rainfall-runoff-method#data), with SPR best obtained through analysis of local event data, but in the absence of records can be obtained from the [Centre for Ecology and Hydrology website](https://fehweb.ceh.ac.uk/) and/or the FEH CD-ROM (if appropriate), and CWI may be entered directly (CWFLAG = 'OBSCW') or, for the design case, the program can calculate it from a digitised version of Figure 3.7 in [FEH Volume 4](/v1/docs/feh-method-feh-rainfall-runoff-method#references). #### Baseflow If using the standard FSSR16/FEH calculation method for baseflow (BFFLAG = F16BF):     $BF = CAREA (33 (CWI - 125) + 3.0 SAAR + 5.5) \times 10^{-5}$ where CAREA, CWI and SAAR are as defined in the [Data section](/v1/docs/feh-method-feh-rainfall-runoff-method#data). Note that several fields within the datafile are left blank (or contain zeros). This is to facilitate easier portability between the FSRBDY and FEHBDY units, i.e. fields describing FSRBDY parameters which are not used by FEHBDY or have no FEHBDY equivalent are ignored and therefore do not need to be overwritten when converting FSRBDY boundary units to FEHBDY boundaries. #### Probable Maximum Flood (PMF) The FEH Rainfall Runoff Method Boundary allows you to obtain a PMF (Probable Maximum Flood) estimation for the catchment using the FEH guidelines. This involves the introduction of the Probable Maximum Precipitation (PMP) design storm hyetograph, including modifications to the CWI and snowmelt contributions, in addition to further modifications to other constituent factors, namely the unit hydrograph and percentage runoff. These modifications are relative to those for the T-year flood estimation. The design storm period remains unchanged. A choice of summer and winter PMPs are available. For further reading and a detailed description of the method, see the [FEH Vol. IV, Chapter 4](/v1/docs/feh-method-feh-rainfall-runoff-method#references) and [Hough and Hollis (2006)](/v1/docs/feh-method-feh-rainfall-runoff-method#references). **Design PMP - EMP** Values for EMP2h and EMP24h, the all-year point estimated maximum precipitations (EMPs) of 2-hour and 24-hour duration for the catchment, must be supplied (in mm). Additionally, the 25-day EMP (EMP25d) is required for design storms of long duration (D > 19.1hrs). Examples of appropriate values can be found in the [FEH Vol. IV, Figures 4.1-4.3](/v1/docs/feh-method-feh-rainfall-runoff-method#references). The hyetograph is then derived from the supplied EMPs, the EM rainfall and seasonal variation factors as given in the FEH, and the areal reduction factor (ARF) at each time interval Dt, 3Dt, ... D. For time intervals of more than 0.5hr, the ARF is calculated by Flood Modeller as a function of STAREA and STDUR from a digitised form of Figure 3.4 in the FEH Volume 4. For time intervals between 0.25hr and 0.5hr, the ARF is interpolated from values in this table. For time intervals of less than 15 minutes, the ARF is extrapolated from this table and may be inaccurate. **Design PMP - Snowmelt** For Winter PMPs, storm event and antecedent period snowmelt contributions may be added, which also may affect the CWI. Snowmelt is invoked by supplying positive values for the snowmelt rates (in mm/day). Alternatively, snowmelt rates can be calculated automatically as: - Melt rate (storm) = 100-yr melt rate + melt rate due to storm rain energy - Antecedent melt rate = 100-yr melt rate + melt rate due to antecedent rain energy The 100-year melt rate and melt rate energy are determined based on the following:      100-yr melt rate = M100:M5 growth rate × 5-yr melt rate Melt rate energy = $0.0125 R_{depth} \Delta T \frac{24}{D}$ where $R_{depth}$ = rain depth (mm) and is that for either the storm event or antecedent period, as appropriate $D$ = storm duration (for storm energy); = antecedent duration (= 2 × storm duration) for antecedent rain energy $\Delta T$ = temperature above freezing (= temperature in ˚C) and the 5-year melt rate and M100:M5 growth rate are obtained by a regression equation based on Northing and Altitude (ALTBAR), and using extreme value analysis based on the nearest available weather station, respectively. When using snowmelt, a value for S100 (100-year snow depth water equivalent, in mm) must also be set. The snowmelt contribution to the precipitation at during each time interval Dt is therefore Dt * melt rate subject to enough snow from the total S100 being available. For Summer PMPs, no snowmelt contribution will be added, irrespective of whether a snowmelt rate is supplied. **Design PMP - CWI** For the Catchment Wetness Index (CWI) for a PMP, the equation used is:     $CWI = 125 + EMa(0.5)^{D/24}$ with the estimated antecedent rainfall given by $EMa = 0.5 \left(ARF_{5D}EMP_{5D} - ARF_{D}EMP_{D} \right)$  where $ARF_{nD}$ is the Areal Reduction Factor related to a time of duration $nD$ $EMP_{nD}$ is the seasonal point Estimated Maximum Precipitation related to a time of duration $nD$  For a Winter PMP, if a snowmelt option is chosen, the CWI is the summation of the CWI due to antecedent rainfall and the CWI due to antecedent melt, and is further incremented by:       $SMa(0.5)^{D/24}$ where $SMa$is the amount of snow melting over the antecedent period $2D$  If the amount of snow remaining from S100 after the contribution from the storm period has been removed is insufficient, then the necessary adjustments are made according to FEH Vol IV, Section 4.3.4. The calculation for baseflow is dependent on CWI, therefore this value will also be different from that for the T-year flood estimation. **Frozen Ground SPR** The calculated percentage runoff method remains the same as that for the FEH T-year flood estimation FEH with the exception that a Winter PMF is assumed to contain the frozen ground adjustment that the SPR cannot be lower than 53%. If this frozen ground adjustment is not required for a Winter PMF, then the user should manually calculate a percentage runoff and select an Observed PR. **Unit hydrograph time-to-peak** The calculated time-to-peak of the PMF instantaneous unit hydrograph is adjusted by a factor of 0.67 from that of the standard FSR unit hydrograph, thus:     $T_{p}(\Delta t)_{PMF} = 0.67 T_{p}(0) + \frac{\Delta t }{2}$  If an Observed time-to-peak is selected, it is assumed that this factor is contained within the value and no adjustment is made. ### General The user is referred to the Flood Estimation Handbook, for full details of the derivation, application and limitations of the rainfall-runoff method. The FEH recommends that where possible, certain of the rainfall-runoff model parameter values obtained from catchment descriptors are replaced with, or revised using, values from observed data. General recommendations on parameter estimation are: - If suitable event data are available then the FEH recommends calculating Tpfrom such data, for example     $T_{p}(0) = 0.879 LAG^{0.951}$ where $LAG$ is the time between the centroid of rainfall and centroid of hydrograph peak - The data interval, t, should be small enough that the flood hydrograph is well defined. A data interval of 10%-20% of the Tp(0) is usually suitable. This should then be rounded off to a convenient value such as 0.25, 0.5, 1 or 2 hours. It is usually important to find the critical values of parameters such as storm duration and storm profile which give the largest flood peak (or largest flood volume in some cases). The sensitivity of the model results to these parameters is described in the FEH and some conclusions are: - For deriving design events the FEH recommends use of the 75% winter profile or 50% summer profiles. On predominantly rural catchments (URBEXT < 0.125) the 75% winter profile is recommended. On more urban catchments (0.125 ≤ URBEXT ≤ 0.5) the 50% summer profile is recommended. These two profiles are included as options in the FEHBDY unit and can be accessed by setting the options RPFLAG to 'WINRP' or 'SUMRP' respectively. - Simulated peak flows can be sensitive to D, the storm duration. For the design case, the recommended storm duration for critical flood peak is:    $D = \left( 1+ \frac{SAAR}{1000} \right) T_{p}$ With D taken as the nearest odd integer multiple of t. Note that with reservoired catchments or semi-distributed catchment modelling it will be necessary to consider a range of durations. The raw value of D, i.e. not adjusted for t, is output under "Critical Storm Duration" in the *.zzb (data summary) file after an unsteady or boundary mode simulation. NB The value of Tp used is the time to peak of the t-hour unit hydrograph. In the FEH procedure the flow return period is not necessarily equal to the storm (rainfall) return period. - In rural or only moderately urbanised catchments (URBEXT < 0.125) the FEH procedure provides a relationship between the storm and flow return periods (e.g. a 140-year storm return period is said to produce a 100-year return period flow). - For urban catchments (0.125 ≤ URBEXT ≤ 0.5), the storm and flow return periods are assumed to be equal. The FEHBDY unit will use this relationship to calculate storm return period if the storm return period in the datafile is set to zero. For highly urbanised catchments (URBEXT > 0.5), FEH procedures are not recommended, and, in such cases, a non-fatal warning message will be displayed. The point rainfall calculation (ERFLAG='FEHER') is derived from a Depth-Duration-Frequency (DDF) model (FEH, Volume 2, Chapter 2). For a valid DDF model, the storm return period must be greater than one year, whereas one greater than 10000 years is not recommended. Furthermore, the parameters were based on observations of return periods of no greater than 1000 years. When importing the parameters for the DDF method from the FEH CD-ROM output, the values used are c,¼,f, rather than c(1km),¼,f(1km). The Flood Modeller implementation of the DDF model is based on a sliding duration rainfall (FEH Volume 2 Section 2.5), rather than a fixed duration. For user-input unit hydrograph ordinates, the user may specify the units of these values from a selection of four commonly used units. These are any combination of 1cm or 1mm of effective rainfall depth covering either the whole catchment area or per 100km2. These units will also be used in the various unit hydrograph output options, either textually (*.zzb and *.zzh files) and graphically, via the user interface. The units specified also apply for output when deriving the unit hydrograph by the FEH (FSR) method. Whilst the FEHBDY is essentially a UK-specific lumped single event rainfall-runoff model, the underlying unit hydrograph approach is applicable to other countries if the various components of the model are defined using local data. ## Data | Field in Data Entry Form | Description | Name in Datafile | | --- | --- | --- | | ALTBAR | Mean catchment altitude (m). Only used in derivation of event and antecedent period snowmelt rates. Range 0.0 to 9968.0 | ALTBAR | | Antecedent Melt Rate | Rate of snowmelt (mm/day) during the antecedent period. Used in CWI calculations. This is converted to the appropriate units and added to the runoff hydrograph. The facility to include snowmelt is not part of the standard rainfall frequency method. | AMRATE | | ARF | Areal reduction factor to relate point rainfall to areal rainfall. If set to zero then the areal reduction factor is calculated by Flood Modeller as a function of STAREA and STDUR from a digitised form of Figure 3.4 in the FEH Volume 4. If STAREA is set to zero then CAREA is used. Range 0.0 to 1.0 | arf | | Baseflow Value | User-defined baseflow (m3/s) | BF | | Minimum flow | Adjusted baseflow (m3/s). Usually zero. Forces the resultant hydrograph to be not less than the specified flow. This may be required for program stability | BFADJS | | Baseflow Calculation Method | Baseflow flag. FEH (bfflag='F16BF'): use FSSR16/FEH calculation method for baseflow, or Observed (bfflag='OBSBF'): use user-defined baseflow value. | BFFLAG | | Simulation Type | Flow values to use during simulation, can be one of: 'base flow' (bfonly) 'peak flow' (pfonly) 'full hydrograph' (; default) For example, if 'base flow' is selected then the baseflow contribution is used for the boundary flow value for the whole simulation. | bfonly | | c | DDF model parameter c | c | | CALCRATES | Override of snowmelt rates provided with calculated values. Override on (CALCRATES = 'CALCTRUE'): snowmelt rates during the storm event (SMRATE) and antecedent period (AMRATE) are recalculated from catchment parameters at simulation runtime. Override off (CALCRATES = 'CALCFALSE'): use user-input values. A blank value is equivalent to CALCFALSE. | CALCRATES | | Tp Coefficient | Calibration factor for Tp (eg CALIB=1.1 increases Tp by 10%). Default value is 1.0 | CALIB | | Catchment Area | Contributing catchment area (km2) | CAREA | | Catchment Wetness Index Flag | Catchment wetness index calculation flag. FEH (cwflag='FSRCW'): use user-input value. For a PMF calculation, selecting FEH uses the PMF CWI calculation (cwflag='PMFCW') | CWFLAG | | Observed CWI | Catchment wetness index (mm) | CWI | | Country Flag | Country location flag. GB for Great Britain or NI for Northern Ireland. | CYFLAG | | Storm Duration | Storm duration (hrs) | D | | d1 | DDF model parameter d1 | d1 | | d2 | DDF model parameter d2 | d2 | | d3 | DDF model parameter d3 | d3 | | DPLBAR | Mean drainage path length (km) | DPLBAR | | DPSBAR | Mean drainage path slope (m/km) | DPSBAR | | e | DDF model parameter e | e | | Easting | Easting coordinate of catchment outflow (Only used in derivation of event and antecedent period snowmelt rates.) | easting | | Em-2h | Estimated maximum 2-hour rainfall (mm). Used in PMP calculation only. | EMP2h | | Em-24h | Estimated maximum 24-hour rainfall (mm). Used in PMP calculation only. | EMP24h | | Em-25d | Estimated maximum 25-day rainfall (mm). Used in PMP calculation only. | EMP25d | | Event Rainfall Flag | Event rainfall flag. FEH (erflag='FEHER'): use FEH method to calculate rainfall depth: PMF (erflag='PMFER'): use the FEH method to calculate a PMP depth or Observed(erflag='OBSER'): use user-input rainfall depth | ERFLAG | | f | DDF model parameter | f | | Use Equal Return Periods | When this field is set to (force= 'FORCE'), the program stops if Tf is not equal to Ts for the SUMMER rainfall profile option. | FORCE | | Boundary Type Flag | Boundary mode flag: 'HYDROGRAPH' (the default if blank) denotes a generated hydrograph (i.e. the unit behaves as a QTBDY); 'HYETOGRAPH' indicates that it behaves as a REBDY-type unit, applying the Rainfall profile to a Rainfall-Only boundary (the latter must be used in conjunction with a lateral inflow node) | hymode | | Node Label | Node label identifier | Label | | Northing | Northing coordinate of catchment outflow (Only used in derivation of event and antecedent period snowmelt rates.) | northing | | n/a | Number of rainfall profile values (maximum of 4600) | nrp | | n/a | Number of unit hydrograph values (maximum of 400) | nuh | | Rainfall depth | Event rainfall precipitation (mm) | P | | PR | Percentage runoff (range 0 - 100) | PR | | Runoff Flag | Percentage runoff calculation flag: FEH (prflag='FEHPR'): use FEH method and SPR to calculate percent runoff; or Observed (prflag='OBSPR'): use user-input PR value. | PRFLAG | | Percentage Runoff Option | Option for variable or fixed percentage run-off. Keyword either 'FIXED' or 'VARIABLE'. If VARIABLE then the method used is to scale the specified PR value (or that calculated from a specified SPR) during the storm based on the method described in FEH Volume 4 Section A.5.2 | PRVAR | | PROPWET | Proportion of time catchment soil moisture deficit (SMD) was below 6mm during the period 1961-1990. (Range 0.0 - 1.0) | PROPWET | | Storm Profile: Rainfall | User-defined rainfall profile starting at tstart with data interval t (mm) | rp | | Storm Profile Flag | Rainfall profile option: - observed - rpflag='OBSRP' - 75% winter or Winter PMP -rpflag= 'WINRP' - 50% summer or Summer PMP - rpflag = 'SUMRP' 75% winter and 50% summer profiles are taken from Figure 3.5 in FEH Volume 4 PMP profiles are used if a PMF event is selected, and are derived using Section 4.3.2 in FEH Volume 4 | RPFLAG | | Use refined rainfall profile | If checked (refine rp=1), uses a finer discretisation of the standard FSR/FEH rainfall profiles - can prevent 'blocky'rainfall profiles for long storm/short data interval events | refine_rp | | SAAR | Standard annual average rainfall (mm) | SAAR | | Hydrograph Scaling Method:Scale | Hydrograph scaling method: 'FULL' (scaling ='FULL'; default) scales the whole hydrograph; Quick Runoff (scaling='RUNOFF'): scaling only applies to the quick runoff component of the hydrograph. Used in conjunction with SCFLAG and scfact. | SCALING | | Method | Hydrograph scaling option. By a factor of ('scflag='SCALE' [default]) - applies the specified factor; To fit peak of (scflag='PEAK') - fits the hydrograph peak to the specified value. | SCFLAG | | Hydrograph Scaling Option | If SCFLAG=PEAK then all hydrograph ordinates are scaled (by a constant value) to achieve a peak flow of scfact If SCFLAG=SCALE then all hydrograph ordinates are multiplied by scfact (default is 1.0) | scfact | | S100 | 100-year snow depth water equivalent (used in PMP calculations only) | S100 | | Snow Melt Rate (Storm) | Rate of snowmelt (mm/day) during the storm event period. This is converted to the appropriate units and added to the runoff hydrograph. The facility to include snowmelt is not part of the standard rainfall frequency method. | SNRATE | | SPR | Standard percentage runoff (range 0 - 100) | SPR | | Storm Area | Storm area (km2) when STAREA is entered as zero then STAREA is set to CAREA | STAREA | | Storm Duration | Storm duration (hrs). Should be the nearest odd integer multiple of t | STDUR | | Data Interval | Time interval used for rainfall profile and unit hydrograph (hrs) | t | | Temperature | Average temperature (oC). Default value 10, minimum 0. Only used in derivation of event and antecedent period snowmelt rates. | TEMP | | Time Delay | Optional delay time (hrs), e.g. if tdelay=2hrs then the hydrograph will start 2hrs after the start time of the simulation | tdelay | | Rainfall Return Period | Storm return period (years). If set to 0.0 then calculated from Figure 3.2 of FEH Volume 4. Must be > 1 year | Ts | | Flood Return Period | Flood return period (years) | Tf | | Tp value | User-defined time to peak of t-hour unit hydrograph (hrs) | Tp | | Tp Calculation Method | Unit hydrograph Tp flag. FEH (tpflag='FEH16P'): use FEH method to calculate unit hydrograph time-to-peak, or Observed (tpflag='OBSTP'): use user-defined value. | TPFLAG | | Storm Return Period | Storm return period (years). If set to 0.0 then calculated from Figure 3.2 of FEH Volume 4. Must be > 1 year | Ts | | UH Ordinate | Unit hydrograph ordinates with data interval t (see units) | uh | | UH Time | unit hydrograph time (hrs) | | | TB Scale Factor | Unit hydrograph time base (TB) adjustment factor. New time base is uhbadj multiplied by the calculated TB value. The UH peak is adjusted accordingly to preserve its unit property. | uhbadj | | n/a | Scaling factor for the unit hydrograph ordinates. Only used if the 'units' keyword is not recognised. If zero or blank then units of m3/s/mm are used. (Not accessible from the user interface.) | uhfctr | | Unit Hydrograph Flag | Unit hydrograph flag. FEH (uhflag='FSRUH'): use standard FSR/FEH triangular hydrograph, or Observed (uhflag='OBSUH'): use user-defined unit hydrograph. | UHFLAG | | Units | Units of the unit hydrograph ordinates. Keyword can be one of: - 'cm100k': 1cm depth over 100km2 (units = m3/s/cm/100km2) - 'cmarea': 1cm depth over catchment area (units = m3/s/cm) - 'mm100k': 1mm depth over 100km2 (units = m3/s/mm/100km2) - 'mmarea': 1mm depth over catchment area (units = m3/s/mm) - The conventional units are m3/s /cm/100km2 | units | | URBEXT | Extent of urban/suburban land cover (range 0.0 - 1.0) | URBEXT | | n/a | Elevation (m AD). Not used in FEHBDY computations | z | ## Datafile Format ### General Line 1 : Keyword 'FEHBDY' [comment] Line 2 : Label Line 3 : z, easting, northing Line 4 : tdelay, t, bfonly, SCFLAG, scfact, hymode Line 5 : CYFLAG Line 6 : CAREA, , , , URBEXT, ALTBAR Line 7 : STAREA, STDUR, , SNRATE, SN100, TEMP, AMRATE, CALCRATES Line 8 : SAAR, , , , FORCE Line 9 : ERFLAG Line 10 : P (if ERFLAG = 'OBSER') Line 10 : Tf, Ts, arf, c, d1, d2, d3, e, f (if ERFLAG = 'FEHER') Line 10 : , , arf (if ERFLAG = 'PMFER') Line 11 : CWFLAG Line 12 : CWI (if CWFLAG = 'OBSCW') Line 12 : (if CWFLAG = 'FSRCW' or 'PMFCW') Line 13 : PRFLAG, PRVAR Line 14 : PR (if PRFLAG = 'OBSPR') Line 14 : SPR (if PRFLAG = 'FEHPR') Line 15 : TPFLAG Line 16 : CALIB, Tp (if TPFLAG = 'OBSTP') Line 16 : CALIB, DPLBAR, DPSBAR, PROPWET (if TPFLAG = 'FEHTP') Line 17 : BFFLAG Line 18 : BFADJS, BF (if BFFLAG = 'OBSBF') Line 18 : BFADJS (if BFFLAG = 'F16BF') Line 19 : UHFLAG Line 20 : 0, units, uhfctr (if UHFLAG = 'FSRUH') Line 20 : nuh, units, uhfctr (if UHFLAG = 'OBSUH') Line 21 to Line 20+nuh : uh (if UHFLAG = 'OBSUH') Line 21+nuh : RPFLAG Line 22+nuh : 0, em2h, em24h, em25d (if ERFLAG = 'PMFER') Line 22+nuh : 0 (if RPFLAG = 'WINRP', ERFLAG¹ 'PMFER') Line 22+nuh : 0 (if RPFLAG = 'SUMRP', ERFLAG¹ 'PMFER') Line 22+nuh : nrp (if RPFLAG = 'OBSRP') Line 23+nuh to Line 22+nuh+nrp : rp (if RPFLAG = 'OBSRP') where items in square brackets are optional, and other parameters are as defined in the [data section](/v1/docs/feh-method-feh-rainfall-runoff-method#data). ### Example ```none FEHBDY M029 0.000 385115 242872 0.000 0.250 SCALEFACT 1.000HYDROGRAPH 1.000 FULL OVERRIDE GB 16.1700 0.000 0.000 0.000 0.005 400 0.000 5.250 0.000 168.471 200.000 10 96.011 CALCTRUE 1175.000 0.000 0.000 0.000 PMFER 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 FSRCW 0.000 FEHPR FIXED 42.760 FEHTP 1.000 5.170 178.100 0.520 F16BF 0.000 FSRUH 0 mmarea WINPMP 1 0 158.000 290.000 500.000 ``` ## References Flood Estimation Handbook (1999)Centre for Ecology and Hydrology, Wallingford, Oxfordshire (formerly the Institute of Hydrology). December 1999*.* ISBN 0 948540 94 X. Vol 1: *Overview*. Vol 2: *Rainfall Frequency Estimation* Vol 3: *Statistical Procedures for Flood Frequency Estimation* Vol 4: *Restatement and Application of the Flood Studies Report Rainfall Runoff Method* Vol 5: *Catchment Descriptors Flood Estimation Handbook at the Centre for Ecology and Hydrology website*. An FEH CD-ROM is also available from the Centre for Ecology and Hydrology containing data for flood frequency estimation covering the entire United Kingdom.Hough, M. N. and Hollis, D (2006)*Rare snowmelt estimation in the United Kingdom*, Meteorological Applications, Science and Technology for Weather and Climate, Volume 5, issue 2.