Modelling pesticides residues

This work is a contribution to the development of a specific method to assess the presence of residues in agricultural commodities. The following objectives are formulated: to identify and describe main processes in environment — plant exchanges, to build of a model to assess the residue concentration at harvest in agricultural commodities, to understand the functioning of the modelled system, to characterise pesticides used in field crops and identify optimisation potentials in phytosanitary measures. The frame for the methodological developments corresponds to the procedure for the evaluation of the toxicity provided for the Life Cycle Impact Assessment methodology and for the method Impact 2002+. In chapter 2, the methodological procedure for the assessment of human toxicity potential is introduced. First the factors of fate and exposure are described, including the notion of harvest fraction, the amount of substance found in harvest per unit of substance emitted initially in the system, the main result of the present study. Then the effect factors and the framework for impact evaluation are introduced. Chapter 3 describes the principles accounted for the building of the fate model. Wheat crop and a restricted list of substances are chosen for these methodological developments. The model is composed by compartments describing the environment and the plant. Its functioning is based on initial amounts of substance in the source compartments, on transfer rates linking the compartments and on a dynamic evolution as a function of time between the treatment and the harvest. Air, soil and formulation deposit on plant are the primer compartments receiving the treated substance. Each transport is described by a transfer rate accounting for the process and for the equilibrium partitioning between the two exchanging compartments. Degradation of substance and plant growth are additional processes considered. Each compartment is described by a linear differential equation for the variation of mass accumulating and dissipating. Their assembly builds the model solved as a function of time. This exact resolution is complemented by additional tools to better understand the system functioning and to provide further approximations of the results: the system is simplified into subsystems describing the source and the receiving plant compartment and analytically solved using interpretable equations. Chapter 4 describes and discusses all transport and dissipation processes determining the fate of the substance in the limits of the system. The recent publications concerning the understanding and the modelling of pesticide transfer from formulation deposit on plant through the cuticular membranes give new possibilities to model pesticide fate and to better account for the direct applications on the plant. In chapter 5, the core model is first applied and its functioning analysed. The low availability and partly unsatisfying quality of data for pesticides description is a main complication for the methodology application: the lack of data for the half-life of the substance in the plant especially leads to a strong extrapolation for this determinant factor. A large difference is observed between early and late applied pesticides with respectively a major release to soil or a release to formulation deposit on plant surface. The initial transport processes quickly distribute the substance in the system. Once each plant compartment has accumulated residues up to a maximum amount, a dissipation phase occurs. The duration of these periods is determinant for the level of residue in harvest. The soil is a determinant source for long term evolutions of the system, for soon applied substances with low degradation rate. The half-life of substance deposited on plant is equal to a few days, but the transfer is fast from formulation deposit to the inner plant, where degradation is generally much slower. The accumulation of substance from the air is mostly negligible. The sum of the subsystems gives an approximation of the total system, useful for interpretation. The possibility to simplify the subsystem by ignoring the transfer back from receiving to the source compartment underlines the low contribution of these transfers in the functioning of the model. An approximated resolution is based on the determination of the maximum accumulated substance and on the subsequent dissipation process. However, an important loss of precision is observed. This approximation is useful for interpretation and for extrapolations. In chapter 6, an evaluation of the model is conducted through a sensitivity and uncertainty analysis. The sensitivity analysis consists of evaluating the effect on the output of a change in an input, on the basis of three complementary approaches: the effect of a fixed change in the input of e.g. 0.1%, the effect of a change specific to the uncertainty of the input and the effect of a change in input value from a minimum to a maximum. The uncertainty of an output is evaluated according to the relative contribution of the confidence factors of the inputs. Results show that the half-lives and the time are the most important factors determining the sensitivity of the system and the propagation of uncertainty. The contribution of the half-life to the confidence factor of the harvest fraction reaches between 30% to 98% of the total uncertainty. The confidence factors of results increase exponentially with the time interval between application and harvest. The role of partition coefficients to the behaviour of substances is highly variable, may be determinant or negligible, with increasing or limiting effect on mobility. Sensitivity and uncertainty for parameters describing the agricultural or environmental system are very variable, but sometimes determinant and so confirmed as essential for the system functioning. Consequently, differences in harvest fractions between substances are only significant if they are high. A first comparison of the computed results with measures of residues obtained by an experiment and with references such as tolerance values lead to a pertinent verification of the overall methodology. Finally, the qualitative comparison with other models underlines the specificities and the originality of the present methodology in particular by comparison with environmental multi-media models running in steady state. In chapter 7, the model is finally applied for an ultimate interpretation. The harvest fractions for more than 100 substances are evaluated. Among all types of substances, low and high levels of residues per treatment are found, representative for the high variability of harvest fractions from 5E-16 for bromoxynil to 7E-03 for tebuconazole sprayed on wheat. The fate process represents the highest source of variation for the toxicity. If the application rate does not explain the high differences in residue level at harvest, the time of application may represent an optimisation potential particularly for late treatments. However, the toxicity needs to account for both fate and effect factors, as only their combination effectively allows to evaluate the toxicity. According to the available list of Human Damage Factors per treatment, problematic substances may be effectively identified and substituted. In chapter 8 answers to questions brought with the objectives bring a conclusion to the study. The appendices include notably the results of harvest fractions and toxicity per unit substance applied, per treatment and per unit cultivated crop area, for the main substances and field crops. A LCA is also presented on the intensity level of wheat production.

[1]  P. Debaeke,et al.  Grain nitrogen content of winter bread wheat (Triticum aestivum L.) as related to crop management and to the previous crop , 1996 .

[2]  Michael Matthies,et al.  Model for uptake of xenobiotics into plants: Validation with bromacil experiments , 1994 .

[3]  D. Kleier Phloem mobility of xenobiotics: I. Mathematical model unifying the weak Acid and intermediate permeability theories. , 1988, Plant physiology.

[4]  M. McLachlan Vegetation-Air Partition Coefficient , 2000 .

[5]  R. Kirkwood Recent developments in our understanding of the plant cuticle as a barrier to the foliar uptake of pesticides , 1999 .

[6]  O. Jolliet,et al.  Harmonisation of Environmental Life Cycle Assessment for Agriculture , 1997 .

[7]  A. Yang,et al.  Study of root uptake and xylem translocation of cinmethylin and related compounds in detopped soybean roots using a pressure chamber technique. , 1990, Plant physiology.

[8]  W. G. Knisel,et al.  GLEAMS: Groundwater Loading Effects of Agricultural Management Systems , 1987 .

[9]  T. Jager,et al.  Modelling the influence of terrestrial vegetation on the environmental fate of xenobiotics , 1998 .

[10]  R. Heijungs,et al.  Environmental life cycle assessment of products , 1992 .

[11]  Bruce Cumings,et al.  A short review , 1983 .

[12]  S. Trapp Dynamic root uptake model for neutral lipophilic organics , 2002, Environmental toxicology and chemistry.

[13]  Frank Brentrup,et al.  Life Cycle Impact assessment of land use based on the hemeroby concept , 2002 .

[14]  Michael S. McLachlan,et al.  Uptake of Airborne Semivolatile Organic Compounds in Agricultural Plants: Field Measurements of Interspecies Variability , 1999 .

[15]  Plant contamination , 1994 .

[16]  D. Rossier L'ECOBILAN, OUTIL DE GESTION ECOLOGIQUE DE L'EXPLOITATION AGRICOLE? , 1999 .

[17]  Peter Baur,et al.  Modelling penetration of plant cuticles by crop protection agents and effects of adjuvants on their rates of penetration , 1994 .

[18]  R. Bromilow,et al.  Principles governing uptake and transport of chemicals , 1995 .

[19]  J. Schnoor,et al.  Predictive Relationships for Uptake of Organic Contaminants by Hybrid Poplar Trees , 1998 .

[20]  J. Schönherr,et al.  Foliar Penetration and Accumulation of Organic Chemicals in Plant Cuticles , 1989 .

[21]  D. J. Strom,et al.  Determining parameters of lognormal distributions from minimal information. , 2000, AIHAJ : a journal for the science of occupational and environmental health and safety.

[22]  K. Racke,et al.  Foliar interception and retention values after pesticide application: a proposal for standardised values for environmental risk assessment , 2002 .

[23]  F. Hsu,et al.  Phloem mobility of xenobiotics VIII. A short review. , 1996, Journal of experimental botany.

[24]  A Newman,et al.  Ranking pesticides by environmental impact. , 1995, Environmental science & technology.

[25]  O. Jolliet,et al.  Assessing Human Health Response in Life Cycle Assessment Using ED10s and DALYs: Part 1—Cancer Effects , 2002, Risk analysis : an official publication of the Society for Risk Analysis.

[26]  M. Knoche,et al.  Spray application factors and plant growth regulator performance: II. Foliar uptake of gibberellic acid and 2,4‐D , 1999 .

[27]  J. Schönherr,et al.  Estimation of path lengths for diffusion of organic compounds through leaf cuticles , 1999 .

[28]  J. Schönherr,et al.  Differences among plant species in cuticular permeabilities and solute mobilities are not caused by differential size selectivities , 1998, Planta.

[29]  R. Schwarzenbach,et al.  Environmental Organic Chemistry , 1993 .

[30]  P. Baur Surfactant effects on cuticular penetration of neutral polar compounds: dependence on humidity and temperature. , 1999, Journal of agricultural and food chemistry.

[31]  Robert S. Boethling,et al.  Handbook of Property Estimation Methods for Chemicals : Environmental Health Sciences , 2000 .

[32]  Jane Bare,et al.  An International Workshop on Life Cycle Impact Assessment Sophistication , 2000 .

[33]  C. T. Chiou,et al.  A partition-limited model for the plant uptake of organic contaminants from soil and water. , 2001, Environmental science & technology.

[34]  O. Jolliet,et al.  Life cycle impact assessment of pesticides on human health and ecosystems , 2002 .

[35]  S. Gyldenkærne,et al.  Ground deposit of pesticides in relation to the cereal canopy density , 1999 .

[36]  M. Abraham,et al.  The use of characteristic volumes to measure cavity terms in reversed phase liquid chromatography , 1987 .

[37]  Michael Matthies,et al.  Chemodynamics and Environmental Modeling , 1998 .

[38]  G. L. Eilrich Tracking the fate of residues from the farm gate to the table , 1991 .

[39]  W. Melander,et al.  Phloem Mobility of Xenobiotics: II. Bioassay Testing of the Unified Mathematical Model. , 1988, Plant physiology.

[40]  Y. Ronen Uncertainty Analysis , 1988 .

[41]  C. Travis,et al.  Bioconcentration of organics in beef, milk, and vegetation. , 1988, Environmental science & technology.

[42]  O. Jolliet,et al.  Defining intake fraction. , 2002, Environmental science & technology.

[43]  Stefan Trapp,et al.  Plant Contamination: Modeling and Simulation of Organic Chemical Processes , 1994 .

[44]  M. Goedkoop,et al.  The Eco-indicator 99, A damage oriented method for Life Cycle Impact Assessment , 1999 .

[45]  F. Wania,et al.  Estimating the influence of forests on the overall fate of semivolatile organic compounds using a multimedia fate model. , 2001, Environmental science & technology.

[46]  D van de Meent,et al.  Priority assessment of toxic substances in life cycle assessment. Part II: assessing parameter uncertainty and human variability in the calculation of toxicity potentials. , 2000, Chemosphere.

[47]  Frank Brentrup,et al.  Application of the Life Cycle Assessment methodology to agricultural production: an example of sugar beet production with different forms of nitrogen fertilisers , 2001 .

[48]  Stefan Schwartz,et al.  Quality Assurance of Exposure Models for Environmental Risk Assessment of Substances , 2000 .

[49]  J. Schönherr,et al.  Accumulation of lipophilic chemicals in plant cuticles: Prediction from octanol/water partition coefficients , 1988 .

[50]  L. Wax,et al.  A Nonlinear Dynamic Simulation Model for Xenobiotic Transport and Whole Plant Allocation Following Foliar Application. I. Conceptual Foundation for Model Development , 2000 .

[51]  Park S. Nobel,et al.  Physicochemical and Environmental Plant Physiology , 1991 .

[52]  O. Jolliet,et al.  Multimedia fate and human intake modeling: spatial versus nonspatial insights for chemical emissions in Western Europe. , 2005, Environmental science & technology.

[53]  D. Beghin,et al.  Genetic differences for nitrogen uptake and nitrogen utilisation efficiencies in winter wheat. , 2000 .

[54]  B. Feil,et al.  Genotypic variation in grain nitrogen content of wheat as affected by mineral nitrogen supply in the soil , 1992 .

[55]  J. Schönherr,et al.  Permeation of lipophilic chemicals across plant cuticles: prediction from partition coefficients and molar volumes , 1988 .

[56]  W. G. Knisel,et al.  Modelling pesticide fate with GLEAMS , 1995 .

[57]  J. Lammel,et al.  Investigations of the energy efficiency of the production of winter wheat and sugar beet in Europe , 1999 .

[58]  L. Schreiber,et al.  Protecting against water loss: analysis of the barrier properties of plant cuticles. , 2001, Journal of experimental botany.

[59]  R. Bromilow,et al.  Relationships between lipophilicity and the distribution of non-ionised chemicals in barley shoots following uptake by the roots , 1983 .

[60]  L. Pussemier,et al.  Model calculations and measurements of uptake and translocation of carbamates by bean plants , 1991 .

[61]  B. T. Grayson,et al.  Partition Coefficients of Active Ingredients between Plant Cuticle and Adjuvants As Related to Rates of Foliar Uptake , 1997 .

[62]  M. Knoche,et al.  Finite dose diffusion studies: I. Characterizing cuticular penetration in a model system using NAA and isolated tomato fruit cuticles. , 2000 .

[63]  M Matthies,et al.  Generic one-compartment model for uptake of organic chemicals by foliar vegetation. , 1995, Environmental science & technology.

[64]  J. Schönherr,et al.  Polydisperse ethoxylated fatty alcohol surfactants as accelerators of cuticular penetration. 2: Separation of effects on driving force and mobility and reversibility of surfactant action , 1999 .

[65]  R. Bromilow,et al.  Relationships between lipophilicity and root uptake and translocation of non-ionised chemicals by barley† , 1982 .

[66]  D. Mackay,et al.  Strategies for including vegetation compartments in multimedia models. , 2001, Chemosphere.

[67]  J. Schönherr,et al.  Modelling Foliar Penetration: Its Role in Optimising Pesticide Delivery , 1999 .

[68]  L. Wax,et al.  A Nonlinear Dynamic Simulation Model for Xenobiotic Transport and Whole Plant Allocation Following Foliar Application. III. Influence of Chemical Properties, Plant Characteristics, and Environmental Parameters on Xenobiotic Absorption and Translocation , 2001 .

[69]  L. Wax,et al.  A nonlinear dynamic simulation model for xenobiotic transport and whole plant allocation following foliar application II. Model validation , 2000 .

[70]  G. H. Willis,et al.  Pesticide persistence on foliage , 1987 .

[71]  O. Jolliet,et al.  Assessing Human Health Response in Life Cycle Assessment Using ED10s and DALYs: Part 2—Noncancer Effects , 2002, Risk analysis : an official publication of the Society for Risk Analysis.

[72]  J. Schönherr,et al.  Relative solubilities of bifenox and 1‐naphthylacetic acid (NAA) in plant cuticles and in selected pure or aqueous glycol additives , 1998 .

[73]  O. Jolliet,et al.  Intake Fraction for Multimedia Pollutants: A Tool for Life Cycle Analysis and Comparative Risk Assessment , 2002, Risk analysis : an official publication of the Society for Risk Analysis.

[74]  Stefan Trapp,et al.  Modelling uptake into roots and subsequent translocation of neutral and ionisable organic compounds , 2000 .

[75]  D. Mackay,et al.  A novel and simple model of the uptake of organic chemicals by vegetation from air and soil. , 1997, Chemosphere.

[76]  Max Henrion,et al.  Uncertainty: A Guide to Dealing with Uncertainty in Quantitative Risk and Policy Analysis , 1990 .

[77]  Gerald Rebitzer,et al.  IMPACT 2002+: A new life cycle impact assessment methodology , 2003 .