Geomorphology-based unit hydrograph models for flood risk management : case study in Brazilian watersheds with contrasting physiographic characteristics

Heavy rainfall in conjunction with an increase in population and intensification of agricultural activities have resulted in countless problems related to flooding in watersheds. Among the techniques available for direct surface runoff (DSR) modeling and flood risk management are the Unit Hydrograph (UH) and Instantaneous Unit Hydrograph (IUH). This study focuses on the evaluation of predictive capability of two conceptual IUH models (Nash and Clark), considering their original (NIUH and CIUH) and geomorphological approaches (NIUHGEO and CIUHGEO), and their advantages over two traditional synthetics UH models – Triangular (TUH) and Dimensionless (DUH), to estimate DSR hydrographs taking as reference two Brazilian watersheds with contrasting geomorphological and climatic characteristics. The main results and conclusions were: i) there was an impact of the differences in physiographical characteristics between watersheds, especially those parameters associated with soil; the dominant rainfall patterns in each watershed had an influence on flood modeling; and ii) CIUH was the most satisfactory model for both watersheds, followed by NIUH, and both models had substantial superiority over synthetic models traditionally employed; iii) although geomorphological approaches for IUH had performances slightly better than TUH and DUH, they should not be considered as standard tools for flood modeling in these watersheds.


INTRODUCTION
Among the different observed natural hazards, Brunda and Shivakumar (2015) reported that flooding is the most frequent.According to Ma et al. (2014), several regions over the world have suffered from problems related to flooding.Its increasing trend is related to the population growth and occupation of higher risk areas, to the economic development and to the effects of climatic change on the hydrological cycle, especially on extreme events (Brunda and Shivakumar 2015, Perumal and Price 2013, Ryu et al. 2016).Floodings frequently result in a loss of human lives and of financial resources (Yadav et al. 2015); therefore, a rapid and accurate forecasting tool is fundamental for rivers prone to flooding.However, monitoring and early warning systems of floodings are not commonly used in developing countries, especially for small and medium-sized watersheds.
Stream flow data sets are fundamental to determine the hydrological response of a watershed to a heavy rainfall event whenever flood risk management or hydrological and hydraulic designs are necessary.The lack of stream flow series has culminated in the development of models, primarily intended for peak stream flow estimation and hydrographs based on a rainfall event (Beskow et al. 2015, Caldeira et al. 2015).
It has been over a century since models applied to hydrology were built to individually represent the hydrological processes (Mello et al. 2016).The first rainfall-runoff model widely used was created in 1851 by Thomas J. Mulvaney, an Irish engineer, which is today called the Rational Method.One of the first studies regarding flood hydrograph modelling began in 1921 and it can be attributed to the Australian engineer Cecil Napier Ross, who developed the concept of a time-area histogram, which seeks to represent the delay of the direct surface runoff (DSR) on each watershed zone and has been applied in various distributed hydrological models available in literature (Beven 2011).
In the case of watersheds with limited data bases, the choice of the model is considered a key step for estimation of DSR hydrographs.Unit Hydrograph (UH) theory, proposed by Le Roy K. Sherman (Sherman 1932), has been largely used.The pioneer study on synthetic UH was conducted by Snyder in 1938 based on data sets of gauged watersheds situated in the region of Appalachian (USA).Other synthetic UH models broadly accepted in literature are the Triangular and Dimensionless, developed in 1970's by Natural Resources Conservation Service (NRCS), formerly Soil Conservation Service (SCS) of United States Department of Agriculture (USDA) (McCuen andBondelid 1983, Folmar et al. 2007).These models are simple and their relations were derived from a large number of UHs that were obtained for American watersheds with different drainage areas and locations (Nourani et al. 2009).
The UH concept was enhanced after the proposition of the Instantaneous Unit Hydrograph (IUH), which corresponds to a DSR hydrograph resulting of a unit effective rainfall (R e ) with an infinitesimal duration (dt  0) applied uniformly over a watershed (Ahmad et al. 2009, Cleveland et al. 2006, Jeng and Coon 2003).Some conceptual models have been applied to derive IUH from R e hyetographs and DSR hydrographs (Jeng and Coon 2003), such as Clark (1945) and Nash (1957) models.Although the Nash's and Clark's IUHs, which are based on the theory of linear systems, have been widely used in some countries for hydrologic designs, their use is limited in Brazil, where traditional synthetic UH models (e.g., Triangular and Dimensionless) are mainly used instead.
UH and IUH models can be employed for gauged and ungauged watersheds (Sule and Alabi 2013).For gauged watersheds, they are derived from observed rainfall and stream flow records those with a geomorphological approach.Because their relations are derived from observed data sets from regions with specific characteristics, their generalization should be analyzed carefully, taking into account the region peculiarities and the quality of results.
In this context, the objective of this study was to evaluate the predictive capability of two conceptual IUH models (Nash and Clark), considering their original and geomorphological approaches, and their advantages over two traditional synthetic UH models (Triangular and Dimensionless), taking as reference the data sets of two Brazilian watersheds (Jaguara creek watershed -JCW and Cadeia river watershed -CRW) with contrasting geomorphological and climatic characteristics.JCW is a representative subwatershed of the Upper Grande river watershed (Southeastern Brazil) with respect to both land-use and soils classes (Oxisols predominantly) and has been a hydrology laboratory in the last 12 years.It drains towards to Camargos Hydropower Plant and has many river stretches vulnerable to flooding (SEMAD 2013).CRW has its drainage area within the municipalities of Canguçu, Morro Redondo and Pelotas, in the southern Brazil, where there have been reported many cases of natural disasters caused by flooding (CEPED 2011).

PHYSIOGRAPHICAL CHARACTERISTICS OF THE STUDY WATERSHEDS
Two Brazilian watersheds (Figure 1) were analyzed in this study including Cadeia river watershed (CRW), situated in Rio Grande do Sul State, and Jaguara creek watershed (JCW), located in Minas Gerais State.Cadeia is one of the main tributaries of the Pelotas River, which has an approximate drainage area of 940 km².Pelotas River is an important affluent to the São Gonçalo canal, which supplies water for human consumption and is a on an event basis making use of either numerical methods or conceptual modeling (e.g., Clark and Nash models).However, in the case of ungauged watersheds, geomorphological characteristics of the watershed can be used to determine necessary parameters for some synthetic models (geomorphological approaches).
Considering that conceptual IUH models depend on observed rainfall and stream flow records and, given the difficulty to have them, geomorphological approaches have been proposed to derive IUHs in watersheds.Rodriguez-Iturbe and Valdez (1979) developed a methodology that considers quantitatively the effect of the watershed geomorphology on IUH.This model can be interpreted as a probability density function of the time spent by the DSR to reach watershed's outlet since that this function depends on watershed morphology.Therefore, geomorphology is added to the IUH by means of indices associated with the drainage network.Bhagwat et al. (2011) stated that geomorphological IUH is important for programs linked to the watershed management when the absence of hydrological information is experienced.Yet, Ghumman et al. (2011) pointed out that this IUH is a useful model for prediction of hydrographs in watersheds due to the recent advances in satellite imagery and tools for data processing.
On evaluating the geomorphological approach proposed by Rodriguez-Iturbe and Valdez (1979) for estimation of Nash's IUH parameters, some scientific studies (Adib et al. 2010, Ghumman et al. 2014, Khaleghi et al. 2011) have been conducted in order to verify its applicability.Nonetheless, the results are quite variable, especially when compared to those produced by other models.Likewise, geomorphological approaches for Clark's IUH have been evaluated and the results are also contrasting among the studies (Kumar et al. 2004, Adib et al. 2010, Sahoo et al. 2006).Constant progress has been made regarding the development and application of UH and IUH models, mainly SAMUEL BESKOW et al. navigable watercourse that joins Patos to Mirim Lagoon.JCW has a drainage area of 29.5 km² and is located in the Upper Grande river region, in the southern Minas Gerais (Southeastern Brazil) and drains directly into the reservoir of the Camargos hydropower plant, owned by Minas Gerais State Energy Company (CEMIG).
The information on relief for CRW was extracted from the cartographic base developed by Hasenack and Weber (2010) on a 1:50,000 scale.
For JCW, the available information was a contour line map with points of known altitude obtained from Brazilian Institute of Geography and Statistics (IBGE).The geomorphological characterization of CRW and JCW was carried out in order to derive all the parameters necessary to be used as inputs for the UH and IUH models employed in this study with the aid of ArcGIS (ESRI 2014).The digital elevation model and land-use map for each watershed are illustrated in Figure 2, whereas, geomorphological characteristics and other important attributes for each watershed are summarized in Table I.

HYDROLOGICAL AND METEOROLOGICAL DATA SETS
Water stage records and the corresponding stream flow values on a thirty-minute time step were used for the outlets of both CRW and JCW.Three monitoring stations in CRW and one in JCW were considered for rainfall information in the same time step.Thiessen Polygons methodology was adopted to delineate the influence area of each rainfall monitoring station in the CRW and to compute aerial average rainfall in the watershed.These events allowed the elaboration of hyetographs and hydrographs, considering a 30-minute time step, which were applied for flood modeling based on UH and IUH theory.Given the necessity of runoff separation for UH and IUH modeling, in this study a methodology based on two inflexion points, defining the start and the end of direct surface runoff, was employed (Agirre et al. 2005, Choi et al. 2011) and the corresponding effective rainfall (R e ) was then determined for each event (R e = DSR depth).The R e hyetographs were determined using the Curve Number method (SCS 1971) following the procedures recommended by Chow et al. (1988) and applied in many flood modeling sudies (Beskow et al. 2009, Choi et al. 2011, López et al. 2012).Initial abstraction (I a ) was considered as the cumulative rainfall until the start time of the DSR observed on the hydrograph and the corresponding Curve Number (CN) was adjusted for each event.

Unit hydrograph (UH) models
In this study, two synthetic Unit Hydrograph models were evaluated: Triangular Unit Hydrograph (TUH) and Dimensionless Unit Hydrograph (DUH).These models were developed by SCS (1971) and are based on the analysis of a number of UHs derived from US watersheds with different drainage areas.The rising time (t a ) for TUH and DUH is considered as the sum of half time of the unit precipitation (P u ) duration and basin lag time (t lag ).The latter parameter accounts for the time from half of the unit precipitation (P u ) to the peak stream flow observed on the hydrograph.In this study, P u was defined as 1 mm occurring uniformly on the watershed during a time D of 30 minutes, while t lag was estimated employing the following equation (SCS, 1971): where t lag is the basin lag time (h), L corresponds to the length of the main river (m), S is the soil potential maximum water storage (mm), which is based on the Curve Number value, and X is the mean watershed slope (%).The CN value adjusted for each event was used in Equation 1 to calculate S.
The peak stream flow (q p ) for TUH and DUH was obtained by the relationship: where q p is in m3 .s -1 , P u corresponds to the unit rainfall (mm), A is watershed area (km2 ), and t a is the rising time (h).
Once t a , q p and base time (t b ) have been determined, the TUH's ordinates are defined at each time of interest.The DUH's ordinates are calculated from the ratio between stream flow and peak stream flow (q/q p ) for a number of ratios of time over rising time (t/t a ) in accordance with the equation: where X corresponds to a Gamma function of the peak factor (FP), which is commonly considered as 484 (Folmar et al. 2007).
The estimated DSR hydrograph (Q), using TUH and DUH, was derived from the convolution equation (Equation 4) in its discrete form for each UH model (TUH and DUH), event and watershed, which was compared to the respective observed hydrograph.For Equation 4 application, it is necessary to define the R e hyetograph in relation to the P u considered (P = R e /P u ), representing a rainfall event occurring uniformly on the watershed of interest and the watershed's UH (q ordinates).
1 1 ( ) where M is the ordinate of the estimated DSR hydrograph and N corresponds to the number of R e intervals.In both watersheds, a 30-min interval was used for the stream flows.

Instantaneous Unit Hydrograph (IUH) models
In this study, Nash's IUH (NIUH), Clark's IUH (CIUH), Nash's Geomorphological IUH (NIUH GEO ) and Clark's Geomorphological IUH (CIUH GEO ) were evaluated for each event and watershed.In all the IUH models assessed, the IUH was converted into an UH, considering ∆t = 30 min, similarly to the procedure employed by López et al. (2012).Afterwards, the estimated DSR hydrograph (Q) was also derived using the convolution equation (Equation 4) in discrete form for each IUH model.

a) Nash's IUH (NIUH)
This model (Nash, 1957) considers that a uniform rainfall occurs throughout the watershed and the surface runoff generated is routed to the respective outlet.The reservoirs account for both attenuations of peak stream flow and translation.The behavior of the reservoir in terms of surface runoff drainage is represented by a parameter k, which reflects the watershed response to the occurrence of one or more effective rainfall events.
The NIUH is described by: where u(t) denotes the NIUH's ordinates, t is the time interval, n and k correspond to the NIUH's parameters and Γ the Gamma function.The parameter k (scale) is equal for all the linear reservoirs and represents a storage constant, having a time unit; the parameter n (shape) is the number of linear reservoirs attenuating the IUH's peak stream flow.
The estimation of n and k for each event was carried out in accordance with the procedures described by Nash (1957), which are based on observed records derived from hyetographs and hydrographs.

b) Clark's IUH (CIUH)
The CIUH considers two important processes related to the transformation of effective rainfall into surface runoff (Clark 1945) -attenuation and translation.The application of the CIUH requires the estimation of two parameters (Ahmad et al. 2009, Che et al. 2014): time of concentration (t c ) and storage coefficient (R).Its mathematical formulation is given by: where u is the CIUH's ordinate, i refers to time, R E corresponds to the effective rainfall uniformly distributed which depends on both the Time-Area Histogram (TAH) and a constant for conversion of units, and C 0 and C 1 are weighting coefficients.
In this study, TAH was obtained with the aid of a tool named as Watershed Isochrones, a computer tool within the ArcGIS environment, developed by the team of the Hydrology and Hydrological Modeling Laboratory, at the Federal University of Pelotas (Brazil).As TAHs depend on the determination of t c , which in turn was adjusted for each rainfall event and watershed, they were elaborated from 30 minutes to 25 hours and between 30 minutes and 12 hours were made for CRW and JCW, respectively, in order to cover all the events in each watershed.The parameter R, in turn, was SAMUEL BESKOW et al.In this study, the NIUH GEO was considered through the formula for n and k proposed by Bhagwat et al. (2011), relating the IUH's scale (k) and shape (n) parameters to the Horton's ratios, based on the results of Rodriguez-Iturbe and Valdez (1979), as follows: where R B and R L correspond to the Horton's laws, represented by bifurcation ratio (R B ) and length ratio (R L ), and R A is defined as the ratio between areas of the watershed.
All the terms of the right side of the Eq.7 can be calculated because they depend on the geomorphological characteristics of the watershed (Table I).The parameter n was iteratively determined in accordance with Eq. 7, whereas, k was computed through the aforementioned dimensionless indices following recommendations of Bhagwat et al. ( 2011): where V is the dynamic velocity (m.s -1 ).Its determination was in accordance with the criterion indicated by Zelazinski (1986) and applied in other scientific studies (e.g., Ghumman et al. 2011, 2014, Kumar et al. 2002).
The CIUH is based on two parameters (t c and R) and on the determination of the watershed's TAH.For the application of the CIUH GEO , t c was estimated considering the empirical equations of Dodge and Ven Te Chow for CRW and JCW, respectively, considering their limitations regarding watershed's size.The storage coefficient (R) was estimated following the recommendation of Mohave County (2009): where t c is the time of concentration (hours), L corresponds to the maximum flow length (km), and A refers to the watershed area (km²).
The TAH for each watershed, corresponding to the t c previously calculated, was obtained with the aid of the "Watershed Isochrones" tool.However, for the CIUH GEO , R and t c were considered invariable in relation to the event, thereby, corresponding to a unique value for each watershed.

PERFORMANCE OF THE MODELS
Performance analysis was carried out by means of two statistical measures: Nash and Sutcliffe efficiency coefficient (C NS ) (Nash and Sutcliffe 1970), and relative error in the estimation of peak stream flow (RE Qp ).Moriasi et al. (2007)

RESULTS AND DISCUSSION
In order to have different scenarios of hydrological response, 12 rainfall-runoff events were evaluated in the CRW (Table II) and 20 events in the JCW (Table III).Some rainfall-runoff events were plotted in Figure 3 and Figure 4 for CRW and JCW, respectively, to emphasize the differences in hydrological behavior between these watersheds.
It can be observed in Tables II and III, that there was a substantial variability in R e values (CRW, 0.4 to 46.6 mm; JCW, 0.3 to 21.6 mm).This wide range of values is desired for UH modeling since it makes the analysis of models' performance broader.It was found that CRW (Figure 3) had rainfall events with duration considerably longer than JCW (Figure 4), which can be attributed to the predominant rainfall climatology of each watershed with respect to the events analyzed (CRW -frontal rainfall and JCWconvective rainfall).
On analyzing CN values (Tables II and III), adjusted for the rainfall-runoff events, one can observe that they were clearly variable among the events, ranging from 11.4 to 87.8 and from 24.8 to 98.7 for CRW and JCW, respectively.Findings similar to these were reported by Agirre et al. (2005) and López et al. (2012).Overall, it was observed a difference between the watersheds with respect to the delay for the beginning of surface runoff, which is related the longer rainfall events observed in CRW, thus culminating in lower maximum rainfall intensities (Figure 3).Such characteristics result in a greater potential for infiltration at the beginning of the rainfall event due to both the lower rainfall intensities and need of a longer time interval to satisfy initial abstractions (I a ).
The assessment of the contrasting rainfall types (Figure 3 and Figure 4) should be analyzed in conjunction with the soil types in both watersheds.According to classification of Embrapa ( 2013), the classes correspond to Red-Yellow Argisol and Grayish Brown Argisol.The argisols are normally deep to very deep soils, varying from well drained to imperfectly drained, and have a sequence of A-Bt-C or A-E-Bt-C horizons in that the B horizon is and increased clay texture in the Bt (Streck et al. 2008).The Grayish Brown Argisols differ from Red-Yellow Argisols basically due to their location in landscape, occurring in lowlands, thereby favoring more accumulation of water and organic matter.However, the predominant soil in JCW is Oxisols, which are very deep, with a high porosity and high clay content but occurring as strong aggregates (Beskow et al. 2013).These differing soil characteristics help elucidate why JCW tends   to experience more infiltration when compared to the CRW.The models' performance can be evaluated in accordance with two statistical measures (C NS and RE Qp ) for each event considered as summarized in Tables IV and V.
The C NS results from Tables IV and V indicate that the hydrographs on average were closer to those observed for both CRW and JCW when deriving from CIUH, followed by NIUH.On the other hand, the use of the DUH and NIUH GEO resulted in the lowest performance for CRW and JCW, respectively.With respect to RE Qp the same behavior was verified.
Yet analyzing C NS values, it was found that all the hydrographs estimated by both CIUH and NIUH to CRW fell into the "Very good fit" category according to the classification proposed by Moriasi et al. (2007).In the case of the JCW, the CIUH and NIUH resulted in hydrographs classified as "Very good fit" for all the events, except for event 16 (NIUH).However, CIUH resulted in a performance markedly greater than NIUH for JCW (see C NS values in Tables IV and V).
Taking the ER Qp statistic as reference, CIUH generated peak stream flows more accurate for a greater number of events in CRW (Table IV) when compared to the other models.The DUH model produced the worst peak stream flow estimation for CRW, while TUH, CIUH GEO and NIUH GEO also had lesser performances with hydrographs estimated for various events framed as "Inadequate" (Van Liew et al. 2007).For JCW (Table V), peak stream flows were better estimated by CIUH and NIUH, while TUH, DUH, NIUH GEO and CIUH GEO produced estimations classified as "Inadequate" for most of the events, thereby indicating that hydrologists should not apply such models, especially for peak stream flow estimation.
One can observe in Table IV that TUH and DUH models had unsatisfactory performance (see C NS values) for CRW, as only three and four events, respectively, out of 12 events presented estimated hydrographs framed as "Good fit" or "Very good fit".Relative to JCW (Table V), the behavior was similar as TUH and DUH models resulted in estimations classified as "Good fit" or "Very good fit" for six and five events, respectively, out of 20 events.The difficulties found in this study with respect to the accuracy of the hydrographs estimated by TUH and DUH go along with results obtained in other studies, such as in Khaleghi et al. (2011) andZakizadeh andMalekinezhad (2015).
The low performance observed in the present study regarding TUH and DUH models (Tables IV  and V) can be partially attributed to the difference between observed and estimated t lag values.Considering the monitored events (hyetograph versus hydrograph), the observed t lag values for CRW (Table I) ranged from 1.9 to 7.3 hours (average value of 4.6 hours), while for JCW (Table II), values between 0.6 and 4.3 hours (average value of 2.9 hours) were found.Otherwise, estimated t lag values were noticeably variable when analyzing a rainfall-runoff event at a time (CRW -2.4 h to 27.6 h; and JCW -0.74 h to 7.5 h).
The parameter t lag varies directly in function of the river length.Also, t lag is directly proportional to S and, consequently, inversely proportional to CN, thus reflecting the watershed's soil, land use and antecedent soil moisture characteristics.In general, t lag was overestimated, especially for CRW; this might have led to underestimates of unit peak stream flows (q p ).This discrepancy in t lag values has been also found in other studies as that carried out by Khaleghi et al. (2011), who evaluated the TUH and DUH to the Kasilian watershed (Iran) based on 21 events and found values from 2.6 to 7.4 hours.It should be emphasized that even using CN values adjusted for each event, t lag values were overestimated, thus indicating that the Equation 1 proposed by SCS (1971) should be carefully applied.
The NIUH model had satisfactory performance for both watersheds (Tables IV and V), with C NS values ranging from 0.65 to 0.98 and between 0.20 and 0.95 for CRW and JCW, respectively.On evaluating NIUH, CIUH and CIUH GEO for seven events in a watershed located in the central region of India, Kumar et al. (2002) noted that NIUH produced the best results for four events.Other studies also presented good results when applying NIUH, as Agirre et al. ( 2005), Saravanan and Manjula (2015), Zakizadeh and Malekinezhad (2015).Even though the NIUH and CIUH are based on observed rainfall-runoff records, Adib et al. (2010) reported that these models may not be capable of estimating DSR hydrographs for specific rainfall events.The results of this study do not corroborate the statement of Adib et al. ( 2010) because all the events were classified as "Very good fit" (C NS values) acording to Moriasi et al. (2010), except when analyzing the potential of NIUH to estimate peak discharges for some events in JCW in which RE Qp values were greater than 25% (Table V) (Van Liew et al. 2007).
On analyzing the applicability of NIUH for the rainfall-runoff events, it was possible to observe that CRW had an average n of 2.6, whereas, for JCW, it was 2.9.Relative to k parameter, it was found an average value of 3.3 h for CRW and 1.7 h for JCW.The amplitude and variability of n and k values observed in the events evaluated in this study corroborate the findings of several other studies, such as Kumar et al. (2002), Adib et al.  2009) and Sahoo et al. (2006).
It should be pointed out the performance of the NIUH GEO for CRW, resulting in estimated hydrographs framed as "Good" or "Very good" (Moriasi et al. 2007) for 8 out of the 12 events.Although C NS statistic has indicated this level of performance, peak discharge did not follow this tendency according to RE Qp values presented in Table IV.On the contrary, NIUH GEO did not have satisfactory accuracy in estimating hydrographs for JCW (Table V).Khaleghi et al. (2011) tested several UH and IUH models for flood estimation in a 68 km² watershed situated in Iran and verified that the geomorphoclimatic IUH, and the geomorphologic IUH initially developed by Rodriguez-Iturbe and Valdez (1979), had performances superior to TUH, Snyder's UH, DUH and Rosso's IUH.The results found in this study are not in total agreement with Khaleghi et al. (2011), who stated that the geomorphologic IUH (NIUH GEO ) can be used as a standard tool for rainfall-runoff transformation in watersheds with scarce data bases.
The results obtained in the present study are in consonance with those reported by Bhaskar et al. (1997), who also found that the NIUH had a performance better than NIUH GEO .These researchers stated that such behavior is expected since the NIUH's parameters are derived from observed data and partially attributed the inferior performance of the NIUH GEO to the use of a geomorphological equation for estimation of the NIUH GEO 's parameters.Bhaskar et al. (1997) highlighted that the equation used is not universal and its transposition for the watershed under analysis is questionable.
Relative to the NIUH GEO 's parameters, n was considered constant (CRW = 3.1; JCW = 2.9) for all the analyzed events.On the other hand, k ranged from 1.25 to 4.93 h (average value = 3.05 h) for CRW, while for JCW, it varied between 0.23 and 0.50 h (average value = 0.30 h).This parameter had a similar behavior to that observed for NIUH.It is worthwhile to highlight the difference in average value of k between NIUH GEO and NIUH (CRW, k = 3.25 h; JCW, k = 1.72 h), especially for JCW.Such variation in k values for event-by-event analysis has also been reported in literature, such as: i) Adib et al. ( 2010), using a similar methology, analysed 13 events monitored in a 67.5 km² watershed in Iran and obtained n equal to 2.83 and k varying from 1.28 and 3.40 h; and ii) Sahoo et al. (2006) verified the applicability of NIUH GEO for 15 rainfall-runoff events monitored for an Indian watershed (about 2,800 km 2 ) and derived geomorphological variables for modelling with the NIUH GEO from topographic maps at the 1:50,000 and 1:250,000 scales.The value of n was 3.74 for the former scale and 3.76 for the latter, whereas, k varied from 1.46 to 2.72 h for the most detailed scale and between 1.57 and 2.93 h for the other.
As aforementioned, CIUH resulted in the best accuracies when estimating hydrographs in that all the simulations were categorized as "Very good fit" (Moriasi et al. 2007) for both watersheds under analysis (Tables IV and V).In addition, its performance was substantially superior than NIUH, which was the other adjusted model assessed in this study.Sarangi et al. (2007) verified the performance of the exponentially distributed IUH (EDIUH), CIUH, and spatially distributed in Canada, and concluded that the EDIUH had the best performance.However, none of the models satisfactorily estimated the DSR hydrograph originated from events with long duration and low intensity, thereby diverging from the findings of this study.The CIUH was capable of appropriately capturing the behavior of hydrographs even for long duration and low intensity rainfall events, as commonly observed in CRW (Table I and Figure 3) due to its climatology.
The results found in this study regarding CIUH agree with those reported by Ghumman et al. (2014), who verified that CIUH had a performance better than NIUH and NIUH GEO for DSR hydrograph modeling in a watershed in Pakistan.Besides, these authors highlighted the importance of the timearea histogram used by CIUH since it accounts for temporal variability in direct surface runoff in the watershed of interest, while this component is not addressed in the NIUH.Contrarily, Ghumman et al. (2014) found that the NIUH provided the best accuracy for estimation of peak stream flows due to its consideration of a larger number of attenuations when compared to the CIUH.
The CIUH's parameters were quite variable among events analyzed for CRW and JCW; time of concentration (t c ) varied from 1.5 to 7.5 h for CRW and from 1 to 8 h for JCW, while R had values ranging between 2.37 and 6.82 h and between 1.00 and 4.38 h for CRW and JCW, respectively.Such amplitude in R and t c values have been observed in other studies concerning the CIUH (Adib et al. 2010, Kumar et al. 2002, Sahoo et al. 2006, Ahmad et al. 2009).Even though CIUH has demonstrated great potential to estimate hydrographs in CRW and JCW, it has not been widely employed due to the difficulties and limitations when: i) deriving R for both gauged and ungauged watersheds; and ii) determining the TAH corresponding to the watershed of interest since this task is quite laborious and the existing techniques are not sufficiently applied.Nevertheless, the development and application of the computer tool named as "Watershed Isochrones", which employs several hydrology operations combined with geoprocessing operations, made such task more straightforward, simple and rapid, given the number of TAHs generated for t c adjusted for each event and watershed.
Overall, CIUH GEO was not as accuracy as its adjusted version (CIUH) for JCW (Table V).It should be emphasized its performance for CRW (Table IV) in which 7 out of 12 events were satisfactorily estimated; however, peak discharges presented satisfactory estimation for only 5 events according to Van Liew et al. (2007) (see RE Qp in Table IV).Because the geomorphologic version of the CIUH used by Kumar et al. (2002) made use of characteristics extracted from the watershed of interest, while CIUH and NIUH were fitted based on the monitored rainfall-runoff events, these authors considered that such geomorphologic model resulted in fair estimated DSR hydrographs.
Adib et al. ( 2010) assessed the geomorphoclimatic CIUH, NIUH, CIUH and NIUH GEO to a 66.7 km 2 watershed in Iran and pointed out that the former and latter models were capable of adequately estimating DSR hydrographs.These researchers also reported that as the geomorphoclimatic CIUH and NIUH GEO do not require rainfall-runoff data sets for calibration, they have been widely employed in hydrology.Comparing the performance of the CIUH, NIUH and two geomorphologic versions of them for hydrograph estimation in an Indian watershed, Sahoo et al. (2006) found that the four models presented fairly accuracy.Nevertheless, according to Sahoo et al. (2006), the fact of not requiring both historical records and calibration makes the geomorphologic versions of CIUH and NIUH advantageous over the conventional CIUH and NIUH.SAMUEL BESKOW et al.
The time of concentration (t c ) along with the maximum DSR flow length (23.68 km for CRW and 10.59 km for JCW) and drainage area, made it possible to estimate R (Equation 9), resulting in 2.12 h and 1.28 h for CRW and JCW, respectively.The values of t c and R used for CIUH GEO were quite different from those adjusted by event for CIUH, especially for JCW, thus exerting influence on CIUH GEO performance (Table IV).As R and t c were not adjusted by event when applying the CIUH GEO , the procedures used to estimate their values may have exerted considerable impact on CIUH GEO modelling.Analogously, evaluating the CIUH GEO with 13 rainfall-runoff events, Adib et al. (2010) verified that the values of t c and R for each event were markedly different from those obtained when modelling CIUH adjusted in the environment of the HEC-HMS through monitored data sets.Sahoo et al. (2006) analyzed the influence of topographic maps' scale on determination of the geomorphological variables necessary for CIUH GEO , considering 15 rainfall-runoff events in an Indian watershed, and found differences in t c and R values.

CONCLUSIONS
Based on the results presented and discussed in this study, it could be concluded that: Not only physiographical characteristics (especially those linked to soil), but also the rainfall pattern, exerted influence on the flood modelling; The CIUH was the model that generated the most accurate surface runoff hydrographs for both watersheds, followed by the NIUH; these models had superiority in relation to the traditionally employed UHs (TUH and DUH) and to NIUH GEO and CIUH GEO ; The geomorphological versions of the CIUH and NIUH presented performances slightly superior when compared to the traditional and simpler models (TUH and DUH), however, they are not recommended to be applied as standard tool for flood modeling in the analyzed watersheds; Finally, this study indicated that it is recommended to evaluate: i) other approaches existing in literature to derive geomorphological parameters for NIUH GEO and CIUH GEO ; ii) the influence of relief information source on determination of parameters of the geomorphological IUH models; iii) the application of these models using other ways of spatial discretization of the watersheds for flood modeling, also considering the rainfall spatial distribution.

Figure 1 -
Figure 1 -Geographical location of CRW and JCW and setup of their hydrological monitoring network.

Figure 2
Figure 2 -a) Hydrologically Consistent Digital Elevation Model for CRW and JCW; and b) Land-use maps for CRW and JCW.
estimated seeking to minimize a given objective function through comparison between observed and estimated stream flows.All the procedures used in this study are in accordance with those applied byChe et al. (2014).c) Nash's Geomorphological IUH (NIUH GEO )

Figure 3 -
Figure 3 -Illustration of some rainfall-runoff events analyzed in CRW.

Figure 4 -
Figure4-Illustration of some rainfall-runoff events analyzed in JCW.