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.

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.

Scale-invariant moving finite elements for nonlinear partial differential equations in two dimensions

A scale-invariant moving finite element method is proposed for the adaptive solution of nonlinear partial differential equations. The mesh movement is based on a finite element discretisation of a scale-invariant conservation principle incorporating a monitor function, while the time discretisation of the resulting system of ordinary differential equations is carried out using a scale-invariant time-stepping which yields uniform local accuracy in time. The accuracy and reliability of the algorithm are successfully tested against exact self-similar solutions where available, and otherwise against a state-of-the-art h-refinement scheme for solutions of a two-dimensional porous medium equation problem with a moving boundary. The monitor functions used are the dependent variable and a monitor related to the surface area of the solution manifold. (c) 2005 IMACS. Published by Elsevier B.V. All rights reserved.

Boundary layer structure and decoupling from synoptic scale flow during NAMBLEX

This paper presents an overview of the meteorology and planetary boundary layer structure observed during the NAMBLEX field campaign to aid interpretation of the chemical and aerosol measurements. The campaign has been separated into five periods corresponding to the prevailing synoptic condition. Comparisons between meteorological measurements ( UHF wind profiler, Doppler sodar, sonic aneometers mounted on a tower at varying heights and a standard anemometer) and the ECMWF analysis at 10 m and 1100 m identified days when the internal boundary layer was decoupled from the synoptic flow aloft. Generally the agreement was remarkably good apart from during period one and on a few days during period four when the diurnal swing in wind direction implies a sea/land breeze circulation near the surface. During these periods the origin of air sampled at Mace Head would not be accurately represented by back trajectories following the winds resolved in ECMWF analyses. The wind profiler observations give a detailed record of boundary layer structure including an indication of its depth, average wind speed and direction. Turbulence statistics have been used to assess the height to which the developing internal boundary layer, caused by the increased surface drag at the coast, reaches the sampling location under a wide range of marine conditions. Sampling conducted below 10 m will be impacted by emission sources at the shoreline in all wind directions and tidal conditions, whereas sampling above 15 m is unlikely to be affected in any of the wind directions and tidal heights sampled during the experiment.