Sub‐grid scale parameterization of hillslope runoff and erosion processes for catchment‐scale models of semi‐arid landscapes

The processes of hillslope runoff and erosion are typically represented at coarse spatial resolution in catchment‐scale models due to computational limitations. Such representation typically fails to incorporate the important effects of topographic heterogeneity on runoff generation, overland flow, and soil erosion. These limitations currently undermine the application of distributed catchment models to understand the importance of thresholds and connectivity on hillslope and catchment‐scale runoff and erosion, particularly in semi‐arid environments. This paper presents a method for incorporating high‐resolution topographic data to improve sub‐grid scale parameterization of hillslope overland flow and erosion models. Results derived from simulations conducted using a kinematic wave overland flow model at 0.5 m spatial resolution are used to parameterize the depth–discharge relationship in the overland flow model when applied at 16 m resolution. The high‐resolution simulations are also used to derive a more realistic parameterization of excess flow shear stress for use in the 16 m resolution erosion model. Incorporating the sub‐grid scale parameterization in the coarse‐resolution model (16 m) leads to improved predictions of overland flow and erosion when evaluated using results derived from high‐resolution (0.5 m) model simulations. The improvement in performance is observed for a range of event magnitudes and is most notable for erosion estimates due to the non‐linear dependency between the rates of erosion and overland flow. Copyright © 2013 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.     

A distributed TOPMODEL for modelling impacts of land‐cover change on river flow in upland peatland catchments

There is global concern about headwater management and associated impacts on river flow. In many wet temperate zones peatlands can be found covering headwater catchments. In the UK there is major concern about how environmental change, driven by human interventions, has altered the surface cover of headwater blanket peatlands. However, the impact of such land‐cover changes on river flow is poorly understood. In particular, there is poor understanding of the impacts of different spatial configurations of bare peat or well‐vegetated, restored peat on river flow peaks in upland catchments. In this paper, a physically based, distributed and continuous catchment hydrological model was developed to explore such impacts. The original TOPMODEL, with its process representation being suitable for blanket peat catchments, was utilized as a prototype acting as the basis for the new model. The equations were downscaled from the catchment level to the cell level. The runoff produced by each cell is divided into subsurface flow and saturation‐excess overland flow before an overland flow calculation takes place. A new overland flow module with a set of detailed stochastic algorithms representing overland flow routing and re‐infiltration mechanisms was created to simulate saturation‐excess overland flow movement. The new model was tested in the Trout Beck catchment of the North Pennines of England and found to work well in this catchment. The influence of land cover on surface roughness could be explicitly represented in the model and the model was found to be sensitive to land cover. Copyright © 2014 John Wiley & Sons, Ltd.     

A physically based distributed subsurface–surface flow dynamics model for forested mountainous catchments

This study was designed to develop a physically based hydrological model to describe the hydrological processes within forested mountainous river basins. The model describes the relationships between hydrological fluxes and catchment characteristics that are influenced by topography and land cover. Hydrological processes representative of temperate basins in steep terrain that are incorporated in the model include intercepted rainfall, evaporation, transpiration, infiltration into macropores, partitioning between preferential flow and soil matrix flow, percolation, capillary rise, surface flow (saturation‐excess and return flow), subsurface flow (preferential subsurface flow and baseflow) and spatial water‐table dynamics. The soil–vegetation–atmosphere transfer scheme used was the single‐layer Penman–Monteith model, although a two‐layer model was also provided. The catchment characteristics include topography (elevation, topographic indices), slope and contributing area, where a digital elevation model provided flow direction on the steepest gradient flow path. The hydrological fluxes and catchment characteristics are modelled based on the variable source‐area concept, which defines the dynamics of the watershed response. Flow generated on land for each sub‐basin is routed to the river channel by a kinematic wave model. In the river channel, the combined flows from sub‐basins are routed by the Muskingum–Cunge model to the river outlet; these comprise inputs to the river downstream. The model was applied to the Hikimi river basin in Japan. Spatial decadal values of the normalized difference vegetation index and leaf area index were used for the yearly simulations. Results were satisfactory, as indicated by model efficiency criteria, and analysis showed that the rainfall input is not representative of the orographic lifting induced rainfall in the mountainous Hikimi river basin. Also, a simple representation of the effects of preferential flow within the soil matrix flow has a slight significance for soil moisture status, but is insignificant for river flow estimations. Copyright © 2005 John Wiley & Sons, Ltd.     

Patch‐scale representation of vegetation within hydraulic models

Submerged aquatic vegetation affects flow, sediment and ecological processes within rivers. Quantifying these effects is key to effective river management. Despite a wealth of research into vegetated flows, the detailed flow characteristics around real plants in natural channels are still poorly understood. Here we present a new methodology for representing vegetation patches within computational fluid dynamics (CFD) models of vegetated channels. Vegetation is represented using a Mass Flux Scaling Algorithm (MFSA) and drag term within the Reynolds‐averaged Navier–Stokes Equations, which account for the mass and momentum effects of the vegetation, respectively. The model is applied using three different grid resolutions (0.2, 0.1 and 0.05 m) using time‐averaged solution methods and compared to field data. The results show that the model reproduces the complex spatial flow heterogeneity within the channel and that increasing the resolution leads to enhanced model accuracy. Future applications of the model to the prediction of channel roughness, sedimentation and key eco‐hydraulic variables are presented, likely to be valuable for informing effective river management. © 2016 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Back‐barrier flooding by storm surges and overland flow

Low‐elevation areas within a sandy barrier island are subject to flooding via saturation overland flow following moderate storm surges and rainfall events. Using a high resolution topographic survey and simple hydrology models, we estimate the discharge and velocities from storm surge return flow and saturation overland flow. Results show that return flow velocities are of the same magnitude as the critical velocity necessary to mobilize sand when a hydraulic connection between the watershed and back‐barrier bay is present. Storms of moderate strength and rainfall intensity may be sufficient to keep the return channels open within the back‐barrier, thus providing natural conduits for water exchange from overwash events during extreme storm surges triggered by hurricanes. Copyright © 2011 John Wiley & Sons, Ltd.     

Development and application of a spatially‐distributed Arctic hydrological and thermal process model (ARHYTHM)

A process‐based, spatially distributed hydrological model was developed to quantitatively simulate the energy and mass transfer processes and their interactions within arctic regions (arctic hydrological and thermal model, ARHYTHM). The model first determines the flow direction in each element, the channel drainage network and the drainage area based upon the digital elevation data. Then it simulates various physical processes: including snow ablation, subsurface flow, overland flow and channel flow routing, soil thawing and evapotranspiration. The kinematic wave method is used for conducting overland flow and channel flow routing. The subsurface flow is simulated using the Darcian approach. The energy balance scheme was the primary approach used in energy‐related process simulations (snowmelt and evapotranspiration), although there are options to model snowmelt by the degree‐day method and evapotranspiration by the Priestley–Taylor equation. This hydrological model simulates the dynamic interactions of each of these processes and can predict spatially distributed snowmelt, soil moisture and evapotranspiration over a watershed at each time step as well as discharge in any specified channel(s). The model was applied to Imnavait watershed (about 2·2 km²) and the Upper Kuparuk River basin (about 146 km²) in northern Alaska. Simulated results of spatially distributed soil moisture content, discharge at gauging stations, snowpack ablations curves and other results yield reasonable agreement, both spatially and temporally, with available data sets such as SAR imagery‐generated soil moisture data and field measurements of snowpack ablation, and discharge data at selected points. The initial timing of simulated discharge does not compare well with the measured data during snowmelt periods mainly because the effect of snow damming on runoff was not considered in the model. Results from the application of this model demonstrate that spatially distributed models have the potential for improving our understanding of hydrology for certain settings. Finally, a critical component that led to the performance of this modelling is the coupling of the mass and energy processes. Copyright © 2000 John Wiley & Sons, Ltd.     

Computational fluid dynamics modelling of boundary roughness in gravel‐bed rivers: an investigation of the effects of random variability in bed elevation

Results from a series of numerical simulations of two‐dimensional open‐channel flow, conducted using the computational fluid dynamics (CFD) code FLUENT, are compared with data quantifying the mean and turbulent characteristics of open‐channel flow over two contrasting gravel beds. Boundary roughness effects are represented using both the conventional wall function approach and a random elevation model that simulates the effects of supra‐grid‐scale roughness elements (e.g. particle clusters and small bedforms). Results obtained using the random elevation model are characterized by a peak in turbulent kinetic energy located well above the bed (typically at y/h = 0·1–0·3). This is consistent with the field data and in contrast to the results obtained using the wall function approach for which maximum turbulent kinetic energy levels occur at the bed. Use of the random elevation model to represent supra‐grid‐scale roughness also allows a reduction in the height of the near‐bed mesh cell and therefore offers some potential to overcome problems experienced by the wall function approach in flows characterized by high relative roughness. Despite these benefits, the results of simulations conducted using the random elevation model are sensitive to the horizontal and vertical mesh resolution. Increasing the horizontal mesh resolution results in an increase in the near‐bed velocity gradient and turbulent kinetic energy, effectively roughening the bed. Varying the vertical resolution of the mesh has little effect on simulated mean velocity profiles, but results in substantial changes to the shape of the turbulent kinetic energy profile. These findings have significant implications for the application of CFD within natural gravel‐bed channels, particularly with regard to issues of topographic data collection, roughness parameterization and the derivation of mesh‐independent solutions. Copyright © 2001 John Wiley & Sons, Ltd.     

Scale invariance and self‐similarity in kinematic wave overland flow in space and time

Fractals are famous for their self‐similar nature at different spatial scales. Similar to fractals, solutions of scale invariant processes are self‐similar at different space–time scales. This unique property of scale‐invariant processes can be utilized to translate the solution of the processes at a much larger or smaller space–time scale (domain) based on the solution calculated on the original space–time scale. This study investigates scale invariance conditions of kinematic wave overland flow process in one‐parameter Lie group of point transformations framework. Scaling (stretching) transformation is one of the one‐parameter Lie group of point transformations and it has a unique importance among the other transformations, as it leads to the scale invariance or scale dependence of a process. Scale invariance of a process yields a self‐similar solution at different space–time scales. However, the conditions for the process to be scale invariant usually dictate various relationships between the scaling coefficients of the dependent and independent variables of the process. Therefore, the scale invariance of a process does not assure a self‐similar solution at any arbitrary space and time scale. The kinematic wave overland flow process is modelled mathematically as initial‐boundary value problem. The conditions to be satisfied by the system of governing equations as well as the initial and boundary conditions of the kinematic wave overland flow process are established in order for the process to be scale invariant. Also, self‐similarity of the solution of the kinematic wave overland flow under the established invariance conditions is demonstrated by various numerical example problems. Copyright © 2011 John Wiley & Sons, Ltd.     

Field investigation and modelling of coupled stream discharge and shallow water‐table dynamics in a small Mediterranean catchment (Sardinia)

In the semi‐arid Mediterranean environment, the rainfall–runoff relationships are complex because of the markedly irregular patterns in rainfall, the seasonal mismatch between evaporation and rainfall, and the spatial heterogeneity in landscape properties. Watersheds often display considerable non‐linear threshold behavior, which still make runoff generation an open research question. Our objectives in this context were: to identify the primary processes of runoff generation in a small natural catchment; to test whether a physically based model, which takes into consideration only the primary processes, is able to predict spatially distributed water‐table and stream discharge dynamics; and to use the hydrological model to increase our understanding of runoff generation mechanisms. The observed seasonal dynamics of soil moisture, water‐table depth, and stream discharge indicated that Hortonian overland‐flow was negligible and the main mechanism of runoff generation was saturated subsurface‐flow. This gives rise to base‐flow, controls the formation of the saturated areas, and contributes to storm‐flow together with saturation overland‐flow. The distributed model, with a 1D scheme for the kinematic surface‐flow, a 2D sub‐horizontal scheme for the saturated subsurface‐flow, and ignoring the unsaturated flow, performed efficiently in years when runoff volume was high and medium, although there was a smoothing effect on the observed water‐table. In dry years, small errors greatly reduced the efficiency of the model. The hydrological model has allowed to relate the runoff generation mechanisms with the land‐use. The forested hillslopes, where the calibrated soil conductivity was high, were never saturated, except at the foot of the slopes, where exfiltration of saturated subsurface‐flow contributed to storm‐flow. Saturation overland‐flow was only found near the streams, except when there were storm‐flow peaks, when it also occurred on hillslopes used for pasture, where soil conductivity was low. The bedrock–soil percolation, simulated by a threshold mechanism, further increased the non‐linearity of the rainfall–runoff processes. Copyright © 2013 John Wiley & Sons, Ltd.     

Morphological signatures of normal faulting in low‐gradient alluvial rivers in south‐eastern Louisiana,USA

We explore the fluvial response to faulting in three low‐gradient, sand‐bed rivers in south‐eastern Louisiana, USA, that flow across active normal faults from footwall (upstream) to hangingwall (downstream). We calculate sinuosity, migration rate and migration direction in order to identify anomalies spatially associated with fault scarps. In two of the rivers we model one‐dimensional steady water flow to identify anomalies in surface water slope, width‐to‐depth ratio, and shear stress. In each of these rivers there is one location where flow modeling suggests potential channel incision through the footwall, as indicated by relatively high surface water slopes and shear stress values. In one of these footwall locations, the river straightens and width‐to‐depth ratios decrease, likely contributing to higher surface water slopes and shear stress. This is in contrast to previous studies that have proposed increased sinuosity across fault footwalls and decreased sinuosity across hangingwalls. However, in two hangingwall locations we also observe relatively less sinuous channels. Other planform changes on the hangingwall include topographic steering of channels along and towards the fault and one example of an avulsion. The most notable anomaly in migration rate occurs on the hangingwall of a fault where a river has cut off a meander loop. Although fluvial response to faulting varies here, comparatively large and small channels exhibit similar responses. Further, Pleistocene fault slip rates are orders of magnitude lower than the channel migration rates, suggesting that faulting should not be a major influence on the fluvial evolution. Nonetheless, notable channel anomalies exist near faults, suggesting that recent fault slip rates are higher than Pleistocene rates, and/or that low‐gradient alluvial channels are more sensitive to faulting than previous studies have suggested. Rivers appear to be influenced by faulting in this setting, however background rates of meander loop cutoff may be just as influential as faulting. Copyright © 2015 John Wiley & Sons, Ltd.     

Simulation of the direct runoff hydrograph at basin outlet

A simple model describing the transformation of effective rainfall to direct runoff through the overland flow mechanism is presented. The model is based on the classical representation of a watershed by a combination of planes and channels. The dynamics of overland flow in each plane is simulated by the non-linear kinematic wave, but the outflow from a given plane is concentrated in the middle of the corresponding drainage channel. The water routing in the channels is carried out by a piece-wise linearized formulation in space of the kinematic wave approximation. Using synthetic events on 10 watersheds, the model was tested by comparing it with results obtained by applying the non-linear kinematic wave to all the elements of the watershed. The model was found to be adequate, even in a form that simplifies the geometric features of the planes through an averaging procedure based on the Horton–Strahler ordering scheme of the watershed. © 1998 John Wiley & Sons, Ltd.     

Hydraulic modelling of fish habitat in urban rivers during high flows

In urban rivers, flow regime and channel morphology are the drivers of physical habitat quality for aquatic species. Peak discharges are increased at high flows as a result of impermeable catchments and channel engineering for flood protection schemes. Hazardous conditions and flashy hydrographs mean that measurement of velocities at high flows is a difficult task. This research uses a three‐dimensional computational fluid dynamics (3D‐CFD) model to simulate hydraulic patterns in two urban river channels. A 3D‐CFD code, called SSIIM, was used to simulate hydraulic conditions in two engineered river reaches of the River Tame, Birmingham, UK. These two sites represent channels with different levels of engineering. Models were calibrated and tested using field measurements. Results show that modelled water surface levels and velocity profiles are well simulated. Calibrated roughness heights are compared with those derived from field measurement of sediment size. Numerical experiments are used to assess the relationship between grid resolution in the vertical dimension and the form of the modelled velocity profiles. Biologists have used laboratory experiments to determine maximum sustainable swimming speeds (MSSS) of fish, often in order to assess what level of a particular pollutant may be tolerable. In this work, simulations of high‐flow hydraulic patterns are used to compare velocity patterns with fish MSSS. Results show that when the water levels rise to fill the first channel of the two‐stage channels at the sites, which occurred 16 times in 2000, MSSS are surpassed in the majority of available habitat, suggesting that excessive velocities at high flows are one factor that limits fish habitat. A comparison between the two reaches shows that there is less available habitat in the more modified reach. Conclusions suggest that an approach that integrates water quality issues and physical channel characteristics must be taken in river rehabilitation schemes, as improvements to water quality alone may not be sufficient to improve habitat quality to the desired level. Copyright © 2002 John Wiley & Sons, Ltd.     

Form roughness and the absence of secondary flow in a large confluence–diffluence,Rio Paraná, Argentina

Confluence–diffluence units are key elements within many river networks, having a major impact upon the routing of flow and sediment, and hence upon channel change. Although much progress has been made in understanding river confluences, and increasing attention is being paid to bifurcations and the important role of bifurcation asymmetry, most studies have been conducted in laboratory flumes or within small rivers with width:depth (aspect) ratios less than 50. This paper presents results of a field‐based study that details the bed morphology and 3D flow structure within a very large confluence–diffluence in the Río Paraná, Argentina, with a width:depth ratio of approximately 200. Flow within the confluence–diffluence is dominated largely by the bed roughness, in the form of sand dunes; coherent, channel‐scale, secondary flow cells, that have been identified as important aspects of the flow field within smaller channels, and assumed to be present within large rivers, are generally absent in this reach. This finding has profound implications for flow mixing rates, sediment transport rates and pathways, and thus the interpretation of confluence–diffluence morphology and sedimentology. Copyright © 2006 John Wiley & Sons, Ltd.     

Modelling long‐term contaminant migration in a catchment at fine spatial and temporal scales using the UP system

A method has been developed to simulate the long‐term migration of radionuclides in the near‐surface of a river catchment, following their release from a deep underground repository for radioactive waste. Previous (30‐year) simulations, conducted using the SHETRAN physically based modelling system, showed that long‐term (many decades) simulations are required to allow the system to reach steady state. Physically based, distributed models, such as SHETRAN, tend to be too computationally expensive for this task. Traditional lumped catchment‐scale models, on the other hand, do not give sufficiently detailed spatially distributed results. An intermediate approach to modelling has therefore been developed which allows flow and transport processes to be simulated with the spatial resolution normally associated with distributed models, whilst being computationally efficient.The approach involves constructing a lumped model in which the catchment is represented by a number of conceptual water storage compartments. The flow rates to and from these compartments are prescribed by functions that summarize the results from physically based distributed models run for a range of characteristic flow regimes. The physically based models used were, SHETRAN for the subsurface compartments, a particle tracking model for overland flow and an analytical model for channel routing. One important advantage of the method used in constructing the lumped model is that it makes down scaling possible, in the sense that fine‐scale information on the distributed hydrological regime, as simulated by the physically based distributed models, can be inferred from the variables in the lumped model that describe the hydrology at the catchment scale. A 250‐year flow simulation has been run and the down scaling process used to infer a 250‐year time‐series of three‐dimensional velocity fields for the subsurface of the catchment. This series was then used to drive a particle tracking simulation of contaminant migration. The concentration and spatial distribution of contaminants simulated by this model for the first 30 years were in close agreement with SHETRAN results. The remaining 220 years highlighted the fact that some of the most important transport pathways to the surface carry contaminants only very slowly so both the magnitude and spatial distribution of concentration in surface soils are not apparent over the shorter SHETRAN simulations. Copyright © 1999 John Wiley & Sons, Ltd.     

A hydrochemical approach to quantify the role of return flow in a surface flow‐dominated catchment

Stormflow generation in headwater catchments dominated by subsurface flow has been studied extensively, yet catchments dominated by surface flow have received less attention. We addressed this by testing whether stormflow chemistry is controlled by either (a) the event‐water signature of overland flow, or (b) the pre‐event water signature of return flow. We used a high‐resolution hydrochemical data set of stormflow and end‐members of multiple storms in an end‐member mixing analysis to determine the number of end‐members needed to explain stormflow, characterize and identify potential end‐members, calculate their contributions to stormflow, and develop a conceptual model of stormflow. The arrangement and relative positioning of end‐members in stormflow mixing space suggest that saturation excess overland flow (26–48%) and return flow from two different subsurface storage pools (17–53%) are both similarly important for stormflow. These results suggest that pipes and fractures are important flow paths to rapidly release stored water and highlight the value of within‐event resolution hydrochemical data to assess the full range and dynamics of flow paths.     

Numerical simulation of three‐dimensional flow hydraulics in a braided channel

Results are presented from a numerical simulation of three‐dimensional flow hydraulics around a mid‐channel bar carried out using the FLUENT/UNS computational fluid dynamics (CFD) software package. FLUENT/UNS solves the three‐dimensional Reynolds‐averaged form of the Navier–Stokes equations. Turbulence closure is achieved using a RNG k–ϵ model. Simulated flow velocities are compared with measured two‐dimensional velocities (downstream and cross‐stream) obtained using an electromagnetic current meter (ECM). The results of the simulation are qualitatively consistent with the flow structures observed in the field. Quantitative comparison of the simulated and measured velocity magnitudes indicates a strong positive correlation between the two (r=0·88) and a mean difference of 0·09 m s⁻¹. Deviations between simulated and measured velocities may be identified that are both random and systematic. The former may reflect a number of factors including subgrid‐scale natural spatial variability in flow velocities associated with local bed structures and measurement uncertainty resulting from problems of ECM orientation. Model mesh configuration, roughness parameterization and inlet boundary condition uncertainty may each contribute to systematic differences between simulated and measured flow velocities. These results illustrate the potential for using CFD software to simulate flow hydraulics in natural channels with complex configurations. They also highlight the need for detailed spatially distributed datasets of three‐dimensional flow variables to establish the accuracy and applicability of CFD software. Copyright © 1999 John Wiley & Sons, Ltd.