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.

Effects of manure and fertilizer on grain yield, soil carbon and phosphorus in a 13-year field trial in semi-arid Kenya

Long-term indicators of soil fertility were assessed by measuring grain yield, soil organic carbon (SOC) and soil Olsen phosphorous for a P-deficient soil. In one set of treatments, goat manure was applied annually for 13 years at 0, 5 and 10 t ha(-1), and intercrops of sorghum/cowpea, millet/green gram and maize/pigeonpea were grown. Yield depended on rainfall and trends with time were not identifiable. Manure caused an upward trend in SOC, but 10 t ha(-1) manure did not give significantly more SOC than 5 t ha(-1). Only 10 t ha(-1) manure increased Olsen P. Measurements of both SOC and Olsen P are recommended. In another set of treatments, manure was applied for four years; the residual effect lasted another seven to eight years when assessed by yield, SOC and Olsen P Treatment with mineral fertilizers provided the same rates of N and P as 5 t hat manure and yields from manure and fertilizer were similar. Fertilizer increased Olsen P but not SOC. Management systems with occasional manure application and intermediate fertilizer applications should be assessed. Inputs and offtakes of C, N and P were measured for three years. Approximately 16, 25 and 11% of C, N and P respectively were stabilized into soil organic matter from 5 t ha(-1) a(-1) manure. The majority of organic P was fixed as soil inorganic P.

Survival, Pb-uptake and behaviour of three species of earthworm in Pb treated soils determined using an OECD-style toxicity test and a soil avoidance test

Mature (clitellate) Eisenia andrei Bouche (ultra epigeic), Lumbricus rubellus Hoffmeister (epigeic), and Aporrectodea caliginosa (Savigny) (endogeic) earthworms were placed in soils treated with Pb(NO3)(2) to have concentrations in the range 1000 to 10 000 mg Pb kg(-1). After 28 days LC50(-95%confidence limit) (+95%confidence limit) values were E. andrei 5824(-361)(+898) mg Pb kg(-1), L. rubellus 2867(-193)(+145) mg Pb kg(-1) and A. caliginosa 2747(-304)(+239) mg Pb kg(-1) and EC50s for weight change were E. andrei 2841(-68)(+150) Pb kg(-1), L. rubellus 1303(-201)(+204) mg Pb kg(-1) and A. caliginosa 1208(-206)(+212) Mg Pb kg(-1). At any given soil Pb concentration, Pb tissue concentrations after 28 days were the same for all three earthworm species. In a soil avoidance test there was no difference between the behaviour of the different species. The lower sensitivity to Pb exhibited by E. andrei is most likely due to physiological adaptations associated with the modes of life of the earthworms, and could have serious implications for the use of this earthworm as the species of choice in standard toxicological testing. (c) 2005 Elsevier Ltd. All rights reserved.

Interactions between earthworms and arsenic in the soil environment: a review

Chemical pollution of the environment has become a major source of concern. In particular, many studies have investigated the impact of pollution on biota in the environment. Studies on metalliferous contaminated mine spoil wastes have shown that some soil organisms have the capability to become resistant to metal/metalloid toxicity. Earthworms are known to inhabit arsenic-rich metalliferous soils and, due to their intimate contact with the soil, in both the solid and aqueous phases, are likely to accumulate contaminants present in mine spoil. Earthworms that inhabit metalliferous contaminated soils must have developed mechanisms of resistance to the toxins found in these soils. The mechanisms of resistance are not fully understood; they may involve physiological adaptation (acclimation) or be genetic. This review discusses the relationships between earthworms and arsenic-rich mine spoil wastes, looking critically at resistance and possible mechanisms of resistance, in relation to soil edaphic factors and possible trophic transfer routes. (C) 2003 Elsevier Science Ltd. 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.

Comparison of EDTA and EDDS as potential soil amendments for enhanced phytoextraction of heavy metals

Phytoextraction has been proposed as an alternative remediation technology for soils polluted with heavy metals or radionuclides, but is generally conceived as too slow working. Enhancing the accumulation of trace pollutants in harvestable plant tissues is a prerequisite for the technology to be practically applicable. The chelating aminopolycarboxylic acid, ethylene diamine tetraacetate (EDTA), has been found to enhance shoot accumulation of heavy metals. However, the use of EDTA in phytoextraction may not be suitable due to its high environmental persistence, which may lead to groundwater contamination. This paper aims to assess whether ethylene diamine disuccinate (EDDS), a biodegradable chelator, can be used for enhanced phytoextraction purposes. A laboratory experiment was conducted to examine mobilisation of Cd, Cu, Cr, Ni, Pb and Zn into the soil solution upon application of EDTA or EDDS. The longevity of the induced mobilisation was monitored for a period of 40 days after application. Estimated effect half lives ranged between 3.8 and 7.5 days for EDDS, depending on the applied dose. The minimum observed effect half life of EDTA was 36 days, while for the highest applied dose no decrease was observed throughout the 40 day period of the mobilisation experiment. Performance of EDTA and EDDS for phytoextraction was evaluated by application to Helianthus annuus. Two other potential chelators, known for their biodegradability in comparison to EDTA, were tested in the plant experiment: nitrilo acetic acid (NTA) and citric acid. Uptake of heavy metals was higher in EDDS-treated pots than in EDTA-treated pots. The effects were still considered insufficiently high to consider efficient remediation. This may be partly due to the choice of timing for application of the soil amendment. Fixing the time of application at an earlier point before harvest may yield better results. NTA and citric acid induced no significant effects on heavy metal uptake. (C) 2004 Elsevier Ltd. All rights reserved.

National and sub-national assessments of soil organic carbon stocks and changes: The GEFSOC modelling system

Soil organic carbon (SOC) plays a vital role in ecosystem function, determining soil fertility, water holding capacity and susceptibility to land degradation. In addition, SOC is related to atmospheric CO, levels with soils having the potential for C release or sequestration, depending on land use, land management and climate. The United Nations Convention on Climate Change and its Kyoto Protocol, and other United Nations Conventions to Combat Desertification and on Biodiversity all recognize the importance of SOC and point to the need for quantification of SOC stocks and changes. An understanding of SOC stocks and changes at the national and regional scale is necessary to further our understanding of the global C cycle, to assess the responses of terrestrial ecosystems to climate change and to aid policy makers in making land use/management decisions. Several studies have considered SOC stocks at the plot scale, but these are site specific and of limited value in making inferences about larger areas. Some studies have used empirical methods to estimate SOC stocks and changes at the regional scale, but such studies are limited in their ability to project future changes, and most have been carried out using temperate data sets. The computational method outlined by the Intergovernmental Panel on Climate Change (IPCC) has been used to estimate SOC stock changes at the regional scale in several studies, including a recent study considering five contrasting eco regions. This 'one step' approach fails to account for the dynamic manner in which SOC changes are likely to occur following changes in land use and land management. A dynamic modelling approach allows estimates to be made in a manner that accounts for the underlying processes leading to SOC change. Ecosystem models, designed for site scale applications can be linked to spatial databases, giving spatially explicit results that allow geographic areas of change in SOC stocks to be identified. Some studies have used variations on this approach to estimate SOC stock changes at the sub-national and national scale for areas of the USA and Europe and at the watershed scale for areas of Mexico and Cuba. However, a need remained for a national and regional scale, spatially explicit system that is generically applicable and can be applied to as wide a range of soil types, climates and land uses as possible. The Global Environment Facility Soil Organic Carbon (GEFSOC) Modelling System was developed in response to this need. The GEFSOC system allows estimates of SOC stocks and changes to be made for diverse conditions, providing essential information for countries wishing to take part in an emerging C market, and bringing us closer to an understanding of the future role of soils in the global C cycle. (C) 2007 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.

Does changing the taxonomical resolution alter the value of soil macroinvertebrates as bioindicators of metal pollution?

Ecological indicators are taxa that are affected by, and indicate effects of, anthropogenic environmental stress or disturbance on ecosystems. There is evidence that some species of soil macrofauna (i.e. diameter > 2 min) constitute valuable biological indicators of certain types of soil perturbations. This study aims to determine which level of taxonomic resolution, (species, family or ecological group) is the best to identify indicator of soil disturbance. Macrofauna were sampled in a set of sites encompassing different land-use systems (e.g. forests, pastures, crops) and different levels of pollution. Indicator taxa were sought using the IndVal index proposed by Dufrene and Legendre [Dufrene, M., Legendre, P., 1997. Species assemblages and indicator species: the need for a flexible asymetrical approach. Ecological Monographs 67, 345-3661. This approach is based on a hierarchical typology of sites. The index value changes along the typology and decreases (increases) for generalist (specialist) faunal units (species, families or ecological groups). Of the 327 morphospecies recorded, 19 were significantly associated with a site type or a group of sites (5.8%). Similarly, species were aggregated to form 59 families among which 17 (28.8%) displayed a significant indicator value. Gathering species into 28 broad ecological assemblages led to 14 indicator groups (50%). Beyond the simple proportion of units having significant association with a given level of the site typology, the proportion of specialist and generalist groups changed dramatically when the level of taxonomic resolution was altered. At the species level 84% of the indicator units were specialist, whereas this proportion decreased to 70 and 43% when families and ecological groups were considered. Because specialist groups are the most interesting type of indicators either in terms of conservation or for management purposes we come to the conclusion that the species level is the most accurate taxonomic level in bioindication studies although it requires a high amount of labour and operator knowledge and is time-consuming. (c) 2005 Published by Elsevier Ltd.

Hydrological controls on soil redox dynamics in a peat-based, restored wetland

Increasing areas of altered wetland are being restored by re-flooding the soil. Evidence in the literature indicates that this practice can induce the redox-mediated release of soil nutrients, thereby increasing the risk of diffuse water pollution. However, for the sake of improving wedand management decisions, there is a need for more detailed studies of the underlying relationship between the hydrological and redox dynamics that explain this risk; this is particularly the case in agricultural peatlands that are commonly targeted for the creation of lowland wet grassland. A 12-month field study was conducted to evaluate the relationship between hydrological fluctuations and soil redox potential (Eh) in a nutrient-rich peat field (32 g N kg(-1) and 1100 mg P kg(-1) in the surface 0-30 cm soil) that had been restored as lowland wet grassland from intensive arable production. Field tensiometers were installed at the 30-, 60- and 90-cm soil depths, and Pt electrodes at the 10-, 30-, 60- and 90-cm depths, for daily logging of soil water tension and Eh, respectively. The values for soil water tension displayed a strong negative relationship (P < 0.001) with monthly dip well observations of water table height. Calculations of soil water potential from the logged tension values were used, therefore, to provide a detailed profile of field water level and, together with precipitation data, explained some of the variation in Eh. For example, during the summer, alternating periods of aerobism (Eh > 330 mV) in the surface, 0-10 cm layer of peat coincided with intense precipitation events. Redox potential throughout the 30-100 cm profile also fluctuated seasonally; indeed, at all depths Eh displayed a strong, negative relationship (P < 0.001) with water table height over the 12-month study period. However, Eh throughout the 30-100 cm profile remained relatively low (< 230 mV), indicating permanently reduced conditions that are associated with denitrification and reductive dissolution of Fe-bound P. The implications of these processes in the N- and P-rich peat for wetland plant diversity and water quality are discussed. (c) 2006 Elsevier B.V. All rights reserved.

Phosphorus dynamics in the ditch system of a restored peat wetland

The restoration of wetlands as bird habitats often involves the maintenance of a fluctuating water regime by careful, localised ditch water management using pumps and sluices. However, there is evidence in the literature to Suggest that alternate flood/drainage cycles can accelerate nutrient cycling and transport within the soil and, therefore, pose a threat to water quality through the process of eutrophication. This study focused on the dynamics and losses of soil P in a recently re-wetted, eutrophic fen peat developed on alluvium ill South west England. During the 2-year Study (2001 and 2002), soil water tensiometry revealed that the field water table (fluctuating annually between +20 and 60 cm relative to ground level) was extensively influenced across most of the 8.4 ha field site by the management of the adjacent ditch water levels. This conservation-led, prescribed water balance was facilitated by the high hydraulic conductivity (1.1 x 10(-s) ms(-1)) of the lower (70-140 cm), degraded layer of peat. However, only during a 7-day period of water table drawdown by intermittent pump drainage, approximately 45 g ha(-1) of dissolved reactive P (DRP) entered the pumped ditch from the field via this degraded layer. Summer rainfall events >35 mm d(-1) also coincided with significant peaks ill ditch water P concentration (up to 200 mu g L-1 DRP). Even larger peaks (Up to 700 mu g L-1 DRP) Occurred With the annual onset of autumn reflooding. These episodic P loss events pose a serious potential threat to biological water quality. (C) 2009 Elsevier B.V. All rights reserved.

Factors influencing the dissolution of phosphate rock in a range of high P-fixing soils from Cameroon

A laboratory incubation experiment was conducted to evaluate the soil factors that influence the dissolution of two phosphate rocks (PRs) of different reactivity (Gafsa, GPR, reactive PR; and Togo-Hahotoe, HPR, low reactivity PR) in seven agricultural soils from Cameroon having variable phosphorus (P)- sorption capacities, organic carbon (C) contents, and exchangeable acidities. Ground PR was mixed with the soils at a rate of 500 mg P kg 21 soil and incubated at 30 degrees C for 85 days. Dissolution of the PRs was determined at various intervals using the Delta NaOH-P method ( the difference of the amount of P extracted by 0.5 M NaOH between the PR-treated soils and the control). Between 4 and 27% of HPR and 33 and 50% of GPR were dissolved in the soils. Calcium (Ca) saturation of cation exchange sites and proton supply strongly affected PR dissolution in these soils. Acid soils with pH-(H2O), < 5 (NKL, ODJ, NSM, MTF) dissolved more phosphate rock than those with pH-(H2O) > 5 (DSC, FGT, BAF). However, the lack of a sufficient Ca sink in the former constrained the dissolution of both PRs. The dissolution of GPR in the slightly acidic soils was limited by increase in Ca saturation and that of HPR was constrained by limited supply in protons. Generally, the dissolution of GPR was higher than that of HPR for each soil. The kinetics of dissolution of PR in the soils was best described by the power function equation P At B. More efficient use of PR in these soils can be achieved by raising the soil cation exchange capacity, thereby increasing the Ca sink size. This could be done by amending such soils with organic materials.