A periodical of the Faculty of Natural and Applied Sciences, UMYU, Katsina
ISSN: 2955 – 1145 (print); 2955 – 1153 (online)
ORIGINAL RESEARCH ARTICLE
Kehinde Raheef Adebayo1*, Buhari Olugbon Yusuf1, Kamil Kayode Katibi1, Azeez Ayinla Adebayo1 and Ademola Lawrence Olaoluwa2
1Department of Food and Agricultural Engineering, Faculty of Engineering and Technology, Kwara State University, Malete, Nigeria
2Department of Agricultural and Bioresources Engineering, College of Engineering, Federal University of Agriculture, Abeokuta, Nigeria
*Corresponding Author: Kehinde Raheef Adebayo
kehinde.adebayo@kwasu.edu.ng
This study estimated unit hydrograph ordinates and developed runoff hydrographs for the Foma River watershed. The Snyder’s and SCS synthetic approaches were used to produce the unit hydrograph ordinates, while the SCS type II curve was used to determine the cumulative excess rainfall values for storm depths with 20-yr and 50-yr return periods. To develop the peak runoff hydrograph, the synthetic unit hydrograph ordinates were convoluted with the cumulative excess rainfall for the 20-yr and 50-yr storm return periods, using the Gumbel Extreme Value Type I (EV-I) probability distribution. Peak runoff values obtained based on the SCS for the watershed ranged from 68.0 m3/s to 112.0 m3/s for the return periods of 20 and 50 years, while peak runoff values obtained based on the Snyders ordinates for the 20- and 50-year return periods varied between 70.0 m3/s and 90.0 m3/s. The findings demonstrate that the SCS approach was the most effective for calculating the flow ordinate required to develop the peak runoff hydrograph for various return periods in the research area.
Keywords: Peak runoff, Return period, Storm, Unit hydrograph, Excess rainfall
-This study demonstrates a novel framework for improved accuracy and reliability in predicting peak runoff rates and hydrograph shapes.
-The application of the Soil Conservation Service (SCS) and Snyder’s Unit hydrograph synthetic methods was used for the study.
- The results reveal that the Soil Conservation Service (SCS) method was the best for the development of peak runoff hydrograph compared to Snyder’s approach.
- This enables more effective flood risk assessment, stormwater management, and water resources planning to enhance resilience and sustainability in hydrological design and management practices in the Foma River watershed, Ilorin, Kwara State.
A unit hydrograph of a basin or watershed is rarely determined using rainfall and runoff data in most Nations worldwide (Yi et al, 2022; Prakash et al., 2025). The absence of gauging stations along most of Nigeria's rivers and streams makes this situation common. Planning and developing water management facilities and other hydraulic infrastructure in undeveloped watersheds often lack fundamental streamflow and rainfall data (Prakash et al. 2025). Nevertheless, methods for creating synthetic unit hydrographs have been developed. These include the Soil Conservation Service (SCS) approach, Snyder's method, Gray's method, and Clark's instantaneous method. Design storm hydrographs are created from unit hydrographs produced using recognized techniques applied to determine the peak discharges of stream flow caused by rainfall.
The 1-hour unit hydrograph, as described by Arora (2004), is the hydrograph that provides a 1 cm depth of direct runoff during a 1-hour storm that occurs evenly throughout the watershed. A hydrograph is a continuous graph that displays the characteristics of stream flow over time. It is typically created using a continuous strip recorder that plots stages against time. A rating curve is then used to convert the hydrograph into a discharge hydrograph. A drainage basin's unit hydrograph is the hydrograph of direct runoff produced uniformly over the basin area at a consistent rate after one unit of effective rainfall of a certain period (WMO, 2009a). The basin's area, slope, orientation, form, altitude, and stream pattern are all considered watershed qualities.
Ogunlela (1996) developed a unit hydrograph for a small agricultural watershed at the University of Ilorin, accounting for the watershed's storage characteristics by routing through an assumed linear reservoir using Clark's approach. At a time to peak of 0.33 hours, he achieved a unit hydrograph peak of 2.97 m3/s. In contrast, he obtained peak flows of 4.53 m3/s (at 0.58 hours) and 6.23 m3/s (at 0.58 hours) for the 25-year, 24-hour and 100-year, 24-hour storm hydrographs. Ayansola and Salami (2009) developed a unit hydrograph for the Awun River basin using the Snyder, SCS, and Grey methods; the peak unit hydrograph values were 299.27 m3/s, 307.28 m3/s, and 2083.40 m3/s, respectively. Salami (2009) assessed storm hydrograph techniques for the Lower Niger River basin watershed downstream of Jebba Dam. With the exception of the Snyder and SCS procedures, which have comparable values, the statistical analysis at the 5% level of significance revealed considerable disparities among the three methods under consideration: Snyder, SCS, and Grey.
Synthetic unit hydrograph methods are widely used for peak discharge estimation in ungauged basins where observed runoff records are unavailable. Recent studies have demonstrated the continued relevance of classical approaches, such as the Soil Conservation Service (SCS) and Snyder methods, particularly when combined with GIS-based extraction of watershed parameters (Patil et al., 2023; Kesgin, 2025). Comparative evaluations in poorly gauged and ungauged basins have shown that these methods remain effective for preliminary flood assessment and hydraulic design applications (Al-Dughairi, 2023; Casado and López, 2025).
Foma River is located approximately 7 kilometres from the Emir’s Palace in Ilorin, Nigeria, at approximately 8.5° N and 4.5° E, as shown in Figure 1 (Remote Sensing in Earth Systems Sciences, 2023). Ilorin has a tropical savanna climate characterized by distinct wet and dry seasons. Mean temperatures range between about 34 °C from November to January and increase to approximately 36 °C between February and April (Ajadi et al., 2016; NiMet, 2020). The total annual rainfall in the area varies from approximately 990.3 mm to 1318 mm (Ajadi et al., 2016; Olaniran, 2002). The river flows freely during rainfall and slows during the dry season (Ayansola and Salami, 2009). The feasibility study indicates that the river has its origin at Wahri, along Gerewu, Ilorin West Local Government Area of Kwara State. It joined the Asa River at a point called Ibu-Afonja, beneath the Sobi Rock, in Ilorin South, Kwara State. The river’s catchment area is calculated as 30.31 km2 (Figure 1). The river, especially during the rainy season, overflows its banks, which, over the years, have resulted in several kilometres of floodplains. Table 1 shows the watershed characteristics of the Foma river.
Table 1: Catchment characteristics of Foma river
| Watershed | L (km) | Lc (km) | A (km2) | Sc (m/m) |
|---|---|---|---|---|
| Foma river | 10.65 | 4.20 | 30.31 | 0.00722 |
Where;
L = The river channel length (Km)
Lc = The river length from the outlet to a point near the centroid (Km)
A = Watershed area (Km2)
Sc = Slope of river channel (m/m)
Figure 1: The Study Area
The Soil Conservation Service (SCS) and Snyder's methodologies were the two approaches used to develop the synthetic peak runoff hydrograph. The SCS and Snyder synthetic unit hydrograph methods were adopted in this study due to their widespread application and proven suitability in ungauged watersheds. Recent applications have demonstrated their robustness in estimating peak runoff hydrographs using basin morphometric characteristics derived from GIS data (Patil et al., 2023; Kesgin, 2025). Similar methodological frameworks have been successfully employed in arid and semi-arid ungauged basins, yielding reliable peak discharge estimates for design storm events (Al-Dughairi, 2023).
In accordance with Viessman et al. (1989), the Soil Conservation Service (SCS) used a dimensionless hydrograph as the basis for their method of creating synthetic unit hydrographs. This method involved determining the catchment slope S, concentration time tc, peak time tp, and peak flow Qp. As previously stated, the catchment properties encompass the basin's area, slope, orientation, shape, altitude, and stream pattern. The annual peak daily rainfall event over the watershed was applied to develop runoff hydrographs using the estimated synthetic unit hydrograph method under the SCS approach. The parameters listed in Table 2 for creating a unit hydrograph using the SCS approach are shown in Table 3.
Table 2: The parameters applied for the generation of the unit hydrograph (SCS Method)
| River Watershed | L (km) | A (km2) | SC (%) | TC (hr) | TL (hr) | Tp (hr) | QP (m3/s) |
|---|---|---|---|---|---|---|---|
| Foma River | 10.65 | 30.31 | 0.00722 | 8.38 | 5.026 | 5.59 | 11.28 |
Table 3: Unit hydrograph by SCS Method
| T (hr) | 0.00 | 2.79 | 5.59 | 8.39 | 11.18 | 13.98 | 16.77 | 19.57 | 20.36 | 25.16 | 27.95 |
|---|---|---|---|---|---|---|---|---|---|---|---|
| Q’p(m3/s) | 0.00 | 4.85 | 11.28 | 7.44 | 3.61 | 1.47 | 0.85 | 0.41 | 0.20 | 0.10 | 0.05 |
Equations 1 to 4 were applied to calculate rainfall excess, and additional parameters employed included rainfall depth at varying return periods and the curve number (CN), with reference to the research area's soils and land use (Plamonia et al. 2025).
\(Q = \frac{{(R^{**} - Ia)}^{2}}{(R^{**} + 0.8S)}\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \) (1)
\(I_{a} = 0.2*S\) \(\ \ \ \ \ \) (2)
\(R^{**}\ = \frac{R^{*}}{24}*R_{T}\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \) (3)
\(S = \frac{25400}{CN} - 254\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \) (4)
Where,
\(R^{**}\)= The accumulated rainfall (mm)
RT = The catchment’s area rainfall recurrence interval (mm)
R*= The precipitation ratio, \(\frac{P_{x}}{P_{24}}\)
S = Total storage volume (mm)
CN = Runoff curve number (Based on the hydrologic soil group B and land use (Agriculture), the CN value of 75 was chosen (Plamonia et al. 2025).
Tables 4 and 5 present the estimations of rainfall excess at varying return periods when\(\ CN\), \(S,\) and \(I_{a}\) are 75, 84.67 and 16.93, respectively.
Table 4: Estimated rainfall excess of Foma River for 20yr, 24-hr storm at P =174.17 mm
| Time (hr) | Precipitation Ratio (P*/24) | Precipitation P** (mm) | Cumulative Rainfall Excess Qd (mm) | Incremental Rainfall Excess (mm) | |
|---|---|---|---|---|---|
| (20-yr., 24-hr) | 174.17 | mm | cm | ||
| 0 | 0.0000 | 0.0000* | 0.0000 | 0.0000 | 0.0000 |
| 3 | 0.0350 | 6.0960* | 0.0000 | 0.0000 | 0.0000 |
| 6 | 0.0800 | 13.9336* | 0.0000 | 0.0000 | 0.0000 |
| 9 | 0.1470 | 25.6030 | 0.0000 | 0.0000 | 0.0000 |
| 12 | 0.6630 | 115.4747 | 52.9930 | 52.9930 | 5.2993 |
| 15 | 0.8540 | 148.7412 | 80.2451 | 27.2521 | 2.7252 |
| 18 | 0.9210 | 160.4106 | 90.2203 | 9.9752 | 0.9980 |
| 21 | 0.9650 | 168.0410 | 96.8623 | 6.6421 | 0.6642 |
| 24 | 1.0000 | 174.1700 | 102.1920 | 5.3287 | 0.5330 |
* P** < Ia
Table 5: Estimated rainfall excess of Foma River for 50yr, 24-hr storm at P =205.36 mm
| Time (hr) | Precipitation Ratio (P*/24) | Precipitation P** (mm) | Cumulative Rainfall Excess Qd (mm) | Incremental Rainfall Excess (mm) | |
|---|---|---|---|---|---|
| (50-yr, 24-hr) | 205.36 | Mm | cm | ||
| 0 | 0.0000 | 0.0000 | 0.0000 | 0.0000 | 0.0000 |
| 3 | 0.0350 | 7.1876 | 0.0000 | 0.0000 | 0.0000 |
| 6 | 0.0800 | 16.4288 | 0.0000 | 0.0000 | 0.0000 |
| 9 | 0.1470 | 30.1879 | 0.0000 | 0.0000 | 0.0000 |
| 12 | 0.6630 | 136.1537 | 69.7025 | 69.7025 | 6.9703 |
| 15 | 0.8540 | 175.3774 | 103.2522 | 33.5493 | 3.3550 |
| 18 | 0.9210 | 189.1366 | 115.4315 | 12.1793 | 1.2179 |
| 21 | 0.9650 | 198.1724 | 123.5188 | 8.0873 | 0.8087 |
| 24 | 1.0000 | 205.3600 | 128.9967 | 6.4779 | 0.6478 |
Snyder's approach uses the watershed's specific features to estimate the peak discharge, lag time, and time to peak. According to Viessman et al. (1989), the hydrograph's features include the lag time (TL), peak discharge (Qp), and effective rainfall duration (Tr). These correlations allow the calculation of the five required unit hydrograph features for a specific effective rainfall duration. The basin lag T'L, the base time T'b, the peak discharge per unit of watershed area Q'p, and the widths W of the unit hydrograph at 50 and 75 % of the peak discharge are the five characteristics. Ayansola and Salami (2009) used these relationships to generate five characteristics of a unit hydrograph required for a given effective rainfall duration using equations 5–12.
Estimation of the Lag time (TL)
The lag time TL calculated through the use of equation (5) as stated below;
\(T_{L\ } = C_{t\ }\left( L*L_{c} \right)^{0.3}\) (5)
Where;
TL is the lag time (hr), and Ct is a coefficient that denotes variations in watershed slope and storage.
The values of Ct range from 1.0 to 2.2 (Arora, 2004) and the mean value of 1.60 was applied.
Estimation of Unit-hydrograph duration, Tr (Storm duration)
The storm's duration was estimated using equation (6).
\(T_{r} = \ \frac{t_{l}}{5.5}\ \) (6)
The new unit hydrograph storm duration (\(t_{r}'\)) and the accompanying basin lag time (\(t_{l}'\))), however, can be calculated from equation (7) in the event that longer storm durations are planned to be achieved for the catchment.
\(t_{r}' = t_{l} + \ \) \(\left( \frac{t_{l} - \ t_{r}}{4} \right)\) (7)
Estimation of peak discharge, Q’p
The peak discharge (\(Q_{p}'\)) can be derived from equation 8 below;
\(Q_{p}' = \ \frac{2.78*\ C_{p}*A}{t_{l}'}\) (8)
Where;
Cp is the coefficient representing the flood wave and storage condition ranges from 0.30 to 0.93 (Arora, 2004) with an average of 0.62.
Estimation of Base time, Tb
Equation (9) was applied for calculating the base time (Tb) in days;
\(T_{b\ \ } = 3\ + 3\ \left( \frac{t_{l}'}{24} \right)\) (9)
Estimation of Time width, W
The time width W50 and W75 of the hydrograph at 50% and 75% of the height of the peak flow ordinate in an hour, estimated from equations 10 and 11, respectively, by the US Army Corps of Engineers (Arora, 2004).
\(W_{50\ } = \ \frac{5.9}{\left( {q'}_{p} \right)\ 1.08}\) (10)
\(W_{75\ } = \ \frac{3.4}{\left( {q'}_{p} \right)\ 1.08}\) (11)
Estimation of Peak discharge per area, q’p
Equation 12 provides the peak discharge per area (cumec/km2);
\({q'}_{p}\ = \frac{{Q'}_{p}}{A}\) (12)
These parameters were substituted to derive the unit hydrograph and results obtained by Snyder’s method were presented in Tables 6 and 7.
Table 6: Parameters for the generation of unit hydrograph (Snyder’s method)
| River Watershed | L (km) |
LC (km) |
A (km2) |
Tl (hr) |
Tr (hr) |
T’l (hr) |
Q’p (m3/s) |
Tb (hr) |
W50 | W75 |
|---|---|---|---|---|---|---|---|---|---|---|
| Foma River | 10.65 | 4.20 | 30.31 | 5.00 | 0.91 | 5.03 | 10.40 | 87.12 | 34.7 | 106.68 |
Table 7: Unit hydrograph by Snyder’s Method
| T (hr) | 0.00 | 2.57 | 5.03 | 8.34 | 12.30 | 15.92 | 17.34 | 21.08 | 28.01 | 30.32 | 30.70 |
|---|---|---|---|---|---|---|---|---|---|---|---|
| Q’p(m3/s) | 0.00 | 5.12 | 10.40 | 6.02 | 5.38 | 4.09 | 3.93 | 2.04 | 1.24 | 0.13 | 0.00 |
The convolution equations developed by using equation 16 to compute the direct runoff (\(Q_{n}\)) by multiplying unit hydrograph ordinates from the catchment (U) and incremental rainfall excess (R), adding and lagging them in sequence as stated in equations 13-16, are used to develop the design runoff hydrographs for selected rainfall events with recurrence intervals of 20 and 50 years.
The direct runoff hydrograph's first ordinate results solely from R1's effective rainfall pulse (Chow et al., 1988; Subramanya, 2013).
\(Q_{1}\ = R_{1}U_{1}\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \) (13)
The second ordinate is calculated from the effective rainfall of R1 and R2
\(Q_{2}\ = \ R_{2}U_{1}\ + R_{1}\ U_{2}\) (14)
The third ordinate is given by
\(Q_{3}\ = \ R_{3}U_{1}\ + R_{2}\ U_{2}\ + R_{1}U_{3}\) (15)
The excess rainfall ordinates are R1, R2, R3, R4 and R5 while the unit hydrograph ordinates are U1, U2, U3, U4, U5, U6, U7, U8, U9, U10 and U11 combined to generate equation (16).
Generally, the peak runoff ordinate is obtained using equation 8;
\(Q_{n}\ = \ R_{1}U_{n}\ + R_{2}\ U_{n - 1}\ + R_{3}U_{n - 2}\ + R_{4}U_{n - 3}\ + R_{5}U_{n - 4}\) (16)
Where;
R = incremental rainfall excess (cm),
U = unit hydrograph ordinate (m3/s/cm).
Rainfall depths corresponding to selected return periods were estimated using the Gumbel Extreme Value Type I (EV-I) probability distribution as depicted in equations 17 and 18. The distribution was applied to annual maximum rainfall series to determine design rainfall magnitudes for different recurrence intervals. The estimated rainfall depths were subsequently used as inputs for the computation of rainfall excess, which served as the basis for runoff and unit hydrograph analysis.
Rainfall depths corresponding to selected return periods were estimated using the Gumbel Extreme Value Type I distribution, which is commonly applied in rainfall and flood frequency analysis for design-oriented applications when data availability is limited (Badou et al., 2021).
\(X_{T} = \ \overline{X} + K_{T\ }\sigma\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \) \((17)\)
where:
\(X_{T}\) - rainfall magnitude corresponding to return period T,
\(\overline{X}\) - mean of the annual maximum rainfall series,
\(\sigma\) - standard deviation of the series,
\(K_{T\ }\)- Gumbel frequency factor.
\(K_{T\ } = \ \frac{\sqrt{6}}{\pi}\left\lbrack \ln\left( \ln\left( \frac{T}{T - 1} \right) \right.\ \right\rbrack\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \) \((18)\)
where:
\(T\) = return period (years).
The peak runoff hydrograph generated for the Foma river watershed for 20-yr and 50- yr return periods using SCS and Snyder’s synthetic unit hydrograph approaches is depicted in Figures 3 and 4. These show the variations between the SCS and Snyder’s method. The Gumbel Extreme Value Type 1 distribution was applied to obtain storm depth values for 20-yr and 50-yr return periods, with values of 174.17 mm and 205.36 mm, respectively.
Figure 3: Runoff hydrographs of different return periods generated from unit hydrograph (SCS method)
Figure 4: Runoff hydrographs of different return periods generated from unit hydrograph (Synder method)
Figure 5: Peak runoff hydrographs for the 20-year return period of SCS and Snyder synthetic unit hydrograph methods
Figure 6: Peak runoff hydrographs for the 50-year return period of SCS and Snyder synthetic unit hydrograph methods
Figure 3 shows that the peak runoff hydrograph estimate occurred over a short duration, ranging from 68.0 m3/s to 112.0 m3/s, for the Foma River using the SCS method. Figure 4 indicates that the peak runoff hydrograph using Snyder’s method ranged from 70.0 m3/s to 90.0 m3/s for the Foma River. From the above, it is shown that the peak runoff hydrograph estimate for the SCS method is close to that of Snyder's method. The results also indicated that peak runoff occurred over a short duration and at a low discharge magnitude for the Foma River.
The 20-yr, 24-hr storm hydrograph discharges are 68.0 m3/s and 70.0 m3/s for the Foma River using both the SCS and Snyder methods. The percentage difference between the two estimates was approximately 2.86%, indicating close agreement between the methods for moderate storm events. This suggests that both approaches yield comparable peak runoff estimates under lower-return-period conditions (Figure 5). The 50-yr, 24-hr storm hydrograph discharges were 90.0 m3/s and 112.0 m3/s for the Foma river for both SCS and Snyder’s methods, also for the same return period; however, the SCS method produced a significantly higher peak discharge of 112.0 m³/s compared to 90.0 m³/s obtained using the Snyder method, representing a percentage difference of approximately 24.4%. This divergence indicates that the SCS method is more sensitive to increases in storm magnitude and return period, producing sharper and higher peak flows than the Snyder approach (Figure 6).
Figures 5 and 6 present schematic runoff hydrographs constructed for comparative purposes using the estimated peak discharges and time-to-peak values derived from the SCS and Snyder methods. These hydrographs are intended to illustrate relative differences in hydrograph shape, peak magnitude, and timing between the two approaches align with WMO (2009b) guidance for ungauged basins.
The stronger response of the SCS method can be attributed to its explicit consideration of watershed land use and soil characteristics through the Curve Number (CN) parameter, which directly influences rainfall excess generation. In contrast, the Snyder method relies on regional coefficients and basin geometry, resulting in a comparatively attenuated hydrograph response, particularly for extreme rainfall events similar to the findings reported by Salami (2009) and Ogunlela and Kasali (2002), who observed that SCS-based hydrographs generally produce higher peak flows than Snyder-based hydrographs for the same watershed conditions.
Overall, the results indicate that the SCS synthetic unit hydrograph method is more responsive to high-magnitude storm events and may therefore be more suitable for flood estimation and hydraulic structure design in the Foma River watershed, while the Snyder method remains useful for preliminary assessments and comparative analysis. The comparative hydrographs generated using the SCS and Snyder methods indicate close agreement for the 20-year return period, while greater divergence is observed for the 50-year storm event. Similar findings have been reported in recent comparative studies, where SCS-based methods tend to produce higher peak discharges under extreme rainfall conditions due to their explicit consideration of land use and soil characteristics (Casado and López, 2025; Kesgin, 2025). In contrast, the Snyder method exhibits a more attenuated runoff response, consistent with observations in other ungauged basin studies (Al-Dughairi, 2023).
The observed consistency of both methods across return periods supports previous findings that SCS and Snyder synthetic unit hydrographs remain suitable for runoff estimation in ungauged basins, with method selection largely dependent on storm magnitude and design objectives (Patil et al., 2023; Casado and López, 2025).
The SCS methodology yielded greater peak discharge values than Snyder's method at 50-yr, 24-hr, indicating a substantial difference between the two approaches (112 m3/s > 90.0 m3/s). In contrast, the SCS method's peak discharge at 20-yr, 24-hr was lower than the Synder method's (68.0 m3/s < 70.0 m3/s) result. In order to construct hydraulic structures, both the SCS and Snyder synthetic unit hydrograph methods demonstrated hydrologically consistent responses across short and long return periods, with peak discharge increasing monotonically with storm severity. The negligible difference (2.86%) in peak discharge estimates for the 20-year storm indicates statistical equivalence of the two methods for moderate events, while the divergence observed at the 50-year return period reflects differences in method sensitivity rather than methodological inadequacy. Consequently, both approaches are suitable for estimating watershed parameters required for unit hydrograph development in the ungauged catchment, with method selection dependent on the design return period and application objective.
I acknowledged the efforts of the students and colleagues that contributed to the study within the Faculty of Engineering and Technology.
Ajadi, M. S., Abiodun, B. J., & Oguntunde, P. G. (2016). Climate variability and change in Ilorin, Nigeria. Journal of Environmental Geography, 9(3–4), 23–33. [Crossref]
Al-Dughairi, A. A. (2023). Runoff hydrographs using Snyder and SCS synthetic unit hydrograph methods: A case study of ungauged watersheds of Imam Turki Bin Abdullah Royal Nature Reserve (ITBA), Saudi Arabia. Journal of Natural Sciences, Life and Applied Sciences, 7(3), 79–94. [Crossref]
Arora, K. R. (2004). Irrigation, water power and water resources engineering. Standard Publishers Distributions.
Ayanshola, A. M., & Salami, A. W. (2009). Evaluation of methods of storm hydrograph development for Asa River in Kwara State, Ilorin, Nigeria. In 1st Annual Civil Engineering Conference, University of Ilorin (pp. 25–34).
Badou, D. F., Adango, A., Hounkpè, J., Bossa, A., Yira, Y., Iboukoun Biao, E., & Afouda, A. A. (2021). Heavy rainfall frequency analysis in the Benin section of the Niger and Volta Rivers basins: Is the Gumbel distribution a one-size-fits-all model? Proceedings of IAHS, 384, 187–194. [Crossref]
Casado, A., & López, N. C. (2025). Comparison of synthetic unit hydrograph methods for flood assessment in a dryland, poorly gauged basin (Napostá Grande, Argentina). AIMS Geosciences, 11(1), 27–46. [Crossref]
Chow, V. T., Maidment, D. R., & Mays, L. W. (1988). Applied hydrology. McGraw-Hill.
Kesgin, E. (2025). Flood discharge estimation in ungauged basins using synthetic unit hydrographs and GIS. Firat University Journal of Experimental and Computational Engineering, 4(2), 375–392. [Crossref]
Nigerian Meteorological Agency (NiMet). (2020). Climate review bulletin for Nigeria. Nigerian Meteorological Agency.
Ogunlela, A. O. (1996). Synthetic hydrograph development for a small agricultural watershed. In Civil Engineering Impact on National Development (pp. 72–83). Faculty of Engineering and Technology, University of Ilorin.
Ogunlela, A. O., & Kasali, M. Y. (2002). Evaluation of four methods of storm hydrograph development for an ungagged watershed. Nigerian Journal of Technological Development, 2(1), 79–106.
Olaniran, O. J. (2002). Rainfall anomalies in Nigeria: The contemporary understanding. Journal of Geography and Regional Planning, 25(1), 1–10.
Patil, S., Gawande, S., & Waghmare, G. (2023). Synthetic unit hydrograph for ungauged basin using Snyder, Taylor-Schwarz model, and SCS method by GIS techniques. Materials Today: Proceedings, 77(Pt. 3), 855–859. [Crossref]
Plamonia, N., Aryani, D., Juniati, A. T., Seknun, N., Putra, A. P., Wahyono, I. B., Zaidan, M., Winarno, B., Kurniawan, B., & Prihartanto. (2025). Analysis of hydrological characteristics and rainfall distribution patterns in the Cimanuk River Basin, West Java, Indonesia. International Journal of Distributed and Parallel Networks (IJDNE), 20(6), 1295–1312. [Crossref]
Prakash, H., Pandey, K. K., & Soni, P. (2025). Peak discharge estimation for ungauged basins: A review. Journal of Water and Climate Change, 16(11), 3483–3502. [Crossref]
Remote Sensing in Earth Systems Sciences. (2023). GIS-assisted flood risk potential mapping of Ilorin and its environs, Kwara State, Nigeria. Remote Sensing in Earth Systems Sciences, 6(2), 245–260. [Crossref]
Salami, A. W. (2009). Evaluation of methods of storm hydrograph development. International Egyptian Engineering Mathematical Society IEEMS, Zagazig University Publications, 6, 17–28.
Subramanya, K. (2013). Engineering hydrology (4th ed.). McGraw-Hill Education.
Viessman, W., Knapp, J. W., & Lewis, G. L. (1989). Introduction to hydrology. Harper and Row Publishers.
World Meteorological Organization (WMO). (2009a). Guide to hydrological practices: Volume II — Management of water resources and application of hydrological practices (WMO-No. 168). [Crossref]
World Meteorological Organization (WMO). (2009b). Guide to hydrological practices, volume I: Hydrology – from measurement to hydrological information (6th ed.). WMO-No. 168.
Yi, B., Chen, L., Zhang, H., Singh, V. P., Jiang, P., Liu, Y., Guo, H., & Qiu, H. (2022). A time-varying distributed unit hydrograph method considering soil moisture. Hydrology and Earth System Sciences, 26, 5269–5289.