An easily implemented agro-hydrological procedure with dynamic root simulation for water transfer in the crop-soil system: validation and application

Models for water transfer in the crop-soil system are key components of agro-hydrological models for irrigation, fertilizer and pesticide practices. Many of the hydrological models for water transfer in the crop-soil system are either too approximate due to oversimplified algorithms or employ complex numerical schemes. In this paper we developed a simple and sufficiently accurate algorithm which can be easily adopted in agro-hydrological models for the simulation of water dynamics. We used a dual crop coefficient approach proposed by the FAO for estimating potential evaporation and transpiration, and a dynamic model for calculating relative root length distribution on a daily basis. In a small time step of 0.001 d, we implemented algorithms separately for actual evaporation, root water uptake and soil water content redistribution by decoupling these processes. The Richards equation describing soil water movement was solved using an integration strategy over the soil layers instead of complex numerical schemes. This drastically simplified the procedures of modeling soil water and led to much shorter computer codes. The validity of the proposed model was tested against data from field experiments on two contrasting soils cropped with wheat. Good agreement was achieved between measurement and simulation of soil water content in various depths collected at intervals during crop growth. This indicates that the model is satisfactory in simulating water transfer in the crop-soil system, and therefore can reliably be adopted in agro-hydrological models. Finally we demonstrated how the developed model could be used to study the effect of changes in the environment such as lowering the groundwater table caused by the construction of a motorway on crop transpiration. (c) 2009 Elsevier B.V. All rights reserved.

Speleothems as sensitive recorders of volcanic eruptions - the Bronze Age Minoan eruption recorded in a stalagmite from Turkey

Mounting evidence exists that variations in sulphur content in stalagmites are closely linked to changes in volcanic or anthropogenic atmospheric sulphur. The strong dependency of sulphur on soil pH and ecosystem storage, however, can result in a delay of several years to decades in the registration of volcanic eruptions and anthropogenic emissions by stalagmites. Here we present synchrotron-radiation based trace element analysis performed on a precisely-dated section of a stalagmite from Sofular Cave in Northern Turkey. As this section covers the time interval of the intensively studied Minoan volcanic eruption between 1600 and 1650 BC, we can test whether this vigorous eruption can be traced in a stalagmite. Of all measured trace elements, only bromine shows a clear short-lived peak at 1621±251621±25 BC, whereas sulphur and molybdenum show peaks later at 1617±251617±25 and 1589±251589±25 respectively. We suggest that all trace element peaks are related to the Minoan eruption, whereas the observed phasing of bromine, molybdenum and sulphur is related to differences in their retention rates in the soil above Sofular Cave. For the first time, we can show that bromine appears to be an ideal volcanic tracer in stalagmites, as it is a prominent volatile component in volcanic eruptions, can be easily leached in soils and rapidly transferred from the atmosphere through the soil and bedrock into the cave and stalagmite respectively. Highly resolved oxygen and carbon isotope profiles indicate that the Minoan eruption had no detectable climatic and environmental impact in Northern Turkey.

Are sulfurous soil amendments (S0, Fe(II)SO4, Fe(III)SO4) an effective tool in the restoration of heathland and acidic grassland after four decades of rock phosphate Fertilization?

This paper deals with the complex issue of reversing long-term improvements of fertility in soils derived from heathlands and acidic grasslands using sulfur-based amendments. The experiment was conducted on a former heathland and acid grassland in the U.K. that was heavily fertilized and limed with rock phosphate, chalk, and marl. The experimental work had three aims. First, to determine whether sulfurous soil amendments are able to lower pH to a level suitable for heathland and acidic grassland re-creation (approximately 3 pH units). Second, to determine what effect the soil amendments have on the available pool of some basic cations and some potentially toxic acidic cations that may affect the plant community. Third, to determine whether the addition of Fe to the soil system would sequester PO4− ions that might be liberated from rock phosphate by the experimental treatments. The application of S0 and Fe(II)SO4− to the soil was able to reduce pH. However, only the highest S0 treatment (2,000 kg/ha S) lowered pH sufficiently for heathland restoration purposes but effectively so. Where pH was lowered, basic cations were lost from the exchangeable pool and replaced by acidic cations. Where Fe was added to the soil, there was no evidence of PO4− sequestration from soil test data (Olsen P), but sequestration was apparent because of lower foliar P in the grass sward. The ability of the forb Rumex acetosella to apparently detoxify Al3+, prevalent in acidified soils, appeared to give it a competitive advantage over other less tolerant species. We would anticipate further changes in plant community structure through time, driven by Al3+ toxicity, leading to the competitive exclusion of less tolerant species. This, we suggest, is a key abiotic driver in the restoration of biotic (acidic plant) communities.

The relationship between diffuse spectral reflectance of the soil and its cation exchange capacity is scale-dependent

Diffuse reflectance spectroscopy (DRS) is increasingly being used to predict numerous soil physical, chemical and biochemical properties. However, soil properties and processes vary at different scales and, as a result, relationships between soil properties often depend on scale. In this paper we report on how the relationship between one such property, cation exchange capacity (CEC), and the DRS of the soil depends on spatial scale. We show this by means of a nested analysis of covariance of soils sampled on a balanced nested design in a 16 km × 16 km area in eastern England. We used principal components analysis on the DRS to obtain a reduced number of variables while retaining key variation. The first principal component accounted for 99.8% of the total variance, the second for 0.14%. Nested analysis of the variation in the CEC and the two principal components showed that the substantial variance components are at the > 2000-m scale. This is probably the result of differences in soil composition due to parent material. We then developed a model to predict CEC from the DRS and used partial least squares (PLS) regression do to so. Leave-one-out cross-validation results suggested a reasonable predictive capability (R2 = 0.71 and RMSE = 0.048 molc kg− 1). However, the results from the independent validation were not as good, with R2 = 0.27, RMSE = 0.056 molc kg− 1 and an overall correlation of 0.52. This would indicate that DRS may not be useful for predictions of CEC. When we applied the analysis of covariance between predicted and observed we found significant scale-dependent correlations at scales of 50 and 500 m (0.82 and 0.73 respectively). DRS measurements can therefore be useful to predict CEC if predictions are required, for example, at the field scale (50 m). This study illustrates that the relationship between DRS and soil properties is scale-dependent and that this scale dependency has important consequences for prediction of soil properties from DRS data

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.

Modelling the movement and survival of the root-feeding clover weevil, Sitona lepidus, in the root-zone of white clover

White clover (Trifolium repens) is an important pasture legume but is often difficult to sustain in a mixed sward because, among other things, of the damage to roots caused by the soil-dwelling larval stages of S. lepidus. Locating the root nodules on the white clover roots is crucial for the survival of the newly hatched larvae. This paper presents a numerical model to simulate the movement of newly hatched S. lepidus larvae towards the root nodules, guided by a chemical signal released by the nodules. The model is based on the diffusion-chemotaxis equation. Experimental observations showed that the average speed of the larvae remained approximately constant, so the diffusion-chernotaxis model was modified so that the larvae respond only to the gradient direction of the chemical signal but not its magnitude. An individual-based lattice Boltzmann method was used to simulate the movement of individual larvae, and the parameters required for the model were estimated from the measurement of larval movement towards nodules in soil scanned using X-ray microtomography. The model was used to investigate the effects of nodule density, the rate of release of chemical signal, the sensitivity of the larvae to the signal, and the random foraging of the larvae on the movement and subsequent survival of the larvae. The simulations showed that the most significant factors for larval survival were nodule density and the sensitivity of the larvae to the signal. The dependence of larval survival rate on nodule density was well fitted by the Michealis-Menten kinetics. (c) 2005 Elsevier B.V All rights reserved.

Genotypic variations in the accumulation of Cd, Cu, Pb and Zn exhibited by six commonly grown vegetables

Metal contaminants in garden and allotment soils could possibly affect human health through a variety of pathways. This study focused on the potential pathway of consumption of vegetables grown on contaminated soil. Five cultivars each of six common vegetables were grown in a control and in a soil spiked with Cd, Cu, Pb and Zn. Highly significant differences in metal content were evident between cultivars of a number of vegetables for several of the contaminants. Carrot and pea cultivars exhibited significant differences in accumulated concentrations of Cd and Cu with carrot cultivars also exhibiting significant differences in Zn. Distinctive differences were also identified when comparing one vegetable to another, legumes (Leguminosae) tending to be low accumulators, root vegetables (Umbelliferae and Liliaceae) tending to be moderate accumulators and leafy vegetables (Compositae and Chenopodiaceae) being high accumulators. (c) 2006 Elsevier Ltd. All rights reserved.

Soil factors associated with zinc deficiency in crops and humans

Zinc deficiency is the most ubiquitous micronutrient deficiency problem in world crops. Zinc is essential for both plants and animals because it is a structural constituent and regulatory co-factor in enzymes and proteins involved in many biochemical pathways. Millions of hectares of cropland are affected by Zn deficiency and approximately one-third of the human population suffers from an inadequate intake of Zn. The main soil factors affecting the availability of Zn to plants are low total Zn contents, high pH, high calcite and organic matter contents and high concentrations of Na, Ca, Mg, bicarbonate and phosphate in the soil solution or in labile forms. Maize is the most susceptible cereal crop, but wheat grown on calcareous soils and lowland rice on flooded soils are also highly prone to Zn deficiency. Zinc fertilizers are used in the prevention of Zn deficiency and in the biofortification of cereal grains.

Evidence of heavy metal movement down the profile of a heavily-sludged soil

Two sites in central England where sewage sludge has been deposited for decades were studied to measure the heavy metal distribution in the soil profiles. The first site (S 1) was a field receiving heavy loads sludge from a nearby wastewater treatment plant, and the second (S2) was a farm applying 'normal' sludge rates of 8 t ha(-1) y(-1) of the same sludge. Soil samples were also taken by a near-by untreated control site. In S I the movement of heavy metals was significant even down to 80 cm depth compared to the control. In S2, the concentrations of lead (Pb) and zinc (Zn) and the organic matter content were higher than the control down to 20 cm, while nickel (Ni) moved significantly down to 80 cm. This underlies. the possibility that the metals bound onto organic surfaces moved along with organic matter down to that depth. The movement of metals in S2 points out the potential risks of applying sewage sludge for a long time.

The influence of mineral solubility and soil solution concentration on the toxicity of copper to Eisenia fetida Savigny

The OECD 14 d earthworm acute toxicity test was used to determine the toxicity of copper added as copper nitrate (Cu(NO3)(2)), copper sulphate (CuSO4) and malachite (Cu-2(OH)(2)(CO3)) to Eisenia fetida Savigny. Cu(NO3)(2), and CuSO4 were applied in both an aqueous (aq) and solid (s) form, Cu-2(OH)(2)(CO3) was added as a solid. Soil solution was extracted by centrifugation, and analysed for copper. Two extractants [0.01 M CaCl2 and 0.005 M diethylenetriminpentaacetic acid (DTPA)] were used as a proxy of the bioavailable copper fraction in the soil. For bulk soil copper content the calculated copper toxicity decreased in the order nitrate > sulphide > carbonate, the same order as decreasing solubility of the metal compounds. For Cu(NO3)(2) and CuSO4, the LC50s obtained were not significantly different when the compound was added in solution or solid form. There was a significant correlation between the soil solution copper concentration and the percentage earthworm mortality for all 3 copper compounds (P less than or equal to 0.05) indicating that the soil pore water copper concentration is important for determining copper availability and toxicity to E. fetida. In soil avoidance tests the earthworms avoided the soils treated with Cu(NO3)(2) (aq and s) and CuSO4 (aq and s), at all concentrations used (110-8750 mug Cu g(-1), and 600-8750 mug Cu g(-1) respectively). In soils treated with Cu-2(OH2)CO3, avoidance behaviour was exhibited at all concentrations greater than or equal to3500 mug Cu g(-1). There was no significant correlation between the copper extracted by either CaCl2 or DTPA and percentage mortality. These two extractants are therefore not useful indicators of copper availability and toxicity to E. fetida.

Spatial covariation of Azotobacter abundance and soil properties: A case study using the wavelet transform

The soil microflora is very heterogeneous in its spatial distribution. The origins of this heterogeneity and its significance for soil function are not well understood. A problem for understanding spatial variation better is the assumption of statistical stationarity that is made in most of the statistical methods used to assess it. These assumptions are made explicit in geostatistical methods that have been increasingly used by soil biologists in recent years. Geostatistical methods are powerful, particularly for local prediction, but they require the assumption that the variability of a property of interest is spatially uniform, which is not always plausible given what is known about the complexity of the soil microflora and the soil environment. We have used the wavelet transform, a relatively new innovation in mathematical analysis, to investigate the spatial variation of abundance of Azotobacter in the soil of a typical agricultural landscape. The wavelet transform entails no assumptions of stationarity and is well suited to the analysis of variables that show intermittent or transient features at different spatial scales. In this study, we computed cross-variograms of Azotobacter abundance with the pH, water content and loss on ignition of the soil. These revealed scale-dependent covariation in all cases. The wavelet transform also showed that the correlation of Azotobacter abundance with all three soil properties depended on spatial scale, the correlation generally increased with spatial scale and was only significantly different from zero at some scales. However, the wavelet analysis also allowed us to show how the correlation changed across the landscape. For example, at one scale Azotobacter abundance was strongly correlated with pH in part of the transect, and not with soil water content, but this was reversed elsewhere on the transect. The results show how scale-dependent variation of potentially limiting environmental factors can induce a complex spatial pattern of abundance in a soil organism. The geostatistical methods that we used here make assumptions that are not consistent with the spatial changes in the covariation of these properties that our wavelet analysis has shown. This suggests that the wavelet transform is a powerful tool for future investigation of the spatial structure and function of soil biota. (c) 2006 Elsevier Ltd. All rights reserved.

Evaluating the century C model using long-term fertilizer trials in the Indo-Gangetic Plains, India

The GEFSOC Project developed a system for estimating soil carbon (C) stocks and changes at the national and sub-national scale. As part of the development of the system, the Century ecosystem model was evaluated for its ability to simulate soil organic C (SOC) changes in environmental conditions in the Indo-Gangetic Plains, India (IGP). Two long-term fertilizer trials (LTFT), with all necessary parameters needed to run Century, were used for this purpose: a jute (Corchorus capsularis L.), rice (Oryza sativa L.) and wheat (Triticum aestivum L.) trial at Barrackpore, West Bengal, and a rice-wheat trial at Ludhiana, Punjab. The trials represent two contrasting climates of the IGP, viz. semi-arid, dry with mean annual rainfall (MAR) of < 800 mm and humid with > 1600 turn. Both trials involved several different treatments with different organic and inorganic fertilizer inputs. In general, the model tended to overestimate treatment effects by approximately 15%. At the semi-arid site, modelled data simulated actual data reasonably well for all treatments, with the control and chemical N + farm yard manure showing the best agreement (RMSE = 7). At the humid site, Century performed less well. This could have been due to a range of factors including site history. During the study, Century was calibrated to simulate crop yields for the two sites considered using data from across the Indian IGP. However, further adjustments may improve model performance at these sites and others in the IGP. The availability of more longterm experimental data sets (especially those involving flooded lowland rice and triple cropping systems from the IGP) for testing and validation is critical to the application of the model's predictive capabilities for this area of the Indian sub-continent. (C) 2007 Elsevier B.V. All rights reserved.

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.