Responses of juvenile mussels to metals in sediment and water of the Tri-State Mining District

<p>The U.S. Geological Survey and collaborators from EcoAnalysts, Inc., completed field and laboratory studies during 2016–19 to evaluate the toxicity of…

Public-domain full text preserved in the Mountain Man Mining Library. Original source: pubs.usgs.gov.

U.S. Department of the Interior U.S. Geological Survey Open-File Report 2023–1024 Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District

Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District By John M. Besser, Chris D. Ivey, James L. Kunz, Nile E. Kemble, Danielle M. Cleveland, Jeffery A. Steevens, Heidi Dunn, and Ryan Foley Open-File Report 2023–1024 U.S. Department of the Interior U.S. Geological Survey

U.S. Geological Survey, Reston, Virginia: 2023 For more information on the USGS—the Federal source for science about the Earth, its natural and living resources, natural hazards, and the environment—visit ://www.usgs.gov or call 1–888–392–8545. For an overview of USGS information products, including maps, imagery, and publications, visit ://store.usgs.gov/ or contact the store at 1–888–275–8747. Any use of trade, firm, or product names is for descriptive purposes only and does not imply endorsement by the U.S. Government. Although this information product, for the most part, is in the public domain, it also may contain copyrighted materials as noted in the text. Permission to reproduce copyrighted items must be secured from the copyright owner. Suggested citation: Besser, J.M., Ivey, C.D., Kunz, J.L., Kemble, N.E., Cleveland, D.M., Steevens, J.A., Dunn, H., and Foley, R., 2023, Responses of juvenile mussels to metals in sediment and water of the Tri-State Mining District: U.S. Geological Survey Open-File Report 2023–1024, 80 p., ://doi.org/​10.3133/​ofr20231024. Associated data for this publication: Ivey, C.D., and Besser, J.M., 2022, Response of juvenile mussels and amphipods to metal concentrations in water and sediment of streams draining the Tri-State Mining District, Missouri, Kansas and Oklahoma, USA: U.S. Geological Survey data release, ://doi.org/​10.5066/​P9R9FHAF. ISSN 2331-1258 (online)

Acknowledgments This project was funded by the U.S. Department of the Interior’s Office of Restoration and Damage Assessment and guided by a team that included Dave Mosby, Suzanne Dunn, Sherry Kircher, and Amy Horner Hanley. Mussel surveys were planned and completed by person­ nel of Ecological Specialists, Inc. (now EcoAnalysts, Inc.), led by Heidi Dunn and Ryan Foley, with assistance from EcoAnalysts, Inc., staff Nathan Badgett, David Ford, Kendall Cranney, and Dan Scoggin and with assistance from U.S. Fish and Wildlife Service personnel David Martinez, Bryan Simmons, and Josh Hundley. It was a pleasure and a great learning experience to work with the malacology team. Larry Tippet and Kyle Foster of the Peoria Tribe of Indians of Oklahoma generously provided a boat, a boat operator, and guidance in the field, allowing us easy access to remote sampling sites. Personnel of the U.S. Geological Survey Columbia Environmental Research Center assisted with fieldwork (Eric Brunson, Rebecca Dorman, Doug Hardesty, and Steve Olson), toxicity testing and water quality analysis (Nile Kemble, Ryan Gettler, and Ning Wang), and chemical analyses (Vanessa Melton, Mike Walther, and Alison Hartman).

Contents Mussel Tests 36

Figures

2.  Boxplot showing zinc concentrations in reconnaissance samples of surface

3.  Graphs showing concentrations and ratios of zinc, lead, and cadmium in pore

5.  Graphs showing characteristics of mussel assemblages at reference sites and

6.  Graphs showing associations of zinc, lead, and cadmium bioaccumulation by Corbicula with metal concentrations in sediment and with species richness of

7.  Graphs showing comparison of mussel survey metrics between low- and

9.  Graphs showing effect of stocking order on mussel endpoints in 2017

10.  Graphs showing mussel endpoints in reference sediments and test sediments

11.  Graphs showing concentration-response plots and 20-percent effect concentrations for effects of sediment zinc on mussel endpoints in 12-week

12.  Graphs showing concentration-response plots and 20-percent effect concentration estimates of EC20s for effects of sediment zinc on amphipod,

13.  Graphs showing comparison of mussel 20-percent effect concentrations

14.  Graphs showing filterable zinc in overlying water of supplemental toxicity

15.  Graphs showing survival and growth of two age/size classes of juvenile mussels in 28-day exposures to 4 composite sediments and 6 metal treatments

16.  Graphs showing relations between mussel endpoints from 2017 sediment toxicity tests and community metrics from quantitative mussel surveys, with

17.  Graph showing estimated 50-percent effect concentration for effects of Tables

2.  Primary sampling sites for mussel community surveys and sediment

3.  Concentrations of metals, major ions, and dissolved organic carbon in samples

4.  Physical and chemical characteristics of sediments from primary sampling

6.  Toxicity of zinc and a three-metal mixture in 12-week water-only tests with

7.  Physical and chemical characteristics of sediment from the Tri-State Mining

8.  Toxicity of sediments to juvenile mussels in 12-week toxicity tests for sediment

9.  Toxicity of Tri-State Mining District sediments in 6-week exposures with

10.  Comparison of toxicity thresholds for Hyalella azteca and Lampsilis siliquoidea

11.  Physical and chemical characteristics of sediment and pore water in

12.  Physical and chemical characteristics of sediment and pore water in

13.  Summary of mussel toxicity trends in 2019 test A: sediment-plus-water toxicity

14.  Summary of mussel toxicity trends in 2019 test B: sediment-plus-water toxicity

Conversion Factors U.S. customary units to International System of Units Multiply By To obtain Flow rate foot per second (ft/s) meter per second (m/s) International System of Units to U.S. customary units Multiply By To obtain Length micrometer (μm) 3.93701×10−5 inch (in.) centimeter (cm) inch (in.) millimeter (mm) inch (in.) meter (m) foot (ft) meter (m) yard (yd) Area square meter (m2) 0.0002471 acre Volume nanoliter (nL) 3.3814×10−8 ounce, fluid (fl. oz) milliliter (mL) ounce, fluid (fl. oz) liter (L) ounce, fluid (fl. oz) liter (L) pint (pt) liter (L) quart (qt) liter (L) gallon (gal) liter (L) cubic inch (in3) Mass milligram (mg) 3.5274×10−5 ounce, avoirdupois (oz) gram (g) ounce, avoirdupois (oz) Temperature in degrees Celsius (°C) may be converted to degrees Fahrenheit (°F) as follows:

°F (1.8 × °C) + 32. Concentrations of chemical constituents in water are given in either milligrams per liter (mg/L), micrograms per liter (µg/L), or grams per liter (g/L). Molar concentrations are given in micromoles per gram (µmol/g). Simultaneously extracted metals concentrations are given in micrograms per gram (μg/g).

Abbreviations ANOVA analysis of variance AVS acid-volatile sulfide CCVS continuing calibration standards CERC Columbia Environmental Research Center CPUE catch per unit effort CRM concentration-response model DOC dissolved organic carbon EC20 20-percent effect concentration EPA U.S. Environmental Protection Agency ESA U.S. Endangered Species Act GM–EC20 geometric mean of 20-percent effect concentrations ICP–MS inductively coupled plasma-mass spectrometer ICVS initial calibration standards LCS laboratory control standard MT metallothionein NRDAR Natural Resource Damage Assessment and Restoration p probability PEQ probable effect quotient PES polyethersulfone PWP-zinc zinc in pore-water peepers QA/QC quality assurance/quality control SE standard error SEM simultaneously extracted metals [SEM–AVS] simultaneously extracted metals normalized to acid-volatile sulfide [ΣSEM–AVS] summed molar concentrations of simultaneously extracted metals relative to acid-volatile sulfide sum-PEQ summed probable effect quotient values

T20 20-percent effect concentration for survival TAC test acceptability criterion TOC total organic carbon TR total recoverable USGS U.S. Geological Survey XRF X-ray fluorescence

greater than ≥ greater than or equal to

less than

Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District By John M. Besser,1 Chris D. Ivey,1 James L. Kunz,1 Nile E. Kemble,1 Danielle M. Cleveland,1 Jeffery A. Steevens,1 Heidi Dunn,2 and Ryan Foley2 Abstract The U.S. Geological Survey and collaborators from EcoAnalysts, Inc., completed field and laboratory studies dur­ ing 2016–19 to evaluate the toxicity of metals to freshwater mussels in streams draining the Tri-State Mining District. This project consisted of (1) sampling and analysis of metals in water and sediment, (2) surveys of mussel assemblages at sites with suitable mussel habitat, (3) toxicity tests with juvenile mussels exposed to zinc or to a mixture of metals (zinc, lead, and cadmium) in water, and (4) toxicity tests to evaluate the contributions of metals in sediment and metals in overly­ ing water to toxic effects on mussels. Field sampling at sites in the Spring River and Neosho River and their tributaries demonstrated wide ranges of metal contamination in water and sediment. Zinc was the predominant toxic metal in water, and concentrations of lead and cadmium were much lower. Mussel areal density and species richness were greater at reference sites with low sediment metal concentrations (for example, zinc, 29–141 micrograms per gram) than at test sites that had higher concentrations of sediment zinc (416–3,420 micro­ grams per gram) as a result of effects of upstream mining activity. Juvenile mussels were highly sensitive to zinc in water in 12-week toxicity tests compared to previous wateronly tests, and adding low levels of waterborne lead and cadmium typical of their occurrence in Tri-State Mining District streams produced greater toxicity. Thresholds for mus­ sel toxicity were at or less than waterborne metal concentra­ tions detected in Tri-State Mining District streams, and sites with waterborne metal concentrations exceeding thresholds had decreased mussel density and decreased mussel species richness. The 12-week toxicity tests with juvenile mussels in Tri-State Mining District sediments also demonstrated nega­ tive mussel responses with metal exposure. Thresholds for reductions in survival, growth, or biomass were at sediment metal concentrations less than thresholds reported for previ­ ous 4-week tests. We documented strong associations between reduced survival in laboratory tests and reduced species rich­ ness in community surveys. Attempts to estimate combined 1U.S. Geological Survey. 2EcoAnalysts, Inc. toxicity thresholds for metals in sediment and overlying water were not successful. These inconclusive results may be attrib­ utable to several factors, including (1) unexpected losses of waterborne metals from solution, (2) differences in sensitivity of different age/size classes of juvenile mussels, (3) disruption of sediment-water equilibria and changes in metal bioavail­ ability, and (4) behavioral or physiological responses allowing juvenile mussels to temporarily reduce or avoid metal expo­ sure. We also observed differences in metal toxicity thresholds between sediment toxicity tests started with different ages/ sizes of test organisms. A followup study that combined expo­ sure to Tri-State Mining District sediments with exposures to multiple levels of waterborne metals demonstrated toxic effects of sediments with low metal concentrations; however, some treatments also indicated unexpected reversals of con­ centration-response trends and reduced toxicity in treatments that had high metal concentrations in overlying water. These unusual responses may reflect development of physiological tolerance to metal toxicity by induction of metal-binding pro­ teins (for example, metallothionein) in response to high metal levels in water. Results of laboratory and field studies indicated strong associations between metal exposure in Tri-State Mining District streams and toxic effects on juvenile freshwater mussels. Mussel community characteristics corresponded to differences in metal concentrations in sediment and water among Tri-State Mining District sampling sites. Responses of juvenile mussels in 12-week water and sediment exposures were strongly correlated with the status of mussel assemblages in Tri-State Mining District streams. The combined results support the hypothesis that exposure to metals from historical mining activities adversely affects freshwater mussel commu­ nities in the Spring River/Neosho River drainage. 1. Introduction and Scope Freshwater mussels (Unionidae) are sensitive to toxic effects of many aquatic contaminants, including metals (Wang and others, 2010, 2013, 2020; Besser and others, 2015). Mussels also are frequently listed as endangered or threatened under the U.S. Endangered Species Act (ESA; Williams and

2    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District others, 1993; Haag and Williams, 2014). Effects of pollu­ tion on ESA-listed mussels can be a legal basis for recovery of natural resource damages by the U.S. Department of the Interior and State and Tribal natural resource trustees through the Natural Resource Damage Assessment and Restoration (NRDAR) program. Two historical mining areas in Missouri, the Southeast Missouri Lead District (Roberts and others, 2016) and the Tri-State Mining District are currently subjects of NRDAR investigations of potential injury to ESA-listed freshwater mussels. Studies of mussel communities in streams draining these areas have documented declines in resident mussel communities, which correspond to metal pollution from historical mining (Angelo and others, 2007; Ingersoll and others, 2008; Besser and others, 2015). In the Tri-State Mining District, Angelo and others (2007) reported reduced mussel taxa richness at sites downstream from mining areas, compared with upstream sites, and determined statistically significant negative correlations between mussel taxa richness and metal concentrations in sediments. Laboratory sedi­ ment toxicity tests completed by the U.S. Geological Survey (USGS) in 2007 with metal-contaminated sediments from streams documented adverse effects on several species of benthic invertebrates (Ingersoll and others, 2008), but toxic effects on juvenile mussels (Lampsilis siliquoidea [Barnes, 1823]; fatmucket) were less severe than effects on other test organisms, notably the freshwater amphipod (Hyalella azteca [Saussure, 1858]; Besser and others, 2015). Lampsilis siliquoi­ dea has been widely used for laboratory toxicity testing and is a reliable surrogate for responses of other unionid mussel spe­ cies, including ESA-listed species (Wang and others, 2007). The purpose of this report is to document metal toxicity thresholds associated with injury to mussels based on longterm exposure to metals in water and in sediment. To accom­ plish this, we designed and completed studies that included four components: 1. Sampling and analysis of metals in water and sediment at selected sites in the Tri-State Mining District, 2. Surveys of mussel assemblages at sites sampled for water and sediment, 3. 12-week toxicity tests to evaluate thresholds for toxic effects of metals in Tri-State Mining District stream water to juvenile mussels and amphipods, and 4. 12-week toxicity tests to evaluate thresholds for toxic effects of metal-contaminated sediment to juvenile mus­ sels, with and without added metals in overlying waters. Sites for collection of water and sediment for chemi­ cal analysis and toxicity testing were colocated with sites of mussel community surveys whenever possible to enable more direct comparisons of mussel responses in the laboratory and in the field. Exposure of mussels to metals in laboratory studies and at field survey sites was characterized by chemi­ cal analyses of metals in sediment, pore water, and overlying water. Low-metals reference sites and test sites potentially affected by metals from mining were characterized to ensure they represented suitable physical habitat for freshwater mussels (EcoAnalysts, Inc., 2018). Toxicity tests with juve­ nile mussels (L. siliquoidea) included exposure to metals in field-collected sediments and laboratory-prepared waterborne mixtures typical of streams in the study area, separately and in combination. Mussel toxicity tests were started with smaller juvenile mussels and had a longer (12-week) exposure period, compared to previous tests with sediments (Ingersoll and oth­ ers, 2008). For comparison with results of the previous 4-week amphipod tests, 6-week sediment toxicity tests also were completed with amphipods (H. azteca; MacDonald and others, 2009; Besser and others, 2015). Data generated during this study are available as a USGS data release (Ivey and Besser, 2022). 2. Metal Concentrations in Water and Sediment and Status of Mussel Community Mining in the Tri-State Mining District took place over a large geographic area (fig. 1) in adjacent parts of Missouri, Kansas, and Oklahoma (Johnson and others, 2016). In south­ western Missouri, large-scale mining took place primarily in lands drained by west-flowing tributaries of the Spring River, including Center Creek, Turkey Creek, and Shoal Creek. Mining-affected areas of southeastern Kansas included small tributaries and the main stem of the Spring River, including an impoundment, Empire Lake. In northeastern Oklahoma, mining-affected streams included highly contaminated Tar Creek and downstream reaches of the Neosho River, as well as the lower Spring River downstream to its confluence with the Neosho River and its terminus in Grand Lake O’ the Cherokees. Subsurface mining throughout the Tri-State Mining District until the 1970s left many large deposits of fine-grained tailings near ore-processing locations. Tailings were distrib­ uted across the landscape by wind and water erosion and by human activities, resulting in metal contamination of stream sediments and leaching of metals into groundwater and surface water. Unlike some other mining areas where Mississippi Valley type ores were mined primarily for lead, ores and mine wastes in the Tri-State Mining District were predominantly enriched with zinc rather than lead, and the greater solubility of zinc led to potentially greater toxicity hazards via exposure to waterborne metals (Besser and others, 2015; Gutiérrez and others, 2020). To date, remediation of contaminated areas in the Tri-State Mining District has primarily been completed under the auspices of the Comprehensive Environmental Response, Compensation, and Liability Act. Remediation efforts have focused on removal of metalcontaminated residential soils, consolidation of mine wastes and capping with uncontaminated soil in upland habitats,

2. Metal Concentrations in Water and Sediment and Status of Mussel Community    3 reburial and capping of mine wastes in subsurface mine workings, and pH neutralization and addition of organic amendments to contaminated subsoils to encourage revegetation (Juracek and Drake, 2016). Additional remedial efforts involving removal of mine waste from ephemeral streams were completed in Missouri and Kansas; however, remedial efforts involving removal or isolation of solid mine wastes from stream channels, floodplains, and impoundments are in the planning stages. Despite these ongoing efforts to control sources of metal contamination, many miles of streams remain contaminated by metal-contaminated sediments and or by dissolved metals leached from mine wastes and contaminated soils. Data from the 2007 USGS investigation of sediment toxicity in the Tri-State Mining District (Ingersoll and others, 2008) were analyzed by MacDonald and others (2009) with the goal of establishing objectives for cleanup based on assessment of residential human health risk and risk of toxicity to aquatic biota. Targets for remediation of stream and reservoir sediments were based mainly on toxicity tests with H. azteca, which was the most sensitive species tested. Despite the evidence of effects of metals on mussel communities of the Spring River drainage reported by Angelo and others (2007), the previous sediment toxicity studies determined that juvenile L. siliquoidea were less sensitive to toxicity of metals than were H. azteca (MacDonald and others, 2009). This study was part of the ongoing NRDAR cases by Trustees Councils in Kansas, Missouri, and Oklahoma and was focused on characterizing toxic effects of metals on freshwater mussels rather than effects of other pollutants or habitat degradation. Therefore, sampling sites for mussel community surveys were limited to sites that currently supported a mussel assemblage or sites judged to be suitable mussel habitat based on site characteristics (EcoAnalysts, Inc., 2018). This restriction meant that sampling would not include some highly contaminated stream reaches, because the habitat of those areas was highly degraded, making it difficult to distinguish between effects of habitat degradation and toxic effects of metals. This study used the “bulk” fraction of sediment particles less than 2 millimeters (<2 mm) in diameter for sediment toxicity testing. This strategy differed from the 2007 USGS study, which sieved sediments so that only the “fine” fraction EXPLANATION Site Primary mussel Primary sediment Primary (all matrices) Reconnaissance N 95° 94°50' 37°10' 36°50' 37° 94°40' 94°30' 94°20' 94°10' Kansas Missouri Oklahoma Arkansas Base from Esri © 1999–2020 and its licensors Figure 1.  Locations of sampling sites. Map modified from EcoAnalysts, Inc. (2018). Inset shows the approximate location of the Tri- State Mining District (hatched area). Inset data from U.S. Environmental Protection Agency (2021a).

4    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District (<0.25-mm diameter) would be used for toxicity testing, to ensure that live small juvenile mussels could be recovered quantitatively from fine sediments after 28 days. Sediment metal concentrations differed between the <0.25-mm fine fraction and the <2.0-mm fraction, however, making it dif­ ficult to compare results of mussel sediment tests to results of tests with other test organisms. The decision to test mus­ sels with <2-mm sediments in the current study also had the advantage of making it easier to collect the large volumes of sediment needed for toxicity testing from sites chosen for mussel surveys. Objectives of Reconnaissance Sampling The objectives of the reconnaissance sampling were as follows: 1. To identify about 25 primary sampling sites widely distributed in the study area with habitat suitable for freshwater mussels (Unionidae); 2. To collect samples of surface water, pore water, sedi­ ments, and tissues of the nonunionid freshwater clam, Corbicula fluminea (O.F. Müller, 1774; Corbicula), at primary sampling sites to characterize metal concentra­ tions and bioavailability and to provide samples for sediment toxicity testing; and 3. To complete qualitative and quantitative surveys of abundance and species richness of unionid mussels at primary sampling sites. Methods Selection of Reconnaissance Sites The selection of primary sampling sites for sediment toxicity testing and mussel community surveys started with a list of about 49 candidate sites, based on published sediment chemistry and toxicity data (for example, Ingersoll and others, 2008; MacDonald and others, 2009; Besser and others, 2015) and mussel community data (Angelo and others, 2007). Sites were also selected using a geographic information system project that delineated historically stable stream reaches that were most likely to support populations of long-lived freshwa­ ter mussels (EcoAnalysts, Inc., 2018; table 1). These original reconnaissance sites were identified by four-digit IDs, consist­ ing of two letters to identify stream plus two numbers to indi­ cate upstream/downstream sequence. During reconnaissance sampling, additional alternate sites were selected near the original sites in order to find suitable habitat; these supplemen­ tal sites are indicated by site IDs with letters added to the end of the original site numbers (table 1). Sites for reconnaissance sampling were characterized for stream order, gradient, and channel stability using geo­ graphic information system data, and the suitability of each candidate site as mussel habitat was assessed using a mussel habitat checklist (appendix 1). This checklist included factors such as distance to bluffs or bedrock exposures; geomorphic stream features (for example, riffles, pools, bars, terraces, and [or] islands); and features that indicate stream stability (for example, evidence of channel migration compared to histori­ cal maps, incised channels, and cut banks), distance from the mouth of the stream, degree of sediment embeddedness, gen­ eral observations of sediment grain size, presence of woody debris, condition of riparian corridor, presence or absence of unionid mussels and Corbicula, and distance to sediment sampling locations. The availability of sediments in the sand-sized and smaller (<2 mm) size fraction in wadeable habitats (<1-meter depth) was assessed at each reconnaissance site. At each of these sites, about 0.5 liter (L) of sediment was collected with a polyvinyl chloride scoop sampler from several potential sampling areas, composited, homogenized, and stored in a precleaned high-density polyethylene container. Sediment samples were dried in an oven at about 60 degrees Celsius (°C). After the total dry weight of the composite sample was determined, the sample was crushed and passed through a 2-mm stainless-steel sieve, and the fraction of the sample in the <2-mm fraction (on a dry mass basis) was determined. A part of each sediment sample was placed in a plastic bag for screening analysis of zinc and lead by X-ray fluorescence (XRF; U.S. Environmental Protection Agency [EPA], 2007). Cadmium concentrations in sediments were typically less than XRF detection limits. Samples were analyzed using a Thermo Niton XL2 GOLDD XRF analyzer (Thermo Scientific, Waltham, Massachusetts) for 1 minute by placing the instru­ ment window directly against a part of the bag that was in full contact with the sediment. Three readings were taken for each sample to ensure a representative sample. Concentrations of metals were expressed as probable effect quotients (PEQs) by dividing the mean concentration of each metal by its probable effect concentration (MacDonald and others, 2000). PEQs for the three metals were assumed to be additive, such that a sedi­ ment with a single metal producing a PEQ of 1.0 (that is, one metal at a concentration equal to its probable effect concentra­ tion) poses a toxicity hazard equal to a sediment with summed PEQ values (sum-PEQ) of 1.0 for multiple metals. These preliminary hazard quotients were used to help select primary sampling sites for sediment toxicity testing. Samples of stream water and pore water also were collected from reconnais­ sance sites for analysis of water quality properties and metals concentrations. Selection of Primary Sampling Sites The 25 primary sampling sites each included suit­ able mussel habitat and had enough fine sediments to allow assessment of sediment toxicity. Within these constraints, site selections sought to include broad geographic coverage of the Spring River and Neosho River drainages basin and broad ranges of metal contamination. Five sites were selected to

2. Metal Concentrations in Water and Sediment and Status of Mussel Community    5 Table 1.  Sites inspected and sampled during reconnaissance, 2016–17. [Dates are given in month/day/year format; NF, North Fork; MO, Missouri; X, site sampled for water, sediment, and mussel surveys; --, no data or not applicable; KS, Kansas; OK, Oklahoma] Stream Site State Water Sediment Mussel Date Latitude (in degrees north) Longitude (in degrees west) N. Fork Spring NF01 MO 7/19/2016 N. Fork Spring NF02 MO 7/19/2016 Spring SP01A MO 7/19/2016 Spring SP01B MO 7/19/2016 Spring SP02 MO 7/20/2016 Spring SP03 KS 10/19/2016 Spring SP05 KS 10/18/2016 Spring SP06 KS 10/18/2016 Spring SP07 KS 10/17/2016 Spring SP08 KS 10/19/2016 Spring SP09 KS 10/5/2016 Spring SP10A KS 10/14/2016 Spring SP10B KS 10/14/2016 Spring SP10C KS 10/14/2016 Spring SP10D KS 10/14/2016 Spring SP10E KS 10/14/2016 Spring SP11A KS 10/3/2016 Spring SP11B KS 10/3/2016 Spring SP11C KS 10/3/2016 Spring SP12A KS 10/5/2016 Spring SP12B KS 10/5/2016 Spring SP13 OK 8/23/2016 Spring SP13A OK 8/23/2016 Spring SP14 OK 8/23/2016 Spring SP15 OK 8/23/2016 Spring SP16 OK 8/24/2016 Spring SP17 OK 8/24/2016 Spring SP18 OK 8/24/2016 Spring SP19 OK 8/24/2016 Spring SP20 OK 8/24/2016 Spring SP21 OK 8/25/2016

6    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District Table 1.  Sites inspected and sampled during reconnaissance, 2016–17.—Continued [Dates are given in month/day/year format; NF, North Fork; MO, Missouri; X, site sampled for water, sediment, and mussel surveys; --, no data or not applicable; KS, Kansas; OK, Oklahoma] Stream Site State Water Sediment Mussel Date Latitude (in degrees north) Longitude (in degrees west) Spring SP22 OK 8/25/2016 Spring SP23 OK 8/25/2016 Spring SP24 OK 8/25/2016 Center MO 10/6/2016 Center MO 7/20/2016 Center CC01G1 MO 7/26/2016 Center CC01G2 MO 7/26/2016 Center CC01A2 MO 7/24/2016 Center CC01A1 MO 7/24/2016 Center CC01B1 MO 7/24/2016 Center CC01B2 MO 7/24/2016 Center MO 7/24/2016 Center CC01D MO 7/24/2016 Center CC01E MO 7/24/2016 Center CC01F MO 7/24/2016 Center MO 7/21/2016 Center MO 7/20/2016 Center CC05D MO 7/20/2016 Center CC05DA MO 7/20/2016 Center CC05DB MO 7/20/2016 Center CC05DC MO 7/21/2016 Center CC05D MO 7/21/2016 Center MO 7/21/2016 Center MO 7/21/2016 Center CCO8 MO 7/21/2016 Center MO 7/21/2016 Center MO 7/21/2016 Center CC15A MO 7/26/2016 Center MO 7/26/2016 Shoal SH01 MO 7/20/2016 Shoal SH01A MO 7/22/2016 Shoal SH02 MO 7/23/2016

2. Metal Concentrations in Water and Sediment and Status of Mussel Community    7 Table 1.  Sites inspected and sampled during reconnaissance, 2016–17.—Continued [Dates are given in month/day/year format; NF, North Fork; MO, Missouri; X, site sampled for water, sediment, and mussel surveys; --, no data or not applicable; KS, Kansas; OK, Oklahoma] Stream Site State Water Sediment Mussel Date Latitude (in degrees north) Longitude (in degrees west) Shoal SH02A MO 7/23/2016 Shoal SH03A MO 7/23/2016 Shoal SH04A MO 7/23/2016 Shoal SH04B MO 7/23/2016 Shoal SC04C MO 7/23/2016 Shoal SH05 MO 7/23/2016 Shoal SH05A MO 7/23/2016 Shoal SH06A MO 7/22/2016 Shoal SH06 MO 7/22/2016 Shoal SH07 MO 7/22/2016 Shoal SH09A MO 7/25/2016 Shoal SH10 MO 7/25/2016 Shoal SH11 MO 7/25/2016 Shoal SH12 MO 7/25/2016 Shoal SH13 MO 7/25/2016 Shoal SH14 MO 7/25/2016 Shoal SH15 MO 7/25/2016 Shoal SH15A MO 7/25/2016 Shoal SH16 KS 10/6/2016 Shoal SH16A KS 10/6/2016 Shoal SH17 KS 10/6/2016 Shoal SH18 KS 10/6/2016 Neosho NRSF OK [missing data] [missing data] [missing data] Neosho NR02 OK 8/26/2016 Neosho NR03 OK 8/26/2016 Neosho NR04 OK 8/26/2016 Neosho NR05 OK 8/26/2016 Neosho NR06 OK 8/27/2016 Tar TC01 OK 8/25/2016 Tar TC02 OK 8/27/2016 Lost OK 8/25/2016

8    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District represent reference conditions, which were defined as sedi­ ments with sum-PEQs of <1.0 for zinc, lead, and cadmium and with habitat characteristics comparable to sites that were potentially affected by metals from mining (test sites). One additional low-metal sample was collected from the Spring River at Waco, Missouri (site SR), to serve as a control sedi­ ment to document performance of the sediment bioassay. Good performance of juvenile mussels in the site SR sediment has been documented in previous studies (Ingersoll and others, 2008; Besser and others, 2015). Mussel Surveys Mussel surveys were completed at reconnaissance sites (qualitative surveys) and at primary sampling sites (quantita­ tive surveys). During the reconnaissance sampling, the status of the mussel community at sites with suitable habitat was assessed by timed searches. The area of the site was delin­ eated with a survey-grade Trimble Global Positioning System, and the site was qualitatively sampled with timed searches to determine unionid distribution, community composition, and catch per unit effort (CPUE). Search effort per site ranged from 30 to 360 person-minutes and averaged 76 minutes. Primary sampling sites were sampled quantitatively for mussels in 2017. At each site, as many as 100 randomly dis­ tributed 0.25-square meter quantitative quadrats were sampled within the delineated unionid habitat area to characterize species richness, density, and age distribution. The number of quadrat samples per site ranged from 15 to 100, depending on the size of the delineated area. Substrate in each quadrat was excavated to a depth of about 15 centimeters (cm), and material was placed into an attached mesh bag. Each sample was washed through a 6-mm sieve and searched for unionids and C. fluminea. All live unionids were identified to species, measured (length in millimeters), and aged (external annu­ lus count). Fresh dead shells were identified to species and counted. Weathered dead and subfossil shells were identified to species and noted as present. Additional details of methods for the mussel community and habitat survey are presented in appendix 1 and in EcoAnalysts, Inc. (2018). Water and Sediment Sampling Metal concentrations and water quality properties affecting metal bioavailability were determined in samples of stream water and in situ pore water collected at reconnais­ sance sites and primary sampling sites. Samples of stream water were collected by subsurface grabs with precleaned polypropylene sample vials. The vials were opened, ultrapure water in the vials was discarded, and the vials were rinsed with and then filled with stream water and sealed below the water surface. Samples of surficial pore water were collected by push-point sampling (Zimmerman and others, 2005) in or near mussel beds. The fritted tip of the sampler was inserted into the surface sediment to a depth of 10 cm, pore water was extracted from the sampler with a syringe, and the first aliquot from each location was used to rinse the sampler tub­ ing and the sample vial and then discarded. The first set of reconnaissance water samples was collected in early summer 2016 (about 40 sites), and a second set of water samples was collected during low flow conditions in late summer and early fall 2016 (25 sites). Metal concentrations from these samples were used to establish metal concentration ranges and ratios of metal mixtures for water-only toxicity tests completed in winter 2017. Sediments for toxicity testing were collected at 25 pri­ mary sampling sites in summer 2017 following methods used for the previous sediment study (Ingersoll and others, 2008; Besser and others, 2015). Sediments were extracted from depositional areas with a 7.5- or 10-cm diameter polyvinyl chloride scoop and were wet sieved through a wash bucket (Wildco Wash Bucket with US no. 10 stainless-steel mesh) to separate sand-sized and smaller particles (<2 mm) for toxicity testing. The use of the <2-mm sediment fraction differed from the <0.25-mm fraction used for testing and analysis in 2007 (Ingersoll and others, 2008). Sediments were sealed in pre­ cleaned high-density polyethylene buckets with site water and substantial head space, transported to the USGS laboratory in a refrigerated trailer, and stored in a walk-in cooler at 4 °C. Sediment was stored 2 months or less before the first sediment test in 2017, then an additional 5 months before the second sediment test in spring 2018. About 1 week before the start of sediment toxicity tests, the sieved (<2 mm) sediments from 24 sites plus the control sediment were homogenized with an electric drill equipped with a stainless-steel auger, and subsamples were collected for analysis of total recoverable (TR) metals, simultaneously extracted metals (SEM), acid-volatile sulfide (AVS), total organic carbon (TOC), and particle-size distribution. Porewater samples were separated from sediments by centrifuga­ tion and by passive diffusion samplers (“peepers”; Brumbaugh and others, 2007; Cleveland and others, 2017). Sediments were added to 250-milliliter (mL) polyethylene conical cen­ trifuge bottles and centrifuged for 30 minutes at 5,000 revo­ lutions per minute or 7,000 times standard gravity (Sorvall LYNX 400 with Bio-flex HC rotor) to generate pore water for analyses of cations, anions, dissolved organic carbon (DOC), and conventional water quality properties (pH, conductivity, alkalinity, hardness, total ammonia). Pore water also was sam­ pled before toxicity tests using large-volume peepers (30 mL; Cleveland and others, 2017) deployed in 1,000-mL beakers containing 100 mL of homogenized sediment and equilibrated for 10 days at 4 °C. Pore water from the large peepers was analyzed for metals and DOC for comparison to pore water obtained from centrifuged samples. Water Analyses Concentrations of filterable metals (zinc, lead, cad­ mium) and major cations (calcium, potassium, magnesium, sodium) were determined in filtered (0.45-micrometer [µm]

2. Metal Concentrations in Water and Sediment and Status of Mussel Community    9 pore diameter) water from grab samples, peeper samples, and push-point samples. Samples were analyzed for metals using an inductively coupled plasma-mass spectrometer (ICP–MS) (ELAN DRC–e; Perkin-Elmer, Waltham, Mass.). Stream water and overlying waters from toxicity tests were filtered with a 0.45-µm polyethersulfone (PES) filter disc and preserved by acidification with subboiling distilled nitric acid (1 percent nitric acid by volume). Concentrations of anions in stream water, in overlying waters of toxicity tests, and in sediment pore waters were filtered (0.45-μm polyproplylene/polyethersulfone) and ana­ lyzed with an ion chromatograph (ICS–1100, Ionpac AS22 anion exchange column, and suppressed conductivity detector; Dionex, Arcade, New York). Anion samples were stored in a laboratory refrigerator at 4 °C for no more than 30 days before analysis. Concentrations of DOC in stream water, overlying waters of sediment toxicity tests, and sediment pore waters (grab, peeper, and [or] push-point samples) were filtered (0.45-μm glass/PES), dispensed into precleaned low-TOC amber glass vials, and acidified with 9 N high-purity sulfuric acid (BDH Aristar Ultra; VWR International, Radnor, Pennsylvania) to a pH of 2 or less within 48 hours of receipt. Acidified samples were refrigerated and held no longer than 28 days before DOC analysis. DOC was measured as nonpurgeable organic carbon by high temperature combustion catalytic oxidationnondispersive infrared spectroscopy using a TOC analyzer (model TOC–L CSH; Shimadzu Scientific Instruments, Inc., Columbia, Maryland). Samples were air sparged to remove inorganic carbon species before the nonpurgeable organic carbon measurements. Sediment and Tissue Analyses Metal analyses were completed on samples of sediment (<2 mm) and Corbicula tissue. Sediment was analyzed for TR metals, SEM and AVS, TOC, and particle-size distribution. Sediments were microwave digested in a mixture of nitric and hydrochloric acids to extract the TR metal fraction, and tissues were digested with nitric acid and hydrogen peroxide. Metal concentrations in these extracts were determined by ICP–MS. Sediment samples for measurement of SEM and AVS were refrigerated and analyzed within 28 days of receipt using deoxygenated reagents in a two-part process. First, the sedi­ ment was mixed with 1 N hydrochloric acid (20:1 mixture of sediment to acid on a wet weight basis) in a nitrogen-purged oxygen-free atmosphere to release sulfide as hydrogen sulfide gas. The hydrogen sulfide was bubbled through an antioxidant trapping solution with a pH of 12 to form free sulfide, which was measured using an ion-selective electrode (Brumbaugh and others, 2011). After sulfide trapping, overlying liquid and the upper 1 cm of sediment were removed and discarded, and the remaining sediment was gently stirred for as much as 1 minute to produce a uniform consistency before sampling the SEM extract. The sediment-hydrochloric acid mixture was settled for a maximum of 15 minutes, and then 40 mL of the extract were filtered through a 0.45-μm PES membrane, diluted in nitric acid, and (or) further digested and (or) volume reduced as needed for ICP–MS analysis. The SEM fraction was analyzed for cadmium, copper, lead, nickel, and zinc. Metal concentrations in Corbicula tissues were used as a surrogate for metal bioaccumulation by unionid mussels. Corbicula collected during quantitative mussel surveys were held in site water at ambient temperatures for 24 hours to allow them to depurate gut contents before they were frozen for storage (appendix 2). In the laboratory, soft tissues were separated, lyophilized, homogenized, and microwave digested in nitric acid and hydrogen peroxide before ICP–MS analyses. TOC content and particle-size distribution of sediments were analyzed at the University of Missouri Soils Laboratory using methods 4H2 and 3A1a from the U.S. Department of Agriculture Kellogg Soil Survey Laboratory (Natural Resources Conservation Service, 2014). The sieved sediment fraction (<2 mm) was characterized for particle-size distribu­ tion using a dispersion and suspension method. The sample was pretreated with peroxide and other reagents to remove organic matter and soluble salts, then ovendried to obtain an initial weight. The particles were then dispersed with a sodium hexametaphosphate solution and mechanically shaken to create a suspension. The sand fraction was removed from the suspension by wet sieving and then fractionated by dry sieving. The clay and fine silt fractions were determined using the suspension remaining from the wet sieving process. This suspension was diluted to 1 L in a sedimentation cylinder and stirred, and then 25-mL aliquots were removed with a pipette at predetermined intervals. The aliquots were dried at 110 °C and weighed. For total carbon analysis, the sieved sediment was packed in foil, weighed, and analyzed for total carbon by an elemental analyzer. The elemental analyzer used catalytic combustion in an oxygenated atmosphere and high tempera­ ture to liberate carbon as carbon dioxide, which was detected by a thermal conductivity detector. Inorganic carbon was determined by treatment of the soil with hydrochloric acid and subsequent manometric measurement of the evolved carbon dioxide, and TOC was estimated by difference of the total and inorganic carbon fractions. Quality Assurance A minimum of three external calibration standards plus a calibration blank were used to calibrate instrument response, and established laboratory quality assurance/quality con­ trol (QA/QC) procedures (for example, laboratory spikes, duplicates, control samples) were used to verify instrument performance throughout the analyses. A summary of QA/QC methods and results is included in appendix 2.

10    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District Results and Discussion Selection of Primary Sampling Sites The USGS and EcoAnalysts, Inc., visited 97 locations in the Spring River-Neosho River drainages between August and October 2016 and completed reconnaissance sampling of water and sediment plus qualitative surveys for freshwater mussels at 42 sites (table 1). Screening analyses of metals in sediment samples by XRF were used to select 25 primary sam­ pling sites for sediment toxicity testing (table 2). These pri­ mary sampling sites included one site for collection of control sediment from the Spring River at Waco, Mo. (Ingersoll and others, 2008; Besser and others, 2015); 5 reference sites (3 in upper Spring River, 2 in upper Center Creek); and 19 test sites (8 in Spring River, 4 in Center Creek, 6 in Shoal Creek, and 1 each in Tar Creek and Neosho River; table 2). Most primary sampling sites were included in both sediment toxicity testing and quantitative mussel community surveys, but several sites in Center Creek (CC09) and Shoal Creek (SH02, SH06, and SH07) were sampled only for sediments for sediment toxicity tests. Physical habitat at these sediment-only sites was judged to be inadequate to support mussel assemblages, and these sites were therefore excluded from quantitative mussel sur­ veys. These sites were included to avoid having long stream reaches without sediment toxicity data and to ensure that the selection of primary sampling sites included the full range of metal concentrations detected in the study area. Only two samples from the Neosho drainage were included as test sites in the sediment study: one from Tar Creek and one from the Neosho River downstream from Tar Creek. Quantitative mus­ sel surveys also included two samples that were not included in the sediment toxicity tests: a reference site in the Neosho drainage (NRSF) and a test site in the lower Spring River (SP19) that was excluded from toxicity testing because it was too close to mussel site SP20. Water Chemistry and Metal Concentrations About one-half of water samples from the reconnaissance sites sampled in 2017 (26 pore-water samples, 23 surface water samples) had measurable zinc concentrations (table 3). Filterable zinc concentrations in surface water and pore water were greater in samples from mining-affected tributaries than in samples from main stem sites in the Spring and Neosho Rivers (fig. 2). Median zinc concentrations in surface water and pore-water samples were 5.6 micrograms per liter (μg/L) as zinc and 8.1 µg/L as zinc, respectively (table 3), and several sites produced samples of surface water or pore water with fil­ terable zinc concentrations greater than 100 µg/L, includ­ ing Center Creek sites CC05D and CC16, Spring River sites SP06 and SP07, Shoal Creek site SH06, and Tar Creek site TC01. Lead (limit of detection, 0.1 μg/L) and cadmium (limit of detection, 0.05 μg/L) were detectable in fewer than onehalf of all surface-water and pore-water samples (44 percent of samples for lead; 28 percent for cadmium). We estimated ratios of zinc:lead (158–310) and zinc:cadmium (451–489) based on water samples in which all three metals were detected (fig. 3). Water chemistry did not differ widely among sites; median DOC concentrations were 3.6 and 6.0 milligrams per liter (mg/L) for surface water and pore water, respectively, and median hardness was 165 mg/L for both sample types. Samples from Tar Creek (TC01) had the greatest concentra­ tions of major ions in surface water and pore water, includ­ ing more than 30 mg/L for sodium, more than 900 mg/L for sulfate, and about 1,000 mg/L for total hardness (table 3). Sediment Metal Concentrations and Toxicity Hazards The 25 sediments selected for toxicity testing included 1 control sediment to characterize the performance of the test organisms and the exposure regime, 5 reference sediments to represent responses of mussels at low-metal sites, and 19 test sediments. Test sediments were further classified as tributary sites, which received direct input of metals from mining areas, and main stem sites on the Spring River and Neosho River, which were downstream from metal-contaminated tributaries (fig. 4). The control sediment and all five reference sediments had sum-PEQ hazard indices (for zinc, lead, and cadmium) of less than 1.0 (table 4), indicating that sediment toxicity was not likely. Sediments from the mining-affected reach of Center Creek had sum-PEQ indices as high as 7.3 (CC07), and the highest sum-PEQ index for Shoal Creek was 4.7 (SH05). Samples from the middle reach of the Spring River, which received inputs from Center Creek (and Turkey Creek—not sampled) had sum-PEQ indices ranging from 1.3 to 1.7 (SP10), and sediments from lower Spring River between Shoal Creek and the confluence of the Neosho River had sum-PEQ indices ranging from 1.7 to 3.4 (SP20). The Tar Creek site (TC01) had the highest sum-PEQ index measured in this study (8.8), but the main stem of the Neosho River below Tar Creek (NR06) had a sum-PEQ index within the range of the refer­ ence sites. Across all sites, zinc was the greatest contributor to the sum-PEQ index (median contribution, 72 percent of total); lead made up almost all of the remainder (median, 28 percent), and cadmium made a minimal contribution (<1 percent). Other measures of metal bioavailability and toxicity hazards indicated similar patterns across the sampling sites. Normalizing summed molar concentrations of five divalent metals (cadmium, copper, nickel, lead, and zinc), expressed as SEM, relative to AVS ([ΣSEM–AVS], expressed as micro­ moles per gram dry weight) did not greatly affect predictions of metal bioavailability, relative to predictions of the simple PEQ index, across the 25 sites (fig. 4B). Although copper and nickel were detected at relatively low concentrations and were not the focus of this study, analysis of concentrations of these metals in SEM extracts was necessary to evaluate the

2. Metal Concentrations in Water and Sediment and Status of Mussel Community    11 Table 2.  Primary sampling sites for mussel community surveys and sediment toxicity testing. [ID, identifier; sum-PEQ, summed probable effect quotient values (see footnote); Hwy, Highway; Rd, Road; Dr, Drive; SP, State park; --, no data or not appli­ cable; Lk, Lake; R, River] Laboratory ID Site Stream or reach Description Site type Mussel survey? Sediment toxicity? 1SumPEQ NF02 Upper Spring North Fork Reference Yes Yes SP01A Upper Spring Hwy 270 Reference Yes Yes SP03 Upper Spring Lawton Reference Yes Yes Center Hwy 110 Reference Yes Yes CC01G2 Center Chapel Reference Yes Yes CC05D Center Above Ben’s Branch Test Yes Yes Center Below Ben’s Branch Test Yes Yes Center Center Creek park Test No Yes Center Below Hwy JJ Test Yes Yes SP05 Middle Spring Below Center Test Yes Yes SP08 Middle Spring Below Turkey Test Yes Yes SP10 Middle Spring Above Empire Test Yes Yes SH02 Shoal Indian Springs Test No Yes SH05 Shoal Cherry Corner Test Yes Yes SH06 Shoal Nighthawk Rd Test No Yes SH07 Shoal Lime Kiln Dr Test No Yes SH10 Shoal Tipton Ford Test Yes Yes SH16A Shoal Below Schermerhorn Park Test Yes Yes SP12B Lower Spring Baxter Springs Test Yes Yes SP13 Lower Spring Blue Springs SP Test Yes Yes SP19 Lower Spring Hatchery main channel Test Yes No SP20 Lower Spring Hatchery side channel Test Yes Yes SP22 Lower Spring Above Grand Lk Test Yes Yes TC01 Neosho Tar Creek Test Yes Yes NRSF Neosho Stepps Ford Reference Yes No NR06 Neosho Below Tar Test Yes Yes SP02 Spring Spring R at Waco Control Yes Yes 1Probable effect quotient=sediment metal concentration/probable effect concentration (MacDonald and others, 2000); summed for zinc, lead, and cadmium. 2Sum-PEQ for mussel-only sites (SP19, NRSF) are X-ray fluorescence estimates from reconnaissance surveys (EcoAnalysts, Inc., 2018).

12    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District Table 3.  Concentrations of metals, major ions, and dissolved organic carbon in samples of surface water and pore water from reconnaissance sites. [ID, identifier; μg/L, microgram per liter; Zn, zinc; Cd, cadmium; Pb, lead; mg/L, milligram per liter; Na, sodium; Mg, magnesium; K, potassium; Ca, calcium; Fe, iron; DOC, dissolved organic carbon; Fl, fluoride; Cl, chlorine; SO4, sulfate; less than; --, no data; %, percent; nd, not determined] Site ID 1Metals (μg/L) 1Cations (mg/L) 1Anions and DOC (mg/L) 2Hardness (mg/L) Zn Cd Pb Na Mg K Ca Fe Fl SO4 DOC Surface water

Cc01G2 <20 <20 Cc05D Nf01 <20 Nf02 <20 Nr02

Nr03

Nr04

Nr05

Nr06

SH01 <20 SH05 SH06 SH07 <20 SH10 <20 SH14 SH15 SH16A SH18 SP01A <20 SP02 <20 SP03

Sp05 Sp06 Sp07

2. Metal Concentrations in Water and Sediment and Status of Mussel Community    13 Table 3.  Concentrations of metals, major ions, and dissolved organic carbon in samples of surface water and pore water from reconnaissance sites.—Continued [ID, identifier; μg/L, microgram per liter; Zn, zinc; Cd, cadmium; Pb, lead; mg/L, milligram per liter; Na, sodium; Mg, magnesium; K, potassium; Ca, calcium; Fe, iron; DOC, dissolved organic carbon; Fl, fluoride; Cl, chlorine; SO4, sulfate; less than; --, no data; %, percent; nd, not determined] Site ID 1Metals (μg/L) 1Cations (mg/L) 1Anions and DOC (mg/L) 2Hardness (mg/L) Zn Cd Pb Na Mg K Ca Fe Fl SO4 DOC Surface water—Continued SP08 SP09 SP10 SP11 SP12A SP12B SP13

Sp15

Sp18

Sp20 Sp21 Sp22 Tc01 1,028 Tc02

% de­ tected 53% 26% 35% 100% 100% 100% 100% 14% 19% 81% 81% 100% 100% Median nd nd nd nd nd nd Minimum Maximum 1,028 Pore water

Cc01G2 <20 <20 Cc05D Nf01 <20 Nf02 <20

14    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District Table 3.  Concentrations of metals, major ions, and dissolved organic carbon in samples of surface water and pore water from reconnaissance sites.—Continued [ID, identifier; μg/L, microgram per liter; Zn, zinc; Cd, cadmium; Pb, lead; mg/L, milligram per liter; Na, sodium; Mg, magnesium; K, potassium; Ca, calcium; Fe, iron; DOC, dissolved organic carbon; Fl, fluoride; Cl, chlorine; SO4, sulfate; less than; --, no data; %, percent; nd, not determined] Site ID 1Metals (μg/L) 1Cations (mg/L) 1Anions and DOC (mg/L) 2Hardness (mg/L) Zn Cd Pb Na Mg K Ca Fe Fl SO4 DOC NR02

Pore water—Continued NR03

Nr04

Nr05 Nr06

SH01 <20 SH05 <20 SH06 SH07 <20 SH10 SH14 <20 SH15 SH16A SH18 SP01A <20 SP02 <20 SP03

Sp05 Sp06 Sp07 Sp08 Sp09 Sp10 Sp11 Sp12A Sp12B Sp13 Sp15 Sp18

2. Metal Concentrations in Water and Sediment and Status of Mussel Community    15 Table 3.  Concentrations of metals, major ions, and dissolved organic carbon in samples of surface water and pore water from reconnaissance sites.—Continued [ID, identifier; μg/L, microgram per liter; Zn, zinc; Cd, cadmium; Pb, lead; mg/L, milligram per liter; Na, sodium; Mg, magnesium; K, potassium; Ca, calcium; Fe, iron; DOC, dissolved organic carbon; Fl, fluoride; Cl, chlorine; SO4, sulfate; less than; --, no data; %, percent; nd, not determined] Site ID 1Metals (μg/L) 1Cations (mg/L) 1Anions and DOC (mg/L) 2Hardness (mg/L) Zn Cd Pb Na Mg K Ca Fe Fl SO4 DOC SP20 Pore water—Continued SP21 SP22 TC01 TC02

% de­ tected 60% 30% 53% 100% 98% 100% 100% 35% 21% 81% 81% 100% 98% Median nd nd nd nd Minimum Maximum 1Samples filtered through 0.45-micrometer filter before analysis. 2Hardness expressed as milligrams per liter equivalent of calcium carbonate.

16    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District 1,000 Filterable zinc, in micrograms per liter Reference (n=12) Tributaries (n=15) Main stem (n=16) EXPLANATION Surface water Pore water Figure 2.  Zinc concentrations in reconnaissance samples of surface water and pore water, 2016–17. Boxes show median and upper and lower quartiles; error bars indicate 5th and 95th percentiles. Red solid and dashed lines are 12-week 20-percent effect concentrations for survival and growth, respectively, from water-only tests with juvenile mussels (Lampsilis siliquoidea). [n, number of samples] A. Pore water B. Surface water 1,000 1,000 Percentage of samples [Zinc:lead=310] [Zinc:cadmium=451] [Zinc:lead=158] [Zinc:cadmium=489] Filterable metal, in micrograms per liter EXPLANATION Zinc Lead Cadmium Figure 3.  Concentrations and ratios of zinc, lead, and cadmium in pore water and surface water from Tri-State Mining District streams, 2016–17. A, pore water; B, surface water. Symbols indicate percentage of samples with each metal at or less than specified concentrations, calculated at 5-percent intervals. Ratios of zinc:lead and zinc:cadmium in brackets are means for intervals where both metals of a pair were detectable. balance between concentrations of AVS and concentrations of divalent metals potentially controlled by AVS (Ankley and others, 1996). Concentrations of filtered zinc in pore-water peepers (PWP-zinc; fig. 4C) varied widely among sampling sites and indicated that the greatest zinc bioavailability was in Center Creek, upper Shoal Creek, and lower Spring River. Mussel Community Surveys Results of the mussel surveys reported by EcoAnalysts, Inc. (2018), are summarized in table 5 and figure 5. Species richness and measures of abundance and areal density of live mussels varied widely among sites but were generally higher at reference sites. The greatest density and species richness of mussels were detected at reference sites in the upper Spring River and North Fork Spring River. The remaining reference sites in Center Creek and Neosho River had mussel densities and species richness that overlapped with test sites, but all reference sites had some live mussels and at least three live species. Mussel assemblages at sites in Center Creek, middle Spring River, and Shoal Creek followed similar patterns, with lower mussel density and species richness in metalcontaminated stream reaches downstream from mining activ­ ity. The percentage of juvenile unionids differed widely among sites, from less than 1 percent at several sites in Center Creek and Shoal Creek to more than 50 percent at several sites in the middle reach of Spring River. The spatial pattern of Corbicula density was nearly opposite that of juvenile unionids, and Corbicula density exceeded 40 per square meter at several tributary sites. The predominance of juvenile mussels in the middle reach of the Spring River (SP08/SP10/SP12B) indi­ cates that recruitment at these sites benefitted from relatively favorable conditions, such as close connectivity with upstream reference sites, whereas repopulation of mussel assemblages in tributary sites may be limited by lack of unionid parental stock upstream and (or) lack of host fishes.

2. Metal Concentrations in Water and Sediment and Status of Mussel Community    17 B A Sum-PEQ (zinc, lead, cadmium) [ΣSEM–AVS], in micromoles per liter Zinc in pore-water peeper samples, in micrograms per liter 1,000 SP02 NF02 SP01 SP03 Site identifier SP05 SP08 SP10 SH02 SH05 SH10 SH07 SH06 SH16 SP12 SP13 SP20 SP22 TC01 NR06 SP02 NF02 SP01 SP03 SP05 SP08 SP10 SH02 SH05 SH10 SH07 SH06 SH16 SP12 SP13 SP20 SP22 TC01 NR06 SP02 NF02 SP01 SP03 SP05 SP08 SP10 SH02 SH05 SH10 SH07 SH06 SH16 SP12 SP13 SP20 SP22 TC01 NR06 Reference EXPLANATION Center Spring/Neosho Shoal Tar Figure 4.  Metal concentrations in sediment and pore water, fall 2017. A, sediment hazard quotients (summed probable effect quotient values [sum-PEQ]) for zinc, lead, and cadmium; B, summed molar concentrations of simultaneously extracted metals relative to acid-volatile sulfide ([ΣSEM–AVS]); C, zinc in pore-water peeper samples. Colors of bars indicate different stream segments. Dashed lines are geometric mean of 20-percent effect concentration estimates for juvenile mussels (table 10).

18    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District Table 4.  Physical and chemical characteristics of sediments from primary sampling sites, fall 2017. [ID, identifier; TOC, total organic carbon; %, percent; SEM, simultaneously extracted metals; μg/g, microgram per gram; Ni, nickel; Cu, copper; Zn, zinc; Cd, cadmium; Pb, lead; SEM–PEQ, simultane­ ously extracted metals probable effect quotient; PEQ, probable effect quotient; sum-PEQ, summed probable effect quotient values (sediment metal concentration/probable effect concentration); ΣSEM, sum of simultaneously extracted metals; μmol/g, micromole per gram; AVS, acid-volatile sulfide; [ΣSEM–AVS], summed molar concentrations of simultaneously extracted metals relative to acid-volatile sulfide; μg/L, microgram per liter] ID Site TOC (%) Clay (%) Silt (%) Total sand (%) Coarse sand (%) SEM (μg/g) SEM–PEQ Sum-PEQ ΣSEM (μmol/g) AVS (μmol/g) [ΣSEM– AVS] (μmol/g) Peeper Zn (μg/L) Ni Cu Zn Cd Pb Zn-PEQ Cd-PEQ PbPEQ (Zn+Cd+Pb) NF02 SP01A SP03 CC01G2 CC05D 1,780 1,310 1,670 SP05 SP08 SP10 SH02 1,570 SH05 1,680 SH06 1,080 SH07 SH10 SH16A SP12B SP13 SP20 1,170 SP22 TC01 3,420 NR06 −0.1 SP02

2. Metal Concentrations in Water and Sediment and Status of Mussel Community    19 Table 5.  Summary of qualitative and quantitative mussel surveys, 2016 and 2017. [CPUE, catch per unit effort; min, minute; SE, standard error of preceding mean; %, percent; --, no data] Site1 Qualitative survey Quantitative survey Live mussels Live species Total species CPUE (per 10 min) Live mussels Live species Total species Mean density per square meter, with standard error Juvenile (%) Mortalities (%) Total SE Adult SE Juvenile SE Corbicula SE NF02 SP01A SP01B SP02 SP03 SP05 SP06 SP08A SP08B SP10C SP12A SP12B SP13 SP19 SP20 SP22 CC01F CC01G2 CC05D CC10B SH02A SH05A SH06A SH07 SH08 SH10B

20    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District Table 5.  Summary of qualitative and quantitative mussel surveys, 2016 and 2017.—Continued [CPUE, catch per unit effort; min, minute; SE, standard error of preceding mean; %, percent; --, no data] Site1 Qualitative survey Quantitative survey Live mussels Live species Total species CPUE (per 10 min) Live mussels Live species Total species Mean density per square meter, with standard error Juvenile (%) Mortalities (%) Total SE Adult SE Juvenile SE Corbicula SE SH15A SH16A SH16B NR03 NRSF NR04 NR05 NR06 1Letters appended to site identifiers indicates that survey was conducted in one or more locations within the main site (as defined in table 1).

2. Metal Concentrations in Water and Sediment and Status of Mussel Community    21 NF02 SP01 SP03 NRSF SP05 SP08 SP10 SH05 SH08 SH10 SH16 SP12 SP13 SP19 SP20 SP22 NR06 Live species Site identifier NF02 SP01 SP03 NRSF SP05 SP08 SP10 SH05 SH08 SH10 SH16 SP12 SP13 SP19 SP20 SP22 NR06 Density (per square meter) A B Figure 5.  Characteristics of mussel assemblages at reference sites (hollow bars) and test sites (filled bars) from quantitative surveys. A, live mussel species; B, mussel density (mean and standard error). Data from EcoAnalysts, Inc. (2018).

22    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District Metals in Corbicula Soft tissues of Asian freshwater clams, C. fluminea (Corbicula), were collected from mussel survey sites to serve as a potential biomonitor of metal exposure in mussels. Angelo and others (2007) reported significant positive correlations of concentrations of zinc, lead, and cadmium in Corbicula tissues with tissues of unionid mussels from the same sites and indicated that Corbicula and unionids experienced similar levels of metal exposure. In the present study, concentrations of zinc, lead, and cadmium in Corbicula were significantly correlated with concentrations of each of these metals in sedi­ ment, although the strength of these associations was stronger for lead and cadmium (probability of type-1 error [p] <0.006) than for zinc (p=0.036; fig. 6). Other studies have had mixed success relating zinc in Corbicula tissue to ambient environ­ mental zinc levels, with weak relations attributed to differ­ ences in zinc bioavailability, especially in sediment (Bonnail and others, 2016), or to physiological regulation of internal zinc concentrations (Besser and others, 2007a). Concentrations of lead and cadmium in Corbicula also had significant nega­ tive associations with unionid species richness (fig. 6B, C). Although the weak correlations of zinc in sediment with zinc in Corbicula tissues limit the usefulness of these data for explaining the overall exposure and effects of metals on Tri-State Mining District mussel assemblages, these data sup­ port the hypothesis that concentrations of lead and cadmium, although much lower than concentrations of zinc, contribute to overall toxic effects on mussels in Tri-State Mining District streams. We evaluated the association of trends in mussel surveys with metal concentrations in sediment and water by comparing responses between groups of sites classified as “low metal” and “high metal.” Low- and high-metal groupings based on sediment metals were divided at a sum-PEQ index of 1.0, the same benchmark used to evaluate reference sites and a level assumed to be equivalent to the probable effect concentration (MacDonald and others, 2000). Waterborne metal concentra­ tions for each site were estimated from geometric means of zinc concentrations in stream water and push-point samples, with low- and high-metal groups separated by the value of 10 μg/L as zinc, which is the lowest concentration that was always detectable in reconnaissance samples and a concentra­ tion that corresponds closely to water-only toxicity thresholds described in the “Toxicity of Waterborne Metals in LongTerm Tests with Juvenile Freshwater Mussels (Spring 2017)” section. Live species richness and CPUE estimated from the qualitative survey and live species richness and mean areal density from the quantitative survey were consistently greater for the low-metal groups, based on either sediment or water (fig. 7), by factors from 1.9 to 8.2. The statistical significance of these differences was evaluated using the Wilcoxon test. For the qualitative survey data, live species richness was significantly greater in the low-metal sites based on zinc in water and sediment (type-1 error p=0.028 for both), and the CPUE was significantly greater in the low-metal sites based on zinc in water (p=0.028). For the quantitative survey data, no significant differences in richness or density were detected between the low- and high-metal sites based on sediment zinc, but richness and density were significantly greater in the lowmetal sites based on zinc in water (p=0.027). These analyses indicate that mussel community metrics are negatively associ­ ated with increases in metal exposure in streams of the study area. Trends in these data were consistent with data presented by Angelo and others (2007), who reported decreases in mus­ sel abundance and richness downstream from mining activities and significant relations between mussel community charac­ teristics and metal concentrations in sediment.

2. Metal Concentrations in Water and Sediment and Status of Mussel Community    23 A B D E F Corbicula zinc, in micrograms per gram Corbicula lead, in micrograms per gram Corbicula cadmium, in micrograms per gram 1,000 r 2=0.22 p=0.036 r 2=0.44 p=0.005 r 2=0.46 p=0.006 r 2=0.01 r 2=0.29 p=0.015 r 2=0.37 p=0.005 Live mussel species Live mussel species Live mussel species Corbicula zinc, in micrograms per gram Corbicula lead, in micrograms per gram Corbicula cadmium, in micrograms per gram SEM-cadmium, in micrograms per gram SEM-lead, in micrograms per gram SEM-zinc, in micrograms per gram 1,000 Figure 6.  Associations of zinc, lead, and cadmium bioaccumulation by Corbicula with metal concentrations in sediment (plots A, C, E) and with species richness of unionid mussels (plots B, D, F). A, Corbicula zinc; B, Corbicula zinc in live mussel species; C, Corbicula lead; D, Corbicula lead in live mussel species; E, Corbicula cadmium; F, Corbicula cadmium in live mussel species. Hollow shapes are site means. Dashed lines are linear regressions. [SEM, simultaneously extracted metals; r 2, coefficient of determination; p, probability; greater than]

24    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District A. Qualitative survey

B. Quantitative survey Sediment Water Sediment Water Live species Live species CPUE (mussels per 10 minutes) Density (mussels per square meter) EXPLANATION Species/low metal Species/high metal CPUE/low metal CPUE/high metal EXPLANATION Species/low metal Species/high metal Density/low metal Density/high metal Figure 7.  Comparison of mussel survey metrics between low- and high-metal sites. A, qualitative survey; B, quantitative survey. White bars are mean species richness; gray bars are mean abundance (catch per unit effort [CPUE] or density). For each pair, the hollow bar represents the low-metal range, and the hatched bar represents the high-metal range, defined by benchmarks for zinc in sediment (summed probable effect quotient value=1) or in water (zinc=10 micrograms per liter). A single asterisk indicates significant difference between the pair (probability [p] less than 0.05); a double asterisk indicates 0.05<p<0.10 (Wilcoxon test).

3. Toxicity of Waterborne Metals to Juvenile Freshwater Mussels    25 3. Toxicity of Waterborne Metals to Juvenile Freshwater Mussels The relatively low sensitivity of juvenile mussels in pre­ vious laboratory toxicity testing with sediments from the TriState Mining District (MacDonald and others, 2009; Besser and others, 2015) is inconsistent with evidence of widespread and severe effects on mussel assemblages (Angelo and others, 2007). These contrasts indicate that laboratory sediment toxic­ ity tests completed to date have not adequately represented the nature and duration of metal exposure experienced by juvenile mussels in the wild. Juvenile mussels primarily live and feed in surficial sediments (Yeager and others, 1994; Cope and others, 2008), but they may be exposed to varying combina­ tions of particulate metals (in sediment and food) and aqueous metals (in pore water and surface water), depending on their burrowing depth and feeding behavior (Kemble and others, 2020). The cumulative metal exposure during the juvenile life stage of freshwater mussels (which may last several years) is greatly underestimated by the standard 28-day sediment toxic­ ity test method (ASTM International, 2020). Recent research has indicated that the length of water-only toxicity tests with juvenile mussels can be extended from 4 weeks to 12 weeks, resulting in greater sensitivity (Wang and others, 2020). Such long-term laboratory exposures may better reflect the responses of juvenile mussels in the wild and thus may better predict changes in mussel assemblages in affected streams (Haag and others, 2019). Objective The objective of this study was to determine toxic­ ity thresholds for sensitive responses of juvenile mussels in 12-week water-only toxicity tests with zinc and with a three-metal mixture (zinc+lead+cadmium) at concentrations detected in streams of the study area. Methods Mussel Culture Gravid adult female fatmucket mussels (L. siliquoidea) were collected from Silver Fork Creek in Boone County, Mo., in fall 2016 and held in the USGS laboratory at reduced temperatures to delay the release of glochidia. In advance of planned toxicity tests, mussel brood stock was transferred to the laboratory of Chris Barnhart at Missouri State University in Springfield, Mo., for transformation of glochidia to the juvenile stage (Barnhart, 2006). Juvenile mussels were then shipped back to the USGS laboratory and reared for testing. Newly transformed juvenile mussels obtained from several females were reared at the USGS Columbia Environmental Research Center (CERC) to a shell length of about 1.0 mm (about 4 weeks). Toxicity Tests Water-only toxicity tests with juvenile mussels were com­ pleted after 12 weeks (Wang and others, 2020). Toxicity tests used intermittent-flow proportional diluters, which delivered five test solutions in a 50-percent dilution series, producing a sixteenfold range of nominal concentrations in each test. Filterable zinc concentrations in these tests broadly overlapped with zinc concentrations measured in stream water and surface (push-point) pore water during reconnaissance sampling (fig. 2). Test 1 was a zinc-only exposure and test 2 was an exposure with a mixture of zinc, lead, and cadmium, based on the ratios of these metals in the stream water and porewater concentrations. Stock solutions with zinc (stock 1) and lead plus cadmium (stock 2) were delivered simultaneously to a diluter by Hamilton syringe pumps (Hamilton, Reno, Nevada). Test 1 had a maximum nominal zinc concentration of 200 µg/L and a low nominal concentration of 12.5 μg/L as zinc. Test 2 included the same series of zinc concentrations plus lead concentrations ranging from 0.125 to 2.0 µg/L and cadmium concentrations ranging from 0.063 to 1.0 µg/L. Water-only toxicity tests used CERC well water diluted to a hardness of about 166 mg/L as calcium carbonate, consistent with the typical hardness of water measured in reconnaissance samples from the Spring River (table 3). Test temperature was maintained with a temperature-controlled water bath at 23 °C with temperature logged every 15 minutes using a HOBO U22 V2. Exposure chambers were 300-mL tall-form glass bea­ kers with a 2.5-cm overflow hole drilled in the side, such that each screened beaker held 200 mL of water. Overflow holes were covered with a stainless-steel screen (50 mesh/279-µm opening) glued with silicone caulk. Splitters on each of the 6 delivery lines fed 2 boxes, each containing 6 replicate bea­ kers (12 beakers per treatment). After day 28, when 4 rep­ licates from each treatment were destructively sampled, the 8 remaining beakers for each treatment were held in a single splitter box. Between days 0 and 28, the diluter cycled every half hour and the 1.0 L of test solution from each line was split 16 ways to deliver 63 mL per beaker per cycle (3.0 L per day), a replacement rate of 15 volumes per day. After day 28, the diluter cycled every hour, delivering 125 mL per beaker per cycle, to produce the same replacement rate. At the start of the test, 10 live mussels (visible foot move­ ment) were impartially transferred into each replicate beaker with 5 mL of silica sand (100- to 400-µm particles; Granusil no. 5010; Unimin Corporation, New Canaan, Connecticut), which was washed with well water overnight and rinsed with deionized water before use. Archive samples of 4 replicates of 10 mussels each were collected for measurements of initial lengths and dry weights. Dry weight of mussels was deter­ mined after being ovendried for 24 hours at 60 °C. During the tests, mussels were fed 3 mL of mixed algae (510-nanoliter cell volume per milliliter) twice per day. The algae mixture was prepared daily by adding 1 mL of Nannochloropsis concentrate and 2 mL of Shellfish Diet concentrate (Reed Mariculture, Campbell, California) into 1.8 L of water. Stocks

26    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District of algae were maintained in aerated cont4ainers at 4 °C and changed daily. Every 2 weeks, surviving mussels in each beaker were transferred into 200-mL glass dishes containing about 100 mL of test solution for determination of survival. Test Endpoints Test endpoints of survival, growth, and biomass were determined after 4 weeks and 12 weeks. Four of 12 repli­ cates per treatment were destructively sampled on day 28 for determination of length, weight, and biomass. Survival in all replicates was assessed by observing foot movement, either spontaneous or in response to gentle probing. Mussels without foot or shell movement, with an empty shell, or with a gaped shell containing swollen or decomposed tissue were recorded as mortalities and removed from test beakers. Surviving mus­ sels from destructively sampled replicates were preserved in 70 percent ethanol for measurements of shell length and dry weight. Shell length was determined by digital photographic analyses. Growth was calculated as the average dry weight of surviving animals in a replicate (total dry weight per number of survivors), and biomass was the total dry weight of surviv­ ing mussels for a replicate. Biomass reflects combined effects of survival and growth and is thus usually more sensitive than survival or growth alone. The biomass endpoint can be con­ sidered an estimate of the contribution of a cohort of juvenile mussels to secondary production at a site. Surviving mussels from the other eight replicates were transferred back into their beakers after the beaker was rinsed with test water and clean sand was added. These eight replicates were sampled for the same endpoints after 12 weeks. The test acceptability criterion (TAC) for the chronic mussel test was a minimum survival of 80 percent in the control after 4 weeks (ASTM International, 2015, 2020). Water Analyses Samples of water from water-only tests were analyzed for concentrations of filterable (<0.45 μm) zinc, lead, and cad­ mium on days 0, 28, 56, and 84 and for major ions and DOC on days 28, 70, and 84, as described in section 2. A minimum of three external calibration standards plus a calibration blank were used to calibrate the instrument response, and continuing calibration blanks and standards, laboratory spikes, duplicates, and other laboratory QA/QC procedures were used to verify instrument performance throughout the analyses (appendix 2). Data Analysis Toxic effects in the zinc-only test and the metal-mixture test were evaluated relative to measured filtered (0.45-µm) zinc concentrations. Chronic effect concentrations, including the no-observed-effect concentration and lowest-observedeffect concentration based on analysis of variance (ANOVA) were determined using SAS/STAT software (version 9.4). The 20-percent effect concentrations (EC20s) were determined based on concentration-response models (CRMs) generated by Toxicity Relationship Analysis Program software, version 1.30 (Erickson, 2010). Results and Discussion Metal Exposure Based on the data presented in section 2, the maximum nominal zinc concentration in the 12-week water-only toxicity tests was 200 µg/L, and the diluter produced nominal zinc concentrations as low as 12.5 µg/L as zinc. Measured concen­ trations of zinc, lead, and cadmium in reconnaissance samples were sorted to align samples by percentiles (fig. 3), and these distributions were used to estimate the ratios of the three met­ als at the concentration ranges where matched percentiles were available. The average zinc:lead ratio was 158 in pore water and 310 in surface water; the average zinc:cadmium ratio was 489 for pore water and 451 for surface water. These ratios were simplified to establish the ratios used in toxicity tests: zinc:lead=​200; zinc:cadmium=​400. Metal concentrations in tests were generally close to nominal, and measured zinc concentrations averaged 94 per­ cent of nominal in the zinc-only test and 92 percent of nominal in the metal-mixture test (table 6). Lead concentrations were detectable in the four highest treatments in the mixture study (nominals, 0.25–2.0 µg/L), and measured lead concentrations averaged 51 percent of nominal in these treatments. Cadmium concentrations were only detectable at the highest concen­ tration of the mixture test (nominal cadmium=1.0 µg/L), and measured cadmium in this treatment was 91 percent of nominal. Lead and cadmium concentrations less than ana­ lytical detection limits in the water-only test were expected, based on the distribution of measured metal concentrations in field-collected samples, but we added these metals at concen­ trations less than detection limits under the assumption that the three metals were at the same ratios across the full range of zinc concentrations. These results did not allow us to precisely measure exposure to lead and cadmium in all treatments of the metal-mixture study, but they indicate that mussels in this test were exposed to these metals at concentrations and ratios consistent with those measured in field-collected samples.

3. Toxicity of Waterborne Metals to Juvenile Freshwater Mussels    27 Table 6.  Toxicity of zinc and a three-metal mixture (zinc, lead, cadmium) in 12-week water-only tests with juvenile mussels. [20-percent effect concentration values and confidence limits are provided in table 10. Toxicity Relationship Analysis Program software models are provided in appendix 3, figure 3.1. Mean and standard errors (SE) are listed for days 28 and 84; μg/L, microgram per liter; n, number of samples; %, percent; mg, milligram; sig, significance; [ ], value less than detection estimated for analysis; ns, mean is not significantly less than the control; mean is significantly less than the control (Dunnett’s test); ANOVA, analysis of variance; p, probability; less than] Zinc (μg/L) n Survival (%) Growth (mean dry weight, mg) Biomass (dry weight, mg) Nominal Measured Mean SE Mean SE Sig Mean SE Sig Zinc test—4 weeks Control [3] ns ns ns ns

Anova

(p=0.078)

(p=0.0001)

(p=0.0001)

Zinc test—12 weeks Control [3] ns ns

Anova

(p=0.0014)

(p<0.0001)

(p<0.0001)

Mixture test—4 weeks Control [3] ns ns ns ns ns ns ns ns

ns

ns

Anova

(p=0.1627)

(p =0.0001)

(p =0.0001)

Mixture test—12 weeks Control [1.2] ns ns ns ns

ns

28    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District Table 6.  Toxicity of zinc and a three-metal mixture (zinc, lead, cadmium) in 12-week water-only tests with juvenile mussels.— Continued [20-percent effect concentration values and confidence limits are provided in table 10. Toxicity Relationship Analysis Program software models are provided in appendix 3, figure 3.1. Mean and standard errors (SE) are listed for days 28 and 84; μg/L, microgram per liter; n, number of samples; %, percent; mg, milligram; sig, significance; [ ], value less than detection estimated for analysis; ns, mean is not significantly less than the control; mean is significantly less than the control (Dunnett’s test); ANOVA, analysis of variance; p, probability; less than] Zinc (μg/L) n Survival (%) Growth (mean dry weight, mg) Biomass (dry weight, mg) Nominal Measured Mean SE Mean SE Sig Mean SE Sig Mixture test—12 weeks—Continued

Anova

(p=0.0001)

(p<0.0001)

(p<0.0001)

Toxicity of Zinc and the Three-Metal Mixture The 4-week and 12-week toxicity tests with zinc only and with the metal mixture had mean survival in controls ranging from 92.5 to 100 percent, exceeding the minimum TAC of 80 percent (table 6). Growth of mussels in the first 28 days was similar between the two tests (control mean dry weight per survivor: 0.563 milligram [mg] in the zinc test and 0.503 mg in the metal-mixture test), but growth was more rapid in the mixture test between days 28 and 84, and final control growth and biomass were more than twice as great in the mixture test. Statistical analyses of differences among treatments (ANOVA and Dunnett’s test on ranked data; table 6) indicated significant reductions in growth and biomass relative to controls after 4 weeks and 12 weeks in the zinc-only test and the metal-mixture test. Survival was not significantly reduced in the zinc-only test on either sampling date or in the metal-mixture test after 4 weeks, but survival at week 12 was significantly reduced in three treatments from the metal-mixture test. Concentration-response plots indicated clear negative trends for growth and biomass, with increased metal concen­ trations in both tests on both sampling dates, but fitting models to the data was difficult because decreased growth and bio­ mass were frequently evident at the lowest metal concentration greater than the control. Without one or more low test con­ centrations with responses close to the control, it is difficult to reliably estimate low effect concentrations, and regression models may predict low-level effects at concentrations less than the control. To avoid underestimating EC20s, growth and biomass data were fitted to two-parameter sigmoid regression models (SigmaPlot, version 14.0). These models estimate only the 50-percent effect concentration and the slope, while fixing the baseline response at the mean control response. Toxic effects of metals on juvenile L. siliquoidea differed between 4-week and 12-week exposures, as well as between zinc-only exposures and exposures to the three-metal mixture. Mussel survival was largely unaffected in the zinc-only test, with differences of 10 percent or less between high-zinc treat­ ments and controls or low-zinc treatments (table 6). Because effects on survival were minimal, thresholds for effects on growth (mean dry weight of individual) and biomass (total dry weight per replicate) were similar. Toxicity thresholds for growth and biomass on day 28 were similar for the zinc-only test (EC20s of 8.3 and 9.0 μg/L as zinc) and the metal-mixture test (EC20s of 9.7 and 8.3 μg/L as zinc; fig. 8). However, toxicity thresholds for the two tests diverged between week 4 and week 12, and EC20s for growth and biomass increased in the zinc-only test (to 16 and 20 μg/L as zinc, respectively) and decreased in the mixture test. By day 84, survival in the metal-mixture test had decreased significantly in the three highest treatments, with an EC20 for survival of 48 μg/L as zinc, and EC20s for growth and biomass had converged at 6.4–6.5 μg/L as zinc (fig. 8). These trends are consistent with the predominance of acute toxicity from zinc. The current EPA water quality criteria for zinc, based on extensive toxic­ ity testing, established a chronic criterion equal to the acute criterion. Based on this evidence, we would expect similar effect concentrations for the zinc-only test after 4 or 12 weeks, but toxicity in the zinc-only test decreased between 4 and 12 weeks. The unusual length of these studies introduces the possibility that surviving animals at the end of the 12-week test had become more tolerant of aqueous zinc than they were during the first 4 weeks of the test. The greater toxicity of the metal mixture after 12 weeks was presumably not caused by zinc alone but included some contributions from lead and cad­ mium, even though they were present at low concentrations in water. This hypothesis is supported by the differences in the chronic toxicity thresholds for these metals. When nominal concentrations of each of the three metals from the mixture test were converted to toxic units based on chronic water qual­ ity criteria (for example, 90 μg/L for zinc, 2.5 μg/L for lead, and 0.72 gram per liter for cadmium at a hardness of 85 mg/L; EPA, 2021b), summed toxic units for lead and cadmium were about equal to those for zinc.

3. Toxicity of Waterborne Metals to Juvenile Freshwater Mussels    29 EC20, in micrograms per liter as zinc EC20, in micrograms per liter as zinc

A. Zinc only Survival Growth Biomass Survival Growth Biomass B. Metal mixture EXPLANATION 4 Weeks 12 Weeks Figure 8.  Effect concentrations for water-only tests. A, with zinc only; B, with zinc in a three-metal mixture (zinc, cadmium, and lead). The 20-percent effect concentrations (EC20s) and 95-percent confidence intervals (error bars) for 4-week and 12-week exposures were estimated with Toxicity Relationship Analysis Program software (appendix 3, fig. 3.1). Symbol indicates the unbounded lower estimate of EC20. Dashed lines indicate zinc concentration of 10 µg/L.

30    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District Results of this study differed in some respects from results of recent zinc chronic tests completed with L. siliquoi­ dea in our laboratory (Wang and others 2020). The most direct comparisons are between the 12-week zinc-only test in the present study and these recent tests, which were not part of the Tri-State Mining District mussel project but were based on the same 12-week water-only test method used for the present study. Wang and others (2020) reported that L. siliquoidea were less sensitive to zinc early in the test, with EC20s of 66 µg/L for growth and biomass after 28 days—about eight­ fold greater (less toxic) than EC20s from the zinc-only test in the present study. However, 12-week EC20s for growth and biomass reported by Wang and others (2020) had decreased to 21–22 µg/L, close to 12-week EC20s from the present zinc-only test (16–20 μg/L as zinc). The primary differences in methodology between the two studies were that the test waters were different dilutions of CERC well water (hardness of 100 mg/L for Wang and others [2020] versus 166 mg/L in the present study) and that the previous test was stocked with larger juvenile L. siliquoidea, with a mean shell diameter of 1.55 mm at the start of the study, compared to 0.62 mm in the present test. The greater toxicity observed in the first 28 days of our zinc-only test may indicate greater sensitivity of the smaller/younger cohort of juvenile mussels to zinc, whereas the convergence of growth and biomass EC20s for the two studies over 12 weeks may represent convergence of zinc tolerance of the two test groups as they became older/larger. In a similar result, Wang and others (2010) reported that 96-hour median lethal concentrations for zinc to juvenile L. siliquoidea increased by a factor of 10 with an increase in the age of juve­ niles from 5 days to 6 months at the start of the test. Implications for Toxic Effects on Mussels These results indicate that the presence of waterborne lead and cadmium in the metal-mixture test resulted in greater chronic toxicity to juvenile mussels compared to the test with zinc alone. In 28-day tests, EC20s were similar for zinc alone and for the zinc-lead-cadmium mixture, with an over­ all range of EC20s from 8.3 to 9.7 μg/L as zinc. By day 84, sensitivity to zinc alone had decreased but sensitivity to the mixture had increased, as indicated by the onset of significant effects on survival and lower EC20s for growth and biomass (6.4–6.5 μg/L as zinc). These differences are not surprising because many studies have demonstrated that zinc is primarily toxic in short-term exposures (for example, up to 28 days in our tests) and that internal regulation of zinc tends to ame­ liorate longer-term toxic effects (Besser and others, 2007b; Mebane and others, 2007; Wang and others, 2010). The increased sensitivity of mussels to the metal mixture in the 12-week study presumably reflects contributions of lead and (or) cadmium to chronic toxicity, even though these metals were present at concentrations two orders of magnitude lower than zinc concentrations. Samples of surface water and pore water collected during summer 2016 frequently exceeded toxicity thresh­ olds for filtered zinc from either the zinc-only test or for the metal-mixture test. The frequency of exceedances of these zinc-based thresholds was about the same for surface-water and pore-water samples and about the same for tributaries and main stem sites. Overall, the frequency of exceedance of zinc thresholds was more than 75 percent for the more sensi­ tive growth thresholds from the metal-mixture test and about 50 percent for thresholds from the zinc-only test. Thresholds for survival—only estimated for the mixture test—were exceeded by an average of more than 25 percent of samples. The frequency of toxicity predicted by exposure to the metal mixture detected in waters of Tri-State Mining District streams was greater than the observed frequency of toxic effects on survival (9 percent of samples) and biomass (12 percent of samples) of juvenile L. siliquoidea in previous 4-week toxic­ ity tests with Tri-State Mining District sediment (Besser and others, 2015). Analyses of metals in Tri-State Mining District sediments and pore waters from the previous tests indicated that ratios of zinc, lead, and cadmium were similar to those in waters analyzed for the present study. The greater frequency of toxicity at relatively low metal concentrations in the wateronly tests may be explained by increased chronic effects in the 12-week exposure, compared to the previous 4-week sediment exposures. However, the lower frequency of toxicity in the previous sediment tests could also reflect low bioavailability of zinc and other metals in sediment pore waters. EC20s esti­ mated from pore-water metal concentrations during 12-week sediment toxicity tests (see section 4) were 4–40 times greater than water-only EC20s, indicating lower metal bioavailability in pore waters. The low thresholds for toxicity of waterborne zinc to juvenile mussels derived from our water-only toxicity tests were consistent with effects on wild mussels in Tri-State Mining District streams. Reconnaissance sampling at 31 sites in 2016 produced qualitative data on mussel assemblages and metal concentrations in stream water and shallow pore water, and 20 of these sites were sampled again in 2017 using quantitative methods (EcoAnalysts, Inc., 2018; table 4). Effects of waterborne metals on mussel assemblages were evaluated by comparing sites with low- and high-metal concentrations in water, with the threshold set at 10 μg/L, a concentration slightly greater than the EC20 for effects on growth in 12-week water-only tests. Sites with average zinc concentrations of 10 μg/L or less had greater average live species richness and greater abundance with density estimated by the CPUE (mussels per time) in the qualitative survey and from density (mussels per area) in the quantitative survey. Averages of these variables were consistently greater for sites in the low-metal group, and means for low-metal sites ranged from 51 to 500 percent greater than high-metal sites. Means for live species richness and CPUE from low-metal sites were significantly greater than means for high-metal sites (fig. 7; Wilcoxon test, p=0.01). Exceptions to this trend included sites SP05 and SP06, in Spring River downstream from Center

4. Chronic Toxicity of Metals in Sediments to Juvenile Mussels    31 Creek, and SH05, in upper Shoal Creek, which had greater richness and abundance of live mussels than other sites in the high-metal group. The relatively robust mussel communities at these high-metal sites may reflect differences in concentra­ tions of waterborne and sediment metals among these sites, for example, relatively low waterborne zinc at SH05 (16 μg/L) and relatively low sediment metals at SP05 (sum-PEQ=1.44). Differences in mussel assemblages among these three sites may also reflect other local conditions, such as differences in the accessibility of these sites for mussel recolonization; for example, sites SP05 and SP06 are a short distance downstream from high-quality communities at reference sites in the upper Spring River. The concordance of results from field sampling and labo­ ratory testing indicates that exposure of freshwater mussels to waterborne metals contributes to adverse effects on freshwater mussel communities in Tri-State Mining District streams. The toxicity predicted by EC20s from water-only tests is greater than that observed in previous exposures with metalcontaminated sediments, but these thresholds are consistent with trends in species richness and density of mussel assem­ blages. It also is possible that effects observed in previous short-term sediment toxicity tests did not adequately represent longer term effects of exposure to metal mixtures. Another fac­ tor that may have affected the relative sensitivity of juvenile mussels in our water-only tests compared to previous sediment tests was the difference in the age/size of juvenile mussels at the start of the tests. Juvenile mussels stocked in the previous sediment tests were larger (mean starting shell lengths of 1.5 and 2.1 mm) than the juveniles stocked into the present wateronly tests (mean starting shell length of 0.62 mm). 4. Chronic Toxicity of Metals in Sediments to Juvenile Mussels In 2007, the USGS characterized sediment toxicity and sediment chemistry at 70 sites across the Tri-State Mining District. This USGS study included laboratory toxicity tests with mussels (L. siliquoidea) and amphipods (H. azteca; Ingersoll and others, 2008; Besser and others, 2015). Toxic effects on both species increased with increasing metal concentrations in sediments, but mussels were less sensitive than amphipods in these tests. Results of mussel community surveys (Angelo and others, 2007) and sediment toxicity tests (Ingersoll and others, 2008) indicated that metals in surface water, sediment, and pore water likely affected mussels in the Tri-State Mining District (MacDonald and others, 2009). At the time of the previous study, options for characterizing toxicity of sediment and water to juvenile freshwater mus­ sels were limited to short-term (28-day) exposures, based on standard methods for water-only tests (for example, Wang and others, 2013; ASTM International, 2015). The 28-day exposure period limited the range of starting sizes of juvenile mussels that could be used for sediment tests because small juveniles do not grow in 28 days to a size that can be readily recovered from sediment, especially if test sediments contain large amounts of coarse particles (that is, medium and coarse sand-sized particles as many as 2 mm in diameter). Recently, our laboratory has successfully completed several water-only toxicity tests for durations of as many as 12 weeks (Wang and others, 2020). The development of this 12-week sediment toxicity method for juvenile mussels provided an opportunity to evaluate the longer term effects of exposure to sediments. Such exposures may be necessary to reliably characterize sensitivity of freshwater mussels to contaminated sediments during their juvenile stage, which may last for many months or years. Effects detected in long-term laboratory exposures of juvenile mussels to contaminants in sediment or water may explain changes observed in wild mussel communities. Recent studies reported that growth of juvenile mussels during 96-day in situ exposures in Kentucky streams was negatively asso­ ciated with chronic chemical stressors and that growth was consistently low in streams that had experienced declines of their mussel faunas (Haag and others, 2019). Objective The objective of this study was to document the fre­ quency and severity of toxic effects experienced by juvenile mussels during 12-week exposures to sediments represent­ ing the range of metal pollution in streams of the study area. Specific questions to be addressed with data from this study included the following: 1. Are responses observed using the 12-week mussel bioassays comparable to those observed in 6-week chronic bioassays with the amphipod H. azteca, the most sensitive species in previous sediment toxicity tests with sediments? 2. Does addition of low levels of metals (zinc, lead, and cadmium) that correspond to ambient concentrations in stream water and pore water result in changes in the responses of mussels that are concurrently exposed to metal-contaminated sediments? 3. What levels of metal exposure in sediments—with and without water column metal exposure—represent thresh­ olds for toxicity to mussels? Methods Sediment Toxicity Tests Two sediment toxicity tests were completed to determine the effects of sediment alone and sediment plus metals added to overlying water. Tests with amphipods and mussels were completed concurrently in fall 2017, and a second mussel test (with metals added to overlying water) was completed in spring 2018. Juvenile mussels were obtained from gravid

32    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District L. siliquoidea collected from Silver Fork Creek in Boone County, Mo., in spring 2017. Glochidia were transformed to juveniles at Missouri State University (Barnhart, 2006), and newly transformed juveniles were reared at the USGS labora­ tory to the appropriate age for testing. Amphipods (H. azteca) were obtained from continuous cultures at the USGS labora­ tory. Cohorts of neonate amphipods produced on the same day were reared to 6–7 days of age before the start of the test. One week before the start of each test, each sediment sample was homogenized and added to exposure beakers, clean overlying water was added, and test chambers were held under static conditions to allow the sediment to equili­ brate with oxygenated overlying water (Ingersoll and others, 2008; Wang and others, 2013). Water renewal was started on day 1. On day 0, groups of 10 test organisms were impartially transferred to each replicate exposure beaker. During stocking, 4 replicate groups of 10 organisms were set aside as archive samples to characterize the average dry weight of starting size for each test. In addition, all replicate groups of mussels were photographed before stocking, and shell lengths were determined by digital image analysis. Mean lengths of juvenile mussels stocked in the two sediment tests were close to the recommended starting size of 1.5 mm, although the mussels stocked in the 2017 test were smaller (1.22 mm average) than those stocked for the 2018 test (1.41 mm). Stocking was not consistent with the random order recommended by ASTM International (2020). Treatments were stocked in order based on laboratory-assigned sample numbers (table 4), with the control (sample 25) stocked first, followed by samples 1–24 in order. The sediment test completed in fall 2017 consisted of whole-sediment toxicity tests with juvenile mussels, L. sili­ quoidea (84 days), and with the amphipod H. azteca (42 days; ASTM International, 2020). This test used sediments (<2 mm) from 25 primary sampling sites. Clean overlying water, consisting of well water diluted to a hardness of 166 mg/L as calcium carbonate, was added at a rate of four volume replace­ ments per day using a proportional diluter system but without any added metals. In the second test in spring 2018, juvenile mussels produced by a different set of females were exposed to the same sediments and overlying water, except that the water added to the system contained a mixture of three metals at nominal concentrations of 10 µg/L for zinc, 0.08 µg/L for lead, and 0.025 µg/L for cadmium. This mixture was intended to represent a typical exposure regime for mussels in streams of the study area, based on results of field sampling of stream water and pore water in 2016 (fig. 1) and on results of wateronly toxicity tests. The mixture was expected to produce about a 10–20-percent reduction in the mussel biomass endpoint, based on 12-week water-only tests with zinc (see section 3) and based on similar tests reported by Wang and others (2020). The metal mixture was prepared in a 700-L reservoir and delivered via diluter delivery lines so that all test chambers would receive the same metal concentration during water renewals. Test conditions for sediment toxicity tests were based on standard and published methods (for example, Wang and others, 2013, 2020; ASTM International, 2015, 2020; Besser and others, 2015). The exposure system was maintained at 23 °C, with ambient laboratory light (about 500 lux) and a 16-hour light:8-​hour dark photoperiod. Test chambers were 300-mL tall-form beakers with screened overflow windows, which contained 100 mL of sediment and 175 mL of overlying water. The duration of sediment tests with freshwater mussels was extended from 28 days to 84 days (12 weeks), based on methods used for 28-day sediment tests (ASTM International, 2020) and 12-week water-only tests (Wang and others, 2020). Each test included a control sediment from the Spring River (SP02; Spring River at Waco, Mo.; Besser and others, 2015) and clean quartz sand (ASTM International, 2020). Mussels were fed 3 mL of algae mixture (510-nanoliter cell volume per milliliter) two times per day. Algae mixtures were prepared daily by adding 1 mL of Nannochloropsis concentrate and 2 mL of Shellfish Diet concentrate (Reed Mariculture, Campbell, Calif.) into 1.8 L of test water. Algae mixtures were kept in a refrigerator at 4 °C in aerated con­ tainers. Algae concentrations in the stock suspension and in selected treatments were measured using an algal cell counter before and after feeding. Amphipods were fed once daily with diatoms (Thalassiosira weissflogii [(Grunow) Fryxell & Hasle], 1200TM; Reed Mariculture, Campbell, Calif.) and with a suspension of flake fish food (Tetramin; Tetra Werke), and feeding rates were increased over time (ASTM International, 2020). Test Endpoints Mussel tests lasted 12 weeks, with endpoints of survival, growth, and biomass. After 4 weeks, mussels in each repli­ cate beaker were separated from sediment by wet sieving, the number of survivors was recorded, and digital images of surviving mussels were obtained for determination of shell length. Survival of juveniles was assessed by observing foot movement, either spontaneous or in response to gentle prob­ ing. Growth was calculated as the average dry weight (after drying at 60 °C for 24 hours) of surviving animals in a repli­ cate (total dry weight per number of survivors), and biomass was the total dry weight of surviving mussels for a replicate. Surviving mussels were then transferred to clean replicate beakers with fresh 100-mL parts of each sediment, which had been equilibrated for 1 week before this transfer as described for the beginning of the test. Mussels were transferred to clean sediment in beakers again after 8 weeks following the same procedure. After 12 weeks, remaining replicates were photographed for length determination and then destructively sampled for determination of survival, growth, and biomass as described previously for the 4-week samples. Amphipod tests were started with eight replicate test chambers for each sediment. After 28 days, sediments from all replicates were sieved to determine survival; survivors from four replicates were sampled destructively and weighed at

4. Chronic Toxicity of Metals in Sediments to Juvenile Mussels    33 60 °C to determine growth (mean dry weight per individual) and biomass (total dry weight per replicate). Survivors in the remaining four replicates were transferred into clean beakers with 5 mL of quartz sand substrate and clean test water and returned to the water delivery system. Starting on day 28, neonates in remaining replicates were counted and removed daily until the test was ended on day 42. On day 42, surviv­ ing amphipods from the remaining replicates were processed as described previously for day-28 samples, except surviv­ ing adults were preserved in formalin to allow determination of the number of males and females in each replicate before determination of dry weights. Reproduction data were used to calculate the total number of young produced per replicate and the total number of young per surviving female. Sampling and Analysis of Water and Sediment During toxicity tests, overlying water was sampled from all 24 sediments and the control during weeks 1, 4, 9, and 12 (2017 test) and weeks 1, 6, and 12 (2018 test) for analysis of zinc, lead, and cadmium. Monitoring of overlying waters also included temperature (daily); dissolved oxygen, pH, and conductivity (weekly); and ammonia, hardness, and alkalinity (expressed as milligrams per liter as calcium carbonate) at the beginning and end of each test. Samples of whole sediment and pore water were collected during mussel toxicity tests from extra replicate chemistry beakers for each sediment, which were treated the same as other replicates. At about day 60, small (2.9-mL) peeper samplers were deployed in all 96 chemistry beakers (2 tests×24 sediments×2 replicates) to sample pore-water metals. After a 10-day equilibration period, peepers were gently removed from sediments, and the contents of two replicate peepers were composited without dilution for metals analysis. After peepers were removed, sediments from chemistry beakers were sampled for analysis of SEM and AVS. Samples collected during tests were analyzed using the same methods for the pretest sampling described in section 2. Quality Assurance for Toxicity Tests Quality assurance in sediment toxicity tests was judged based on performance of the control sediment and the sand control (ASTM International, 2020, annex 5). The TAC for mussel survival after 12 weeks was assumed to be 80 per­ cent, the same TAC established for the 28-day sediment test. No TAC has been established for mussel growth in sediment tests. The sand control was used to assess the quality of test water and the adequacy of diets fed to mussels, independent of nutrients supplied by the control sediment. For tests using a reference envelope, growth for reference samples should not greatly exceed that in the sand controls to avoid setting too high a standard for performance of test sediments. Test acceptability for the amphipod test was judged by performance of the site SR control sediment with the following TACs: survival greater than or equal to (≥) 80 percent on days 28 and 42; mean dry weight ≥0.35 gram on day 28 and ≥0.5 gram on day 42; and reproduction ≥6 young per female (ASTM International, 2020, annex 2). Data Analysis Statistical analyses of survival, weight, biomass, and reproduction data were completed in accordance with require­ ments outlined by the EPA (2000) and ASTM International (2015, 2020). Responses of animals in test sediments were evaluated by comparison to the range of responses in reference sediments using a “reference envelope” approach (for exam­ ple, Wang and others, 2013; Besser and others, 2015; Steevens and others, 2020). Although various statistical methods have been used to compare data to a reference envelope (for exam­ ple, Hunt and others, 2001), we analyzed our sediment toxicity data by ANOVA with rank-transformed data, and we used Fisher’s least significant difference test to determine which test means were significantly less than the lowest reference mean. Relations of metal concentrations in sediment and water to responses of mussels in laboratory tests or in the field were examined using Toxicity Relationship Analysis Program (ver­ sion 1.30; Erickson, 2010) or SigmaPlot (version 14; SPSS Inc.). These models were used to estimate concentrations of metals or metal mixtures associated with various levels of biological effects (appendix 3). We used the EC20, the con­ centration associated with a 20-percent reduction of a speci­ fied endpoint, as the primary basis for comparing sensitivity among species and endpoints with the rationale that 20 percent represents a biologically meaningful level of effect that should be statistically separable from control or reference responses. Results and Discussion Characterization of Sediments and Pore Water Sediments collected from primary sampling sites in summer 2017 were characterized for physical and chemical characteristics that could affect the toxicity and bioavailabil­ ity of metals, and most of these analyses were repeated with sediments from chemistry beakers sampled during each test (table 7). Most characteristics of sediments tended to indicate high metal bioavailability and high toxicity risks. Zinc concen­ trations in the SEM fraction exceeded the zinc probable effect concentration (that is, PEQ index exceeded 1) in almost all test sediments, and the sum-PEQ index for zinc+lead+cadmium exceeded 5.0 in six of the test sites. Sediment characteristics that tend to ameliorate metal toxicity hazards were generally low, and most samples had low TOC (median of 0.77 percent) and had concentrations of SEM that exceeded the molar AVS binding capacity (median [ΣSEM–AVS] of 7.7 micromoles per gram). Sediment texture was dominated by sand-sized particles that tend to have low metal-binding capacity. Given these conditions, it is not surprising that zinc concentra­ tions in pore water were high, and 8 of 19 test sediments

34    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District Table 7.  Physical and chemical characteristics of sediment from the Tri-State Mining District during toxicity tests, 2017 and 2018. [ID, identifier; SEM, simultaneously extracted metals; μg/g, microgram per gram; Zn, zinc; Cd, cadmium; Pb, lead; PEQ, probable effect quotient (=SEM concentration/probable effect concentration); sumPEQ, summed probable effect quotient values; ΣSEM, sum of simultaneously extracted metals; μmol/g, micromole per gram; AVS, acid-volatile sulfide; [ΣSEM–AVS], summed molar concentrations of simul­ taneously extracted metals relative to acid-volatile sulfide; μg/L, microgram per liter; --, no data or not applicable; less than] Site SEM (μg/g) PEQ Sum-PEQ ΣSEM (μmol/g) AVS (μmol/g) [ΣSEM–AVS] (μmol/g) Peeper metals (μg/L) Zn Cd Pb Zn Cd Pb Zn Cd Pb 2017 test NF02 SP01A SP03 CC01G2 CC05D 2,060 1,430 1,550 SP05 SP08 SP10 SH02 1,400 SH05 1,670 SH06 1,100 SH07 SH10 SH16A SP12B SP13 SP20 1,100 SP22 TC01 3,460 NR06 SP02 Median

4. Chronic Toxicity of Metals in Sediments to Juvenile Mussels    35 Table 7.  Physical and chemical characteristics of sediment from the Tri-State Mining District during toxicity tests, 2017 and 2018.—Continued [ID, identifier; SEM, simultaneously extracted metals; μg/g, microgram per gram; Zn, zinc; Cd, cadmium; Pb, lead; PEQ, probable effect quotient (=SEM concentration/probable effect concentration); sumPEQ, summed probable effect quotient values; ΣSEM, sum of simultaneously extracted metals; μmol/g, micromole per gram; AVS, acid-volatile sulfide; [ΣSEM–AVS], summed molar concentrations of simul­ taneously extracted metals relative to acid-volatile sulfide; μg/L, microgram per liter; --, no data or not applicable; less than] Site SEM (μg/g) PEQ Sum-PEQ ΣSEM (μmol/g) AVS (μmol/g) [ΣSEM–AVS] (μmol/g) Peeper metals (μg/L) Zn Cd Pb Zn Cd Pb Zn Cd Pb 2018 test NF02 SP01A SP03 CC01G2 CC05D 2,100 1,100 1,500 SP05 SP08 SP10 SH02 1,500 SH05 1,700 SH06 1,000 SH07 SH10 SH16A SP12B SP13 SP20 1,100 SP22 TC01 3,100 NR06 SR Median

36    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District had pore-water zinc concentrations greater than 50 µg/L. Concentrations of zinc or metal mixtures in sediment (zincPEQ, sum-PEQ, and [ΣSEM–AVS]) did not differ consistently in test sediments between 2017 and 2018, but concentrations of PWP-zinc were less in the 2018 test than in the 2017 test. A total of 11 test sediments demonstrated decreases in PWP-zinc in the 2018 test (median decrease of 36 μg/L), and 8 sediments demonstrated increases (median increase of 19 μg/L; table 7). Mussel Tests The overall average length of mussels stocked in the 2017 test was 1.22 mm; however, when digitized images of the mussels stocked for these tests were measured, we observed a trend for decreasing mean lengths of mussels across the 25 sediments. Sediments stocked last had mean lengths that averaged 15–20 percent less than those stocked first (fig. 9A). Because the control and reference sediments were the first sediments stocked, this gradient in starting length could have affected the comparisons made to evaluate toxicity of test sediments, especially for sediments that were stocked with smaller mussels late in the sequence. The 2018 mussel test used a cohort of larger test organisms (mean day-0 length of 1.41 mm). This test was stocked in two cycles in an attempt to reduce selection bias. First, all test chambers were stocked with five mussels, then all chambers were stocked with a sec­ ond group of five mussels. Despite this change, we observed a similar decreasing trend for day-0 shell lengths in the order the sediments were stocked. Mean day-0 lengths in the 2018 test ranged from about 1.6 mm in the first sediment stocked to about 1.3 mm in the last sediments stocked. Controls in 2017 and 2018 mussel tests had a mean 12-week survival of 82.5 percent, exceeding the TAC mini­ mum of 80 percent. Growth of mussels at 12 weeks in the 2017 and 2018 tests was similar for the controls (SR; 3.9- and 3.8-mg mean dry weight, respectively), but mean growth in the sand control was lower in 2017 than in 2018 (2.7 mg versus 3.8 mg). Poor growth in the sand control could indicate that mussels were underfed, but this hypothesis is not sup­ ported by the consistency of mussel growth in the reference sites (lowest reference mean of 3.6 mg for both tests). Sediment Toxicity Survival and growth of juvenile mussels in the 2017 whole-sediment test differed substantially among reference and test sites. After 12 weeks, 14 of 19 test sites had a mean survival less than the lowest reference mean, and 4 test sites (CC07, CC09, CC10, SH16A) had a mean survival signifi­ cantly less than all means in the reference envelope (ANOVA/ Fisher’s least significant difference; table 8). Growth in both length and dry weight indicated fewer differences among sedi­ ments than survival, but most test sites had growth means at day 83 that were less than the range for the reference envelope (table 8). Overall, mean growth in dry weight was less than the reference envelope in 13 of 19 test sediments and means for 2 sediments (CC07 and SP20) were significantly reduced relative to all reference means (table 8). All but 1 of the 19 test sites had a mean biomass less than the reference envelope, and 12 test sites had means significantly less than all refer­ ence means. Results of the 2017 mussel sediment test may have been affected by variation in the starting size of juvenile mus­ sels because of nonrandom stocking. A significant trend for decreased shell lengths with increased stocking order was observed on day 0 (fig. 9A) and was still significant on day 28 (fig. 9C), so differences in growth among treatments during this exposure period cannot be attributed solely to metal expo­ sure. However, growth expressed as length increments (net increase in length) for the 12-week exposure period was not significantly correlated with stocking order (fig. 9D), indicat­ ing that effects of nonrandom stocking on mussel size were transient. Stocking treatments with smaller mussels could also have led to reduced survival because of either inherent greater sensitivity of smaller mussels or the technical difficulty of recovering smaller mussels at day 28; however, survival at day 28 was not significantly correlated with stocking order (fig. 9B), indicating no effect of starting size of mussels on survival. These results indicate that a significant bias in start­ ing size because of nonrandom stocking skewed treatments that started with smaller mussels toward lesser growth on day 28 but not at the end of the 12-week test. At the end of the 12-week exposure, sediments with reduced survival, growth, and biomass were clustered in several stream reaches most affected by mining activities: Center Creek (CC05D to CC10), Middle Spring River (SP05 and SP08), lower Shoal Creek (SH10 TO SH18A), lower Spring River (SP12B and SP13), and Tar Creek (TC01; fig. 10; table 8). Results of the 2018 mussel test indicated less variation of endpoints among test sediments (table 8). This difference was most noticeable for survival, and no 2018 test sediments had a mean survival less than the reference envelope. Mussel growth in dry weight differed among sediments, and the mean for 8 of 19 test sites was less than the reference envelope; the range of these differences in dry weight was relatively narrow, and none of the means for test sites were significantly less than all the reference means. Seven of the eight sediments with reduced growth in dry weight in the 2018 test also had reduced growth in dry weight in the 2017 test. The lesser toxicity observed in the 2018 mussel test was unexpected, considering the addition of the metal mixture to the overlying water of this test. We expected that addition of aqueous metals would produce an increase in toxicity rela­ tive to the sediment-only tests completed in 2017; however, metal analyses of overlying waters indicated that we had underestimated the potential for loss of dissolved metals to the reservoir, to the water delivery system, and (or) to the test chambers and sediments. Although nominal zinc concentra­ tions of 10 µg/L were automatically added to overlying water, measured zinc concentrations in overlying water were no greater than those measured in the 2017 sediment test (median

4. Chronic Toxicity of Metals in Sediments to Juvenile Mussels    37 B r 2=0.71 p<0.0001 r 2=0.67 p=0.0001 r 2=0.003 p=0.99 r 2=0.03 p=0.5 Day-0 length, in millimeters Day-28 survival, in percent D Length increment, days 0–83, in millimeters Day-28 length, in millimeters A 1 2 3 4 5 6 7 8 9 10 1112 131415 16 171819 20 2122 232425 1 2 3 4 5 6 7 8 9 10 1112 131415 16 171819 20 2122 232425 1 2 3 4 5 6 7 8 9 10 1112 131415 16 171819 20 2122 232425 1 2 3 4 5 6 7 8 9 10 1112 131415 16 171819 20 2122 232425 Stocking order Figure 9. Effect of stocking order (1=first stocked, 25=last stocked) on mussel endpoints in 2017 sediment test. A, length of mussels at stocking; B, survival of mussels at day 28; C, length of mussels at day 28; D, increment of length increase between days 0 and 83. Means (number of samples=4) are shown with linear regressions (dashed lines). r 2, coefficient of determination; p, probability; less than] zinc concentration of 4.6 µg/L in 2017 versus 3.2 µg/L in 2018). As a result, the 2018 study did not provide a test of our hypothesis about the contribution of metals in overlying water to toxic effects on juvenile mussels; however, the similar pat­ tern of growth results between the two tests indicated that the responses of mussels in the two tests were not simply random. The amphipod test met the TAC for survival (≥80 per­ cent) in the control and in all reference sediments on days 28 and 42. Amphipod growth did not meet the 28-day TAC of 0.35-mg dry weight per survivor in either the control or the reference sediments, but growth in the controls recovered to approach the 42-day TAC (0.47 mg versus 0.50 mg).

38    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District Table 8.  Toxicity of sediments to juvenile mussels in 12-week toxicity tests for sediment (fall 2017) and sediment-plus-water (spring 2018) tests. [Means with standard error (SE) and results of reference envelope (RE) analysis. Fisher’s least significant difference test (on ranks) was used for comparisons between test and reference sites. n, number of samples; mg; milligram; --, no comparison with reference; mean less than all reference means; nd, mean not different from reference; x, mean significantly less than all reference means; p, probability of type-1 error] Site Whole-sediment tests—Fall 2017 Site Sediment-plus-water test—Spring 2018 n Survival (of 10) Growth (mean dry weight, mg) Biomass (total dry weight, mg) n Survival (of 10) Mean dry weight (mg) Biomass (dry weight, mg) Mean SE RE Mean SE RE Mean SE RE Mean SE RE Mean SE RE Mean SE RE CTRL CTRL NF02 NF02 SP01A SP01 SP03 SP03 CC01G2 CC05D

nd

CC05D nd nd nd

nd

nd

nd

nd nd nd

nd

nd SP05

nd

SP05 nd

nd SP08

nd

SP08 nd nd nd SP10 nd nd SP10 nd nd nd SHO2 nd

SH02 nd

nd SHO5 nd

SH05 nd nd nd SH06 nd

SH06 nd nd nd SH07

SH07 nd nd nd SH10

nd

SH10 nd nd nd SH16A

SH16 nd

nd SP12B

SP12 nd nd nd SP13

SP13 nd

nd SP20

SP20 nd

nd SP22

SP22 nd

nd TC01

TC01 nd nd nd NR06 nd

TC06 nd

p<0.00001 p<0.00001 p<0.00001 p<0.00001 p=0.0002 p=0.0026

4. Chronic Toxicity of Metals in Sediments to Juvenile Mussels    39 Site identifier Site identifier Site identifier SR NF02 SP01 SP03 SP05 SP08 SP10 SH02 SH05 SH06 SH07 SH10 SH16 SP12 SP13 SP20 SP22 TC01 NR06

SR NF02 SP01 SP03 SP05 SP08 SP10 SH02 SH05 SH06 SH07 SH10 SH16 SP12 SP13 SP20 SP22 TC01 NR06 Survival, day 83, in percent A Biomass, day 83, total, in milligrams B SR NF02 SP01 SP03 SP05 SP08 SP10 SH02 SH05 SH06 SH07 SH10 SH16 SP12 SP13 SP20 SP22 TC01 NR06

Growth, day 83, mean dry weight, in milligrams Figure 10.  Mussel endpoints in reference sediments (hollow bars) and test sediments (filled bars) during sediment toxicity tests in 2017. A, survival; B, growth expressed as mean dry weight, in milligrams; C, biomass, in milligrams. Bars are means and error bars are standard errors (number of samples=4). Dashed horizontal line indicates lowest reference mean. Asterisk indicates mean is significantly less than the lowest reference mean (rank analysis of variance with Fisher’s least significant difference test, probability of type-1 error less than 0.05). Reproduction in the control sediment met the TAC of 6.0 young per female, but two of the five reference sediments had reproduction less than the TAC. Toxicity of sediments to H. azteca (table 9) followed trends similar to those in the 2017 mussel test. Sediments from the three sites with highest metal concentrations (Center Creek sites CC07 and CC09 and Tar Creek site TC01) indicated sta­ tistically significant decreases from reference sites for amphi­ pod survival on days 28 and (or) 42, and sediments from sites in middle Spring River (SP05, SP08) and Shoal Creek (SH07 and SH10) also had one or more endpoints below the refer­ ence envelope. All these sites also had reductions in one or more endpoints in the 2017 mussel test. Amphipod endpoints with the greatest number of reductions relative to the reference envelope were day-28 biomass (7 of 19 test sites) and day-42 survival (6 sites). Day-42 growth and biomass indicated minimal effects of sediments. Reproduction was low in several of the same sediments that had low day-42 survival, but only one site (CC07, which had no reproduction) fell to less than the wide reference range for reproduction. Reference sediment SP03 had unexplained low values for all day-42 sublethal endpoints, but even if this sample was excluded from the ref­ erence envelope, the (slightly) narrower reference range would identify toxic effects only in a few additional test sediments.

40    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District Table 9.  Toxicity of Tri-State Mining District sediments in 6-week exposures with amphipods, Hyalella azteca, fall 2017. [Means with standard error (SE) and results of reference envelope (RE) analysis. Fisher’s least significant difference test (on ranks) was used for comparisons between test and reference sites. ID, identifier; n, number of samples; %, percent; mg; milligram; --, no comparison with reference; ns, not significant; mean less than all reference means; x, mean significantly less than all reference means ANOVA, analysis of variance; p, probability of type-1 error] Site 4-week endpoints 6-week endpoints n Day-28 survival (%) Growth (mean dry weight, mg) Biomass (total dry weight, mg) n Day-42 survival (%) Growth (mean dry weight, mg) Biomass (total dry weight, mg) Young per female Mean SE RE Mean SE RE Mean SE RE Mean SE RE Mean SE RE Mean SE RE Mean SE RE CTRL NF02 SP01A SP03 CC01G2 CC05D ns ns ns ns ns ns ns

ns ns ns ns

ns ns ns ns

ns

ns

Sp05

ns

ns ns ns ns SP08

ns ns

ns ns ns SP10 ns ns ns ns ns ns ns SH02 ns ns ns ns ns ns ns SH05 ns ns ns ns ns ns ns SH06 ns ns ns ns ns ns ns SH07 ns ns

ns ns ns ns SH10 ns ns ns

ns ns ns SH16A ns ns ns ns ns ns ns SP12B ns ns ns ns ns ns ns SP13 ns ns ns ns ns ns ns SP20 ns ns ns ns ns ns ns SP22 ns

ns ns ns ns TC01

ns ns ns NR06 ns

ns ns ns ns ANOVA (p) p<0.0001 p<0.0001 p<0.0001 p<0.0001 p=0.06 p=0.0007 p=0.0005

4. Chronic Toxicity of Metals in Sediments to Juvenile Mussels    41 Results of the amphipod test are consistent with the trends seen in the other tests completed for the present study, with relatively severe short-term responses followed by moderating effects during longer exposures. Sediment Toxicity Thresholds We estimated EC20s for multiple endpoints from 12-week mussel tests and 6-week amphipod tests from CRMs based on several metrics of metal exposure. Survival, growth (mean dry weight), and biomass of juvenile mussels in the 2017 sediment test indicated decreasing trends with increasing metal concentrations in sediment or pore water. Concentrationresponse plots and EC20 estimates for 2017 test endpoints based on SEM-zinc concentrations are presented in figure 11. Survival and growth EC20s had relatively weak regression fits (r2=0.22), and biomass had a better fit (r compile2=0.48) and greater sensitivity. Regression models for the 2017 test that focused on endpoint responses over days 28–83 (weeks 4–12), to minimize effects of nonrandom stocking, had small differ­ ences in model fit (better for survival, worse for growth) and endpoint sensitivity (better for survival, worse for growth; appendix 3, fig. 3.3). Across all endpoints, the geometric mean of EC20s (GM–EC20) for SEM-zinc were 2.47 PEQ A B Survival, day 83 (of 10) Growth, day 83 mean dry weight, in milligrams Biomass, day 83, in milligrams EC20=2.71 r 2=0.22 EC20=3.82 r 2=0.22 EC20=1.5 r 2=0.43 SEM-zinc (PEQ) SEM-zinc (PEQ) SEM-zinc (PEQ) Figure 11.  Concentration-response plots and 20-percent effect concentrations (EC20s) for effects of sediment zinc on mussel endpoints in 12-week sediment toxicity tests in 2017. A, survival; B, growth expressed as mean dry weight, in milligrams; C, biomass, in milligrams. EC20s were estimated from three-parameter sigmoid regressions (SigmaPlot; appendix 3). Circles are site means. Solid lines are regressions. [SEM, simultaneously extracted metals; PEQ, probable effect quotient; r 2, coefficient of determination]

42    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District Table 10.  Comparison of toxicity thresholds for Hyalella azteca and Lampsilis siliquoidea estimated from tests with Tri-State Mining District sediment, 2007 and 2017. [EC20, 20-percent effect concentration; r2, coefficient of determination; %, percent; GM–EC20, geometric mean of EC20s for 2017 sediment tests; SEM, simul­ taneously extracted metals; Zn, zinc; PEQ, probable effect quotient (sediment metal concentration/probable effect concentration); --, no data or not applicable; [SEM–AVS], simultaneously extracted metals normalized to acid-volatile sulfide; μmol/g, micromole per gram; PWP–Zn, zinc in pore-water peepers; μg/L, microgram per liter; OW–Zn, zinc on overlying water; greater than; values in parentheses indicate r2; values in brackets indicate 95% confidence limits] Exposure Units 2017 EC20 with (r 2) or [95% confidence limits] threshold (GM– EC20) threshold 2007/2017 ratio Survival (12 week) Growth (12 week) Biomass (12 week) Reproduction 12-week mussel test (2017) SEM–Zn (0–12 week) PEQ (0.22) (0.22) (0.43) SEM–Zn (4–12 week) PEQ (0.39) (0.10) (0.43) SEM-mix Sum-PEQ −0.32 (0.26) (0.48) [SEM–AVS] μmol/g (0.44) (0.24) (0.68 PWP–Zn μg/L (0.18) (0.25) (0.46) OW–Zn μg/L (0.44) (0.25) Water Zn (Zn only) μg/L [12–21] Water Zn (metal mix) μg/L [45–51] [5–8] 6-week amphipod test (2017) SEM–Zn PEQ 3(0.13) 4(0.02) 5(0.68) 1T20=20-percent toxicity threshold from MacDonald and others (2009). 2Denotes 4-week period for survival. 3Denotes 4-week period for growth. 4Denotes 6-week period for biomass. 5Denotes 6-week period for reproduction. 6Denotes 2007/2017 ratio. (1,134 µg/g) for the full 12-week test models and 3.15 PEQ (1,446 µg/g) for the interval between 4 and 12 weeks (table 10), indicating that the sensitivity of the mussels was less during the later stages of the exposure but that the associa­ tion of responses with metals remained equally strong. CRMs based on SEM-zinc for the 42-day tests with H. azteca in 2017 are presented in figure 12 (for zinc PEQs) and appendix 3, figure 3.8 (for sum-PEQ, [ΣSEM–AVS], and PWP-zinc). Mussel growth, expressed as shell length, indicated a narrow range of differences among sediments in the 2017 test, and these data did not produce well-defined CRMs (results not shown). Additional CRMs for mussel endpoints in the 2017 test were based on other measures of sediment metal exposure: (1) the mixture of zinc, lead, and cadmium in the SEM fraction, expressed as sum-PEQ, and (2) SEM normalized to AVS ([SEM–AVS]; table 10; appendix 3, fig. 3.7). These models indicated consistent patterns across the different endpoints, with the lowest EC20s for biomass and highest EC20s for growth (fig. 13A). Although CRMs for sediment metals indicate low sensitivity of the growth endpoint, the consistently lower EC20s for biomass compared to survival indicate that reduced growth contributed to the sensitive biomass endpoint. Previous mussel sediment toxicity tests with sediments from the Tri-State Mining District in 2007 (Ingersoll and others, 2008) used sieved fine (<0.25-mm) sediments, which tend to have greater metal concentrations (and may therefore

4. Chronic Toxicity of Metals in Sediments to Juvenile Mussels    43 EC20=3.6 r 2=0.41 EC20=6.0 r 2=0.13 EC20=2.85 r 2=0.21 EC20=2.9 r 2=0.22 D. Reproduction, day 42 C. Biomass, day 42 A. Survival, day 28 B. Growth, day 28 Survival, in percent SEM-zinc (PEQ) Mean dry weight, in milligrams Biomass, in milligrams Young per female Figure 12.  Concentration-response plots and 20-percent effect concentration (EC20) estimates of EC20s for effects of sediment zinc on amphipod, Hyalella azteca, endpoints in 2017 test. A, day-28 survival; B, day-28 growth expressed as dry weight, in milligrams; C, day-42 biomass, in milligrams; D, day-42 reproduction. Circles indicate site means. Curved line is three-parameter sigmoid regression (SigmaPlot; appendix 3). [SEM, simultaneously extracted metals; PEQ, probable effect quotient; r 2, coefficient of determination] tend to be more toxic) than the bulk (<2-mm) sediment frac­ tion used for 2007 amphipod tests and for tests with species in 2017 and 2018. Thresholds estimated from the 2007 tests indicate that the sensitivity of juvenile mussels in fine sedi­ ments was an order of magnitude less than the sensitivity of amphipods in bulk sediments (6.4 PEQ for amphipods versus 64 PEQ for mussels; table 10). Sediments in the <0.25-mm fraction had greater average metal concentrations than those in the <2-mm fraction, but these differences were not great, with bulk sediments having concentrations of cadmium, lead, and zinc that averaged 75–84 percent of the concentrations in the fine fraction (Ingersoll and others, 2008). Mussel effect thresh­ olds (20-percent effect concentrations for survival [T20s]) [ΣSEM–AVS] in 2007 tests were calculated separately for the <0.25- and <2-mm sediment fractions, and toxicity thresholds estimated from these data also were similar for fine and bulk sediments (53 and 64 micromoles per gram, respectively). Given the relatively small differences in metal concentrations

44    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District A. Sediment EC20s B. Water EC20s SEM (PEQ) Survival Growth Biomass Survival Growth Biomass [SEM–AVS], in micrograms per liter 1,000 Filterable zinc, in micrograms per liter Zinc-PEQ Sum–PEQ [SEM]–AVS PWP-zinc Water zinc Water mix OW-zinc EXPLANATION EXPLANATION Figure 13.  Comparison of mussel 20-percent effect concentrations (EC20) among toxicity endpoints and among exposure metrics. A, in sediment; B, in water. [SEM, simultaneously extracted metals; PEQ, probable effect quotient; [SEM–AVS], simultaneously extracted metals normalized to acid-volatile sulfide; sum-PEQ, summed PEQ values for zinc, lead, and cadmium; PWP-zinc, zinc in pore-water peepers; water zinc, zinc in water-only test; water mix, zinc in water-only test with three-metal mixture; OW-zinc, zinc in overlying water]

4. Chronic Toxicity of Metals in Sediments to Juvenile Mussels    45 in these different sediment size fractions, it is likely that the greater effect concentrations for mussels in 2007 tests reflected lesser sensitivity of mussels to metals in shorter sediment exposures rather than differences in metal concentrations between sediment size fractions. Toxicity thresholds derived from measured zinc concen­ trations in pore water or overlying water indicated similarities and differences with sediment thresholds (table 10; fig. 13B). In general, thresholds for zinc in water, like sediment thresh­ olds, were greatest for survival and lowest for biomass; however, this pattern was sometimes obscured by differences among sample types, notably the roughly twentyfold differ­ ences in thresholds for pore waters (peepers) and overlying waters in 2017 sediment tests. These large differences indicate that one or both sample types did not accurately represent metal exposure by juvenile mussels. Differences in aqueous metal exposure may be mediated by burrowing behavior, which can either increase or decrease exposure to aqueous metals. Different sample types may also contain different mix­ tures of toxic metals, and results of the two water-only tests indicate lower thresholds for the three-metal mixture than for zinc alone, indicating a meaningful contribution of lead and cadmium to toxicity of the mixture. The 2018 mussel sediment-plus-water test had consis­ tent baseline responses for survival, growth, and biomass endpoints after 12 weeks, but 2018 endpoints did not indicate consistent decreasing trends at higher metal concentrations. Growth data from the 2018 test (table 8) indicated decreases in nine test sites relative to the reference envelope, and several sites approached 20-percent reductions in growth from the reference mean. All these sites also demonstrated growth less than the reference envelope in the 2017 test, but several sites with high metal exposure demonstrated greater growth in the 2018 test, preventing us from formulating a useful CRM. Toxicity thresholds from long-term tests with mussels and amphipods were generally lower than those estimated by MacDonald and others (2009) from 28-day toxicity tests. Thresholds for SEM-zinc are compared in table 10 among species, endpoints, and years. The GM–EC20 for SEM-zinc toxicity to mussel endpoints in 12-week tests in 2017 was more than twentyfold lower than the T20 values for 4-week mussel tests in 2007 (MacDonald and others, 2009), and the mussel GM–EC20 for SEM-zinc also was less than the GM– EC20 for amphipods (2.6 PEQ and 3.9 PEQ, respectively) in the 2017 tests. MacDonald and others (2009) recommended using thresholds based on toxicity of sediment zinc to H. azteca for ecological risk assessment of metals in Tri-State Mining District sediments, but results of our 12-week mussel tests produced GM–EC20s for mussels that were lower than either T20s for amphipods from 28-day tests in 2007 or GM– EC20s for amphipod 6-week tests in 2017 (table 10). These results support our hypothesis that thresholds for toxicity to L. siliquoidea estimated from 12-week exposures are lower than those estimated from short-term tests and that juvenile mussels are equally sensitive or more sensitive than amphipods to toxic effects of metals in sediments. Differences Between 2017 and 2018 Mussel Tests The two sediment tests demonstrated unexpected differ­ ences in sensitivity despite testing the same sediment samples and the same test organism. These differences could be due to differences in the quality or sensitivity of test organisms or to differences in metal exposure between the two tests. The use of smaller mussels to stock the 2017 test also could have contributed to the greater sensitivity observed in that test. Differences in starting size of mussels among treatments because of nonrandom stocking could have affected differ­ ences in growth early in the 2017 study, but trends for reduced survival, growth, and biomass among treatments after the full 12-week exposure were not closely correlated with the starting size of mussels. Differences in the starting size of juvenile mussels among 2018 treatments were similar to those observed in the 2017 test, but these differences did not have any effect on toxic effects in long-term tests. Mussels exposed to sediments in the 2017 and 2018 tests indicated reductions in growth, relative to reference sediments, in many sediments in the same range of sediment metal exposure. Many of the same sediments that were toxic in the 2017 mussel test also were the most toxic sediments in the amphipod test, indicating that differences in toxicity among treatments in the 2017 tests were predominantly related to metal exposure, not to effects of stocking bias. The absence of clear concentration-response trends in the 2018 study may reflect differences in metal exposure, notably decreases in PWP-zinc during the 6-month storage period between tests for several sediments that had indicated high toxicity in the 2017 test. The apparent differences in sensitivity of mussels between the 2017 and 2018 tests did not correspond to differences in sediment metal concentra­ tions because the sum-PEQ index averaged about 25 percent greater for sediment from the 2018 test compared to 2017 and the [ΣSEM–AVS] value was nearly identical across the two tests. However, PWP-zinc averaged about 30 percent lower in the 2018 test. This difference was strongly affected by several sites (CC07, CC10, and SP20) that had high PWP-zinc concentrations and had significant reductions in growth and (or) survival in the 2017 test. Thus, it seems that the lack of sensitivity to sediments in the 2018 test may reflect reduced metal sensitivity because of the greater starting size of the mussels and (or) changes in bioavailability of sediment metals during storage. Another factor that could have affected differences in sen­ sitivity between the 2017 and 2018 tests was the addition of metals into overlying water during the 2018 test. Our expecta­ tion was that these added metals would be physically and bio­ logically available, for example, as free dissolved metals that would be bioavailable to juvenile mussels via the gill or other body surfaces. However, it is possible that these added metals were either not physically available (for example, if juvenile mussels remained buried in sediment and were not exposed to additional metals in overlying water) or not bioavailable

46    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District (for example, because of sorption to surfaces of sediments or test chamber surface or partitioning to resistant phases such as metal sulfides). Either of these scenarios could account for the reduced toxicity observed in the 2018 test, although the amount of metals added to the overlying water (nominal, 10 µg/L as zinc) seems too small to stimulate changes in mussel burrowing behavior or to substantially shift chemical equilibria. The amount of metals added to overlying water in the 2018 test was much less than the added metal concentra­ tions (20–80 μg/L as zinc) associated with reversals of toxicity in the 2019 supplemental test (see section 5). Toxicity of Metals in Water and Sediment Results of whole-sediment toxicity tests with and wateronly toxicity tests that simulated exposure to waterborne met­ als in Tri-State Mining District streams (see section 3) were consistent with the hypothesis that mussels in Tri-State Mining District streams are adversely affected by long-term exposures to metals. Our attempts to better understand the relative con­ tributions of exposure to metals in sediment versus overlying water to adverse effects on juvenile mussels were unsuccess­ ful, and we have only indirect evidence of the relative contri­ butions of metal exposure via sediment and water exposure routes. However, the nature of the effects we observed in separate tests with metals in water and sediment give us some clues about the similarity of the exposure routes and mecha­ nisms of toxicity. Exposure to waterborne zinc alone pro­ duced short-term effects on growth that decreased in severity between 4-week and 12-week exposures. Water-only expo­ sures that combined zinc with low levels of lead and cadmium, based on concentrations measured in streams, also produced short-term effects on mussel growth, but these effects became more severe during the 12-week tests; mussels in the metalmixture tests also experienced reduced survival and more severe reductions in biomass than those observed in the zinconly test. Toxic effects in our 12-week sediment tests consisted of reduced survival, especially during the first 4 weeks of exposure, and reduced growth that persisted over the duration of the study. The limited severity of the toxic effects observed in sediment tests may reflect our intentional selection of sam­ pling sites that supported mussel assemblages and tended to have relatively low sediment metal concentrations. We were unable to document increased toxicity because of added metals in overlying water in the 2018 sediment test, but zinc concentrations in overlying water of sediment tests exceeded EC20s for mussel growth and biomass (6.4–6.5 μg/L as zinc) determined from the water-only test with the threemetal mixture. The lowest threshold based on zinc in pore water peepers (80 μg/L as zinc for effects on biomass) was an order of magnitude greater than the water-only thresholds, but the frequency of exceedance of this threshold was the same as for the water-only threshold (5 of 19 test sites). These similar patterns of exceedance of water-based thresholds, despite differences in the zinc concentration in pore water versus overlying water, indicate that contaminated sediments are the dominant exposure route and that water-based threshold tests simply reflect aqueous concentrations in equilibrium with sediment thresholds. Sensitivity of Different Sizes of Juvenile Mussels Differences in sensitivity of different cohorts of juvenile L. siliquoidea to metal toxicity may be attributable to differ­ ences in starting age/size. Decreasing sensitivity with increas­ ing age/size is a common phenomenon in ecotoxicology, and Wang and others (2010) reported tenfold increases in acute zinc median lethal concentrations for L. siliquoidea between tests with newly transformed (5-day-old; about 0.4 mm) juveniles and tests with 6-month-old (2.07 mm) juveniles. The water-only tests described in section 3 were started with smaller juvenile mussels than are currently recommended for sediment testing because these tests included only a clean sand substrate, and small size did not cause difficulties for recovering either live mussels or shells of dead mussels. The EC20s estimated from these tests were as low as any previous 12-week water-only tests with juvenile mussels (Wang and others, 2020). Also, differences in sensitivity between the 2017 and 2018 sediment tests may reflect the use of two size ranges of juvenile mussels for these tests. It also is possible that unin­ tentional stocking of different-sized mussels across different treatments of the sediment tests affected their response to sedi­ ment metals, but we see little evidence of this phenomenon. The unanswered questions that remained after these tests led us to design a supplemental study to address more directly the question of water-plus-sediment exposures (see section 5). These studies provided another opportunity to address the question of differences in sensitivity with size and age of juve­ nile mussels. Finally, the selection of sediments for the supple­ mental studies allowed us to evaluate if other characteristics of sediments other than toxic metals might affect the survival and growth of juvenile mussels. 5. Toxicity of Metals in Sediment and Water to Size Classes of Juvenile Mussels Background This objective of this supplemental study was to investi­ gate the relative contribution of waterborne and sediment met­ als to toxic effects on juvenile mussels, which was one of the original research questions of this project that was left unan­ swered because of an unsuccessful attempt to enrich overlying waters with metals during the 2018 sediment toxicity tests. The design for the supplemental experiments included 5 levels of an aqueous metal mixture and 4 sediment treatments. Waterborne metals were added in a 50-percent dilution series (plus a control) that delivered a sixteenfold range of nominal

5. Toxicity of Metals in Sediment and Water to Size Classes of Juvenile Mussels    47 metal concentrations. Sediments were prepared by composit­ ing parts of Tri-State Mining District sediments remaining after about 2 years of storage. Sediment composites included sediments with low, medium, and high sediment metal con­ centrations in combination with low and high percentages of coarse sand (particle diameters of 0.5–2.0 mm). Objectives The supplemental studies were designed to address the following research questions raised by results of sediment toxicity tests in 2017 and 2018. 1. What levels of metals added to overlying water will sub­ stantially increase toxicity in combination with metals in sediment and pore water? 2. Does the starting age/size of juvenile mussels affect their sensitivity to metals in sediment and overlying water? 3. Do juvenile mussels respond to sediment physical characteristics (for example, high coarse sand content) in addition to high metal concentrations? Methods Experimental Design In fall 2019, four sediment composites were pre­ pared from stored sediments collected in summer/fall 2017 (tables 11 and 12). Sediment 1 was a reference low-metal sediment dominated by fine sediment particles (sum-PEQ=1.4; coarse sand=4.3 percent). Sediments 2 and 3 had intermediate metals concentrations (sum-PEQ of about 3) but differed in sediment texture (sediment 2, 9.3 percent coarse sand; sedi­ ment 3, 58 percent coarse sand). Sediment 4 had high metals (sum-PEQ=7.8) and a coarse texture (43 percent coarse sand). Two tests with these sediments were completed in fall 2019: test A was stocked with juvenile mussels with average starting shell lengths of 0.91 mm (about 1 month old), and test B was stocked with larger juveniles from the same maternal cohort with average starting lengths of 1.68 mm (about 2 months old). Diluters delivered water containing five levels of a metals mixture (zinc/lead/cadmium) in the same ratios used in the water-only test (see section. 2) plus a control water. Waterborne metal treatments 1 (low) to 5 (high) were sepa­ rated by factors of 0.5, and the nominal high zinc concentra­ tion was 80 µg/L for test A and 160 µg/L for test B (table 12). These solutions were automatically prepared in diluted well water (hardness of 166 mg/L) and delivered to test chambers at the rate of four volume replacements per day. Toxicity and Chemistry Methods Short-term (28-day) tests were based on ASTM International (2015) and the methods from previous stud­ ies of effects of mining-derived metals on juvenile mussels (MacDonald and others, 2009; Besser and others, 2015). This approach is similar to terminating the 12-week sediment tests described in section 4 after 4 weeks. Longer exposures could not be completed because they would have required replace­ ment of sediment in test beakers, and the supply of sedi­ ments was limited. Samples of sediment and pore water were collected from separate chemistry beakers during the test and were analyzed for SEM, AVS, and pore-water (peeper) zinc as described in section 4. Data Analysis Data from tests A and B were analyzed by separate twoway ANOVAs using rank-transformed data. These analyses evaluated the significance of overall differences among the 4 sediment treatments and among the 6 water treatments, as well as the interactive effects of sediment and water treat­ ments. Means for toxicity endpoints in sediment treatments were compared with sediment 1 and means of water treat­ ments were compared with the water control. Results and Discussion Chemistry Metal concentrations in sediment composites and overly­ ing waters of 2019 tests were close to targets (tables 11 and 12). Sediment 4 had concentrations of zinc-PEQ, sum-PEQ, and [ΣSEM–AVS] about 2.5 times greater than those in sediments 2 and 3. Zinc concentrations in pore waters cor­ responded closely to concentrations in sediment, with little or no increase of pore-water metals in treatments with high metal concentrations added to overlying water. Zinc concentrations in pore water of sediments 2 and 3 varied by more than a factor of two among water treatments, but this variation was independent of overlying water treatments. In contrast to our previous attempts to control metal concentrations in overly­ ing water, overlying water in test chambers remained close to nominal concentrations delivered by the diluters (fig. 14). Zinc concentrations in overlying water indicated only a minor effect of pore-water zinc, which was noticeable only in treatments with high pore-water zinc and low zinc in overlying water.

48    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District Table 11.  Physical and chemical characteristics of sediment and pore water in supplemental tests, fall 2019, for pretest chemistry. [TOC, total organic carbon; %, percent; SEM, simultaneously extracted metals; Zn, zinc; Pb, lead; PEQ, probable effect quotient (sediment metal concentration/probable effect concentration); Cd, cadmium; sum-PEQ, summed probable effect quotient values; ΣSEM, sum of simultaneously extracted metals; μmol/g, micromole per gram; AVS, acid-volatile sulfide; [ΣSEM–AVS], summed molar concentrations of simultaneously extracted metals relative to acid-volatile sulfide; DOC, dissolved organic carbon; mg/L, milligram per liter; Na, sodium; Mg, magnesium; K, potassium; Ca, calcium; Mn, manganese; Cl, chlo­ rine; SO4, sulfate] Sediment characteristics Sediment metals Pore-water chemistry Sediment TOC (%) Clay (%) Silt (%) Total sand (%) Coarse sand (%) SEM– Zn SEM– Pb (PEQ) SEM– Cd (PEQ) SumPEQ (SEM) ΣSEM (μmol/g) AVS (μmol/g) [ΣSEM– AVS] (μmol/g) DOC Hardness (mg/L) Na (mg/L) Mg (mg/L) K (mg/L) Ca (mg/L) Mn (mg/L) (mg/L) SO4 (mg/L)

5. Toxicity of Metals in Sediment and Water to Size Classes of Juvenile Mussels    49 Table 12.  Physical and chemical characteristics of sediment and pore water in supplemental tests, fall 2019, for midtest chemistry. [Zn, zinc; μg/L, microgram per liter; nom., nominal; sed., sediment; less than; greater than] Zn in pore water (μg/L) Zn in overlying water (μg/L) Water Test A Test A Test B Nom. Sed. Sed. Sed. Sed. Nom. Sed. Sed. Sed. Sed. Nom. Sed. Sed. Sed. Sed.

Toxicity Survival of small juvenile mussels in test A was unaf­ fected by sediment or water treatments (fig. 15A). In test B, survival in sediment 2 (medium metals/fine texture) was significantly reduced relative to the reference sediment (sedi­ ment 1); this reduction was most evident in water treatment 1 (nominal, 10 μg/L as zinc), and survival gradually increased at higher water treatments. Growth in test A was significantly reduced in sediment 2 compared to the reference sediment across all six water treatments (fig. 15C). Growth also was reduced in sediment 4 (high metals, coarse texture) at the three lowest water treatments of test A, but growth increased at the two highest water treatments, reaching levels similar to sediments 1 and 3. Growth effects were less affected in test B, but overall growth in sediment 2 was significantly less than growth in the reference sediment. Growth in test B followed similar complex patterns among all sediment treatments: level or decreasing growth in water treatments as much as 40 µg/L, increased growth in the 80-µg/L treatment, and sharply decreased growth in the highest water treatment (160 µg/L). Concentration-Response Relations These results indicate that metals in sediment and overly­ ing water affected toxicity to juvenile mussels. The clearest example of the effect of sediment metals was the growth data for sediment 2 in test A, where growth was reduced relative to most other sediments but indicated no differences across the six water treatments (0–80 μg/L as zinc; fig. 15C). Variation of toxicity with metal concentrations in water did not follow simple concentration-response trends. Data for all endpoints that indicated toxicity also indicated “reversals”—decreases in toxicity with increasing metal concentrations in water. This pattern was most prominent in test A for effects of growth in sediment 4, which were as low as sediment 2 at low water concentrations but increased to near the reference sediment at higher water metals treatments (fig. 15C). Similar reversals of toxicity were evident for growth data from both tests and for survival data for sediment 2 in test B. These reversals were across similar ranges of waterborne metals concentrations (20–40 μg/L as zinc in test A; 40–80 μg/L as zinc in test B), and toxicity (especially decreased growth) resumed at concen­ trations greater than these concentration ranges. These unusual reversals of concentration-response trends indicate that exposure to increasing metal concentrations in overlying water elicited a protective response in the juvenile mussels. One protective mechanism known to be active in freshwater mussels is the induction of metallothionein (MT), a low molecular weight protein that has a strong affinity for binding divalent metals. Induction of MT can exist in unionid mussels and a wide range of other aquatic taxa, consistent with its role in regulation of essential metals (Amiard and others, 2006). Studies with the unionid Anodonta grandis (Say, 1829) documented that MT can be induced in mussels by exposure to toxic metals (Couillard and others, 1993). Conditions lead­ ing to MT induction vary widely among published studies, but induction of MT has often been associated with cadmium exposure (Amiard and others, 2006). Once MT is induced, its high cysteine content can bind strongly to cationic metals including cadmium, zinc, and lead, a process that can lead to reduced toxicity (Couillard and others, 1993). The rapid­ ity of MT induction varies widely among published labora­ tory and field studies, but induction over periods as short as several days has been reported (Amiard and others, 2006). In studies with C. fluminea (Correia and others, 2002), MT was induced by exposure to waterborne copper, but not by copper-spiked sediments, resulting in concentration-response trends like those observed in the present study. Aqueous cop­ per was toxic at low concentrations, with no induction of MT. At intermediate copper concentrations, toxicity was detected within 1 day but decreased to control levels after several days, concurrent with maximum induction of MT. At the highest copper concentrations tested, toxicity was detected without increases in MT. This response at high metal exposure levels has been described as “spillover” because MT binding capac­ ity and biosynthetic capacity are exceeded and toxic effects ensue (Baudrimont and others, 1999, Correia and others,

50    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District A B Filterable zinc, in percent of nominal Filterable zinc in control, in micrograms per liter Sediment composite EXPLANATION Test 1 Test 2 EXPLANATION Test 1 Test 2 Figure 14.  Filterable zinc in overlying water of supplemental (sediment-plus-water) toxicity test, 2019. A, zinc as a percentage of nominal (bars means and error bars standard deviations for the three highest water treatments); B, mean zinc concentrations in water controls (no waterborne zinc added) of four sediment treatments.

5. Toxicity of Metals in Sediment and Water to Size Classes of Juvenile Mussels    51 Water treatment (nominal zinc), in micrograms per liter A. Survival of small mussels (test A) B. Survival of large mussels (test B) C. Growth of small mussels (test A) D. Growth of large mussels (test B) Mean survival (of 10) Growth, in milligrams Growth, in milligrams Mean survival (of 10) ANOVA Sediment (none) Water (none) ANOVA Sediment: 2 (−) Water: (none) ANOVA Sediment: 2, 4 (−) Water: 4 (+) ANOVA Sediment: 2 (−) Water: 4 (+) Figure 15.  Survival and growth of two age/size classes of juvenile mussels in 28-day exposures to 4 composite sediments (shown as numbered and colored symbols) and 6 metal treatments in overlying water, fall 2019. A, survival of small mussels; B, survival of large mussels; C, growth of small mussels; D, growth of large mussels. Colored numbers and error bars indicate means and standard error for water treatments (4 samples per treatment). Text boxes indicate which responses in sediment treatments were significantly different from reference sediment 1 and (or) which water treatments differed significantly from the control (two-way analysis of variance [ANOVA] with mean comparisons by Fisher’s least significant difference). [−, decrease; +, increase]

52    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District 2002). Although we did not attempt to verify induction of MT analytically, the sequence of responses observed in these published studies was consistent with the sequence of events observed in tests A and B: (1) initial toxicity, (2) reversal/ reduction of toxicity, and (3) spillover toxicity observed in tests A and B. Trends observed in toxicity data from the 2019 tests are presented in tables 13 and 14 in a simplified matrix of color-coded cells that indicate which treatments seem to have experienced toxicity, reversals, and spillover. Mussel Age/Size and Behavior Physiological responses to metal exposure may depend on several aspects of the juvenile mussel and the environment. The sensitivity of mussels to metal toxicity may depend on their size and age, as indicated by the greater growth response of the smaller (test A) mussels in sediment 2, compared to the larger (test B) mussels. This differential sensitivity also is reflected in the thresholds for reversal and restoration of growth effects—about 40–80 μg/L as zinc in test A and 80–160 μg/L as zinc in test B. The results of tests A and B also indicate that behavioral responses to environmental conditions may modify the metal exposure of juvenile mussels. The uniformly low growth of Table 13.  Summary of mussel toxicity trends in 2019 test A: sediment-plus-water toxicity test with small mussels. [Data for this table are shown in figure 15. Sediments: 1, low metals/fine (reference); 2, medium metals/fine; 3, medium metals/coarse; 4, high metals/coarse. Responses: NT (no color), not toxic; T (yellow), toxic; R (green), reversal; TT (orange), spillover (toxicity increase). μg/L as zinc, microgram per liter as zinc] Sediment Water (μg/L as zinc) Survival NT NT NT NT NT NT NT NT NT NT NT NT NT NT NT NT NT NT NT NT NT NT NT NT Growth NT NT NT NT R TT T T T T T T NT NT NT NT R TT T T T R R TT Table 14.  Summary of mussel toxicity trends in 2019 test B: sediment-plus-water toxicity test with large mussels. [Data for this table are shown in figure 15. Sediments: 1, low metals/fine (reference); 2, medium metals/fine; 3, medium metals/coarse; 4, high metals/coarse. Responses: NT (no color), not toxic; T (yellow), toxic; R (green), reversal; TT (orange), spillover (toxicity increase). μg/L as zinc, microgram per liter as zinc] Sediment Water (μg/L as zinc) Survival NT NT NT NT NT NT NT T T T R NT NT NT NT NT NT NT NT NT NT NT NT NT Growth NT NT NT NT R TT NT T T T R TT NT NT NT NT R TT NT NT NT R R TT

6. Summary and Conclusions    53 mussels in sediment 2 in test A included low growth in the water control treatment, which can only be attributed to expo­ sure to metals in sediments. In addition, the lack of any rever­ sal and restoration of toxicity in this treatment indicates that these mussels were not substantially affected by levels of met­ als in overlying waters. One plausible explanation for these results is that mussels in this treatment remained burrowed in sediment 2 throughout the test. Such burrowing behavior is consistent with the consensus that juvenile mussels live and feed within the sediment for the first few years of their life (Yeager and others, 1994; Cope and others, 2008). However, a recent study by Kemble and others (2020) reported deeper burrowing of L. siliquoidea in soft, fine-textured sediment compared to coarse sand, indicating that burrowing behav­ ior may vary depending on habitat conditions. Unlike finetextured sediments 1 and 2, the predominance of coarse sand in sediments 3 and 4 may have limited the ability of small juvenile mussels to burrow to avoid high metal concentrations in overlying water. The resulting exposure to waterborne met­ als could have induced MT-mediated resistance to physiologic detoxification at low levels of metals in overlying water, reversal of toxic effects at intermediate levels of water metals, and increased toxicity at the highest waterborne exposure concentrations. Toxicity of Sediment Versus Overlying Water The strong effect of sediment metals on toxicity made it difficult to isolate the effect of metals in overlying water. Treatments that indicated toxicity at control or low metal concentrations in water (for example, sediments 2 and 4 in test A) provide strong evidence of toxicity of sediment metals. At the other extreme, only the toxicity in the water treatments with the highest metals, with nominal zinc concentrations of 80–160 μg/L as zinc, could be confidently attributed to toxic effects of overlying water. Overall, toxicity in these tests was consistent with sediment toxicity thresholds estimated from the 2017 sediment tests. Composite sediment 2 caused severe reductions in growth of the small (1-mm) mussels in test A at SEM-zinc of 2.0 PEQ and sum-PEQ of 3.3 PEQ. This same sediment was associated with toxic effects on survival and growth of larger (2-mm) mussels in test B. In contrast, no tox­ icity was observed in sediment 3 in either test, despite sediment metal concentrations (zinc-PEQ=1.6; sum-PEQ=2.7) that were close to those in sediment 2. The contrast between the consistent toxicity of low-sand sediment 2 and the minimal toxicity in high-sand sediment 3 indicates that sediment texture affected sediment toxicity indirectly by forcing alternative routes of exposure. Juvenile mussels may have been able to burrow more deeply into fine sediments (Kemble and others, 2020), where they would have been primarily exposed to sediment metals, whereas mussels in coarse sediment may have been unable to borrow and thus had greater exposure to waterborne metals. The combination of coarse sand and high metals in sediment 4 indicated a mixture of these responses, with severe toxicity (presumably because of sediment metals) in the control water and low water treatments, reversal of growth effects at intermediate water treatments, and resumption of growth effects at the highest water treatment. Based on this limited dataset, it seems that coarse sediments may limit burrowing behavior of juvenile mussels—perhaps only the smallest size classes—and leave juvenile mussels more vulnerable to toxicity of overlying water. However, the evidence that waterborne metals may stimulate tolerance to metal exposure indicates that toxicity thresholds estimated from our water-only tests may overestimate the sensitivity of juvenile mussels to metals in overlying water. 6. Summary and Conclusions Species richness and abundance of mussel assemblages were negatively associated with metal exposure—Long-term exposure to metals in sediment, water, or diet has been pro­ posed as a plausible cause for reduced mussel abundance and reduced mussel species richness in streams of the study area (Angelo and others, 2007). In the present study, analyses of metals in tissues of Corbicula collected during mussel surveys demonstrated that metals were biologically available for uptake by bivalve mollusks (fig. 6) and that metal concentrations in sediment and water were negatively associated with mussel community metrics (fig. 16). In the qualitative mussel survey, CPUE (mussels/time) and species richness differed signifi­ cantly between sites with low versus high metals in water or sediment (p<0.02). In the quantitative survey, species richness and density (mussels/area) also differed significantly between low- and high-metal sediments (p<0.03), but associations between these metrics and metal concentrations in water were less significant (p<0.10). Although mussel community met­ rics did not generally follow the clear concentration-response trends typical of controlled laboratory tests, species richness (expressed as a percentage of maximum species detected in each stream reach) was reduced by 50 percent at sites with SEM-zinc concentrations of 0.65 PEQ (298 µg/g) or greater (fig. 17). Water-only toxicity tests and whole-sediment toxicity tests documented significant toxic effects on survival, growth, and biomass of juvenile mussels.—Effect concentrations for juve­ nile mussels exposed to waterborne zinc in 12-week tests were consistent with results of other recent chronic zinc tests with L. siliquoidea (Wang and others, 2020). In tests with water­ borne zinc only, toxic effects were limited to reduced growth and biomass, and thresholds for these endpoints were greater (indicating lesser sensitivity) after 12 weeks than after 4 weeks. Despite this reduced toxicity with increasing age/size, chronic effects of zinc only on juvenile mussels in these 12-week tests were detected at filtered zinc concentrations far less than the hardness-based water quality criteria for zinc, for example, 180 µg/L at the hardness of the water tested (165 mg/L as cal­ cium carbonate; EPA, 2021b).

54    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District A B r 2=0.29 (p=0.015) r 2=0.38 (p=0.005) Density of mussels per square meter Live mussel species Mussel biomass, day 83, in milligrams Mussel survival, day 83 (of 10) EXPLANATION Reference site Test site Figure 16.  Relations between mussel endpoints from 2017 sediment toxicity tests and community metrics from quantitative mussel surveys, with linear regressions. A, mussel biomass versus density; B, mussel survival versus species richness. Hollow symbols are reference sites; solid symbols are test sites. Solid line is linear regression. [r 2, coefficient of determination; p, probability]

6. Summary and Conclusions    55 EC50=0.65 PEQ r 2=0.41 SEM-zinc (PEQ) Species richness, in percent of maximum by stream EXPLANATION Spring River Center Creek Shoal Creek Neosho River Figure 17.  Estimated 50-percent effect concentration (EC50) for effects of sediment zinc on mussel species richness. Species richness values are expressed as a percentage of the maximum number of species for each stream. Curved line is three-parameter sigmoid regression. [PEQ, probable effect quotient, r 2, coefficient of determination; SEM, simultaneously extracted metals] Mussels exposed to mixtures of zinc, lead, and cadmium in water for 12 weeks not only indicated greater sensitivity of effects on growth and biomass but also indicated significant effects on survival that were not observed after 4 weeks. The range of water-only EC20s for effects of the metal mixture on growth and biomass of juvenile mussels (6.4–20 µg/L as zinc) was similar to waterborne zinc concentrations in stream water (>10 μg/L) that were associated with decreases in abundance and species richness in mussel assemblages (fig. 16). Whole-sediment toxicity tests completed in fall 2017 also demonstrated significant toxic effects on juvenile mus­ sels. After 12 weeks, mussels in most of the 19 test sediments had mean survival, growth, and biomass the ranges of means for reference sediments, and two or more test sediments had means for each endpoint that were significantly less than all reference means (table 4). Toxicity thresholds for juvenile mussels estimated from several metrics of metal exposure (SEM-zinc, sum-PEQ for the zinc/lead/cadmium mixture, [ΣSEM–AVS], and pore-water zinc) in the 12-week sediment test were less than thresholds from 4-week mussel tests with sediments sampled in 2007 (MacDonald and others, 2009; table 10). For example, the GM–EC20 for SEM-zinc was 2.49 PEQ or 1,143 µg/g, which is much lower than the T20 of 23,700 µg/g determined for sediment zinc (table 10). In con­ trast, toxicity thresholds for amphipods (H. azteca) in 6-week tests in 2017 were within a factor of two of thresholds from 4-week tests in 2007. Concurrent exposure of juvenile mussels to metals in water and sediment did not always cause increased toxic effects.—The first sediment-plus-water test (spring 2018) indi­ cated little toxicity, even though this test included exposure to sediments that were toxic in the 2017 whole-sediment test plus a water-phase exposure with a nominal zinc concentration of 10 μg/L (in a mixture with lead and cadmium) that was consis­ tent with EC20s derived from the water-only tests (fig. 8). The reduced toxicity observed in the 2018 test was contrary to our expectations that addition of metals to overlying water would increase toxicity, relative to sediments without added metals in overlying water. Several factors could help explain the lack of toxicity in the 2018 water-plus-sediment test: (1) the mean starting size of mussels in the 2018 test was greater than the 2017 test (1.4 mm versus 1.2 mm; these older, larger mus­ sels could have been more tolerant to the metal exposure); (2) PWP-zinc concentrations averaged lower in the 2018 test, consistent with a possible decrease in metal bioavailability; and (3) addition of low levels of waterborne metals in the 2018 test could have induced metal tolerance as we observed in the 2019 study. Results of the 2019 water-plus-sediment test indicate that mussels developed tolerance to metal toxicity.—This test used composite sediment samples, which were selected to produce four sediment treatments with combinations of low and high metal concentrations and low and high coarse sand content. The water-phase exposure was delivered by a proportional

56    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District diluter to produce ranges of the three-metal mixture of as much as 80 μg/L as zinc (test A) and 160 μg/L as zinc (test B). Results of this test indicated that significant toxic effects at low waterborne metal concentrations (as much as 20–40 μg/L) were reversed at higher waterborne metal concentrations (tables 13 and 14). We hypothesize that these reversals of toxicity at zinc concentrations greater than water-only toxicity thresholds (20–80 μg/L as zinc) were attributable to induction of metal-binding protein(s), such as MT (Amiard and oth­ ers, 2006). The pattern of toxicity we observed—toxicity at low waterborne metal concentrations, reduced toxicity with increasing metal concentrations, and a return to toxicity at the highest concentrations tested—is consistent with the pattern of MT effects in published studies (for example, Correia and others, 2002). Burrowing behavior of juvenile mussels may have affected their development of tolerance to metal toxicity.—It is widely accepted that juvenile mussels in the wild remain buried in fine-textured sediments (Yeager and others, 1994). In 2019 tests with sediments dominated by coarse sand (com­ posites 3 and 4), toxic effects corresponded more closely to metal concentrations in overlying water than to metals in sediment, perhaps because juvenile mussels could not burrow into these coarse sediments to avoid metals in overlying water (Kemble and others, 2020). Reversals of toxicity were most evident in sediment 2 (tables 11, 12, 13, and 14), which had metal concentrations close to sediment effects thresholds, and at low added zinc concentrations in water. At higher added metal concentrations in water, we hypothesize that toxicity was reduced because of induction of MT, but this protection was overwhelmed at the highest added metal concentrations (typically >40 μg/L). Toxic effects on mussels in the 2017 study were detected at sediment metal concentrations less than thresholds for toxicity in 2007 tests.—EC20s estimated from the 2017 sedi­ ment test were not as robust as the T20s derived by Mac­ Donald and others (2009) because the previous thresholds were derived from larger sample sets (70 sediments tested with amphipods and 43 with mussels, compared to 25 sedi­ ments tested in 2017) and at wider ranges of sediment metal concentrations. Maximum zinc, lead, and cadmium concentra­ tions were 3–5 times greater in 2007 samples (Ingersoll and others, 2008). To reduce the influence of individual EC20s estimated from the 2017 test, we calculated composite EC20s for mussels and amphipods based on GM–EC20s for multiple endpoints. The zinc GM–EC20s for mussels and amphipods were 2.8 PEQ and 3.4 PEQ, respectively. These 2017 results are lower than the T20s for both species estimated by Mac­ Donald and others (2009), consistent with longer exposures (12 weeks for mussels, 6 weeks for amphipods) in the 2017 tests. The GM–EC20s from 2017 tests were less than the T20s by factors from 5 to 23 for mussels and by factors of about 0.7 to 2.4 for amphipods. The GM–EC20s from the 2017 mussel test (for example, the EC20 for SEM-zinc of 2.5 PEQ (or 1,148 μg/g) were more protective than either 2007 mus­ sel thresholds (T20 for SEM-zinc=6.4 PEQ; MacDonald and others, 2009) or 2017 amphipod thresholds (EC20 for SEMzinc=4.0 PEQ). Thresholds based on other metrics of metal exposure, including sum-PEQ and [ΣSEM–AVS], indicated similar trends between species and years, and GM–EC20s for these variables from 2017 tests were less than corresponding T20 values (table 10). Adverse effects on mussel community metrics and toxic effects in water-only toxicity tests were detected at waterborne metal concentrations typical of Tri-State Mining District streams.—Waterborne zinc concentrations seemed to explain the similarity of toxic effects on juvenile mussels across dif­ ferent exposure types. Similar pore water (push-point) zinc concentrations were associated with thresholds for effects in the mussel survey and in the water-only toxicity tests. The zinc EC20s for the water-only (metal-mixture) test ranged from 6.4 μg/L (growth and biomass) to 48 μg/L (survival; table 6). This range is consistent with the significant reductions in com­ munity metrics at sites with zinc in surface water and pushpoint pore water greater than 10 μg/L (fig. 7). In addition, sediment 2 in the 2019 water-plus-sediment tests (tests A and B) indicated toxic effects on survival and growth associated with overlying water at zinc concentrations in a similar range, 14–43 μg/L (tables 11 and 12). Both of these results indicate that waterborne zinc concentrations in the range of 6.3 to 43 μg/L, co-occurring with proportional levels of lead and Cd, were associated with effects on juvenile mussels. However, toxicity thresholds based on Zn in peeper samples (EC20s from 86 to 215 μg/L; table 10) were greater than those for the water-only tests. Toxic effects in the 2019 water-plus-sediment test (sediment 2) were also associated with higher Zn con­ centrations in peeper samples (58–200 μg/L). These greater thresholds for peeper metals may indicate lower bioavailabil­ ity of pore water in contact with deep sediment, or they may indicate that mussels inhabiting surface sediments are exposed to lower metal concentrations, for example, concentrations closer to EC20s for overlying water or water-only tests. Toxic effects on mussels in sediment toxicity tests were consistent with effects on mussel community metrics.— Survival and biomass of mussels in the 2017 sediment toxicity test were significantly correlated with density and species richness of mussel assemblages (fig. 16). These associations indicate that reductions in survival and (or) growth could contribute to reduced abundance of mussels in the field. These associations also are consistent with a recent study of mussel assemblages in defaunated Tennessee streams, which con­ cluded that long-term reduction in mussel growth was associ­ ated with failure of juvenile recruitment and loss of mussel species (Haag and others, 2019). Community metrics indicated significant differences between high- and low-metal sites (fig. 7), but these metrics were more variable than endpoints of laboratory tests and they generally did not produce good fits to CRMs. However, we used a sigmoidal model to estimate that a 50-percent effect concentration for reduction in mussel spe­ cies richness—normalized to the maximum richness for each stream segment—was 0.64 PEQ (or 294 µg/g as SEM-zinc (fig. 17).

6. Summary and Conclusions    57 Table 15.  Frequency of mussel community effects and toxic effects in laboratory tests. [qual., qualitative; CPUE, catch per unit effort; quant., quantitative; X, less than reference; O, not less than reference; --, no data] Site Sites with responses less than reference Percentage of sites with effect Qual. species Qual. CPUE Quant. species Quant. density Quant. survey Survey mean Toxicity Toxicity mean Center Creek CC01G2 O CC05D CC10B Middle Spring River SP05 O O SP06 O O SP08A O SP08B SP10C O Shoal Creek SH02A SH05A O O SH06A SH07 SH08 O O O SH10B SH15A SH16A SH16B ─ Lower Spring River SP12B SP12C O O SP13 O SP19 O O O O SP20 ─ O O SP22 O Neosho River NR03 NR04 NR05 NR06 O O O

58    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District Toxicity tests and mussel surveys identified similar stream reaches with effects on mussel communities.—Rankings of severity of effects based on toxicity tests and mussel surveys were generally in agreement, but some responses differed in severity. Results of the qualitative mussel surveys indicated severe effects relative to reference conditions across almost all test sites. Richness and CPUE for 78 percent of test sites had responses less than the 5th percentile of the reference enve­ lope (table 15). Results of quantitative surveys and sediment toxicity tests had wider ranges of sensitivity, and effects were detected at 35–58 percent of sites for the quantitative sur­ vey and at 4–48 percent of sites for sediment toxicity. These graded responses made quantitative survey data and toxicity data more useful for distinguishing different levels of effects among sites (table 15). The quantitative survey and sedi­ ment toxicity tests indicated that the greatest effects of metals toxicity were detected in Center Creek, where 67 percent of test sediments were toxic and 88 percent of sites had reduced quantitative survey metrics. Upper Shoal Creek indicated little or no toxicity and no effects in either quantitative surveys or toxicity tests, despite two sites that exceeded sediment toxicity thresholds (SH02 and SH05); however, toxic effects were detected at several sites in the downstream reach (SH07 through SH16B). The middle Spring River (downstream from Center Creek; sites SP08 and SP10) indicated no toxicity but indicated effects in mussel surveys. Several sites in lower Spring River below Empire Lake (SP12B, SP13, and SP22) had a high frequency of toxicity and moderate effects in the quantitative survey. We did not have enough colocated survey and toxicity data to evaluate the severity of effects in the Neosho River drainage. References Cited Amiard, J., Amiardtriquet, C., Barka, S., Pellerin, J., and Rainbow, P., 2006, Metallothioneins in aquatic invertebrates—Their role in metal detoxification and their use as biomarkers: Aquatic Toxicology (Amsterdam, Netherlands), v. 76, no. 2, p. 160–202. [Also available at ://doi.org/​10.1016/​j​.aquatox.2​005.08.015.] Angelo, R.T., Cringan, M.S., Chamberlain, D.L., Stahl, A.J., Haslouer, S.G., and Goodrich, C.A., 2007, Residual effects of lead and zinc mining on freshwater mussels in the Spring River Basin (Kansas, Missouri, and Oklahoma, USA): Science of the Total Environment, v. 384, no. 1–3, p. 467–496. [Also available at ://doi.org/​10.1016/​j.s​ citotenv.2​007.05.045.] Ankley, G.T., Di Toro, D.M., Hansen, D.J., and Berry, W.J., 1996, Technical basis and proposal for deriving sediment quality criteria for metals: Environmental Toxicology and Chemistry, v. 15, no. 12, p. 2056–2066. [Also available at ://doi.org/​10.1002/​etc.5620151202.] ASTM International, 2015, Standard guide for conducting lab­ oratory toxicity tests with freshwater mussels, E2455–06: West Conshohocken, Pa., ASTM International, Annual Book of ASTM Standards, v. 11.06, p. 1378–1429. ASTM International, 2020, Test method for measuring the toxicity of sediment-associated contaminants with fresh­ water invertebrates, E1706–05: West Conshohocken, Pa., ASTM International, Annual Book of ASTM Standards, v. 11.06, 118 p. Barnhart, M.C., 2006, Buckets of muckets—A compact sys­ tem for rearing juvenile freshwater mussels: Aquaculture (Amsterdam, Netherlands), v. 254, no. 1–4, p. 227–233. [Also available at ://doi.org/​10.1016/​j.aqu​aculture.2​ 005.08.028.] Baudrimont, M., Andrès, S., Metivaud, J., Lapaquellerie, Y., Ribeyre, F., Maillet, N., Latouche, C., and Boudou, A., 1999, Field transplantation of the freshwater bivalve Corbicula fluminea along a polymetallic contamina­ tion gradient (river Lot, France)—II. Metallothionein response to metal exposure: Environmental Toxicology and Chemistry, v. 18, no. 11, p. 2472–2477. [Also available at ://doi.org/​10.1002/​etc.5620181113.] Besser, J.M., Brumbaugh, W.G., May, T.W., and Schmitt, C.J., 2007a, Biomonitoring of lead, zinc, and cadmium in streams draining lead-mining and non-mining areas, southeast Missouri, USA: Environmental Monitoring and Assessment, v. 129, no. 1–3, p. 227–241. [Also available at ://doi.org/​10.1007/​s10661-​006-​9356-​9.] Besser, J.M., Ingersoll, C.G., Brumbaugh, W.G., Kemble, N.E., May, T.W., Wang, N., MacDonald, D.D., and Roberts, A.D., 2015, Toxicity of sediments from lead-zinc mining areas to juvenile freshwater mussels (Lampsilis siliquoi­ dea) compared to standard test organisms: Environmental Toxicology and Chemistry, v. 34, no. 3, p. 626–639. [Also available at ://doi.org/​10.1002/​etc.2849.] Besser, J.M., Mebane, C.A., Mount, D.R., Ivey, C.D., Kunz, J.L., Greer, I.E., May, T.W., and Ingersoll, C.G., 2007b, Sensitivity of mottled sculpins (Cottus bairdi) and rainbow trout (Onchorhynchus mykiss) to acute and chronic toxicity of cadmium, copper, and zinc: Environmental Toxicology and Chemistry, v. 26, no. 8, p. 1657–1665. [Also available at ://doi.org/​10.1897/​06-​571R.1.] Bonnail, E., Sarmiento, A., DelValls, T.A., Nieto, and J.M., Riba, 2016, Assessment of metal contamination, bioavail­ ability, toxicity and bioaccumulation in extreme metal­ lic environments (Iberian Pyrite Belt) using Corbicula fluminea: Science of the Total Environment, v. 544, p. 1031–1044. [Also available at ://doi.org/​10.1016/​j.s​ citotenv.2​015.11.131.]

References Cited    59 Brumbaugh, W.G., Hammerschmidt, C.R., Zanella, L., Rogevich, E., Salata, G., and Bolek, R., 2011, Interlaboratory comparison of measurements of acidvolatile sulfide and simultaneously extracted nickel in spiked sediments: Environmental Toxicology and Chemistry, v. 30, no. 6, p. 1306–1309. [Also available at ://doi.org/​10.1002/​etc.506.] Brumbaugh, W.G., May, T.W., Besser, J.M., Allert, A.L., and Schmitt, C.J., 2007, Assessment of elemental concentrations in streams of the New Lead Belt in southeastern Missouri, 2002–05: U.S. Geological Survey Scientific Investigations Report 2007–5057, 60 p. [Also available at ://doi.org/​ 10.3133/​sir20075057.] Cleveland, D., Brumbaugh, W.G., and MacDonald, D.D., 2017, A comparison of four porewater sampling methods for metal mixtures and dissolved organic carbon and the impli­ cations for sediment toxicity evaluations: Environmental Toxicology and Chemistry, v. 36, no. 11, p. 2906–2915. [Also available at ://doi.org/​10.1002/​etc.3884.] Cope, W.G., Bringolf, R.B., Buchwalter, D.B., Newton, T.J., Ingersoll, C.G., Wang, N., Augspurger, T., Dwyer, F.J., Barnhart, M.C., Neves, R.J., and Hammer, E., 2008, Differential exposure, duration, and sensitivity of unionoi­ dean bivalve life stages to environmental contaminants: Journal of the North American Benthological Society, v. 27, no. 2, p. 451–462. [Also available at ://doi.org/​ 10.1899/​07-​094.1.] Correia, A.D., Lima, G., Costa, M.H., and Livingstone, D.R., 2002, Studies on biomarkers of copper exposure and toxicity in the marine amphipod Gammarus locusta (Crustacea)—I. Induction of metallothionein and lipid peroxidation: Biomarkers, v. 7, no. 5, p. 422–437. [Also available at ://doi.org/​10.1080/​135475002760413516.] Couillard, Y., Campbell, P.G.C., and Tessier, A., 1993, Response of metallothionein concentrations in a freshwa­ ter bivalve (Anodonta grandis) along an environmental cadmium gradient: Limnology and Oceanography, v. 38, no. 2, p. 299–313. [Also available at ://doi.org/​10.4319/​ lo.1993.38.2.0299.] EcoAnalysts, Inc., 2018, Final report—Tri-State Mining District unionid assessment, Missouri, Kansas, and Oklahoma, 2016–2018: Columbia, Mo., U.S. Fish and Wildlife Service, prepared by EcoAnalysts, Inc., 101 p. [Also available at :​//www.cerc​.usgs.gov/​orda_​docs/​ DocHandler.ashx?​task=​get&ID=​8373.] Erickson, R.J., 2010, Toxicity Relationship Analysis Program (TRAP), version 1.21: Washington, D.C., U.S. Environmental Protection Agency, EPA/600/C–11/002, 36 p. Gutiérrez, M., Qiu, X., Collette, Z.J., and Lurvey, Z.T., 2020, Metal content of stream sediments as a tool to assess reme­ diation in an area recovering from historic mining contami­ nation: Minerals (Basel), v. 10, no. 3, art. 247, 12 p. [Also available at ://doi.org/​10.3390/​min10030247.] Haag, W.R., Culp, J.J., McGregor, M.A., Bringolf, R., and Stoeckel, J.A., 2019, Growth and survival of juvenile freshwater mussels in streams—Implications for under­ standing enigmatic mussel declines: Freshwater Science, v. 38, no. 4, p. 753–770. [Also available at ://doi.org/​ 10.1086/​705919.] Haag, W.R., and Williams, J.D., 2014, Biodiversity on the brink—An assessment of conservation strategies for North American freshwater mussels: Hydrobiologia, v. 735, no. 1, p. 45–60. [Also available at ://doi.org/​10.1007/​s10750-​ 013-​1524-​7.] Hunt, J.W., Anderson, B.S., Phillips, B.M., Newman, J., Tjeerdema, R.S., Fairey, R., Puckett, H.M., Stephenson, M., Smith, R.W., Wilson, C.J., and Taberski, K.M., 2001, Evaluation and use of sediment toxicity reference sites for statistical comparisons in regional assessments: Environmental Toxicology and Chemistry, v. 20, no. 6, p. 1266–1275. [Also available at ://doi.org/​10.1002/​ etc.5620200615.] Ingersoll, C.G., Besser, J.M., Brumbaugh, W.G., Ivey, C.D., Kemble, N.E., Kunz, J.L., May, T.W., Wang, N., MacDon­ ald, D.D., and Smorong, D.E., 2008, Sediment chemistry, toxicity, and bioaccumulation data report for the U.S. Environmental Protection Agency—Department of the Interior sampling of metal-contaminated sediment in the Tri-State Mining District in Missouri, Oklahoma, and Kansas: Columbia, Mo., and Nanaimo, British Columbia, U.S. Geological Survey and MacDonald Environmental Sciences, 287 p. Ivey, C.D., and Besser, J.M., 2022, Response of juvenile mus­ sels and amphipods to metal concentrations in water and sediment of streams draining the Tri-State Mining District, Missouri, Kansas and Oklahoma, USA: U.S. Geological Survey data release, ://doi.org/​10.5066/​P9R9FHAF. Johnson, A.W., Gutiérrez, M., Gouzie, D., and McAliley, L.R., 2016, State of remediation and metal toxicity in the Tri-State Mining District, USA: Chemosphere, v. 144, p. 1132–1141. [Also available at ://doi.org/​10.1016/​j.che​ mosphere.2​015.09.080.] Juracek, K.E., and Drake, K.D., 2016, Mining-related sedi­ ment and soil contamination in a large Superfund site— Characterization, habitat implications, and remediation: Environmental Management, v. 58, no. 4, p. 721–740. [Also available at ://doi.org/​10.1007/​s00267-​016-​0729-​8.]

60    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District Kemble, N.E., Besser, J.M., Steevens, J., and Hughes, J.P., 2020, Assessment of burrowing behavior of freshwater juvenile mussels in sediment: Freshwater Mollusk Biology and Conservation, v. 23, no. 2, p. 69–81. [Also available at ://doi.org/​10.31931/​.v23i2.2020.69-​81.] MacDonald, D.D., Ingersoll, C.G., and Berger, T.A., 2000, Development and evaluation of consensus-based sediment quality guidelines for freshwater ecosystems: Archives of Environmental Contamination and Toxicology, v. 39, no. 1, p. 20–31. [Also available at ://doi.org/​10.1007/​ s002440010075.] MacDonald, D.D., Smorong, D.E., Ingersoll, C.G., Besser, J.M., Brumbaugh, W.G., Kemble, N.E., May, T.W., Ivey, C.D., Irving, S., and O’Hare, M., 2009, Development and evaluation of sediment and pore-water toxicity thresholds to support sediment quality assessments in the Tri-State Mining District (TSMD), Missouri, Oklahoma, and Kansas: Nanaimo, British Columbia, MacDonald Environmental Sciences Ltd., prepared for U.S. Environmental Protection Agency, [variously paged]. [Also available at ​ s://archiv​e.epa.gov/​region07/​cleanup/​npl-​archive/​web/​pdf/​ 000030284900.pdf.] Mebane, C.A., Hennessy, D.P., and Dillon, F.S., 2007, Developing acute-to-chronic toxicity ratios for lead, cadmium, and zinc using rainbow trout, a mayfly, and a midge: Water, Air, and Soil Pollution, v. 188, no. 1–4, p. 41–66. [Also available at ://doi.org/​10.1007/​s11270-​ 007-​9524-​8.] Natural Resources Conservation Service, 2014, Kellogg soil survey laboratory methods manual: Lincoln, Nebr., U.S. Department of Agriculture, Natural Resources Conservation Service Soil Survey Investigations Report no. 42, v. 5.0, 1,001 p. Roberts, A.D., Hundley, J., Rosenberger, A., Bouska, K., Mosby, D., Simmons, B., and Lindner, G., 2016, Draft report—Quantitative survey of freshwater mussels (Unionidae) and assessment of sediment contamina­ tion in the Big River, Missouri: Columbia, Mo., U.S. Fish and Wildlife Service, 63 p. [Also available at ://www.fws.gov/​sites/​default/​files/​documents/​nrdar-​ final-​study%20report-​SEMO%​20Big%20Ri​ver%20muss​ el%20and%2​0sediment-​2016.pdf.] Steevens, J.A., Besser, J.M., Dorman, R.A., and Sparks, D.W., 2020, Influence of remediation on sediment toxicity within the Grand Calumet River, Indiana, USA: Chemosphere, v. 249, art. 126056, 11 p. [Also available at ://doi.org/​ 10.1016/​j.che​mosphere.2​020.126056.] U.S. Environmental Protection Agency [EPA], 2000, Methods for measuring the toxicity and bioaccumulation of sediment-associated contaminants with freshwater inver­ tebrates: Duluth, Minn., U.S. Environmental Protection Agency, EPA/600/R–99/064, 133 p. U.S. Environmental Protection Agency [EPA], 2001, Methods for collection, storage and manipulation of sediments for chemical and toxicological analyses—Technical manual: Washington, D.C., U.S. Environmental Protection Agency, EPA–823–B–01–002, [variously paged]. U.S. Environmental Protection Agency [EPA], 2007, Method 6200—Field portable X-ray fluorescence spectrometry for the determination of elemental concentrations in soil and sediment (revision 0): U.S. Environmental Protection Agency, SW–846, 32 p. U.S. Environmental Protection Agency [EPA], 2021a, Abandoned mine lands: U.S. Environmental Protection Agency web page, accessed November 8, 2021, at ://www.epa.gov/​superfund/​abandoned-​mine-​lands. U.S. Environmental Protection Agency [EPA], 2021b, National recommended water quality criteria—Aquatic life criteria table: EPA web page, accessed November 30, 2021, at ://www.epa.gov/​wqc/​national-​recommended-​water-​ quality-​criteria-​aquatic-​life-​criteria-​table. Wang, N., Ingersoll, C.G., Hardesty, D.K., Ivey, C.D., Kunz, J.L., May, T.W., Dwyer, F.J., Roberts, A.D., Augspurger, T., Kane, C.M., Neves, R.J., and Barnhart, M.C., 2007, Acute toxicity of copper, ammonia, and chlorine to glochidia and juveniles of freshwater mussels (Unionidae): Environmental Toxicology and Chemistry, v. 26, no. 10, p. 2036–2047. [Also available at ://doi.org/​10.1897/​06-​523R.1.] Wang, N., Ingersoll, C.G., Ivey, C.D., Hardesty, D.K., May, T.W., Augspurger, T., Roberts, A.D., van Genderen, E., and Barnhart, M.C., 2010, Sensitivity of early life stages of freshwater mussels (Unionidae) to acute and chronic toxicity of lead, cadmium, and zinc in water: Environmental Toxicology and Chemistry, v. 29, no. 9, p. 2053–2063. [Also available at ://doi.org/​10.1002/​etc.250.] Wang, N., Ingersoll, C.G., Kunz, J.L., Brumbaugh, W.G., Kane, C.M., Evans, R.B., Alexander, S., Walker, C., and Bakaletz, S., 2013, Toxicity of sediments potentially contaminated by coal mining and natural gas extraction to unionid mussels and commonly tested benthic invertebrates: Environmental Toxicology and Chemistry, v. 32, no. 1, p. 207–221. [Also available at ://doi.org/​10.1002/​ etc.2032.]

References Cited    61 Wang, N., Kunz, J.L., Cleveland, D.M., Steevens, J.A., Hammer, E.J., Van Genderen, E., Ryan, A.C., and Schlekat, C.E., 2020, Evaluation of acute and chronic toxicity of nickel and zinc to 2 sensitive freshwater benthic inver­ tebrates using refined testing methods: Environmental Toxicology and Chemistry, v. 39, no. 11, p. 2256–2268. [Also available at ://doi.org/​10.1002/​etc.4841.] Williams, J.D., Warren, M.L., Jr., Cummings, K.S., Harris, J.L., and Neves, R.J., 1993, Conservation status of fresh­ water mussels of the United States and Canada: Fisheries (Bethesda, Md.), v. 18, no. 9, p. 6–22. [Also available at ://doi.org/​10.1577/​1548-​8446(19​93)018%3C0​ 006:CSOFMO​%3E2.0.CO;​2.] Yeager, M.M., Cherry, D.S., and Neves, R.J., 1994, Feeding and burrowing behaviors of juvenile rainbow mussels, Villosa iris (Bivalvia:Unionidae): Journal of the North American Benthological Society, v. 13, no. 2, p. 217–222. [Also available at ://doi.org/​10.2307/​1467240.] Zimmerman, M.J., Massey, A.M., and Campo, K.W., 2005, Pushpoint sampling for defining spatial and temporal variations in contaminant concentrations in sediment pore water near the ground-water/surface-water interface: U.S. Geological Survey Scientific Investigations Report 2005–5036, 70 p. [Also available at ://doi.org/​10.3133/​ sir20055036.]

Appendixes 1–3

64    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District Appendix 1.  Methods for Mussel Surveys Reconnaissance Site Selection Reconnaissance sites were selected based on the follow­ ing criteria: accessibility, published sediment chemistry and toxicity data, available unionid community data, and geo­ graphic information system data. Geographic information sys­ tem data were used to (1) delineate historically stable stream reaches; (2) complete a broad-scale geographic overview for identification of habitat features generally associated with freshwater mussel communities such as gravel bars, riffles, and shoals; and (3) identify access points. Because unionids generally require permanent, flowing water above stable, gravel-dominated substrates with a compo­ nent of finer grained particles such as sand, sites were chosen with these criteria in mind and pool habitats were avoided. Areas where the overlay of 1963 topographic maps and cur­ rent maps indicated shifts in the channel also were avoided because unionid communities are generally detected in stable river channels. Reconnaissance sites were chosen to represent the range of potential contamination (that is, reference, low contamina­ tion, medium contamination, and high contamination) and were distributed as evenly as possible across States, basins, and streams while meeting other site selection criteria. Sites also were placed upstream and downstream from mining input streams to document any change in sediment/water contami­ nation and unionid community characteristics. As parts of the Spring River, Center Creek, Shoal Creek, and Neosho River were traversed, additional sites were investigated based on field observations such as the presence of shell material on the banks or in the shallows and (or) seemingly suitable unionid habitat (stable but permeable substrate with some current velocity). Because these conditions vary with the river system, reference reaches were investigated first to facilitate understanding of mussel habitat within the Spring and Neosho River systems. An initial 56 sites were selected in 2016, with additional sites added during the 2016 reconnaissance survey. Additional reconnaissance surveys were completed in 2017 and 2018 to fill in geographic gaps and to further investigate river reaches of interest. A total of 102 sites were investigated: 2 on the North Fork Spring River, 34 on the Spring River, 27 on Center Creek, 30 on Shoal Creek, 1 on Lost Creek, 2 on Tar Creek, and 6 on the Neosho River (EcoAnalysts, Inc., 2018). 2016–17 Reconnaissance Survey Accessibility varied by site and by State and included streamside, boat, and canoe access. Long stretches of river were accessible in Missouri, whereas only point access was available in much of Kansas and Oklahoma. In Missouri, most sites were accessed through bridge right-of-ways or other public lands, and sites upstream and downstream from the access were investigated while walking in the stream or from a canoe or boat. Because Kansas and Oklahoma required land­ owner permission to access and sample sites, only sites with landowner access permission to the river were investigated. Permissions were coordinated before fieldwork when possible and in the field when needed. A unionid habitat checklist was developed by Eco­ Analysts, Inc.; the U.S. Fish and Wildlife Service; and the U.S. Geological Survey Columbia Environmental Research Center to incorporate unionid habitat factors that might be qualitatively assessed. Unionid habitat requirements can vary by species and by stream, and efforts to describe suitable unionid habitat have ranged from simple variables such as substrate composition and current velocity to more complex hydraulic parameters such as shear stress, boundary Reynold’s number, and relative substrate stability (Ecological Specialists, Inc., 2014). For the reconnaissance survey, simple attributes that could be qualitatively assessed were used. Habitat attri­ butes that allow unionids to persist include current veloci­ ties slow enough to allow juveniles to settle to the substrate; substrate that is sufficiently stable so that it is not scoured during high flow events; available food, dissolved oxygen, and minerals; favorable water temperatures and water quality; and habitat for host fish (Strayer, 2008; Haag, 2012). Because the Asian clam, Corbicula fluminea (O.F. Müller, 1774), has been used as a unionid surrogate in the past (Angelo and others, 2007), the presence of C. fluminea also was documented, indi­ cating that sediment and water quality were at least sufficient for this tolerant species. A map was created at each site, and habitat observations recorded included the following information: Geomorphic stream features (riffle, pool, run), Bank stability, Stream stability, Substrate composition (Wentworth scale; Wentworth, 1922), Substrate stability/embeddedness/burrowing capacity, Instream depositional features, and Distance to tributaries.

Appendix 1.  Methods for Mussel Surveys    65 Each site was subjected to an initial reconnaissance search consisting of field team members dispersing throughout the area searching for the presence of unionids, C. fluminea, and (or) suitable unionid habitat. Methods included a shore­ line search for shell material, snorkeling throughout the area and fanning the substrate, and scuba diving in deeper areas. If a unionid community was discovered and (or) field biolo­ gists agreed that the site harbored suitable unionid habitat, the area was delineated with a survey-grade Trimble Global Positioning System and was qualitatively sampled with timed searches to determine unionid distribution throughout the area, unionid community composition, and catch per unit effort (qualitative sites). Qualitative search effort per site ranged from 30 to 360 person-minutes and averaged 76 minutes. All unionid data were recorded on standard data sheets, and habitat characterizations were logged within field notes and within the habitat checklist. All data sheets were scanned upon return from the field, and data were entered into a database and underwent a line-by-line quality assurance/qual­ ity control procedure. All Global Positioning System points and features were downloaded and placed into an ArcGIS catalog. Digital photographs taken in the field were stored on a project-specific memory card with a unique identifier, and photographs were compiled into an archive set. This archive is unedited and will not be opened. The project-specific memory card was tracked through a chain of custody form and is cur­ rently housed at EcoAnalysts, Inc. 2017 Quantitative Survey After the reconnaissance survey, qualitative unionid data, habitat characteristics, and results of sediment and pore-water sampling were analyzed for unionid community character­ istics (number of live species, total number of species, catch per unit effort); habitat suitability; and level of lead, zinc, and cadmium contamination. Using this information, 22 sites were selected for quantitative unionid sampling (EcoAnalysts, Inc., 2018). Sites were distributed over a range of contamination levels on the North Fork Spring River, Spring River, Center Creek, Shoal Creek, and Neosho River. At each site, as many as 100 randomly distributed 0.25-square meter quantitative quadrat samples were collected within the previously delin­ eated unionid habitat area/unionid community to characterize species richness, density, and age distribution. The number of quantitative samples depended on the size of the delineated area and ranged from 15 to 100 (EcoAnalysts, Inc., 2018). Each quantitative sample was excavated to a depth of about 15 centimeters, and material was placed into an attached mesh bag. Each sample was washed through a 6-millimeter sieve and searched for unionids and C. fluminea. All live unionids were identified to species, measured (length in millimeters), and aged (external annuli count). Fresh dead shells were iden­ tified to species and counted. Weathered dead and subfossil shells were identified to species and noted as present. All live unionids were returned to the river near their collection point, and at least one individual of each species was photographed. Live C. fluminea detected during quantitative sampling were enumerated and retained to obtain soft tissue as unionid sur­ rogates for determination of tissue metals concentration. Live C. fluminea were placed in a plastic sample bag (one com­ posite sample per site) filled with water and allowed to expel gut contents over a 24-hour period, after which they were transferred to a clean plastic jar, placed in a cooler filled with dry ice, and frozen. All C. fluminea samples were transferred to the U.S. Geological Survey Columbia Environmental Research Center after field sampling. The following quantitative metrics were used to compare sites within each stream: Total number of live unionids, Total number of live species (species richness), Total number of species (live and dead), Mean total unionid density (per square meter) and standard error (SE), Mean adult unionid density (per square meter) and SE, Mean juvenile unionid density (per square meter) and SE, Percentage of juvenile unionids, Percentage of unionid mortality, and Mean C. fluminea density (per square meter) and SE. Quantitative microscale habitat features (depth [in cen­ timeters], current velocity [in feet per second], and substrate constituents) were recorded at each quantitative sample point to compare physical habitat characteristics between sites in each stream. For each quantitative sample, a Gravelometer was used to characterize particle sizes (in millimeters) in each of the four corners and the middle of the quadrat frame (modified pebble count). Pebble count data were entered into a program (Reference Reach Spreadsheet, version 4.2, Ohio Department of Natural Resources) to calculate the size greater than 50 percent and 84 percent of particles, respectively (known as the D50 and D84, respectively). Additionally, qualitative habitat characteristics, including the following, were recorded: Dominant land use (forest/natural, agriculture, and so on), General site habitat (primary or secondary channel; run, riffle, and [or] pool), Riparian corridor size (less than or equal to 10 meters or greater than 10 meters), Canopy cover (shaded, mostly shaded, mostly open, open),

66    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District Substrate stability (loose, somewhat loose, somewhat stable, stable), and Instream depositional features (gravel bars, detritus). References Cited Angelo, R.T., Cringan, M.S., Chamberlain, D.L., Stahl, A.J., Haslouer, S.G., and Goodrich, C.A., 2007, Residual effects of lead and zinc mining on freshwater mussels in the Spring River Basin (Kansas, Missouri, and Oklahoma, USA): Science of the Total Environment, v. 384, no. 1–3, p. 467–496. [Also available at ://doi.org/​10.1016/​j.s​ citotenv.2​007.05.045.] EcoAnalysts, Inc., 2018, Final report—Tri-State Mining District unionid assessment, Missouri, Kansas, and Oklahoma, 2016–2018: Columbia Mo., U.S. Fish and Wildlife Service, prepared by EcoAnalysts, Inc., 101 p. Ecological Specialists, Inc., 2014, Final report—Unionid mus­ sel habitat construction/creation summary: St. Louis, Mo., U.S. Army Corps of Engineers, St. Louis District, prepared by Ecological Specialists, Inc., 64 p. plus appendices. Haag, W.R., 2012, North American freshwater mussels— Natural history, ecology, and conservation: New York, Cambridge University Press, 505 p. Strayer, D.L., 2008, Freshwater mussel ecology, a multifactor approach to distribution and abundance: Berkeley, Calif., University of California Press, 204 p. Wentworth, C.K., 1922, A scale of grade and class terms for clastic sediments: Journal of Geology, v. 30, no. 5, p. 377–392.

Appendix 2.  Quality Assurance/Quality Control for Chemical Analyses    67 Appendix 2.  Quality Assurance/Quality Control for Chemical Analyses This appendix describes quality assurance and quality control measures, which were not described in the main text for the purpose of conciseness, used during the quantifica­ tion of elements (zinc, cadmium, lead, nickel, copper), major ions (fluoride, chloride, sulfate, dissolved organic carbon (DOC), sodium, magnesium, potassium, calcium, iron), and acid-volatile sulfide as appropriate in surface and pore water, sediment, and Corbicula fluminea (O.F. Müller, 1774) tissues. As described in the main text, analyses of elements and major cations (sodium, magnesium, potassium, calcium, iron) were completed using inductively coupled plasma-mass spectrom­ etry (ICP–MS). Anions were quantified using ion exchange chromatography with suppressed conductivity detection, and concentrations of DOC were quantified using a combustion analyzer. The total number of samples and the ranges of detec­ tion or quantification limits are provided in table 2.1. Quality control measures for ICP–MS analyses included the use of at least three National Institute of Standards and Technology-traceable calibration standards and a calibra­ tion blank, laboratory control standards (LCSs), interference check solutions, analysis duplicates, analysis spikes, and dilution analyses. Recoveries of the analytes in the calibration standards were 90 to 110 percent, and the internal standard recoveries were 60 to 120 percent throughout the analyses. Second source initial calibration verification standards (ICVS) and continuing calibration verification standards (CCVS) and blanks were analyzed every 10 samples; recoveries of the analytes from the ICVS and CCVS were 90 to 110 percent for nickel, copper, lead, zinc, and cadmium and 80 to 120 percent for the major cations. Analyte recoveries from the LCS and interference check solution and recoveries of the analysis spikes are summarized by analyte in table 2.2 below for ICP– MS analyses. The relative percentage differences between analysis duplicates and the percentage difference between dilution analyses also are listed in table 2.2. Additional quality control samples for sediment and C. fluminea tissues were the analysis of certified reference materials, method replicates, and method spikes (table 2.2). The quality control results indicated that the ICP–MS measurements were generally in control. Quality control measures for anion and DOC analy­ ses included the use of at least three National Institute of Standards and Technology-traceable calibration standards, analysis duplicates, and analysis spikes. A laboratory control solution was included in the anion analyses for some, but not all, analytical sequences. Recoveries of the analytes in the calibration standards for anions and DOC were generally 80 to 120 percent. ICVS and continuing calibration standards (CCVS and blanks) were analyzed every 10 samples; recov­ eries of the analytes from the ICVS and CCVS were 80 to 120 percent. Analyte recoveries from the LCSs and recoveries of the analysis spikes are summarized by analyte in table 2.2 for the anion and DOC analyses. The relative percentage dif­ ferences between analysis duplicates and the percentage differ­ ence between dilution analyses also are listed in table 2.2. The quality control results indicated that the ICP–MS measure­ ments were generally within control ranges.

68    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District Table 2.1.  Summary of reporting limits, limits of quantification, or limits of detection; sample numbers; and numbers of censored sample results. [Concentration results for major anions and cations were censored at a predetermined reporting limit (RL); dissolved organic carbon results were censored at the limit of detection (LOD), which was estimated on a batch-by-batch basis. Results for lead, cadmium, zinc, nickel, and copper were censored at the limit of quan­ tification (LOQ); the LOQs also were estimated on a batch-by-batch basis. Results for sediments and Corbicula fluminea tissues are reported on a dry-weight basis. The total number of samples analyzed for the analyte of interest, and the number of samples that were censored (percentage of total number of samples that were censored) are provided. TR, total recoverable; SEM, simultaneously extracted metals] Matrix, analyte Unit of measurement Range of RLs, LODs, or LOQs Total number of samples Number of samples censored (percentage of total) Surface water Fluoride Milligrams per liter 35 (64) Chloride Milligrams per liter 0 (0) Sulfate Milligrams per liter 0 (0) Dissolved organic carbon Milligrams per liter 0.1–0.6 2 (4) Zinc Micrograms per liter 1–40 188 (34) Cadmium Micrograms per liter 0.02–2 313 (72) Lead Micrograms per liter 0.03–1 304 (70) Nickel Micrograms per liter 0.2–0.3 0 (0) Copper Micrograms per liter 0.4–0.6 0 (0) Sodium Milligrams per liter 0 (0) Magnesium Milligrams per liter 0 (0) Potassium Milligrams per liter 0 (0) Calcium Milligrams per liter 0 (0) Iron Milligrams per liter 49 (89) Pore water Fluoride Milligrams per liter 80 (82) Chloride Milligrams per liter 1 (1) Sulfate Milligrams per liter 31 (32) Dissolved organic carbon Milligrams per liter 0.1–0.6 0 (0) Zinc Micrograms per liter 1–40 23 (13) Cadmium Micrograms per liter 0.02–2 145 (80) Lead Micrograms per liter 0.03–1 64 (35) Nickel Micrograms per liter 0.2–0.3 0 (0) Copper Micrograms per liter 0.4–0.6 21 (84) Sodium Milligrams per liter 0 (0) Magnesium Milligrams per liter 0 (0) Potassium Milligrams per liter 0 (0) Calcium Milligrams per liter 0 (0) Iron Milligrams per liter 41 (42) Corbicula fluminea tissue Zinc Micrograms per gram 1–4 0 (0) Cadmium Micrograms per gram 0 (0) Lead Micrograms per gram 0 (0) Sediment TR-nickel Micrograms per gram 0 (0) TR-copper Micrograms per gram 0 (0) TR-zinc Micrograms per gram 0 (0)

Appendix 2.  Quality Assurance/Quality Control for Chemical Analyses    69 Table 2.1.  Summary of reporting limits, limits of quantification, or limits of detection; sample numbers; and numbers of censored sample results.—Continued [Concentration results for major anions and cations were censored at a predetermined reporting limit (RL); dissolved organic carbon results were censored at the limit of detection (LOD), which was estimated on a batch-by-batch basis. Results for lead, cadmium, zinc, nickel, and copper were censored at the limit of quan­ tification (LOQ); the LOQs also were estimated on a batch-by-batch basis. Results for sediments and Corbicula fluminea tissues are reported on a dry-weight basis. The total number of samples analyzed for the analyte of interest, and the number of samples that were censored (percentage of total number of samples that were censored) are provided. TR, total recoverable; SEM, simultaneously extracted metals] Matrix, analyte Unit of measurement Range of RLs, LODs, or LOQs Total number of samples Number of samples censored (percentage of total) Sediment—Continued TR-cadmium Micrograms per gram 0.02–0.03 0 (0) TR-lead Micrograms per gram 0.1–0.3 0 (0) Acid-volatile sulfide Micromoles per gram 0.1–0.7 5 (6) SEM-nickel Micrograms per gram 0.5–3 1 (2) SEM-copper Micrograms per gram 0.2–0.3 0 (0) SEM-zinc Micrograms per gram 1–5 0 (0) SEM-cadmium Micrograms per gram 0.1–0.2 0 (0) SEM-lead Micrograms per gram 0 (0)

70    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District Table 2.2.  Summary of quality control results. [Target recoveries for quality control samples are generally 90 to 110 percent for analytes in laboratory control standard (LCS), analysis spikes, and method spikes, with the exception of major ions and dis­ solved organic carbon, which have a target recovery range of 80 to 120 percent. The target recovery range for metals in the interference check solution is 80 to 120 percent, and recoveries of analytes from certi­ fied reference materials (CRMs) depend on the quantification method. For example, a total recoverable metal in a tissue will have a much more comparable recovery (80 to 120 percent) relative to a sediment material. Sediments are certified based on total extractions rather than total recoverable or simultaneous extractions with acid-volatile sulfide; therefore, it is reasonable to expect a much lower recovery of ana­ lytes from sediment CRMs. Target differences between analysis duplicates, method replicates, and dilution analysis samples (undiluted versus 5x diluted) are generally less than a 10-percent relative standard deviation or relative percentage difference. Few outliers from the target ranges were detected; outliers are not included in the ranges below (see footnotes). %, percent; RPD, relative percentage difference; RSD, relative standard deviation; NA, not applicable, this type of quality control (QC) sample was not analyzed for this analysis; TR, total recoverable; SEM, simultaneously extracted metals] Matrix, analytes LCS (% recovery) Interference check (% recovery) Analysis spike (% recovery) Analysis duplicates (RPD) Dilution analyses (% difference) CRMs (% recovery) Method replicates (RSD or RPD) Method spikes (% recovery) Surface and pore water Fluoride NA 95.0–106.5 NA NA NA NA Chloride NA 90.0–122.0 0–8.2 NA NA NA NA Sulfate NA 88.7–116.0 0–8.0 NA NA NA NA Dissolved organic carbon NA NA 85.8–106.2 10–9.5 NA NA NA NA Zinc 93.7–110.9 83.0–119.0 90.6–107.3 0–5.9 0–7.3 NA NA NA Cadmium 95.4–105.6 286.6–128.2 91.8–106.4 0–5.2 0–5.0 NA NA NA Lead 90.7–105.0 85.0–99.8 90.2–105.2 0–4.9 0–6.3 NA NA NA Nickel 97.0–106.6 89.5–91.5 95.7–104.5 0–0.6 1.3–3.7 NA NA NA Copper 91.0–107.5 88.8–89.3 94.7–97.6 2.7–2.8 0.3–0.7 NA NA NA Sodium 91.7–101.7 96.0–106.0 96.0–104.0 0–4.5 0–8.2 NA NA NA Magnesium 97.0–106.7 98.0–104.0 94.0–104.0 0–5.0 0–5.1 NA NA NA Potassium 93.1–104.0 96.0–108.0 100.0–114.0 0–6.9 0–9.5 NA NA NA Calcium 391.4–120.0 94.0–106.0 96.0–120.0 1.9–6.9 0–3.1 NA NA NA Iron 96.0–106.0 95.6–104.0 0–4.1 0–4.1 NA NA NA Corbicula fluminea tissue Zinc 99.3–102.3 92.5–98.8 96.7–101.7 0–2.2 0.5–6.5 97.2–100.0 9.4–31.4 496.8–101.3 Cadmium 100.8–106.6 116.4–119.0 93.6–100.8 1.2–2.2 1.0–5.4 96.2–100.0 9.7–22.7 593.8–104.9 Lead 98.4–100 95.6–98.0 1.0–1.5 0.4–2.7 686.8–100.0 9.3–11.6 95.8–108.6 Sediment TR-nickel 97.6–99.0 90.4–95.1 0–1.9 0.5–1.7 83.2–86.1 78.3–95.1 TR-copper 100.9–101.1 95.1–101.1 0.9–7.2 1.2–1.8 88.7–100.0 96.7–98.9 TR-zinc 100.8–103.3 86.9–89.3 91.4–97.2 0.8–4.9 1.4–4.3 84.0–91.8 12.2–19.3 94.8–106.4 TR-cadmium 97.9–105.3 104.0–114.2 96.3–98.4 0.6–1.4 1.2–1.9 91.2–100.0 14.8–17.5 96.6–103.9 TR-lead 94.7–101.1 90.8–93.4 96.6–98.7 0–0.9 0.9–1.5 69.0–100.0 2.9–35.7 97.8–113.6 Acid-volatile sulfide NA NA 103.0–134.4 0–2.2 NA NA 0–18.2 103.4–136.5

Appendix 2.  Quality Assurance/Quality Control for Chemical Analyses    71 Table 2.2.  Summary of quality control results.—Continued [Target recoveries for quality control samples are generally 90 to 110 percent for analytes in laboratory control standard (LCS), analysis spikes, and method spikes, with the exception of major ions and dis­ solved organic carbon, which have a target recovery range of 80 to 120 percent. The target recovery range for metals in the interference check solution is 80 to 120 percent, and recoveries of analytes from certi­ fied reference materials (CRMs) depend on the quantification method. For example, a total recoverable metal in a tissue will have a much more comparable recovery (80 to 120 percent) relative to a sediment material. Sediments are certified based on total extractions rather than total recoverable or simultaneous extractions with acid-volatile sulfide; therefore, it is reasonable to expect a much lower recovery of ana­ lytes from sediment CRMs. Target differences between analysis duplicates, method replicates, and dilution analysis samples (undiluted versus 5x diluted) are generally less than a 10-percent relative standard deviation or relative percentage difference. Few outliers from the target ranges were detected; outliers are not included in the ranges below (see footnotes). %, percent; RPD, relative percentage difference; RSD, relative standard deviation; NA, not applicable, this type of quality control (QC) sample was not analyzed for this analysis; TR, total recoverable; SEM, simultaneously extracted metals] Matrix, analytes LCS (% recovery) Interference check (% recovery) Analysis spike (% recovery) Analysis duplicates (RPD) Dilution analyses (% difference) CRMs (% recovery) Method replicates (RSD or RPD) Method spikes (% recovery) Sediment—Continued SEM-nickel 95.2–103.5 88.5–103.0 99.3–100.7 1.5–2.4 22.6–44.1 88.7–97.7 SEM-copper 93.3–102.0 91.9–113.0 100.0–101.9 1.9–3.9 0–1.0 48.1–58.8 6.1–10.7 89.7–99.3 SEM-zinc 94.7–103.9 88.1–108.0 97.0–101.3 0–4.6 0–2.6 67.7–83.9 6.1–8.3 90.2–110.0 SEM-cadmium 94.4–102.4 111.0–120.2 95.2–102.0 0.7–2.2 0.7–2.7 66.7–89.8 0–8.8 89.6–105.0 SEM-lead 93.6–103.3 84.3–91.3 94.2–96.1 0–1.9 0.5–2.8 67.8–79.4 5.2–18.5 85.7–104.9 1One analysis duplicate pair for dissolved organic carbon had a 28.6-percent RPD; all other QC results indicated that the analyses were in control. 2Two cadmium results in the interference check solution had greater than expected recoveries (173.4 percent and 189.4 percent), which was likely due to an uncorrected interference in the presence of sig­ nificantly elevated concentrations of molybdenum in the check solution. Results for cadmium in the other QC samples indicated that this was isolated to the interference check solution and that the results were ultimately in control. 3Calcium in one laboratory control solution had a recovery of 140.0 percent; all other QC results were within the acceptable range. 4One method spike for zinc had a recovery of 68.2 percent. This was likely due to underspiking of the sample. 5One method spike for cadmium had a recovery of 71.9 percent. This was likely due to underspiking of the sample. 6One certified reference material had a lead recovery of 135.1 percent; all other QC results indicated that the analyses were in control.

72    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District Appendix 3.  Estimation of Effect Concentrations for Metals in Water and Sediment The estimation of effect concentrations for toxic effects of metals in water and sediment on freshwater mus­ sels (Lampsilis siliquoidea; Barnes, 1823) and amphipods (Hyalella azteca; Saussure, 1858) are shown in figures 3.1 to 3.8. In figures 3.1 and 3.2, the models were produced using Toxicity Relationship Analysis Program software (version 1.30a). When filtered zinc in the control was less than the detection limit, the background zinc concentration was estimated (for example, control=3 micrograms per liter). In figures 3.3 to 3.8, the models were derived using threeparameter logistic regression in SigmaPlot (version 14.5).

Appendix 3.  Estimation of Effect Concentrations for Metals in Water and Sediment    73 F. Water-only zinc, in micrograms per liter, versus day-84 biomass (dry weight, in milligrams) Log (exposure variable) Log (exposure variable) Log (exposure variable) Log (exposure variable) Log (exposure variable) Log (exposure variable) Effect (variable) Effect (variable) Effect (variable) Effect (variable) Effect (variable) Effect (variable) A. Water-only zinc, in micrograms per liter, versus day-28 survival B. Water-only zinc, in micrograms per liter, versus day-84 survival (of 10) C. Water-only zinc, in micrograms per liter, versus day-28 growth (average dry weight, in milligrams) D. Water-only zinc, in micrograms per liter, versus day-84 growth (average dry weight, in milligrams) E. Water-only zinc, in micrograms per liter, versus day-28 biomass (dry weight, in milligrams) Figure 3.1.  Toxicity Relationship Analysis Program estimates of mussel effect concentration values for filterable zinc (in micrograms per liter) in 12-week water-only tests with zinc, spring 2017. A, water-only zinc versus day-28 survival; B, water-only zinc versus day-84 survival; C, water-only zinc versus day-28 growth; D, water-only zinc versus day-84 growth; E, water-only zinc versus day-28 biomass; F, water-only zinc versus day-84 biomass. Circles indicate treatment means; lines indicate regressions.

74    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District A. Water-only mixture (normal zinc), in micrograms per liter, versus day-28 survival (of 10) B. Water-only mixture (normal zinc), in micrograms per liter, versus day-84 survival (of 10) C. Water-only mixture (normal zinc), in micrograms per liter, versus day-28 growth (average dry weight, in milligrams) E. Water-only mixture (normal zinc), in micrograms per liter, versus day-28 biomass (average dry weight, in milligrams) F. Water-only mixture (normal zinc), in micrograms per liter, versus day-84 biomass (average dry weight, in milligrams) D. Water-only mixture (normal zinc), in micrograms per liter, versus day-84 growth (average dry weight, in milligrams) Effect variable Effect variable Log (exposure variable) Log (exposure variable) Log (exposure variable) Effect variable Effect variable Effect variable Effect variable Log (exposure variable) Log (exposure variable) Log (exposure variable) Figure 3.2.  Toxicity Relationship Analysis Program estimates of mussel effect concentration values for filterable zinc (in micrograms per liter) in 12-week water-only tests with three-metal mixture, spring 2017. A, water-only mixture versus day-28 survival; B, water-only mixture versus day-84 survival; C, water-only mixture versus day-28 growth; D, water-only mixture versus day-84 growth; E, water-only mixture versus day-28 biomass; F, water-only mixture versus day-84 biomass. Circles indicate treatment means; lines indicate regressions.

Appendix 3.  Estimation of Effect Concentrations for Metals in Water and Sediment    75 A. Survival, day 83 B. Growth, day 83 E. Growth, days 28–83 C. Biomass, day 83 D. Survival, days 28–83 F. Biomass, days 28–83 Survival (of 10) Mean dry weight, in milligrams SEM-zinc (PEQ) SEM-zinc (PEQ) Zinc-PEQ SEM-zinc (PEQ) Biomass, in milligrams Biomass, in milligrams Mean dry weight, days 28–83, in milligrams Survival (of 10) Zinc-PEQ Zinc-PEQ Figure 3.3.  SigmaPlot estimates of mussel 20-percent effect concentrations for simultaneously extracted metals- (SEM-) zinc (as probable effect quotient [PEQ]) in 12-week sediment toxicity tests, spring 2017. A, survival, day 83; B, growth, day 83; C, biomass, day 83; D, survival, days 28–83; E, growth, days 28–83; F, biomass, days 28–83. Circles indicate treatment means; lines indicate regressions.

76    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District A. Survival, day 83 C. Biomass, day 83 SEM-zinc, lead, cadmium (sum-PEQ) Sum-PEQ Sum-PEQ Survival (of 10) Growth, in milligrams dry weight Biomass, in milligrams B. Growth, day 83 Figure 3.4.  SigmaPlot estimates of mussel 20-percent effect concentrations for sum of zinc+lead+cadmium (as probable effect quotient [PEQ]) in 12-week sediment toxicity tests, spring 2017. A, survival, day 83; B, growth, day 83; C, biomass, day 83. Circles indicate treatment means; lines indicate regressions. [SEM, simultaneously extracted metals].

Appendix 3.  Estimation of Effect Concentrations for Metals in Water and Sediment    77 C. Biomass, day 83 [ΣSEM–AVS], in micromoles per gram [ΣSEM–AVS], in micromoles per gram [ΣSEM–AVS], in micromoles per gram Survival (of 10) Mean dry weight, in milligrams Biomass, in milligrams A. Survival, day 83 B. Growth, day 83 Figure 3.5.  SigmaPlot estimates of mussel 20-percent effect concentrations for summed molar concentrations of simultaneously extracted metals relative to acid-volatile sulfide ([ΣSEM–AVS]; in micromoles per gram) in 12-week sediment toxicity tests, spring 2017. A, survival, day 83; B, growth, day 83; C, biomass, day 83. Circles indicate treatment means. Solid lines indicate regressions.

78    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District C. Biomass, day 83 PWP-zinc 1,000 PWP-zinc 1,000 PWP-zinc 1,000 Survival (of 10) Growth, in milligrams dry weight Biomass, in milligrams A. Survival, day 83 B. Growth, day 83 Figure 3.6.  SigmaPlot estimates of mussel 20-percent effect concentrations for zinc in pore-water peepers (PWP-zinc; in micrograms per liter as zinc) in 12-week sediment toxicity tests, spring 2017. A, survival, day 83; B, growth, day 83; C, biomass, day 83. Circles indicate treatment means. Solid lines indicate regressions.

Appendix 3.  Estimation of Effect Concentrations for Metals in Water and Sediment    79 B. Growth versus zinc in overlying water A. Survival versus zinc in overlying water C. Biomass, day 83 Overlying water zinc, in micrograms per liter Overlying water zinc, in micrograms per liter Overlying water zinc, in micrograms per liter Survival, day 83, in percent Growth, day 83, in milligrams dry weight Biomass, in milligrams per reproduction Figure 3.7.  SigmaPlot estimates of mussel 20-percent effect concentrations for zinc in overlying water (in micrograms per liter as zinc) in 12-week sediment toxicity tests, spring 2017. A, survival versus zinc in overlying water; B, growth versus zinc in overlying water; C, biomass, day 83. Circles indicate treatment means. Solid lines indicate regressions.

80    Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District A. Survival, day 28 B. Growth, day 28 C. Biomass, day 42 D. Reproduction, day 42 Survival, in percent Mean dry weight, in milligrams Dry weight per replicate, in milligrams Young per female SEM-zinc (PEQ) SEM-zinc (PEQ) SEM-zinc (PEQ) SEM-zinc (PEQ) Figure 3.8.  SigmaPlot estimates of amphipod 20-percent effect concentrations for simultaneously extracted metals- (SEM-) zinc (as probable effect quotient [PEQ]) in 12-week sediment toxicity tests, spring 2017. A, survival, day 28; B, growth, day 28; C, biomass, day 42; D, reproduction, day 42. Circles indicate treatment means. Solid lines indicate regressions.

For more information about this publication, contact: Director, USGS Columbia Environmental Research Center 4200 New Haven Road Columbia, MO 65201 573–875–5399 For additional information, visit: h​://www​.usgs.gov/​centers/​cerc Publishing support provided by the Rolla and Sacramento Publishing Service Centers

Besser and others—Responses of Juvenile Mussels to Metals in Sediment and Water of the Tri-State Mining District—OFR 2023–1024 ISSN 2331-1258 (online) ://doi.org/​10.3133/​ofr20231024