Modelling the impact of land use change on water quality in agricultural catchments

Export coefficient modelling was used to model the impact of agriculture on nitrogen and phosphorus loading on the surface waters of two contrasting agricultural catchments. The model was originally developed for the Windrush catchment where the highly reactive Jurassic limestone aquifer underlying the catchment is well connected to the surface drainage network, allowing the system to be modelled using uniform export coefficients for each nutrient source in the catchment, regardless of proximity to the surface drainage network. In the Slapton catchment, the hydrological path-ways are dominated by surface and lateral shallow subsurface flow, requiring modification of the export coefficient model to incorporate a distance-decay component in the export coefficients. The modified model was calibrated against observed total nitrogen and total phosphorus loads delivered to Slapton Ley from inflowing streams in its catchment. Sensitivity analysis was conducted to isolate the key controls on nutrient export in the modified model. The model was validated against long-term records of water quality, and was found to be accurate in its predictions and sensitive to both temporal and spatial changes in agricultural practice in the catchment. The model was then used to forecast the potential reduction in nutrient loading on Slapton Ley associated with a range of catchment management strategies. The best practicable environmental option (BPEO) was found to be spatial redistribution of high nutrient export risk sources to areas of the catchment with the greatest intrinsic nutrient retention capacity.

Land use for England and Wales: evaluation of management options to support 'good ecological status' in surface freshwaters

This paper analyses historic records of agricultural land use and management for England and Wales from 1931 and 1991 and uses export coefficient modelling to hindcast the impact of these practices on the rates of diffuse nitrogen (N) and phosphorus (P) export to water bodies for each of the major geo-climatic regions of England and Wales. Key trends indicate the importance of animal agriculture as a contributor to the total diffuse agricultural nutrient loading on waters, and the need to bring these sources under control if conditions suitable for sustaining 'Good Ecological Status' under the Water Framework Directive are to be generated. The analysis highlights the importance of measuring changes in nutrient loading in relation to the catchment-specific baseline state for different water bodies. The approach is also used to forecast the likely impact of broad regional scale scenarios on nutrient export to waters and highlights the need to take sensitive land out of production, introduce ceilings on fertilizer use and stocking densities, and controls on agricultural practice in higher risk areas where intensive agriculture is combined with a low intrinsic nutrient retention capacity, although the uncertainties associated with the modelling applied at this scale should be taken into account in the interpretation of model output. The paper advocates the need for a two-tiered approach to nutrient management, combining broad regional policies with targeted management in high risk areas at the catchment and farm scale.

Developing a reliable strategy to infer the effective soil hydraulic properties from field evaporation experiments for agro-hydrological models

The Richards equation has been widely used for simulating soil water movement. However, the take-up of agro-hydrological models using the basic theory of soil water flow for optimizing irrigation, fertilizer and pesticide practices is still low. This is partly due to the difficulties in obtaining accurate values for soil hydraulic properties at a field scale. Here, we use an inverse technique to deduce the effective soil hydraulic properties, based on measuring the changes in the distribution of soil water with depth in a fallow field over a long period, subject to natural rainfall and evaporation using a robust micro Genetic Algorithm. A new optimized function was constructed from the soil water contents at different depths, and the soil water at field capacity. The deduced soil water retention curve was approximately parallel but higher than that derived from published pedo-tranfer functions for a given soil pressure head. The water contents calculated from the deduced soil hydraulic properties were in good agreement with the measured values. The reliability of the deduced soil hydraulic properties was tested in reproducing data measured from an independent experiment on the same soil cropped with leek. The calculation of root water uptake took account for both soil water potential and root density distribution. Results show that the predictions of soil water contents at various depths agree fairly well with the measurements, indicating that the inverse analysis is an effective and reliable approach to estimate soil hydraulic properties, and thus permits the simulation of soil water dynamics in both cropped and fallow soils in the field accurately. (C) 2009 Elsevier B.V. All rights reserved.

Development of the Surface Urban Energy and Water balance Scheme (SUEWS) for cold climate cities

The Surface Urban Energy and Water Balance Scheme (SUEWS) is developed to include snow. The processes addressed include accumulation of snow on the different urban surface types: snow albedo and density aging, snow melting and re-freezing of meltwater. Individual model parameters are assessed and independently evaluated using long-term observations in the two cold climate cities of Helsinki and Montreal. Eddy covariance sensible and latent heat fluxes and snow depth observations are available for two sites in Montreal and one in Helsinki. Surface runoff from two catchments (24 and 45 ha) in Helsinki and snow properties (albedo and density) from two sites in Montreal are also analysed. As multiple observation sites with different land-cover characteristics are available in both cities, model development is conducted independent of evaluation. The developed model simulates snowmelt related runoff well (within 19% and 3% for the two catchments in Helsinki when there is snow on the ground), with the springtime peak estimated correctly. However, the observed runoff peaks tend to be smoother than the simulated ones, likely due to the water holding capacity of the catchments and the missing time lag between the catchment and the observation point in the model. For all three sites the model simulates the timing of the snow accumulation and melt events well, but underestimates the total snow depth by 18–20% in Helsinki and 29–33% in Montreal. The model is able to reproduce the diurnal pattern of net radiation and turbulent fluxes of sensible and latent heat during cold snow, melting snow and snow-free periods. The largest model uncertainties are related to the timing of the melting period and the parameterization of the snowmelt. The results show that the enhanced model can simulate correctly the exchange of energy and water in cold climate cities at sites with varying surface cover.

Modeling plant transpiration under limited soil water: comparison of different plant and soil hydraulic parameterizations and preliminary implications for their use in land surface models

Accurate estimates of how soil water stress affects plant transpiration are crucial for reliable land surface model (LSM) predictions. Current LSMs generally use a water stress factor, β, dependent on soil moisture content, θ, that ranges linearly between β = 1 for unstressed vegetation and β = 0 when wilting point is reached. This paper explores the feasibility of replacing the current approach with equations that use soil water potential as their independent variable, or with a set of equations that involve hydraulic and chemical signaling, thereby ensuring feedbacks between the entire soil–root–xylem–leaf system. A comparison with the original linear θ-based water stress parameterization, and with its improved curvi-linear version, was conducted. Assessment of model suitability was focused on their ability to simulate the correct (as derived from experimental data) curve shape of relative transpiration versus fraction of transpirable soil water. We used model sensitivity analyses under progressive soil drying conditions, employing two commonly used approaches to calculate water retention and hydraulic conductivity curves. Furthermore, for each of these hydraulic parameterizations we used two different parameter sets, for 3 soil texture types; a total of 12 soil hydraulic permutations. Results showed that the resulting transpiration reduction functions (TRFs) varied considerably among the models. The fact that soil hydraulic conductivity played a major role in the model that involved hydraulic and chemical signaling led to unrealistic values of β, and hence TRF, for many soil hydraulic parameter sets. However, this model is much better equipped to simulate the behavior of different plant species. Based on these findings, we only recommend implementation of this approach into LSMs if great care with choice of soil hydraulic parameters is taken

Sensitivity of a coupled climate model to canopy interception capacity

The canopy interception capacity is a small but key part of the surface hydrology, which affects the amount of water intercepted by vegetation and therefore the partitioning of evaporation and transpiration. However, little research with climate models has been done to understand the effects of a range of possible canopy interception capacity parameter values. This is in part due to the assumption that it does not significantly affect climate. Near global evapotranspiration products now make evaluation of canopy interception capacity parameterisations possible. We use a range of canopy water interception capacity values from the literature to investigate the effect on climate within the climate model HadCM3. We find that the global mean temperature is affected by up to -0.64 K globally and -1.9 K regionally. These temperature impacts are predominantly due to changes in the evaporative fraction and top of atmosphere albedo. In the tropics, the variations in evapotranspiration affect precipitation, significantly enhancing rainfall. Comparing the model output to measurements, we find that the default canopy interception capacity parameterisation overestimates canopy interception loss (i.e. canopy evaporation) and underestimates transpiration. Overall, decreasing canopy interception capacity improves the evapotranspiration partitioning in HadCM3, though the measurement literature more strongly supports an increase. The high sensitivity of climate to the parameterisation of canopy interception capacity is partially due to the high number of light rain-days in the climate model that means that interception is overestimated. This work highlights the hitherto underestimated importance of canopy interception capacity in climate model hydroclimatology and the need to acknowledge the role of precipitation representation limitations in determining parameterisations.

Evaluation of water properties in HEA–HEMA hydrogels swollen in aqueous-PEG solutions using thermoanalytical techniques

Hydrogels are polymeric materials used in many pharmaceutical and biomedical applications due to their ability to form 3D hydrophilic polymeric networks, which can absorb large amounts of water. In the present work, polyethylene glycols (PEG) were introduced into the hydrogel liquid phase in order to improve the mechanical properties of hydrogels composed of 2-hydroxyethylacrylate and 2-hydroxyethylmethacrylate (HEA–HEMA) synthesized with different co-monomer compositions and equilibrated in water or in 20 % water–PEG 400 and 600 solutions. The thermoanalytical techniques [differential scanning calorimetry (DSC) and thermogravimetry (TG)] were used to evaluate the amount and properties of free and bound water in HEA–HEMA hydrogels. The internal structure and the mechanical properties of hydrogels were studied using scanning electron microscopy and friability assay. TG “loss-on-drying” experiments were applied to study the water-retention properties of hydrogels, whereas the combination of TG and DSC allowed estimating the total amount of freezable and non-freezing water in hydrogels. The results show that the addition of viscous co-solvent (PEG) to the liquid medium results in significant improvement of the mechanical properties of HEA–HEMA hydrogels and also slightly retards the water loss from the hydrogels. A redistribution of free and bound water in the hydrogels equilibrated in mixed solutions containing 20 vol% of PEGs takes place.

Leaf trapping and retention of particles by holm oak and other common tree species in Mediterranean urban environments

Holm oak (Quercus ilex), a widespread urban street tree in the Mediterranean region, is widely used as biomonitor of persistent atmospheric pollutants, especially particulate-bound metals. By using lab- and field-based experimental approaches, we compared the leaf-level capacity for particles’ capture and retention between Q. ilex and other common Mediterranean urban trees: Quercus cerris, Platanus × hispanica, Tilia cordata and Olea europaea. All applied methods were effective in quantifying particulate capture and retention, although not univocal in ranking species performances. Distinctive morphological features of leaves led to differences in species’ ability to trap and retain particles of different size classes and to accumulate metals after exposure to traffic in an urban street. Overall, P. × hispanica and T. cordata showed the largest capture potential per unit leaf area for most model particles (Na+ and powder particles), and street-level Cu and Pb, while Q. ilex acted intermediately. After wash-off experiments, P. × hispanica leaves had the greatest retention capacity among the tested species and O. europaea the lowest. We concluded that the Platanus planting could be considered in Mediterranean urban environments due to its efficiency in accumulating and retaining airborne particulates; however, with atmospheric pollution being typically higher in winter, the evergreen Q. ilex represents a better year-round choice to mitigate the impact of airborne particulate pollutants.

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

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

The impact of fallowing and green manuring on soil conditions and the growth of sugarcane

There are currently concerns within some sugar industries that long-term monoculture has led to soil degradation and consequent yield decline. An investigation was conducted in Swaziland to assess the effects of fallowing and green manuring practices, over a seven-month period, on sugarcane yields and the physical properties of a poorly draining clay soil. In the subsequent first sugarcane crop after planting, yields were improved from 129 t ha(-1) under continuous sugarcane to 141-144 t ha(-1) after fallowing and green manuring, but there were no significant responses in the first and second ratoon crops. Also, in the first crop after planting, root length index increased from 3.5 km m(-2) under continuous sugarcane to 5.2-6.8 km m(-2) after fallowing, and improved rooting was still evident in the first ratoon crop where there had been soil drying during the fallow period. Soil bulk density, total porosity and water-holding capacity were not affected by the fallowing practices. However, air-filled porosity increased from 11% under continuous sugarcane to 16% after fallowing, and steady state ponded infiltration rates were increased from 0.61 mm h(-1) to 1.34 mm h(-1), but these improvements were no longer evident after a year back under sugarcane. Levels of soil organic matter were reduced in all cases, probably as a result of the tillage operations involved. In the plant crop, root length was well correlated with air-filled porosity, indicating the importance of improving belowground air supply for crop production on poorly draining clay soils.

Macrophyte community structure and species occurrence in relation to environmental determinants in the ephemeral aquatic habitats of Gavdos, Greece

The aims of this study were to explore the environmental factors that determine the distribution of plant communities in temporary rock pools and provide a quantitative analysis of vegetation-environment relationships for five study sites on the island of Gavdos, southwest of Crete, Greece. Data from 99 rock pools were collected and analysed using Two-Way Indicator Species Analysis (TWINSPAN), Detrended Correspondence Analysis (DCA) and Canonical Correspondence Analysis (CCA) to identify the principal communities and environmental gradients that are linked to community distribution. A total of 46 species belonging to 21 families were recorded within the study area. The dominant families were Labiatae, Gramineae and Compositae while therophytes and chamaephytes were the most frequent life forms. The samples were classified into six community types using TWINSPAN, which were also corroborated by CCA analysis. The principal gradients for vegetation distribution, identified by CCA, were associated with water storage and water retention ability, as expressed by pool perimeter and water depth. Generalised Additive Models (GAMs) were employed to identify responses of four dominant rock pool species to water depth. The resulting species response curves showed niche differentiation in the cases of Callitriche pulchra and Tillaea vaillantii and revealed competition between Zannichellia pedunculata and Chara vulgaris. The use of classification in combination with ordination techniques resulted in a good discrimination between plant communities. Generalised Additive Models are a powerful tool in investigating species response curves to environmental gradients. The methodology adopted can be employed for improving baseline information on plant community ecology and distribution in Mediterranean ephemeral pools.

Assessment of the quality of cottage cheese produced from standard and protein-fortified skim milk

The effects of milk protein fortification on the texture and microstructure of cottage cheese curd were evaluated. Protein powder (92.6% protein) was added to the skim milk at a level of 0.4% (w/w) to produce curds. Control curds with no protein powder addition were also produced. These curds were analysed for differences in yield, total solids, curd size, texture and structure. It was found that the addition of protein powder contributed to a significant yield increase, which can be attributed to increased water retention, with better curd size distribution. Control curds were firmer than the fortified curds and the structure showed less open-pore structure as revealed by electron microscopy. However, the addition of dressing masked the textural differences, and a sensory panel was unable to distinguish between cheeses produced from fortified milk and controls.

Production and dissolution rates of earthworm-secreted calcium carbonate

Earthworms secrete granules of calcium carbonate. These are potentially important in soil biogeochemical cycles and are routinely recorded in archaeological studies of Quaternary soils. Production rates of calcium carbonate granules by the earthworm Lumbricus terrestris L. were determined over 27 days in a range of soils with differing chemical properties (pH, organic matter content, water holding capacity, bulk composition, cation exchange capacity and exchangeable cations). Production rate varied between soils, lay in the range 0–0.043 mmolCaCO3 (0–4.3 mg) earthworm−1 d−1 with an average rate of 8 × 10−3 mmolCaCO3 (0.8 mg) earthworm−1 d−1 and was significantly correlated (r = 0.68, P ≤ 0.01) with soil pH. In a second experiment lasting 315 days earthworms repeatedly (over periods of 39–57 days) produced comparable masses of granules. Converting individual earthworm granule production rates into fluxes expressed on per hectare of land per year basis depends heavily on estimates of earthworm numbers. Using values of 10–20 L. terrestris m−2 suggests a rate of 18– 3139 molCaCO3 ha−1 yr−1. Data obtained from flow-through dissolution experiments suggest that at near neutral pH, granule geometric surface areanormalised dissolution rates are similar to those for other biogenic and inorganic calcium carbonate. Fits of the data to the dissolution relationship r = k(1 − ˝)n where r = dissolution rate, k = a rate constant, ˝ = relative saturation and n = the reaction order gave values of k = 1.72 × 10−10 mol cm−2 s−1 and n = 1.8 for the geometric surface area-normalised rates and k = 3.51 × 10−13 mol cm−2 s−1 and n = 1.8 for the BET surface area-normalised rates. In 196 day leaching column experiments trends in granule dissolution rate referenced to soil chemistry corresponded to predictions made by the SLIM model for dissolution of limestone in soil. If soil solution approaches saturation with respect to calcium carbonate granule dissolution will slow or even stop and granules be preserved indefinitely. Granules have the potential to be a small but significant component of the biogeochemical cycling of C and Ca in soil.

Effect of drainage and wildfire on peat hydrophysical properties

Consecutive multiple disturbances to northern peatlands can dramatically impact peat hydrophysical properties. We examine the impact of a double disturbance (drainage and wildfire) on the hydrophysical and moisture retention properties of peat, a key regulator of peatland ecohydrological resilience, and compare this with the impact to each individual disturbance (drainage and wildfire). The compound effect of drainage and wildfire resulted in a shift of the surface datum down the peat profile, revealing denser peat. Less-dense near-surface peats that regulate water-table position and near-surface moisture content, both favourable to Sphagnum recolonization, were lost. At a drained peatland that was then subject to wildfire, peat bulk density increased by 14.1%, von Post humification class increased by two categories and water retention increased by 15.6%, compared with an adjacent burned but undrained (single disturbance) portion of the fen. We discuss the key hydrophysical metrics of peatland vulnerability and outline how they are affected by the isolated impacts of drainage and wildfire, as well as their combined effects. We demonstrate that multiple peatland disturbances have likely led to an increase in hydrological limitations to Sphagnum recovery, which may impact peatland ecohydrological resilience.