# Toward continental hydrologic-hydrodynamic modeling in South America

Vinícius A. Siqueira, Rodrigo C. D. Paiva, Ayan S. Fleischmann, Fernando M.

Fan, Anderson L. Ruhoff, Paulo R. M. Pontes, Adrien Paris, Stéphane

Calmant, Walter Collischonn

► To cite this version:

Vinícius A. Siqueira, Rodrigo C. D. Paiva, Ayan S. Fleischmann, Fernando M. Fan, Anderson L. Ruhoff, et al.. Toward continental hydrologic-hydrodynamic modeling in South America. Hydrology and Earth System Sciences, 2018, 22, pp.4815-4842. 10.5194/hess-22-4815-2018 . insu-03670736

HAL Id: insu-03670736

https://insu.hal.science/insu-03670736

Submitted on 17 May 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Toward continental hydrologic–hydrodynamic modeling in South America

Vinícius A. Siqueira¹, Rodrigo C. D. Paiva¹, Ayan S. Fleischmann¹, Fernando M. Fan¹, Anderson L. Ruhoff¹, Paulo R. M. Pontes², Adrien Paris³,⁴,ª, Stéphane Calmant³, and Walter Collischonn¹

¹Instituto de Pesquisas Hidráulicas, Universidade Federal do Rio Grande
do Sul (UFRGS), Porto Alegre, 91501-970, Brazil

²Instituto Tecnológico Vale (ITV), Belém, 66055-090, Brazil

³LEGOS, Université de Toulouse, CNRS, CNES, IRD, UPS, Toulouse, France

⁴GET, Université de Toulouse, UPS, CNRS, IRD, Toulouse, France

ªnow at: Collecte Localisation Satellite (CLS), Ramonville-Saint-Agne, 31520, France

Correspondence: Vinícius A. Siqueira (vinisiquera@gmail.com)

Received: 25 April 2018 – Discussion started: 2 May 2018

Revised: 21 August 2018 – Accepted: 28 August 2018 – Published: 18 September 2018

**Abstract.** Providing reliable estimates of streamflow and hydrological fluxes is a major challenge for water resources management over national and transnational basins in South America. Global hydrological models and land surface models are a possible solution to simulate the terrestrial water cycle at the continental scale, but issues about parameterization and limitations in representing lowland river systems can place constraints on these models to meet local needs. In an attempt to overcome such limitations, we extended a regional, fully coupled hydrologic-hydrodynamic model (MGB; *Modelo hidrológico de Grandes Bacias*) to the continental domain of South America and assessed its performance using daily river discharge, water levels from independent sources (in situ, satellite altimetry), estimates of terrestrial water storage (TWS) and evaporation (ET) from remote sensing and other available global datasets. In addition, river discharge was compared with outputs from global models acquired through the eartH2Observe project (HTESSEL/CaMa-Flood, LIS-FLOOD and WaterGAP3), providing the first cross-scale assessment (regional/continental × global models) that makes use of spatially distributed, daily discharge data. A satisfactory representation of discharge and water levels was obtained (Nash–Sutcliffe efficiency, NSE > 0.6 in 55 % of the cases) and the continental model was able to capture patterns of seasonality and magnitude of TWS and ET, especially over the largest basins of South America. After the comparison

with global models, we found that it is possible to obtain considerable improvement on daily river discharge, even by using current global forcing data, just by combining parameterization and better routing physics based on regional experience. Issues about the potential sources of errors related to both global- and continental-scale modeling are discussed, as well as future directions for improving large-scale model applications in this continent. We hope that our study provides important insights to reduce the gap between global and regional hydrological modeling communities.

## 1 Introduction

Reliable simulations of streamflow dynamics and related processes are vital to support water resources management regarding water security, natural hazards, navigation, agriculture and energy production. Therefore, improved predictions of the hydrological system can aid policymakers and stakeholders in making better decisions, also fostering actions to reduce risk and impacts on water resources under current and future conditions. In South America, recent important floods (e.g., Marengo et al., 2012; Hoyos et al., 2013; Ovando et al., 2016) and droughts (Melo et al., 2016; Erfanian et al., 2017), together with uncertainties about the potential effects of climate change (Marengo et al., 2009), are encouraging new strategies for meeting social, economic and environmen-
tal water demands in large river basins all over the continent, some of them extending beyond political borders.

In this context, large-scale hydrological models arise as important tools for simulating the terrestrial phase of the water cycle. Despite limitations related to observed (in situ) data, especially in developing countries, advances in computational resources and remote sensing technologies are pushing such models toward continental and global scales with increasing resolution (Wood et al., 2011; Bierkens, 2015; Bierkens et al., 2015; Sood and Smathkin, 2015). Currently, estimates of water fluxes at these scales are usually obtained using two modeling frameworks, namely global hydrological models (GHMs) and land surface models (LSMs) (Haddeland et al., 2011). While GHMs are more concerned with water resources assessment and lateral transfer of water, thus enabling quantification of human impacts and water scarcity/stress (e.g., Döll et al., 2009; Wada et al., 2011), LSMs were primarily developed to provide lower boundary conditions for atmospheric circulation models, i.e., having a particular focus on vertical fluxes of heat and water (e.g., Balsamo et al., 2009). The latter are often coupled (i.e., in offline mode) to global river routing models designed for transporting water along drainage networks (e.g., Decharme et al., 2008; Zaitchik et al., 2010; Yamazaki et al., 2011; Getirana et al., 2012), which also enables the conversion from surface and groundwater runoff into river discharge and other surface water variables (e.g., flood extent, water level).

Although global-scale models can provide valuable spatiotemporal estimates of water fluxes and projections of those estimates (Sood and Smakthin, 2015), their ability to reproduce discharge observations at basin scale and to address practical water management issues is still limited (Archfield et al., 2015; Hattermann et al., 2018). Inaccuracies in runoff estimation from GHMs and LSMs may be first attributed to the uncertainty in global satellite precipitation products (Tian and Peters-Lidard, 2010; Sperna Weiland et al., 2015), but several studies have shown considerable differences between model outputs even when using the same meteorological forcing, given the lack of knowledge about runoff generation processes and deficiencies in parameter estimation (e.g., Haddeland et al., 2011; Gudmundsson et al., 2012; Zhou et al., 2012; Beck et al., 2017a). In particular, calibration has been found to have the largest impact on storage fluxes, evapotranspiration and discharge in comparison to variations in model structure and forcing data (Müller Schmied et al., 2014), which is a reason to call for efforts on this exercise as many of the GHMs and LSMs are not calibrated (Sood and Smathkin, 2015; Zhang et al., 2016; Beck et al., 2017a).

An alternative to overcome some limitations of GHMs and LSMs is to expand the spatial domain of hydrological models that were initially developed for catchment to regional scales. Applying these models at national (e.g., Crooks et al., 2014) to continental domains (e.g., Abbaspour et al., 2015; Pechlivanidis and Arheimer, 2015; Donnelly et al., 2016) translates into a better use of local expert knowledge and

country-specific datasets that may be difficult to reach globally. At the same time, it is possible to focus on regionally relevant processes that are usually not included or not well resolved in global models. In South America, for example, several previous studies suggested that lateral water fluxes in large lowland rivers should be resolved using hydrodynamic routing (e.g., Paiva et al., 2011, 2013; Paz et al., 2011, 2014; Yamazaki et al., 2011; Pontes et al., 2017; Zhao et al., 2017), while GHMs generally apply methods based on constant/variable velocity or a kinematic simplification of the St. Venant equations (see the overview by Kauffeldt et al., 2016 and Bierkens, 2015). Even if LSMs can be coupled offline to more physically based global river routing models (e.g., Yamazaki et al., 2011; Getirana et al., 2017b), calibration in the latter is likely to compensate for errors in runoff generation (Pappenberger et al., 2010; Getirana et al., 2013; Hodges, 2013) and lack of relevant vertical hydrological processes linked to river-floodplain dynamics (e.g., Pedinotti et al., 2012; Paz et al., 2014; Fleischmann et al., 2018). In turn, fully coupled large-scale hydrologic-hydrodynamic models (e.g., Paiva et al., 2013) can handle the above interactions while using one single modeling framework, and are now feasible for using in continental domains because recent routing schemes (e.g., Bates et al., 2010) have proved to be computationally efficient for both regional (Getirana et al., 2017b; Pontes et al., 2017; Fleischmann et al., 2018) and global simulations (Yamazaki et al., 2013).

Over the past decades, skill in streamflow prediction has been emphasized in catchment- to regional-scale modeling (Archfield et al., 2015), but there is a growing opportunity to perform further spatial analyses rather than just relying on point measurements. Currently, a wide range of remote sensing products can be used to assess other variables than discharge, such as terrestrial water storage (e.g., Tapley et al., 2004; Watkins et al., 2015), evapotranspiration (e.g., Miralles et al., 2011; Zhang et al., 2018), soil moisture (e.g., Kerr et al., 2012) or water surface elevation derived from satellite altimetry (e.g., Santos da Silva et al., 2010). Previous studies have shown the utility of the aforementioned datasets not only to validate hydrological/routing models (e.g., Alkama et al., 2010; Getirana et al., 2012; Paiva et al., 2013) but also as an interesting tool to constrain and estimate model parameters (Getirana, 2010; Werth and Güntner, 2010; López López et al., 2017). Therefore, remote sensing products can be helpful for continental-scale modeling in assessing regions where streamflow data are scarce, as well as to outline areas in which future model improvements are potentially needed.

In parallel, the interest in building catchment/regional models up to continental domains, together with global models trying to be locally relevant through hyper-resolution (Wood et al., 2011), fosters the need to reduce the gap between these two modeling communities (Bierkens et al., 2015; Archfield et al., 2015). If the primary goal of a continental model is to provide estimates of river discharge to support regional water management practices, the results can
also be compared with outputs of global models to assess both the performance and potential shortcomings of these models under a regional perspective. In recent years, there has been an increasing number of studies assessing outputs of LSMs through multimodel intercomparison (e.g., Zaitchik et al., 2010; Gudmundsson et al., 2012; Zhou et al., 2012), sometimes in conjunction with GHMs (e.g., Haddeland et al., 2011; Beck et al., 2017a). Other recent studies focused on intercomparison of regional models in large basins around the world, usually relating overall performance to a single gauge station and having a particular interest in monthly statistics (e.g., Huang et al., 2016; Eisner et al., 2017; Krysanova et al., 2017). Moreover, little attention has been given to the intercomparison of global- and regional-scale models, and the existing studies only focused on monthly to annual flows (e.g., Zhang et al., 2016) and projection of climate change impacts using a small number of gauge stations (e.g., Gosling et al., 2011; Hattermann et al., 2018). As streamflow is highly variable over space and at short timescales (i.e., daily), model performance should be assessed with spatially distributed data within large basins and at sub-monthly intervals (e.g., Wu et al., 2014; Beck et al., 2017a; Zhao et al., 2017), but to our knowledge, no cross-scale (i.e., regional/continental × global models) intercomparison with a comprehensive evaluation of daily river discharge has been carried out over South America.

In this paper, we aim to start bridging this gap by (i) extending a regional-scale, fully coupled hydrologic–hydrodynamic model to the entirety of South America, assessing its ability to represent discharge and other hydrological variables across the continent; (ii) exploring how discharge estimates from continental-scale modeling perform when compared to those estimated with state-of-the-art global models and (iii) identifying the issues that should be addressed for improving continental-/global-scale modeling in this continent. The next sections provide a brief description of (i) the major river systems of South America, (ii) modeling approaches, (iii) datasets selected for validation, (iv) calibration procedures, (v) global models selected for discharge comparison and (vi) metrics used for assessment of results.

## 2 Overview of the major South American river systems

South America is one of the most freshwater-abundant regions on Earth, contributing around 30% of the global runoff to the oceans (Clark et al., 2015) despite only having 12% of the total land area. Because of a combination of wide latitudinal extent (10° N–55° S), major orographic features and strong oceanic influences (Garcia and Mechoso, 2005; Vera et al., 2006; Garreaud et al., 2009), the continent is subject to a diverse climate that feeds 6 out of the 10 largest basins in the world in terms of mean annual discharge, 4 of them only within the Amazon (Latrubesse et al., 2005). In partic-

ular, the Amazon is probably the most relevant hydrological system of the world, which attracts great scientific attention due to its ecological importance and role in local to global climate (Werth and Avissar, 2002; Vera et al., 2006).

Figure 1 presents general information about South America, including hydroclimatic characteristics, major wetlands and hydrological regions. According to classifications of the Food Agriculture Organization (FAO) and Brazilian National Water Agency (ANA), the continent is partitioned into 27 major hydrological divisions, with the largest ones (Amazon, La Plata and Orinoco basins) sharing much of their water between different countries. Because almost 80% of the territory lies between the tropics, the continent is dominated by an equatorial (tropical) climate that drives high amounts of rainfall especially near the equator, but extreme conditions also exist and range from very arid (such as the Atacama desert in northern Chile) to polar areas over the Andes Cordillera and south of Patagonia (southern Argentina and Chile). As a result, water availability is subject to be highly variable with large discrepancies in mean annual flows, for example, $8.7 \text{ mm yr}^{-1}$ at Desaguadero River in the Colorado basin (Canalejas gauge station, $1.8 \times 10^5 \text{ km}^2$) compared with $2100 \text{ mm yr}^{-1}$ at Japura River (Puerto Cordoba gauge station, $1.5 \times 10^5 \text{ km}^2$), the latter located in the Amazon basin.

Regarding river dynamics, flows in the Amazon are largely affected by floodplains over extensive flat terrains, causing significant flood peak delay and attenuation (Richey et al., 1989; Alsdorf et al., 2007; Yamazaki et al., 2011; Paiva et al., 2013). Because flood waves have travel times in the order of a few months, sometimes out of phase because of the seasonal differences in precipitation (Richey et al., 1989), rivers are subject to strong backwater effects that extend for several hundred kilometers upstream of the river mouth or the confluence of its tributaries (Meade et al., 1991; Getirana and Paiva, 2013; Paiva et al., 2013), existing in both high- and low-water periods (Trigg et al., 2009). To the north, the Orinoco basin shares some characteristics of the Amazon, such as a unimodal flood pulse and low interannual variability of floodplain inundation, especially in the Llanos region (Hamilton et al., 2002). In addition, the high amplitude of the mainstem water level – 14 to 16 m – produces backwater effects in its tributaries with strong hydrological gradients (Rosales et al., 2002).

At the heart of South America, the La Plata basin plays a major role in terms of agriculture, hydroelectricity and the economy in general, corresponding to almost 70% of the combined gross domestic product in countries such as Argentina, Paraguay, Bolivia and Brazil (Vera et al., 2006; Barros et al., 2006). La Plata is mainly composed of the Paraná and Paraguay rivers and, to a lesser extent, the Uruguay River, which are completely different with respect to flow conditions. For example, the upper Paraná is largely regulated by reservoirs (about 50% of mean annual flow, according to Su and Lettenmaier, 2009) and provides around 75%
<figure><img src="image_3.png" /><figcaption>**Figure 1.** South America maps showing (a) countries and major hydrological regions according to FAO and ANA classifications, (b) Major wetlands and lowland regions, adapted from Lehner and Döll (2004), (c) mean annual precipitation derived from the MSWEP dataset (Beck et al., 2017b), (d) Köppen–Geiger updated climate classification from Kottek et al. (2006), (e) relief map based on the Bare-Earth SRTM (O'Loughlin et al., 2016), including main rivers and (f) mean annual flow at discharge gauge stations used in this study.</figcaption></figure>

of discharge up to the confluence with Paraguay River, despite the similar drainage area of both basins (~10⁶ km²) (Barros et al., 2006). The Paraguay basin, on the other hand, is largely influenced by one of the largest wetlands in the world – the Pantanal – a complex anabranching river–floodplain system characterized by very gentle slopes that can be less than 1.5 cm km⁻¹ (Tucci and Clarke, 1998; Berbery and Barros, 2002; Paz et al., 2011; Bravo et al., 2012). Much of the water stored in the floodplain at high water does not return to the channels during the drying phase and become available for evaporation and infiltration

(Paz et al., 2014), so flood waves are lagged by about 4–6 months (Tucci and Clarke, 1998; Hamilton et al., 2002), up to 60% water volume can be lost (Gonçalves et al., 2011), and the shape of downstream hydrographs is strongly modified (Bravo et al., 2012). Moreover, huge areas on the right overbank of the Paraguay River (Chaco Plain) have poorly defined drainage networks associated with alluvial megafans (Latrubesse, 2015), which makes this basin one of the most challenging regions for hydrological modeling in South America.
The headwaters of La Plata also border important hydrological systems such as the Tocantins–Araguaia and São Francisco in a tropical wet–dry biome called Brazilian Cerrado (i.e., Brazilian savanna). The former is composed of Tocantins and Araguaia rivers, flowing parallel northwards until joining approximately 550 km upstream of the basin mouth near the Amazon delta. While the Tocantins is marked by a cascade of large dams, the Araguaia River is much less altered and hosts the huge Bananal plain, which contributes up to 30% of reduction in peak discharge due to floodplain inundation (Lininger and Latrubesse, 2016). Regarding São Francisco, two-thirds of the runoff is generated at the upper part of the basin (Allasia et al., 2006) and the mainstem crosses a semiarid region known as the “drought polygon”, which affects several parts of Northeast Brazil including the Parnaíba basin. At the south of the continent, rivers flowing to the Atlantic Ocean correspond to less than 40% of the area (i.e., > 60% is related to endorheic basins) (Pasquini and Depetris, 2007) and their annual cycles usually show two maxima, one associated with the winter rainfalls and another with snowmelt during spring and early summer (Rivera et al., 2018).

# 3 Methods

## 3.1 MGB model

### 3.1.1 Model description

The MGB, *Modelo hidrológico de Grandes Bacias* (Large-Scale Hydrological Model), is a conceptual, semi-distributed, large-scale hydrological model first presented by Collischonn et al. (2007). The choice of MGB for this study was motivated by several past applications in South America (Allasia et al., 2006), which encompassed rapid response (e.g., Collischonn et al., 2005; Siqueira et al., 2016a) to markedly seasonal and often slow response basins (e.g., Bravo et al., 2012; Paiva et al., 2013; Fan et al., 2016; Pontes et al., 2017). In its most recent version, basins are divided into unit catchments (Paiva et al., 2011; Pontes et al., 2017), each one containing a single river reach with an associated floodplain and hydrological vertical water balance. Combinations of soil type and land use within each unit catchment are categorized as Hydrological Response Units (HRUs). Water and energy budgets are computed independently for each HRU of each unit catchment, and soil is depicted as a bucket model with a single layer. Canopy interception is represented in terms of leaf area index (LAI) and evapotranspiration (soil, plant transpiration and open water evaporation) is calculated using the Penman–Monteith equation. Surface runoff is produced using the variable contribution area concept following the Arno model (Todini, 1996), while groundwater and subsurface flows are computed, respectively, with linear and nonlinear functions according to water availability

in the soil layer. Runoff from each one of the components (surface, subsurface and groundwater) is propagated to the stream network using linear reservoirs (i.e., hillslope routing). As MGB was primarily developed for tropical regions, snow processes are not represented in the current model version.

Flow routing in river channels can be computed using the Muskingum–Cunge method (Collischonn et al., 2007), a one-dimensional full hydrodynamic method (Paiva et al., 2013) or the local inertial method (Pontes et al., 2017). In this work, the MGB was applied with the inertial routing as described by Pontes et al. (2017), which uses the 1-D version of the explicit local inertial approximation proposed by Bates et al. (2010). The routing structure of MGB is similar to that described by Yamazaki et al. (2011, 2013); i.e., the volume of water stored in a given unit catchment is the only prognostic variable, while other variables such as flow depth and flooded area are diagnosed from the stored volume using floodplain profiles derived using sub-grid topography. The floodplain is treated as a simple storage model, and the water level for a given time step is assumed to be constant along the entire unit catchment. In addition, the model accounts for evaporation in floodplains and infiltration from flooded areas to the unsaturated soil (Fleischmann et al., 2018); thus feedbacks between hydrological and hydrodynamic modules can also be represented (i.e., a two-way coupling approach). Further details on model water balance and flow routing equations are presented in the Supplement Sect. S1.

### 3.1.2 GIS processing

All geoprocessing steps were conducted using an adapted version of the IPH-Hydro Tools package (Siqueira et al., 2016b), using the 15 arcsec HydroSHEDS flow direction map (Lehner et al., 2008) as the main input. We chose the latter because it has received extensive corrections to address topological problems in flat areas and endorheic regions (Lehner et al., 2008) and has been successfully applied for river routing in other studies (e.g., Yamazaki et al., 2013; Zhao et al., 2017). An upstream area threshold of 1000 km² was adopted for the onset of drainage networks, while unit catchments and river reaches were delineated using a fixed-length vector-based discretization of $\Delta x = 15$ km (see Sect. S1.3 for details). This length threshold was selected to ensure a balance between model stability and efficiency, resulting in an improved resolution when compared with configurations used by Yamazaki et al. (2013) and Geitirana et al. (2017b) for river routing with the local inertial method at global and Amazon domains (grids with 0.25° resolution), respectively.

To estimate sub-grid floodplain topography, we first computed the Height Above Nearest Drainage (HAND) map (Rennó et al., 2008) using flow directions and drainage networks derived from HydroSHEDS together with the Bare-Earth SRTM v.1 digital elevation model (DEM) (O'Loughlin
et al., 2016), and further, a floodplain profile was created at each unit catchment relating the HAND value, flooded area and water volume, similarly to that done by Yamazaki et al. (2013). The Bare-Earth SRTM, resampled from 3 to 15 arcsec to match the HydroSHEDS resolution, was adopted to account for vegetation biases in floodplains since the C-band radar used by the original SRTM is not able to penetrate fully through the canopy (Carabajal and Hardling, 2005; Berry et al., 2007). Channel bed elevation at a given unit catchment was estimated, subtracting channel bankfull depth from river bank height (i.e., the elevation at bankfull depth) (Yamazaki et al., 2011; Paiva et al., 2013), the latter also derived from Bare-Earth SRTM. However, one of the drawbacks of using an unconditioned DEM such as the Bare-Earth SRTM is the high level of noise affecting channel bank elevations, which need to be attenuated to avoid excessive inundation in low-relief areas. Instead of applying smoothing algorithms that modify the original DEM values (e.g., Paiva et al., 2011), a simple linear regression was fitted to DEM pixels located over drainage networks within each unit catchment (river reach). Channel bank heights were set as the smoothed elevation associated with the center pixel of each river reach, while the original DEM values remained unchanged, for example, when computing the HAND model and associated floodplain profiles (see Sects. S1.5 and S1.6 for more details).

### 3.1.3 River hydraulic geometry

Because flow routing is very sensitive to river geometry (Yamazaki et al., 2011; Getirana et al., 2012; Paiva et al., 2013), channel parameters such as bankfull width and depth must be properly defined. However, detailed information about channel geometry is usually not available for large-scale basins and a very common approach is to adopt classic hydraulic geometry relationships (HGs) (Leopold and Maddock, 1953) for specific sites according to drainage area or discharge (Decharme et al., 2008; Yamazaki et al., 2011; Getirana et al., 2012; Paiva et al., 2013; Luo et al., 2017; Pontes et al., 2017). Here, the global database of Andreadis et al. (2013) was used to set initial values of bankfull widths and depths, which were derived from 2-year return period flows using the Global Runoff Database Center (GRDC) data and universal HGs obtained from several rivers around the world. In addition, regional HGs (Beighley and Gummaldi, 2011; Paiva et al., 2011, 2013; Pontes, 2016) and width estimates based on satellite imagery from Pontes (2016) were included to improve the global channel geometries of Andreadis et al. (2013) for Amazon and La Plata basins.

### 3.1.4 Model forcing

The Multi-Source Weighted Ensemble Precipitation – MSWEP v1.1 (Beck et al., 2017b) – was used as precipitation input to the rainfall-runoff module of the MGB model. This is a 3-hourly, global-scale dataset (0.25° resolution) that optimally combines satellite, reanalysis and daily gauge data, and it has been evaluated with satisfactory results in a recent comparison of several precipitation datasets (Beck et al., 2017c). Regarding climate variables used to compute evapotranspiration, mean monthly data for the period 1961–1990 were retrieved from the Climate Research Unit (CRU) Global Climate v.2 (New et al., 2002), which provides long-term climatologies of temperature, pressure, radiation and wind speed for all land areas at 10' resolution.

### 3.1.5 Land use and soil data

Herein, we used the 400 m resolution HRU merged product (soil + land use) for South America (Fan et al., 2015), which is available at https://www.ufrgs.br/ish (last access: 24 April 2018). Basically, the soil map is a combination of the Brazilian database RADAMBrasil and the FAO Digitized Soil Map of the World and Derived Soil Properties, the latter included to account for areas lying outside of Brazil. Land use classification was retrieved from the Global Land Cover map, which was generated using Envisat MERIS fine-resolution (300 m) satellite imagery over the year 2009. Regional land use maps of some Brazilian states were further included in the HRU merged product to improve level of detail.

## 3.2 Validation datasets

### 3.2.1 Discharge and water level data

Daily records of discharge were collected from several national hydrological services including Agência Nacional de Águas (ANA/Brazil: http://www.snirh.gov.br/hidroweb/, last access: 24 April 2018), Operador Nacional do Serviço Elétrico (ONS/Brazil, Reservoir naturalized flows: http://ons.org.br/, last access: 24 April 2018), Instituto Nacional del Agua (INA/Argentina: http://bdhi.hidricosargentina.gob.ar/, last access: 24 April 2018), Instituto de Hidrología, Meteorología y Estudios Ambientales (IDEAM/Colombia: http://www.ideam.gov.co/solicitud-de-informacion, last access: 24 April 2018), Servicio Nacional de Meteorología y Hidrología (SENAMHI/Peru and Bolivia), Dirección General de Aguas (DGA/Chile: http://snia.dga.cl/BNAConsultas/, last access: 24 April 2018) and other databases such as the Environmental Research Observatory for geodynamical, Hydrological and Biogeochemical control of erosion/alteration and material transport in the Amazon (ORE-HyBam: http://www.ore-hybam.org, last access: 24 April 2018) and the GRDC (http://www.bafg.de/GRDC/, last access: 24 April 2018, for Ecuador and NE South America). Based on expert knowledge, we masked out gauges heavily influenced by up-
stream reservoirs and only included gauges with more than 10 000 km² of drainage area. In a few cases, however, this threshold was lowered to include, at least, a small number of gauges due to the lack of available data. Short time series with less than 5 years of records were also excluded from analysis. The complete list of discharge gauge stations can be found in the Supplement Table S.4.1.

Satellite altimetry data were obtained from the THEIA/Hydroweb website (http://hydroweb.theia-land.fr/, last access: 24 April 2018). Within this database, time series of water surface elevation (WSE) were extracted manually using the methodology presented in Santos da Silva et al. (2010) and are provided at virtual stations (VSs) where the satellite ground track forms a crossover with the river network (~10–40 cm of water level accuracy). A total of 841 VSs were found over the Amazon basin, 10 over the Orinoco basin and 29 over the La Plata basin. Data are derived from the observations of Envisat and Jason-2 for the period of 2002–2010 (35-day repeat orbit) and 2008–2010 (10-day repeat orbit), respectively. In situ stage data from ANA gauge stations were also obtained for the Brazilian territory and were filtered using the same criteria as discharge data.

### 3.2.2 Terrestrial water storage (TWS)

Launched in 2002, the Gravity Recovery and Climate Experiment (GRACE) measures temporal changes in the Earth's gravity field (Tapley et al., 2004). Several studies have shown the ability of GRACE to detect continental water storage variations at large spatial scales (e.g., Wahr et al., 1998; Ramillien et al., 2004; Tapley et al., 2004), which can provide insights into hydrological modeling about potential deficiencies in process description, parameters and input data (Schmidt et al., 2008). Here, we used the Release 05 JPL RL05M v2 mass concentration (mascon) estimates available on the GRACE Tellus website (https://grace.jpl.nasa.gov, last access: 24 April 2018). The JPL RL05M mascon solution solves for monthly gravity anomalies in terms of 3° × 3° equal-area spherical cap mascons, while using a Coastline Resolution Improvement (CRI) filter to discriminate between land and ocean mass portions of each mascon that spans coastlines (Wiese et al., 2016). Regarding the traditional spherical harmonic (SH) approach, which has been widely used in the last decade for global studies (Wahr et al., 1998; Landerer and Swenson, 2012), mascon solutions can also be applied to regional scales (Scanlon et al., 2016) and do not require the user to apply any postprocessing filters to the data, lowering the dependence on using scale factors to recover signal loss (Watkins et al., 2015). Uncertainties in GRACE mascon solutions (3° × 3°) over South America are around 10–15 mm of equivalent water thickness, and have been found to be similar or slightly lower in relation to SH solutions (Scanlon et al., 2016).

Despite the native resolution (3° × 3°), mascon grids are provided with a spatial sampling of 0.5° × 0.5°. We kept the original resolution, computing a simple average of 0.5° grid pixels located inside 3° × 3° mascon locations, as signals at sub-mascon resolution cannot be considered independent of each other. Time series of simulated TWS were first derived by summing water stored in all hydrological compartments, including rivers, floodplains, soil, groundwater and vegetation canopy, in each time step. Similar to Paiva et al. (2013), the modeled TWS was then resampled as the weighted mean of TWS of all unit catchments within each 3° × 3° equal-area mascon cell, using the former drainage area as weight. To ensure agreement with GRACE data, anomalies of simulated TWS were obtained by subtracting the long-term mean computed for the period between 2004 and 2009.

### 3.2.3 Evapotranspiration (ET)

Reference values of ET were extracted from the Climate Data Record (CDR) (Zhang et al., 2018), which is available at http://stream.princeton.edu:8080/opendap/MEaSURES/WC_MULTISOURCES_WB_050/ (last access: 24 April 2018). Within this dataset, 10 gridded global ET products estimated from satellite (five), reanalysis (two) and LSMs (three) were optimally combined at 0.5° resolution using weighted averaging and a Bayesian merging technique. The weight of each product is related to the inverse of the ensemble spread, and the deviation from the ensemble mean is assumed as a proxy of the uncertainty/error in individual products. Together with other optimally merged variables provided by the CDR dataset (precipitation, runoff and TWS), estimates of ET were further adjusted with a Constrained Kalman Filter to ensure terrestrial water budget closure at each 0.5° grid cell (Zhang et al., 2018). For comparison purposes, the modeled ET was spatially aggregated into cells of 0.5° resolution using the unit catchment drainage area as weight.

## 3.3 Model adjustment

Model calibration is commonly performed to improve agreement between observations and model results. However, the traditional gauge-by-gauge calibration used in regional hydrological modeling is not very common in continental to global domains (Archfield et al., 2015; Bierkens, 2015, Samaniego et al., 2017; Mizukami et al. 2017) because it can lead to spatial discontinuities of parameters (i.e., patchwork patterns) and overfitting to account for limitations in data and model structure. In other words, good results of discharge may not reflect a suitable depiction of the underlying hydrological processes, so modelers are more likely to “get the right answers for the wrong reasons” (Kirchner, 2006)

In an attempt to reduce the only dependency of river gauges during calibration, regions of parameter sets were derived by intersecting the global map of lithol-
ogy/geology of Durr et al. (2005) with large South American basins/hydrological regions, the latter shown in Fig. 1a. The Amazon and La Plata basins were further divided into their main tributaries prior to intersection due to their large spatial extent. As parameter sets do not correspond to a single gauge station, but rather to regions defined by geological characteristics, multiple gauges were calibrated at the same time using the same parameter set. It is worth mentioning that calibration still remains a challenge for hydrological modeling with respect to large-scale domains (Mizukami et al., 2017), and assessing the suitability of emerging parameter regionalization techniques (e.g., Samaniego et al., 2017) may be investigated in the future because it is beyond the scope of this paper.

The set of parameters used for model calibration is listed in Table 1, including their respective typical ranges of values. Other parameters used to compute energy balance and evapotranspiration (e.g., LAI, superficial resistance, albedo and canopy height) were defined a priori for each HRU vegetation type according to Collischonn (2001, and references within). Automatic calibration was not used herein to keep coherent values according to soil type and land cover, thus aiding to reduce model overparameterization. Therefore, sensitivity analysis of rainfall-runoff parameters was continuously performed as part of the manual calibration process. Regarding the hydrodynamic module, downstream boundary condition at oceans and lakes (endorheic basins) was set to a constant water level for simplification. Manning coefficient values were globally set to 0.03, with adjustments in specific rivers of the Amazon basin according to Paiva et al. (2013). Infiltration from floodplains to the soil column was considered and only calibrated for the Pantanal region ($K_{inf} = 10 \text{ mm day}^{-1}$) since previous studies showed that vertical hydrological processes largely influence model results in this area (e.g., Paz et al., 2014). It is worth mentioning that model sensitivity to river geometry and infiltration parameters were previously assessed by Paiva et al. (2013) and Fleischmann et al. (2018), respectively.

As a result of the calibration procedure, several model parameter sets were manually adjusted and can be summarized into the following median values and percentile ranges (p5–p95): $W_m = 500 (50–1500)$ mm; $b = 0.2 (0.02–1.5)$; $K_{bas} = 0.2 (0.01–3.0) \text{ mm day}^{-1}$, $K_{int} = 2 (0.1–50) \text{ mm day}^{-1}$, $X_L = 0.67 (0.1–0.67)$, $C_s = 15 (5–35)$, $C_i = 120 (20–200)$ and $C_b = 1200 (800–6000)$ h.

## 3.4 River discharge from GHMs and LSMs

Discharge outputs from state-of-the-art global models were acquired through the eartH2Observe Water Cycle Integrator (WCI; ftp://wci.earth2observe.eu, last access: 24 April 2018). The WCI hosts multidecadal global water resources reanalysis datasets produced by 10 GHMs and LSMs, providing multi-scale (regional, continental and global) estimates of meteorological and hydrological water balance variables. We

selected outputs from the 0.25° resolution Water Resources Reanalysis run 2 (WRR-2) baseline, which is an improved dataset over the WRR-1 (0.5°) produced by the initial project run (Schellekens et al., 2017). Models in the WRR-2 baseline are forced with MSWEP precipitation (1979–2014) and bias-corrected ERA-Interim data using the WFDEI correction methodology (see Dutra et al., 2017). Among the global models in WRR-2, river discharge at 0.25° resolution was only available for one LSM, the HTESSEL offline coupled to CaMa-Flood (Balsamo et al., 2009; Yamazaki et al., 2011), and two GHMs, namely, LISFLOOD (van der Knijf et al., 2010) and WaterGAP3 (Döll et al., 2009). The latter two have some degree of calibration and performed relatively well in terms of runoff in a recent model intercomparison (Beck et al., 2017a), while the former is uncalibrated but uses a state-of-the-art hydrodynamic routing model. Within the eartH2Observe project, LISFLOOD and WaterGAP3 were run at, respectively, 0.1° and 0.08333° resolutions, and discharge was then resampled to 0.25° for WRR-2 (Dutra et al., 2017). A brief overview of the structure of these models is shown in Table S2.1.

Because these models are grid-based, we followed a similar procedure to that in Zhao et al. (2017) to match grid cells to corresponding river gauge stations. First, we applied an automatic routine to find the cell coordinates nearest to the gauge locations. Cells were selected when the difference in the upstream area was within 5%; otherwise, the surrounding cell with minimum upstream area difference was selected. Gauges associated with cells whose drainage area differed more than 15% were excluded from the analysis. This procedure was performed separately for each global model to deal with differences between their respective drainage networks. Moreover, due to the spatial resolution mismatch of LISFLOOD and WaterGAP3, flow accumulation grids were recomputed using their respective flow direction maps (at 0.1° and 0.08333°) and were resampled to the same resolution of discharge grids (0.25°). The corresponding cells were then extensively validated with a thorough, GIS-assisted visual inspection, supported by long-term mean annual discharge grids (derived from each global model) to minimize errors of gauge mislocation.

## 3.5 Metrics for assessment of results

MGB simulation was carried out between 1 January 1990 and 31 December 2009 using a daily time step and a warmup period of 2 years to eliminate the influence of initial conditions. Model results were assessed in terms of discharge, water levels, ET and TWS, while simulated river discharge was further compared with the output of global models. Table 2 lists all efficiency metrics used for the assessment of model results. Statistics such as the Kling–Gupta efficiency and delay index were the same as used in Kling et al. (2012) and Paiva et al. (2013), respectively.
<table><caption>**Table 1.** Parameters used to calibrate the rainfall-runoff module of MGB, including their respective (typical) range of values.</caption><thead><tr><th>Parameter</th><th>Description</th><th>Unit</th><th>Min</th><th>Max</th></tr></thead><tbody><tr><td>W<sub>m</sub></td><td>Maximum water storage</td><td>mm</td><td>50</td><td>2000</td></tr><tr><td>b</td><td>Controls the distribution of water storage capacity of the soil</td><td>–</td><td>0.01</td><td>1.6</td></tr><tr><td>K<sub>bas</sub></td><td>Percolation rate from soil to groundwater</td><td>mm Δt<sup>-1</sup></td><td>0.1</td><td>4</td></tr><tr><td>K<sub>int</sub></td><td>Saturated soil hydraulic conductivity</td><td>mm Δt<sup>-1</sup></td><td>4</td><td>40</td></tr><tr><td>K<sub>inf</sub></td><td>Infiltration rate from floodplains when soil is completely dry</td><td>mm Δt<sup>-1</sup></td><td>–</td><td>–</td></tr><tr><td>XL</td><td>Soil porosity index</td><td>–</td><td>Default = 0.67</td><td>–</td></tr><tr><td>C<sub>B</sub></td><td>Groundwater reservoir residence time</td><td>h</td><td>800</td><td>8000</td></tr><tr><td>C<sub>I</sub></td><td>Adjustment factor for subsurface reservoir residence time</td><td>–</td><td>50</td><td>200</td></tr><tr><td>C<sub>S</sub></td><td>Adjustment factor for superficial reservoir residence time</td><td>–</td><td>1</td><td>30</td></tr></tbody></table>

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Discharge performance over South America in terms of (a) correlation (r), (b) Kling–Gupta efficiency (KGE) and (c) Nash–Sutcliffe (NSE).</figcaption></figure>

# 4 Results and discussion

## 4.1 Model validation

### 4.1.1 River discharge

Simulated daily discharge was compared with in situ observations and results were presented in maps of performance metrics (r, KGE and NSE) at each gauge station (Fig. 2). In addition, the runoff coefficient (RC = Q<sub>mean</sub>/P<sub>mean</sub>) was calculated for each gauge station and was plotted against its respective KGE and drainage area (Fig. 3). There is a good agreement between simulated and observed flows in several regions of South America, as NSE and KGE values are larger than 0.6 in 55% and 70% of the cases, respectively. Model performance is clearly higher in the southern and southeastern regions of Brazil, including central Ama-

zon. On the other hand, performance decreases in regions marked by semiarid to arid climates, such as in Northeast Brazil, west and southwest of the La Plata basin, most parts of Argentina and northern Chile. For instance, a poor correlation (r < 0.2) is observed in a semiarid region covered by the Colorado basin, where snowmelt/glacier melt contribute a large amount to total runoff and corresponds to the main source of water for human activities (Rivera et al., 2017). Other locations with lower performance (NSE < 0.2) usually refer to regions strongly influenced by orography (around Andes Cordillera), which are expected due to larger uncertainties of satellite-derived precipitation in these areas (Tian and Peters-Lidard, 2010; Paiva et al., 2013). River discharge at gauge stations with the RC ranging between 0.3 and 0.6 is generally well represented in all spatial scales, while performance tends to be lower for RC < 0.3 and highly variable for rivers with lower drainage areas.
<table><caption>**Table 2.** Efficiency metrics used in this study.</caption><thead><tr><th>Metric</th><th>Abbreviation</th><th>Assessment</th><th>Variables</th><th>Equation</th></tr></thead><tbody><tr><td>Pearson's correlation coefficient</td><td>r</td><td>Linear correlation</td><td>Discharge, water level, TWS</td><td>$$ \frac{\sum[(x_{sim} - \mu_{sim})(x_{obs} - \mu_{obs})]}{\sqrt{\sum[(x_{sim} - \mu_{sim})^2 \sum[(x_{obs} - \mu_{obs})^2]}} $$</td></tr><tr><td>Modified Kling-Gupta efficiency</td><td>KGE</td><td>Overall performance</td><td>Discharge</td><td>$$ 1 - \sqrt{(1 - r)^2 + \left(1 - \frac{\mu_{sim}}{\mu_{obs}}\right)^2 + \left(1 - \frac{CV_{sim}}{CV_{obs}}\right)^2} $$</td></tr><tr><td>Nash–Sutcliffe efficiency</td><td>NSE</td><td>High flows</td><td>Discharge, water level</td><td>$$ 1 - \sum \frac{(x_{sim} - x_{obs})^2}{(x_{obs} - \mu_{obs})^2} $$</td></tr><tr><td>Log-transformed Nash–Sutcliffe efficiency</td><td>NSElog</td><td>Low flows</td><td>Discharge</td><td>$$ 1 - \sum \frac{[\log(x_{sim}) - \log(x_{obs})]^2}{[\log(x_{obs}) - \log(\mu_{obs})]^2} $$</td></tr><tr><td>Overall BIAS (%)</td><td>Bias</td><td>Under- and overestimation (volume)</td><td>Discharge</td><td>$$ \sum x_{sim} - \sum x_{obs} \frac{(\sigma_{sim} - \sigma_{obs}) \cdot 100}{\sigma_{obs}} $$</td></tr><tr><td>BIAS in standard deviation (%)</td><td>σBIAS</td><td>Under- and overestimation (anomalies)</td><td>Water levels, TWS</td><td>$$ \max_{r_{XY}} [x_{sim} - x_{obs}], \quad \text{lag} [-100; +100] $$</td></tr><tr><td>Delay index (days)</td><td>DI</td><td>Timing errors</td><td>Discharge</td><td></td></tr><tr><td>Root mean square error</td><td>RMSE</td><td>Deviation of predicted values</td><td>ET, TWS</td><td></td></tr></tbody></table>

where $x_{sim}$ is the simulated variable; $x_{obs}$ is the observed variable; $\mu_{sim}$ is the mean of simulated variable; $\mu_{obs}$ is the mean of observed variable; CV is the coefficient of variation, equal to $\sigma/\mu$; $r_{XY}$ is the cross-correlation.

<figure><img src="image_5.png" /><figcaption>**Figure 3.** MGB model performance (KGE) versus RC and drainage area over South America.</figcaption></figure>

Figure 4 shows daily simulated discharge for some of the large South American rivers. The agreement between simulated and observed discharge is notable for both high and low flows in most of the cases, which indicates the model's ability to simulate regional- to continental-scale rivers ($10^5$ km$^2$ to $4.7 \times 10^6$ km$^2$) with different flow regimes. Results in the Amazon basin (e.g., Obidos, $NSE_{HD} = 0.89$) are comparable to other regional studies (e.g., Getirana et al., 2012; Paiva et al., 2013; Luo et al., 2017), while better performance is found over several of its tributaries (e.g., Purus, Madeira and Japura rivers). Figure 4 also highlights the improvements of MGB using a hydrodynamic (HD) model over a non-hydrodynamic (noHD) routing method. In the Paraguay River, peak flows are dramatically reduced at the Amolar gauge station when using the HD routing (up to -75%), and a similar behavior can be seen at Puerto Bermejo ($NSE_{noHD} = -5.8$ to $NSE_{HD} = 0.42$), located about 1600 km downstream near the confluence with the Paraná River. Previous attempts of regional hydrological modeling in this basin that did not account for the floodplain inundation in the Pantanal (e.g., using the calibrated VIC model; see Su and Lettenmaier, 2009) reported negative NSE values for Puerto Bermejo, even at the monthly timescale. Differences in performance between noHD and HD are also quite remarkable (especially in terms of NSE) at gauge stations of Conceicao do Araguaia and Calamar in the lower Magdalena, where a pronounced attenuation effect is observed. On the other hand, in some rivers, such as the Uruguay at Garruchos ($NSE_{noHD} = 0.85$; $NSE_{HD} = 0.82$), Paraná at Itaipu ($NSE_{noHD} = 0.91$; $NSE_{HD} = 0.87$) and Toçantins at Descarreto ($NSE_{noHD} = 0.72$; $NSE_{HD} = 0.70$), the routing method has a minor impact. In the case of Orinoco at Ciudad Bolivar, both hydrographs look similar, but the NSE suggests that results are improved when HD routing is used ($NSE_{noHD} = 0.83$; $NSE_{HD} = 0.9$).
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Comparison between observed (black) and simulated discharge for major South American rivers. Model results are shown considering both hydrodynamic (HD; red color) and non-hydrodynamic (noHD; gray color) routing methods.</figcaption></figure>

### 4.1.2 Water levels

Results of performance metrics regarding water levels are presented in Fig. 5. For a suitable comparison, observed data and modeled WSE were first converted into anomalies (i.e., by subtracting their respective long-term mean: $h_{\text{new}} = h - \bar{h}$) to keep values with the same reference. In addition, Fig. 6 shows time series of simulated water level anomalies (hereafter referred to as water levels) for some of the large rivers

of South America, which were plotted against in situ water levels and satellite altimetry. In general, the results obtained for the assessed gauges and VSs are considered satisfactory in terms of correlation ($r \ge 0.8$ in 80% of cases) and Nash–Sutcliffe (NSE $\ge 0.6$ in 60% of cases), with a reasonable performance for amplitudes ($-30\% < \sigma$ BIAS $< 30\%$ in 50% of the cases). Like prior studies in the same region (e.g., Geitirana et al., 2012, 2017b; Paiva et al., 2013), water levels are well represented in central Amazon, where a good per-
formance is observed at tributaries and along the main river
down to its mouth in the Atlantic Ocean. Amplitudes are
overestimated (σBIAS > 30%) in southeast tributaries such
as Madeira, Xingu and Tapajos, as well as in headwaters lo-
cated in the northwest part of the basin. Outside the Ama-
zon, there are acceptable results in the Orinoco (e.g., lower
Meta), Uruguay and Tocantins–Araguaia basins, where the
model generally performs well in all assessed metrics. Large
overestimations in the standard deviation (σBIAS > +50 %)
are systematically found over the São Francisco main stem,
which are reflected by very low values of NSE (<0.2). On the
Paraguay River, a reasonable agreement between observed
and simulated water levels is observed at Amolar, but perfor-
mance significantly reduces in both correlation and σBIAS
for downstream regions (e.g., at Porto Murtinho). In the lat-
ter case, model results are clearly advanced in time and do
not capture rapid variations of water level originating from
lateral contributions of tributaries.

Performance in water levels is directly related to the agree-
ment between simulated and observed discharge. On the
other hand, even if discharge is well represented, there are
uncertainties related to Manning values and also to river
widths and depths derived from HGs, which do not reflect
singularities of cross sections such as narrowing or widen-
ing of rivers at both gauge and VS locations. Previous stud-
ies have demonstrated the large influence of channel geom-
etry and roughness on both amplitude and timing of water
levels, especially over the Amazon basin (e.g., Yamazaki
et al., 2011; Paiva et al., 2013; Paris et al., 2016; Luo et
al., 2017). Moreover, riverbed profiles are subject to DEM
errors that can hardly be reduced through simple profile-
smoothing procedures. For example, datasets used to remove
the vegetation bias in Bare-Earth SRTM (IceSAT, vegeta-
tion height maps, uncorrected SRTM) have different spatial
resolutions that lead to artifacts around the edges of vegeta-
tion patches (O'Loughlin et al., 2016), producing additional
noise on riverbed elevations. In addition to vegetation, other
SRTM error sources such as stripe noise (Rodriguez et al.,
2006) significantly affect large flat areas on the La Plata basin
(Yamazaki et al., 2017), which can ultimately impact model
results. Model resolution and the ability to route discharge
in downstream multi-directions (e.g., rivers with bifurcations
and anabranching networks) can affect simulated water lev-
els and flooded areas (e.g., Mateo et al., 2017), which has
been taken into account in recent studies with MGB (e.g.,
Pontes et al., 2017; Fleischmann et al., 2018) but not in this
continental model application. Other simplifications in the
model may also affect results, like the assumptions of rect-
angular cross-sections and zero velocity of water stored on
floodplains.

4.1.3 Evapotranspiration (ET)

Figure 7 shows the magnitude and seasonality of ET aver-
aged for major basins in South America, as well as the mag-

nitude of errors (RMSE) comparing modeled values to the
optimal estimate of the CDR dataset. The ratio of the RMSE
to the CDR uncertainty (RMSE<sub>unc</sub>) was calculated to outline
regions where simulated ET tends to deviate from the optimal
CDR value; i.e., values above unity indicate that the model
error is larger than the mean deviation of all datasets (used in
CDR) from their ensemble mean.

Results show that MGB can capture patterns of ET
over the South America region. Simulated ET values are
within the CDR uncertainty range in most of the conti-
nent (RMSE<sub>unc</sub> < 1), with errors varying between 10 and
30 mm month<sup>-1</sup>. A good agreement in terms of magni-
tude and seasonality of ET is observed for the Amazon,
Tocantins–Araguaia and mainly for La Plata basin, but the
model also performs reasonably well in São Francisco. Con-
versely, larger deviations (RMSE<sub>unc</sub> > 1) are found at low
latitudes (20° S–10° N), where RMSE values reach up to
50 mm month<sup>-1</sup>. ET is underestimated during the dry season
in basins such as Orinoco (DJF), Amazon and Tocantins–
Araguaia (JJA), while it is largely overestimated from the
onset to the end of the dry (wet) season in São Francisco
and Parnaíba (Orinoco). The latter two are clearly affected
by a temporal lag in ET seasonality, where simulated values
are delayed by approximately 1 month with respect to CDR
estimations. At midlatitudes ( > 20° S), large RMSE values
(above 50 mm month<sup>-1</sup> and RMSE<sub>unc</sub> > 2) are only observed
in a narrow N–S range over the southern Andes.

Regarding issues about the timing and magnitude of sim-
ulated ET, meteorological forcing probably has an influence
on model performance because long-term mean climate data
are used for ET computations. Another possible reason is the
lack of spatial variability of moisture in the MGB soil col-
umn. Guswa et al. (2002) compared a simple bucket (one
layer) to a physically based Richards model, showing that
large discrepancies can occur with respect to the relation-
ship between ET and average root-zone saturation, as well as
in timing and intensity of transpiration, especially for water-
limited conditions. Moreover, Wang et al. (2006) found that
timescales of evapotranspiration can differ significantly be-
tween a single-layer and multi-layer soil scheme due to non-
linear interactions that occur in the latter. Indeed, ET is ex-
pected to respond quickly at the beginning of the rainy sea-
son due to an increase of water availability at the soil sur-
face layer, which cannot be well represented with a single-
layer, bucket-type model like MGB. In contrast to some of
the datasets used in CDR (e.g., reanalysis and LSMs), the
MGB does not account for snow processes, which may ex-
plain the large RMSE values over the southern Andean re-
gion.

Although it is beyond the scope of this study for a full as-
sessment of ET estimations derived from hydrological mod-
els and other sources, it is important to note that errors pre-
sented here correspond to the difference between both esti-
mates (MGB and CDR). ET is one of the most uncertain wa-
ter balance variables due to its high spatial and temporal vari-
<figure><img src="image_7.png" /><figcaption>**Figure 5.** MGB performance for simulated water levels over South America in terms of (a) correlation (r), (b) Nash–Sutcliffe (NSE) and (c) bias in standard deviation (σBIAS). In situ and satellite altimetry locations are shown by circle and square symbols, respectively.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 6.** Comparison between simulated (blue) and observed (black) water level anomalies in major South American rivers for in situ gauges (continuous lines) and satellite altimetry (circles) at virtual stations (VSs).</figcaption></figure>
<figure><img src="image_9.png" /><figcaption>**Figure 7.** Comparison between MGB and CDR ET estimates in terms of RMSE (a) and seasonality for major South American basins (b). The light gray area represents the proxy of the CDR uncertainty, i.e., the mean deviation of all datasets (within CDR) from the ensemble mean (Zhang et al., 2018).</figcaption></figure>

ability; thus it is difficult to validate given the lack of ground
observations (Miralles et al., 2011; Zhang et al., 2018). Even
accounting for ~70% of the weight in CDR dataset in com-
parison to LSM and reanalysis (Zhang et al., 2018), remote
sensing products of ET are based on Penman–Monteith or
Priestley–Taylor equations that depend on vegetation indices
and meteorological forcing derived from satellite/reanalysis
data, which are associated with many uncertainties (Mi-
ralles et al., 2011; Vinukollu et al., 2011). Christoffersen
et al. (2014) and Maeda et al. (2017) showed that most re-
mote sensing and land surface models are unable to consis-
tently reproduce ET seasonal cycles in tropical areas (across
the Amazon basin) when compared with eddy covariance
measurements and ET estimates from water balance. In the
Amazon, for example, ET seasonality is regulated by radi-
ation, rainfall and how vegetation assimilates water and en-
ergy (Restrepo-Coupe et al., 2013; Maeda et al., 2017). Other
limitations are also associated with the vegetation cover frac-

tion and how ET is partitioned between transpiration, soil and
canopy evaporation (Miralles et al., 2011).

4.1.4 Terrestrial water storage (TWS)

Figure 8 shows the performance of simulated TWS anomalies in comparison to observations from GRACE mascon solutions. To evaluate the ability of MGB to reproduce monthly variations of TWS, both simulated and observed time series of TWS were averaged to the scale of large basins in South America and are presented in Fig. 9.

In general, the results show that MGB has the ability to
represent TWS anomalies over the continent. There is a good
temporal correlation in most parts of tropical South America
(r > 0.75), as well as in temperate regions with dry summer
between latitudes of 30 and 40° S. Amplitudes of TWS are
reasonably well simulated (−20% < σ BIAS < 20%), mainly
in central Brazil, parts of the northeast, south of La Plata and
areas of southern Chile. On the other hand, performance typ-
<figure><img src="image_10.png" /><figcaption>**Figure 8.** Comparison between MGB and GRACE (JPL RL05 v2 mascon solution) TWS anomalies in terms of (a) bias in standard deviation, (b) correlation and (c) RMSE.</figcaption></figure>

ically decreases in semiarid to arid climates, as can be seen in regions such as northern Chile, Colorado basin, west of La Plata and southern Argentina. High negative bias with large RMSE (> 150 mm) is observed in the northeast Amazon and west of the Orinoco, whereas large overestimations are found mainly over coastlines at low latitudes in the Southern Hemisphere (0–20° S). Moreover, in some regions characterized by polar climate, at the extreme south or in areas over the Andes, modeled TWS anomalies are markedly underestimated ($\sigma$BIAS < −80 %).

Modeled TWS is in good agreement with GRACE observations over the Amazon, Tocantins–Araguaia, São Francisco and Parnaíba basins ($r > 0.9$, $|\sigma$BIAS| < 15 % and RMSE < 45 mm), capturing both the interannual variability and amplitude of TWS anomalies for the analyzed period. MGB is also successful in representing annual changes of TWS in the La Plata basin, but with an overestimation ($\sigma$BIAS = 22 %) probably caused by high positive $\sigma$BIAS in the Paraguay and Chaco regions. Errors in the La Plata basin (RMSE = 24 mm) are in the same order as those in the Amazon (RMSE = 26 mm). In addition, larger amplitude differences clearly occur in the Orinoco ($\sigma$BIAS = −32 %), but with a pronounced RMSE (> 60 mm) that mainly originates from the eastern part of the basin.

The good agreement found in the Amazon basin can be attributed to the explicit representation of the surface water reservoir (channels + floodplains), which has been demonstrated to play an important role in both the magnitude and timing of TWS (Alkama et al., 2010; Paiva et al., 2013; Getirana et al., 2017a). Other authors have pointed out that the contribution of surface storage to TWS is also potentially high in the Orinoco (~ 45 %) (e.g., Frappart et al., 2014), suggesting a large underestimation of the soil storage (in the eastern part of the basin) because anomalies of water level were reasonably well simulated. Indeed, surface storage has

been understood as a major component of TWS variability over tropical regions of South America, and may also be relevant for large rivers crossing semiarid areas such as the São Francisco (Getirana et al., 2017a). In the case of La Plata, the TWS amplitude is likely to be amplified if surface water is anticipated in time (Getirana et al., 2017a), which probably occurs due to the low correlation of water levels previously simulated for the Paraguay River. In addition, the absence of a well-defined river system due to very flat terrains (e.g., Chaco region, in the west part of La Plata) potentially favors the dominance of the groundwater dynamics over TWS, as already reported by Kuppel et al. (2015) in the western Pam pas more in the south.

It is worth mentioning that, for regions of South America located in midlatitudes, TWS is dominated by interannual variability rather than the seasonal cycle (Humphrey et al., 2016), where TWS amplitudes are generally lower and errors more apparent. Previous studies showed a strong negative trend in GRACE mascon solutions in the Colorado basin (e.g., Scanlon et al., 2016) that can be associated with a decrease in snow water equivalent over the dry Central Andes (Rivera et al., 2017). Moreover, negative variations in glacier mass have been reported in southern Argentina/Chile over the Patagonia Icefields (Chen et al., 2007), which is probably the main factor responsible for the large RMSE observed in the extreme south of the continent. TWS in nearby semiarid areas is potentially affected by snowmelt/glacier melt because the latter is an important water source of Patagonian rivers flowing to the Atlantic Ocean (Pasquini and Depetris, 2007; Rivera et al., 2018). Inconsistencies along coastlines are also expected because of the smaller size of land masses that increase uncertainty in GRACE estimates (Wiese et al., 2016). Finally, in addition to issues related to model parameterization and depiction of hydrological processes (e.g., snowmelt), artificial reservoirs (dams) and lakes are not in-
<figure><img src="image_11.png" /><figcaption>**Figure 9.** Comparison between MGB and GRACE (JPL RL05 v2 mascon solution) monthly TWS anomalies for major South American basins.</figcaption></figure>

cluded in the current version of the South America MGB model, leading to additional uncertainties in TWS estimation.

## 4.2 Cross-scale comparison of river discharge from continental x global models

This section presents an assessment of MGB-simulated discharge in comparison to the outputs from HTESSEL/CaMa-Flood, LISFLOOD and WaterGAP3 models, extracted from WRR-2 in the context of the eartH2Observe project (Schellekens et al., 2017; Dutra et al., 2017). This offers an interesting opportunity to evaluate to what extent discharge estimates can be improved at the continental scale of South America, as well as to identify the major shortcomings that should be addressed. To provide a concise spatial analysis, discharge from global models was reduced to its ensemble mean (Ensemble GM), and results are presented in terms of the difference of each metric (indicated by “dMetric”), i.e., by subtracting the performance of MGB from the performance of the Ensemble GM (Fig. 10). Bias and DI values are given in terms of absolute differences (d_Abs(metric)) to make both under- and overestimations comparable. Therefore, positive values indicate that MGB outperforms the ensemble mean of global models and vice versa. Detailed performance metrics of each model can be found in the figures in Sect. S3.

The continental model presents improvements for all metrics over most of the South America regions when compared with the global ensemble mean. In relative terms, a better agreement of simulated and observed discharge (d_KGE > 0.8) is observed over semiarid regions (e.g., eastern/northeastern Brazil and most parts of Argentina), which

are strongly impacted by bias in the Ensemble GM (d_Abs(BIAS) > 60%). In tropical regions with marked seasonality and dry winter (e.g., upper Paraná headwaters), differences in bias are lower (–10% < d_Abs(BIAS) < 20%), which indicates that KGE performance mainly depends on the variability of flows that is not captured by the Ensemble GM. Correlation is considerably higher over the Paraguay River (d_r > 0.4), highlighting the strong influence of hydrodynamic effects and complex processes in the Pantanal and Chaco regions, as documented by regional studies (e.g., Paz et al., 2011, 2014; Bravo et al., 2012; Pontes, 2016). There is also a clear correlation improvement in rivers such as the Araguaia, Amazonas and lower Magdalena, which are also affected by river–floodplain interactions with consequent flood peak attenuation (e.g., Paiva et al., 2013; Lininger and Latrubesse, 2016; Angarita et al., 2018; Pontes et al., 2017). A similar performance is observed for timing (d_Abs(DI)) with absolute differences being larger than 20 days, which also occur in the main stem of the Orinoco basin. In terms of NSE, the largest differences in performance previously observed for the KGE now extend to the main Amazon River, to its tributaries in the eastern region (i.e., Tapajos and Xingu) and also to both Magdalena and Tocantins–Araguaia basins, with values of d_NSE ≥ 0.8. With respect to low to medium flows (d_NSElog), there is a similar pattern to d_KGE (except for eastern Brazil), although with more pronounced differences in the Amazon and Magdalena regions.

The Ensemble GM performs relatively well in all statistics over temperate regions with the absence of lowland rivers (e.g., southern Brazil and southern Chile) and outlines specific locations where the continental approach may be somewhat limited. For example, correlation is slightly reduced for
<figure><img src="image_12.png" /><figcaption>**Figure 10.** Difference between performances of MGB and the Ensemble GM for discharge metrics. Values considered as not significant (gray) are within the ranges -0.05 to +0.05 (d_r, d_KGE, d_NSE and d_NSElog), -5% to 5% (d_Abs(BIAS)) and -2 to +2 days (d_Abs(DI)).</figcaption></figure>

the continental model (-0.1 < d_r < -0.2) in areas over the Parnaíba basin and Chile, while a marked decrease in timing performance (d_Abs (DI) < -20 days) is observed in dry Argentinian rivers like Salado (southwest of La Plata basin) and Desaguadero (Colorado basin). Poor estimates of river geometry and large overestimation of flows in these regions may be causing excessive flooded areas and consequent peak attenuation. Regarding intermediate to low flows, considerable differences in model performance are mainly observed over specific rivers in eastern Brazil and parts of the Amazon

basin near the Andes Cordillera (d_NSElog < -0.8), as well as in regions over southern Chile (-0.2 < d_NSElog < -0.6) that are potentially affected by snowmelt.

Table 3 shows differences in median discharge statistics for each global model and also for the Ensemble GM in comparison to the MGB continental model. Because LIS-FLOOD and WaterGAP3 account for reservoir impacts in their model structure, gauge stations with naturalized discharge data (n = 98) were excluded from the analysis to provide a fair assessment. Except for the Ensemble GM, dif-
<figure><img src="image_13.png" /><figcaption>**Figure 11.** Box plots of global (above center line) and continental (below center line) model performances for different South American regions. MGB model configurations: hydrodynamic routing with calibrated rainfall-runoff parameters (MGB_HD Calib), hydrodynamic routing with uncalibrated rainfall-runoff parameters (MGB_HD_noCalib), non-hydrodynamic routing with calibrated rainfall-runoff parameters (MGB_noHD Calib) and non-hydrodynamic routing with uncalibrated rainfall-runoff parameters (MGB_noHD_noCalib). Gauge stations with naturalized flows were removed from the analysis to provide a fair comparison.</figcaption></figure>

ferences in performance regarding each pair of models (i.e., difference between continental and global) are quite similar for KGE (~0.45) and NSElog (~0.5), while being highly variable for both NSE (~1 to ~1.8) and bias (~4% to ~30%). Differences in median DI are between 1 and 2 days, which can be important for cases where flood timing is around this order of magnitude. Among the estimates from global models only, the Ensemble GM outperforms four out of the six metrics analyzed (KGE, NSE, NSElog and DI) with correlation (d_r=0.03) equivalent to the best of global models for this metric (LISFLOOD, d_r=0.02). A reduction in performance only occurs when bias is evaluated, where 50% of the gauge stations have an absolute difference equal to or greater than 11% compared with differences in HTESSEL/CaMa-Flood (d_Abs (BIAS) ≈ 8%) and WaterGAP3 (d_Abs (BIAS) ≈ 4%). In the assessment by Beck

et al. (2017a) for basins < 10 000 km² around the world, LISFLOOD also had an advantage in correlation when compared with other global models, while WaterGAP3 demonstrated problems related to baseflow index, which may be indicated here by the largest difference of NSElog (d_NSElog ≈ 0.6).

The set of box plots shown in Fig. 11 summarizes the individual performance of continental and global models. Results are presented for some of the representative South America basins and for the entire continental region, using a subset of metrics (KGE, NSE, BIAS and DI). In addition, a further analysis of the continental model performance was carried out using a few degraded configurations: hydrodynamic routing with calibrated rainfall-runoff parameters, i.e., the reference simulation (MGB_HD Calib), hydrodynamic routing with uncalibrated rainfall-runoff parameters (MGB_HD_noCalib), non-hydrodynamic routing with cal-
<table><caption>**Table 3.** Median values of discharge metrics for South America, computed as the performance difference between continental and global models. Lower values show better performance for a given global model when benchmarked against the MGB continental hydrologic-hydrodynamic model. Gauge stations with naturalized flows were removed from the analysis to provide a fair comparison.</caption><thead><tr><th>Model difference</th><th>d_r</th><th>d_KGE</th><th>d_NSE</th><th>d_NSElog</th><th>d_Abs(BIAS) (%)</th><th>d_Abs(DI) (days)</th></tr></thead><tbody><tr><td>MGB – HTESSEL/CaMa-Flood</td><td>0.11</td><td>0.48</td><td>1.35</td><td>0.53</td><td>8.3</td><td>1</td></tr><tr><td>MGB – LISFLOOD</td><td>0.02</td><td>0.44</td><td>1.86</td><td>0.48</td><td>32.5</td><td>1.5</td></tr><tr><td>MGB – WaterGAP3</td><td>0.16</td><td>0.44</td><td>1.05</td><td>0.57</td><td>3.8</td><td>2</td></tr><tr><td>MGB – Ensemble GM</td><td>0.03</td><td>0.26</td><td>0.91</td><td>0.27</td><td>11.0</td><td>0</td></tr></tbody></table>

ibrated rainfall-runoff parameters (MGB_noHD Calib) and non-hydrodynamic routing with uncalibrated rainfall-runoff parameters (MGB_noHD_noCalib). For the uncalibrated MGB versions (noCalib), a single set of parameters was adopted corresponding to the median values resulting from model adjustment (as shown in Sect. 3.3). It is important to note that only rainfall-runoff parameters were reduced to their median values, while river routing parameters (Manning coefficient and river geometries) remained unchanged.

Results indicate that global models have important limitations in representing daily discharge in South America. In absolute terms, more than 40% (60%) of the gauge stations show KGE (NSE) values that are negative or close to zero. These models tend to overestimate discharge in the continent, with median bias ranging between +10% and +50%. In general, the performance among global models is variable according to the analyzed region and metric, which is supported by the large box plot ranges. None of the models has a clear advantage with respect to all statistics, and this is especially valid for NSE and KGE. In the Amazon, KGE values present a more uniform pattern than in other regions with a median value close to 0.5, while models agree in a reasonable number of positive KGE and NSE values. Performance in the La Plata basin is highly variable between models, and this is the only region in which both systematic underestimation (HTESSEL/CaMa-Flood with median BIAS ≈ −20%) and overestimation (LISFLOOD with median BIAS ≈ +20%) is observed. Moreover, performance of the global models in basins with semiarid regions (e.g., São Francisco and Parnaíba) is extremely poor for KGE and NSE (median < −1 and < −2 respectively), which is probably associated with a dramatic overestimation of flows in these regions. WaterGAP3 shows a lower bias for all basins, but simulated peak flows occur too early according to DI for South America (percentile 25% of DI ≈ −10 days). On the other hand, LISFLOOD appears to have a systematic delay in flow timing with more pronounced values over the Amazon (median DI ≈ +10 days), and a strong wet bias. For instance, median values of bias in LISFLOOD are larger than 40% for the entire continent and exceed 100% in basins such as São Francisco and Parnaíba. Absolute DI values are generally lower for HTESSEL/CaMa-Flood (median DI closer to 0) and this model usually shows an intermediate performance with respect to other metrics in comparison to LISFLOOD and WaterGAP3.

Furthermore, the Ensemble GM shows a better overall performance when compared with each of these models alone, but still produces a similar number of negative KGE and NSE values (33% and 60% of the gauges, respectively).

Simulated discharge after setting MGB with a single set of rainfall-runoff parameters (MGB_noHD_noCalib) results in positive values for both KGE and NSE median values, varying between 0.3 and 0.6 for the entire continent. The uncalibrated version of the continental model outperforms global models in South America except for basins with semiarid regions (e.g., São Francisco and Parnaíba), where performances seem to be very dependent on parameter adjustment. The introduction of hydrodynamic routing (MGB_HD_noCalib) causes a slight improvement in NSE and KGE but this effect is more evident in the Amazon and especially over the La Plata basin (percentile 25% of NSE changes from −1.5 to 0). Improvements in flow timing (DI) for both Amazon and La Plata are also observed after including the HD routing method, although excessive delays occur in São Francisco and Parnaíba because of the large bias that leads to an excess of floodplain attenuation (see MGB_HD_noCalib × MGB_HD_Calib). Furthermore, box plot ranges are considerably smaller for KGE, NSE and bias in the default MGB simulation (MGB_HD_Calib) with respect to global models, and this reduction can be mostly seen in both MGB calibrated versions.

Our results agree with other studies from the literature, highlighting the large influence of model structure and parameterization in addition to meteorological forcing (e.g., Haddeland et al., 2011; Gudmundsson et al., 2012; Zhou et al., 2012; Beck et al., 2017a). Regarding global models, other studies also stress the large number of negative NSE values resulting from LSMs and GHMs in many basins around the world (e.g., Zhang et al., 2016; Beck et al., 2017a), including South America (e.g., Wu et al., 2014). In particular, WaterGAP3 is expected to produce a lower bias because it is calibrated in terms of mean annual flow (Döll et al., 2009; Müller Schmied et al., 2014), and the systematic advance in timing is probably caused by the simple variable velocity equation (based on Manning) used for computing flow routing. In the case of LISFLOOD, large overestimation of
flows in comparison to other global models has already been reported in the context of the eartH2Observe project (Beck et al., 2017a), where it showed the lowest estimate of potential evaporation (Schellekens et al., 2017). Indeed, this excessive wet bias is one of the possible reasons for the observed delay in flows, but this is not the only factor because large overestimations are concomitantly found in regions where DI is negative. This is the case for eastern Amazon tributaries located downstream of Obidos, such as the Tapajos and Xingu (see figures in Sect. S3). The interplay between the Manning coefficient and groundwater parameters and their influence on flow timing of LISFLOOD has been shown in recent studies (e.g., Revilla-Romero et al., 2015; Zajac et al., 2017), and here may compensate for limitations of the double kinematic wave (channel + floodplain) used for river routing, especially in the Amazon. This suggests that calibrating large basins with lowland river systems using few downstream stations (such as Obidos) should be taken with care if hydrodynamic routing is not accounted for in the model structure.

Although some authors pointed to a clear (and general) underestimation of HTESSEL (e.g., Haddeland et al., 2011; Beck et al., 2017a), our results showed that it only occurs in the La Plata Basin and not in other regions, which may be related to the precipitation forcing used in WRR-2 (MSWEP). With respect to model performances, the relatively better flow timing of HTESSEL can be attributed to CaMa-Flood routing, but the advantages of this coupling were below the expected ones when looking at other statistics. It is worth mentioning that default parameters for river routing were used within WRR-2 simulations; i.e., the CaMa-Flood model was not calibrated (Dutra et al., 2017). This could be one of the reasons that low values of NSE are found over the Amazon main stem for this model (see figures in Sect. S3). Characteristics such as timing and magnitude of flood waves in the hydrodynamic routing are very sensitive to channel geometry and roughness (Yamazaki et al., 2011; Paiva et al., 2013), but also to DEM vegetation effects (Baugh et al., 2013) that can impact the subgrid floodplain profiles (Paiva et al., 2011; Yamazaki et al., 2011). Furthermore, Zhao et al. (2017) emphasize that the benefit of CaMa-Flood highly depends on the runoff fields simulated by the coupled LSM. Our results showed that discharge estimates of an uncalibrated model are improved over Amazon and (mainly) La Plata basins after inclusion of hydrodynamic routing (MGB_HD_noCalib), provided that channel geometry and floodplain topography are reasonably well estimated. Apart from the particular issues of WaterGAP3, LISFLOOD and HTESSEL/CaMa-Flood over South America, our findings reinforce the conclusion of other authors who recommend the ensemble mean of global models as the most reliable estimate (e.g., Haddeland et al., 2011; Gudmundsson et al., 2012; Schellekens et al., 2017; Hattermann et al., 2018), and it occurs even when discharge of a small number of models is averaged (three in the present case).

As outlined by many studies (e.g., Haddeland et al., 2011; Gudmundsson et al., 2012; Pechlivanidis and Arheimer, 2015; Zhang et al., 2016; Krysanova et al., 2017), performance of both regional and global models generally reduces when there is a transition from wet to dry conditions. In semi-arid regions, satellites have several limitations in capturing rainfall intensities due to the local, convective nature of the precipitation, and they often overestimate the occurrence of rainfall because raindrops are likely to evaporate (i.e., subcloud evaporation) before reaching the surface (Dinku et al., 2010; Sunilkumar et al., 2015; Beck et al., 2017c). In addition, runoff generation mechanisms are strongly nonlinear and depend too much on storage processes, which are parameterized with large uncertainty (Gudmundsson et al., 2012). For instance, there is little knowledge about the influence of transmission losses, their partitioning between its main components (e.g., infiltration/evaporation from channels or floodplains) (Jarihani et al., 2015) and the dominant mechanisms of losing/gaining water according to different periods of the wet season (Costa et al., 2013). Processes such as reinfiltration of surface runoff, lateral redistribution of subsurface runoff and hydraulic-connected stream-aquifer interactions have been shown to be necessary for hydrological modeling in Northeast Brazil (Güntner and Bronstert, 2004; Costa et al., 2012, 2013), but are not explicitly accounted for in any structure of the assessed models. Therefore, a systematic underestimation of continental ET and consequent overestimation of flows is expected in dry regions (e.g., Alkama et al., 2010; Haddeland et al., 2011; Zhang et al., 2016). Uncertainties about human interferences in water resources (i.e., small ponds and reservoirs, water abstractions) may also play an important role (Hanasaki et al., 2018), especially in regions where data are scarce. Nevertheless, the complexity of the global models assessed herein makes it difficult to explain the real factors that impact discharge estimates.

Model resolution can partially explain performances in headwater catchments or in areas with complex orography. In addition to issues related to the quality of satellite products (Tian and Peters-Lidard, 2010), aspects that potentially affect these regions are the shape (grid or unit catchment) and size of computational elements, as well as the downscaling method of rainfall fields to force the hydrological models (Rahman et al., 2009), which is in our case a simple inverse distance weighting interpolation. It is worth noting that we have used a drainage area threshold of 1000 km² for headwater catchments (i.e., the onset of drainage networks), while this same area for global models varies between ~100 and ~625 km² according to their respective grid resolutions (between 0.08333° and 0.25°). In addition, results of discharge were evaluated at gauges monitoring at least 10 000 km², meaning that at these points the continental model has at least ~20 unit catchments forced by ~16 MSWEP pixels. On the other hand, response to precipitation at these smaller basins may occur at an hourly scale, while model forcing and analyses were performed at a daily scale.
# 5 Summary and conclusions

For the first time, a regional-scale, fully coupled hydrologic–hydrodynamic model (MGB) was applied to a continental domain (South America). Model results were assessed using observed discharge and water levels from both in situ/satellite altimetry in an unprecedented gauge network over the continent, together with estimates of TWS and ET from remote sensing and other data sources. In addition, a cross-scale assessment (i.e., regional/continental × global models, the latter acquired from the eartH2Observe project) was conducted with the novelty of using spatially distributed, daily discharge data for model comparison.

Regarding continental modeling, analyses showed a satisfactory agreement between simulated and observed discharge, with NSE (KGE) > 0.6 in more than 55% (70%) of the gauges. The performance was generally better in large rivers and humid regions, but worse in areas with semiarid to arid climates, influence of snowmelt or draining complex orography such as the Andes Cordillera. Similar results were found for water levels (both in situ and satellite altimetry), despite 50% of the gauges showing large under- and overestimation of amplitudes (> |30%|). The model was able to capture patterns of seasonality and magnitude of ET and TWS in many parts of the continent, especially when results were averaged to the scale of large South American basins (e.g., Amazon, La Plata, Orinoco, Tocantins–Araguaia, São Francisco, Parnaíba). In addition, model errors in simulating discharge were also found in other hydrological variables, which demonstrates the importance of assessing model results using multiple data sources. Uncertainties were attributed to deficiencies in process representation and simplifications in parameterization, as well as to limitations of the datasets used as model input and validation.

The cross-scale comparison shed light on the extent to which it is possible to improve discharge estimates in South America. Global models (HTESSEL/CaMa-Flood, LISFLOOD and WaterGAP3) presented negative NSE values in a large number of gauges (> 60%) and showed highly variable performances when evaluated over multiple gauges within large basins. A considerable improvement in performance was found when the continental model was compared with individual global models, reaching median differences of around 0.45 for both KGE and NSElog, being larger than unity for NSE. By using the ensemble mean of global models as their best estimate, large differences in absolute bias (> 60%) were detected mainly in eastern/northeastern Brazil and regions over Argentina, as well as in São Francisco, Parnaíba and Magdalena basins. Differences in timing of more than 20 days were found in rivers with floodplain effects, such as the Amazon, La Plata, Tocantins–Araguaia, Orinoco and lower Magdalena. Nevertheless, global models demonstrated a good ability to predict daily discharge over temperate, humid regions with the absence of lowland rivers (e.g.,

southern Brazil and southern Chile), while performing reasonably in the Amazon basin.

The analyses also showed that model calibration and hydrodynamic routing cannot be neglected if simulation of daily discharge in this continent is desired. Calibration was found to be a key factor to model performance in most regions but mainly in drier basins (e.g., Parnaíba and São Francisco), where models generally fail to represent the underlying hydrological processes. In addition, a hydrodynamic routing module was essential to achieve a suitable representation of both magnitude and timing in major river systems, especially in cases where flows are dramatically attenuated by floodplains (e.g., the Paraguay River). However, the expected benefit of coupling hydrodynamic and hydrological models occurs when river geometries are reasonably well represented and calibration of rainfall-runoff parameters is performed together. This must be conducted by looking at different flow signatures (low and high flows, bias, timing of hydrographs) in a spatially distributed way, i.e., not considering only a single downstream gauge of large basins to reduce potential issues related to parameter compensation. As many of the approaches used in this study are applicable to global models, our findings suggest that large improvements on estimated discharge can be achieved by the latter even without a significant increase in the number of computational elements.

Regardless of the scale (global or continental), limitations still remain in some regions of South America and can be explored as “stress tests” in model evaluation studies. For instance, characteristics of the La Plata basin such as complex floodplains, extensive rivers with mild slopes, significant reservoir regulation and existence of several climatic zones make it a unique test bed for model assessments. In addition, large basins with semiarid conditions located in Northeast Brazil (e.g., Parnaíba) or also with snowmelt-driven regimes in the south of the continent (e.g., Colorado) are interesting examples of stress model performance. Despite the cross-scale assessments conducted herein only encompassing river discharge, we recommend that other studies should also include in situ water levels and observation-inferred variables (e.g., TWS from GRACE, ET from remote sensing/global datasets and satellite altimetry) when possible. To facilitate access to in situ data, we call for cooperation among South American countries to produce a continental dataset that can be used by a broader audience, thus contributing to reduce the gap between regional and global modeling communities.

Finally, this study shows that extending a regional, fully coupled hydrologic-hydrodynamic model to the whole South American continent is feasible. This underscores the importance of regional expert knowledge, which can indicate relevant hydrological processes and datasets to be included in continental/global model simulations. We hope that moving from regional toward continental hydrologic-hydrodynamic modeling will bring new opportunities for operational practices such as real-time hydrological forecasting, which is the
topic of ongoing research. Nevertheless, several improve-
ments should be carried out in the model structure, not only
to achieve a better understanding of the underlying processes,
but also to provide further insights about human impacts on
South American water resources. This includes the represen-
tation of reservoirs, lakes and water abstractions. It is also
important to address uncertainties in model parameters, and
these should be investigated further.

*Data availability.* Results from the MGB model are available to the public at http://www.ufrgs.br/lsh (last access: 11 September 2018). All other datasets used in the present study can be accessed using the websites cited in this paper.

**The Supplement related to this article is available**
**online at https://doi.org/10.5194/hess-22-4815-2018-**
**supplement.**

*Author contributions.* VAS worked on the MGB model for South America, designed the study, performed model calibration/simulations, managed all datasets (in situ, remote sensing and global models) and produced the manuscript; RCDP helped with the study design, MGB code optimization and interpretation of results, contributed to writing of the manuscript and conducted a paper review; ASF and FMF assisted with MGB calibration and conducted a paper review; ALR helped with the interpretation of results. PRMP provided support with GIS routines and MGB simulations over the La Plata basin; AP and SC provided support with satellite altimetry data and associated analyses; WC was the research supervisor, contributed to writing of the manuscript and conducted a paper review.

***Competing interests.*** The authors declare that they have no conflict of interest.

***Special issue statement.*** This article is part of the special issue “Integration of Earth observations and models for global water resource assessment”. It is not associated with a conference.

***Acknowledgements.*** The first author would like to acknowledge the Brazilian National Council for Scientific Research (CNPq) for the financial support (Project: “South America Flood Awareness System – SAFAS”, under grant number 422422/2016-9); the institutions ANA, SENAMHIs, IDEAM, DGA, INA, HyBAM, GRDC and ONS for providing in situ discharge and water stage data; the LEGOS-HydroWeb for both Envisat and Jason-2 satellite altimetry data; the Jet Propulsion Laboratory (JPL) for GRACE mascon data; Princeton University for the CDR-ET data; the eartH2Observe project for WRR-2 discharge from global models (WaterGAP3, LISFLOOD and HTESSEL/CaMa-Flood) and, in particular, Martina Floerke, Gabriel Fink, Ad de Roo, Emanuel Dutra and Gianpaolo Balsamo for their kind support with the global models used herein.

Edited by: Albrecht Weerts

Reviewed by: Guy J.-P. Schumann and one anonymous referee

References

Abbaspour, K. C., Rouholahnejad, E., Vaghefi, S., Srinivasan, R., Yang, H., and Kløve, B.: A continental-scale hydrology and water quality model for Europe: Calibration and uncertainty of a high-resolution large-scale SWAT model, J. Hydrol., 524, 733–752, https://doi.org/10.1016/j.jhydrol.2015.03.027, 2015.

Alkama, R., Decharme, B., Douville, H., Becker, M., Cazenave, A., Sheffield, J., Voldoire, A., Tyteca, S., and Le Moigne, P.: Global evaluation of the ISBA-TRIP continental hydrological system. Part I: Comparison to GRACE terrestrial water storage estimates and in situ river discharges, J. Hydrometeorol., 11, 583–600, https://doi.org/10.1175/2010JHM1211.1, 2010.

Allasia, D. G., Da Silva, B. C., Collischonn, W., and Tucci, C. E. M.: Large basin simulation experience in South America, IAHS-AISH Publication, 360–370, 2006.

Alsdorf, D., Bates, P., Melack, J., Wilson, M., and Dunne, T.: Spatial and temporal complexity of the Amazon flood measured from space, Geophys. Res. Lett., 34, L08402, https://doi.org/10.1029/2007GL029447, 2007.

Andreadis, K. M., Schumann, G. J. P., and Pavelsky, T.: A simple global river bankfull width and depth database, Water Resour. Res., 49, 7164–7168, https://doi.org/10.1002/wrcr.20440, 2013.

Angarita, H., Wickel, A. J., Sieber, J., Chavarro, J., Maldonado-Ocampo, J. A., Herrera-R., G. A., Delgado, J., and Purkey, D.: Basin-scale impacts of hydropower development on the Mompós Depression wetlands, Colombia, Hydrol. Earth Syst. Sci., 22, 2839–2865, https://doi.org/10.5194/hess-22-2839-2018, 2018.

Archfield, S. A., Clark, M., Arheimer, B., Hay, L. E., McMillan, H., Kiang, J. E., Seibert, J., Hakala, K., Bock, A., Wagener, T., Farmer, W. H., Andréassian, V., Attinger, S., Viglione, A., Knight, R., Markstrom, S., and Over, T.: Accelerating advances in continental domain hydrologic modeling, Water Resour. Res., 51, 10078–10091, https://doi.org/10.1002/2015wr017498, 2015.

Balsamo, G., Beljaars, A., Scipal, K., Viterbo, P., van den Hurk, B., Hirschi, M., and Betts, A. K.: A Revised Hydrology for the ECMWF Model: Verification from Field Site to Terrestrial Water Storage and Impact in the Integrated Forecast System, J. Hydrometeorol., 10, 623–643, https://doi.org/10.1175/2008jhm1068.1, 2009.

Barros, V., Clarke, R., and Dias, P. S.: Climate change in the La Plata basin, Publication of the Inter-American Institute for Global Change Research (IAI), São José dos Campos, Brazil, 2006.

Bates, P. D., Horritt, M. S., and Fewtrell, T. J.: A simple inertial formulation of the shallow water equations for efficient two-dimensional flood inundation modelling, J. Hydrol., 387, 33–45, https://doi.org/10.1016/j.jhydrol.2010.03.027, 2010.

Baugh, C. A., Bates, P. D., Schumann, G., and Trigg, M. A.: SRTM vegetation removal and hydrodynamic modeling accuracy, Water Resour. Res., 49, 5276–5289, https://doi.org/10.1002/wrcr.20412, 2013.

Beck, H. E., van Dijk, A. I. J. M., de Roo, A., Dutra, E., Fink, G., Orth, R., and Schellekens, J.: Global evaluation of runoff from 10
state-of-the-art hydrological models, Hydrol. Earth Syst. Sci., 21, 2881-2903, https://doi.org/10.5194/hess-21-2881-2017, 2017a.

Beck, H. E., van Dijk, A. I. J. M., Levizzani, V., Schellekens, J., Miralles, D. G., Martens, B., and de Roo, A.: MSWEP: 3-hourly 0.25° global gridded precipitation (1979-2015) by merging gauge, satellite, and reanalysis data, Hydrol. Earth Syst. Sci., 21, 589-615, https://doi.org/10.5194/hess-21-589-2017, 2017b.

Beck, H. E., Vergopolan, N., Pan, M., Levizzani, V., van Dijk, A. I. J. M., Weedon, G. P., Brocca, L., Pappenberger, F., Huffman, G. J., and Wood, E. F.: Global-scale evaluation of 22 precipitation datasets using gauge observations and hydrological modeling, Hydrol. Earth Syst. Sci., 21, 6201-6217, https://doi.org/10.5194/hess-21-6201-2017, 2017c.

Beighley, R. E. and Gummadi, V.: Developing channel and floodplain dimensions with limited data: A case study in the Amazon Basin, Earth Surf. Proc. Land., 36, 1059-1071, https://doi.org/10.1002/esp.2132, 2011.

Berbery, E. H. and Barros, V. R.: The hydrologic cycle of the La Plata basin in South America, J. Hydrometeorol., 3, 630-645, https://doi.org/10.1175/1525-7541(2002)003<0630:THCOTL>2.0.CO;2, 2002.

Berry, P. A. M., Garlick, J. D., and Smith, R. G.: Near-global validation of the SRTM DEM using satellite radar altimetry, Remote Sens. Environ., 106, 17-27, https://doi.org/10.1016/j.rse.2006.07.011, 2007.

Bierkens, M. F. P.: Global hydrology 2015: State, trends, and directions, Water Resour. Res., 51, 4923-4947, https://doi.org/10.1002/2015wr017173, 2015.

Bierkens, M. F. P., Bell, V. A., Burek, P., Chaney, N., Condon, L. E., David, C. H., de Roo, A., Döll, P., Drost, N., Famiglietti, J. S., Flörke, M., Gochis, D. J., Houser, P., Hut, R., Keune, J., Kollet, S., Maxwell, R. M., Reager, J. T., Samaniego, L., Sudicky, E., Sutanudjaja, E. H., van de Giesen, N., Winsemius, H., and Wood, E. F.: Hyper-resolution global hydrological modelling: what is next?, Hydrol. Process., 29, 310-320, https://doi.org/10.1002/hyp.10391, 2015.

Bravo, J. M., Allasia, D., Paz, A. R., Collischonn, W., and Tucci, C. E. M.: Coupled Hydrologic-Hydraulic Modeling of the Upper Paraguay River Basin, J. Hydrol. Eng., 17, 635-646, https://doi.org/10.1061/(ASCE)HE.1943-5584.0000494, 2012.

Carabajal, C. C. and Harding, D. J.: ICESat validation of SRTM C-band digital elevation models, Geophys. Res. Lett., 32, L22S01, https://doi.org/10.1029/2005gl023957, 2005.

Chen, J. L., Wilson, C. R., Tapley, B. D., Blankenship, D. D., and Ivins, E. R.: Patagonia icefield melting observed by gravity recovery and climate experiment (GRACE), Geophys. Res. Lett., 34, L22501, https://doi.org/10.1029/2007gl031871, 2007.

Christoffersen, B. O., Restrepo-Coupe, N., Arain, M. A., Baker, I. T., Cestaro, B. P., Ciais, P., Fisher, J. B., Galbraith, D., Guan, X. D., Gulden, L., van den Hurk, B., Ichii, K., Imbuzeiro, H., Jain, A., Levine, N., Miguez-Machor, G., Poulter, B., Roberti, D. R., Sakaguchi, K., Sahoo, A., Schaefer, K., Shi, M. J., Verbeeck, H., Yang, Z. L., Araujo, A. C., Kruijt, B., Manzi, A. O., da Rocha, H. R., von Randow, C., Muza, M. N., Borak, J., Costa, M. H., de Gonçalves, L. G. G., Zeng, X. B., and Saleska, S. R.: Mechanisms of water supply and vegetation demand govern the seasonality and magnitude of evapotranspiration in Amazonia and Cerrado, Agr. Forest Meteorol., 191, 33-50, https://doi.org/10.1016/j.agrformet.2014.02.008, 2014.

Clark, E. A., Sheffield, J., van Vliet, M. T. H., Nijssen, B., and Lettenmaier, D. P.: Continental Runoff into the Oceans (1950-2008), J. Hydrometeorol., 16, 1502-1520, https://doi.org/10.1175/jhm-d-14-0183.1, 2015.

Collischonn, W.: Simulação Hidrológica de Grandes Bacias, Phd Dissertation, Doutorado em Recursos Hídricos e Saneamento Ambiental, Instituto de Pesquisas Hidráulicas, Universidade Federal do Rio Grande do Sul, Porto Alegre, 270 pp., 2001 (in Portuguese).

Collischonn, W., Haas, R., Andreolli, I., and Tucci, C. E. M.: Forecasting River Uruguay flow using rainfall forecasts from a regional weather-prediction model, J. Hydrol., 305, 87-98, https://doi.org/10.1016/j.jhydrol.2004.08.028, 2005.

Collischonn, W., Allasia, D., Da Silva, B. C., and Tucci, C. E. M.: The MGB-IPH model for large-scale rainfall-runoff modelling, Hydrolog. Sci. J.-Journal Des Sciences Hydrologiques, 52, 878-895, https://doi.org/10.1623/hysj.52.5.878, 2007.

Costa, A. C., Bronstert, A., and de Araújo, J. C.: A channel transmission losses model for different dryland rivers, Hydrol. Earth Syst. Sci., 16, 1111-1135, https://doi.org/10.5194/hess-16-1111-2012, 2012.

Costa, A. C., Foerster, S., de Araújo, J. C., and Bronstert, A.: Analysis of channel transmission losses in a dryland river reach in north-eastern Brazil using streamflow series, groundwater level series and multi-temporal satellite data, Hydrol. Process., 27, 1046-1060, https://doi.org/10.1002/hyp.9243, 2013.

Crooks, S. M., Kay, A. L., Davies, H. N., and Bell, V. A.: From Catchment to National Scale Rainfall-Runoff Modelling: Demonstration of a Hydrological Modelling Framework, Hydrology, 1, 63-88, https://doi.org/10.3390/hydrology1010063, 2014.

Decharme, B., Douville, H., Prigent, C., Papa, F., and Aires, F.: A new river flooding scheme for global climate applications: Offline evaluation over South America, J. Geophys. Res.-Atmos., 113, D11110, https://doi.org/10.1029/2007jd009376, 2008.

Dinku, T., Ceccato, P., Cressman, K., and Connor, S. J.: Evaluating detection skills of satellite rainfall estimates over desert lo-cust recession regions, J. Appl. Meteorol. Clim., 49, 1322-1332, https://doi.org/10.1175/2010JAMC2281.1, 2010.

Döll, P., Fiedler, K., and Zhang, J.: Global-scale analysis of river flow alterations due to water withdrawals and reservoirs, Hydrol. Earth Syst. Sci., 13, 2413-2432, https://doi.org/10.5194/hess-13-2413-2009, 2009.

Donnelly, C., Andersson, J. C. M., and Arheimer, B.: Using flow signatures and catchment similarities to evaluate the E-HYPE multi-basin model across Europe, Hydrolog. Sci. J.-Journal Des Sciences Hydrologiques, 61, 255-273, https://doi.org/10.1080/02626667.2015.1027710, 2016.

Durr, H. H., Meybeck, M., and Durr, S. H.: Lithologic composition of the Earth's continental surfaces derived from a new digital map emphasizing riverine material transfer, Global Biogeochem. Cy., 19, GB4S10, https://doi.org/10.1029/2005gb002515, 2005.

Dutra, E., Balsamo, G., Calvet, J.-C., Munier, S., Burke, S., Fink, G., van Dijk, A., Martinez-de la Torre, A., van Beek, R., de Roo, A., and Polcher, J.: Report on the improved Water Resources Reanalysis. Earth2Observe, Report No.: 5.2., 94 pp., https://doi.org/10.13140/RG.2.2.14523.67369, 2017.

Eisner, S., Florke, M., Chamorro, A., Daggupati, P., Donnelly, C., Huang, J., Hundecha, Y., Koch, H., Kalugin, A., Krylenko, I.,
Mishra, V., Piniewski, M., Samaniego, L., Seidou, O., Wallner, M., and Krysanova, V.: An ensemble analysis of climate change impacts on streamflow seasonality across 11 large river basins, Climatic Change, 141, 401–417, https://doi.org/10.1007/s10584-016-1844-5, 2017.

Erfanian, A., Wang, G., and Fomenko, L.: Unprecedented drought over tropical South America in 2016: Significantly under-predicted by tropical SST, Sci. Rep., 7, 5811, https://doi.org/10.1038/s41598-017-05373-2, 2017.

Fan, F. M., Buarque, D. C., Pontes, P. R. M., and Collischonn, W.: Um mapa de unidades de resposta hidrológica para a América do Sul, XXI Simpósio Brasileiro de Recursos Hídricos, Brasília, PAP019919, 2015.

Fan, F. M., Collischonn, W., Quiroz, K. J., Sorribas, M. V., Buarque, D. C., and Siqueira, V. A.: Flood forecasting on the Tocantins River using ensemble rainfall forecasts and real-time satellite rainfall estimates, J. Flood Risk Manag., 9, 278–288, https://doi.org/10.1111/jfr3.12177, 2016.

Fleischmann, A. S., Siqueira, V. A., Paris, A., Collischonn, W., Paiva, R. C. D., Pontes, P. R. M., Crétaux, J.-F., Bergé-Nguyen, M., Biancamaria, S., Gosset, M., Calmant, S., and Tanimoun, B. A.: Modeling hydrologic and hydrodynamic processes in basins with semi-arid wetlands, J. Hydrol., 561, 943–959, 2018.

Frappart, F., Papa, F., Malbeteau, Y., León, J., Ramillien, G., Prigent, C., Seoane, L., Seyler, F., and Calmant, S.: Surface Freshwater Storage Variations in the Orinoco Floodplains Using Multi-Satellite Observations, Remote Sensing, 7, 89–110, https://doi.org/10.3390/rs70100089, 2014.

Garcia, N. O. and Mechoso, C. R.: Variability in the discharge of South American rivers and in climate, Hydrolog. Sci. J.-Journal Des Sciences Hydrologiques, 50, 459–478, https://doi.org/10.1623/hysj.50.3.459.65030, 2005.

Garreaud, R. D., Vuille, M., Compagnucci, R., and Marengo, J.: Present-day South American climate, Palaeogeogr. Palaeocl., 281, 180–195, https://doi.org/10.1016/j.palaeo.2007.10.032, 2009.

Getirana, A. C. V.: Integrating spatial altimetry data into the automatic calibration of hydrological models, J. Hydrol., 387, 244–255, https://doi.org/10.1016/j.jhydrol.2010.04.013, 2010.

Getirana, A. C. V., Boone, A., Yamazaki, D., Decharme, B., Papa, F., and Mognard, N.: The Hydrological Modeling and Analysis Platform (HyMAP): Evaluation in the Amazon Basin, J. Hydrometeorol., 13, 1641–1665, https://doi.org/10.1175/jhm-d-12-021.1, 2012.

Getirana, A. C. V., Boone, A., Yamazaki, D., and Mognard, N.: Automatic parameterization of a flow routing scheme driven by radar altimetry data: Evaluation in the Amazon basin, Water Resour. Res., 49, 614–629, https://doi.org/10.1002/wrcr.20077, 2013.

Getirana, A. C. V. and Paiva, R. C. D.: Mapping large-scale river flow hydraulics in the Amazon Basin, Water Resour. Res., 49, 2437–2445, https://doi.org/10.1002/wrcr.20212, 2013.

Getirana, A., Kumar, S., Girotto, M., and Rodell, M.: Rivers and Floodplains as Key Components of Global Terrestrial Water Storage Variability, Geophys. Res. Lett., 44, 10359–10368, https://doi.org/10.1002/2017GL074684, 2017a.

Getirana, A., Peters-Lidard, C., Rodell, M., and Bates, P. D.: Trade-off between cost and accuracy in large-scale surface wa-

ter dynamic modeling, Water Resour. Res., 53, 4942–4955, https://doi.org/10.1002/2017wr020519, 2017b.

Gonçalves, H. C., Mercante, M. A., and Santos, E. T.: Hydrological cycle, Brazilian Journal of Biology, 71, 241–253, https://doi.org/10.1590/s1519-69842011000200003, 2011.

Gosling, S. N., Taylor, R. G., Arnell, N. W., and Todd, M. C.: A comparative analysis of projected impacts of climate change on river runoff from global and catchment-scale hydrological models, Hydrol. Earth Syst. Sci., 15, 279–294, https://doi.org/10.5194/hess-15-279-2011, 2011.

Gudmundsson, L., Tallaksen, L. M., Stahl, K., Clark, D. B., Dumont, E., Hagemann, S., Bertrand, N., Gerten, D., Heinke, J., Hanasaki, N., Voss, F., and Koirala, S.: Comparing Large-Scale Hydrological Model Simulations to Observed Runoff Percentiles in Europe, J. Hydrometeorol., 13, 604–620, https://doi.org/10.1175/jhm-d-11-083.1, 2012.

Güntner, A. and Bronstert, A.: Representation of landscape variability and lateral redistribution processes for large-scale hydrological modelling in semi-arid areas, J. Hydrol., 297, 136–161, https://doi.org/10.1016/j.jhydrol.2004.04.008, 2004.

Guswa, A. J., Celia, M. A., and Rodriguez-Iturbe, I.: Models of soil moisture dynamics in ecohydrology: A comparative study, Water Resour. Res., 38, 1166, https://doi.org/10.1029/2001wr000826, 2002.

Haddeland, I., Clark, D. B., Franssen, W., Ludwig, F., Voß, F., Arnell, N. W., Bertrand, N., Best, M., Folwell, S., Gerten, D., Gomes, S., Gosling, S. N., Hagemann, S., Hanasaki, N., Harding, R., Heinke, J., Kabat, P., Koirala, S., Oki, T., Polcher, J., Stacke, T., Viterbo, P., Weedon, G. P., and Yeh, P.: Multimodel Estimate of the Global Terrestrial Water Balance: Setup and First Results, J. Hydrometeorol., 12, 869–884, https://doi.org/10.1175/2011jhm1324.1, 2011.

Hamilton, S. K., Sippel, S. J., and Melack, J. M.: Comparison of inundation patterns among major South American floodplains, J. Geophys. Res.-Atmos., 107, 8038, https://doi.org/10.1029/2000jd000306, 2002.

Hanasaki, N., Yoshikawa, S., Pokhrel, Y., and Kanae, S.: A global hydrological simulation to specify the sources of water used by humans, Hydrol. Earth Syst. Sci., 22, 789–817, https://doi.org/10.5194/hess-22-789-2018, 2018.

Hattermann, F. F., Vetter, T., Breuer, L., Su, B. D., Daggupati, P., Donnelly, C., Fekete, B., Florke, F., Gosling, S. N., Hoffmann, P., Liersch, S., Masaki, Y., Motovilov, Y., Muller, C., Samaniego, L., Stacke, T., Wada, Y., Yang, T., and Krysanova, V.: Sources of uncertainty in hydrological climate impact assessment: a cross-scale study, Environ. Res. Lett., 13, 015006, https://doi.org/10.1088/1748-9326/aa9938, 2018.

Hodges, B. R.: Challenges in Continental River Dynamics, Environ. Modell. Soft., 50, 16–20, https://doi.org/10.1016/j.envsoft.2013.08.010, 2013.

Hoyos, N., Escobar, J., Restrepo, J. C., Arango, A. M., and Ortiz, J. C.: Impact of the 2010–2011 La Nina phenomenon in Colombia, South America: The human toll of an extreme weather event, Appl. Geogr., 39, 16–25, https://doi.org/10.1016/j.apgeog.2012.11.018, 2013.

Huang, S., Kumar, R., Flörke, M., Yang, T., Hundecha, Y., Kraft, P., Gao, C., Gelfan, A., Liersch, S., Lobanova, A., Strauch, M., van Ogtrop, F., Reinhardt, J., Haberlandt, U., and Krysanova, V.: Evaluation of an ensemble of regional hydrological models in 12
large-scale river basins worldwide, Climatic Change, 141, 381–397, https://doi.org/10.1007/s10584-016-1841-8, 2016.

Humphrey, V., Gudmundsson, L., and Seneviratne, S. I.: Assessing Global Water Storage Variability from GRACE: Trends, Seasonal Cycle, Subseasonal Anomalies and Extremes, Surv. Geophys., 37, 357–395, https://doi.org/10.1007/s10712-016-9367-1, 2016.

Jarihani, A. A., Larsen, J. R., Callow, J. N., McVicar, T. R., and Johansen, K.: Where does all the water go? Partitioning water transmission losses in a data-sparse, multi-channel and low-gradient dryland river system using modelling and remote sensing, J. Hydrol., 529, 1511–1529, https://doi.org/10.1016/j.jhydrol.2015.08.030, 2015.

Kauffeldt, A., Wetterhall, F., Pappenberger, F., Salamon, P., and Thielen, J.: Technical review of large-scale hydrological models for implementation in operational flood forecasting schemes on continental level, Environ. Modell. Soft., 75, 68–76, https://doi.org/10.1016/j.envsoft.2015.09.009, 2016.

Kerr, Y. H., Waldteufel, P., Richaume, P., Wigneron, J. P., Ferrazzoli, P., Mahmoodi, A., Al Bitar, A., Cabot, F., Gruhier, C., Julea, S. E., Leroux, D., Mialon, A., and Delwart, S.: The SMOS soil moisture retrieval algorithm, IEEE T. Geosci. Remote, 50, 1384–1403, https://doi.org/10.1109/TGRS.2012.2184548, 2012.

Kirchner, J. W.: Getting the right answers for the right reasons: Linking measurements, analyses, and models to advance the science of hydrology, Water Resour. Res., 42, W03S04, https://doi.org/10.1029/2005wr004362, 2006.

Kling, H., Fuchs, M., and Paulin, M.: Runoff conditions in the upper Danube basin under an ensemble of climate change scenarios, J. Hydrol., 424–425, 264–277, https://doi.org/10.1016/j.jhydrol.2012.01.011, 2012.

Kottek, M., Grieser, J., Beck, C., Rudolf, B., and Rubel, F.: World map of the Köppen-Geiger climate classification updated, Meteorol. Z., 15, 259–263, https://doi.org/10.1127/0941-2948/2006/0130, 2006.

Krysanova, V., Vetter, T., Eisner, S., Huang, S., Pechlivanidis, I., Strauch, M., Gelfan, A., Kumar, R., Aich, V., Arheimer, B., Chamorro, A., Van Griensven, A., Kundu, D., Lobanova, A., Mishra, V., Plötner, S., Reinhardt, J., Seidou, O., Wang, X., Wortmann, M., Zeng, X., and Hattermann, F. F.: Intercomparison of regional-scale hydrological models and climate change impacts projected for 12 large river basins worldwide – A synthesis, Environ. Res. Lett., 12, 105002, https://doi.org/10.1088/1748-9326/aa8359, 2017.

Kuppel, S., Houspanossian, J., Nosetto, M. D., and Jobbagy, E. G.: What does it take to flood the Pampas?: Lessons from a decade of strong hydrological fluctuations, Water Resour. Res., 51, 2937–2950, https://doi.org/10.1002/2015wr016966, 2015.

Landerer, F. W. and Swenson, S. C.: Accuracy of scaled GRACE terrestrial water storage estimates, Water Resour. Res., 48, W04531, https://doi.org/10.1029/2011wr011453, 2012.

Latrubesse, E. M., Stevaux, J. C., and Sinha, R.: Tropical rivers, Geomorphology, 70, 187–206, https://doi.org/10.1016/j geomorph.2005.02.005, 2005.

Latrubesse, E. M.: Large rivers, megafans and other Quaternary avulsive fluvial systems: A potential “who’s who” in the geological record, Earth-Sci. Rev., 146, 1–30, https://doi.org/10.1016/j.earscirev.2015.03.004, 2015.

Lehner, B. and Döll, P.: Development and validation of a global database of lakes, reservoirs and wetlands, J. Hydrol., 296, 1–22, https://doi.org/10.1016/j.jhydrol.2004.03.028, 2004.

Lehner, B., Verdin, K., and Jarvis, A.: New global hydrography derived from spaceborne elevation data, EOS, 89, 93–94, https://doi.org/10.1029/2008EO100001, 2008.

Leopold, L. B. and Maddock, T. J.: The hydraulic geometry of stream channels and some physiographic implications, U.S. Geological Survey Professional Paper, 252, 56 pp., 1953.

Lininger, K. B. and Latrubesse, E. M.: Flooding hydrology and peak discharge attenuation along the middle Araguaia River in central Brazil, Catena, 143, 90–101, https://doi.org/10.1016/j.catena.2016.03.043, 2016.

López López, P., Sutanudjaja, E. H., Schellekens, J., Sterk, G., and Bierkens, M. F. P.: Calibration of a large-scale hydrological model using satellite-based soil moisture and evapotranspiration products, Hydrol. Earth Syst. Sci., 21, 3125–3144, https://doi.org/10.5194/hess-21-3125-2017, 2017.

Luo, X., Li, H.-Y., Leung, L. R., Tesfa, T. K., Getirana, A., Papa, F., and Hess, L. L.: Modeling surface water dynamics in the Amazon Basin using MOSART-Inundation v1.0: impacts of geomorphological parameters and river flow representation, Geosci. Model Dev., 10, 1233–1259, https://doi.org/10.5194/gmd-10-1233-2017, 2017.

Maeda, E. E., Ma, X., Wagner, F. H., Kim, H., Oki, T., Eamus, D., and Huete, A.: Evapotranspiration seasonality across the Amazon Basin, Earth Syst. Dynam., 8, 439–454, https://doi.org/10.5194/esd-8-439-2017, 2017.

Marengo, J. A., Jones, R., Alves, L. M., and Valverde, M. C.: Future change of temperature and precipitation extremes in south america as derived from the precis regional climate modeling system, Int. J. Climatol., 29, 2241–2255, https://doi.org/10.1002/joc.1863, 2009.

Marengo, J. A., Tomasella, J., Soares, W. R., Alves, L. M., and Nobre, C. A.: Extreme climatic events in the Amazon basin, Theor. Appl. Climatol., 107, 73–85, https://doi.org/10.1007/s00704-011-0465-1, 2012.

Mateo, C. M. R., Yamazaki, D., Kim, H., Champathong, A., Vaze, J., and Oki, T.: Impacts of spatial resolution and representation of flow connectivity on large-scale simulation of floods, Hydrol. Earth Syst. Sci., 21, 5143–5163, https://doi.org/10.5194/hess-21-5143-2017, 2017.

Meade, R. H., Rayol, J. M., Daconceicao, S. C., and Natividade, J. R. G.: Backwater effects in the Amazon river basin of Brazil, Environ. Geol. Water S., 18, 105–114, https://doi.org/10.1007/bf01704664, 1991.

Melo, D. D. C. D., Scanlon, B. R., Zhang, Z., Wendland, E., and Yin, L.: Reservoir storage and hydrologic responses to droughts in the Paraná River basin, south-eastern Brazil, Hydrol. Earth Syst. Sci., 20, 4673–4688, https://doi.org/10.5194/hess-20-4673-2016, 2016.

Miralles, D. G., De Jeu, R. A. M., Gash, J. H., Holmes, T. R. H., and Dolman, A. J.: Magnitude and variability of land evaporation and its components at the global scale, Hydrol. Earth Syst. Sci., 15, 967–981, https://doi.org/10.5194/hess-15-967-2011, 2011.

Mizukami, N., Clark, M. P., Newman, A. J., Wood, A. W., Gutmann, E. D., Nijssen, B., Rakovec, O., and Samaniego, L.: Towards seamless large-domain parameter estimation for
hydrologic models, Water Resour. Res., 53, 8020–8040, https://doi.org/10.1002/2017WR020401, 2017.

Müller Schmied, H., Eisner, S., Franz, D., Wattenbach, M., Portmann, F. T., Flörke, M., and Döll, P.: Sensitivity of simulated global-scale freshwater fluxes and storages to in-put data, hydrological model structure, human water use and calibration, Hydrol. Earth Syst. Sci., 18, 3511–3538, https://doi.org/10.5194/hess-18-3511-2014, 2014.

New, M., Lister, D., Hulme, M., and Makin, I.: A high-resolution data set of surface climate over global land areas, Clim. Res., 21, 1–25, 2002.

O'Loughlin, F. E., Paiva, R. C. D., Durand, M., Alsdorf, D. E., and Bates, P. D.: A multi-sensor approach towards a global vegetation corrected SRTM DEM product, Remote Sens. Environ., 182, 49–59, https://doi.org/10.1016/j.rse.2016.04.018, 2016.

Ovando, A., Tomasella, J., Rodriguez, D. A., Martinez, J. M., Siqueira-Junior, J. L., Pinto, G. L. N., Passy, P., Vauchel, P., Nor-iega, L., and von Randow, C.: Extreme flood events in the Boli-vian Amazon wetlands, J. Hydrol.-Regional Studies, 5, 293–308, https://doi.org/10.1016/j.ejrh.2015.11.004, 2016.

Paiva, R. C. D., Collischonn, W., and Tucci, C. E. M.: Large scale hydrologic and hydrodynamic modeling using limited data and a GIS based approach, J. Hydrol., 406, 170–181, https://doi.org/10.1016/j.jhydrol.2011.06.007, 2011.

Paiva, R. C. D., Buarque, D. C., Collischonn, W., Bonnet, M. P., Frappart, F., Calmant, S., and Bulhões Mendes, C. A.: Large-scale hydrologic and hydrodynamic modeling of the Amazon River basin, Water Resour. Res., 49, 1226–1243, https://doi.org/10.1002/wcr.20067, 2013.

Pappenberger, F., Cloke, H. L., Balsamo, G., Ngo-Duc, T., and Oki, T.: Global runoff routing with the hydrological component of the ECMWF NWP system, Int. J. Climatol., 30, 2155–2174, https://doi.org/10.1002/joc.2028, 2010.

Paris, A., Dias de Paiva, R., Santos da Silva, J., Medeiros Moreira, D., Calmant, S., Garambois, P. A., Collischonn, W., Bonnet, M. P., and Seyler, F.: Stage-discharge rating curves based on satellite altimetry and modeled discharge in the Amazon basin, Water Resour. Res., 52, 3787–3814, https://doi.org/10.1002/2014WR016618, 2016.

Pasquini, A. I. and Depetris, P. J.: Discharge trends and flow dynamics of South American rivers draining the southern Atlantic seaboard: An overview, J. Hydrol., 333, 385–399, https://doi.org/10.1016/j.jhydrol.2006.09.005, 2007.

Paz, A. R., Collischonn, W., Tucci, C. E. M., and Padovani, C. R.: Large-scale modelling of channel flow and floodplain inundation dynamics and its application to the Pantanal (Brazil), Hydrol. Process., 25, 1498–1516, https://doi.org/10.1002/hyp.7926, 2011.

Paz, A. R., Collischonn, W., Bravo, J. M., Bates, P. D., and Baugh, C.: The influence of vertical water balance on modelling Pantanal (Brazil) spatio-temporal inundation dynamics, Hydrol. Process., 28, 3539–3553, https://doi.org/10.1002/hyp.9897, 2014.

Pechlivanidis, I. G. and Arheimer, B.: Large-scale hydrologi-cal modelling by using modified PUB recommendations: the India-HYPE case, Hydrol. Earth Syst. Sci., 19, 4559–4579, https://doi.org/10.5194/hess-19-4559-2015, 2015.

Pedinotti, V., Boone, A., Decharme, B., Crétaux, J. F., Mognard, N., Panthou, G., Papa, F., and Tanimoun, B. A.: Evaluation of the ISBA-TRIP continental hydrologic system over the Niger basin

using in situ and satellite derived datasets, Hydrol. Earth Syst. Sci., 16, 1745–1773, https://doi.org/10.5194/hess-16-1745-2012, 2012.

Pontes, P. R. M.: Modelagem hidrológica e hidrodinâmica integrada da bacia do Prata, Phd Dissertation, Instituto de Pesquisas Hidráulicas, Universidade Federal do Rio Grande do Sul, Porto Alegre, 194 pp., 2016 (in Portuguese).

Pontes, P. R. M., Fan, F. M., Fleischmann, A. S., de Paiva, R. C. D., Buarque, D. C., Siqueira, V. A., Jardim, P. F., Sorribas, M. V., and Collischonn, W.: MGB-IPH model for hydrological and hydraulic simulation of large floodplain river systems coupled with open source GIS, Environ. Modell. Softw., 94, 1–20, https://doi.org/10.1016/j.envsoft.2017.03.029, 2017.

Rahman, S., Bagtzoglou, A. C., Hossain, F., Tang, L., Yarbrough, L. D., and Easson, G.: Investigating Spatial Downscaling of Satellite Rainfall Data for Streamflow Simulation in a Medium-Sized Basin, J. Hydrometeorol., 10, 1063–1079, https://doi.org/10.1175/2009JHM1072.1, 2009.

Ramillien, G., Cazenave, A., and Brunau, O.: Global time varia-tions of hydrological signals from GRACE satellite gravimetry, Geophys. J. Int., 158, 813–826, https://doi.org/10.1111/j.1365-246X.2004.02328.x, 2004.

Rennó, C. D., Nobre, A. D., Cuartas, L. A., Soares, J. V., Hod-nett, M. G., Tomasella, J., and Waterloo, M. J.: HAND, a new terrain descriptor using SRTM-DEM: Mapping terra-firme rain-forest environments in Amazonia, Remote Sens. Environ., 112, 3469–3481, https://doi.org/10.1016/j.rse.2008.03.018, 2008.

Restrepo-Coupe, N., da Rocha, H. R., Hutyra, L. R., da Araujo, A. C., Borma, L. S., Christoffersen, B., Cabral, O. M. R., de Ca-margo, P. B., Cardoso, F. L., da Costa, A. C. L., Fitzjarrald, D. R., Goulden, M. L., Kruijrt, B., Maia, J. M. F., Malhi, Y. S., Manzi, A. O., Miller, S. D., Nobre, A. D., von Randow, C., Sa, L. D. A., Sakai, R. K., Tota, J., Wofsy, S. C., Zanchi, F. B., and Saleska, S. R.: What drives the seasonality of photosynthesis across the Amazon basin? A cross-site analysis of eddy flux tower mea-surements from the Brasil flux network, Agr. Forest Meteorol., 182, 128–144, https://doi.org/10.1016/j.agrformet.2013.04.031, 2013.

Revilla-Romero, B., Beck, H. E., Burek, P., Salamon, P., de Roo, A., and Thielen, J.: Filling the gaps: Cali-brating a rainfall-runoff model using satellite-derived sur-face water extent, Remote Sens. Environ., 171, 118–131, https://doi.org/10.1016/j.rse.2015.10.022, 2015.

Richey, J. E., Mertes, L. A. K., Dunne, T., Victoria, R. L., Forsberg, B. R., Tancredi, A., and Oliveira, E.: Sources and routing of the Amazon River flood wave, Global Biogeochem. Cy., 3, 191–204, https://doi.org/10.1029/GB003i003p00191, 1989.

Rivera, J. A., Penalba, O. C., Villalba, R., and Araneo, D. C.: Spatio-temporal patterns of the 2010–2015 extreme hydrologi-cal drought across the Central Andes, Argentina, Water, 9, 652, https://doi.org/10.3390/w9090652, 2017.

Rivera, J. A., Araneo, D. C., Penalba, O. C., and Villalba, R.: Regional aspects of streamflow droughts in the An-dean rivers of Patagonia, Argentina. Links with large-scale climatic oscillations, Hydrol. Res., 49, 134–149, https://doi.org/10.2166/nh.2017.207, 2018.

Rodriguez, E., Morris, C. S. and Belz, J. E.: A global assessment of the SRTM performance, Photogramm. Eng. Rem. S., 72, 249–260, https://doi.org/10.14358/pers.72.3.249, 2006.
Rosales, J., Vispo, C., Dezzeo, N., Blanco-Belmonte, L., Knab-Vispo, C., González, N., Daza, F., Bradley, C., Gilvear, D., and Escalante, G.: Ecohydrology of riparian forests in the Orinoco River Basin, The ecohydrology of South American rivers and wetlands, Special Publication No. 6, IAHS Press, Wallingford, UK, 93–110, 2002.

Samaniego, L., Kumar, R., Thober, S., Rakovec, O., Zink, M., Wanders, N., Eisner, S., Müller Schmied, H., Sutanudjaja, E. H., Warrach-Sagi, K., and Attinger, S.: Toward seamless hydrologic predictions across spatial scales, Hydrol. Earth Syst. Sci., 21, 4323–4346, https://doi.org/10.5194/hess-21-4323-2017, 2017.

Santos da Silva, J., Calmant, S., Seyler, F., Rotunno Filho, O. C., Cochonneau, G., and Mansur, W. J.: Water levels in the Amazon basin derived from the ERS 2 and ENVISAT radar altimetry missions, Remote Sens. Environ., 114, 2160–2181, https://doi.org/10.1016/j.rse.2010.04.020, 2010.

Scanlon, B. R., Zhang, Z., Save, H., Wiese, D. N., Landerer, F. W., Long, D., Longuevergne, L., and Chen, J.: Global evaluation of new GRACE mascon products for hydrologic applications, Water Resour. Res., 52, 9412–9429, https://doi.org/10.1002/2016wr019494, 2016.

Schellekens, J., Dutra, E., Martínez-de la Torre, A., Balsamo, G., van Dijk, A., Sperna Weiland, F., Minvielle, M., Calvet, J.-C., Decharme, B., Eisner, S., Fink, G., Flörke, M., Peßenteiner, S., van Beek, R., Polcher, J., Beck, H., Orth, R., Calton, B., Burke, S., Dorigo, W., and Weedon, G. P.: A global water resources ensemble of hydrological models: the eartH2Observe Tier-1 dataset, Earth Syst. Sci. Data, 9, 389–413, https://doi.org/10.5194/essd-9-389-2017, 2017.

Schmidt, R., Flechtner, F., Meyer, U., Neumayer, K. H., Dahle, C., Konig, R., and Kusche, J.: Hydrological Signals Observed by the GRACE Satellites, Surv. Geophys., 29, 319–334, https://doi.org/10.1007/s10712-008-9033-3, 2008.

Siqueira, V. A., Collischonn, W., Fan, F. M., and Chou, S. C.: Ensemble flood forecasting based on operational forecasts of the regional Eta EPS in the Taquari-Antas basin, RBRH, 21, 587–602, https://doi.org/10.1590/2318-0331.011616004, 2016a.

Siqueira, V. A., Fleischmann, A., Jardim, P. F., Fan, F. M., and Collischonn, W.: IPH-Hydro Tools: A GIS coupled tool for watershed topology acquisition in an open-source environment, RBRH, 21, 274–287, https://doi.org/10.21168/rbrh.v21n1.p274-287, 2016b.

Sood, A. and Smakhtin, V.: Global hydrological models: a review, Hydrol. Sci. J., 60, 549–565, https://doi.org/10.1080/02626667.2014.950580, 2015.

Sperna Weiland, F. C., Vrugt, J. A., van Beek, R. P. H., Weerts, A. H., and Bierkens, M. F. P.: Significant uncertainty in global scale hydrological modeling from precipitation data errors, J. Hydrol., 529, 1095–1115, https://doi.org/10.1016/j.jhydrol.2015.08.061, 2015.

Su, F. G. and Lettenmaier, D. P.: Estimation of the Surface Water Budget of the La Plata Basin, J. Hydrometeorol., 10, 981–998, https://doi.org/10.1175/2009jhm1100.1, 2009.

Sunilkumar, K., Narayana Rao, T., Saikranthi, K., and Purnachandra Rao, M.: Comprehensive evaluation of multisatellite precipitation estimates over India using grided rainfall data, J. Geophys. Res., 120, 8987–9005, https://doi.org/10.1002/2015JD023437, 2015.

Tapley, B. D., Bettadpur, S., Watkins, M., and Reigber, C.: The gravity recovery and climate experiment: Mission overview and early results, Geophys. Res. Lett., 31, L09607, https://doi.org/10.1029/2004gl019920, 2004.

Tian, Y. and Peters-Lidard, C. D.: A global map of uncertainties in satellite-based precipitation measurements, Geophys. Res. Lett., 37, L24407, https://doi.org/10.1029/2010GL046008, 2010.

Todini, E.: The ARNO rainfall-runoff model, J. Hydrol., 175, 339–382, 1996.

Trigg, M. A., Wilson, M. D., Bates, P. D., Horritt, M. S., Alsdorf, D. E., Forsberg, B. R., and Vega, M. C.: Amazon flood wave hydraulics, J. Hydrol., 374, 92–105, https://doi.org/10.1016/j.jhydrol.2009.06.004, 2009.

Tucci, C. E. M. and Clarke, R. T.: Environmental issues in the la Plata Basin, Int. J. Water Resour. D., 14, 157–173, https://doi.org/10.1080/07900629849376, 1998.

Vera, C., Baez, J., Douglas, M., Emmanuel, C. B., Marengo, J., Meitin, J., Nicolini, M., Nogues-Paegle, J., Paegle, J., Penalba, O., Salio, P., Saulo, C., Dias, M. A. S., Dias, P. S., and Zipser, E.: The South American low-level jet experiment, B. Am. Meteorol. Soc., 87, 63–77, https://doi.org/10.1175/bams-87-1-63, 2006.

Vinukollu, R. K., Meynadier, R., Sheffield, J., and Wood, E. F.: Multi-model, multi-sensor estimates of global evapotranspiration: climatology, uncertainties and trends, Hydrol. Process., 25, 3993–4010, https://doi.org/10.1002/hyp.8393, 2011.

Wada, Y., van Beek, L. P. H., and Bierkens, M. F. P.: Modelling global water stress of the recent past: on the relative importance of trends in water demand and climate variability, Hydrol. Earth Syst. Sci., 15, 3785–3808, https://doi.org/10.5194/hess-15-3785-2011, 2011.

Wahr, J., Molenaar, M., and Bryan, F.: Time variability of the Earth's gravity field: Hydrological and oceanic effects and their possible detection using GRACE, J. Geophys. Res.-Sol. Ea., 103, 30205–30229, 1998.

Wang, A., Zeng, X., Shen, S. S. P., Zeng, Q. C., and Dickinson, R. E.: Time scales of land surface hydrology, J. Hydrometeorol., 7, 868–879, https://doi.org/10.1175/JHM527.1, 2006.

Watkins, M. M., Wiese, D. N., Yuan, D. N., Boening, C., and Landerer, F. W.: Improved methods for observing Earth's time variable mass distribution with GRACE using spherical cap mascons, J. Geophys. Res.-Sol. Ea., 120, 2648–2671, https://doi.org/10.1002/2014JB011547, 2015.

Werth, D. and Avissar, R.: The local and global effects of Amazon deforestation, J. Geophys. Res.-Atmos., 107, 8087, https://doi.org/10.1029/2001jd000717, 2002.

Werth, S. and Güntner, A.: Calibration analysis for water storage variability of the global hydrological model WGHM, Hydrol. Earth Syst. Sci., 14, 59–78, https://doi.org/10.5194/hess-14-59-2010, 2010.

Wiese, D. N., Landerer, F. W., and Watkins, M. M.: Quantifying and reducing leakage errors in the JPL RL05M GRACE mascon solution, Water Resour. Res., 52, 7490–7502, https://doi.org/10.1002/2016wr019344, 2016.

Wood, E. F., Roundy, J. K., Troy, T. J., van Beek, L. P. H., Bierkens, M. F. P., Blyth, E., de Roo, A., Doll, P., Ek, M., Famiglietti, J., Gochis, D., van de Giesen, N., Houser, P., Jaffe, P. R., Kollet, S., Lehner, B., Lettenmaier, D. P., Peters-Lidard, C., Sivapalan, M., Sheffield, J., Wade, A., and Whitehead, P.: Hyperresolution global land surface modeling: Meeting a grand challenge
for monitoring Earth's terrestrial water, Water Resour. Res., 47, W05301, https://doi.org/10.1029/2010wr010090, 2011.

Wu, H., Adler, R. F., Tian, Y. D., Huffman, G. J., Li, H. Y., and Wang, J. J.: Real-time global flood estimation using satellite-based precipitation and a coupled land surface and routing model, Water Resour. Res., 50, 2693–2717, https://doi.org/10.1002/2013wr014710, 2014.

Yamazaki, D., Kanae, S., Kim, H., and Oki, T.: A physically based description of floodplain inundation dynamics in a global river routing model, Water Resour. Res., 47, W04501, https://doi.org/10.1029/2010WR009726, 2011.

Yamazaki, D., De Almeida, G. A. M., and Bates, P. D.: Improving computational efficiency in global river models by implementing the local inertial flow equation and a vector-based river network map, Water Resour. Res., 49, 7221–7235, https://doi.org/10.1002/wrcr.20552, 2013.

Yamazaki, D., Ikeshima, D., Tawatari, R., Yamaguchi, T., O'Loughlin, F., Neal, J. C., Sampson, C. C., Kanae, S., and Bates, P. D.: A high-accuracy map of global terrain elevations, Geophys. Res. Lett., 44, 5844–5853, https://doi.org/10.1002/2017GL072874, 2017.

Zaitchik, B. F., Rodell, M., and Olivera, F.: Evaluation of the Global Land Data Assimilation System using global river discharge data and a source-to-sink routing scheme, Water Resour. Res., 46, W06507, https://doi.org/10.1029/2009wr007811, 2010.

Zajac, Z., Revilla-Romero, B., Salamon, P., Burek, P., Hirpa, F., and Beck, H.: The impact of lake and reservoir parameterization on global streamflow simulation, J. Hydrol., 548, 552–568, https://doi.org/10.1016/j.jhydrol.2017.03.022, 2017.

Zhang, Y., Zheng, H., Chiew, F. H. S., Arancibia, J. P., and Zhou, X.: Evaluating Regional and Global Hydrological Models against Streamflow and Evapotranspiration Measurements, J. Hydrometeorol., 17, 995–1010, https://doi.org/10.1175/jhm-d-15-0107.1, 2016.

Zhang, Y., Pan, M., Sheffield, J., Siemann, A. L., Fisher, C. K., Liang, M., Beck, H. E., Wanders, N., MacCracken, R. F., Houser, P. R., Zhou, T., Lettenmaier, D. P., Pinker, R. T., Bytheway, J., Kummerow, C. D., and Wood, E. F.: A Climate Data Record (CDR) for the global terrestrial water budget: 1984–2010, Hydrol. Earth Syst. Sci., 22, 241–263, https://doi.org/10.5194/hess-22-241-2018, 2018.

Zhao, F., Veldkamp, T. I. E., Frieler, K., Schewe, J., Ostberg, S., Willner, S., Schaubberger, B., Gosling, S. N., Schmied, H. M., Portmann, F. T., Leng, G., Huang, M., Liu, X., Tang, Q., Hanasaki, N., Biemans, H., Gerten, D., Satoh, Y., Pokhrel, Y., Stacke, T., Ciais, P., Chang, J., Ducharne, A., Guimberteau, M., Wada, Y., Kim, H., and Yamazaki, D.: The critical role of the routing scheme in simulating peak river discharge in global hydrological models, Environ. Res. Lett., 12, W06507, https://doi.org/10.1088/1748-9326/aa7250, 2017.

Zhou, X., Zhang, Y., Wang, Y., Zhang, H., Vaze, J., Zhang, L., Yang, Y., and Zhou, Y.: Benchmarking global land surface models against the observed mean annual runoff from 150 large basins, J. Hydrol., 470–471, 269–279, https://doi.org/10.1016/j.jhydrol.2012.09.002, 2012.