River2D

Quick summary

Figure 1: River 2D layout images: (a) R2D_Bed; (b) R2D_Mesh; and (c) River2D.

Figure 2: Water depth (m) output from River 2D. Ocreza River, Portugal.

Figure 3: Water velocity (m/s) output from River 2D. Ocreza River, Portugal.

Figure 4: Weighted Usable Area (WUA) (m2) for the European Eel. Ocreza River, Portugal.

File:River2d wua discharge.jpg

Figure 5: Weighted Usable AreaWUA (m2) and discharge in the Ocreza River for Iberian barbel (adult and juvenile), European eel (adult and juvenile) and Gudgeon.

Developed by: collaboratively between the Freshwater Institute in Winnipeg, the Civil and Environmental Department of the University of Alberta in Edmonton, the Midcontinent Ecological Science Center of the U.S. Geological Survey in Ft. Collins, and the Fisheries Division of the Alberta Government in Cochrane.

Date: 2002

Type: Tool

Suitable for the following [[::Category:Measures|measures]]:

Introduction

The River 2D is a two dimensional depth-averaged finite element hydrodynamic model that has been customized for fish habitat evaluation studies (Steffler, 2000). The River2D model simulates hydraulic conditions in natural rivers from topographic data input, and uses the habitat suitability curves containing known biological preference data, to calculate the potential habitat for specific species life-history stages by obtaining the Weighted Usable Area (WUA).

The model suite consists of four programs: R2D_Bed, R2D_Ice, R2D_Mesh and River2D (Figure 1). R2D_Bed was designed for editing bed topography data while R2D_Ice is intended for developing ice topographies to be used in the modelling of ice-covered domains. The R2D_Mesh program is used for the development of computational meshes that will ultimately be input for River2D. The latter is then used to solve for the water depths and velocities using the 2D shallow water equations throughout the discretized domain. The habitat model to calculate the WUA is incorporated in the River2D program.

Application

The River2D can be used to calculate the water depths and velocities for any point in the river reach. The final goal is to calculate the WUA for a given fish species and life-stages by using Habitat Suitability Curves.

As input data, the model requires channel bed topography, roughness and transverse eddy viscosity distributions, boundary conditions, and initial flow conditions. Boundary conditions usually take the form of a specified total discharge at inflow sections and fixed water surface elevations at outflow sections. Locating flow boundaries some distance away from areas of interest is important to reduce the effect of boundary conditions uncertainties. The 2D finite element model works with a triangular irregular mesh.

Topography forms the template of the modelling mesh which is draped over the boundary surface and which forms the basis for the numerical solution to the governing flow equations. The challenge of building a mesh is to distribute the nodes in such a way that the most accurate solution is obtained for a particular purpose. Closely spaced nodes in areas of high interest, gradual changes in node spacing, and regularity of elements are important considerations.

Outputs from the model are two (horizontal) velocity components and a depth at each point or node (Figure 2 and Figure 3). Velocity distributions in the vertical are assumed to be uniform and pressure distributions are assumed to be hydrostatic.

The fish habitat component of River2D is based on the Weighted Usable Area (WUA) (Bovee et al., 1998) concept used in the PHABSIM family of fish habitat models. The WUA is calculated as an aggregate of the product of a Composite Suitability Index (CSI, range 0.0 - 1.0) evaluated at every node in the domain and the "tributary area" associated with that node in the finite element mesh. The Suitability Index (SI) for each parameter is evaluated by linear interpolation from an appropriate Habitat Suitability Curves (HSC) to be supplied separately. Velocities and depths are taken directly from the hydrodynamic component of the model. The channel index values may depend on channel substrate or cover for different fish species and life stages. These values are interpolated from a separate channel index file to the computational nodes. The interpolation may be linear (continuous) or nearest neighbour (discrete). At the last stage, the River 2D computes the WUA (Figure 4) as the product, harmonic mean, or minimum value.

For habitat studies, the assumption is that the habitat is optimal when it reaches a maximum value of WUA (Figure 5).

The River 2D has been widely used to access the environmental flows (Rivaes et al., 2017), to study fish movements due to hydropeaking (Boavida et al., 2017; Yarnell et al., 2012), to calculate the suitability of spawning grounds (Gard, 2009), and for river restoration studies (Boavida et al., 2010; García de Jalón and Gortázar, 2007; Lacey and Millar, 2004).

Other information

In addition to the use on regular Windows and Linux workstations, the upcoming version 3.0 of BASEMENT (March 2019) is designed to run on graphical processing units (GPUs) and distributed memory computer clusters allowing engineers to tackle problems with very large computational domains or long simulation time.

Relevant literature

• Vetsch D., Siviglia A., Caponi F., Ehrbar D., Gerke E., Kammerer S., Koch A., Peter S., Vanzo D., Vonwiller L., Facchini M., Gerber M., Volz C., Farshi D., Mueller R., Rousselot P., Veprek R., Faeh R. 2018. System Manuals of BASEMENT, Version 2.8. Laboratory of Hydraulics, Glaciology and Hydrology (VAW). ETH Zurich. Available from http://www.basement.ethz.ch.

Contact information

Website (download, tutorials, etc.): http://www.basement.ethz.ch

Email: basement@ethz.ch