Ammonia volatilization in soil treated with tannery sludge

The utilization of tannery sludge in agricultural areas can be an alternative for its disposal and recycling. Despite this procedure may cause the loss of nitrogen by ammonia volatilization, there is no information about this process in tropical soils. For two years a field experiment was carried out in Rolandia (Parana State, Brazil), to evaluate the amount of NH(3) volatilization due to tannery sludge application on agricultural soil. The doses of total N applied varied from zero to 1200 kg ha(-1), maintained at the surface for 89 days, as usual in this region. The alkalinity of the tannery sludge used was equivalent to between 262 and 361 g CaCO(3) per kg. Michaelis-Menten equation was adequate to estimate NH(3)-N volatilization kinetics. The relation between total nitrogen applied as tannery sludge and the potentially volatilized NH(3)-N, calculated by the chemical-kinetics equation resulted in an average determination coefficient of 0.87 (P > 0.01). In this period, the amount of volatilized NH(3) was more intense during the first 30 days; the time to reach half of the maximum NH(3) volatilization (K(m)) was 13 an 9 days for the first and second experiments, respectively. The total loss as ammonia in the whole period corresponded in average to 17.5% of the total N applied and to 35% of the NH(4)(+)-N present in the sludge. If tannery sludge is to be surface applied to supply N for crops, the amounts lost as NH(3) must be taken into consideration. (C) 2010 Elsevier Ltd. All rights reserved.

Contribution of soil organic carbon to the ion exchange capacity of tropical soils

Highly weathered soils represent about 3 billion ha of the tropical region. Oxisols represent about 60% of the Brazilian territory (more than 5 million km 2), in areas of great agricultural importance. Soil organic carbon (SOC) can be responsible for more than 80% of the cation exchange capacity (CEC) of highly weathered soils, such as Oxisols and Ultisols. The objective of this study was to estimate the contribution of the SOC to the CEC of Brazilian soils from different orders. Surface samples (0.0 to 0.2 m) of 30 uncultivated soils (13 Oxisols, 6 Ultisols, 5 Alfisols, 3 Entisols, I Histosol, 1 Inceptisol. and I Molisol), under native forests and from reforestation sites from Sao Paulo State, Brazil, were collected in order to obtain a large variation of (electro)chemical, physical, and mineralogical soil attributes. Total content of SOC was quantified by titulometric and colorimetric methods. Effective cation exchange capacity (ECEC) was obtained by two methods: the indirect method-summation-estimated the ECECi from the sum of basic cations (Ca+ Mg+ K+ Na) and exchangeable Al; and the direct ECECd obtained by the compulsive exchange method, using unbuffered BaCl2 solution. The contribution of SOC to the soil CEC was estimated by the Bennema statistical method. The amount of SOC var ied from 6.6 g kg(-1) to 213.4 g kg(-1). while clay contents varied from 40 g kg(-1) to 716 g kg(-1). Soil organic carbon contents were strongly associated to the clay contents, suggesting that clay content was the primary variable in controling the variability of SOC contents in the samples. Cation exchange capacity varied from 7.0 mmol(c) kg(-1) to 137.8 mmol(c) kg(-1) and had a positive Correlation with SOC. The mean contribution (per grain) of the SOC (1.64 mmol(c)) for the soil CEC was more than 44 times higher than the contribution of the clay fraction (0.04 mmol(c),). A regression model that considered the SOC content as the only significant variable explained 60% of the variation in the soil total CEC. The importance of SOC was related to soil pedogenetic process, since its contribution to the soil CEC was more evident in Oxisols with predominance of Fe and Al (oxihydr)oxides in the mineral fraction or in Ultisols, that presented illuviated clay. The influence of SOC in the sign and in the magnitude of the net charge of soils reinforce the importance of agricultural management systems that preserve high levels of SOC, in order to improve their sustainability.

Deciphering the weathering processes using environmental mineralogy and geochemistry: Towards an integrated model of laterite and podzol genesis in the Upper Amazon Basin

This article reports major results from collaborative research between France and Brazil on soil and water systems, carried out in the Upper Amazon Basin. It reveals the weathering processes acting in the partly inundated, low elevation plateaus of the Basin, mostly covered by evergreen forest. Our findings are based on geochemical data and mineral spectroscopy that probe the crystal chemistry of Fe and Al in mineral phases (mainly kaolinite, Al- and Fe-(hydr)oxides) of tropical soils (laterites). These techniques reveal crystal alterations in mineral populations of different ages and changes of metal speciation associated with mineral or organic phases. These results provide an integrated model of soil formation and changes (from laterites to podzols) in distinct hydrological compartments of the Amazon landscapes and under altered water regimes. (C) 2010 Academie des sciences. Published by Elsevier Masson SAS. All rights reserved.

Leaf area index estimation with MODIS reflectance time series and model inversion during full rotations of Eucalyptus plantations

The leaf area index (LAI) of fast-growing Eucalyptus plantations is highly dynamic both seasonally and interannually, and is spatially variable depending on pedo-climatic conditions. LAI is very important in determining the carbon and water balance of a stand, but is difficult to measure during a complete stand rotation and at large scales. Remote-sensing methods allowing the retrieval of LAI time series with accuracy and precision are therefore necessary. Here, we tested two methods for LAI estimation from MODIS 250m resolution red and near-infrared (NIR) reflectance time series. The first method involved the inversion of a coupled model of leaf reflectance and transmittance (PROSPECT4), soil reflectance (SOILSPECT) and canopy radiative transfer (4SAIL2). Model parameters other than the LAI were either fixed to measured constant values, or allowed to vary seasonally and/or with stand age according to trends observed in field measurements. The LAI was assumed to vary throughout the rotation following a series of alternately increasing and decreasing sigmoid curves. The parameters of each sigmoid curve that allowed the best fit of simulated canopy reflectance to MODIS red and NIR reflectance data were obtained by minimization techniques. The second method was based on a linear relationship between the LAI and values of the GEneralized Soil Adjusted Vegetation Index (GESAVI), which was calibrated using destructive LAI measurements made at two seasons, on Eucalyptus stands of different ages and productivity levels. The ability of each approach to reproduce field-measured LAI values was assessed, and uncertainty on results and parameter sensitivities were examined. Both methods offered a good fit between measured and estimated LAI (R(2) = 0.80 and R(2) = 0.62 for model inversion and GESAVI-based methods, respectively), but the GESAVI-based method overestimated the LAI at young ages. (C) 2010 Elsevier Inc. All rights reserved.

A split-pot experiment with sorghum to test a root water uptake partitioning model

Correct modeling of root water uptake partitioning over depth is an important issue in hydrological and crop growth models. Recently a physically based model to describe root water uptake was developed at single root scale and upscaled to the root system scale considering a homogeneous distribution of roots per soil layer. Root water uptake partitioning is calculated over soil layers or compartments as a function of respective soil hydraulic conditions, specifically the soil matric flux potential, root characteristics and a root system efficiency factor to compensate for within-layer root system heterogeneities. The performance of this model was tested in an experiment performed in two-compartment split-pot lysimeters with sorghum plants. The compartments were submitted to different irrigation cycles resulting in contrasting water contents over time. The root system efficiency factor was determined to be about 0.05. Release of water from roots to soil was predicted and observed on several occasions during the experiment; however, model predictions suggested root water release to occur more often and at a higher rate than observed. This may be due to not considering internal root system resistances, thus overestimating the ease with which roots can act as conductors of water. Excluding these erroneous predictions from the dataset, statistical indices show model performance to be of good quality.

Combining GIS and Expert Opinion to Model Landscapes for a Smallholder Forest Extension Program in Leyte, the Philippines

A Geographic Information System (GIS) was used to model datasets of Leyte Island, the Philippines, to identify land which was suitable for a forest extension program on the island. The datasets were modelled to provide maps of the distance of land from cities and towns, land which was a suitable elevation and slope for smallholder forestry and land of various soil types. An expert group was used to assign numeric site suitabilities to the soil types and maps of site suitability were used to assist the selection of municipalities for the provision of extension assistance to smallholders. Modelling of the datasets was facilitated by recent developments of the ArcGIS® suite of computer programs and derivation of elevation and slope was assisted by the availability of digital elevation models (DEM) produced by the Shuttle Radar Topography (SRTM) mission. The usefulness of GIS software as a decision support tool for small-scale forestry extension programs is discussed.

Improving wheat simulation capabilities in Australia from a cropping systems perspective - III. The integrated wheat model (I_WHEAT)

Previous work has identified several short-comings in the ability of four spring wheat and one barley model to simulate crop processes and resource utilization. This can have important implications when such models are used within systems models where final soil water and nitrogen conditions of one crop define the starting conditions of the following crop. In an attempt to overcome these limitations and to reconcile a range of modelling approaches, existing model components that worked demonstrably well were combined with new components for aspects where existing capabilities were inadequate. This resulted in the Integrated Wheat Model (I_WHEAT), which was developed as a module of the cropping systems model APSIM. To increase predictive capability of the model, process detail was reduced, where possible, by replacing groups of processes with conservative, biologically meaningful parameters. I_WHEAT does not contain a soil water or soil nitrogen balance. These are present as other modules of APSIM. In I_WHEAT, yield is simulated using a linear increase in harvest index whereby nitrogen or water limitations can lead to early termination of grainfilling and hence cessation of harvest index increase. Dry matter increase is calculated either from the amount of intercepted radiation and radiation conversion efficiency or from the amount of water transpired and transpiration efficiency, depending on the most limiting resource. Leaf area and tiller formation are calculated from thermal time and a cultivar specific phyllochron interval. Nitrogen limitation first reduces leaf area and then affects radiation conversion efficiency as it becomes more severe. Water or nitrogen limitations result in reduced leaf expansion, accelerated leaf senescence or tiller death. This reduces the radiation load on the crop canopy (i.e. demand for water) and can make nitrogen available for translocation to other organs. Sensitive feedbacks between light interception and dry matter accumulation are avoided by having environmental effects acting directly on leaf area development, rather than via biomass production. This makes the model more stable across environments without losing the interactions between the different external influences. When comparing model output with models tested previously using data from a wide range of agro-climatic conditions, yield and biomass predictions were equal to the best of those models, but improvements could be demonstrated for simulating leaf area dynamics in response to water and nitrogen supply, kernel nitrogen content, and total water and nitrogen use. I_WHEAT does not require calibration for any of the environments tested. Further model improvement should concentrate on improving phenology simulations, a more thorough derivation of coefficients to describe leaf area development and a better quantification of some processes related to nitrogen dynamics. (C) 1998 Elsevier Science B.V.

Subantarctic Macquarie Island - a model ecosystem for studying animal-derived nitrogen sources using N-15 natural abundance

Plants collected from diverse sites on subantarctic Macquarie Island varied by up to 30 parts per thousand in their leaf delta(15)N values. N-15 natural abundance of plants, soils, animal excrement and atmospheric ammonia suggest that the majority of nitrogen utilised by plants growing in the vicinity of animal colonies or burrows is animal-derived. Plants growing near scavengers and animal higher in the food chain had highly enriched delta(15)N values (mean = 12.9 parts per thousand), reflecting the highly enriched signature of these animals' excrement, while plants growing near nesting penguins and albatross, which have an intermediate food chain position, had less enriched delta(15)N values (> 6 parts per thousand). Vegetation in areas affected by rabbits had lower delta(15)N values (mean = 1.2 parts per thousand), while the highly depleted delta(15)N values (below -5 parts per thousand) of plants at upland plateau sites inland of penguin colonies, suggested that a portion of their nitrogen is derived from ammonia (mean N-15 = -10 parts per thousand) lost during the degradation of penguin guano. Vegetation in a remote area had delta(15)N values near -2 parts per thousand. These results contrast with arctic and subarctic studies that attribute large variations in plant N-15 values to nitrogen partitioning in nitrogen-limited environments. Here, plant N-15 reflects the N-15 Of the likely nitrogen sources utilised by plants.

A mathematical model for estimating the extent of solute- and water-flux heterogeneity in multiple sample percolation experiments

Multiple sampling is widely used in vadose zone percolation experiments to investigate the extent in which soil structure heterogeneities influence the spatial and temporal distributions of water and solutes. In this note, a simple, robust, mathematical model, based on the beta-statistical distribution, is proposed as a method of quantifying the magnitude of heterogeneity in such experiments. The model relies on fitting two parameters, alpha and zeta to the cumulative elution curves generated in multiple-sample percolation experiments. The model does not require knowledge of the soil structure. A homogeneous or uniform distribution of a solute and/or soil-water is indicated by alpha = zeta = 1, Using these parameters, a heterogeneity index (HI) is defined as root 3 times the ratio of the standard deviation and mean. Uniform or homogeneous flow of water or solutes is indicated by HI = 1 and heterogeneity is indicated by HI > 1. A large value for this index may indicate preferential flow. The heterogeneity index relies only on knowledge of the elution curves generated from multiple sample percolation experiments and is, therefore, easily calculated. The index may also be used to describe and compare the differences in solute and soil-water percolation from different experiments. The use of this index is discussed for several different leaching experiments. (C) 1999 Elsevier Science B.V. All rights reserved.

Commentary: Using the convection-dispersion model and transit time density functions in the analysis of organ distribution kinetics

The convection-dispersion model and its extended form have been used to describe solute disposition in organs and to predict hepatic availabilities. A range of empirical transit-time density functions has also been used for a similar purpose. The use of the dispersion model with mixed boundary conditions and transit-time density functions has been queried recently by Hisaka and Sugiyanaa in this journal. We suggest that, consistent with soil science and chemical engineering literature, the mixed boundary conditions are appropriate providing concentrations are defined in terms of flux to ensure continuity at the boundaries and mass balance. It is suggested that the use of the inverse Gaussian or other functions as empirical transit-time densities is independent of any boundary condition consideration. The mixed boundary condition solutions of the convection-dispersion model are the easiest to use when linear kinetics applies. In contrast, the closed conditions are easier to apply in a numerical analysis of nonlinear disposition of solutes in organs. We therefore argue that the use of hepatic elimination models should be based on pragmatic considerations, giving emphasis to using the simplest or easiest solution that will give a sufficiently accurate prediction of hepatic pharmacokinetics for a particular application. (C) 2000 Wiley-Liss Inc. and the American Pharmaceutical Association J Pharm Sci 89:1579-1586, 2000.

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.

Three-dimensional model for multi-component reactive transport with variable density groundwater flow

PHWAT is a new model that couples a geochemical reaction model (PHREEQC-2) with a density-dependent groundwater flow and solute transport model (SEAWAT) using the split-operator approach. PHWAT was developed to simulate multi-component reactive transport in variable density groundwater flow. Fluid density in PHWAT depends not on only the concentration of a single species as in SEAWAT, but also the concentrations of other dissolved chemicals that can be subject to reactive processes. Simulation results of PHWAT and PHREEQC-2 were compared in their predictions of effluent concentration from a column experiment. Both models produced identical results, showing that PHWAT has correctly coupled the sub-packages. PHWAT was then applied to the simulation of a tank experiment in which seawater intrusion was accompanied by cation exchange. The density dependence of the intrusion and the snow-plough effect in the breakthrough curves were reflected in the model simulations, which were in good agreement with the measured breakthrough data. Comparison simulations that, in turn, excluded density effects and reactions allowed us to quantify the marked effect of ignoring these processes. Next, we explored numerical issues involved in the practical application of PHWAT using the example of a dense plume flowing into a tank containing fresh water. It was shown that PHWAT could model physically unstable flow and that numerical instabilities were suppressed. Physical instability developed in the model in accordance with the increase of the modified Rayleigh number for density-dependent flow, in agreement with previous research. (c) 2004 Elsevier Ltd. All rights reserved.

A cost-benefit analysis of hedgerow intercropping in the Philippine uplands using the SCUAF model

Soil erosion in the Philippine uplands is severe. Hedgerow intercropping is widely advocated as an effective means of controlling soil erosion from annual cropping systems in the uplands. However, few farmers adopt hedgerow intercropping even in areas where it has been vigorously promoted. This may be because farmers find hedgerow intercropping to be uneconomic compared to traditional methods of farming. This paper reports a cost-benefit analysis comparing the economic returns from traditional maize farming with those from hedgerow intercropping in an upland community with no past adoption of hedgerows. A simple erosion/productivity model, Soil Changes Under Agroforestry (SCUAF), is used to predict maize yields over 25 years. Economic data were collected through key informant surveys with experienced maize farmers in an upland community. Traditional methods of open-field farming of maize are economically attractive to farmers in the Philippine uplands. In the short term, establishment costs are a major disincentive to the adoption of hedgerow intercropping. In the long term, higher economic returns from hedgerow intercropping compared to open-field farming are realised, but these lie beyond farmers' limited planning horizons.