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.

Host plant recognition by the root feeding clover weevil, Sitona lepidus (Coleoptera : Curculionidae)

This study investigated the ability of neonatal larvae of the root-feeding weevil, Sitona lepidus Gyllenhal, to locate white clover Trifolium repens L. (Fabaceae) roots growing in soil and to distinguish them from the roots of other species of clover and a co-occurring grass species. Choice experiments used a combination of invasive techniques and the novel technique of high resolution X-ray microtomography to non-invasively track larval movement in the soil towards plant roots. Burrowing distances towards roots of different plant species were also examined. Newly hatched S. lepidus recognized T. repens roots and moved preferentially towards them when given a choice of roots of subterranean clover, Trifolium subterraneum L. (Fabaceae), strawberry clover Trifolium fragiferum L. (Fabaceae), or perennial ryegrass Lolium perenne L. (Poaceae). Larvae recognized T. repens roots, whether released in groups of five or singly, when released 25 mm (meso-scale recognition) or 60 mm (macro-scale recognition) away from plant roots. There was no statistically significant difference in movement rates of larvae.

Causes of the apparent scale independence of fractal indices associated with forest fragmentation in Bolivia

In a recent investigation, Landsat TM and ETM+ data were used to simulate different resolutions of remotely-sensed images (from 30 to 1100 m) and to analyze the effect of resolution on a range of landscape metrics associated with spatial patterns of forest fragmentation in Chapare, Bolivia since the mid-1980s. Whereas most metrics were found to be highly dependent on pixel size, several fractal metrics (DLFD, MPFD, and AWMPFD) were apparently independent of image resolution, in contradiction with a sizeable body of literature indicating that fractal dimensions of natural objects depend strongly on image characteristics. The present re-analysis of the Chapare images, using two alternative algorithms routinely used for the evaluation of fractal dimensions, shows that the values of the box-counting and information fractal dimensions are systematically larger, sometimes by as much as 85%, than the "fractal" indices DLFD, MPFD, and AWMFD for the same images. In addition, the geometrical fractal features of the forest and non-forest patches in the Chapare region strongly depend on the resolution of images used in the analysis. The largest dependency on resolution occurs for the box-counting fractal dimension in the case of the non-forest patches in 1993, where the difference between the 30 and I 100 m-resolution images corresponds to 24% of the full theoretical range (1.0 to 2.0) of the mass fractal dimension. The observation that the indices DLFD, MPFD, and AWMPFD, unlike the classical fractal dimensions, appear relatively unaffected by resolution in the case of the Chapare images seems due essentially to the fact that these indices are based on a heuristic, "non-geometric" approach to fractals. Because of their lack of a foundation in fractal geometry, nothing guarantees that these indices will be resolution-independent in general. (C) 2006 International Society for Photogrammetry and Remote Sensing, Inc. (ISPRS). Published by Elsevier B.V. All rights reserved.

Sensitivity analysis of a catchment-scale nitrogen model

There are now considerable expectations that semi-distributed models are useful tools for supporting catchment water quality management. However, insufficient attention has been given to evaluating the uncertainties inherent to this type of model, especially those associated with the spatial disaggregation of the catchment. The Integrated Nitrogen in Catchments model (INCA) is subjected to an extensive regionalised sensitivity analysis in application to the River Kennet, part of the groundwater-dominated upper Thames catchment, UK The main results are: (1) model output was generally insensitive to land-phase parameters, very sensitive to groundwater parameters, including initial conditions, and significantly sensitive to in-river parameters; (2) INCA was able to produce good fits simultaneously to the available flow, nitrate and ammonium in-river data sets; (3) representing parameters as heterogeneous over the catchment (206 calibrated parameters) rather than homogeneous (24 calibrated parameters) produced a significant improvement in fit to nitrate but no significant improvement to flow and caused a deterioration in ammonium performance; (4) the analysis indicated that calibrating the flow-related parameters first, then calibrating the remaining parameters (as opposed to calibrating all parameters together) was not a sensible strategy in this case; (5) even the parameters to which the model output was most sensitive suffered from high uncertainty due to spatial inconsistencies in the estimated optimum values, parameter equifinality and the sampling error associated with the calibration method; (6) soil and groundwater nutrient and flow data are needed to reduce. uncertainty in initial conditions, residence times and nitrogen transformation parameters, and long-term historic data are needed so that key responses to changes in land-use management can be assimilated. The results indicate the general, difficulty of reconciling the questions which catchment nutrient models are expected to answer with typically limited data sets and limited knowledge about suitable model structures. The results demonstrate the importance of analysing semi-distributed model uncertainties prior to model application, and illustrate the value and limitations of using Monte Carlo-based methods for doing so. (c) 2005 Elsevier B.V. All rights reserved.

Why Bt cotton pays for small-scale producers in South Africa

There is much debate about the potential benefits (and costs) of genetically modified (GM) crop technology for developing countries. Studies have been carried out in Argentina, China, Indonesia and most recently India1 to assess the impact of Bacillus thuringiensis (Bt) cotton on farmers in those regions.

Bt-cotton boosts the gross margin of small-scale cotton producers in South Africa

This paper explores some of the issues involved in the Genetic Modification (GM) debate by focusing on one crop that has been modified for pest resistance, cotton (Gossypium hirsutum), and commercially released to small-scale farmers in the Makhathini Flats, KwaZulu Natal, the Republic of South Africa. This was the first commercial release of a GM variety (Bt-cotton) in Sub-Saharan Africa, and thus provides valuable and timely insights into some of the potential advantages and disadvantages of the technology for small-scale farmers in Africa. Even though there are wider concerns regarding the vulnerability of small-scale farmers in the area, the survey results suggest that Bt-cotton generated higher yields and gross margins than non-Bt-cotton. In addition, Bt-cotton significantly reduced the use of pesticide with consequent potential benefits to human health and the environment.

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.

Over-parameterised, uncertain 'mathematical marionettes' - How can we best use catchment water quality models? An example of an 80-year catchment-scale nutrient balance

The Integrated Catchment Model of Nitrogen (INCA-N) was applied to the River Lambourn, a Chalk river-system in southern England. The model's abilities to simulate the long-term trend and seasonal patterns in observed stream water nitrate concentrations from 1920 to 2003 were tested. This is the first time a semi-distributed, daily time-step model has been applied to simulate such a long time period and then used to calculate detailed catchment nutrient budgets which span the conversion of pasture to arable during the late 1930s and 1940s. Thus, this work goes beyond source apportionment and looks to demonstrate how such simulations can be used to assess the state of the catchment and develop an understanding of system behaviour. The mass-balance results from 1921, 1922, 1991, 2001 and 2002 are presented and those for 1991 are compared to other modelled and literature values of loads associated with nitrogen soil processes and export. The variations highlighted the problem of comparing modelled fluxes with point measurements but proved useful for identifying the most poorly understood inputs and processes thereby providing an assessment of input data and model structural uncertainty. The modelled terrestrial and instream mass-balances also highlight the importance of the hydrological conditions in pollutant transport. Between 1922 and 2002, increased inputs of nitrogen from fertiliser, livestock and deposition have altered the nitrogen balance with a shift from possible reduction in soil fertility but little environmental impact in 1922, to a situation of nitrogen accumulation in the soil, groundwater and instream biota in 2002. In 1922 and 2002 it was estimated that approximately 2 and 18 kg N ha(-1) yr(-1) respectively were exported from the land to the stream. The utility of the approach and further considerations for the best use of models are discussed. (C) 2008 Elsevier B.V. All rights reserved.

A contour-advective semi-lagrangian numerical algorithm for simulating fine-scale conservative dynamical fields

This paper describes a novel numerical algorithm for simulating the evolution of fine-scale conservative fields in layer-wise two-dimensional flows, the most important examples of which are the earth's atmosphere and oceans. the algorithm combines two radically different algorithms, one Lagrangian and the other Eulerian, to achieve an unexpected gain in computational efficiency. The algorithm is demonstrated for multi-layer quasi-geostrophic flow, and results are presented for a simulation of a tilted stratospheric polar vortex and of nearly-inviscid quasi-geostrophic turbulence. the turbulence results contradict previous arguments and simulation results that have suggested an ultimate two-dimensional, vertically-coherent character of the flow. Ongoing extensions of the algorithm to the generally ageostrophic flows characteristic of planetary fluid dynamics are outlined.

Intake, growth and meat quality of steers given diets based on varying proportions of maize silage and grass silage

Simmental × Holstein-Friesian steers were offered four forage diets. These comprised grass silage (G); proportionately 0·67 grass silage, proportionately 0·33 maize silage (GGM); 0·33 grass silage, 0·67 maize silage ( MMG); maize silage ( M) from 424 (s.d. = 11·5) kg to slaughter at a minimum weight of 560 kg. Forages were mixed and offered ad libitum. Steers were offered 2 kg of a concentrate daily, the concentrate being formulated such that all steers had similar crude protein intakes across dietary treatments. A sample of steers was slaughtered at the beginning of the experimental period to allow the calculation of the rate of gain of the carcass and its components. Carcass dissection of a sample of steers allowed the development of a prediction equation of carcass composition based on thoracic limb dissection of all carcasses. Forage dry matter intake and live-weight gain increased linearly as maize silage replaced grass silage in the forage mixture, resulting in improvements in food conversion ratio (all P = 0·001). Killing-out proportion increased with maize silage inclusion ( P < 0·001) but fat and conformation scores did not differ significantly between diets. However, increasing maize inclusion in the diet resulted in a greater weight ( P = 0·05) and proportion ( P = 0·008) of fat in the carcass, and significant increases in internal fat deposition. The inclusion of maize led to a progressive increase in the daily gains of carcass ( P < 0·001), and significant increases in the daily gains of both fat ( P < 0·001) and lean tissue ( P < 0·001). Fat colour was more yellow in cattle given diets G and GGM than diets MMG and M ( P < 0·001) and colour intensity was lower on diet M than the other three diets ( P < 0·001). There were no significant differences in any aspects of eating quality between diets. Therefore, maize silage has the potential to reduce the time taken for finishing beef animals to achieve slaughter weight with no apparent detrimental effects on subsequent meat quality.

Digestion, rumen fermentation and circulating concentrations of insulin, growth hormone and IGF-1 in steers fed diets based on different proportions of maize silage and grass silage

Replacing grass silage with maize silage results in a fundamental change in the ratio of structural to non-structural carbohydrates with commensurate changes in rumen fermentation patterns and nutrient utilisation. This study investigated the effects of feeding four forage mixtures, namely grass silage (G); 67 g/100 g grass silage133 g/100 g maize silage (GGM); 67 g/100 g maize silage133/100 g grass silage (MMG); maize silage (M) to four ruminally and duodenally canulated Holstein Friesian steers. All diets were formulated to be isonitrogenous (22.4 g N/kg DM) using a concentrate mixture. Dietary dry matter (DM) and organic matter (OM) digestibility increased with ascending maize silage inclusion (P,0.1) whereas starch and neutral detergent fibre digestibility declined (P,0.05). Ratio of non-glucogenic to glucogenic precursors in the rumen fluid increased with maize silage inclusion (P,0.01) with a commensurate reduction in rumen pH (P,0.05). Mean circulating concentrations of insulin were greatest and similar in diets MMG and GGM, lower in diet M and lowest in diet G (P,0.01). There were no effects of diet on the mean circulating concentration of growth hormone (GH), or the frequency, amplitude and duration of GH pulses, or the mean circulating concentrations of IGF-1. Increasing levels of DM, OM and starch intakes with the substitution of grass silage with maize silage affected overall digestion, nutrient partitioning and subsequent circulating concentrations of insulin.

The spacing of orographic rainbands triggered by small-scale topography

A combination of idealized numerical simulations and analytical theory is used to investigate the spacing between convective orographic rainbands over the Coastal Range of western Oregon. The simulations, which are idealized from an observed banded precipitation event over the Coastal Range, indicate that the atmospheric response to conditionally unstable flow over the mountain ridge depends strongly on the subridge-scale topographic forcing on the windward side of the ridge. When this small-scale terrain contains only a single scale (l) of terrain variability, the band spacing is identical to l, but when a spectrum of terrain scales are simultaneously present, the band spacing ranges between 5 and 10 km, a value that is consistent with observations. Based on the simulations, an inviscid linear model is developed to provide a physical basis for understanding the scale selection of the rainbands. This analytical model, which captures the transition from lee waves upstream of the orographic cloud to moist convection within it, reveals that the spacing of orographic rainbands depends on both the projection of lee-wave energy onto the unstable cap cloud and the growth rate of unstable perturbations within the cloud. The linear model is used in tandem with numerical simulations to determine the sensitivity of the band spacing to a number of environmental and terrain-related parameters.

The triggering of orographic rainbands by small-scale topography

The triggering of convective orographic rainbands by small-scale topographic features is investigated through observations of a banded precipitation event over the Oregon Coastal Range and simulations using a cloud-resolving numerical model. A quasi-idealized simulation of the observed event reproduces the bands in the radar observations, indicating the model’s ability to capture the physics of the band-formation process. Additional idealized simulations reinforce that the bands are triggered by lee waves past small-scale topographic obstacles just upstream of the nominal leading edge of the orographic cloud. Whether a topographic obstacle in this region is able to trigger a strong rainband depends on the phase of its lee wave at cloud entry. Convective growth only occurs downstream of obstacles that give rise to lee-wave-induced displacements that create positive vertical velocity anomalies w_c and nearly zero buoyancy anomalies b_c as air parcels undergo saturation. This relationship is quantified through a simple analytic condition involving w_c, b_c, and the static stability N_m^2 of the cloud mass. Once convection is triggered, horizontal buoyancy gradients in the cross-flow direction generate circulations that align the bands parallel to the flow direction.

Hybridization between Brassica napus and B. rapa on a national scale in the United Kingdom

Measures blocking hybridization would prevent or reduce biotic or environmental change caused by gene flow from genetically modified (GM) crops to wild relatives. The efficacy of any such measure depends on hybrid numbers within the legislative region over the life-span of the GM cultivar. We present a national assessment of hybridization between rapeseed (Brassica napus) and B. rapa from a combination of sources, including population surveys, remote sensing, pollen dispersal profiles, herbarium data, local Floras, and other floristic databases. Across the United Kingdom, we estimate that 32,000 hybrids form annually in waterside B. rapa populations, whereas the less abundant weedy populations contain 17,000 hybrids. These findings set targets for strategies to eliminate hybridization and represent the first step toward quantitative risk assessment on a national scale.