Phosphorus sorption and availability in soils amended with animal manures and sewage sludge

Soils that receive large applications of animal wastes and sewage sludge are vulnerable to releasing environmentally significant concentrations of dissolved P available to subsurface flow owing to the gradual saturation of the soil's P sorption capacity. This study evaluated P sorption (calculated from Langmuir isotherms) and availability of P (as CaCl2-P and resin P) in soils incubated for 20 d with poultry litter, poultry manure, cattle slurry, municipal sewage sludge, or KH2PO4, added on a P-equivalent basis (100 mg P kg(-1)). All the P sources had a marked negative effect on P sorption and a positive effect on P availability in all soils. In the cattle slurry- and KH2PO4- treated soils, the decreases in P sorption maximum (19-66%) and binding energy (25-89%) were consistently larger than the corresponding decreases (7-41% and 11-30%) in poultry litter-, poultry manure-, and sewage sludge-treated soils. The effects of cattle slurry and KH2PO4 on P availability were, in most cases, larger than those of the other P sources. In the poultry litter, poultry manure, and sewage sludge treatments, the increase in soil solution P was inversely related (R-2 = 0.75) to the input of Ca from these relatively high Ca (13.5-42 g kg(-1)) sources. Correlation analyses implied that the magnitude of the changes in P sorption and availability was not related to the water-extractable P content of the P sources. Future research on the sustainable application of organic wastes to agricultural soils needs to consider the non-P- as well as P-containing components of the waste.

Contribution of N Species and P Fractions to Stream Water Quality in Agricultural Catchments

The contribution from agricultural catchments to stream nitrogen and phosphorus concentrations was assessed by evaluation of the chemical composition of these nutrients in agricultural runoff for both surface and subsurface flow pathways. A range of land uses (grazed and ungrazed grassland, cereals, roots) in intensive agricultural systems was studied at scales from hillslope plots (0.5m2) to large catchment (>300km2). By fractionating the total nutrient load it was possible to establish that most of the phosphorus was transported in the unreactive (particulate and organic) fraction via surface runoff. This was true regardless of the scale of measurement. The form of the nitrogen load varied with land use and grazing intensity. High loads of dissolved inorganic nitrogen (with >90% transported as NH4-N) were recorded in surface runoff from heavily grazed land. In subsurface flow from small (2km2) subcatchments and in larger (>300 km2) catchments, organic nitrogen was found to be an important secondary constituent of the total nitrogen load, comprising 40% of the total annual load.

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.

Influence of UV radiation on chlorophyll, and antioxidant enzymes of wetland plants in different types of constructed wetland

A surface- and vertical subsurface-flow-constructed wetland were designed to study the response of chlorophyll and antioxidant enzymes to elevated UV radiation in three types of wetland plants (Canna indica, Phragmites austrail, and Typha augustifolia). Results showed that (1) chlorophyll content of C. indica, P. austrail, and T. augustifolia in the constructed wetland was significantly lower where UV radiation was increased by 10 and 20 % above ambient solar level than in treatment with ambient solar UV radiation (p < 0.05). (2) The malondialdehyde (MDA) content, guaiacol peroxidase (POD), superoxide dismutase (SOD), and catalase (CAT) activities of wetland plants increased with elevated UV radiation intensity. (3) The increased rate of MDA, SOD, POD, and CAT activities of C. indica, P. australis, and T. angustifolia by elevated UV radiation of 10 % was higher in vertical subsurface-flow-constructed wetland than in surface-flow-constructed wetland. The sensitivity of MDA, SOD, POD, and CAT activities of C. indica, P. austrail, and T. augustifolia to the elevated UV radiation was lower in surface-flow-constructed wetland than in the vertical subsurface-flow-constructed wetland, which was related to a reduction in UV radiation intensity through the dissolved organic carbon and suspended matter in the water. C. indica had the highest SOD and POD activities, which implied it is more sensitive to enhanced UV radiation. Therefore, different wetland plants had different antioxidant enzymes by elevated UV radiation, which were more sensitive in vertical subsurface-flow-constructed wetland than in surface-flow-constructed wetland.

Multi‐method global sensitivity analysis (MMGSA) for modelling floodplain hydrological processes

When studying hydrological processes with a numerical model, global sensitivity analysis (GSA) is essential if one is to understand the impact of model parameters and model formulation on results. However, different definitions of sensitivity can lead to a difference in the ranking of importance of the different model factors. Here we combine a fuzzy performance function with different methods of calculating global sensitivity to perform a multi-method global sensitivity analysis (MMGSA). We use an application of a finite element subsurface flow model (ESTEL-2D) on a flood inundation event on a floodplain of the River Severn to illustrate this new methodology. We demonstrate the utility of the method for model understanding and show how the prediction of state variables, such as Darcian velocity vectors, can be affected by such a MMGSA. This paper is a first attempt to use GSA with a numerically intensive hydrological model.

Multi-method global sensitivity analysis (MMGSA) for modelling floodplain hydrological processes

When studying hydrological processes with a numerical model, global sensitivity analysis (GSA) is essential if one is to understand the impact of model parameters and model formulation on results. However, different definitions of sensitivity can lead to a difference in the ranking of importance of the different model factors. Here we combine a fuzzy performance function with different methods of calculating global sensitivity to perform a multi-method global sensitivity analysis (MMGSA). We use an application of a finite element subsurface flow model (ESTEL-2D) on a flood inundation event on a floodplain of the River Severn to illustrate this new methodology. We demonstrate the utility of the method for model understanding and show how the prediction of state variables, such as Darcian velocity vectors, can be affected by such a MMGSA. This paper is a first attempt to use GSA with a numerically intensive hydrological model

Catchment-scale modelling of flow and nutrient transport in the Chalk unsaturated zone

The unsaturated zone exerts a major control on the delivery of nutrients to Chalk streams, yet flow and transport processes in this complex, dual-porosity medium have remained controversial. A major challenge arises in characterising these processes, both at the detailed mechanistic level and at an appropriate level for inclusion within catchment-scale models for nutrient management. The lowland catchment research (LOCAR) programme in the UK has provided a unique set of comprehensively instrumented groundwater-dominated catchments. Of these, the Pang and Lambourn, tributaries of the Thames near Reading, have been a particular focus for research into subsurface processes and surface water-groundwater interactions. Data from LOCAR and other sources, along with a new dual permeability numerical model of the Chalk, have been used to explore the relative roles of matrix and fracture flow within the unsaturated zone and resolve conflicting hypotheses of response. From the improved understanding gained through these explorations, a parsimonious conceptualisation of the general response of flow and transport within the Chalk unsaturated zone was formulated. This paper summarises the modelling and data findings of these explorations, and describes the integration of the new simplified unsaturated zone representation with a catchment-scale model of nutrients (INCA), resulting in a new model for catchment-scale flow and transport within Chalk systems: INCA-Chalk. This model is applied to the Lambourn, and results, including hindcast and forecast simulations, are presented. These clearly illustrate the decadal time-scales that need to be considered in the context of nutrient management and the EU Water Framework Directive. (C) 2007 Elsevier B.V. All rights reserved.

The influence of viscous and latent heating on crystal-rich magma flow in a conduit

The flow dynamics of crystal-rich high-viscosity magma is likely to be strongly influenced by viscous and latent heat release. Viscous heating is observed to play an important role in the dynamics of fluids with temperature-dependent viscosities. The growth of microlite crystals and the accompanying release of latent heat should play a similar role in raising fluid temperatures. Earlier models of viscous heating in magmas have shown the potential for unstable (thermal runaway) flow as described by a Gruntfest number, using an Arrhenius temperature dependence for the viscosity, but have not considered crystal growth or latent heating. We present a theoretical model for magma flow in an axisymmetric conduit and consider both heating effects using Finite Element Method techniques. We consider a constant mass flux in a 1-D infinitesimal conduit segment with isothermal and adiabatic boundary conditions and Newtonian and non-Newtonian magma flow properties. We find that the growth of crystals acts to stabilize the flow field and make the magma less likely to experience a thermal runaway. The additional heating influences crystal growth and can counteract supercooling from degassing-induced crystallization and drive the residual melt composition back towards the liquidus temperature. We illustrate the models with results generated using parameters appropriate for the andesite lava dome-forming eruption at Soufriere Hills Volcano, Montserrat. These results emphasize the radial variability of the magma. Both viscous and latent heating effects are shown to be capable of playing a significant role in the eruption dynamics of Soufriere Hills Volcano. Latent heating is a factor in the top two kilometres of the conduit and may be responsible for relatively short-term (days) transients. Viscous heating is less restricted spatially, but because thermal runaway requires periods of hundreds of days to be achieved, the process is likely to be interrupted. Our models show that thermal evolution of the conduit walls could lead to an increase in the effective diameter of flow and an increase in flux at constant magma pressure.

Impact of Tibetan Orography and Heating on the Summer Flow over Asia

The influence on the summer flow over Asia of both the orographic and thermal forcing of the Tibetan Plateau is investigated using a sequence of idealised experiments with a global primitive equation model. The zonally averaged flow is prescribed and both realistic and idealised orography and heating are used. There is some similarity between the responses to the two forcings when applied separately. The upper tropospheric Tibetan anticyclone is predominantly forced by the heating but also weakly by the orography. Below this, both forcings tend to give air descending in an equatorward anticyclonic circulation down the isentropes to the west and rising in a similar poleward circulation to the east. However the heating-only response has a strong ascending southwesterly flow that is guided around the south and south-east of the orography when it is included. On the northern side, the westerly flow over the orography gives ascent on the upslope and descent on the downslope. It is found that heating over the Plateau leads to a potential vorticity (PV) minimum and that if it is sufficiently strong the flow is unstable, producing a quasi-biweekly oscillation. During this oscillation the Tibetan anticyclone changes between a single centre over the southwestern side of the Plateau and a split/double structure with centres over China and the Middle East. These characteristics are similar to observed variability in the region. Associated with this quasi-biweekly oscillation are significant variations in the strength of the ascent over the Plateau and the Rossby wave pattern over the North Pacific. The origin of the variability is instability associated with the zonally extended potential vorticity PV minimum on a θ-surface, as proposed by Hsu and Plumb (2000). This minimum is due to the tendency to reduce the PV above the heating over the Plateau and to advection by the consequent anticyclone of high PV around from the east and low PV to the west. The deep convection to the south and southeast of the Plateau tends to suppress the quasi-biweekly oscillation because the low PV produced above it acts to reduce the meridional PV gradient reversal. The occurrence of the oscillation depends on the relative magnitude of the heating in the two regions.

Spatial variability of flow statistics within regular building arrays

Turbulence statistics obtained by direct numerical simulations are analysed to investigate spatial heterogeneity within regular arrays of building-like cubical obstacles. Two different array layouts are studied, staggered and square, both at a packing density of λp=0.25 . The flow statistics analysed are mean streamwise velocity ( u− ), shear stress ( u′w′−−−− ), turbulent kinetic energy (k) and dispersive stress fraction ( u˜w˜ ). The spatial flow patterns and spatial distribution of these statistics in the two arrays are found to be very different. Local regions of high spatial variability are identified. The overall spatial variances of the statistics are shown to be generally very significant in comparison with their spatial averages within the arrays. Above the arrays the spatial variances as well as dispersive stresses decay rapidly to zero. The heterogeneity is explored further by separately considering six different flow regimes identified within the arrays, described here as: channelling region, constricted region, intersection region, building wake region, canyon region and front-recirculation region. It is found that the flow in the first three regions is relatively homogeneous, but that spatial variances in the latter three regions are large, especially in the building wake and canyon regions. The implication is that, in general, the flow immediately behind (and, to a lesser extent, in front of) a building is much more heterogeneous than elsewhere, even in the relatively dense arrays considered here. Most of the dispersive stress is concentrated in these regions. Considering the experimental difficulties of obtaining enough point measurements to form a representative spatial average, the error incurred by degrading the sampling resolution is investigated. It is found that a good estimate for both area and line averages can be obtained using a relatively small number of strategically located sampling points.

Structure of turbulent flow over arrays of cubical roughness

The structure of turbulent flow over large roughness consisting of regular arrays of cubical obstacles is investigated numerically under constant pressure gradient conditions. Results are analysed in terms of first- and second-order statistics, by visualization of instantaneous flow fields and by conditional averaging. The accuracy of the simulations is established by detailed comparisons of first- and second-order statistics with wind-tunnel measurements. Coherent structures in the log region are investigated. Structure angles are computed from two-point correlations, and quadrant analysis is performed to determine the relative importance of Q2 and Q4 events (ejections and sweeps) as a function of height above the roughness. Flow visualization shows the existence of low-momentum regions (LMRs) as well as vortical structures throughout the log layer. Filtering techniques are used to reveal instantaneous examples of the association of the vortices with the LMRs, and linear stochastic estimation and conditional averaging are employed to deduce their statistical properties. The conditional averaging results reveal the presence of LMRs and regions of Q2 and Q4 events that appear to be associated with hairpin-like vortices, but a quantitative correspondence between the sizes of the vortices and those of the LMRs is difficult to establish; a simple estimate of the ratio of the vortex width to the LMR width gives a value that is several times larger than the corresponding ratio over smooth walls. The shape and inclination of the vortices and their spatial organization are compared to recent findings over smooth walls. Characteristic length scales are shown to scale linearly with height in the log region. Whilst there are striking qualitative similarities with smooth walls, there are also important differences in detail regarding: (i) structure angles and sizes and their dependence on distance from the rough surface; (ii) the flow structure close to the roughness; (iii) the roles of inflows into and outflows from cavities within the roughness; (iv) larger vortices on the rough wall compared to the smooth wall; (v) the effect of the different generation mechanism at the wall in setting the scales of structures.

The effect of flow field & geometry on the dynamic contact angle

A number of recent experiments suggest that, at a given wetting speed, the dynamic contact angle formed by an advancing liquid-gas interface with a solid substrate depends on the flow field and geometry near the moving contact line. In the present work, this effect is investigated in the framework of an earlier developed theory that was based on the fact that dynamic wetting is, by its very name, a process of formation of a new liquid-solid interface (newly “wetted” solid surface) and hence should be considered not as a singular problem but as a particular case from a general class of flows with forming or/and disappearing interfaces. The results demonstrate that, in the flow configuration of curtain coating, where a liquid sheet (“curtain”) impinges onto a moving solid substrate, the actual dynamic contact angle indeed depends not only on the wetting speed and material constants of the contacting media, as in the so-called slip models, but also on the inlet velocity of the curtain, its height, and the angle between the falling curtain and the solid surface. In other words, for the same wetting speed the dynamic contact angle can be varied by manipulating the flow field and geometry near the moving contact line. The obtained results have important experimental implications: given that the dynamic contact angle is determined by the values of the surface tensions at the contact line and hence depends on the distributions of the surface parameters along the interfaces, which can be influenced by the flow field, one can use the overall flow conditions and the contact angle as a macroscopic multiparametric signal-response pair that probes the dynamics of the liquid-solid interface. This approach would allow one to investigate experimentally such properties of the interface as, for example, its equation of state and the rheological properties involved in the interface’s response to an external torque, and would help to measure its parameters, such as the coefficient of sliding friction, the surface-tension relaxation time, and so on.