A DEPTH – DURATION – FREQUENCY ANALYSIS FOR SHORT - DURATION RAINFALL EVENTS IN E NGLAND AND W ALES

.


INTRODUCTION
Estimates of the magnitude of rainfall events of a given duration with an expected annual exceedance probability p, are an important component of current methods of flood frequency estimation, used in the design and assessment of flood defence schemes, bridges and reservoir spillways, as well as urban drainage systems.Rainfall frequency estimates are also a key input to mapping studies of the risk of surface water flooding.The estimates can be obtained from depth-duration frequency (DDF) models, in which the relationship between the rainfall depth, event duration and event rarity is integrated in a unique framework.In a DDF model, it is required that frequency curves for different durations do not cross, meaning that the rainfall depth that is exceeded with probability p should increase monotonically with increasing event duration.The probability p is typically expressed as a return period T, with p=1/T, as events larger than those corresponding to the quantile that is expected to be exceeded with probability p should happen on average every T years.
DDF models, which are often referred to as Intensity Duration Frequency (IDF) models, can then serve two purposes: to estimate the rainfall depth of a hypothetical event with a given duration and rarity, and to assess the rarity of a storm event with known rainfall depth and duration.Svensson and Jones (2010) give an overview of different DDF models used in several countries, showing the large array of possible approaches to rainfall frequency estimation.Many of the countries included in the review use some form of index rainfall approach combined with regional estimation of growth curves for different durations, although some countries were reported to use the linear regression approach.The idea behind the latter approach is to fit a statistical distribution separately to the single series of block maxima of different accumulation periods and then to fit regression models across the different durations or frequencies, so that increasing rainfall depths are estimated for increasing durations given a certain frequency.See Koutsoyiannis et al. (1998) for a discussion of the mathematical formulation of the relationship between the duration and frequency of rainfall events, and a general discussion of DDF modelling.Although the relationship between rainfall depths and frequencies has been studied for several decades, there is still much interest in identifying methods to derive DDF curves (e.g.Overeem et al., 2008) and in the actual derivation of DDF curves to be used at different sites of interest (e.g.Jiang and Tung, 2013).
One interesting finding of the review in Svensson and Jones (2010) is that, in several contries, different models are used depending on the duration and rarity of the rainfall events of interest.The need for different models for different durations and frequencies stems from the diffculty of developing models that can provide reliable results across several rainfall durations and frequencies.One country where several DDF models are currently in use is the UK: the main models are presented below and are the main focus of this study.
In the UK, the most widely used DDF models are those presented in the Volume II of the Flood Studies Report (FSR, Natural Environment Research Council, 1975) and in Volume 2 of the Flood Estimation Handbook (FEH99, Faulkner, 1999), which mostly superseded the FSR methods.Recently, a new model (FEH13, Stewart et al., 2013) has been developed, with the specific aim of overcoming the issues encountered when the original FEH99 model is used to estimate rare events.Estimates from the FEH13 model have only been available to practitioners since November 2015, and have therefore not yet been widely used in practice.Furthermore, the performance of the FEH13 model for short duration events (i.e.under 1 hour) is still being assessed, since most of the model evaluation focussed on the estimation of the frequency of long-duration events.Considering that the FEH13 model aimed to improve rainfall frequency estimates for rare events with durations longer than 1 hour, it is not yet clear how it will perform for the more frequent events of very short duration which are of interest in this study.
The FSR and FEH99 DDF models are based on an index-rainfall approach and were developed with the scope of providing nationwide rainfall frequency estimates.The FEH99 method was calibrated on a larger network of stations with longer records than the FSR method and, unlike the FSR method, incorporated a spatial model in which data from nearby stations were used for rainfall frequency estimation at a given location.On the other hand, the FEH99 method was calibrated using data with an accumulation period of at least 1 hour while, in the development of the FSR method, some data with an accumulation period of 1 minute were also used.Compared to the FSR method, the FEH99 method has been found to give much larger estimates of rainfall depth for the very long return periods required for reservoir safety assessment (MacDonald and Scott, 2001;Babtie Group, 2000).As a result, the FSR and FEH99 methods are both still used, but for different cases that depend on the duration and rarity of the design event to be estimated.As Svensson and Jones (2010) report, the FSR method can be used to estimate return periods of rainfall events with accumulation periods between 1 minute and 25 days and return periods longer than 1000 years, and is recommended for the estimation of rainfall depths associated with return periods up to 10,000 years The FEH99 method provides estimates of rainfall accumulations between 1 hour and 8 days, with return periods shorter than 1,000 years and, although rainfall frequencies up to return periods of 10,000 years can technically be estimated, their use is not recommended.The newly developed FEH13 might replace the FSR and the FEH99 as the recommended model to use to estimate the magnitude of very rare events, but the official guidelines have not yet been amended.The FEH99 method can also be extended to estimate the frequency of rainfall events with accumulation periods shorter than 1 hour, although, as no sub-hourly data were used in the calibration of the method, extrapolation to durations below 30 minutes is strongly discouraged.The coexistence of the FEH99 and FSR methods results in uncertainty when estimates are needed for sub-hourly rainfall events.These cases go beyond the range of reliable estimates for the FSR, a relatively old model that was calibrated on fairly short records with very limited sub-hourly data, and beyond the intended use of the FEH99, a more complex and structured model that was calibrated using a dense network of stations but no sub-hourly data at all.Small catchments (i.e.smaller than 25km 2 ) and plot-sized areas are expected to be particularly vulnerable to short, intense cloudbursts, due to their short response times.As Faulkner et al. (2012) emphasise, reliable estimates of sub-hourly design rainfalls are therefore needed to allow credible flow and hydrograph estimates for the smallest catchments using rainfall-runoff techniques.The suggestions in Faulkner et al. (2012) motivated the second phase of the Environment Agency's Estimating Flood Peaks and Hydrographs for Small Catchments project.The project aims to improve the estimation of flood frequencies in small catchments and encompasses, among other things, an assessment of the most appropriate methods to estimate the frequency of very short duration rainfall, which this study is concerned with.A novel at-site DDF modelling strategy is discussed and an application of the proposed model is presented using data series available at selected sites that give a reasonable geographical coverage of England and Wales, for which relatively long records of subhourly rainfall are available.The proposed model does not follow the traditional approaches and uses instead the data across all durations to fit a unique model.Rainfall frequency curves estimated with the proposed method are compared to those estimated with the FSR and FEH99 DDF models, and to empirical return level estimates.
The stations and datasets used in the study are introduced in the next section.Subsequently, a unified GEV model is proposed and its performance for the stations under study is discussed.The performance of the unified GEV, FSR and FEH99 models for short-duration rainfall frequency estimation are compared in the section Comparisons of the unified GEV results to current methods.The final section of the paper contains the conclusions and final remarks.

DATA
From the large number of tipping bucket rain gauges managed by the Environment Agency (EA) and Natural Resources Wales (NRW) and providing sub-hourly rainfall data, a subset with sufficiently long records was identified that could allow for good spatial coverage of the area.Sub-hourly data for the rainfall stations are available as time of tip (ToT) at some sites series and as aggregated 15-minute accumulation series at other sites.In the selection of stations to be included in this scoping study, priority was given to those for which ToT data were available, to allow very short durations to be investigated.It appears that long ToT series are more readily available in some regions (the English Midlands and Wales), hence the final subset of stations included in the study is a compromise between the competing needs of having long series and maintaining a good coverage of England and Wales (E&W).In particular, the sites were chosen to be at least 35 km apart.The final selected stations are shown in Figure 1, which also shows the Standard-period Annual Average Rainfall (SAAR) for the years 1961-1990(Spackman, 1993)).The shortest series in the dataset is 15 years long; the longest two are each 46 years long.A total of nine ToT series and ten 15-minute series are included in the study dataset.The analysis was performed on the annual and seasonal maxima of the different accumulations, with two six-month seasons included in the study.The final dataset was compiled from the ToT and 15-minute series, following two slightly different workflows as outlined below.
For all series, a month was considered complete if at least 75% of the data in the month were nonmissing.Finally, the annual and seasonal maxima series were constructed from the monthly maxima series.A year or season was considered complete if no more than one monthly record within that year or season was incomplete.Approximately 89% of the station-seasons have at least 99% of valid data points, 98% of the station-seasons have at least 90% of valid data points and in just one instance is the percentage of valid data points in a season lower than 80% (summer rainfall series of 1995 at Victoria Park, which has a total of 79.3% valid data points).Overall, for all stations, for the series across all years and seasons more than 99% of the total number of data points are recorded as valid, giving reasonable confidence in the quality of the available data and confidence that the maxima were captured.Annual maxima were extracted as the maximum single value recorded in each calendar year.Summer maxima were extracted as the maximum value recorded in the months from May to October inclusive.Winter maxima were extracted as the maximum value recorded in the months from November to April inclusive.The availability of the raw time of tip information for the tipping bucket stations allows for the extraction of series at a 1-minute resolution and additionally at coarser or even finer resolutions.
However, the level of precision that can be reached in high resolution series depends greatly on the tip volume of the instrument, a property that might change slightly in time (e.g.due to sediment collecting in the bucket) or more significantly over time (e.g. if the specific gauge used at a station is replaced by a different model).Furthermore, the tip volume might be different at different stations, thus creating inconsistencies in the precision across different stations.The discrete nature of the tipping bucket measurements is also the underlying reason why, in a number of months, the recorded 1-minute and 2-minute maxima have the same value, and why several annual and seasonal maxima are identical across a number of years.These issues are more common in the earlier years of the record, during which time the data were measured at a coarser resolution.The issues connected to systematic errors in tipping buckets are known (Molini et al., 2005, and references therein).In particular, lower intensities tend to be overestimated and higher intensities tend to be underestimated.Methods to quantify the systematic error of each station are beyond the scope of this study, and the data extracted from the original series are used in all subsequent analysis without further adjustment.The issues connected with the original data series should nevertheless be acknowledged as they can have an impact on the estimation procedures discussed in the section Results for the at-site analysis and in the comparisons discussed in the section Comparisons of the unified GEV results to current method.Due to differences in the underlying data collection methods, the series of maxima extracted from the ToT and the 15-minute series do not provide the same information for accumulations of 15 minutes or greater.The time of tip maxima are computed using a sliding window, so the 15-minute annual maximum value (for example) corresponds to the actual largest amount of rainfall recorded in any 15minute interval in the year.However, the maximum obtained from the 15-minute records instead corresponds to the maximum amount recorded in one predefined 15-minute interval, which is likely to be lower than the actual maximum amount of rainfall that could have been recorded in a 15-minute interval without a fixed start time.The true maximum rainfall is most likely to be under-recorded when its duration is the same as the fixed-duration recording unit, as the rainfall event is very unlikely to align neatly with the station clock.However, when longer durations are considered, the alignment between the rainfall event and the station clock is less important, as the depths of rainfall at the tail ends of the storm, which are difficult to capture exactly, become less and less important to the storm depth as a whole.
To adjust the maxima extracted from the 15-minute stations so that they are closer to the higher values that would be attained using sliding windows, correction factors were introduced.For each time of tip record, fixed-period (15-minute) annual and seasonal maxima were extracted for durations of 15, 30, 45, 60, 90 and 120 minutes.These series correspond to the maxima that would be obtained if the data for the time of tip stations were stored as 15-minute series (fixed window) rather than time of tip series (sliding window).The average ratio between the sliding window maxima and the fixed window maxima at each duration, shown in Table 1, is used as a sliding window correction factor for that duration.In the rest of this work, the maxima extracted from the 15-minute series are multiplied by the appropriate correction factor to give estimates of the equivalent sliding window maxima.Due to the different ranges of time resolution present in the two different data sources, two separate analyses are carried out: one which uses only the series extracted from the ToT stations and covers the range of The FSR, FEH99 and FEH13 DDF models build on a large set of available gauges and allow the estimation of frequency curves for a number of durations across the whole UK.In particular, the FEH99 and the FEH13 have complex spatial model components so that estimates for rainfall frequencies at one point are built incorporating information from nearby gauges.Such complex spatial structures are unattainable with the subset of stations available in this study.Given the exploratory scope of this work, a simpler model is proposed: the model allows the estimation of a station's DDF curves based solely on the data series available for that station; it does not have a component to include information from nearby stations.
The proposed model builds on extreme value theory (Coles, 2001), assuming that block (e.g. annual or seasonal) maxima follow a Generalised Extreme Value (GEV) distribution: X~GEV(ξ, α, κ) where X indicates the random variable that describes rainfall block maxima and ξ, α and κ are the location, scale and skewness parameters of the GEV distribution respectively.The cumulative distribution function of a GEV distributed random variable X~GEV(ξ, α, κ) is defined as: The set on which the variable X is defined, e.g. the values that might be observed in a sample from a population with underlying distribution X, is governed by the skewness parameter as follows: A DDF model for short-duration rainfall events in England and Wales The distribution is bounded for the case in which κ ≠ 0, with the lower and upper bound being a linear combination of the distribution parameters.The skewness parameter therefore defines whether an upper or lower bound for the values of X exists.
The quantile function for the GEV distribution, which is used to build frequency curves, is derived as: where F is the non-exceedance probability, corresponding to  = 1 − 1  ⁄ for the T-year event.The desired property of a DDF model is that the quantile functions for increasing durations of rainfall accumulation, D, do not cross.This means that, denoting by x(F,D) the rainfall depths of durations D associated to a certain non-exceedance probability F, for d0 < d1 one should have x(F,d0) < x(F,d1).
The proposed model uses the relationship between the GEV parameters shown in Equation ( 2) and stems from some empirical properties observed via visual explorations of the estimated parameters for the different durations at each station (see the section Results for the at-site analysis).The GEV distribution can be shown to be the asymptotic distribution of sample maxima (see Coles, 2001) and has often been used as an underlying distribution in the development of DDF models (among others Overeem et al., 2008, Jiang andTung, 2013).According to the godness of fit test presented in Kjeldsen and Prosdocimi (2015), the GEV distribution was deemed acceptable for a large majority of the series analysed in the study.When estimating frequency curves, it is expected that no upper limit should be attainable by the rainfall values at any duration, so the skewness parameter is constrained in the proposed model to be negative.The model development is presented below only for the case in which κ < 0, although similar ideas would apply for κ > 0. It is assumed that the skewness parameter κ is constant across all durations, while the location and scale parameters are dependent on the duration D: ξ() and α().Taking ℓ to be the lower bound of the distribution, and assuming this to be the same for all durations, the following relationship is obtained from the inequality in Equation (2): The quantile function shown in Equation (3) can then be updated to a quantile function   (), which depends on the event duration D via the location parameter ξ(): Provided that ξ() is monotonically increasing, the function   () is a monotonically increasing function of D, so that the estimated frequency curves give consistent results for increasing durations.
For the case of the British rain gauges under study, the following relationship is proposed to model the location as a function of the event duration, based on the observed properties of the location parameter for a GEV distribution fitted separately for each different duration across all stations (see Figure 2 in the next Section): which is an increasing function of D provided that its first derivative is positive: The scale function is determined by a combination of the lower bound (ℓ), the skewness parameter (κ) and the location function (ξ()) according to Equation (4).The proposed unified GEV model then requires the estimation of a total of six parameters (a, b, c, g, ℓ, ), a relatively parsimonious model which, given some constraints in the location function, allows for consistent frequency estimates for different durations.It is possible that an even simpler formulation could be used for Equation ( 6), but the suggested function originates from the methods discussed in Stewart et al. (2013) and seems to give reasonable results.
The proposed unified GEV model uses a different strategy to obtain consistent estimates for rainfall frquencies than many published works, which use approaches based on linear regression across estimates for the different durations.The unified GEV model presented in this paper instead seeks to fit a unique model to all series at once, so that all available information is used to estimate the DDF curves.The development of the model is inspired by some of the discussion in Stewart et al. (2013), on the development of the statistical framework used in the FEH13 model.
The basic novel idea behind the proposed model is to ensure that monotonic quantile functions are obtained by constraining some of the parameters of the rainfall distribution to have common properties across different durations.It is possible that for a different set of durations, or a new set of gauging stations that exhibit different properties, the assumptions of which common distributional properties are be shared across durations might be different.Furthermore, the functional relationship between the location and the duration shown in equation ( 6) might not be valid.Nevertheless the building blocks of the proposed model could be adapted to accommodate different data behaviours: the unified GEV is an addition to the possible modelling approaches used for at-site estimation of DDF curves.

RESULTS FOR THE AT-SITE ANALYSIS
For each station separately, the parameters of the unified GEV model (a, b, c, g, ℓ, ) are estimated via maximum likelihood, which ensures some optimal properties for parameter estimates (Coles, 2001).The unified GEV model is fitted to the data from all the ToT stations and to all the series with accumulations of at least 15-minutes, in two different fitting procedures.The location function, shown in Equation ( 6), and the relationship between the scale and other parameters, shown in Equation ( 4), are used in the two fitting procedures.Estimated rainfall depth-duration-frequency curves for the annual, winter and summer series for the ToT station Dowdeswell are shown in Figure 3, together with the block maxima extracted from the original series plotted using Gringorten plotting positions.The frequency curves seem to fit the data reasonably well.Due to the constraints in the model structure that ensures that the location function is monotonically increasing for increasing durations, the frequency curves computed from the formula in Equation ( 5) tend to fan out.A noticeable feature of the data is that the winter maxima tend to be much smaller than the summer maxima, which also appear to be the annual maxima.Results for the other ToT stations have similar properties to the ones shown in Figure 3   Figure 4 shows the estimated location and scale functions, together with the skewness and lower bound parameters of the unified GEV model, for annual data at all 19 stations, considering accumulations of 15 to 120 minutes.As in Figure 2, the original estimates for the GEV parameters obtained from an Lmoments estimation procedure fitted to each duration separately are also shown.Again, the fitted location functions seem to be mostly in agreement with the estimates obtained from the different durations, while more variability can be seen in the estimation of the scale function in the top right panel.In particular, the scale functions for Victoria Park and Otterbourne are very flat, as a result of the estimates for the skewness parameters at these stations being very close to zero.The estimated lower bounds for these two stations are also very small: -37.7 at Otterbourne and -124.6 at Victoria Park (censored in Figure 4).The fact that the skewness parameter for these stations is estimated to be very close to zero in the unified GEV model is likely to be connected to the fact that some series in these stations appear to have a finite upper bound (e.g.positive skewness) for some durations.In the unified GEV model, the skewness parameter is required to be negative and to be unique for all durations, so that the final estimate is a summary of the properties of all durations.If the behaviour of the series at a station differs across durations, the final estimates need to be a compromise between the different tendencies of each series.Nevertheless, the final fit of the estimated frequency curves compared to the annual maxima shown in Figure 5 seem to indicate that overall an acceptable fit is obtained for the series at Otterbourne.The estimated frequency curves shown in Figure 5 have similar properties to those shown in Figure 3 the curves have a tendency to fan out and the annual extremes appear to mostly driven by summer rather than winter events.
Seasonal differences are not explored further in this analysis, but the estimates obtained from the different stations could be employed in the future to develop correction factors to obtain seasonal estimates from sub-hourly annual estimates, similarly to Kjeldsen et al. (2006).The unified GEV proved to be a flexible and reliable modelling approach which could give reasonable estimates across different seasons.

COMPARISONS OF THE UNIFIED GEV RESULTS TO CURRENT METHODS
The estimated depths obtained with the methods currently in use (FSR and FEH99) and the proposed unified GEV, corresponding to some pre-specified frequencies, are compared to the empirical estimates obtained from the recorded data series at each station.Since reliable estimates of very rare events cannot be obtained from the relatively short records available (the median record length for the observed series is 24 years), the comparison is limited to the 2-, 5-and 10-year return periods.The empirical estimates are obtained as the median (50 th percentile), 80 th percentile and 90 th percentile of the recorded data series.For some series, the record length would be less than 2T years long when estimating the 10-year event: these empirical estimates might be less precise.The comparison is performed on every station for durations of at least 15-minutes, and the fitted unified GEV models shown in Figure 4 are used to estimate the rainfall depths.
Figure 6, Figure 7 and Figure 8 display the relative differences between the rainfall depths, as estimated with the different methods, and the empirical quantile corresponding to the specific frequencies for the A DDF model for short-duration rainfall events in England and Wales 2-, 5-and 10-year return periods respectively.For example, the left panel of Figure 6 shows, for each station and each duration, the value (R2FSR -R2Observed)/ R2Observed, where R2FSR and R2Observed indicate the estimated 2-year rainfall of the given duration at a station and the empirical 2-year event estimated from the observed data, respectively.The unified GEV model is the only method directly fitted to the observed data, which explains the much better performance of that model in comparison to the FSR and FEH99 models for the 2-year events.In particular, the FSR appears to give consistently positively biased estimates for the 2-year events (Figure 6), with lower variabilities in the error for longer durations.The FEH99 seems to perform well on average, although the results are more variable than the unified GEV.The results for the longer return periods show more variation, with the unified GEV performing slightly better in terms of the variability of the error.The unified GEV, an at-site model fitted directly to the observed data, performs quite well for most stations.Among the models currently used in the UK for rainfall frequency estimation, the FEH99 seems to give acceptable results, across all return periods, with an error variability comparable to the FSR estimate.These comparisons are based only on empirical estimates of events with a relatively short return period, and it is not clear how the different models differ for the estimation of rare events, for which no reliable empirical estimates can be obtained from the observed series.An assessment of the accuracy of the A DDF model for short-duration rainfall events in England and Wales different estimation methods for longer return periods would in fact require reliable information on the real frequency of short-duration rainfall events, which cannot be easily retrieved.The overall relative difference between the FEH99 and FSR, which were developed with the purpose of allowing DDF estimation for the whole UK, and the estimates obtained from the unified GEV model, estimated using only at-site data is investigated in Figure 9.The Figure shows, for a large range of return periods, the relative difference between the design events estimated by FSR and the unified GEV (left panel) and the difference between the design events estimated by FEH99 and the unified GEV (right panel), for all stations and all durations.The average relative differences across all stations for each duration are also shown.It should be noted that a large difference between the standard methods and the unified GEV estimates does not necessarily indicate poor performance of the standard methods: the unified GEV models are fitted to the recorded data series, which are at most 46 years long.It is therefore very likely that unified GEV estimates would be more accurate for shorter rather than longer return periods.
Nevertheless, what is visible in the plots is that the variability is much larger for the longer return periods for all durations.Furthermore, the FSR seems to give consistently larger results than the unified GEV for short return periods, but the difference between the two estimates become smaller for return periods longer than 10 years.For the 15-minute events the difference is more marked and the FSR seems to give much smaller estimates than the unified GEV for longer return periods.The difference between the FEH99 and the unified GEV results instead appear to increase for longer return periods, although for shorter return periods (up to 5 years) the difference in the two estimates is on average very small.At very long return periods, it appears that the average difference between the unified GEV and the FEH99 estimates is smaller for events of long duration.between the location and scale parameter and the event duration.The effect of the assumptions, enforced to ensure the consistency of the frequency curves, is that curves for different durations diverge more and more as return period increases.The model might therefore give extremely large rainfall depth estimates for very long return periods.The model is designed to be fitted to block maximum series of sub-hourly data at single stations and does not have a procedure to integrate information from nearby stations in the rainfall frequency estimation.The estimation of the model parameters was carried out by maximum likelihood estimation, a procedure that attains some optimal properties when applied to large samples.The final model parameter estimates are influenced by properties of the observed data series and issues might arise when the actual properties of the observed series do not match well with the properties that are assumed during model building.Nevertheless, the proposed model gives overall satisfactory results and fits the empirical data quite well, using a relatively small number of parameters.The new estimated frequency curves are compared to those obtained using the FSR and FEH99 methods currently employed in the UK.Although no sub-hourly data were used in the model calibration, the FEH99 method seems to give acceptable results for all of the sub-hourly durations under study, at least for the return periods for which reliable empirical rainfall frequencies can be estimated.The FSR estimates seem to overestimate the rainfall depths for short return periods, although the bias is less marked for longer return periods.In addition, the difference between the FEH99 and the FSR estimates becomes larger for rarer events.However, the comparisons could only be carried out on a small set of stations, and a more in-depth analysis would be needed to give a robust indication of the behaviour of the different models.Potentially, it could be useful to develop a full DDF model for short duration rainfall events at a national scale, in which information from different stations could be used in a unique framework.The relative scarcity of long and precise records of subhourly data would be the major obstacle to overcome in the potential development of a DDF model for the whole UK.Most of the available ToT records are fairly short and most are located in only a few areas of the UK.Due to the nature of tipping buckets, the measurement of very short duration events is likely to be biased, especially in less recent years, which would undermine the quality of any estimation procedure.

Figure 1 .
Figure 1.Location of the 19 stations included in the study.The record length of the annual maxima series is indicated in the location of each station: red numbers indicate ToT stations, blue numbers indicate 15-minute stations.The underlying colouring of the map corresponds to the Standard-period Average Annual Rainfall (SAAR 61-90) for England and Wales.

Figure 2 .
Figure 2.Estimated parameters obtained from an L-moment estimation for each duration separately (dots) and from the proposed unified GEV distribution (lines)annual series.Colours indicate the different ToT stations.Distribution parameters in each panel are (clockwise from top left): the location, scale, lower bound and skewness.To make the figure readable, lower bound estimates below -40 are not shown.
and are shown in Prosdocimi et al. (2016).

Figure 3 .
Figure 3.Estimated frequency curves for the station at Dowdeswell for the annual (left panel), winter (central panel) and summer (right panel) series, superimposed on the Gringorten plotting positions for each duration, starting from 1 minute

Figure 4 .
Figure 4.Estimated parameters obtained from an L-moment estimation for each duration separately (dots) and from the proposed unified GEV distribution (lines)annual series.Colours indicate the different stations with series of at least 15minute accumulations.Parameters in each panel are (clockwise from top left): the location, scale, lower bound and skewness.To make the figure readable, lower bound estimates below -40 are not shown.

Figure 5 .
Figure 5.Estimated frequency curves for the station at Otterbourne for the annual (left panel), winter (central panel) and summer (right panel) series, superimposed on the Gringorten plotting positions, for durations from 15 to 120 minutes.

Figure 6 .
Figure 6.Relative difference between the 2-year rainfall depths estimated by different methods and the 50 th percentile of the recorded series.Larger dots correspond to stations with longer records.

Figure 7 .
Figure 7. Relative difference between the 5-year rainfall depths estimated by different methods and the 80 th percentile of the recorded series.Larger dots correspond to stations with longer records.

Figure 8 .
Figure 8. Relative difference between the 10-year rainfall depths estimated by different methods and the 90 th percentile of the recorded series.Larger dots correspond to stations with longer records.

Figure 9 :
Figure 9: Relative differences between the FSR (left panel) and FEH99 (right panel) estimates and the unified GEV estimates for all stations and all durations.Thick lines indicate the average differences for all durations across all stations.

Table 1 .
The correction factors applied to the maxima obtained from 15-minute series, for different seasons and event