Estimating percolation and lateral water flow on sloping land in rainfed lowland rice ecosystem

Quantifying water losses in paddy fields assists estimation of water availability in rainfed lowland rice ecosystem. Little information is available on water balance in different toposequence positions of sloped rainfed lowland. Therefore, the aim of this work was to quantify percolation and the lateral water flow with special reference to the toposequential variation. Data used for the analysis was collected in Laos and northeast Thailand. Percolation and water tables were measured on a daily basis using a steel cylindrical tube with a lid and perforated PVC tubes, respectively. Percolation rate was determined using linear regression analysis of cumulative percolation. Assuming that the total amount of evaporation and transpiration was equivalent to potential evapotranspiration, the lateral water flow was estimated using the water balance equation. Separate perched water and groundwater tables were observed in paddy fields on coarse-textured soils. The percolation rate varied between 0 and 3 mm/day across locations, and the maximum water loss by lateral movement was more than 20 mm/day. Our results are in agreement with the previously reported findings, and the methodology of estimating water balance components appears reasonably acceptable. With regard to the toposequential variation, the higher the position in the topoesquence, the greater potential for water loss because of higher percolation and lateral flow rates.

Application of a coupled ground-surface water flow model to simulate periodic groundwater flow influenced by a sloping boundary, capillarity and vertical flows

Theoretical developments as well as field and laboratory data have shown the influence of the capillary fringe on water table fluctuations to increase with the fluctuation frequency. The numerical solution of a full, partially saturated flow equation can be computationally expensive. In this paper, the influence of the capillary fringe on water table fluctuations is simplified through its parameterisation into the storage coefficient of a fully-saturated groundwater flow model using the complex effective porosity concept [Nielsen, P., Perrochet, P., 2000. Water table dynamics under capillary fringes: experiments and modelling. Advances in Water Resources 23 (1), 503-515; Nielsen, P., Perrochet, P., 2000. ERRATA: water table dynamics under capillary fringes: experiments and modelling (Advances in Water Resources 23 (2000) 503-515). Advances in Water Resources 23, 907-908]. The model is applied to sand flume observations of periodic water table fluctuations induced by simple harmonic forcing across a sloping boundary, analogous to many beach groundwater systems. While not providing information on the moisture distribution within the aquifer, this approach can reasonably predict the water table fluctuations in response to periodic forcing across a sloping boundary. Furthermore, he coupled ground-surface water model accurately predicts the extent of the seepage face formed at the sloping boundary. (C) 2005 Elsevier Ltd. All rights reserved.

A new analytical solution for water table fluctuations in coastal aquifers with sloping beaches

The Boussinesq equation appears as the zeroth-order term in the shallow water flow expansion of the non-linear equation describing the flow of fluid in an unconfined aquifer. One-dimensional models based on the Boussinesq equation have been used to analyse tide-induced water table fluctuations in coastal aquifers. Previous analytical solutions for a sloping beach are based on the perturbation parameter, epsilon(N) = alphaepsilon cot beta (in which beta is the beach slope, alpha is the amplitude parameter and epsilon is the shallow water parameter) and are limited to tan(-1) (alphaepsilon) much less than beta less than or equal to pi/2. In this paper, a new higher-order solution to the non-linear boundary value problem is derived. The results demonstrate the significant influence of the higher-order components and beach slope on the water table fluctuations. The relative difference between the linear solution and the present solution increases as 6 and a increase, and reaches 7% of the linear solution. (C) 2003 Elsevier Ltd. All rights reserved.

Effect of water depth and water velocity upon the surfacing frequency of the bimodally respiring freshwater turtle, Rheodytes leukops

This study examines the effect of increasing water depth and water velocity upon the surfacing behaviour of the bimodally respiring turtle, Rheodytes leukops. Surfacing frequency was recorded for R. leukops at varying water depths (50, 100, 150 cm) and water velocities (5, 15, 30 cm s(-1)) during independent trials to provide an indirect cost-benefit analysis of aquatic versus pulmonary respiration. With increasing water velocity, R. leukops decreased its surfacing frequency twentyfold, thus suggesting a heightened reliance upon aquatic gas exchange. An elevated reliance upon aquatic respiration, which presumably translates into a decreased air-breathing frequency, may be metabolically more efficient for R. leukops compared to the expenditure (i.e. time and energy) associated with air-breathing within fast-flowing riffle zones. Additionally, R. leukops at higher water velocities preferentially selected low-velocity microhabitats, presumably to avoid the metabolic expenditure associated with high water flow. Alternatively, increasing water depth had no effect upon the surfacing frequency of R. leukops, suggesting little to no change in the respiratory partitioning of the species across treatment settings. Routinely long dives (>90 min) recorded for R. leukops indicate a high reliance upon aquatic O-2 uptake regardless of water depth. Moreover, metabolic and temporal costs attributed to pulmonary gas exchange within a pool-like environment were likely minimal for R. leukops, irrespective of water depth.

Effect of particle hydrophobicity on particle and water transport across a flotation froth

Froth recovery measurements have been conducted in both the presence (three-phase froth) and absence (two-phase froth) of particles of different contact angles in a specially modified laboratory flotation column. Increasing the particle hydrophobicity increased the flow rate of particles entering the froth, while the recovery of particles across the froth phase itself also increased for particle contact angles to 63 and at all vertical heights of the froth column. However, a further increase in the contact angle to 69 resulted in lower particle recovery across the froth phase. The reduced froth recovery for particles of 69 contact angle was linked to significant bubble coalescence within the froth phase. The reduced froth recovery occurred uniformly across the entire particle size range, and was, presumably, a result of particle detachment from coalescing bubbles. Water flow rates across the froth phase also varied with particle contact angle. The general trend was a decrease in the concentrate flow rate of water with increasing particle contact angle. An inverse relationship between water flow rate and bubble radius was also observed, possibly allowing prediction of water flow rate from bubble size measurements in the froth. Comparison of the froth structure, defined by bubble size, gas hold-up and bubble layer thickness, for two- and three-phase froths, at the same frother concentration, showed there was a relationship between water flow rate and froth structure. (c) 2005 Elsevier B.V. All rights reserved.

The water balance on sloping land in rainfed lowland rice ecosystem

Quantifying water losses in paddy fields assists estimation of water availability in rainfed lowland rice ecosystem. There is, however, no definite method for determining the water losses, and little information is available on water balance in different toposequence positions of a sloped rainfed lowland. Therefore, the aim of this work was to quantify percolation and the lateral water flow with special reference to the toposequential variation. Data used for the analysis was collected in Laos and northeast Thailand. Percolation and water tables were measured on a daily basis. The percolator is a steel cylindrical tube with a lid to prevent water loss from evapotranspiration. The water table meter is a short PVC tube for determination of perched water table and a long PVC tube for groundwater table, and the side is perforated with 5-mm diameter holes at 20-mm distance. Percolation rate was determined using linear regression analysis of cumulative percolation. Assuming that the total amount of evaporation and transpiration was equivalent to potential evapotranspiration, the lateral water flow was estimated using the water balance equation. Our results are in agreement with the previously reported findings, and the methodology of estimating water balance components appears reasonably acceptable. With regard to the toposequential variation, the higher position in the topoesquence, the greater potential of the water losses because of higher percolation and lateral flow rates.

On sensor network segmentation for urban water distribution monitoring

Identifying water wastage in forms of leaks in a water distribution network of any city becomes essential as droughts are presenting serious threats to few major cities. In this paper, we propose a deployment of sensor network for monitoring water flow in any water distribution network. We cover the issues related with designing such a dedicated sensor network by considering types of sensors required, sensors' functionality, data collection, and providing computation serving as leak detection mechanism. The main focus of this paper is on appropriate network segmentation that provides the base for hierarchical approach to pipes' failure detection. We show a method for sensors allocation to the network in order to facilitate effective pipes monitoring. In general, the identified computational problem belongs to hard problems. The paper shows a heuristic method to build effective hierarchy of the network segmentation.

SimuMap: a computational system for spatial modelling

The paper presents a computational system based upon formal principles to run spatial models for environmental processes. The simulator is named SimuMap because it is typically used to simulate spatial processes over a mapped representation of terrain. A model is formally represented in SimuMap as a set of coupled sub-models. The paper considers the situation where spatial processes operate at different time levels, but are still integrated. An example of such a situation commonly occurs in watershed hydrology where overland flow and stream channel flow have very different flow rates but are highly related as they are subject to the same terrain runoff processes. SimuMap is able to run a network of sub-models that express different time-space derivatives for water flow processes. Sub-models may be coded generically with a map algebra programming language that uses a surface data model. To address the problem of differing time levels in simulation, the paper: (i) reviews general approaches for numerical solvers, (ii) considers the constraints that need to be enforced to use more adaptive time steps in discrete time specified simulations, and (iii) scaling transfer rates in equations that use different time bases for time-space derivatives. A multistep scheme is proposed for SimuMap. This is presented along with a description of its visual programming interface, its modelling formalisms and future plans. (C) 2003 Elsevier Ltd. All rights reserved.

Numerical investigation of seawater intrusion at Gooburrum, Bundaberg, Queensland, Australia

Seawater intrusion in coastal agricultural areas due to groundwater abstraction is a major environmental problem along the northeastern coast of Australia. Management options are being explored using numerical modelling, however, questions remain concerning the appropriate level of sophistication in models, choice of seaward boundary conditions, and how to accommodate heterogeneity and data uncertainty. The choice of seaward boundary condition is important since it affects the amount of salt transported into the aquifers and forms the focus of the present study. The impact of this boundary condition is illustrated for the seawater-intrusion problem in the Gooburrum aquifers, which occur within Tertiary sedimentary strata. A two-dimensional variable-density groundwater and solute-transport model was constructed using the computer code 2DFEMFAT (Cheng et al. 1998). The code was tested against an experiment for a steady-state freshwater-saltwater interface and against the Elder (Elder 1967) free-convection problem. Numerical simulations show that the imposition of the commonly-used equivalent hydrostatic freshwater heads, combined with a constant salt concentration at the seaward boundary, results in overestimated seawater intrusion in the lower Gooburrum aquifer. Since the imposition of this boundary condition allows water flow across the boundary, which subsequently takes salt into the aquifer, a careful check is essential to estimate whether too much mass of salt is introduced.

Chloride as a signature indicator of soil textural and hydrologic stratigraphies in variable charge deep profiles

Soil properties that influence water movement through profiles are important for determining flow paths, reactions between soil and solute, and the ultimate destination of solutes. This is particularly important in high rainfall environments. For highly weathered deep profiles, we hypothesize that abrupt changes in the distribution of the quotient [QT = (silt + sand)/clay] reflect the boundaries between textural units or textural (TS) and hydrologic (HS) stratigraphies. As a result, QT can be used as a parameter to characterize TS and as a surrogate for HS. Secondly, we propose that if chloride distributions were correlated with QT, under non-limiting anion exchange, then chloride distributions can be used as a signature indicator of TS and HS. Soil cores to a depth of 12.5 in were taken from 16 locations in the wet tropical Johnstone River catchment of northeast Queensland, Australia. The cores belong to nine variable charge soil types and were under sugarcane (Saccharun officinarum-S) production, which included the use of potassium chloride, for several decades. The cores were segmented at I m depth increments and subsamples were analysed for chloride, pH, soil water content (theta), clay, silt and sand contents. Selected bores were capped to serve as piezometers to monitor groundwater dynamics. Depth incremented QT, theta and chloride correlated, each individually, significantly with the corresponding profile depth increments, indicating the presence of textural, hydrologic and chloride gradients in profiles. However, rapid increases in QT down the profile indicated abrupt changes in TS, suggesting that QT can be used as a parameter to characterize TS and as a surrogate for HS. Abrupt changes in chloride distributions were similar to QT, suggesting that chloride distributions can be used as a signature indicator of QT (TS) and HS. Groundwater data indicated that chloride distributions depended, at least partially, on groundwater dynamics, providing further support to our hypothesis that chloride distribution can be used as a signature indicator of HS. Copyright (c) 2005 John Wiley & Sons, Ltd.

Effect of rainfall as a component of climate change on estuarine fish production in Queensland, Australia

The speculation that climate change may impact on sustainable fish production suggests a need to understand how these effects influence fish catch on a broad scale. With a gross annual value of A$ 2.2 billion, the fishing industry is a significant primary industry in Australia. Many commercially important fish species use estuarine habitats such as mangroves, tidal flats and seagrass beds as nurseries or breeding grounds and have lifecycles correlated to rainfall and temperature patterns. Correlation of catches of mullet (e.g. Mugil cephalus) and barramundi (Lates calcarifer) with rainfall suggests that fisheries may be sensitive to effects of climate change. This work reviews key commercial fish and crustacean species and their link to estuaries and climate parameters. A conceptual model demonstrates ecological and biophysical links of estuarine habitats that influences capture fisheries production. The difficulty involved in explaining the effect of climate change on fisheries arising from the lack of ecological knowledge may be overcome by relating climate parameters with long-term fish catch data. Catch per unit effort (CPUE), rainfall, the Southern Oscillation Index (SOI) and catch time series for specific combinations of climate seasons and regions have been explored and surplus production models applied to Queensland's commercial fish catch data with the program CLIMPROD. Results indicate that up to 30% of Queensland's total fish catch and up to 80% of the barramundi catch variation for specific regions can be explained by rainfall often with a lagged response to rainfall events. Our approach allows an evaluation of the economic consequences of climate parameters on estuarine fisheries. thus highlighting the need to develop forecast models and manage estuaries for future climate chan e impact by adjusting the quota for climate change sensitive species. Different modelling approaches are discussed with respect to their forecast ability. (c) 2006 Elsevier Ltd. All rights reserved.

Modelling trickle irrigation: Comparison of analytical and numerical models for estimation of wetting front position with time

Irrigation practices that are profligate in their use of water have come under closer scrutiny by water managers and the public. Trickle irrigation has the propensity to increase water use efficiency but only if the system is designed to meet the soil and plant conditions. Recently we have provided a software tool, WetUp (http://www.clw.csiro.au/products/wetup/), to calculate the wetting patterns from trickle irrigation emitters. WetUp uses an analytical solution to calculate the wetted perimeter for both buried and surface emitters. This analytical solution has a number of assumptions, two of which are that the wetting front is defined by water content at which the hydraulic conductivity (K) is I mm day(-1) and that the flow occurs from a point source. Here we compare the wetting patterns calculated with a 2-dimensional numerical model, HYDRUS2D, for solving the water flow into typical soils with the analytical solution. The results show that the wetting patterns are similar, except when the soil properties result in the assumption of a point source no longer being a good description of the flow regime. Difficulties were also experienced with getting stable solutions with HYDRUS2D for soils with low hydraulic conductivities. (c) 2005 Elsevier Ltd. All rights reserved.

Incorporating level set methods in Geographical Information Systems (GIS) for land-surface process modelling

Land-surface processes include a broad class of models that operate at a landscape scale. Current modelling approaches tend to be specialised towards one type of process, yet it is the interaction of processes that is increasing seen as important to obtain a more integrated approach to land management. This paper presents a technique and a tool that may be applied generically to landscape processes. The technique tracks moving interfaces across landscapes for processes such as water flow, biochemical diffusion, and plant dispersal. Its theoretical development applies a Lagrangian approach to motion over a Eulerian grid space by tracking quantities across a landscape as an evolving front. An algorithm for this technique, called level set method, is implemented in a geographical information system (GIS). It fits with a field data model in GIS and is implemented as operators in map algebra. The paper describes an implementation of the level set methods in a map algebra programming language, called MapScript, and gives example program scripts for applications in ecology and hydrology.