Solute transport by surface runoff from low-angle slopes: theory and application

The removal of chemicals in solution by overland how from agricultural land has the potential to be a significant source of chemical loss where chemicals are applied to the soil surface, as in zero tillage and surface-mulched farming systems. Currently, we lack detailed understanding of the transfer mechanism between the soil solution and overland flow, particularly under field conditions. A model of solute transfer from soil solution to overland flow was developed. The model is based on the hypothesis that a solute is initially distributed uniformly throughout the soil pore space in a thin layer at the soil surface. A fundamental assumption of the model is that at the time runoff commences, any solute at the soil surface that could be transported into the soil with the infiltrating water will already have been convected away from the area of potential exchange. Solute remaining at the soil surface is therefore not subject to further infiltration and may be approximated as a layer of tracer on a plane impermeable surface. The model fitted experimental data very well in all but one trial. The model in its present form focuses on the exchange of solute between the soil solution and surface water after the commencement of runoff. Future model development requires the relationship between the mass transfer parameters of the model and the time to runoff: to be defined. This would enable the model to be used for extrapolation beyond the specific experimental results of this study. The close agreement between experimental results and model simulations shows that the simple transfer equation proposed in this study has promise for estimating solute loss to surface runoff. Copyright (C) 2000 John Wiley & Sons, Ltd.

X-ray absorption and phase contrast imaging to study the interplay between plant roots and soil structure

Plant performance is, at least partly, linked to the location of roots with respect to soil structure features and the micro-environment surrounding roots. Measurements of root distributions from intact samples, using optical microscopy and field tracings have been partially successful but are imprecise and labour-intensive. Theoretically, X-ray computed micro-tomography represents an ideal solution for non-invasive imaging of plant roots and soil structure. However, before it becomes fast enough and affordable or easily accessible, there is still a need for a diagnostic tool to investigate root/soil interplay. Here, a method for detection of undisturbed plant roots and their immediate physical environment is presented. X-ray absorption and phase contrast imaging are combined to produce projection images of soil sections from which root distributions and soil structure can be analyzed. The clarity of roots on the X-ray film is sufficient to allow manual tracing on an acetate sheet fixed over the film. In its current version, the method suffers limitations mainly related to (i) the degree of subjectivity associated with manual tracing and (ii) the difficulty of separating live and dead roots. The method represents a simple and relatively inexpensive way to detect and quantify roots from intact samples and has scope for further improvements. In this paper, the main steps of the method, sampling, image acquisition and image processing are documented. The potential use of the method in an agronomic perspective is illustrated using surface and sub-surface soil samples from a controlled wheat trial. Quantitative characterization of root attributes, e.g. radius, length density, branching intensity and the complex interplay between roots and soil structure, is presented and discussed.

Soil carbon determination by high temperature combustion - A comparison with dichromate oxidation procedures and the influence of charcoal and carbonate carbon on the measured value

The measurement of organic carbon in soils has traditionally used dichromate oxidation procedures including the Wakley and Black and the Heanes methods. The measurement of carbon in soils by high temperature combustion is now widely used providing a rapid automated procedure without the use of toxic chemicals. This procedure however measures total carbon thus requiring some means of correction for soil samples containing carbonate and charcoal forms of carbon. This paper examines the effects of known additions of charcoal to a range of soil types on the results obtained by the Walkley and Black, Heanes and combustion methods. The results show, that while the charcoal carbon does not react under Walkley and Black conditions, some proportion does so with the Heanes method. A comparison of six Australian Soil and Plant Analysis Council reference soil samples by the three methods showed good agreement between the Heanes method, the combustion method and only slightly lower recoveries by the Walkley and Black procedure. Carbonate carbon will cause an overestimation of soil organic carbon by the combustion method thus requiring a separate determination of carbonate carbon to be applied as a correction. This work shows that a suitable acid pre-treatment of alkaline soils in the sample boats followed by a drying step eliminates the carbonate carbon prior to combustion and the need for an additional measurement. The measurement of carbon in soils by high temperature combustion in an oxygen atmosphere has been shown to be a rapid and reliable method capable of producing results in good agreement with one of the established dichromate oxidation procedures.

Numerical modelling of single-layer folding: clarification of an issue regarding the possible effect of computer codes and the influence of initial irregularities

The influence of initial perturbation geometry and material propel-ties on final fold geometry has been investigated using finite-difference (FLAC) and finite-element (MARC) numerical models. Previous studies using these two different codes reported very different folding behaviour although the material properties, boundary conditions and initial perturbation geometries were similar. The current results establish that the discrepancy was not due to the different computer codes but due to the different strain rates employed in the two previous studies (i.e. 10(-6) s(-1) in the FLAC models and 10(-14) s(-1) in the MARC models). As a result, different parts of the elasto-viscous rheological field were bring investigated. For the same material properties, strain rate and boundary conditions, the present results using the two different codes are consistent. A transition in Folding behaviour, from a situation where the geometry of initial perturbation determines final fold shape to a situation where material properties control the final geometry, is produced using both models. This transition takes place with increasing strain rate, decreasing elastic moduli or increasing viscosity (reflecting in each case the increasing influence of the elastic component in the Maxwell elastoviscous rheology). The transition described here is mechanically feasible but is associated with very high stresses in the competent layer (on the order of GPa), which is improbable under natural conditions. (C) 2000 Elsevier Science Ltd. All rights reserved.

Returns to farm-level soil conservation on tropical steep slopes: The case of the Eritrean highlands

This study conducts an economic analysis of investment in simple soil conservation technologies in the highlands of Eritrea. The data used in the analysis were obtained from a farm survey and supplemented with data from secondary sources. Risk analysis techniques are used to take account of the uncertainties regarding the relationship between soil erosion and crop yield. The financial analysis reveals negative net present values (NPVs) and internal rates of return (IRRs) below 12 per cent for various slope categories. On the other hand, the economic analysis returns positive NPVs and IRRs of over 20 per cent. The results clearly indicate that in-vestment in soil conservation technology may not be a viable short-term proposition from the farmer's point of view and yet the net social benefits are positive. There is a strong case for government to provide incentives for soil conservation in view of the economic benefits.

Disaggregation of polygons of surficial geology and soil maps using spatial modelling and legacy data

Examples from the Murray-Darling basin in Australia are used to illustrate different methods of disaggregation of reconnaissance-scale maps. One approach for disaggregation revolves around the de-convolution of the soil-landscape paradigm elaborated during a soil survey. The descriptions of soil ma units and block diagrams in a soil survey report detail soil-landscape relationships or soil toposequences that can be used to disaggregate map units into component landscape elements. Toposequences can be visualised on a computer by combining soil maps with digital elevation data. Expert knowledge or statistics can be used to implement the disaggregation. Use of a restructuring element and k-means clustering are illustrated. Another approach to disaggregation uses training areas to develop rules to extrapolate detailed mapping into other, larger areas where detailed mapping is unavailable. A two-level decision tree example is presented. At one level, the decision tree method is used to capture mapping rules from the training area; at another level, it is used to define the domain over which those rules can be extrapolated. (C) 2001 Elsevier Science B.V. All rights reserved.

Finite element modelling of heat transfer through permeable cracks in hydrothermal systems with upward throughflow

We use the finite element method to model the heat transfer phenomenon through permeable cracks in hydrothermal systems with upward throughflow. Since the finite element method is an approximate numerical method, the method must be validated before it is used to soh,e any new, kind of problem. However, the analytical solution, which can be used to validate the finite element method and other numerical methods, is rather limited in the literature, especially, for the problem considered here. Keeping this in mind, we have derived analytical solutions for the temperature distribution along the vertical axis of a crack in a fluid-saturated porous layer. After the finite element method is validated by comparing the numerical solution with the analytical solution for the same benchmark problem, it is used to investigate the pore-fluid flow and heat transfer in layered hydrothermal systems with vertical permeable cracks. The related analytical and numerical results have demonstrated that vertical cracks are effective and efficient members to transfer heat energy from the bottom section to the top section in hydrothermal systems with upward throughflow.

Spatial data mining for enhanced soil map modelling

The principle of using induction rules based on spatial environmental data to model a soil map has previously been demonstrated Whilst the general pattern of classes of large spatial extent and those with close association with geology were delineated small classes and the detailed spatial pattern of the map were less well rendered Here we examine several strategies to improve the quality of the soil map models generated by rule induction Terrain attributes that are better suited to landscape description at a resolution of 250 m are introduced as predictors of soil type A map sampling strategy is developed Classification error is reduced by using boosting rather than cross validation to improve the model Further the benefit of incorporating the local spatial context for each environmental variable into the rule induction is examined The best model was achieved by sampling in proportion to the spatial extent of the mapped classes boosting the decision trees and using spatial contextual information extracted from the environmental variables.

A strategy to fill gaps in soil survey over large spatial extents: an example from the Murray-Darling basin of Australia

We re-mapped the soils of the Murray-Darling Basin (MDB) in 1995-1998 with a minimum of new fieldwork, making the most out of existing data. We collated existing digital soil maps and used inductive spatial modelling to predict soil types from those maps combined with environmental predictor variables. Lithology, Landsat Multi Spectral Scanner (Landsat MSS), the 9-s digital elevation model (DEM) of Australia and derived terrain attributes, all gridded to 250-m pixels, were the predictor variables. Because the basin-wide datasets were very large data mining software was used for modelling. Rule induction by data mining was also used to define the spatial domain of extrapolation for the extension of soil-landscape models from existing soil maps. Procedures to estimate the uncertainty associated with the predictions and quality of information for the new soil-landforms map of the MDB are described. (C) 2002 Elsevier Science B.V. All rights reserved.

Scope for further analytical solutions for constant flux infiltration into a semi-infinite soil profile or redistribution in a finite soil profile

[1] We attempt to generate new solutions for the moisture content form of the one-dimensional Richards' [1931] equation using the Lisle [1992] equivalence mapping. This mapping is used as no more general set of transformations exists for mapping the one-dimensional Richards' equation into itself. Starting from a given solution, the mapping has the potential to generate an infinite number of new solutions for a series of nonlinear diffusivity and hydraulic conductivity functions. We first seek new analytical solutions satisfying Richards' equation subject to a constant flux surface boundary condition for a semi-infinite dry soil, starting with the Burgers model. The first iteration produces an existing solution, while subsequent iterations are shown to endlessly reproduce this same solution. Next, we briefly consider the problem of redistribution in a finite-length soil. In this case, Lisle's equivalence mapping is generalized to account for arbitrary initial conditions. As was the case for infiltration, however, it is found that new analytical solutions are not generated using the equivalence mapping, although existing solutions are recovered.

Analysis of soil wetting and solute transport in subsurface trickle irrigation

The increased use of trickle or drip irrigation is seen as one way of helping to improve the sustainability of irrigation systems around the world. However, soil water and solute transport properties and soil profile characteristics are often not adequately incorporated in the design and management of trickle systems. In this paper, we describe results of a simulation study designed to highlight the impacts of soil properties on water and solute transport from buried trickle emitters. The analysis addresses the influence of soil hydraulic properties, soil layering, trickle discharge rate, irrigation frequency, and timing of nutrient application on wetting patterns and solute distribution. We show that (1) trickle irrigation can improve plant water availability in medium and low permeability fine-textured soils, providing that design and management are adapted to account for their soil hydraulic properties, (2) in highly permeable coarse-textured soils, water and nutrients move quickly downwards from the emitter, making it difficult to wet the near surface zone if emitters are buried too deep, and (3) changing the fertigation strategy for highly permeable coarse-textured soils to apply nutrients at the beginning of an irrigation cycle can maintain larger amounts of nutrient near to and above the emitter, thereby making them less susceptible to leaching losses. The results demonstrate the need to account for differences in soil hydraulic properties and solute transport when designing irrigation and fertigation management strategies. Failure to do this will result in inefficient systems and lost opportunities for reducing the negative environmental impacts of irrigation.

Use of bulk density profiles from X-radiography to examine structural crust models

Exponential and sigmoidal functions have been suggested to describe the bulk density profiles of crusts. The present work aims to evaluate these conceptual models using high resolution X-radiography. Repacked seedbeds from two soil materials, air-dried or prewetted by capillary rise, were subjected to simulated rain, which resulted in three types of structural crusts, namely, slaking, infilling, and coalescing. Bulk density distributions with depth were generated using high-resolution (70 mum), calibrated X-ray images of slices from the resin-impregnated crusted seedbeds. The bulk density decreased progressively with depth, which supports the suggestion that a crust should be considered as a nonuniform layer. For the slaking and the coalescing crusts, the exponential function underestimated the strong change in bulk density across the morphologically defined transition between the crust and the underlying material; the sigmoidal function provided a better description. Neither of these crust models effectively described the shape of the bulk density profiles through the whole seedbed. Below the infilling and slaking crusts, bulk density increased linearly with depth as a result of slumping. In the coalescing crusted seedbed, the whole seedbed uniformly collapsed and most of the bulk density change within the crust could be ascribed to slumping (0.33 g cm(-3)) rather than to crusting (0.12 g cm(-3)). Finally, (i) X-radiography appears as a unique tool to generate high resolution bulk density profiles and (ii) in structural crusts, bulk density profiles could be modeled using the existing exponential and sigmoidal crusting models, provided a slumping model would be coupled.

Platinum nanoparticles embledded in layer-by-layer films from SnO2/polyallylamine for ethanol electrooxidation

Self-assembled films from SnO2 and polyallylamine (PAH) were deposited on gold via ionic attraction by the layer-by-layer(LbL) method. The modified electrodes were immersed into a H2PtCl6 solution, a current of 100 mu A was applied, and different electrodeposition times were used. The SnO2/PAH layers served as templates to yield metallic platinum with different particle sizes. The scanning tunnel microscopy images show that the particle size increases as a function of electrodeposition time. The potentiodynamic profile of the electrodes changes as a function of the electrodeposition time in 0.5 mol L-1 H2SO4, at a sweeping rate of 50mVs(-1). Oxygen-like species are formed at less positive potentials for the Pt-SnO2/PAH film in the case of the smallest platinum particles. Electrochemical impedance spectroscopy measurements in acid medium at 0.7 V show that the charge transfer resistance normalized by the exposed platinum area is 750 times greater for platinum electrode (300 k Omega cm(2)) compared with the Pt-SnO2/PAH film with 1 min of electrodeposition (0.4 k Omega cm(2)). According to the Langmuir-Hinshelwood bifunctional mechanism, the high degree of coverage with oxygen-like species on the platinum nanoparticles is responsible for the electrocatalytic activity of the Pt-SnO2/PAH concerning ethanol electrooxidation. With these features, this Pt-SnO2/PAH film may be grown on a proton exchange membrane (PEM) in direct ethanol fuel cells (DEFC). (c) 2008 Elsevier B.V. All rights reserved.

Impedance characteristics of the diamond/carbon fiber electrodes for electrical double-layer capacitor

The electrochemical performance of carbon fibers (CF) and boron-doped diamond electrodes grown on carbon fiber substrate (BDD/CF) was studied. CF substrates were obtained from polyacrylonitrile precursor heat treated at two different temperatures of 1000 and 2000 degrees C to produce the desirable CF carbon graphitization index. This graphitization process influenced the CF conductivity and its chemical surface, also analyzed from X-ray photoelectron spectroscopy measurements. These three-dimensional CF structures allowed a high incorporation of diamond films compared to other carbon substrates such as glass carbon or HOPG. The electrochemical responses, from these four classes of electrodes, were evaluated focusing their application as electrical double-layer capacitors using cyclic voltammetry and impedance measurements. Cyclic voltammetry results revealed that the electrode formed from BDD grown on CF-2000 presented a typical capacitor behavior with the best rectangular shape, compared to those electrodes of CF or BDD/CF-1000. Furthermore, the BDD/CF-2000 electrode presented the lowest impedance, associated to its significant capacitance value of 1940 mu F/cm(2) taking into account the BDD films. This behavior was attributed to the strong dependence between diamond coating texture and the CF graphitization temperature. The largest surface area of BDD/CF-2000 was promoted by its singular film growth mechanism associated to the substrate chemical surface. (c) 2008 Elsevier B.V. All rights reserved.

Surface modification of metals by calcium carbonate thin films on a layer-by-layer polyelectrolyte matrix

A multilayer organic film containing poly(acrylic acid) and chitosan was fabricated on a metallic support by means of the layer-by-layer technique. This film was used as a template for calcium carbonate crystallization and presents two possible binding sites where the nucleation may be initiated, either calcium ions acting as counterions of the polyelectrolyte or those trapped in the template gel network formed by the polyelectrolyte chains. Calcium carbonate formation was carried out by carbon dioxide diffusion, where CO, was generated from ammonium carbonate decomposition. The CaCO3 nanocrystals obtained, formed a dense, homogeneous, and continuous film. Vaterite and calcite CaCO3 crystalline forms were detected. (c) 2007 Elsevier B.V All rights reserved.