Spatial analysis of the error in a model of soil nitrogen

Models of the dynamics of nitrogen in soil (soil-N) can be used to aid the fertilizer management of a crop. The predictions of soil-N models can be validated by comparison with observed data. Validation generally involves calculating non-spatial statistics of the observations and predictions, such as their means, their mean squared-difference, and their correlation. However, when the model predictions are spatially distributed across a landscape the model requires validation with spatial statistics. There are three reasons for this: (i) the model may be more or less successful at reproducing the variance of the observations at different spatial scales; (ii) the correlation of the predictions with the observations may be different at different spatial scales; (iii) the spatial pattern of model error may be informative. In this study we used a model, parameterized with spatially variable input information about the soil, to predict the mineral-N content of soil in an arable field, and compared the results with observed data. We validated the performance of the N model spatially with a linear mixed model of the observations and model predictions, estimated by residual maximum likelihood. This novel approach allowed us to describe the joint variation of the observations and predictions as: (i) independent random variation that occurred at a fine spatial scale; (ii) correlated random variation that occurred at a coarse spatial scale; (iii) systematic variation associated with a spatial trend. The linear mixed model revealed that, in general, the performance of the N model changed depending on the spatial scale of interest. At the scales associated with random variation, the N model underestimated the variance of the observations, and the predictions were correlated poorly with the observations. At the scale of the trend, the predictions and observations shared a common surface. The spatial pattern of the error of the N model suggested that the observations were affected by the local soil condition, but this was not accounted for by the N model. In summary, the N model would be well-suited to field-scale management of soil nitrogen, but suited poorly to management at finer spatial scales. This information was not apparent with a non-spatial validation. (c),2007 Elsevier B.V. All rights reserved.

A simple model for predicting soil temperature in snow-covered and seasonally frozen soil: model description and testing

Microbial processes in soil are moisture, nutrient and temperature dependent and, consequently, accurate calculation of soil temperature is important for modelling nitrogen processes. Microbial activity in soil occurs even at sub-zero temperatures so that, in northern latitudes, a method to calculate soil temperature under snow cover and in frozen soils is required. This paper describes a new and simple model to calculate daily values for soil temperature at various depths in both frozen and unfrozen soils. The model requires four parameters average soil thermal conductivity, specific beat capacity of soil, specific heat capacity due to freezing and thawing and an empirical snow parameter. Precipitation, air temperature and snow depth (measured or calculated) are needed as input variables. The proposed model was applied to five sites in different parts of Finland representing different climates and soil types. Observed soil temperatures at depths of 20 and 50 cm (September 1981-August 1990) were used for model calibration. The calibrated model was then tested using observed soil temperatures from September 1990 to August 2001. R-2-values of the calibration period varied between 0.87 and 0.96 at a depth of 20 cm and between 0.78 and 0.97 at 50 cm. R-2 -values of the testing period were between 0.87 and 0.94 at a depth of 20cm. and between 0.80 and 0.98 at 50cm. Thus, despite the simplifications made, the model was able to simulate soil temperature at these study sites. This simple model simulates soil temperature well in the uppermost soil layers where most of the nitrogen processes occur. The small number of parameters required means, that the model is suitable for addition to catchment scale models.

Quantifying uncertainty in critical loads: (A) literature review

Critical loads are the basis for policies controlling emissions of acidic substances in Europe. The implementation of these policies involves large expenditures, and it is reasonable for policymakers to ask what degree of certainty can be attached to the underlying critical load and exceedance estimates. This paper is a literature review of studies which attempt to estimate the uncertainty attached to critical loads. Critical load models and uncertainty analysis are briefly outlined. Most studies have used Monte Carlo analysis of some form to investigate the propagation of uncertainties in the definition of the input parameters through to uncertainties in critical loads. Though the input parameters are often poorly known, the critical load uncertainties are typically surprisingly small because of a "compensation of errors" mechanism. These results depend on the quality of the uncertainty estimates of the input parameters, and a "pedigree" classification for these is proposed. Sensitivity analysis shows that some input parameters are more important in influencing critical load uncertainty than others, but there have not been enough studies to form a general picture. Methods used for dealing with spatial variation are briefly discussed. Application of alternative models to the same site or modifications of existing models can lead to widely differing critical loads, indicating that research into the underlying science needs to continue.

Quantifying uncertainty in critical loads: (B) acidity mass balance critical loads on a sensitive site

This paper reports an uncertainty analysis of critical loads for acid deposition for a site in southern England, using the Steady State Mass Balance Model. The uncertainty bounds, distribution type and correlation structure for each of the 18 input parameters was considered explicitly, and overall uncertainty estimated by Monte Carlo methods. Estimates of deposition uncertainty were made from measured data and an atmospheric dispersion model, and hence the uncertainty in exceedance could also be calculated. The uncertainties of the calculated critical loads were generally much lower than those of the input parameters due to a "compensation of errors" mechanism - coefficients of variation ranged from 13% for CLmaxN to 37% for CL(A). With 1990 deposition, the probability that the critical load was exceeded was > 0.99; to reduce this probability to 0.50, a 63% reduction in deposition is required; to 0.05, an 82% reduction. With 1997 deposition, which was lower than that in 1990, exceedance probabilities declined and uncertainties in exceedance narrowed as deposition uncertainty had less effect. The parameters contributing most to the uncertainty in critical loads were weathering rates, base cation uptake rates, and choice of critical chemical value, indicating possible research priorities. However, the different critical load parameters were to some extent sensitive to different input parameters. The application of such probabilistic results to environmental regulation is discussed.

Field trial using bone meal amendments to remediate mine waste derived soil contaminated with zinc, lead and cadmium

Bone meal amendments are being considered as a remediation method for metal-contaminated wastes. In various forms (biogenic, geogenic or synthetic), apatite, the principal mineral constituent of bone, has shown promise as an amendment to remediate metal-contaminated soils via the formation of insoluble phosphates of Pb and possibly other metals. The efficacy of commercially available bovine bone meal in this role was investigated in a field trial at Nenthead, Cumbria with a mine waste derived soil contaminated with Zn, Pb and Cd. Two 5 m(2) plots were set up: the first as a control and the second, a treatment plot where the soil was thoroughly mixed with bone meal to a depth of 50 cm at a soil to amendment ratio of 25:1 by weight. An array of soil solution samplers (Rhizon SMS (TM)) were installed in both plots and the soil pore water was collected and analysed for Ca, Cd, Zn and Pb regularly over a period of 2 a. Concurrently with the field trial, a laboratory trial with 800 mm high and 100 mm wide leaching Columns Was conducted using identical samplers and with soil from the held site. A substantial release of Zn, Pb, Cd and Ca was observed associated with the bone meal treatment. This release was transient in the case of the leaching columns, and showed seasonal variation in the case of the field trial. It is proposed that this effect resulted from metal complexation with organic acids released during breakdown of the bone meal organic fraction and was facilitated by the relatively high soil pH of 7.6-8.0. Even after this transient release effect had subsided or when incinerated bone meal was substituted in order to eliminate the organic fraction, no detectable decrease in dissolved metals was observed and no P was detected in solution, in contrast with an earlier small column laboratory study. It is concluded that due to the relative insolubility of apatite at above-neutral pH, the rate of supply of phosphate to soil solution was insufficient to result in significant precipitation of metal phosphates and that this may limit the effectiveness of the method to more acidic soils. (c) 2008 Elsevier Ltd. All rights reserved.

Influences of space, soil, nematodes and plants on microbial community composition of chalk grassland soils

Microbial communities respond to a variety of environmental factors related to resources (e.g. plant and soil organic matter), habitat (e.g. soil characteristics) and predation (e.g. nematodes, protozoa and viruses). However, the relative contribution of these factors on microbial community composition is poorly understood. Here, we sampled soils from 30 chalk grassland fields located in three different chalk hill ridges of Southern England, using a spatially explicit sampling scheme. We assessed microbial communities via phospholipid fatty acid (PLFA) analyses and PCR-denaturing gradient gel electrophoresis (DGGE) and measured soil characteristics, as well as nematode and plant community composition. The relative influences of space, soil, vegetation and nematodes on soil microorganisms were contrasted using variation partitioning and path analysis. Results indicate that soil characteristics and plant community composition, representing habitat and resources, shape soil microbial community composition, whereas the influence of nematodes, a potential predation factor, appears to be relatively small. Spatial variation in microbial community structure was detected at broad (between fields) and fine (within fields) scales, suggesting that microbial communities exhibit biogeographic patterns at different scales. Although our analysis included several relevant explanatory data sets, a large part of the variation in microbial communities remained unexplained (up to 92% in some analyses). However, in several analyses, significant parts of the variation in microbial community structure could be explained. The results of this study contribute to our understanding of the relative importance of different environmental and spatial factors in driving the composition of soil-borne microbial communities.

Soil moisture mediates association between the winter North Atlantic Oscillation and summer growth in the Park Grass Experiment

Several aspects of terrestrial ecosystems are known to be associated with the North Atlantic Oscillation (NAO) through effects of the NAO on winter climate, but recently the winter NAO has also been shown to be correlated with the following summer climate, including drought. Since drought is a major factor determining grassland primary productivity, the hypothesis was tested that the winter NAO is associated with summer herbage growth through soil moisture availability, using data from the Park Grass Experiment at Rothamsted, UK between 1960 and 1999. The herbage growth rate, mean daily rainfall, mean daily potential evapotranspiration (PE) and the mean and maximum potential soil moisture deficit (PSMD) were calculated between the two annual cuts in early summer and autumn for the unlimed, unfertilized plots. Mean and maximum PSMD were more highly correlated than rainfall or PE with herbage growth rate. Regression analysis showed that the natural logarithm of the herbage growth rate approximately halved for a 250 mm increase in maximum PSMD over the range 50-485 mm. The maximum PSMD was moderately correlated with the preceding winter NAO, with a positive winter NAO index associated with greater maximum PSMD. A positive winter NAO index was also associated with low herbage growth rate, accounting for 22% of the interannual variation in the growth rate. It was concluded that the association between the winter NAO and summer herbage growth rate is mediated by the PSMD in summer.

Simulating current global river runoff with a global hydrological model: model revisions, validation, and sensitivity analysis

Global hydrological models (GHMs) model the land surface hydrologic dynamics of continental-scale river basins. Here we describe one such GHM, the Macro-scale - Probability-Distributed Moisture model.09 (Mac-PDM.09). The model has undergone a number of revisions since it was last applied in the hydrological literature. This paper serves to provide a detailed description of the latest version of the model. The main revisions include the following: (1) the ability for the model to be run for n repetitions, which provides more robust estimates of extreme hydrological behaviour, (2) the ability of the model to use a gridded field of coefficient of variation (CV) of daily rainfall for the stochastic disaggregation of monthly precipitation to daily precipitation, and (3) the model can now be forced with daily input climate data as well as monthly input climate data. We demonstrate the effects that each of these three revisions has on simulated runoff relative to before the revisions were applied. Importantly, we show that when Mac-PDM.09 is forced with monthly input data, it results in a negative runoff bias relative to when daily forcings are applied, for regions of the globe where the day-to-day variability in relative humidity is high. The runoff bias can be up to - 80% for a small selection of catchments but the absolute magnitude of the bias may be small. As such, we recommend future applications of Mac-PDM.09 that use monthly climate forcings acknowledge the bias as a limitation of the model. The performance of Mac-PDM.09 is evaluated by validating simulated runoff against observed runoff for 50 catchments. We also present a sensitivity analysis that demonstrates that simulated runoff is considerably more sensitive to method of PE calculation than to perturbations in soil moisture and field capacity parameters.

Climate vs. soil factors in local adaptation of two common plant species

Evolutionary theory suggests that divergent natural selection in heterogeneous environments can result in locally adapted plant genotypes. To understand local adaptation it is important to study the ecological factors responsible for divergent selection. At a continental scale, variation in climate can be important while at a local scale soil properties could also play a role. We designed an experiment aimed to disentangle the role of climate and ( abiotic and biotic) soil properties in local adaptation of two common plant species. A grass (Holcus lanatus) and a legume ( Lotus corniculatus), as well as their local soils, were reciprocally transplanted between three sites across an Atlantic-Continental gradient in Europe and grown in common gardens in either their home soil or foreign soils. Growth and reproductive traits were measured over two growing seasons. In both species, we found significant environmental and genetic effects on most of the growth and reproductive traits and a significant interaction between the two environmental effects of soil and climate. The grass species showed significant home site advantage in most of the fitness components, which indicated adaptation to climate. We found no indication that the grass was adapted to local soil conditions. The legume showed a significant home soil advantage for number of fruits only and thus a weak indication of adaptation to soil and no adaptation to climate. Our results show that the importance of climate and soil factors as drivers of local adaptation is species-dependent. This could be related to differences in interactions between plant species and soil biota.

Long-term response of the nematode community to elevated atmospheric CO2 in a temperate dry grassland soil

Long-term effects of the elevated atmospheric CO2 on biosphere have been in focus of research since the last few decades. In this experiment undisturbed soil monoliths of loess grassland were exposed to an elevated CO2 environment (two-times the ambient CO2 level) for a period of six years with the aid of the open top chamber method. Control without a chamber and CO2 elevation was applied as well. Elevated CO2 level had very little impact oil soil food web. It did not influence either root and microbial biomass or microbial and nematode community structure. The only significant response was that density of the bacterial feeder genus Heterocephalobus increased in the chamber with elevated CO2 concentration. Application of the open top chambers initiated more changes on nematodes than the elevated CO2 level. Open top chamber (OTC) method decreased nematode density (total and plant feeder as well) to less than half of the original level. Negative effect was found on the genus level in the case of fungal feeder Aphelenchoides, plant feeder Helicotylenchus and Paratylenchus. It is very likely that the significantly lower belowground root biomass and partly its decreased quality reflected by the increased C/N ratio are the main responsible factors for the lower density of the plant feeder nematodes in the plots of chambers. According to diversity profiles, MI and MI(2-15) parameters, nematode communities in the open top chambers (both on ambient and elevated CO2 level) seem to be more structured than those under normal circumstances six years after start of the experiment.

Inter-laboratory variation of in vitro cumulative gas production profiles of feeds using manual and automated methods

A study was conducted to estimate variation among laboratories and between manual and automated techniques of measuring pressure on the resulting gas production profiles (GPP). Eight feeds (molassed sugarbeet feed, grass silage, maize silage, soyabean hulls, maize gluten feed, whole crop wheat silage, wheat, glucose) were milled to pass a I mm screen and sent to three laboratories (ADAS Nutritional Sciences Research Unit, UK; Institute of Grassland and Environmental Research (IGER), UK; Wageningen University, The Netherlands). Each laboratory measured GPP over 144 h using standardised procedures with manual pressure transducers (MPT) and automated pressure systems (APS). The APS at ADAS used a pressure transducer and bottles in a shaking water bath, while the APS at Wageningen and IGER used a pressure sensor and bottles held in a stationary rack. Apparent dry matter degradability (ADDM) was estimated at the end of the incubation. GPP were fitted to a modified Michaelis-Menten model assuming a single phase of gas production, and GPP were described in terms of the asymptotic volume of gas produced (A), the time to half A (B), the time of maximum gas production rate (t(RM) (gas)) and maximum gas production rate (R-M (gas)). There were effects (P<0.001) of substrate on all parameters. However, MPT produced more (P<0.001) gas, but with longer (P<0.001) B and t(RM gas) (P<0.05) and lower (P<0.001) R-M gas compared to APS. There was no difference between apparatus in ADDM estimates. Interactions occurred between substrate and apparatus, substrate and laboratory, and laboratory and apparatus. However, when mean values for MPT were regressed from the individual laboratories, relationships were good (i.e., adjusted R-2 = 0.827 or higher). Good relationships were also observed with APS, although they were weaker than for MPT (i.e., adjusted R-2 = 0.723 or higher). The relationships between mean MPT and mean APS data were also good (i.e., adjusted R 2 = 0. 844 or higher). Data suggest that, although laboratory and method of measuring pressure are sources of variation in GPP estimation, it should be possible using appropriate mathematical models to standardise data among laboratories so that data from one laboratory could be extrapolated to others. This would allow development of a database of GPP data from many diverse feeds. (c) 2005 Published by Elsevier B.V.

Influence of management type on Diptera communities of coniferous plantations and deciduous woodlands

The contribution of four types of secondary woodlands to Scottish invertebrate biodiversity was investigated for coniferous plantation forestry, riparian ash-alder woodlands, early successional deciduous woodlands and climax deciduous woodlands. Considerable variation in the type and intensity of management within these four woodland types existed. Adult Diptera from 21 families, representing diverse trophic and ecological guilds, were sampled from 31 woodlands in the Aberdeenshire region of northeast Scotland, between June and August 2001. Environmental differences between woodlands were recorded at each site using environmental parameters such as pH and organic matter content, vegetation characteristics, including percentage canopy cover and dominant field layer plant species. Multivariate ordination techniques detected significant responses in the Dipteran communities to soil type, organic matter content, soil pH, field layer plant species richness, dominant field layer plant species and percentage cover of Pteridium aquilinum. Responses in terms of Dipteran abundance, species richness, diversity and evenness were observed to soil type and dominant species of the field layer vegetation. The role of woodland type and management in diversifying Diptera communities is discussed with a view to maintain and possibly enhance Dipteran and other invertebrate communities in Scottish secondary woodlands. (C) 2002 Elsevier Science B.V. All rights reserved.

Population performance of collembolans feeding on soil fungi from different ecological niches

The potential reproductive value of arbuscular mycorrhizal fungi (Gloinus intraradices and Glomus invermaium), root pathogenic fungi (Rhizoctonia solani and Fusarium culmorum) and saprotrophic fungi (Penicillium hordei and Trichoderma harzianum) were examined for the collembolans Folsomia candida Willem and Folsomia fimetaria L. Dried baker's yeast (Saccharomyces cerevisiae) was used as a reference standard food in laboratory cultures. Collembolan performance was determined as final size, fecundity and population growth rate after when fed the fungal food sources for 31 days. The mycorrhizal fungi gave the least growth and fecundity compared with the other fungi, but G. intraradices gave good fecundity for F. candida. In terms of growth, Baker's yeast was a high-quality food for both adults and juveniles of both species, but it was a poorer food in terms of fecundity of F. candida. Preference of the fungi in all possible pairwise combinations showed that although F. fimetaria did not perform well on Glomus spp. and F. candida did not grow well on Glomus spp. their preference for these fungi did not reflect this. The highest fecundity was seen with the root pathogen F. culmorum. Different quality indicators such as the C:N ratio of the fungal food sources as well as other biological parameters are discussed in relation to their reproductive value and Collembola preferential feeding. (c) 2007 Elsevier Ltd. All rights reserved.

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

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

The influence of organic acids in relation to acid deposition in controlling the acidity of soil and stream waters on a seasonal basis

Much uncertainty still exists regarding the relative importance of organic acids in relation to acid deposition in controlling the acidity of soil and surface waters. This paper contributes to this debate by presenting analysis of seasonal variations in atmospheric deposition, soil solution and stream water chemistry for two UK headwater catchments with contrasting soils. Acid neutralising capacity (ANC), dissolved organic carbon (DOC) concentrations and the Na:Cl ratio of soil and stream waters displayed strong seasonal patterns with little seasonal variation observed in soil water pH. These patterns, plus the strong relationships between ANC, Cl and DOC, suggest that cation exchange and seasonal changes in the production of DOC and seasalt deposition are driving a shift in the proportion of acidity attributable to strong acid anions, from atmospheric deposition, during winter to predominantly organic acids in summer.