SCALING LIMIT FOR A DRAINAGE NETWORK MODEL

We consider the two-dimensional version of a drainage network model introduced ill Gangopadhyay, Roy and Sarkar (2004), and show that the appropriately rescaled family of its paths converges in distribution to the Brownian web. We do so by verifying the convergence criteria proposed in Fontes, Isopi, Newman and Ravishankar (2002).

Lithologic Control on the Scaling Properties of the First-Order Streams of Drainage Networks: A Monofractal Analysis

The movement of water through the landscape can be investigated at different scales. This study dealt with the interrelation between bedrock lithology and the geometry of the overlying drainage systems. Parameters of fractal analysis, such as fractal dimension and lacunarity, were used to measure and quantify this relationship. The interrelation between bedrock lithology and the geometry of the drainage systems has been widely studied in the last decades. The quantification of this linkage has not yet been clearly established. Several studies have selected river basins or regularly shaped areas as study units, assuming them to be lithologically homogeneous. This study considered irregular distributions of rock types, establishing areas of the soil map (1:25,000) with the same lithologic information as study units. The tectonic stability and the low climatic variability of the study region allowed effective investigation of the lithologic controls on the drainage networks developed on the plutonic rocks, the metamorphic rocks, and the sedimentary materials existing in the study area. To exclude the effect of multiple in- and outflows in the lithologically homogeneous units, we focused this study on the first-order streams of the drainage networks. The geometry of the hydrologic features was quantified through traditional metrics of fluvial geomorphology and scaling parameters of fractal analysis, such as the fractal dimension, the reference density, and the lacunarity. The results demonstrate the scale invariance of both the drainage networks and the set of first-order streams at the study scale and a relationship between scaling in the lithology and the drainage network.

Influence of rainfall and catchment characteristics on urban stormwater quality

This thesis presents the outcomes of a comprehensive research study undertaken to investigate the influence of rainfall and catchment characteristics on urban stormwater quality. The knowledge created is expected to contribute to a greater understanding of urban stormwater quality and thereby enhance the design of stormwater quality treatment systems. The research study was undertaken based on selected urban catchments in Gold Coast, Australia. The research methodology included field investigations, laboratory testing, computer modelling and data analysis. Both univariate and multivariate data analysis techniques were used to investigate the influence of rainfall and catchment characteristics on urban stormwater quality. The rainfall characteristics investigated included average rainfall intensity and rainfall duration whilst catchment characteristics included land use, impervious area percentage, urban form and pervious area location. The catchment scale data for the analysis was obtained from four residential catchments, including rainfall-runoff records, drainage network data, stormwater quality data and land use and land cover data. Pollutants build-up samples were collected from twelve road surfaces in residential, commercial and industrial land use areas. The relationships between rainfall characteristics, catchment characteristics and urban stormwater quality were investigated based on residential catchments and then extended to other land uses. Based on the influence rainfall characteristics exert on urban stormwater quality, rainfall events can be classified into three different types, namely, high average intensity-short duration (Type 1), high average intensity-long duration (Type 2) and low average intensity-long duration (Type 3). This provides an innovative approach to conventional modelling which does not commonly relate stormwater quality to rainfall characteristics. Additionally, it was found that the threshold intensity for pollutant wash-off from urban catchments is much less than for rural catchments. High average intensity-short duration rainfall events are cumulatively responsible for the generation of a major fraction of the annual pollutants load compared to the other rainfall event types. Additionally, rainfall events less than 1 year ARI such as 6- month ARI should be considered for treatment design as they generate a significant fraction of the annual runoff volume and by implication a significant fraction of the pollutants load. This implies that stormwater treatment designs based on larger rainfall events would not be feasible in the context of cost-effectiveness, efficiency in treatment performance and possible savings in land area needed. This also suggests that the simulation of long-term continuous rainfall events for stormwater treatment design may not be needed and that event based simulations would be adequate. The investigations into the relationship between catchment characteristics and urban stormwater quality found that other than conventional catchment characteristics such as land use and impervious area percentage, other catchment characteristics such as urban form and pervious area location also play important roles in influencing urban stormwater quality. These outcomes point to the fact that the conventional modelling approach in the design of stormwater quality treatment systems which is commonly based on land use and impervious area percentage would be inadequate. It was also noted that the small uniformly urbanised areas within a larger mixed catchment produce relatively lower variations in stormwater quality and as expected lower runoff volume with the opposite being the case for large mixed use urbanised catchments. Therefore, a decentralised approach to water quality treatment would be more effective rather than an "end-of-pipe" approach. The investigation of pollutants build-up on different land uses showed that pollutant build-up characteristics vary even within the same land use. Therefore, the conventional approach in stormwater quality modelling, which is based solely on land use, may prove to be inappropriate. Industrial land use has relatively higher variability in maximum pollutant build-up, build-up rate and particle size distribution than the other two land uses. However, commercial and residential land uses had relatively higher variations of nutrients and organic carbon build-up. Additionally, it was found that particle size distribution had a relatively higher variability for all three land uses compared to the other build-up parameters. The high variability in particle size distribution for all land uses illustrate the dissimilarities associated with the fine and coarse particle size fractions even within the same land use and hence the variations in stormwater quality in relation to pollutants adsorbing to different sizes of particles.

First flush behaviour in urban residential catchments

The current state of knowledge in relation to first flush does not provide a clear understanding of the role of rainfall and catchment characteristics in influencing this phenomenon. This is attributed to the inconsistent findings from research studies due to the unsatisfactory selection of first flush indicators and how first flush is defined. The research study discussed in this thesis provides the outcomes of a comprehensive analysis on the influence of rainfall and catchment characteristics on first flush behaviour in residential catchments. Two sets of first flush indicators are introduced in this study. These indicators were selected such that they are representative in explaining in a systematic manner the characteristics associated with first flush. Stormwater samples and rainfall-runoff data were collected and recorded from stormwater monitoring stations established at three urban catchments at Coomera Waters, Gold Coast, Australia. In addition, historical data were also used to support the data analysis. Three water quality parameters were analysed, namely, total suspended solids (TSS), total phosphorus (TP) and total nitrogen (TN). The data analyses were primarily undertaken using multi criteria decision making methods, PROMETHEE and GAIA. Based on the data obtained, the pollutant load distribution curve (LV) was determined for the individual rainfall events and pollutant types. Accordingly, two sets of first flush indicators were derived from the curve, namely, cumulative load wash-off for every 10% of runoff volume interval (interval first flush indicators or LV) from the beginning of the event and the actual pollutant load wash-off during a 10% increment in runoff volume (section first flush indicators or P). First flush behaviour showed significant variation with pollutant types. TSS and TP showed consistent first flush behaviour. However, the dissolved fraction of TN showed significant differences to TSS and TP first flush while particulate TN showed similarities. Wash-off of TSS, TP and particulate TN during the first 10% of the runoff volume showed no influence from corresponding rainfall intensity. This was attributed to the wash-off of weakly adhered solids on the catchment surface referred to as "short term pollutants" or "weakly adhered solids" load. However, wash-off after 10% of the runoff volume showed dependency on the rainfall intensity. This is attributed to the wash-off of strongly adhered solids being exposed when the weakly adhered solids diminish. The wash-off process was also found to depend on rainfall depth at the end part of the event as the strongly adhered solids are loosened due to impact of rainfall in the earlier part of the event. Events with high intensity rainfall bursts after 70% of the runoff volume did not demonstrate first flush behaviour. This suggests that rainfall pattern plays a critical role in the occurrence of first flush. Rainfall intensity (with respect to the rest of the event) that produces 10% to 20% runoff volume play an important role in defining the magnitude of the first flush. Events can demonstrate high magnitude first flush when the rainfall intensity occurring between 10% and 20% of the runoff volume is comparatively high while low rainfall intensities during this period produces low magnitude first flush. For events with first flush, the phenomenon is clearly visible up to 40% of the runoff volume. This contradicts the common definition that first flush only exists, if for example, 80% of the pollutant mass is transported in the first 30% of runoff volume. First flush behaviour for TN is different compared to TSS and TP. Apart from rainfall characteristics, the composition and the availability of TN on the catchment also play an important role in first flush. The analysis confirmed that events with low rainfall intensity can produce high magnitude first flush for the dissolved fraction of TN, while high rainfall intensity produce low dissolved TN first flush. This is attributed to the source limiting behaviour of dissolved TN wash-off where there is high wash-off during the initial part of a rainfall event irrespective of the intensity. However, for particulate TN, the influence of rainfall intensity on first flush characteristics is similar to TSS and TP. The data analysis also confirmed that first flush can occur as high magnitude first flush, low magnitude first flush or non existence of first flush. Investigation of the influence of catchment characteristics on first flush found that the key factors that influence the phenomenon are the location of the pollutant source, spatial distribution of the pervious and impervious surfaces in the catchment, drainage network layout and slope of the catchment. This confirms that first flush phenomenon cannot be evaluated based on a single or a limited set of parameters as a number of catchment characteristics should be taken into account. Catchments where the pollutant source is located close to the outlet, a high fraction of road surfaces, short travel time to the outlet, with steep slopes can produce high wash-off load during the first 50% of the runoff volume. Rainfall characteristics have a comparatively dominant impact on the wash-off process compared to the catchment characteristics. In addition, the pollutant characteristics also should be taken into account in designing stormwater treatment systems due to different wash-off behaviour. Analysis outcomes confirmed that there is a high TSS load during the first 20% of the runoff volume followed by TN which can extend up to 30% of the runoff volume. In contrast, high TP load can exist during the initial and at the end part of a rainfall event. This is related to the composition of TP available for the wash-off.

Using a structural approach to identify relationships between soil and erosion in a semi-humid forested area, South India

Biogeochemical and hydrological cycles are currently studied on a small experimental forested watershed (4.5 km(2)) in the semi-humid South India. This paper presents one of the first data referring to the distribution and dynamics of a widespread red soil (Ferralsols and Chromic Luvisols) and black soil (Vertisols and Vertic intergrades) cover, and its possible relationship with the recent development of the erosion process. The soil map was established from the observation of isolated soil profiles and toposequences, and surveys of soil electromagnetic conductivity (EM31, Geonics Ltd), lithology and vegetation. The distribution of the different parts of the soil cover in relation to each other was used to establish the dynamics and chronological order of formation. Results indicate that both topography and lithology (gneiss and amphibolite) have influenced the distribution of the soils. At the downslope, the following parts of the soil covers were distinguished: i) red soil system, ii) black soil system, iii) bleached horizon at the top of the black soil and iv) bleached sandy saprolite at the base of the black soil. The red soil is currently transforming into black soil and the transformation front is moving upslope. In the bottom part of the slope, the chronology appears to be the following: black soil > bleached horizon at the top of the black soil > streambed > bleached horizon below the black soil. It appears that the development of the drainage network is a recent process, which was guided by the presence of thin black soil with a vertic horizon less than 2 in deep. Three distinctive types of erosional landforms have been identified: 1. rotational slips (Type 1); 2. a seepage erosion (Type 2) at the top of the black soil profile; 3. A combination of earthflow and sliding in the non-cohesive saprolite of the gneiss occurs at midslope (Type 3). Types 1 and 2 erosion are mainly occurring downslope and are always located at the intersection between the streambed and the red soil-black soil contact. Neutron probe monitoring, along an area vulnerable to erosion types 1 and 2, indicates that rotational slips are caused by a temporary watertable at the base of the black soil and within the sandy bleached saprolite, which behaves as a plane of weakness. The watertable is induced by the ephemeral watercourse. Erosion type 2 is caused by seepage of a perched watertable, which occurs after swelling and closing of the cracks of the vertic clay horizon and within a light textured and bleached horizon at the top of black soil. Type 3 erosion is not related to the red soil-black soil system but is caused by the seasonal seepage of saturated throughflow in the sandy saprolite of the gneiss occurring at midslope. (c) 2006 Elsevier B.V. All rights reserved.

A general geomorphological recession flow model for river basins

Recession flows in a basin are controlled by the temporal evolution of its active drainage network (ADN). The geomorphological recession flow model (GRFM) assumes that both the rate of flow generation per unit ADN length (q) and the speed at which ADN heads move downstream (c) remain constant during a recession event. Thereby, it connects the power law exponent of -dQ/dt versus Q (discharge at the outlet at time t) curve, , with the structure of the drainage network, a fixed entity. In this study, we first reformulate the GRFM for Horton-Strahler networks and show that the geomorphic ((g)) is equal to D/(D-1), where D is the fractal dimension of the drainage network. We then propose a more general recession flow model by expressing both q and c as functions of Horton-Strahler stream order. We show that it is possible to have = (g) for a recession event even when q and c do not remain constant. The modified GRFM suggests that is controlled by the spatial distribution of subsurface storage within the basin. By analyzing streamflow data from 39 U.S. Geological Survey basins, we show that is having a power law relationship with recession curve peak, which indicates that the spatial distribution of subsurface storage varies across recession events. Key Points The GRFM is reformulated for Horton-Strahler networks. The GRFM is modified by allowing its parameters to vary along streams. Sub-surface storage distribution controls recession flow characteristics.

Study of dynamic behaviour of recession curves

Since Brutsaert and Neiber (1977), recession curves are widely used to analyse subsurface systems of river basins by expressing -dQ/dt as a function of Q, which typically take a power law form: -dQ/dt=kQ, where Q is the discharge at a basin outlet at time t. Traditionally recession flows are modelled by single reservoir models that assume a unique relationship between -dQ/dt and Q for a basin. However, recent observations indicate that -dQ/dt-Q relationship of a basin varies greatly across recession events, indicating the limitation of such models. In this study, the dynamic relationship between -dQ/dt and Q of a basin is investigated through the geomorphological recession flow model which models recession flows by considering the temporal evolution of its active drainage network (the part of the stream network of the basin draining water at time t). Two primary factors responsible for the dynamic relationship are identified: (i) degree of aquifer recharge (ii) spatial variation of rainfall. Degree of aquifer recharge, which is likely to be controlled by (effective) rainfall patterns, influences the power law coefficient, k. It is found that k has correlation with past average streamflow, which confirms the notion that dynamic -dQ/dt-Q relationship is caused by the degree of aquifer recharge. Spatial variation of rainfall is found to have control on both the exponent, , and the power law coefficient, k. It is noticed that that even with same and k, recession curves can be different, possibly due to their different (recession) peak values. This may also happen due to spatial variation of rainfall. Copyright (c) 2012 John Wiley & Sons, Ltd.

What mainly controls recession flows in river basins?

The ubiquity of the power law relationship between dQ/dt and Q for recession periods (-dQ/dt kQ(alpha); Q being discharge at the basin outlet at time t) clearly hints at the existence of a dominant recession flow process that is common to all real basins. It is commonly assumed that a basin, during recession events, functions as a single phreatic aquifer resting on a impermeable horizontal bed or the Dupuit-Boussinesq (DB) aquifer, and with time different aquifer geometric conditions arise that give different values of alpha and k. The recently proposed alternative model, geomorphological recession flow model, however, suggests that recession flows are controlled primarily by the dynamics of the active drainage network (ADN). In this study we use data for several basins and compare the above two contrasting recession flow models in order to understand which of the above two factors dominates during recession periods in steep basins. Particularly, we do the comparison by selecting three key recession flow properties: (1) power law exponent alpha, (2) dynamic dQ/dt-Q relationship (characterized by k) and (3) recession timescale (time period for which a recession event lasts). Our observations suggest that neither drainage from phreatic aquifers nor evapotranspiration significantly controls recession flows. Results show that the value of a and recession timescale are not modeled well by DB aquifer model. However, the above mentioned three recession curve properties can be captured satisfactorily by considering the dynamics of the ADN as described by geomorphological recession flow model, possibly indicating that the ADN represents not just phreatic aquifers but the organization of various sub-surface storage systems within the basin. (C) 2014 Elsevier Ltd. All rights reserved.

`Universal' recession curves and their geomorphological interpretation

The study of recession flows offers fundamental insights into basin hydrological processes and, in particular, into the collective behavior of the governing dominant subsurface flows and properties. We use here an existing geomorphological interpretation of recession dynamics, which links the exponent in the classic recession curve -dQ/dt - kQ(alpha) to the geometric properties of the time-varying drainage network to study the general properties of recession curves across a wide variety of river basins. In particular, we show how the parameter k depends on the initial soil moisture state of the basin and can be made to explicitly depend on an index discharge, representative of initial sub-subsurface storage. Through this framework we obtain a non-dimensional, event-independent, recession curve. We subsequently quantify the variability of k across different basins on the basis of their geometry, and, by rescaling, collapse curves from different events and basins to obtain a generalized, or `universal', recession curve. Finally, we analyze the resulting normalized recession curves and explain their universal characteristics, lending further support to the notion that the statistical properties of observed recession curves bear the signature of the geomorphological structure of the networks producing them. (C) 2014 Elsevier Ltd. All rights reserved.

Relevo e drenagem do centro da cidade do Rio de Janeiro: enchentes e políticas públicas no passado e no presente

As enchentes se constituem em um dos mais freqüentes e complexos problemas nas áreas urbanas em todo o mundo, causando danos à população e às atividades econômicas a estas associadas. A cidade do Rio de Janeiro possui um longo histórico de inundações, outrora ligadas sobretudo às suas características topográficas e climáticas, que, em virtude dos processo de ocupação e urbanização, iniciados no século XVI, foram potencializas pelas interferências no meio físico, através das políticas públicas. Esta dissertação visa analisar as causas das inundações, bem como suas implicações na organização do espaço, na área central da cidade do Rio de Janeiro, evidenciando formas e processos pretéritos e presentes. As enchentes no centro da cidade foram mapeadas, com diferenciação entre os setores análogos, nos quais destacam-se as áreas mais críticas de ocorrências, em virtude dos danos decorrentes das inundações. Para tanto foram feitos levantamentos bibliográficos e trabalhos de campo. São feitas também recomendações que se originaram dos resultados dos levantamentos e das análises realizadas, visando contribuir para iniciativas que busquem soluções efetivas para o velho problema de enchentes na área central do Rio de Janeiro. Como resultado, pode-se concluir que uma série de fatores conjugados contribuem para as ocorrências atuais, àqueles inseridos em uma escala mais local, como nivelamentos de ruas, que ainda hoje remontam antigas feições da cidade ou a ineficiência da rede de drenagem atual, como também fatores mais abrangentes, ligados aos divisores topográficos e drenagem associada.

Limnology in Sicily: A history of fresh waters surrounded by the sea

The aim of this short article is to trace the history of limnology in Sicily, from the end of the last century up to the present, and pay a little homage to a scientist to whom limnology is deeply indebted: G.E. Hutchinson. Due to its insular and climatic conditions, Sicily is characterised by a drainage network formed by numerous short, torrent-like rivers, and by a few small, natural lakes. The geological characteristics of the island strongly condition the quality of these small waterbodies which generally have surface areas of less than 0.2 km super(2). The earliest observations on Sicilian lentic waters go back more than a century; in particular, at the end of the 19th century, some of the issues regarding the natural lakes had already been brought into focus. Subsequent studies on natural waterbodies concentrated on the geomorphology of landslide lakes or lakes created by the dissolution of the gypsum tableland. However, many of the waterbodies no longer exist because of land reclamation which took place up to the first half of the 1950s. During the last 35 years, there has been a notable increase in limnological publications. In addition, these studies show a more careful and integrated approach to the limnological aspects of waterbodies, compared with the early studies.

Modelling soil erosion and transport in the Burrishoole catchment, Newport, Co. Mayo, Ireland

The Burrishoole catchment is situated in County Mayo, on the northwest coast of the Republic of Ireland. Much of the catchment is covered by blanket peat that, in many areas, has become heavily eroded in recent years. This is thought to be due, primarily, to the adverse effects of forestry and agricultural activities in the area. Such activities include ploughing, drainage, the planting and harvesting of trees, and sheep farming, all of which are potentially damaging to such a sensitive landscape if not managed carefully. This article examines the sediment yield and hydrology of the Burrishoole catchment. Flow and sediment concentrations were measured at 8-hourly intervals from 5 February 2001 to 8 November 2001 with an automatic sampler and separate flow gauge, and hourly averages were recorded between 4 July 2002 and 6 September 2002 using an automatic river monitoring system [ARMS]. The authors describe the GIS-based model of soil erosion and transport that was applied to the Burrishoole catchment during this study. The results of these analyses were compared, in a qualitative manner, with the aerial photography available for the Burrishoole catchment to see whether areas that were predicted to contribute large proportions of eroded material to the drainage network corresponded with areas where peat erosion could be identified through photo-interpretation.

On the multiple ecological roles of water in river networks

The distribution and movement of water can influence the state and dynamics of terrestrial and aquatic ecosystems through a diversity of mechanisms. These mechanisms can be organized into three general categories wherein water acts as (1) a resource or habitat for biota, (2) a vector for connectivity and exchange of energy, materials, and organisms, and (3) as an agent of geomorphic change and disturbance. These latter two roles are highlighted in current models, which emphasize hydrologic connectivity and geomorphic change as determinants of the spatial and temporal distributions of species and processes in river systems. Water availability, on the other hand, has received less attention as a driver of ecological pattern, despite the prevalence of intermittent streams, and strong potential for environmental change to alter the spatial extent of drying in many regions. Here we summarize long-term research from a Sonoran Desert watershed to illustrate how spatial patterns of ecosystem structure and functioning reflect shifts in the relative importance of different 'roles of water' across scales of drainage size. These roles are distributed and interact hierarchically in the landscape, and for the bulk of the drainage network it is the duration of water availability that represents the primary determinant of ecological processes. Only for the largest catchments, with the most permanent flow regimes, do flood-associated disturbances and hydrologic exchange emerge as important drivers of local dynamics. While desert basins represent an extreme case, the diversity of mechanisms by which the availability and flow of water influence ecosystem structure and functioning are general. Predicting how river ecosystems may respond to future environmental pressures will require clear understanding of how changes in the spatial extent and relative overlap of these different roles of water shape ecological patterns. © 2013 Sponseller et al.

Identifying the influence of geology, land use, and anthropogenic activities on riverine sulfate on a watershed scale by combining hydrometric, chemical and isotopic approaches

Throughout the last few decades, sulfate concentrations in streamwater have received considerable attention due to their dominant role in anthropogenic acidification of surface waters. The objectives of this study conducted in the Oldman River Basin in Alberta (Canada) were to determine the influence of geology, land use and anthropogenic activities on sources, concentrations and fluxes of riverine sulfate on a watershed scale. This was achieved by combining hydrological, chemical and isotopic techniques. Surface water samples were collected from the main stem and tributaries of the Oldman River on a monthly basis between December 2000 and March 2003 and analyzed for chemical and isotopic compositions. At a given sampling site, sulfate sources were primarily dependent on geology and did not vary with time or flow condition. With increasing flow distance a gradual shift from ?34S values > 10 ‰ and ?18O values > 0 ‰ of riverine sulfate indicating evaporite dissolution and soil-derived sulfate in the predominantly forested headwaters, to negative ?34S and ?18O values suggested that sulfide oxidation was the predominant sulfate source in the agriculturally used downstream part of the watershed. Significant increases in sulfate concentrations and fluxes with downstream distance were observed, and were attributed to anthropogenically enhanced sulfide oxidation due to the presence of an extensive irrigation drainage network with seasonally varying water levels. Sulfate-S exports in an artificially drained subbasin (64 kg S/ha/yr) were found to exceed those in a naturally drained subbasin (4 kg S/ha/yr) by an order of magnitude. Our dataset suggests that the naturally occurring process of sulfide oxidation has been enhanced in the Oldman River Basin by the presence of an extensive network of drainage and irrigation canals.

Correlating Alpine glaciation with Adriatic sea-level changes through lake and alluvial stratigraphy

We compare lake and alluvial stratigraphy along a frame connecting the southern Alpine foothills and the Adriatic Sea, with the aim of matching the effects of Alpine glaciation and sea-level changes on sedimentation during the last glacial cycle. The palynostratigraphy of Lake Fimon provided proxies for regional vegetation and climate change and was coupled with sediment petrography, loss on ignition and magnetic susceptibility, disentangling alluvial phases from fluvioglacial activity related to culminations of the southeastern Alpine glaciers. The Fimon area was not reached by alluvial fans during the penultimate glacial maximum, nor by the sea transgression during the last interglacial, but a closed lake soon developed at the Eemian onset due to enhanced rainfall. Sea-level fall at glacial inception triggered the entrenchment of the drainage network in the plain reaching the outer Fimon Basin. Slow aggradation, but no sign of fluvioglacial activity, lasted to 38.2 +/- 1.45 cal. ka BP, when a major forest withdrawal took place, coeval to the spread of alluvial fans. By 27.5 perpendicular to 0.5 cal. ka BP the Fimon Basin was dammed by the Brenta outwash system. The main step of of forest recovery commenced at around (15.8) cal. ka BP, when apex trenching of the outwash fans was triggered by the glacier's decay. Copyright (C) 2011 John Wiley & Sons, Ltd.