Exploring the spatial relations between soil physical properties and apparent electrical conductivity

Relations between the apparent electrical conductivity of the soil (ECa) and top- and sub-soil physical properties were examined for two arable fields in southern England (Crowmarsh Battle Farms and the Yattendon Estate). The spatial variation of ECa and the soil properties was explored geostatistically. The variogram ranges showed that ECa varied on a similar spatial scale to many of the soil physical properties in both fields. Several features in the map of kriged predictions of ECa were also evident in maps of the soil properties. In addition, the correlation coefficients showed a strong relation between ECa and several soil properties. A moving correlation analysis enabled differences in the relations between ECa and the soil properties to be examined within the fields. The results indicated that relations were inconsistent; they were stronger in some areas than others. A regression of ECa on the principal component scores of the leading components for both fields showed that the first two components accounted for a large proportion of the variance in ECa, whereas the others accounted for little or none. For Crowmarsh topsoil sand and clay, loss on ignition and volumetric water measured in the autumn had large correlations on the first component, and for Yattendon they were large for topsoil sand and clay, and autumn and spring volumetric water. The cross-variograms suggested strong coregionalization between ECa and several soil physical properties; in particular subsoil sand and silt at Crowmarsh, and subsoil sand and clay at Yattendon. The structural correlations from the linear model of coregionalization confirmed the strength of the relations between ECa and the subsoil properties. Nevertheless, no one property was consistently important for both fields. Although a map of ECa can indicate the general patterns of spatial variation in the soil, it is not a substitute for information on soil properties obtained by sampling and analysing the soil. Nevertheless, it could be used to guide further sampling. (c) 2005 Elsevier B.V. All rights reserved.

Fate of prions in soil: Longevity and migration of recPrP in soil columns

Research into transmissible spongiform encephalopathy (TSE) diseases has become a high priority worldwide in recent years yet remarkably little is known about the behaviour of TSE infectivity in the environment. The resilience and stability of prion proteins could lead to soils becoming a potential reservoir of TSE infectivity as a result of contamination from activities such as infected carcass burial or the dispersion of effluents from slaughter houses, or by contamination of pastures by infected animals, (e.g. scrapie in sheep). Knowledge of the fate of prion proteins in soils, and associated physico-chemical conditions which favour migration, can be used to help prevent re-infection of animals through grazing, to protect watercourses and develop good management practices. In two consecutive experiments of 9 and 6 months, the migration of recombinant ovine PrP (recPrP) in soil columns was followed under contrasting levels of microbial activity (normal versus reduced), under varying regimes of soil water content and redox potential, and in two different soil types (loamy sand and clay loam). At each analysis time (1, 3, 6 or 9 months), in both soil types, full-length recPrP was detected in the original contaminated layer, indicating the resilience and stability of recPrP under varied soil conditions, even in the presence of active soil microbial populations. Evidence of protein migration was found in every soil column at the earliest analysis time (1 or 3 months), but was restricted to a maximum distance of 1 cm, indicative of limited initial mobility in soils followed by strong adsorption over the following days to weeks. The survival of recPrP in the soil over a period of at least 9 months was demonstrated. In this study, recPrP was used as an indicator for potential TSE infectivity, although infectivity tests should be carried out before conclusions can be drawn regarding the infection risk posed by prions in soil. However, it has been demonstrated that soil is likely to act as a significant barrier to the dispersion of contaminated material at storage or burial sites. (c) 2007 Elsevier Ltd. All rights reserved.

Estimating detailed soil water profile records from point measurements

Temporal and spatial patterns of soil water content affect many soil processes including evaporation, infiltration, ground water recharge, erosion and vegetation distribution. This paper describes the analysis of a soil moisture dataset comprising a combination of continuous time series of measurements at a few depths and locations, and occasional roving measurements at a large number of depths and locations. The objectives of the paper are: (i) to develop a technique for combining continuous measurements of soil water contents at a limited number of depths within a soil profile with occasional measurements at a large number of depths, to enable accurate estimation of the soil moisture vertical pattern and the integrated profile water content; and (ii) to estimate time series of soil moisture content at locations where there are just occasional soil water measurements available and some continuous records from nearby locations. The vertical interpolation technique presented here can strongly reduce errors in the estimation of profile soil water and its changes with time. On the other hand, the temporal interpolation technique is tested for different sampling strategies in space and time, and the errors generated in each case are compared.

Modelling stream and soil water nitrate dynamics during experimentally increased nitrogen deposition in a coniferous forest catchment at Gardsjon, Sweden

Increased atmospheric deposition of inorganic nitrogen (N) may lead to increased leaching of nitrate (NO3-) to surface waters. The mechanisms responsible for, and controls on, this leaching are matters of debate. An experimental N addition has been conducted at Gardsjon, Sweden to determine the magnitude and identify the mechanisms of N leaching from forested catchments within the EU funded project NITREX. The ability of INCA-N, a simple process-based model of catchment N dynamics, to simulate catchment-scale inorganic N dynamics in soil and stream water during the course of the experimental addition is evaluated. Simulations were performed for 1990-2002. Experimental N addition began in 1991. INCA-N was able to successfully reproduce stream and soil water dynamics before and during the experiment. While INCA-N did not correctly simulate the lag between the start of N addition and NO 2 3 breakthrough, the model was able to simulate the state change resulting from increased N deposition. Sensitivity analysis showed that model behaviour was controlled primarily by parameters related to hydrology and vegetation dynamics and secondarily by in-soil processes.

Joined-up archaeology at Old Scatness, Shetland: thin section analysis of the site and hinterland

Excavations on the multi-period settlement at Old Scatness, Shetland have uncovered a number of Iron Age structures with compacted, floor-like layers. Thin section analysis was undertaken in order to investigate and compare the characteristics of these layers. The investigation also draws on earlier analyses of the Iron Age agricultural soil around the settlement and the midden deposits that accumulated within the settlement, to create a 'joined-up' analysis which considers the way material from the settlement was used and then recycled as fertiliser for the fields. Peat was collected from the nearby uplands and was used for fuel and possibly also for flooring. It is suggested that organic-rich floors from the structures were periodically removed and the material was spread onto the fields as fertilisers. More organic-rich material may have been used selectively for fertiliser, while the less organic peat ash was allowed to accumulate in middens. Several of the structures may have functioned as byres, which suggests a prehistoric plaggen system.

Measuring root traits in barley (Hordeum vulgare ssp vulgare and ssp spontaneum) seedlings using gel chambers, soil sacs and X-ray microtomography

Root characteristics of seedlings of five different barley genotypes were analysed in 2D using gel chambers, and in 3D using soil sacs that were destructively harvested and pots of soil that were assessed non-invasively using X-ray microtomography. After 5 days, Chime produced the greatest number of root axes (similar to 6) and Mehola significantly less (similar to 4) in all growing methods. Total root length was longest in GSH01915 and shortest in Mehola for all methods, but both total length and average root diameter were significantly larger for plants grown in gel chambers than those grown in soil. The ranking of particular growth traits (root number, root angular spread) of plants grown in gel plates, soil sacs and X-ray pots was similar, but plants grown in the gel chambers had a different order of ranking for root length to the soil-grown plants. Analysis of angles in soil-grown plants showed that Tadmore had the most even spread of individual roots and Chime had a propensity for non-uniform distribution and root clumping. The roots of Mehola were less well spread than the barley cultivars supporting the suggestion that wild and landrace barleys tend to have a narrower angular spread than modern cultivars. The three dimensional analysis of root systems carried out in this study provides insights into the limitations of screening methods for root traits and useful data for modelling root architecture.

Predicted soil organic carbon stocks and changes in Kenya between 1990 and 2030

Under the United Nations Framework Convention on Climate Change (UNFCCC), Non-Annex 1 countries such as Kenya are obliged to report green house gas (GHG) emissions from all sources where possible, including those from soils as a result of changes in land use or land management. At present, the convention encourages countries to estimate emissions using the most advanced methods possible, given the country circumstances and resources. Estimates of soil organic carbon (SOC) stocks and changes were made for Kenya using the Global Environment Facility Soil Organic Carbon (GEFSOC) Modelling System. The tool conducts analysis using three methods: (1) the Century general ecosystem model; (2) the RothC soil C decomposition model; and (3) the Intergovernmental Panel on Climate Change (IPCC) method for assessing soil C at regional scales. The required datasets included: land use history, monthly mean precipitation, monthly mean minimum and maximum temperatures for all the agro-climatic zones of Kenya and historical vegetation cover. Soil C stocks of 1.4-2.0 Pg (0-20 cm), compared well with a Soil and Terrain (SOTER) based approach that estimated similar to .8-2.0 Pg (0-30 cm). In 1990 48% of the country had SOC stocks of < 18 t C ha(-1) and 20% of the country had SOC stocks of 18-30 t C ha(-1), whereas in 2000 56% of the country had SOC stocks of < 18 t C ha(-1) and 31% of the country had SOC stocks of 18-30 t C ha(-1). Conversion of natural vegetation to annual crops led to the greatest soil C losses. Simulations suggest that soil C losses remain substantial throughout the modelling period of 1990-2030. All three methods involved in the GEFSOC System estimated that there would be a net loss of soil C between 2000 and 2030 in Kenya. The decline was more marked with RothC than with Century or the IPCC method. In non-hydric soils the SOC change rates were more pronounced in high sandy soils compared to high clay soils in most land use systems. (C) 2007 Elsevier B.V. All rights reserved.

Variograms of ancillary data to aid sampling for soil surveys

To provide reliable estimates for mapping soil properties for precision agriculture requires intensive sampling and costly laboratory analyses. If the spatial structure of ancillary data, such as yield, digital information from aerial photographs, and soil electrical conductivity (EC) measurements, relates to that of soil properties they could be used to guide the sampling intensity for soil surveys. Variograins of permanent soil properties at two study sites on different parent materials were compared with each other and with those for ancillary data. The ranges of spatial dependence identified by the variograms of both sets of properties are of similar orders of magnitude for each study site, Maps of the ancillary data appear to show similar patterns of variation and these seem to relate to those of the permanent properties of the soil. Correlation analysis has confirmed these relations. Maps of kriged estimates from sub-sampled data and the original variograrns showed that the main patterns of variation were preserved when a sampling interval of less than half the average variogram range of ancillary data was used. Digital data from aerial photographs for different years and EC appear to show a more consistent relation with the soil properties than does yield. Aerial photographs, in particular those of bare soil, seem to be the most useful ancillary data and they are often cheaper to obtain than yield and EC data.

Determining nugget: sill ratios of standardized variograms from aerial photographs to krige sparse soil data

Maps of kriged soil properties for precision agriculture are often based on a variogram estimated from too few data because the costs of sampling and analysis are often prohibitive. If the variogram has been computed by the usual method of moments, it is likely to be unstable when there are fewer than 100 data. The scale of variation in soil properties should be investigated prior to sampling by computing a variogram from ancillary data, such as an aerial photograph of the bare soil. If the sampling interval suggested by this is large in relation to the size of the field there will be too few data to estimate a reliable variogram for kriging. Standardized variograms from aerial photographs can be used with standardized soil data that are sparse, provided the data are spatially structured and the nugget:sill ratio is similar to that of a reliable variogram of the property. The problem remains of how to set this ratio in the absence of an accurate variogram. Several methods of estimating the nugget:sill ratio for selected soil properties are proposed and evaluated. Standardized variograms with nugget:sill ratios set by these methods are more similar to those computed from intensive soil data than are variograms computed from sparse soil data. The results of cross-validation and mapping show that the standardized variograms provide more accurate estimates, and preserve the main patterns of variation better than those computed from sparse data.

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.

Sampling the soil in long-term forest plots: The implications of spatial variation

Long-term monitoring of forest soils as part of a pan-European network to detect environmental change depends on an accurate determination of the mean of the soil properties at each monitoring event. Forest soil is known to be very variable spatially, however. A study was undertaken to explore and quantify this variability at three forest monitoring plots in Britain. Detailed soil sampling was carried out, and the data from the chemical analyses were analysed by classical statistics and geostatistics. An analysis of variance showed that there were no consistent effects from the sample sites in relation to the position of the trees. The variogram analysis showed that there was spatial dependence at each site for several variables and some varied in an apparently periodic way. An optimal sampling analysis based on the multivariate variogram for each site suggested that a bulked sample from 36 cores would reduce error to an acceptable level. Future sampling should be designed so that it neither targets nor avoids trees and disturbed ground. This can be achieved best by using a stratified random sampling design.

Sensitivity analysis of a catchment-scale nitrogen model

There are now considerable expectations that semi-distributed models are useful tools for supporting catchment water quality management. However, insufficient attention has been given to evaluating the uncertainties inherent to this type of model, especially those associated with the spatial disaggregation of the catchment. The Integrated Nitrogen in Catchments model (INCA) is subjected to an extensive regionalised sensitivity analysis in application to the River Kennet, part of the groundwater-dominated upper Thames catchment, UK The main results are: (1) model output was generally insensitive to land-phase parameters, very sensitive to groundwater parameters, including initial conditions, and significantly sensitive to in-river parameters; (2) INCA was able to produce good fits simultaneously to the available flow, nitrate and ammonium in-river data sets; (3) representing parameters as heterogeneous over the catchment (206 calibrated parameters) rather than homogeneous (24 calibrated parameters) produced a significant improvement in fit to nitrate but no significant improvement to flow and caused a deterioration in ammonium performance; (4) the analysis indicated that calibrating the flow-related parameters first, then calibrating the remaining parameters (as opposed to calibrating all parameters together) was not a sensible strategy in this case; (5) even the parameters to which the model output was most sensitive suffered from high uncertainty due to spatial inconsistencies in the estimated optimum values, parameter equifinality and the sampling error associated with the calibration method; (6) soil and groundwater nutrient and flow data are needed to reduce. uncertainty in initial conditions, residence times and nitrogen transformation parameters, and long-term historic data are needed so that key responses to changes in land-use management can be assimilated. The results indicate the general, difficulty of reconciling the questions which catchment nutrient models are expected to answer with typically limited data sets and limited knowledge about suitable model structures. The results demonstrate the importance of analysing semi-distributed model uncertainties prior to model application, and illustrate the value and limitations of using Monte Carlo-based methods for doing so. (c) 2005 Elsevier B.V. All rights reserved.

Moving towards a more mechanistic approach in the determination of soil heat flux from remote measurements - II. Diurnal shape of soil heat flux

For vegetated surfaces, calculation of soil heat flux, G, with the Exact or Analytical method requires a harmonic analysis of below-canopy soil surface temperature, to obtain the shape of the diurnal course of G. When determining G with remote sensing methods, only composite (vegetation plus soil) radiometric brightness temperature is available. This paper presents a simple equation that relates the sum of the harmonic terms derived for the composite radiometric surface temperature to that of belowcanopy soil surface temperature. The thermal inertia, Gamma(,) for which a simple equation has been presented in a companion paper, paper I, is used to set the magnitude of G. To assess the success of the method proposed in this paper for the estimation of the diurnal shape of G, a comparison was made between 'remote' and in situ calculated values from described field sites. This indicated that the proposed method was suitable for the estimation of the shape of G for a variety of vegetation types and densities. The approach outlined in paper I, to obtain Gamma, was then combined with the estimated harmonic terms to predict estimates of G, which were compared to values predicted by empirical remote methods found in the literature. This indicated that the method proposed in the combination of papers I and II gave reliable estimates of G, which, in comparison to the other methods, resulted in more realistic predictions for vegetated surfaces. This set of equations can also be used for bare and sparsely vegetated soils, making it a universally applicable method. (C) 2007 Elsevier B.V. All rights reserved.

Effects of metal pollution on soil macroinvertebrate burrow systems

A reduction in the numbers of macroinvertebrates present in soil may have a negative effect on soil structure, infiltration rates, and gas exchanges. Soil pollution by metal is known to have a detrimental effect on soil macrofauna. The aim of the present study was to evaluate (1) direct and indirect effects of soil pollution on soil macroinvertebrate bioturbation and (2) effects of the two macroinvertebrate communities found in a polluted and a nonpolluted area (one supposed sensitive, the other tolerant to metals) on burrow systems parameters. Macroinvertebrate porosity was studied using X-ray tomography. Three-dimensional reconstructions and characterisation of the burrow system were obtained using image analysis. Results showed that metal pollution principally affected the spatial distribution of macropores (more macropores were found near the soil surface) and the shape of the burrow system (branching rate was higher in the polluted soil), whereas soil macroinvertebrate composition principally affects burrow density parameters (the number of burrows was higher for the sensitive macroinvertebrate community).

Effect of mycorrhizae and pH change at the root-soil interface on phosphorus uptake by Sorghum using a rhizocylinder technique

Sorghum (Sorghum bicolor) was grown for 40 days in. rhizocylinder (a growth container which permitted access to rh zosphere and nonrhizosphere soil), in two soils of low P status. Soils were fertilized with different rates of ammonium and nitrate and supplemented with 40 mg phosphorus (P) kg(-1) and inoculated with either Glomus mosseae (Nicol. and Gerd.) or nonmycorrhizal root inoculum.. N-serve (2 mg kg(-1)) was added to prevent nitrification. At harvest, soil from around the roots was collected at distances of 0-5, 5-10, and 10-20 mm from the root core which was 35 mm diameter. Sorghum plants, with and without mycorrhiza, grew larger with NH4+ than with NO3- application. After measuring soil pH, 4 3 suspensions of the same sample were titrated against 0.01 M HCl or 0.01 M NaOH until soil pH reached the nonplanted pH level. The acid or base requirement for each sample was calculated as mmol H+ or OFF kg(-1) soil. The magnitude of liberated acid or base depended on the form and rate of nitrogen and soil type. When the plant root was either uninfected or infected with mycorrhiza., soil pH changes extended up to 5 mm from the root core surface. In both soils, ammonium as an N source resulted in lower soil pH than nitrate. Mycorrhizal (VAM) inoculation did not enhance this difference. In mycorrhizal inoculated soil, P depletion extended tip to 20 mm from the root surface. In non-VAM inoculated soil P depletion extended up to 10 mm from the root surface and remained unchanged at greater distances. In the mycorrhizal inoculated soils, the contribution of the 0-5 mm soil zone to P uptake was greater than the core soil, which reflects the hyphal contribution to P supply. Nitrogen (N) applications that caused acidification increased P uptake because of increased demand; there is no direct evidence that the increased uptake was due to acidity increasing the solubility of P although this may have been a minor effect.