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.

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.

A Network Model of Static Foam Drainage

On the basis of a more realistic tetrakaidecahedral structure of foam bubbles, a network model of static foam drainage has been developed. The model considers the foam to be made up of films and Plateau borders. The films drain into the adjacent Plateau borders, which in turn form a network through which the liquid moves from the foam to the liquid pool. From the structure, a unit flow cell was found, which constitutes the foam when stacked together both horizontally and vertically. Symmetry in the unit flow cell indicates that the flow analysis of a part of it can be employed to obtain the drainage for the whole foam. Material balance equations have been written for each segment of this subsection, ensuring connectivity, and solved with the appropriate boundary and initial conditions. The calculated rates of drainage, when compared with the available experimental results, indicate that the model predicts the experimental results well.

Ferrolysis induced soil transformation by natural drainage in Vertisols of sub-humid South India

In sub-humid South India, recent studies have shown that black soil areas (Vertisols and vertic Intergrades), located on flat valley bottoms, have been rejuvenated through the incision of streambeds, inducing changes in the pedoclimate and soil transformation. Joint pedological, geochemical and geophysical investigations were performed in order to better understand the ongoing processes and their contribution to the chemistry of local rivers. The seasonal rainfall causes cycles of oxidation and reduction in a perched watertable at the base of the black soil, while the reduced solutions are exported through a loamy sand network. This framework favours a ferrolysis process, which causes low base saturation and protonation of clay, leading to the weathering of 2:1 then 1:1 clay minerals. Maximum weathering conditions occur at the very end of the wet season, just before disappearance of the perched watertable. Therefore, the by-products of soil transformation are partially drained off and calcareous nodules, then further downslope, amorphous silica precipitate upon soil dehydration. The ferrolysed area is fringing the drainage system indicating that its development has been induced by the streambed incision. The distribution of C-14 ages of CaCO3 nodules suggests that the ferrolysis process started during the late Holocene, only about 2 kyr B.P. at the studied site and about 5 kyr B.P. at the watershed outlet. The results of this study are applied to an assessment of the physical erosion rate (4.8x10(-3) m/kyr) since the recent reactivation of the erosion process. (C) 2010 Elsevier B.V. All rights reserved.

Studies on removal of metal ions and sulphate reduction using rice husk and Desulfotomaculum nigrificans with reference to remediation of acid mine drainage

The utility of rice husk as an adsorbent for metal ions such as iron, zinc and copper from acid mine water was assessed. The adsorption isotherms exhibited Langmuirian behavior and were endothermic in nature. The free energy values for adsorption of the chosen metal ions onto rice husk were found to be highly negative attesting to favorable interaction. Over 99% Fe3+, 98% of Fe2+ and Zn2+ and 95% Cu2+ uptake was achieved from acid mine water, with a concomitant increase in the pH value by two units using rice husk. The remediation studies carried out on acid mine water and simulated acid mine water pretreated with rice husk indicated successful growth of Desulfotomaculum nigrificans (D. nigrificans). The amount of sulphate bioreduction in acid mine water at an initial pH of 5.3 was enhanced by D. nigrificans from 21% to 40% in the presence of rice husk filtrate supplemented with carbon and nitrogen. In simulated acid mine water with fortified husk filtrate, the sulphate reduction was even more extensive, with an enhancement to 73%. Concurrently, almost 90% Fe2+, 89% Zn2+ and 75% Cu2+ bioremoval was attained from simulated acid mine water. Metal adsorption by rice husk was confirmed in desorption experiments in which almost complete removal of metal ions from the rice husk was achieved after two elutions using 1 M HCl. The possible mechanisms of metal ion adsorption onto rice husk and sulphate reduction using D. nigrificans are discussed.

A model for static foam drainage

A model for static foam drainage, based on the pentagonal dodecahedral shape of bubbles, that takes into account the surface mobility of both films and Plateau border walls has been developed. The model divides the Plateau borders into nearly horizontal and nearly vertical categories and assigns different roles to them. The films are assumed to drain into all the adjacent Plateau borders equally. The horizontal Plateau borders are assumed to receive liquid from films and drain into vertical Plateau borders, which in turn form the main component for gravity drainage. The model yields the liquid holdup values for films, horizontal Plateau borders and vertical Plateau borders as functions of height and time. The model has been tested on static foams whose cumulative drainage was measured as a function of time. The experimental data on the effect of foam height, initial holdup, surface viscosity, etc. can be explained by the model quantitatively.

Drainage and separation factors for static foams containing agglomerates of microbial cells

The presence of cell agglomerates has been found to influence significantly the rates of liquid drainage from static foams. The process of drainage has been modelled by considering the foam to be made of pentagonal dodecahedral bubbles yielding films, nearly horizontal and nearly vertical Plateau borders. The films are assumed to drain into both kinds of Plateau borders equally. The horizontal Plateau borders are assumed to receive liquid from the films and drain into the vertical Plateau borders, which, in turn, form the main flow paths for gravity drainage. The drainage process is assumed to be similar to that for pure liquid until a stage is reached where the size of the cell agglomerates become equivalent to those of films and Plateau borders. Thereafter, a squeezing flow mechanism has been formulated where the aggromerates deform and flow. The model based on the above assumptions has been verified against experimental results and has been found to predict not only drainage data but also the separation of cell agglomerates from broths.

A new model for drainage of static foams

Existing theories of foam drainage assume bubbles as pentagonal dodecahedrons, though a close-packed structure built with cells of this shape is not space-filling. The present work develops a theory for calculating drainage rates based on the more realistic beta-tetrakaidecahedral shape for the bubbles. In contrast with the earlier works, three types of films, and Plateau borders had to be considered in view of the more complex shape used in the present work. The exchange of liquid between Plateau borders was treated in a way different From earlier theories, using the idea that the volume of junctions of Plateau borders is negligible. For foams made of large bubble sizes, the present model performs as well as the previous models, but when bubble size is small, its predictions of drainage rates from static foams are in better agreement with the experimental observations.

Consolidation behavior of clayey soils under radial drainage

For Barren's degree of consolidation, U-r, versus the time factor, T-r, relationship for soils undergoing consolidation with radial drainage for the equal vertical strain condition, a simple method has been developed to determine the value of the coefficient of consolidation with radial drainage c(r). Theoretical log(10)(d(e)(2)/t) versus U-r curves where d(e) is the diameter of influence and r is the real time for the different known value of c(r) have been generated. A method has been developed wherein both the theoretical and experimental behaviors of soils undergoing consolidation with radial drainage can be simultaneously compared and studied on the same plot. The experimental log(10)(d(e)(2)/t) versus U-r curves have been compared with the theoretical curves. Effects of initial compression, secondary compression, and duration of load increment are studied. Simple procedures are presented for calculating the values of c(r) using the experimental log(10)(d(e)(2)/t) versus U-r curves. A comparative study of the coefficient of consolidation and the coefficient of permeability between the cases of vertical and radial drainage has been done.

Influence of land use changes in river basins on diversity and distribution of amphibians

Land-use changes influence local biodiversity directly, and also cumulatively, contribute to regional and global changes in natural systems and quality of life. Consequent to these, direct impacts on the natural resources that support the health and integrity of living beings are evident in recent times. The Western Ghats being one of the global biodiversity hotspots, is reeling under a tremendous pressure from human induced changes in terms of developmental projects like hydel or thermal power plants, big dams, mining activities, unplanned agricultural practices,monoculture plantations, illegal timber logging, etc. This has led to the once contiguous forest habitats to be fragmented in patches, which in turn has led to the shrinkage of original habitat for the wildlife, change in the hydrological regime of the catchment, decreased inflow in streams,human-animal conflicts, etc. Under such circumstances, a proper management practice is called for requiring suitable biological indicators to show the impact of these changes, set priority regions and in developing models for conservation planning. Amphibians are regarded as one of the best biological indicators due to their sensitivity to even the slightest changes in the environment and hence they could be used as surrogates in conservation and management practices. They are the predominating vertebrates with a high degree of endemism (78%) in Western Ghats. The present study is an attempt to bring in the impacts of various land-uses on anuran distribution in three river basins. Sampling was carried out for amphibians during all seasons of 2003-2006 in basins of Sharavathi, Aghanashini and Bedthi. There are as many as 46 species in the region, one of which is new to science and nearly 59% of them are endemic to the Western Ghats. They belong to nine families, Dicroglossidae being represented by 14 species,followed by Rhacophoridae (9 species) and Ranidae (5 species). Species richness is high in Sharavathi river basin, with 36 species, followed by Bedthi 33 and Aghanashini 27. The impact of land-use changes, was investigated in the upper catchment of Sharavathi river basin. Species diversity indices, relative abundance values, percentage endemics gave clear indication of differences in each sub-catchment. Karl Pearson’s correlation coefficient (r) was calculated between species richness, endemics, environmental descriptors, land-use classes and fragmentation metrics. Principal component analysis was performed to depict the influence of these variables. Results show that sub-catchments with lesser percentage of forest, low canopy cover, higher amount of agricultural area, low rainfall have low species richness, less endemic species and abundant non-endemic species, whereas endemism, species richness and abundance of endemic species are more in the sub-catchments with high tree density, endemic trees, canopy cover, rainfall and lower amount of agriculture fields. This analysis aided in prioritising regions in the Sharavathi river basin for further conservation measures.

A comparative study of deglaciation in two neighbouring basins (Warwan and Bhut) of Western Himalaya

Glaciers of the Himalaya contribute significantly in the processes linking atmosphere, biosphere and hydrosphere, thus need to be monitored in view of the climatic variations. In this direction, many studies have been carried out during the last two decades and satellite-based multispectral data have been used extensively for this purpose throughout the world. The present study is aimed at mapping of glaciers in two adjacent basins (Warwan and Bhut) of the Western Himalaya with almost similar altitude and latitude and comparing the changes in the two time-frames with respect to three parameters, i.e. area, debris cover and area altitude distribution of glaciers. The two time-frames are topographical maps of 1962 and IRS LISS III images of 2001/02. Deglaciation was observed in both the basins with 19% and 9% loss in the glaciated area in Warwan and Bhut respectively. This difference may be due to: (i) the smaller size of the glaciers of the Warwan Basin (e.g. 164 glaciers having <1 sq. km area in comparison to 101 glaciers in the Bhut Basin), (ii) lower percentage of moraine cover in Warwan (18) than in the Bhut Basin (30) and (iii) higher percentage of glaciated area lying below 5100 m (80) in Warwan than in the Bhut Basin (70).

Variation of Snowline and Mass Balance of Glaciers of Warwan and Bhut Basins of Western Himalaya Using Remote Sensing Technique

Glaciers are natural reservoirs of fresh water in frozen state and sensitive indicators of climate change. Among all the mountainous glaciated regions, glaciers of Himalayas form one of the largest concentrations of ice outside the Polar Regions. Almost all the major rivers of northern India originate from these glaciers and sustain perennial flow. Therefore, in view of the importance and role of the glaciers in sustaining the life on the Earth, monitoring the health of glaciers is necessary. Glacier's health is monitored in two ways (i) by mapping the change in extent of glaciers (ii) by finding variation in the annual mass balance. This paper has been discussed the later approach for monitoring the health of glaciers of Warwan and Bhut basins. Mass balance of glaciers of these two basins was determined based on the extraction of snow line at the end of ablation season. A series of satellite images of AWiFS sensor were analysed for extraction of snowline on the glaciers for the period of 2005, 2006 and 2007. The snow line at the end of ablation season is used to compute accumulation area ratio (AAR = Accumulation area/Glacier area) for each glacier of basins. An approach based on relationship of AAR to specific mass balance (computed in field) for glaciers of Basapa basin was employed in the present study. Mean of specific mass balance of individual glacier for the year 2005, 2006 and 2007 of Warwan basin was found to be -ve 0.19 m, -ve 0.27 m and -ve 0.2 m respectively. It is 0.05 m, -ve 0.11 m and -ve 0.19 m for Bhut basin. The analysis suggests a loss of 4.3 and 0.83 kmA(3) of glacier in the monitoring period of 3 years for Warwan and Bhut basins respectively. The overall results suggest that the glaciers of Warwan basin and Bhut basins have suffered more loss of ice than gain in the monitoring period of 3 years.

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.

`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.

Estimation of `drainable' storage - A geomorphological approach

Storage of water within a river basin is often estimated by analyzing recession flow curves as it cannot be `instantly' estimated with the aid of available technologies. In this study we explicitly deal with the issue of estimation of `drainable' storage, which is equal to the area under the `complete' recession flow curve (i.e. a discharge vs. time curve where discharge continuously decreases till it approaches zero). But a major challenge in this regard is that recession curves are rarely `complete' due to short inter-storm time intervals. Therefore, it is essential to analyze and model recession flows meaningfully. We adopt the wellknown Brutsaert and Nieber analytical method that expresses time derivative of discharge (dQ/dt) as a power law function of Q : -dQ/dt = kQ(alpha). However, the problem with dQ/dt-Q analysis is that it is not suitable for late recession flows. Traditional studies often compute alpha considering early recession flows and assume that its value is constant for the whole recession event. But this approach gives unrealistic results when alpha >= 2, a common case. We address this issue here by using the recently proposed geomorphological recession flow model (GRFM) that exploits the dynamics of active drainage networks. According to the model, alpha is close to 2 for early recession flows and 0 for late recession flows. We then derive a simple expression for drainable storage in terms the power law coefficient k, obtained by considering early recession flows only, and basin area. Using 121 complete recession curves from 27 USGS basins we show that predicted drainable storage matches well with observed drainable storage, indicating that the model can also reliably estimate drainable storage for `incomplete' recession events to address many challenges related to water resources. (C) 2014 Elsevier Ltd. All rights reserved.