Microscale patterns of habitat fragmentation and disturbance events as a result of chemical applications: effects on Folsomia candida (Collembola) populations Mattia Meli Annemette Palmqvist, Valery E. Forbes 1
Introduction In assessing soil quality and the ecological risk of contamination by chemicals, homogeneous contamination is often assumed However, soils are a very heterogeneous environment Both contamination and the response of soil organisms can be assumed to be heterogeneous, too Consequences for the exposure of soil organisms and for the extrapolation of risk from the individual to the population level 2
Approach Model organism: Folsomia candida (Willem, 1902), Collembola - Standard test arthropod (OECD-guideline 232, 2009) - Common & widely distributed - Parthenogenetic - Life cycle parameters well known Spatially explicit individual-based population model - individuals sense and avoid contaminated habitat with defined probability - Rules and parameters based on data from the literature - Parameterized and evaluated using Pattern Oriented Modeling (POM) Meli M., Auclerc A., Palmqvist A., Forbes V.E. and Grimm V. Population-level consequences of spatially heterogeneous exposure to heavy metals in soil: an individual-based model of springtails. Ecological Modelling 250, 2013 3
Model contaminant Copper sulphate (CuSO 4 ) - Persistent - Extensively used as fungicide - Most widely distributed pollutant among metals Individual-level toxicity data - Reproduction Fecundity reduction (normalized to control) = 0.2189 ln concentration 0.8743 R 2 = 0.919 - Survival Survival reduction (normalized to control) = 0.0824 ln concentration 0.1366 Sandifer & Hopkin, Chemosphere, 1996: Effects of ph on the toxicity of cadmium, copper, lead and zinc to Folsomia candida Willem, 1902 (Collembola) in a standard laboratory test system R 2 = 0.847 - Egg viability % eggs hatching (normalized to control) = 0.2243 ln concentration + 1.8893 R 2 = 0.932 Xu et al., Insect Science, 2009: Effects of copper, lead and zinc in soil on egg development and hatching of Folsomia candida - Avoidance % avoidance = 5.7475 ln concentration 1.4235 Boiteau et al., Environmental Pollution, 2011: Avoidance tests with Folsomia candida for the assessment of copper contamination in agricultural soils R 2 = 0.926 4
Conceptual model Energy < threshold? FORAGE Survive? DIE AGIG DESIT-DEPEDET OVIPOSITIO Age maturation time? Time since last reproduction = Time between broods? REPRODUCE Eggs ewborn juveniles HATCH Age hatching time 5
Movement & Avoidance behaviour Food heterogeneously distributed: typically, 10 % of grid-cells are food sources Sense food? Turn random angle Turn towards food cell Cell in front of you contaminated? Avoid? Move 1 step Turn toward an uncontaminated cell Energy spent for moving > threshold? OR Found food? Move 1 step and increase exposure counter EXIT 6
Sensitivity analysis Sensitivity analysis of parameters not directly determined from literature Directly determined Sensitive (inversely determined) Insensitive (initial values maintained) 7
Pattern oriented parameterization Excess food Slightly limiting food Limiting food Log population size 1000 100 10 1 0 20 40 60 80 100 120 0 20 40 60 80 100 120 Days 0 20 40 60 80 100 120 Observed data from Usher et al.. Oecologia, 1971: Studies on populations of Folsomia candida (Insecta: Collembola). The productivity of populations in relation to food and exploitation observed simulated 95% CL Resulting mean density Range Estimated r Range Observed 463,21 207,62 774,67 0,178 0,166 0,199 Simulated 548,6 442 670 0,163 0,158 0,175 Observed data from Seifert et al., J Kans Entomol Soc, 1979: estimation of the intrinsic rate of increase (r) and density-dependent population size of F. candida cultures 8
Model performance 28-days exposure in artificial soil 20 ± 1 C temperature Mean ± Standard Error of mean for 4 replicates Instantaneous population growth rate 0,10 0,08 0,06 0,04 0,02 0,00-0,02-0,04-0,06-0,08-0,10 0 12,5 50 200 800 3200 12800 Concentration (mg kg -1 ) observed simulated Observed data from Herbert et al., Ecotoxicology and Environmental Safety, 2004: Comparison of instantaneous rate of population increase and critical-effect estimates in Folsomia candida exposed to four toxicants 9
Simulation scenarios - Midpoint displacement algorithm used to generate fractal landscapes - 3 degrees of spatial autocorrelation of contaminated habitat (H) - Varying % of uncontaminated habitat (0-30%) H=0 H=0.72 H=1 - Simulation of a 30 day drought period - Progressive reduction of soil relative humidity - Effects of drought on F. candida derived from empirical observations Waagner et al, Soil biology & biochemistry, 2011: Recovery of reproduction after drought in the soil living Folsomia candida (Collembola) Relative humidity (%) 100 99,5 99 98,5 98 97,5 97 96,5 96 95,5 0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 10 Time since beginning of drought period (days)
Results - Constant temperature - Initial population: 100 individuals (adults and juveniles) - Varying % of uncontaminated habitat (0-30%) - Same concentration (4500 mg Cu kg -1 ) 6000 With avoidance Without avoidance Final abundance (individuals m -2 ) 5000 4000 3000 2000 1000 0 0-2 2-4 4-6 6-8 H=0 H=0.72 H=1 8-10 10-12 12-14 14-16 16-18 18-20 % clean habitat 20-22 22-24 24-26 26-28 28-30 Error bars: Standard Error of mean for 10 replicates 0-2 2-4 4-6 6-8 H=0 H=0.72 H=1 8-10 10-12 12-14 % clean habitat 14-16 16-18 18-20 20-22 22-24 24-26 26-28 28-30 11
Results -20% clean habitat - Varying degrees of spatial autocorrelation of contaminated habitat - Same concentration (4500 mg Cu kg -1 ) - Drought period from day 200 to day 230 Abundance (individuals m- 2 ) 8000 6000 4000 2000 With avoidance H=0 H=0.72 H=1 Without avoidance 0 0 50 100 150 200 250 300 0 50 100 150 200 250 300 Mean ± Standard Error of mean for 10 replicates Time (days) 12
Conclusions Application method of a pesticide influences its distribution in soil. This can affect the impact on soil-dwelling organisms: scattered or clumped distributions of toxicant in soil do not have the same impact on F. candida populations. Avoidance should also be taken into account Modelling can add ecological realism to ERA by incorporating effects of natural stress that occur in the field, such as drought or food scarcity, and that may delay population recovery Models can be used in an exploratory way to predict the effects of different risk management options 13
Thank you For your attention! mattia@ruc.dk ACKOWLEDGEMETS: This research has been financially supported by the European Union under the 7 th Framework Programme (project acronym CREAM, contract number PIT-GA-2009-238148)