An assessment of the fine sediment dynamics in an upland river system: INCA-Sed modifications and implications for fisheries

There is a need for better links between hydrology and ecology, specifically between landscapes and riverscapes to understand how processes and factors controlling the transport and storage of environmental pollution have affected or will affect the freshwater biota. Here we show how the INCA modelling framework, specifically INCA-Sed (the Integrated Catchments model for Sediments) can be used to link sediment delivery from the landscape to sediment changes in-stream. INCA-Sed is a dynamic, process-based, daily time step model. The first complete description of the equations used in the INCA-Sed software (version 1.9.11) is presented. This is followed by an application of INCA-Sed made to the River Lugg (1077 km2) in Wales. Excess suspended sediment can negatively affect salmonid health. The Lugg has a large and potentially threatened population of both Atlantic salmon (Salmo salar) and Brown Trout (Salmo trutta). With the exception of the extreme sediment transport processes, the model satisfactorily simulated both the hydrology and the sediment dynamics in the catchment. Model results indicate that diffuse soil loss is the most important sediment generation process in the catchment. In the River Lugg, the mean annual Guideline Standard for suspended sediment concentration, proposed by UKTAG, of 25 mg l− 1 is only slightly exceeded during the simulation period (1995–2000), indicating only minimal effect on the Atlantic salmon population. However, the daily time step simulation of INCA-Sed also allows the investigation of the critical spawning period. It shows that the sediment may have a significant negative effect on the fish population in years with high sediment runoff. It is proposed that the fine settled particles probably do not affect the salmonid egg incubation process, though suspended particles may damage the gills of fish and make the area unfavourable for spawning if the conditions do not improve.

Moving towards a more mechanistic approach in the determination of soil heat flux from remote measurements - I. A universal approach to calculate thermal inertia

In this paper we pledge that physically based equations should be combined with remote sensing techniques to enable a more theoretically rigorous estimation of area-average soil heat flux, G. A standard physical equation (i.e. the analytical or exact method) for the estimation of G, in combination with a simple, but theoretically derived, equation for soil thermal inertia (F), provides the basis for a more transparent and readily interpretable method for the estimation of G; without the requirement for in situ instrumentation. Moreover, such an approach ensures a more universally applicable method than those derived from purely empirical studies (employing vegetation indices and albedo, for example). Hence, a new equation for the estimation of Gamma(for homogeneous soils) is discussed in this paper which only requires knowledge of soil type, which is readily obtainable from extant soil databases and surveys, in combination with a coarse estimate of moisture status. This approach can be used to obtain area-averaged estimates of Gamma(and thus G, as explained in paper II) which is important for large-scale energy balance studies that employ aircraft or satellite data. Furthermore, this method also relaxes the instrumental demand for studies at the plot and field scale (no requirement for in situ soil temperature sensors, soil heat flux plates and/or thermal conductivity sensors). In addition, this equation can be incorporated in soil-vegetation-atmosphere-transfer models that use the force restore method to update surface temperatures (such as the well-known ISBA model), to replace the thermal inertia coefficient.

A universal agro-hydrological model for water and nitrogen cycles in the soil-crop system SMCR_N: critical update and further validation

Agro-hydrological models have widely been used for optimizing resources use and minimizing environmental consequences in agriculture. SMCRN is a recently developed sophisticated model which simulates crop response to nitrogen fertilizer for a wide range of crops, and the associated leaching of nitrate from arable soils. In this paper, we describe the improvements of this model by replacing the existing approximate hydrological cascade algorithm with a new simple and explicit algorithm for the basic soil water flow equation, which not only enhanced the model performance in hydrological simulation, but also was essential to extend the model application to the situations where the capillary flow is important. As a result, the updated SMCRN model could be used for more accurate study of water dynamics in the soil-crop system. The success of the model update was demonstrated by the simulated results that the updated model consistently out-performed the original model in drainage simulations and in predicting time course soil water content in different layers in the soil-wheat system. Tests of the updated SMCRN model against data from 4 field crop experiments showed that crop nitrogen offtakes and soil mineral nitrogen in the top 90 cm were in a good agreement with the measured values, indicating that the model could make more reliable predictions of nitrogen fate in the crop-soil system, and thus provides a useful platform to assess the impacts of nitrogen fertilizer on crop yield and nitrogen leaching from different production systems. (C) 2010 Elsevier B.V. All rights reserved.

Reduction of arsenic uptake by lettuce with ferrous sulfate applied to contaminated soil

Soil contamination by arsenic (As) presents a hazard in many countries and there is a need for techniques to minimize As uptake by plants. A proposed in situ remediation method was tested by growing lettuce (Lactuca sativa L. cv. Kermit) in a greenhouse pot experiment on soil that contained 577 mg As kg(-1), taken from a former As smelter site. All combinations of iron (Fe) oxides, at concentrations of 0.00, 0.22, 0.54, and 1.09% (w/w), and lime, at concentrations of 0.00, 0.27, 0.68, and 1.36% (w/w), were tested in a factorial design. To create the treatments, field-moist soil, commercial-grade FeSO4, and ground agricultural lime were mixed and stored for one week, allowing Fe oxides to precipitate. Iron oxides gave highly significant (P < 0.001) reductions in lettuce As concentrations, down to 11% of the lettuce As concentration for untreated soil. For the Fe oxides and lime treatment combinations where soil pH was maintained nearly constant, the lettuce As concentration declined in an exponential relationship with increasing FeSO4 application rate and lettuce yield was almost unchanged. Iron oxides applied at a concentration of 1.09% did not give significantly lower lettuce As concentrations than the 0.54% treatment. Simultaneous addition of lime with FeSO4 was essential. Ferrous sulfate with insufficient lime lowered soil pH and caused mobilization of Al, Ba, Co, Cr, Cu, Fe, K, Mg, Mn, Na, Ni, Pb, Sr, and Zn. At the highest Fe oxide to lime ratios, Mn toxicity caused severe yield loss.

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.

A geostatistical analysis of the spatial variation of soil mineral nitrogen and potentially available nitrogen within an arable field.

The technology for site-specific applications of nitrogen (N) fertilizer has exposed a gap in our knowledge about the spatial variation of soil mineral N, and that which will become available during the growing season within arable fields. Spring mineral N and potentially available N were measured in an arable field together with gravimetric water content, loss on ignition, crop yield, percentages of sand, silt, and clay, and elevation to describe their spatial variation geostatistically. The areas with a larger clay content had larger values of mineral N, potentially available N, loss on ignition and gravimetric water content, and the converse was true for the areas with more sandy soil. The results suggest that the spatial relations between mineral N and loss on ignition, gravimetric water content, soil texture, elevation and crop yield, and between potentially available N and loss on ignition and silt content could be used to indicate their spatial patterns. Variable-rate nitrogen fertilizer application would be feasible in this field because of the spatial structure and the magnitude of variation of mineral N and potentially available N.

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.

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.

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

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

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.

Estimating the spatial scale of herbicide and soil interactions by nested sampling, hierarchical analysis of variance and residual maximum likelihood

An unbalanced nested sampling design was used to investigate the spatial scale of soil and herbicide interactions at the field scale. A hierarchical analysis of variance based on residual maximum likelihood (REML) was used to analyse the data and provide a first estimate of the variogram. Soil samples were taken at 108 locations at a range of separating distances in a 9 ha field to explore small and medium scale spatial variation. Soil organic matter content, pH, particle size distribution, microbial biomass and the degradation and sorption of the herbicide, isoproturon, were determined for each soil sample. A large proportion of the spatial variation in isoproturon degradation and sorption occurred at sampling intervals less than 60 m, however, the sampling design did not resolve the variation present at scales greater than this. A sampling interval of 20-25 m should ensure that the main spatial structures are identified for isoproturon degradation rate and sorption without too great a loss of information in this field.

Adsorption of water vapor by bare soil in an olive grove in southern Spain

Data for water vapor adsorption and evaporation are presented for a bare soil (sandy loam, clay content 15%) in a southern Spanish olive grove. Water losses and gains were measured using eight high-precision minilysimeters, placed around an olive tree, which had been irrigated until the soil reached field capacity (similar to 0.22 m(3) m(-3)). They were subsequently left to dry for 10 days. A pair of lysimeters was situated at each of the main points of the compass (N, E, S, W), at a distance of 1 m (the inner set of lysimeters; ILS) and 2 m (the outer set of lysimeters; OLS), respectively, from the tree trunk. Distinct periods of moisture loss (evaporation) and moisture gain (vapor adsorption) could be distinguished for each day. Vapor adsorption often started just after noon and generally lasted until the (early) evening. Values of up to 0.7 mm of adsorbed water per day were measured. Adsorption was generally largest for the OLS (up to 100% more on a daily basis), and increased during the dry down. This was mainly the result of lower OLS surface soil moisture contents (period-average absolute difference similar to 0.005 m(3) m(-3)), as illustrated using various analyses employing a set of micrometeorological equations describing the exchange of water vapor between bare soil and the atmosphere. These analyses also showed that the amount of water vapor adsorbed by soils is very sensitive to changes in atmospheric forcing and surface variables. The use of empirical equations to estimate vapor adsorption is therefore not recommended.

The diurnal course of soil moisture as measured by various dielectric sensors: Effects of soil temperature and the implications for evaporation estimates

Soil moisture content, theta, of a bare and vegetated UK gravelly sandy loam soil (in situ and repacked in small lysimeters) was measured using various dielectric instruments (single-sensor ThetaProbes, multi-sensor Profile Probes, and Aquaflex Sensors), at depths ranging between 0.03 and I m, during the summers of 2001 (in situ soil) and 2002 (mini-lysimeters). Half-hourly values of evaporation, E, were calculated from diurnal changes in total soil profile water content, using the soil water balance equation. For the bare soil field, Profile Probes and ML2x ThetaProbes indicated a diurnal course of theta that did not concur with typical soil physical observations: surface layer soil moisture content increased from early morning until about midday, after which theta declined, generally until the early evening. The unexpected course of theta was positively correlated to soil temperature, T-s, also at deeper depths. Aquaflex and ML1 ThetaProbe (older models) outputs, however, reflected common observations: 0 increased slightly during the night (capillary rise) and decreased from the morning until late afternoon (as a result of evaporation). For the vegetated plot, the spurious diurnal theta fluctuations were less obvious, because canopy shading resulted in lower amplitudes of T-s. The unrealistic theta profiles measured for the bare and vegetated field sites caused diurnal estimates of E to attain downward daytime and upward night-time values. In the mini-lysimeters, at medium to high moisture contents, theta values measured by (ML2x) ThetaProbes followed a relatively realistic course, and predictions of E from diurnal changes in vertically integrated theta generally compared well with lysimeter estimates of E. However, time courses of theta and E became comparable to those observed for the field plots when the soil in the lysimeters reached relatively low values of theta. Attempts to correct measured theta for fluctuations in T, revealed that no generally applicable formula could be derived. (c) 2005 Elsevier B.V. All rights reserved.

Reduction of arsenic uptake by lettuce with ferrous sulfate applied to contaminated soil

Soil contamination by arsenic (As) presents a hazard in many countries and there is a need for techniques to minimize As uptake by plants. A proposed in situ remediation method was tested by growing lettuce (Lactuca sativa L. cv. Kermit) in a greenhouse pot experiment on soil that contained 577 mg As kg(-1), taken from a former As smelter site. All combinations of iron (Fe) oxides, at concentrations of 0.00, 0.22, 0.54, and 1.09% (w/w), and lime, at concentrations of 0.00, 0.27, 0.68, and 1.36% (w/w), were tested in a factorial design. To create the treatments, field-moist soil, commercial-grade FeSO4, and ground agricultural lime were mixed and stored for one week, allowing Fe oxides to precipitate. Iron oxides gave highly significant (P < 0.001) reductions in lettuce As concentrations, down to 11% of the lettuce As concentration for untreated soil. For the Fe oxides and lime treatment combinations where soil pH was maintained nearly constant, the lettuce As concentration declined in an exponential relationship with increasing FeSO4 application rate and lettuce yield was almost unchanged. Iron oxides applied at a concentration of 1.09% did not give significantly lower lettuce As concentrations than the 0.54% treatment. Simultaneous addition of lime with FeSO4 was essential. Ferrous sulfate with insufficient lime lowered soil pH and caused mobilization of Al, Ba, Co, Cr, Cu, Fe, K, Mg, Mn, Na, Ni, Pb, Sr, and Zn. At the highest Fe oxide to lime ratios, Mn toxicity caused severe yield loss.

Consumer welfare and the loss induced by withholding information: The case of BSE in Italy

The paper develops a measure of consumer welfare losses associated with withholding it formation about a possible link between BSE and vCJD. The Cost of Ignorance (COI) is measured by comparing the utility of the informed choice with the utility of the uninformed choice, under conditions of improved information. Unlike previous work that is largely based on a single equation demand model, the measure is obtained retrieving a cost,function from a dynamic Almost Ideal Demand System. The estimated perceived loss for Italian consumers due to delayed information ranges from 12 percent to 54 percent of total meat expenditure, depending on the month assumed to embody correct beliefs about the safety level of beef.