Water-Quality, Bed-Sediment, and Invertebrate Tissue Trace-Element Concentrations for Tributaries in the Clark Fork Basin, Montana, October 2019–September 2020

Open-File Report 2022-1090
Prepared in cooperation with the U.S. Environmental Protection Agency
By: , and 

Links

Acknowledgments

Prepared in cooperation with the U.S. Environmental Protection Agency.

Abstract

Water, bed sediment, and invertebrate tissue were sampled in streams from Butte to near Missoula, Montana, as part of a monitoring program in the Clark Fork Basin. The sampling program was completed by the U.S. Geological Survey, in cooperation with the U.S. Environmental Protection Agency, to characterize aquatic resources in the Clark Fork Basin and monitor trace elements associated with historical mining and smelting activities. Sampling sites were on the Clark Fork River and a subset of its tributaries. Water samples were collected periodically at 22 sites from October 2019 through September 2020. Bed-sediment and tissue samples were collected once at 12 sites during July 2020.

Water-quality data included concentrations of major ions, dissolved organic carbon, nitrogen (nitrate plus nitrite), trace elements, and suspended sediment. Daily values of turbidity were determined at four sites. Bed-sediment data included trace-element concentrations in the fine-grained (less than 0.063 millimeter) fraction. Biological data included trace-element concentrations in whole-body tissue of selected aquatic benthic invertebrates. Statistical summaries of water-quality, bed-sediment, and invertebrate tissue trace-element data for sites in the Clark Fork Basin were provided for the period of record: March 1985–September 2020.

Introduction

The Clark Fork originates near the town of Warm Springs in western Montana at the confluence of Silver Bow and Warm Springs Creeks (fig. 1). Along the 148-mile reach from Silver Bow Creek in Butte to the Clark Fork near Missoula, six major tributaries enter: Blacktail Creek, Warm Springs Creek, Little Blackfoot River, Flint Creek, Rock Creek, and Blackfoot River. Principal surface-water uses in the 6,000-square-mile Clark Fork Basin above Missoula include irrigation, stock watering, small-scale industry (Cannon and Johnson, 2004), and habitat for trout fisheries. Current (as of 2020) land uses are primarily cattle production, logging, mining, residential development, and recreation. Large-scale mining and smelting were prevalent in the basin for more than 100 years but are now either discontinued or reduced in scale.

Map showing creeks and rivers and sites, labeled by type, in the study area.
Figure 1.

Location of the study area in the Clark Fork Basin, Montana.

Copper, gold, silver, and lead ores were extensively mined, milled, and smelted in the drainages of Silver Bow and Warm Springs Creeks from about the 1860s to the 1980s. Moderate- and small-scale mining also took place in the basins of most of the major tributaries to the upper Clark Fork Basin. Tailings produced during past mineral processing commonly contained large quantities of trace elements such as arsenic, cadmium, copper, lead, and zinc. Eroded tailings mixed with stream sediment and deposited downstream in stream channels, on flood plains, in the Warm Springs Ponds, and at the former Milltown Reservoir (fig. 1; Andrews, 1987), whose dam (Milltown Dam, not shown) was breached on March 28, 2008.

Concern about the toxicity of trace elements to the aquatic ecosystem and human health has resulted in a comprehensive effort by State, Federal, Tribal, and private entities to characterize the aquatic resources in the Clark Fork Basin. This effort was designed to guide and monitor remedial activities and to evaluate the effectiveness of remediation and cleanup. Water-quality data have been collected by the U.S. Geological Survey (USGS) in the Clark Fork Basin since 1985 (Clark and others [2021] and references therein). Trace-element data for bed sediment and aquatic benthic invertebrates have been collected intermittently since 1986 as part of studies on the contamination of bed-sediment quality and bioaccumulation of metals led by the USGS (Axtmann and Luoma, 1991; Cain and others, 1992, 1995; Axtmann and others, 1997; Hornberger and others, 1997). In March 1993, these efforts merged into a long-term (March 1985–September 2019) monitoring program for water, bed sediment, and tissue in the basin. The program was implemented by the USGS, in cooperation with the U.S. Environmental Protection Agency, to systematically quantify the seasonal and annual variability of selected constituents.

The purpose of this report is to present water-quality data from samples collected at 22 sites and bed-sediment and biological data from samples collected at 12 sites in the Clark Fork Basin from October 2019 through September 2020 (fig. 1). Associated quality-assurance data also are presented. Statistical summaries are provided for all data collected for the period of record in an accompanying USGS data release (Clark and others, 2022).

Sampling Locations and Data Types

Sampling sites for the monitoring program in the Clark Fork Basin are from Butte to near Missoula (fig. 1). The sampling sites include sites on the Clark Fork main stem, including Silver Bow Creek, and three major tributaries (Blacktail Creek, Warm Springs Creek, and Blackfoot River) and three smaller tributaries (Mill Creek, Willow Creek, and Lost Creek). The sites, types of data collected, and period of record for each data type are listed in table 1. Main-stem sampling sites were selected to divide the Clark Fork Basin into reaches of approximate uniform length. Each reach encompasses either a major tributary or a depositional environment (that is, Warm Springs Ponds or the former Milltown Reservoir). Major tributaries were sampled to describe water-quality, bed-sediment, and biological characteristics of important hydrologic sources in the Clark Fork Basin. Major tributaries also provided reference comparisons to the main stem. The three smaller tributaries were sampled to gain better spatial resolution on sources of trace elements entering the Clark Fork Basin in an area of historical metal-processing activities near Anaconda, Montana. Water-quality samples were collected at 22 sites, 6 to 8 times per year on a schedule designed to represent seasonal and hydrological variability. Daily turbidity values were computed using continuous turbidity monitors recording data every 15 minutes at four sites. Bed-sediment and invertebrate tissue samples were collected at 12 sites. Continuous streamflow data measured every 15 minutes were collected at 21 sites.

Table 1.    

Type and period of data collection at sampling sites in the Clark Fork Basin, Montana, through September 2020.

[Dates are listed in month/year format. USGS, U.S. Geological Survey; --, no data; P, present (September 2020); D, discontinued]

USGS site number (fig. 1) USGS site name Continuous-record streamflow Periodic
water quality1
Daily suspended sediment Daily turbidity Fine-grained bed sediment2 Tissue2
12323230 Blacktail Creek at Harrison Avenue, at Butte -- 03/1993–08/1995, 12/1996–08/2003, 12/2004–06/2020 -- -- -- --
12323233 Blacktail Creek above Grove Gulch, at Butte 06/2020–P 06/2020–P -- -- -- --
12323242 Silver Bow Creek at Montana Street, at Butte 06/2020–P 06/2020–P -- -- -- --
12323250 Silver Bow Creek below Blacktail Creek, at Butte 10/1983–P 03/1993–08/1995,
12/1996–P
-- -- -- --
12323600 Silver Bow Creek at Opportunity 07/1988–P 03/1993–08/1995,
12/1996–P
03/1993–09/1995, D -- 07/1992–P 07/1992, 08/1994–08/1995, 08/1997–P
12323670 Mill Creek near Anaconda 10/2004–P 12/2004–P -- 06/2006–09/2012, D -- --
12323700 Mill Creek at Opportunity 04/2003–P 03/2003–P -- 04/2013–P -- --
12323710 Willow Creek near Anaconda 03/2005–P 12/2004–P -- 06/2006–09/2012, D -- --
12323720 Willow Creek at Opportunity 04/2003–P 03/2003–P -- 04/2013–P -- --
12323750 Silver Bow Creek at Warm Springs 03/1972–09/1979,
04/1993–P
03/1993–P 04/1993–09/1995, D -- 07/1992–P 07/1992–P
12323760 Warm Springs Creek near Anaconda 10/1997–P 10/2005–P -- 05/2006–09/2012, D -- --
12323770 Warm Springs Creek at Warm Springs 10/1983–P 03/1993–P -- 04/2013–P 08/1995, 08/1997, 08/1999, 08/2002, 08/2005, 08/2008, 08/2011, 08/2014, 08/2017, 07/2020 08/1995, 08/1997, 08/1999, 08/2002, 08/2005, 08/2008, 08/2011, 08/2014, 08/2017, 07/2020
12323800 Clark Fork near Galen 07/1988–P 07/1988–P -- -- 08/1987, 08/1991–P 08/1987, 08/1991–P
12323840 Lost Creek near Anaconda 10/2004–P 12/2004–P -- 05/2006–P -- --
12323850 Lost Creek near Galen 04/2003–P 03/2003–P -- -- -- --
461415112450801 Clark Fork below Lost Creek, near Galen -- -- -- -- 08/1996–P 08/1996–P
461559112443301 Clark Fork at county bridge, near Racetrack -- -- -- -- 08/1996–P 08/1996–P
461903112440701 Clark Fork at Dempsey Creek diversion, near Racetrack -- -- -- -- 08/1996–P 08/1996–P
12324200 Clark Fork at Deer Lodge 10/1978–P 03/1985–P 03/1985–08/1986, 04/1987–03/2003, 08/2003–2014, D 03/2016–09/2016, D 08/1986–08/1987,
08/1990–P
08/1986–08/1987,
08/1990–P
12324400 Clark Fork above Little Blackfoot River, near Garrison 02/2009–P 03/2009–P -- -- 08/2009–08/2019 08/2009–08/2019
12324680 Clark Fork at Goldcreek 10/1977–P 03/1993–P -- -- 07/1992–P 07/1992–P
12331800 Clark Fork near Drummond 04/1993–P 03/1993–P -- -- 08/1986, 08/1987,
08/1991–08/2019
08/1986, 08/1991–08/2019
12334550 Clark Fork at Turah Bridge, near Bonner 03/1985–P 03/1985–P 03/1985–03/2003,
08/2003–09/2016, D
-- 08/1986, 08/1991–P 08/1986, 08/1991–P
12340000 Blackfoot River near Bonner 10/1939–P 03/1985–P 07/1986–04/1987,
06/1988–09/1995,
10/2005–09/216, D
-- 08/1986, 8/1987, 08/1991, 08/1993–08/1996, 08/1998–08/2001, 09/2003,
08/2006–P
08/1986, 08/1987, 08/1991, 08/1993, 08/1996, 08/1998, 09/2000, 09/2003,
08/2006–P
12340500 Clark Fork above Missoula 03/1929–P 07/1986–P3 07/1986–04/1987, 06/1988–01/1996, 03/1996–03/2003, 08/2003–09/216, D 04/2007–09/2007, D 08/1997–P 08/1997–P
Table 1.    Type and period of data collection at sampling sites in the Clark Fork Basin, Montana, through September 2020.
1

Onsite measurements of physical properties and laboratory analyses for selected major ions, trace elements, and suspended sediment. Before March 1993, laboratory analyses included only trace elements and suspended sediment. Beginning in 2012, dissolved organic carbon and turbidity analyses were included at selected sites. Beginning in 2013, nutrient sample analyses were included for two sites (12323230 and 12323250) near Butte, Montana.

2

Once annual laboratory analyses of fine-grained bed sediment and aquatic benthic insect tissue for trace elements.

3

Before October 1989, water-quality data for Clark Fork above Missoula included only suspended-sediment data.

Properties measured onsite and constituents for which water, bed-sediment, and biota samples were analyzed are listed in table 2. Data-quality objectives for analyses of water samples are listed in table 3. Results of onsite measurements of stream properties; laboratory analyses of water-quality, bed-sediment, and biota samples; quality-assurance data for water year 2020 (October 1, 2019, through September 30, 2020); and statistical summaries of long-term water-quality, bed-sediment, and tissue data collected between March 1985 and September 2020 are provided in the accompanying data release (Clark and others, 2022).

Table 2.    

Properties and constituents measured onsite or analyzed in water, bed-sediment, and biota samples from the Clark Fork Basin, Montana.
Water Bed sediment Tissue
Property Constituent Constituent Constituent
Streamflow Hardness (calculated) Arsenic Arsenic
pH Calcium Cadmium Cadmium
Specific conductance Magnesium Chromium Chromium
Temperature Potassium Copper Copper
Turbidity Sodium Iron Iron
Alkalinity Lead Lead
Chloride Manganese Manganese
Fluoride Nickel Nickel
Silica Zinc Zinc
Sulfate
Nitrate plus nitrite
Cadmium
Copper
Iron
Lead
Manganese
Zinc
Arsenic
Dissolved organic carbon
Suspended sediment
Table 2.    Properties and constituents measured onsite or analyzed in water, bed-sediment, and biota samples from the Clark Fork Basin, Montana.

Table 3.    

Data-quality objectives for analyses of water samples collected in the Clark Fork Basin, Montana.

[lab, laboratory; NTRU, nephelometric turbidity ratio unit; --, not determined; mg/L, milligram per liter; µg/L, microgram per liter; mm, millimeter]

Constituent Data-quality objectives
Detectability Precision Bias
Laboratory reporting
level
Maximum relative standard deviation of replicate analyses (percent) Maximum deviation of spike recovery (percent)
Turbidity, unfiltered, lab 2.0 NTRU 20 --
Calcium, filtered 0.04 mg/L 20 --
Magnesium, filtered 0.02 mg/L 20 --
Potassium, filtered 0.6 mg/L 20 --
Sodium, filtered 0.8 mg/L 20 --
Alkalinity, filtered, lab 8.0 mg/L 20 --
Chloride, filtered 0.04 mg/L 20 --
Fluoride, filtered 0.02 mg/L 20 --
Silica, filtered 0.1 mg/L 20 --
Sulfate, filtered 0.04 mg/L 20 --
Nitrate plus nitrite, filtered 0.02 mg/L 20 --
Cadmium, filtered 0.06 µg/L 20 25
Cadmium, unfiltered recoverable 0.06 µg/L 20 25
Copper, filtered 0.8 µg/L 20 25
Copper, unfiltered recoverable 0.8 µg/L 20 25
Iron, filtered 10 µg/L 20 25
Iron, unfiltered recoverable 10 µg/L 20 25
Lead, filtered 0.04 µg/L 20 25
Lead, unfiltered recoverable 0.12 µg/L 20 25
Manganese, filtered 0.8 µg/L 20 25
Manganese, unfiltered recoverable 0.8 µg/L 20 25
Zinc, filtered 4 µg/L 20 25
Zinc, unfiltered recoverable 4 µg/L 20 25
Arsenic, filtered 0.2 µg/L 20 25
Arsenic, unfiltered recoverable 0.2 µg/L 20 25
Organic carbon, filtered 0.46 mg/L 20 --
Sediment, suspended, percentage finer than 0.062 mm 1 percent 20 --
Sediment, suspended 1 mg/L 20 --
Table 3.    Data-quality objectives for analyses of water samples collected in the Clark Fork Basin, Montana.

Quality assurance of data was maintained using documented procedures described in the following sections. These quality-assurance data were designed to provide environmentally representative data. Acceptable results of the procedures were verified with quality-assurance samples that were collected systematically to provide a measure of the accuracy, precision, and bias of the environmental data and to identify variability associated with sampling, processing, or analysis.

Trace-Element Concentrations and Physical Properties of Surface-Water Samples

Water-quality data consist of onsite stream properties and laboratory determination of concentrations of chemical and physical constituents (listed in table 2) in stream samples. Water samples were collected at 22 sites in the Clark Fork Basin 6 to 8 times per year on a schedule designed to describe seasonal and hydrologic variability. Turbidity monitors were operated seasonally (March to September 2020) at four sites; continuous turbidity data (recorded every 15 minutes) were used to compute daily mean turbidity values. Water-quality data are available through the USGS National Water Information System (NWIS) database at https://waterdata.usgs.gov/nwis/qw (U.S. Geological Survey, 2022).

Methods

Water samples were collected and composited from vertical transits throughout the stream depth at multiple locations along the channel cross section using depth- and width-integration methods described by Ward and Harr (1990), Edwards and Glysson (1999), and the USGS (U.S. Geological Survey, variously dated). These methods provided vertically and laterally discharge-weighted composite samples that represented the flow passing through the cross section of the stream. Samples were collected with isokinetic depth-integrating water-quality samplers (Davis and the Federal Interagency Sedimentation Project, 2005) constructed of plastic or coated with a nonmetallic rubber-coating paint and equipped with polytetrafluoroethylene nozzles.

Instantaneous streamflow was determined at the time of water sampling either by direct measurement or from stage-discharge rating tables (Rantz and others, 1982). Daily mean streamflow values during ice periods were labeled as estimated because backwater affected the stage-discharge relation. Onsite measurements of pH, specific conductance, and water temperature were measured during water sample collection. Sample filtration and preservation were completed according to procedures described by Ward and Harr (1990), Horowitz and others (1994), and the USGS (U.S. Geological Survey, variously dated).

Concentrations of arsenic, cadmium, copper, lead, manganese, and zinc in filtered samples (0.45-micrometer [µm] pore size) were measured using inductively coupled plasma-mass spectrometry (ICP–MS; Garbarino and others, 2006). Concentrations of calcium, magnesium, and iron in filtered samples were measured using inductively coupled plasma-atomic emission spectrometry (Fishman and Friedman, 1989). Calcium and magnesium concentrations were used to calculate water hardness. Concentrations of potassium, sodium, alkalinity, chloride, fluoride, silica, sulfate, nitrogen (nitrate plus nitrite), and dissolved organic carbon were measured in filtered samples collected at selected sites in the upper Clark Fork Basin. Potassium was measured by Standard Method 3120 (American Public Health Association, 1998); sodium and silica were analyzed by inductively coupled plasma-atomic emission spectrometry (Fishman, 1993); alkalinity was measured by electrometric titration (Fishman and Friedman, 1989); chloride, fluoride, and sulfate were measured by ion chromatography (Fishman and Friedman, 1989); nitrogen (nitrate plus nitrite) was measured by colorimetric enzymatic reduction (Patton and Kryskalla, 2011); and dissolved organic carbon was measured by persulfate oxidation and infrared spectrometry (Brenton and Arnett, 1993). Recoverable concentrations of trace elements were measured in unfiltered samples that were first digested with dilute hydrochloric acid (Hoffman and others, 1996). For cadmium, iron, lead, and manganese, the digested samples were analyzed by ICP–MS as described by Garbarino and Struzeski (1998). For arsenic, copper, and zinc, the digested samples were analyzed by ICP–MS as described by Garbarino and others (2006). Turbidity was measured using Standard Method 2130 (American Public Health Association, 1998) in selected unfiltered samples. All samples were analyzed at the USGS National Water Quality Laboratory (NWQL) in Denver, Colorado.

Water samples for analysis of suspended sediment were collected from multiple vertical transits when periodic water samples were collected. Water samples were analyzed for suspended-sediment concentration and the percentage of suspended-sediment mass finer than a 0.062-millimeter (mm) diameter (silt size and smaller) by the USGS Wyoming-Montana Water Science Center Sediment Laboratory (hereinafter referred to as the “Wyoming-Montana Sediment Laboratory”) in Helena, Mont., according to methods described by Guy (1969) and Dodge and Lambing (2006).

Continuous turbidity was measured from early spring (after ice breakup) to early fall (before stream freezeup) using model 6136 turbidity sensors (YSI, Inc, Yellow Springs, Ohio) at four tributary sites in the upper Clark Fork Basin near Anaconda (table 1). Turbidity values were recorded every 15 minutes and are stored in the USGS NWIS database (U.S. Geological Survey, 2022) at https://waterdata.usgs.gov/mt/nwis/current?type=quality. The in-field values differ from the results of laboratory-measured turbidity in discrete water-quality samples because of differences in the instrumentation and sampling procedures. Continuous recordings enable determination of the minimum and maximum turbidity values for each day and a daily mean turbidity, which is based on the average of all values in a 24-hour period. Procedures for the operation of continuous turbidity monitors and for daily record computations are described by Wagner and others (2006).

Results

Water-quality data from samples collected periodically during water year 2020 (October 1, 2019, through September 30, 2020) are listed in the accompanying data release (“table_4_water_quality_data_Clark_Fork_Basin.xlsx” in Clark and others [2022]). Daily maximum, minimum, and mean turbidity at four sites are listed in the accompanying data release along with monthly summary statistics (“table_5_daily_turbidity_at_Mill_Creek_at_Opportunity.xlsx,” “table_6_daily_turbidity_at_Willow_Creek_at_Opportunity.xlsx,” “table_7_daily_turbidity_at_Warm_Springs_Creek_at_Warm_Springs.xlsx,” and “table_8_daily_turbidity_at_Lost_Creek_near_Anaconda.xlsx” in Clark and others [2022]). Water-quality data are available in the USGS NWIS database (U.S. Geological Survey, 2022).

Quality Assurance

Quality-assurance procedures used for the collection and field processing of water samples are described by Ward and Harr (1990), Horowitz and others (1994), Edwards and Glysson (1999), Lambing (2006), and the USGS (U.S. Geological Survey, variously dated). Standard procedures used by the NWQL for internal sample handling and quality assurance are described by Friedman and Erdmann (1982), Jones (1987), and Pritt and Raese (1995). Quality-assurance procedures used by the Wyoming-Montana Sediment Laboratory are described by Dodge and Lambing (2006). Standard procedures used for the calibration, measurement, and quality assurance of turbidity monitors are described by Anderson (2005).

The quality of analytical results reported for water samples was evaluated using quality-control samples that were sampled and analyzed concurrently with primary environmental samples. These quality-control samples consisted of replicates, spikes, and blanks that provided quantitative information on the precision and bias of the overall field and laboratory processes. The number of quality-control samples represented about 15 percent of the total number of water samples.

Replicate data provided an assessment of the precision (reproducibility) of analytical results and collection variability (U.S. Geological Survey, variously dated). Replicate samples were collected in the field (field replicate) either by repeating the collection process (sequential or concurrent replicate) to obtain two or more independent composite samples or by splitting a single composite sample into two or more subsamples (split replicate). The individual replicate samples were analyzed separately. Likewise, a single sample can be analyzed two or more times in the laboratory to obtain a measure of analytical precision (laboratory replicate).

Precision of analytical results for field replicates can be affected by numerous sources of variability within the field and laboratory environments, including sample collection, processing, and analysis. Overall precision for samples exposed to field and laboratory sources of variability were provided by obtaining replicate stream samples for chemical analysis by splitting a composite stream sample. Replicate stream samples for suspended-sediment analysis were obtained in the field by collecting two independent cross-sectional samples. Analyses of field replicate samples indicated the reproducibility of environmental data that were affected by the combined potential variability introduced by field and laboratory processes.

In addition to analyzing quality-control samples submitted from the field, internal quality-assurance practices are completed systematically by the NWQL to provide quality control of analytical procedures (D.L. Stevenson, U.S. Geological Survey, written commun., 2012). These internal practices include analyses of quality-control samples such as calibration standard samples, standard reference water samples, replicate samples, deionized-water blank samples, or spiked samples at a proportion equivalent to at least 10 percent of the samples loaded. The NWQL participates in a blind-sample program in which standard reference water samples prepared by the USGS Quality Systems Branch are routinely inserted into the sample line for each analytical method at a frequency proportional to the sample load. The laboratory also participates in external evaluation studies and audits with the National Environmental Laboratory Accreditation Program, the U.S. Environmental Protection Agency, Environment and Climate Change Canada, and the USGS Quality Systems Branch to assess analytical performance.

Precision of analytical results for laboratory replicates, which exclude field sources of variability, was determined using two independent chemical analyses of aliquots from a single sample selected from the group of samples constituting each analytical run. A separate analysis of the sample was completed at the beginning and end of each analytical run to provide information on the reproducibility of laboratory analytical results independent of variability caused by field sample collection and processing. Laboratory replicates of suspended-sediment samples were not obtainable because the samples were consumed during the analysis.

Spiked samples were used to evaluate bias, which measures the ability of an analytical method to accurately quantify a known amount of analyte added to a sample. In the laboratory, deionized-water blank samples and aliquots of stream samples were spiked with known amounts of the same trace elements for which water samples were being analyzed. Analyses of spiked blanks indicated if the spiking procedure and analytical method are within control for water that is presumably free of chemical interferences. Analyses of spiked aliquots of stream samples indicated if the chemical matrix of the stream water interferes with the analytical measurement and if these interferences could contribute substantial bias to reported trace-element concentrations for stream samples.

Field blank samples were submitted for every field trip and analyzed to identify the presence and magnitude of contamination that could bias analytical results. Field blanks consisted of deionized water that is certified as constituent free and is processed in the field through clean sampling equipment used to collect stream samples. These blanks were subjected to the same processing (sample splitting, filtration, preservation, transportation, and laboratory handling) as stream samples. Blank samples were analyzed for the same constituents as stream samples to detect contamination.

All water samples were handled in accordance with chain-of-custody procedures that provide documentation of sample identity, shipment, receipt, and laboratory handling (Driscoll and Hatcher, 2010). All environmental and quality-control samples submitted from a sampling episode were stored in a secure area of the NWQL and analyzed as a discrete sample group, independent of other samples submitted to the NWQL; therefore, the quality-control data apply solely to the analytical results for stream samples reported herein and provide a direct measure of data quality for this study.

Data-quality objectives (table 3) were established by the USGS and U.S. Environmental Protection Agency as part of the study plan for the expanded monitoring program started in 1993. The objectives identify the analytical requirements of detectability and serve as a guide for identifying questionable data by establishing acceptable limits for precision and bias of laboratory results. Comparisons of quality-control data to data-quality objectives were used to evaluate if sampling and analytical procedures produced environmentally representative data in a consistent manner. Data that did not meet the objectives were evaluated for acceptability; if necessary, additional quality-control samples were submitted and corrective action was taken.

The precision of analytical results was determined by calculating the standard deviation of the differences in concentrations between replicate analyses. These replicate analyses consisted of pairs of field replicate samples and laboratory replicates. Standard deviations were calculated according to the following equation (Taylor, 1987):

S = d 2 2 k
(1)
where

S

is the standard deviation of the difference in concentration between replicate analyses,

d

is the difference in concentration between each pair of replicate analyses, and

k

is the number of pairs of replicate analyses.

Precision also was expressed as a relative standard deviation (RSD), in percent, which was computed from the standard deviation and the mean concentration for all the replicate analyses. Expressing precision relative to a mean concentration standardized the comparison of precision among individual constituents. The RSD was calculated according to the following equation (Taylor, 1987):

R S D = S x ¯ × 100
(2)
where

RSD

is the relative standard deviation,

S

is the standard deviation, and

̅x

is the mean concentration for all replicate analyses.

Sample results and the corresponding field replicate data are listed in the accompanying data release (“table_9_results_field_replicates_Clark_Fork_Basin.xlsx” in Clark and others [2022]). The precision for each constituent estimated from field replicate analyses and from laboratory replicate analyses is listed in the accompanying data release (“table_10_precision_field_replicates_Clark_Fork_Basin.xlsx” and “table_11_precision_laboratory_replicates_Clark_Fork_Basin.xlsx” in Clark and others [2022]). Statistics summarizing the precision of analytical results for field and laboratory replicates were calculated using unrounded values. Precision estimates for field and laboratory replicates were within the acceptable 20-percent RSD limit for all constituents except for suspended-sediment concentration, which resulted from one sample pair having concentrations exceeding 350 milligrams per liter (mg/L; “table_10_precision_field_replicates_Clark_Fork_Basin.xlsx” in Clark and others [2022]). No adjustments were made to analytical data based on replicate analysis precision.

Recovery efficiencies for constituent analyses were determined by comparing a sample and a spiked aliquot of the same sample. The data-quality objective for acceptable spike recoveries of trace elements in water samples was a maximum deviation of 25 percent from a theoretical 100-percent recovery of an added constituent. At the NWQL, a spiked deionized-water blank sample and a spiked aliquot of a stream sample were prepared and analyzed along with the original unspiked sample. The differences between the spiked and unspiked sample concentrations were determined and used to compute recovery, in percent, according to equation 3:

R = D C × 100
(3)
where

R

is the spike recovery, in percent;

D

is the difference between the spiked and unspiked sample concentrations; and

C

is the concentration of material used to spike the sample.

If the spike recovery of a trace element was outside a range of 75 to 125 percent, the instrument was recalibrated, and the sample set and all spiked samples were reanalyzed for that element until recoveries were improved to the extent possible. Recovery efficiency for individual trace elements in spiked blank samples and in spiked stream samples is listed in the accompanying data release (“table_12_recovery_efficiency_deionized_water_blank_samples.xlsx” and “table_13_recovery_efficiency_stream_samples.xlsx” in Clark and others [2022]). The mean recovery for spiked blank samples ranged from 97.6 to 104 percent (“table_12_recovery_efficiency_deionized_water_blank_samples.xlsx” in Clark and others [2022]). The 95-percent confidence intervals (Taylor, 1987) for the mean recovery for each constituent for which spiked blank samples were analyzed did not exceed a 25-percent deviation from an expected 100-percent recovery. The mean recovery for spiked stream samples (“table_13_recovery_efficiency_stream_samples.xlsx” in Clark and others [2022]) ranged from 80.4 to 111 percent. The 95-percent confidence intervals for the mean spike recovery for each constituent for which stream water samples were analyzed did not exceed a 25-percent deviation from an expected 100-percent recovery except for copper, unfiltered recoverable (33.3 and 128 percent); iron, unfiltered recoverable (66.7 percent); lead, unfiltered recoverable (37.4 and 126 percent); manganese, unfiltered recoverable (29.6 and 132 percent); zinc, unfiltered recoverable (35.1 and 138 percent); arsenic, filtered (60.1 and 125 percent); and arsenic, unfiltered recoverable (129 percent). No adjustments were made to analytical data based on the mean spike recovery.

High or low bias is indicated if the 95-percent confidence interval does not include 100-percent recovery, thereby indicating a consistent deviation or bias, either high or low. Confidence intervals for percentage recovery include 100 percent for all laboratory-spiked blank samples (“table_12_recovery_efficiency_deionized_water_blank_samples.xlsx” in Clark and others [2022]) except for cadmium, unfiltered recoverable (102–105 percent); lead, filtered (96.5–99.5 percent); and manganese, unfiltered recoverable (96.2–99.0). Confidence intervals for percentage recovery include 100 percent for all laboratory-spiked stream samples except for iron, filtered (86.8–99.1 percent; “table_13_recovery_efficiency_stream_samples.xlsx” in Clark and others [2022]). No adjustments were made to analytical results for stream samples based on spike recoveries.

Analytical results for field blanks are listed in the accompanying data release (“table_14_analyses_of_field_blanks_for_water_samples.xlsx” in Clark and others [2022]). Field blanks with constituent concentrations less than or equal to the laboratory reporting level (LRL) indicated that the process of sample collection, field processing, and laboratory analysis is presumably free of contamination. If detectable concentrations of trace elements in field blanks were greater than or equal to twice the LRL, the concentrations were noted during data review. Analytical results from the field blank were evaluated for evidence of a consistent trend that could indicate systematic contamination. Sporadic, infrequent, nonconsecutive exceedances of twice the LRL most likely represented random contamination or laboratory-instrument calibration error that was not persistent in the process and was not likely to cause positive bias in a long-term record of analytical results; however, if concentrations for a particular constituent exceeded twice the LRL in field blanks from two consecutive field trips, additional blank samples were collected and submitted for analysis to identify the source of contamination.

Constituent concentrations in field blanks almost always were less than the LRL. Two sample concentrations of chloride, filtered (0.04 microgram per liter [µg/L] and 0.06 µg/L), exceeded the LRL of 0.02 µg/L. Two sample concentrations of silica, filtered (0.10 mg/L and 0.08 mg/L), exceeded the LRL of 0.05 mg/L. One sample concentration of copper, filtered (1.4 µg/L), exceeded the LRL of 0.4 µg/L. Two sample concentrations of copper, unfiltered recoverable (1.1 µg/L, and 2.6 µg/L), exceeded the LRL of 0.4 µg/L. One sample concentration of iron, unfiltered recoverable (5.2 µg/L), exceeded the LRL of 5.0 µg/L. Two sample concentrations of organic carbon, filtered (0.26 mg/L and 0.30 mg/L), exceeded the LRL of 0.23 mg/L. No adjustments were made to water-quality sample data based on a review of these results.

Bed-Sediment Data

Bed-sediment data for the long-term monitoring program in the Clark Fork Basin consisted of trace-element concentrations in the fine-grained (less than 0.063 mm) fraction of bed-sediment samples. Bed-sediment samples were collected once annually at 12 sites (fig. 1; table 1) during low, stable flow conditions at the same time of year as previous samples to facilitate data comparisons among years. Travel restrictions during the 2020 pandemic limited field crew to two personnel and Clark Fork above Little Blackfoot River, near Garrison (12324400), and Clark Fork near Drummond (12331800) were not sampled because of site access challenges with a limited crew. Water year 2020 bed-sediment data are available in Clark and others (2022). Warm Springs Creek at Warm Springs (12323770) is sampled once every 3 years and was sampled during water year 2020.

Methods

Fine-grained bed-sediment samples were collected in July 2020 using protocols described by Axtmann and Luoma (1991). Samples were collected from the surfaces of streambed deposits in areas near the edge of the stream using an acid-washed polypropylene scoop. Whenever possible, samples were collected from both sides of the stream.

Individual samples of bed sediment were collected by scooping material from the surfaces of three to five randomly selected deposits along pools or low-velocity areas. The three to five individual samples were combined to form a single composite sample. This collection process was repeated three times to obtain three composite samples. Each composite sample was wet sieved onsite through a 0.063-mm polyester-mesh sieve using ambient stream water. The fraction of bed sediment in each composite sample that was finer than 0.063 mm was collected in an acid-washed 500-milliliter (mL) polyethylene bottle and transported to the laboratory on ice.

Bed-sediment samples were processed and analyzed at the USGS Metal Bioavailability Laboratory in Menlo Park, California. Bed-sediment samples were oven dried at 60 degrees Celsius and ground into smaller particle sizes using an acid-washed, ceramic mortar and pestle. Single aliquots of about 0.5–0.6 gram of sediment from each of the three composite bed-sediment samples were digested using a hot, concentrated nitric acid reflux according to methods described by Luoma and Bryan (1981). Laboratory replicates were analyzed by taking an aliquot from one of the three sieved replicate samples at each site. After a 2-week digestion period, the aliquots were evaporated to dryness on a hot plate. The dry residue was reconstituted in 10 mL of 0.6 N (normal) hydrochloric acid. The reconstituted aliquots were then filtered through a 0.45-µm pore-size filter using a syringe and an in-line disposable filter cartridge. The filtrate was diluted to a 1 to 10 ratio with 0.6 N hydrochloric acid. These final solutions were analyzed for arsenic, cadmium, chromium, copper, iron, lead, manganese, nickel, and zinc using inductively coupled plasma-optical emission spectrometry (Hornberger and others, 1997). The smallest concentration of a constituent that can be reliably reported for analyses of bed sediment is termed the minimum reporting level (MRL).

Results

Solid-phase concentrations of trace elements measured in samples of fine-grained bed sediment collected during July 2020 are listed in the accompanying data release (“table_15_bed_sediment_data_Clark_Fork_Basin.xlsx” in Clark and others [2022]). Liquid-phase concentrations, measured in micrograms per milliliter, were analyzed in the reconstituted aliquots of digested bed sediment. Solid-phase concentrations, measured in micrograms per gram, were calculated using the following equation:

μg/g = μg/mL volume of digested sample, in mL dilution ratio dry weight of sample, in grams
(4)
where

µg/g

is micrograms of trace element per gram of sediment, by dry weight;

µg/mL

is micrograms per milliliter of liquid-phase trace element; and

mL

is milliliters.

The reported solid-phase concentrations (Clark and others, 2022) are the means of all analyses for replicate aliquots from each composite bed-sediment sample collected at each site. Because the conversion from liquid-phase to solid-phase concentration is dependent on the dilution ratio and the dry weight of the sample, MRLs for some trace elements might differ among sites and among years.

Quality Assurance

The USGS protocols for field collection and processing of bed-sediment samples are designed to prevent contamination from metal sources. Nonmetallic sampling and processing equipment (white plastic scoop, funnel-frame apparatus, and 500-mL sample bottles) was acid washed and rinsed with deionized water before the collection of the first sample. Polyester-mesh sieves were washed in laboratory-grade detergent and rinsed with deionized water. All equipment received a field rinse onsite with native water. Sampling equipment used at more than one site was rinsed thoroughly between sites with site-specific stream water. Separate sieves were used at each site and, therefore, did not require decontamination between sites.

Quality assurance of analytical results for bed-sediment samples included laboratory-instrument calibration with standard solutions and analysis of quality-control samples designed to identify the presence and magnitude of bias (Ellen V. Axtmann, U.S. Geological Survey, written commun., 1994). Quality-control samples consisted of standard reference materials (SRMs) issued by the National Institute of Standards and Technology (NIST) and procedural blanks. In total, 12 low-concentration SRMs, 12 high-concentration SRMs, and 13 procedural blanks were analyzed.

SRMs are commercially prepared materials that have certified concentrations of trace elements. Analyses of SRMs were used to indicate the ability of the method to accurately measure a known quantity of a constituent. Multiple analyses of SRMs were made to derive a mean and 95-percent confidence interval for recovery. Recovery efficiency for trace-element analyses of SRMs for bed sediment is listed in the accompanying data release (“table_16_recovery_efficiency_bed_sediment_samples.xlsx” in Clark and others [2022]). Two SRMs, consisting of agricultural soils and representing low (SRM sample 2709, NIST San Joaquin soils; hereinafter referred to as “SRM sample 2709”) and high concentrations (SRM sample 2711, NIST Montana soil II; hereinafter referred to as “SRM sample 2711”) of trace elements, were analyzed to test recovery efficiency for a range of concentrations that are likely detected in the bed sediment in streams in the upper Clark Fork Basin. The digestion used for the bed-sediment samples was not a “total” digestion (does not liberate elements associated with crystalline lattices) but rather a “near-total” digestion. As a result, 100-percent recoveries were unlikely for elements strongly associated with crystalline phases. Recoveries of trace elements from SRMs processed by “near-total” digestion are useful to identify trace elements weakly or strongly associated with sediment particles and to evaluate the consistency of recoveries among multiple sets of analyses.

Mean recoveries of arsenic, cadmium, chromium, copper, iron, lead, manganese, nickel, and zinc from SRM sample 2709 ranged from 49.0 to 92.1 percent (“table_16_recovery_efficiency_bed_sediment_samples.xlsx” in Clark and others [2022]). Elements with low certified concentrations (cadmium, 0.371 microgram per gram) or with a strong association with crystalline lattices (arsenic, chromium, lead) had lower percentage recoveries (mean recoveries between 49.0 and 76.2 percent). Mean recoveries in SRM sample 2711 ranged from 70.6 to 97.0 percent (“table_16_recovery_efficiency_bed_sediment_samples.xlsx” in Clark and others [2022]). Because of the strong binding nature of chromium’s crystalline lattice, chromium had the lowest mean recovery at 70.6 percent. No adjustments were made to trace-element concentrations in bed-sediment samples based on recovery efficiencies.

Procedural blanks for bed-sediment samples consisted of the analysis of the same reagents used for sample digestion and reconstitution. Concentrated nitric acid used for sample digestion was heated and evaporated to dryness. After evaporation, 0.6 N hydrochloric acid was added to reconstitute the dry residue. Analytical results of procedural blanks for bed sediment were reported as a liquid-phase concentration, in micrograms per milliliter. A procedural blank was prepared and analyzed concurrently with bed-sediment samples for each site. Concentrations of trace elements in all procedural blanks were less than the MRL for all elements (“table_17_analyses_of_procedural_blanks_for_bed_sediment_samples.xlsx” in Clark and others [2022]). No adjustments were made to analytical data based on procedural blanks.

Tissue Concentrations

Tissue data for the long-term monitoring program in the Clark Fork Basin consist of trace-element concentrations in the whole-body tissue of aquatic benthic invertebrates. Invertebrate samples were collected once annually at the same 12 sites and on the same dates as bed-sediment samples (fig. 1; table 1), 8 of which have USGS streamgages, allowing for a direct comparison of biological data with bed-sediment and water-quality data. Travel restrictions during the 2020 pandemic limited field crew to two personnel and Clark Fork above Little Blackfoot River, near Garrison (12324400), and Clark Fork near Drummond (12331800) were not sampled because of the site access challenges. Biological data for water year 2020 are available in a USGS data release (“table_18_biological_data_Clark_Fork_Basin.xlsx” in Clark and others [2022]). Warm Springs Creek at Warm Springs (12323770) is sampled once every 3 years and was sampled during water year 2020.

Methods

Insect samples were collected using protocols described in Hornberger and others (1997). Benthic insects at immature stages were collected with a large nylon-mesh kick net. A single riffle at each site was sampled repeatedly until an adequate number of individual insects was collected to provide enough mass for analysis (for example, about 10 stoneflies, about 100–200 caddisflies). Two caddisfly species of the genus Hydropsyche (Hydropsyche cockerelli and Hydropsyche occidentalis) were targeted for collection in this study because they were detected at most sites. Hydropsyche species (spp.) that could not be positively identified were categorized as Hydropsyche spp. or Hydropsyche morosa group. On the few occasions when Hydropsyche were not present, other caddisflies, including Brachycentrus spp. and Rhyacophila spp., were collected. The caddisfly Arctopsyche grandis and the stoneflies Claassenia sabulosa and Hesperoperla spp. were collected where available to represent additional insect taxa commonly distributed in the Clark Fork Basin but with different feeding traits.

Samples of each taxon were sorted by genus in the field and placed in acid-washed plastic containers. Samples were frozen in a small amount of ambient stream water on dry ice within 30 minutes of collection. Between 1986 and 1998, macroinvertebrate containers were kept on ice to allow the insects to evacuate their gut contents (depurate) for 6 to 8 hours. Excess water was drained, and insects were frozen for transport to the laboratory. Since 1999, samples were immediately frozen on dry ice in the field to reduce the possibility of metal loss through intracellular breakdown during depuration. A comparison of immediately frozen to depurated samples indicated that although no substantial difference was detected for most metals, concentrations of copper were about 20 percent lower in the depurated samples than in the samples that were immediately frozen. The data were not adjusted for this difference.

Invertebrate samples were processed and analyzed at the USGS Metal Bioavailability Laboratory in Menlo Park, Calif. Insects were thawed, rinsed with ultrapure deionized water to remove particulate matter, and then sorted to their lowest possible taxonomic level. If large numbers of specimens were collected at a site, similar-sized individuals were composited into replicate subsamples. Subsamples were placed in tared scintillation vials and oven dried at 70 degrees Celsius. Subsamples were weighed to obtain a final dry weight and digested by reflux using concentrated nitric acid (Cain and others, 1992). After digestion, insect samples were evaporated to dryness on a hot plate. The dry residue was reconstituted in 0.6 N hydrochloric acid, filtered through a 0.45-µm pore-size filter, and analyzed undiluted using inductively coupled plasma-optical emission spectrometry (Hornberger and others, 1997) for arsenic, cadmium, chromium, copper, iron, lead, manganese, nickel, and zinc. The smallest concentration of a constituent that can be reliably reported for analyses of tissue samples is termed the MRL.

Results

Concentrations of trace elements in whole-body tissue of aquatic invertebrates collected during July 2020 are listed in the accompanying data release (“table_18_biological_data_Clark_Fork_Basin.xlsx” in Clark and others [2022]). The variability in the number of composite samples among species and among sites reflects differences in invertebrate abundance (the number of composite samples increases with the relative abundance of insects). Liquid-phase concentrations, in micrograms per milliliter, analyzed in the reconstituted samples were converted to solid-phase concentrations, in micrograms per gram, using equation 4 (used earlier in this report to calculate solid-phase concentrations of trace elements in bed sediment). All tissue samples were analyzed undiluted (the dilution ratio was 1 to 1).

Quality Assurance

The USGS protocols for field collection and processing of tissue samples were designed to prevent contamination from metal sources. Nonmetallic nets, sampling equipment, and processing equipment were used in all sample collection. Equipment was acid washed and rinsed in ultrapure deionized water before the first sample collection. Nets and equipment were thoroughly rinsed in stream water at each main-stem site. Clean nets were used at each tributary site.

Quality control of analytical results for tissue samples included laboratory-instrument calibration with standard solutions and analyses of quality-control samples designed to quantify precision and to identify the presence and magnitude of bias. Quality-control samples consisted of 12 replicates of the certified reference material (CRM) TORT–3 (lobster hepatopancreas), which was purchased from the National Research Council Canada. Quality-control samples were analyzed in a proportion equivalent to about 20 percent of the total number of biota samples. Recovery efficiencies for trace-element analyses of the TORT–3 CRM are listed in the accompanying data release (“table_19_recovery_efficiency_biota_samples.xlsx” in Clark and others [2022]). Data-quality objectives have not been established for analytical recovery in tissue, but percentage recoveries indicate analytical performance.

Mean CRM recoveries for TORT–3 ranged from 85.3 to 112 percent for all constituents except lead (“table_19_recovery_efficiency_biota_samples.xlsx” in Clark and others [2022]). The mean recoveries were within 15 percent of theoretical 100-percent recovery (based on the 95-percent confidence interval) for arsenic, cadmium, chromium, copper, iron, manganese, nickel, and zinc. Recoveries for lead were less than the analytical detection limits for the CRM and are not reported here; however, lead concentrations in the tissue samples are much greater than detection limits and are included in the accompanying USGS data release (Clark and others, 2022). No adjustments were made to trace-element concentrations in biota samples based on recovery efficiencies.

Procedural blanks for biota consisted of undiluted aliquots of the same reagents used to digest and reconstitute tissue of aquatic insects. Analytical results of procedural blanks for biota (“table_20_analyses_of_procedural_blanks_for_tissue_samples.xlsx” in Clark and others [2022]) were reported as a liquid-phase concentration, in micrograms per milliliter. A procedural blank was prepared and analyzed concurrently with tissue samples for each site. Concentrations of trace elements in all procedural blanks were less than the MRL; therefore, no adjustments to the data were necessary.

Statistical Summaries of Data

Statistical summaries of long-term water-quality, bed-sediment, and biological data for the Clark Fork Basin are provided in an accompanying USGS data release (“table_21_long_term_water_quality_data_Clark_Fork_Basin.xlsx,” “table_22_long_term_bed_sediment_data_Clark_Fork_Basin.xlsx,” and “table_23_long_term_biological_data_Clark_Fork_Basin.xlsx” in Clark and others [2022]) for the period of record at each site. The summaries include the period of record; number of samples; and maximum, minimum, mean, and median concentrations.

The summaries do not include data for supplemental samples collected at selected sites that targeted high-flow conditions or maintenance drawdowns of Milltown Reservoir, which might disproportionately skew the long-term statistics relative to the other sites in the network. Sample results at sites that have been sampled for other projects can be accessed in the USGS NWIS database at https://waterdata.usgs.gov/mt/nwis/qw. Statistical summaries of long-term bed-sediment and biological data are based on results of samples collected once each year during the indicated years. Because not all sites were sampled for bed sediment and biota every year, the data for some sites do not represent a consecutive annual record. Statistical summaries are not presented for discontinued sites.

Statistics for long-term bed-sediment data (“table_22_long_term_bed_sediment_data_Clark_Fork_Basin.xlsx” in Clark and others [2022]) were based on the mean trace-element concentrations determined for each year from the mean of the analyses of composite samples; therefore, the number of samples for bed sediment represents the number of years that the constituent was analyzed. The number of arsenic samples for bed sediment was smaller than the number for other trace elements because sampling for arsenic began in September 2003. In addition, the number of samples analyzed for silver in bed sediment is smaller because analysis for this constituent was discontinued in 2004.

In contrast, statistics for long-term tissue data (“table_23_long_term_biological_data_Clark_Fork_Basin.xlsx” in Clark and others [2022]) were based on individual analyses for each composite sample collected rather than on a single mean concentration for each year. Differences in the number of composited invertebrate samples among species reflected differences in species abundance within and between sites and among years. As a result, the statistics for insects described a wider range of variation in trace-element concentrations than would be evident if results from individual composite samples were averaged. Also, the number of samples for arsenic in insect samples is smaller than the number for other trace elements because sampling for arsenic began in September 2003. The abundance of aquatic invertebrates at a site in a given year limits the biomass of the sample, which in turn, may result in varied MRLs. When MRLs vary among years, differences in concentration with time are difficult to determine, especially when a large percentage of the samples have concentrations less than MRLs.

The presence or absence of insect species at a given site can vary among years and may result in different taxa being analyzed in the long-term period of record. Because Hydropsyche insects were not sorted to the species level during 1986–89, statistics for sites sampled during those years were based on the results of all Hydropsyche species combined. At some sites, statistics for the Hydropsyche morosa group were based on the combined results for two or more species because these samples could not clearly be identified to the species level, but the individual insects had morosa characteristics.

References Cited

American Public Health Association, 1998, Standard methods for the examination of water and wastewater (20th ed.): Washington, D.C., American Public Health Association, American Water Works Association, and Water Environment Federation, 1,325 p.

Anderson, C.W., 2005, Turbidity (ver. 2.1): U.S. Geological Survey Techniques of Water-Resources Investigations, book 9, chap. A6, sec. 6.7, 64 p. [Also available at https://doi.org/10.3133/twri09A6.7.]

Andrews, E.D., 1987, Longitudinal dispersion of trace metals in the Clark Fork River, Montana, in Averett, R.C., and McKnight, D.M., eds., Chemical quality of water and the hydrologic cycle: Chelsea, Mich., Lewis Publishers, p. 179–191.

Axtmann, E.V., Cain, D.J., and Luoma, S.N., 1997, Effect of tributary inflows on the distribution of trace metals in fine-grained sediment and benthic insects of the Clark Fork River, Montana: Environmental Science & Technology, v. 31, no. 3, p. 750–758. [Also available at https://doi.org/10.1021/es960417y.]

Axtmann, E.V., and Luoma, S.N., 1991, Large scale distribution of metal contamination in the fine-grained sediment of the Clark Fork River, Montana: Applied Geochemistry, v. 6, no. 1, p. 75–88. [Also available at https://doi.org/10.1016/0883-2927(91)90064-V.]

Brenton, R.W., and Arnett, T.L., 1993, Methods of analysis by the U.S. Geological Survey National Water Quality Laboratory—Determination of dissolved organic carbon by UV-promoted persulfate oxidation and infrared spectrometry: U.S. Geological Survey Open-File Report 92–480, 12 p. [Also available at https://doi.org/10.3133/ofr92480.]

Cain, D.J., Luoma, S.N., and Axtmann, E.V., 1995, Influence of gut content in immature aquatic insects on assessments of environmental metal contamination: Canadian Journal of Fisheries and Aquatic Sciences, v. 52, no. 12, p. 2736–2746. [Also available at https://doi.org/10.1139/f95-862.]

Cain, D.J., Luoma, S.N., Carter, J.L., and Fend, S.V., 1992, Aquatic insects as bioindicators of trace element contamination in cobble-bottom rivers and streams: Canadian Journal of Fisheries and Aquatic Sciences, v. 49, no. 10, p. 2141–2154. [Also available at https://doi.org/10.1139/f92-237.]

Cannon, M.R., and Johnson, D.R., 2004, Estimated water use in Montana in 2000: U.S. Geological Survey Scientific Investigations Report 2004–5223, 61 p. [Also available at https://doi.org/10.3133/sir20045223.]

Clark, G.D., Hornberger, M.I., Hepler, E.J., Cleasby, T.E., and Heinert, T.L., 2021, Water-quality, bed-sediment, and invertebrate tissue trace-element concentrations for tributaries in the Clark Fork Basin, Montana, October 2018–September 2019: U.S. Geological Survey Open-File Report 2021–1027, 16 p., accessed November 16, 2021, at https://doi.org/10.3133/ofr20211027.

Clark, G.D., Hornberger, M.I., Hepler, E.J., and Heinert, T.L., 2022, Results of water-quality, bed-sediment, and invertebrate tissue trace-element concentrations for tributaries in the Clark Fork Basin, Montana, October 2019–September 2020: U.S. Geological Survey data release, https://doi.org/10.5066/P93BP9P8.

Davis, B.E., and the Federal Interagency Sedimentation Project, 2005, A guide to the proper selection and use of federally approved sediment and water-quality samplers: U.S. Geological Survey Open-File Report 2005–1087, 20 p. [Also available at https://doi.org/10.3133/ofr20051087.]

Dodge, K.A., and Lambing, J.H., 2006, Quality-assurance plan for the analysis of suspended sediment by the U.S. Geological Survey in Montana: U.S. Geological Survey Open-File Report 2006–1242, 25 p. [Also available at https://doi.org/10.3133/ofr20061242.]

Driscoll, A., and Hatcher, J., 2010, Chain of custody: U.S. Geological Survey National Water Quality Laboratory, SOP QUAX0030.4, effective May 5, 2010, 17 p.

Edwards, T.K., and Glysson, G.D., 1999, Field methods for measurement of fluvial sediment: U.S. Geological Survey Techniques of Water-Resources Investigations, book 3, chap. C2, 89 p. [Also available at https://doi.org/10.3133/twri03C2.]

Fishman, M.J., ed., 1993, Methods of analysis by the U.S. Geological Survey National Water Quality Laboratory—Determination of inorganic and organic constituents in water and fluvial sediments: U.S. Geological Survey Open-File Report 93–125, 217 p. [Also available at https://doi.org/10.3133/ofr93125.]

Fishman, M.J., and Friedman, L.C., 1989, Methods for determination of inorganic substances in water and fluvial sediments: U.S. Geological Survey Techniques of Water-Resources Investigations, book 5, chap. A1, 545 p. [Also available at https://doi.org/10.3133/twri05A1.]

Friedman, L.C., and Erdmann, D.E., 1982, Quality assurance practices for the chemical and biological analyses of water and fluvial sediments: U.S. Geological Survey Techniques of Water-Resources Investigations, book 5, chap. A6, 181 p. [Also available at https://pubs.usgs.gov/twri/twri5a6/.]

Garbarino, J.R., Kanagy, L.K., and Cree, M.E., 2006, Determination of elements in natural-water, biota, sediment, and soil samples using collision/reaction cell inductively coupled plasma-mass spectrometry: U.S. Geological Survey Techniques and Methods, book 5, chap. B1, 88 p. [Also available at https://doi.org/10.3133/tm5B1.]

Garbarino, J.R., and Struzeski, T.M., 1998, Methods of analysis by the U.S. Geological Survey National Water Quality Laboratory—Determination of elements in whole-water digests using inductively coupled plasma-optical emission spectrometry and inductively coupled plasma-mass spectrometry: U.S. Geological Survey Open-File Report 98–165, 101 p. [Also available at https://doi.org/10.3133/ofr98165.]

Guy, H.P., 1969, Laboratory theory and methods for sediment analysis: U.S. Geological Survey Techniques of Water-Resources Investigations, book 5, chap. C1, 58 p. [Also available at https://doi.org/10.3133/twri05C1.]

Hoffman, G.L., Fishman, M.J., and Garbarino, J.R., 1996, Methods of analysis by the U.S. Geological Survey National Water Quality Laboratory—In-bottle acid digestion of whole-water samples: U.S. Geological Survey Open-File Report 96–225, 28 p. [Also available at https://doi.org/10.3133/ofr96225.]

Hornberger, M.I., Lambing, J.H., Luoma, S.N., and Axtmann, E.V., 1997, Spatial and temporal trends of trace metals in surface water, bed sediment, and biota of the upper Clark Fork Basin, Montana, 1985–95: U.S. Geological Survey Open-File Report 97–669, 84 p. [Also available at https://doi.org/10.3133/ofr97669.]

Horowitz, A.J., Demas, C.R., Fitzgerald, K.K., Miller, T.L., and Rickert, D.A., 1994, U.S. Geological Survey protocol for the collection and processing of surface-water samples for the subsequent determination of inorganic constituents in filtered water: U.S. Geological Survey Open-File Report 94–539, 57 p. [Also available at https://doi.org/10.3133/ofr94539.]

Jones, B.E., 1987, Quality control manual of the U.S. Geological Survey’s National Water Quality Laboratory: U.S. Geological Survey Open-File Report 87–457, 17 p. [Also available at https://doi.org/10.3133/ofr87457.]

Lambing, J.H., comp., 2006, Quality-assurance plan for water-quality activities of the U.S. Geological Survey Montana Water Science Center: U.S. Geological Survey Open-File Report 2006–1275, 39 p. [Also available at https://doi.org/10.3133/ofr20061275.]

Luoma, S.N., and Bryan, G.W., 1981, A statistical assessment of the form of trace metals in oxidized estuarine sediments employing chemical extractants: Science of the Total Environment, v. 17, no. 2, p. 165–196. [Also available at https://doi.org/10.1016/0048-9697(81)90182-0.]

Patton, C.J., and Kryskalla, J.R., 2011, Colorimetric determination of nitrate plus nitrite in water by enzymatic reduction, automated discrete analyzer methods: U.S. Geological Survey Techniques and Methods, book 5, chap. B8, 34 p. [Also available at https://doi.org/10.3133/tm5B8.]

Pritt, J.W., and Raese, J.W., eds., 1995, Quality assurance/quality control manual—National Water Quality Laboratory: U.S. Geological Survey Open-File Report 95–443, 35 p. [Also available at https://doi.org/10.3133/ofr95443.]

Rantz, S.E., and others, 1982, Measurement and computation of streamflow: U.S. Geological Survey Water-Supply Paper 2175, v. 2, 631 p., accessed November 2019 at https://doi.org/10.3133/wsp2175.

Taylor, J.K., 1987, Quality assurance of chemical measurements: Chelsea, Mich., Lewis Publishers, 328 p.

U.S. Geological Survey, 2022, USGS water-quality data for the Nation, in USGS water data for the Nation: U.S. Geological Survey National Water Information System database, accessed January 2022, at https://doi.org/10.5066/F7P55KJN. [Water-quality data directly accessible at https://waterdata.usgs.gov/nwis/qw.]

U.S. Geological Survey, [variously dated], National field manual for the collection of water-quality data: U.S. Geological Survey Techniques of Water-Resources Investigations, book 9, chaps. A1–A10, accessed November 2019 at https://water.usgs.gov/owq/FieldManual/.

Wagner, R.J., Boulger, R.W., Jr., Oblinger, C.J., and Smith, B.A., 2006, Guidelines and standard procedures for continuous water-quality monitors—Station operation, record computation, and data reporting: U.S. Geological Survey Techniques and Methods, book 1, chap. D3, 51 p., 8 attachments. [Also available at https://doi.org/10.3133/tm1D3.]

Ward, J.R., and Harr, C.A., eds., 1990, Methods for collection and processing of surface-water and bed-material samples for physical and chemical analyses, U.S. Geological Survey Open-File Report 90–140, 71 p. [Also available at https://doi.org/10.3133/ofr90140.]

Conversion Factors

U.S. customary units to International System of Units

Multiply By To obtain
mile (mi) 1.609 kilometer (km)
square mile (mi2) 2.590 square kilometer (km2)

Temperature in degrees Celsius (°C) may be converted to degrees Fahrenheit (°F) as follows:

°F=(1.8×°C)+32

Datum

Horizontal coordinate information is referenced to the North American Datum of 1927 (NAD 27).

Supplemental Information

Concentrations of chemical constituents in water are given either in milligrams per liter (mg/L) or micrograms per liter (µg/L).

Pore size is given in micrometers (µm).

Suspended-sediment sizes are given in millimeters (mm) and weights are given in grams (g) of dried material.

Bottle capacities or liquid measurements are given in milliliters (mL).

Liquid-phase trace-element concentrations are given in micrograms per milliliter (µg/mL), and solid-phase concentrations are given in micrograms per gram (µg/g).

A water year is the 12-month period from October 1 through September 30 and is designated by the calendar year in which it ends. For example, water year 2020 is the period from October 1, 2019, through September 30, 2020.

Abbreviations

CRM

certified reference material

ICP–MS

inductively coupled plasma-mass spectrometry

LRL

laboratory reporting level

MRL

minimum reporting level

NIST

National Institute of Standards and Technology

NWIS

National Water Information System

NWQL

National Water Quality Laboratory

RSD

relative standard deviation

spp.

species

SRM

standard reference material

USGS

U.S. Geological Survey

For more information about this publication, contact:

Director, USGS Wyoming-Montana Water Science Center

3162 Bozeman Avenue

Helena, MT 59601

406–457–5900

For additional information, visit: https://www.usgs.gov/centers/wy-mt-water/

Publishing support provided by the

Rolla Publishing Service Center

Disclaimers

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

Clark, G.D., Hornberger, M.I., Hepler, E.J., and Heinert, T.L., 2022, Water-quality, bed-sediment, and invertebrate tissue trace-element concentrations for tributaries in the Clark Fork Basin, Montana, October 2019–September 2020: U.S. Geological Survey Open-File Report 2022–1090, 17 p., https://doi.org/10.3133/ofr20221090.

ISSN: 2331-1258 (online)

Study Area

Publication type Report
Publication Subtype USGS Numbered Series
Title Water-quality, bed-sediment, and invertebrate tissue trace-element concentrations for tributaries in the Clark Fork Basin, Montana, October 2019–September 2020
Series title Open-File Report
Series number 2022-1090
DOI 10.3133/ofr20221090
Year Published 2022
Language English
Publisher U.S. Geological Survey
Publisher location Reston, VA
Contributing office(s) Wyoming-Montana Water Science Center
Description Report: vii, 17 p.; Data Release; Dataset
Country United States
State Montana
Other Geospatial Clark Fork Basin
Online Only (Y/N) Y
Google Analytic Metrics Metrics page
Additional publication details