Co-kriging when the soil and ancillary data are not co-located

Data such as digitized aerial photographs, electrical conductivity and yield are intensive and relatively inexpensive to obtain compared with collecting soil data by sampling. If such ancillary data are co-regionalized with the soil data they should be suitable for co-kriging. The latter requires that information for both variables is co-located at several locations; this is rarely so for soil and ancillary data. To solve this problem, we have derived values for the ancillary variable at the soil sampling locations by averaging the values within a radius of 15 m, taking the nearest-neighbour value, kriging over 5 m blocks, and punctual kriging. The cross-variograms from these data with clay content and also the pseudo cross-variogram were used to co-krige to validation points and the root mean squared errors (RMSEs) were calculated. In general, the data averaged within 15m and the punctually kriged values resulted in more accurate predictions.

Comparing sampling needs for variograms of soil properties computed by the method of moments and residual maximum likelihood

It has been generally accepted that the method of moments (MoM) variogram, which has been widely applied in soil science, requires about 100 sites at an appropriate interval apart to describe the variation adequately. This sample size is often larger than can be afforded for soil surveys of agricultural fields or contaminated sites. Furthermore, it might be a much larger sample size than is needed where the scale of variation is large. A possible alternative in such situations is the residual maximum likelihood (REML) variogram because fewer data appear to be required. The REML method is parametric and is considered reliable where there is trend in the data because it is based on generalized increments that filter trend out and only the covariance parameters are estimated. Previous research has suggested that fewer data are needed to compute a reliable variogram using a maximum likelihood approach such as REML, however, the results can vary according to the nature of the spatial variation. There remain issues to examine: how many fewer data can be used, how should the sampling sites be distributed over the site of interest, and how do different degrees of spatial variation affect the data requirements? The soil of four field sites of different size, physiography, parent material and soil type was sampled intensively, and MoM and REML variograms were calculated for clay content. The data were then sub-sampled to give different sample sizes and distributions of sites and the variograms were computed again. The model parameters for the sets of variograms for each site were used for cross-validation. Predictions based on REML variograms were generally more accurate than those from MoM variograms with fewer than 100 sampling sites. A sample size of around 50 sites at an appropriate distance apart, possibly determined from variograms of ancillary data, appears adequate to compute REML variograms for kriging soil properties for precision agriculture and contaminated sites. (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.

Exploring the enzymatic landscape: distribution and kinetics of hydrolytic enzymes in soil particle-size fractions

The location of extracellular enzymes within the soil architecture and their association with the various soil components affects their catalytic potential. A soil fractionation study was carried out to investigate: (a) the distribution of a range of hydrolytic enzymes involved in C, N and P transformations, (b) the effect of the location on their respective kinetics, (c) the effect of long-term N fertilizer management on enzyme distribution and kinetic parameters. Soil (silty clay loam) from grassland which had received 0 or 200 kg N ha(-1) yr(-1) was fractionated, and four particle-size fractions (> 200, 200-63, 63-2 and 0. 1-2 mum) were obtained by a combination of wet-sieving and centrifugation, after low-energy ultrasonication. All fractions were assayed for four carbohydrases (beta-cellobiohydrolase, N-acetyl-beta-glucosammidase, beta-glucosidase and beta-xylosidase), acid phosphatase and leucine-aminopeptidase using a microplate fluorimetric assay based on MUB-substrates. Enzyme kinetics (V-max and K-m) were estimated in three particle-size fractions and the unfractionated soil. The results showed that not all particle-size fractions were equally enzymatically active and that the distribution of enzymes between fractions depended on the enzyme. Carbohydrases predominated in the coarser fractions while phosphatase and leucine-aminopeptidase were predominant in the clay-size fraction. The Michaelis constant (K.) varied among fractions, indicating that the association of the same enzyme with different particle-size fractions affected its substrate affinity. The same values of Km were found in the same fractions from the soil under two contrasting fertilizer management regimes, indicating that the Michaelis constant was unaffected by soil changes caused by N fertilizer management. (C) 2004 Elsevier Ltd. All rights reserved.

Adsorption of water vapor by bare soil in an olive grove in southern Spain

Data for water vapor adsorption and evaporation are presented for a bare soil (sandy loam, clay content 15%) in a southern Spanish olive grove. Water losses and gains were measured using eight high-precision minilysimeters, placed around an olive tree, which had been irrigated until the soil reached field capacity (similar to 0.22 m(3) m(-3)). They were subsequently left to dry for 10 days. A pair of lysimeters was situated at each of the main points of the compass (N, E, S, W), at a distance of 1 m (the inner set of lysimeters; ILS) and 2 m (the outer set of lysimeters; OLS), respectively, from the tree trunk. Distinct periods of moisture loss (evaporation) and moisture gain (vapor adsorption) could be distinguished for each day. Vapor adsorption often started just after noon and generally lasted until the (early) evening. Values of up to 0.7 mm of adsorbed water per day were measured. Adsorption was generally largest for the OLS (up to 100% more on a daily basis), and increased during the dry down. This was mainly the result of lower OLS surface soil moisture contents (period-average absolute difference similar to 0.005 m(3) m(-3)), as illustrated using various analyses employing a set of micrometeorological equations describing the exchange of water vapor between bare soil and the atmosphere. These analyses also showed that the amount of water vapor adsorbed by soils is very sensitive to changes in atmospheric forcing and surface variables. The use of empirical equations to estimate vapor adsorption is therefore not recommended.

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.

The effect of soil phosphorus on particulate phosphorus in land runoff

Accumulation of surplus phosphorus (P) in the soil and the resulting increased transport of P in land runoff contribute to freshwater eutrophication. The effects of increasing soil P (19–194 mg Olsen-P (OP) kg−1) on the concentrations of particulate P (PP), and sorption properties (Qmax, k and EPCo) of suspended solids (SS) in overland flow from 15 unreplicated field plots established on a dispersive arable soil were measured over three monitoring periods under natural rainfall. Concentrations of PP in plot runoff increased linearly at a rate of 2.6 μg litre−1 per mg OP kg−1 of soil, but this rate was approximately 50% of the rate of increase in dissolved P (< 0.45 μm). Concentrations of SS in runoff were similar across all plots and contained a greater P sorption capacity (mean + 57%) than the soil because of enrichment with fine silt and clay (0.45–20 μm). As soil P increased, the P enrichment ratio of the SS declined exponentially, and the values of P saturation (Psat; 15–42%) and equilibrium P concentration (EPCo; 0.7–5.5 mg litre−1) in the SS fell within narrower ranges compared with the soils (6–74% and 0.1–10 mg litre−1, respectively). When OP was < 100 mg kg−1, Psat and EPCo values in the SS were smaller than those in the soil and vice-versa, suggesting that eroding particles from soils with both average and high P fertility would release P on entering the local (Rosemaund) stream. Increasing soil OP from average to high P fertility increased the P content of the SS by approximately 10%, but had no significant (P > 0.05) effect on the Psat, or EPCo, of the SS. Management options to reduce soil P status as a means of reducing P losses in land runoff and minimizing eutrophication risk may therefore have more limited effect than is currently assumed in catchment management.

Soil disturbed using a strip tillage implement on a range of soil types and the effects on sugar beet establishment

The area of soil disturbed using a single tine is well documented. However, modern strip tillage implements using a tine and disc design have not been assessed in the UK or in mainland Europe. Using a strip tillage implement has potential benefits for European agriculture where economic returns and sustainability are key issues. Using a strip tillage system a narrow zone is cultivated leaving most of the straw residue on the soil surface. Small field plot experiments were undertaken on three soil types and the operating parameters of forward speed, tine depth and tine design were investigated together with measurements of seedbed tilth and crop emergence. The type of tine used was found to be the primary factor in achieving the required volume of disturbance within a narrow zone whilst maintaining an area of undisturbed soil with straw residue on the surface. The winged tine produced greater disturbance at a given depth compared with the knife tine. Increasing forward speed did not consistently increase the volume of disturbance. In a sandy clay loam the tilth created and emergence of sugar beet by strip tillage and ploughing were similar but on a sandy loam the strip tillage treatments generally gave a finer tilth but poorer emergence particularly at greater working depth.

Carbon fractions associated with silt-size particles in surface and subsurface soil horizons

Information on the distribution and behavior of C fractions in soil particle sizes is crucial for understanding C dynamics in soil. At present little is known about the behavior of the C associated with silt-size particles. We quantified the concentrations, distribution, and enrichment of total C (TC), readily oxidizable C (ROC), hotwater- extractable C (HWC), and cold-water-extractable C (CWC) fractions in coarse (63–20-mm), medium (20–6.3-mm), and fine (6.3–2-mm) silt-size subfractions and in coarse (2000–250 mm) and fine (250–63 mm) sand and clay (<2-mm) soil fractions isolated from bulk soil (<2 mm), and 2- to 4-mm aggregate-size fraction of surface (0–25 cm) and subsurface (25–55 cm) soils under different land uses. All measured C fractions varied significantly across all soil particle-size fractions. The highest C concentrations were associated with the <20-mm soil fractions and peaked in the medium (20–6.3-mm) and fine (6.3–2-mm) silt subfractions in most treatments. Carbon enrichment ratios (ERC) revealed the dual behavior of the C fractions associated with the medium silt-size fraction, demonstrating the simultaneous enrichment of TC and ROC, and the depletion of HWC and CWC fractions. The medium silt (20–6.3-mm) subfraction was identified in this study as a zone where the associated C fractions exhibit transitory qualities. Our results show that investigating subfractions within the silt-size particle fraction provides better understanding of the behavior of C fractions in this soil fraction.

Determination of the gas diffusion coefficient of a peat grassland soil

Peatland habitats are important carbon stocks that also have the potential to be significant sources of greenhouse gases, particularly when subject to changes such as artificial drainage and application of fertilizer. Models aiming to estimate greenhouse gas release from peatlands require an accurate estimate of the diffusion coefficient of gas transport through soil (Ds). The availability of specific measurements for peatland soils is currently limited. This study measured Ds for a peat soil with an overlying clay horizon and compared values with those from widely available models. The Ds value of a sandy loam reference soil was measured for comparison. Using the Currie (1960) method, Ds was measured between an air-filled porosity (ϵ) range of 0 and 0.5 cm3 cm−3. Values of Ds for the peat cores ranged between 3.2 × 10−4 and 4.4 × 10−3 m2 hour−1, for loamy clay cores between 0 and 4.7 × 10−3 m2 hour−1 and for the sandy reference soil they were between 5.4 × 10−4 and 3.4 × 10−3 m2 hour−1. The agreement of measured and modelled values of relative diffusivity (Ds/D0, with D0 the diffusion coefficient through free air) varied with soil type; however, the Campbell (1985) model provided the best replication of measured values for all soils. This research therefore suggests that the use of the Campbell model in the absence of accurately measured Ds and porosity values for a study soil would be appropriate. Future research into methods to reduce shrinkage of peat during measurement and therefore allow measurement of Ds for a greater range of ϵ would be beneficial.

A new data set of soil mineralogy for dust-cycle modeling

The mineralogy of airborne dust affects the impact of dust particles on direct and indirect radiative forcing, on atmospheric chemistry and on biogeochemical cycling. It is determined partly by the mineralogy of the dust-source regions and partly by size-dependent fractionation during erosion and transport. Here we present a data set that characterizes the clay and silt-sized fractions of global soil units in terms of the abundance of 12 minerals that are important for dust–climate interactions: quartz, feldspars, illite, smectite, kaolinite, chlorite, vermiculite, mica, calcite, gypsum, hematite and goethite. The basic mineralogical information is derived from the literature, and is then expanded following explicit rules, in order to characterize as many soil units as possible. We present three alternative realizations of the mineralogical maps, taking the uncertainties in the mineralogical data into account. We examine the implications of the new database for calculations of the single scattering albedo of airborne dust and thus for dust radiative forcing.

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.

Designing a sampling scheme to reveal correlations between weeds and soil properties at multiple spatial scales

Weeds tend to aggregate in patches within fields and there is evidence that this is partly owing to variation in soil properties. Because the processes driving soil heterogeneity operate at different scales, the strength of the relationships between soil properties and weed density would also be expected to be scale-dependent. Quantifying these effects of scale on weed patch dynamics is essential to guide the design of discrete sampling protocols for mapping weed distribution. We have developed a general method that uses novel within-field nested sampling and residual maximum likelihood (REML) estimation to explore scale-dependent relationships between weeds and soil properties. We have validated the method using a case study of Alopecurus myosuroides in winter wheat. Using REML, we partitioned the variance and covariance into scale-specific components and estimated the correlations between the weed counts and soil properties at each scale. We used variograms to quantify the spatial structure in the data and to map variables by kriging. Our methodology successfully captured the effect of scale on a number of edaphic drivers of weed patchiness. The overall Pearson correlations between A. myosuroides and soil organic matter and clay content were weak and masked the stronger correlations at >50 m. Knowing how the variance was partitioned across the spatial scales we optimized the sampling design to focus sampling effort at those scales that contributed most to the total variance. The methods have the potential to guide patch spraying of weeds by identifying areas of the field that are vulnerable to weed establishment.

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.