Erosion/productivity modelling of maize farming in the Philippine uplands part III: Economic analysis of alternative farming methods

Two previous papers in this series (Nelson et al., this issue) described the use of the Agricultural Production Systems Simulator (APSIM) to simulate the effect of erosion on maize yields from open-field farming and hedgerow intercropping in the Philippine uplands. In this paper, maize yields simulated with APSIM are used to compare the economic viability of intercropping maize between leguminous shrub hedgerows with that of continuous and fallow open-field farming of maize. The analysis focuses on the economic incentives of upland farmers to adopt hedgerow intercropping, discussing farmers' planning horizons, access to credit and security of land tenure, as well as maize pricing in the Philippines. Insecure land tenure has limited the planning horizons of upland farmers, and high establishment costs reduce the economic viability of hedgerow intercropping relative to continuous and fallow open-field farming in the short term, In the long term, high discount rates and share-tenancy arrangements in which landlords do not contribute to establishment costs reduce the economic viability of hedgerow intercropping relative to fallow open-field farming, (C) 1998 Elsevier Science Ltd. 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.

Differentiation of soil properties related to the spatial association of wheat roots and soil macropores

Under certain soil conditions, e.g. hardsetting clay B-horizons of South-Eastern Australia, wheat plants do not perform as well as would be expected given measurements of bulk soil attributes. In such soils, measurement indicates that a large proportion (80%) of roots are preferentially located in the soil within 1 mm of macropores. This paper addresses the question of whether there are biological and soil chemical effects concomitant with this observed spatial relationship. The properties of soil manually dissected from the 1-3 mm wide region surrounding macropores, the macropore sheath, were compared to those that are measured in a conventional manner on the bulk soil. Field specimens of two different soil materials were dissected to examine biological differentiation. To ascertain whether the macropore sheath soil differs from rhizosphere soil, wheat was grown in structured and repacked cores under laboratory conditions. The macropore sheath soil contained more microbial biomass per unit mass than both the bulk soil and the rhizosphere. The bacterial population in the macropore sheath was able to utilise a wider range of carbon substrates and to a greater extent than the bacterial population in the corresponding bulk soil. These differences between the macropore sheath and bulk soil were almost non-existent in the repacked cores. Evidence for larger numbers of propagules of the broad host range fungus Pythium in the macropore sheath soil were also obtained.

Macropore sheath: quantification of plant root and soil macropore association

Plants require roots to supply water, nutrients and oxygen for growth. The spatial distribution of roots in relation to the macropore structure of the soil in which they are growing influences how effective they are at accessing these resources. A method for quantifying root-macropore associations from horizontal soil sections is illustrated using two black vertisols from the Darling Downs, Queensland, Australia. Two-dimensional digital images were obtained of the macropore structure and root distribution for an area 55 x 55 mm at a resolution of 64 mu m. The spatial distribution of roots was quantified over a range of distances using the K-function. In all specimens, roots were shown to be clustered at short distances (1-10 mm) becoming more random at longer distances. Root location in relation to macropores was estimated using the function describing the distance of each root to the nearest macropore. From this function, a summary variable, termed the macropore sheath, was defined. The macropore sheath is the distance from macropores within which 80% of roots are located. Measured root locations were compared to random simulations of root distribution to establish if there was a preferential association between roots and macropores. More roots were found in and around macropores than expected at random.

Empirical catchment-wide rainfall erosivity models for two rivers in the humid tropics of Australia

The concept of rainfall erosivity is extended to the estimation of catchment sediment yield and its variation over time. Five different formulations of rainfall erosivity indices, using annual, monthly and daily rainfall data, are proposed and tested on two catchments in the humid tropics of Australia. Rainfall erosivity indices, using simple power functions of annual and daily rainfall amounts, were found to be adequate in describing the interannual and seasonal variation of catchment sediment yield. The parameter values of these rainfall erosivity indices for catchment sediment yield are broadly similar to those for rainfall erosivity models in relation to the R-factor in the Universal Soil Loss Equation.

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.

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.

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.

Plant growth limitation and nutrient loss following piled burning in slash and burn agriculture

In the forest zone of Cameroon, small-scale family farmers practicing traditional slash and burn practices achieve a clear field by piled burning of the branches and trunks of cleared vegetation. Plant growth inhibition on ash patches, and the risk of nutrient loss from these areas, was evaluated on field plots on which 0.5 t m(-2) or 1.0 t m(-2) of wood was piled and burnt, and in laboratory studies. The ash produced by burning was strongly alkaline, and laboratory bio-assessment studies showed that the saline, high pH conditions produced in ash patches prevented germination and plant growth for up to two wet seasons, as is observed in the field. Field and laboratory studies demonstrated rapid release (1 wet season) of K and S from the ash and the loss of a substantial portion of these nutrients from the soil profile by leaching. In contrast, leaching carries Mg from the ash gradually (3 to 4 wet seasons), while Ca, Cu, Zn and P are leached slowly. The nutrients contained in ash patches are considered at risk of loss both through leaching (K and S) and by erosion of ash (Ca, Mg, Cu, Zn and P). Farmers should be encouraged to spread ash patches prior to cultivation in order to exploit the nutrient content of ash and to lessen the risk of nutrient loss.

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.

Australia-wide predictions of soil properties using decision trees

This paper describes the construction of Australia-wide soil property predictions from a compiled national soils point database. Those properties considered include pH, organic carbon, total phosphorus, total nitrogen, thickness. texture, and clay content. Many of these soil properties are used directly in environmental process modelling including global climate change models. Models are constructed at the 250-m resolution using decision trees. These relate the soil property to the environment through a suite of environmental predictors at the locations where measurements are observed. These models are then used to extend predictions to the continental extent by applying the rules derived to the exhaustively available environmental predictors. The methodology and performance is described in detail for pH and summarized for other properties. Environmental variables are found to be important predictors, even at the 250-m resolution at which they are available here as they can describe the broad changes in soil property.