Abstracts of communications: Proceedings of the thirty-nineth meeting of the Agricultural Research Modellers' Group

This group, which is concerned with the applications of mathematics to agricultural science, was formed in 1970 and has since met at approximately yearly intervals in London for one-day meetings. The thirty-ninth meeting of the group, chaired by Professor N. Crout of the University of Nottingham, was held in the Kohn Centre at the Royal Society, 6 Carlton House Terrace, London on Friday, 30 March 2007 when the following papers were read.

Uncertainty in the estimation of potential evapotranspiration under climate change.

21st century climate change is projected to result in an intensification of the global hydrological cycle, but there is substantial uncertainty in how this will impact freshwater availability. A relatively overlooked aspect of this uncertainty pertains to how different methods of estimating potential evapotranspiration (PET) respond to changing climate. Here we investigate the global response of six different PET methods to a 2 °C rise in global mean temperature. All methods suggest an increase in PET associated with a warming climate. However, differences in PET climate change signal of over 100% are found between methods. Analysis of a precipitation/PET aridity index and regional water surplus indicates that for certain regions and GCMs, choice of PET method can actually determine the direction of projections of future water resources. As such, method dependence of the PET climate change signal is an important source of uncertainty in projections of future freshwater availability.

Enhanced phosphopeptide isolation by Fe(III)-IMAC using 1,1,1,3,3,3-hexafluoroisopropanol

IMAC can be used to selectively enrich phosphopeptides from complex peptide mixtures, but co-retention of acidic peptides together with the failure to retain some phosphopeptides restricts the general utility of the method. In this study Fe(III)-IMAC was qualitatively and quantitatively assessed using a panel of phosphopeptides, both synthetic and derived from proteolysis of known phosphoproteins, to identify the causes of success and failure in the application of this technique. Here we demonstrate that, as expected, peptides with a more acidic amino acid content are generally more efficiently purified and detected by MALDI-MS after Fe(III)-IMAC than those with a more basic content. Modulating the loading buffer used for Fe(III)-IMAC significantly affects phosphopeptide binding and suggests that conformational factors that lead to steric hindrance and reduced accessibility to the phosphate are important. The use of 1,1,1,3,3,3-hexafluoroisopropanol is shown here to significantly improve Fe(III)-IMAC enrichment and subsequent detection of phosphopeptides by MALDI-MS.

Enhanced phosphopeptide isolation by Fe(III)-IMAC using 1,1,1,3,3,3-hexafluoroisopropanol

IMAC can be used to selectively enrich phosphopeptides from complex peptide mixtures, but co-retention of acidic peptides together with the failure to retain some phosphopeptides restricts the general utility of the method. In this study Fe(III)-IMAC was qualitatively and quantitatively assessed using a panel of phosphopeptides, both synthetic and derived from proteolysis of known phosphoproteins, to identify the causes of success and failure in the application of this technique. Here we demonstrate that, as expected, peptides with a more acidic amino acid content are generally more efficiently purified and detected by MALDI-MS after Fe(III)-IMAC than those with a more basic content. Modulating the loading buffer used for Fe(III)-IMAC significantly affects phosphopeptide binding and suggests that conformational factors that lead to steric hindrance and reduced accessibility to the phosphate are important. The use of 1,1,1,3,3,3-hexa-fluoroisopropanol is shown here to significantly improve Fe(III)-IMAC enrichment and subsequent detection of phosphopeptides by MALDI-MS.

Oxidation of oleic acid at the air-water interface and its potential effects on cloud critical supersaturations

The oxidation of organic films on cloud condensation nuclei has the potential to affect climate and precipitation events. In this work we present a study of the oxidation of a monolayer of deuterated oleic acid (cis-9-octadecenoic acid) at the air-water interface by ozone to determine if oxidation removes the organic film or replaces it with a product film. A range of different aqueous sub-phases were studied. The surface excess of deuterated material was followed by neutron reflection whilst the surface pressure was followed using a Wilhelmy plate. The neutron reflection data reveal that approximately half the organic material remains at the air-water interface following the oxidation of oleic acid by ozone, thus cleavage of the double bond by ozone creates one surface active species and one species that partitions to the bulk (or gas) phase. The most probable products, produced with a yield of similar to(87 +/- 14)%, are nonanoic acid, which remains at the interface, and azelaic acid (nonanedioic acid), which dissolves into the bulk solution. We also report a surface bimolecular rate constant for the reaction between ozone and oleic acid of (7.3 +/- 0.9) x 10(-11) cm(2) molecule s(-1). The rate constant and product yield are not affected by the solution sub-phase. An uptake coefficient of ozone on the oleic acid monolayer of similar to 4 x 10(-6) is estimated from our results. A simple Kohler analysis demonstrates that the oxidation of oleic acid by ozone on an atmospheric aerosol will lower the critical supersaturation needed for cloud droplet formation. We calculate an atmospheric chemical lifetime of oleic acid of 1.3 hours, significantly longer than laboratory studies on pure oleic acid particles suggest, but more consistent with field studies reporting oleic acid present in aged atmospheric aerosol.

Oxidation of biogenic and water-soluble compounds in aqueous and organic aerosol droplets by ozone: a kinetic and product analysis approach using laser Raman tweezers

The results of an experimental study into the oxidative degradation of proxies for atmospheric aerosol are presented. We demonstrate that the laser Raman tweezers method can be used successfully to obtain uptake coeffcients for gaseous oxidants on individual aqueous and organic droplets, whilst the size and composition of the droplets is simultaneously followed. A laser tweezers system was used to trap individual droplets containing an unsaturated organic compound in either an aqueous or organic ( alkane) solvent. The droplet was exposed to gas- phase ozone and the reaction kinetics and products followed using Raman spectroscopy. The reactions of three different organic compounds with ozone were studied: fumarate anions, benzoate anions and alpha pinene. The fumarate and benzoate anions in aqueous solution were used to represent components of humic- like substances, HULIS; a alpha- pinene in an alkane solvent was studied as a proxy for biogenic aerosol. The kinetic analysis shows that for these systems the diffusive transport and mass accommodation of ozone is relatively fast, and that liquid- phase di. ffusion and reaction are the rate determining steps. Uptake coe. ffcients, g, were found to be ( 1.1 +/- 0.7) x 10(-5), ( 1.5 +/- 0.7) x 10 (-5) and ( 3.0 - 7.5) x 10 (-3) for the reactions of ozone with the fumarate, benzoate and a- pinene containing droplets, respectively. Liquid- phase bimolecular rate coe. cients for reactions of dissolved ozone molecules with fumarate, benzoate and a- pinene were also obtained: k(fumarate) = ( 2.7 +/- 2) x 10 (5), k(benzoate) = ( 3.5 +/- 3) x 10 (5) and k(alpha-pinene) = ( 1-3) x 10(7) dm(3) mol (-1) s (- 1). The droplet size was found to remain stable over the course of the oxidation process for the HULIS- proxies and for the oxidation of a- pinene in pentadecane. The study of the alpha- pinene/ ozone system is the first using organic seed particles to show that the hygroscopicity of the particle does not increase dramatically over the course of the oxidation. No products were detected by Raman spectroscopy for the reaction of benzoate ions with ozone. One product peak, consistent with aqueous carbonate anions, was observed when following the oxidation of fumarate ions by ozone. Product peaks observed in the reaction of ozone with alpha- pinene suggest the formation of new species containing carbonyl groups.

Metal cluster expansion by addition of low valent group 14 reagents: III. Reaction of bis[bis(trimethylsilyl)methyl]tin with M3(CO)12 or its derivatives (M = Fe, Ru, or Os): the crystal and molecular structures of M3(μ-SnR2)2(CO)10(R = CH(SiMe3)2; M = Ru or Os)

The stannylene [SnR2] (R = CH(SiMe3)2) reacts in different ways with the three dodecacarbonyls of the iron triad: [Fe3(CO)12] gives [Fe2(CO)8(μ-SnR2)], [Ru3(CO)12] gives the planar pentametallic cluster [Ru3(CO)10(μ-SnR2)2], for which a full structural analysis is reported, while [Os3(CO)12] fails to react. Different products are also obtained from three nitrile derivatives: [Fe3-(CO)11(MeCN)] gives [Fe2(CO)6(μ-SnR2)2], which has a structure significantly different from that of known Fe2Sn2 clusters, [Ru3(CO)10(MeCN)2] gives the pentametallic cluster described above, while [Os3(CO)10(MeCN)2] gives the isostructural osmium analogue, which shows the unusual feature of a CO group bridging two osmium atoms.

New beetles (Insecta : Coleoptera) from the Berriasian Purbeck Limestone Group, Dorset, UK

Four new beetle species are described from the Lower Cretaceous Purbeck Limestone Group of southern England: Mesogyrus anglicus sp. nov. (Gyrinidae), Coptoclavella purbeckensis sp. nov. (Coptoclavidae), Palaeodytes incompleta sp. nov. (Dytiscidae) and Cretorabus suleatus sp. nov. (Carabidae). The first three taxa were aquatic; the last is terrestrial but may have frequented the margins of water bodies. © 2005 Elsevier Ltd. All rights reserved.

Self-activation and out-group contrast

It has been suggested that higher in-group identifiers primed with an out-group stereotype show contrastive behavioral responses because they activate the in-group, social-self. However, priming the personal-self can lead to contrastive judgments. We investigated whether personal self-activation was also evident for higher identifiers primed with an out-group. An experiment demonstrated that higher identifiers primed with an out-group showed faster responses to self-words than higher identifiers primed with the in-group. This findings suggest that the personal-self is also activated for higher identifiers primed with an out-group, and this self-activation may underlie their contrastive responding.

Plant species and functional group effects on abiotic and microbial soil properties and plant-soil feedback responses in two grasslands

1 Plant species differ in their capacity to influence soil organic matter, soil nutrient availability and the composition of soil microbial communities. Their influences on soil properties result in net positive or negative feedback effects, which influence plant performance and plant community composition. 2 For two grassland systems, one on a sandy soil in the Netherlands and one on a chalk soil in the United Kingdom, we investigated how individual plant species grown in monocultures changed abiotic and biotic soil conditions. Then, we determined feedback effects of these soils to plants of the same or different species. Feedback effects were analysed at the level of plant species and plant taxonomic groups (grasses vs. forbs). 3 In the sandy soils, plant species differed in their effects on soil chemical properties, in particular potassium levels, but PLFA (phospholipid fatty acid) signatures of the soil microbial community did not differ between plant species. The effects of soil chemical properties were even greater when grasses and forbs were compared, especially because potassium levels were lower in grass monocultures. 4 In the chalk soil, there were no effects of plant species on soil chemical properties, but PLFA profiles differed significantly between soils from different monocultures. PLFA profiles differed between species, rather than between grasses and forbs. 5 In the feedback experiment, all plant species in sandy soils grew less vigorously in soils conditioned by grasses than in soils conditioned by forbs. These effects correlated significantly with soil chemical properties. None of the seven plant species showed significant differences between performance in soil conditioned by the same vs. other plant species. 6 In the chalk soil, Sanguisorba minor and in particular Briza media performed best in soil collected from conspecifics, while Bromus erectus performed best in soil from heterospecifics. There was no distinctive pattern between soils collected from forb and grass monocultures, and plant performance could not be related to soil chemical properties or PLFA signatures. 7 Our study shows that mechanisms of plant-soil feedback can depend on plant species, plant taxonomic (or functional) groups and site-specific differences in abiotic and biotic soil properties. Understanding how plant species can influence their rhizosphere, and how other plant species respond to these changes, will greatly enhance our understanding of the functioning and stability of ecosystems.

An adaptive approach to implementing bivariate group sequential clinical trial designs

In clinical trials, situations often arise where more than one response from each patient is of interest; and it is required that any decision to stop the study be based upon some or all of these measures simultaneously. Theory for the design of sequential experiments with simultaneous bivariate responses is described by Jennison and Turnbull (Jennison, C., Turnbull, B. W. (1993). Group sequential tests for bivariate response: interim analyses of clinical trials with both efficacy and safety endpoints. Biometrics 49:741-752) and Cook and Farewell (Cook, R. J., Farewell, V. T. (1994). Guidelines for monitoring efficacy and toxicity responses in clinical trials. Biometrics 50:1146-1152) in the context of one efficacy and one safety response. These expositions are in terms of normally distributed data with known covariance. The methods proposed require specification of the correlation, ρ between test statistics monitored as part of the sequential test. It can be difficult to quantify ρ and previous authors have suggested simply taking the lowest plausible value, as this will guarantee power. This paper begins with an illustration of the effect that inappropriate specification of ρ can have on the preservation of trial error rates. It is shown that both the type I error and the power can be adversely affected. As a possible solution to this problem, formulas are provided for the calculation of correlation from data collected as part of the trial. An adaptive approach is proposed and evaluated that makes use of these formulas and an example is provided to illustrate the method. Attention is restricted to the bivariate case for ease of computation, although the formulas derived are applicable in the general multivariate case.

A practical comparison of group-sequential and adaptive designs

Sequential methods provide a formal framework by which clinical trial data can be monitored as they accumulate. The results from interim analyses can be used either to modify the design of the remainder of the trial or to stop the trial as soon as sufficient evidence of either the presence or absence of a treatment effect is available. The circumstances under which the trial will be stopped with a claim of superiority for the experimental treatment, must, however, be determined in advance so as to control the overall type I error rate. One approach to calculating the stopping rule is the group-sequential method. A relatively recent alternative to group-sequential approaches is the adaptive design method. This latter approach provides considerable flexibility in changes to the design of a clinical trial at an interim point. However, a criticism is that the method by which evidence from different parts of the trial is combined means that a final comparison of treatments is not based on a sufficient statistic for the treatment difference, suggesting that the method may lack power. The aim of this paper is to compare two adaptive design approaches with the group-sequential approach. We first compare the form of the stopping boundaries obtained using the different methods. We then focus on a comparison of the power of the different trials when they are designed so as to be as similar as possible. We conclude that all methods acceptably control type I error rate and power when the sample size is modified based on a variance estimate, provided no interim analysis is so small that the asymptotic properties of the test statistic no longer hold. In the latter case, the group-sequential approach is to be preferred. Provided that asymptotic assumptions hold, the adaptive design approaches control the type I error rate even if the sample size is adjusted on the basis of an estimate of the treatment effect, showing that the adaptive designs allow more modifications than the group-sequential method.