Watershed-scale runoff routing and solute transport in a spatially aggregated hydrological framework

Dissertation submitted in partial fulfillment of the requirements for the Degree of Master of Science in Geospatial Technologies

Runoff and sediment transport in a degraded area

Gully erosion occurs by the combined action of splash, sheetwash and rill-wash (interrill and rill erosion). These erosion processes have a great capacity for both sediment production and sediment transport. The objectives of this experiment were to evaluate hydrological and sediment transport in a degraded area, severely dissected by gullies; to assess the hydraulic flow characteristics and their aggregate transport capacity; and to measure the initial splash erosion rate. In the study area in Guarapuava, State of Paraná, Brazil (lat 25º 24' S; long 51º24' W; 1034 m asl), the soil was classified as Cambissolo Húmico alumínico, with the following particle-size composition: sand 0.116 kg kg-1; silt 0.180 kg kg-1; and clay 0.704 kg kg-1. The approach of this research was based on microcatchments formed in the ground, to study the hydrological response and sediment transport. A total of eight rill systems were simulated with dry and wet soil. An average rainfall of 33.7 ± 4.0 mm was produced for 35 to 54 min by a rainfall simulator. The equipment was installed, and a trough was placed at the end of the rill to collect sediments and water. During the simulation, the following variables were measured: time to runoff, time to ponding, time of recession, flow velocity, depth, ratio of the initial splash and grain size. The rainsplash of dry topsoil was more than twice as high as under moist conditions (5 g m-2 min-1 and 2 g m-2 min-1, respectively). The characteristics of the flow hydraulics indicate transition from laminar to turbulent flow [Re (Reynolds number) 1000-2000]. In addition, it was observed that a flow velocity of 0.12 m s-1 was the threshold for turbulent flow (Re > 2000), especially at the end of the rainfall simulation. The rill flow tended to be subcritical [Fr (Froude Number) < 1.0]. The variation in hydrological attributes (infiltration and runoff) was lower, while the sediment yield was variable. The erosion in the rill systems was characterized as limited transport, although the degraded area generated an average of 394 g m-2 of sediment in each simulation.

Hydrological processes obtained on the plot scale under four simulated rainfall tests during the cycle of different crop systems

The cropping system influences the interception of water by plants, water storage in depressions on the soil surface, water infiltration into the soil and runoff. The aim of this study was to quantify some hydrological processes under no tillage cropping systems at the edge of a slope, in 2009 and 2010, in a Humic Dystrudept soil, with the following treatments: corn, soybeans, and common beans alone; and intercropped corn and common bean. Treatments consisted of four simulated rainfall tests at different times, with a planned intensity of 64 mm h-1 and 90 min duration. The first test was applied 18 days after sowing, and the others at 39, 75 and 120 days after the first test. Different times of the simulated rainfall and stages of the crop cycle affected soil water content prior to the rain, and the time runoff began and its peak flow and, thus, the surface hydrological processes. The depth of the runoff and the depth of the water intercepted by the crop + soil infiltration + soil surface storage were affected by the crop systems and the rainfall applied at different times. The corn crop was the most effective treatment for controlling runoff, with a water loss ratio of 0.38, equivalent to 75 % of the water loss ratio exhibited by common bean (0.51), the least effective treatment in relation to the others. Total water loss by runoff decreased linearly with an increase in the time that runoff began, regardless of the treatment; however, soil water content on the gravimetric basis increased linearly from the beginning to the end of the rainfall.

Statistical methods for river runoff prediction

Methods used to analyze one type of nonstationary stochastic processes?the periodically correlated process?are considered. Two methods of one-step-forward prediction of periodically correlated time series are examined. One-step-forward predictions made in accordance with an autoregression model and a model of an artificial neural network with one latent neuron layer and with an adaptation mechanism of network parameters in a moving time window were compared in terms of efficiency. The comparison showed that, in the case of prediction for one time step for time series of mean monthly water discharge, the simpler autoregression model is more efficient.

Impacts of upland open drains upon runoff generation: a numerical assessment of catchment-scale impacts

Shallow upland drains, grips, have been hypothesized as responsible for increased downstream flow magnitudes. Observations provide counterfactual evidence, often relating to the difficulty of inferring conclusions from statistical correlation and paired catchment comparisons, and the complexity of designing field experiments to test grip impacts at the catchment scale. Drainage should provide drier antecedent moisture conditions, providing more storage at the start of an event; however, grips have higher flow velocities than overland flow, thus potentially delivering flow more rapidly to the drainage network. We develop and apply a model for assessing the impacts of grips on flow hydrographs. The model was calibrated on the gripped case, and then the gripped case was compared with the intact case by removing all grips. This comparison showed that even given parameter uncertainty, the intact case had significantly higher flood peaks and lower baseflows, mirroring field observations of the hydrological response of intact peat. The simulations suggest that this is because delivery effects may not translate into catchment-scale impacts for three reasons. First, in our case, the proportions of flow path lengths that were hillslope were not changed significantly by gripping. Second, the structure of the grip network as compared with the structure of the drainage basin mitigated against grip-related increases in the concentration of runoff in the drainage network, although it did marginally reduce the mean timing of that concentration at the catchment outlet. Third, the effect of the latter upon downstream flow magnitudes can only be assessed by reference to the peak timing of other tributary basins, emphasizing that drain effects are both relative and scale dependent. However, given the importance of hillslope flow paths, we show that if upland drainage causes significant changes in surface roughness on hillslopes, then critical and important feedbacks may impact upon the speed of hydrological response. Copyright (c) 2012 John Wiley & Sons, Ltd.

Biogeochemical processes controlling oxygen and carbon isotopes of diatom silica in Late Glacial to Holocene lacustrine rhythmites

Biogeochemical cycles and sedimentary records in lakes are related to climate controls on hydrology and catchment processes. Changes in the isotopic imposition of the diatom frustules (δ 18 O diatom and δ 13 C diatom ) in lacustrine sediments can be used to reconstruct palaeoclimatic and palaeoenvironmental changes. The Lago Chungará (Andean Altiplano, 18°15 ′ S, 69°10 ′ W, 4520 masl) diatomaceous laminated sediments are made up of white and green multiannual rhythmites. White laminae were formed during short-term diatom super-blooms, and are composed almost exclusively of large-sized Cyclostephanos andinus.These diatoms bloom during mixing events when recycled nutrients from the bottom waters are brought to the surface and/or when nutrients are introduced from the catchment during periods of strong runoff. Conversely, the green laminae are thought to have been deposited over several years and are composed of a mixture of diatoms (mainly smaller valves of C. andinus and Discostella stelligera ) and organic matter. These green laminae reflect the lake's hydrological recovery from a status favouring the diatom super-blooms (white laminae) towards baseline conditions. δ 18 O diatom and δ 13 C diatom from 11,990 to 11,530 cal years BP allow us to reconstruct shifts in the precipitation/evaporation ratio and changes in the lake water dissolved carbon concentration, respectively. δ 18 O diatom values indicate that white laminae formation occurred mainly during low lake level stages, whereas green laminae formation generally occurred during high lake level stages. The isotope and chronostratigraphical data together suggest that white laminae deposition is caused by extraordinary environmental events. El Niño-Southern Oscillation and changes in solar activity are the most likely climate forcing mechanisms that could trigger such events, favouring hydrological changes at interannual-to-decadal scale. This study demonstrates the potential for laminated lake sediments to document extreme pluriannual events.

Rainfall and runoff water quality of the Pang and Lambourn, tributaries of the River Thames, south-eastern England

The water quality of rainfall and runoff is described for two catchments of two tributaries of the River Thames, the Pang and Lambourn. Rainfall chemistry is variable and concentrations of most determinands decrease with increasing volume of catch probably due to 'wash out' processes. Two rainfall sites have been monitored, one for each catchment. The rainfall site on the Lambourn shows higher chemical concentrations than the one for the Pang which probably reflects higher amounts of local inputs from agricultural activity, Rainfall quality data at a long-term rainfall site on the Pang (UK National Air Quality Archive) shows chemistries similar to that for the Lambourn site. but with some clear differences. Rainfall chemistries show considerable variation on an event-to-event basis. Average water quality concentrations and flow-weighted concentrations as well as fluxes vary across the sites, typically by about 30%. Stream chemistry is much less variable due to the main Source of water coming from aquifer sources of high storage. The relationship between rainfall and runoff chemistry at the catchment outlet is described in terms of the relative proportions of atmospheric and within-catchment sources. Remarkably, in view of the quantity of agricultural and sewage inputs to the streams, the catchments appear to be retaining both P and N.

Simulating current global river runoff with a global hydrological model: model revisions, validation, and sensitivity analysis

Global hydrological models (GHMs) model the land surface hydrologic dynamics of continental-scale river basins. Here we describe one such GHM, the Macro-scale - Probability-Distributed Moisture model.09 (Mac-PDM.09). The model has undergone a number of revisions since it was last applied in the hydrological literature. This paper serves to provide a detailed description of the latest version of the model. The main revisions include the following: (1) the ability for the model to be run for n repetitions, which provides more robust estimates of extreme hydrological behaviour, (2) the ability of the model to use a gridded field of coefficient of variation (CV) of daily rainfall for the stochastic disaggregation of monthly precipitation to daily precipitation, and (3) the model can now be forced with daily input climate data as well as monthly input climate data. We demonstrate the effects that each of these three revisions has on simulated runoff relative to before the revisions were applied. Importantly, we show that when Mac-PDM.09 is forced with monthly input data, it results in a negative runoff bias relative to when daily forcings are applied, for regions of the globe where the day-to-day variability in relative humidity is high. The runoff bias can be up to - 80% for a small selection of catchments but the absolute magnitude of the bias may be small. As such, we recommend future applications of Mac-PDM.09 that use monthly climate forcings acknowledge the bias as a limitation of the model. The performance of Mac-PDM.09 is evaluated by validating simulated runoff against observed runoff for 50 catchments. We also present a sensitivity analysis that demonstrates that simulated runoff is considerably more sensitive to method of PE calculation than to perturbations in soil moisture and field capacity parameters.

Hillslope scale surface runoff, sediment, and nutrient losses associated with tramline wheelings

Research on arable sandy loam and silty clay loam soils on 4° slopes in England has shown that tramlines (i.e. the unseeded wheeling areas used to facilitate spraying operations in cereal crops) can represent the most important pathway for phosphorus and sediment loss from moderately sloping fields. Detailed monitoring over the October–March period in winters 2005–2006 and 2006–2007 included event-based sampling of surface runoff, suspended and particulate sediment, and dissolved and particulate phosphorus from hillslope segments (each ∼300–800 m2) established in a randomized block design with four replicates of each treatment at each of two sites on lighter and heavier soils. Experimental treatments assessed losses from the cropped area without tramlines, and from the uncropped tramline area, and were compared to losses from tramlines which had been disrupted once in the autumn with a shallow tine. On the lighter soil, the effects of removal or shallow incorporation of straw residues was also determined. Research on both sandy and silty clay loam soils across two winters showed that tramline wheelings represented the dominant pathway for surface runoff and transport of sediment, phosphorus and nitrogen from cereal crops on moderate slopes. Results indicated 5·5–15·8% of rainfall lost as runoff, and losses of 0·8–2·9 kg TP ha−1 and 0·3–4·8 t ha−1 sediment in tramline treatments, compared to only 0·2–1·7% rainfall lost as runoff, and losses of 0·0–0·2 kg TP ha−1 and 0·003–0·3 t ha−1 sediment from treatments without tramlines or those where tramlines had been disrupted. The novel shallow disruption of tramline wheelings using a tine once following the autumn spray operation consistently and dramatically reduced (p < 0·001) surface runoff and loads of sediment, total nitrogen and total phosphorus to levels similar to those measured in cropped areas between tramlines. Results suggest that options for managing tramline wheelings warrant further refinement and evaluation with a view to incorporating them into spatially-targeted farm-level management planning using national or catchment-based agri-environment policy instruments aimed at reducing diffuse pollution from land to surface water systems. Copyright © 2010 John Wiley & Sons, Ltd.

PERSiST: a flexible rainfall-runoff modelling toolkit for use with the INCA family of models

Runoff generation processes and pathways vary widely between catchments. Credible simulations of solute and pollutant transport in surface waters are dependent on models which facilitate appropriate, catchment-specific representations of perceptual models of the runoff generation process. Here, we present a flexible, semi-distributed landscape-scale rainfall-runoff modelling toolkit suitable for simulating a broad range of user-specified perceptual models of runoff generation and stream flow occurring in different climatic regions and landscape types. PERSiST (the Precipitation, Evapotranspiration and Runoff Simulator for Solute Transport) is designed for simulating present-day hydrology; projecting possible future effects of climate or land use change on runoff and catchment water storage; and generating hydrologic inputs for the Integrated Catchments (INCA) family of models. PERSiST has limited data requirements and is calibrated using observed time series of precipitation, air temperature and runoff at one or more points in a river network. Here, we apply PERSiST to the river Thames in the UK and describe a Monte Carlo tool for model calibration, sensitivity and uncertainty analysis

Investigating hydrological regimes and processes in a set of catchments with temporary waters in Mediterranean Europe

Seven catchments of diverse size in Mediterranean Europe were investigated in order to understand the main aspects of their hydrological functioning. The methods included the analysis of daily and monthly precipitation, monthly potential evapotranspiration rates, flow duration curves, rainfall runoff relationships and catchment internal data for the smaller and more instrumented catchments. The results showed that the catchments were less dry than initially considered. Only one of them was really semi-arid throughout the year. All the remaining catchments showed wet seasons when precipitation exceeded potential evapotrans-piration, allowing aquifer recharge, wet runoff generation mechanisms and relevant baseflow contribution. Nevertheless, local infiltration excess (Hortonian) overland flow was inferred during summer storms in some catchments and urban overland flow in some others. The roles of karstic groundwater, human disturbance and low winter temperatures were identified as having an important impact on the hydrological regime in some of the catchments.

Connectivity processes and riparian vegetation of the upper Paraná River, Brazil

In fluvial systems, the relationship between a dominant variable (e.g. flood pulse) and its dependent ones (e.g. riparian vegetation) is called connectivity. This paper analyzes the connectivity elements and processes controlling riparian vegetation for a reach of the upper Paraná River (Brazil) and estimates the future changes in channel-vegetation relationship as a consequence of the managing of a large dam. The studied reach is situated 30km downstream from the Porto Primavera Dam (construction finished in 1999). Through aerial photography (1:25,000, 1996), RGB-CBERS satellite imagery and a previous field botany survey it was possible to elaborate a map with the five major morpho-vegetation units: 1) Tree-dominated natural levee, 2) Shrubby upper floodplain, 3) Shrub-herbaceous mid floodplain, 4) Grass-herbaceous lower floodplain and 5) Shrub-herbaceous flood runoff channel units. By use of a detailed topographic survey and statistical tools each morpho-vegetation type was analyzed according to its connectivity parameters (frequency, recurrence, permanence, seasonality, potamophase, limnophase and FCQ index) in the pre- and post-dam closure periods of the historical series. Data showed that most of the morpho-vegetation units were predicted to present changes in connectivity parameters values after dam closing and the new regime could affect, in different intensity, the river ecology and particularly the riparian vegetation. The methods used in this study can be useful for dam impact studies in other South American tropical rivers. © 2012 Elsevier Ltd.

Spatial distribution of runoff depth from a watershed located in a cuesta relief area of São Paulo State, Brazil

This study was undertaken in a 1566 ha drainage basin situated in an area with cuesta relief in the state of São Paulo, Brazil. The objectives were: 1) to map the maximum potential soil water retention capacity, and 2) to simulate the depth of surface runoff in each geographical position of the area based on a typical rainfall event. The database required for the development of this research was generated in the environment of the geographical information system ArcInfo v.10.1. Undeformed soil samples were collected at 69 points. The ordinary kriging method was used in the interpolation of the values of soil density and maximum potential soil water retention capacity. The spherical model allowed for better adjustment of the semivariograms corresponding to the two soil attributes for the depth of 0 to 20 cm, while the Gaussian model enabled a better fit of the spatial behavior of the two variables for the depth of 20 to 40 cm. The simulation of the spatial distribution revealed a gradual increase in the depth of surface runoff for the rainfall event taken as example (25 mm) from the reverse to the peripheral depression of the cuesta (from west to east). There is a positive aspect observed in the gradient, since the sites of highest declivity, especially those at the front of the cuesta, are closer to the western boundary of the watershed where the lowest depths of runoff occur. This behavior, in conjunction with certain values of erodibility and depending on the land use and cover, can help mitigate the soil erosion processes in these areas.

Potential impacts of polar fronts on sedimentation processes at Abrolhos coral reef (South-West Atlantic Ocean/Brazil)

Several reefs of the world have undergone changes in community due to sedimentation processes. It has been suggested that Abrolhos reefs (Brazil/South-West Atlantic) are subjected to a steady coastal influence, although there is still little information regarding this assumption. In this work, we have analyzed a set of environmental parameters concerning sedimentation characteristics at the Abrolhos reefs, near 18 degrees S-39 degrees W. The analysis included remote sensing, model and in situ data to provide a three-dimensional quantitative description of the processes that influence sediment apportionment to the reefs. Mineralogy and natural radioactivity of sediment trapped at three reef sites in a transect perpendicular to the coastline were used in conjunction with numerical weather prediction model and remote sensing databases. We have observed an increase of around 100% of sediment flux during the summer compared to the winter season. A comparison of regional rainfall regime, sediment plume dynamics and a year-around monitoring of polar fronts trajectories and surface wind showed that the wind-driven resedimentation due to polar front activity is the major contributor to the intensification of sedimentation processes at the offshore area of Abrolhos reefs, despite river runoff from mainland. (C) 2007 Elsevier Ltd. All rights reserved.