Runoff sources and land cover change in the Amazon: an end-member mixing analysis from small watersheds

The flowpaths by which water moves from watersheds to streams has important consequences for the runoff dynamics and biogeochemistry of surface waters in the Amazon Basin. The clearing of Amazon forest to cattle pasture has the potential to change runoff sources to streams by shifting runoff to more surficial flow pathways. We applied end-member mixing analysis (EMMA) to 10 small watersheds throughout the Amazon in which solute composition of streamwater and groundwater, overland flow, soil solution, throughfall and rainwater were measured, largely as part of the Large-Scale Biosphere-Atmosphere Experiment in Amazonia. We found a range in the extent to which streamwater samples fell within the mixing space determined by potential flowpath end-members, suggesting that some water sources to streams were not sampled. The contribution of overland flow as a source of stream flow was greater in pasture watersheds than in forest watersheds of comparable size. Increases in overland flow contribution to pasture streams ranged in some cases from 0% in forest to 27-28% in pasture and were broadly consistent with results from hydrometric sampling of Amazon forest and pasture watersheds that indicate 17- to 18-fold increase in the overland flow contribution to stream flow in pastures. In forest, overland flow was an important contribution to stream flow (45-57%) in ephemeral streams where flows were dominated by stormflow. Overland flow contribution to stream flow decreased in importance with increasing watershed area, from 21 to 57% in forest and 60-89% in pasture watersheds of less than 10 ha to 0% in forest and 27-28% in pastures in watersheds greater than 100 ha. Soil solution contributions to stream flow were similar across watershed area and groundwater inputs generally increased in proportion to decreases in overland flow. Application of EMMA across multiple watersheds indicated patterns across gradients of stream size and land cover that were consistent with patterns determined by detailed hydrometric sampling.

Collapse behavior of soil in a Brazilian region affected by a rising water table

Collapsible soils are usually nonsaturated, low density, and metastable-structured soils that are known to exhibit a volume reduction following an episode of moisture increase or suction reduction. This paper describes the collapsible behavior of clayey sand based on controlled soil suction tests carried out on undisturbed samples from the city of Pereira Barreto, in the State of Sao Paulo, Brazil. Foundation settlements due to soil collapse are common in this region and occurred during the filling of the reservoir of the Tres Irmaos Dam, which induced the elevation of the groundwater table in different parts of Pereira Barreto. This paper shows that collapse strains depend on the stress and soil suction acting in the sample and that saturation is not necessary for a collapse to occur. The influence of soil suction, gradual wetting, and the wetting and drying cycle on the collapsible behavior of the soil is also shown and discussed.

A comparison of dynamic and quasi-static results for the flow towards a well

The present analysis takes into account the acceleration term in the differential equation of motion to obtain exact dynamic solutions concerning the groundwater flow towards a well in a confined aquifer. The results show that the error contained in the traditional quasi-static solution is very small in typical situations.

A novel continuous time representation for the scheduling of pipeline systems with pumping yield rate constraints

Pipeline systems play a key role in the petroleum business. These operational systems provide connection between ports and/or oil fields and refineries (upstream), as well as between these and consumer markets (downstream). The purpose of this work is to propose a novel MINLP formulation based on a continuous time representation for the scheduling of multiproduct pipeline systems that must supply multiple consumer markets. Moreover, it also considers that the pipeline operates intermittently and that the pumping costs depend on the booster stations yield rates, which in turn may generate different flow rates. The proposed continuous time representation is compared with a previously developed discrete time representation [Rejowski, R., Jr., & Pinto, J. M. (2004). Efficient MILP formulations and valid cuts for multiproduct pipeline scheduling. Computers and Chemical Engineering, 28, 1511] in terms of solution quality and computational performance. The influence of the number of time intervals that represents the transfer operation is studied and several configurations for the booster stations are tested. Finally, the proposed formulation is applied to a larger case, in which several booster configurations with different numbers of stages are tested. (C) 2007 Elsevier Ltd. All rights reserved.

Carbon and nitrogen cycling in a tropical Brazilian soil cropped with sugarcane and irrigated with wastewater

Carbon (C) and nitrogen (N) dynamics in agro-systems can be altered as a consequence of treated sewage effluent (TSE) irrigation. The present study evaluated the effects of TSE irrigation over 16 months on N concentrations in sugarcane (leaves, stalks and juice), total soil carbon (TC), total soil nitrogen (TN), NO(3)(-)-N in soil and nitrate (NO(3)(-)) and dissolved organic carbon (DOC) in soil solution. The soil was classified as an Oxisol and samplings were carried out during the first productive crop cycle, from February 2005 (before planting) to September 2006 (after sugarcane harvest and 16 months of TSE irrigation). The experiment was arranged in a complete block design with five treatments and four replicates. Irrigated plots received 50% of the recommended mineral N fertilization and 100% (T100), 125% (T125), 150% (T150) and 200% (T200) of crop water demand. No mineral N and irrigation were applied to the control plots. TSE irrigation enhanced sugarcane yield but resulted in total-N inputs(804-1622 kg N ha(-1)) greater than exported N (463-597 kg N ha(-1)). Hence, throughout the irrigation period, high NO(3)(-) concentrations (up to 388 mg L(-1) at T200) and DOC (up to 142 mg L(-1) at T100) were measured in soil solution below the root zone, indicating the potential of groundwater contamination. TSE irrigation did not change soil TC and TN. (C) 2009 Elsevier B.V. All rights reserved.

Effectiveness of two systems for the spraying of pesticides in citrus trees

Citriculture normally uses high application volumes in pesticide solutions (of 2.000 to 5.000 L ha(-1)) to control pests and diseases that affect the crop, which generates an increase in operational costs. For this reason, diverse systems of application are being developed to reduce application volumes and improve the uniformity of pesticide deposition. The goal of this work was to evaluate the efficiency of two application systems of pesticides in citrus trees. One system used a prototype for terrestrial application with rotary disc atomizers that are widely used in agricultural aviation, and the other system used hollow cone tip hydraulics. For the treatment of the trees the insecticide Metidation was used at the dose of 180 gr per hectare. To study the droplet spectrum, water-sensitive papers were installed at different positions in the trees canopy, and for the study of insecticide deposition leaves of the treated plants were collected. The water-sensitive papers were collected and analyzed using a computerized image analysis system (e-Sprinkle, EMBRAPA, Sao Paulo, Brazil), and the leaves analyzed by the technique of gas chromatography. Pesticide deposition was similar in both application system, although the solution volume used by the application system equipped with rotary disc atomizers was one quarter of the volume used by the application system equipped with hydraulic tips, reducing considerably the cost of the phytosanitary treatments.

Characterization of the microtube emitters used in a novel micro-sprinkler

The microtube is a simple and cheap emitter that was widely used throughout the world in the early days of drip irrigation. Its length can be adjusted according to the pressure distribution along the lateral line and the discharge from the microtube can be adjusted by its length. This not only counters the pressure loss due to pipe friction but also makes it suitable for undulating and hilly conditions, where pressure in the lateral line varies considerably according to the differences in elevation. This is the major problem facing the designer, i.e., emitter flow changes as the acting pressure head changes. In this study, a novel micro-sprinkler system is proposed that uses microtube as the emitter and where the length of the microtube can be varied in response to pressure changes along the lateral to give uniformity of emitter discharges. The objective of this work is to develop and validate empirical and semi-theoretical equations for the emitter hydraulics. Laboratory testing of two microtube emitters of different diameter over a range of pressures and discharges was used in the development of the equations relating pressure and discharge, and pressure and length for these emitters. The equations proposed will be used in the design of the micro-sprinkler system, to determine the length of microtube required to give the nominal discharge for any given pressure. The semi-theoretical approach underlined the importance of accurate measurements of the microtube diameter and the uncertainty in the estimation of the friction factor for these tubes.

Podzolization as a deferralitization process: dynamics and chemistry of ground and surface waters in an Acrisol-Podzol sequence of the upper Amazon Basin

Hydrochemical processes involved in the development of hydromorphic Podzols are a major concern for the upper Amazon Basin because of the extent of the areas affected by such processes and the large amounts of organic carbon and associated metals exported to the rivers. The dynamics and chemical composition of ground and surface waters were studied along an Acrisol-Podzol sequence lying in an open depression of a plateau. Water levels were monitored along the sequence over a period of 2 years by means of piezometers. Water was sampled in zero-tension lysimeters for groundwater and for surface water in the drainage network of the depression. The pH and concentrations of organic carbon and major elements (Si, Fe and Al) were determined. The contrasted changes reported for concentrations of Si, organic carbon and metals (Fe, Al) mainly reflect the dynamics of the groundwater and the weathering conditions that prevail in the soils. Iron is released by the reductive dissolution of Fe oxides, mostly in the Bg horizons of the upslope Acrisols. It moves laterally under the control of hydraulic gradients and migrates through the iron-depleted Podzols where it is exported to the river network. Aluminium is released from the dissolution of Al-bearing minerals (gibbsite and kaolinite) at the margin of the podzolic area but is immobilized as organo-Al complexes in spodic horizons. In downslope positions, the quick recharge of the groundwater and large release of organic compounds lead to acidification and a loss of metals (mainly Al), previously stored in the Podzols.

Experimental Analysis of Froude Number Effect on Air Entrainment in the Hydraulic Jump

A hydraulic jump is characterized by strong energy dissipation and mixing, large-scale turbulence, air entrainment, waves and spray. Despite recent pertinent studies, the interaction between air bubbles diffusion and momentum transfer is not completely understood. The objective of this paper is to present experimental results from new measurements performed in rectangular horizontal flume with partially-developed inflow conditions. The vertical distributions of void fraction and air bubbles count rate were recorded for inflow Froude number Fr1 in the range from 5.2 to 14.3. Rapid detrainment process was observed near the jump toe, whereas the structure of the air diffusion layer was clearly observed over longer distances. These new data were compared with previous data generally collected at lower Froude numbers. The comparison demonstrated that, at a fixed distance from the jump toe, the maximum void fraction Cmax increases with the increasing Fr1. The vertical locations of the maximum void fraction and bubble count rate were consistent with previous studies. Finally, an empirical correlation between the upper boundary of the air diffusion layer and the distance from the impingement point was provided.

Geochemical mass-balance and oxygen-isotope constrains on silcrete formation and its paleoclimatic implications in Southern Australia

Two geographically distinct silcrete associations are present in southern Australia, inland and eastern; these were sampled in central South Australia and central Victoria, respectively, At each site, both silicified and immediately adjacent unsilicified parent material were collected. Analytical data from these pairs were used to construct isocons, assuming Zr immobility, and to calculate the volume change and amount of silica introduced during silicification, These results, together with whole-rock oxygen isotope compositions, were used to determine the delta(18)O of th, introduced silica, The results show that the eastern silcretes in central Victoria are probably linked genetically to the associated basalts, weathering of which supplied the introduced silica, This conclusion is based on the close spatial connection between the two, as well as the substantial amount of introduced silica in the silcretes (greater than in the inland silcretes), resulting in volume increases in some eastern silcretes, Oxygen isotopic calculations for the silcretes indicate that the silica precipitated from groundwaters at temperatures slightly higher than present conditions. Silcrete formation apparently occurred during the Miocene and Pliocene (basalts in Victoria younger than Pliocene lack associated silcrete) and may reflect the much wetter climate in southeastern Australia at that time. The inland silcretes of central South Australia can be divided into pedogenic (the most common) and groundwater varieties. The pedogenic silcretes, which show typical soil features like columnar and nodular textures, contain moderate amounts of introduced silica that precipitated by evaporation from saline groundwaters, For the groundwater silcretes, which have massive textures and formed at or close to the water table, insufficient data are available to determine the mode of formation. The inland pedogenic silcretes have probably been farming from the Eocene-Miocene to the present, implying that conditions of seasonally high evaporation have occurred in central Australia during this time period. Thus silcrete formation depends on a complex interplay between climate and silica supply, and it is impossible to generalize that the presence of silcrete is indicative of a particular climate. Likewise, the elemental composition of silcretes, particularly Ti content, is not necessarily of climatic significance, Nevertheless, detailed geochemical and oxygen isotopic studies of a silcrete and its parent material can elucidate the mechanisms of silcrete formation, and if evaporation is indicated as a major factor in silcrete formation, then the climate at the time was likely to have been at least seasonally arid.

Finite element modelling of temperature gradient driven rock alteration and mineralization in porous rock masses

We present finite element simulations of temperature gradient driven rock alteration and mineralization in fluid saturated porous rock masses. In particular, we explore the significance of production/annihilation terms in the mass balance equations and the dependence of the spatial patterns of rock alteration upon the ratio of the roll over time of large scale convection cells to the relaxation time of the chemical reactions. Special concepts such as the gradient reaction criterion or rock alteration index (RAI) are discussed in light of the present, more general theory. In order to validate the finite element simulation, we derive an analytical solution for the rock alteration index of a benchmark problem on a two-dimensional rectangular domain. Since the geometry and boundary conditions of the benchmark problem can be easily and exactly modelled, the analytical solution is also useful for validating other numerical methods, such as the finite difference method and the boundary element method, when they are used to dear with this kind of problem. Finally, the potential of the theory is illustrated by means of finite element studies related to coupled flow problems in materially homogeneous and inhomogeneous porous rock masses. (C) 1998 Elsevier Science S.A. All rights reserved.

An investigation of localised soil heterogeneities on solute transport using a multisegment percolation system

A multisegment percolation system (MSPS) consisting of 25 individual collection wells was constructed to study the effects of localised soil heterogeneities on the transport of solutes in the vadose zone. In particular, this paper discusses the transport of water and nutrients (NO3-, Cl-, PO43-) through structurally stable, free-draining agricultural soil from Victoria, Australia. A solution of nutrients was irrigated onto the surface of a large undisturbed soil core over a 12-h period. This was followed by a continuous irrigation of distilled water at a fate which did not cause pending for a further 18 days. During this time, the volume of leachate and the concentration of nutrients in the leachate of each well were measured. Very significant variation in drainage patterns across a small spatial scale was observed. Leaching of nitrate-nitrogen and chloride from the core occurred two days after initial application. However, less than 1% of the total applied phosphate-phosphorus leached from the soil during the 18-day experiment, indicating strong adsorption. Our experiments indicate considerable heterogeneity in water flow patterns and solute leaching on a small spatial scale. These results have significant ramifications for modelling solute transport and predicting nutrient loadings on a larger scale.

Tidal effects on sea water intrusion in unconfined aquifers

A variable-density groundwater model is used to analyse the effects of tidal fluctuations on sea-water intrusion in an unconfined aquifer. It is shown that the tidal activity forces the sea-water to intrude further inland and it also creates a thicker interface than would occur without tidal effects. Moreover, the configuration of the interface is radically changed when the tidal fluctuations are included. This is because of the dramatic changes in the flow pattern and velocity of the groundwater near the shoreline. For aquifer depths much larger than tidal amplitudes, the tidal fluctuation does not have much effect on how far the sea-water intrudes into the aquifer; nevertheless, a significant change in the configuration of concentration contours because of the effect of tidal fluctuations is observed. This change is more noticeable at the top of the aquifer, near the water table, than at the bottom of the aquifer, and is caused by the infiltration of salt water into the top of the aquifer at higher tidal levels. A flatter beach slope, therefore, intensifies this phenomenon. The interface configurations do not change noticeably over the course of a tidal cycle. Neglecting tidal fluctuation effects results in an inaccurate evaluation of the water table elevation at the land end of the aquifer, although no distinguishable difference is seen between the water tables near the shoreline. Where the landward boundary condition is a constant head, the effects of tidal fluctuations on sea-water intrusion are more pronounced than for cases where the landward boundary condition is a specified flux. Also it is shown that the effects of tidal fluctuations are more significant for a sloping beach than for a vertical shoreline and the salt water intrudes further inland for the sloping case. (C) 1999 Elsevier Science B.V. All rights reserved.

Finite element modelling of reactive mass transport problems in fluid-saturated porous media

We use the finite element method to solve reactive mass transport problems in fluid-saturated porous media. In particular, we discuss the mathematical expression of the chemical reaction terms involved in the mass transport equations for an isothermal, non-equilibrium chemical reaction. It has turned out that the Arrhenius law in chemistry is a good mathematical expression for such non-equilibrium chemical reactions especially from the computational point of view. Using the finite element method and the Arrhenius law, we investigate the distributions of PH (i.e. the concentration of H+) and the relevant reactive species in a groundwater system. Although the main focus of this study is on the contaminant transport problems in groundwater systems, the related numerical techniques and principles are equally applicable to the orebody formation problems in the geosciences. Copyright (C) 1999 John Wiley & Sons, Ltd.

Numerical and experimental study of seepage in unconfined aquifers with a periodic boundary condition

The assessment of groundwater conditions within an unconfined aquifer with a periodic boundary condition is of interest in many hydrological and environmental problems. A two-dimensional numerical model for density dependent variably saturated groundwater flow, SUTRA (Voss, C.I., 1984. SUTRA: a finite element simulation model for saturated-unsaturated, fluid-density dependent ground-water flow with energy transport or chemically reactive single species solute transport. US Geological Survey, National Center, Reston, VA) is modified in order to be able to simulate the groundwater flow in unconfined aquifers affected by a periodic boundary condition. The basic flow equation is changed from pressure-form to mixed-form. The model is also adjusted to handle a seepage-face boundary condition. Experiments are conducted to provide data for the groundwater response to the periodic boundary condition for aquifers with both vertical and sloping faces. The performance of the numerical model is assessed using those data. The results of pressure- and mixed-form approximations are compared and the improvement achieved through the mixed-form of the equation is demonstrated. The ability of the numerical model to simulate the water table and seepage-face is tested by modelling some published experimental data. Finally the numerical model is successfully verified against present experimental results to confirm its ability to simulate complex boundary conditions like the periodic head and the seepage-face boundary condition on the sloping face. (C) 1999 Elsevier Science B.V. All rights reserved.