logo CRAS
Comptes Rendus. Géoscience
Original Article - External Geophysics, Climate
Variable backwater and channel roughness: effects on Solimões River discharge
Comptes Rendus. Géoscience, Volume 352 (2020) no. 3, pp. 185-198.


Historically, streamflow daily values are obtained indirectly by applying a rating curve on water level data. In several Amazon basin stations, a backwater effect impacts the flow, thereby invalidating the rating curve method. To account for this effect, we generated a daily discharge data series by applying the Manning equation to daily water level data from the Manacapuru station located on the Solimões River. Each Manning parameter had to be established based on its relationship to the water level. The roughness coefficient was obtained by inverting the Manning equation for each discharge measurement, and a “coefficient–stage” relationship was established. The slope of the energy line was estimated by the water-surface slope between two fluviometric stations. Discretizing the equations for each parameter enhanced the discharge estimation for the Manacapuru station, thereby generating an accurate discharge data series.

Historiquement les valeurs journalières de débit sont obtenues de manière indirecte en appliquant une courbe de tarage aux valeurs correspondantes de hauteur d’eau. Mais pour beaucoup de stations du bassin amazonien le débit mesuré est affecté par un effet de barrage hydraulique dû à l’influence variable du débit entrant d’un affluent et l’on ne peut plus appliquer une simple courbe hauteur/débit. Pour prendre en compte cet effet de mascaret, nous avons créé pour la station de Manacapuru sur le Río Solimoes une série de débits journaliers en appliquant l’équation de Manning aux hauteurs journalières correspondantes. Chaque paramêtre de Manning a été établi en fonction du niveau d’eau. Le coefficient de rugosité a été calculé en inversant les termes de l’équation de Manning pour chaque mesure de débit et une courbe “coefficient de rugosité–niveau” a pu être établie. La pente de la ligne d’énergie a été estimée à partir de la pente de la ligne d’eau avec la station hydrométrique la plus proche. En discrétisant les équations pour chaque paramètre nous avons pu améliorer l’estimation du débit à la station de Manacapuru et générer ainsi une série de données de débits plus précise.

Online First:
Published online:
DOI: 10.5802/crgeos.35
Keywords: Backwater effect, Manning equation, Roughness coefficient, Rating curve, Maximum flow, Solimões River
Luna Gripp Simões Alves 1; Demetrius David da Silva 2; Philippe Vauchel 3; Pascal Fraizy 3; Naziano Pantoja Filizola 4

1 Geological Survey of Brazil (CPRM), Federal University of Viçosa (UFV), Manaus-AM, Brazil
2 Federal University of Viçosa (UFV), Viçosa-MG, Brazil
3 French National Research Institute for Sustainable Development, France
4 Federal University of Amazonas, Manaus-AM, Brazil
     author = {Luna Gripp Sim\~oes Alves and Demetrius David da Silva and Philippe Vauchel and Pascal Fraizy and Naziano Pantoja Filizola},
     title = {Variable backwater and channel roughness: effects on {Solim\~oes} {River} discharge},
     journal = {Comptes Rendus. G\'eoscience},
     pages = {185--198},
     publisher = {Acad\'emie des sciences, Paris},
     volume = {352},
     number = {3},
     year = {2020},
     doi = {10.5802/crgeos.35},
     language = {en},
TI  - Variable backwater and channel roughness: effects on Solimões River discharge
JO  - Comptes Rendus. Géoscience
PY  - 2020
DA  - 2020///
SP  - 185
EP  - 198
VL  - 352
IS  - 3
PB  - Académie des sciences, Paris
UR  - https://doi.org/10.5802/crgeos.35
DO  - 10.5802/crgeos.35
LA  - en
ID  - CRGEOS_2020__352_3_185_0
ER  - 
%0 Journal Article
%T Variable backwater and channel roughness: effects on Solimões River discharge
%J Comptes Rendus. Géoscience
%D 2020
%P 185-198
%V 352
%N 3
%I Académie des sciences, Paris
%U https://doi.org/10.5802/crgeos.35
%R 10.5802/crgeos.35
%G en
%F CRGEOS_2020__352_3_185_0
Luna Gripp Simões Alves; Demetrius David da Silva; Philippe Vauchel; Pascal Fraizy; Naziano Pantoja Filizola. Variable backwater and channel roughness: effects on Solimões River discharge. Comptes Rendus. Géoscience, Volume 352 (2020) no. 3, pp. 185-198. doi : 10.5802/crgeos.35. https://comptes-rendus.academie-sciences.fr/geoscience/articles/10.5802/crgeos.35/

Full text

1. Introduction

Streamflow data provide information on water availability across time and space and constitute the main source of information used in hydrological studies such as water resource planning and management [WMO 2010], hydroelectric energy generation [Stickler et al. 2013], flood forecasting and control [MoussaBocquillon 2009], and understanding the hydrological cycle from a global perspective [Huntington 2006].

Population centers in the Amazon basin are located along the main river channels, where water is used for both consumptive purposes and transport of people and consumer goods. In this region, the modeling of hydrological systems focused on understanding the occurrence of extreme events is essential for mitigating the related negative impacts [Collischonn et al. 2007; Getirana et al. 2009; Guimberteau et al. 2012; Paiva et al. 2011; Yamazaki et al. 2011].

Discharge series provide key information for calibrating and validating these hydrological models at several scales, and they are pivotal to understanding the impacts caused by both local changes in land cover and climate change on a larger scale [Mahmoud et al. 2009; Peña-Arancibia et al. 2015].

Historically, the most common method used to obtain a continuous discharge data series is the conversion of a water stage data series via the application of a curve that describes the relationship between the stage and the discharge in a specific cross section of the river. This approach is called the “discharge rating” or “rating curve” method [WMO 2010].

The process typically consists of fitting a parametric function such as Q(t) = a(h(t) − h0)c, where Q(t) represents the discharge at time t, h(t) is the concurrent water level, h0 is the water level corresponding to zero flow, and a and c are coefficients related to cross-sectional characteristics and are calibrated for each specific station [Rantz 1982]. Using this approximation is essential since the discharge field measurements are complex and expensive, and in most cases, they remain impracticable for continuous monitoring [Durand et al. 2016; Muste et al. 2011].

However, for this approximation to be valid, certain hydraulic conditions must be observed such as permanent flow and stable channel. Ideally, the stream gauge must be installed in a channel where there is a feature at the downstream section that forces the flow through a critical depth such as a rock riffle or artificially constructed weir. The downstream section feature “controls” the flow, creating a unique relationship between the stage and the discharge [Holmes et al. 2016]. These conditions, however, are not always feasible in the Amazon basin. In this poorly inhabited basin, the choice of gauge sites depends first on the proximity of inhabitants who can perform the hydrological observations.

Thus, gauges are frequently installed next to the confluence of large rivers, where an ambiguous relationship can occur between the stage and the discharge, under a backwater effect induced by changes in downstream water levels [Mansanarez et al. 2016; WMO 2010]. In large Amazon rivers, even in the main channel several kilometers from the mouth, this backwater effect also exists and can be caused by water changes between the main channel and the floodplains [Meade et al. 1991]. In hydraulic terms, this variable backwater affects the flow energy line and renders it non-parallel to the channel bottom. Thus, it must be considered as a complementary variable when estimating discharge [Jaccon and Cudo 1989].

The stage–fall–discharge method and the constant-fall method are the most commonly applied approaches for stage–discharge relations affected by variable backwater [ISO 9123 2001]. Their application requires the consideration of the slope of the energy line in the stretch, which is estimated by the difference between the water levels at the main station and at an auxiliary station located on the same river [ISO 9123 2001; Petersen-Øverleir and Reitan 2009].

These processes are supported by the method in Hall [1916], which is based on the Chezy–Manning equation. Specifically, the “real” discharge Qr has been proposed, which considers the backwater effect, and it can be calculated as Qr=QnSrSn, where Qn is the “reference discharge”, Sr is the “real” slope, and Sn is the “normal” slope of the energy line.

Figure 1.

Overview maps of the Solimões catchment and the Solimões River between the Manacapuru and Careiro stations. Elaborated by authors. Source of the satellite image: Google Earth.

The application of these methods can produce acceptable discharge estimates for flows experiencing a backwater effect. Nevertheless, there is an additional factor that must be considered when determining an accurate discharge series: the variation in roughness along the channel [Mansanarez et al. 2016].

In natural channels, resistive forces exerted by the banks and bed over the flow can vary as the water level changes. The resistance caused by friction from deeper flows can be less expressive than that from shallow flows [Pan et al. 2016]. However, when the flow reaches the floodplain, which is common in the Amazon plain [Rudorff et al. 2014], the resistive forces exerted by more densely vegetated surfaces can represent an increase in the roughness for that water level. Therefore, in situations where these variations are significant, it is important to consider their effects when calculating discharge [Le Coz et al. 2014].

Although the Manning equation was adapted from the Chezy equation with the aim of forcing the roughness coefficient to be constant in a specific section, several studies reported that for natural courses, the roughness coefficient varies as a function of the water level [Dingman 2009; Ferguson 2010; Le Coz et al. 2014; Pan et al. 2016]. Therefore, based on the knowledge of the slope of the energy line and the relationship between the roughness coefficient and the water level, it is possible to estimate the discharge using the Manning equation while still considering the effects of both the variable roughness and the backwater.

The aim of this work was to generate an accurate daily discharge data series for the Solimões River, which experiences flow under the backwater effect. Variations in the energy slope and the roughness of the channel had to be considered in the computations. Thus, it was necessary to establish the relationship between each Manning equation parameter and the water level in the channel.

2. Materials and methods

2.1. Study area

The proposed method for computing streamflow was tested for the Manacapuru station (Code 14100000). The Manacapuru station is located on the Amazon plain, 80 km upstream from the confluence with the Negro River. It is the most downstream station of the Solimões River. This confluence marks the start of the Amazonas River. The Manacapuru station represents a drainage basin of approximately 2.2 million km2 or 36% of the total Amazon basin area (Figure 1).

Figure 2.

Cross-sectional profile of the Manacapuru station obtained from acoustic discharge measurements performed in June 2015 (2077 cm stage).

In this Amazon plain region, the flow dynamics of big rivers induce an important effect of backwater, which can be observed several kilometers from the river mouths. This effect is a result of the seasonal storage of water on floodplains and is associated with timing differences in tributary inflow [Meade et al. 1991; Trigg et al. 2009].

The downstream backwater effect on the Manacapuru station affects the relationship between the stage and the discharge in this station, and a unique rating curve is not able to accurately calculate the discharge data series [Jaccon 1986]. In hydraulic terms, the downstream backwater effect can be quantified as the variation in the dimensionless hydrostatic pressure gradient [Tsai 2005].

2.2. Modeling methods and data

To consider the downstream backwater effect in the discharge calculation, we propose to estimate the discharge by the Manning equation:

where Q(m3s1) is the streamflow, A (m2) is the wetted area, Rh (m) is the hydraulic radius, Sf (m⋅m−1) is the slope of the energy line of the flow, and n (s⋅m−1∕3) is the Manning roughness coefficient.

The variation in the hydrostatic pressure gradient over time was represented by the variation in the “energy slope” parameter of the Manning equation. The other parameters of the Manning equation were established as a function of the water level, which is explained in Section 3. Then, the Manning equation was applied to the daily stage series to generate the discharge data series.

The geometric parameters of (1) (the area A and hydraulic radius Rh parameters) were estimated as functions of the water level based on a cross-sectional profile obtained from an acoustic profiler discharge measurement performed in high waters (Figure 2). The established profile could be used to simulate the geometric parameters at every daily water level.

The “slope of the energy line” parameter was estimated as the “water-surface slope” between the Manacapuru and an auxiliary station. This estimate assumed that variations in the velocity head are negligible in the stretch and the flow gradually varies. The auxiliary station was the Careiro station (Code 14040000), which is located 90 km downstream from the Manacapuru station at a branch of the Amazonas River. The water levels of the Careiro station were corrected to the Manacapuru datum based on information provided by Moreira et al. [2010]. The water-surface slope was calculated as the difference between the Manacapuru and Careiro water levels divided by the distance between the two stations (Figure 1).

To establish the relationship between the Manning roughness coefficient n and the water level, we used discharge field measurement data and the corresponding water level. With the discharge value available and all the other parameters of the Manning equation calculated for the associated water level, the Manning roughness coefficient n was computed by inverting (1). Then, for each measurement, we obtained the n coefficient, which was plotted with the corresponding water level, and fitted the equation that best described the relationship. With the relationship established, the roughness coefficient could be inversely estimated for daily data at each water level.

Then, all the parameters of the Manning equation were calculated as a function of daily water levels, and the equation was applied to generate a daily discharge data series. The software Hydraccess [Vauchel 2005] was used to facilitate the estimation and calculation of the Manning parameters.

The water stage, discharge measurements, and daily discharge series generated by traditional rating curves were obtained from Brazilian National Water Agency (ANA) databases (available at http://www.snirh.gov.br/hidroweb). For the present analyses, only the discharge measurements collected using acoustic equipment were considered because previous results indicated that the method by which the measurements are performed significantly affects the relationship between the water level and the discharge at the Manacapuru station. The statistical analyses and the comparison of hydrographs that supported these findings are described in Alves et al. [2017b].

Figure 3.

Relationship between the stage (cm) and discharge measurements (m3⋅s−1) at the Manacapuru station. Data are discretized by the average velocity.

Usually, discharge field measurements were performed once every 4 months in the Amazon basin. For this study, we complemented the traditional database with information from two projects by the Geological Survey of Brazil (SGB), “Sistemas de Alerta Hidrológico” and “Dinâmica Fluvial do Sistema Negro-Solimões-Amazonas”, which are supported by the SGB and the ANA (i.e., the Brazilian National Water Agency). In these contexts, discharge measurements were performed every month for 8 years (from 2009 to 2016) using acoustic Doppler current profilers (ADCPs) at 300 or 600 kHz coupled to Global Positioning System (GPS) equipment on a large vessel. An example of an ADCP measurement result is shown in Figure 2. The use of GPS was important because this station presents a moving bed during certain periods of the year, which can make it difficult to obtain the position by the bottom tracking of ADCP [Mueller et al. 2009]. This complemented database generated an unprecedented amount of detailed information for the Amazon river.

The 85 discharge measurements available varied between 30 500 and 163 300 m3⋅s−1, with an average of 107 889 ± 33 893 m3⋅s−1 (average ± standard deviation). They were obtained between the stages of 4.71 and 20.77 m, which basically encompass the entire range of stages observed in the historical series of the Manacapuru station. From the 16 245 daily stage data points available in the Manacapuru series, only 15 (0.09%) were lower than the minimum measurement stage and only 1 (0.01%) was higher than the stage observed for the maximum discharge measurement.

2.3. Model evaluation

The proposed methodology for discharge estimates was evaluated by analyzing statistical indexes that are typically used in hydrological modeling, that is, the mean absolute error (MAE) (2), the root mean square error (RMSE) (3), and the Nash–Sutcliffe efficiency (NSE) index (3) [Moriasi et al. 2007].

where (Q1obs,,Qkobs) are the observed discharges (from discharge field measurements), (Q1calc,,Qkcalc) are the associated calculated discharges, k is the number of discharge measurements, and Qobs¯ is the average discharge of all observed data. To compare the proposed methodology with the traditional approach, the rating curve discharge data series available in Hidroweb was also submitted to the statistical index calculations and the results were compared.

Figure 4.

Monthly box plots of the slopes of the surface lines for the complete water level series from Manacapuru and Careiro stations (1977–2016). The dashed line represents the average slope value. Upper and lower vertical line extremes indicate the maximum and minimum values, respectively, and box boundaries represent first and third quarters.

3. Results

The relationship between the discharge field measurements and stages is presented in Figure 3. Data are discretized by color as a function of the average velocity (measured discharge divided by section area) obtained in the measurements. As can be observed in the figure, the relationship between discharge and stage is not unique along the total range of observed data.

For example, at an intermediary stage of 16.00 m, the discharge data varied between 99 160 and 124 300 m3⋅s−1, presenting a difference of 25% in the discharge values. For lower stages, dispersion also occurs. The discharge associated with a stage of 7.65 m (40 000 m3⋅s−1) was 31% lower than the discharge measured for a stage of 6.65 m (56 439 m3⋅s−1). The full dispersion of stage–discharge along the plot seems to indicate that the backwater effect is active at all stages of variation observed in the database.

Analyzing the different patterns of average velocities shows that the dispersion between points is caused by different flow velocities associated with a determined stage value. For example, for the same discharge of 100 000 m3⋅s−1, the velocity varies between 1.15 and 1.34 m⋅s−1 and is associated with stages between 16.42 and 13.23 m (Figure 3).

In hydraulic terms, these different velocity patterns can be explained by the variation in the slope of the energy line of the flow, which is induced by the variable backwater effect. In this case, the flow is controlled by the downstream level.

3.1. Slope of the energy line

The variation in the slope of the energy line was estimated as the variation in the water-surface slope between the Manacapuru and Careiro stations. This variation was visible throughout the year, mainly between August and March (Figure 4).

For the months of April, May, and June, that is, the period when the maximum stages occur at Solimões River, the variability of the slope (shown by the height of the box plots) is notably lower than that observed in other months. However, the water-surface slope is shown to be higher in magnitude in these months than the annual average. From August to October, the slopes reach the minimum values (with a higher variability), and between November and January, they tend to increase in magnitude.

Variations in the water-surface slope along this stretch are mainly caused by the storage effect of floodplains, which is present in a major part of the Solimões and Amazonas systems [Meade et al. 1991; Rudorff et al. 2014]. When the flow affluent to the Manacapuru station starts to decrease at the end of the full period, which occurs between June and July, the water level at this station also starts to decrease. As the water level in the channel falls, the water that was stored in the floodplains during the full period starts to return to the channel. Then, the water level at the downstream station is reduced by a lower rate than that at the Manacapuru station, thus causing a softer downstream recession curve [Alsdorf et al. 2000].

Figure 5 presents the stage records observed at Manacapuru and Careiro stations in 2015 (the year when the maximum stage occurred at the Manacapuru station) to illustrate this process. From August, when the downstream level falls at a smaller rate relative to the upstream station, the water-surface slope tends to be reduced as seen in the figure.

Figure 5.

Stage (cm) and slope (m⋅m−1) variations in 2015. Careiro stages converted to Manacapuru datum based on data provided by Moreira et al. [2010].

An inverse pattern is observed when the flow affluent to the Manacapuru station begins to increase. The water level increases in the main channel, and the water starts to flow from the channel to the floodplains. The water that overflows to the floodplains delays the rate at which the downstream level increases. In this case, the water-surface slope starts to increase. In Figure 5, this process starts between October and November. In Figure 4, the process can be observed from the end of October to December.

Changes in slope magnitude are associated with the rising and falling processes of the river, and the monthly variation in the slope (illustrated by the height of box plots in Figure 4) can be explained based on the time when the inversion of rising and falling limbs occurs.

For the rising limb, the maximum annual stages occurred in June and July in the Manacapuru station in 97% of the years in the database. Therefore, in May, June, and July in almost all the years, the rivers and the floodplains are full and the slope magnitudes are always similar. This results in an almost constant value of the water-surface slope of approximately 4.4 × 10−5 m⋅m−1 in the stretch.

Figure 6.

Area and hydraulic radius of the Manacapuru section calculated as a function of the water level.

For the falling limb, the minimum annual stage has been recorded between the months of October and January. For these months, depending on when the processes have changed, the slopes can be characteristic of the rising or falling limbs and vary between 3.0 and 5.0 × 10−5 m⋅m−1.

It is worth noting that the changes between the channel and floodplains occur not only by the river overflowing its banks but also via breaks in the lateral levees in the Amazon system [Rudorff et al. 2014; Trigg et al. 2012]. Therefore, the processes of filling and draining of floodplains take place even before the occurrence of the bankfull level, and such dynamics impact the whole level range throughout the year.

Figure 7.

Manning roughness coefficients (s⋅m−1∕3) estimated for each field discharge measurement as a function of the measured water level (cm).

The effects of floodplain storage on the slope variability and, consequently, on the stage–discharge relationship have already been investigated by Meade et al. [1991] for several stations along the Solimões River. However, these authors did not have access to data to correct water levels to a unique datum. Thus, they recommended discretizing the stage–discharge data in rising and falling periods to improve the discharge estimates. The present results indicate that rising/falling discretization is not clear for this station and the consideration of slope parameters is paramount to obtaining an accurate data series for the Manacapuru station.

3.2. Geometric parameters

The relationships between the wetted area and hydraulic radius parameters and water level are presented in Figure 6. The relationship between the wetted area and the water level presents a unique trend and is approximately linear throughout the whole range of stage variations. For the hydraulic radius, however, a break can be observed as the right bank becomes flatter next to the 13.00 m stage (see Figure 2). For this flat region, a small increase in water level causes a large increase in the wetted perimeter and a small increase in the wetted area. The hydraulic radius can thus decrease even under a rising water level.

3.3. Manning roughness coefficient

The estimated Manning coefficient varied between 0.035 and 0.063 s⋅m−1∕3 with an average of 0.045 ± 0.004 s⋅m−1∕3 (average ± standard deviation) as presented in Figure 7. The plot of the Manning roughness coefficient and the water level shows that they are inversely related for low and medium stages (Figure 7). This pattern has been already observed in several other studies [Dingman 2009; Ferguson 2010].

The inverse relationship can be explained by the variation in the impact exerted by the resistive forces of the bed and bank friction in natural channel flow [Mansanarez et al. 2016]. For a shallow flow, the impact of the resistive forces exerted by the riverbed is great, resulting in high values of the Manning coefficient. For a deeper flow, this impact is softer; thus, the Manning coefficient decreases.

Such a relationship is also observed from a wider point of view when comparing the roughness coefficients estimated for different water courses of variable magnitude [Ferguson 2010]. By analyzing a compilation of published data in which discharges varied from 0.02 to 3000.00 m3⋅s−1, Ferguson [2010] observed a significantly inverse correlation between n and discharge for several of the analyzed water courses. In their analysis, however, they excluded data where the flow overtopped the bank and was affected by the floodplain vegetation.

When considering the flow over the floodplain, the Manning roughness coefficient–stage relationship is modified and can tend asymptotically toward a constant value [Domeneghetti et al. 2012; Ferguson 2010; Moramarco and Singh 2010]. In the floodplain, which is more densely vegetated, the roughness and Manning coefficients can be higher than those in the main channel, and they increase as the flow progresses toward the plain. For the Manacapuru section, the floodplain and the higher parts of the banks present denser vegetation (grass and large trees with exposed roots; see Figure 8). Therefore, trends in the coefficients for the section as a whole depend on the effects of the roughness variations in the floodplains relative to the variations in the main channel [Le Coz et al. 2014].

Figure 8.

Images from the left bank of the Manacapuru section. (a) Image taken on 08/27/2017 at the 1580 cm stage illustrating the density of exposed tree roots, which affects the surface roughness. (b) Image taken on 07/27/2017 at the 1857 cm stage when the roots were flooded, thus highlighting the considerable presence of grasses on the banks.

Figure 9.

Hydrographs obtained by applying the Manning equation (ME) and rating curve (RC) methodologies and field measurements. Bars represent the uncertainty of discharge measurements (±5.3%) for a 95% confidence interval [ISO 748 1997].

In this work, the function that best fitted the relationship between the roughness coefficient and the water level at the Manacapuru station was a quadratic polynomial with an inflection point near the 17.00 m stage. Researchers studying this relationship for other sections have observed similar results [Le Coz et al. 2014; Pan et al. 2016].

It is worth noting that although the Manning coefficient is a physical parameter related to bed and bank roughness, it has been used as a calibration coefficient in a substantial number of hydraulic and hydrological models. In such cases, the observed values are not always physically interpretable or justifiable [Di Baldassarre et al. 2010]. In the present work, these coefficients were obtained by the direct application of the Manning equation using field discharge measurements performed using the best available technology at this site (ADCP equipment associated with GPS). In fact, the magnitude of the roughness coefficients that was obtained was within the range recommended by the classical literature [Chow 1959]. Therefore, the presented results can be used as reference coefficients for related future studies at this station.

3.4. Discharge calculation

Once all the relationships between the Manning parameters and water levels were established, it was possible to apply this equation to daily stage data, thereby enabling the generation of the daily discharge data series. Hydrographs obtained by applying the Manning equation as well as the data generated by using simple rating curves and discharge field measurements are presented in Figure 9. The results of the model evaluation indexes for discharge estimations are presented in Table 1.

Table 1.

Model evaluation indexes: mean absolute error (MAE), root mean square error (RMSE), and Nash–Sutcliffe efficiency (NSE) for both methodologies (RC—rating curve and ME—Manning equation). The “observed values” are relative to discharge field measurements

MAE (m3⋅s−1)RMSE (m3⋅s−1)NSE

Table 2.

Maximum annual discharge obtained by both methodologies (RC—rating curve and ME—Manning equation)

Maximum annual discharge (m3⋅s−1)
Year RC ME Difference
2009160 615143 90016 715
2010137 967131 3006 667
2011142 567131 20011 367
2012163 467150 30013 167
2013156 932147 0009 932
2014158 953147 50011 453
2015164 780154 10010 680
2016134 695136 000−1 305
Average9 834

An analysis of the statistical indexes shows that both methods can be considered adequate to estimate the discharge at the Manacapuru station. Considering that the average measured discharge is 110 014 m3⋅s−1, the values of both the MAE and the RMSE errors are lower than 10% of the average discharge. For the NSE index, values closer to 1 indicate that the two methods are acceptable for estimating discharge.

In a comparative analysis, estimating discharges using the Manning equation led to better results in terms of the MAE, RMSE, and NSE. The hydrographs presented in Figure 9 support the results of the evaluation indexes, indicating that the main differences are related to the maximum discharge estimates. In the upper peaks, the results obtained by the Manning equation seem to have a better fit to the discharge measurements than the rating curve, which overestimates the maximum peak discharges.

Table 2 presents the maximum annual discharge obtained by each methodology. The differences can reach approximately 17 000 m3⋅s−1 such as for 2009. On average, the differences between the annual discharges estimated by the different methods are approximately 10 000 m3⋅s−1.

For hydrological modeling, this magnitude of variation could be reasonable because the differences for both methodologies are approximately 10% of the average discharge. However, it is of paramount importance to note that this work focuses on the generation of an accurate discharge database that can be used to calibrate and validate several kinds of other models.

As an example of the broad application of this database, discharge data series of the Manacapuru station were used to estimate the solid discharge of the Amazonas River to the ocean [Park and Latrubesse 2014], create a multidecadal hydrological retrospective [Correa et al. 2017], and forecast streamflow [Paiva et al. 2013]. Thus, it is essential for the discharge series to be as accurate as possible.

4. Discussion

Applying the Manning equation has two main advantages relative to the traditional rating curve. Each parameter is considered separately, which allows for a consideration of variations, and the energy slope is included in the computations.

For example, the discontinuity of the “stage–hydraulic radius” curve is not considered in the rating curve calculations although this condition of an irregular cross section can be incorporated into the estimation of the rating curves by using different curves for each range of stage variation [Jaccon and Cudo 1989]. Nevertheless, under the analyzed conditions of the Manacapuru station, this break is not clear in the stage–discharge representation (see Figure 3) since the dispersion caused by the variable energy slope is visually more expressive.

The observed data dispersion in the stage–discharge relationship is mainly caused by a variable backwater effect associated with the changes between channels and floodplains. To consider this factor in calculations, the slope of the energy flow line estimated by the water-surface slope was considered in the discharge estimates. According to WMO recommendations [WMO 2010], to estimate the water-surface slope between two stations, the distance between the main and the auxiliary station should be enough to establish a difference of 15 cm in observed stages. Nevertheless, we have to consider that the distance between them is sufficiently short to have a linear water profile between the two gauges. The stretch between the stations should be as uniform as possible because the contractions and obstructions in the channel cannot cause significant breaks in the water-surface slope between the gauges [Petersen-Øverleir and Reitan 2009]. The auxiliary station should be sufficiently far from the backwater source [WMO 2010]. However, in practical terms, these assumptions are difficult to achieve under Amazonian conditions.

In the present work, the observed falls between the stations varied between 2.70 m and 4.71 m. The gauges are separated by a distance of 90 km in a stretch defined as “anabranching” [Latrubesse 2015], which includes levees and large islands as shown in Figure 1. From a theoretical point of view, this stretch is far from ideal for considering the kinetic energy loss constant and estimating the slope of the energy line as the slope of the surface between the two stations. Nevertheless, the fine adjustments of the estimated discharge series to the field discharge measurements indicate that it is fully feasible to consider these assumptions in a region with such a flat relief such as the Solimões basin.

Because the roughness coefficient decreases as a function of the water level and tends to reach a constant value, these results also reflect the importance of detailing these relationships, which is especially significant when considering the fact that in discharge estimates, the calculated discharge is inversely proportional to the coefficient at a one-to-one rate. Considering an average coefficient instead of using a variable coefficient could lead to an overestimation of the minimum discharge and an underestimation of the maximum flows. Such errors are undesirable when considering that the discharge data are used to calibrate and validate hydrological models for extreme events. These results were described by Alves et al. [2017a] for the Solimões River at the Manacapuru station and corroborated by Mansanarez et al. [2016] for the Madeira River.

Finally, it is worth noting that the present analysis was possible only due to the availability of a robust field discharge database with temporal details sufficiently fine to represent the dominant hydrological processes. A database with this amount of detail is not available for most Brazilian stations.

To apply the methodology for stations without such detailed data, it is likely that the relationships obtained between the parameters and the water level will have to be extrapolated, which is the standard procedure for several other methodologies. However, the advantage of applying the Manning equation is related to the possibility of analyzing each parameter separately. For the present station, it is reasonable to understand that by analyzing each relationship obtained for the parameters (area, hydraulic radius, and roughness coefficient), similar results could be obtained even when extreme data are not available.

For the roughness coefficient, the best fit was found for a quadratic equation. However, from a practical perspective, an initial stretch that is inversely proportional to the water level should be considered, followed by a second stretch that tends toward a constant value. Such considerations can be applied for other stations without sufficient data to cover the entire range of water level variations.

At the Manacapuru station and several other stations located on large Amazon rivers, the hydrographs are stable over time and present a unique period of rising and a unique period of falling. Therefore, the discharge measurements performed once a month were sufficient to adequately represent the evolution of hydrological processes. For stations located in smaller catchments with faster responses to precipitation events, the frequency of field measurements must be enhanced by focusing on the measurements of extreme events.

Hence, it is important to highlight the importance of dedicating efforts toward obtaining satisfactory field databases for establishing an adequate stage–discharge relationship and determining the best method to estimate discharge. This information improves the generation of continuous discharge data series, especially for flows affected by variable backwater effects.

5. Conclusions

The results verified the importance of considering the variable backwater effect on discharge calculations for water level data series. The slope of the energy line, which is estimated by the water-surface slope between two gauges, can be adequately applied for this purpose. Individualizing each of the Manning equation parameters as a function of the water level is also extremely important for enhancing the quality of the generated data series. When these data are available, the Manning equation can be applied for each daily water level, thereby generating an accurate discharge data series. Applying the Manning equation appeared to enhance the results compared with the traditional methodology of using rating curves, especially in terms of the maximum discharge peaks. To obtain this accurate discharge data series, it is crucial to have a robust field database to support the development and evaluation of the best method for establishing stage–discharge relationships.


[Alsdorf et al., 2000] D. E. Alsdorf; J. M. Melack; T. Dunne; L. A. K. Mertes; L. L. Hess; L. C. Smith Interferometric radar measurements of water level changes on the Amazon flood plain, Nature, Volume 404 (2000) no. 6774, pp. 174-177 | Article

[Alves et al., 017a] L. G. S. Alves; D. D. Silva; N. P. Filizola; F. F. Pruski Estimativa do coeficiente de Manning para cálculo de vazão em regime sob efeito de remanso hidráulico na bacia Amazônica, XX Simpósio Brasileiro de Recursos Hídricos, Associação Brasileira de Recursos Hídricos, Porto Alegre-RS, Brazil, 2017a

[Alves et al., 017b] L. G. S. Alves; D. D. Silva; N. P. Filizola; F. F. Pruski Stage-discharge relation in non-uniform flow based on Strickler–Manning equation on Amazon basin, ASABE Annual International Meeting, American Society of Agricultural and Biological Engineers, Spokane, Washington, 2017b

[Chow, 1959] V. T. Chow Open Channel Hydraulics, McGraw-Hill, New York, New York, USA, 1959

[Collischonn et al., 2007] W. Collischonn; D. Allasia; B. C. Da Silva; C. E. M. Tucci The MGB-IPH model for large-scale rainfall—runoff modelling, Hydrol. Sci. J., Volume 52 (2007) no. 5, pp. 878-895 | Article

[Correa et al., 2017] S. W. Correa; R. C. D. de Paiva; J. C. Espinoza; W. Collischonn Multi-decadal hydrological retrospective: case study of Amazon floods and droughts, J. Hydrol., Volume 549 (2017) no. Supplement C, pp. 667-684 | Article

[Dingman, 2009] S. L. Dingman Fluvial Hydraulics, Oxford University Press, New York, New York, USA, 2009

[Di Baldassarre et al., 2010] G. Di Baldassarre; G. Schumann; P. D. Bates; J. E. Freer; K. J. Beven Flood-plain mapping: a critical discussion of deterministic and probabilistic approaches, Hydrol. Sci. J.–J. Des Sci. Hydrol., Volume 55 (2010) no. 3, pp. 364-376 | Article

[Domeneghetti et al., 2012] A. Domeneghetti; A. Castellarin; A. Brath Assessing rating-curve uncertainty and its effects on hydraulic model calibration, Hydrol. Earth Syst. Sci., Volume 16 (2012) no. 4, pp. 1191-1202 | Article

[Durand et al., 2016] M. Durand; C. J. Gleason; P.-A. Garambois; D. Bjerklie; L. C. Smith; H. Roux; E. Rodriguez; P. D. Bates; T. M. Pavelsky; J. Monnier An intercomparison of remote sensing river discharge estimation algorithms from measurements of river height, width, and slope, Water Resour. Res., Volume 52 (2016) no. 6, pp. 4527-4549 | Article

[Ferguson, 2010] R. Ferguson Time to abandon the Manning equation?, Earth Surf. Process. Landf., Volume 35 (2010) no. 15, pp. 1873-1876 | Article

[Getirana et al., 2009] A. C. V Getirana; M.-P. Bonnet; S. Calmant; E. Roux; O. C. Rotunno Filho; W. J. Mansur Hydrological monitoring of poorly gauged basins based on rainfall–runoff modeling and spatial altimetry, J. Hydrol., Volume 379 (2009) no. 3, pp. 205-219 | Article

[Guimberteau et al., 2012] M. Guimberteau; G. Drapeau; J. Ronchail; B. Sultan; J. Polcher; J. M. Martinez; C. Prigent; J.-L. Guyot; G. Cochonneau; J. C. Espinoza Discharge simulation in the sub-basins of the Amazon using ORCHIDEE forced by new datasets, Hydrol. Earth Syst. Sci., Volume 16 (2012) no. 3, pp. 911-935 | Article

[Hall, 1916] M. R. Hall A Method of Determining the Daily Discharge of Rivers if Variable Slope, USGPO, Washington D.C., Washington, USA, 1916

[Holmes Jr., 2016] R. R. Holmes Jr. River rating complexity, River Flow (G. Constantinescu; M. Garcia; D. Hanes, eds.), CRC Press, 2016, pp. 10-14

[Huntington, 2006] T. G. Huntington Evidence for intensification of the global water cycle: review and synthesis, J. Hydrol., Volume 319 (2006) no. 1–4, pp. 83-95 | Article

[ISO 748, 1997] ISO 748 Measurement of Liquid Flow in Open Channels—Velocity-area Methods, 1997 (Geneva, Switzerland)

[ISO 9123, 2001] ISO 9123 Measurement of Liquid Flow in Open Channels—Stage-fall-discharge Relationships (V. Methods, ed.), 2001 (Geneva, Switzerland)

[Jaccon and Cudo, 1989] G. Jaccon; K. J. Cudo Hidrologia-curva-chave: analise e traçado, 1989

[Jaccon, 1986] G. Jaccon Estudo da curva-chave do posto fluviometrico de Manacapuru: no rio Solimoes, 1986

[Latrubesse, 2015] E. M. Latrubesse Large rivers, megafans and other Quaternary avulsive fluvial systems: a potential “who’s who” in the geological record, Earth-Science Rev., Volume 146 (2015), pp. 1-30 | Article

[Le Coz et al., 2014] J. Le Coz; B. Renard; L. Bonnifait; F. Branger; R. Le Boursicaud Combining hydraulic knowledge and uncertain gaugings in the estimation of hydrometric rating curves: a Bayesian approach, J. Hydrol., Volume 509 (2014), pp. 573-587 | Article

[Mahmoud et al., 2009] M. Mahmoud; Y. Liu; H. Hartmann; S. Stewart; T. Wagener; D. Semmens; R. Stewart; H. Gupta; D. Dominguez; F. Dominguez A formal framework for scenario development in support of environmental decision-making, Environ. Model. Softw., Volume 24 (2009) no. 7, pp. 798-808 | Article

[Mansanarez et al., 2016] V. Mansanarez; J. Le Coz; B. Renard; M. Lang; G. Pierrefeu; P. Vauchel Bayesian analysis of stage- fall- discharge rating curves and their uncertainties, Water Resour. Res., Volume 52 (2016) no. 9, pp. 7424-7443 | Article

[Meade et al., 1991] R. H. Meade; J. M. Rayol; S. C. Da Conceicão; J. R. G. Natividade Backwater effects in the Amazon River basin of Brazil, Environ. Geol. Water Sci., Volume 18 (1991) no. 2, pp. 105-114 | Article

[Moramarco and Singh, 2010] T. Moramarco; V. P. Singh Formulation of the entropy parameter based on hydraulic and geometric characteristics of river cross sections, J. Hydrol. Eng., Volume 15 (2010) no. 10, pp. 852-858 | Article

[Moreira et al., 2010] D. M. Moreira; O. C. Rotunno Filho; S. Calmant Rede de referência altimétrica para avaliação da altimetria por satélites e estudos hidrológicos na região amazônica (Vol. 1, Issue 1) [Universidade Federal do Rio de Janeiro], 2010 (doi:10.1017/CBO9781107415324.004)

[Moriasi et al., 2007] D. N. Moriasi; J. G. Arnold; M. W. Van Liew; R. L. Bingner; R. D. Harmel; T. L. Veith Model evaluation guidelines for systematic quantification of accuracy in watershed simulations, Trans. ASABE, Volume 50 (2007) no. 3, pp. 885-900 | Article

[Moussa and Bocquillon, 2009] R. Moussa; C. Bocquillon On the use of the diffusive wave for modelling extreme flood events with overbank flow in the floodplain, J. Hydrol., Volume 374 (2009) no. 1–2, pp. 116-135 | Article

[Mueller et al., 2009] D. S. Mueller; C. R. Wagner; M. S. Rehmel; K. A. Oberg; F. Rainville Measuring Discharge with Acoustic Doppler Current Profilers from a Moving Boat, US Department of the Interior, US Geological Survey, 2009

[Muste et al., 2011] M. Muste; H.-C. Ho; D. Kim Considerations on direct stream flow measurements using video imagery: outlook and research needs, J. Hydro-Environ. Res., Volume 5 (2011) no. 4, pp. 289-300 | Article

[Paiva et al., 2011] R. C. D. Paiva; D. C. Buarque; W. Collischonn; M. V Sorribas; D. G. S. da Allasia; C. A. B. Mendes Using TRMM rainfall estimates in hydrologic and hydrodynamic modelling of the Amazon Basin, Grace, Remote Sensing and Ground-based Methods in Multi-scale Hydrology, Volume 343, IASH-Press, Wallingford, Oxfordshire, UK, 2011, pp. 72-77

[Paiva et al., 2013] R. C. D. Paiva; W. Collischonn; M.-P. Bonnet; L. G. G. De Goncalves; S. Calmant; A. Getirana; J. S. Da Silva Assimilating in situ and radar altimetry data into a large-scale hydrologic-hydrodynamic model for streamflow forecast in the Amazon, Hydrol. Earth Syst. Sci., Volume 17 (2013) no. 7, pp. 2929-2946 | Article

[Pan et al., 2016] F. Pan; C. Wang; X. Xi Constructing river stage-discharge rating curves using remotely sensed river cross-sectional inundation areas and river bathymetry, J. Hydrol., Volume 540 (2016), pp. 670-687 | Article

[Park and Latrubesse, 2014] E. Park; E. M. Latrubesse Modeling suspended sediment distribution patterns of the Amazon River using MODIS data, Remote Sens. Environ., Volume 147 (2014) no. 0, pp. 232-242 | Article

[Petersen-Øverleir and Reitan, 2009] A. Petersen-Øverleir; T. Reitan Bayesian analysis of stage–fall–discharge models for gauging stations affected by variable backwater, Hydrol. Process., Volume 23 (2009) no. 21, pp. 3057-3074 | Article

[Peña Arancibia et al., 2015] J. L. Peña-Arancibia; Y. Zhang; D. E. Pagendam; N. R. Viney; J. Lerat; A. I. J. M. van Dijk; J. Vaze; A. J. Frost Streamflow rating uncertainty: characterisation and impacts on model calibration and performance, Environ. Model. Softw., Volume 63 (2015), pp. 32-44 | Article

[Rantz, 1982] S. E. Rantz Measurement and Computation of Streamflow: Volume 2, Computation of Discharge, USGPO, Washington D.C., Washington, USA, 1982

[Rudorff et al., 2014] C. M. Rudorff; J. M. Melack; P. D. Bates Flooding dynamics on the lower Amazon floodplain: 1. Hydraulic controls on water elevation, inundation extent, and river- floodplain discharge, Water Resour. Res., Volume 50 (2014) no. 1, pp. 619-634 | Article

[Stickler et al., 2013] C. M. Stickler; M. T. Coe; M. H. Costa; D. C. Nepstad; D. G. McGrath; L. C. P. Dias; H. O. Rodrigues; B. S. Soares-Filho Dependence of hydropower energy generation on forests in the Amazon Basin at local and regional scales, Proc. Natl Acad. Sci. USA, Volume 110 (2013) no. 23, pp. 9601-9606 | Article

[Trigg et al., 2009] M. A. Trigg; M. D. Wilson; P. D. Bates; M. S. Horritt; D. E. Alsdorf; B. R. Forsberg; M. C. Vega Amazon flood wave hydraulics, J. Hydrol., Volume 374 (2009) no. 1–2, pp. 92-105 | Article

[Trigg et al., 2012] M. A. Trigg; P. D. Bates; M. D. Wilson; G. Schumann; C. Baugh Floodplain channel morphology and networks of the middle Amazon River, Water Resour. Res., Volume 48 (2012) no. 10, p. W10504 | Article

[Tsai, 2005] C. W. Tsai Flood routing in mild-sloped rivers—Wave characteristics and downstream backwater effect, J. Hydrol., Volume 308 (2005) no. 1–4, pp. 151-167 | Article

[Vauchel, 2005] P. Vauchel HYDRACCESS: Software for Management and Processing of Hydro-meteorological Data, Institut de Recherche pour le Développement, Paris, France, 2005

[WMO, 2010] WMO Manual on Stream Gauging (Vols. 1-Field), Chairperson Publications Board, Geneva, Switzerland, 2010

[Yamazaki et al., 2011] D. Yamazaki; S. Kanae; H. Kim; T. Oki A physically based description of floodplain inundation dynamics in a global river routing model, Water Resour. Res., Volume 47 (2011) no. 4, p. W04501 | Article