The analysis of ranked observations of soil structure using indicator geostatistics

Structure is an important physical feature of the soil that is associated with water movement, the soil atmosphere, microorganism activity and nutrient uptake. A soil without any obvious organisation of its components is known as apedal and this state can have marked effects on several soil processes. Accurate maps of topsoil and subsoil structure are desirable for a wide range of models that aim to predict erosion, solute transport, or flow of water through the soil. Also such maps would be useful to precision farmers when deciding how to apply nutrients and pesticides in a site-specific way, and to target subsoiling and soil structure stabilization procedures. Typically, soil structure is inferred from bulk density or penetrometer resistance measurements and more recently from soil resistivity and conductivity surveys. To measure the former is both time-consuming and costly, whereas observations made by the latter methods can be made automatically and swiftly using a vehicle-mounted penetrometer or resistivity and conductivity sensors. The results of each of these methods, however, are affected by other soil properties, in particular moisture content at the time of sampling, texture, and the presence of stones. Traditional methods of observing soil structure identify the type of ped and its degree of development. Methods of ranking such observations from good to poor for different soil textures have been developed. Indicator variograms can be computed for each category or rank of structure and these can be summed to give the sum of indicator variograms (SIV). Observations of the topsoil and subsoil structure were made at four field sites where the soil had developed on different parent materials. The observations were ranked by four methods and indicator and the sum of indicator variograms were computed and modelled for each method of ranking. The individual indicators were then kriged with the parameters of the appropriate indicator variogram model to map the probability of encountering soil with the structure represented by that indicator. The model parameters of the SIVs for each ranking system were used with the data to krige the soil structure classes, and the results are compared with those for the individual indicators. The relations between maps of soil structure and selected wavebands from aerial photographs are examined as basis for planning surveys of soil structure. (C) 2007 Elsevier B.V. All rights reserved.

Problems with determining the particle size distribution of chalk soil and some of their implications

Particle size distribution (psd) is one of the most important features of the soil because it affects many of its other properties, and it determines how soil should be managed. To understand the properties of chalk soil, psd analyses should be based on the original material (including carbonates), and not just the acid-resistant fraction. Laser-based methods rather than traditional sedimentation methods are being used increasingly to determine particle size to reduce the cost of analysis. We give an overview of both approaches and the problems associated with them for analyzing the psd of chalk soil. In particular, we show that it is not appropriate to use the widely adopted 8 pm boundary between the clay and silt size fractions for samples determined by laser to estimate proportions of these size fractions that are equivalent to those based on sedimentation. We present data from field and national-scale surveys of soil derived from chalk in England. Results from both types of survey showed that laser methods tend to over-estimate the clay-size fraction compared to sedimentation for the 8 mu m clay/silt boundary, and we suggest reasons for this. For soil derived from chalk, either the sedimentation methods need to be modified or it would be more appropriate to use a 4 pm threshold as an interim solution for laser methods. Correlations between the proportions of sand- and clay-sized fractions, and other properties such as organic matter and volumetric water content, were the opposite of what one would expect for soil dominated by silicate minerals. For water content, this appeared to be due to the predominance of porous, chalk fragments in the sand-sized fraction rather than quartz grains, and the abundance of fine (<2 mu m) calcite crystals rather than phyllosilicates in the clay-sized fraction. This was confirmed by scanning electron microscope (SEM) analyses. "Of all the rocks with which 1 am acquainted, there is none whose formation seems to tax the ingenuity of theorists so severely, as the chalk, in whatever respect we may think fit to consider it". Thomas Allan, FRS Edinburgh 1823, Transactions of the Royal Society of Edinburgh. (C) 2009 Natural Environment Research Council (NERC) Published by Elsevier B.V. All rights reserved.

Spatial analysis of the error in a model of soil nitrogen

Models of the dynamics of nitrogen in soil (soil-N) can be used to aid the fertilizer management of a crop. The predictions of soil-N models can be validated by comparison with observed data. Validation generally involves calculating non-spatial statistics of the observations and predictions, such as their means, their mean squared-difference, and their correlation. However, when the model predictions are spatially distributed across a landscape the model requires validation with spatial statistics. There are three reasons for this: (i) the model may be more or less successful at reproducing the variance of the observations at different spatial scales; (ii) the correlation of the predictions with the observations may be different at different spatial scales; (iii) the spatial pattern of model error may be informative. In this study we used a model, parameterized with spatially variable input information about the soil, to predict the mineral-N content of soil in an arable field, and compared the results with observed data. We validated the performance of the N model spatially with a linear mixed model of the observations and model predictions, estimated by residual maximum likelihood. This novel approach allowed us to describe the joint variation of the observations and predictions as: (i) independent random variation that occurred at a fine spatial scale; (ii) correlated random variation that occurred at a coarse spatial scale; (iii) systematic variation associated with a spatial trend. The linear mixed model revealed that, in general, the performance of the N model changed depending on the spatial scale of interest. At the scales associated with random variation, the N model underestimated the variance of the observations, and the predictions were correlated poorly with the observations. At the scale of the trend, the predictions and observations shared a common surface. The spatial pattern of the error of the N model suggested that the observations were affected by the local soil condition, but this was not accounted for by the N model. In summary, the N model would be well-suited to field-scale management of soil nitrogen, but suited poorly to management at finer spatial scales. This information was not apparent with a non-spatial validation. (c),2007 Elsevier B.V. All rights reserved.

Plant roots release phospholipid surfactants that modify the physical and chemical properties of soil

Plant root mucilages contain powerful surfactants that will alter the interaction of soil solids with water and ions, and the rates of microbial processes. The lipid composition of maize, lupin and wheat root mucilages was analysed by thin layer chromatography and gas chromatography-mass spectrometry. A commercially available phosphatidylcholine (lecithin), chemically similar to the phospholipid surfactants identified in the mucilages, was then used to evaluate its effects on selected soil properties. The lipids found in the mucilages were principally phosphatidylcholines, composed mainly of saturated fatty acids, in contrast to the lipids extracted from root tissues. In soil at low tension, lecithin reduced the water content at any particular tension by as much as 10 and 50% in soil and acid-washed sand, respectively. Lecithin decreased the amount of phosphate adsorption in soil and increased the phosphate concentration in solution by 10%. The surfactant also reduced net rates of ammonium consumption and nitrate production in soil. These experiments provide the first evidence we are aware of that plant-released surfactants will significantly modify the biophysical environment of the rhizosphere.

Mineral and geochemical characterization of a leptic aluandic soil and a thapto aluandic-ferralsol developed on trachytes in Mount Bambouto (Cameroon volcanic line)

Mineral and geochemical investigations were carried out on soil samples and fresh rock (trachytes) from two selected soil profiles (TM profile on leptic aluandic soils and TL profile on thapto aluandic-ferralsols) from Mount Bambouto to better understand geochemical processes and mineral paragenesis involved in the development of soils in this environment. In TM profile, the hydrated halloysites and goethite occur in the weathered saprolite boulders of BC horizon while dehydrated halloysite, gibbsite and goethite dominate the soils matrices of BC and A horizons. In TL profile, the dehydrated halloysites and goethite are the most abundant secondary minerals in the weathered saprolites of C and BC horizons while gibbsite, hematite and kaolinite occur in the soil matrices of BC, B and A horizons. The highest gibbsite content is in the platy nodules of B horizon. In both soil profiles, organo-metal complexes (most likely of AI and Fe) are present in the surface A horizon. Geochemically, between the fresh rock and the weathered saprolites in both soils, SiO2, K2O, CaO, Na2O and MgO contents decrease strongly while Fe2O3 and Al2O3 tend to accumulate. The molar ratio of SiO2/Al2O3 (Ki) and the sum of Ca, Mg, K and Na ions (TRB) also decreases abruptly between fresh rocks and the weathered saprolites, but increases significantly at the soil surface. The TM profile shows intense Al enrichment whereas the TL profile highlights enrichment in both AI and Fe as the weathering progresses upwards. Both soil profiles are enriched in Ni, Cu, Ba and Co and depleted in U, Th, Ta, Hf, Y, Sr, Pb, Zr and Zn relative to fresh rock. They also show a relatively low fractionation of the rare earth elements (REE: La, Nd, Sm, Eu, Tb, Yb and Lu), except for Ce which tends to be enriched in soils compared to CI chondrite. All these results give evidence of intense hydrolysis at soil deep in Mount Bambouto resulting in the formation of halloysite which progressively transforms into gibbsite and/or dehydrated halloysite. At the soil surface, the prominent pedogenetic process refers to andosolization with formation of organo-metal complexes. In TL profile, the presence of kaolinite in soil matrices BC and B horizons is consistent with ferralitization at soil deep. In conclusion, soil forming processes in Mount Bambouto are strongly influenced by local climate: (i) in the upper mountain (>2000 m), the fresh, misty and humid climate favors andosolization; whereas (ii) in the middle lands (1700-2000 m) with a relatively dry climate, both andosolization at the soil surface and ferralitization at soil deep act together. (C) 2009 Elsevier B.V. All rights reserved.

Remote estimation of thermal inertia and soil heat flux for bare soil

A method is presented which allows thermal inertia (the soil heat capacity times the square root of the soil thermal diffusivity, C(h)rootD(h)), to be estimated remotely from micrometeorological observations. The method uses the drop in surface temperature, T-s, between sunset and sunrise, and the average night-time net radiation during that period, for clear, still nights. A Fourier series analysis was applied to analyse the time series of T-s . The Fourier series constants, together with the remote estimate of thermal inertia, were used in an analytical expression to calculate diurnal estimates of the soil heat flux, G. These remote estimates of C(h)rootD(h) and G compared well with values derived from in situ sensors. The remote and in situ estimates of C(h)rootD(h) both correlated well with topsoil moisture content. This method potentially allows area-average estimates of thermal inertia and soil heat flux to be derived from remote sensing, e.g. METEOSAT Second Generation, where the area is determined by the sensor's height and viewing angle. (C) 2003 Elsevier B.V. All rights reserved.

Field trials to assess the uptake of arsenic by vegetables from contaminated soils and soil remediation with iron oxides

The uptake of arsenic (As) by plants from contaminated soils presents a health hazard that may affect the use of agricultural and former industrial land. Methods for limiting the hazard are desirable. A proposed remediation treatment comprises the precipitation of iron (Fe) oxides in the contaminated soil by adding ferrous sulfate and lime. The effects on As bioavailability were assessed using a range of vegetable crops grown in the field. Four UK locations were used, where soil was contaminated by As from different sources. At the most contaminated site, a clay loam containing a mean of 748 mg As kg(-1) soil, beetroot, calabrese, cauliflower, lettuce, potato, radish and spinach were grown. For all crops except spinach, ferrous sulfate treatment caused a significant reduction in the bioavailability of As in some part of the crop. Application of ferrous sulfate in solution, providing 0.2% Fe oxides in the soil (0-10 cm), reduced As uptake by a mean of 22%. Solid ferrous sulfate was applied to give concentrations of 0.5% and 1% Fe oxides: the 0.5% concentration reduced As uptake by a mean of 32% and the 1% concentration gave no significant additional benefit. On a sandy loam containing 65 mg As kg(-1) soil, there was tentative evidence that ferrous sulfate treatment up to 2% Fe oxides caused a significant reduction in lettuce As, but calabrese did not respond. At the other two sites, the effects of ferrous sulfate treatment were not significant, but the uptake of soil As was low in treated and untreated soils. Differences between sites in the bioavailable fraction of soil As may be related to the soil texture or the source of As. The highest bioavailability was found on the soil which had been contaminated by aerial deposition and had a high sand content. (C) 2003 Elsevier Science B.V. All rights reserved.

Moisture can be the dominant environmental parameter governing cadaver decomposition in soil

Forensic taphonomy involves the use of decomposition to estimate postmortem interval (PMI) or locate clandestine graves. Yet, cadaver decomposition remains poorly understood, particularly following burial in soil. Presently, we do not know how most edaphic and environmental parameters, including soil moisture, influence the breakdown of cadavers following burial and alter the processes that are used to estimate PMI and locate clandestine graves. To address this, we buried juvenile rat (Rattus rattus) cadavers (∼18 g wet weight) in three contrasting soils from tropical savanna ecosystems located in Pallarenda (sand), Wambiana (medium clay), or Yabulu (loamy sand), Queensland, Australia. These soils were sieved (2 mm), weighed (500 g dry weight), calibrated to a matric potential of -0.01 megapascals (MPa), -0.05 MPa, or -0.3 MPa (wettest to driest) and incubated at 22 °C. Measurements of cadaver decomposition included cadaver mass loss, carbon dioxide-carbon (CO2-C) evolution, microbial biomass carbon (MBC), protease activity, phosphodiesterase activity, ninhydrin-reactive nitrogen (NRN) and soil pH. Cadaver burial resulted in a significant increase in CO2-C evolution, MBC, enzyme activities, NRN and soil pH. Cadaver decomposition in loamy sand and sandy soil was greater at lower matric potentials (wetter soil). However, optimal matric potential for cadaver decomposition in medium clay was exceeded, which resulted in a slower rate of cadaver decomposition in the wettest soil. Slower cadaver decomposition was also observed at high matric potential (-0.3 MPa). Furthermore, wet sandy soil was associated with greater cadaver decomposition than wet fine-textured soil. We conclude that gravesoil moisture content can modify the relationship between temperature and cadaver decomposition and that soil microorganisms can play a significant role in cadaver breakdown. We also conclude that soil NRN is a more reliable indicator of gravesoil than soil pH.

Microbial decomposition of skeletal muscle tissue (Ovis aries) in a sandy loam soil at different temperatures

A laboratory experiment was conducted to determine the effect of temperature (2, 12, 22 °C) on the rate of aerobic decomposition of skeletal muscle tissue (Ovis aries) in a sandy loam soil incubated for a period of 42 days. Measurements of decomposition processes included skeletal muscle tissue mass loss, carbon dioxide (CO2) evolution, microbial biomass, soil pH, skeletal muscle tissue carbon (C) and nitrogen (N) content and the calculation of metabolic quotient (qCO2). Incubation temperature and skeletal muscle tissue quality had a significant effect on all of the measured process rates with 2 °C usually much lower than 12 and 22 °C. Cumulative CO2 evolution at 2, 12 and 22 °C equaled 252, 619 and 905 mg CO2, respectively. A significant correlation (P<0.001) was detected between cumulative CO2 evolution and tissue mass loss at all temperatures. Q10s for mass loss and CO2 evolution, which ranged from 1.19 to 3.95, were higher for the lower temperature range (Q10(2– 12 °C)>Q10(12–22 °C)) in the Ovis samples and lower for the low temperature range (Q10(2–12 °C)C)) in the control samples. Metabolic quotient and the positive relationship between skeletal muscle tissue mass loss and cumulative CO2 evolution suggest that tissue decomposition was most efficient at 2 °C. These phenomena may be due to lower microbial catabolic requirements at lower temperature.

A laboratory incubation method for determining the rate of microbiological degradation of skeletal muscle tissue in soil

A controlled laboratory experiment is described, in principle and practice, which can be used for the of determination the rate of tissue decomposition in soil. By way of example, an experiment was conducted to determine the effect of temperature (12°C, 22°C) on the aerobic decomposition of skeletal muscle tissue (Organic Texel × Suffolk lamb (Ovis aries)) in a sandy loam soil. Measurements of decomposition processes included muscle tissue mass loss, microbial CO2 respiration, and muscle tissue carbon (C) and nitrogen (N). Muscle tissue mass loss at 22°C always was greater than at 12°C (p < 0.001). Microbial respiration was greater in samples incubated at 22°C for the initial 21 days of burial (p < 0.01). All buried muscle tissue samples demonstrated changes in C and N content at the end of the experiment. A significant correlation (p < 0.001) was demonstrated between the loss of muscle tissue-derived C (C1) and microbially-respired C (Cm) demonstrating CO2 respiration may be used to predict mass loss and hence biodegradation. In this experiment Q10 (12°C - 22°C) = 2.0. This method is recommended as a useful tool in determining the effect of environmental variables on the rate of decomposition of various tissues and associated materials.

Soil quality assessment based on carbon stratification index in different olive grove management practices in Mediterranean areas

In Mediterranean areas, conventional tillage increases soil organic matter losses, reduces soil quality, and contributes to climate change due to increased CO2 emissions. CO2 sequestration rates in soil may be enhanced by appropriate agricultural soil management and increasing soil organic matter content. This study analyzes the stratification ratio (SR) index of soil organic carbon (SOC), nitrogen (N) and C:N ratio under different management practices in an olive grove (OG) in Mediterranean areas (Andalusia, southern Spain). Management practices considered in this study are conventional tillage (CT) and no tillage (NT). In the first case, CT treatments included addition of alperujo (A) and olive leaves (L). A control plot with no addition of olive mill waste was considered (CP). In the second case, NT treatments included addition of chipped pruned branches (NT1) and chipped pruned branches and weeds (NT2). The SRs of SOC increased with depth for all treatments. The SR of SOC was always higher in NT compared to CT treatments, with the highest SR of SOC observed under NT2. The SR of N increased with depth in all cases, ranging between 0.89 (L-SR1) and 39.11 (L-SR3 and L-SR4).The SR of C:N ratio was characterized by low values, ranging from 0.08 (L-SR3) to 1.58 (NT1-SR2) and generally showing higher values in SR1 and SR2 compared to those obtained in SR3 and SR4. This study has evaluated several limitations to the SR index such as the fact that it is descriptive but does not analyze the behavior of the variable over time. In addition, basing the assessment of soil quality on a single variable could lead to an oversimplification of the assessment. Some of these limitations were experienced in the assessment of L, where SR1 of SOC was the lowest of the studied soils. In this case, the higher content in the second depth interval compared to the first was caused by the intrinsic characteristics of this soil's formation process rather than by degradation. Despite the limitations obtained SRs demonstrate that NT with the addition of organic material improves soil quality.

Effect of air-drying pre-treatment on the characterization of forest soil carbon pools

Archived soils could represent a valuable resource for the spatio-temporal inventory of soil carbon stability. However, archived soils are usually air-dried before storage and the impact of a drying pretreatment on physically and chemically-defined C fractions has not yet been fully assessed. Through the comparison of field-moist and corresponding air-dried (at 25oC for 2 weeks) forest soil samples, we examined the effect of air-drying on: a) the quantity and the quality of cold- (CWEC) and hot-water (HWEC) extractable C and b) the concentration of C in physically isolated fractions (free- and intra-aggregate light and organo-mineral). Soil samples were collected from the organic (O) and mineral (A and B) horizons of three different forest soils from southeastern England: (i) Cambisol under Pine (Pinus nigra); (ii) Cambisol under Beech (Fagus sylvatica) and (iii) Gleysol under oak (Quercus robur). CWEC concentrations for dry samples were up to 2 times greater than for corresponding field moist samples and had significantly (p < 0.001) higher phenolic content. However, the effect of drying pretreatment on HWEC, its phenolic content was not significant (p > 0.05) for most samples. Dried soils had significantly (p < 0.001) higher concentrations of free light-C while having lower concentrations of intra-aggregate-C when compared to moist samples (p < 0.001). However, fine silt and clay fractions were not significantly affected by the drying pretreatment (p=0.789). Therefore, based on the results obtained from gleysol and cambisol forest soils studied here, C contents in hot-water extractions and fine particle size physical fractions (< 25µm) seem to be robust measurements for evaluating C fractions in dried stored forest soils. Further soil types should be tested to evaluate the wider generality of these findings.

Organic contaminant content and physico-chemical characteristics of waste materials recycled in agriculture

A range of wastes representative of materials currently applied, or with future potential to be applied, to agricultural land in the UK as fertilisers and soil improvers or used as animal bedding in livestock production, were investigated. In addition to full physico-chemical characterization, the materials were analysed for a suite of priority organic contaminants. In general, contaminants were present at relatively low concentrations. For example, polychlorinated dibenzo-p-dioxins/dibenzofurans and polychlorinated biphenyls in biosolids and compost-like-outputs (CLOs) were, in most cases, between 5-50 times lower than proposed and implemented European limit values for biosolids or composts applied to agricultural land. However, the technical basis for these limits may need to be re-evaluated. Polybrominated, and mixed halogenated, dibenzo-p-dioxins/dibenzofurans are not currently considered in risk assessments of dioxins and dioxin-like chemicals, but were detected in the biosolids and compost-like-outputs and their potential contribution to the overall toxic equivalency will be assessed. Other, ‘emerging’ contaminants such as perfluoralkyl compounds (PFCs) and organophosphate flame retardants were detected in several of the waste materials, and their potential significance is discussed. The study is part of a wider research programme that will provide evidence to improve confidence in the use of waste-derived materials in agriculture and establish guidelines to protect the food chain where necessary.