Continuous Modeling of the Mkurumudzi River Catchment in Kenya Using the HEC-HMS Conceptual Model : Calibration , Validation , Model Performance Evaluation and Sensitivity Analysis

The Mkurumudzi River originates in the Shimba hills and runs through Kwale County on the Kenyan Coast. Study on this river has been informed by the many economic activities that the river supports, which include sugarcane plantations, mining, tourism and subsistence farming. The main objective of this study was to use the soil moisture accounting (SMA) model specified in the Hydrologic Engineering Center-Hydrologic Modeling System (HEC-HMS) settings for the continuous modeling of stream flow in the Mkurumudzi catchment. Data from past years were compared with observed stream flow data in order to evaluate whether the model can be used for further prediction. The calibration was performed using data from 1988 to 1991 and validation for the period from 1992 to 1995 at a daily time step. The model performance was evaluated based on computed statistical parameters and visual checking of plotted hydrographs. For the calibration period of the continuous modeling, the performance of the model was very good, with a coefficient of determination R2 = 0.80, Nash-Sutcliffe Efficiency NSE = 0.80, index of agreement d = 0.94, and a Root Mean Squared Error (RMSE)/observations’ standard deviation ratio—RSR = 0.46. Similarly, the continuous model performance for the validation period was good, with R2 = 0.67, NSE = 0.65, RSR = 0.62 and d = 0.88. Based on these performance results, the SMA model in the HEC-HMS was found to give a satisfactory prediction of stream flow in the Mkurumudzi Catchment. The sensitivity analysis of the model parameters was performed, and the different parameters were ranked according to their sensitivity in terms of percent change in simulated runoff volume, peaks, Nash-Efficiency, seven-day low flow and base flow index. Sensitivity analysis helped to understand the relationships between the key model parameters and the variables.


Introduction
The hydrology of any catchments depends highly on climatic conditions, especially rainfall.Hydrological modeling helps in the understanding of all the processes of rainfall transformation into runoff.The Hydrologic Modeling System (HEC-HMS), developed by the US Army Corps of Engineers Hydrologic Engineering Center in 1998, is designed for both continuous and event-based hydrologic modeling [1].It offers different options for modeling various components of the hydrologic Hydrology 2018, 5, 44; doi:10.3390/hydrology5030044www.mdpi.com/journal/hydrologycycle.For watershed modeling, the HEC-HMS model contains four components: (1) Basin component, (2) Meteorological component, (3) Control specification, and (4) Input data component (time series, paired data and gridded data).For effective use of hydrological models, calibration, validation, performance evaluation and sensitivity analysis need to be undertaken.After calibration is done using one part of the data, the validation process follows in order to analyze the performance of the simulation using the remaining part of the data [2].Then sensitivity analysis is done to find the most sensitive parameters, so as to understand how great the impact of each model parameter is on the results.This paper addresses three specific objectives: (1) calibration and validation of the HEC-HMS hydrologic model in the Mkurumudzi catchment, (2) Model Performance Evaluation by statistical measures, and (3) sensitivity analysis of the model's parameters.

Materials and Methods
The study site, the Mkurumudzi River basin, measures approximately 230 km 2 , is located 50 km south of Mombasa, Kenya, and stretches 30 km inland from the Indian Ocean.The river, which is approximately 40 km long, is the main waterway flowing from the Shimba Hills and drains into the Indian Ocean.The river basin experiences a sub-humid climate, with average rainfall ranging from 800 mm to 1300 mm in the short (October-December) and long (March-May) rainfall periods, respectively, based on records from 1959 to 2012.Mean annual evaporation is about 2170 mm/year, giving an aridity index of approximately 0.55 [3].The mean annual minimum and maximum temperatures are 22.8 and 30.0 • C, respectively.The warmest months are between November and April, with mean temperatures of 26.0 to 28.0 • C, while the cooler months have a temperature ranging between 24.0 and 26.0 • C. The mean daily evapotranspiration rate is 4.4 mm/day.The Mkurumudzi River has a gauging site (RGS-3KD06) near the Shimba Hills.The waterway has stable dry season base flows, which are sustained by groundwater.The catchment land is mainly used for agriculture (maize, beans, cowpeas, millet and sorghum, okra, cassava) (Figure 1) and livestock husbandry, commercial mining (Base titanium Ltd., Ukunda, Kenya), commercial farming of sugarcane (by Kwale International Sugar Company KISCOL, Kwale, Kenya), tourism associated with the sea and the Shimba Hills National Reserve, and for fishing, mainly in the Indian Ocean.For effective use of hydrological models, calibration, validation, performance evaluation and sensitivity analysis need to be undertaken.After calibration is done using one part of the data, the validation process follows in order to analyze the performance of the simulation using the remaining part of the data [2].Then sensitivity analysis is done to find the most sensitive parameters, so as to understand how great the impact of each model parameter is on the results.This paper addresses three specific objectives: (1) calibration and validation of the HEC-HMS hydrologic model in the Mkurumudzi catchment, (2) Model Performance Evaluation by statistical measures, and (3) sensitivity analysis of the model's parameters.

Materials and Methods
The study site, the Mkurumudzi River basin, measures approximately 230 km 2 , is located 50 km south of Mombasa, Kenya, and stretches 30 km inland from the Indian Ocean.The river, which is approximately 40 km long, is the main waterway flowing from the Shimba Hills and drains into the Indian Ocean.The river basin experiences a sub-humid climate, with average rainfall ranging from 800 mm to 1300 mm in the short (October-December) and long (March-May) rainfall periods, respectively, based on records from 1959 to 2012.Mean annual evaporation is about 2170 mm/year, giving an aridity index of approximately 0.55 [3].The mean annual minimum and maximum temperatures are 22.8 and 30.0 °C, respectively.The warmest months are between November and April, with mean temperatures of 26.0 to 28.0 °C, while the cooler months have a temperature ranging between 24.0 and 26.0 °C.The mean daily evapotranspiration rate is 4.4 mm/day.The Mkurumudzi River has a gauging site (RGS-3KD06) near the Shimba Hills.The waterway has stable dry season base flows, which are sustained by groundwater.The catchment land is mainly used for agriculture (maize, beans, cowpeas, millet and sorghum, okra, cassava) (Figure 1) and livestock husbandry, commercial mining (Base titanium Ltd., Ukunda, Kenya), commercial farming of sugarcane (by Kwale International Sugar Company KISCOL, Kwale, Kenya), tourism associated with the sea and the Shimba Hills National Reserve, and for fishing, mainly in the Indian Ocean.

Datasets
The data used in this study includes: rainfall, maximum and minimum temperature, insolation, wind speed and solar radiation, data for evapotranspiration calculation (Penman-Monteith method [4]) and flow data.Rainfall data are from the Kidongo Park Gate station, which is coded 9439065 .Maximum and minimum temperature, insolation, wind speed and solar radiation data are from the NASA POWER dataset.Flow data are from the gauging station 3KD06 for the period 1956 to 2002.The study period is chosen as the period where data are not having missing gaps (Table 1).

Software
ArcGIS of ESRI was used to visualize and edit soil and land use data, while Microsoft Excel was used for data formatting and for evapotranspiration calculation using the Penman-Monteith

Datasets
The data used in this study includes: rainfall, maximum and minimum temperature, insolation, wind speed and solar radiation, data for evapotranspiration calculation (Penman-Monteith method [4]) and flow data.Rainfall data are from the Kidongo Park Gate station, which is coded 9439065 .Maximum and minimum temperature, insolation, wind speed and solar radiation data are from the NASA POWER dataset.Flow data are from the gauging station 3KD06 for the period 1956 to 2002.The study period is chosen as the period where data are not having missing gaps (Table 1).

Software
ArcGIS of ESRI was used to visualize and edit soil and land use data, while Microsoft Excel was used for data formatting and for evapotranspiration calculation using the Penman-Monteith equation.The Arc Hydro tool in ArcGIS [5] was used for watershed delineation based on a Digital Elevation Model (DEM) and river maps.The Hydrologic Modeling System (HEC-HMS) was used to simulate the hydrological processes in the basin.

HEC-HMS Model Application
Before the HEC-HMS project can be run, we must have the following four components: basin, meteorological, input data and control specification components.The basin component is created based on a background map containing the delineated sub-basins derived from the Arc Hydro tool 10.3 from ESRI.The meteorological component contains the observed precipitation and discharge data; the control specification component determines the simulation period and the time step to be used.The observed historical data of the Kidongo Park Gate rain station representing each sub-basin and one stream gauge station (3KD06) were used for model calibration and validation.Due to the location of the 3KD06 station, only the two (02) upstream sub-basins M1 and M2 (Figure 3) were used.A daily time step was used for the simulation, given the time interval of the available observed data.equation.The Arc Hydro tool in ArcGIS [5] was used for watershed delineation based on a Digital Elevation Model (DEM) and river maps.The Hydrologic Modeling System (HEC-HMS) was used to simulate the hydrological processes in the basin.

HEC-HMS Model Application
Before the HEC-HMS project can be run, we must have the following four components: basin, meteorological, input data and control specification components.The basin component is created based on a background map containing the delineated sub-basins derived from the Arc Hydro tool 10.The Soil Moisture Accounting (SMA) Loss Method included in HEC-HMS was employed to model infiltration losses combined with canopy and surface methods.The canopy is a component of the sub-basin which is intended to represent the presence of plants in the area.The canopy method specified in the HEC-HMS settings is mainly used for continuous simulations.All rainfall is trapped until the canopy storage capacity is filled.Subsequently, the excess precipitation falls on the soil surface after passing through the canopy.The surface method specified in the HEC-HMS settings represents the depression of the ground surface where water accumulates after the pores in the soil are filled to the soil's field capacity.The water on the surface then accumulates in the depression storage until runoff begins to occur as part of the precipitation percolating deep into the groundwater zone.Runoff begins when the precipitation rate exceeds the infiltration rate of the soil (Figure 4) and the storage in the depressions is full.Values for the canopy and surface storage were obtained from the analysis of Land Use and the DEM maps as derived from Tables 2 and 3, respectively.The Soil Moisture Accounting (SMA) Loss Method included in HEC-HMS was employed to model infiltration losses combined with canopy and surface methods.The canopy is a component of the sub-basin which is intended to represent the presence of plants in the area.The canopy method specified in the HEC-HMS settings is mainly used for continuous simulations.All rainfall is trapped until the canopy storage capacity is filled.Subsequently, the excess precipitation falls on the soil surface after passing through the canopy.The surface method specified in the HEC-HMS settings represents the depression of the ground surface where water accumulates after the pores in the soil are filled to the soil's field capacity.The water on the surface then accumulates in the depression storage until runoff begins to occur as part of the precipitation percolating deep into the groundwater zone.Runoff begins when the precipitation rate exceeds the infiltration rate of the soil (Figure 4) and the storage in the depressions is full.Values for the canopy and surface storage were obtained from the analysis of Land Use and the DEM maps as derived from Tables 2 and 3, respectively.The required inputs for the SMA model are presented in Table 4.The maximum infiltration rate was determined as the upper limit of the rate of water entry from surface storage into the soil [8].The  The required inputs for the SMA model are presented in Table 4.The maximum infiltration rate was determined as the upper limit of the rate of water entry from surface storage into the soil [8].The values for maximum infiltration rate were obtained based on the soil analysis in the catchment and represent the saturated hydraulic conductivity (Table 5).Impervious area was defined as the percentage of the area under urban civilization for the sub-basin using Google Earth with the aid of the Land Use map.Soil water storage was defined as the porosity in Table 5, which is the available space that water can occupy in the soil.Tension storage was determined from Soil-Plant-Air-Water (SPAW) computer software (6.02.75, United States Department of Agriculture-USDA, Washington, DC, USA) [8] by considering it as the field capacity of the soil based on the soil texture values.The soil percolation rate and the first groundwater layer (GW1) percolation rate were chosen as the average hydraulic conductivity of all sub-basins as obtained from SPAW software based on soil texture [9].Storage coefficients and depths of GW1 and GW2 were determined based on a stream flow recession analysis of historical flow data.The values of the percolation rate of GW2 were obtained during the calibration process.The SCS (Soil Conservation Service) unit hydrograph method was used to model the transformation of precipitation excess into direct surface runoff.The other important parameter is the Lag time (T LAG ), an important concept in linear modelling of basin response which was calculated using Equation (1) [11]: T LAG = L 0.8 (S + 1) 0. The Curve Number (CN) was manually estimated for the sub-basins, based on the hydrologic soil group (SHG) and the land cover type.After determining the required soil and land cover characteristics, the CN was estimated for each unit of the sub-basin, followed by area-weighting for the whole sub-basin.The tables used for computation are found in the Technical Release Number 55 [12].The linear reservoir method was employed to model base flow with the following required parameters:

•
Groundwater 1 initial (m 3 /s): initial base flow at the beginning of the simulation for the first layer of groundwater.

•
Groundwater 1 coefficient (h): the response time of the sub-basin as specified in the SMA model.

•
Groundwater 1 reservoir is used so that the base flow is routed through several sequential reservoirs.The base flow is attenuated when the number of reservoirs is increased.
The same parameters are also defined for the second layer of groundwater.The Lag routing model was used to model the reaches through the lag time, which is the time that falls between the rainfall and flow peaks.

Calibration and Validation
Before a hydrological model can be considered to have outputs that are reliable, it needs to be calibrated and validated using observed stream flow.The simulated stream flow must be compared to the observed stream flow to evaluate the goodness of fit and conclude whether the model is able to predict and present credible results.In this work, the model was calibrated using the identified parameters to achieve good fit between the simulated and observed data.The auto-calibration (through optimization trials) tool available in the HEC-HMS model was used for optimizing the estimates of the model parameters.The Nelder-Mead [14] method was used as the search method for optimization.However, the auto-calibration process in HEC-HMS does not necessarily converge to the desired optimum results, so in addition, a manual calibration was performed.The input parameters can be determined by estimating, observing and measuring the watershed and stream characteristics, except for the percolation rate of the groundwater through the second layer (GW2), which is estimated during calibration and optimized to obtain the best fit of simulated flow to observed flow.
The validation process entailed using the optimized parameters of a different period and confirming the goodness of fit for the observed and simulated stream flow.

Model Performance Evaluation
The HEC-HMS model performance evaluation involved assessing the goodness of fit in the observed and simulated stream flow using: 1.
The Percentage Error in Volume (PEV) where Vol o , Vol s are the observed and simulated volumes, respectively.2.
The percentage Error in Peak Flow (PEPF) where Q o(peak) , Q s(peak) are the observed and simulated flows, respectively.

3.
The Coefficient of correlation (R 2 ) where O i , S i are the observed and simulated flows at time i, respectively; and Ō, S are the average observed and simulated flows during the calibration period, respectively.R 2 is a standard regression criterion which is limited because it only evaluates the linear relationships between observed and simulated variables.Legates and Mccabe [15] and Moriasi et al. [16] recommended that a good model efficiency criterion have at least three important components, which are improvements over the coefficient R 2 : one dimensionless statistic, one absolute error index statistic and one graphical technique.That means that none of the statistics should be used alone; applied together, they form a set of model selection criteria which offsets each other's limitations [17,18].According to the categorization of the most used statistics, we added the following criteria:
The absolute error index represented by the Root Mean Squared Error (RMSE)-standard deviation ratio (RSR) of observations given by: RSR ranges from 0 (optimal value) to 0.5 (very good performance) ratings for both calibration and validation periods.Lower values of RSR indicate a lower root mean square error normalized by the standard deviation of the observations, which indicates the appropriateness of the model simulation [16,20].
To interpret the results, Tables 6 and 7 from [16] were used as a guide.

Sensitivity Analysis
Sensitivity analysis is a critical component of rainfall-runoff modeling that helps to identify influential parameters.Each of the optimized parameters must be used as inputs into the model to get the desired simulated results.Some of these parameters are more sensitive than others.That is, a minor change in the value can lead to a big difference in the simulated flow.Thus, the most sensitive parameters of the model need to be precisely estimated.
To do so, the model was first run with the optimized model parameters obtained after calibration and validation.Next, one parameter at a time method was applied: the value of each parameter was varied from −30% to +30% in increments of 10%, keeping all other parameters constant.The output values (simulated runoff data, volume, peaks, seven-day low flow (7-MALF) and baseflow index (BFI)) were analyzed to determine variation with respect to the initial estimates of the parameters.The elasticity ratio (e) [21] was used to rank the parameters in descending order from most to the least sensitive.Also called the relative sensitivity, e expresses the relative change in the dependent variable with respect to the independent variable.The elasticity ratio is invariant to the dimensions of the variables and is given by Equation ( 9) [22,23].A greater elasticity ratio indicates a more highly sensitive variable.

e =
∆O/O ∆I/I = %change output %change input (9) where O, I are the output and the input variables, respectively.

Calibration and Validation
The manual calibration of thirteen HEC-HMS parameters was conducted using observed stream flow data from the 3KD06 gauging station for the period 1988 to 1991.The calibration was started using the basic values that were also used for sensitivity analysis.These values were manually modified until a good fit between the simulated and observed stream flows was obtained.The goodness of fit was evaluated using hydrograph visualization and computed statistics values.The observed and simulated stream flow comparison graphs are shown in Figure 5

Calibration and Validation
The manual calibration of thirteen HEC-HMS parameters was conducted using observed stream flow data from the 3KD06 gauging station for the period 1988 to 1991.The calibration was started using the basic values that were also used for sensitivity analysis.These values were manually modified until a good fit between the simulated and observed stream flows was obtained.The goodness of fit was evaluated using hydrograph visualization and computed statistics values.The observed and simulated stream flow comparison graphs are shown in Figure 5 for the calibration period (1988 to 1991).The comparison shows a close agreement between simulated stream flow and observed stream flow in terms of peak values and an acceptable agreement in terms of stream flow distribution.The parameters were optimized to obtain a closer agreement between the simulated and observed stream flow.Table 8 shows the optimized values of the calibration parameters for the study area.Model validation involved running the model using the same input parameters used for the calibration process.Based on the optimized parameters, the model was run for the validation period (1992 to 1995) to check the capability of the model to predict runoff at the 3KD06 gaging station.The observed and simulated stream flow comparison graphs are shown in Figure 6 for the validation period (1992 to 1995).The comparison shows an acceptable agreement between the simulated stream flow and the observed stream flow in terms of peak value, as well as the stream flow distribution.Table 8 shows the optimized values of the calibration parameters for the study area.Model validation involved running the model using the same input parameters used for the calibration process.Based on the optimized parameters, the model was run for the validation period (1992 to 1995) to check the capability of the model to predict runoff at the 3KD06 gaging station.The observed and simulated stream flow comparison graphs are shown in Figure 6 for the validation period (1992 to 1995).The comparison shows an acceptable agreement between the simulated stream flow and the observed stream flow in terms of peak value, as well as the stream flow distribution.

Model Performance Evaluation
Continuous model performance evaluation was conducted for each year and for the calibration and validation periods.The time series of simulated and observed flows from the results of the simulation run in the HEC-HMS model were analyzed in Microsoft Excel to compute the statistics used for performance evaluation.The statistics that were used are described in Section 2.5, presented earlier.The performance ratings of these parameters are presented in Tables 6 and 7.During calibration for the 3KD06 gauging station, the values of PEV, PEPF, R 2 , d, NSE and RSR were found to be 8.2%, 10.8%, 0.80, 0.94, 0.80 and 0.46, respectively.Similarly, during the validation, the model evaluation criteria for PEV, PEPF, R 2 , d, NSE and RSR were found to be 21.7%, 46.1%, 0.67, 0.88, 0.65 and 0.62, respectively.Table 9 displays the performance evaluation of the SMA model specified in the HEC-HMS settings for individual years and for the calibration and validation periods.According to Tables 6 and 7, the model performance ranges from satisfactory to very good, except for the year 1995, which may be due to various uncertainties in the observed stream flow data.According to the results, and based on the values of NSE, PEPF, PEV and R 2 calculated, the model performance is satisfactory.Furthermore, the index of agreement (d) ranges from 0.81 to 0.96, and the RSR from 0.37 to 0.74, which are also indicators of good model performance.

Model Performance Evaluation
Continuous model performance evaluation was conducted for each year and for the calibration and validation periods.The time series of simulated and observed flows from the results of the simulation run in the HEC-HMS model were analyzed in Microsoft Excel to compute the statistics used for performance evaluation.The statistics that were used are described in Section 2.5, presented earlier.The performance ratings of these parameters are presented in Tables 6 and 7.During calibration for the 3KD06 gauging station, the values of PEV, PEPF, R 2 , d, NSE and RSR were found to be 8.2%, 10.8%, 0.80, 0.94, 0.80 and 0.46, respectively.Similarly, during the validation, the model evaluation criteria for PEV, PEPF, R 2 , d, NSE and RSR were found to be 21.7%, 46.1%, 0.67, 0.88, 0.65 and 0.62, respectively.Table 9 displays the performance evaluation of the SMA model specified in the HEC-HMS settings for individual years and for the calibration and validation periods.According to Tables 6 and 7, the model performance ranges from satisfactory to very good, except for the year 1995, which may be due to various uncertainties in the observed stream flow data.According to the results, and based on the values of NSE, PEPF, PEV and R 2 calculated, the model performance is satisfactory.Furthermore, the index of agreement (d) ranges from 0.81 to 0.96, and the RSR from 0.37 to 0.74, which are also indicators of good model performance.

Sensitivity Analysis
The sensitivities of thirteen SMA parameters (except the five initial conditions for the five storage layers in the SMA model; see Table 4) were analyzed by separately varying each parameter from −30% to +30% in increments of 10%.A one-parameter-at-a-time method was used for this sensitivity analysis; changing one parameter while keeping the others constant.The percentage of variation in simulated volume, peak, seven day low flow (7-MALF), baseflow index (BFI) and Nash-Sutcliffe Efficiency (NSE) were plotted against the percentage of variation of each parameter, as shown in Figures 7-11  Legend Very Good Good Satisfactory Unsatisfactory

Sensitivity Analysis
The sensitivities of thirteen SMA parameters (except the five initial conditions for the five storage layers in the SMA model; see Table 4) were analyzed by separately varying each parameter from −30% to +30% in increments of 10%.A one-parameter-at-a-time method was used for this sensitivity analysis; changing one parameter while keeping the others constant.The percentage of variation in simulated volume, peak, seven day low flow (7-MALF), baseflow index (BFI) and Nash-Sutcliffe Efficiency (NSE) were plotted against the percentage of variation of each parameter, as shown in Figures 7-11           After analyzing the elasticity of the different parameters, they were ranked from the most to the least sensitive.The sensitivity analysis was done to determine the sensitivity of the computed runoff volume, the computed peak, the Nash-Sutcliffe, the seven day low flow (7-MALF) and the baseflow index (BFI) to the SMA model parameters.The runoff volume was found to be more sensitive to the GW1 percolation rate, the GW1 storage coefficient and the GW1 water storage (Table 10), while the peak volume was found to be more sensitive to the impervious area, the GW1 water storage, the GW1 percolation rate and the GW1 storage coefficient (Table 11).At the same time, the Nash-Sutcliffe Efficiency was found to be more sensitive to GW1 water storage, GW1 coefficient storage and GW1 percolation rate respectively (Table 12).The seven-day low flow was found to be more sensitive to the GW2 storage coefficient, the GW2 water storage and the GW2 percolation rate (Table 13).Finally, the baseflow index was found to be more sensitive to the impervious area, the GW2 percolation rate and the GW2 storage coefficient (Table 14).After analyzing the elasticity of the different parameters, they were ranked from the most to the least sensitive.The sensitivity analysis was done to determine the sensitivity of the computed runoff volume, the computed peak, the Nash-Sutcliffe, the seven day low flow (7-MALF) and the baseflow index (BFI) to the SMA model parameters.The runoff volume was found to be more sensitive to the GW1 percolation rate, the GW1 storage coefficient and the GW1 water storage (Table 10), while the peak volume was found to be more sensitive to the impervious area, the GW1 water storage, the GW1 percolation rate and the GW1 storage coefficient (Table 11).At the same time, the Nash-Sutcliffe Efficiency was found to be more sensitive to GW1 water storage, GW1 coefficient storage and GW1

Discussion
Before using the HEC-HMS model to get an accurate prediction of runoff in the Mkurumudzi catchment, the model needed to be well calibrated using the SMA parameters.Most of these parameters are related in nature to the soil properties and require careful observation and some field investigations to reach a given level of accuracy.However, in the Mkurumudzi catchment, such records are lacking, and no investigations were conducted for the present study; instead, the parameters were obtained from the literature and other secondary sources.With this kind of data estimation, the results obtained are highly satisfactory.Evapotranspiration is one of the most important factors for continuous modeling and was estimated using the Penman-Monteith equation [4].
The percentage error in volume (PEV) for the continuous model ranges from 1.8% to 49.4%, with satisfactory global model performance, except in 1995.The percentage error in peaks ranges from 10.8% to 46.3%, indicating very good and good model performance, respectively, except for 1994.The Nash-Sutcliffe efficiencies of the model range from 0.45 (1995) to 0.80, indicating satisfactory to very good model performances, except for the year 1995.The coefficient of correlation R 2 ranges from 0.52 to 0.87, which, according to Table 6, indicates a satisfactory to very good performance.These results are similar to the NSE and the RSR results.As noted, the years 1994 and 1995 exhibit very low model efficiencies and low correlation between observation and simulation; these differences can be attributed to uncertainties in the stream flow dataset.Probably, the use of data from multiple gauging stations would have improved the model performance.
During the sensitivity analysis of the soil moisture accounting (SMA) parameters, the groundwater layer parameters and the impervious area were found to be some of the most sensitive parameters for runoff simulation.Fleming and Neary [24] performed a similar sensitivity analysis of a continuous HEC-HMS model for the Dale Hollow basin in Kentucky and Tennessee.They found the maximum infiltration rate, the maximum soil depth, and the tension zone depth to be the most sensitive parameters.However, Ouedraogo et al. [25], in 2017, also used HEC-HMS to perform runoff simulation in Ruiru reservoir catchment.They found soil storage to be the most sensitive parameter, followed by the groundwater storage coefficient and the soil tension storage capacity.These findings are different from the current study results.Singh and Jain [9] also conducted continuous hydrological modeling in Vamsadhara River Basin (India) using the SMA model and found the soil storage to be the most sensitive parameter.

Figure 1 .
Figure 1.Land use in the Mkurumudzi basin.

Figure 1 .
Figure 1.Land use in the Mkurumudzi basin.
3 from ESRI.The meteorological component contains the observed precipitation and discharge data; the control specification component determines the simulation period and the time step to be used.The observed historical data of the Kidongo Park Gate rain station representing each sub-basin and one stream gauge station (3KD06) were used for model calibration and validation.Due to the location of the 3KD06 station, only the two (02) upstream sub-basins M1 and M2 (Figure 3) were used.A daily time step was used for the simulation, given the time interval of the available observed data.
lag time (h).L = hydraulic length of the watershed (ft.).Y = watershed slope (%).S = maximum retention in the watershed (mm) as defined by: S = 25400/CN − 254 (2) CN = SCS curve number for the watershed for the calibration period (1988 to 1991).The comparison shows a close agreement between simulated stream flow and observed stream flow in terms of peak values and an acceptable agreement in terms of stream flow distribution.The parameters were optimized to obtain a closer agreement between the simulated and observed stream flow.Hydrology 2018, 5, x FOR PEER REVIEW 11 of 19

Figure 7 .
Figure 7. Percentage changes in simulated volume plotted against the percentage variation of each parameter.Figure 7. Percentage changes in simulated volume plotted against the percentage variation of each parameter.

Figure 7 .
Figure 7. Percentage changes in simulated volume plotted against the percentage variation of each parameter.Figure 7. Percentage changes in simulated volume plotted against the percentage variation of each parameter.

Figure 7 .
Figure 7. Percentage changes in simulated volume plotted against the percentage variation of each parameter.

Figure 8 .
Figure 8. Percentage changes in simulated peak plotted against the percentage variation of each parameter.Figure 8. Percentage changes in simulated peak plotted against the percentage variation of each parameter.

Figure 8 . 19 Figure 9 .
Figure 8. Percentage changes in simulated peak plotted against the percentage variation of each parameter.Figure 8. Percentage changes in simulated peak plotted against the percentage variation of each parameter.Hydrology 2018, 5, x FOR PEER REVIEW 14 of 19

Figure 9 .
Figure 9. Percentage changes in simulated NSE plotted against the percentage variation of each parameter.Figure 9. Percentage changes in simulated NSE plotted against the percentage variation of each parameter.

Figure 9 .
Figure 9. Percentage changes in simulated NSE plotted against the percentage variation of each parameter.

Figure 10 .
Figure 10.Percentage changes in simulated 7-MALF plotted against the percentage variation of each parameter.Figure 10.Percentage changes in simulated 7-MALF plotted against the percentage variation of each parameter.

Figure 10 .
Figure 10.Percentage changes in simulated 7-MALF plotted against the percentage variation of each parameter.Figure 10.Percentage changes in simulated 7-MALF plotted against the percentage variation of each parameter.Hydrology 2018, 5, x FOR PEER REVIEW 15 of 19

Figure 11 .
Figure 11.Percentage changes in simulated BFI plotted against the percentage variation of each parameter.

Figure 11 .
Figure 11.Percentage changes in simulated BFI plotted against the percentage variation of each parameter.

Table 1 .
Summary of annual hydrological data of the Mkurumudzi river basin and gauging station 3KD06.

Table 9 .
Performance evaluation of the continuous HEC-HMS model.

Table 10 .
SMA parameters sensitivity ranking for runoff volume.

Table 14 .
SMA parameters sensitivity ranking for BFI.