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 food web properties explain ecosystem services across European land use systems

Intensive land use reduces the diversity and abundance of many soil biota, with consequences for the processes that they govern and the ecosystem services that these processes underpin. Relationships between soil biota and ecosystem processes have mostly been found in laboratory experiments and rarely are found in the field. Here, we quantified, across four countries of contrasting climatic and soil conditions in Europe, how differences in soil food web composition resulting from land use systems (intensive wheat rotation, extensive rotation, and permanent grassland) influence the functioning of soils and the ecosystem services that they deliver. Intensive wheat rotation consistently reduced the biomass of all components of the soil food web across all countries. Soil food web properties strongly and consistently predicted processes of C and N cycling across land use systems and geographic locations, and they were a better predictor of these processes than land use. Processes of carbon loss increased with soil food web properties that correlated with soil C content, such as earthworm biomass and fungal/bacterial energy channel ratio, and were greatest in permanent grassland. In contrast, processes of N cycling were explained by soil food web properties independent of land use, such as arbuscular mycorrhizal fungi and bacterial channel biomass. Our quantification of the contribution of soil organisms to processes of C and N cycling across land use systems and geographic locations shows that soil biota need to be included in C and N cycling models and highlights the need to map and conserve soil biodiversity across the world.

Improved evaporative flux partitioning and carbon flux in the land surface model JULES: impact on the simulation of land surface processes in temperate Europe

The primary role of land surface models embedded in climate models is to partition surface available energy into upwards, radiative, sensible and latent heat fluxes. Partitioning of evapotranspiration, ET, is of fundamental importance: as a major component of the total surface latent heat flux, ET affects the simulated surface water balance, and related energy balance, and consequently the feedbacks with the atmosphere. In this context it is also crucial to credibly represent the CO2 exchange between ecosystems and their environment. In this study, JULES, the land surface model used in UK weather and climate models, has been evaluated for temperate Europe. Compared to eddy covariance flux measurements, the CO2 uptake by the ecosystem is underestimated and the ET overestimated. In addition, the contribution to ET from soil and intercepted water evaporation far outweighs the contribution of plant transpiration. To alleviate these biases, adaptations have been implemented in JULES, based on key literature references. These adaptations have improved the simulation of the spatio-temporal variability of the fluxes and the accuracy of the simulated GPP and ET, including its partitioning. This resulted in a shift of the seasonal soil moisture cycle. These adaptations are expected to increase the fidelity of climate simulations over Europe. Finally, the extreme summer of 2003 was used as evaluation benchmark for the use of the model in climate change studies. The improved model captures the impact of the 2003 drought on the carbon assimilation and the water use efficiency of the plants. It, however, underestimates the 2003 GPP anomalies. The simulations showed that a reduction of evaporation from the interception and soil reservoirs, albeit not of transpiration, largely explained the good correlation between the carbon and the water fluxes anomalies that was observed during 2003. This demonstrates the importance of being able to discriminate the response of individual component of the ET flux to environmental forcing.

Pre-Columbian raised-field agriculture and land use in the Bolivian Amazon

We present an integrated palaeoecological and archaeobotanical study of pre-Columbian raised-field agriculture in the Llanos de Moxos, a vast seasonally inundated forest–savanna mosaic in the Bolivian Amazon. Phytoliths from excavated raised-field soil units, together with pollen and charcoal in sediment cores from two oxbow lakes, were analysed to provide a history of land use and agriculture at the El Cerro raised-field site. The construction of raised fields involved the removal of savanna trees, and gallery forest was cleared from the area by AD 310. Despite the low fertility of Llanos de Moxos soils, we determined that pre-Columbian raised-field agriculture sufficiently improved soil conditions for maize cultivation. Fire was used as a common management practice until AD 1300, at which point, the land-use strategy shifted towards less frequent burning of savannas and raised fields. Alongside a reduction in the use of fire, sweet potato cultivation and the exploitation of Inga fruits formed part of a mixed resource strategy from AD 1300 to 1450. The pre-Columbian impact on the landscape began to lessen around AD 1450, as shown by an increase in savanna trees and gallery forest. Although agriculture at the site began to decline prior to European arrival, the abandonment of raised fields was protracted, with evidence of sweet potato cultivation occurring as late as AD 1800.

Improved seasonal prediction of the hot summer of 2003 over Europe through better representation of uncertainty in the land surface

Methods to explicitly represent uncertainties in weather and climate models have been developed and refined over the past decade, and have reduced biases and improved forecast skill when implemented in the atmospheric component of models. These methods have not yet been applied to the land surface component of models. Since the land surface is strongly coupled to the atmospheric state at certain times and in certain places (such as the European summer of 2003), improvements in the representation of land surface uncertainty may potentially lead to improvements in atmospheric forecasts for such events. Here we analyse seasonal retrospective forecasts for 1981–2012 performed with the European Centre for Medium-Range Weather Forecasts’ (ECMWF) coupled ensemble forecast model. We consider two methods of incorporating uncertainty into the land surface model (H-TESSEL): stochastic perturbation of tendencies, and static perturbation of key soil parameters. We find that the perturbed parameter approach considerably improves the forecast of extreme air temperature for summer 2003, through better representation of negative soil moisture anomalies and upward sensible heat flux. Averaged across all the reforecasts the perturbed parameter experiment shows relatively little impact on the mean bias, suggesting perturbations of at least this magnitude can be applied to the land surface without any degradation of model climate. There is also little impact on skill averaged across all reforecasts and some evidence of overdispersion for soil moisture. The stochastic tendency experiments show a large overdispersion for the soil temperature fields, indicating that the perturbation here is too strong. There is also some indication that the forecast of the 2003 warm event is improved for the stochastic experiments, however the improvement is not as large as observed for the perturbed parameter experiment.

Heathland restoration on former agricultural land: effects of artificial acidification on the availability and uptake of toxic metal cations

Lowland heath is an internationally important habitat type that has greatly declined in abundance throughout Western Europe. In recent years this has led to a growing interest in the restoration of heathland on agricultural land. This generally requires the use of chemical treatments to return soil chemical conditions to those appropriate for the support of heathland ecosystems. However, the potential for negative impacts on the environment due to the potential of these treatments to increase the availability of trace metals via raised soil acidity requires investigation. A large-scale field study investigated the effect of two chemical treatments used in heathland restoration, elemental sulphur and ferrous sulphate, on soil acidity and whether it is possible to predict the effect of the treatments on availability of two potentially toxic cations (Al and Cd) in the soil along with their subsequent accumulation in the shoots of the grass Agrostis capillaris. Results showed that both treatments decreased soil pH, but that only elemental sulphur produced a pH similar to heathland soil. The availability of Al, measured by extraction with 1 M ammonium nitrate, could not be predicted by soil pH, depth in the soil and total Al concentration in the soil. By contrast, availability of Cd could be predicted from these three variables. Concentrations of both Al and Cd in the shoots of A. capillaris showed no significant relationship with the extractable concentration in the soil. Results are discussed in light of the possible environmental impacts of the chemical restoration techniques.

Fine-scale temporal characterization of trends in soil water dissolved organic carbon and potential drivers

Long-term monitoring of surface water quality has shown increasing concentrations of Dissolved Organic Carbon (DOC) across a large part of the Northern Hemisphere. Several drivers have been implicated including climate change, land management change, nitrogen and sulphur deposition and CO2 enrichment. Analysis of stream water data, supported by evidence from laboratory studies, indicates that an effect of declining sulphur deposition on catchment soil chemistry is likely to be the primary mechanism, but there are relatively few long term soil water chemistry records in the UK with which to investigate this, and other, hypotheses directly. In this paper, we assess temporal relationships between soil solution chemistry and parameters that have been argued to regulate DOC production and, using a unique set of co-located measurements of weather and bulk deposition and soil solution chemistry provided by the UK Environmental Change Network and the Intensive Forest Monitoring Level II Network . We used statistical non-linear trend analysis to investigate these relationships at 5 forested and 4 non-forested sites from 1993 to 2011. Most trends in soil solution DOC concentration were found to be non-linear. Significant increases in DOC occurred mostly prior to 2005. The magnitude and sign of the trends was associated qualitatively with changes in acid deposition, the presence/absence of a forest canopy, soil depth and soil properties. The strongest increases in DOC were seen in acidic forest soils and were most clearly linked to declining anthropogenic acid deposition, while DOC trends at some sites with westerly locations appeared to have been influenced by shorter-term hydrological variation. The results indicate that widespread DOC increases in surface waters observed elsewhere, are most likely dominated by enhanced mobilization of DOC in surficial organic horizons, rather than changes in the soil water chemistry of deeper horizons. While trends in DOC concentrations in surface horizons have flattened out in recent years, further increases may be expected as soil chemistry continues to adjust to declining inputs of acidity.

Complexation of copper by sewage sludge-derived dissolved organic matter: Effects on soil sorption behaviour and plant uptake

The complexation of Cu by sewage sludge-derived dissolved organic matter (SSDOM) is a process by which the environmental significance of the element may become enhanced due to reduced soil sorption and, hence, increased mobility. The work described in this paper used an ion selective electrode procedure to show that SSDOM complexation of Cu was greatest at intermediate pH values because competition between hydrogen ions and Cu for SSDOM binding sites, and between hydroxyl ions and SSDOM as Cu ligands, was lowest at such values. Batch sorption experiments further showed that the process of Cu complexation by SSDOM provided an explanation for enhanced desorption of Cu from the solid phase of a contaminated, organic matter-rich, clay loam soil, and reduced adsorption of Cu onto the solid phase of a sandy loam soil. Complexation of Cu by SSDOM did not affect uptake of Cu by spring barley plants, when compared to free ionic Cu, in a sand-culture pot experiment. However, it did appear to lead to greater biomass yields of the plant; perhaps indicating that the Cu-SSDOM complex had a lower toxicity towards the plant than the free Cu ion.

Preparation of consistent soil data sets for modelling purposes: Secondary SOTER data for four case study areas

The common GIS-based approach to regional analyses of soil organic carbon (SOC) stocks and changes is to define geographic layers for which unique sets of driving variables are derived, which include land use, climate, and soils. These GIS layers, with their associated attribute data, can then be fed into a range of empirical and dynamic models. Common methodologies for collating and formatting regional data sets on land use, climate, and soils were adopted for the project Assessment of Soil Organic Carbon Stocks and Changes at National Scale (GEFSOC). This permitted the development of a uniform protocol for handling the various input for the dynamic GEFSOC Modelling System. Consistent soil data sets for Amazon-Brazil, the Indo-Gangetic Plains (IGP) of India, Jordan and Kenya, the case study areas considered in the GEFSOC project, were prepared using methodologies developed for the World Soils and Terrain Database (SOTER). The approach involved three main stages: (1) compiling new soil geographic and attribute data in SOTER format; (2) using expert estimates and common sense to fill selected gaps in the measured or primary data; (3) using a scheme of taxonomy-based pedotransfer rules and expert-rules to derive soil parameter estimates for similar soil units with missing soil analytical data. The most appropriate approach varied from country to country, depending largely on the overall accessibility and quality of the primary soil data available in the case study areas. The secondary SOTER data sets discussed here are appropriate for a wide range of environmental applications at national scale. These include agro-ecological zoning, land evaluation, modelling of soil C stocks and changes, and studies of soil vulnerability to pollution. Estimates of national-scale stocks of SOC, calculated using SOTER methods, are presented as a first example of database application. Independent estimates of SOC stocks are needed to evaluate the outcome of the GEFSOC Modelling System for current conditions of land use and climate. (C) 2007 Elsevier B.V. All rights reserved.

Remediation of soils contaminated with petrol and diesel with lime

Lime treatment of hydrocarbon-contaminated soils offers the potential to stabilize and solidify these materials, with a consequent reduction in the risks associated with the leachate emanating from them. This can aid the disposal of contaminated soils or enable their on-site treatment. In this study, the addition of hydrated lime and quicklime significantly reduced the leaching of total petroleum hydrocarbons (TPH) from soils polluted with a 50:50 petrol/diesel mixture. Treatment with quicklime was slightly more effective, but hydrated lime may be better in the field because of its ease of handling. It is proposed that this occurs as a consequence of pozzolanic reactions retaining the hydrocarbons within the soil matrix. There was some evidence that this may be a temporary effect, as leaching increased between seven and 21 days after treatment, but the TPH concentrations in the leachate of treated soils were still one order of magnitude below those of the control soil, offering significant protection to groundwater. The reduction in leaching following treatment was observed in both aliphatic and aromatic fractions, but the latter were more affected because of their higher solubilty. The results are discussed in the context of risk assessment, and recommendations for future research are made.

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.

The Phosphorus Indicators Tool: a simple model of diffuse P loss from agricultural land to water

The Phosphorus Indicators Tool provides a catchment-scale estimation of diffuse phosphorus (P) loss from agricultural land to surface waters using the most appropriate indicators of P loss. The Tool provides a framework that may be applied across the UK to estimate P loss, which is sensitive not only to land use and management but also to environmental factors such as climate, soil type and topography. The model complexity incorporated in the P Indicators Tool has been adapted to the level of detail in the available data and the need to reflect the impact of changes in agriculture. Currently, the Tool runs on an annual timestep and at a 1 km(2) grid scale. We demonstrate that the P Indicators Tool works in principle and that its modular structure provides a means of accounting for P loss from one layer to the next, and ultimately to receiving waters. Trial runs of the Tool suggest that modelled P delivery to water approximates measured water quality records. The transparency of the structure of the P Indicators Tool means that identification of poorly performing coefficients is possible, and further refinements of the Tool can be made to ensure it is better calibrated and subsequently validated against empirical data, as it becomes available.

Fine sediment delivery and transfer in lowland catchments: modelling suspended sediment concentrations in response to hydrological forcing

Fine sediment delivery to and storage in stream channel reaches can disrupt aquatic habitats, impact river hydromorphology, and transfer adsorbed nutrients and pollutants from catchment slopes to the fluvial system. This paper presents a modelling toot for simulating the time-dependent response of the fine sediment system in catchments, using an integrated approach that incorporates both land phase and in-stream processes of sediment generation, storage and transfer. The performance of the model is demonstrated by applying it to simulate in-stream suspended sediment concentrations in two lowland catchments in southern England, the Enborne and the Lambourn, which exhibit contrasting hydrological and sediment responses due to differences in substrate permeability. The sediment model performs well in the Enborne catchment, where direct runoff events are frequent and peak suspended sediment concentrations can exceed 600 mg l(-1). The general trends in the in-stream concentrations in the Lambourn catchment are also reproduced by the model, although the observed concentrations are low (rarely exceeding 50 mg l(-1)) and the background variability in the concentrations is not fully characterized by the model. Direct runoff events are rare in this highly permeable catchment, resulting in a weak coupling between the sediment delivery system and the catchment hydrology. The generic performance of the model is also assessed using a generalized sensitivity analysis based on the parameter bounds identified in the catchment applications. Results indicate that the hydrological parameters contributing to the sediment response include those controlling (1) the partitioning of runoff between surface and soil zone flows and (2) the fractional loss of direct runoff volume prior to channel delivery. The principal sediment processes controlling model behaviour in the simulations are the transport capacity of direct runoff and the in-stream generation, storage and release of the fine sediment fraction. The in-stream processes appear to be important in maintaining the suspended sediment concentrations during low flows in the River Enborne and throughout much of the year in the River Lambourn. Copyright (c) 2007 John Wiley & Sons, Ltd.

Evaluation of two soil carbon models using two Kenyan long term experimental datasets

RothC and Century are two of the most widely used soil organic matter (SOM) models. However there are few examples of specific parameterisation of these models for environmental conditions in East Africa. The aim of this study was therefore, to evaluate the ability of RothC and the Century to estimate changes in soil organic carbon (SOC) resulting from varying land use/management practices for the climate and soil conditions found in Kenya. The study used climate, soils and crop data from a long term experiment (1976-2001) carried out at The Kabete site at The Kenya National Agricultural Research Laboratories (NARL, located in a semi-humid region) and data from a 13 year experiment carried out in Machang'a (Embu District, located in a semi-arid region). The NARL experiment included various fertiliser (0, 60 and 120 kg of N and P2O5 ha(-1)), farmyard manure (FYM - 5 and 10 t ha(-1)) and plant residue treatments, in a variety of combinations. The Machang'a experiment involved a fertiliser (51 kg N ha(-1)) and a FYM (0, 5 and 10 t ha(-1)) treatment with both monocropping and intercropping. At Kabete both models showed a fair to good fit to measured data, although Century simulations for treatments with high levels of FYM were better than those without. At the Machang'a site with monocrops, both models showed a fair to good fit to measured data for all treatments. However, the fit of both models (especially RothC) to measured data for intercropping treatments at Machang'a was much poorer. Further model development for intercrop systems is recommended. Both models can be useful tools in soil C Predictions, provided time series of measured soil C and crop production data are available for validating model performance against local or regional agricultural crops. (C) 2007 Elsevier B.V. All rights reserved.

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.