Estimating suspended sediment concentrations from turbidity measurements and the calibration problem

In situ turbidity meters are being increasingly used to generate continuous records of suspended sediment concentration in rivers. However, the usefulness of the information obtained depends heavily on the existence of a close relationship between fluctuations in suspended sediment concentration and turbidity and the calibration procedure that relates suspended sediment concentration to the turbidity meter's signal. This study assesses the relationship between suspended sediment concentration and turbidity for a small (1·19 km²) rural catchment in southern Brazil and evaluates two calibration methods by comparing the estimates of suspended sediment concentration obtained from the calibrated turbidity readings with direct measurements obtained using a USDH 48 suspended sediment sampler. With the first calibration method, the calibration relationship is derived by relating the turbidity readings to simultaneous measurements of concentration obtained from suspended sediment samples collected from the vicinity of the turbidity probe during flood events. With the second method, the calibration is based on the readings obtained from the turbidity meter when the probe immersed in samples of known concentration prepared using soils collected from the catchment. Overall, there was a close link between fluctuations in suspended sediment concentration and turbidity in the stream at the outlet of the catchment, and the estimates of sediment concentration obtained using the first calibration method corresponded closely with the conventionally measured sediment concentrations. However, use of the second calibration method introduced appreciable errors. When the estimated sediment concentrations were compared with the measured values, the mean errors were ± 122 mg l^?1 and + 601 mg l^?1 for the first and second calibration procedures respectively. Copyright © 2007 John Wiley & Sons, Ltd.     

Reducing the impacts of flood-induced reservoir turbidity on a regional water supply system

This paper proposed an integrated simulation model to incorporate the impact of flood-induced reservoir turbidity into water supply. The integrated model includes a regional water allocation model and a one-dimensional settling model of cohesive particles based on Kynch’s theory. It simulates the settling of sediment flocculation in a turbid reservoir. The restrictions of water supply during floods is mimicked by simulating turbidity profiles for control points and then quantifying the associated treatment capability of raw water in the regional water allocation model for each time step. This framework can simulate shortages caused by flood-induced high turbidity as well as extended droughts, thus provide a basis for comprehensive evaluations of emergent and regular water supply facilities. A case study of evaluating different measures to mitigate the impact of turbid reservoir on water supply in northern Taiwan is presented to demonstrate the efficacy of the proposed approach.     

Experimental and numerical modelling of sedimentation in a rectangular shallow basin

Numerical simulation of flows in shallow reservoirs has to be checked for its consistency in predicting real flow conditions and sedimentation patterns. Typical flow patterns may exhibit flow separation at the inlet, accompanied by several recirculation and stagnation areas all over the reservoir surface. The aim of the present research project is to study the influence of the geometry of a reservoir on sediment transport and deposition numerically and experimentally, focusing on a prototype reservoir depth between 5 and 15 m as well as suspended sediment transport.
A series of numerical simulations is presented and compared with scaled laboratory experiments, with the objective of testing the sensitivity to different flow and sediment parameters and different turbulence closure schemes. Different scenarios are analyzed and a detailed comparison of preliminary laboratory tests and some selected simulations are presented.
The laboratory experiments show that suspended sediment transport and deposition are determined by the initial flow pattern and by the upstream and downstream boundary conditions. In the experiments, deposition in the rectangular basin systematically developed along the left bank, although inflow and outflow were positioned symmetrically along the centre of the basin. Three major horizontal eddies developed influencing the sediment deposition pattern. Although asymmetric flow patterns are privileged, a symmetric pattern can appear from time to time. This particular behaviour could also be reproduced by a two-dimensional depth-averaged flow and sediment transport model （CCHE2D）. The paper presents numerical simulations using different turbulence closure schemes （k-ε and eddy viscosity models）. In spite of the symmetric setup, these generally produced an asymmetric flow pattern that can easily switch sides depending on the assumptions made for the initial and boundary conditions. When using the laboratory experiment as a reference, the most reliable numerical results have been obtai     

Evaluating suspended sediment and turbidity reduction projects in a glacially conditioned catchment through dynamic regression and fluvial process-based modelling

Elevated turbidity (T_n) and suspended sediment concentrations (SSC) during and following flood events can degrade water supply quality and aquatic ecosystem integrity. Streams draining glacially conditioned mountainous terrain, such as those in the Catskill Mountains of New York State, are particularly susceptible to high levels of T_n and SSC sourced from erosional contact with glacial-related sediment. This study forwards a novel approach to evaluate the effectiveness of stream restoration best management practices (BMPs) meant to reduce stream T_n and SSC, and demonstrates the approach within the Stony Clove sub-basin of the Catskills, a water supply source for New York City. The proposed approach is designed to isolate BMP effects from natural trends in T_n and SSC caused by trends in discharge and shifts in average T_n or SSC per unit discharge (Q) following large flood events. We develop Dynamic Linear Models (DLMs) to quantify how T_n-Q and SSC-Q relationships change over time at monitoring stations upstream and downstream of BMPs within the Stony Clove and in three other sub-basins without BMPs, providing observational evidence of BMP effectiveness. A process-based model, the River Erosion Model, is then developed to simulate natural, hydrology-driven SSC-Q dynamics in the Stony Clove sub-basin (absent of BMP effects). We use DLMs to compare the modelled and observed SSC-Q dynamics and isolate the influence of the BMPs. Results suggest that observed reductions in SSC and T_n in the Stony Clove sub-basin have been driven by a combination of declining streamflow and the installed BMPs, confirming the utility of the BMPs for the monitored hydrologic conditions.     

Targeting of intervention areas to reduce reservoir sedimentation in the Tana catchment (Kenya) using SWAT

Abstract

A measurement campaign was carried out in the Upper Tana basin (Kenya) to quantify soil erosion and reservoir sedimentation rates, including a bathymetric reservoir survey and sediment load sampling during one year. Then, distributed soil erosion modelling was performed to study sediment budgets throughout the basin and to evaluate the potential of upstream erosion control through vegetated contour strips and check dams. Finally, the areas where these measures would be most effective were identified and local stakeholder associations to implement them were prioritized. The influence of the scale of implementation was evaluated by using the model to consider three adoption scenarios. This study illustrates the relevance of distributed erosion models to target erosion control measures when sufficient information on the eroding areas is available from field surveys. Bathymetric surveys were fundamental to validate the long-term model response, while point measurements were valuable to verify the spatial variability of model predictions.

Editor Z.W. Kundzewicz; Associate editor G. Mahé

Citation Hunink, J.E., Niadas, I.A., Antonaropoulos, P., Droogers, P., and de Vente, J., 2013. Targeting of intervention areas to reduce reservoir sedimentation in the Tana catchment (Kenya) using SWAT. Hydrological Sciences Journal, 58 (3), 600–614.     

Significance of particle interaction to the modelling of cohesive sediment transport in rivers

Aggregation processes of fine sediments have rarely been integrated in numerical simulations of cohesive sediment transport in riverine systems. These processes, however, can significantly alter the hydrodynamic characteristics of suspended particulate matter (SPM), modifying the particle settling velocity, which is one of the most important parameters in modelling suspended sediment dynamics. The present paper presents data from field measurements and an approach to integrate particle aggregation in a hydrodynamic sediment transport model. The aggregation term used represents the interaction of multiple sediment classes (fractions) with corresponding multiple deposition behaviour. The k–ε–turbulence model was used to calculate the coefficient of vertical turbulent mixing needed for the two‐dimensional vertical‐plane simulations. The model has been applied to transport and deposition of tracer particles and natural SPM in a lake‐outlet lowland river (Spree River, Germany). The results of simulations were evaluated by comparison with field data obtained for two levels of river discharge. Experimental data for both discharge levels showed that under the prevailing uniform hydraulic conditions along the river reach, the settling velocity distribution did not change significantly downstream, whereas the amount of SPM declined. It was also shown that higher flow velocities (higher fluid shear) resulted in higher proportions of fast settling SPM fractions. We conclude that in accordance with the respective prevailing turbulence structures, typical aggregation mechanisms occur that continuously generate similar distribution patterns, including particles that settle toward the river bed and thus mainly contribute to the observed decline in the total SPM concentration. In order to determine time‐scales of aggregation and related mass fluxes between the settling velocity fractions, results of model simulations were fitted to experimental data for total SPM concentration and of settling velocity frequency distributions. The comparison with simulations for the case of non‐interacting fractions clearly demonstrated the practical significance of particle interaction for a more realistic modelling of cohesive sediment and contaminant transport. Copyright © 2004 John Wiley & Sons, Ltd.     

Assessment of reservoir sedimentation using remote sensing and recommendations for desilting Patratu Reservoir,India

ABSTRACT

This study describes the assessment of reservoir sedimentation of the Patratu Reservoir using Satellite Remote Sensing (SRS). The sedimentation assessment was carried out using satellite data and reservoir water level data from 2006 to 2012. Water spread area was analysed from satellite data. The Normalized Difference Water Index (NDWI) has been used to delineate open water features and to enhance the presence of water surface in satellite imagery of the Patratu Reservoir. Water spread area of the reservoir at a particular elevation on the date of the passing of the satellite was used to develop an elevation-area curve. For the present case, fluctuation of water level was found to vary from 387.096 to 406.152 m. The linear interpolation/extrapolation technique has been employed to assess the water spread area of Patratu Reservoir at different elevations. Further, these areas were used to compute the live storage capacity of the reservoir between two elevations by the Prismoidal formula. From the study, it was found that due to sedimentation, the live storage capacity of Patratu Reservoir has reduced from 101.95 to 89.96 hm³, thus showing capacity loss of 11.76% in a span of 44 years. To increase the live storage capacity of the reservoir it is proposed to adopt manual and mechanical digging combined with flushing for desilting of the deposited sediment.

EDITOR D. Koutsoyiannis ASSOCIATE EDITOR S. Kanae     

Development of a real-time forecasting model for turbidity current arrival time to improve reservoir desilting operation

ABSTRACT

Among various strategies for sediment reduction, venting turbidity currents through dam outlets can be an efficient way to reduce suspended sediment deposition. The accuracy of turbidity current arrival time forecasts is crucial for the operation of reservoir desiltation. A turbidity current arrival time (TCAT) model is proposed. A multi-objective genetic algorithm (MOGA), a support vector machine (SVM) and a two-stage forecasting technique are integrated to obtain more effective long lead-time forecasts of inflow discharge and inflow sediment concentration. The multi-objective genetic algorithm (MOGA) is applied for determining the optimal inputs of the forecasting model, support vector machine (SVM). The two-stage forecasting technique is implemented by adding the forecasted values to candidate inputs for improving the long lead-time forecasting. Then, the turbidity current arrival time from the inflow boundary to the reservoir outlet is calculated. To demonstrate the effectiveness of the TCAT model, it is applied to Shihmen Reservoir in northern Taiwan. The results confirm that the TCAT model forecasts are in good agreement with the observed data. The proposed TCAT model can provide useful information for reservoir sedimentation management during desilting operations.     

Fine sediment delivery and transfer in lowland catchments: modelling suspended sediment concentrations in response to hydrological forcing

Fine sediment delivery to and storage in stream channel reaches can disrupt aquatic habitats, impact river hydromorphology, and transfer adsorbed nutrients and pollutants from catchment slopes to the fluvial system. This paper presents a modelling tool for simulating the time‐dependent response of the fine sediment system in catchments, using an integrated approach that incorporates both land phase and in‐stream processes of sediment generation, storage and transfer. The performance of the model is demonstrated by applying it to simulate in‐stream suspended sediment concentrations in two lowland catchments in southern England, the Enborne and the Lambourn, which exhibit contrasting hydrological and sediment responses due to differences in substrate permeability. The sediment model performs well in the Enborne catchment, where direct runoff events are frequent and peak suspended sediment concentrations can exceed 600 mg l^?1. The general trends in the in‐stream concentrations in the Lambourn catchment are also reproduced by the model, although the observed concentrations are low (rarely exceeding 50 mg l^?1) and the background variability in the concentrations is not fully characterized by the model. Direct runoff events are rare in this highly permeable catchment, resulting in a weak coupling between the sediment delivery system and the catchment hydrology. The generic performance of the model is also assessed using a generalized sensitivity analysis based on the parameter bounds identified in the catchment applications. Results indicate that the hydrological parameters contributing to the sediment response include those controlling (1) the partitioning of runoff between surface and soil zone flows and (2) the fractional loss of direct runoff volume prior to channel delivery. The principal sediment processes controlling model behaviour in the simulations are the transport capacity of direct runoff and the in‐stream generation, storage and release of the fine sediment fraction. The in‐stream processes appear to be important in maintaining the suspended sediment concentrations during low flows in the River Enborne and throughout much of the year in the River Lambourn. Copyright © 2007 John Wiley & Sons, Ltd.     

A revised soil erosion budget for India: role of reservoir sedimentation and land‐use protection measures

Among the different controls of erosion budget at basin level, the relative impact of dams and land management is yet to be investigated. In this paper, the impact of dams on sediment yield has been assessed by using a conceptual modelling framework which considers the gross erosion and the cascade of dams constructed on a river network. The sediment budget has been estimated based on the gross erosion, deposition of sediment in reservoirs, and sediment yields of 23 mainland river basins of India. The gross erosion of the country is estimated as 5.11 ± 0.4 Gt yr^?1 or 1559 t km^?2 yr^?1, out of which 34.1 ± 12% of the total eroded soil is deposited in the reservoirs, 22.9 ± 29% is discharged outside the country (mainly to oceans), and the remaining 43.0 ± 41% is displaced within the river basins. The river basins of northern India contribute about 81% of the total sediment yield from landmass while the share of southern river basins is 19%. The components of revised sediment budget for India are prominently influenced by the sediment trapped in reservoirs and the treatment of catchment areas by soil and water conservation measures. Analysis of sediment deposition in 4937 reservoirs indicated the average annual percentage capacity loss as 1.04% though it varies from 0.8% to >2% per year in smaller dams (1–50 Mm³ capacity) and from <0.5% to 0.8% per year in larger dams (51 to >1000 Mm³ capacity). Siltation of smaller dams poses a serious threat to their ecosystem services as they cater to a wider population for domestic, agricultural, and industrial purposes. Amongst the environment controls, land use significantly impacts the gross erosion rate and specific sediment yield as compared to climatic and topographic parameters. However, to analyse their integrated effect on the complex processes of sediment fluxes in a basin, further research efforts are needed. Copyright © 2016 John Wiley & Sons, Ltd.     

Challenges to the representation of suspended sediment transfer using a depth‐averaged flux

The sediment saturation recovery process (i.e. the adaptation of suspended sediment concentration [SSC] to local forcing) is the main feature of the non‐equilibrium suspended sediment transport (SST) frequently occurring in fluvial, estuarine and coastal waters. In order to quantitatively describe this phenomenon, a series solution is analytically derived, including the evolution of both vertical SSC profile and near‐bed sediment flux (NBSF), and is verified by net erosion and net deposition experiments, respectively. The results suggest that the sediment saturation recovery process involves vertically varying fluxes that are not represented correctly by depth‐averaging. Consequently, a vertical two‐dimensional (2D) combined scheme is established and applied respectively in to a dredged trench and to a sand wave feature to demonstrate this argument. By analyzing the variations of the calculated depth‐averaged SSC and NBSF we reveal that the equilibrium state presented by the sediment carrying capacity (SCC) form of the NBSF, which is usually applied in depth‐integrated SST models, lags behind the actual dynamic bed equilibrium state. Moreover, the key factor α, the so‐called saturation recovery coefficient within this form, is not only a function of local Rouse number but also is influenced by the local SSC profile. Finally, a three‐dimensional (3D) non‐orthogonal curvilinear body‐fitted SST model is developed and validated in the Yangtze estuary, China, combined with the in situ hourly hydrographic data from August 14–15, 2007 during spring tide in the wet season. Model results confirm that the vertically varying sediment saturation recovery process, the discrepancies between the actual and SCC form of NBSF and non‐constant value of α are significant in actual real geomorphic cases. The quantitative morphological change resulting from variations in environmental conditions may not be correctly represented by uncorrected depth‐integrated SST models if they do not treat the effects of vertical motion on the sediment saturation recovery process. Copyright © 2016 John Wiley & Sons, Ltd.     

An examination of geochemical modelling approaches to tracing sediment sources incorporating distribution mixing and elemental correlations

The identification of sediment sources is fundamental to the management of increasingly scarce water resources. Tracing the origin of sediment with elemental geochemistry is a well‐established approach to determining sediment provenance. Fundamental to the confident apportionment of sediment to their lithogenic sources is the modelling process. Recent approaches have incorporated distributions throughout the modelling process including source contribution terms for two end‐member sources. The shift from modelling source samples to modelling samples drawn from distributions has removed relationships, including potential correlations between elemental concentrations, from the modelling process. Here, we present a novel modelling approach that re‐incorporates correlations between elemental concentrations and models distributions for source contribution terms for multiple source end members. Artificial mixtures, based on catchment sources samples, were created to test the accuracy of this correlated distribution model and also examine modelling approaches used in the literature. The most accurate model incorporates correlations between elements, uses the absolute mixing model difference and does not use any weighting. This model was then applied to identify the sources of sediment in three South East Queensland catchments and demonstrated that Quaternary Alluvium is the most dominant source of sediment in these catchments (μ 44%, σ 12%). This study demonstrates that it is important to understand how different weightings may impact modelling results. Copyright © 2014 John Wiley & Sons, Ltd.     

Stream turbidity responses to storm events in a pristine rainforest watershed on the Coral Coast of southern Fiji

On the Coral Coast of Viti Levu Island in Fiji, inadequate knowledge of suspended sediment delivery patterns in small pristine coastal watersheds hinders any future assessment of accelerated erosion in disturbed areas nearby. This study adopts a rainfall–stream turbidity monitoring approach in the Votua Creek, which drains a small, steep but minimally-disturbed coastal rainforest catchment. Storm rainfall characteristics, stream depth and water turbidity were continuously monitored over one complete Fiji wet season from October 2009 to April 2010. The aim was to evaluate whether these parameters provide sufficient information to illustrate basic features of storm–sediment transport responses, in the case of limited stream gauging and very simple sediment rating curves. This is important because Pacific Island nations like Fiji do not have the resources to initiate long-term gauging and sediment sampling pro-grammes across numerous small catchments.
A significant power function demonstrates that turbidity (T) is a suitable proxy for total suspended solids (TSS) for turbidity measurements above 5 NTU, with TSS?0.930T1.111 (r?0.98, Po0.001). Over the study period, 10 individual storms 11.2–120.1 mm in size produced a‘significant turbidity response’ (STR) in the Votua Creek. Rainfall parameters (totals and intensities) showed positive linear relationships (r ? 0.72–0.94) with stream turbidity parameters (mean, maximum, duration), whilst relationships of similar strength (r ? 0.76–0.98) were also derived between stream flow depth and turbidity. This implies that for small rainforest watersheds in Fiji, rainfall parameters offer no substantial disadvantage over flow as predictors of stream sediment responses to major storms. Event-based analysis revealed that negative (anticlockwise) hysteresis is a typical flow–turbidity pattern for STR events. Negative hysteresis is produced when secondary episodes of renewed (heavy) rainfall occur after maximum intensity, in the later phase of storm events. Tropical Cyclone Mick in December 2009 generated the largest flood and the greatest turbidity response (Tmax ? 1021 NTU, Tmean ? 207 NTU). This concurs with earlier work confirming that tropical cyclones are the most important events for sediment transport in Fiji stream networks.     

Erosion modelling for land management in the Tahoe basin,USA: scaling from plots to forest catchments

Abstract

Recent developments in hydrological modelling of river basins are focused on prediction in ungauged basins, which implies the need to improve relationships between model parameters and easily-obtainable information, such as satellite images, and to test the transferability of model parameters. A large-scale distributed hydrological model is described, which has been used in several large river basins in Brazil. The model parameters are related to classes of physical characteristics, such as soil type, land use, geology and vegetation. The model uses two basin space units: square grids for flow direction along the basin and GRU—group response units—which are hydrological classes of the basin physical characteristics for water balance. Expected ranges of parameter values are associated with each of these classes during the model calibration. Results are presented of the model fitting in the Taquari-Antas River basin in Brazil (26 000 km² and 11 flow gauges). Based on this fitting, the model was then applied to the Upper Uruguay River basin (52 000 km²), having similar physical conditions, without any further calibration, in order to test the transferability of the model. The results in the Uruguay basin were compared with recorded flow data and showed relatively small errors, although a tendency to underestimate mean flows was found.     

Hydro‐meteorological drivers and sources of suspended sediment flux in the pro‐glacial zone of the retreating Castle Creek Glacier,Cariboo Mountains,British Columbia,Canada

Glaciers are major agents of erosion that increase sediment load to the downstream fluvial system. The Castle Creek Glacier, British Columbia, Canada, has retreated ~1.0 km in the past 70 years. Suspended sediment concentration (SSC) and streamflow (Q) were monitored independently at five sites within its pro‐glacial zone over a 60 day period from July to September 2011, representing part of the ablation season. Meteorological data were collected from two automatic weather stations proximal to the glacier. The time‐series were divided into hydrologic days and the shape and magnitude of the SSC response to hydro‐meteorological conditions (‘cold and wet’, ‘hot and dry’, ‘warm and damp’, and ‘storm’) were categorized using principal component analysis (PCA) and cluster analysis (CA). Suspended sediment load (SSL) was computed and summarized for the categories. The distribution of monitoring sites and results of the multivariate statistical analyses describe the temporal and spatial variability of suspended sediment flux and the relative importance of glacial and para‐glacial sediment sources in the pro‐glacial zone. During the 2011 study period, ~ 60% of the total SSL was derived from the glacial stream and sediment deposits proximal to the terminus of the glacier; during ‘storm’ events, that contribution dropped to ~40% as the contribution from diffuse and point sources of sediment throughout the pro‐glacial zone and within the meltwater channels increased. While ‘storm’ events accounted for just 3% of the study period, SSL was ~600% higher than the average over the monitoring period, and ~20% of the total SSL was generated in that time. Determining how hydro‐meteorological conditions and sediment sources control sediment fluxes will assist attempts to predict how pro‐glacial zones respond to future climate changes. Copyright © 2015 John Wiley & Sons, Ltd.     

Application of the honey-bees mating programming (HBMP) algorithm to sediment concentration modelling

Abstract

The honey-bees mating programming (HBMP) algorithm is introduced as a novel tool for predicting suspended sediment concentration for the Mad River catchment near Arcata, USA. The paper also applies gene expression programming (GEP) as a comparison and shows that these two approaches can the produce transparent, nonlinear relationships between the independent and dependent variables. Some modifications have been made to the HBMP algorithm to improve its capability and efficiency. The results achieved from this method and GEP are compared with two different sediment rating curves based on regression techniques. The findings show that the results from both the HBMP and GEP methods are promising and outperform the results obtained from the sediment rating curves.

Editor D. Koutsoyiannis; Associate editor L. See     

Relationship between texture and fractions of inorganic phosphorus in the surface sediment of a reservoir

Inorganic P was fractionated into three categories (NaOH-, BD- and HCI-extractable) for 30 surface sediment samples in La Minilla Reservoir (Sevilla, Spain). The amount of reactive P extracted with NaOH (NaOH-RP) and with HCI (HCI-RP) correlated in a multivariate regression with the clay and sand content of the sediment. This multivariate function should aid in predicting the amount of phosphorus available to the sediment organisms, and it can also contribute to the knowledge of the phosphorus budget of the reservoir.Supported by a Grant from CSIC-CONICET.     

Magnitude and variation in the contribution of bank erosion to the suspended sediment load of the River Severn,UK

The wide range of studies describing the role of bank erosion in fluvial sediment supply have mostly lumped amounts of bank erosion into coarse temporal units, such as years. This paper investigates sediment yields from individual bank erosion events within the upper River Severn, UK (basin area 380 km²). Manual erosion pins and photo-electronic erosion pins were used to estimate bank erosion, and turbidity meters were used to determine suspended sediment transport. At the annual time-scale, the silt-clay fraction of bank-derived sediment accounted for an equivalent of 17 per cent of the suspended load, increasing to an average of 38 per cent at the monthly timescale, and then to an average of 64 per cent at the event timescale. This research highlighted that for an upland catchment, bank erosion was an important supply of suspended sediment, and that for some flood events bank erosion can supply more sediment than is transported. © 1997 John Wiley & Sons, Ltd.     

Assessing the performance of a physically based hydrological model using a proxy‐catchment approach in an agricultural environment

Physically based models are useful frameworks for testing intervention strategies designed to reduce elevated sediment loads in agricultural catchments. Evaluating the success of these strategies depends on model accuracy, generally established by a calibration and evaluation process. In this contribution, the physically based SHETRAN model was assessed in two similar U.K. agricultural catchments. The model was calibrated on the Blackwater catchment (18 km²) and evaluated in the adjacent Kit Brook catchment (22 km²) using 4 years of 15 min discharge and suspended sediment flux data. Model sensitivity to changes in single and multiple combinations of parameters and sensitivity to changes in digital elevation model resolution were assessed. Model flow performance was reasonably accurate with a Nash–Sutcliffe efficiency coefficient of 0.78 in Blackwater and 0.60 in Kit Brook. In terms of event prediction, the mean of the absolute percentage of difference (μAbs_diff) between measured and simulated flow volume (Qv), peak discharge (Qp), sediment yield (Sy), and peak sediment flux (Sp) showed larger values in Kit Brook (48% [Qv], 66% [Qp], 298% [Sy], and 438% [Sp]) compared with the Blackwater catchment (30% [Qv], 41% [Qp], 106% [Sy], and 86% [Sp]). Results indicate that SHETRAN can produce reasonable flow prediction but performs less well in estimation of sediment flux, despite reasonably similar hydrosedimentary behaviour between catchments. The sensitivity index showed flow volume sensitive to saturated hydraulic conductivity and peak discharge to the Strickler coefficient; sediment yield was sensitive to the overland flow erodibility coefficient and peak sediment flux to raindrop/leaf soil erodibility coefficient. The multiparameter sensitivity analysis showed that different combinations of parameters produced similar model responses. Model sensitivity to grid resolution presented similar flow volumes for different digital elevation model resolutions, whereas event peak and duration (for both flow and sediment flux) were highly sensitive to changes in grid size.