Rodrigo Cauduro Dias de Paiva

MGB-IPH model for hydrological and hydraulic simulation of large floodplain river systems coupled with open source GIS

Abstract Large-scale hydrological models are useful tools for water resources studies, however, river network flow routing is generally represented using simplified methods, which may lead to simulation errors in flat regions. We present recent improvements to the large-scale hydrological model MGB-IPH to improve its capability of simulating large river basins with extensive floodplains. We also describe the coupling of MGB-IPH to an open source GIS and a large set of developed pre-processing tools with a user-friendly interface for remote sensing data preparation and output visualization. The new features implemented are demonstrated applying the model to the whole Araguaia river basin (380,000 km2). Results are compared to the previous MGB-IPH routing method, observed flow and water level data and remote sensing imagery, showing improvement in the representation of floodplain inundation dynamics. The test case also shows that the proposed model software framework amplifies possibilities of large-scale simulation of ungauged basins.

Improved error estimates of a discharge algorithm for remotely sensed river measurements: Test cases on Sacramento and Garonne Rivers

We present an improvement to a previously presented algorithm that used a Bayesian Markov Chain Monte Carlo method for estimating river discharge from remotely sensed observations of river height, width, and slope. We also present an error budget for discharge calculations from the algorithm. The algorithm may be utilized by the upcoming Surface Water and Ocean Topography (SWOT) mission. We present a detailed evaluation of the method using synthetic SWOT-like observations (i.e., SWOT and AirSWOT, an airborne version of SWOT). The algorithm is evaluated using simulated AirSWOT observations over the Sacramento and Garonne Rivers that have differing hydraulic characteristics. The algorithm is also explored using SWOT observations over the Sacramento River. SWOT and AirSWOT height, width, and slope observations are simulated by corrupting the “true” hydraulic modeling results with instrument error. Algorithm discharge root mean square error (RMSE) was 9% for the Sacramento River and 15% for the Garonne River for the AirSWOT case using expected observation error. The discharge uncertainty calculated from Mannings equation was 16.2% and 17.1%, respectively. For the SWOT scenario, the RMSE and uncertainty of the discharge estimate for the Sacramento River were 15% and 16.2%, respectively. A method based on the Kalman filter to correct errors of discharge estimates was shown to improve algorithm performance. From the error budget, the primary source of uncertainty was the a priori uncertainty of bathymetry and roughness parameters. Sensitivity to measurement errors was found to be a function of river characteristics. For example, Steeper Garonne River is less sensitive to slope errors than the flatter Sacramento River.

The potential impact of new Andean dams on Amazon fluvial ecosystems

Increased energy demand has led to plans for building many new dams in the western Amazon, mostly in the Andean region. Historical data and mechanistic scenarios are used to examine potential impacts above and below six of the largest dams planned for the region, including reductions in downstream sediment and nutrient supplies, changes in downstream flood pulse, changes in upstream and downstream fish yields, reservoir siltation, greenhouse gas emissions and mercury contamination. Together, these six dams are predicted to reduce the supply of sediments, phosphorus and nitrogen from the Andean region by 69, 67 and 57% and to the entire Amazon basin by 64, 51 and 23%, respectively. These large reductions in sediment and nutrient supplies will have major impacts on channel geomorphology, floodplain fertility and aquatic productivity. These effects will be greatest near the dams and extend to the lowland floodplains. Attenuation of the downstream flood pulse is expected to alter the survival, phenology and growth of floodplain vegetation and reduce fish yields below the dams. Reservoir filling times due to siltation are predicted to vary from 106–6240 years, affecting the storage performance of some dams. Total CO2 equivalent carbon emission from 4 Andean dams was expected to average 10 Tg y-1 during the first 30 years of operation, resulting in a MegaWatt weighted Carbon Emission Factor of 0.139 tons C MWhr-1. Mercury contamination in fish and local human populations is expected to increase both above and below the dams creating significant health risks. Reservoir fish yields will compensate some downstream losses, but increased mercury contamination could offset these benefits.

On river‐floodplain interaction and hydrograph skewness

Understanding hydrological processes occurring within a basin by looking at its outlet hydrograph can improve and foster comprehension of ungauged regions. In this context, we present an extensive examination of the roles that floodplains play on driving hydrograph shapes. Observations of many river hydrographs with large floodplain influence are carried out and indicate that a negative skewness of the hydrographs is present among many of them. Through a series of numerical experiments and analytical reasoning, we show how the relationship between flood wave celerity and discharge in such systems is responsible for determining the hydrograph shapes. The more water inundates the floodplains upstream of the observed point, the more negatively skewed is the observed hydrograph. A case study is performed in the Amazon River Basin, where major rivers with large floodplain attenuation (e.g., Purus, Madeira, and Jurua) are identified with higher negative skewness in the respective hydrographs. Finally, different wetland types could be distinguished by using this feature, e.g., wetlands maintained by endogenous processes, from wetlands governed by overbank flow (along river floodplains). A metric of hydrograph skewness was developed to quantify this effect, based on the time derivative of discharge. Together with the skewness concept, it may be used in other studies concerning the relevance of floodplain attenuation in large, ungauged rivers, where remote sensing data (e.g., satellite altimetry) can be very useful.

Comparison of numerical schemes of river flood routing with an inertial approximation of the Saint Venant equations

The one-dimensional flow routing inertial model, formulated as an explicit solution, has advantages over other explicit models used in hydrological models that simplify the Saint-Venant equations. The main advantage is a simple formulation with good results. However, the inertial model is restricted to a small time step to avoid numerical instability. This paper proposes six numerical schemes that modify the one-dimensional inertial model in order to increase the numerical stability of the solution. The proposed numerical schemes were compared to the original scheme in four situations of river’s slope (normal, low, high and very high) and in two situations where the river is subject to downstream effects (dam backwater and tides). The results are discussed in terms of stability, peak flow, processing time, volume conservation error and RMSE (Root Mean Square Error). In general, the schemes showed improvement relative to each type of application. In particular, the numerical scheme here called Prog Q(k+1)xQ(k+1) stood out presenting advantages with greater numerical stability in relation to the original scheme. However, this scheme was not successful in the tide simulation situation. In addition, it was observed that the inclusion of the hydraulic radius calculation without simplification in the numerical schemes improved the results without increasing the computational time.

Hydrological Forecasting Practices in Brazil

Abstract This chapter presents the current hydrological and flood forecasting practices in Brazil, including the main forecast applications, the different kinds of techniques that are currently being employed and the institutions involved in forecast generation. A brief overview of Brazil is provided, including aspects related to its geography, climate, hydrology, and flood hazards. A general discussion about the Brazilian practices on hydrological short- and medium-range forecasting is presented. The chapter also offers detailed examples of some hydrological forecasting systems that are operational or in a research/preoperational phase. Finally, some suggestions are given about how the forecasting practices in Brazil can be understood nowadays, and what the perspectives are for the future.

Uncertainty assessment in hydrodynamic modeling of floods generated by dam break

The objective of the present work was to evaluate sources of uncertainty in dam break simulations for areas with contrasting characteristics: one in a narrow valley with high slopes, and another an open valley with low slopes. A sensitivity analysis of the hydrodynamic model Hec-Ras 5.03 was performed, varying the input data of the model (Manning coefficient, breach configurations, reservoir volume, type of valley topography and equation considered). These variations cause different changes in peak flow, peak time, maximum depth and a maximum speed for different sections of the downstream watercourse. It was concluded that there are uncertainties in determining the input data that impacts in many ways at the generated flood wave, considering both a hydraulic variable of interest, the distance of the section from the dam and the mean geomorphological characteristic of the downstream valley. Topography is not always the most important input data, which allows the possibility of use of low resolution topographies to estimate the peak time in some sections, depending on the valley. Finally, safety coefficients for dam break studies are suggested, aiming to represent uncertainties of the input data in the generated results.

Hydraulic Causes for Basin Hydrograph Skewness

It has been suggested that hydrograph skewness depends on the relative dominance between hillslope and channel network transport processes, where the former ones make the hydrograph positively skewed while the latter tends to make it negatively skewed. More recently, however, the role of river hydraulics in shaping the hydrograph has been highlighted. We present a set of numerical modelling experiments using a hydrodynamic model of river networks in which we investigate how channel hydraulics influence the shape of the hydrograph, particularly its skewness. We further investigate the influence of baseflow, rainfall intensity, cross section geometry and basin scale on unit hydrograph response. We show that river hydraulics has a decisive role because positively skewed hydrographs may occur even when water inputs to the river network are negatively skewed, and in catchments whose width function is also negatively skewed. We show additional results related to the effect of the degree of non-linearity in the relationship between celerity and discharge, the effects of baseflow and rainfall intensity. These further confirm that hydraulic factors may be decisive in determining hydrograph shape.

Efficient number of calibrated cross sections bottom levels on a hydrodynamic model using the SCE-UA algorithm. Case study: Madeira River

Hydrodynamic models are important tools for simulating river water level and flow. A considerable fraction of the hydrodynamic model errors are related to parameters uncertainties. As cross sections bottom levels considerably affect water level simulation, this parameter has to be well estimated for flood studies. Automatic calibration performance and processing time depend on the search space dimension, which is related to the number of calibrated parameters. This paper shows the application of the Shuffled Complex Evolution (SCE-UA) optimization algorithm to assess the number of cross sections bottom levels used in calibration. Also was evaluated the extent of algorithm exploration regarding computational processing time and accuracy. It was tested the calibration of 2, 4, 7 and 10 cross sections bottom levels (2PAR, 4PAR, 7PAR and 10PAR calibration configurations) of a 1,100 km reach of the Madeira River. 7PAR and 10PAR representation had better fitness (lower objective function value) on cross sections used for calibration; however, the error on other cross sections (2 validation gauging stations) was higher than 2PAR and 4PAR calibration. The short number (5) of gauging stations used in calibration has limited the number of calibrated parameters to represent adequately the river level profile. Finally, this paper shows a contribution for the parsimonious selection of parameters regarding the spatial distribution of observation sites used in calibration.