Development of a General Protocol to Enhance the Hydrological Analysis Techniques for Urban Catchments in Ireland

The hydrological analysis of urban catchments is a critical consideration for all major civil engineering projects. The purpose of this article is to develop a general protocol to enhance the established techniques for estimating runoff from Irish urban catchments. The chosen case study is the Rye Water catchment, which is covering two significant urban centres, Maynooth and Leixlip. The findings were compared with the flow data available on the Office of Public Works (OPW) website, to verify the reliability of the proposed approach. The study showed that the methodology described in this paper is applicable to Irish catchments and the protocol, with minor improvements can be adopted in Ireland to support designers in their approach to extreme rainfall events.


Introduction
Floods became a serious topic worldwide in the last ten years.The European Union (EU) issued a number of directives and regulations to face the problem, however, the diversity of climates and coping capacities of the different countries means the adoption of a consistent transnational approach is challenging.European countries accepted the regulations that were proposed by the EU as general guidelines adapting them depending on the population coping capacities, the climate of the country, the differences in terms of legislation, and many other factors.
One of the most relevant documents issued by the European Commission is the Flood Directive (FD) 2007/60/EC [1], which requires the member states: • to undertake Preliminary Flood Risk Assessments (PFRA); • to produce Flood Hazard Maps and Flood Risk Maps; and, • to develop Flood Risk Management Plans.
In order to produce the above-mentioned deliverables a high level of detail is needed, especially for the hydrological analyses.The FD 2007/60/EC does not give any indication regarding the procedure to use to develop the hydrological analysis, leaving the member countries free to follow the most appropriate approach.
Ireland can be considered in the average as regards flood risk data according to Hans de Moel's article "Flood maps in Europe-methods, availability and use" [2].
Despite the good material produced, Irish practitioners complain of the absence of a general procedure to follow in developing their hydrological analyses, particularly in the case of urban catchments.O'Sullivan et al. [3] assessed the difficulties that were encountered in selecting a suitable approach to develop hydrological analysis in Ireland.The study, undertaken as part of the Office of Public Works (OPW) Flood Studies Update Programme, comprised both quantitative (posted questionnaires) and qualitative research (focus groups) targeting a range of organisations that are involved in the planning and design of urban drainage infrastructure.This study highlighted, amongst other things, the relatively limited data available, particularly regarding soil parameters.
The main goals of this articles are: • to develop a general protocol that will be able to help users in developing hydrological analyses of urban catchments; and, • to compare the data available online and in literature with data estimated through the proposed methodology for two identified catchments.

Current Scenario
Many publications have been issued in order to guide users on the management of the flooding problem before and after any engineering projects.One of the key organisations regarding the Flood Management Strategies is the OPW a leading agency that is aimed to minimize the impacts of flooding through sustainable planning.The most relevant documents that are issued by the OPW in November 2009 are, as follows:

•
The Planning System and Flood Risk Management-Guidelines for Planning Authorities [4] • The Planning System and Flood Risk Management Guidelines for Planning Authorities-Technical Appendices [5] Both of the documents provide useful information for designers and local authorities to scope all possible impacts on the environment of future developments and indication on how to implement possible mitigation measures to reduce the associated risks.
The OPW has been also active in developing some websites containing flood information, interactive maps, and helpful publications to help end users in developing their hydrological analyses, including the Catchment Flood Risk Assessment and Management (CFRAM) Programme, the Flood Studies Update (FSU) Web Portal, and the Flood Info website.

The CFRAM Programme
The CFRAM Programme is a national programme commenced in 2011 and is central to the medium to long-term strategy for the reduction and management of flood risk in Ireland.As part of this programme, a website was released containing all of the information regarding the Catchment Flood Risk Assessment and Management in Ireland.
The Programme delivers on core components of the National Flood Policy [6], adopted in 2004, and on the requirements of the FD 2007/60/EC [1].The Irish CFRAM programme is being carried out in parallel with similar programmes across the European Union.
OPW works in close partnership with all local authorities in delivering the objectives of the CFRAM Programme.Implementation of the requirements EU Floods Directive is being coordinated with the requirements of the EU Water Framework Directive [7] and the current River Basin Management Plans [8].
According to the objectives set for the CFRAM Programme, the (FD) 2007/60/EC will comprise three stages:

•
The Preliminary Flood Risk Assessment (PFRA): 2011 The FSU Web Portal is an online website providing design rainfall and flood estimation tools for Ireland.
The FSU Research Programme was commissioned, managed and funded by OPW and is a substantial update of the Flood Studies Report (NERC, 1975) [9], which provided methodologies for flood estimation in Ireland and the United Kingdom.The FSU Web Portal is available for the OPW, local authorities, third level institutions, and private sector practitioners, and it represents a key instrument for professionals working in the area of flood risk assessment and management in Ireland.

Flood Info
Flood Info is a recently released website, which was launched on 3 May 2018 and providing access to the published Flood Plans along with the Flood Maps, as developed by the OPW and information about flood risk management in Ireland.The project is part of a more general context, in which the Minister of State Kevin 'Boxer' Moran T.D., in collaboration with the OPW and Flood Relief, published 29 Flood Risk Management Plans to address flood risk in Ireland.The Government has committed almost €1bn to flood relief measures as part of the National Development Plan 2018-2027 [10].
The main goal of the project is the completion of the largest ever study of flood risk in Ireland and the key outputs of the CFRAM Programme are as follows: • 40,000 Flood Maps showing the flood risk for 300 communities that support planning decisions and emergency response; and, • 29 Flood Risk Management Plans to cover the whole country, with the proposed flood relief measures-informed by costs, benefits and environmental factors-to address the flood risk in each community and nationwide.
Part of the project, has been produced in the report "Implementing the National Flood Risk Policy" [11], published on the website to accompany the publication of the Flood Plans.
Those initiatives, which were promoted mainly by the OPW, seem to sign a mark-up boundary in going towards a new approach of the flood risk management.
On the other side, there is a hidden risk in producing an over-abundance of documents.O'Sullivan et al. [3], after making a ranking the most common used methods (Modified Rational Method, Rational Method and Wallingford Procedure the three most commonly used) provides a clear picture of whose are the greatest difficulties and shortcomings that were encountered by the users in dealing with hydrological analyses, including: The emerged trend is not positive considering the number of engineering developments that are influenced by those analyses.The most complained aspect seems to be the lack of a national structured approach, which may help users in choosing the correct path in their professional applications.
A good example of national guidance is represented by "Delivering benefits through evidence-Rainfall runoff management for developments" [12], a report that was published in UK in October 2013 and commissioned by the Environment Agency.The purpose of the report is to provide advice to regulators, developers, and local authorities on the management of stormwater drainage for developments and in particular to assist in sizing of storage elements for the control and treatment of stormwater runoff.
The production of a similar document could help Irish users in the future in having only one reference document to follow in dealing with different type of developments.That document would need, in any case, as a starting point, a preliminary research programme to assess the most suitable approach depending on the types of developments.

Study Area
The study area (Figure 1) lies in between County Meath and Kildare, situated in north Leinster, along the east coast of Ireland.The catchments chosen for the study lie along the Rye Water, which rises in County Meath and flows South-East for 19 miles, and it is one of the main tributary of the river Liffey.developments and in particular to assist in sizing of storage elements for the control and treatment of stormwater runoff.The production of a similar document could help Irish users in the future in having only one reference document to follow in dealing with different type of developments.That document would need, in any case, as a starting point, a preliminary research programme to assess the most suitable approach depending on the types of developments.

Study Area
The study area (Figure 1) lies in between County Meath and Kildare, situated in north Leinster, along the east coast of Ireland.The catchments chosen for the study lie along the Rye Water, which rises in County Meath and flows South-East for 19 miles, and it is one of the main tributary of the river Liffey.

Proposal Methodology
The goal of the article was the delineation of a general protocol to simplify the methods that are adopted for the hydrological analyses of Irish urban catchments.The methodology has been applied on two catchments, as mentioned in Section 2.2, Catchment 1, and Catchment 2.
The methodology used to pursue the objective can be summarised in the phases as follows: 1. rainfall data analysis; 2. catchment definition in Geographical Information Systems (GIS); 3. evaluation of the time of concentration; and, 4. flow estimation These four simple stages can be expanded and detailed, but they give to the users a first idea on how to approach the study of the catchment of interest.The catchments considered in the research belong both to the Rye Water.The first one, of larger dimensions (200 square Kilometres approximately), named Catchment 1, and a sub-catchment of the first, named Catchment 2 and derived from an upstream closure section of the Rye Water.

Proposal Methodology
The goal of the article was the delineation of a general protocol to simplify the methods that are adopted for the hydrological analyses of Irish urban catchments.The methodology has been applied on two catchments, as mentioned in Section 2.2, Catchment 1, and Catchment 2.
The methodology used to pursue the objective can be summarised in the phases as follows: 1. rainfall data analysis; 2.
evaluation of the time of concentration; and, 4.
flow estimation These four simple stages can be expanded and detailed, but they give to the users a first idea on how to approach the study of the catchment of interest.

Rainfall Data Analysis
An indirect method was chosen to evaluate the flow for different return periods using rainfall data as input.
The first step was the collection of the rainfall data for the proposed catchment.The purpose of the study was, from the analysis of the historical data of extreme rainfall events, deriving the rainfall depth at 1, 3, 6, 12, and 24 h immediately subsequent to the event.The use of this extreme rainfall event data is quite common in Italy, where usually the regional organisation responsible for the collection and the publication of rainfall data also provides the distribution of the rainfall depth at 1, 3, 6, 12, and 24 h measured at the start of the event.Evidence of the approach can be found in literature in Crisci et al. [13] and Pinto et al. [14].
The challenge, when considering the unavailability of this type of data in Ireland, was the estimation of them starting by historical data available.The information have been supplied by the web portal of Met Eirann, the Irish Metereological Service web portal responsible for monitoring, analysing, and predicting Ireland's weather and climate.
The gauging stations considered relevant for the analysed catchments were initially four: A gauging station very close to the analysed catchment has been excluded from the analysis because of its lack of hourly rainfall distribution data.The hourly rainfall data in fact were essential to define Peaks Over Threshold (POT) rainfall events for durations of 1, 3, 6, 12, and 24 h.
The first stage was the identification of the daily maximum rainfall depth per year, through the daily distribution, and once identified the day of heaviest rainfall, it was possible to reconstruct the rainfall event through the hourly distribution of the rainfall depth.After cumulating at 1, 3, 6, 12 and 24 h the rainfall depths, it was possible producing the long series of the rainfall, the values for a number of years variable from 13 to 20 years.
The rainfall depths for 1, 3, 6, 12, and 24 h were the input for the application of the Generalised Extreme Values (GEV) probability distribution.The intent was the estimation of the Probable Maximum Precipitation and the calculation of the formula describing the temporal distribution of the rainfall depth [15] (p.137), as follows: where a represents the rainfall depth of duration one hour and fixed return period and n is a rainfall factor.This procedure has been applied to all of the four gauging stations above mentioned.

Catchment Definition in GIS
The second important stage was the definition of the catchment and the representation of the stream network in GIS.
The hydrological modelling has been performed in GRASS (Geographic Resources Analysis Support System) GIS, which is commonly referred to as GRASS, a free and open source GIS software suite used for geospatial data management and analysis, image processing, graphics and maps production, spatial modelling, and visualization.The software is part of Quantum GIS Package (QGIS), whicb is a professional GIS application that is built as Free and Open Source Software (FOSS).QGIS has been used in that research for a few conversions of raster and vector files.
The input data elaborated by GRASS was a Digital Elevation Model (DEM) of 30 m resolution of Ireland, imported in Tiff format and downloaded from the European Environment Agency database.
The elaboration of the DEM was functional to the representation of the final stream network of the area and the identification of the catchment of interest.

Evaluation of the Time of Concentration
Obtained from the shape of the catchment in GRASS, the vector file has been imported in QGIS and the Polygons of Thiessen have been traced to assess how much each gauging station was relevant in terms of rainfall characteristics for the proposed catchments.
In this method, the catchment area is divided into subareas using rain gages as hubs of polygons.The subareas are used to determine ratios that are multiplied by the subarea rainfall and summed to get the catchment average depth.The ratios are the percentages of area in the basin, represented by each rain gage [16].
The result, as shown in Figure 2, has been the exclusion of Dublin Airport and Phoenix Park from the following calculations.The purpose of the polygons of Thiessen, was, in fact, the evaluation of the features of a probable rainfall event within the catchment boundaries; considering that, the only two polygons included within were the ones of Dunsany and Casement their rainfall regimes can be considered as similar to the regime of the catchment.The proportion phase is to quantify how "similar" are these rainfalls.The final product was a temporal distribution of an extreme rainfall event for different return periods, respectful of the effect of each gauging station on the catchment of interest.

Evaluation of the Time of Concentration
Obtained from the shape of the catchment in GRASS, the vector file has been imported in QGIS and the Polygons of Thiessen have been traced to assess how much each gauging station was relevant in terms of rainfall characteristics for the proposed catchments.
In this method, the catchment area is divided into subareas using rain gages as hubs of polygons.The subareas are used to determine ratios that are multiplied by the subarea rainfall and summed to get the catchment average depth.The ratios are the percentages of area in the basin, represented by each rain gage [16].
The result, as shown in Figure 2, has been the exclusion of Dublin Airport and Phoenix Park from the following calculations.The purpose of the polygons of Thiessen, was, in fact, the evaluation of the features of a probable rainfall event within the catchment boundaries; considering that, the only two polygons included within were the ones of Dunsany and Casement their rainfall regimes can be considered as similar to the regime of the catchment.The proportion phase is to quantify how "similar" are these rainfalls.The final product was a temporal distribution of an extreme rainfall event for different return periods, respectful of the effect of each gauging station on the catchment of interest.It is important specify that the weighting of rainfall data stage was necessary for Catchment 1 only, once Catchment 2 was located in the area of influence of Dunsany, used to progress Phases 3 and 4.
The Polygons of Thiessen procedure was functional to the estimation of the concentration time, the time for a drop of water to flow from the remotest point in the catchment to a closure section of the catchment [17] (p.499).
The formulas to estimate the concentration time are dependent on the catchment characteristics, and this explains why the catchment definition in GIS has been put as the antecedent stage.
The estimation of the concentration time is a critical task, and considering the empirical nature of formulas that are available in literature, the results obtained have been averaged.
The formulas used are: • Kinematic Wave Formulas-Morgali and Linsley (1965) [18] (p.250) It is important specify that the weighting of rainfall data stage was necessary for Catchment 1 only, once Catchment 2 was located in the area of influence of Dunsany, used to progress Phases 3 and 4.
The Polygons of Thiessen procedure was functional to the estimation of the concentration time, the time for a drop of water to flow from the remotest point in the catchment to a closure section of the catchment [17] (p.499).
The formulas to estimate the concentration time are dependent on the catchment characteristics, and this explains why the catchment definition in GIS has been put as the antecedent stage.
The estimation of the concentration time is a critical task, and considering the empirical nature of formulas that are available in literature, the results obtained have been averaged.
The formulas used are: • Kinematic Wave Formulas-Morgali and Linsley (1965) [18] (p.250) (2) where To be solved in t c substituting the intensity with the following equation: where: -L is the length of channel/ditch from headwater to outlet where: where: The formulas that are mentioned are different, however most of them consider the geometric characteristics of the basin as input.The values measured in GRASS were the area, the height of the catchment above the closure section, the length of overland flow, the length of channel/ditch from headwater to outlet, and the average overland slope.The Gauss-Strickler roughness has been derived by values suggested in literature for similar context and the intensity rainfall on the base of the calculations performed in the rainfall data phase.
The value that was obtained with Kirpich equation has been discarded and the remaining values averaged.It was expected that Kirpich formula could provide unreliable values when considering that the equation has been developed from SCS data for seven rural basins in Tennessee with well-defined channel and steep slopes (3% to 10%), while the catchments of interest were both quite flat.

Flow Estimation Methods
The methods and formulas used in the third phase of the proposed general protocol have been different for Catchment 1 and Catchment 2, therefore the following sections will be distinguished.
Catchment 1 The methods used to estimate the flow for the Catchment 1 were:

Soil Conservation Service Curve Number
The first method used to evaluate the flow has been the Soil Conservation Service Curve Number (SCS-CN), an approach that was widely used as a simple method for predicting direct runoff volume for a given rainfall event [19].The CN parameter values corresponding to various soil, land cover, and land management conditions can be selected from tables, but it is preferable to estimate the CN value from the measured rainfall runoff data if available.
The SCS-CN method is based on the following basic form calculating runoff from rainfall depth, where -P is the total rainfall -I a is the initial abstraction -Q is the direct runoff -S is the potential maximum retention Based on a second assumption, that the amount of initial abstraction is a fraction of the potential maximum retention Equation ( 7) becomes The potential retention S is expressed in terms of the dimensionless curve number (CN) through the relationship taking values from 0, when S → ∞, to 100, when S = 0.This definition was originally applied to the English metric system (with S in inches).In the SI units (with S in mm), : According to [20], the Equation ( 9) has been changed as: The Equation ( 13) applied to Equation (10) provides: The input data have been prepared in Excel and the analysis have been performed in HEC-HMS.

Flood Studies Update Programme
The second method used was the Flood Studies Update (FSU) through the OPW Hydro NET portal.The portal aims to estimate Q med in ungauged catchments, or at sites for which there is no data to use or transfer data from a nearby site, which was preferably upstream or downstream from the site of interest.
All of the information regarding the theoretical background of the FSU are included in the Final Report-Work Package 2.3 Flood Estimation in Ungauged Catchments [20], developed by the Department of Geography of the Irish Climate Analysis and Research Units (ICARUS) in June 2009.
The use of the FSU was necessary to make a comparison with the data obtained through the SCS-CN.The comparison was helpful, not only to get a first confirmation of the exactness of the delineation of the catchment, but also to calculate a number of flow values of different return periods.

Rational Method
The rational method [20] has been in use for over 150 years and it remains the most widely used method to estimate peak flows from urban and small rural ungauged catchments [21].
The formula that is used links the peak flow distribution to the catchment area, rainfall intensity (mm/h), and runoff coefficient, as follows: where: - The method is based on the assumptions that rainfall intensity and storm duration are uniform over the area of study.In addition, the other two hypotheses are considered: duration of the storm equal to the time of concentration of the catchment and runoff coefficient constant during a storm [22].The equation above mentioned must be divided by 360 to work with International System units.
Catchment 2 For Catchment 2, when considering an area of approximately 20 Km 2 , have been applied three different methods to estimate the flow:

The Institute of Hydrology Report 124
The Institute of Hydrology Report 124 (IH 124) Report was a research to examine the response of small catchments, less than 25 km 2 , to rainfall, and to derive an improved flood estimation equation [23].A total of 84 sites were used to validate the method.
The report worked out also an equation to estimate mean annual flood, QBAR, for small rural and urban catchments: It also worked out an equation to estimate mean annual flood, QBAR, for small rural and urban catchments.

Rainfall Data Analysis
In an Indirect Method the Rainfall Data Analysis stage represents a key phase in describing the rainfall features of a catchment.
Below have been reported the most relevant tables and charts (Figure 3 and Table 1), according to the theoretical background, as detailed in Section 2.3.1.At first sight, it is quite clear how much the two pluviometric regimes are different each other.Casement is in fact characterised by higher rainfall depths, especially for return periods that are above 100 years.On the other hand, Dunsany shows similar values to Casement for rainfall of medium duration (within the six hours) and low return periods (5 to 50 years).We can assume that for small events and not long return period design, approximating one gauging stations with the other could be acceptable.For the purpose of the article and considering that we wanted to be as detailed as we can, we proceeded for Catchment 1 with the Thiessen Polygons.
Below have been listed in Table 1 the value of a and n for the two gauging stations.The parameters are needed for the formula (1), as presented in Section 2.3.1.At first sight, it is quite clear how much the two pluviometric regimes are different each other.Casement is in fact characterised by higher rainfall depths, especially for return periods that are above 100 years.On the other hand, Dunsany shows similar values to Casement for rainfall of medium duration (within the six hours) and low return periods (5 to 50 years).We can assume that for small events and not long return period design, approximating one gauging stations with the other could be acceptable.For the purpose of the article and considering that we wanted to be as detailed as we can, we proceeded for Catchment 1 with the Thiessen Polygons.
Below have been listed in Table 1 the value of a and n for the two gauging stations.The parameters are needed for the formula (1), as presented in Section 2.3.1.The rainfall data results that were obtained by Casement and Dunsany have been used both for Catchment 1 and Dunsany only for Catchment 2, according to the Thiessen Method, as mentioned in Section 2.3.3.

Catchment Definition in GIS
The input datum to define the catchment in GRASS has been the 30 m resolution DEM, downloaded in Tiff format from the European Environment Agency portal, which is an open access website providing DEM to be used for any working use.
The images below enclosed have been produced through the cartographic composer of GRASS and represent some important deliverables in a catchment definition.The process followed to obtain Figure 4, Figure 5a

Catchment Definition in GIS
The input datum to define the catchment in GRASS has been the 30 m resolution DEM, downloaded in Tiff format from the European Environment Agency portal, which is an open access website providing DEM to be used for any working use.
The images below enclosed have been produced through the cartographic composer of GRASS and represent some important deliverables in a catchment definition.The process followed to obtain Figures 4, 5a

Catchment Definition in GIS
The input datum to define the catchment in GRASS has been the 30 m resolution DEM, downloaded in Tiff format from the European Environment Agency portal, which is an open access website providing DEM to be used for any working use.
The images below enclosed have been produced through the cartographic composer of GRASS and represent some important deliverables in a catchment definition.The process followed to obtain Figures 4, 5a   The catchment created can be then converted into a vector file or a dxf and imported in CAD or in Google Earth.The procedure was applied for both catchments, reiterating the last point for the creation of Catchment 2.
The software GRASS has been very helpful in the catchment definition stage.GRASS is an user friendly software and allows for the creation of a mapset folder where all of the rasters created can be automatically saved during the process.In addition, the software allows for the 3D View of all the raster, including the DEM, a very important tool for professionals working with terrain models.

Evaluation of the Time of Concentration
The four gauging stations that are considered in the Rainfall Data Analysis have been weighted according to the Thiessen Method [16], as explained in detail in Section 2.3.3.
Once evaluated the rainfall parameters for Catchment 1, and used Dunsany parameters for Catchment 2, according to the results that were obtained during the second stage, it was possible to then summarise the geometric parameters of the catchments in Table 2 and use them as input for the time of concentration phase.
The formulas used for this stage have been in total four, with one value discarded from the average.It is always good practice to have more than one methods to evaluate the time of concentration, when considering the empirical nature of the formulas and the different measures involved.The procedure was applied separately on Catchment 1 and Catchment 2.
Below listed the main results that were obtained for Catchment 1 and Catchment 2 and the emerged observations.
The great advantage of the first method is the versatility in the evaluation of the time of concentration, function of the return period.
Looking purely at numbers, it can be seen that the time of concentration is approximately halved going from five years return period to 1000 years, a sign that time of concentration and return period are not linked by a close inverse correlation (increasing the return periods 200 times the time of concentration does not become 200 times smaller).
The values obtained through the Giandotti and Pezzoli methods are very similar to the values obtained for five years return period with Morgali and Linsley, therefore, for low scale designs we can assume that those methods would be applicable.In terms of average values, as expected, the highest return period values were the most affected by the use of the additional formulas.
In general, the final values could be considered reliable, when considering the nature of the surface and the dimensions of the catchment, however, it can be appreciated a reduction of time of concentration in the 1000 years return period scenario approximately around the 16%.
The times of concentration obtained have been implemented with the parameters of Table 3 to derive the temporal distribution of the rainfall, which plotted against the return periods provide the curve presented in Figure 7.The catchment created can be then converted into a vector file or a dxf and imported in CAD or in Google Earth.The procedure was applied for both catchments, reiterating the last point for the creation of Catchment 2.
The software GRASS has been very helpful in the catchment definition stage.GRASS is an user friendly software and allows for the creation of a mapset folder where all of the rasters created can be automatically saved during the process.In addition, the software allows for the 3D View of all the raster, including the DEM, a very important tool for professionals working with terrain models.

Evaluation of the Time of Concentration
The four gauging stations that are considered in the Rainfall Data Analysis have been weighted according to the Thiessen Method [16], as explained in detail in Section 2.3.3.
Once evaluated the rainfall parameters for Catchment 1, and used Dunsany parameters for Catchment 2, according to the results that were obtained during the second stage, it was possible to then summarise the geometric parameters of the catchments in Table 2 and use them as input for the time of concentration phase.
The formulas used for this stage have been in total four, with one value discarded from the average.It is always good practice to have more than one methods to evaluate the time of concentration, when considering the empirical nature of the formulas and the different measures involved.The procedure was applied separately on Catchment 1 and Catchment 2.
Below listed the main results that were obtained for Catchment 1 and Catchment 2 and the emerged observations.
The great advantage of the first method is the versatility in the evaluation of the time of concentration, function of the return period.
Looking purely at numbers, it can be seen that the time of concentration is approximately halved going from five years return period to 1000 years, a sign that time of concentration and return period are not linked by a close inverse correlation (increasing the return periods 200 times the time of concentration does not become 200 times smaller).
The values obtained through the Giandotti and Pezzoli methods are very similar to the values obtained for five years return period with Morgali and Linsley, therefore, for low scale designs we can assume that those methods would be applicable.In terms of average values, as expected, the highest return period values were the most affected by the use of the additional formulas.
In general, the final values could be considered reliable, when considering the nature of the surface and the dimensions of the catchment, however, it can be appreciated a reduction of time of concentration in the 1000 years return period scenario approximately around the 16%.
The times of concentration obtained have been implemented with the parameters of Table 3 to derive the temporal distribution of the rainfall, which plotted against the return periods provide the curve presented in Figure 7.A similar procedure has been applied also to Catchment 2, however, as above-mentioned, the Thiessen Method has been bypassed as initial stage, once the catchment was entirely within Dunsany.The rainfall parameters adopted were therefore the values of Dunsany, as shown in Table 1.
Implementing the plan dimensions of Catchment 2 in the time of concentration formulas have obtained the following results, summarised in Tables 4 and 5.
In Catchment 2, we can appreciate an inverse trend, comparing Morgali and Linsley and Giandotti and Pezzoli methods, in fact, in that case, the highest return period values obtained through the first method are the most similar to the times that were evaluated with Giandotti and Pezzoli.That relevance lead us to the conclusion that the application of those methods could be well implementable for small catchments and high return periods.
As regards the time of concentration obtained with Kirpich, which is usually applicable to small catchments, we can appreciate an improvement, but not sufficient anyway to promote the value in the average phase, and therefore it has been discarded.A similar procedure has been applied also to Catchment 2, however, as above-mentioned, the Thiessen Method has been bypassed as initial stage, once the catchment was entirely within Dunsany.The rainfall parameters adopted were therefore the values of Dunsany, as shown in Table 1.
Implementing the plan dimensions of Catchment 2 in the time of concentration formulas have obtained the following results, summarised in Tables 4 and 5.
In Catchment 2, we can appreciate an inverse trend, comparing Morgali and Linsley and Giandotti and Pezzoli methods, in fact, in that case, the highest return period values obtained through the first method are the most similar to the times that were evaluated with Giandotti and Pezzoli.That relevance lead us to the conclusion that the application of those methods could be well implementable for small catchments and high return periods.
As regards the time of concentration obtained with Kirpich, which is usually applicable to small catchments, we can appreciate an improvement, but not sufficient anyway to promote the value in the average phase, and therefore it has been discarded.

Flow Estimation
As anticipated in Section 2.3.4,the methods that were applied to evaluate the flow values for Catchment 1 have been different from the methods adopted for Catchment 2, therefore the following sections have been divided.

Soil Conservation Service Method
The Soil Conservation Service method (SCS-CN) has been applied through Excel to prepare the effective rainfall depth, and the data have been imported in the Hydrologic Modeling System (HEC-HMS) to run the analysis.HEC-HMS is an open source software, as provided by US Army Corps of Engineers, to simulate complete hydrologic processes of watershed systems.
An example of the performed analysis have been reported below in Figures 8 and 9.The unit of measurement used have been years for the return period and m 3 /s for the flow values.

Flow Estimation
As anticipated in Section 2.3.4,the methods that were applied to evaluate the flow values for Catchment 1 have been different from the methods adopted for Catchment 2, therefore the following sections have been divided.

Soil Conservation Service Method
The Soil Conservation Service method (SCS-CN) has been applied through Excel to prepare the effective rainfall depth, and the data have been imported in the Hydrologic Modeling System (HEC-HMS) to run the analysis.HEC-HMS is an open source software, as provided by US Army Corps of Engineers, to simulate complete hydrologic processes of watershed systems.
An example of the performed analysis have been reported below in Figures 8 and 9.The unit of measurement used have been years for the return period and m 3 /s for the flow values.The great utility of the software is in particular the possibility to run the model when considering many return period scenarios within the same file, helpful for several analyses running contemporary.
Below have been enclosed a summary table (Table 6) showing all of the values obtained through the SCS-CN Method applied through the HEC-HMS software.Return Period (Rp) 3 Flow (Q) 3  5 27.3 10 36.6 50 66.6 100 71.3 300 89.9 1000 109.8 3 Unit of measurements adopted: The flow values obtained with 5 years return period are almost double within the first 50 years return period to reach a flatter shape from 50 to 300 with a low increase until 1000 years return periods.The great utility of the software is in particular the possibility to run the model when considering many return period scenarios within the same file, helpful for several analyses running contemporary.

Flood Studies Update Programme
Below have been enclosed a summary table (Table 6) showing all of the values obtained through the SCS-CN Method applied through the HEC-HMS software.The flow values obtained with 5 years return period are almost double within the first 50 years return period to reach a flatter shape from 50 to 300 with a low increase until 1000 years return periods.

Flood Studies Update Programme
The Flood Studies Update Programme offer to its users the possibility to perform a series of rainfall and flood applications, as listed:

•
Rainfall Depth Duration Frequency (DDF) The tool used to perform the analyses has been the "Flood frequencies", which allow for the users on the base of a sequence of options to perform an analysis producing flood values and other flood information.
The stages required by the portal to evaluate the flood values are: Within each section, the designer can pick the most appropriate site and deliver a flood frequency analysis based on a series of input information he should know or decide on the base of the project.
Below included the main deliverables of the flow data estimation through the FSU Programme for Catchment 1 (Figures 10 and 11), as reported in the "Review Analysis".
The Flood Studies Update Programme offer to its users the possibility to perform a series of rainfall and flood applications, as listed:

•
Rainfall Depth Duration Frequency (DDF) The tool used to perform the analyses has been the "Flood frequencies", which allow for the users on the base of a sequence of options to perform an analysis producing flood values and other flood information.
The stages required by the portal to evaluate the flood values are: Within each section, the designer can pick the most appropriate site and deliver a flood frequency analysis based on a series of input information he should know or decide on the base of the project.
Below included the main deliverables of the flow data estimation through the FSU Programme for Catchment 1 (Figures 10 and 11), as reported in the "Review Analysis".Unfortunately, the OPW portal does not have a dedicated tool for maps creation and downloading a report, which is probably helpful to be added in a future version of the portal.

Rational Method
The Rational Method is one of the most common used methods worldwide and in Ireland, in particular, according to O'Sullivan et al. [3].
Below included in Table 7 the values obtained through the Rational Method applied on Catchment 1.The flow values obtained are very similar to the ones evaluated through the SCS-CN method.The most important consideration that can be made on the base of what obtained is that it is not important how much sophisticated the calculation method is, if the input values are correct and the method can be considered to be applicable to the catchment of interest, then the result should always be comparable.The results confirmed this hypothesis.Unfortunately, the OPW portal does not have a dedicated tool for maps creation and downloading a report, which is probably helpful to be added in a future version of the portal.

Rational Method
The Rational Method is one of the most common used methods worldwide and in Ireland, in particular, according to O'Sullivan et al. [3].
Below included in Table 7 the values obtained through the Rational Method applied on Catchment 1.The flow values obtained are very similar to the ones evaluated through the SCS-CN method.The most important consideration that can be made on the base of what obtained is that it is not important how much sophisticated the calculation method is, if the input values are correct and the method can be considered to be applicable to the catchment of interest, then the result should always be comparable.The results confirmed this hypothesis.

The Institute of Hydrology Report 124
The first method used for Catchment 2 to estimate the flow values has been the IH124, an often used method in Ireland and United Kingdom (UK) for small catchments.Below included a summary table (Table 8) of the results.At first sight, we can immediately appreciate a general decrease of flow values, when comparing the results with the data obtained for Catchment 1, with reduced variations going from low to high return periods.It can be then made the assumption that for small catchments the dependence on the return period is not as close, as seen for Catchment 1, ten times bigger than Catchment 2.

Flood Studies Update Programme
The procedure used through the OPW Portal is the same used for Catchment 1, however the results that were obtained are slightly different.Below presented the deliverables that were obtained (Figures 12 and 13) and immediately after the comments related.

The Institute of Hydrology Report 124
The first method used for Catchment 2 to estimate the flow values has been the IH124, an often used method in Ireland and United Kingdom (UK) for small catchments.Below included a summary table (Table 8) of the results.Return Period (Rp) 5 Flow (Q) 5  5 At first sight, we can immediately appreciate a general decrease of flow values, when comparing the results with the data obtained for Catchment 1, with reduced variations going from low to high return periods.It can be then made the assumption that for small catchments the dependence on the return period is not as close, as seen for Catchment 1, ten times bigger than Catchment 2.

Flood Studies Update Programme
The procedure used through the OPW Portal is the same used for Catchment 1, however the results that were obtained are slightly different.Below presented the deliverables that were obtained (Figures 12 and 13) and immediately after the comments related.An element worth being mentioned is the slight difference in the shape of Catchment 2, when compared to the one obtained in GRASS.That is probably due to the uncertainty of the OPW model derived for the closure section.The portal approximate, in fact, the catchment using a number of points that are available for users, is therefore very difficult to get exactly the same results in terms of shape

Rational Method
The rational method is probably for Catchment 2 the method more complete and respectful of the background information provided in phase 1, 2, and 3 of the proposed protocol.The results have been listed in Table 9.

Discussion and Conclusions
The findings of the research opened the possibility for users to apply an alternative method to develop hydrological analyses on urban catchments.
The different methods used to estimate flow data for Catchment 1 and 2 provided encouraging results in terms of values.
For Catchment 1, the data obtained have been summarised in Table 10.An element worth being mentioned is the slight difference in the shape of Catchment 2, when compared to the one obtained in GRASS.That is probably due to the uncertainty of the OPW model derived for the closure section.The portal approximate, in fact, the catchment using a number of points that are available for users, is therefore very difficult to get exactly the same results in terms of shape

Rational Method
The rational method is probably for Catchment 2 the method more complete and respectful of the background information provided in phase 1, 2, and 3 of the proposed protocol.The results have been listed in Table 9.

Discussion and Conclusions
The findings of the research opened the possibility for users to apply an alternative method to develop hydrological analyses on urban catchments.
The different methods used to estimate flow data for Catchment 1 and 2 provided encouraging results in terms of values.
For Catchment 1, the data obtained have been summarised in Table 10.It is evident that the greatest discrepancy is between the SCS-CN and the FSU method, providing in any case differences in percentage under the 20% with the two exceptions of flows evaluated for 5 and 10 years return period.The different values could be explained by the size of catchment, usually in fact, when hydrological analyses are performed for these catchments, the purpose is the designing of an important engineering structure and the return periods of design are above the 30 years.The deriving percentage is therefore acceptable if considered to be that for smaller catchments lower return periods would be chosen.
Looking at Catchment 2, the flow values that were obtained have been listed in Table 11.The results for Catchment 2 are definitely more similar between the three different methods.In fact, differences in terms of percentage error definitely lower than the 10% can be appreciated.
The main goal of the article was the creation of a general protocol, and the identified main steps emerged in the analysis have been reported in Figure 14.It is evident that the greatest discrepancy is between the SCS-CN and the FSU method, providing in any case differences in percentage under the 20% with the two exceptions of flows evaluated for 5 and 10 years return period.The different values could be explained by the size of the catchment, usually in fact, when hydrological analyses are performed for these catchments, the purpose is the designing of an important engineering structure and the return periods of design are above the 30 years.The deriving percentage is therefore acceptable if considered to be that for smaller catchments lower return periods would be chosen.
Looking at Catchment 2, the flow values that were obtained have been listed in Table 11.The results for Catchment 2 are definitely more similar between the three different methods.In fact, differences in terms of percentage error definitely lower than the 10% can be appreciated.
The main goal of the article was the creation of a general protocol, and the identified main steps emerged in the analysis have been reported in Figure 14.The first element of novelty added to the research has been the different approach to the rainfall data.The attention that was dedicated to the maximum rainfall event within the subsequent 24 h is in fact very different approach from the use of the Standard Annual Average Rainfall (SAAR) [24], which is one of the most commonly used parameters in Ireland and in UK, which is useful, but at the same time, very generic.
Another element of modernisation proposed in this study has been shaping of the catchment in GRASS.The creation of the catchment is in fact an important action in any hydrological analysis because it helps the users in understanding the flow network dynamics, but also the real conformation of the ground.
Nevertheless, one of the most innovative element introduced in the research was the application of the SCS-CN method to an Irish catchment.The method has not been used often, probably, when considering the lack of maps showing the Curve Number values for Ireland, however the method proved to be effectively applicable also to Irish catchments.The lack of a CN map has been solved through the estimation of a reliable value through the Curve Number assessment table [16].
The positive results emerged in the research encourage users in trying the alternative method proposed.The main advantages of the procedure are indeed:

•
A more basic and coherent subdivision of the working phases • A more engaged participation of users in making decisions relevant for their analyses • A better understanding of the catchment and the rainfall regime characterising it • A more defined view of the hydrological parameters that shapes the data • An help for users in becoming critical reviewer of their work

•
A different approach to the instruments available in literature and online The advantages mentioned are an incentive for users and institutions to think about an alternative procedure that is more user friendly and easy to approach.
Key insights for a possible future research could be: • the assessment of the methods available in literature in order to create a user manual for designers; and, • the creation of a unified database of information containing high quality data and maps, in particular, There are, in fact, several applications and GIS map tools that are available online, however the disaggregation of sources and the lack of guidance give the users the perception of being abandoned in the development of their analysis.
In conclusion, the general protocol identified can be represented as a viable alternative procedure to develop hydrological analysis in urban catchments and to support users in dealing with data and resource available in literature and online.

Figure 1 .
Figure 1.Google Earth view of the study area.The catchments considered in the research belong both to the Rye Water.The first one, of larger dimensions (200 square Kilometres approximately), named Catchment 1, and a sub-catchment of the first, named Catchment 2 and derived from an upstream closure section of the Rye Water.

Figure 1 .
Figure 1.Google Earth view of the study area.

Figure 2 .
Figure 2. Thiessen Polygons for the Catchment 1 and 2 (image created in Quantum Geographical Information Systems (QGIS)).

Figure 2 .
Figure 2. Thiessen Polygons for the Catchment 1 and 2 (image created in Quantum Geographical Information Systems (QGIS)).

Figure 3 .
Figure 3. (a) Probable Maximum Precipitation Curve for the gauging station Casement and (b) Dunsany.

Figure 3 .
Figure 3. (a) Probable Maximum Precipitation Curve for the gauging station Casement and (b) Dunsany.
,b and Figure 6a,b can be summarised as follows: • importing the DEM into the working area; • sinking the DEM too fill any possible gaps; • creation of the flow direction raster; • creation of the flow accumulation raster (the raster calculating the water flowing from the upstream areas into each pixel); • extraction of the stream network (creating a buffering zone of to the network); and, • definition of the catchment on the base of the closure section chosen on the network.Geosciences 2018, 8, x FOR PEER REVIEW 11 of 22 ,b and 6a,b can be summarised as follows: • importing the DEM into the working area; • sinking the DEM too fill any possible gaps; • creation of the flow direction raster; • creation of the flow accumulation raster (the raster calculating the water flowing from the upstream areas into each pixel); • extraction of the stream network (creating a buffering zone of to the network); and, • definition of the catchment on the base of the closure section chosen on the network.

Figure 4 .
Figure 4. Two-dimensional (2D) Geographic Resources Analysis Support System (GRASS) View of Catchment 1 and 2 visualised on the Digital Elevation Model (DEM).
,b and 6a,b can be summarised as follows: • importing the DEM into the working area; • sinking the DEM too fill any possible gaps; • creation of the flow direction raster; • creation of the flow accumulation raster (the raster calculating the water flowing from the upstream areas into each pixel); • extraction of the stream network (creating a buffering zone of to the network); and, • definition of the catchment on the base of the closure section chosen on the network.

Figure 5 .Figure 6 .
Figure 5. (a) 2D GRASS View of Catchment 1 visualised on the DEM (b) and showing flow network.

Figure 6 .
Figure 6.(a) 2D GRASS View of Catchment 2 visualised on the DEM (b) and showing flow network.

Figure 8 .
Figure 8. Hydrologic Modeling System (HEC-HMS) Summary Results table obtained for the five years return period scenario.

Figure 8 .
Figure 8. Hydrologic Modeling System (HEC-HMS) Summary Results table obtained for the five years return period scenario.

Figure 10 .
Figure 10.Flood Studies Update (FSU) Web Portal Screen Capture of Catchment 1.Figure 10.Flood Studies Update (FSU) Web Portal Screen Capture of Catchment 1.

Figure 10 .
Figure 10.Flood Studies Update (FSU) Web Portal Screen Capture of Catchment 1.Figure 10.Flood Studies Update (FSU) Web Portal Screen Capture of Catchment 1.

Figure 12 .
Figure 12.FSU Web Portal Screen Capture of Catchment 2.Figure 12. FSU Web Portal Screen Capture of Catchment 2.

Figure 12 .
Figure 12.FSU Web Portal Screen Capture of Catchment 2.Figure 12. FSU Web Portal Screen Capture of Catchment 2.

Table 1 .
Value characteristic of the Probable Maximum Rainfall distribution for Casement and Dunsany.The rainfall data results that were obtained by Casement and Dunsany have been used both for Catchment 1 and Dunsany only for Catchment 2, according to the Thiessen Method, as mentioned in Section 2.3.3.

Table 1 .
Value characteristic of the Probable Maximum Rainfall distribution for Casement and Dunsany.

Table 2 .
Plan Dimensions of Catchment 1 and Catchment 2.

Table 2 .
Plan Dimensions of Catchment 1 and Catchment 2.

Table 3 .
Value Characteristic of the catchment weighted between Dunsany and Casement.

Table 3 .
Value Characteristic of the catchment weighted between Dunsany and Casement.

Table 4 .
Estimation of the Time of Concentration Tc for Catchment 1.

Table 5 .
Estimation of the Time of Concentration T c for Catchment 2.

Table 4 .
Estimation of the Time of Concentration Tc for Catchment 1.

Table 5 .
Estimation of the Time of Concentration Tc for Catchment 2.

Table 6 .
Flow values estimated through the Soil Conservation Service Curve Number (SCS-CN) for Catchment 1.

Table 6 .
Flow values estimated through the Soil Conservation Service Curve Number (SCS-CN) for Catchment 1.

Table 7 .
Flow Values estimated through the Rational Method for Catchment 1.

Table 7 .
Flow Values estimated through the Rational Method for Catchment 1.

Table 8 .
Flow values estimated through the IH124 Method for Catchment 2.

Table 8 .
Flow values estimated through the IH124 Method for Catchment 2.

Table 9 .
Flow values estimated through the Rational Method for Catchment 2.

Table 9 .
Flow values estimated through the Rational Method for Catchment 2.

Table 10 .
Flow data obtained for Catchment 1 through the proposed methods.

Table 11 .
Flow data obtained for Catchment 2 through the proposed methods.

Table 10 .
Flow data obtained for Catchment 1 through the proposed methods.Unit of measurements adopted: Rp [years] Q [m 3 /s]. 7

Table 11 .
Flow data obtained for Catchment 2 through the proposed methods.