Four small watersheds, Kamwamphula, Luelo, Kanyanga and Mphiri, near Emusa (Lundazi District) in Eastern Province, Zambia, were studied to document transitions in land use over time and to project the impacts of land use and topography on runoff, erosion and sediment delivery. Landuse was delineated from 2007 IKONOS image (one meter pixel), and Landsat imagery was used to depict the historic changes in landuse between the period of 1989 and 2007. The GWLF model was used to predict the impact of the landuses on the hydrology of the area. There had been an increase in clearing of forest area mainly due to the expansion of cropland area. The highest rate of clearing was predicted for the Kamwamphula watershed where the forest cover decreased from 95% to 71% over the 18 year period. The GWLF model was used to predict the impact of the landuse on the hydrology and sediment delivery. In comparison with the limited field data available from the four watersheds, the GWLF model gave poor prediction of streamflow, probably because the hydrology of the area is poorly understood and dambo function in the landscape is not well represented in the model. Highest runoff, erosion and sediment yields came from Luelo watershed which has steeper slopes and, less vegetative cover and poor permeability of soils. The GWLF gave poor prediction of streamflow, probably, because the hydrology of the area is poorly understood. / Master of Science

In 2005 an assessment of existing Total Maximum Daily Load studies by the U.S. Environmental Protection Agency showed sediment as the fourth leading cause of water quality impairment. A source assessment is important in developing a successful TMDL. Past research efforts have focused on controlling erosion sources in agricultural and urban land areas. New research suggests major contributions to overall sediment loads may be due to stream channel degradation. Monitoring and modeling techniques to assess the contribution of channel sediment to overall sediment load are needed to determine the reductions necessary to meet water quality standards. This research focused on testing the ability of watershed and reach-scale models to predict stream channel degradation. Model predictions were compared to estimates developed from a system of erosion pins and scour chains. A 500-m experimental reach in Blacksburg, VA, USA, was selected as the focus of channel degradation monitoring and modeling efforts. A series of over 250 erosion pins and seven scour chains were installed systematically throughout the experimental reach. A monthly monitoring program measured channel degradation for the period from July 2005 - June 2006. Point data were interpolated across individual bank segments to produce an estimate of soil erosion volume. Measured soil bulk densities were then used to calculate the estimated mass loading to Stroubles Creek from channel degradation. Two watershed models and one reach-scale model were developed to predict sediment loading to the stream channel from channel degradation. The Generalized Watershed Loading Function (GWLF) was selected to represent watershed models with limited channel degradation process detail; the Soil and Water Assessment Tool (SWAT) represented the level of channel degradation detail seen in the majority of watershed models; and the CONservation Channel Evolution and Pollutant Transport System (CONCEPTS) reach-scale model was used to evaluate the effectiveness of a detailed process model. Monthly model predictions were compared to retreat rates measured using the erosion pin network. Sediment loading to the stream from bank retreat was estimated as 41 tonnes/yr, based on erosion pin measurements. GWLF, SWAT, and CONCEPTS predicted stream channel sediment contributions of 8 tonnes/yr, 1500 tonnes/yr and 4 tonnes/yr, respectively. Theil-Sen non-parametric simple linear regression was used to test agreement between monthly model predictions and erosion pin estimates. No significant agreement was found between any model predictions and measured retreat, using a conservative a-value of 0.2. GWLF model predictions underpredicted measured channel degradation, but most closely approximated observed data. This result is likely due to similarities in climate and watershed characteristics for the Stroubles Creek watershed and the Pennsylvania watershed used in the empirical model development. SWAT predicted retreat rates exceeded measured values by two orders of magnitude. This result is explained by the inability of SWAT to predict daily flow and sediment discharge. Highly sensitive channel degradation parameters and the lack of calibration data also contributed to SWAT simulation error. CONCEPTS simulation predicted monthly retreat rates slightly less than GWLF. The lack of agreement between CONCEPTS simulation and observed data was mainly the result of limited input data availability. SWAT daily discharge predictions were used as CONCEPTS input data and likely contributed to poor model agreement. Poor estimation of sensitive sediment input parameters may have also contributed to underpredictions by CONCEPTS. Results showed the potential of screening-level watershed models in channel degradation prediction and the importance of flow and sediment time series discharge data in detailed process-based simulation. The limited flexibility of the GWLF channel degradation algorithm makes it unsuitable for evaluating the effects of stream restoration. SWAT and CONCEPTS should only be used for evaluation if appropriate input data are available. Future research will focus on the development of a long-term flow and sediment monitoring data set. Few long-term data sets of this nature exist, making channel degradation modeling difficult. Development of long-term data will allow more accurate modeling and better assessment of channel restoration impacts on channel degradation. Further modeling with GWLF in geographic regions outside the Eastern United States is also needed to determine the scope of applicability of the GWLF channel degradation empirical relationship. Additional research should also focus on the significance of subaerial processes for watersheds of various sizes and on the development of algorithms to simulate these processes. / Master of Science

A key question related to climate change projections is how will aquatic systems respond to changes in variables such as temperature and precipitation? This thesis uses GWLF, a simple catchment scale model to explore potential impacts of climate change on water quantity and quality. River discharge and nutrient loads were modelled for several warmer world scenarios. For one catchment in New York, USA changes in annual dissolved phosphorus loads decreased along with annual streamflow, and particulate phosphorus loads increased for a single future climate scenario. For the Galten catchment of Lake Mälaren, Sweden, the spring melt peak observed historically was reduced for six future scenarios. Peak runoff and dissolved phosphorus and nitrogen load maxima occurred in winter rather than early spring. A synthesis of model results for dissolved inorganic nitrogen (DIN) loadings for five European catchments showed changes in the timing and magnitude of peak DIN load for several future scenarios. In northern Europe, changes were largely due to increased winter streamflow and reduced snow pack and spring melt runoff. In western Europe, DIN loads increased in winter or early spring due to increased precipitation. The biological response for a warmer future scenario was modelled for the Galten basin of Lake Mälaren using GWLF coupled to a lake phytoplankton model and a physical lake model. An increase in cyanobacteria biomass accompanied by a decline in diatom biomass resulted from changes in the timing of nutrient export from the catchment. The projected increase in lake temperature favored an overall increase in total phytoplankton biomass. Lastly, a method based on hourly measurements of colored dissolved organic matter (CDOM) fluorescence provides the appropriate data for dissolved organic carbon (DOC) model parameterization and may also be used for surveillance of organic matter inputs to lakes.

<p>Övergödning utgör ett av de allvarligaste hoten mot Östersjöns miljö. Storleken av näringsbelastningen till havet behöver därför bestämmas med hjälp av tillgängliga matematiska modeller. Modellen ”Generalised Watershed Loading Functions” (GWLF), en ickedistribuerad parametermodell som uppskattar hydrologi och månatlig näringsbelastning, tillämpades på avrinningsområdet till Daugava som mynnar i Östersjön. Syftet med studien var att genom modellering av historisk transport av näringsämnen till Östersjön ta fram parametrar och indata som sedan kan användas vid applicering av GWLF på omkringliggande avrinningsområden. Data från 1990-talet användes för kalibrering av modellen och data från 1980-talet för validering. Årlig kvävebelastning modellerades med R2värdet 0,78 för kalibreringsperioden. Modellerad årlig kvävebelastning för valideringsperioden underskattades med ungefär 30 % vilket troligen beror på att kvävekoncentrationer i grundvatten och ytavrinning minskade mellan 1980- och 1990-talen.</p><p>Fosforbelastningen underskattades jämfört med rapporterade värden vilket troligen beror på att enskilda avlopp inte inkluderades och att rapporterade punktutsläpp är för låga.</p><p>Modifikationer av modellen föreslås för prediktion av näringsbelastningar under lång tid och behovet av harmoniserad, uppdaterad och lättillgänglig data för näringstransportsmodellering diskuteras.</p> / <p>Eutrophication is one of the most serious threats to the Baltic Sea environment. Nutrient loading into the sea therefore needs to be quantified by available mathematical models. The Generalised Watershed Loading Functions (GWLF), a lumpedparameter model that predicts hydrology and monthly nutrient loads, was applied to the Daugava River Basin, discharging into the Baltic Sea. The aim of the study was to model historic transport of nutrients into the Baltic Sea and thereby produce estimates of parameters and input data needed for a spatial extension of the GWLF to surrounding river basins.</p><p>Calibration data were taken from the 1990’s and validation data from the 1980’s. Yearly nitrogen loads were modelled with an R2 value of 0.78 for the calibration period. Predicted yearly nitrogen loads for the validation period were about 30 % lower than reported values, probably depending on decreasing groundwater and runoff concentrations between the 1980’s and 1990’s. Phosphorus loads were underestimated compared to reported values, the main reason probably being the exclusion of septic systems and too low reported point sources.</p><p>Modifications of the model are suggested for longterm predictions of nutrient loads and the need for harmonised, uptodate and generally accessible data for nutrient transport modelling discussed.</p>

A escassez de dados hidrológicos no Brasil é um problema recorrente em muitas regiões, principalmente em se tratando de dados hidrométricos, produção de sedimentos e qualidade da água. A pesquisa por modelos de bacias hidrográficas tem aumentado nas últimas décadas, porém, a estimativa de dados hidrossedimentológicos a partir de modelos mais sofisticados demanda de grande número de variáveis, que devem ser ajustadas para cada sistema natural, o que dificulta a sua aplicação. O objetivo principal desta tese foi avaliar diferentes ferramentas de modelagem utilizadas para a estimativa da vazão, produção de sedimentos e qualidade da água e, em particular, comparar os resultados obtidos de um modelo hidrológico físico semi-distribuído, o Soil Water Assessment Tool (SWAT) com os resultados obtidos a partir de modelos hidrológicos concentrados, com base na metodologia do número da curva de escoamento do Soil Conservation Service (SCS-CN) e no modelo Generalized Watershed Loading Function (GWLF). Buscou-se avaliar e apresentar em quais condições o uso de cada modelo deve ser recomendado, ou seja, quando o esforço necessário para executar o modelo semi-distribuído leva a melhores resultados efetivos. Em relação à simulação da vazão, os resultados dos dois modelos foram altamente influenciados pelos dados de precipitação, indicando que existem, possivelmente, falhas ou erros de medição que poderiam ter influenciado negativamente os resultados. Portanto, foi proposto aplicar o modelo semi-distribuído com dados de precipitação interpolados (DPI) de alta resolução para verificar a eficiência de seus resultados em comparação com os resultados obtidos com a utilização dos dados de precipitação observados (DPO). Para simulação da produção de sedimentos, e das concentrações de nitrogênio e fósforo, o SWAT realiza uma simulação hidrológica mais detalhada, portanto, fornece resultados ligeiramente melhores para parâmetros de qualidade da água. O uso do modelo semi-distribuído também foi ampliado para simular uma bacia hidrográfica sob a influência do reservatório, a fim de verificar a potencialidade do modelo para esse propósito. Os modelos também foram aplicados para identificar quais os impactos potenciais das mudanças no uso do solo previstas e em andamento. Os cenários estudados foram: I &#8211; cenário atual, II &#8211; cenário tendencial, com o aumento da mancha urbana e substituição do solo exposto e de parte da mata nativa por uso agrícola; III &#8211; cenário desejável, complementa o crescimento urbano tendencial com aumento de áreas de reflorestamento. As metodologias foram aplicadas em duas bacias hidrográficas localizadas no Sudeste do Brasil. A primeira é a bacia do rio Jacaré-Guaçu, incluída na Unidade de Gerenciamento de Recursos Hídricos 13 (UGRHI-13), a montante da confluência do rio das Cruzes, com uma área de 1934 km2. O segundo caso de estudo, é a bacia do rio Atibaia, inserida na UGRHI-5, tem uma área de 2817,88 km2 e abrange municípios dos estados de São Paulo e Minas Gerais. Como principal conclusão, o desempenho do modelo semi-distribuído para estimar a produção de sedimentos, e as concentrações de nitrogênio e fósforo foi ligeiramente melhor do que as simulações do modelo concentrado SCS-CN e GWLF, mas essa vantagem pode não compensar o esforço adicional de calibrá-lo e validá-lo. / The lack of hydrological data in Brazil is a recurrent problem in many regions, especially in hydrometric data, sediment yield and water quality. The research by simplified models has increased in the last decades, however, the estimation of hydrossedimentological data from these more sophisticated models demands many variables, which must be adjusted for each natural system, which makes it difficult to apply. At times it is necessary to respond quickly without much precision in the results, in these situations, simpler models with few parameters can be the solution. The objective of this research is to evaluate different modelling tools used estimate streamflow, sediments yield and nutrients loads values, and namely to compare the results obtained from a physically-based distributed hydrological model (SWAT) with the results from a lumped hydrological, the Soil Conservation Service (SCS-CN) and the Generalized Watershed Loading Function (GWLF) model. Both models use the curve number (CN) concept, determined from land use, soil hydrologic group and antecedent soil moisture conditions and were run with a daily time step. We are particularly interested in understanding under which conditions the use of each model is to be recommended, namely when does the addition effort required to run the distributed model leads to effective better results. The input variables and parameters of the lumped model are assumed constant throughout the watershed, while the SWAT model performs the hydrological analysis at a small unit level, designated as hydrological response units (HRUs), and integrates the results at a sub-basin level. In relation to the flow simulation, the results of the two models were highly influenced by the rainfall data, indicating that, possibly, faults or measurement errors could have negatively influenced the results. Therefore, it was proposed to apply the distributed model with high-resolution grids of daily precipitation to verify the efficiency of its results when compared to rainfall data. For simulation of sediment, nitrogen and phosphorus, SWAT performs a more detailed simulation and thus provides slightly better results. The use of the SWAT was also extended to simulate the influence of reservoir, in order to verify the potentiality of the model, in relation to the simulation. The models also were used to identify which are potential impacts of the ongoing land use changes. The scenarios were: I - Current scenario, II - trend scenario, with the increase of urban land and replacement of the exposed soil and part of the native forest by agricultural use; III - desirable scenario complements the trend urban growth with the replacement of exposed soil and part of the agricultural use by reforestation. The methodologies were applied on two watersheds located in the Southeast of Brazil. The first one is the Jacaré-Guaçu river basin, included in the Water Resources Management Unit 13 (UGRHI-13), upstream of Cruzes river confluence, with an area of 1934 km2. The second watershed is the Atibaia River Basin, a part of Water Resources Management Unit 5 (UGRHI-5). It has an area of 2817.88 km2 and covers municipalities of the states of São Paulo and Minas Gerais.

A escassez de dados hidrológicos no Brasil é um problema recorrente em muitas regiões, principalmente em se tratando de dados hidrométricos, produção de sedimentos e qualidade da água. A pesquisa por modelos de bacias hidrográficas tem aumentado nas últimas décadas, porém, a estimativa de dados hidrossedimentológicos a partir de modelos mais sofisticados demanda de grande número de variáveis, que devem ser ajustadas para cada sistema natural, o que dificulta a sua aplicação. O objetivo principal desta tese foi avaliar diferentes ferramentas de modelagem utilizadas para a estimativa da vazão, produção de sedimentos e qualidade da água e, em particular, comparar os resultados obtidos de um modelo hidrológico físico semi-distribuído, o Soil Water Assessment Tool (SWAT) com os resultados obtidos a partir de modelos hidrológicos concentrados, com base na metodologia do número da curva de escoamento do Soil Conservation Service (SCS-CN) e no modelo Generalized Watershed Loading Function (GWLF). Buscou-se avaliar e apresentar em quais condições o uso de cada modelo deve ser recomendado, ou seja, quando o esforço necessário para executar o modelo semi-distribuído leva a melhores resultados efetivos. Em relação à simulação da vazão, os resultados dos dois modelos foram altamente influenciados pelos dados de precipitação, indicando que existem, possivelmente, falhas ou erros de medição que poderiam ter influenciado negativamente os resultados. Portanto, foi proposto aplicar o modelo semi-distribuído com dados de precipitação interpolados (DPI) de alta resolução para verificar a eficiência de seus resultados em comparação com os resultados obtidos com a utilização dos dados de precipitação observados (DPO). Para simulação da produção de sedimentos, e das concentrações de nitrogênio e fósforo, o SWAT realiza uma simulação hidrológica mais detalhada, portanto, fornece resultados ligeiramente melhores para parâmetros de qualidade da água. O uso do modelo semi-distribuído também foi ampliado para simular uma bacia hidrográfica sob a influência do reservatório, a fim de verificar a potencialidade do modelo para esse propósito. Os modelos também foram aplicados para identificar quais os impactos potenciais das mudanças no uso do solo previstas e em andamento. Os cenários estudados foram: I &#8211; cenário atual, II &#8211; cenário tendencial, com o aumento da mancha urbana e substituição do solo exposto e de parte da mata nativa por uso agrícola; III &#8211; cenário desejável, complementa o crescimento urbano tendencial com aumento de áreas de reflorestamento. As metodologias foram aplicadas em duas bacias hidrográficas localizadas no Sudeste do Brasil. A primeira é a bacia do rio Jacaré-Guaçu, incluída na Unidade de Gerenciamento de Recursos Hídricos 13 (UGRHI-13), a montante da confluência do rio das Cruzes, com uma área de 1934 km2. O segundo caso de estudo, é a bacia do rio Atibaia, inserida na UGRHI-5, tem uma área de 2817,88 km2 e abrange municípios dos estados de São Paulo e Minas Gerais. Como principal conclusão, o desempenho do modelo semi-distribuído para estimar a produção de sedimentos, e as concentrações de nitrogênio e fósforo foi ligeiramente melhor do que as simulações do modelo concentrado SCS-CN e GWLF, mas essa vantagem pode não compensar o esforço adicional de calibrá-lo e validá-lo. / The lack of hydrological data in Brazil is a recurrent problem in many regions, especially in hydrometric data, sediment yield and water quality. The research by simplified models has increased in the last decades, however, the estimation of hydrossedimentological data from these more sophisticated models demands many variables, which must be adjusted for each natural system, which makes it difficult to apply. At times it is necessary to respond quickly without much precision in the results, in these situations, simpler models with few parameters can be the solution. The objective of this research is to evaluate different modelling tools used estimate streamflow, sediments yield and nutrients loads values, and namely to compare the results obtained from a physically-based distributed hydrological model (SWAT) with the results from a lumped hydrological, the Soil Conservation Service (SCS-CN) and the Generalized Watershed Loading Function (GWLF) model. Both models use the curve number (CN) concept, determined from land use, soil hydrologic group and antecedent soil moisture conditions and were run with a daily time step. We are particularly interested in understanding under which conditions the use of each model is to be recommended, namely when does the addition effort required to run the distributed model leads to effective better results. The input variables and parameters of the lumped model are assumed constant throughout the watershed, while the SWAT model performs the hydrological analysis at a small unit level, designated as hydrological response units (HRUs), and integrates the results at a sub-basin level. In relation to the flow simulation, the results of the two models were highly influenced by the rainfall data, indicating that, possibly, faults or measurement errors could have negatively influenced the results. Therefore, it was proposed to apply the distributed model with high-resolution grids of daily precipitation to verify the efficiency of its results when compared to rainfall data. For simulation of sediment, nitrogen and phosphorus, SWAT performs a more detailed simulation and thus provides slightly better results. The use of the SWAT was also extended to simulate the influence of reservoir, in order to verify the potentiality of the model, in relation to the simulation. The models also were used to identify which are potential impacts of the ongoing land use changes. The scenarios were: I - Current scenario, II - trend scenario, with the increase of urban land and replacement of the exposed soil and part of the native forest by agricultural use; III - desirable scenario complements the trend urban growth with the replacement of exposed soil and part of the agricultural use by reforestation. The methodologies were applied on two watersheds located in the Southeast of Brazil. The first one is the Jacaré-Guaçu river basin, included in the Water Resources Management Unit 13 (UGRHI-13), upstream of Cruzes river confluence, with an area of 1934 km2. The second watershed is the Atibaia River Basin, a part of Water Resources Management Unit 5 (UGRHI-5). It has an area of 2817.88 km2 and covers municipalities of the states of São Paulo and Minas Gerais.

Övergödning utgör ett av de allvarligaste hoten mot Östersjöns miljö. Storleken av näringsbelastningen till havet behöver därför bestämmas med hjälp av tillgängliga matematiska modeller. Modellen ”Generalised Watershed Loading Functions” (GWLF), en ickedistribuerad parametermodell som uppskattar hydrologi och månatlig näringsbelastning, tillämpades på avrinningsområdet till Daugava som mynnar i Östersjön. Syftet med studien var att genom modellering av historisk transport av näringsämnen till Östersjön ta fram parametrar och indata som sedan kan användas vid applicering av GWLF på omkringliggande avrinningsområden. Data från 1990-talet användes för kalibrering av modellen och data från 1980-talet för validering. Årlig kvävebelastning modellerades med R2värdet 0,78 för kalibreringsperioden. Modellerad årlig kvävebelastning för valideringsperioden underskattades med ungefär 30 % vilket troligen beror på att kvävekoncentrationer i grundvatten och ytavrinning minskade mellan 1980- och 1990-talen. Fosforbelastningen underskattades jämfört med rapporterade värden vilket troligen beror på att enskilda avlopp inte inkluderades och att rapporterade punktutsläpp är för låga. Modifikationer av modellen föreslås för prediktion av näringsbelastningar under lång tid och behovet av harmoniserad, uppdaterad och lättillgänglig data för näringstransportsmodellering diskuteras. / Eutrophication is one of the most serious threats to the Baltic Sea environment. Nutrient loading into the sea therefore needs to be quantified by available mathematical models. The Generalised Watershed Loading Functions (GWLF), a lumpedparameter model that predicts hydrology and monthly nutrient loads, was applied to the Daugava River Basin, discharging into the Baltic Sea. The aim of the study was to model historic transport of nutrients into the Baltic Sea and thereby produce estimates of parameters and input data needed for a spatial extension of the GWLF to surrounding river basins. Calibration data were taken from the 1990’s and validation data from the 1980’s. Yearly nitrogen loads were modelled with an R2 value of 0.78 for the calibration period. Predicted yearly nitrogen loads for the validation period were about 30 % lower than reported values, probably depending on decreasing groundwater and runoff concentrations between the 1980’s and 1990’s. Phosphorus loads were underestimated compared to reported values, the main reason probably being the exclusion of septic systems and too low reported point sources. Modifications of the model are suggested for longterm predictions of nutrient loads and the need for harmonised, uptodate and generally accessible data for nutrient transport modelling discussed.