Linking pesticide exposure and spatial dynamics: An individual-based model of wood mouse (Apodemus sylvaticus) populations in agricultural landscapes

The wood mouse is a common and abundant species in agricultural landscape and is a focal species in pesticide risk assessment. Empirical studies on the ecology of the wood mouse have provided sufficient information for the species to be modelled mechanistically. An individual-based model was constructed to explicitly represent the locations and movement patterns of individual mice. This together with the schedule of pesticide application allows prediction of the risk to the population from pesticide exposure. The model included life-history traits of wood mice as well as typical landscape dynamics in agricultural farmland in the UK. The model obtains a good fit to the available population data and is fit for risk assessment purposes. It can help identify spatio-temporal situations with the largest potential risk of exposure and enables extrapolation from individual-level endpoints to population-level effects. Largest risk of exposure to pesticides was found when good crop growth in the “sink” fields coincided with high “source” population densities in the hedgerows. Keywords: Population dynamics, Pesticides, Ecological risk assessment, Habitat choice, Agent-based model, NetLogo

Incorporating toxicokinetics into an individual-based model for more realistic pesticide exposure estimates: A case study of the wood mouse

The potential risk of agricultural pesticides to mammals typically depends on internal concentrations within individuals, and these are determined by the amount ingested and by absorption, distribution, metabolism, and excretion (ADME). Pesticide residues ingested depend, amongst other things, on individual spatial choices which determine how much and when feeding sites and areas of pesticide application overlap, and can be calculated using individual-based models (IBMs). Internal concentrations can be calculated using toxicokinetic (TK) models, which are quantitative representations of ADME processes. Here we provide a population model for the wood mouse (Apodemus sylvaticus) in which TK submodels were incorporated into an IBM representation of individuals making choices about where to feed. This allows us to estimate the contribution of individual spatial choice and TK processes to risk. We compared the risk predicted by four IBMs: (i) “AllExposed-NonTK”: assuming no spatial choice so all mice have 100% exposure, no TK, (ii) “AllExposed-TK”: identical to (i) except that the TK processes are included where individuals vary because they have different temporal patterns of ingestion in the IBM, (iii) “Spatial-NonTK”: individual spatial choice, no TK, and (iv) “Spatial-TK”: individual spatial choice and with TK. The TK parameters for hypothetical pesticides used in this study were selected such that a conventional risk assessment would fail. Exposures were standardised using risk quotients (RQ; exposure divided by LD50 or LC50). We found that for the exposed sub-population including either spatial choice or TK reduced the RQ by 37–85%, and for the total population the reduction was 37–94%. However spatial choice and TK together had little further effect in reducing RQ. The reasons for this are that when the proportion of time spent in treated crop (PT) approaches 1, TK processes dominate and spatial choice has very little effect, and conversely if PT is small spatial choice dominates and TK makes little contribution to exposure reduction. The latter situation means that a short time spent in the pesticide-treated field mimics exposure from a small gavage dose, but TK only makes a substantial difference when the dose was consumed over a longer period. We concluded that a combined TK-IBM is most likely to bring added value to the risk assessment process when the temporal pattern of feeding, time spent in exposed area and TK parameters are at an intermediate level; for instance wood mice in foliar spray scenarios spending more time in crop fields because of better plant cover.

Neonicotinoid pesticide exposure impairs crop pollination services provided by bumblebees

Recent concern over global pollinator declines has led to considerable research on the effects of pesticides on bees1, 2, 3, 4, 5. Although pesticides are typically not encountered at lethal levels in the field, there is growing evidence indicating that exposure to field-realistic levels can have sublethal effects on bees, affecting their foraging behaviour1, 6, 7, homing ability8, 9 and reproductive success2, 5. Bees are essential for the pollination of a wide variety of crops and the majority of wild flowering plants10, 11, 12, but until now research on pesticide effects has been limited to direct effects on bees themselves and not on the pollination services they provide. Here we show the first evidence to our knowledge that pesticide exposure can reduce the pollination services bumblebees deliver to apples, a crop of global economic importance. Bumblebee colonies exposed to a neonicotinoid pesticide provided lower visitation rates to apple trees and collected pollen less often. Most importantly, these pesticide-exposed colonies produced apples containing fewer seeds, demonstrating a reduced delivery of pollination services. Our results also indicate that reduced pollination service delivery is not due to pesticide-induced changes in individual bee behaviour, but most likely due to effects at the colony level. These findings show that pesticide exposure can impair the ability of bees to provide pollination services, with important implications for both the sustained delivery of stable crop yields and the functioning of natural ecosystems.

Sublethal neonicotinoid insecticide exposure reduces solitary bee reproductive success

1. Pollinating insects provide crucial and economically important ecosystem services to crops and wild plants, but pollinators, particularly bees, are globally declining as a result of various driving factors, including the prevalent use of pesticides for crop protection. Sublethal pesticide exposure negatively impacts numerous pollinator lifehistory traits, but its influence on reproductive success remains largely unknown. Such information is pivotal, however, to our understanding of the long-term effects on population dynamics. 2 We investigated the influence of field-realistic trace residues of the routinely used neonicotinoid insecticides thiamethoxam and clothianidin in nectar substitutes on the entire life-time fitness performance of the red mason bee Osmia bicornis. 3 We show that chronic, dietary neonicotinoid exposure has severe detrimental effects on solitary bee reproductive output. Neonicotinoids did not affect adult bee mortality; however, monitoring of fully controlled experimental populations revealed that sublethal exposure resulted in almost 50% reduced total offspring production and a significantly male-biased offspring sex ratio. 4 Our data add to the accumulating evidence indicating that sublethal neonicotinoid effects on non-Apis pollinators are expressed most strongly in a rather complex, fitness-related context. Consequently, to fully mitigate long-term impacts on pollinator population dynamics, present pesticide risk assessments need to be expanded to include whole life-cycle fitness estimates, as demonstrated in the present study using O. bicornis as a model.

Plastic films for polytunnels can prolong the effective residual life of cypermethrin to over 6 months

Synthetic pyrethroid insecticides are degraded almost entirely by ultraviolet (UV)-catalysed oxidation. A bioassay using the beetle Tribolium confusum duVal caged on bandages soaked in 0.04% a.i. cypermethrin showed large differences in residual insecticide-life under three plastic films available for cladding polytunnels. Cypermethrin exposed to a UV film that transmitted 70% of UVB and 80% of UVA killed all beetles for 8 weeks, compared to only 3 weeks for cypermethrin exposed in a clear plastic envelope. Cypermethrin under a UV-absorbing film that reduced the transmission of UVB and UVA to 14% and 50%, respectively, gave a complete kill for 17 weeks. Reducing the transmission of UVB to virtually zero, and that of UVA to only 3%, using a UV-opaque film prolonged the effective life of the cypermethrin residue to 26 weeks, and some beetles were still killed for a further 11 weeks. Even after this time, beetles exposed to cypermethrin from the UV-opaque treatment were still affected by the insecticide, and only showed near-normal mobility after 24 months of pesticide exposure to the UV-opaque film. These results have implications for the recommended intervals between cypermethrin treatment and crop harvest, and on the time of introduction of insect-based biological control agents, when UV-opaque films are used in commercial horticulture.

A novel method for assessing risks to pollinators from plant protection products using honeybees as a model species

Pollination is one of the most important ecosystem services in agroecosystems and supports food production. Pollinators are potentially at risk being exposed to pesticides and the main route of exposure is direct contact, in some cases ingestion, of contaminated materials such as pollen, nectar, flowers and foliage. To date there are no suitable methods for predicting pesticide exposure for pollinators, therefore official procedures to assess pesticide risk are based on a Hazard Quotient. Here we develop a procedure to assess exposure and risk for pollinators based on the foraging behaviour of honeybees (Apis mellifera) and using this species as indicator representative of pollinating insects. The method was applied in 13 European field sites with different climatic, landscape and land use characteristics. The level of risk during the crop growing season was evaluated as a function of the active ingredients used and application regime. Risk levels were primarily determined by the agronomic practices employed (i.e. crop type, pest control method, pesticide use), and there was a clear temporal partitioning of risks through time. Generally the risk was higher in sites cultivated with permanent crops, such as vineyard and olive, than in annual crops, such as cereals and oil seed rape. The greatest level of risk is generally found at the beginning of the growing season for annual crops and later in June–July for permanent crops.

An energy budget agent-based model of earthworm populations and its application to study the effects of pesticides

Earthworms are important organisms in soil communities and so are used as model organisms in environmental risk assessments of chemicals. However current risk assessments of soil invertebrates are based on short-term laboratory studies, of limited ecological relevance, supplemented if necessary by site-specific field trials, which sometimes are challenging to apply across the whole agricultural landscape. Here, we investigate whether population responses to environmental stressors and pesticide exposure can be accurately predicted by combining energy budget and agent-based models (ABMs), based on knowledge of how individuals respond to their local circumstances. A simple energy budget model was implemented within each earthworm Eisenia fetida in the ABM, based on a priori parameter estimates. From broadly accepted physiological principles, simple algorithms specify how energy acquisition and expenditure drive life cycle processes. Each individual allocates energy between maintenance, growth and/or reproduction under varying conditions of food density, soil temperature and soil moisture. When simulating published experiments, good model fits were obtained to experimental data on individual growth, reproduction and starvation. Using the energy budget model as a platform we developed methods to identify which of the physiological parameters in the energy budget model (rates of ingestion, maintenance, growth or reproduction) are primarily affected by pesticide applications, producing four hypotheses about how toxicity acts. We tested these hypotheses by comparing model outputs with published toxicity data on the effects of copper oxychloride and chlorpyrifos on E. fetida. Both growth and reproduction were directly affected in experiments in which sufficient food was provided, whilst maintenance was targeted under food limitation. Although we only incorporate toxic effects at the individual level we show how ABMs can readily extrapolate to larger scales by providing good model fits to field population data. The ability of the presented model to fit the available field and laboratory data for E. fetida demonstrates the promise of the agent-based approach in ecology, by showing how biological knowledge can be used to make ecological inferences. Further work is required to extend the approach to populations of more ecologically relevant species studied at the field scale. Such a model could help extrapolate from laboratory to field conditions and from one set of field conditions to another or from species to species.

Simulating human and environmental exposure from hand-held knapsack pesticide application: Be-WetSpa-Pest, an integrative, spatially explicit modeling approach

This paper presents an integrative and spatially explicit modeling approach for analyzing human and environmental exposure from pesticide application of smallholders in the potato producing Andean region in Colombia. The modeling approach fulfills the following criteria: (i) it includes environmental and human compartments; (ii) it contains a behavioral decision-making model for estimating the effect of policies on pesticide flows to humans and the environment; (iii) it is spatially explicit; and (iv) it is modular and easily expandable to include additional modules, crops or technologies. The model was calibrated and validated for the Vereda La Hoya and was used to explore the effect of different policy measures in the region. The model has moderate data requirements and can be adapted relatively easy to other regions in developing countries with similar conditions.

Gene transcription in Daphnia magna: effects of acute exposure to a carbamate insecticide and an acetanilide herbicide.

Daphnia magna is a key invertebrate in the freshwater environment and is used widely as a model in ecotoxicological measurements and risk assessment. Understanding the genomic responses of D. magna to chemical challenges will be of value to regulatory authorities worldwide. Here we exposed D. magna to the insecticide methomyl and the herbicide propanil to compare phenotypic effects with changes in mRNA expression levels. Both pesticides are found in drainage ditches and surface water bodies standing adjacent to crops. Methomyl, a carbamate insecticide widely used in agriculture, inhibits acetylcholinesterase, a key enzyme in nerve transmission. Propanil, an acetanilide herbicide, is used to control grass and broad-leaf weeds. The phenotypic effects of single doses of each chemical were evaluated using a standard immobilisation assay. Immobilisation was linked to global mRNA expression levels using the previously estimated 48h-EC(1)s, followed by hybridization to a cDNA microarray with more than 13,000 redundant cDNA clones representing >5000 unique genes. Following exposure to methomyl and propanil, differential expression was found for 624 and 551 cDNAs, respectively (one-way ANOVA with Bonferroni correction, Ppesticides promoted transcriptional changes in energy metabolism (e.g., mitochondrial proteins, ATP synthesis-related proteins), moulting (e.g., chitin-binding proteins, cuticular proteins) and protein biosynthesis (e.g., ribosomal proteins, transcription factors). Methomyl induced the transcription of genes involved in specific processes such as ion homeostasis and xenobiotic metabolism. Propanil highly promoted haemoglobin synthesis and up-regulated genes specifically related to defence mechanisms (e.g., innate immunity response systems) and neuronal pathways. Pesticide-specific toxic responses were found but there is little evidence for transcriptional responses purely restricted to genes associated with the pesticide target site or mechanism of toxicity.

Accumulated metal speciation in earthworm populations with multi-generational exposure to metalliferous soils: cell fractionation and high energy synchrotron analyses

Predicting metal bioaccumulation and toxicity in soil organisms is complicated by site-specific biotic and abiotic parameters. In this study we exploited tissue fractionation and digestion techniques, combined with X-ray absorption spectroscopy (XAS), to investigate the whole-body and subcellular distributions, ligand affinities, and coordination chemistry of accumulated Pb and Zn in field populations of the epigeic earthworm Lumbricus rubellus inhabiting three contrasting metalliferous and two unpolluted soils. Our main findings were (i) earthworms were resident in soils with concentrations of Pb and Zn ranging from 1200 to 27 000 mg kg(-1) and 200 to 34 000 mg kg(-1), respectively; (ii) Pb and Zn primarily accumulated in the posterior alimentary canal in nonsoluble subcellular fractions of earthworms; (iii) site-specific differences in the tissue and subcellular partitioning profiles of populations were observed, with earthworms from a calcareous site partitioning proportionally more Pb to their anterior body segments and Zn to the chloragosome-rich subcellular fraction than their acidic-soil inhabiting counterparts; (iv) XAS indicated that the interpopulation differences in metal partitioning between organs were not accompanied by qualitative differences in ligand-binding speciation, because crystalline phosphate-containing pyromorphite was a predominant chemical species in the whole-worm tissues of all mine soil residents. Differences in metal (Pb, Zn) partitioning at both organ and cellular levels displayed by field populations with protracted histories of metal exposures may reflect their innate ecophysiological responses to essential edaphic variables, such as Ca2+ status. These observations are highly significant in the challenging exercise of interpreting holistic biomarker data delivered by "omic" technologies.

Human exposure modelling for chemical risk assessment: a review of current approaches and research and policy implications

A wide variety of exposure models are currently employed for health risk assessments. Individual models have been developed to meet the chemical exposure assessment needs of Government, industry and academia. These existing exposure models can be broadly categorised according to the following types of exposure source: environmental, dietary, consumer product, occupational, and aggregate and cumulative. Aggregate exposure models consider multiple exposure pathways, while cumulative models consider multiple chemicals. In this paper each of these basic types of exposure model are briefly described, along with any inherent strengths or weaknesses, with the UK as a case study. Examples are given of specific exposure models that are currently used, or that have the potential for future use, and key differences in modelling approaches adopted are discussed. The use of exposure models is currently fragmentary in nature. Specific organisations with exposure assessment responsibilities tend to use a limited range of models. The modelling techniques adopted in current exposure models have evolved along distinct lines for the various types of source. In fact different organisations may be using different models for very similar exposure assessment situations. This lack of consistency between exposure modelling practices can make understanding the exposure assessment process more complex, can lead to inconsistency between organisations in how critical modelling issues are addressed (e.g. variability and uncertainty), and has the potential to communicate mixed messages to the general public. Further work should be conducted to integrate the various approaches and models, where possible and regulatory remits allow, to get a coherent and consistent exposure modelling process. We recommend the development of an overall framework for exposure and risk assessment with common approaches and methodology, a screening tool for exposure assessment, collection of better input data, probabilistic modelling, validation of model input and output and a closer working relationship between scientists and policy makers and staff from different Government departments. A much increased effort is required is required in the UK to address these issues. The result will be a more robust, transparent, valid and more comparable exposure and risk assessment process. (C) 2006 Elsevier Ltd. All rights reserved.

Production risk, pesticide use and GM crop technology in South Africa

Technology involving genetic modification of crops has the potential to make a contribution to rural poverty reduction in many developing countries. Thus far, pesticide-producing Bacillus thuringensis (Bt) varieties of cotton have been the main GM crops under cultivation in developing nations. Several studies have evaluated the farm-level performance of Bt varieties in comparison to conventional ones by estimating production technology, and have mostly found Bt technology to be very successful in raising output and/or reducing pesticide input. However, the production risk properties of this technology have not been studied, although they are likely to be important to risk-averse smallholders. This study investigates the output risk aspects of Bt technology by estimating two 'flexible risk' production function models allowing technology to independently affect the mean and higher moments of output. The first is the popular Just-Pope model and the second is a more general 'damage control' flexible risk model. The models are applied to cross-sectional data on South African smallholders, some of whom used Bt varieties. The results show no evidence that a 'risk-reduction' claim can be made for Bt technology. Indeed, there is some evidence to support the notion that the technology increases output risk, implying that simple (expected) profit computations used in past evaluations may overstate true benefits.

Impacts of a pesticide on pollinator species richness at different spatial scales

Pesticides are an important potential cause of biodiversity and pollinator decline. Little is known about the impacts of pesticides on wild pollinators in the field. Insect pollinators were sampled in an agricultural system in Italy with the aim of detecting the impacts of pesticide use. The insecticide fenitrothion was over 150 times greater in toxicity than other pesticides used in the area, so sampling was set up around its application. Species richness of wild bees, bumblebees and butterflies were sampled at three spatial scales to assess responses to pesticide application: (i) the ‘field’ scale along pesticide drift gradients; (ii) the ‘landscape’ scale sampling in different crops within the area and (iii) the ‘regional’ scale comparing two river basins with contrasting agricultural intensity. At the field scale, the interaction between the application regime of the insecticide and the point in the season was important for species richness. Wild bee species richness appeared to be unaffected by one insecticide application, but declined after two and three applications. At the landscape scale, the species richness of wild bees declined in vine fields where the insecticide was applied, but did not decline in maize or uncultivated fields. At the regional scale, lower bumblebee and butterfly species richness was found in the more intensively farmed basin with higher pesticide loads. Our results suggest that wild bees are an insect pollinator group at particular risk from pesticide use. Further investigation is needed on how the type, quantity and timing of pesticide application impacts pollinators.

Metal bioaccumulation and cellular fractionation in an epigeic earthworm (Lumbricus rubellus): the interactive influences of population exposure histories, site-specific geochemistry and mitochondrial genotype

Subcellular fractionation techniques were used to describe temporal changes (at intervals from T0 to T70 days) in the Pb, Zn and P partitioning profiles of Lumbricus rubellus populations from one calcareous (MDH) and one acidic (MCS) geographically isolated Pb/Zn-mine sites and one reference site (CPF). MDH and MCS individuals were laboratory maintained on their native field soils; CPF worms were exposed to both MDH and MCS soils. Site-specific differences in metal partitioning were found: notably, the putatively metal-adapted populations, MDH and MCS, preferentially partitioned higher proportions of their accumulated tissue metal burdens into insoluble CaPO4-rich organelles compared with naive counterparts, CPF. Thus, it is plausible that efficient metal immobilization is a phenotypic trait characterising metal tolerant ecotypes. Mitochondrial cytochrome oxidase II (COII) genotyping revealed that the populations indigenous to mine and reference soils belong to distinct genetic lineages, differentiated by 13%, with 7 haplotypes within the reference site lineage but fewer (3 and 4, respectively) in the lineage common to the two mine sites. Collectively, these observations raise the possibility that site-related genotype differences could influence the toxico-availability of metals and, thus, represent a potential confounding variable in field-based eco-toxicological assessments.