Hydraulic model calibration for extreme floods in bedrock‐confined channels: case study from northern Thailand

Palaeoflood reconstructions based on stage evidence are typically conducted in data‐poor field settings. Few opportunities exist to calibrate the hydraulic models used to estimate discharge from this evidence. Consequently, an important hydraulic model parameter, the roughness coefficient (e.g. Manning's n), is typically estimated by a range of approximate techniques, such as ‘visual estimation’ and semi‐empirical equations. These techniques contribute uncertainty to resulting discharge estimates, especially where the study reach exhibits sensitivity in the discharge–Manning's n relation. We study this uncertainty within a hydraulic model for a large flood of known discharge on the Mae Chaem River, northern Thailand. Comparison of the ‘calibrated’ Manning's n with that obtained from semi‐empirical equations indicates that these underestimate roughness. Substantial roughness elements in the extra‐channel zone, inundated during large events, contribute significant additional sources of flow resistance that are captured neither by the semi‐empirical equations, nor by existing models predicting stage–roughness variations. This bedrock channel exhibits a complex discharge–Manning's n relation, and reliable estimates of the former are dependent upon realistic assignment of the latter. Our study demonstrates that a large recent flood can provide a valuable opportunity to constrain this parameter, and this is illustrated when we model a palaeoflood event in the same reach, and subsequently examine the magnitude–return period consequences of discharge uncertainty within a flood frequency analysis, which contributes its own source of uncertainty. Copyright © 2005 John Wiley & Sons, Ltd.     

Hydrological characteristics of vegetated river flows: a laboratory flume study

Abstract

Laboratory flume experiments were undertaken to measure the vertical profiles of mean flow velocity for three different flow discharges and four different stem densities of Hydrilla verticillata. The data were used to calculate three parameters, namely Manning's roughness coefficient, the Reynolds number and the Froude number. In addition, empirical equations were obtained for the vertical distribution of measured flow velocity within the transitional zone and above the plant canopy. The results show that: (a) the vertical distribution of measured flow velocity exhibits three zone profiles; (b) Manning's roughness coefficient decreases with increasing depth-averaged flow velocity; (c) the relationship between Manning's roughness coefficient and the depth-averaged flow velocity is within the smooth left inverse curve; (d) Manning's roughness coefficient significantly changes with increasing density of Hydrilla; (e) the Froude number is independent of the density of Hydrilla; and (f) both the Reynolds number and the Froude number increase with increasing depth-averaged flow velocity.

Citation Shi, J.Z., Li, Y.-H., Hughes, J.M.R., and Zhao, M., 2013. Hydrological characteristics of vegetated river flows: a laboratory flume study. Hydrological Sciences Journal, 58 (5), 1047–1058.

Editor Z.W. Kundzewicz     

Modelling the bulk flow of a bedrock‐constrained,multi‐channel reach of the Mekong River,Siphandone, southern Laos

The general nature of bulk flow within bedrock single‐channel reaches has been considered by several studies recently. However, the flow structure of a bedrock‐constrained, large river with a multiple channel network has not been investigated previously. The multiple channel network of the Siphandone wetlands in Laos, a section of the Mekong River, was modelled using a steady one‐dimensional hydraulic model. The river network is characterized by a spatially‐varying channel‐form leading to significant changes in the bulk flow properties between and along the channels. The challenge to model the bulk flow in such a remote region was the lack of ideal boundary conditions. The flow models considered both low flow, high inbank and overbank flows and were calibrated using SPOT satellite sensor imagery and limited field data concerning water levels. The application of the model highlighted flow characteristics of a large multi‐channel network and also further indicated the field data that would be required to properly characterize the flow field empirically. Important results included the observation that adjacent channels within the network had different water surface slopes for the same moments in time; thus calibration data for modelling similar systems needs to account for these significant local differences. Further, the in‐channel hydraulic roughness coefficient strongly varied from one cross‐section to the next (Manning's ‘n’ range: 0·01 to 0·10). These differences were amplified during low flow but persisted in muted form during high discharges. Copyright © 2012 John Wiley & Sons, Ltd.     

Supraglacial melt channel networks in the Jakobshavn Isbræ region during the 2007 melt season

Supraglacial channels are an important mechanism for surface water transport over the ablation zone of western Greenland. The first assessment of the spatio‐temporal distribution of surface melt channels and their relationship to supraglacial lakes over the Jakobshavn Isbræ region of Western Greenland was analysed using Landsat Enhanced Thematic Mapper Plus panchromatic images during the 2007 melt season. A total of 1188 melt channels were delineated and show an increase in the number of melt channels throughout the season, reaching a peak on 9 August. Water‐filled melt channels advanced to a maximum elevation of 1647 m on 9 August and attained a minimum average slope of 0.009 on 8 July. The ablation zone demonstrates two hydrologic modes, where crevasse and moulin terminating channels dominate at elevations <800 m and higher‐order channel networks >800 m. Development of higher‐order networks is interrupted by flow divergence due to partitioning of melt water into vertical infiltration through moulins and crevasse fields prevalent at lower elevations. Tributary and connector networks between 800 and 1200 m in elevation are correlated with fewer lake occurrences, lower surface velocities (~50 m a^?1), and ice flow dominated by internal deformation over basal sliding. High‐order channels are associated with lake basins that exceed melt water storage capacity. Evolution of channel networks is coupled to changes in melt water production, runoff, and ice dynamics with implication for the englacial and subglacial environments. © 2013 The Authors. Hydrological Processes Published by John Wiley & Sons, Ltd.     

Near‐surface hydraulic conductivity of northern hemisphere glaciers

The hydrology of near‐surface glacier ice remains a neglected aspect of glacier hydrology despite its role in modulating meltwater delivery to downstream environments. To elucidate the hydrological characteristics of this near‐surface glacial weathering crust, we describe the design and operation of a capacitance‐based piezometer that enables rapid, economical deployment across multiple sites and provides an accurate, high‐resolution record of near‐surface water‐level fluctuations. Piezometers were employed at 10 northern hemisphere glaciers, and through the application of standard bail–recharge techniques, we derive hydraulic conductivity (K) values from 0.003 to 3.519 m day^?1, with a mean of 0.185 ± 0.019 m day^?1. These results are comparable to those obtained in other discrete studies of glacier near‐surface ice, and for firn, and indicate that the weathering crust represents a hydrologically inefficient aquifer. Hydraulic conductivity correlated positively with water table height but negatively with altitude and cumulative short‐wave radiation since the last synoptic period of either negative air temperatures or turbulent energy flux dominance. The large range of K observed suggests complex interactions between meteorological influences and differences arising from variability in ice structure and crystallography. Our data demonstrate a greater complexity of near‐surface ice hydrology than hitherto appreciated and support the notion that the weathering crust can regulate the supraglacial discharge response to melt production. The conductivities reported here, coupled with typical supraglacial channel spacing, suggest that meltwater can be retained within the weathering crust for at least several days. Not only does this have implications for the accuracy of predictive meltwater run‐off models, but we also argue for biogeochemical processes and transfers that are strongly conditioned by water residence time and the efficacy of the cascade of sediments, impurities, microbes, and nutrients to downstream ecosystems. Because continued atmospheric warming will incur rising snowline elevations and glacier thinning, the supraglacial hydrological system may assume greater importance in many mountainous regions, and consequently, detailing weathering crust hydraulics represents a research priority because the flow path it represents remains poorly constrained.     

Impacts of sediment load on Manning coefficient in supercritical shallow flow on steep slopes

The Manning equation is one of the most widely used formulae for calculating the velocity of shallow overland flow in hydrological and erosion models. Precise estimation of the Manning's friction coefficient (n) is critical to determining overland flow and soil erosion processes. Few studies have been conducted to quantify the effects of sediment load on Manning's n on steep slopes. This study was conducted to investigate the potential effects of sediment load on Manning's n in a flume with a fixed bed, under wide ranges of hydraulics and sediment loads. Slope gradient varied from 8·7 to 34·2%, unit flow rate from 0·66 to 5·26 × 10^?3 m² s^?1, and sediment load from 0 to 6·95 kg m^?1 s^?1. The Reynolds number ranged from 350 to 5899. Results showed that Manning's n varied in both sediment‐free and sediment‐laden flows ranging from 0·012 to 0·055. The apparent Manning's coefficients of sediment‐laden flow were much greater than those of sediment‐free flow. The mean Manning coefficient of sediment‐laden flow was 51·27% greater than the mean value of sediment‐free flow. For sediment‐laden flow, Manning's n could be estimated with a power function of unit flow discharge and sediment content. Further studies are needed to quantify the potential effects of sediment load on the Manning's n on erodible beds and in fields. Copyright © 2010 John Wiley & Sons, Ltd.     

Improvement of surface run‐off in the hydrological model ParFlow by a scale‐consistent river parameterization

We propose an improvement of the overland‐flow parameterization in a distributed hydrological model, which uses a constant horizontal grid resolution and employs the kinematic wave approximation for both hillslope and river channel flow. The standard parameterization lacks any channel flow characteristics for rivers, which results in reduced river flow velocities for streams narrower than the horizontal grid resolution. Moreover, the surface areas, through which these wider model rivers may exchange water with the subsurface, are larger than the real river channels potentially leading to unrealistic vertical flows. We propose an approximation of the subscale channel flow by scaling Manning's roughness in the kinematic wave formulation via a relationship between river width and grid cell size, following a simplified version of the Barré de Saint‐Venant equations (Manning–Strickler equations). The too large exchange areas between model rivers and the subsurface are compensated by a grid resolution‐dependent scaling of the infiltration/exfiltration rate across river beds. We test both scaling approaches in the integrated hydrological model ParFlow. An empirical relation is used for estimating the true river width from the mean annual discharge. Our simulations show that the scaling of the roughness coefficient and the hydraulic conductivity effectively corrects overland flow velocities calculated on the coarse grid leading to a better representation of flood waves in the river channels.     

The topographic data source of digital terrain models as a key element in the accuracy of hydraulic flood modelling

The effects of the topographic data source and resolution on the hydraulic modelling of floods were analysed. Seven digital terrain models (DTMs) were generated from three different altimetric sources: a global positioning system (GPS) survey and bathymetry; high‐resolution laser altimetry data LiDAR (light detection and ranging); and vectorial cartography (1:5000). Hydraulic results were obtained, using the HEC‐RAS one‐dimensional model, for all seven DTMs. The importance of the DTM's accuracy on the hydraulic modelling results was analysed within three different hydraulic contexts: (1) the discharge and water surface elevation results from the hydraulic model; (2) the delineation of the flooded area; and (3) the relative sensitivity of the hydraulic model to changes in the Manning's n roughness coefficient. The contour‐based DTM was the least accurate with a root mean square error (RMSE) of 4·5 m in the determination of the water level and a variation of up to 50 per cent in the estimation of the inundated area of the floodplain. The GPS‐based DTM produced more realistic water surface elevation results and variations of up to 8 per cent in terms of the flooded area. The laser‐based model's RMSE for water level was 0·3 m, with the flooded area varying by less than 1 per cent. The LiDAR data also showed the greatest sensitivity to changes in the Manning's roughness coefficient. An analysis of the effect of mesh resolution indicated an influence on the delineation of the flooded area with variations of up to 7·3 per cent. In addition to determining the accuracy of the hydraulic modelling results produced from each DTM, an analysis of the time–cost ratio of each topographic data source illustrates that airborne laser scanning is a cost‐effective means of developing a DTM of sufficient accuracy, especially over large areas. Copyright © 2006 John Wiley & Sons, Ltd.     

Hydraulic properties affected by litter and stem cover under overland flow

Vegetation stems and litter cover have different effects on sediment transport capacity under the same experimental conditions, which in essence, may be due to differences in their hydraulic properties, but the availability of comparative studies is limited. This study aimed to compare the hydraulic properties affected by litter and stem cover, compare differences in the drag forces exerted by litter and stems on overland flow, and develop new Manning's n and flow velocity equations for litter cover. Two series of flume experiments were conducted with the same slope gradients (8.8%, 17.6%, 26.8%) and flow discharge rates (0.5, 1.0 × 10⁻³ m³ s⁻¹). Artificial Gramineae stems with a 0%–30% cover level and Pinus tabulaeformis litter with a 0%–70% cover level were used in series 1 and series 2, respectively. The flow velocity and depth were measured. The results showed that the Froude number and flow velocity affected by stem cover were much lower than those affected by litter cover, while the opposite trend was observed in the relative magnitude of the Reynolds number, flow depth and shear stress. The form resistance caused by stems was 22–57 times greater than that caused by litter for the same cover level, which suggests that stem cover contributes more than litter cover to increasing the flow resistance and reducing the flow's ability for sediment detachment and transport. Two new equations for calculating Manning's n and flow velocity under the influence of litter cover were developed, with R² and NSE values of 0.96. The results of this study contribute to revealing the mechanisms of the differences of the effects of stem and litter cover on soil erosion.     

Flow resistance in steep mountain streams

Resistance to flow at low to moderate stream discharge was examined in five small (12–77 km² drainage area) tributaries of Chilliwack River, British Columbia, more than half of which exhibit planar bed morphology. The resulting data set is composed of eight to 12 individual estimates of the total resistance to flow at 61 cross sections located in 13 separate reaches of five tributaries to the main river. This new data set includes 625 individual estimates of resistance to flow at low to moderate river stage. Resistance to flow in these conditions is high, highly variable and strongly dependent on stage. The Darcy–Weisbach resistance factor (ff) varies over six orders of magnitude (0·29–12 700) and Manning's n varies over three orders of magnitude (0·047–7·95). Despite this extreme range, both power equations at the individual cross sections and Keulegan equations for reach‐averaged values describe the hydraulic relations well. Roughness is divided into grain and form (considered as all non‐grain sources) components. Form roughness is the dominant component, accounting for about 90% of the total roughness of the system (i.e., form roughness is on average 8.6 times as great as grain roughness). Of the various quantitative and qualitative form‐roughness indicators observed, only the sorting coefficient (σ = D₈₄/D₅₀) correlates well with form roughness. Copyright © 2008 John Wiley & Sons, Ltd.     

Assessment of vegetation effects on hydraulics and of feedbacks on plant survival and zonation in ephemeral channels

The interaction of vegetation and flow in channels is important for understanding the influences of forces in channels and effects on erosion, sediment flux and deposition; it has implications for channel habitats, channel instability and restoration schemes. Methods are needed for calculating forces on plants and data are required on thresholds for plant destruction and survival. A simple method of calculating the effect of hydraulics on vegetation and its zonation within ephemeral channels is described. Detailed cross section surveys of channel morphology, vegetation and estimates of Manning's n are input into the software program WinXSPRO to calculate the hydraulics of flows across the channel for a given event or flow level, incorporating subdivision into zones of differing morphology and vegetation across the section. This was applied to a number of cross sections on ephemeral channels in SE Spain and typical roughness values for Mediterranean vegetation types in channels were assessed. The method is demonstrated with reference to two well‐documented floods in SE Spain, in September 1997 on the Torrealvilla and in October 2003 along the Salada. These flows led to the mortality of herbs, reed and smaller shrub species. Some damage to larger shrubs and trees occurred, but trees such as Tamarisk (Tamarix canariensis) were shown to withstand high forces. Some grasses were highly resistant to removal and induced sedimentation. Significant erosion was limited to areas with little vegetation covering the channel floor. Further quantification of resistance of vegetation to flows and upper threshold values for removal is continuing by relating calculated hydraulic conditions using the methods outlined to measurements of vegetation responses in events at monitoring sites. Copyright © 2009 John Wiley & Sons, Ltd.     

The formation of supraglacial debris covers by primary dispersal from transverse englacial debris bands

Glaciological controls on debris cover formation are investigated from the perspective of primary dispersal of supraglacial debris across a melting ice surface. This involves the migration of angled debris septa outcrops across a melting, thinning glacier ablation zone. Three measures of a glacier's ability to evacuate supraglacial debris are outlined: (1) a concentration factor describing the focusing of englacial debris into specific supraglacial mass loads; (2) the rate of migration of a septum outcrop relative to the local ice surface; and (3) a downstream velocity differential between a slower septum outcrop and the faster ice surface velocity. Measures (1) and (2) are inversely related, while measure (3) increases down‐glacier to explain why slow‐moving, thinning ice rapidly becomes debris covered. Data from Glacier d'Estelette (Italian Alps) are used to illustrate these processes, and to explore the potential for debris cover formation and growth in different glaciological environments. The transition from a ‘clean’, transport‐dominated to a debris‐covered ablation‐dominated glacier is explained by the melting out of more closely‐spaced debris septa, in combination with the geometric interactions of angled septa and ice surface in a field of reducing flow and increasing ablation. The growth and shrinkage of debris covers are most sensitive to glaciological changes at glaciers with gently‐dipping debris‐bearing foliation, but less sensitive at high‐compression glaciers whose termini are constrained by moraine dams and other forms of obstruction. These findings show that a variety of debris‐covered glacier types will show a spectrum of response characteristics to negative mass balance. Copyright © 2013 John Wiley & Sons, Ltd.     

Measuring the flow resistance of submerged grass

Through laboratory experiments conducted in a grass‐lined flume, the hydraulic resistance of grass is measured and quantified. For the grass examined, it is found that Manning's n value is greater than those recommended in well‐established texts such as Chow (1959. Open Channel Flow. McGraw‐Hill: Singapore), relatively lower than those predicted by n‐UR methods, but corresponds well with the value found from calibration studies of two‐ and three‐dimensional numerical models. The assumption of a uniform Manning's n value with flow depth, which is often made in numerical modelling, may be invalid depending on the relative submergence of the vegetation. Drag coefficients are evaluated for a method applicable to three‐dimensional numerical models. Further detailed experimental investigation and application of these approaches within a numerical modelling framework is now recommended. Copyright © 2002 John Wiley & Sons, Ltd.     

VELOCITY–DISCHARGE RELATIONSHIPS DERIVED FROM DYE TRACER EXPERIMENTS IN GLACIAL MELTWATERS: IMPLICATIONS FOR SUBGLACIAL FLOW CONDITIONS

Repeated dye tracer tests were undertaken from individual moulins at Haut Glacier d'Arolla, Switzerland, over a number of diurnal discharge cycles during the summers of 1989–1991. It was hoped to use the concepts of at-a-station hydraulic geometry to infer flow conditions in subglacial channels from the form of the velocity–discharge relationships derived from these tests. The results, however, displayed both clockwise and anticlockwise velocity–discharge hysteresis, in addition to the simple power function relationship assumed in the hydraulic geometry approach. Clockwise hysteresis seems to indicate that a moulin drains into a small tributary channel rather than directly into an arterial channel, and that discharges in the two channels vary out of phase with each other. Anticlockwise hysteresis is accompanied by strong diurnal variations in the value of dispersivity derived from the dye breakthrough curve, and is best explained by hydraulic damming of moulins or sub/englacial passageways. Despite the complex velocity–discharge relationships observed, some indication of subglacial flow conditions may be obtained if tributary channels comprise only a small fraction of the drainage path and power function velocity–discharge relationships are derived from dye injections conducted during periods when the supraglacial discharge entering the moulin and the bulk discharge vary in phase. Analyses based on this premise suggest that both open and closed channel flow occur beneath Haut Glacier d'Arolla, and that flow conditions are highly variable at and between sites.     

Comparing the performance of a 2‐D finite element and a 2‐D finite volume model of floodplain inundation using airborne SAR imagery

The performances of a finite volume model (SFV) and finite element model (TELEMAC‐2D) in reproducing inundation on a 16 km reach of the river Severn, United Kingdom, are compared. Predicted inundation extents are compared with 4 airborne synthetic aperture radar images of a major flood event in November 2000, and these are used to calibrate 2 values of Manning's n for the channel and floodplain. The four images are shown to have different capacities to constrain roughness parameters, with the image acquired at low flow rate doing better in determining these parameters than the image acquired at approximately peak flow. This is assigned to the valley filling nature of the flood and the associated insensitivity of flood extent to changes in water level. The level of skill demonstrated by the models, when compared with inundation derived using a horizontal water free surface, also increases as flow rate drops. The two models show markedly different behaviours to the calibration process, with TELEMAC showing less sensitivity and lower optimum values for Manning's n than SFV. When the models are used in predictive mode, calibrated against one image and predicting another, SFV performs better than TELEMAC. Copyright © 2007 John Wiley & Sons, Ltd.     

A method for estimating spatially variable seepage and hydraulic conductivity in channels with very mild slopes

Infiltration along ephemeral channels plays an important role in groundwater recharge in arid regions. A model is presented for estimating spatial variability of seepage due to streambed heterogeneity along channels based on measurements of streamflow‐front velocities in initially dry channels. The diffusion‐wave approximation to the Saint‐Venant equations, coupled with Philip's equation for infiltration, is connected to the groundwater model MODFLOW and is calibrated by adjusting the saturated hydraulic conductivity of the channel bed. The model is applied to portions of two large water delivery canals, which serve as proxies for natural ephemeral streams. Estimated seepage rates compare well with previously published values. Possible sources of error stem from uncertainty in Manning's roughness coefficients, soil hydraulic properties and channel geometry. Model performance would be most improved through more frequent longitudinal estimates of channel geometry and thalweg elevation, and with measurements of stream stage over time to constrain wave timing and shape. This model is a potentially valuable tool for estimating spatial variability in longitudinal seepage along intermittent and ephemeral channels over a wide range of bed slopes and the influence of seepage rates on groundwater levels. Copyright © 2012 John Wiley & Sons, Ltd.     

The impact of flow discharge on the hydraulic characteristics of headcut erosion processes in the gully region of the Loess Plateau

Headcut erosion is associated with major hydraulic changes induced by the gully head of concentrated flow. However, the variation in the hydraulic characteristics of the headcut erosion process is still not clear in the gully region of the Loess Plateau. A series of rainfall combined scouring experiments (flow discharges ranging from 3.6 to 7.2 m³ hr⁻¹, with 0.8 mm min⁻¹ rainfall intensity) were conducted on experimental plots to clarify the variation in the hydraulic parameters induced by gully head and erosion processes under different flow discharges. The results showed that concentrated flows in the catchment area and gully bed were turbulent (Reynolds number ranging from 1,876 to 6,693) and transformed between supercritical and subcritical (Froude number ranging from 0.96 to 3.73). The hydraulic parameters, such as the flow velocity, Reynolds number, shear stress, stream power, Darcy–Weisbach friction factor, and unit stream power in the catchment area were 0.45–0.59 m s⁻¹, 2086–6693, 1.96–5.33 Pa, 0.89–2.86 W m⁻², 0.08–0.16, and 0.023–0.031 m s⁻¹, respectively. When the concentrated flows dropped from the gully head, the hydraulic parameters in the gully bed decreased by 3.39–26.07%, 1.49–29.99%, 65.19–67.14%, 67.25–74.96%, 28.53–61.31%, and 67.82–77.14%, respectively, which contributed to the flow energy consumption at the gully head. As flow discharge increased, Reynolds number, shear stress, and stream power increased, while flow velocity, Froude number, unit stream power, and Darcy–Weisbach friction factor did not. The flow energy consumption at the gully head was 9.66–10.13, 13.25–13.74, 15.68–16.41, and 19.28–20.25 J s⁻¹, respectively, under different flow discharges and accounted for 60.58–68.50% of the flow energy consumption of the experimental plots. Generally, the sediment discharges increased rapidly at the initial stage, then increased slowly, and finally reached a steady state condition, which showed a significant declining logarithmic trend with experimental duration (P<.01) and increased with increasing flow discharge. Accordingly, the flow energy consumption was significantly correlated with the sediment yield. These findings could improve our understanding of the hydraulic properties and flow energy characteristics of headcut erosion.     

The frictional resistance characteristics of a bedrock‐influenced river channel

Much of the published literature relating to reach‐scale total ?ow resistance concentrates on alluvial rivers and little mention is made of the more extreme resistance found in bedrock‐in?uenced river systems. This paper presents the results of a detailed investigation of total ?ow resistance variation for ?ve channel types on the Sabie River, South Africa (cohesive mixed anastomosed, uncohesive mixed anastomosed, mixed pool–rapid, alluvial braided and alluvial single thread) over a range of discharges. Manning's ‘n’ and Darcy‐Weisbach f resistance coef?cients are calculated using Barnes' methodology (H. H. Barnes, US Geological Survey Water Supply Paper 1849, 1967). These data are analysed at the scales of channel type and morphologic unit and the resultant resistance values are compared with the literature. Generally the low‐?ow friction values estimated for the bedrock‐in?uenced channel types are higher than any reported. Flow resistance quanti?cation of the alluvial channel type study sites on the Sabie River produced values slightly higher than earlier studies of similar morphologies. In situations where signi?cant areas of irregular bedrock are exposed (cohesive mixed anastomosed and mixed pool–rapid channel types) the high ?ow resistance coef?cients at low discharges probably re?ect energy dissipation due to hydraulic jumps and internal distortions; however, these values fall as discharge increases to values comparable to those recorded for step–pool rivers. It is suggested that the slightly elevated resistance coef?cients recorded in the alluvial channel types are the result of the vegetation roughness component, operating when the macro‐channel bar features are inundated. Copyright © 2004 John Wiley & Sons, Ltd.     

A kinematic‐wave‐based distributed watershed model using FEM,GIS and remotely sensed data

For the appropriate management of water resources in a watershed, it is essential to calculate the time distribution of runoff for the given rainfall event. In this paper, a kinematic‐wave‐based distributed watershed model using finite element method (FEM), geographical information systems (GIS) and remote‐sensing‐based approach is presented for the runoff simulation of small watersheds. The kinematic wave equations are solved using FEM for overland and channel flow to generate runoff at the outlet of the watershed concerned. The interception loss is calculated by an empirical model based on leaf area index (LAI). The Green‐Ampt Mein Larson (GAML) model is used for the estimation of infiltration. Remotely sensed data has been used to extract land use (LU)/land cover (LC). GIS have been used to prepare finite element grid and input files such as Manning's roughness and slope. The developed overland flow model has been checked with an analytical solution for a hypothetical watershed. The model has been applied to a gauged watershed and an ungauged watershed. From the results, it is seen that the model is able to simulate the hydrographs reasonably well. A sensitivity analysis of the model is carried out with the calibrated infiltration parameters, overland flow Manning's roughness, channel flow Manning's roughness, time step and grid size. The present model is useful in predicting the hydrograph in small, ungauged watersheds. Copyright © 2007 John Wiley & Sons, Ltd.     

Sensitivity analysis of hydraulic conductivity and Manning's n parameters lead to new method to scale effective hydraulic conductivity across model resolutions

Hydrological modelling is an important tool for research, policy, and management, but uncertainty remains about parameters transferability from field observations made at small scale to models at the catchment scale and larger. This uncertainty compels the need to develop parameter relationships that are translatable across scale. In this study, we compare the changes to modelled processes as resolution is coarsened from 100‐m to 1‐km in a topographically complex, 255‐km² Colorado River headwater catchment. We conducted a sensitivity analysis for hydraulic conductivity (K) and Manning's n parameters across four orders of magnitude. Results showed that K acts as a moderator between surface and subsurface contributions to streamflow, whereas n moderates the duration of high intensity, infiltration‐excess flow. The parametric sensitivity analysis informed development of a new method to scale effective hydraulic conductivity across modelling resolutions in order to compensate for the loss of topographic gradients as resolution is coarsened. A similar mathematical relationship between n and lateral resolution changes was not found, possibly because n is also sensitive to time discretization. This research provides an approach to translate hydraulic conductivity parameters from a calibrated coarse model to higher resolutions where the number of simulations are limited by computational demand.