Search

O modelo inercial unidimensional de propagação de vazão, formulado através de uma solução explícita, possui resultados comparáveis aos modelos que consideram as equações completas de Saint-Venant. A principal vantagem é a simplicidade na formulação com bons resultados. No entanto, o modelo inercial está restrito à um passo de tempo pequeno para evitar instabilidade numérica. Este trabalho propõe seis esquemas numéricos que modifcam o modelo inercial unidimensional buscando aumentar a estabilidade numérica da solução. O desempenho dos esquemas numéricos propostos foi comparado em quatro situações de declividade do rio (normal, baixa, alta e muito alta) e em duas situações onde o rio está sujeito ao efeito de jusante (remanso e maré). O modelo hidrodinâmico HEC-RAS foi utilizado para validação das soluções. Os resultados são discutidos em termos de estabilidade, vazão de pico, tempo de processamento, erro de conservação de volume e raiz do erro médio quadrático da vazão (Root Mean Square Error – RMSE). De forma geral os esquemas mostraram melhora relativa a cada tipo de aplicação. Em particular, o esquema numérico aqui denominado Prog Q(k+1)xQ(k+1) se destacou apresentando vantagens com maior estabilidade numérica em relação ao esquema original. Porém este esquema não foi bem sucedido na situação de simulação com maré. Além disso, observou-se que a inclusão do cálculo do raio hidráulico sem simplifcação nos esquemas numéricos melhorou os resultados sem aumentar o tempo computacional. The one-dimensional flow routing inertial model, formulated as an explicit solution, has advantages over other explicit models used in hydrological models that simplify the Saint-Venant equations. The main advantage is a simple formulation with good results. However, the inertial model is restricted to a small time step to avoid numerical instability. This paper proposes six numerical schemes that modify the one-dimensional inertial model in order to increase the numerical stability of the solution. The proposed numerical schemes were compared to the original scheme in four situations of river’s slope (normal, low, high and very high) and in two situations where the river is subject to downstream effects (dam backwater and tides). The results are discussed in terms of stability, peak flow, processing time, volume conservation error and RMSE (Root Mean Square Error). In general, the schemes showed improvement relative to each type of application. In particular, the numerical scheme here called Prog Q(k+1)xQ(k+1) stood out presenting advantages with greater numerical stability in relation to the original scheme. However, this scheme was not successful in the tide simulation situation. In addition, it was observed that the inclusion of the hydraulic radius calculation without simplification in the numerical schemes improved the results without increasing the computational time.

System for hydrological forecasting and alert running in an operational way are important tools for floods impacts reduction. The present study describes the development and results evaluation of an operational discharge forecasting system of the upper Uruguay River basin, sited in Southern Brazil. Developed system was operated every day to provide experimental forecasts with special interest for Barra Grande and Campos Novos hydroelectric power plants reservoirs inflow, with 10 days in advance. We present results of inflow forecasted for floods occurred between July 2013 to July 2016, the period which the system was operated. Forecasts results by visual and performance metrics analysis showed a good fit with observations in most cases, with possibility of floods occurrence being well predicted with antecedence of 2 to 3 days. Comparing the locations, it was noted that the sub-basin of Campos Novos, being slower in rainfall-runoff transformation, is easier forecasted. The difference in predictability between the two basins can be observed by the coefficient of persistence, which is positive from 12h in Barra Grande and from 24h to Campos Novos. These coefficient values also show the value of the rainfall-runoff modeling for forecast horizons of more than one day in the basins. RESUMO Sistemas operacionais de previsão hidrológica e alerta são ferramentas úteis para a redução de impactos das inundações. Este estudo apresenta o desenvolvimento e uma avaliação dos resultados de um sistema de previsão de vazão operacional para a bacia do alto rio Uruguai, localizada no sul do Brasil. O sistema foi operado diariamente para fornecer previsões experimentais com interesse especial nas afluências dos reservatórios das usinas hidrelétricas de Barra Grande e Campos Novos, com antecedência de 10 dias. Aqui são apresentados resultados das previsões de afluências para as principais cheias ocorridas entre julho de 2013 e julho de 2016, período no qual o sistema foi operado. Os resultados das previsões por análise visuais e de métricas de desempenho mostraram um bom ajuste com observações na maioria dos casos, prevendo a possível ocorrência de inundações com 2 a 3 dias de antecedência. Na comparação entre as duas bacias notou-se que um a sub-bacia de Campos Novos, por ser mais lenta na transformação chuva-vazão, é mais fácil de ter suas vazões previstas. A diferença na previsibilidade entre as duas bacias pode ser observada pelo coeficiente de persistência, que é positivo a partir de 12h em Barra Grande e a partir de 24h para Campos Novos. Estes valores de coeficiente também mostram a importância da transformação chuva-vazão para a previsão de horizontes de mais de um dia nas bacias.

RESUMO A delimitação de bacias hidrográficas, geração da rede de drenagem e determinação de características hidráulicas de um rio de interesse são partes importantes de estudos na área de hidrologia. Atualmente muitas dessas informações são obtidas com o processamento de modelos digitais de elevação (MDEs) em softwares comerciais de SIG, como o ArcGIS e o IDRISI. Por outro lado, pacotes de SIG para uso livre, ou seja, gratuitos e de código aberto, têm aumentado significativamente nos últimos anos, e as vantagens desses pacotes incluem ampla distribuição e customização, desenvolvimento continuado pela comunidade de usuários e atendimento a necessidades específicas. Este trabalho apresenta o pacote livre (open-source) denominado IPH-Hydro Tools, um conjunto de ferramentas acoplado ao software livre MapWindow GIS criado para facilitar a aquisição de informações topológicas em bacias hidrográficas, bem como realização de etapas de pré-processamento em modelos hidrológicos a exemplo do MGB-IPH. Para avaliar a aplicabilidade e o desempenho da ferramenta desenvolvida foram realizados testes específicos, através da comparação dos resultados do IPH-Hydro Tools em relação a outros pacotes de SIG (ArcGIS, IDRISI, WhiteBox) disponíveis para esta finalidade. O IPH-Hydro Tools apresentou qualidade de rede de drenagem geralmente superior aos demais pacotes e menor tempo de processamento necessário para delimitação de bacias, apesar de algumas limitações como incompatibilidade em relação a matrizes muito grandes e dificuldade na representação da rede de drenagem em áreas extensas de mesma cota, a exemplo de reservatórios e rios muito largos. ABSTRACT Watershed delineation, drainage network generation and determination of river hydraulic characteristics are important issues in hydrological sciences. In general, this information can be obtained from Digital Elevation Models (DEM) processing within GIS commercial softwares, such as ArcGIS and IDRISI. On the other hand, the use of open source GIS tools has increased significantly, and their advantages include free distribution, continuous development by user communities and full customization for specific requirements. Herein, we present the IPH-Hydro Tools, an open source tool coupled to MapWindow GIS software designed for watershed topology acquisition, including preprocessing steps in hydrological models such as MGB-IPH. In addition, several tests were carried out assessing the performance and applicability of the developed tool, given by a comparison with available GIS packages (ArcGIS, IDRISI, WhiteBox) for similar purposes. The IPH-Hydro Tools provided satisfactory results on tested applications, allowing for better drainage network and less processing time for catchment delineation. Regarding its limitations, the developed tool was incompatible with huge terrain data and showed some difficulties to represent drainage networks in extensive flat areas, which can occur in reservoirs and large rivers.

International audience; In this study, rating curves (RCs) were determined by applying satellite altimetry to a poorly gauged basin. This study demonstrates the synergistic application of remote sensing and watershed modeling to capture the dynamics and quantity of flow in the Amazon River Basin, respectively. Three major advancements for estimating basin‐scale patterns in river discharge are described. The first advancement is the preservation of the hydrological meanings of the parameters expressed by Manning's equation to obtain a data set containing the elevations of the river beds throughout the basin. The second advancement is the provision of parameter uncertainties and, therefore, the uncertainties in the rated discharge. The third advancement concerns estimating the discharge while considering backwater effects. We analyzed the Amazon Basin using nearly one thousand series that were obtained from ENVISAT and Jason‐2 altimetry for more than 100 tributaries. Discharge values and related uncertainties were obtained from the rain‐discharge MGB‐IPH model. We used a global optimization algorithm based on the Monte Carlo Markov Chain and Bayesian framework to determine the rating curves. The data were randomly allocated into 80% calibration and 20% validation subsets. A comparison with the validation samples produced a Nash‐Sutcliffe efficiency ( urn:x-wiley:00431397:media:wrcr22024:wrcr22024-math-0001) of 0.68. When the MGB discharge uncertainties were less than 5%, the urn:x-wiley:00431397:media:wrcr22024:wrcr22024-math-0002 value increased to 0.81 (mean). A comparison with the in situ discharge resulted in an urn:x-wiley:00431397:media:wrcr22024:wrcr22024-math-0003 value of 0.71 for the validation samples (and 0.77 for calibration). The urn:x-wiley:00431397:media:wrcr22024:wrcr22024-math-0004 values at the mouths of the rivers that experienced backwater effects significantly improved when the mean monthly slope was included in the RC. Our RCs were not mission‐dependent, and the urn:x-wiley:00431397:media:wrcr22024:wrcr22024-math-0005 value was preserved when applying ENVISAT rating curves to Jason‐2 altimetry at crossovers. The cease‐to‐flow parameter of our RCs provided a good proxy for determining river bed elevation. This proxy was validated against Acoustic Doppler current profiler (ADCP) cross sections with an accuracy of more than 90%. Altimetry measurements are routinely delivered within a few days, and this RC data set provides a simple and cost‐effective tool for predicting discharge throughout the basin in nearly real time.

Modelos hidrológicos são capazes de representar processos do ciclo hidrológico através de equações matemáticas. Dentre esses processos, a propagação de vazão nos rios é de grande importância, especialmente em bacias de grande escala. Em geral, o módulo de propagação de vazão implementado nos modelos hidrológicos é simplificado, sendo adotada uma versão da metodologia de onda cinemática, na maior parte dos casos. Apesar disso, modelos que utilizam essas metodologias produzem resultados aceitáveis na maioria das aplicações. Entretanto, bacias com rios de baixa declividade, com trechos sujeitos a remanso, ou com extensas planícies de inundação, podem não ser bem representados por metodologias de propagação simplificadas. Uma alternativa é o uso de modelos hidrodinâmicos para propagação de vazão. Entretanto, esses modelos necessitam de mais dados de entrada e possuem equacionamentos complexos. Uma opção de complexi- dade intermediária que pode ser adotada é utilizar uma aproximação das equações de Saint-Venant que despreza apenas o termo de inércia advectiva na equação dinâmica. Esta aproximação vem sendo chamada de modelo inercial, ou método inercial. Esse artigo apresenta uma proposta de uma nova versão do modelo hidrológico MGB-IPH, em que o módulo de propagação de vazões na rede de drenagem, originalmente baseado no método Muskingum-Cunge, é substituído pelo método inercial. Além disso, o artigo também apresenta uma forma de representar redes de drenagem com ilhas ou laços, e extensas áreas inundáveis. A bacia do rio Araguaia foi escolhida para testar a metodologia, por possuir rios com alta e baixa declividade e regiões de planície de inundação, como a Ilha do Bananal. Os resultados do modelo MGB-IPH com propagação Inercial foram comparados com a propagação por Muskingum-Cunge através de medidas de desempenho entre as vazões calculadas e observadase mapas de inundação. Os resultados mostraram que utilizando o método inercial, o modelo MGB-IPH representa melhor os efeitos de atenuação da onda de cheia do que utilizando o métodoMuskingum-Cunge. Os mapas de inundação calculados pelo modelo MGB-IPH com o método de propagação inercial foram comparados com imagens do sensor MODIS mostrando que é possível representar, em escala regional, grandes regiões inundáveis como a Ilha do Bananal, mesmo sem dados detalhados de seções transversais e de topografia. Hydrological models are capable of representing basin scale processes through mathematic equations. Among them, flow routing is important, especially for large-scale basins. Flow routing methods used in hydrologic models are commonly simplified, and show satisfactory results in most cases. However, low slope rivers, rivers with backwater effects, or rivers with large flood plains cannot be well represented by simplified flow routing methods. An alternative for them is the full Saint-Venant equations model (hydrodynamic model). Nevertheless, this model demands more input data and its equations are more complex. An easier option than the full Saint-Venant equations model is the Inertial model, which ignores the inertial term from the Momentum equation. This paper presents an improvement of the MGB-IPH hydrologic model, where the simplified flow routing method is replaced by the Inertial model. Furthermore, this paper shows a methodology to represent island sand/or bifurcations in rivers. The Araguaia river basin was selected as a case study due to its particular physical characteristics. Results of the MGB-IPH with Inertial flow routing model were compared to conventional MGB-IPH Muskingum-Cunge method results, and to observed flow data. General outcomes show that observed flow data are better represented by the new improved MGB-IPH model. Furthermore, the simulated flood plain areas were visually well represented in comparison with MODIS imagery.

Previsões hidrológicas de curto e médio prazo normalmente estão associadas a muita incerteza, principalmente relacionadas com o uso de informações de previsões meteorológicas. Uma forma de considerar esta incerteza no processo de previsão e de tomada de decisão é o uso de previsões probabilísticas, obtidas na forma de previsões por conjunto (ou por ensemble). Nesta modalidade de previsões geralmente são gerados diversos cenários futuros possíveis, com o objetivo de quantificar as incertezas no processo de geração das previsões. Estudos recentes têm demonstrado benefícios promissores no uso de previsões por conjunto em detrimento de previsões determinísticas. O presente trabalho propõe uma abordagem complementar relacionada com o tema. Primeiramente é feita uma breve introdução sobre como previsões de vazão são utilizadas atualmente no Brasil, e como previsões por conjunto poderiam ser adotadas neste contexto. Em seguida é apresentada a análise dos resultados de um jogo realizado durante o XX Simpósio Brasileiro de Recursos Hídricos em 2013 na cidade de Bento Gonçalves (RS), que contou com a participação do público presente no evento no papel tomadores de decisão utilizando informações de previsão de vazão, usando e não utilizando informações de previsões probabilísticas. Os resultados do jogo apresentaram, de forma geral, que previsões com probabilidades leva- ram na maior parte dos casos a decisões melhores e mais coerentes entre múltiplas pessoas. Acredita-se que este argumento ajude a motivar desenvolvimentos de mais pesquisas e sistemas operacionais de previsão por conjunto no cenário brasileiro. Short and medium term hydrological forecasts are usually associated with large uncertainty, especially concerning the use of information from numerical wea- ther forecasts. One way to consider this uncertainty in the forecasting and decision-making process is the use of probabilistic predictions, obtained as ensemble forecasts. Usually at these kinds of forecasts several possible future scenarios are generated in order to quantify the uncertainties in the process. Recent studies have shown promising benefits on using ensemble forecasts rather than deterministic ones. This paper proposes a complementary approach related to the topic. First a brief introduction to streamflow forecasts current use in Brazil is made, and discussed how ensemble forecasts could be taken in this context. Then we present results of a game held during the XX Brazilian Symposium on Water Resources in 2013 in the city of Bento Gonçalves (RS), which had the public participation at the event playing the whole of decision makers that uses forecast information, using and not using information from probabilistic forecasts. The game results showed, in general, that probabilistic forecasts led in most cases to better and more consistent decisions among multiple persons. We believe that this argument helps to motivate the development of more research and more operational ensemble forecasting systems in the Brazilian scenario.

Na simulação numérica de escoamento, a simplificação da equação dinâmica pode ser útil, tanto em relação à estabilidade da solução numérica quanto ao ganho computacional. Tendo em vista esta preocupação, este trabalho avalia o uso de um modelo obtido ao se desconsiderar o termo advectivo da equação de Saint-Venant, denominado de modelo inercial local, ou simplesmente modelo inercial, na simulação de escoamentos influenciados por fechamentos de compor- tas, com enfoque na onda superficial provocada por estes fechamentos. Os resultados do modelo inercial são comparados com os resultados do modelo HEC -RAS, que utiliza uma solução da equação de Saint-Venant completa, a fim de se verificar no que a simplificação do termo convectivo acarreta. Os resultados obtidos com o modelo inercial se assemelharam aos do modelo HEC-RAS dependendo do tipo de onda formada, da vazão, e da altura do escoamento. Para escoamentos com a mesma vazão e altura, a validade do modelo inercial pode ser identificada utilizando adimensionais que relacionam características da onda com a altura do escoamento. Entretanto, para escoamentos expressivamente diferentes, estes adimensionais não indicaram a validade da simulação pelo modelo inercial, devendo-se analisar o peso dos termos da equação dinâmica. Observa-se que o código do programa com a aproximação inercial possui menos oscilações numéricas, quando se adota uma malha grosseira, do que o código do programa com as equações completas. Com o aumento da declividade do canal não foi identificado se importância do termo advectivo da equação dinâmica é maior ou menor. A partir deste trabalho, acredita-se que os testes realizados para ondas geradas pelo fechamento de comportas podem se estender para ondas superficiais formadas por outras fontes. In a numerical flow simulation, the simplification of partial equations can be useful for numerical solution stability and computational effort. This study evaluates the model obtained from the Saint-Venant equation simplification, which ignores the advective term, called inertial model in flow simulation influen- ced by closing gates, focusing on the surface wave caused by this closure. Inertial model results are compared with HEC-RAS model results, which use full Saint-Venant equation. For different problems, an inertial model can be appropriately used depending on wave type, flow rate and flow height. For similar flows (similar flow rate and height) the simulation validity can be verified by wave steepness, but for different flows, wave steepness should not be used and analysis of dynamic equation terms should be done. It is observed that inertial approximation code has less numerical oscillations that the HEC-RAS code. Channel slope does not indicate if the advective term is or not important. It is thought that waves generated by closing gate tests may be extend to surface waves generated by other sources.

A estimativa de recargas aquíferas é um fator determinante para a gestão sustentável das águas subterrâneas, uma vez que permite a manutenção da quantidade e qualidade da água dos aquíferos. Para a estimativa de recargas aquíferas existem diversos métodos que podem ser empregados e que variam de acordo com as características da área estudada e as escalas a serem trabalhadas. No presente estudo foram avaliadas as estimativas de recarga na Bacia Hidrográfica do Rio Ibicuí utilizando os métodos de Separação de Escoamento de Base e Variação dos Níveis d’Água. Levando em consideração as limitações de escala, os resultados da estimativa de recarga aquífera pelos dois métodos convergiram.

In the present study we investigate the use of medium-range streamflow forecasts in the Doce River basin (Brazil), at the reservoir of Aimorés Hydro Power Plant (HPP). During daily operations this reservoir acts as a “trap” to the sediments that originate from the upstream basin of the Doce River. This motivates a cleaning process called “pass through” to periodically remove the sediments from the reservoir. The “pass through” or “sediments flushing” process consists of a decrease of the reservoir’s water level to a certain flushing level when a determined reservoir inflow threshold is forecasted. Then, the water in the approaching inflow is used to flush the sediments from the reservoir through the spillway and to recover the original reservoir storage. To be triggered, the sediments flushing operation requires an inflow larger than 3000m3/s in a forecast horizon of 7 days. This lead-time of 7 days is far beyond the basin’s concentration time (around 2 days), meaning that the forecasts for the pass through procedure highly depends on Numerical Weather Predictions (NWP) models that generate Quantitative Precipitation Forecasts (QPF). This dependency creates an environment with a high amount of uncertainty to the operator. To support the decision making at Aimorés HPP we developed a fully operational hydrological forecasting system to the basin. The system is capable of generating ensemble streamflow forecasts scenarios when driven by QPF data from meteorological Ensemble Prediction Systems (EPS). This approach allows accounting for uncertainties in the NWP at a decision making level. This system is starting to be used operationally by CEMIG and is the one shown in the present study, including a hindcasting analysis to assess the performance of the system for the specific flushing problem. The QPF data used in the hindcasting study was derived from the TIGGE (THORPEX Interactive Grand Global Ensemble) database. Among all EPS available on TIGGE, three were selected: ECMWF, GEFS, and CPTEC. As a deterministic reference forecast, we adopt the high resolution ECMWF forecast for comparison. The experiment consisted on running retrospective forecasts for a full five-year period. To verify the proposed objectives of the study, we use different metrics to evaluate the forecast: ROC Curves, Exceedance Diagrams, Forecast Convergence Score (FCS). Metrics results enabled to understand the benefits of the hydrological ensemble prediction system as a decision making tool for the HPP operation. The ROC scores indicate that the use of the lower percentiles of the ensemble scenarios issues for a true alarm rate around 0,5 to 0,8 (depending on the model and on the percentile), for the lead time of seven days. While the false alarm rate is between 0 and 0,3. Those rates were better than the ones resulting from the deterministic reference forecast. Exceedance diagrams and forecast convergence scores indicate that the ensemble scenarios provide an early signal about the threshold crossing. Furthermore, the ensemble forecasts are more consistent between two subsequent forecasts in comparison to the deterministic forecast. The assessments results also give more credibility to CEMIG in the realization and communication of flushing operation with the stakeholders involved.

International audience; In this paper, a hydrologic/hydrodynamic modeling of the Amazon River basin iscpresented using the MGB-IPH model with a validation using remotely sensed observations. Moreover, the sources of model errors by means of the validation and sensitivity tests are investigated, and the physical functioning of the Amazon basin is also explored. The MGBIPH is a physically based model resolving all land hydrological processes and here using a full 1-D river hydrodynamic module with a simple floodplain storage model. Riverfloodplain geometry parameters were extracted from the SRTM digital elevation model, and the model was forced using satellite-derived rainfall from TRMM3B42. Model results agree with observed in situ daily river discharges and water levels and with three complementary satellite-based products: (1) water levels derived from ENVISAT altimetry data; (2) a global data set of monthly inundation extent; and (3) monthly terrestrial water storage (TWS) anomalies derived from the Gravity Recovery and Climate Experimental mission. However, the model is sensitive to precipitation forcing and river-floodplain parameters. Most of the errors occur in westerly regions, possibly due to the poor quality of TRMM 3B42 rainfall data set in these mountainous and/or poorly monitored areas. In addition, uncertainty in river-floodplain geometry causes errors in simulated water levels and inundation extent, suggesting the need for improvement of parameter estimation methods. Finally, analyses of Amazon hydrological processes demonstrate that surface waters govern most of the Amazon TWS changes (56%), followed by soil water (27%) and ground water (8%). Moreover, floodplains play a major role in stream flow routing, although backwater effects are also important to delay and attenuate flood waves.

In this work, a first-order upwind and a high-order flux-limiter schemes for solving the advection–diffusion equation on unstructured grids were evaluated. The numerical schemes were implemented as a module of an unstructured two-dimensional depth-averaged circulation model for shallow lakes (IPH-UnTRIM2D), and they were applied to the Guaíba River in Brazil. Their performances were evaluated by comparing mass conservation balance errors for two scenarios of a passive tracer released into the Guaíba River. The circulation model showed good agreement with observed data collected at four water level stations along the Guaíba River, where correlation coefficients achieved values up to 0.93. In addition, volume conservation errors were lower than 1% of the total volume of the Guaíba River. For all scenarios, the higher order flux-limiter scheme has been shown to be less diffusive than a first-order upwind scheme. Accumulated conservation mass balance errors calculated for the flux limiter reached 8%, whereas for a first-order upwind scheme, they were close to 18% over a 15-day period. Although both schemes have presented mass conservation errors, these errors are assumed negligible compared with kinetic processes such as erosion, sedimentation or decay rates.

Os procedimentos para decidir tecnicamente sobre pedidos de outorga de direito de uso da água envolvem etapas de comparação entre a demanda e a disponibilidade de água, as quais têm forte relação com os atributos espaciais da rede de drenagem e da localização dos diferentes e água em uma bacia hidrográfica. Idealmente, portanto, modelos hidrológicos utilizados para analisar a disponibilidade de água nos locais onde a outorga é solicitada deveriam ser integrados a um Sistema de Informação Geográfica (SIG). Neste trabalho apresentamos uma metodologia integração de modelos hidrológicos e SIG para auxílio à tomada de decisão em pedidos de outorga de direito de uso da água. A metodologia envolve: 1) a utilização de SIG para pré-processamento de dados visando a aplicação de um modelo hidrológico chuva-vazão; 2) a aplicação do chuva-vazão para gerar estimativas da vazão de referência adotada como disponibilidade hídrica em cada trecho de em que é dividida a bacia; 3) a transferência dos resultados do modelo hidrológico de volta ao SIG; e 4) a aplicação de um Sistema de Suporte à Decisão no qual cálculos simples Grande do Sul, de balanço hídrico são realizados internamente em um SIG. É apresentado um teste da metodologia na bacia do rio dos Sinos, no Rio Grande do Sul-Brasil, região em as demandas ultrapassam a disponibilidade de água. Os resultados no estudo de caso demonstraram que é possível aprimorar a metodologia de análise de outorga ao SIG e modelagem hidrológica para auxílio à tomada de decisão, facilitando assim o aproveitamento sustentável dos recursos hídricos numa bacia baseados no conhecimento da variabilidade sazonal e espacial da sua disponibilidade. The procedures for technical deciding about water use permits, involve steps of comparison between de‑ mand and availability of water which has a strong relation‑ ship between the spatial attributes of the drainage network and the location of different water users in a watershed. Therefore, ideally, the models used to make these analyses should be integrated into a Geographic Information System ‑ GIS. We propose a methodology for integrating GIS and hydrological models to support decision making in analysis of water use permits requests. The methodology involves 1) the use of GIS for pre‑processing of data for the application of a rainfall runoff hydrologic model; 2) the application of rainfall‑runoff model to generate estimates of the flow of reference adopted as water availability in every stretch of river; 3) transfer the results of the hydrological model back to the GIS; and 4) the application of a Decision Support System in which simple calculations of water balance are done entirely in GIS environment. The presented meth‑ odology was tested at the Dos Sinos River basin, at Rio Grande do Sul, Brazil, region where the demands exceed the water availability. The results of this case study showed that it is possible to improve the methodology used to analyze the granting of water rights for use permits when GIS and hydrological models are integrated.

É crescente a disponibilidade e utilização de campos de chuva estimados por sensoriamento remoto ou calculados por modelos de circulação da atmosfera, os quais são freqüentemente utilizados como entrada para modelos hidrológicos distribuídos. A distribuição espacial dos campos de chuva estimados é altamente relevante e deve ser avaliada frente aos campos de chuva observados. Este artigo propõe um método de comparação espaço-temporal entre campos de chuva observados e estimados baseado na comparação pixel a pixel e na construção de tabelas de contingência. Duas abordagens são utilizadas: (i) a análise integrada no espaço gera índices de performance que retratam a qualidade do campo de chuva estimada em reproduzir a ocorrência de chuva observada ao longo do tempo; (ii) a análise integrada no tempo produz mapas dos índices de performance que resumem a destreza das estimativas de ocorrência de chuva em cada pixel. Como exemplo de aplicação, é analisada a chuva estimada na climatologia do modelo global de circulação da atmosfera CPTEC/COLA sobre a bacia do Rio Grande. Utilizando-se cinco índices de performance, o método proposto permitiu identificar variações sazonais e padrões espaciais na performance das estimativas de chuva em relação a campos de chuva derivados de observações em pluviômetros. There is an increasing availability and application of precipitation fields estimated by remote sensing or calculated by atmospheric circulation models, which are frequently used as input for distributed hydrological models. The spatial distribution of the estimated precipitation fields is extremely important and must be verified against observed precipitation fields. This paper proposes a method for spatiotemporal comparison between observed and estimated precipitation fields based on a pixel by pixel comparison and on contingency tables. Two distinct approaches are carried out: (i) the spatial integrated analysis produces skill scores denoting the ability of the estimated precipitation field in reproducing the occurrence of observed precipitation along the time; (ii) the time integrated analysis generates maps of skill scores showing the reliability of the precipitation estimates in each pixel. As an example of application, the estimated precipitation climatology of the CPTEC/COLA global circulation model over the Rio Grande basin is assessed. Using five different skill scores, the proposed method identified seasonal variations and spatial patterns in the performance of the estimated precipitation fields in relation to precipitation fields derived from pluviometers measurement.

En los últimos años los cálculos de lluvia obtenidos por satélite son una fuente de datos potencialmente útil para estudios hidrológicos y pueden ser aplicados en modelos de transformación lluvia-escorrentía. La propuesta es de mayor interés para regiones remotas, donde los datos de lluvia obtenidos por métodos convencionales son limitados, debido a la baja densidad o escasez de estaciones pluviométricas. El objetivo de este trabajo es evaluar la calidad de los productos de lluvia por satélite como datos de entrada a un modelo hidrológico en la cuenca del río Huallaga en Perú. Se evaluaron dos productos de lluvia por satélite: 3B42-TMPA y CMORPH, además, lluvia obtenida por método convencional y una última fuente que resulta de las correcciones del producto CMORPH con lluvia convencional denotado como MIXTO. El periodo de simulación está comprendido entre octubre/2003 y febrero/2009. Los caudales simulados obtenidos por el modelo hidrológico para las cuatro fuentes de lluvia se asemejan a los regímenes hidrológicos, valores máximos para épocas húmedas y mínimos para épocas secas; sin embargo, con tendencias a subestimar los caudales observados. Las expectativas de considerar la lluvia obtenida por satélite como dato para complementar la medición de pluviómetros en la cuenca del río Huallaga son poco alentadoras. Los coeficientes de desempeño mostraron valores muy bajos de lo normal (coeficiente de Nash-Sutcliffe menor a 0,50) concluyendo que los productos de lluvia por satélite evaluados en este estudio, como datos de entrada al modelo hidrológico no resultan confiables. In recent years satellite rainfall estimates are a source of potentially useful data for hydrological studies and can be applied in rainfall-runoff transformation models. This study is of greater interest in remote regions where rainfall data obtained by conventional methods is limited due to low or poor density of surface rainfall stations. The aim of this study is to evaluate the quality of satellite rainfall products as input into hydrological model for the Huallaga River basin in Peru. Thus, two satellite rainfall products were evaluated: 3B42-TMPA y CMORPH and also rainfall obtained by conventional method; and a final source resulting from the corrections to the CMORPH product with conventional rainfall, referred to as MIXTO. The simulation period is from October/2003 to February/2009. The simulated discharge obtained from the hydrological model for the four sources of rain water regimes could be reproduced, maximum values for wet seasons and minimum values for dry seasons, however, it showed a trend to underestimate the observed discharge. Considering the satellite rainfall as a complement of raingauge measuring in the Huallaga River basin is not encouraging. Statistical skill showed very low values, under normal (Nash-Sutcliffe coefficient is less than 0.50) and concluded that the satellite rainfall products evaluated in this study as inputs into the hydrological model are not reliable.

A simulação matemática do escoamento em rios tem sido realizada principalmente com modelos unidimensionais (1D) acoplados ou não a células de armazenamento na planície. Para o caso de rios com grandes planícies de inundação, o extravasamento de água do canal para a planície e a propagação da inundação na planície governam a passagem da onda de cheia. Caso o interesse do estudo seja representar esses processos, métodos tradicionais não são adequados e uma abordagem recente tem sido proposta, baseada no acoplamento de um modelo 1D para simular o escoamento na calha principal e um modelo tipo raster para simular o escoamento bidimensional na planície. Este artigo apresenta um sistema computacional de simulação desenvolvido com esse propósito, com exemplo de aplicação à Bacia do Alto Paraguai. Os resultados obtidos indicam excelente reprodução do regime de vazões nas calhas dos rios e a simulação do padrão espaço-temporal de inundações na planície (Pantanal) bastante coerente. As áreas de inundação permanente são ditadas pela água extravasada dos canais principais, enquanto a ocorrência de precipitação é responsável pelo pulso de inundação sazonal, o que ressalta a importância da incorporação dos processos hidrológicos verticais na modelagem da planície. The mathematical simulation of river flow routing has been performed mostly using one-dimensional (1D) models coupled or not with floodplain storage cells. For rivers with large floodplains, the overflow of water from the main channel to the floodplains and the propagation of inundation are the main governing factors of flood wave routing. If the modeling study aims at representing these processes, traditional methods are not adequate. A more recent approach has been proposed, based on coupling a 1D model for channel flow routing and a raster-based model for two-dimensional floodplain inundation routing. This paper presents a computational simulation system developed for this purpose, and its application to the Upper Paraguay river basin. The results obtained show an excellent reproduction of the flow regime in the main channels and the coherent simulation of inundation patterns in space and time. Permanently flooded areas are governed by water exchanged from the main channels, while precipitation plays the key role for the seasonal flood pulse, highlighting the importance of incorporating vertical hydrological processes to floodplain modeling.

A evapotranspiração é uma das variáveis mais importantes do ciclo hidrológico, consistindo na ligação entre energia, clima e disponibilidade hídrica. É considerado um fenômeno bastante complexo, pois depende da interação entre diversas variáveis climáticas e do tipo e estágio de desenvolvimento da vegetação. O objetivo deste trabalho foi a análise de um conjunto amostral de variáveis climáticas, coletadas em dois sítios de monitoramento de vórtices turbulentos instalados em áreas de cultivos de cana-de-açúcar e remanescentes florestais de cerrado sensu stricto na bacia do rio Grande — SP, para identificar o comportamento de cada uma dessas variáveis climáticas e verificar quais apresentam uma maior influência no processo de evapotranspiração. Para tanto, foram aplicadas técnicas de estatística multivariada de análise de componentes principais e agrupamento de variáveis. Concluiu-se que o processo de evapotranspiração, para a área estudada, é comandado por três componentes, que explicam 75 e 82% da variância dos dados em áreas de cerrado e cultivos de cana-de-açúcar, respectivamente: (1) variação da oferta de energia, representada principalmente pelas variáveis Ki e Rn, (2) processo de circulação atmosférica local, representada pelas variáveis de velocidade do vento e pressão atmosférica, que regulam o sistema de precipitação e as condições de umidade do ar e do solo, e, (3) condições da vegetação, que apresentam modificações significativas com a alternância de estações secas e úmidas. One of the most important variables in the hydrological cycle is evapotranspiration, which contributes to the connection between energy, climate and water availability. It is considered a very complex phenomenon, for it depends on the interaction between several climate variables to identify the behavior of each climate variable and to assess factors which highly influence the evapotranspiration process. A set of samples of climate variables was analyzed. They were collected at two turbulent vortex monitoring sites installed in sugar cane fields and remaining “cerrado” sensu stricto in the Grande River basin, São Paulo, Brazil. Therefore, multivariate statistical techniques for analyzing main components and for grouping variables were applied. It was concluded that the evapotranspiration process consists of three components which explain 75 and 82% of variance data in “cerrado” and sugar cane fields, respectively: (1) variation of energy supply, represented mainly by Ki and Rn variables; (2) local atmospheric circulation process, represented by wind speed and atmospheric pressure variables, ruling the system of precipitation and air humidity and soil moisture, and; (3) conditions of the vegetation presenting significant changes in dry and humid seasons.

Os complexos processos do ciclo hidrológico podem ser representados por meio da modelagem hidrológica, sendo os modelos que simulam o processo de transformação da chuva em vazão os mais utilizados. Esses modelos são baseados em equações matemáticas que descrevem, de forma simplificada, o comportamento hidrológico da bacia e possuem parâmetros que devem ser definidos através de um processo de calibração. A calibração manual, por tentativa e erro pode ser uma tarefa tediosa, sobretudo quando o usuário do modelo é inexperiente. A calibração automática, por sua vez, utiliza técnicas numéricas de otimização baseadas no uso intensivo de computadores. Esse trabalho apresenta um algoritmo evolucionário multiobjetivo de otimização desenvolvido por Vrugt et al. (2003) e aplicado na calibração automática do modelo hidrológico IPH II. Os resultados obtidos são promissores, o algoritmo conseguiu uma aproximação uniforme do frente de Pareto nos diferentes testes realizados, mantendo os extremos da mesma bem representados. Ainda mostrou algumas vantagens sobre outro algoritmo evolucionário multi-objetivo atualmente utilizado na calibração automática do modelo hidrológico IPH II. The complex processes of the hydrological cycle can be represented through hydrological modeling, being the models that simulate the rainfall-runoff process the most used of them. These models are based in mathematical equations that describe, in a simplified way, the hydrological behavior of the basin and possess parameters that must be defined through a process of calibration. The manual calibration, by trial and error, can be a tedious task, especially when the model's user is inexperienced. The automatic calibration, however, utilizes numerical optimization techniques based in the intensive use of computers. This study presents a multi-objective evolutionary algorithm of optimization developed by Vrugt et al. (2003) and applied in the automatic calibration of the IPH II hydrological model. The obtained results are encouraging: the algorithm produced a uniform approach of the Pareto Front in all the different tests carried out, keeping well represented its extremities. Additionally this method displayed some advantages over another multi-objective evolutionary algorithm currently used for the automatic calibration of the IPH II hydrological model.

Este artigo apresenta uma discussão quanto à viabilidade de implementação do hidrograma ecológico ao avaliar: seus métodos de determinação; práticas que concorrem para sua implementação, e; a realidade da política brasileira de gerenciamento de recursos hídricos e meio ambiente, no tocante à conservação de rios. Como resultado, uma estrutura mínima para sua implementação consistiria em: legitimar o ecossistema como usuário de água, com prioridade inferior aos humanos básicos apenas; ajustar o gerenciamento de recursos hídricos para considerar hidrogramas ecológicos; definir hidrogramas regionais típicos, como parâmetro para definições locais; classificar rios com relação às suas características prévias à intervenção humana na bacia, incluindo características geomorfológicas, ecológicas e hidrológicas; e elaborar programas de pesquisa e monitoramento de índices que contribuam para refinar estimativas iniciais de hidrogramas ecológicos ao longo de sua operação. This paper presents some debate over the feasibility of applying environmental flows which mimic the natural regime through the assessment of: its prescription methods; water management techniques which would benefit its implementation, and; the Brazilian Water Resources and Environmental Management, especially fluvial systems conservation. As a result, the minimum framework to its implementation would consist of: legitimizing fluvial ecosystems as water resources users; setting management tools to consider environmental flows; establishing regional typical hydrographs, to serve as parameters to local definitions; classifying rivers based on its modification from undisturbed state, including geomorphological, ecological and hydrological aspects; developping research and monitoring programs to refine and evaluate the river status along its management.

A gestão de um ecossistema aquático sujeito a alterações antrópicas envolve analisar diferentes cenários de desenvolvimento econômico-social e ambiental no entorno e no próprio sistema hídrico. Ecossistemas aquáticos rasos tais como lagos e banhados estão suscetíveis à mudança de seu estado trófico quando submetidos a altas cargas de nutrientes e matéria orgânica. A produção primária responde a partir da disponibilidade de nutrientes, sugerindo que alterações na base da cadeia alimentar podem influenciar níveis tróficos acima. A eutrofização está altamente correlacionada com a variação na composição das populações de fitoplâncton, a qual pode ser caracterizada como uma forte indicadora biológica do nível trófico destes ecossistemas. Com o objetivo de avaliar o potencial de eutrofização de ecossistemas aquáticos rasos foi desenvolvido um modelo que representa os principais processos abióticos e bióticos destes sistemas. O modelo proposto é constituído de três módulos: (a) hidrodinâmico; (b) de transporte de nutrientes e (c) biológico relacionado à produtividade do fitoplâncton. Neste artigo foi realizada uma análise de sensibilidade dos parâmetros do modelo biológico, calibração dos parâmetros do módulo hidrodinâmico e simulações de eutrofização no Sistema Hidrológico do Taim (SHT) visando identificar regiões com alto potencial de eutrofização, utilizando a estimativa de clorofila a como indicador. A análise de sensibilidade dos parâmetros biológicos mostrou que a concentração de clorofila a é mais sensível às alterações dos parâmetros relacionados aos efeitos de temperatura da água e às perdas por respiração e consumo por zooplâncton e outros organismos. As simulações de eutrofização mostraram que as regiões mais rasas (litorâneas) no SHT apresentam condições mais favoráveis à eutrofização, enquanto que no canal principal (região pelágica) e na região de Banhado, onde existe redução da radiação solar devido a barreira de macrófitas aquáticas presente, a produtividade é baixa. Water resource management involves analyzing the different economic-social scenarios and environmental development. A shallow aquatic ecosystem such as wetlands or lakes can become eutrophic due affluent conditions, as well as physical, chemical and biological elements. Primary production responds to nutrient availability, suggesting that alterations at the top of the food-web can influence the highest trophic levels. Eutrophication is correlated with the variation in the composition of the phytoplankton population, which can be characterized as a strong biological indicator of the trophic level of those ecosystems. Aiming to evaluate the potential of eutrophication in shallow aquatic ecosystems a model was developed to represent the main biotic and biotic processes of these systems. The considered model is comprised of three modules: (a) hydrodynamic; (b) nutrient transport and (c) biological which deals with phytoplankton productivity. This paper presents the sensitivity analysis of the biological parameters model, parameter calibration of the hydrodynamic model and eutrophication simulation of the Taim Hydrologic System (SHT-Sistema Hidrólogico do Taim), whose aim is to identify regions with a high eutrophication potential. The sensitivity analysis of the biological parameters has shown that the chlorophyll a concentration is most sensitive to changes in the parameters involving water temperature effect and the losses due to respiration and grazing. The eutrophication simulations have shown that the shallow regions present the most favorable conditions for eutrophication, while in the deep region (pelagic region) and in the wetland region where solar radiation is decreased due to the aquatic macrophyte barrier found there, productivity is lower.

O comportamento hidrológico de grandes bacias envolve a integração da variabilidade espacial e temporal de um grande número de processos. No passado, o desenvolvimento de modelos matemáticos precipitação – vazão, para representar este comportamento de forma simplificada, permitiu dar resposta às questões básicas de engenharia. No entanto, estes modelos não permitiram avaliar os efeitos de modificações de uso do solo e a variabilidade da resposta em grandes bacias. Este trabalho apresenta o desenvolvimento e a validação de um modelo hidrológico distribuído utilizado para representar os processos de transformação de chuva em vazão em grandes bacias hidrográficas (maiores do que 10.000 km2). Uma grade regular de células de algumas dezenas ou centenas de km2 é utilizada pelo modelo para representar os processos de balanço de água no solo; evapotranspiração; escoamentos: superficial, sub-superficial e subterrâneo na célula; e o escoamento na rede de drenagem em toda a bacia hidrográfica. A variabilidade espacial é representada pela distribuição das características da bacia em células regulares ao longo de toda a bacia, e pela heterogeneidade das características no interior de cada célula. O modelo foi aplicado na bacia do rio Taquari Antas, no Rio Grande do Sul, na bacia do rio Taquari, no Mato Grosso do Sul, e na bacia do rio Uruguai, entre Rio Grande do Sul e Santa Catarina. O tamanho destas bacias variou entre, aproximadamente, 30.000 km2 e 75.000 km2. Os parâmetros do modelo foram calibrados de forma manual e automática, utilizando uma metodologia de calibração automática multi-objetivo baseada em um algoritmo genético. O modelo foi validado pela aplicação em períodos de verificação diferentes do período de calibração, em postos fluviométricos não considerados na calibração e pela aplicação em bacias próximas entre si, com características físicas semelhantes. Os resultados são bons, considerando a capacidade do modelo de reproduzir os hidrogramas observados, porém indicam que novas fontes de dados, como os fluxos de evapotranspiração para diferentes coberturas vegetais, serão necessários para a plena utilização do modelo na análise de mudanças de uso do solo.

As barragense os reservatórios são projetados considerando um risco de falha. As consequências de uma falha devem ser avaliadas, como parte do projeto de uma barragem e também para barragens já existentes, com o objetivo de diminuir os riscos a que está exposta a população. Este estudo buscou revisar aspectos relevantes do impacto do rompimento de uma barragem, avaliar as incertezas envolvidas quando se utilizam modelos matemáticos como ferramenta básica, e desenvolver a análise para o reservatório de Ernestina, no rio Jacuí, no Rio Grande do Sul. O modelo simplificado SMPDBK e o modelo hidrodinâmico DAMBRK, do National Weather Service (Serviço meteorológico), dos E.U.A., utilizados como metodologia básica, são descritos em todos os aspectos relevantes à simulação do rompimento de barragens. Os modelos SMPDBK e, mais profundamente, DAMIBRK foram comparados, testados e aplicados ao caso do rompimento hipotético da barragem de Ernestina. Foi realizada uma análise de sensibilidade dos resultados do modelo DAMBRK a diversos parâmetros e hipóteses de representação matemática. O rompimento hipotético da barragem de Ernestina foi analisado, desde a probabilidade do seu acontecimento, até a estimativa preliminar da população atingida e dos seus prquízos. O modelo SMPDBK apresentou resultados razoáveis frente ao modelo DAMBRK. Os resultados do modelo DAMBRK mostraram-se especialmente sensíveis ao tamanho máximo que a brecha atinge no seu desenvolvimentoe ao nível da água no interior do reservatório quando se inicia o rompimento. O tempo de formação da brecha e a rugosidade do leito do rio também se revelaram importantes. Os resultados preliminares mostraram que o rompimento hipotético de Ernestina deve atingir diretamente as residências de cerca de 800 habitantes. Aproximadamente 600 hectares de terra cultivada devem ser atingidos e pelo menos quatro pontes sobre o rio Jaçuí devem ser destruídas. As barragens de Passo Real e Maia Filho, bem como todo o sistema de geração de energia do rio Jacuí, não devem sofrer conseqüências. Dams and reservoirs are designed assuming a level of risk. Consequentes of the failure of a new dam must be taken into account as part of its design and, for existing dams, as a means to lower the risk it poses to the population. This work reviewed relevant aspects of dam break analysis, availed the uncertainty when using a mathematical modeling technique and developed a dam break analysis for the Ernestina dam, on Jacuí River, Rio Grande do Sul. The simplified mathematical modem SMPDBK and the hydrodynamic model DAMBRK, both from National Weather Service, U.S.A ., were used as the basic methodology, and are described in all aspects relevant to dam break analysis. Both models, especially DAMBRK, were compared, tested ad applied to the hypothetical ruptura of Ernestina dam. The sensitivity of DAMBRK to many parameters and mathematical representation analysis was analyzed. The hypothetical failure of Ernestina dam was analyzed, from the failure risk to a preliminary estimate of the endangered population. Results of the SMPDBK model were reasonable as compared to DAMBRK's. DAMBRK results showed great sensitivity to maximum breach width and initial reservoir water level. A little lower, but also important, was the sensitivity of DAMBRK to breach formation time interval and river boundary friction. According to the preliminary results, the hypothetical Ernestina dam failure will cause direct damages to the houses of 800 people. Near 600 hectares of agriculture land will be submerged and almost tour bridges will be destroyed. The downstream located dams, Passo Real and Maia Filho, and the whole Jacuí Generation System, will not be damaged.

A América do Sul (AS) possui uma área de ~17,8 milhões de km² e contribui com cerca de 30 % do volume total de água escoado para os oceanos. Gerenciar os recursos hídricos continentais na AS é uma tarefa desafiadora devido à heterogeneidade das características climáticas, suscetibilidade a cheias e secas, complexidade dos sistemas de rios e existência de usos múltiplos da água, em particular a produção de energia onde a hidroeletricidade corresponde por ~50 % da capacidade instalada na AS. Previsões de vazão, por exemplo, são essenciais para a operação de reservatórios e para melhorar a preparação contra eventos extremos, os quais podem ocorrer em múltiplas escalas espaciais e afetar extensas áreas inclusive para além de fronteiras políticas. Neste sentido, há necessidade de ferramentas que possibilitem quantificar fluxos hidrológicos com consistência espaço-temporal na AS, bem como prever a condição futura dos rios com estimativa de incerteza para dar suporte à tomada de decisão. O objetivo desta tese foi estabelecer uma base técnica para modelagem hidrológica com abrangência da América do Sul e dar contribuições para previsão de vazões em escala continental, no horizonte de médio prazo (até 15 dias). O modelo hidrológico MGB foi adaptado para simulação hidrológica–hidrodinâmica em múltiplas bacias sul americanas a partir de técnicas de geoprocessamento, produtos de sensoriamento remoto e outras bases de dados globais. O modelo desenvolvido, denominado MGB–SA, teve representação satisfatória de vazões em intervalo de tempo diário, apresentando acurácia superior à de modelos globais estado-da-arte. Também foram representados de forma adequada padrões de sazonalidade e magnitude de evapotranspiração e armazenamento terrestre de água nas principais bacias da AS, com desempenho razoável para os níveis simulados. Através de um experimento de previsão que assume o modelo hidrológico como livre de erros, foi identificado que previsões de chuva por conjunto (ensemble) com horizonte de médio prazo podem contribuir para a destreza das previsões de vazão com até 15 dias à frente. Na avaliação de desempenho geral das previsões de vazão, observou-se maior destreza para regiões situadas ao leste e sudeste da América do Sul, incluindo bacias como Uruguai, Alto Paraná, São Francisco, Parnaíba e Tocantins-Araguaia. Em relação à discriminação entre ocorrência e não ocorrência de eventos de vazão elevada, a destreza da previsão foi maior nos rios do sul do Brasil. Além disso, as previsões exibiram maior destreza para acurácia geral e discriminação em regiões de clima temperado e úmido (sub-tropical). Quando considerados ambos os erros do modelo hidrológico e da precipitação prevista, as previsões de vazão com horizonte de médio prazo demonstraram baixa ou nenhuma destreza na maioria das localidades analisadas além de um grande viés probabilístico, com a maior parte das observações situando-se fora dos limites do conjunto mesmo para antecedências mais longas após 1 semana à frente. Com aplicação de pós-processamento estatístico a destreza das previsões de vazão aumentou substancialmente, melhorando a consistência estatística e proporcionando largura de incerteza geralmente inferior ao da climatologia das vazões observadas. O desempenho de séries temporais de previsão por conjunto obtidas de distribuições preditivas pós-processadas também foi avaliado, demonstrando limitações principalmente em condições de baixo espalhamento do conjunto bruto. Por fim, os resultados deste trabalho sugerem que o modelo desenvolvido pode servir como uma ferramenta para dar suporte a pesquisas que contemplem os recursos hídricos da AS de forma integrada, e que previsões de vazão com destreza no horizonte de médio prazo podem ser obtidas na escala continental mesmo com dados in situ limitados. South America (SA) has an area of ~17.8 million km² and contributes around 30 % of the global runoff to the oceans. Managing the continental waters in SA is a challenge given the wide range of climatic characteristics, susceptibility to floods and droughts, complexity of river systems and existence of multiple water uses, in particular the energy production where hydropower accounts for ~50 % of the total installed capacity. Streamflow forecasts, for instance, are essential for reservoir operation and to improve preparedness against extreme hydrometeorological events, which can occur on multiple spatial scales and affect extensive regions including areas beyond the political boundaries. In this context, there is a need for tools that enable to quantify hydrological fluxes in SA with spatiotemporal consistency, as well as to predict future river conditions with uncertainty estimates to support decision making. The objective of this thesis was to establish a technical basis for hydrological modeling in the South America domain and contribute to medium-range (up to 15 days) streamflow forecasting at the continental scale. The MGB hydrological model was adapted for simulation of multiple South American basins by using hydrologic–hydrodynamic methods, geoprocessing techniques, remote sensing products and other global scale datasets. The developed model, named MGB–SA, showed a satisfactory representation of daily discharges and exhibited higher accuracy compared with state-of-the-art global models. Seasonality and magnitude of both evapotranspiration and terrestrial water storage were adequately represented at the scale of main SA basins, with a reasonable performance regarding the simulated water levels. Through a forecasting experiment that assumes the hydrological model as error-free, it was found that medium-range ensemble precipitation forecasts can contribute to the skill of streamflow forecasts up to 15 days ahead. Overall skill was higher in the eastern and southeastern SA regions, which include large basins such as Uruguay, Upper Paraná, São Francisco, Parnaíba and Tocantins-Araguaia. Skill in discrimination between the occurrence and non-occurrence of high flow events was mainly observed in rivers of southern Brazil. In addition, forecasts showed better skill in terms of overall quality and discrimination in warm temperate, fully humid regions. When the errors of both hydrological model and predicted precipitation were taken into account, medium-range streamflow forecasts exhibited low or no skill in most of the analyzed locations in addition to a large probabilistic bias, with most of the observations falling outside the ensemble limits even at longer lead times, after 1 week ahead. Statistical postprocessing substantially improved streamflow forecast skill, enhancing statistical consistency and generally providing lower uncertainty width compared with that of observed climatology. The performance of streamflow ensemble traces derived from postprocessed predictive distributions was also evaluated, showing limitations mainly in conditions for which the raw ensemble spread is low. Finally, the results of this work suggest that the developed model can be used to simulate the South American water resources in an integrated manner, and that skillful, medium-range streamflow forecasts could be obtained at the continental scale of SA even with limited in situ data.

A região semiárida do Brasil apresenta um regime de precipitação com alta variabilidade espacial e temporal, além de elevada evaporação potencial, o que limita a disponibilidade de água para as atividades econômicas. Uma solução historicamente adotada para este problema, tanto pelo Estado como por empreendedores individuais, foi a construção de reservatórios, estruturas capazes de transpor a disponibilidade hídrica no tempo, de um período chuvoso para um seco. No Nordeste brasileiro, o processo de construção de açudes se iniciou ainda no século XIX e hoje há mais de 105.000 barramentos apenas no estado do Ceará, com predominância dos de pequeno porte. Essa densa rede de reservatórios é de difícil gestão, pois pouco se sabe sobre as características de cada estrutura e sobre as interações e efeitos sobre a hidrologia local. Modelos hidrológicos surgem como uma solução para este problema, porém, em sua maioria, não são idealizados para incorporar estruturas antrópicas em seu escopo, a exemplo do Modelo de Grandes Bacias (MGB-IPH). Assim, este trabalho tem o objetivo de desenvolver um módulo de reservatórios para incorporação ao MGB-IPH visando à modelagem de bacias hidrográficas de grande escala considerando densas redes de pequenos reservatórios sobre as quais existe pouco ou nenhum dado. Para isso, foram desenvolvidas metodologias de análise e geração de dados sobre pequenos reservatórios distribuídos sobre a sub-bacia do Banabuiú, bacia adotada para estudo de caso. Primeiramente, foram investigadas as relações e a influência de diferentes formas de representação dos açudes dentro do modelo, variando-se a quantidade de açudes representados explicitamente (on-grid) e a quantidade utilizada para gerar os açudes equivalentes (aqueles em sub-grid) através da utilização de diferentes densidades de drenagem. Em relação ao reservatório equivalente, foi feita uma comparação entre a metodologia proposta da tabela incremental que considera as respectivas áreas de contribuição de cada açude, bem como as interações entre eles; com uma versão alternativa que considerava os reservatórios equivalentes no exutório de cada minibacia. Por fim, foram analisados diferentes cenários de interação entre pequenos e grandes reservatórios através de curvas de permanência de atendimento a demandas. A partir dos resultados obtidos, pode-se inferir que a utilização de maiores densidades de drenagem foi a abordagem mais adequada, apresentando os melhores resultados globais do coeficiente de Nash-Sutcliffe (NS) para a bacia. Por sua vez, na comparação entre os reservatórios equivalentes, a metodologia proposta apresentou-se a mais adequada, sendo capaz de melhorar globalmente os coeficientes do modelo. Quanto às curvas de permanência, observou-se que a presença dos pequenos reservatórios reduz a entrada de água nos reservatórios estratégicos, diminuindo a capacidade destes açudes de atender às demandas da bacia. Dessa forma, embora existam limitações quanto aos resultados apresentados, observa-se que houve um considerável avanço na modelagem desses tipos de bacias dentro do modelo MGB-IPH. The semi-arid region of Brazil presents a precipitation regime with high spatial and temporal variability, in addition to high potential evaporation ratings, which limits the availability of water for economic activities. A historical solution adopted to solve this problem, both by the State and by individual entrepreneurs, was the construction of reservoirs. These structures are capable of transferring water availability over time, from a wet period to a dry one. In the Brazilian Northeast, the process of construction of these dams date back to the 19th century and nowadays there are more than 105.000 dams in the state of Ceará alone, being the small ones predominant. This dense network of reservoirs is difficult to manage, as little is known about the characteristics of each structure and the interactions and effects on local hydrology. Hydrological models emerge as a solution to this problem. However, for the most part, they are not designed to incorporate anthropic structures in their scope, such as the Modelo de Grandes Bacias (MGB-IPH). Hence, the goal of this paper is to develop a reservoir module for incorporation into the MGB-IPH aiming at modeling large-scale hydrographic basins considering dense networks of small reservoirs on which little or no data exists. On this context, a methodology for analysis and data generation was developed and applied on the group of small reservoirs distributed over the Banabuiú sub-basin, adopted as a case study. Firstly, the relationships and influence of different forms of representation of the dams within the model were investigated, varying the amount of dams represented explicitly (on-grid) and the amount used to generate the equivalent dams (those in sub-grid). These analyses were made by using different drainage densities. In relation to the equivalent reservoir, an incremental table methodology was proposed, considering the respective areas of contribution of each dam, as well as the interactions between them. Then, a comparison was made between this methodology and an alternative version that considered equivalent reservoirs at the outlet of each mini-basin. Finally, different scenarios of interaction between small and large reservoirs were analyzed through demand permanence curves. From the results obtained, it can be inferred that the use of higher drainage densities was the most appropriate approach, presenting the best overall results of the Nash-Sutcliffe coefficient. In the comparison between equivalent reservoirs, the proposed methodology was the most adequate, being able to globally improve the coefficients of the model. As for the permanence curves, it was observed that the presence of small reservoirs reduces the entry of water into strategic reservoirs, reducing their capacity to supply the basin needs. Thus, although there are limitations regarding the results presented, it is observed that there has been a considerable advance in the modeling of these types of basins within the MGB-IPH model.

As reservas de água subterrânea estão entre os componentes do balanço hídrico com maior incerteza de quantificação. O entendimento da sua disponibilidade ainda é limitado em comparação com outras reservas hídricas como rios, atmosfera, solo e lagos. A busca por esse recurso tem aumentado tanto nas áreas úmidas quanto nas áreas secas. Explorar os processos de interação entre as águas subterrâneas e os rios em aquíferos de grande escala (como o comportamento das reservas em períodos secos, eventos de recarga episódica e variabilidade anual), é de grande interesse para a indústria, a economia e a qualidade de vida da população. O desenvolvimento e a avaliação de ferramentas é um passo inicial para o uso sustentável das águas subterrâneas. O objetivo desta tese é avançar no entendimento de processos hidrogeológicos de grande escala que ocorrem no Brasil em diferentes climas e formações aquíferas a partir de ferramentas inovadoras e complementares de monitoramento in situ, dados de sensoriamento remoto e modelagem hidrológica. Verificou-se que o uso do GRACE para detectar variações nas reservas de água subterrânea para um aquífero sedimentar no semiárido brasileiro apresentou resultados adequados. As reservas de água subterrânea estão entre os componentes do balanço hídrico com maior incerteza de quantificação. O entendimento da sua disponibilidade ainda é limitado em comparação com outras reservas hídricas como rios, atmosfera, solo e lagos. A busca por esse recurso tem aumentado tanto nas áreas úmidas quanto nas áreas secas. Explorar os processos de interação entre as águas subterrâneas e os rios em aquíferos de grande escala (como o comportamento das reservas em períodos secos, eventos de recarga episódica e variabilidade anual), é de grande interesse para a indústria, a economia e a qualidade de vida da população. O desenvolvimento e a avaliação de ferramentas é um passo inicial para o uso sustentável das águas subterrâneas. O objetivo desta tese é avançar no entendimento de processos hidrogeológicos de grande escala que ocorrem no Brasil em diferentes climas e formações aquíferas a partir de ferramentas inovadoras e complementares de monitoramento in situ, dados de sensoriamento remoto e modelagem hidrológica. Verificou-se que o uso do GRACE para detectar variações nas reservas de água subterrânea para um aquífero sedimentar no semiárido brasileiro apresentou resultados adequados. Groundwater volumes are among the water balance components with the greatest uncertainty of quantification, the understanding of its availability is still limited compared to other water reserves such as rivers, atmosphere, soil, and lakes. The search for this continuous supply resource throughout the year has increased in wet and dry areas. Exploring hydrological, hydrogeological, and surface-groundwater interaction processes among these large-scale South American aquifers, such as the dynamics associated with dry periods response, recharge events, and interannual variability, is of great interest to the industry, the economy, and the quality of life of the regional population. And the development and testing of tools for researching these aspects is a primer step for sustainable usage of groundwaters in the South American domain. The main objective of this thesis is to advance in the determination of large-scale hydrogeological process in Brazil, South America in different climates and aquifer formations from innovative and complementary tooling of intensive field monitoring, remote sensing data, and hydrological modeling. We found that the use of GRACE to detect variations in groundwater reserves showed adequate results for a small-scale sedimentary aquifer in the Brazilian semi-arid region. The findings were promising to improve the understanding of droughts at different scales in those areas. GRACE data also showed itself as an essential tool for monitoring groundwater volumes in the other two aquifers in humid subtropical areas and investigated in this thesis (Caiuá Aquifer and SASG). We also found that in those humid subtropical areas, the high soil moisture storage has an important role in the occurrence of large episodic recharge events. Atypical rainfall in winter periods was responsible for the increase in soil moisture that explains the larger events. The changes in aquifer storage caused by episodic recharge events are long-lasting and directly affect low flows in rivers with implications on hydro-climatic variability. We also brought important findings related to groundwater variations in fractured aquifer systems, which are complex to predict. Significant contributions related to groundwater level variations, transit times, regional trends, and interaction with the rivers in the SASG were presented.

A América do Sul abriga alguns dos maiores sistemas hídricos do planeta, frequentemente associados a grandes planícies de inundação, como o Pantanal e várias áreas da Amazônia. Áreas úmidas (AU’s) interfluviais são também encontrados no continente, com características geomorfológicas particulares, e vegetações de savana e gramíneas únicas. As AU’s da América do Sul provêm diversos serviços ecossistêmicos, como suporte à biodiversidade, provisão de alimento e atenuação de cheias. Humanos têm se estabelecido ao redor de AU’s por milênios, se beneficiando dos recursos providos por elas. Eles se adaptaram ao seu regime de inundação, e adaptaram sua paisagem, definindo o que tem sido chamado de sistemas sociedade-água. Por outro lado, um número crescente de pessoas têm sido negativamente afetado por cheias extremas. Da escala continental à local, esta tese convida o leitor a uma jornada através de importantes AU’s da América do Sul e suas particulares dinâmicas de inundação, sob a luz da era dos satélites e dos grandes avanços em modelagem hidrológica-hidrodinâmica das últimas décadas. Este trabalho é baseado na proposta de uma escala continental de pesquisa sobre AU’s, e é baseado em uma abordagem de hidrologia comparativa. Inundações são estudadas em múltiplas dimensões, de processos de AU’s naturais à questão do perigo para humanos. A primeira parte apresenta uma série de estudos sobre as AU’s da bacia amazônica, desde o desenvolvimento de modelos 1D e 2D para simular processos hidrológicos em tipos contrastantes de AU’s na bacia do Rio Negro, até a intercomparação de 29 produtos de inundação e avaliação de tendências de inundações de longo prazo para a escala da bacia amazônica. Enquanto a maioria dos estudos de AU’s foi conduzida nas várzeas do rio Amazonas, importantes lacunas do conhecimento permanecem para a compreensão da dinâmica hidrológica de áreas interfluviais como Llanos de Moxos e as savanas do rio Negro, onde a inundação é menos previsível e mais rasa. A segunda parte da tese utiliza dados oriundos de satélites relacionados a múltiplas variáveis hidrológicas (níveis d’água, armazenamento total de água, extensão de áreas inundadas, precipitação e evapotranspiração) para estudar a hidrologia de 12 grandes sistemas de AU’s do continente. São destacadas as grandes diferenças entre planícies de inundação e AU’s interfluviais em termos de amplitude anual de níveis d’água, defasagem entre precipitação e inundação, e dinâmica de evapotranspiração. Por fim, a última parte da tese aborda o componente de perigo de inundação das interações sociedade-água através de avaliações em grande escala da dinâmica de inundação e dos efeitos de infraestruturas construídas (como barragens) na atenuação de cheias. A dinâmica das grandes cheias de 1983, um dos anos mais extremos já registrados no continente, é avaliada com um modelo hidrológico continental. Depois, a capacidade de modelos continentais para simular o continuum entre rios, planícies de inundação e reservatórios que existe em grandes bacias hidrográficas é avaliada com estudos de casos para importantes bacias afetadas pela intervenção humana (bacia dos rios Paraná e Itajaí-Açu). Enquanto esta tese avança a compreensão de relevantes processos hidrológicos relacionados a inundações na América do Sul em múltiplas escalas, bem como seus efeitos positivos e negativos nas sociedades humanas e ecossistemas em geral, importantes lacunas do conhecimento persistem e fomentam importantes oportunidades de pesquisa futuras. O lançamento de várias missões satelitais orientadas a hidrologia, e uma cada vez mais crescente capacidade computacional, faz da agenda continental de hidrologia relacionada a AU’s e inundações um grande tópico de pesquisa para os próximos anos. South America hosts some of the major river systems on Earth, often associated with large floodplains that are inundated every year, such as the Pantanal and many Amazon wetlands. Interfluvial wetland complexes are also found across the continent, with particular geomorphic settings and unique savanna or grassland vegetation. South American wetlands can provide distinctive ecosystem services such as biodiversity supporting, food provision and flood attenuation. On the other hand, humans have settled around wetlands for millennia, benefiting from all resources they provide, and have adapted to its flood regime as well adapted its landscape, defining what has been called human-water systems. Yet, an increasing number of South American people have been negatively affected by extreme floods. Moving from continental to local scales, this thesis invites the readers to a journey across major South American wetland systems and their unique hydrological dynamics, under the light of the satellite era and the breakthrough advances on hydrologic-hydrodynamic modeling in the last decades. This work is founded on the proposition of a continental wetland research agenda, and based on a comparative hydrology approach. Floods are studied through both natural wetland processes and hazard dimensions. The first part presents a set of studies on the Amazon basin wetlands, from the development of 1D and 2D models to simulate hydrological processes in contrasting wetland types in the Negro river basin to the basin-wide intercomparison of 29 inundation products and assessment of long-term inundation trends. While most wetland studies have been conducted over the central Amazon floodplains, major knowledge gaps remain for understanding the hydrological dynamics of interfluvial areas such as the Llanos de Moxos and Negro savannas, where the inundation is less predictable and shallower. The second part of the thesis leverages satellite-based datasets of multiple hydrological variables (water levels, total water storage, inundation extent, precipitation and evapotranspiration) to address the hydrology of 12 large wetland systems in the continent. It shows the major differences among river floodplains and interfluvial wetlands on the water level annual amplitude, time lag between precipitation and inundation, and evapotranspiration dynamics. Finally, the third part addresses the flood hazard component of human-wetland interactions through large-scale assessments of flood hazard dynamics and effects of built infrastructure (dams) on flood attenuation. The dynamics of the great 1983 floods, one of the most extreme years ever recorded in the continent, is assessed with a continental hydrological model. Then, the capabilities of continental models to simulate the river-floodplain-reservoir continuum that exists across large river basins are assessed with case studies for major river basins affected by human intervention (Itajaí-Açu and upper Paraná river basins in Brazil). While this thesis enlightens some relevant hydrological processes regarding South American floods and their positive and negative effects to human societies and ecosystems in general, major knowledge gaps persist and provide great research opportunities for the near future. The launching of many hydrology-oriented satellite missions, and an ever-growing computational capacity, make the continental hydrology agenda related to wetlands and floods a great research topic for the upcoming years.

Riparian zones (RZs) have a clear distinct behaviour than the rest of the landscape. Particularly in water-limited regions, such as the Brazilian Savannah (Cerrado biome), where dry season may extend 5 months, the difference between riparian and upland zones is highly pronounced due to vegetation water access to groundwater, and this can have implications on the climatic and hydrological cycles. In order to quantify this difference at large-scale, it was herein proposed to (1) map RZs using topographical information, (2) investigate how land cover is distributed among topographic gradients and (3) investigate vegetation behaviour through remote sensing vegetation measurements and evapotranspiration (ET) estimation. A 140,000 km² upland region inside the Cerrado biome, called the Urucuia aquifer system, was chosen as study site. The region has seen a huge agricultural expansion during the last decades, with mechanized and irrigated crops increasingly using water from its underground reserves, which associated with climate change can have a big impact on the ecosystem, and understanding the role of RZs can be essential to quantify this impact. The height above nearest drainage (HAND) index was used to map RZs, by visually assessing bellow which values the index provided a reasonable RZ buffer comparing with Google Earth imagery. We also used HAND to quantify across its values the historical land cover distribution obtained by the MapBiomas database, and analyse vegetation behaviour in RZs and upland zones (UZs) using remote sensing vegetation measurements of normalized difference vegetation index (NDVI) and normalized difference moisture index (NDMI) and ET estimation from the surface energy balance algorithm for land (SEBAL). A necessary step for HAND computation is a defined stream network, for which the main challenge is identifying channel heads. Herein it was developed an algorithm that produced a varying draining area threshold (vDAT) map for channel initiation, using the topographic position index (TPI) as an auxiliary variable. This algorithm was tested using MERIT-DEM. With the stream network, HAND values bellow 5 m provided the best RZ buffer. As for land cover distribution, we captured that forests naturally occur more densely in the extreme values of HAND (very shallow and very deep) and that farmland historical occupation in the Urucuia region occur more in the upper portions of the terrain, possibly due to soil conditions stablished during landscape formation and evolution. As for vegetation activity, the land cover class seems to have more influence on vegetation behaviour than topographic position, for all indicators computed. Yet, NDMI values in Riparian Forests are greater than in Upland Forests, particularly towards drier conditions, in terms of both seasonality (drier months) and inter-annual variability (drier years). Despite this indication of more water available in RZs than UZs, the ET estimation could not capture these differences, possibly due to difficulties in estimating this variable in natural vegetation with high degree of water stress.

Os processos de tomada de decisão que envolvem recursos hídricos no Brasil vêm se baseando unicamente em análises do comportamento histórico de variáveis hidrológicas e climatológicas. Entretanto, isto pode levar a estratégias equivocadas relacionadas ao uso dos recursos naturais, pois o impacto das mudanças climáticas no comportamento das vazões futuras pode alterar consideravelmente a disponibilidade de recursos hídricos. Na presente pesquisa é proposta uma análise da variação das disponibilidades hídricas de bacias hidrográficas contidas na área de contribuição à Laguna dos Patos, Rio Grande do Sul, Brasil, que pode ser consequência das mudanças climáticas previstas por modelos globais do CMIP5 (Coupled Model Intercomparison Project Phase 5), publicados no mais recente relatório do IPCC (International Panel on Climate Change). Dois períodos futuros de 30 anos foram simulados, de 2006 a 2035 e de 2051 a 2080, utilizando o modelo hidrológico MGB-IPH, considerando dois ctnários extremos de emissões de gases de efeito estufa e vinte modelos climáticos. Foi possível verificar uma tendência de aumento das vazões na bacia hidrográfica simulada, pois mais de 60% das simulações realizadas indicou algum percentual de aumento de vazões médias em todas as unidades de discretização modeladas. A análise dos resultados das simulações indicou que os dados dos modelos climáticos HadGEM2-ES e GFDL-CM3 utilizados como dados de entrada no modelo hidrológico são os que fornecem, respectivamente, limiares superiores e inferiores de vazão do conjunto de resultado simulados. A avaliação conjunta dos resultados gerados por estes dois modelos, associados ao cenário de altas emissões de gases de efeito estufa, é capaz de gerar cenários extremos de vazões que são projetados para o futuro considerando as mudanças climáticas na bacia hidrográfica da Laguna dos Patos. Já os resultados fornecidos pelos modelos bcc-csm1-1, BNU-ESM e CNRM-CM5 são similares à mediana do conjunto das simulações geradas por todos os modelos avaliados nesta pesquisa. Além disso, a região mais ao norte da área de estudo foi identificada como as que apresenta maior sensibilidade às mudanças climáticas projetadas pelos modelos globais do CMIP5. The decision-making processes involving water resources in Brazil have been based on the historical behavior of hydroclimatological variables. However, this may lead to inappropriate strategies related to the use of natural resources, since the impact of future climate may change the water resources availability. This research proposes an analysis of the water availability variation of hydrographic basins contained in the contribution area of Patos Lagoon, which can be a consequence of climatic changes predicted by climate models and scenarios of CMIP5 (Coupled Model Intercomparison Project Phase 5), published in the latest International Panel on Climate Change (IPCC) report. Two 30-year future periods were simulated, from 2006 to 2035 and from 2051 to 2080, using the MGB-IPH hydrological model, considering two extreme scenarios of greenhouse gas emissions and twenty climate models. It was possible to verify a trend of increase of the flows in the simulated hydrographic basin, since more than 60% of the simulations performed indicated some percentage of increase of average flows in all the discretized modeling units. The analysis of the simulations results indicated that the data from climatic models HadGEM2-ES and GFDL-CM3 used as input data in the hydrological model are the ones that respectively provide upper and lower flow thresholds of the simulations’ ensemble. The joint evaluation of the results generated by these two models, associated to the scenario of high greenhouse gas emissions, is capable of covering extreme flow scenarios predicted for the future considering climate change in Patos Lagoon basin. The results provided by the models bcc-csm1-1, BNU-ESM and CNRM-CM5 are similar to the median of the set of simulations generated by all the models evaluated in this research. In addition, the northernmost region of the study area was identified as having the highest sensitivity to climate change projected by global models of CMIP5 published in AR5.

Cheias são hoje os processos naturais que mais impactam pessoas no mundo todo, causando mortes e prejuízos estimados em US$104 bilhões por ano. Os seus efeitos são cada vez mais importantes na medida em que a sociedade se desenvolve e se instala em locais passíveis de inundação. Para combater os danos das cheias os sistemas de previsão e de infraestrutura tem aumentado sua complexidade, criando assim demanda para metodologias relativamente complexas de modelagem hidrológica-hidrodinâmica que dependem por sua vez da caracterização adequada das bacias hidrográficas e dos rios. Em muitos casos a aquisição de características da calha dos rios é obtida através de relações geomorfológicas, que para vários casos não apresenta resultados satisfatórios. O presente estudo propõe uma técnica de obtenção de características da seção transversal de rios a partir da análise da relação entre a vazão (Q) e a celeridade (c) de ondas de cheia observadas em trechos de rios, apoiada pela estimativa de largura a partir de imagens de sensoriamento remoto. Esta se baseou na divisão da área superficial do rio pelo comprimento do mesmo. Para o desenvolvimento do primeiro produto (seção transversal) propõe-se a utilização do modelo inercial de propagação de vazões associado ao algoritmo de otimização SCE-UA para replicar as celeridades observadas alterando apenas os parâmetros da calha dos rios. Foram feitos estudos de caso para os rios Piquiri e Iguaçu, no estado do Paraná, Cuiabá no Mato Grosso, Miranda no Mato Grosso do Sul e Piracicaba em São Paulo. Nestes existem dados telemétricos de vazão que permitem extrair a curva cxQ empiricamente. Os resultados gerados mostraram que não apenas foi possível replicar o comportamento da curva cxQ característica nos locais onde ocorre inundação mas também retirar informações de profundidade, largura e coeficiente de manning muito próximos aos observados. Ainda, a largura obtida foi muito próxima àquela resultante da metodologia proposta para sua aproximação via sensoriamento remoto. Estes resultados só não foram observados para o trecho do rio Iguaçu, onde a relação cxQ não apresentou o comportamento de áreas alagadas, ou seja, a diminuição da celeridade a partir de uma vazão de extravasamento. Desta maneira este trabalho demonstrou que é possível a obtenção de um conjunto de parâmetros efetivos da calha de rios a partir apenas de dados de vazão para regiões alagáveis, além de apresentar uma metodologia simples e eficaz para obtenção da largura em qualquer rio. Flooding is the natural processe that impact most people in the world today, causing deaths and losses estimated at $ 104 billion a year. Its effects are increasingly important as the society develops and settles in flood-prone places. To combat flood damage, forecasting and infrastructure systems has increased its complexity, thus creating demand for relatively complex hydrological-hydrodynamic modeling methodologies which in turn depend on the proper characterization of river basins. In many cases, the acquisition of characteristics of river channels is obtained through geomorphological relations but for several cases it does not present satisfactory results. The present study proposes a technique to obtain characteristics of the cross section of rivers from the analysis of the relation between flow (Q) and the wave speed (c) observed in stretches of rivers, supported by the estimation of width from remote sensing images. For the development of the first product (cross-section) it is proposed the use of the inertial flow propagation model associated with the SCE-UA genetic algorithm to replicate observed velocities by changing only the river channel parameters. Case studies were carried out for the Piquiri and Iguaçu rivers in the state of Paraná, Cuiabá in Mato Grosso, Miranda in Mato Grosso do Sul and Piracicaba in São Paulo. In these, there were flow telemetry data that allowed to extract the cxQ curve. The results showed that not only was it possible to replicate the behavior of the characteristic cxQ curve in the places where flooding occurs, but also to extract information of depth, width and manning very close to those observed. Also, the width obtained was very close to the one obtained by the methodological proposal for its approach through remote sensing. This is based on the division of the surface area of the river by its length. These results were not observed for the stretch of the Iguaçu River, where the relation cxQ is not does not behave as in flooded areas, that is, there isn't a decrease of celerity from an extravasation flow. In this way, this work demonstrated that it would be possible to obtain a set of effective parameters of the river channel from flow data for floodable regions. In addition it presents a simple and effective methodology to obtain the width in any river.

A estimativa de valores de vazão de referência máximos, médios e mínimos em locais sem dados constitui um desafio constante nos estudos de hidrologia. Em muitos casos este desafio é enfrentado utilizando as técnicas tradicionais de regionalização de vazões com suas diversas variantes e técnicas. Em outros casos são adotadas as técnicas de modelagem hidrológica chuva-vazão. No presente trabalho é apresentada uma ampla revisão bibliográfica referente aos métodos que são utilizados a nível global em pós da obtenção de dados de vazões em locais sem monitoramento. A metodologia para estimativa de vazões em locais sem dados proposta aqui é baseada na modelagem hidrológica com um modelo chuva-vazão distribuído, combinada com uma técnica de assimilação de dados de vazão observados. São apresentados testes da metodologia utilizando o modelo MGB-IPH para todas as bacias contidas total ou parcialmente no estado do Rio Grande do Sul. Os resultados são avaliados através da comparação entre as vazões de referência obtidas pela metodologia proposta e as vazões de referência calculadas com base nos dados observados em postos fluviométricos que não são utilizados no processo de assimilação de dados e cujos locais são chamados de Pseudo Não Monitorados. As vazões de referência estimadas com base na metodologia proposta também são comparadas com os valores encontrados utilizando equações de regionalização de vazões desenvolvidas para a mesma região e um método simples de interpolação linear. A avaliação dos métodos foi realizada através da quantificação dos erros relativos médios resultantes dos 4 métodos, comparando Vazões Calculadas vs. Vazões Observadas em aqueles Locais Pseudo Não Monitorados. Em geral observa-se que a assimilação de dados contribui na melhora da estimativa das vazões de referência em locais sem dados quando comparados com os outros métodos testados a exceção do caso das vazões Q50 e Q10 onde a diminuição do ERM é mínima. Os resultados obtidos nesta pesquisa encorajam o avanço nas pesquisas dentro dessa linha seja aprimorando o próprio método testado aqui ou utilizando outras metodologias de assimilação no acoplamento ao modelo hidrológico MGB-IPH e outros. The estimate of maximum, average and minimum reference flow in places without observed data is a constant challenge in studies of hydrological engineering. In many cases this challenge is faced using traditional techniques of regionalization of flow with its various variants and techniques. In other cases, the techniques of hydrologicalrainfall–runoff modelling are adopted. In the present work a broad bibliographic review is presented referring to the methods that are used at global level in post of obtaining flow data in ungauged catchments. This research presents an empirical data assimilation method applied within the MGB-IPH model to minimize the errors in the estimation of streamflow, performed through hydrologic simulations in the Rio Grande do Sul State of Brazil. The results are evaluated by comparing of the reference flow rates obtained by the proposed methodology and reference streamflows calculated based on data observed at fluviometric stations, called Pseudo-Ungauged Locations, not used in data assimilation process. The reference streamflow estimated based on the proposed methodology are also compared with the values found using regionalization equations flows developed for the same region and a simplelinear interpolation method. The evaluation of the methods was carried out by quantifying the Average Relative Errors resulting from the three methods, comparing calculated flows vs. observed flows in those sites called Pseudo-Ungauged Locations. In general it is observed that data assimilation contributes to the improvement of the estimation of the reference flows in places with no data when compared with the other methods tested, except for the case of the flows Q50 and Q10 where the decrease of ERM is minimal. The results obtained in this research encourage the advancement of research within this line either by improving the method used here or by using other assimilation methodologies in the MGB-IPH or other hydrological model.

Neste estudo, trata-se de avaliar fontes de precipitação baseadas em estimativas por satélite e técnicas de assimilação de dados para previsão de vazões por meio do modelo hidrológico distribuído MGB-IPH. A insuficiente representatividade espacial dos pluviômetros torna difícil a correta representação dos campos de precipitações. Por outro lado, as estimativas de satélite, embora forneçam uma descrição espacial mais consistente, são potencialmente menos acuradas. Sendo assim, procura-se utilizar métodos que combinem os dados de ambas as fontes para gerar um campo de precipitação mais consistente. Neste trabalho, implementaramse dois modelos de combinação pluviômetro-satélite, CHUVSAT e MERGEHQ, através de uma metodologia de interpolação. Por outro lado, as técnicas de assimilação de dados acoplados aos modelos de previsão hidrológica são também de interesse neste estudo, pois minimizam as incertezas associadas ao processo de calibração de parâmetros, às variáveis de estado e dados de entrada do modelo hidrológico. Para esse propósito, escolheu-se a bacia do rio Tocantins e implementou-se particularmente a técnica de assimilação de dados de tipo sequencial chamado na literatura de filtro de partículas, conjuntamente com o método de filtro Kalman por conjunto e o método de assimilação AsMGB atualmente acoplado ao modelo MGB-IPH. O estudo mostra que a precipitação combinada utilizada como dado de entrada na simulação hidrológica permitiu reproduzir adequadamente os hidrogramas observados para o período de calibração e validação. Já para o caso das vazões resultantes, durante a etapa de previsão, a precipitação combinada mostrou-se com melhor desempenho em termos estatísticos que os métodos sem combinar, sobretudo após 24 horas de antecedência. Finalmente, a técnica de assimilação de dados por filtro de partículas conseguiu absorver os erros da simulação melhorando as medidas de desempenho na etapa de previsão sendo superior ao modelo de previsão sem considerar assimilação. The objective of this study is to evaluate precipitation sources based on satellite estimates and data assimilation techniques for prediction of flows by means of the distributed hydrological model MGB-IPH. The insufficient spatial availability of rain gauges makes difficult to represent precipitation fields appropriately. In contrast, satellite estimates, although providing a more consistent spatial description, are potentially less accurate. Thus, raingauge satellite merging methods that combine data from both sources to generate a more consistent precipitation field are used herein. For this purpose, two models namely CHUVSAT and MERGEHQ were implemented using an interpolation technique. On the other hand, data assimilation techniques coupled with hydrological forecasting models are also assessed in this study. The assimilation process minimizes the uncertainties associated with the parameter calibration procedure, variable state and hydrological input data. In this manner, the sequential data assimilation technique namely particle filter in conjunction with the Kalman filter method and the assimilation method AsMGB, which is currently coupled to the MGBIPH model, were implemented and applied to the Tocantis basin. The obtained results showed that the combined precipitation used as input data in the hydrological simulation allowed reproducing adequately the observed hydrograms for the periods of calibration and validation. In the case of the resulting flows during the forecast stage, the merging precipitation was shown to perform better in statistical terms than the uncombined methods, especially after 24 hours in advance. Finally, the data assimilation technique by particle filter was able to absorb all simulation errors, improving the performance measures in the forecasting stage, thus being superior to the forecasting model without considering assimilation.

Grandes áreas inundáveis prestam importantes serviços ambientais, como manutenção da biodiversidade e regulação de clima, cheias e ciclos biogeoquímicos, e estão presentes nas mais variadas regiões geográficas, desde úmidas áreas tropicais como a Amazônia, até locais áridos como os grandes deltas internos da África. Complexos processos hidrológicos e hidrodinâmicos ocorrem nestas áreas, onde existe intensa interação entre áreas inundadas, biosfera e atmosfera, e onde ocorrem alterações significativas das ondas de cheia que são por elas propagadas sazonalmente. Neste contexto, o presente trabalho apresenta, através de dois métodos distintos, formas com que podemos aprofundar nosso conhecimento da dinâmica destas grandes áreas inundáveis. Inicialmente, é abordado o efeito causado pela interação rio-planície de inundação na forma dos hidrogramas. Constatou-se, através de observação de diversos hidrogramas defluentes de grandes sistemas inundáveis, que a ascensão destes é tipicamente mais lenta que a recessão, caracterizando uma assimetria negativa. Então, a partir de uma série de abordagens numéricas, analíticas e empíricas, mostrou-se que a este fenômeno é decorrente da inversa relação entre velocidade da onda de cheia (celeridade) e vazão que ocorre nestes ambientes devido ao armazenamento nas planícies. Um estudo de caso com os principais rios da Amazônica evidenciou a relevância da avaliação de assimetria de hidrogramas para compreensão da dinâmica hidrológica em regiões inundáveis. A segunda etapa deste trabalho consistiu no aprimoramento e validação de um modelo acoplado hidrológico-hidrodinâmico distribuído para simulação de bacias com grandes áreas inundáveis, capaz de modelar matematicamente a dinâmica existente e a aumentar nossa capacidade preditiva da inundação destas regiões. Foi realizada uma aplicação na bacia do Alto Rio Níger, no Oeste da África, onde está localizado o Delta Interno do Níger, uma vasta planície no Deserto do Sahel inundada sazonalmente, e onde importantes perdas por evaporação ocorrem. O modelo desenvolvido representou de forma satisfatória vazões, níveis e áreas inundadas ao longo da região, permitindo a identificação de processos relevantes na inundação do Delta. O acoplamento entre uma representação explícita da hidrodinâmica dos canais da planície (armazenamento, bifurcações, remanso) e o balanço vertical hidrológico (infiltração de água da planície para o solo e evapotranspiração) mostrou-se necessário para o satisfatório resultado do modelo. Assim, com esta pesquisa, acredita-se que novos elementos foram contribuídos para a compreensão da hidrologia de grandes áreas inundáveis, a partir da identificação do importante efeito da interação rio-planície na assimetria de hidrogramas, e do desenvolvimento de técnicas de modelagem matemática que evidenciaram a relevância da representação de interações entre hidrologia e hidrodinâmica em modelos de grandes escala. Large floodable areas provide important environment services, such as biodiversity maintenance and regulation of climate, floods and biogeochemical cycles, and are present throughout different geographical settings, from tropical, wet areas like the Amazon, to arid ones such as the large African Inner deltas. Complex hydrologic and hydrodynamic processes occur in these areas, such as an intense interaction between flooded areas, biosphere and atmosphere, and alterations of flood waves that propagate through it. In this context, this study presents, through two distinct methods, ways with which we can deepen our understanding of the dynamics of such floodable areas. Firstly, the effect of river-floodplain interaction on hydrograph shape is addressed. It was noted, by observing a series of hydrographs routed through floodable systems, that hydrograph ascension is typically slower than the recession, characterizing a negative skewness. Then, through a series of numerical, analytical and empirical approaches, it was shown that this phenomenon exists due to the inverse relationship between flood wave speed (celerity) and discharge that occurs in these environments because of floodplain storage. A case study in main Amazon Basin tributaries was performed and indicated the relevance of hydrograph skewness evaluation for comprehension of hydrological dynamics in floodable areas. The study second step consisted in the development and validation of a distributed, coupled hydrologic-hydrodynamic model for simulation of basins with large floodable areas, able to mathematically model and to improve our prediction capacity of flooding in these areas. An application was performed in Upper Niger River, West Africa, where the Niger Inland Delta is located, which is a vast floodplain in the Sahel Desert that is seasonally flooded, and where important evaporation losses occur. The developed model satisfactorily represented discharges, levels and flooded areas across the basin, allowing the identification of relevant flooding processes in the Delta. The coupling between explicit representation of floodplain channel hydrodynamics (storage, bifurcations, backwater) and vertical hydrological balance (floodplain water infiltration into soil column, and evapotranspiration) proved a necessary model structure for the satisfactorily obtained results. Finally, with this research, we understand that new elements were contributed to the comprehension of hydrology of large floodable areas, from the identification of the river-floodplain interaction effects on hydrograph skewness to the development and application of mathematical modelling techniques, showing the relevance of the interaction between hydrology and hydrodynamic in such areas.

A bacia hidrográfica do rio da Prata (BHP) é a quinta maior bacia hidrográfica do mundo em termos de área de drenagem e a segunda maior da América do Sul. Os recursos hídricos da bacia têm um papel fundamental em termos de economia, produção de energia hidrelétrica, navegação, agropecuária dos países que a compõe e está fortemente relacionado com problemas relacionados à eventos extremos. Apesar da importância hidrológica da BHP para a América do Sul, poucas pesquisas estudaram o seu comportamento hidrológico de forma integrada e utilizando representações mais físicas dos processos hidrológicos. Essas pesquisas, por exemplo, possuem limitações relacionadas a representação do processo de propagação de vazão, escoamento na planície de inundação do modelo hidrológico utilizado e ao fato de não considerar o efeito da operação dos reservatórios, que é um aspecto importante na hidrologia da BHP. Diante disso, o objetivo dessa tese foi desenvolver bases técnicas para a modelagem hidrológica integrada da bacia do rio da Prata, e mais especificamente, o aperfeiçoamento de um modelo hidrológico para esta bacia que possa servir como ferramenta comum para o apoio à gestão integrada de recursos hídricos. Nesse contexto, foi utilizado o Modelo Hidrológico de Grandes Bacias, MGB-IPH, com um aperfeiçoamento do modulo de propagação de vazões, que corresponde a utilização do modelo Inercial e a representação da planície de inundação. Além disso, a tese também apresentou um avanço nas técnicas de discretização da bacia hidrográfica em unidades menores, visando a possibilidade de representação de regiões muito planas e com drenagem complexa. Essa tese também apresentou uma análise acerca do impacto da operação dos reservatórios do Alto Paraná em algumas vazões de referência do baixo rio Paraná. Os resultados da tese foram divididos em cinco partes: Calibração e verificação; efeito do modelo de propagação utilizado nas vazões; verificação das áreas inundadas; Avaliação do impacto dos reservatórios nas vazões do Baixo Paraná; análise de sensibilidade. Os resultados dessa Tese mostraram que é possível, através de uma modelagem integrada, utilizando a versão aperfeiçoada do modelo MGB-IPH, representar os processos hidrológicos da Bacia do Prata. Apesar das incertezas apresentadas, o modelo conseguiu representar vazões, níveis e áreas inundadas de forma satisfatória na maioria dos pontos analisados. Os resultados da tese também mostraram que o impacto da operação dos reservatórios pode ser significativo, principalmente nas vazões extremas. The La Plata River Basin (LPB) is the fifth largest basin in the world and the second one largest in South America in terms of drainage area. The water resources of LPB have an important role for the economy, hydroelectric power, navigation, agriculture of the countries that compose the basin and they are strongly related with hydrologic extreme events. Despite the importance of LPB’s hydrology to South America, there are few researches that have studied the hydrologic behavior in an integrated manner and using more physical representations of hydrologic processes. These researches have limitations related to representation of the flow propagation processes, floodplain representation and no consideration of the effect of reservoir operation, which is an important aspect of LPB. Thus, the aim of this research was developing technical basis for the integrated hydrologic modeling of the LPB improving a hydrologic and hydraulic model for this basin. This model can be useful as common tool to support the integrated management of LPB’s water resources. We used the Large Scale Hydrologic Model, MGB-IPH, with an improvement of the flow propagation module, that include the use of Inertial flow propagation model and a simplified scheme to floodplain representation. Moreover, the research also introduced an advance in techniques of river basin discretization into smaller units, aiming at the possibility of representation of very flat regions with complex drainage. This research also presented an analysis of the impact of reservoirs’ operation of Alto Paraná river basin. The results were divided into five parts: calibration and verification; effect of flow propagation model; check the flooded areas; Impact of reservoirs in the Lower Paraná flows and sensitivity analysis. The results showed that it is possible through an integrated modeling, using the improved version of the MGB-IPH model, represent the hydrological processes of the Plata Basin. Despite the uncertainties presented, the model could represent flows, flooded areas and levels satisfactorily in most analyzed points. The results also showed that the impact of the reservoirs operation can be significant, especially in extreme flows.

Mudanças climáticas podem afetar a distribuição espacial e temporal das variáveis hidrológicas, tendo como consequências alterações nos regimes de precipitação e vazão dos rios. Aumentos ou reduções no volume de escoamento de uma bacia hidrográfica podem, por exemplo, produzir danos aos ecossistemas, afetar a produção de alimentos, abastecimento de água, navegação e geração de energia. Atualmente buscam-se relações que permitam entender os processos de mudanças climáticas a fim avaliar os impactos e mitigá-los, assim como avaliar as incertezas inerentes ao processo de modelagem hidrológica de tais mudanças. Neste contexto este trabalho apresenta uma metodologia de quantificação e análise de incertezas para estudos de mudanças climáticas, tomando como estudo de caso a bacia hidrográfica do Rio Uruguai (BHRU) com área aproximada de 110,000 Km². Para tanto três fontes de incerteza foram analisadas e comparadas: o modelo hidrológico, técnicas de remoção de viés e modelos climáticos. O modelo hidrológico MGB-IPH foi avaliado quanto ao processo de parametrização, utilizando diferentes períodos de simulação para calibração: (i) Período de calibração 1 – MGB/P1: representando a série completa de observações de 1960-1990 com verificação no período de 1992-1999; (ii) Período de calibração 2 - MGB/P2: calibração em período seco e verificação de período de cheias (iii) Período de calibração 3 – MGB/P3: calibração em período característico de cheias e verificação de período de estiagem. Três diferentes técnicas de remoção de viés foram aplicadas para analisar o grau de incerteza que a escolha de um determinado método de correção pode agregar ao resultado final: (i) RV1 - Técnica de Mapeamento Quantil-Quantil; (ii) RV2 - Técnica de Escalonamento Linear e (iii) RV3 - Técnica Delta change. Os modelos climáticos globais (GMC’s) foram analisados quanto a sua estrutura, comparando projeções de cinco diferentes modelos: MPEH5 (ECHAM5/MPIOM), GFCM21 (GFDL-CM2.1), MRCGCM (MRI-CGCM2.3.2), HADCM3 (UKMO-HadCM3) e NCCCSM (CCSM3). Adicionalmente, também foram analisadas as projeções climáticas de cinco diferentes versões do modelo climático regional (RCM) ETA/CPTEC: CT20, CT40, LOW, MID e HIGH. Inicialmente os resultados das simulações provenientes de cada uma destas fontes foram comparados de maneira isolada e em seguida de maneira combinada. Portanto, a metodologia foi dividida em Etapa (1) e Etapa (2). A Etapa (1) teve por objetivo responder a seguinte pergunta: Qual dentre as fontes de incerteza selecionadas agrega maior variação ao resultado final? Ou seja, qual destas fontes propaga maior incerteza em termos de impactos de mudanças climáticas na BHRU? Os resultados obtidos por cada uma das fontes de incerteza foram comparados em termos de anomalias de vazões médias de longo período (QMLP), máximas e mínimas anuais. Na Etapa (2) foi realizada a análise total de incerteza, ou seja, a análise combinada dos resultados obtidos na Etapa (1). As anomalias de vazões foram apresentadas utilizando as curvas de distribuição acumulada (CDF’s) e a incerteza total expressa pela diferença entre os percentis 5% e 95%. Considerando os resultados obtidos para as vazões médias de longo período (QMLP), as fontes podem ser ordenadas de forma decrescente, em relção ao grau de incerteza que propagam: modelos climáticos globais > modelos climáticos regionais > técnicas de remoção de viés > modelo hidrológico. Para as vazões extremas os RCM’s apresentam as maiores variações de anomalias se comparadas às dos modelos hidrológicos e técnicas de remoção de viés, inclusive para ambos os extremos, máximos e mínimos. Esta variação se dá principalmente, pelos resultados de LOW e MID. Estas informações podem ajudar os gestores e tomadores de decisão no adequado gerenciamento e planejamento dos recursos hídricos sob condições de mudanças climáticas, assim como o entendimento da incerteza associada. Climate change can affect the spatial and temporal distribution of hydrological variables, with the consequences of changes in precipitation regimes and river flows. Increase or decrease the flow of rivers, for example, can cause damage to ecosystems, affecting food production, water supply, navigation and power generation. Currently seeking to relationships that allow understand climate change processes in order to assess the impacts and mitigate them, and assess the uncertainties inherent in hydrologic modeling process of such climate change. This thesis aimed at the development of a methodology for quantification and analysis of uncertainties for climate change studies in hydrology , taking as a case study the basin of the Uruguay River (BHRU) with a drainage area near 110,000 km². For that three sources of uncertainty were analyzed and compared: the hydrologic model, bias removal techniques and climate models. The hydrological model MGB-IPH was evaluated for parameterization, using different simulation periods for calibration: (i) MGB /P1: full range with calibration period (1960-1990) and validate (1992-1999); (ii) MGB / P2: calibrated in the period of dry and validated in the flood season (iii) MGB/P3: calibrated in the period of floods and validated in the dry season. Three different bias correction methods were applied to analyze the degree of uncertainty that the choice of a particular method of correction can add to the final result: (i) RV1 - Quantil-Quantil Mapping; (ii) RV2 - Linear Scaling, and (iii) RV3 - Delta Change Technique. Global climate models (GMC's) were analyzed for their structure, comparing projections of five different models: MPEH5 (ECHAM5/MPI-OM), GFCM21 (GFDLCM2.1), MRCGCM (MRI-CGCM2.3.2), HADCM3 (UKMO-HadCM3) e NCCCSM (CCSM3). Additionaly, climate projections from five different versions of the regional climate model (RCM) ETA / CPTEC were also analyzed: CT20, CT40, LOW, MID e HIGH. Initially the simulation results from each of the sources of uncertainty were compared individually (single propagation) and then in a combined way (multi propagation). Therefore, the methodology was divided in Step (1) and Step (2). Step (1) aimed to answer: Which of the selected sources of uncertainty adds more variation to the final result? Which of these sources propagates greater uncertainty in terms of impacts of climate change on BHRU? The results for each of the sources of uncertainty were compared in terms of long-term mean flow (QMLP), maximum and minimum annual flow. In Step (2) total uncertainty analysis was performed, therefore the combined analysis of the results obtained in Step (1). The anomalies in discharge were presented using the cumulative distribution function (CDF's) and the total uncertainty expressed by the difference between the percentiles 5% and 95%. Throughout the application of the proposed methodology it was concluded that: (i) for the extremes (maximum and minimum) annual discharges the largest source of uncertainty are the projections of the RCM's, followed by the the bias removal technique and finally the hydrological model; (ii) for the QMLP the largest source of uncertainty are followed global climate models, then the regional climate models. This information can help managers and decision makers in the proper management and planning of water resources under climate change conditions, as well as the understanding of the associated uncertainty.

O tempo que a água fica retida nos rios e planícies de inundação afeta diretamente os processos biogeoquímicos e ecológicos relacionados com os ciclos hidrológicos. Apesar de existirem métodos simplificados e sofisticados para estimar o tempo de viagem e o tempo de residência, a análise de grandes sistemas hidrológicos (e.g. maior do que 104 km2) e de sua variabilidade no tempo e no espaço ainda é um desafio. Neste trabalho foi desenvolvido um modelo numérico de rastreamento de águas, numa abordagem probabilística, capaz de fornecer informações associadas à origem, trajetória e dinâmica das escalas de tempo de renovação de água superficial, com enfoque em bacia hidrográficas. Experimentos numéricos foram utilizados para estudar a dinâmica do escoamento de água na bacia Amazônica, com ênfase nas interações e trocas de água entre rio e planícies de inundação. Devido à importância e magnitude dessa bacia, a presença de grandes áreas inundáveis (17% da bacia) e águas de características distintas associadas a sua origem (i.e. águas brancas, claras e negras), a bacia Amazônica configura um estudo de caso interessante para simulações baseadas em rastreamento. Os resultados demonstraram que a planície de inundação tem um efeito de armazenamento importante para a dinâmica da bacia. As simulações sugerem que: (i) o tempo de viagem da água no canal principal, em sua maioria, é menor do que 70 dias; (ii) cerca de 20% das águas da bacia passam pela planície de inundação; (iii) a vazão de troca entre rio e planície de inundação pode chegar a ordem de 105 m3/s, para a bacia inteira; (iv) se as planícies de inundação atuam como áreas de armazenamento durante a cheia, elas podem acarretar um aumento substancial no tempo de viagem das águas (i.e. maior do que 200 dias); (v) e, nesse caso, águas retidas na planície de inundação podem representar uma grande parcela do escoamento da bacia durante o período de águas baixas. Por fim, outros aspectos sobre a variabilidade espacial e temporal do tempo de viagem das águas superficiais e da composição da água, em relação a sua origem, foram estudados. The time that water is retained in rivers and floodplains directly affects the biogeochemical and ecological processes associated to the hydrological cycles. Although there are both simplified and sophisticated methods to estimate travel time and residence time, the analysis of large hydrological systems (e.g. larger than 104 km2) and their variability in time and space is still a challenge. In this work, we developed a numerical model for water-packet tracking, with a probabilistic approach, capable to provide information associated to the origin, trajectory and dynamics of the surface water renewal time scales, focusing application to watersheds. Numerical experiments were used to study the dynamics of water flow in the Amazon basin, with emphasis on interaction and water exchanges between river and floodplains. Due to the importance and magnitude of this basin, the presence of large floodable areas (17% of the basin) and waters of distinct characteristics associated with its origin (e.g. white waters, clear waters and black waters), the Amazon basin constitutes an interesting case study for simulations Based on tracking. The results showed that the floodplain has an important storage effect for the dynamics of the basin. The simulations suggest that: (i) the travel time of the water in the main channel is less than 70 days; (ii) about 20% of the waters of the basin pass through the floodplain; (iii) the flow between river and floodplain can reach the order of 105 m3/s at basin scale; (iv) if the floodplains act as storage areas during the flooding period, they could increase substantially the travel time of the basin waters (i.e. larger than 200 days); (v) also, in this case, water retained in the floodplains could represent a large portion of the river discharges during the low water period. In addition, other aspects on the spatial and temporal variability of the travel time and composition of surface waters were also explored.

Fundação de Amparo a Pesquisa no Estado do Rio Grande do Sul The sediment yield causes a number of problems, and among them stands out silting of rivers and reservoirs, and contamination of surface water. Sediments result from the processes of erosion, transport, storage and compression. They depend on the physical characteristics of the watershed, as well as soil use and land cover. To study the sediment yield of a river is important to estimate the amount and origin of the sediments transported and deposited in a watershed. Therefore, it is possible to obtain knowledge about the unit of study, which leads to better planning and use of water resources. From liquid and solid discharge measurements in a river section, and conducting a sedimentometrical balance in Vacacaí Mirim reservoir, located at Santa Maria, Rio Grande do Sul, Brazil, the objective of this research is to characterize the sediment yield of a small watershed. The methodology was based on two stages: the first one was about the sedimentometrical characterization of the study area, and it was carried out between January and November of 2014, with the continuous monitoring of flow and sediments in fluviossedimentometrical station João Goulart, whose laboratory samples were also characterized and had their sediment concentration determined; the second stage was about the execution of sedimentometrical balance in the reservoir with the help of HEC-HMS 4.0 modeling program. Sediment entries, storage and output were checked in Vacacaí Mirim reservoir, which influences the study area. The results for the hydrossedimentometrical monitoring were satisfactory, getting a key curve with determination coefficient of 0.96; sedimentometrical curves ranged from 0.93 to 0.98; concentration ratio between suspended sediment and turbidity was 0.84, being considered acceptable; particle size characteristics were predominantly silt, for suspended material, and sand, for dragged material. Balance calibration for sedimentometrical hydrologic model gave a Nasch-Sutcliffe between 0.13 to 0.53 and 0.44 in validation; about reservoir sediments balance, it was noted and therefore considered that all sediment that enters the reservoir is retained and more sediment is produced downstream thereof mouth of the river it reaches the watershed. A produção de sedimentos ocasiona uma série de problemas, entre eles destacam-se assoreamento dos rios e reservatórios e contaminação das águas superficiais. Os sedimentos são resultado dos processos de erosão, transporte, depósito e compactação. Dependem das características físicas da bacia hidrográfica e do uso e ocupação do solo. Estudar a produção de sedimentos em um rio é importante para estimar a quantidade e a origem dos sedimentos transportados e depositados em uma bacia hidrográfica, podendo assim, obter conhecimento da unidade de estudo e realizar um melhor planejamento e aproveitamento dos recursos hídricos. O objetivo desta pesquisa é de caracterizar a produção de sedimentos de uma pequena bacia hidrográfica, a partir de medições de descarga líquida e sólida na seção do rio e realizar um balanço sedimentométrico no reservatório do Vacacaí Mirim no município de Santa Maria. A metodologia se baseou em duas etapas: primeira etapa foi à realização da caracterização sedimentométrica da área de estudo com o monitoramento contínuo das vazões e sedimentos na estação fluviossedimentométrica João Goulart no período de janeiro a novembro de 2014 e na determinação da concentração de sedimentos e caracterização das amostras em laboratório; a segunda etapa foi à realização do balanço sedimentométrico no reservatório com o auxílio do programa de modelagem HEC-HMS 4.0 onde foram verificados as entradas, armazenamento e saída de sedimentos no reservatório Vacacaí Mirim que influencia a área de estudo. Os resultados obtidos para o monitoramento hidrossedimentométrico foram satisfatórios, obtendo uma curva-chave com coeficiente de determinação de 0,96; as curvas sedimentométricas variando de 0,93 a 0,98; a relação da concentração de sedimentos em suspensão e a turbidez foi de 0,84 sendo consideradas aceitáveis; as características granulométricas do material em suspensão foi predominante o silte e para o material em arraste foi à areia. Para o balanço sedimentométrico a calibração hidrológica do modelo obteve-se um Nasch-Sutcliffe entre 0,13 a 0,53 e na validação de 0,44; as considerações observadas no balanço de sedimentos no reservatório foram que todo o sedimento que entra fica retido no reservatório e é produzido mais sedimentos a jusante do mesmo que chega até o exutório da bacia.

A previsão hidrológica possibilita a identificação antecipada de eventos de cheia potencialmente causadores de inundação, o que é de grande importância para a atuação de entidades como a Defesa Civil. Quando se deseja estender a antecedência no tempo em relação a estes eventos, principalmente nos casos onde a bacia de interesse é relativamente rápida, torna-se necessária a incorporação de previsões quantitativas de precipitação (QPF) na modelagem hidrológica, as quais podem ser obtidas a partir de modelos numéricos de previsão do tempo. Entretanto, a falta de acurácia atribuída a estas previsões de chuva, dadas de forma determinística, vem promovendo sua substituição por sistemas de previsão meteorológica por conjunto (EPS - Ensemble Prediction Systems), cuja finalidade é a geração de possíveis estados futuros da atmosfera para considerar as incertezas associadas ao seu estado inicial e às deficiências na representação física dos modelos de previsão do tempo. Neste contexto, o presente estudo teve por objetivo avaliar uma metodologia de previsão de cheias por conjunto na bacia do Taquari-Antas/RS até a cidade de Encantado (19.000 km²), localizada na região Sul do Brasil. Para tanto, foi utilizado o modelo hidrológico MGB-IPH acoplado a diferentes sistemas de previsão, sendo eles: (i) EPS Regional ETA, de curto prazo (até 72 horas) com 5 membros de diferentes parametrizações; (ii) EPS Global ECMWF de médio prazo (até 10 dias) com 50 membros de condições iniciais perturbadas, incluindo perturbação estocástica nos parâmetros de ajuste do modelo e; (iii) Previsão Determinística do Modelo Regional ETA (até 7 dias). A avaliação das previsões consistiu em dois hindcastings distintos, envolvendo uma análise visual de eventos singulares ocorridos em 06/06/2014 e 21/07/2011 além de uma análise estatística no período de Mar/2014 - Nov/2014. Durante a análise visual foi possível identificar, a partir de antecedências de 5 a 6 dias, uma persistência na previsão dada pelo crescente número de membros do conjunto de médio prazo (ECMWF) com superação dos limiares de referência, na medida em que se aproximavam os eventos de cheia. Apesar da grande incerteza na magnitude das previsões hidrológicas para o conjunto de curto prazo, a vazão máxima foi relativamente bem prevista por pelo menos 1 membro em quase todas as antecedências, enquanto que a previsão do timing dos eventos foi considerada de boa confiabilidade. Durante a avaliação estatística foi possível notar uma falta de espalhamento nos conjuntos, com tendência de subestimativa de acordo com o aumento da antecedência. Em uma comparação com previsões determinísticas, as previsões por conjunto demonstraram maior acurácia principalmente até 72 horas de antecedência, com destaque para a maior probabilidade de detecção dos limiares de referência e manutenção de falso alarme a níveis reduzidos. Além disso, verificou-se também que a agregação de previsões efetuadas em tempo anterior naquelas atuais acarreta em ligeira ampliação do espalhamento do conjunto e maiores probabilidades de detecção dos limiares de alerta para os membros mais elevados, apesar da redução no desempenho em termos de acurácia e viés. De forma geral, as previsões por conjunto apresentam potencial para servir como uma informação complementar em sistemas de alerta contra cheias, possibilitando uma melhor preparação dos agentes envolvidos durante a ocorrência destes eventos. Hydrological forecasting plays an important role for issuing flood warnings, allowing for anticipation and better preparation of authorities at the occurrence of such events. In order to extend lead time in a flood forecast, especially when the catchment response time is relatively fast, it may be useful to couple a hydrological model to quantitative precipitation forecasts (QPF), usually obtained directly from numerical weather prediction (NWP). However, deterministic (i.e. single) QPF are usually referred to many errors and lack of accuracy, mainly caused by uncertainties on initial state of the atmosphere and on physical representation of weather forecasting models. To address these shortcomings, it becomes necessary to take into account the uncertainties associated to rain forecasts, which can be represented by Ensemble Prediction Systems (EPS). The purpose of such systems is to provide different trajectories of the atmosphere by perturbations on its initial condition and on parameterization schemes of the models, generating an ensemble of forecasts that can be used as input to hydrological modelling (HEPS). In this context, the present study aimed to assess a methodology of ensemble flood forecasting on Taquari-Antas basin up to the city of Encantado/RS (19.000 km²), located in southern Brazil. Therefore, the hydrological model MGB-IPH was coupled to different forecasting systems: (i) Short Range EPS ETA (up to 72 hours), a regional model with 5 members of different parameterization schemes; (ii) Medium Range EPS ECMWF (up to 10 days), a global model with 50 members of perturbed initial conditions and stochastic perturbation in the model parameters; (iii) Deterministic ETA Model (up to 7 days). The forecasts were evaluated by two different hindcastings, which includes a visual assessment of singular events occurred on 2011 and 2014 and a statistical analysis for the period between Mar/2014 and Nov/2014. It was possible to identify a forecast persistence on medium-range for the selected events, by the increasing number of members exceeding the reference thresholds from lead times up to 5 - 6 days. On the short range, although large uncertainties in the magnitude of hydrological forecasts were found, the peak discharge was well forecasted - at least for a single member - in nearly all lead times, whereas the prediction of the peak timing was considered reliable. Regarding to statistical evaluation, an inadequate spread in the ensemble was observed from short- to medium-range, with a tendency of underestimation for longer lead times. In a comparison with deterministic forecasts, the ensemble forecasts showed higher accuracy especially up to 72 hours in advance, including highlights on greater probability of detection (POD) above the reference thresholds even with low false alarm rates. It also was found that the a combination of previous forecasts on the recent ones leads to a slight increase of ensemble spread and POD for higher members, despite the performance reduction in terms of accuracy and bias. In summary, the hydrological ensemble forecasts demonstrated a good potential to serve as an additional information within a Flood Alert System.

Foi desenvolvida uma ferramenta para simular a geração, o transporte e a redistribuição de sedimentos em bacias de grande escala, através da integração de uma metodologia simplificada para estimativa da distribuição espacial e temporal da perda de solo na bacia com uma metodologia de transporte de sedimentos em rios com planícies de inundação, e com o modelo hidrológico distribuído MGB-IPH. O MGB-IPH é um modelo hidrológico distribuído de grandes escalas, baseado em processos, que adota uma discretização da bacia em minibacias e o conceito de Unidades de Respostas Hidrológicas (URH). O modelo usa equações de base física para simular os processos hidrológicos, como o modelo de evapotranspiração de Pennan Monteith, e utiliza a abordagem de Muskingum-Cunge e um modelo hidrodinâmico 1D para propagar as vazões nos rios, incluindo efeitos de remanso e inundações sazonais. A integração considera a possibilidade de aplicação em bacias com escassez de dados e o uso de utilização de Sistemas de Informações Geográficas e dados espacializados. Buscou-se representar a dinâmica dos sedimentos em grandes bacias, tendo como estudo de caso a bacia do rio Madeira, caracterizada pela sua grande extensão territorial (~1,4.106 km2) e elevada produção de sedimentos. A MUSLE é utilizada para estimar a perda de solo em cada minibacia considerando um fator topográfico LS bidimensional extraído do SRTM/MDE e os demais parâmetros obtidos da literatura ou por formulações específicas. O transporte de sedimentos nos rios é realizado utilizando uma equação de advecção para partículas de silte e argila (materiais em suspensão) e uma equação da continuidade para partículas de areia (carga de fundo). Uma troca lateral de sedimentos em suspensão também é permitida entre os rios e suas planícies de inundação. As metodologias acopladas ao MGB forneceram resultados satisfatórios na bacia do rio Madeira, possibilitando uma visão distribuída da localização das principais fontes de geração de sedimentos na bacia e dos trechos de rios com maiores cargas de sedimentos transportadas. A dinâmica geral dos sedimentos, desde os principais locais de geração de sedimentos na bacia até a sazonalidade do transporte de sedimentos em suspensão nos principais rios, foi representada de forma distribuída. As concentrações e as cargas de sedimentos em suspensão simuladas foram comparadas com dados observados, nas escalas anual, mensal e diária, com ajustes satisfatórios entre elas ao longo dos principais rios da bacia. A transferência lateral de sedimentos entre rio e planície indicaram que 35% de toda a carga de sedimentos provenientes dos Andes é depositada nas planícies. Essa representação também permitiu o amortecimento dos picos e o retardo da propagação das concentrações. As concentrações e as cargas médias mensais de sedimentos em suspensão apresentaram comportamento semelhante aqueles descritos por outros estudos realizados neste bacia, os quais a avaliaram a partir de dados de campo. Apesar dos resultados satisfatórios, foi feita uma avaliação das potencialidades e limitações do modelo e fornecidas recomendações para trabalhos futuros. It was developed a tool for modelling the production, transport and redistribution of sediments for large scales river basins, by integrating a simplified methodology to predict the spatial and temporal distribution of the soil loss over the basin with a methodology to transport sediments along river with floodplains, and couple them to the MGB-IPH hydrological distributed model. The MGB-IPH is a large-scale, distributed and process based hydrological model that uses a catchment based discretization and the Hydrological Response Units (HRU) approach. It uses physical based equations to simulate the hydrological processes, such as the Penman Monteith model for evapotranspiration, and uses the Muskingum Cunge approach and a full 1D hydrodynamic model for river routing; including backwater effects and seasonal flooding. The integration considers possible applications of the model in basins with data scarcity and the use of geographic information systems and spatially distributed information. It aimed to represent the sediments dynamics of large basin using the Madeira river basin as a case study, which is a large basin (~1,4.106 km2) with high amount of sediment transportation. The MUSLE model is used to predicts basin sediment production et each catchment considering a bi-dimensional topographic LS factor extract from SRTM/DEM and values of the other parameters derived from literature or specific equations. Sediment transport along the river are performed using an advection equation for silt and clay (suspended materials) and a sediment continuity equation for sand (bed load). A lateral exchange of suspended sediments is also allowed between rivers and its floodplains. The methodologies coupled with the MGB-IPH model provided satisfactory results in the Madeira river basin, allowing a distributed view of sediment production sites over the basin and of the rivers reaches transporting large amount of sediments. General sediment dynamic was represented in a spatialized way, including the main sources of sediment productions and the suspended sediment transport seasonality. Concentrations and loads of simulated suspended sediment were compared with observed data at annual, monthly and daily scales, shown satisfactory adjustments between them along the main rivers of the basin. Lateral exchanges of suspended sediment between rivers and floodplains indicated that near 35% of all Andeans loads are deposited into floodplains. This representation also allowed the damping of the peaks and the routing delay of the concentrations and of the loads of suspended sediments. Mean monthly concentrations and loads of suspended sediments showed similar behaviors of those described by other researches at this basin using observed data. Despite the satisfactory results, an assessment was made of the potential and limitations of the model and provided recommendations for future work.

A previsão com antecedência de curto e médio prazo da vazão em diferentes locais de bacias hidrográficas geralmente é benéfica ao permitir uma resposta antecipada a eventos hidrológicos como cheias, e a operação mais eficiente de obras hidráulicas como barragens. Atualmente, cada vez mais se tem reconhecida a importância da inclusão das incertezas na geração de previsões hidrológicas, feita através de previsões por conjunto (ou ensemble). Neste tipo de previsão são feitas inferências sobre cenários possíveis futuros através da consideração de, por exemplo, múltiplas trajetórias possíveis dos estados da atmosfera, que ao serem aplicadas em um modelo hidrológico resultam em distribuições de trajetórias de vazões. Várias aplicações recentes tem sugerido a possibilidade da tomada de melhores decisões para o futuro quando fundamentadas neste conhecimento das incertezas. No Brasil, um uso predominante de previsões hidrológicas é na operação de reservatórios de usinas hidroelétricas, que constituem a maior fonte de energia do País. As previsões nestes casos são utilizadas tanto para a operação normal do sistema nacional, feita de forma centralizada, como para a operação local das usinas em casos de cheia, onde é necessário velar pela segurança da barragem e pela atenuação de impactos a jusante e/ou a montante dos barramentos. Contudo, a forma como as previsões de vazão são geradas e usadas no cenário nacional não são baseadas em técnicas de previsão por conjunto, onde a própria pesquisa local sobre os potenciais benefícios destas formas de geração de previsões pode ser classificada como incipiente. Assim, o objetivo principal deste estudo é investigar benefícios em termos de qualidade e persistência do uso de previsões de afluência por conjunto para reservatórios em grandes bacias hidrográficas brasileiras em curto e médio prazo. Para cumprir com estes objetivos foram propostos ensaios de previsão de vazão por conjunto para três bacias hidrográficas brasileiras: Alto São Francisco, Doce, e Tocantins. O modelo hidrológico MGBIPH foi aplicado para a execução de previsões retroativas (hindcastings) alimentado por dados de chuva provindos de três diferentes sistemas de previsão meteorológica por conjunto (ECMWF-pf, GEFS, e CPTEC-pf) e mais uma previsão determinística de referência (ECMWF-fc), todos disponíveis na base de dados denominada TIGGE. De uma forma geral, as previsões por conjunto, principalmente dos modelos ECMWF-pf e GEFS, se mostraram superiores em termos de qualidade e persistência na comparação com a previsão determinística. E o uso do Super Ensemble, formado pela combinação dos três modelos mostrou-se uma alternativa entre as melhores testadas, principalmente por ser também uma estratégia robusta. Para uma estratégia de defesa contra cheias, as análises indicam benefícios para a consideração dos membros superiores dos conjuntos, e já para uma estratégia de operação de reservatórios essa visão pode ser mais focada em vazões médias, as quais podem conter algum viés. Já a comparação entre as bacias mostrou que resultados não podem ser transportados de um local para outro, apesar de estarem no mesmo clima. Em relação às incertezas, notou-se que a modelagem hidrológica amplifica as incertezas na previsão na medida em que os estados do modelo da grande bacia evoluem. De qualquer forma, acreditase que os resultados mostram que mais investimentos em técnicas de previsão por conjunto e suas aplicações são um caminho a ser seguido para ampliar os benefícios do uso de previsões hidrológicas. Short to medium-term streamflow forecasts at different locations in a watershed are generally beneficial to allow an early response to hydrological events such as floods, and more efficient operation of hydraulic structures such as dams. Currently, an increasingly recognition has been given to the including of uncertainties in the generation of hydrological forecasts, what is usually made producting the so called Ensemble Forecast. In this kind of forecast inferences about possible future scenarios are made by considering, for example, multiple possible trajectories of the atmospheric states, which when applied to a hydrological model results in streamflow trajectories distributions. Several recent applications suggested the possibility of better decisions making based on this uncertainties knowledge. In Brazil, a predominant use of hydrological forecasts is for hydropower reservoirs operation, which are the largest source of energy for the country. Future inflows estimates in these cases are used either for normal operation of the national system, done centrally, as for local operation of the dams in cases of floods, where it is necessary to ensure the dam safety and the mitigation of impacts downstream and/or upstream. However, the currently technique used to generate the operational forecasts is not based on ensembles, and the Brazilian local research on the potential benefits of these forms of forecasts production can be classified as incipient. Thus, the aim of this Thesis was to investigate benefits in terms of quality and persistence of using short to medium-term ensemble inflow forecast for reservoirs located on large Brazilian river basins. To fulfill these objectives streamflow forecast tests have been proposed for three Brazilian river basins: Alto San Francisco, Doce, and Tocantins. The hydrological model MGB-IPH was applied to perform retroactive forecasts (hindcastings) within a period of tests forced by rainfall data from three different ensemble weather forecasting systems (ECMWFpf, GEFs, and CPTEC-pf) and a deterministic prediction reference (ECMWF-fc), all available in the TIGGE archive. In general, the ensemble predictions, especially from ECMWF-pf and GEFs models, were superior in quality and persistence in comparison to the deterministic reference. And the use of the Super ensemble composed by the combination of the three ensemble models was shown to be among the results, and also a robust strategy. For a flood protection strategy, the analyzes indicate benefits in the consideration of the upper bounds of the ensembles, and for a reservoir operation strategy that vision could be more focused on average flow rates, which may contain some verified bias. The comparison between the basins results showed that one can not transport results and considerations from one location to another, despite being in the same climate region. Regarding uncertainties, it was noted that hydrological modeling amplifies the uncertainty in the forecasts, in some extent due to the large basin evolution of state variables. Anyway, it is believed that more investments in ensemble forecasting techniques and its applications shown to be a good way to make better use of forecasts.

Os efeitos das alterações de uso da terra na hidrologia das bacias hidrográficas têm sido, ao longo dos anos, amplamente discutidos para pequenas bacias experimentais. O uso de modelos hidrológicos para descrever esse tipo de processo tem se tornado promissor no que tange o estudo de grandes bacias. Sendo assim, o objetivo deste trabalho é avaliar uma metodologia para analisar o impacto das mudanças de uso da terra sobre o regime hidrológico em bacias hidrográficas de grande escala. Este estudo foi realizado por meio de simulação. Para isso foi utilizado um modelo hidrológico distribuído de base física, o Modelo de Grandes Bacias (MGB-IPH). Para realizar simulações de impactos das mudanças de uso da terra, inicialmente avaliou-se a capacidade do modelo MGB-IPH de captar as alterações produzidas pelas mudanças de uso. A análise de sensibilidade mostrou que a vazão gerada na bacia é inversamente proporcional aos parâmetros de índice de área foliar (IAF) e de altura média da vegetação (h), uma vez que quanto maior o IAF, mais interceptação ocorrerá na bacia e, quanto maior o valor de h, mais facilidade haverá nas trocas aerodinâmicas, ambos favorecendo o processo de evaporação. Já em relação aos parâmetros albedo (a) e resistência superficial (rs), a vazão gerada é diretamente proporcional, pois quanto maior o a maior a energia refletida pela superfície e, consequentemente menos energia estará disponível para o processo de evaporação. Em relação à rs, o seu aumento implica em uma menor perda por transpiração, disponibilizando mais água para gerar escoamento. As simulações de análise de sensibilidade mostraram que o MGB-IPH reproduz, de forma coerente, o comportamento médio das bacias submetidas ao desmatamento em um grande número de estudos experimentais. Assim, avaliou-se os impactos das mudanças de uso da terra no regime hidrológico em bacias de grande escala através de um estudo de caso. Este estudo de caso foi realizado na Bacia Hidrográfica do Rio Paraná. Os resultados das simulações mostraram que o desmatamento possui forte influência sobre as variáveis hidrológicas, como interceptação, evapotranspiração e vazão. Devido ao desmatamento foram verificadas reduções na interceptação e na evapotranspiração ao longo dos anos. Em relação a vazão, verificou-se que as simulações foram capaz de captar a não estacionariedade existente nas vazões observadas. Além disso, foi possível verificar que a vegetação pode ser associada a 39% dos incrementos observados na vazão média de longo termo, quando comparados os períodos anterior e posterior a 1970 na Bacia Hidrográfica do Rio Paraná. Já as variações na precipitação são responsáveis por 61% dos impactos na vazão média. Em relação às vazões extremas foram verificados incrementos, sendo eles mais expressivos nas vazões mínimas do que nas vazões máxima. Assim, pode-se afirmar que a metodologia proposta é consistente com os efeitos observados na bacia e com os resultados de estudos experimentais realizados em diferentes regiões do mundo. The effects of land use change on the hydrology of watersheds have been widely discussed during several years in small basins. The hydrological models have been promising for describe the effects of land use change in hydrology in large basins. Thus, the aim of this study was to evaluate a methodology to study the impact of land use changes on the hydrological processes in large basins. This study was conducted by hydrological simulation. For this, a physically based distributed hydrological model, Modelo de Grandes Bacias (MGB-IPH), was used. The ability of MGB-IPH model for study the effects of land-use change in hydrology processes was first analyzed. The results of sensibility analysis showed that the flow generated in the basin is inversely proportional to the parameters of leaf area index (LAI) and average vegetation height (h), since the higher the LAI more interception, and the higher the h easier will be the aerodynamic changes, both favoring the evaporation process . The flow is directly proportional to the parameters albedo (a) and surface resistance (rs), since the greater the a, more energy reflected by the surface, and consequently less energy is available for the evaporation process. The increasing in rs parameter implies a smaller loss by transpiration, and then more runoff is generated. The sensibility analysis showed that the MGB-IPH reproduces consistently, the average behavior of the basins subjected to deforestation on a large number of experimental studies. Thus, we evaluated the impacts of land use changes on the hydrological regime in large-scale basins by a case study. This case study was conducted in the Paraná River Basin. The simulation results showed that deforestation has a strong influence by hydrological variables, such as interception, evapotranspiration and flow. Due to deforestation were observed interception and evapotranspiration reduction over the years. It was found that the simulations were able to capture the non-stationarity in the observed streamflows. Furthermore, we found that 39% of the observed increase in average flow of long term is due to vegetation change, when comparing the periods before and after 1970 in the Paraná River Basin, while the 61% is due precipitation changes. Increments were observed in extreme flows. These increments were more evident at low flows than the maximum flow. Thus, it can be stated that the proposed methodology is consistent with the effects observed in the Rio Paraná basin and with the results of experimental studies conducted in different regions of the world.

Os reservatórios se constituem numa das principais formas estabelecidas pelo homem de modificar o ciclo hidrológico ao armazenar água para sua utilização no futuro. Em algumas regiões do país este mecanismo é de fundamental importância para o atendimento das demandas consultivas, seja pelas condições adversas do clima, como é o caso do semiárido nordestino, como também em função da alta demanda de água requerida para o atendimento da irrigação, como ocorre em algumas áreas do Rio Grande do Sul. Este trabalho apresenta uma metodologia de balanço hídrico em rios e reservatórios de pequeno porte, cujo foco é o atendimento aos instrumentos da Política Nacional de Recursos Hídricos, como a outorga de direito de recursos hídricos e as políticas de planejamento. A metodologia está baseada na integração entre o modelo hidrológico chuva-vazão MGB-IPH como recurso para obtenção das séries afluentes de cada reservatório e trecho de rio do sistema hídrico. Associado ao modelo hidrológico, está a integração com os Sistemas de Informações Geográficas, que possuem a vantagem de obtenção das características físicas da bacia de forma automática, facilitando também a associação com os objetos que representam o sistema hídrico (i.e. demandas e reservatórios). O modelo de simulação realiza a contabilização do movimento de água através de um sistema de reservatórios e trechos de rio, cuja progressão de cálculo é realizada de montante para jusante. Como técnica de otimização, foi utilizado o algoritmo SCE-UA, baseado na teoria dos algoritmos evolucionários. A metodologia desenvolvida foi aplicada na Bacia do Rio Quaraí, caracterizada pelo intenso uso da água para o atendimento das demandas da irrigação, e também pelo estabelecimento de centenas de pequenas estruturas de reservação de água, destinadas ao abastecimento destas demandas. Os resultados indicaram que não há água suficiente para o atendimento das demandas da irrigação em 100% do tempo. Em média, apenas 75% do volume total de água requerido para o atendimento das demandas atuais poderia ser atendido, considerando o período simulado de 20 anos. Nos anos mais críticos pode ocorrer o esgotamento do curso principal do Quaraí e outros locais, decorrência das retiradas de água, além do efeito da reservação de água pelos açudes, pois admitiu-se a hipótese de não existirem descarregadores de fundo nessas estruturas. Como alternativas para a gestão da água na bacia, apontam-se a definição de critérios de outorga para retiradas e liberação de água em reservatórios, além de implementação de novas estruturas de reservação e regularização. Reservoirs are the most important form established by the man to modify the hydrological cycle to store water for use in the future. In some regions of Brazil this mechanism is crucial to meet the advisory needs, either by adverse weather conditions, such as the semi-arid northeast, but also due to the high demand of water required to meet the irrigation, as occurs in some areas of Rio Grande do Sul. This work presents a methodology of water balance in rivers and small reservoirs, whose focus is to achieve the goals of the instruments of the National Water Resources Policy, as the water rights and the planning policies. The methodology is based on the integration between the hydrological rainfall-runoff model MGB-IPH as a resource for obtaining streamflow series from upstreams of each reservoir and for all sub-watershed defined on the water system. Associated with the hydrological model is the integration with Geographic Information Systems, which have the advantage of obtaining the physical characteristics of the basin automatically, it also facilitates association with objects that represent the water system (i.e. demands and reservoirs). The simulation model performs accounting for the movement of water through a system of reservoirs and river stretches whose progression calculation is performed from upstream to downstream. As optimization technique, the SCE-UA algorithm, based on the theory of evolutionary algorithms was used. The procedure was applied in the Quaraí River Basin, characterized by intense use of water to meet the irrigation demands, and also by the establishment of hundreds of small structures reservation of water, destined to supply these demands. The results indicated that there is not water enough to meet the demands of irrigation at 100% of the time. On average, only 75% of the total volume of water required to meet the current demands could be met, considering the simulated period of 20 years. In the most critical years can occur exhaustion of the main course of the Quaraí and other sites, due to withdrawals plus the effect of the reservation by dams, it admitted the hypothesis that there are no releases in this structures. As alternatives to water management in the basin, point the definition of licensing criteria for withdrawals and releases in dams, as well as implementing new reservoir with regularization structures.

A modelagem matemática é uma valiosa ferramenta para avaliação de alterações da qualidade da água causadas por descargas diversas. Para a simulação de bacias hidrográficas a integração com sistemas de informações geográficas (SIG) é muito útil para a descrição detalhada e automatizada da região de interesse, e para a interpretação dos resultados. Neste contexto de desenvolvimento de soluções integradas que o presente trabalho tem como objetivo contribuir para a melhoria no gerenciamento da qualidade da água em bacias hidrográficas, onde é apresentado o desenvolvimento de um modelo acoplado com SIG, aplicável para condições de baixa disponibilidade de dados, denominado SIAQUA-IPH. Adotou-se para o modelo a utilização de soluções analíticas da equação de advecção-dispersão para simular o impacto de lançamentos com duração variável em cursos de água, incluindo lançamentos instantâneos, ou de duração curta. No modelo é adotado um esquema de convolução dos polutogramas para a representação de lançamentos múltiplos e de confluências existentes na bacia. Ainda, para a simulação de constituintes comumente encontrados em esgotos (DBO, OD, N, P e Coliformes) foi adotada a utilização de soluções conhecidas, como a equação de Streeter-Phelps. Soluções analíticas foram priorizadas visando diminuir restrições relacionadas à capacidade computacional, tempo de simulação e outros erros numéricos. A topologia adotada na concepção do modelo consiste na subdivisão da bacia hidrográfica em pequenos trechos de rio, que pode ser obtida automaticamente através de ferramentas de geoprocessamento, o que caracteriza o primeiro fator de integração com SIG desenvolvida. Esta é a mesma topologia utilizada pelo modelo MGB-IPH, por isto os resultados do modelo hidrológico podem ser utilizados diretamente pelo modelo de qualidade. Como segundo fator de integração com SIG, o SIAQUA-IPH foi desenvolvido plenamente acoplado à plataforma MapWindow GIS®, de uso livre, onde o SIG é a forma como o usuário interage com o modelo, inserido os dados de entrada e visualizando os resultados. Para a calibração do SIAQUA-IPH adotou-se a possibilidade do uso do algoritmo de calibração automática SIMPLEX de Nelder-Mead. Para o cálculo do coeficiente de dispersão longitudinal foi adotado o uso de formulações práticas, selecionadas a partir de uma matriz de seleção desenvolvida. Testes hipotéticos de funcionamento do modelo nas bacias do rio Ijuí (RS) e do rio São Francisco (MG) evidenciaram um desempenho pleno de todas as soluções adotadas. Em um terceiro ensaio de representação de dados de traçadores, obtidos de Rigo (1992), foram obtidos resultados satisfatórios em que as diferenças computadas servem de referência para usuários sobre as incertezas dos métodos utilizados. Espera-se que as soluções adotadas no desenvolvimento do SIAQUA-IPH sirvam para várias aplicações, incluindo estudos prévios de impacto, avaliações posteriores de despejos ocorridos, e até suporte à decisão em tempo real, como no caso de derramamentos acidentais. Mathematical modeling is a valuable tool to assess water quality changes produced by different kinds of discharges. For watershed simulations, the coupling with Geographic Information Systems (GIS) is very useful for automated and detailed description of the studied region and for helping results interpretation. On this context of integrated solutions development that this research aims to contribute for river basins water quality management improving, where is exposed the development of a model coupled with GIS and applicable to low data availability conditions, called SIAQUA-IPH. In the model was adopted the use of advection-dispersion equation analytical solutions to simulate the impacts of various duration discharges into rivers. A polutograph convolution scheme is adopted for representing multiple spills and confluences in the basin. For constituents commonly found in sewage simulation (BOD, OD, N, P and Coliforms) the model also allows the use of classical solutions, such as the Streeter-Phelps equation. The model is designed to use primarily analytical solutions aiming to reduce constraints related to computational capacity, simulation time and others numerical issues. The topology adopted in the model consists in the hydrography discretization into small river reaches, which can be obtained automatically through geoprocessing tools, which is the primary factor of developed GIS coupling. Also, as this is the same topology used by the MGB-IPH model, the hydrological model results can be used directly by the water quality model. As a second GIS integration factor, SIAQUA-IPH was developed fully coupled to MapWindow GIS ® platform, free distributed, where the GIS is the way the user interacts with the model, enters the input information and access results. Calibration of SIAQUA-IPH can be performed through Nelder-Mead SIMPLEX automatic calibration algorithm. For the longitudinal dispersion coefficients estimation was adopted the use of empirical formulations, selected through a selection matrix developed. Tests of hypothetical releases in the Ijuí (RS) and in the São Francisco (MG) river basins were carried out, which showed a full performance of all the adopted solutions. Also, a tracer data representation test was conducted using data showed by Rigo (1992) in the Paraíba do Sul River (RJ), where obtained results were considered satisfactory and computed differences can be used as a guide for uncertainty of the adopted methods. It is expected that the SIAQUA-IPH serves for various applications, including prior studies of pollutants releases, later evaluations of discharges and to support real-time decision in case of accidental spills.

Le bassin Amazonien est connu comme le plus grand système hydrologique du monde et pour son rôle important sur le système terre, influençant le cycle du carbone et le climat global. Les pressions anthropiques récentes, telles que la déforestation, les changements climatiques, la construction de barrage hydro-électriques, ainsi que l’augmentation des crues et sécheresse extrêmes qui se produisent dans cette région, motivent l’étude de l’hydrologie du bassin Amazonien. Dans le même temps, des méthodes hydrologiques de modélisation et de surveillance par observation satellitaire ont été développées qui peuvent fournir les bases techniques à cette fin. Ce travail a eu pour objectif la compréhension et la prévision du régime hydrologique du bassin Amazonien. Nous avons développé et évaluer des techniques de modélisation hydrologique-hydrodynamique de grande échelle, d’assimilation de données in situ et spatiales et de prévision hydrologique. L’ensemble de ces techniques nous a permis d’explorer le fonctionnement du bassin Amazonien en terme de processus physiques et de prévisibilité hydrologique. Nous avons utilisé le modèle hydrologique-hydrodynamique de grande échelle MGB-IPH pour simuler le bassin, le forçage précipitation étant fourni par l’observation spatiale. Les résultats de la modélisation sont satisfaisants lorsque validés à partir de données in situ de débit et de hauteurs d’eau mais également de données dérivées de l’observation spatiale incluant les niveaux d’eau déduits de l’altimétrie radar, le contenu en eau total issu de la gravimétrie satellitaire, l’extension des zones inondées. Nous avons montré que les eaux superficielles sont responsables en grande partie de la variation du stock total d’eau, l’influence des grands plans d’eau sur la variabilité spatiale des précipitations et l’influence des plaines d’inondation et des effets de remous sur la propagation des ondes de crues. Nos analyses ont montré le rôle prépondérant des conditions initiales, en particulier des eaux superficielles, pour la prévisibilité des grands fleuves Amazoniens, la connaissance des précipitations futures n’ayant qu’une influence secondaire. Ainsi, pour améliorer l’estimation des variables d’état hydrologiques, nous avons développé, pour la première fois, un schéma d’assimilation de données pour un modèle hydrologique-hydrodynamique de grande échelle, pour l’assimilation de données de jaugeages in situ et dérivées de l’altimétrie radar (débit et hauteur d’eau), dont les résultats se sont montrés satisfaisants. Nous avons également développé un prototype de système de prévision des débits pour le bassin Amazonien, basé sur le modèle initialisé avec les conditions initiales optimales fournies par le schéma d’assimilation de données, et en utilisant la pluie estimée par satellite disponible en temps réel. Les résultats ont été prometteurs, le modèle étant capable de prévoir les débits dans les principaux fleuves Amazoniens avec une antécédence importante (entre 1 et 3 mois), permettant d’anticiper, par exemple, la sècheresse extrême de 2005. Ces résultats démontrent le potentiel de la modélisation hydrologique appuyé par l’observation spatiale pour la prévision des débits avec une grande antécédence dans les grands bassins versant mondiaux. A bacia Amazônica se destaca como o principal sistema hidrológico do mundo e pelo seu importante papel no sistema terrestre, influenciando o ciclo de carbono e o clima global. Recentes pressões antrópicas, como o desflorestamento, mudanças climáticas e a construção de barragens hidroelétricas, somados às crescentes cheias e secas extremas ocorridas nesta região, motivam o estudo da hidrologia da bacia Amazônica. Ao mesmo tempo, têm se desenvolvido métodos hidrológicos de modelagem e monitoramento via sensoriamento remoto que podem fornecer as bases técnicas para este fim. Este trabalho objetivou a compreensão e previsão da hidrologia da bacia Amazônica. Foram desenvolvidas e avaliadas diversas técnicas, incluindo de modelagem hidrológica-hidrodinâmica de larga escala, de assimilação de dados in situ e de sensoriamento remoto, e de previsão hidrológica. Este conjunto de técnicas foi utilizado para compreender o funcionamento da bacia Amazônica em termos de seus processos hidrológicos e sua previsibilidade hidrológica. O modelo hidrológico-hidrodinâmico de larga escala MGB-IPH foi utilizado para simular a bacia, sendo forçado com dados de chuva estimados por satélite. O modelo mostrou bom desempenho em uma validação detalhada contra observações de vazões e cotas in situ além de dados oriundos de sensoriamento remoto, incluindo níveis d’água de altimetria por radar, armazenamento d’água de gravimetria espacial e extensão de áreas alagadas. Mostrou-se a dominância das águas superficiais nas variações do armazenamento de água, a influência dos grandes corpos d’água sobre a variabilidade espacial da precipitação, além da importância das várzeas da inundação e efeitos de remanso sobre a propagação das ondas de cheia Amazônicas. As condições hidrológicas iniciais, com destaque para as águas superficiais, mostraram dominar a previsibilidade hidrológica nos grandes rios amazônicos, tendo assim a precipitação no futuro um papel secundário. Portanto, afim de melhor estimar os estados hidrológicos, de forma pioneira, foi desenvolvido um esquema de assimilação de dados para um modelo hidrológicohidrodinâmico de larga escala para assimilar informações in situ e de altimetria por radar, cujo desempenho se mostrou satisfatório. Desenvolveu-se também um protótipo de sistema de previsão de vazões para a bacia Amazônica, baseado no modelo inicializado com condições iniciais ótimas do esquema de assimilação de dados e utilizando precipitação estimada por satélite disponível em tempo real. Os resultados foram promissores e o modelo foi capaz de prever vazões nos principais rios amazônicos com grande antecedência (~1 a 3 meses), antecipando, por exemplo, a grande seca de 2005. Estes resultados mostram o potencial da modelagem hidrológica de larga escala apoiada por informação de sensoriamento remoto na previsão de vazões com alta antecedência nas grandes bacias hidrográficas do mundo. The Amazon basin is known as the world’s main hydrological system and by its important role in the earth system, carbon cycle and global climate. Recent anthropogenic pressure, such as deforestation, climate change and the construction of hydropower dams, together with increasing extreme floods and droughts, encourage the research on the hydrology of the Amazon basin. On the other hand, hydrological methods for modeling and remotely sensed observation are being developed, and can be used for this goal. This work aimed at understanding and forecasting the hydrology of the Amazon River basin. We developed and evaluated techniques for large scale hydrologic-hydrodynamic modeling, data assimilation of both in situ and remote sensing data and hydrological forecasting. By means of these techniques, we explored the functioning of the Amazon River basin, in terms of its physical processes and its hydrological predictability. We used the MGB-IPH large scale hydrologichydrodynamic model forced by satellite-based precipitation. The model had a good performance when extensively validated against in situ discharge and stage measurements and also remotely sensed data, including radar altimetry-based water levels, gravimetric-based terrestrial water storage and flood inundation extent. We showed that surface waters governs most of the terrestrial water storage changes, the influence of large water bodies on precipitation spatial variability and the importance of the floodplains and backwater effects on the routing of the Amazon floodwaves. Analyses showed the dominant role of hydrological initial conditions, mainly surface waters, on hydrological predictability on the main Amazon Rivers, while the knowledge of future precipitation may be secondary. Aiming at the optimal estimation of these hydrological states, we developed, for the first time, a data assimilation scheme for both gauged and satellite altimetry-based discharge and water levels into a large scale hydrologic-hydrodynamic model, and it showed a good performance. We also developed a forecast system prototype, where the model is based on initial conditions gathered by the data assimilation scheme and forced by satellite-based precipitation. Results are promising and the model was able to provide accurate discharge forecasts in the main Amazon rivers even for very large lead times (~1 to 3 months), predicting, for example, the historical 2005 drought. These results point to the potential of large scale hydrological models supported with remote sensing information for providing hydrological forecasts well in advance at world’s large rivers and poorly monitored regions.

A previsão e emissão de alertas antecipados constituem um dos principais elementos na prevenção dos impactos ocasionados por eventos de cheias. Uma das formas utilizadas para se obter uma ampliação do horizonte de previsão é através do uso da modelagem chuva-vazão associada à previsão de precipitação, tipicamente derivada de modelos meteorológicos. A precipitação, no entanto, é uma das variáveis que impõe maior dificuldade na previsão meteorológica, sendo considerada uma das principais fontes de incerteza nos resultados da previsão de cheias. A previsão por conjunto é uma técnica originalmente desenvolvida nas ciências atmosféricas e procura explorar as incertezas associadas às condições iniciais e/ou deficiências na estrutura dos modelos meteorológicos com intuito de melhorar sua previsibilidade. A partir de diferentes modelos meteorológicos ou de diferentes condições iniciais de um único modelo, são gerados um conjunto de previsões que representam possíveis trajetórias dos processos atmosféricos ao longo do horizonte de previsão. Pesquisas recentes, principalmente na Europa e Estados Unidos, têm mostrado resultados promissores do acoplamento de previsões meteorológicas por conjunto à modelos hidrológicos para realizar previsões de cheia. Essa pesquisa trata da avaliação do benefício da previsão de cheias por conjunto em curto prazo, em uma bacia de médio porte, utilizando dados e de ferramentas para previsão de vazões disponíveis em modo operacional no Brasil. Como estudo de caso foi utilizada a bacia do Rio Paraopeba (12.150km²), de clima tipicamente tropical, localizada na região sudeste do Brasil. A metodologia proposta para geração das previsões hidrológicas utilizou o modelo hidrológico MGB-IPH alimentado por um conjunto previsões de precipitação de diferentes modelos, com diferentes condições iniciais e parametrizações, dando origem a distintos cenários de previsão de vazões. Como parâmetro de referência na avaliação do desempenho das previsões por conjunto foi utilizada uma previsão hidrológica determinística única, baseada em uma previsão de precipitação obtida da combinação ótima de saídas de diversos modelos meteorológicos. As previsões foram realizadas retrospectivamente no período entre ago/2008 e mai/2011, sendo analisadas durante o período chuvoso dos anos hidrológicos (out-abr). Os resultados das previsões de cheia por conjunto foram avaliados através de uma representação determinística, considerando a média dos membros do conjunto, assim como através de uma representação probabilística, considerando todos os membros, através de medidas de desempenho específicas para esse fim. Na avaliação determinística, a média do conjunto hidrológico apresentou resultados similares aos obtido com a previsão determinística de referência, embora tenha apresentado benefício significativo em relação à maior parte dos membros do conjunto. A avaliação das previsões de cheia por conjunto, por sua vez, mostrou a existência de uma superestimativa e de um subespalhamento dos membros em relação às observações, sobretudo nos primeiros intervalos de tempo da previsão. Na comparação dos resultados das previsões de eventos do tipo dicótomos, que consideram a superação ou não de vazões limites de alerta, o 9º decil das previsões por conjunto mostrou superioridade em relação à previsão determinística de referência e mesmo a média do conjunto, sendo possível obter, na maior parte dos casos analisados, um aumento significativo na proporção de eventos corretamente previstos mantendo as taxas de alarmes falsos em níveis reduzidos. Esse benefício foi, de modo geral, maior para maiores antecedências e vazões limites, situações mais importantes num contexto de prevenção de cheias. Os resultados mostraram ainda que, em média, uma diminuição do número de membros do conjunto diminui seu desempenho nas previsões. The forecasting and issuing of early warnings represent a key element to prevent the impacts of flood events. An alternative to extend forecasting horizon is the use of rainfall-runoff modeling coupled with precipitation forecasts derived from numerical weather prediction (NWP) models. However, NWP models have difficulty to accurately predict precipitation due to the extremely sensitivity of the initial conditions. Therefore, this variable represents one of the major sources of uncertainties in flood forecasting. A probabilistic or ensemble forecasting approach was originally developed in the atmospheric sciences and then applied to other research areas. This procedure explores the uncertainties related to initial conditions and deficiencies in the structure of NWP models intending to improve its predictability. Using different NWP models or different initial conditions of a single model, an ensemble forecast showing possible trajectories of atmospheric processes over the forecast horizon are produced. Recent studies developed in Europe and the United States have shown promising results in flood forecasting using hydrological models fed by NWP ensemble outputs. The present research assess the performance of short term ensemble flood forecasting in a medium size tropical basin, based on data and streamflow forecasting tools available in operational mode in Brazil. The Paraopeba River basin (12,150 km²), located in the upper portion of the São Francisco River basin, in Southeastern Brazil, was selected as a case study. The proposed methodology used the MGB-IPH hydrological coupled to an ensemble of precipitation forecasts generated by several models with different initial conditions and parameterizations. The results are several scenarios of streamflow forecasts. A single deterministic streamflow forecast, based on a quantitative precipitation forecast derived from the optimal combination of several outputs of NWP models, was used as a reference to assess the performance of the streamflow ensemble forecasts. The streamflow forecasts were performed between aug/2008 and may/2011 and were analyzed during the rainy seasons (austral summer). The results from the ensemble flood forecasting were assessed by deterministic and probabilistic performance measures, with the ensemble mean being used by the former, and specific assessment measure by the later. Based on the deterministic assessment, the ensemble mean showed similar results to those obtained by the deterministic reference forecast, although showing better performance over most of the ensemble members. Based on the probabilistic performance measures, however, results showed the existence of an ensemble overforecasting and underspread of the members in regard to observed values, especially during the first lead times. The results for predictions of dichotomous events, which mean exceeding or not flood warning thresholds, showed that the 9th decile of the ensemble over performed the deterministic forecast and even the ensemble mean. In most cases, it was observed an increase in the proportion of correctly forecasted events while keeping false alarm rates at low levels. This benefit was generally higher for higher flow thresholds and for longer lead times, which are the most important situations for flood mitigation. The results show, also, that, in average, a reduction in the number of ensemble members decreases the performance of ensemble flood forecasts.

Conselho Nacional de Desenvolvimento Científico e Tecnológico ABSTRACT: In regions of natural or planted forests, rainfall interception by vegetation can be an important variable in the hydrological cycle and should not be neglected. This research analyzes the process of interception in different forest types (native and planted) in the region of Santa Maria - RS, at the head of Vacacaí Mirim. The variables Rainfall (R), stemflow (Sf), and throughfall (Tf) were monitored in the period from Dec/09 to Fev/11. Rainfall was monitored in two tipping bucket rain gauges and meteorological station installed in the subbasin on July 14th, 2010. Throughfall was monitored in three plots (areas I, II and III) with dimensions of 12 m x 50 m, 21 m x 12 m and 10 m x 20 m, respectively. The sampling points were installed 40 or 60 cm from the ground and spaced 1 m in straight line. In area I (native forest), 101 stakes were allocated to 40 instruments; where 20 collectors remained in fixed locations and 20 were rotated randomly. In the areas II (native forest) and III (eucalyptus) 60 sampling points were allocated to 20 collectors; at each new collection they were placed one sampling point ahead the former. Stemflow was based on values suggested in literature and by monitoring the flow in three selected trees within each plot, and the interception was calculated as the difference between precipitation and the sum of throughfall and stemflow. The results showed that interception is a key component of the hydrological cycle and, in general, underestimated by the hydrological models.The highest interception value was checked for native forest and the largest stemflow for eucalyptus. The stemflow was higher than the values observed by other authors, revealing that this variable cannot be ruled out in the study of interception. On the other hand, the ranges of variation of throughfall and interception were close to those reported in other studies. The results seem also to suggest that the depth parameter adopted for the calculation of interception using the leaf area index is underestimated and should be 0.4. Statistical analysis of the number of collectors used and necessary to ensure error of 5 mm of throughfall in relation to its average and 90% of confidence, showed that for vegetation such as Eucalyptus, the use of 20 rotating collectors were sufficient in 95.45% of the collections; which may be a reference to works in which there is a pattern of similarity of species. For native forest we observed that when the vegetation has a low density and the similar pattern of canopy, stem diameter and height, the use of 20 rotating collectors would be enough to ensure statistical requirements set out in 71.43% of the samples monitored. When the native vegetation is dense, with standard range of height, stem diameter and crown area, we observed the need for a higher number of collectors equipment of throughfall. The use of 40 collectors (20 fixed and 20 rotational) assured that the statistical requirements stipulated were guaranteed in 60.98% of the samples monitored. RESUMO: Nas regiões de floretas naturais ou plantadas, a interceptação da precipitação pela vegetação pode ser uma importante variável do ciclo hidrológico, não devendo ser negligenciada. Esta pesquisa analisou o processo de interceptação em diferentes formações florestais (nativa e plantada) na região de Santa Maria - RS, na cabeceira do Vacacaí Mirim. As variáveis monitoradas no período de dez/09 a fev/11 foram: precipitação inicdente (P), escoamento pelo tronco (Esct) e precipitação interna (Pi). A precipitação incidente foi monitorada em dois pluviômetros tipo cubas-basculantes e pela estação climatológica instalada na sub-bacia em 14/7/2010. A precipitação interna foi monitorada em três parcelas (área I, II e III) com dimensões de 12 m x 50 m, 21 m x 12 m e 10 m x 20 m, respectivamente. As estacas foram instaladas a 40 ou 60 cm do solo, sendo espaçadas de 1 m e em linha reta. Na área I (mata nativa), foram alocadas 101 estacas para disposição de 40 instrumentos, sendo que 20 permaneceram fixos e 20 foram rotacionados. Já nas áreas II (mata nativa) e III (eucalipto) foram alocadas 60 estacas para disposição de 20 coletores, sendo que a cada nova coleta eles eram realocados uma estaca a frente da anterior. O escoamento pelo tronco foi determinado com base em valores sugeridos em literatura e por meio de monitoramento do escoamento em três árvores escolhidas dentro de cada parcela e a interceptação foi calculada pela diferença entre P e a soma de Pi e Esct. Os resultados mostraram que a interceptação é uma componente fundamental do ciclo hidrológico e, em geral, subestimada pelos modelos hidrológicos. O maior valor de interceptação foi verificado para mata nativa e o maior Esct para eucalipto. O escoamento pelo tronco foi superior aos valores verificados por outros autores, revelando que essa variável não pode ser descartada no estudo da interceptação. Já as faixas de variação de precipitação interna e da interceptação ficaram próximas aquelas registradas em outros estudos. Os resultados parecem sugerir, ainda, que o parâmetro de lâmina adotado para o cálculo de interceptação com utilização do índice de área foliar está subestimado e deveria ser 0,4. A análise estatística do número de aparelhos utilizados e necessários para garantir 5 mm de erro de Pi em relação a sua média e 90% de confiança, mostrou que para vegetação do tipo eucalipto, o uso de 20 instrumentos rotacionáveis seriam suficientes em 95,45% das coletas realizadas. O que pode ser uma referência para trabalhos em que existe um padrão de semelhança de espécies. Já para mata nativa observou-se que, quando a vegetação apresenta baixa densidade e mesmo padrão de copa, diâmetro do caule e altura, o uso de 20 instrumentos rotacionáveis seria suficiente para garantir os requisitos estatísticos estipulados em 71,43% das coletas monitoradas. Quando a vegetação nativa é densa, com padrão variado de altura, diâmetro do caule e área da copa, percebeu-se a necessidade de um número mais elevado de equipamentos coletores de Pi. O uso de 40 instrumentos, (20 aparelhos rotacionáveis e 20 fixos), garantiu que em 60,98% das coletas monitoradas os requisitos estatísticos estipulados fossem assegurados.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior Areas under intensive tobacco crop cultivation have been impacting the water balance and increasing soil erosion in Southern Brazil. Application of large amounts of mineral fertilizers, higher than the tobacco requirement, contributes to excessive phosphorus loads in soil and water bodies. The assessment of potential environmental impacts due to changes in land use and management practices is necessary to achieve the sustainable management of natural resources such as soil and water. In this study the hydrological processes, the sediment yield and the phosphorus transfer in the Arroio Lino watershed were evaluated by using the Soil and Water Assessment Tool (SWAT) model. The small watershed (4.18 km2 ) is located in Southern Brazil and its main land use is tobacco crop under conventional tillage. Measured streamflow and sediment yield at the watershed outlet were used for model streamflow sensitivity analysis, calibration and validation. Phosphorus loads at four sub-watersheds (A1, A2, B, C) and at the watershed outlet were used for model sensitivity analysis and calibration. A Latin Hypercube (LH) and One-factor- At-a-Time (OAT) sensitivity analysis was performed on input variables. Model calibration was performed with the Shuffled Complex Evolution Algorithm-Uncertainty Analysis (SCE-UA). Time series plots and standard statistical measures were used to verify model predictions. Adaptations of the model parameters for the reality of the region in which the watershed is located were made during the calibration and validation of the model. The predicted monthly streamflow matched the observed values, with a Nash Sutcliffe coefficient (NSE) of 0.87 for calibration and 0.76 for validation. In an attempt to account for sediment transport and deposition processes across the landscape, the SWAT model was modified to simulate landscape sediment transport capacity. The new deposition routine performed better during calibration (NSE of 0.70) than SWAT standard version (NSE of -0.14) in the studied watershed. The modified model provided reasonable simulations of sediment transport across the landscape positions. The results suggest that the integration of the sediment deposition routine in SWAT increases model predictions accuracy in steeper areas, while at the same time significantly improves the ability to predict spatial distribution of sediment deposition areas. The predicted P loads are in the order of magnitude of the measured ones, however, the model failed to predict the P loads in three sub-watersheds (A1, A2 and B). Although occupying only 29% of the total land cover, cropland is the primary source of nutrients in the watershed (80%). After calibration of streamflow, sediment yield and phosphorus loads, the simulation of distinct management scenarios was done for the Arroio Lino watershed. Three scenarios of management practices: conventional tillage (CT), minimum tillage (MT) and no-tillage cultivation (NT) with reduction of 50% of fertilizer rate application were tested over a period of 30 years. No-tillage practices did not significantly affect water yield, but greatly affected sediment yield due to reduction of soil erosion. The soluble phosphorus losses decreased mainly when the fertilizer doses decreased. The simulation results suggest that conventional tillage practices should be replaced by less intensive tillage practices in order to minimize environmental impacts caused by a single land use. Áreas com culturas de tabaco cultivadas intensivamente, no sul do Brasil, vêm afetando a qualidade da água e aumentando a erosão do solo nesta região. A aplicação de grandes quantidades de fertilizantes minerais, maiores do que a exigência da cultura, contribui para a excessiva carga de fósforo nos solos e corpos d'água. A avaliação dos impactos potenciais das mudanças no uso e manejo do solo é fundamental para propiciar uma gestão sustentável dos recursos naturais como solo e água. Neste estudo foram avaliados os processos hidrológicos, a produção de sedimentos e a transferência de fósforo na bacia hidrográfica do Arroio Lino, por meio do modelo Soil and Water Assessment Tool (SWAT). A pequena bacia (4,18 km2 ) localiza-se no sul do Brasil e seu principal uso da terra é a cultura de fumo sob plantio convencional. As vazões e a produção de sedimentos medidos no exutório da bacia foram utilizadas para análise de sensibilidade, calibração e validação dos parâmetros do modelo. Cargas de fósforo em quatro sub-bacias (A1, A2, B, C) e no exutório foram utilizados para a análise de sensibilidade e calibração do modelo. A análise de sensibilidade foi feita com o uso de um algoritmo que combina as técnicas de Hipercubo Latino (LH) e One-factor-At-a-Time (OAT). A calibração foi realizada com o algoritmo Shuffled Complex Evolution (SCE-UA). Análises gráficas e medidas estatísticas foram utilizadas para verificar as previsões do modelo. Adaptações nos parâmetros do modelo foram feitas durante as etapas de calibração e validação, tendo em vista a realidade da região em que a bacia está localizada. Na simulação de vazões mensais obtiveram-se valores do índice de eficiência de Nash e Sutcliffe (NSE) de 0,87 na calibração e 0,76 na validação. Com o intuito de representar melhor os processos de deposição de sedimentos nas encostas, o modelo SWAT foi modificado com a incorporação de um termo de capacidade de transporte de sedimentos na paisagem. Com a nova rotina de sedimentos obteve-se um melhor desempenho durante a calibração (NSE de 0,70) em relação ao modelo SWAT original (NSE de -0,14). As simulações com o modelo modificado foram satisfatórias para o transporte de sedimentos em diferentes posições da paisagem. Os resultados sugerem que a integração da rotina de deposição de sedimentos no SWAT aumenta a precisão do modelo de previsões em áreas mais íngremes e, ao mesmo tempo, melhora significativamente a capacidade de prever a distribuição espacial das áreas de deposição de sedimentos. As previsões de cargas de fósforo foram na ordem de grandeza das cargas medidas; no entanto, o modelo não conseguiu prever satisfatoriamente as cargas de fósforo em três sub-bacias (A1, A2 e B). Apesar das lavouras ocuparem apenas 29% da área total da bacia, de acordo com a simulação elas são a principal fonte de nutrientes na bacia hidrográfica (80%). Após a calibração da vazão, da produção de sedimentos e de fósforo, foi realizada a simulação de cenários de manejo do solo na bacia hidrográfica do Arroio Lino. Três cenários de práticas de manejo: preparo convencional (PC), cultivo mínimo (CM) e plantio direto de cultivo (NT) com redução de 50% da taxa de aplicação de fertilizantes foram testados durante um período de 30 anos. A prática de plantio direto não afetou significativamente o escoamento, no entanto afetou fortemente a produção de sedimentos devido à redução da erosão do solo. Houve redução das perdas de fósforo principalmente devido à redução das doses de fertilizantes. Os resultados indicam que as práticas de preparo convencional do solo deveriam ser substituídas por práticas de cultivo mínimo ou direto, a fim de minimizar os impactos ambientais causados por um determinado uso do solo.

Neste trabalho avaliou-se o impacto de possíveis cenários de mudanças climáticas sobre os recursos hídricos na bacia do Rio Ibicuí (35,158.43 Km².) afluente da margem esquerda do rio Uruguai em território brasileiro. Os impactos de mudanças climáticas na bacia do Rio Ibicuí foram avaliados através do uso de um modelo hidrológico, analises de sensibilidade de vazão e cenários de mudanças climáticas. As analises de sensibilidade mostraram que as vazões mínimas são menos sensíveis a mudanças de precipitação do às médias e máximas. Cenários de mudanças climáticas para o ano de 2050 (janela de 30 anos centrada em 2050) foram obtidos através do gerador de cenários MAGICC/SCENGEN, que considera resultados de vinte modelos climáticos globais (GCM’s) e três cenários de emissão de gases de efeito estufa. A totalidade dos resultados indicaram aumento da temperatura, porém divergências quanto aos resultados de referentes a precipitação. Alguns dos GCM’s indicaram diminuição da precipitação média anual enquanto outros indicaram o aumento desta variável. Os resultados mensais de modificação na temperatura e precipitação média previstos para 2050 na região da bacia do Ibicuí, obtidos pelos vinte GCMs para cada um dos três cenários, foram utilizados para perturbar as séries observadas e simular o impacto sobre o regime hidrológico através do modelo MGB-IPH. As incertezas relacionadas às projeções de anomalias de precipitação e temperatura, previstas pelo conjunto de 20 GCM’s, foram amplificadas ao avaliar o impacto de mudanças climáticas no regime de vazões. The impact of possible climate change scenarios on water resources in the Ibicui River basin was assessment in this work. The Ibicui River (35,158.43 Km²) is tributary of Uruguay River and it is located on border of southern of Brazil. The impacts of climate change on the Ibicui River basin were evaluated using a large scale distribuited hydrological model, sensitivity analyses and climate change scenarios. Sensitivity analyses showed that minimum discharges in the Ibicui River are less sensitive to change in mean precipitation when compared to maximum and mean discharge. Climate change scenarios for 2050 (30-years time interval centered in 2050) was obtained through a climate scenario generator, MAGICC/SCENGEN, considering the results from twenty general circulation models (GCM’s) and three emission scenarios of greenhouse gas emissions. Results indicate an increase in the annual mean temperature, but divergences about precipitation. Some GCM’s indicated increase in the annual mean precipitation while others indicated decrease. The hydrological impact of all projections of climate change for 2050 was simulated using the MGB-IPH model, by perturbing the observed times series of temperature and precipitation according to relative changes projected by GMC’s. The uncertainties in projected climate changes are amplified in the hydrological results.