Applying SHETRAN in a Tropical West African Catchment ( Dano , Burkina Faso ) — Calibration , Validation , Uncertainty Assessment

This study presents the calibration and validation of the physically based spatially distributed hydrological and soil erosion model SHETRAN for the Dano catchment, Burkina Faso. A sensitivity analysis of six model parameters was performed to assess the model response and to reduce the number of parameters for calibration. The hydrological component was calibrated and validated using observed discharge data of two years. Statistical quality measures (R2, NSE, KGE) ranged from 0.79 to 0.66 during calibration and validation. The calibrated hydrological component was used to feed the erosion modeling. The simulated suspended sediment load (SSL) was compared with turbidity-based measurements of SSL of two years. Achieved quality measures are comparable to other SHETRAN studies. Uncertainties of measured discharge and suspended sediment concentration were determined to assess the propagated uncertainty of SSL. The comparison of measurement uncertainties of discharge and SSL with parameter uncertainty of the corresponding model output showed that simulated discharge and SSL were frequently outside the large measured uncertainty bands. A modified NSE was used to incorporate measurement and parameter uncertainty into the efficiency evaluation of the model. The analyses of simulated erosion sources and spatial patterns showed the importance of river erosion contributing more than 60% to the total simulated sediment loss.


Introduction
Soil degradation by water-related soil erosion is a major environmental problem threatening food security, income, and environmental health especially in tropical and subtropical countries [1][2][3][4].On the one hand, systemic and natural reasons such as unfavorable climatic conditions and the structural instability of soils resulting from low soil organic carbon (SOC) content are responsible for a high erosion risk in these regions [1].On the other hand, socio-economic factors contribute to the problem and are often responsible for the severity of soil erosion.Among those factors, increased pressure on land resources through population growth is highlighted as a major reason for accelerated soil erosion [4][5][6][7].
Soil erosion strongly varies in space and time.Thus measuring soil erosion requires a large personnel and financial effort and despite advances in measurement technology it is often impossible to perform measurements over the required spatial and temporal scales.This especially applies to data scarce regions such as West Africa.To overcome these drawbacks soil erosion models have been frequently used [8,9].Erosion models have also been implemented to predict the effect of land use and climate change on soil erosion and to identify areas where erosion control measures are necessary [9].In general, three types of erosion models are differentiated: Empirical erosion models, conceptual models, and physically based erosion models.Physically based erosion models are based on physical principles such as the conservation of mass and momentum [10].After an evaluation of the available erosion models the physically based spatially distributed soil erosion model SHETRAN [11] was chosen to simulate hydrological and soil erosion processes in a tropical West African catchment.The two main reasons for using SHETRAN in this study are its ability to simulate the dominant erosion processes [12,13] and the continuity in simulation necessary for the prediction of land use and climate change.The study considers soil erosion by surface runoff.
Physically based models, such as SHETRAN, need to be calibrated and validated.Although parameters of physically based models theoretically do not need calibration, adjustments are necessary to account for an unrealistic representation of environmental properties such as grid size or channel geometry.In the present context, calibration is therefore considered as the adjustment of parameter values to overcome the unrealistic representation of environmental properties.The calibration process requires knowledge on the sensitivity of results to model parameters (i) to better assess the model response; (ii) to reduce the number of parameters for calibration; and (iii) to define the parameter uncertainty [14,15].Outputs of environmental models are subject to uncertainty that is related to the simplification of the model, the uncertainty of the model parameter, and the measurement uncertainty [8,16,17].Measurement uncertainty refers to the uncertainty of measured data used as input (e.g., precipitation) and to calibrate and validate the model.Studying the uncertainty associated with the modeling output is important since decision making is often based on the output of environmental models [18].The present study focuses on the assessment of two sources of uncertainty in environmental modeling, the parameter and the measurement uncertainty.Information on and measurements of the required model parameters are difficult if possible at all and often not available for the study area or on the required spatial or temporal scale.Quantifying parameter uncertainty is necessary for the interpretation of model outcomes and its application in environmental planning [9,15].
Measured variables, such as water discharge and suspended sediment concentration (SSC), are subject to uncertainties.These uncertainties are propagated if different measurements are combined to calculate variables that are used to calibrate and validate the model [18][19][20][21].For instance, the suspended sediment load (SSL), which is frequently used to validate erosion models, is calculated using the SSC and water discharge.Defining ranges of measured uncertainties is important to contextualize the simulated output.
Studies on soil erosion in West Africa often use USLE-based modelling approaches since data required to run and validate complex, physically based and spatially distributed hydrological soil erosion models are rarely available (see Table 1).Empirically based erosion models are relatively easy to use and only a few input data sets are required to run these models.However, their process representation and applicability to complex conditions of land use and climate change is limited [9].Most of the studies listed in Table 1 are not comparable with the present study since they use different modeling approaches regarding catchment size, continuity, and considered output.Despite numerous studies that applied SHETRAN in different regions and with multiple objectives (see Table 2) the model has not been tested in the West African environment.The environmental properties of the study catchment such as the low slope angles and the peculiar rainfall pattern may provide a new challenge to SHETRAN.Therefore, the present study may also serve as a model test of the particular environmental conditions found in the study region.Furthermore, the present study can also be considered as an independent check on the ease of use of the model as most of the studies that used SHETRAN before were conducted by members of the team that developed the model.Thus, this study aims to: (1) assess the uncertainty of measured discharge and SSL used to calibrate and validate the hydrological and erosion components of SHETRAN; (2) perform a detailed sensitivity analysis to define parameter ranges and to reduce the number of calibration parameters; (3) use a Latin Hypercube Sampling approach to calibrate the model and to define uncertainty bounds of simulated discharge and SSL; (4) evaluate model performance considering the uncertainty of measured data used to compare the model output and parameter uncertainty.

Study Area
The investigated Dano catchment in the Ioba province covers an area of 126 km 2 and is located in the Southwest of Burkina Faso (Figure 1).The study area is in a focal watershed of the WASCAL program (West African Science Service Center on Climate Change and Adapted Land Use).The multidisciplinary program is designed to study the influence of climate and land use/land cover change on human and environmental systems and to enhance their resilience.
Agricultural land use is the most important land use category in the region (Figure 1d).The agricultural area has expanded in recent decades due to a growing demographic pressure indicated by an annual population growth of 3%.It has gradually been intensified accompanied by reduced fallow periods and expansion to marginal land areas with adverse effects on soil fertility [4,[46][47][48][49].Since 1990 each year on average 2% of the savanna in the study area was converted to agricultural land [49].The general appearance of the vegetation in the Sudano savannah is dominated by open forests and wide arborous and shrubby areas.The main staple food crops cultivated in the region are sorghum (Sorghum bicolor), millet (Pennisetum glaucum), maize (Zea mays), cowpeas (Vigna unguiculata), and groundnut (Arachidis hypogaea).Cotton (Gossypium hirsutum) is the most important cash crop.During the rainy season between 40% and 70% of the soil is covered by vegetation.
The catchment is dominated by a flat and slightly undulating landscape characterized by low slope gradients (average and maximum gradients are 3.1% and 38%, respectively, Figure 1b) and an elevation ranging from 236 to 565 m above sea level (m a.s.l).The annual mean temperature is 28.6 • C and annual precipitation ranges from 800 to 1200 mm/a for the period 1951-2005 [50].The rainfall pattern is uni-modal and characterized by a distinct rainy season from May to October and a dry season from November to April.Eighty percent of the rain falls between July and September with high rainfall intensities.As an example from the Dano catchment, 60 mm/h were measured as maximum in 2014.The flow regime is ephemeral and the channel geometry is divers ranging from strongly incised (3-4 m) clearly defined channels to broader inland valleys.Information on the ranges of measured parameters is given in Table 3.
Most of the soils (73%) are plinthosols according to the World Reference Base for soil resources (WRB) [51] characterized by a high content of coarse particles and a plinthitic subsurface layer in the first meter of the profile.Other soils that were formed in the region are gleysols, cambisols, lixisols, leptosols, and stagnosols (Figure 1e).

Data Sources
In order to calibrate and validate SHETRAN, multiple data sets are required (Table 4).Existing data were complemented by a measurement network consisting of five automatic climate stations and pluviometers (tipping bucket type) and one discharge and turbidity station.Additionally, physical and chemical analyses of soil samples were done to retrieve necessary soil parameters for the erosion model (Table 4).Notes: 1 α and n are van Genuchten empirical parameters; 2 Ksat refers to the saturated hydraulic conductivity; 3 θsat to the saturated water content; 4 θres to the residual water content; 5 LAI to the leaf area index; and 6 ETp/ETa ratio to the ratio of potential evapotranspiration to actual evapotranspiration.

Model Description
Modeling of hydrological and erosion processes was performed using the physically based, spatially distributed and raster-based model SHETRAN [11,12].SHETRAN is based on SHE (Système Hydrologique Europeen) which was jointly developed by the British Institute of Hydrology, the Danish Hydraulic Institute and the French consulting company SOGREAH [54].

Model Description
Modeling of hydrological and erosion processes was performed using the physically based, spatially distributed and raster-based model SHETRAN [11,12].SHETRAN is based on SHE (Système Hydrologique Europeen) which was jointly developed by the British Institute of Hydrology, the Danish Hydraulic Institute and the French consulting company SOGREAH [54].During the last thirty years SHETRAN has been continuously improved and equipped with new components that include e.g., the sediment component [12,55] and a fully 3D subsurface water flow component [56].A summary of SHETRAN applications with various objectives and in different regions is given in Table 2.
Detailed information about the model is given in Bathurst [57].A short overview of the most important hydrological process descriptions of the model is summarized in the following list:

•
Fully 3D subsurface flow simulation based on Richards' equation.

•
Infiltration is calculated using Richards' equation.

•
Overland and channel flow is calculated using the diffusive wave approximations of the full Saint-Venant equation.

•
Potential evapotranspiration (ETp): Potential plant transpiration, evaporation from intercepting surfaces and from bare soil as well as water bodies was calculated externally based on the Penman-Monteith equation [58] and added as input into SHETRAN.

•
Actual evapotranspiration (ETa) is estimated based on the approach introduced by Feddes et al. [59] where the ratio ETa/ETp is a function of soil moisture tension.The ratio ETa/ETp at field capacity is the input parameter and the reduction of ETa with decreasing soil moisture tension is calculated based on this parameter.

•
Interception is calculated based on the approach by Rutter et al. [60,61] who relates interception to the leaf area index, the vegetation cover, and the maximum depth of water on leaves.
The parameterization and calibration of land use and soil properties was done based on data obtained from literature and measurements (see Tables 4 and 5).The parameters (θ sat , θ res , α, n) used to describe the soil water retention curve after van Genuchten [62] were determined from soil texture and organic matter content following Rawls and Brakensiek [63].Measured saturated hydraulic conductivity (K sat ) was used for the top soil horizon.For the remaining horizons K sat was calculated using soil texture and organic matter content following Brakensiek and Rawls [64].
SHETRAN requires different types of input data.Spatially distributed data, including digital elevation model (DEM), the soil and land use map, were used in a raster format with a grid resolution of 200 m × 200 m.The applied resolution is relatively coarse compared with other applications of SHETRAN with resolutions typically below 100 m (Table 2).Nevertheless, the topography of the study area is characterized by long straight slopes which are well represented in this resolution.Zhang [33] applied a resolution of 2 km to a larger catchment (705 km 2 ) and compared it with resolutions of 0.5 and 1 km.The performance measure using the Nash-Sutcliff-Efficiency (NSE) decreased by 3.7% with decreasing resolution (from 1 to 2 km) as a result of information loss as land use and soil type maps become coarser.
Precipitation and potential evapotranspiration (ETp) are given as time series over two years for each of the five stations considered in the modeled catchment.The area that is represented by each station is determined by Thiessen polygons.A pre-processing software uses the DEM to determine the river geometry and produce the input files [65].The temporal resolution of 1 hour used here is the standard timestep of SHETRAN and commonly used in other studies (see Table 2).The precipitation input has an hourly timestep.
A short summary of erosion processes simulated by SHETRAN is given below.Soil detachment is accounted for by three separate equations describing detachment by raindrop/leaf drip (Equation ( 1)) [66], by overland flow (Equation ( 2)) [67] and by channel flow (Equation ( 3)) [68]: where D r is the rate of soil detachment (kg/m 2 /s), F w (-) accounts for the protection against drop detachment by surface water, k r is the raindrop impact erodibility coefficient (J −1 ), C g is the proportion of ground covered by near ground vegetation (%), C r is the rock cover (-), M r /M d is the momentum squared of raindrops/leaf drips reaching the ground per unit time and area (kg 2 /s 3 ), where D q is the rate of soil detachment per unit area (kg/m 2 /s), k f is the overland flow erodibility coefficient (kg/m 2 /s), C r is the proportion of ground shielded by rock cover (-), τ is the shear stress exerted by overland flow (N/m 2 ), τ ec is the critical shear stress for the initiation of motion (N/m 2 ), where E b is the detachment rate of bank material per unit area (kg/m 2 /s), BKB is the bank erodibility coefficient (kg/m 2 /s), τ bc is the critical shear stress for the initiation of motion of bank material (N/m 2 ) and τ b is the shear stress acting on the bank (N/m 2 ).
Sediment is transported based on the transport capacity of overland (Equation ( 4)) [69] and channel flow (Equation ( 5)) [70]: where G tot is the transport capacity rate for overland flow (m 3 /s), τ is the shear stress (N/m 2 ), p is the water density (kg/m 3 ), l is the width of flow (m), Q is the water discharge (m 3 /s), D 50 is the median sediment diameter, δ and a are parameters, where G i is the transport capacity rate of particle size in group i (m 3 /s), D i is the particle diameter in size group i (m), H is the water flow depth (m), U is the mean water flow velocity (m/s), u * is the shear velocity (m/s), G gr,i is the dimensionless sediment transport rate for sediment size group i.Further details are given in Morgan and Nearing [71] and Wicks [55].

Model Sensitivity, Calibration, and Validation
Several parameters of SHETRAN need to be calibrated by comparing simulated and observed variables.Prior to the calibration of parameters to which the model output is most sensitive the corresponding initial values were identified based on previous studies that used SHETRAN and based on sensitivity analyses.The sensitivity analyses were done based on the "one factor at a time" (OFAT) method using Equation (6) [72]: where SI 90 is the sensitivity index, O 90 and O -90 the model output resulting from a parameter value increased or decreased by 90%, and O 0 the model output from the base run.A list of parameters to which the model output responds sensitively and the corresponding calibration ranges used in this study are given in Table 5.The parameter range of k f is quite low compared with that indicated in literature.However, as this parameter is considered to be a calibration parameter [55], we assume that the range is representative for the soil properties in the study area.
Latin Hypercube Sampling (LHS) [73] was used to generate 300 parameter sets within the defined value ranges.This is considered as a reasonable compromise between the necessary model executions which are dependent on the number of parameters used and the run time.The hydrological component of SHETRAN was calibrated based on the observed hydrograph in 2014.The soil erosion component was calibrated based on the observed SSL in 2014.The model performance was statistically evaluated by the coefficient of determination (R 2 ), the Nash-Sutcliff efficiency (NSE) [74], and the Kling-Gupta efficiency (KGE) [75,76].The model was validated using data from the year 2015.

Measurement Uncertainty
A power regression model was used to describe the relation between measured water level and water discharge.The relation between measured suspended sediment concentration (SSC) and turbidity was defined by a linear regression model.Polynomial and power regression equations of confidence intervals calculated by Equation ( 7) [82] were used to express the measured uncertainty of both regression models.
where ∆X i is the confidence interval of the predictor, Y i is the response variable, t α is the t-value at α = 0.05 significance level, SE is the standard error, DF are the degrees of freedom, X m is the mean of X, and SS xx is the sum of squared differences.
The uncertainty of water discharge and the combined uncertainty of SSL were compared with the simulated discharge and SSL through visual inspection and the calculation of R-and P-factors.The P-factor gives the proportion of the variable in question which is within the corresponding uncertainty bands in percent.The R-factor is defined as the mean width of the uncertainty band divided by the standard deviation of the variable in question [83].
Uncertainties related to field measurements of discharge and SSC sampling as well as laboratory work were not explicitly accounted for in this study.

Parameter Uncertainty
The ten best simulations were chosen out of the 300 parameter sets based on the sum of NSE, R 2 , and KGE calculated separately for discharge and SSL.The ten best hydrological simulations were combined with the ten best sediment simulations to define the parameter uncertainty bounds.Among these, the parameter set with the highest sum of performance measures and reasonable parameter values is considered to give the best representation of measured discharge and SSL.

Uncertainty Based Modification of the NSE
Based on the work by Harmel and Smith [18] and Harmel et al. [84] a modified error term e i (O i − P i ) was incorporated into the traditional calculation of the NSE that considers the uncertainty of both measured and parameter uncertainty at each observation.A correction factor (CF i ) is calculated based on the degree of overlap of the assumed distributions of the observed and predicted values (see Equation (8)).In the present study measured and predicted values are assumed to be normally distributed.
CF i ranges from 0 (total overlap) to 1 (no overlap): 8)   where prob(o i < P imax ) and prob(o i < P imin ) are the probability distributions of the observed value o i limited by the maximum (P imax ) and minimum (P imin ) predicted value and prob(p i < O imax ) and prob(p i < O imin ) are the probability distributions of the predicted value p i limited by the maximum (O imax ) and minimum (O imin ) observed 95%-uncertainty bounds.
The error term e i (O i − P i ) is then multiplied by CF i and substituted into the NSE to get the modified NSE m (Equation ( 9)).

Measurement Uncertainty
Figure 2a,b show the stage-discharge regression line and the corresponding uncertainty ranges of the years 2014 and 2015.Two different rating curves were used in order to consider the changes of channel morphology and the human intervention in the channel system.However, data collection in the focus area is challenging and the limited number of observations is also reflected in the uncertainty bands.Uncertainties increase with increasing water level and discharge due to the chosen power regression equation and the sample properties.Overbank flow was observed following intense rainfall events but could not be measured due to inaccessibility.Therefore, extrapolation beyond the measured range was done despite the increasing uncertainties during peak flows.
Figure 2c shows the linear regression line and the uncertainty band of suspended sediment concentration (SSC) and turbidity.The uncertainty ranges are almost parallel due to the linear regression equation.As a result from the error propagation the combined uncertainty of SSL is quite large.

Measurement Uncertainty
Figure 2a,b show the stage-discharge regression line and the corresponding uncertainty ranges of the years 2014 and 2015.Two different rating curves were used in order to consider the changes of channel morphology and the human intervention in the channel system.However, data collection in the focus area is challenging and the limited number of observations is also reflected in the uncertainty bands.Uncertainties increase with increasing water level and discharge due to the chosen power regression equation and the sample properties.Overbank flow was observed following intense rainfall events but could not be measured due to inaccessibility.Therefore, extrapolation beyond the measured range was done despite the increasing uncertainties during peak flows.
Figure 2c shows the linear regression line and the uncertainty band of suspended sediment concentration (SSC) and turbidity.The uncertainty ranges are almost parallel due to the linear regression equation.As a result from the error propagation the combined uncertainty of SSL is quite large.

Model Sensitivity
Figure 3 shows the model sensitivity to the investigated soil, land use and erosion parameters (Table 5) following the OFAT method.Total and maximum water discharge as well as suspended sediment load is used for comparison.
Changing parameter values of the ratio Eta/ETp (Figure 3a) specified for each land use strongly influence total runoff which is indicated by a relatively high SI90 of 1.8.As this parameter directly affects the actual evapotranspiration (ETa), an increase of e.g., 90% leads to 60% less surface runoff due to higher ETa.

Model Sensitivity
Figure 3 shows the model sensitivity to the investigated soil, land use and erosion parameters (Table 5) following the OFAT method.Total and maximum water discharge as well as suspended sediment load is used for comparison.
Changing parameter values of the ratio Eta/ETp (Figure 3a) specified for each land use strongly influence total runoff which is indicated by a relatively high SI 90 of 1.8.As this parameter directly affects the actual evapotranspiration (ETa), an increase of e.g., 90% leads to 60% less surface runoff due to higher ETa.
The roughness coefficient KSTR (Figure 3b) is given for each land use type and controls the surface roughness.Larger KSTR results in faster surface runoff and therefore especially influences the maximum runoff.However, interactions between surface roughness, infiltration, and evapotranspiration also lead to a change of total catchment runoff.The total runoff responds less sensitively (SI 90 = 0.4).An increase of 90% results in 6% higher total runoff but increases the maximum runoff by 52%.The higher sensitivity of the maximum discharge is also shown by SI 90max of 1.2.
The roughness coefficient KSTR (Figure 3b) is given for each land use type and controls the surface roughness.Larger KSTR results in faster surface runoff and therefore especially influences the maximum runoff.However, interactions between surface roughness, infiltration, and evapotranspiration also lead to a change of total catchment runoff.The total runoff responds less sensitively (SI90 = 0.4).An increase of 90% results in 6% higher total runoff but increases the maximum runoff by 52%.The higher sensitivity of the maximum discharge is also shown by SI90max of 1.2. Figure 3c,d show the effect of the changing erodibility coefficients kf and kr on the sediment yield respectively.An increase leads to higher erosion in both cases.However, the results are more sensitive to kr in comparison with kf as indicated by the higher SI90.An increase of kr by 90% leads to 10% higher total sediment yield while an increase of Kf by 90% results in only 2.8% higher sediment yield. Figure 3c,d show the effect of the changing erodibility coefficients k f and k r on the sediment yield respectively.An increase leads to higher erosion in both cases.However, the results are more sensitive to k r in comparison with k f as indicated by the higher SI 90 .An increase of k r by 90% leads to 10% higher total sediment yield while an increase of K f by 90% results in only 2.8% higher sediment yield.
Figure 3e shows the model's response to a changing river bank erodibility coefficient (BKB).The SI 90 is the highest in comparison with the other parameter indicating the importance of bank erosion for the catchment sediment yield.An increase of 90% leads to a 40% higher sediment yield.
DLSMAX can be considered as maximum sediment storage depth above which the soil is protected against erosion [33].In other words, the soil material that cannot be transported due to an insufficient transport capacity is considered to be available and stored as loose sediment till it reaches DLSMAX.Figure 3f shows that a 90% increase of DLSMAX leads to a 14.5% higher SSL since more soil material can be stored as available sediment.

Hydrological Modelling
Two parameters (ETa/ETp, KSTR) of SHETRAN to which the results are sensitive were used to calibrate the hydrological component.Although other parameters such as K sat or soil hydrological parameters are reported to be important, a choice was made based on previous studies [34,39,40,85,86] in order to limit run time.Based on the sum of R 2 , KGE, and NSE, several parameter sets gave satisfactory to good quality measures according to the equifinality concept introduced by Beven and Freer [16].The hydrograph of the simulation having the highest sum of the performance indices and a reasonable parameter setting is shown in Figure 4.For the calibration period the parameter uncertainty is based on the ten best parameter sets.
Figure 3e shows the model's response to a changing river bank erodibility coefficient (BKB).The SI90 is the highest in comparison with the other parameter indicating the importance of bank erosion for the catchment sediment yield.An increase of 90% leads to a 40% higher sediment yield.
DLSMAX can be considered as maximum sediment storage depth above which the soil is protected against erosion [33].In other words, the soil material that cannot be transported due to an insufficient transport capacity is considered to be available and stored as loose sediment till it reaches DLSMAX.Figure 3f shows that a 90% increase of DLSMAX leads to a 14.5% higher SSL since more soil material can be stored as available sediment.

Hydrological Modelling
Two parameters (ETa/ETp, KSTR) of SHETRAN to which the results are sensitive were used to calibrate the hydrological component.Although other parameters such as Ksat or soil hydrological parameters are reported to be important, a choice was made based on previous studies [34,39,40,85,86] in order to limit run time.Based on the sum of R 2 , KGE, and NSE, several parameter sets gave satisfactory to good quality measures according to the equifinality concept introduced by Beven and Freer [16].The hydrograph of the simulation having the highest sum of the performance indices and a reasonable parameter setting is shown in Figure 4.For the calibration period the parameter uncertainty is based on the ten best parameter sets.The calculated NSE for the best hydrological simulation is 0.7 and 0.66, the KGE 0.79 and 0.76 and the R 2 0.72 and 0.7 for calibration (2014) and validation ( 2015) respectively.The model performance is good and in the range of other studies that used SHETRAN (see Table 2).Among these studies R 2 and NSE values above 0.5 are frequently reported.
Measurement uncertainty, as presented by the 95%-confidence interval, is large for discharge especially during peak flows.The maximum uncertainty ranges from 17.3 to 40.3 m 3 /s.However, P-factors show that the model is often not able to simulate discharge within the measured uncertainty bounds even if the bounds are wide (Figure 4).This is supported by the small difference between NSE and NSE m .A higher NSE m would signify a greater overlap between the simulated and observed distributions.However, overlapping areas are observed only during peak flows and that does not change the NSE m substantially.Figure 4c shows in more detail that the rising base flow during the rainy season is not well represented by the model.The simulated low flow is frequently below the measured uncertainty band.During base flow conditions even the parameter uncertainty range does not overlap with the measured uncertainty range.This is not surprising since low flow was not in the focus of this study and therefore parameters that control low flow were not considered in the calibration process.The comparison of the hyetographs from all climate stations suggests that overestimated peaks during the rainy season are attributed to the spatial assignment of climate stations which was done using Thiessen polygons.As the polygon sizes are unevenly distributed (range from 0.3 to 82 km 2 ) this method may not be appropriate to account for localized precipitation events.Hence, local precipitation events may result in errors if this method is applied.Other interpolation methods as inverse distance are not implemented in the model code.

Erosion Modelling
The total simulated suspended sediment load responds sensitively to four model parameters (k f , k r , BKB, and DLSMAX) that were included in the LHS.
Figure 5 shows the simulated SSL having the highest sum of the performance indices for the calibration and validation period.The NSE is 0.4 and 0.2, the KGE 0.3 and 0.01 and the R 2 is 0.47 and 0.37 for calibration and validation respectively.The NSE is in the range of other SHETRAN studies (see Table 2).However, few studies used the erosion component and if so performance indices are sometimes not reported.Other physically based erosion models may perform worse or better (−0.75 ≤ NSE ≤ 0.94) as shown by model comparison studies of de Vente et al. [8] and Pandey et al. [9] but a comparison with other models is difficult due to differences regarding the model setup (spatial/temporal scale) and the chosen output variable [8].Given the various sources of measurement uncertainty a NSE of larger than 0.7 cannot be expected [8].Nevertheless, it has to be noted that especially for the validation period, quality measures are not satisfactory.This may be related to the large differences between 2014 and 2015 regarding the observed SSL.The measured annual erosion rate in 2015 is almost three times higher than in 2014 whereas the simulated erosion rate is only 1.6 times higher in 2015.The simulated and measured annual sediment yield is 700 and 970 t for 2014, 1045 and 2725 t for 2015.Reasons for the differences between simulated and observed SSL may be the setting of the erosion parameters as some values are markedly different from those found in literature.However, an adjustment of erosion parameters to different conditions as discussed by Bathurst [57] is not possible during continuous simulations.Another reason may be the hourly timestep of the precipitation input which may be too long since erosion often occurs at sub-hourly periods.
Parameter and measurement uncertainty are shown in Figure 5b.The assessment of the measured uncertainties of discharge and SSC results in relatively large combined uncertainty bands of SSL especially during peak flow periods.Possible measurement errors may be attributed to the discharge rating curve that does not cover the full range of recorded water levels and the calibration of the turbidity sensor that is also subject to uncertainty as shown in Section 2.5.1 [21,87].The large measured uncertainty is also reflected in higher R-factors and results in higher P-factors since simulated SSL and its uncertainty band are more often within the large measured uncertainty.This is also supported by the NSE m which is higher than the classical version.

Simulated Erosion Sources
Table 6 shows the relative contribution and the sediment yield of each source as simulated by SHETRAN.The interpretation of the results is associated with large uncertainties since results from fingerprinting analyses necessary to validate the simulated results are not yet available.Furthermore, knowledge on the erosion parameters is limited but the model output is strongly controlled by the parameterization.Among the sediment sources listed in Table 6 water contributes up to 1% to the sediment yield of the catchment.As the erosion parameters in SHETRAN are linked to the soil types an additional soil type would have been necessary to account for the conditions of areas covered by water.The number of soil types is limited in SHETRAN.Therefore, an additional soil type could not be implemented.
The simulated range of river erosion (including bank erosion and incision) dominates contributions from hillslope erosion.Between 68% and 89% of the simulated sediment loss are supplied from river erosion whereas 11% to 32% is eroded on the hillslopes.Simulated erosion rates for the calibration period range from 0.008 to 0.081 t/ha/year, for the entire catchment.The simulated values are low but in the range of the measured SSL (0.04-0.13 t/ha for 2014).However, the parameterization of a model to well simulate very small erosion rates is also challenging as discussed by Nearing [88].Hillslope erosion rates derived from 137 Cs measurement on hillslopes in the same area [25] are three orders of magnitude higher compared with the simulated hillslope erosion rates (0.005-0.022 t/ha).In the study by Roose [31] soil erosion rates measured under different experimental conditions in West Africa range from 0.01 to 90 t/ha/year.Walling et al. [89] assessed the sediment budget of a catchment in Zambia and measured 0.2 t/ha/year.The same study

Simulated Erosion Sources
Table 6 shows the relative contribution and the sediment yield of each source as simulated by SHETRAN.The interpretation of the results is associated with large uncertainties since results from fingerprinting analyses necessary to validate the simulated results are not yet available.Furthermore, knowledge on the erosion parameters is limited but the model output is strongly controlled by the parameterization.Among the sediment sources listed in Table 6 water contributes up to 1% to the sediment yield of the catchment.As the erosion parameters in SHETRAN are linked to the soil types an additional soil type would have been necessary to account for the conditions of areas covered by water.The number of soil types is limited in SHETRAN.Therefore, an additional soil type could not be implemented.
The simulated range of river erosion (including bank erosion and incision) dominates contributions from hillslope erosion.Between 68% and 89% of the simulated sediment loss are supplied from river erosion whereas 11% to 32% is eroded on the hillslopes.Simulated erosion rates for the calibration period range from 0.008 to 0.081 t/ha/year, for the entire catchment.The simulated values are low but in the range of the measured SSL (0.04-0.13 t/ha for 2014).However, the parameterization of a model to well simulate very small erosion rates is also challenging as discussed by Nearing [88].Hillslope erosion rates derived from 137 Cs measurement on hillslopes in the same area [25] are three orders of magnitude higher compared with the simulated hillslope erosion rates (0.005-0.022 t/ha).In the study by Roose [31] soil erosion rates measured under different experimental conditions in West Africa range from 0.01 to 90 t/ha/year.Walling et al. [89] assessed the sediment budget of a catchment in Zambia and measured 0.2 t/ha/year.The same study assessed the channel bank and gully contribution by using the fingerprint method to be in the order of 17%.Data collected in 2015 from plot measurements close to the study site indicate a range between 0.05 and 0.6 t/ha/year.Based on the comparison with measured ranges, hillslope erosion seems to be underestimated whereas the contribution of river bank and bed erosion is overestimated by the model.Knowledge on the relative contribution of hillslopes and rivers to the total catchment erosion is limited.The link between the sediment mobilization in the source area and the sediment yield measured at the outlet is difficult to study due to a lack of knowledge regarding the magnitude and residence time of sediment in storage [89,90].Additionally, information on the erodibility parameters of SHETRAN is often obtained from previous studies and has rarely been validated against measured plot data.Hence, setting the parameter range in order to reflect the catchment conditions is difficult and mainly based on the modelers' perception of the main erosion processes taking place in the catchment.

Catchment Distributed Erosion
Figure 6 shows the spatial pattern of soil erosion and deposition on hillslopes as simulated by the best SHETRAN simulation for the calibration period.Erosion ranges from 1.6 to almost 0 t/ha/year and deposition from 0 to 1.8 t/ha/year.Regarding erosion 55% of the grid cells are within the range of −0.03 to <0 t/ha/year.However, given the spatial heterogeneity of the controlling factors that are considered (slope, land use, soil properties, hydrological conditions) and the complex model approach it is difficult to clearly identify erosion hot spots and to explicitly attribute these hot spots to a single factor or reason.One example in this context is the hilly area in the western part of the catchment: We assumed the highest erosion rates would be found here due to the steep slopes and the partly practiced agriculture but simulated erosion is not especially high compared with other parts of the catchment.The small erosion rates are attributed to the low simulated surface runoff due to a small drainage area of these cells.Hence, rain drop detachment may be high but surface runoff is insufficient to transport the available sediment.High erosion rates simulated close to the river channel are a result of higher surface runoff simulated in cells close to the valley bottoms as a result of larger drainage areas and overbank flow.Field observations confirm overbank flow during large events but the areas are rather characterized by deposition of fine material and not necessarily by erosion as simulated.To get a more realistic representation of the vegetation close to the rivers (gallery forests), a fifth land use type with the same properties as savannah vegetation was introduced around the river links.Though this led to lower erosion rates, it negatively affected the performance of the simulated hydrograph and was therefore discarded.

Conclusions
This study applied the hydrological soil erosion model SHETRAN in a tropical West African catchment and investigated measurement and parameter uncertainty.From the investigations we can draw the following conclusions: (1) The performed uncertainty analyses of observed discharge reveals large uncertainty bands especially during peak flows (max.uncertainty from 17.3 (−34.1% of measured value) to 40.3 m 3 /s (+53.1% of measured value)) which was attributed on the one hand to the power law chosen for the rating curves and on the other hand to the sample properties.As a result of the intrinsic measurement errors and the error propagation the combined uncertainty of SSL is quite large (max.uncertainty from 2.8 to 8 kg/s).(2) Two hydrological parameters were tested regarding the sensitivity of the model response.
Whereas the ratio ETa/ETp affects total catchment runoff, the roughness coefficient KSTR has greater effect on the maximum runoff.Among the four tested erosion related model parameters the river bank erodibility coefficient BKR had the largest impact on the model response.
Parameter ranges of the overland flow erodibility coefficient kf and DLSMAX were quite low which was explained by the higher soil erodibility of the soils found in the study area.(3) The performance indices of simulated discharge are good (≥0.66) and comparable with other studies that used SHETRAN.Among these studies R 2 and NSE values above 0.5 are frequently reported.However, SHETRAN often underestimates base flow which could be explained by the missing calibration of hydrological subsurface parameters.Some peaks were not well represented due to the differences between real and model spatial representation of rainfall.The performance indices of the simulated SSL are comparable with the few studies that used SHETRAN to simulate soil erosion and that indicated model performance.As the calculation of SSC is based on the relation between turbidity and sediment concentration, input of organic material into the river following the burning of crop residues and grassland may lead to high turbidity readings although the measured weight is low [91,92].Thus, the mismatch between observed and simulated SSL at the start of the rainy season may also be explained by the method used to obtain the sedigraph.(4) The combined uncertainty assessment of measured and simulated discharge showed that SHETRAN frequently underestimates base flow despite large measured uncertainty bounds.

Conclusions
This study applied the hydrological soil erosion model SHETRAN in a tropical West African catchment and investigated measurement and parameter uncertainty.From the investigations we can draw the following conclusions: (1) The performed uncertainty analyses of observed discharge reveals large uncertainty bands especially during peak flows (max.uncertainty from 17.3 (−34.1% of measured value) to 40.3 m 3 /s (+53.1% of measured value)) which was attributed on the one hand to the power law chosen for the rating curves and on the other hand to the sample properties.As a result of the intrinsic measurement errors and the error propagation the combined uncertainty of SSL is quite large (max.uncertainty from 2.8 to 8 kg/s).(2) Two hydrological parameters were tested regarding the sensitivity of the model response.
Whereas the ratio ETa/ETp affects total catchment runoff, the roughness coefficient KSTR has greater effect on the maximum runoff.Among the four tested erosion related model parameters the river bank erodibility coefficient BKR had the largest impact on the model response.Parameter ranges of the overland flow erodibility coefficient k f and DLSMAX were quite low which was explained by the higher soil erodibility of the soils found in the study area.(3) The performance indices of simulated discharge are good (≥0.66) and comparable with other studies that used SHETRAN.Among these studies R 2 and NSE values above 0.5 are frequently reported.However, SHETRAN often underestimates base flow which could be explained by the missing calibration of hydrological subsurface parameters.Some peaks were not well represented due to the differences between real and model spatial representation of rainfall.The performance indices of the simulated SSL are comparable with the few studies that used SHETRAN to simulate soil erosion and that indicated model performance.As the calculation of SSC is based on the relation between turbidity and sediment concentration, input of organic material into the river following the burning of crop residues and grassland may lead to high turbidity readings although the measured weight is low [91,92].Thus, the mismatch between observed and simulated SSL at the start of the rainy season may also be explained by the method used to obtain the sedigraph.
(4) The combined uncertainty assessment of measured and simulated discharge showed that SHETRAN frequently underestimates base flow despite large measured uncertainty bounds.
The modified NSE m used to include both uncertainties in the quality assessment showed that the overlapping areas of distributions are rarely observed and small.As a result of the large uncertainty of observed SSL the model uncertainty is almost always within the range of measured uncertainty bounds.This is also reflected by a slightly higher NSE m in comparison with the traditional NSE.The erosion sources simulated by SHETRAN do not correspond with the sources reported in the literature.The contribution of river bank and bed erosion may be too high and the erosion on hillslopes too low.However, knowledge on this point is limited.Results from fingerprint analyses may help to validate the simulated output.
We showed that the physically based spatially distributed erosion model SHETRAN offers chances and challenges.SHETRAN provides a better representation of erosion processes, especially in environments that are characterized by low frequency and high magnitude erosion events, than low-complexity models that focus on mean annual erosion rates (such as the USLE).So far, the application to larger river catchments is limited due to the incomplete knowledge on the model parameters and the limited availability of model input data at a large scale and with appropriate resolution.Yet, the modeling results obtained here help to improve the parameterization of large-scale, low-complexity erosion models and to improve the representation of the strong temporal variability of erosion rates in the Sudano savanna zone of Africa.Consequently, modeling erosion on a small scale is needed to optimize parameter estimation and process understanding which in turn helps to improve large scale modeling [93].Erosion modeling with SHETRAN should therefore focus on better and broader knowledge on the erosion parameters, including the definition of parameter ranges for different environmental conditions.The application of SHETRAN to erosion plots may be an opportunity to better assess these parameter values.

Figure 1 .
Figure 1.Location map of the Dano catchment: (a) location of the catchment and Burkina Faso in West Africa; (b) slope of the catchment; (c) model catchment; (d) land use map [53]; (e) soil map (data base: soil survey done by Ozias Hounkpatin, University of Bonn, Institute of Crop Science and Resource Conservation, Soil Science and Soil Ecology).

Figure 1 .
Figure 1.Location map of the Dano catchment: (a) location of the catchment and Burkina Faso in West Africa; (b) slope of the catchment; (c) model catchment; (d) land use map [53]; (e) soil map (data base: soil survey done by Ozias Hounkpatin, University of Bonn, Institute of Crop Science and Resource Conservation, Soil Science and Soil Ecology).

Figure 3 .
Figure 3. Scatter plots showing the sensitivity of water discharge Q to the ratio of actual to potential evapotranspiration (ETa/ETp) (a); and surface roughness (KSTR) (b); and the sensitivity of the catchment suspended sediment load (SSL) to the raindrop erodibility coefficient (kr) (c); the overland flow erodibility (kf) (d); the bank erodibility (BKB) (e); and the depth of loose sediment (DLSMAX) (f).Red dashed lines indicate the base run used for comparison.SI90max indicates the SI90 for the maximum discharge.

Figure 3 .
Figure 3. Scatter plots showing the sensitivity of water discharge Q to the ratio of actual to potential evapotranspiration (ETa/ETp) (a); and surface roughness (KSTR) (b); and the sensitivity of the catchment suspended sediment load (SSL) to the raindrop erodibility coefficient (k r ) (c); the overland flow erodibility (k f ) (d); the bank erodibility (BKB) (e); and the depth of loose sediment (DLSMAX) (f).Red dashed lines indicate the base run used for comparison.SI 90max indicates the SI 90 for the maximum discharge.

Figure 4 .
Figure 4. Observed and simulated daily discharge (Q) over the calibration and validation period (a); and over a selected period for the calibration (b); and validation (c) period.The measured 95%-confidence interval is given for both periods, simulated parameter uncertainty only for the calibration year.

Figure 4 .
Figure 4. Observed and simulated daily discharge (Q) over the calibration and validation period (a); and over a selected period for the calibration (b); and validation (c) period.The measured 95%-confidence interval is given for both periods, simulated parameter uncertainty only for the calibration year.

Water 2017, 9 , 101 14 of 21 Figure 5 .
Figure 5. Observed and simulated daily suspended sediment load (SSL) over the calibration and validation period (a); and a selected period for the calibration (b); and validation (c) period.Measurement uncertainty is given for both periods, parameter uncertainty only for the calibration.

Figure 5 .
Figure 5. Observed and simulated daily suspended sediment load (SSL) over the calibration and validation period (a); and a selected period for the calibration (b); and validation (c) period.Measurement uncertainty is given for both periods, parameter uncertainty only for the calibration.

Figure 6 .
Figure 6.Annual soil erosion (-) and deposition rate (+) on hillslopes as simulated for the calibration period.

Figure 6 .
Figure 6.Annual soil erosion (-) and deposition rate (+) on hillslopes as simulated for the calibration period.

Table 1 .
Selected studies on soil erosion modelling in West Africa.NSE is the Nash-Sutcliff efficiency.

Table 2 .
Selected studies using SHETRAN for water flow and/or sediment flow simulations.

Table 3 .
Range of the measured parameters for the years 2014 and 2015.SSC refers to suspended sediment concentration and SSL to suspended sediment load.

Table 4 .
Applied datasets and required inputs for SHETRAN.α and n are van Genuchten empirical parameters; 2 K sat refers to the saturated hydraulic conductivity; 3 θ sat to the saturated water content; 4 θ res to the residual water content; 5 LAI to the leaf area index; and6ETp/ETa ratio to the ratio of potential evapotranspiration to actual evapotranspiration.

Table 4 .
Applied datasets and required inputs for SHETRAN.

Table 5 .
Soil, land use and erosion parameters in SHETRAN.

Table 6 .
Relative contribution and specific sediment yield of the different erosion sources as simulated by the best SHETRAN run for the year 2014 (calibration) and 2015 (validation).Min. and max.indicate the ranges of the ten considered simulations.
Notes:1River bank and bed erosion rates are calculated in reference to the total catchment size.