River#
River storage and related fluxes are calculated based on section 4.7 of Müller Schmied et al 2021 [1].
Water balance#
River storage \(S_r\) \([m^3]\) is computated as
where \({Q}_{r,in}\) is inflow into the river compartment [\({m}^{3} {d}^{-1}\)], \({Q}_{r,out}\) is the streamflow [\({m}^{3} {d}^{-1}\)], and \({NA}_{s,r}\) is the net abstractions of surface water from the river [\({m}^{3} {d}^{-1}\)].
Inflows#
The inflow \({Q}_{r,in}\) into the river compartment (if there are no surface water bodies) is the sum of soil surface runoff \(({R}_{s})\), groundwater discharge \(({Q}_{g})\), and upstream streamflow. Otherwise fraction of \({R}_{s}\) and \({Q}_{g}\) (in humid cells) is routed through the surface water bodies (See: Model Schematic). The outflow from the surface water body preceding the river compartment then becomes part of \({Q}_{r,in}\). In addition, negative net abstractions \(({NA}_{s})\) values due to high return flows from irrigation with groundwater lead to a net increase in storage. Thus, if no surface water bodies exist in the cell, negative \({NA}_{s}\) is added to \({Q}_{r,in}\).
Outflows#
\({Q}_{r,out}\) is defined as the streamflow or river discharge that leaves a cell and is transported to the next cell downstream.
This is calculated as:
where \({v}\) \([m/s]\) is the velocity of the river flow, and \({l}\) is the river length \([m]\). \({l}\) is calculated as the product of the cell’s river segment length, derived from the HydroSHEDS drainage direction map [2], and a meandering ratio specific to that cell described in Verzano et al., 2012 [3].
The velocity \({v}\) is calculated using the Manning–Strickler equation:
where \({n}\) is river bed roughness \([–]\), \({R}_{h}\) is the hydraulic radius of the river channel \([m]\) and \({s}\) is river bed slope \([{m}*{m}^{-1}]\) [1] .
Daily varying \({R}_{h}\) is calculated assuming a trapezoidal river cross section with a slope of 0.5. \({R}_{h}\) then can be calculated as a function of daily varying river depth \({D}_{r}\) and temporally constant bottom width \({W}_{r,bottom}\) [3]. WaterGAP implements a consistent method for determining daily width and depth as a function of river water storage. Bankfull flow conditions are assumed to occur at the initial time step and the initial volume of water stored in the river is calculated as:
where \({S}_{r,max}\) is the maximum volume of water that can be stored in the river at bankfull depth \([{m}^3]\), \({D}_{r,bf}\) \([{m}]\) and \({W}_{r,bf}\) \([{m}]\) are river depth and top width at bankfull conditions, respectively, and \({W}_{r,bottom}\) is river bottom width \([{m}]\). River water depth \({D}_{r}\) \([{m}]\) is simulated to change at each time step with actual \({S}_{r}\) as:
WaterGap also computes net cell runoff \({R}_{n,c}\) \([{m}*{m*d}^{-1}]\), which is the part of the cell precipitation that has neither been evapotranspirated nor stored. It is calculated as:
where \({A}_{cont}\) is the continental area (0.5° times 0.5° grid cell area minus ocean area) of the grid cell (\([m^2]\)). Renewable water resources are calculated as long-term annual mean \({R}_{nc}\) under naturalized conditions (without human impact). It is possible for renewable water resources to be negative if evapotranspiration in a grid cell is higher than precipitation due to evapotranspiration from global lakes, reservoirs or wetlands that receive water from upstream cells.