Flood Estimation and Control in a Micro-Watershed Using GIS-Based Integrated Approach

: Flood analyses when using a GIS-based integrated approach have been successfully applied around the world in large-sized watersheds. This study employed hydrological-hydraulic modeling to analyze ﬂash ﬂoods by integrating HEC-HMS, HEC-RAS, and ArcGIS software for ﬂood evaluation and control in a micro-watershed in the Samaru River, Nigeria. The watershed boundaries, its characteristics (soil and land use), the topographical survey, and the intensity duration frequency curve (IDF) of the study area were produced using data-driven techniques. The HEC-HMS model was used to derive the peak discharges for 2-, 5, 10, 25, 50-, 100-, and 200-year return periods with the frequency storm method. Afterward, the water surface proﬁles for the respective return periods were estimated using the HEC-RAS hydrodynamic model. The simulated design ﬂood for the 2-, 5, 10, 25, 50-, 100-, and 200-year return periods at the reference location (the NUGA gate culvert) were 3.5, 6.8, 9.1, 12.1, 14.3, 16.6, and 19.0 m 3 /s, respectively, while those at the watershed outlet for the respective return periods were 7.5, 14.9, 20.3, 27.3, 32.6, 38.0, and 43.5 m 3 /s, respectively (with a water height of 0.9 m, 1.1 m, 1.3 m, 1.33 m, 1.38 m, 1.5 3m, and 1.8 m, respectively), at the NUGA gate culvert cross-section. The maximum water depths of about 0.9 m and 1.0 m were recorded in the right and left overbanks, which were similar to the simulated water depth for the 2-and 5-year return periods. Hence, for the smart control of ﬂoods passing through the river and major hydraulic structures, a minimum design height of 1.50 m is recommended. For the most economic trapezoidal channel section, a normal depth of 1.50 m, a bottom width of 1.73 m, a top width of 3.50 m, and a free board of 0.30 m is proposed to curb the overtopping of ﬂoods along the channel sub-sections. The ﬁndings of this study could help hydraulic engineers minimize ﬂooding in streams and rivers overbanks in a micro-watershed.


Introduction
As a result of climate change, flooding is an increasing natural disaster that is bedeviling every nation worldwide.It causes catastrophic damage to ecosystems, life losses, and economic losses [1][2][3][4][5][6][7][8].In Nigeria, flooding effects are similar to those in the rest of the world [6, [9][10][11][12][13][14].According to Cirella and Iyalomhe, [15] flooding displaces more people than any other disaster in Nigeria because about 20% of the population are at risk of flooding.Recently, floods have displaced about 2.3 million people, and this has resulted in the unfortunate loss of 363 lives.Additionally, approximately 16 million individuals were affected in various ways, thus leading to severe economic losses.Consequently, a state of emergency was recently declared in nine states of the nation by the Nigerian government [10,11,13,16].
Storm water flows through an open channel with high velocity as runoff; this flow may be unsteady, non-uniform, or backflow-dependent upon the channel geometry, obstructions in the channel, catchment characteristics, and the amount of runoff from the catchment.High floods tends to erode channel bottoms, which may lead to channel bank collapse, as well as to modifications to the channel geometry [17][18][19].This usually results in spilling out of a reasonable quantity of water from the channel [18,20].It is not possible for a natural or artificial waterway to remain unchanged during the urban development process because land use changes have a great influence on runoff generation [18,21,22].Residential and industrial development have been the main activities taking place in micro-watersheds.In most cases, the stormwater from an upper catchment is conveyed through an unlined natural open channel to the watershed outlets.Thus, some sections of these channels collapse, and several diversion channels are formed by the flowing water; as such, a reasonable quantity of water can be lost via these channels, and this can be coupled with evaporation and infiltration losses.It is anticipated that these losses could be so high as to result in inadequate discharges, which are unable to reach the outlet.Moreover, a flow conveyance structure along the channel is normally constructed to allow access, but it might be in danger of washing away or collapsing as water overflows its top when there is excess rainfall.Having said this, owing to the above concerns, it has become imperative to model the flow and optimize channel geometry in order to design preventive and control measures through which to limit water loss and the collapse of conveyance structures.This can only be achieved by carrying out detailed hydrologic and hydraulic modeling.
Floodplain mapping has, in the past, been carried out using conventional methods [23,24].These methods are tedious, time-consuming, and relatively expensive.In addition, they require considerable fieldwork and the maintenance of long-term river and stream records, which are not readily available for catchments (which are mostly ungauged) in a large developing country such as Nigeria [9,25,26].These challenges could be addressed using an integrated approach of hydrologic and hydrodynamic modeling along with Arc-GIS analysis, which is less tedious and faster, more effective, and reliable, and it makes the immediate assessment of natural resources possible and more accurate [9,21,27].Hydrological modeling aids in determining the catchment response to rainfall events [21].Floodplain modeling is mainly concerned with the hydrological behavior of watersheds for a particular rainfall event.The development of flood management strategies and flood damage reduction and mitigation for an ungauged catchment has remained a big challenge to engineers and policymakers [6,28].This is genuinely true because detailed hydrological modeling is highly necessary for the estimation of flood volume and its recurrence time, which are the key parameters on which floods solutions are based [29].Hydraulic modeling results help to provide a set of actions and preventive measures that can be undertaken to curtail flood occurrences and flood aftermaths [30][31][32][33].
In this study, the Hydrologic Engineering Center-Hydrological Model System (HEC-HMS) and the Hydrologic Engineering Center-River Analysis System (HEC-RAS) were coupled alongside Arc-GIS to perform hydrologic and hydrodynamic modeling in a microwatershed.These packages were produced by the Hydrologic Engineering Center (HEC) of the US Army Corps of Engineers.The latter was used for simulating rainfall runoff, and it has the capability of modeling a wide range of geographic expanses, including flood hydrology [18,21,34].The former was used for the simulation of design floods and water surface profiles, which relies on I-D energy equation results [35,36].Abdessamed and Abderrazak, [37] who coupled their work with HEC-RAS and HEC-HMS, performed rainfall runoff modeling for floodplain inundation mapping in a data-sparse environment.Laouacheria and Mansouri [38] employed the frequency storm approach in the HEC-HMS model to simulate runoff in a small urban catchment in the northeast region of Algeria.An integrated hydrologic and hydraulic modeling was performed for a mediumsized, ungagged urban watershed using HEC-HMS and HEC-RAS by Natarajan and Radhakrishnan [35].Bruno et al. [39] applied hydrological and hydraulic modeling via integrating HEC-HMS and HEC-RAS, and they analyzed flash flood events in a small urban stream.
In this study, we conducted a comprehensive analysis of a data-sparse micro-watershed in Nigeria by employing a coupled hydrological-hydraulic modeling methodology that uses HEC-HMS and HEC-RAS.The goal was to provide a framework for other researchers to assess and validate their own flash flood modeling methodologies within a similar microwatershed context.The study analyzed flood frequencies, peak discharges, and their effects on the flood frequencies.The specific objective of the study was to determine water surface profiles for different return periods and to propose smart flood control structures that will safely convey the excess floods to the watershed outlet.It is hoped that the information from this study will help inform the policy and decision makers within relevant authorities to be aware of flood-prone areas along streams, which will, in turn, help in planning for future flooding events.

Study Area and Data Collection
The study area is within Zaria, Kaduna State, Nigeria.Zaria is one of the provinces that makes up the Central High Plains of Northern Nigeria, and it is approximately 670 m above mean sea level (Figure 1).It is located about 950 km away from Lagos; it covers a total landmass of about 61 km 2 , thus making it a nodal point in terms of road and rail transport because of its locational factors [40,41].[40,41].This stream flows in the southeast direction of Samaru through Ganga Uku and Ahmadu Bello University, Zaria to the Kubanni reservoir.For the land use characteristics of the watershed (Figure 2a), it has six main classes, which are as follows: urban build up (57.5%), open space (4.2%), agricultural land (8.1%), vegetation (26.2%), wet vegetation (3%), and Water bodies (1%).The highest elevation in the watershed is 716 m, while the lowest elevation is 640 m.Also, the catchment is very sloped at the downstream part of the Samaru River, as shown in Figure 2b.The mean slope of the watershed is 2.8%, while the maximum and minimum slopes are 15% and 0%, respectively.Analysis of the soil hydrologic characteristics shows that the dominant hydrological soil group is C, which covers 3.76 km 2 (90.4%) of the watershed and is at the upstream part (Figure 2c); meanwhile, group D covers 0.39 km 2 (9.6%) at the far downstream near the outlet.The curve number (CN) grid of the watershed shows (Figure 2d) that the maximum curve number is 100 while the minimum is 70, and the mean curve number in the watershed is 86.2.The data used in this study and their sources are summarized in Table 1.The channel cross-section data was obtained through a topographical survey that was constructed using Total Station (Tables S1 and S2).Climate and geospatial data were utilized to drive the HEC-HMS and HEC-RAS models.The Landsat 8 Operational Land Imager (OLI) with spatial resolution of 30 m, wasobtained from the United States Geological Survey (USGS) Global Visualization Viewer GLOVIS.The land use classification was achieved using the supervised classification method, and it was evaluated using onsite data.The accuracy of the LULC map was assessed using ground truthing and by comparing the data from the field.A confusion matrix was constructed to estimate the producer, user, and kappa accuracy (with a kappa index of 96%).The Shuttle Radar Topography Mission (SRTM) DEM data with a spatial resolution of 30 m were downloaded from GLOVIS, and the source of the data for the soil map was the Food and Agricultural Organization (FOA).A time series of the climate data (rainfall data from the Nigerian College of Aviation Technology, Zaria station) was obtained from the Nigerian Meteorological Agency (NiMet), which covered the period from 1967 to 2017 (i.e., 50 years of data).

Intensity Duration Frequency Curve IDF
Intensity duration frequency curve (IDF) development is the first step for HEC-HMS hydrologic simulation when using the frequency storm method.The peak discharges (Qmax) were determined for different return periods.The maximum daily (24 h) precipitation data for 1968-2017 in Samaru was reduced to a shorter time scale of 10, 20, 30, 60, 120, 360, 720, and 1440 min using the formula recommended by the Indian Meteorological Department [42,43].The formula has been widely used and has gained acceptance worldwide [44][45][46][47][48].
The required precipitation depth of less than 24 h is given by Equation (1). , where P t = the required precipitation depth for a duration less than 24 h in mm, P 24 = the daily precipitation depth in mm, and t = the required time duration in hours.
To determine the intensity, duration, and frequency curves of the Samaru stream watershed, vital input data for the HEC-HMS simulation and Gumbel's statistical distribution for frequency analysis were utilized.This method is often used to predict extreme hydrological events such as floods [49,50].The cumulative distribution function (cdf) of the Gumbel extreme value distribution (maximum) is given by the following [50]: where F(X) = the probability distribution function of a random variable x, and y = the reduced variate.
The probability of exceedance is related to the return period (T).The probability of exceedance of the excluded events is thus where F 1 (X) = the probability of exceedance at a return period (T), and T = the return period.
Combining Equations ( 2) and ( 3) yields where T = the return period, and y = the reduced variate.
The reduced variate from Equation ( 4) may be transformed into the following form The Gumbel's reduced variate is represented by a linear relationship with Gumbel's reduced mean variable y n , standard deviation σ n , and frequency factor K. Thus, we have Hence, K is then The frequency factor K for a Gumbel distribution is also given by Equation ( 8) where T = the return periods, and K = the frequency factor.The Gumbel equation for the total maximum precipitation for a given return period is given as where P T = the magnitude of precipitation for a T return period, P = the mean maximum precipitation for n years, and σ x = the standard deviation of the sample size.The mean maximum precipitation for n years P is given by the expression below: where P = the mean maximum precipitation for n years, P i = the maximum precipitation of the ith year, n = the number of sample size/number of years of record.The standard deviation of the sample size σ x is given by Equation ( 11) where σ x = the standard deviation of the sample size, P = the mean maximum precipitation for n years, P i = the maximum precipitation of the ith year, and n = the sample size/number of years of record.The intensity of precipitation (mm/h) for any time duration t for all the respective return periods was calculated using Equation ( 12) where I = the intensity of rainfall in mm/h, P T = the precipitation depth in mm, and t = the time duration of rainfall in hours.
For the estimation of the maximum precipitation when using the Gumbel distribution in this study, 50 years of maximum annual precipitation data of the Samaru stream watershed from 1968 to 2017 were considered.

Model Setup and Evaluation
The methodological sequence of the model connections, evaluation, and development used in this study, from the onset to completion, is presented in Figure 3.

HEC-HMS Model
The HEC-HMS basin model contains information relevant to the physical attributes of a watershed, such as the basin area, river reach connectivity, or reservoir data.The model of the Samaru basin contains four of these elements, and they were developed in Arc-GIS using the HEC-GeoHMS extension.A total of 18 hydrologic elements were generated for the Samaru model, which consists of 9 sub-basins, 4 river reaches, 4 junctions, and 1 outlet (at the point where the Samaru stream drains into the Malmo River).Although, the meteorologic model was developed based on the climatic data of the study area, afterward, the control specifications were defined.In this study, hourly time intervals were used and simulated for only two days.This period was selected based on the temporal size of the observed streamflow.The HEC-HMS model for the Samaru watershed was calibrated and validated with the August and September 2014 daily discharge data, which were obtained from a gauge station located at the outlet of the watershed.The choice of the August and September data was made because we are interested in peak discharge, and in those periods of the year the watershed usually experiences flood events.An auto-optimization operation and the peak weighted root mean square error (PWRMSE) were selected as objective functions in the HEC-HMS model during the calibration of the model parameters.This approach was selected because of its simplicity and high performance [21,37,51].A sensitivity analysis of the model parameter was performed by varying the values of the respective model parameters by ±20% at 5% intervals.The parameters checked were the curve number, lag time, percentage imperviousness, and Muskingum K.The percentage change in the simulated peak flow and volume as a result of variations in the respective model parameters was plotted.The model's performance was evaluated using statistical hydrologic indices of the agreement between the observed and simulated runoff values.The indices used were the root mean square error (RMSE), the Nash-Sutcliffe efficiency (NSE), and percentage bias (PBIAS).
where Q obs = the observed storm runoff (mm), Q comp = the computed runoff (mm), Q abs = the mean observed storm runoff (mm), N = the total number of rainfall runoff events, and i = an integer varying from 1 to N.

HEC-RAS Model
The development of the HEC-RAS model was a stepwise process that requires three main components: geometric data, flow data, and plan data.The objective of hydraulic modeling was to convert flow values into water surface elevation along the stream.HEC RAS solves the energy equation by using the standard step-backwater method to calculate the water surface profiles of successive channel cross-sections [52].Equation ( 16) presents the energy equation.The geometric data involves the establishment of a river system schematic (connectivity of the river system).It was achieved by entering cross-section data, which involved defining all of the necessary junction information, adding hydraulic structures, and interpolating cross-sections.The cross-section data were entered after the river schematic process was completed.The cross-section data gave the stream a geometric boundary.The cross-sections were measured at representative locations from Ganga Uku to the ABU detention basin.Manning's roughness coefficient selection involves judgment, skills, and subjectivity.Barnes [53] presented detailed photographs of natural rivers and creeks, as well as their respective roughness values.Manning's "n" value was assigned based on field observations and via a comparison of the Samaru stream main channel and the Barnes standard photographs.After the comparison, the roughness value was found to be between 0.03 and 0.035, and it was selected for the Samaru stream main channel.The expansion and contraction coefficients were entered manually into the HEC RAS cross-section data editor.USACE [54] gives typical values of expansion and contraction coefficient-based flow transition types.Based on the knowledge of the physical characteristics of the stream, the transition was gradual, and the default values of 0.3 and 0.1 were assigned for expansion and contraction coefficients, respectively.Seven pieces of discharge data were entered for the 2-year, 5-year, 10-year, 25-year, 50-year, 100-year, and 200-year return periods, and these were obtained from HEC HMS model outputs.The critical depth was selected as the upstream boundary condition, whereas the channel slope was selected as the downstream boundary condition (after which a steady flow analysis was executed).
where y 1 and y 2 are the water depths in the two cross sections, Z 1 and Z 2 are the heights of the main channel above the datum, v 1 and v 2 are the average velocities, α 1 and α 2 are the velocity weighing coefficients, g is acceleration due to gravity, and h e is the head loss of energy level.

Hydraulic Evaluation
The floods for various return periods simulated by HEC-HMS software were utilized as flow data in HEC-RAS for surface water profiles and velocity computations; in addition, basic parameters were required for hydraulic evaluation.The most economic channel section was designed for the Samaru stream using the HEC-HMS and HC-RAS outputs.According to Chow [55] the most economic channel section is the one where the maximum discharge gives the least wetted perimeter.Manning's equation was used for the computation of the channel dimension, and a free board of 20% of the normal flow depth was selected.The channel stability was checked based on the procedure recommended by Muhammad et al. [56].The mean velocity and tractive stress were first calculated.Afterward, an assessment of the existing stability was achieved by comparing the estimated local and instantaneous shear velocities to the velocity values.If the existing conditions were found to be stable and were in agreement with other project objectives, the channel was considered stable and no further analysis was required.

Intensity Duration Frequency Curve (IDF)
The IDF of the Samaru watershed was developed using daily rainfall data for a span of 50 years.The maximum daily rainfall data for each year was discretized into short time steps.The Gumbel distribution was used for flood frequency analysis.Table 2 presents the rainfall depth for different return periods at different time steps, while Figure 4 shows the IDF curve of the Samaru stream watershed.From both Table 2 and Figure 4, it was observed that the rainfall depth increased with an increase in the return period.The IDF and the return period were used as the input in the meteorological model of HEC-HMS for peak flow simulation at various return periods.

HEC-HMS Parameters Sensitivity Analysis
Model sensitivity analysis was conducted to determine which parameters of the model had a significant influence on the model output.The contribution of model parameters to an overall error in the model predictions was ranked with the aid of sensitivity analysis.Four ( 4) parameters of the model were subjected to sensitivity analysis.The curve number, impervious area, lag time, and Muskingum K were checked on a ±20% increment with a 5% interval.The curve number was determined to be the most sensitive parameter, as evident from Figure 5.An increase in the curve number by 20% made the peak discharge increase by 13.2%, and, when the curve number was decreased by 20%, the peak discharge was observed to decrease by 12.1%.The reason for it being the most sensitive parameter is that it depends on land use, hydrological soil groups, and antecedent moisture conditions, which are known to be the major influential factors in runoff generation.Lag time was also found to be sensitive as it tended to increase the peak discharge by 7.2% when increased by 20%, and it decreased the peak discharge by 6.8% when shortened by 20%. Figure 6 shows how these parameters deviated from the origin.

HEC-HMS Model Calibration, Validation, and Performance Evaluation
The model was found to be sensitive to curve number (CN), lag time, and Muskingum K.The HEC-HMS model calibration was achieved by auto-optimizing the sensitive parameter of the model.The parameters optimized were the curve number, lag time, and Muskingum K.The optimization operation was conducted on parameters that were sensitive to the model output.Table 3 gives the initial and optimized values of the sensitive parameters of the model.The discharge data from August 2014 were used for the calibration, which spanned from 1 August to 31 August 2014 and showed good agreement between the simulated outflows at all junctions, including the outlet.The observed outflows that produced hydrograph shapes and the time to peak were relatively the same, as shown in Figure 6.The root mean square error standard deviation (RMSE), Nash-Sutcliffe coefficient (NSE), and percentage bias (PBIAS) before and after optimization are presented in Table 4, while the simulated peak discharge and the total volume at the outlet were 5.3 m 3 /s and 583.1 mm before optimization and 4.1 m 3 /s and 401.7 mm after optimization.However, the performance of the model improved significantly when the RMSE, Nash-Sutcliffe coefficient, and percentage biases were within the acceptable threshold limit.The model was validated with the daily discharge data for September 2014, which spanned from September 1 to September 30.The validation was achieved by running the calibrated model with the September 2014 (Figure 6) events.The validation shows a very close match between the observed and simulated outflows at the watershed outlet, but it slightly underestimated the peak discharge.Based on model performance indices, the result was closely similar to the result obtained in the calibration session.The RMSE, NSE, and PBIAS of the validation result were 0.50, 0.78, and 20.8%, respectively, indicating very good performance (Table 4).These results are in good agreement with those obtained by Derdour et al. [51] (who worked on the Ain Sefra watershed and the Ksour Mountains (SW Algeria)) and Khaddor et al. [57] (who modeled the Boukhalef watershed in Morocco) when using the HEC-HMS model.Moreover, other researchers used the SCS-CN method for rainfall runoff simulations, and their model performance results matched with the results obtained in this research.

Simulation of Peak Runoff
The simulation results of the HEC-HMS model when using the frequency storm method for 2-year, 5-year, 10-year, 25-year, 50-year, 100-year, and 200-year periods at the upstream (Ganga Uku culvert), the NUGA gate culvert, towards the outlet (Dan Fodio culvert), and the watershed outlet in m 3 /s are presented in Table 5.It was observed from the results that the discharge increases with increased flow length.This is because of the contributions the river receives from adjacent sub-basins.For instance, the peak discharge at the NUGA gate culvert for a 50-year return period was 14.3 m 3 /s, which is less than the 27.3 m 3 /s at the Dan Fodio culvert, which is at its downstream.This is due to the fact that the stream receives runoff from Dan Fodio's impervious surfaces through the tributary upstream of the Dan Fodio culvert.The result also follows a normal pattern of hydrological analysis, as the peak discharge was observed to increase with an increase in the return period.The flood hydrograph for the 100-year Tr at the NUGA gate culvert and the watershed outlet are shown in Figure 7.As seen from the figure, the flood hydrographs have a shape similar to the ones obtained by Derdour et al. [51].

Hydraulic Evaluation of the River System of Watershed
The results for the hydraulic simulation with HE-RAS and the water surface profiles of some of the sections (where overtopping and spilling occur for flows of different return periods) were reviewed.Manning's coefficient was not calibrated due to the lack of highquality aerial imagery data.However, Manning's values were carefully assigned to the main channel and the overbanks based on rigorous field observation and via consultations from literature sources.The culvert height was measured to be 1.24 m during a topographical survey, but the simulation results for the 2-year, 5-year, 10-year, 25-year, 50-year, 100-year, and 200-year return periods' floods yielded water heights of 0.9 m, 1.1 m, 1.3m, 1.33 m, 1.38 m, 1.53 m, and 1.8 m, respectively.The flow characteristics at the 26 (Ganga Uku), 22 (NUGA gate culvert), 8 (Dan Fodio culvert), and 1 (Detention basin) stations for the 100-year Tr profile are shown in Table 6.The longitudinal profile shows the variation of the water surface profile for each cross-section in the channel.Figure 9 shows the longitudinal water surface profiles for the floods in a 100-year period.The profile varies significantly from section to section, as seen from the figures, and this is evident because the natural streams usually have varying cross-sections.The longitudinal profile shows the variation of the water surface profile for each cross-section in the channel.Figure 9 shows the longitudinal water surface profiles for the floods in a 100-year period.The profile varies significantly from section to section, as seen from the figures, and this is evident because the natural streams usually have varying cross-sections.

(c) Water Depth in Overbanks
The water depth in the right and left overbanks for the different return periods is presented in Figure 10.The maximum water depth of the right overbank is about 0.9 m, while it is 1.0 m in the left overbank.Furthermore, it was observed that the water depth in both overbanks increased with an increase in the return period.

(c) Water Depth in Overbanks
The water depth in the right and left overbanks for the different return periods is presented in Figure 10.The maximum water depth of the right overbank is about 0.9 m, while it is 1.0 m in the left overbank.Furthermore, it was observed that the water depth in both overbanks increased with an increase in the return period.

Smart Flood Control Structures
The 100-year design flood was selected along with the other basic data required for hydraulic design, such as the channel cross-section and flow characteristics.To safely convey the peak runoff from the Samaru stream watershed to the stream outlet, a hydraulically economical trapezoidal channel of a constant cross-section was designed, which was achieved by considering the peak runoff of a 100-year return period at the Dan Fodio culvert junction.The HEC-HMS model output gave the design peak discharge at the Dan Fodio culvert as 31.8m 3 /s.The mean slope of the channel was obtained from the topographic survey as 0.0154.Chow [55] proposed a Manning's coefficient for concrete as 0.012.
The channel was found to be stable after checking the value of the mean velocity to be 2.86 m/s, which falls below the allowable velocity limit for a concrete lining.Hence, the channel was deemed stable against erosion for all flow durations, which meant that the velocity was lower than the allowable velocity.The summary of the design channel is presented in Table 7.

Conclusions
GIS-based techniques were employed to estimate flood events in the micro-watershed of Samaru via multi-model inference.The HEC-HMS model was developed to evaluate the effect of rainfall on the peak discharges of different return periods, while the HEC-RAS model was used to analyze the water surface profiles of different return periods and their characteristics so as to develop flash flood adaptation and mitigation measures.The following conclusions were drawn from the study.
(a) The hydrological modeling of the watershed was achieved using HEC-HMS.The intensity duration frequency curve (IDF) of the watershed was successfully developed and used for peak discharge estimation.The calibration and validation of the HEC-HMS results showed a close match between the observed and simulated flows in August and September 2014.(b) The calibration performance evaluation gave an RMSE, NSE, and PBIAS of 0.6, 0.67, and 20.54%, respectively, while the validation performance evaluation produces were 0.50, 0.783, and 20.77%, respectively.(c) A parameter sensitivity analysis of the HEC-HMS model revealed that the model outputs were more sensitive to the curve number.This may be attributed to the land use changes in the watershed, which led to increasing urbanization and paved surfaces; thus, land development in the watershed should be monitored.(d) The peak runoff rates from the Samaru watershed for 2-, 5-, 10-, 25-, 50-, 100-, and 200-year return periods at all junctions and reaches along the Samaru stream, including the watershed outlet, were successfully estimated.The peak discharges at the NUGA gate culvert for 2-, 5-, 10-, 25-, 50-, 100-, and 200-year return periods were estimated as 3.7, 6.8, 9.1, 12.1, 14.3, 16.6, and 19.0 m 3 /s, respectively, while those at the watershed outlet were found to be 5.7, 11.0, 14.9, 20.0, 23.9, 27.8, and 31.8 m 3 /s, respectively.(e) The hydraulic modeling of the Samaru stream was performed using the HEC-RAS model with river cross-sections; the flow data from the 2-, 5-, 10-, 25-, 50-, 100-, and 200-year return periods; and Manning's coefficient as the input data.The modeling analysis yielded water depth, water velocity, and other flow characteristics at all of the cross-sections that were entered for the respective return periods.(f) The floods of the 25-year return period and above produced water surface profiles that inundated the overbanks at cross-sections 22, 24, 9, and 8. Furthermore, the NUGA gate culvert was also found to be overtopped by floods in 25-year return periods and above.(g) The maximum water depth in the right and left overbanks was found to be about 0.9 m and 1.0 m, respectively.This was used for channel design, which was achieved by considering the floods in a 100-year return period.The most economical channel section design gave approximate channel dimensions of 1.50 m as the normal flow depth, 1.73 m as the bottom width, 3.50 m as the channel top width, and a free board of 0.30 m.(h) These findings enable an improved understanding of parameter behavior in the HEC-HMS model application, as well as aids in the identification of the dominant processes under different hydrological scenarios that will allow for the accurate estimation of floods in micro-watersheds.

Figure 1 .
Figure 1.Location of the Samaru stream watershed on a Nigerian map.The sub-catchment is a 1st-order minor tributary of the Kubanni River within the Kubanni drainage basin, lying between 110 • 08 32 and 110 • 09 38 N, and 70 • 38 36 and

Figure 2 .
Figure 2. Characteristics of the Samaru watershed.(a) Land use of the study area, (b) slope, (c) hydrological soil map, and (d) curve number grid.

Figure 3 .
Figure 3. Flow chart of the methodology.

Figure 4 .
Figure 4.The IDF curve of the Samaru watershed.

Figure 5 .
Figure 5. Sensitivity analysis on a 100−year return period simulation.

Figure 7 .
Figure 7. Flood hydrograph of the 100−year Tr at (a) the NUGA gate culvert and (b) the watershed outlet.
(a) Cross-sectional Water Surface ProfileThe cross-sectional water surface profile at the NUGA gate culvert, Dan Fodio culvert, downstream of the NUGA gate culvert, and the detention basin are shown in Figure8.As seen from the cross-section at the NUGA gate culvert (XS-22), virtually all of the flows of different return periods produced a water surface profile that overtopped the main channel and the overbanks.Hence, this signified a flash flood or overtopping of the culvert.

Figure 8 .
Figure 8. Cross-sectional profile at the NUGA gate culvert, Dan Fodio culvert, downstream of the NUGA gate Culvert, and the detention basin.

Water 2023 ,
15, x FOR PEER REVIEW 15 of 20 (b) Longitudinal Water Surface Profile

Figure 9 .
Figure 9. Longitudinal water surface profile for floods in a 100-year period.

1 Figure 9 .
Figure 9. Longitudinal water surface profile for floods in a 100-year period.

Figure 10 .
Figure 10.Water depth in (a) the right overbank and (b) the left overbank.

Table 1 .
Data type and sources of the data.

Table 2 .
Intensity-frequency and duration values for 24 h events.

Table 3 .
Initial and optimized parameters of the HEC-HMS model.

Table 4 .
Performance indices during verification.

Table 5 .
Peak discharges at a location of interest in m 3 /s.

Table 6 .
Flow characteristics at some selected river stations.

Table 7 .
Summary of the design channel parameters.