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.

High-resolution bulk density images, using calibrated X-ray radiography of impregnated soil slices

Bulk density of undisturbed soil samples can be measured using computed tomography (CT) techniques with a spatial resolution of about 1 mm. However, this technique may not be readily accessible. On the other hand, x-ray radiographs have only been considered as qualitative images to describe morphological features. A calibration procedure was set up to generate two-dimensional, high-resolution bulk density images from x-ray radiographs made with a conventional x-ray diffraction apparatus. Test bricks were made to assess the accuracy of the method. Slices of impregnated soil samples were made using hardsetting seedbeds that had been gamma scanned at 5-mm depth increments in a previous study. The calibration procedure involved three stages: (i) calibration of the image grey levels in terms of glass thickness using a staircase made from glass cover slips, (ii) measurement of ratio between the soil and resin mass attenuation coefficients and the glass mass attenuation coefficient, using compacted bricks of known thickness and bulk density, and (iii) image correction accounting for the heterogeneity of the irradiation field. The procedure was simple, rapid, and the equipment was easily accessible. The accuracy of the bulk density determination was good (mean relative error 0.015), The bulk density images showed a good spatial resolution, so that many structural details could be observed. The depth functions were consistent with both the global shrinkage and the gamma probe data previously obtained. The suggested method would be easily applied to the new fuzzy set approach of soil structure, which requires generation of bulk density images. Also, it would be an invaluable tool for studies requiring high-resolution bulk density measurement, such as studies on soil surface crusts.

A technique to assess small-scale heterogeneity of chemical properties in soil aggregates

A new method is presented which allows the separation of the soil aggregate exterior from the aggregate core. The method employs a combination of aggregate freezing with rapid separation of aggregate exteriors using ultrasonic energy. The factors influencing the thickness of the removed aggregate surface layer include water content of the aggregate prior to freezing, temperature difference between that of the frozen aggregate and that of the liquid it is submerged in during sonification, sonification time and energy, and the type of the immersion liquid. The success of the method and the thickness of the removed aggregate surface were examined using barium ( Ba2+) as a tracer. Barium ( as BaCl2) is rapidly absorbed by soil and is present at only very low levels in natural soils. Surface layers of 0.2 - 0.4 cm thickness were successfully removed from aggregates of 1 - 4 cm diameter. Two examples are given from soils in northern NSW to demonstrate the occurrence of small- scale heterogeneity in soil chemical properties. Compared with the surface fraction, a 4 - 7% higher calcium concentration was found in the core fraction of a clay loam soil ( Dermosol). Conversely, on a cracking clay soil ( Vertosol), atrazine concentration was around 15 times greater in the aggregate surface fractions compared with core fractions. Compared with the traditional estimation of soil chemical properties on homogenised bulk soil samples, it is suggested that separate analysis of aggregate surface and core fractions could provide useful additional information on the relationships between soil properties and environmental responses.

Structure repair of a compacted vertisol with wet-dry cycles and crops

We hypothesized that the four rotation crops: wheat (Triticum aestivum L.), sorghum [Sorghum bicolor (L.) Merr.], lablab [Lablab purpureus (L.) Sweet] and mung bean [ Vigna radiata (L.) R. Wilczek] differ in their ability to repair soil structure. The study was conducted on a Typic Haplustert, Queensland, Australia, locally termed a Black Earth and considered a prime cropping soil. Large (0.5-m depth by 0.3-m diam.) soil cores, collected from compacted wheel furrows in an irrigated cotton (Gossypium hirsutum L.) field, were subjected to three, six, or nine wet-dry cycles that simulated local flood irrigation practices. After each cycle, soil profiles were sampled for clod bulk density, image analysis of soil structure, and evapotranspiration. Generally, all crops improved soil structure over the initial field condition but lablab and mung bean gave improvements to greater depths and more rapidly than wheat and sorghum. Mung bean and lablab caused up to a threefold increase in clod porosity in the 0.1- to 0.4-m soil layer after only three wet-dry cycles, whereas sorghum required nine wet-dry cycles to increase clod porosity in only the 0.2- to 0.3-m layer, and wheat gave no improvement even after nine wet-dry cycles. Image analysis of soil structure showed that lablab and mung bean rapidly (by three wet-dry cycles) produced smaller peds with more interconnected pore space than wheat and sorghum. By nine wet-dry cycles, sorghum achieved deep cracking of the soil but the material between the cracks remained large and dense. Evapotranspiration was double under lablab and mung bean compared with wheat and sorghum. Our results indicate greater cycles of wetting and drying under lablab and mung bean than wheat and sorghum that have led to rapid repair of soil compaction.

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.

Experimental examination of solute transport by surface runoff from low-angle slopes

The removal of chemicals in solution, by overland flow from agricultural land has the potential to be a significant source of chemical loss from zero-till and surface mulched farming systems. The objective of this study was to determine the magnitude of solute loss by surface runoff from agricultural systems. Previous experiments have enhanced the understanding of the exchange process, but the initial soil conditions together with the tracer application method in these experiments have meant that in some cases the results have limited applicability to field situations. In this study, two different sets of experiments were carried out to determine the magnitude of solute loss by surface runoff. These experiments entailed the surface application of bromide to (1) field scale plots 18 m long by 2 m wide and (2) repacked soil cores 236 mm in diameter; followed by the application of simulated rainfall in both cases. The most substantial finding of the field experiments was that the quantities of solute in surface runoff varied greatly with soil type and structure (0.07-14.9% of the applied bromide). Also, on some soils, large quantities of tracer were measured in the surface runoff even after several hours of infiltration. The experiments on soil cores showed that soil structure plays an important role in the quantity of chemical that may be transported in the surface runoff. These field results showed that, in certain systems, solute movement by overland flow is an important transport mechanism, which should be considered when budgeting for chemical loss. (C) 2000 Elsevier Science B.V. All rights reserved.

Modelling solute transport in structured soils: Performance evaluation of the ADR and TRM models

The movement of chemicals through the soil to the groundwater or discharged to surface waters represents a degradation of these resources. In many cases, serious human and stock health implications are associated with this form of pollution. The chemicals of interest include nutrients, pesticides, salts, and industrial wastes. Recent studies have shown that current models and methods do not adequately describe the leaching of nutrients through soil, often underestimating the risk of groundwater contamination by surface-applied chemicals, and overestimating the concentration of resident solutes. This inaccuracy results primarily from ignoring soil structure and nonequilibrium between soil constituents, water, and solutes. A multiple sample percolation system (MSPS), consisting of 25 individual collection wells, was constructed to study the effects of localized soil heterogeneities on the transport of nutrients (NO3-, Cl-, PO43-) in the vadose zone of an agricultural soil predominantly dominated by clay. Very significant variations in drainage patterns across a small spatial scale were observed tone-way ANOVA, p < 0.001) indicating considerable heterogeneity in water flow patterns and nutrient leaching. Using data collected from the multiple sample percolation experiments, this paper compares the performance of two mathematical models for predicting solute transport, the advective-dispersion model with a reaction term (ADR), and a two-region preferential flow model (TRM) suitable for modelling nonequilibrium transport. These results have implications for modelling solute transport and predicting nutrient loading on a larger scale. (C) 2001 Elsevier Science Ltd. All rights reserved.

3D reconstruction and quantification of macropores using X-ray computed tomography and image analysis

Axial X-ray Computed tomography (CT) scanning provides a convenient means of recording the three-dimensional form of soil structure. The technique has been used for nearly two decades, but initial development has concentrated on qualitative description of images. More recently, increasing effort has been put into quantifying the geometry and topology of macropores likely to contribute to preferential now in soils. Here we describe a novel technique for tracing connected macropores in the CT scans. After object extraction, three-dimensional mathematical morphological filters are applied to quantify the reconstructed structure. These filters consist of sequences of so-called erosions and/or dilations of a 32-face structuring element to describe object distances and volumes of influence. The tracing and quantification methodologies were tested on a set of undisturbed soil cores collected in a Swiss pre-alpine meadow, where a new earthworm species (Aporrectodea nocturna) was accidentally introduced. Given the reduced number of samples analysed in this study, the results presented only illustrate the potential of the method to reconstruct and quantify macropores. Our results suggest that the introduction of the new species induced very limited chance to the soil structured for example, no difference in total macropore length or mean diameter was observed. However. in the zone colonised by, the new species. individual macropores tended to have a longer average length. be more vertical and be further apart at some depth. Overall, the approach proved well suited to the analysis of the three-dimensional architecture of macropores. It provides a framework for the analysis of complex structures, which are less satisfactorily observed and described using 2D imaging. (C) 2002 Elsevier Science B.V. All rights reserved.

Simultaneous X-ray imaging of plant root growth and water uptake in thin-slab systems

Direct and simultaneous observation of root growth and plant water uptake is difficult because soils are opaque. X-ray imaging techniques such as projection radiography or Computer Tomography (CT) offer a partial alternative to such limitations. Nevertheless, there is a trade-off between resolution, large field-of-view and 3-dimensionality: With the current state of the technology, it is possible to have any two. In this study, we used X-ray transmission through thin-slab systems to monitor transient saturation fields that develop around roots as plants grow. Although restricted to 2-dimensions, this approach offers a large field-of-view together with high spatial and dynamic resolutions. To illustrate the potential of this technology, we grew peas in 1 cm thick containers filled with soil and imaged them at regular intervals. The dynamics of both the root growth and the water content field that developed around the roots could be conveniently monitored. Compared to other techniques such as X-ray CT, our system is relatively inexpensive and easy to implement. It can potentially be applied to study many agronomic problems, such as issues related to the impact of soil constraints (physical, chemical or biological) on root development.

Slaking, dispersion, and crust formation

Soil structure is generally defined as the arrangement, orientation, and organization of the primary particles of sand, silt, and clay into compound aggregates, which exhibit properties that are unequal to the properties of a mass of nonaggregated material with a similar texture.6 Therefore the nature of soil structure is that it conveys specific properties to the soil and any alteration, i.e., breakdown or structural development, to the soil structural units will affect the physical properties of the soil. The aggregation and organization of the soil particles tend to form a hierarchical order4, 5 where the lower orders tend to have higher densities and greater internal strength than the higher orders. A schematic diagram of the hierarchical nature of soil structural elements in a clay soil is given in Fig. 1.4 Clay particles tend to form domains (packets of parallel clay sheets, generally consisting of 5-7 sheets), in turn several domains form clusters, followed by several orders of clusters, micro- and macroaggregates. The hierarchical nature implies that the destruction of a lower order will result in the destruction of all higher hierarchical orders. An example is the dispersion of sodic clay domains which results in the destruction of all higher orders, resulting in a dense soil with low hydraulic conductivity. Hence the clay domains are the fundamental building blocks of the soil and its integrity may determine the soil's physical properties and behavior.

Rice growth and post-rice mungbean in relation to two puddling intensities under glasshouse conditions

The effect of soil puddling on growth of lowland rice (Oryza sativa) and post-rice mungbean (Vigna radiata) was investigated using mini rice beds under controlled glasshouse conditions. Each mini rice bed was approximately 1 m(3) in size. Three different soil types were used: a well-drained, permeable loam; a hardsetting, structurally unstable silty loam; and a medium clay. Rice yields were reduced by low puddling compared with high puddling intensity on the loam but not affected on the heavier textured soils (silty loam and clay). Yield of mungbean was reduced on highly puddle, structurally unstable soil, indicating that puddling should be reduced on structurally unstable soils. Under glasshouse condition where crop establishment was not a limiting factor and plant available water in 0.65 m of soil was 100 mm, mungbean yields of >1 t/ha were achieved. However, under conditions where subsoil water reserves were depleted for the production of vegetative biomass during initial optimal growing condition, grain yield remained well below 1 t/ha.

Conventional detection methodology is limiting our ability to understand the roles and functions of fine roots

We lack a thorough conceptual and functional understanding of fine roots. Studies that have focused on estimating the quantity of fine roots provide evidence that they dominate overall plant root length. We need a standard procedure to quantify root length/biomass that takes proper account of fine roots. Here we investigated the extent to which root length/biomass may be underestimated using conventional methodology, and examined the technical reasons that could explain such underestimation. Our discussion is based on original X-ray-based measurements and on a literature review spanning more than six decades. We present evidence that root-length recovery depends strongly on the observation scale/spatial resolution at which measurements are carried out; and that observation scales/resolutions adequate for fine root detection have an adverse impact on the processing times required to obtain precise estimates. We conclude that fine roots are the major component of root systems of most (if not all) annual and perennial plants. Hence plant root systems could be much longer, and probably include more biomass, than is widely accepted.

A two-dimensional fuzzy random model of soil pore structure

A new conceptual model for soil pore-solid structure is formalized. Soil pore-solid structure is proposed to comprise spatially abutting elements each with a value which is its membership to the fuzzy set ''pore,'' termed porosity. These values have a range between zero (all solid) and unity (all pore). Images are used to represent structures in which the elements are pixels and the value of each is a porosity. Two-dimensional random fields are generated by allocating each pixel a porosity by independently sampling a statistical distribution. These random fields are reorganized into other pore-solid structural types by selecting parent points which have a specified local region of influence. Pixels of larger or smaller porosity are aggregated about the parent points and within the region of interest by controlled swapping of pixels in the image. This creates local regions of homogeneity within the random field. This is similar to the process known as simulated annealing. The resulting structures are characterized using one-and two-dimensional variograms and functions describing their connectivity. A variety of examples of structures created by the model is presented and compared. Extension to three dimensions presents no theoretical difficulties and is currently under development.

Bubbly Flow Structure in Hydraulic Jump

In an open channel, a hydraulic jump is the rapid transition from super- to sub-critical flow associated with strong turbulence and air bubble entrainment in the mixing layer. New experiments were performed at relatively large Reynolds numbers using phase-detection probes. Some new signal analysis provided characteristic air-water time and length scales of the vortical structures advecting the air bubbles in the developing shear flow. An analysis of the longitudinal air-water flow structure suggested little bubble clustering in the mixing layer, although an interparticle arrival time analysis showed some preferential bubble clustering for small bubbles with chord times below 3 ms. Correlation analyses yielded longitudinal air-water time scales Txx*V1/d1 of about 0.8 in average. The transverse integral length scale Z/d1 of the eddies advecting entrained bubbles was typically between 0.25 and 0.4, irrespective of the inflow conditions within the range of the investigations. Overall the findings highlighted the complicated nature of the air-water flow