Arable crop rotations under abiotic stress: a dynamic economic model

Break crops and multi-crop rotations are common in arable farm management, and the soil quality inherited from a previous crop is one of the parameters that determine the gross margin that is achieved with a given crop from a given parcel of land. In previous work we developed a dynamic economic model to calculate the potential yield and gross margin of a set of crops grown in a selection of typical rotation scenarios, and we reported use of the model to calculate coexistence costs for GM maize grown in a crop rotation. The model predicts economic effects of pest and weed pressures in monthly time steps. Validation of the model in respect of specific traits is proceeding as data from trials with novel crop varieties is published. Alongside this aspect of the validation process, we are able to incorporate data representing the economic impact of abiotic stresses on conventional crops, and then use the model to predict the cumulative gross margin achievable from a sequence of conventional crops grown at varying levels of abiotic stress. We report new progress with this aspect of model validation. In this paper, we report the further development of the model to take account of abiotic stress arising from drought, flood, heat or frost; such stresses being introduced in addition to variable pest and weed pressure. The main purpose is to assess the economic incentive for arable farmers to adopt novel crop varieties having multiple ‘stacked’ traits introduced by means of various biotechnological tools available to crop breeders.

Microbial and biochemical soil quality indicators and their potential for differentiating areas under contrasting agricultural management regimes

The aim of this study was to examine interrelationships between functional biochemical and microbial indicators of soil quality, and their suitability to differentiate areas under contrasting agricultural management regimes. The study included five 0.8 ha areas on a sandy-loam soil which had received contrasting fertility and cropping regimes over a 5 year period. These were organically managed vegetable, vegetable -cereal and arable rotations, an organically managed grass clover ley, and a conventional cereal rotation. The organic areas had been converted from conventional cereal production 5 years prior to the start of the study. All of the biochemical analyses, including light fraction organic matter (LFOM) C and N, labile organic N (LON), dissolved organic N and water-soluble carbohydrates showed significant differences between the areas, although the nature of the relationships between the areas varied between the different parameters, and were not related to differences in total soil organic matter content. The clearest differences were seen in LFOM C and N and LON, which were higher in the organic arable area relative to the other areas. In the case of the biological parameters, there were differences between the areas for biomass-N, ATP, chitin content, and the ratios of ATP: biomass and basal respiration: biomass. For these parameters, the precise relationships between the areas varied. However, relative to the conventionally managed area, areas under organic management generally had lower biomass-N and higher ATP contents. Arbuscular mycorrhizal fungus colonization potential was extremely low in the conventional area relative to the organic areas. Further, metabolic diversity and microbial community level physiological profiles, determined by analysis of microbial community metabolism using Biolog GN plates and the activities of eight key nutrient cycling enzymes, grouped the organic areas together, but separated them from the conventional area. We conclude that microbial parameters are more effective and consistent indicators of management induced changes to soil quality than biochemical parameters, and that a variety of biochemical and microbial analyses should be used when considering the impact of management on soil quality. (C) 2004 Elsevier Ltd. All rights reserved.

The accuracy of using the Ulysses result of the spatial invariance of the radial heliospheric field to compute the open solar flux

We survey observations of the radial magnetic field in the heliosphere as a function of position, sunspot number, and sunspot cycle phase. We show that most of the differences between pairs of simultaneous observations, normalized using the square of the heliocentric distance and averaged over solar rotations, are consistent with the kinematic "flux excess" effect whereby the radial component of the frozen-in heliospheric field is increased by longitudinal solar wind speed structure. In particular, the survey shows that, as expected, the flux excess effect at high latitudes is almost completely absent during sunspot minimum but is almost the same as within the streamer belt at sunspot maximum. We study the uncertainty inherent in the use of the Ulysses result that the radial field is independent of heliographic latitude in the computation of the total open solar flux: we show that after the kinematic correction for the excess flux effect has been made it causes errors that are smaller than 4.5%, with a most likely value of 2.5%. The importance of this result for understanding temporal evolution of the open solar flux is reviewed.

A comparison of energy use in conventional and integrated arable farming systems in the UK

The LINK Integrated Farming Systems (LINK-IFS) Project (1992-1997) was setup to compare conventional and integrated arable farming systems (IAFS), concentrating on practical feasibility and economic viability, but also taking into account the level of inputs used and environmental impact. As part of this, an examination into energy use within the two systems was also undertaken. This paper presents the results from that analysis. The data used is from the six sites within the LINK-IFS Project, spread through the arable production areas of England and from the one site in Scotland, covering the 5 years of the project. The comparison of the energy used is based on the equipment and inputs used to produce I kg of each crop within the conventional and integrated rotations, and thereby the overall energy used for each system. The results suggest that, in terms of total energy used, the integrated system appears to be the most efficient. However, in terms of energy efficiency, energy use per kilogram of output, the results are less conclusive. (C) 2003 Elsevier Science B.V. All rights reserved.

The role of weeds in supporting biological diversity within crop fields

Weeds are major constraints on crop production, yet as part of the primary producers within farming systems, they may be important components of the agroecosystem. Using published literature, the role of weeds in arable systems for other above-ground trophic levels are examined. In the UK, there is evidence that weed flora have changed over the past century, with some species declining in abundance, whereas others have increased. There is also some evidence for a decline in the size of arable weed seedbanks. Some of these changes reflect improved agricultural efficiency, changes to more winter-sown crops in arable rotations and the use of more broad-spectrum herbicide combinations. Interrogation of a database of records of phytophagous insects associated with plant species in the UK reveals that many arable weed species support a high diversity of insect species. Reductions in abundances of host plants may affect associated insects and other taxa. A number of insect groups and farmland birds have shown marked population declines over the past 30 years. Correlational studies indicate that many of these declines are associated with changes in agricultural practices. Certainly reductions in food availability in winter and for nestling birds in spring are implicated in the declines of several bird species, notably the grey partridge, Perdix perdix . Thus weeds have a role within agroecosystems in supporting biodiversity more generally. An understanding of weed competitivity and the importance of weeds for insects and birds may allow the identification of the most important weed species. This may form the first step in balancing the needs for weed control with the requirements for biodiversity and more sustainable production methods.

Evaluation of data for atmospheric models: master equation/RRKM calculations on the combination reaction, BrO+NO2 -> BrONO2, a conundrum

Experimental data for the title reaction were modeled using master equation (ME)/RRKM methods based on the Multiwell suite of programs. The starting point for the exercise was the empirical fitting provided by the NASA (Sander, S. P.; Finlayson-Pitts, B. J.; Friedl, R. R.; Golden, D. M.; Huie, R. E.; Kolb, C. E.; Kurylo, M. J.; Molina, M. J.; Moortgat, G. K.; Orkin, V. L.; Ravishankara, A. R. Chemical Kinetics and Photochemical Data for Use in Atmospheric Studies, Evaluation Number 15; Jet Propulsion Laboratory: Pasadena, California, 2006)(1) and IUPAC (Atkinson, R.; Baulch, D. L.; Cox, R. A.: R. F. Hampson, J.; Kerr, J. A.; Rossi, M. J.; Troe, J. J. Phys. Chem. Ref. Data. 2000, 29, 167) 2 data evaluation panels, which represents the data in the experimental pressure ranges rather well. Despite the availability of quite reliable parameters for these calculations (molecular vibrational frequencies (Parthiban, S.; Lee, T. J. J. Chem. Phys. 2000, 113, 145)3 and a. value (Orlando, J. J.; Tyndall, G. S. J. Phys. Chem. 1996, 100,. 19398)4 of the bond dissociation energy, D-298(BrO-NO2) = 118 kJ mol(-1), corresponding to Delta H-0(circle) = 114.3 kJ mol(-1) at 0 K) and the use of RRKM/ME methods, fitting calculations to the reported data or the empirical equations was anything but straightforward. Using these molecular parameters resulted in a discrepancy between the calculations and the database of rate constants of a factor of ca. 4 at, or close to, the low-pressure limit. Agreement between calculation and experiment could be achieved in two ways, either by increasing Delta H-0(circle) to an unrealistically high value (149.3 kJ mol(-1)) or by increasing , the average energy transferred in a downward collision, to an unusually large value (> 5000 cm(-1)). The discrepancy could also be reduced by making all overall rotations fully active. The system was relatively insensitive to changing the moments of inertia in the transition state to increase the centrifugal effect. The possibility of involvement of BrOONO was tested and cannot account for the difficulties of fitting the data.

Modelling the cabling of rope systems

This paper presents a theoretical model of the torsional characteristics of parallel multi-part rope systems. In such systems, the ropes may cable, or wrap around each other, depending on the combination of applied torque, rope tension, length and spacing between the rope parts. Cabling constitutes a failure that might be retrievable but as such can seriously affect the performance of the rope system. The torsional characteristics of the system are very different before and after cabling, and theoretical models are given for both situations. Laboratory tests were performed on both two and four rope systems, with measurements being made of torque at rotations from 0 to 360 deg. Tests were run with different rope spacings, tensions and lengths and the results compared with predictions from the theoretical model. The conclusion from the test results was that the theoretical model predicts both the pre- and post-cabling torsional behaviour with an acceptable level of accuracy.

Modelling the torsional interaction of wire and polyester fibre ropes used for offshore moorings

Oil rig mooring lines have traditionally consisted of chain and wire rope. As production has moved into deeper water it has proved advantageous to incorporate sections of fibre rope into the mooring lines. However, this has highlighted torsional interaction problems that can occur when ropes of different types are joined together. This paper describes a method by which the torsional properties of ropes can be modelled and can then be used to calculate the rotation and torque for two ropes connected in series. The method uses numerical representations of the torsional characteristics of both the ropes, and equates the torque generated in each rope under load to determine the rotation at the connection point. Data from rope torsional characterization tests have been analysed to derive constants used in the numerical model. Constants are presented for: a six-strand wire rope; a torque-balanced fibre rope; and a fibre rope that has been designed to be torque-matched to stranded wire rope. The calculation method has been verified by comparing predicted rotations with measured test values. Worked examples are given for a six-strand wire rope connected, firstly, to a torque-balanced fibre rope that offers little rotational restraint, and, secondly, to a fibre rope whose torsional properties are matched to that of the wire rope.

Marker placement to describe the wrist movements during activities of daily living in cyclical tasks

Objective. To describe the wrist kinematics during movement through free range of motion and activities of daily living using a cyclical task. Design. The wrist angles were initially calculated in a calibration trial and then in two selected activities of daily living (jar opening and carton pouring). Background. Existing studies which describe the wrist movement do not address the specific application of daily activities. Moreover, the data presented from subject to subject may differ simply because of the non-cyclical nature of the upper limbs movements. Methods. The coordinates of external markers attached to bone references on the forearm and dorsal side of the hand were obtained using an optical motion capture system. The wrist angles were derived from free motion trials and successively calculated in four healthy subjects for two specific cyclical daily activities (opening a jar and pouring from a carton). Results. The free motions trial highlighted the interaction between the wrist angles. Both the jar opening and the carton pouring activity showed a repetitive pattern for the three angles within the cycle length. In the jar-opening task, the standard deviation for the whole population was 10.8degrees for flexion-extension, 5.3degrees for radial-ulnar deviation and 10.4degrees for pronation-supination. In the carton-pouring task, the standard deviation for the whole population was 16.0degrees for flexion-extension, 3.4degrees for radial-ulnar deviation and 10.7degrees for pro nation-supination. Conclusion. Wrist kinematics in healthy subjects can be successfully described by the rotations about the axes of marker-defined coordinates systems during free range of motion and daily activities using cyclical tasks.

An analytical framework for appraising Integrated Arable Farming Systems

Integrated Arable Farming Systems are examined from the perspective of the farmer considering the use of such techniques, and data are presented which suggest that the uptake of the approach may expose the manager to a greater degree of risk. Observations are made about the possible uptake of such systems in the UK and the implications this may have for agricultural and environmental policy in general.

Three dimensional picture of the local structure of 1,4-Polybutadiene from a complete atomistic model and neutron scattering data

An efficient method of combining neutron diffraction data over an extended Q range with detailed atomistic models is presented. A quantitative and qualitative mapping of the organization of the chain conformation in both glass and liquid phase has been performed. The proposed structural refinement method is based on the exploitation of the intrachain features of the diffraction pattern by the use of internal coordinates for bond lengths, valence angles and torsion rotations. Models are built stochastically by assignment of these internal coordinates from probability distributions with limited variable parameters. Variation of these parameters is used in the construction of models that minimize the differences between the observed and calculated structure factors. A series of neutron scattering data of 1,4-polybutadiene at the region 20320 K is presented. Analysis of the experimental data yield bond lengths for C-C and C=C of 1.54 and 1.35 Å respectively. Valence angles of the backbone were found to be at 112 and 122.8 for the CCC and CC=C respectively. Three torsion angles corresponding to the double bond and the adjacent R and β bonds were found to occupy cis and trans, s(, trans and g( and trans states, respectively. We compare our results with theoretical predictions, computer simulations, RIS models, and previously reported experimental results.

Riluzole neuroprotection in a parkinson’s disease model involves suppression of reactive astrocytosis but not GLT-1 regulation

Background Riluzole is a neuroprotective drug used in the treatment of motor neurone disease. Recent evidence suggests that riluzole can up-regulate the expression and activity of the astrocyte glutamate transporter, GLT-1. Given that regulation of glutamate transport is predicted to be neuroprotective in Parkinson's disease, we tested the effect of riluzole in parkinsonian rats which had received a unilateral 6-hydroxydopamine injection into the median forebrain bundle. Results Rats were treated with intraperitoneal riluzole (4 mg/kg or 8 mg/kg), 1 hour before the lesion then once daily for seven days. Riluzole produced a modest but significant attenuation of dopamine neurone degeneration, assessed by suppression of amphetamine-induced rotations, preservation of tyrosine hydroxylase positive neuronal cell bodies in the substantia nigra pars compacta and attenuation of striatal tyrosine hydroxylase protein loss. Seven days after 6-hydroxydopamine lesion, reactive astrocytosis was observed in the striatum, as determined by increases in expression of glial fibrillary acidic protein, however the glutamate transporter, GLT-1, which is also expressed in astrocytes was not regulated by the lesion. Conclusions The results confirm that riluzole is a neuroprotective agent in a rodent model of parkinson’s disease. Riluzole administration did not regulate GLT-1 levels but significantly reduced GFAP levels, in the lesioned striatum. Riluzole suppression of reactive astrocytosis is an intriguing finding which might contribute to the neuroprotective effects of this drug.

Advances in sequential data assimilation and numerical weather forecasting: an Ensemble Transform Kalman-Bucy Filter, a study on clustering in deterministic Ensemble Square Root Filters, and a test of a new time stepping scheme in an atmospheric model

This dissertation deals with aspects of sequential data assimilation (in particular ensemble Kalman filtering) and numerical weather forecasting. In the first part, the recently formulated Ensemble Kalman-Bucy (EnKBF) filter is revisited. It is shown that the previously used numerical integration scheme fails when the magnitude of the background error covariance grows beyond that of the observational error covariance in the forecast window. Therefore, we present a suitable integration scheme that handles the stiffening of the differential equations involved and doesn’t represent further computational expense. Moreover, a transform-based alternative to the EnKBF is developed: under this scheme, the operations are performed in the ensemble space instead of in the state space. Advantages of this formulation are explained. For the first time, the EnKBF is implemented in an atmospheric model. The second part of this work deals with ensemble clustering, a phenomenon that arises when performing data assimilation using of deterministic ensemble square root filters in highly nonlinear forecast models. Namely, an M-member ensemble detaches into an outlier and a cluster of M-1 members. Previous works may suggest that this issue represents a failure of EnSRFs; this work dispels that notion. It is shown that ensemble clustering can be reverted also due to nonlinear processes, in particular the alternation between nonlinear expansion and compression of the ensemble for different regions of the attractor. Some EnSRFs that use random rotations have been developed to overcome this issue; these formulations are analyzed and their advantages and disadvantages with respect to common EnSRFs are discussed. The third and last part contains the implementation of the Robert-Asselin-Williams (RAW) filter in an atmospheric model. The RAW filter is an improvement to the widely popular Robert-Asselin filter that successfully suppresses spurious computational waves while avoiding any distortion in the mean value of the function. Using statistical significance tests both at the local and field level, it is shown that the climatology of the SPEEDY model is not modified by the changed time stepping scheme; hence, no retuning of the parameterizations is required. It is found the accuracy of the medium-term forecasts is increased by using the RAW filter.

Ensemble clustering in deterministic ensemble Kalman filters

Ensemble clustering (EC) can arise in data assimilation with ensemble square root filters (EnSRFs) using non-linear models: an M-member ensemble splits into a single outlier and a cluster of M−1 members. The stochastic Ensemble Kalman Filter does not present this problem. Modifications to the EnSRFs by a periodic resampling of the ensemble through random rotations have been proposed to address it. We introduce a metric to quantify the presence of EC and present evidence to dispel the notion that EC leads to filter failure. Starting from a univariate model, we show that EC is not a permanent but transient phenomenon; it occurs intermittently in non-linear models. We perform a series of data assimilation experiments using a standard EnSRF and a modified EnSRF by a resampling though random rotations. The modified EnSRF thus alleviates issues associated with EC at the cost of traceability of individual ensemble trajectories and cannot use some of algorithms that enhance performance of standard EnSRF. In the non-linear regimes of low-dimensional models, the analysis root mean square error of the standard EnSRF slowly grows with ensemble size if the size is larger than the dimension of the model state. However, we do not observe this problem in a more complex model that uses an ensemble size much smaller than the dimension of the model state, along with inflation and localisation. Overall, we find that transient EC does not handicap the performance of the standard EnSRF.