A probabilistic approach to chronic effects assessments for listed species in a vernal pool case study.

Ecological risk assessments for potential pesticide impacts on species listed as threatened or endangered must ensure that decisions to grant registration or establish water quality standards will not jeopardize species or their critical habitats. Pesticides are designed to affect pest species via physiological pathways that may be shared by some nontarget species for which toxicity data are usually unavailable, creating a need for robust methods to estimate acute and chronic toxicity with minimal data. We used a unique probabilistic approach to estimate the risk of chronic effects of two organophosphate (OP) pesticides on the vernal pool fairy shrimp Branchinecta lynchi. Acute toxicity estimates were derived from Monte Carlo (MC) sampling of acute toxicity distributions developed from interspecies relationships using surrogate species. Within each MC draw, acute values were divided by an acute to chronic ratio (ACR) sampled from a distribution of ACRs for OP pesticides and invertebrates, producing a distribution of chronic effects concentrations. The estimated exposure concentrations (EECs) were sampled from distributions representing different environmental conditions. Risk was characterized using probability distributions of acute toxicity, ACRs, and EECs in a probabilistic analysis, as well as partial probabilistic variations that used only some distributions whereas some variables were used deterministically. A deterministic risk quotient (RQ) was compared with the results of probabilistic methods to compare the approaches. Risk varied across exposure scenarios and the number of variables that were handled probabilistically, increasing as the number of variables drawn from distributions increased. The magnitude of RQs was not correlated with the probability that EECs would exceed chronic thresholds, and comparison of the two approaches demonstrates the limited interpretability of RQs. Our novel probabilistic approach to estimating chronic risk with minimal data incorporates uncertainty underlying both exposure and effects assessments for listed species. Integr Environ Assess Manag 2024;00:1-13. Published 2024. This article is a U.S. Government work and is in the public domain in the USA.

Estimating dermal contact soil exposure for amphibians.

Chemical exposure estimation through the dermal route is an underemphasized area of ecological risk assessment for terrestrial animals. Currently, there are efforts to create exposure models to estimate doses from this pathway for use in ecological risk assessment. One significant limitation has been insufficient published data to characterize exposure and to support the selection and parameterization of appropriate models, particularly for amphibians in terrestrial habitats. Recent publications measuring pesticide doses to terrestrial-phase amphibians have begun to rectify this situation. We collated and summarized available measurements of terrestrial amphibian dermal exposure to pesticides from 11 studies in which researchers measured tissue concentrations associated with known pesticide experimental application rates. This data set included tissue concentrations in 11 amphibian species and 14 different pesticides. We then compared the results of two screening exposure models that differed based on surface area scaling approaches as a function of body weight (one based on birds as surrogates for amphibians and another amphibian-specific) to the measured tissue residue concentrations. We define a false-negative rate for each screening model as the proportion of amphibians for which the predicted concentration is less than the observed concentration (i.e., underestimate), contrary to the intent of screening models, which are intended to have a bias for higher exposure concentrations. The screening model that uses birds as surrogates did not have any instances where estimated expected avian doses were less than measured amphibian body burdens. When using the amphibian-specific exposure model that corrected for differences between avian and amphibian surface area, measured concentrations were greater than model estimates for 11.3% of the 1158 comparisons. The database of measured pesticide concentrations in terrestrial amphibians is provided for use in calculating bioconcentration factors and for future amphibian dermal exposure model development. Integr Environ Assess Manag 2023;19:9-16. © 2022 SETAC. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

A multi-scale approach for identification of potential pesticide use sites impacting vernal pool critical habitat in California.

Spatially explicit ecological risk assessment (ERA) requires estimating the overlap between chemical and receptor distribution to evaluate the potential impacts of exposure on nontarget organisms. Pesticide use estimation at field level is prone to error due to inconsistencies between ground-reporting and geospatial data coverage; attempts to rectify these inconsistencies have been limited in approach and rarely scaled to multiple crop types. We built upon a previously developed Bayesian approach to combine multiple crop types for a probabilistic determination of field-crop assignments and to examine co-occurrence of critical vernal pool habitats and bifenthrin application within a 5-county area in California (Madera, Merced, Sacramento, San Joaquin, and Stanislaus counties). We incorporated a multi-scale repeated sampling approach with an area constraint to improve the delineation of field boundaries and better capture variability in crop assignments and rotation schemes. After comparing the accuracy of the spatial probabilistic approach to USDA Census of Agriculture crop acreage data, we found our approach allows more specificity in the combination of crop types represented by the potential application area and improves acreage estimates when compared to traditional deterministic approaches. In addition, our multi-scale sampling scheme improved estimates of bifenthrin acreage variability for co-occurrence analysis and allowed for estimates of crop rotations that were previously uncaptured. Our approach could be leveraged for more realistic, spatially resolved exposure and effects models both in and outside of California.

Investigating Vernal Pool Fairy Shrimp Exposure to Organophosphate Pesticides: Implications for Population-Level Risk Assessment.

Vernal pool fairy shrimp, Branchinecta lynchi, is a freshwater crustacean endemic to California and Oregon, including California's Central Valley. B. lynchi is listed as a Federally Threatened species under the US Endangered Species Act, and as a vulnerable species on the IUCN Red List. Threats that may negatively impact vernal pool fairy shrimp populations include pesticide applications to agricultural land use (e.g., agrochemicals such as organophosphate pesticides) and climate changes that impact vernal pool hydrology. Pop-GUIDE (Population model Guidance, Use, Interpretation, and Development for Ecological risk assessment) is a comprehensive tool that facilitates development and implementation of population models for ecological risk assessment and can be used to document the model derivation process. We employed Pop-GUIDE to document and facilitate the development of a population model for investigating impacts of organophosphate pesticides on vernal pool fairy shrimp populations in California's Central Valley. The resulting model could be applied in combination with field assessment and laboratory-based chemical analysis to link effects from pesticide exposure to adverse outcomes in populations across their range. B. lynchi has a unique intra-annual life cycle that is largely dependent upon environmental conditions. Future deployment of this population model should include complex scenarios consisting of multiple stressors, whereby the model is used to examine scenarios that combine chemical stress resulting from exposure to pesticides and climate changes.

Demonstration of an online web services tool incorporating automatic retrieval and comparison of precipitation data.

Input data acquisition and preprocessing is time-consuming and difficult to handle and can have major implications on environmental modeling results. US EPA's Hydrological Micro Services Precipitation Comparison and Analysis Tool (HMS-PCAT) provides a publicly available tool to accomplish this critical task. We present HMS-PCAT's software design and its use in gathering, preprocessing, and evaluating precipitation data through web services. This tool simplifies catchment and point-based data retrieval by automating temporal and spatial aggregations. In a demonstration of the tool, four gridded precipitation datasets (NLDAS, GLDAS, DAYMET, PRISM) and one set of gauge data (NCEI) were retrieved for 17 regions in the United States and evaluated on 1) how well each dataset captured extreme events and 2) how datasets varied by region. HMS-PCAT facilitates data visualizations, comparisons, and statistics by showing the variability between datasets and allows users to explore the data when selecting precipitation datasets for an environmental modeling application.

Comparison and Evaluation of Gridded Precipitation Datasets in a Kansas Agricultural Watershed Using SWAT.

Gridded precipitation datasets are becoming a convenient substitute for gauge measurements in hydrological modeling; however, these data have not been fully evaluated across a range of conditions. We compared four gridded datasets (Daily Surface Weather and Climatological Summaries [DAYMET], North American Land Data Assimilation System [NLDAS], Global Land Data Assimilation System [GLDAS], and Parameter-elevation Regressions on Independent Slopes Model [PRISM]) as precipitation data sources and evaluated how they affected hydrologic model performance when compared with a gauged dataset, Global Historical Climatology Network-Daily (GHCN-D). Analyses were performed for the Delaware Watershed at Perry Lake in eastern Kansas. Precipitation indices for DAYMET and PRISM precipitation closely matched GHCN-D, whereas NLDAS and GLDAS showed weaker correlations. We also used these precipitation data as input to the Soil and Water Assessment Tool (SWAT) model that confirmed similar trends in streamflow simulation. For stations with complete data, GHCN-D based SWAT-simulated streamflow variability better than gridded precipitation data. During low flow periods we found PRISM performed better, whereas both DAYMET and NLDAS performed better in high flow years. Our results demonstrate that combining gridded precipitation sources with gauge-based measurements can improve hydrologic model performance, especially for extreme events.

A sensitivity analysis of pesticide concentrations in California Central Valley vernal pools.

Vernal pools are ephemeral wetlands that provide critical habitat to many listed species. Pesticide fate in vernal pools is poorly understood because of uncertainties in the amount of pesticide entering these ecosystems and their bioavailability throughout cycles of wet and dry periods. The Pesticide Water Calculator (PWC), a model used for the regulation of pesticides in the US, was used to predict surface water and sediment pore water pesticide concentrations in vernal pool habitats. The PWC model (version 1.59) was implemented with deterministic and probabilistic approaches and parameterized for three agricultural vernal pool watersheds located in the San Joaquin River basin in the Central Valley of California. Exposure concentrations for chlorpyrifos, diazinon and malathion were simulated. The deterministic approach used default values and professional judgment to calculate point values of estimated concentrations. In the probabilistic approach, Monte Carlo (MC) simulations were conducted across the full input parameter space with a sensitivity analysis that quantified the parameter contribution to model prediction uncertainty. Partial correlation coefficients were used as the primary sensitivity metric for analyzing model outputs. Conditioned daily sensitivity analysis indicates curve number (CN) and the universal soil loss equation (USLE) parameters as the most important environmental parameters. Therefore, exposure estimation can be improved efficiently by focusing parameterization efforts on these driving processes, and agricultural pesticide inputs in these critical habitats can be reduced by best management practices focused on runoff and sediment reductions.

A unified approach for protecting listed species and ecosystem services in isolated wetlands using community-level protection goals.

The protection of listed species through the Ecological Risk Assessment (ERA) process is encumbered by the number and diversity of species that need protection and the limited data available to inform assessments. Ecological communities within isolated ecosystems often contain a number of biologically diverse endemic, endangered, and threatened species, as well as providing numerous ecosystem services (ES). We propose an approach that develops community-level protection goals using isolated wetlands that includes both listed species and Service Providing Units (SPUs) that drive ES for ecological risk assessments (ERAs). Community-level protection goals are achieved by developing a protection community and weighing lines of evidence to determine a set of focal species within that community upon which to base the assessment. Lines of evidence include chemical mechanism of action, likely routes of exposure, and taxa susceptibility, as well as relationships among species, and other ecological factors. We demonstrate the process using case studies of chlorpyrifos in California vernal pools and coal ash effluent in Carolina bays. In the California vernal pool case study, listed species were the primary SPUs for the ES provided by the critical habitat. The weight of evidence demonstrated the honey bee as the focal species for the terrestrial environment and the vernal pool fairy shrimp as the focal species for the aquatic environment. The protection community within the Carolina bay case study was more taxonomically diverse than vernal pools for both listed species and SPUs, with amphibians identified as the focal species for which to target mitigation goals and hazard levels. The approach presented here will reduce the time and resource investment required for assessment of risk to listed species and adds an ES perspective to demonstrate value of assessments beyond listed species concerns.

Ecohydrological Index, Native Fish, and Climate Trends and Relationships in the Kansas River Basin.

This study quantified climatological and hydrological trends and relationships to presence and distribution of two native aquatic species in the Kansas River Basin over the past half-century. Trend analyses were applied to indicators of hydrologic alteration (IHAs) at 34 streamgages over a 50-year period (1962-2012). Results showed a significant negative trend in annual streamflow for 10 of 12 western streamgages (up to -7.65 mm/50 yr) and smaller negative trends for most other streamgages. Significant negative trends in western Basin streamflow were more widespread in summer (12 stations) than winter or spring (6 stations). The negative-trend magnitude and significance decreased from west to east for maximum-flow IHAs. Minimum- flow IHAs, however, significantly decreased at High Plains streamgages but significantly increased at Central Great Plains streamgages. Number of zero-flow days showed positive trends in the High Plains. Most streamgages showed negative trends in low- and high-flow pulse frequency and high-flow pulse duration, and positive trends in low-flow pulse duration. These results were consistent with increasing occurrence of drought. Shift in occurrence from present (1860-1950) to absent (2000-2012) was significantly related (p<0.10) to negative trends of 1-day maximum flows (both species) and indices associated with reduced spawning-season flows for Plains Minnow and shifting annual-flow timing and increased flow intermittency for Common Shiner. Both species were absent for all western Basin sites and had different responses to hydrological index trends at eastern Basin sites. These results demonstrate ecohydrological index changes impact distributions of native fish and suggest target factors for assessment or restoration activities.

Pasture BMP effectiveness using an HRU-based subarea approach in SWAT.

Many conservation programs have been established to motivate producers to adopt best management practices (BMP) to minimize pasture runoff and nutrient loads, but a process is needed to assess BMP effectiveness to help target implementation efforts. A study was conducted to develop and demonstrate a method to evaluate water-quality impacts and the effectiveness of two widely used BMPs on a livestock pasture: off-stream watering site and stream fencing. The Soil and Water Assessment Tool (SWAT) model was built for the Pottawatomie Creek Watershed in eastern Kansas, independently calibrated at the watershed outlet for streamflow and at a pasture site for nutrients and sediment runoff, and also employed to simulate pollutant loads in a synthetic pasture. The pasture was divided into several subareas including stream, riparian zone, and two grazing zones. Five scenarios applied to both a synthetic pasture and a whole watershed were simulated to assess various combinations of widely used pasture BMPs: (1) baseline conditions with an open stream access, (2) an off-stream watering site installed in individual subareas in the pasture, and (3) stream or riparian zone fencing with an off-stream watering site. Results indicated that pollutant loads increase with increasing stocking rates whereas off-stream watering site and/or stream fencing reduce time cattle spend in the stream and nutrient loads. These two BMPs lowered organic P and N loads by more than 59% and nitrate loads by 19%, but TSS and sediment-attached P loads remained practically unchanged. An effectiveness index (EI) quantified impacts from the various combinations of off-stream watering sites and fencing in all scenarios. Stream bank contribution to pollutant loads was not accounted in the methodology due to limitations of the SWAT model, but can be incorporated in the approach if an amount of bank soil loss is known for various stocking rates. The proposed methodology provides an adaptable framework for pasture BMP assessment and was utilized to represent a consistent, defensible process to quantify the effectiveness of BMP proposals in a BMP auction in eastern Kansas.