A modeled comparison of direct and food web-mediated impacts of common pesticides on Pacific salmon

Pesticide residues in juvenile Chinook salmon and prey items of the Sacramento River watershed, California - A comparison of riverine and floodplain habitats

Juvenile Chinook salmon (Oncorhynchus tshawytscha) of the Sacramento River system encounter many anthropogenically-induced stressors while rearing and migrating to the Pacific Ocean. Located in a prominent agricultural region, the watershed serves as a source of notable contaminants including pesticides. Salmon rearing in riverine and floodplain areas are potentially exposed to these compounds via dietary exposure, which can vary based on selected food webs. Previous studies have suggested that juvenile Chinook salmon rearing in riverine and floodplain environments of the Sacramento River watershed are characterized by different dietary preferences, with potential for contrasting pesticide exposure between habitats. To examine the potential for pesticide exposure, juvenile Chinook salmon and known dietary items were collected in the mainstem Sacramento River and an adjacent floodplain, the Yolo Bypass, in 2019 and 2020, and analyzed for 33 pesticides, including degradates and isomers. Organochlorine pesticides including the DDX group (p,p'-DDT, p,p'-DDD and p,p'-DDE) were prevalent in all examined biota. There was a significantly greater number of total pesticide detections across all classes in zooplankton compared to macroinvertebrates, coupled with higher bifenthrin concentrations in zooplankton across regions and years, which may indicate different exposure potential depending on fish dietary preferences. Detection frequencies and concentrations of organochlorines were higher in prey items during flooding than in drought conditions, suggesting resuspension of legacy compounds. Significantly higher concentrations of organochlorines were recorded in floodplain rearing fish compared to the Sacramento River. These findings suggest that within these habitats, juvenile Chinook salmon feeding primarily on zooplankton within the water column may be exposed to a greater range of pesticides than those feeding on benthic macroinvertebrates, and that the benefits of floodplain rearing may come at a cost of increased organochlorine exposure.

Passive sampling and ecohydrologic modeling to investigate pesticide surface water loading in the Zollner Creek watershed, Oregon, USA

In the U.S. Pacific Northwest and California contaminants entering surface water may harm Endangered Species Act (ESA) listed salmonid species and consequently there is ongoing concern regarding agricultural practices and resulting pesticide surface water loading may adversely impact salmonid species, their food web, and habitat. Characterizing pesticide exposure in surface water at the watershed scale and beyond is challenging due to uncertainty regarding pesticide use practices and sparse monitoring data. We report here a 2-year continuous deployment of passive sampling devices (PSDs) for monitoring of pesticides in surface water at the outflow of the Zollner Creek watershed located within the Willamette Basin, Oregon, USA. This watershed is predominately agricultural and within the geographic range of two ESA listed Pacific salmonid species. Grab and passive sampling monitoring data were used to evaluate the performance of a probabilistic application of the Soil and Water Assessment Tool (SWAT), a physically based process model which integrates institutional and local knowledge and expertise to investigate the relationship between land use practices and pesticide surface water loading at the watershed scale. SWAT estimated pesticide surface water concentrations for the pesticides chlorpyrifos and trifluralin followed temporal trend in PSD monitoring results and the 5th to 95th percentile range of estimated pesticide concentrations based on the probabilistic assessment encompassed 65-76% of the observed PSD concentrations. Evaluation of model estimates for metolachlor in surface water was challenged by insufficient publicly available grab sample monitoring data. A process to estimate pesticide surface water concentrations on biologically relevant time scales and comparison to screening level aquatic life benchmarks is presented. Additionally, model estimates were used to characterize the variance in surface water concentrations in this small hydrologically responsive watershed to determine grab sampling frequency adequate for model evaluation.

The effects of temperature and salinity on the endocrinology in two life stages of juvenile rainbow/steelhead trout (Oncorhynchus mykiss)

The San Francisco Bay Delta is experiencing seasonally warmer waters and saltwater intrusion into historically freshwater ecosystems due to climate change. Steelhead/rainbow trout (Oncorhynchus mykiss) are resident in the Bay-Delta from juvenile development through the smoltification process. Due to increases in sea level, premature seawater (SW) acclimation may co-occur with increased temperatures on pre-smolt juveniles. To evaluate the interactive effects of salinity and temperature on juvenile life stages of salmonids, rainbow trout alevin (3 days post-hatching) were exposed to 13, 16.4 and 19°C for 10 days and then challenged for 24 h to 18 parts per thousand SW. Similarly, fry (4 weeks post-hatching) were exposed to 13, 16.4 and 19°C for 2 weeks (14 days) and then challenged to SW. Estradiol-17β (E₂ ), cortisol, triiodothyronine (T₃ ) and thyroxine (T₄ ) were measured in whole animal homogenates and muscle tissue using enzyme-linked immunosorbent assays. Transcripts of gill Na⁺ /K⁺ ATPase β (NKAα1b), brain growth hormone I (gh1) and brain gonadotropin-releasing hormone receptor 2 (gnrh2) were also measured. Alevin exhibited a significant temperature-dependent decrease in survival, and fry showed a temperature-dependent decrease in condition factor. The gene expression of NKAα1b, gh1 and gnrh2 significantly decreased in all SW-challenged alevin, and a significant decrease in gnrh2 expression was observed in fry with temperature. Alevin T₃ and T₄ concentrations significantly increased with increasing temperature. There was a temperature-dependent increase in E₂ of fry but not of alevin. The results of this study demonstrate that increasing temperature and SW exposure may adversely affect the survival and SW acclimation of alevin and fry stages of salmonids and that the tolerances of younger juvenile stages should be considered when assessing the response of salmonid populations to climate change stressors.

The effects of combinations of limited ration and diazinon exposure on acetylcholinesterase activity, growth and reproduction in Oryzias latipes, the Japanese medaka

Environmental contamination can negatively impact fish populations. In addition to acute toxicity leading to death, toxicants can reduce fish growth and lower reproduction. The potential for adverse population level effects of environmental contaminants are estimated to conduct risk assessments from laboratory toxicity tests that most often measure apical endpoints related to growth, survival and reproduction. The relationships between these effect endpoints are being evaluated to predict shifts in fish population demography better after exposure to environmental toxicants. Environmental contaminants can also affect fish populations indirectly by reducing prey biomass. However, estimating the magnitude of the combined effects of prey reduction and direct toxicity is difficult and rarely attempted. Here we describe a toxicity test designed to estimate the effect on Japanese medaka of both reduced food and chronic exposure to diazinon, an acetylcholinesterase inhibiting organophosphate pesticide. Fish were reared with limited food ration and/or diazinon exposure through a full life cycle to assess possible interactions between the two stressors in their effects on growth and reproduction. Diazinon exposure (10 or 20 μg/L), reduced ration (50% and 25% of ad libitum), or combinations of both lowered growth rates and reproductive output of Japanese medaka. In addition, growth and reproduction alone were modeled, and then various relationships between the two stressors (diazinon and ration) and how they relate to growth and reproduction were modeled.

Legacy habitat contamination as a limiting factor for Chinook salmon recovery in the Willamette Basin, Oregon, USA

In the western United States, the long-term recovery of many Pacific salmon populations is inextricably linked to freshwater habitat quality. Industrial activities from the past century have left a legacy of pollutants that persist, particularly near working waterfronts. The adverse impacts of these contaminants on salmon health have been studied for decades, but the population-scale consequences of chemical exposure for salmonids are still poorly understood. We estimated acute and delayed mortality rates for seaward migrating juvenile Chinook salmon that feed and grow in a Superfund-designated area in the Lower Willamette River in Portland, Oregon. We combined previous, field-collected exposure data for juvenile Chinook salmon together with reduced growth and disease resistance data from earlier field and laboratory studies. Estimates of mortality were then incorporated into a life cycle model to explore chemical habitat-related fish loss. We found that 54% improved juvenile survival—potentially as a result of future remediation activities—could increase adult Chinook salmon population abundance by more than 20%. This study provides a framework for evaluating pollution remediation as a positive driver for species recovery.

Effects of glyphosate formulations on the population dynamics of two freshwater cladoceran species

The general objective of this work is to experimentally assess the effects of acute glyphosate pollution on two freshwater cladoceran species (Daphnia magna and Ceriodaphnia dubia) and to use this information to predict the population dynamics and the potential for recovery of exposed organisms. Five to six concentrations of four formulations of glyphosate (4-Gly) (Eskoba®, Panzer Gold®, Roundup Ultramax® and Sulfosato Touchdown®) were evaluated in both cladoceran species through acute tests and 15-day recovery tests in order to estimate the population dynamics of microcrustaceans. The endpoints of the recovery test were: survival, growth (number of molts), fecundity, and the intrinsic population growth rate (r). A matrix population model (MPM) was applied to r of the survivor individuals of the acute tests, followed by a Monte Carlo simulation study. Among the 4-Gly tested, Sulfosato Touchdown® was the one that showed higher toxicity, and C. dubia was the most sensitive species. The Monte Carlo simulation study showed an average value of λ always <1 for D. magna, indicating that its populations would not be able to survive under natural environmental conditions after an acute Gly exposure between 0.25 and 35 a.e. mg L-1. The average value of λ for C. dubia was also <1 after exposure to Roundup Ultramax®: 1.30 and 1.20 for 1.21 and 2.5 mg a.e. L-1,respectively. The combined methodology-recovery tests and the later analysis through MPM with a Monte Carlo simulation study-is proposed to integrate key demographic parameters and predict the possible fate of microcrustacean populations after being exposed to acute 4-Gly contamination events.

Complex mixtures of dissolved pesticides show potential aquatic toxicity in a synoptic study of Midwestern U.S. streams

Aquatic organisms in streams are exposed to pesticide mixtures that vary in composition over time in response to changes in flow conditions, pesticide inputs to the stream, and pesticide fate and degradation within the stream. To characterize mixtures of dissolved-phase pesticides and degradates in Midwestern streams, a synoptic study was conducted at 100 streams during May-August 2013. In weekly water samples, 94 pesticides and 89 degradates were detected, with a median of 25 compounds detected per sample and 54 detected per site. In a screening-level assessment using aquatic-life benchmarks and the Pesticide Toxicity Index (PTI), potential effects on fish were unlikely in most streams. For invertebrates, potential chronic toxicity was predicted in 53% of streams, punctuated in 12% of streams by acutely toxic exposures. For aquatic plants, acute but likely reversible effects on biomass were predicted in 75% of streams, with potential longer-term effects on plant communities in 9% of streams. Relatively few pesticides in water-atrazine, acetochlor, metolachlor, imidacloprid, fipronil, organophosphate insecticides, and carbendazim-were predicted to be major contributors to potential toxicity. Agricultural streams had the highest potential for effects on plants, especially in May-June, corresponding to high spring-flush herbicide concentrations. Urban streams had higher detection frequencies and concentrations of insecticides and most fungicides than in agricultural streams, and higher potential for invertebrate toxicity, which peaked during July-August. Toxicity-screening predictions for invertebrates were supported by quantile regressions showing significant associations for the Benthic Invertebrate-PTI and imidacloprid concentrations with invertebrate community metrics for MSQA streams, and by mesocosm toxicity testing with imidacloprid showing effects on invertebrate communities at environmentally relevant concentrations. This study documents the most complex pesticide mixtures yet reported in discrete water samples in the U.S. and, using multiple lines of evidence, predicts that pesticides were potentially toxic to nontarget aquatic life in about half of the sampled streams.

Mechanistic modeling of insecticide risks to breeding birds in North American agroecosystems

Insecticide usage in the United States is ubiquitous in urban, suburban, and rural environments. There is accumulating evidence that insecticides adversely affect non-target wildlife species, including birds, causing mortality, reproductive impairment, and indirect effects through loss of prey base, and the type and magnitude of such effects differs by chemical class, or mode of action. In evaluating data for an insecticide registration application and for registration review, scientists at the United States Environmental Protection Agency (USEPA) assess the fate of the insecticide and the risk the insecticide poses to the environment and non-target wildlife. Current USEPA risk assessments for pesticides generally rely on endpoints from laboratory based toxicity studies focused on groups of individuals and do not directly assess population-level endpoints. In this paper, we present a mechanistic model, which allows risk assessors to estimate the effects of insecticide exposure on the survival and seasonal productivity of birds known to forage in agricultural fields during their breeding season. This model relies on individual-based toxicity data and translates effects into endpoints meaningful at the population level (i.e., magnitude of mortality and reproductive impairment). The model was created from two existing USEPA avian risk assessment models, the Terrestrial Investigation Model (TIM v.3.0) and the Markov Chain Nest Productivity model (MCnest). The integrated TIM/MCnest model was used to assess the relative risk of 12 insecticides applied via aerial spray to control corn pests on a suite of 31 avian species known to forage in cornfields in agroecosystems of the Midwest, USA. We found extensive differences in risk to birds among insecticides, with chlorpyrifos and malathion (organophosphates) generally posing the greatest risk, and bifenthrin and λ-cyhalothrin (pyrethroids) posing the least risk. Comparative sensitivity analysis across the 31 species showed that ecological trait parameters related to the timing of breeding and reproductive output per nest attempt offered the greatest explanatory power for predicting the magnitude of risk. An important advantage of TIM/MCnest is that it allows risk assessors to rationally combine both acute (lethal) and chronic (reproductive) effects into a single unified measure of risk.

Population modeling for pesticide risk assessment of threatened species-A case study of a terrestrial plant, Boltonia decurrens

Although population models are recognized as necessary tools in the ecological risk assessment of pesticides, particularly for species listed under the Endangered Species Act, their application in this context is currently limited to very few cases. The authors developed a detailed, individual-based population model for a threatened plant species, the decurrent false aster (Boltonia decurrens), for application in pesticide risk assessment. Floods and competition with other plant species are known factors that drive the species' population dynamics and were included in the model approach. The authors use the model to compare the population-level effects of 5 toxicity surrogates applied to B. decurrens under varying environmental conditions. The model results suggest that the environmental conditions under which herbicide applications occur may have a higher impact on populations than organism-level sensitivities to an herbicide within a realistic range. Indirect effects may be as important as the direct effects of herbicide applications by shifting competition strength if competing species have different sensitivities to the herbicide. The model approach provides a case study for population-level risk assessments of listed species. Population-level effects of herbicides can be assessed in a realistic and species-specific context, and uncertainties can be addressed explicitly. The authors discuss how their approach can inform the future development and application of modeling for population-level risk assessments of listed species, and ecological risk assessment in general. Environ Toxicol Chem 2017;36:480-491. © 2016 SETAC.

Assessing the risks of pesticides to threatened and endangered species using population modeling: A critical review and recommendations for future work

United States legislation requires the US Environmental Protection Agency to ensure that pesticide use does not cause unreasonable adverse effects on the environment, including species listed under the Endangered Species Act (ESA; hereafter referred to as listed species). Despite a long history of population models used in conservation biology and resource management and a 2013 report from the US National Research Council recommending their use, application of population models for pesticide risk assessments under the ESA has been minimal. The pertinent literature published from 2004 to 2014 was reviewed to explore the availability of population models and their frequency of use in listed species risk assessments. The models were categorized in terms of structure, taxonomic coverage, purpose, inputs and outputs, and whether the models included density dependence, stochasticity, or risk estimates, or were spatially explicit. Despite the widespread availability of models and an extensive literature documenting their use in other management contexts, only 2 of the approximately 400 studies reviewed used population models to assess the risks of pesticides to listed species. This result suggests that there is an untapped potential to adapt existing models for pesticide risk assessments under the ESA, but also that there are some challenges to do so for listed species. Key conclusions from the analysis are summarized, and priorities are recommended for future work to increase the usefulness of population models as tools for pesticide risk assessments. Environ Toxicol Chem 2016;35:1904-1913. © 2016 SETAC.

Year-Round Monitoring of Contaminants in Neal and Rogers Creeks, Hood River Basin, Oregon, 2011-12, and Assessment of Risks to Salmonids

Pesticide presence in streams is a potential threat to Endangered Species Act listed salmonids in the Hood River basin, Oregon, a primarily forested and agricultural basin. Two types of passive samplers, polar organic chemical integrative samplers (POCIS) and semipermeable membrane devices (SPMDs), were simultaneously deployed at four sites in the basin during Mar. 2011-Mar. 2012 to measure the presence of pesticides, polybrominated diphenyl ethers (PBDEs), and polychlorinated biphenyls (PCBs). The year-round use of passive samplers is a novel approach and offers several new insights. Currently used pesticides and legacy contaminants, including many chlorinated pesticides and PBDEs, were present throughout the year in the basin's streams. PCBs were not detected. Time-weighted average water concentrations for the 2-month deployment periods were estimated from concentrations of chemicals measured in the passive samplers. Currently used pesticide concentrations peaked during spring and were detected beyond their seasons of expected use. Summed concentrations of legacy contaminants in Neal Creek were highest during July-Sept., the period with the lowest streamflows. Endosulfan was the only pesticide detected in passive samplers at concentrations exceeding Oregon or U.S. Environmental Protection Agency water-quality thresholds. A Sensitive Pesticide Toxicity Index (SPTI) was used to estimate the relative acute potential toxicity among sample mixtures. The acute potential toxicity of the detected mixtures was likely greater for invertebrates than for fish and for all samples in Neal Creek compared to Rogers Creek, but the indices appear to be low overall (<0.1). Endosulfans and pyrethroid insecticides were the largest contributors to the SPTIs for both sites. SPTIs of some discrete (grab) samples from the basin that were used for comparison exceeded 0.1 when some insecticides (azinphos methyl, chlorpyrifos, malathion) were detected at concentrations near or exceeding acute water-quality thresholds. Early life stages and adults of several sensitive fish species, including salmonids, are present in surface waters of the basin throughout the year, including during periods of peak estimated potential toxicity. Based on these data, direct toxicity to salmonids from in-stream pesticide exposure is unlikely, but indirect impacts (reduced fitness due to cumulative exposures or negative impacts to invertebrate prey populations) are unknown.

Targeted Gene Expression in Zebrafish Exposed to Chlorpyrifos-Oxon Confirms Phenotype-Specific Mechanisms Leading to Adverse Outcomes

Zebrafish models for mild, moderate, and severe acute organophosphorus poisoning were previously developed by exposing zebrafish larvae to chlopyrifos-oxon. The phenotype of these models was characterized at several levels of biological organization. Oxidative stress and mitochondrial dysfunction were found to be involved in the development of the more severe phenotype. Here we used targeted gene expression to understand the dose-responsiveness of those two pathways and their involvement on generating the different zebrafish models. As the severe phenotype is irreversible after only 3 h of exposure, we also analyzed the response of the oxidative stress pathway at 3 and 24 h. Some of the genes related to oxidative stress were already differentially expressed at 3 h. There was an increase in differentially expressed genes related to both oxidative stress and mitochondrial function from the more mild to the more severe phenotype, suggesting the involvement of these mechanisms in increasing phenotype severity. Temporal data suggest that peroxynitrite leading to lipid peroxidation might be involved in phenotype transition and irreversibility.