Abstract

Long-term mining activities in the mountains around Creede, Colorado have resulted in significant contamination in soils and water in the Willow Creek floodplain. Total major and trace were determined for soils and water and sequential chemical extraction for soils. Objectives were to determine concentrations and potential reactivity of trace elements and investigate their relationship with other soil and water properties. Water trace elements showed significant variability among sites, ranging from 347 to 12108 μg/L. Relative trend showed (Zn > Sr > Ba) > (Mn > W > Cd) > (Sn > V Ni Cu > Co) > (Ag). Soil trace elements showed significant short-range spatial variability, ranging from 2819 to 19274 mg/kg. Relative trend showed (Pb Zn > Mn > Ba > P) > (As > Cu > Sr > V > Cd > Sb Ag) > (Co Cr > Mo Sn Ni) > (Be W > Se Hg). Predominant fractions were oxide, specifically-sorbed/carbonate bound, and residual. Water soluble and exchangeable fractions showed (Zn Cd) > Pb and Cd > Zn > Pb, respectively. Mobility factors for highly contaminated soils showed Cd Zn > Pb > Cu > As.

1. Introduction

The measurement of the total extractable pool of trace elements has been commonly used to assess the environmental levels or background amounts of trace elements in soils and water [17]. Selective fractionations or “geochemical partitioning” have been used to evaluate the potential reactivity of these elements [811] and are widely used in soil pollution studies, providing qualitative evidence about trace element reactivity and indirectly of their bioavailability [1215]. Trace element bioavailability in soils is a complex phenomenon, affected by many factors such as total concentration, pH, organic matter, clay, and redox conditions [7, 1620]. Trace element toxicity in waters is strongly affected by site-specific water quality factors such as pH, hardness, and other dissolved constituents [21, 22]. An assessment of trace element fate, bioavailability, and transport (e.g., surface and groundwater) is required in order to predict potential contamination and impact upon soil and water quality.

Historic mining activities in the mountains around Creede, Colorado, began around 1889 and continued until 1985. Underground mining of silver and base metals has resulted in Zn and Cd contamination of ground and surface water in and along the broad floodplain of Willow Creek below Creede [23]. Willow Creek, a tributary to the Rio Grande, is polluted from drainage from various mine adits and rock piles upstream of Creede and by leachates from a gravel-capped tailings pile below [23]. Ores common to this area contain Cd which occurs mainly in the Zn sulfides sphalerite and wurtzite and is recovered with Zn usually from polymetallic ores containing Pb and Cu [24].

Major and trace element concentrations (Al, Cd, Cu, Fe, Mn, Pb, and Zn) in river systems are of interest because of their many natural and mine-related sources and their impact on aquatic life through a variety of mechanisms, including both acute and chronic toxic effects of aqueous metals as well as from metal exposure in soils and sediments [22, 23, 2537].

To assess, monitor, and revitalize the Willow Creek watershed, the Willow Creek Reclamation Committee (WCRC) partnered with numerous state and Federal agencies, including the US Environmental Protection Agency (USEPA); US Forest Service (USFS); US Geological Service (USGS); Colorado Department of Public Health and Environment (CDPHE); Colorado Division of Minerals and Geology; the US Department of Agriculture, Natural Resources Conservation Service (USDA-NRCS), and others. Assessment reports produced by WCRC and/or its partners include [22, 3843].

The USDA-NRCS designed a sampling scheme to help assess contamination of Willow Creek from nonpoint sources of metals from historical mining activities. In August 2005, USDA-NRCS collected and analyzed soil and water samples from Willow Creek. Initial intent of sampling was to provide background survey levels to help guide potential future samplings for more intensive investigations and later expanded to provide comment on potential soil metal/metalloid bioavailability. This paper reports the results of this investigation. The objectives of this study were as follows: (1) to determine the concentrations and potential reactivity of trace elements using total analysis and sequential chemical fractionation so as to (a) determine the extent and variability of trace element dispersal along the creek and (b) investigate the relationship between trace element concentrations and their potential reactivity with other soil and water properties as well as vegetation or lack thereof along the stream. Knowledge and understanding of these relationships are important to researchers when evaluating trace element data for use in predictive models for environmental purposes. In addition to the geochemistry, detailed site descriptions and other basic soil and water characterization were considered necessary to aid in the interpretation and application of these elemental data.

2. Materials and Methods

2.1. Study Area

The Willow Creek watershed is an area of 103 km2 in Mineral County, Colorado in the northern San Juan Mountains. Willow Creek is formed by the confluence of the East and West Willow Creeks and is a tributary of the Rio Grande River. The City of Creede (Figure 1) and most of the Creede Mining District are located within the watershed. The Creede Mining District is in the central part of the San Juan Volcanic Field. Principal metals produced in the Creede District included Ag, Pb, Zn, Cu, and Au. Production of these metals was from veins along faults that formed during subsidence of the Creede caldera.

Bedrock geology is dominated by igneous rocks related to extensive tertiary volcanism, with ashflow tuffs the common rock in the watershed above Creede and the dominant rock for mineral deposits in the mining areas along the East and West Willow Creeks [41]. The area beneath Creede and southward to the confluence of Willow Creek and the Rio Grande is predominately the Creede Formation which is mostly lake and river deposits of reworked ash, sand and gravel, and travertine deposits from numerous mineral hot springs forming after the eruption of the Creede Caldera [41]. Much of the older rocks in the lower part of the watershed, especially along the waterways, are covered by Quaternary alluvium from Willow Creek and the Rio Grande [44].

The Willow Creek watershed has a 1.6 km vertical relief and is characterized by steep slopes and distinct vegetative communities. Due to the substantial elevation variability, the climate within the watershed is diverse, affecting the vegetative communities, stream flow magnitude and timing, water temperature, ground-water recharge, and other watershed characteristics [41]. The Willow Creek watershed crosses four major life zones from the high mountains to the Rio Grande Valley, with the origin of Willow Creek in the Alpine Zone (tufted hairgrass-sedge and willow), in the sub-Alpine Zone (Engleman spruce, sub-Alpine fir, Thurber fescue grassland), in Montane Zone (aspen, Douglas Fir, and bristlecone pine), and in the Foothills Zone (Arizona fescue grasslands) [41]. In Creede, the annual average precipitation, temperature, and snowfall are 335 mm, 16°C, and 122 cm, respectively, with the highest and lowest precipitation in August and December and the highest and lowest temperature in July and January, respectively [41]. The largest volume of snow or rain is at the highest elevations. Aquatic resources within the watershed are primarily streams but also include wetlands and riparian areas.

2.2. Sampling Strategy: Soils and Water

Soils are naturally variable, their properties changing horizontally across the landscape as well as vertically down the soil profile. Soil map units derived from changes in topography, underlying geology, and dominant vegetation type can be used for horizontal subdivisions, with soil horizons as subdivisions of vertical change [45]. Environmental and/or anthropogenic disturbances introduce additional variation into natural landscapes. In the field, the US soil survey has routinely used pedon sampling to obtain representative samples [46, 47]. In environmental studies as with soil survey, intuitive sampling often forms the basis of exploratory sampling which may be used for qualitative assessment of soils where impact or damage is visible or anticipated [45, 47, 48]. Pedon sampling has been used in studies when the area affected is known and/or usually visible, and the types of contaminants are either known or unknown [11, 4951]. Replicate sampling is often included in sampling plans for water quality, the primary purpose of which is to identify and/or quantify the variability in all or part of the sampling and analysis system [52]. Amount and composition of water samples vary strongly with small changes in location along water body [53]. The primary objective of water sampling is the same as with soil sampling, that is, to obtain a representative sample. For the reasons stated above, the approach to the USDA-NRCS exploratory soil/water sampling at Willow Creek (Aug 9–11, 2005) utilized pedon sampling in the collection of soil samples and water samples replicated on different day (sequential) and same day (concurrent) for laboratory characterization.

All water and soil samples were collected from sites located downstream of Creede on the Willow Creek floodplain. Water sampling was performed in a unidirectional manner with the water flow. Willow Creek has been substantially diverged from its presettlement form due to mining activities [43]. At the time of sampling, Willow Creek was in a braided form. Water sampling was conducted during low-flow conditions when the metal-loading pattern not only shows the metal sources entering the streams on a continual basis but also those sources contributing to high concentrations during the low-flow winter months when mine drainage is less diluted by other water sources and the most toxic conditions likely occur [22, 35, 38]. It was observed that flow conditions in some areas of Willow Creek were free flowing, whereas in other areas, flow was physically obstructed causing water to pool and stagnate, suggesting stream waters had both oxidizing and reducing conditions at the time of sampling.

Water samples were collected from seven sites over a three-day period. Two sites (Westfork and Eastfork) were collected at higher elevations than other sites and are referred to as “baseline” or “reference” samples. Some water samples were collected in conjunction with the pedon sampling [47] and as such bear the same name as the soil. Table 1 provides the names for the water samples and their soil survey identification numbers, with embedded codes indicating concurrent and/or sequential sampling. For example Creek water shows the soil survey numbers S05CO-079-003A-1 through S05CO-079-003C-3, with the A through C representative of sequential samples and 1 through 3 representative of concurrent samples.

Study soils have aquic conditions, undergoing periods of saturation and reduction [54]. Sampling in August during low-flow conditions provided a good representation of more oxidizing conditions in study soils compared to other times of the year with high-flow rates (high precipitation or snowmelt). Soil surface materials were collected from ten sites, including two sites (Spruce and Railroad), where additional underlying materials were collected. Figure 2 shows sites of these sampled soils. Table 2 provides longitude and latitude for sites and also shows soil name, soil survey identification number, and sampling depth (cm).

2.3. Field and Laboratory Characterization: Soils and Water
2.3.1. Soil and Water Sampling

Soils were described, sampled, and analyzed using standard USDA-NRCS soil survey methods [47, 55]. Water samples were collected using standard USDA-NRCS [47] and USGS [5658] sampling protocols, expedited for transport to the laboratory, and analyzed using standard USDA-NRCS soil survey laboratory methods [47].

2.3.2. Soil Analyses

Bulk soil samples were air-dried, homogenized, and sieved to <2 mm in the laboratory to remove rock fragments. Samples received standard laboratory characterization performed on air-dried <2 mm soil with resulting data reported on an oven-dry basis. Particle-size analysis was determined by sieve and pipette methods, following pretreatments to remove organic matter (OM) and soluble salts, and chemical dispersion with sodium hexametaphosphate. Total C, N, and S were determined by dry combustion. Soil pH 1 : 1 soil:water was determined. Exchangeable cations (Ca, Mg, K, and Na) were extracted with buffered (pH 7.0) NH4OAc and measured by atomic absorption spectrometry (AAS). Cation exchange capacity (CEC-7) was determined by Kjeldahl titration of the NH4-saturated soil. Base saturation by CEC-7 was determined by dividing the sum of NH4OAc extractable bases by CEC-7 and multiplying by 100. Acidity was extracted by BaCl2-TEA, pH 8.2. Exchangeable Al was extracted with 1 M KCl and determined by inductively couple plasma emission spectrometry (ICP-AES) on samples with pH . If electrical conductivity (EC) for 1 : 2 soil water extract was  dS/m, a saturated paste extract was prepared and analyzed for cations (Ca, Mg, Na, and K) by AAS and anions (F, Cl, PO4, Br, OAc, SO4, NO2, and NO3) by ion chromatography and carbonate and bicarbonate (CO3 and HCO3) by acid titration. Dithionite-citrate extracts were analyzed for Fe, Al, and Mn (Fed, Ald, and Mnd, resp.) by AAS. Sodium pyrophosphate extracts were analyzed for Fe, Al, Si, and Mn (Fep, Alp, Sip, and Mnp, resp.) by AAS. Acid oxalate extracts were analyzed for Fe, Al, and P (Feo Alo, and Po, resp.) by ICP-AES. Optical grain counts were conducted on the coarse-silt (0.02–0.05 mm), very fine sand (0.05–0.10 mm), or fine-sand (0.10–0.25 mm) fraction. Crystalline clay minerals (<2 μm) were identified by X-ray diffraction analysis.

2.3.3. Soil Elemental Analyses

Soil samples for major and trace element analysis were prepared metal-free by grinding approximately 10 g of <2 mm soil in a silicon nitride ball mill to <74 μm. Major elements (Al, Ca, Fe, K, Mg, Mn, Na, P, Si, Sr, Ti, and Zr) were determined on soil samples (0.25 g) by microwave digestion (180°C for 9.5 min) combined with 4 mL HF + 9 mL HNO3 + 3 mL HCl. Following digestion, boric acid was added to neutralize HF (2.5% in 100 mL final volume) and extracts analyzed by ICP-AES. Total extractable (TE) trace elements (Ag, As, Ba, Be, Cd, Cr, Co, Cu, Hg, Mn, Mo, Ni, P, Pb, Sb, Se, Sn, Sr, V, W, and Zn) were determined on soil samples (0.5 g) by microwave digestion (180°C for 7 min) using aqua regia (9 mL HNO3 + 3 mL HCl). Elements extracted by this method are designated as TE, and the sum of all elements extracted by this method are designated SUMTE. In order to simplify presentation of study results, Mn, P, and Sr are included as trace elements.

2.3.4. Soil Sequential Fractionations

Sequential chemical fractionations were based on the procedure of Tessier et al. [10] and Keller and Vedy [8]. Sequential fractionations were determined on separate samples (0.1 g) in 50 mL centrifuge tubes. Extracts were replicated, with mean values reported. Residues were separated from supernatants by 80-min centrifugation at 4000 rpm. Soils were extracted twice on fractions 3 and 4 with designated reagents (10 mL), and both extracts combined. Washing steps (10 mL H2O, shake 15 min) were performed between fractionations 2 through 5, with resulting supernatants separated from residues by 80-min centrifugation and solutions combined with initial extracts. Concentrated (16 M) HNO3 was added to fractions 1 and 2 to acidify final extracts. Trace element concentrations in all extracts were analyzed by inductively coupled plasma mass spectrometry (ICP-MS). Percent recovery was determined by summing the trace elements from the sequential-extraction process (SUMSEQ), divided by the amount determined by the TE elemental analysis by microwave digestion (SUMTE), and multiplied by 100.

Total extractable (TE) and sequential fractions 1 through 6 are designated herein as water soluble (WS), exchangeable (EX), specifically sorbed/carbonate bound (SS/CAR), oxide bound (OX), organic/sulfide bound (OM/S), and residual (RES). The sums of total extractable and sequential fractions are designated herein as SUMTE and SUMSEQ, respectively (Table 3). This designated chemical partitioning of trace elements is based on the evaluations of the sequential fractionation method used herein as applied in various studies [6, 812, 59]. In this study, the CH3COONa/CH3COOH reagent is used for both calcareous and noncalcareous soils and is designated as SS/CAR. Tessier et al. [10] indicated that the pH 5.0 used in the SS/CAR fraction solubilizes not only the carbonate-bound fraction if present but also a large fraction of specifically adsorbed trace elements from various solid substrates. Similarly, the H2O2 extraction (fraction 5) is designated as OM/S, recognizing that not only OM but also sulfide minerals may be extracted to a large extent [10].

2.3.5. Water Analyses

Water samples were analyzed for cations (Ca, Mg, Na, and K) by AAS and anions (F, Cl, PO4, Br, OAc, SO4, NO2, and NO3) by ion chromatography and carbonate and bicarbonate (CO3 and HCO3) by acid titration. Electrical conductivity and pH were determined on water samples. Water samples were filtered (0.45 μm), acidified with HCl, and analyzed for major (Al, Ca, Fe, K, Mg, and Sr) and trace elements (Ag, As, Ba, Be, Cd, Co, Cr, Cu, Mn, Mo, Ni, P, Pb, Sb, Sn, Sr, V, W, and Zn) by ICP-AES. For purposes of this study, Sr is included as trace element.

2.3.6. Quality Control

Quality control (QC) included blanks, duplicates, certified reference material (CRM), and standards from the National Institute of Standards and Technology (NIST). Reference material for soil samples included CRM Loam Soil C, high-purity standards, Charleston, SC and NIST 2709a San Joaquin Soil. The high purity standard CRM-TMDW was used for water samples. Elemental recoveries and relative standard deviation (RSD) for soil and water QC standards ranged from 90% to 115% and ≤11%, respectively. Percent recovery for QC soil standards, determined by dividing SUMSEQ by SUMTE and multiplying by 100, ranged from 124 to 130 percent.

3. Results and Discussion

3.1. Water: Chemical Properties

Water samples were less acidic than most of the sampled soils, with pH ranging from 6.5 to 7.5 (slightly acid to slightly alkaline) (Table 3). Of the sampled waters, the Channel and Creek water samples were generally less acidic than those waters sampled in conjunction with the soils (Gravelbar, Cottonwood, and Spruce). The relatively lower pH for the water associated with the sampled soils and for the water from the second day sampling of the Creek water may be due in part to the capturing of more sediment in the sample compared to the other samplings. In general, there are relatively small differences in water pH within sites, that is, same-day samples (concurrent) and different-day samples (sequential). For example, Channel water showed a mean and standard deviation of over a three-day period. These data suggest minor daily fluctuations in the stream or channel pH may be anticipated under relatively similar environmental conditions. This lack of high daily variability within the time frame of sampling may also indicate that the water system has a reasonably good buffering system. Poorly buffered waters and high daily variability in pH are often attributed to biological processes that affect the carbonate buffering system.

The water pH for the two locations at the higher elevations designated as “reference” or “baseline” sites (Westfork and Eastfork) were 7.3 and 7.4, respectively. These data suggest the pH of the water at the upper elevations was not dramatically different from the water at the lower elevations (excluding those associated with sampled soils) during the time frame of sampling. As water moves through soils and rocks, the pH may increase or decrease as additional chemical reactions occur. The pH within streams can impact toxic materials, with high acidity or alkalinity tending to convert insoluble metal sulfides to soluble forms, thereby increasing the concentrations of toxic metals.

Other data [43] showed significant pH variability in groundwater across the Willow Creek floodplain, ranging from 2.6 to 7.0. Yochum et al. [43] suggests pH can be a good indicator for groundwater contamination, that is, the link between low pH and groundwater contamination can be used for early detection for ground water contamination plume movement, minimizing the expense of water quality analyses.

Electrical conductivity (EC) showed variability among sites, with the lowest EC values in the “reference” sites, with 0.06 and 0.11 dS/m in the Eastfork and Westfork waters, respectively (Table 3). Highest EC values were found in those waters sampled in conjunction with the soils, ranging from 0.16 to 0.22 dS/m. Data showed good agreement within sites, with the Channel and Creek waters showing the greatest consistency with means and standard deviations of and  dS/m, respectively, over a three-day period. Similar to pH, these data suggest minor fluctuations in the EC of the stream or channels may be anticipated under relatively similar environmental conditions. Additionally, these data may indicate that the salt load may increase in transport by stream waters from higher elevations to the valley floor and that these salts would be more significant in areas where water stagnates or pools with soils and sediments.

Predominant ions in water samples were SO4 > Ca > HCO3 (Table 3). Ion concentrations showed good agreement across all samples (among and within sites), with means and standard deviations of , , and  mmol/L for SO4, Ca, and HCO3, respectively. Sulfate is a product of pyrite oxidation as well as the oxidation and dissolution of other metal-bearing sulfide minerals. Carbonate buffering results from chemical equilibrium between calcium, carbonate, bicarbonate (HCO3), carbon dioxide, and hydrogen ions in the water and carbon dioxide in the atmosphere. Buffering causes waters to resist changes in pH [60]. The amount of buffering is primarily determined by carbonate and bicarbonate concentration which are introduced into the water from dissolved calcium carbonate (i.e., limestone) and similar minerals in the watershed [61]. Glucose is the primary product of photosynthesis, with acetate (OAc) forming from the glucose as plants decompose. In comparison, Kimball et al. [38] found the stream water at the injection point in West Willow Creek a CaHCO3-CaSO4 type water, reflecting the chemical weathering of bedrock in the watershed [44]. Kimball et al. [38] further determined these ions and pH remained constant until the inflow of the Nelson Tunnel, at which point the chemical character of water changed to a CaSO4 type, and pH was lower. The Nelson Tunnel discharge is located on West Willow Creek about 2.4 km upstream of Creede and is considered responsible for 75% of the total inorganics (heavy metals) contamination in Willow Creek [42]. With the inflow of East willow creek, pH became more basic again, Ca and SO4 diluted, and alkalinity remained nearly constant [38]. For this study, soil and water samples were collected below Creede on the southern part of the Willow creek floodplain.

3.2. Water: Major and Trace Elements

Metals in natural waters undergo continual changes between dissolved, precipitated, and sorbed-to-sediment forms, with the rate of these adsorption, desorption, and precipitation processes depending on pH, redox potential, water chemistry, and composition of bottom and suspended sediments [62]. These forms have traditionally been subdivided into two fractions, “dissolved” and “particulate”, according to an operationally defined limit (0.45 μm) and separated by filtration [63]. Using this partitioning, study waters were filtered with a 0.45-μm filter, with the filtrate analyzed for major and trace elements. In all samples, Fe and Al concentrations were below detection limits (Table 4). The pattern for Mn concentrations was not similar to any other trace element and showed spatial variability, ranging from below detection limits to 173.2 g/L (Table 5). Other water pollution studies [37, 38, 63, 64] have examined additional fractions (dissolved, colloidal, and particulate), with the dissolved fraction generally characterized as the ultrafiltered (<10 kD). Using this partitioning, Kimball et al. [38] similarly found Al and Fe mostly near the limits of detection and Mn measurable in streams near Creede. Differences between concentrations of Fe and Al compared to Mn in study waters may due in part to the geochemical behavior of these solutes, that is, Fe and Al are among the most reactive constituents likely precipitating and forming colloidal materials, with Mn less reactive and more likely remaining in the aqueous phase.

Total concentrations of major elements showed slight variability among sites, ranging from 12 mg/L in Eastfork to 35 mg/L in Spruce (S05CO-079-004B). Overall, Eastfork had lower levels of all major elements compared to other sites. Total major element concentrations were in relatively good agreement across all samples (among and within sites), with mean and standard deviation of  mg/L. The relative trend among major elements showed Ca > Na > Mg > K, with good agreement across all samples (among and within sites) showing means and standard deviations of , , , and  mg/L, respectively.

Total concentrations of trace elements showed significant variability among sites, ranging from 345 to 12108 g/L, with much of this disparity explained by Zn levels (Table 5). Total trace element concentrations within sites were in relatively good agreement during the August sampling. It is expected that these concentrations may vary seasonally or year-round. In a study of the upper Animas River watershed in Colorado, affected by acid drainage, Besser et al. [35] found the annual stream discharge (with high runoff in spring and considerably lower discharge in winter) drives the trace element concentrations in the water column, resulting in aqueous concentrations of dissolved Zn and Cu typically near annual minima during summer and reaching annual maxima in late winter when these concentrations may be greater by a factor of three or more [65].

Water geochemistry showed a wide suite of trace elements present across the sampling area. The relative trend among trace elements showed . Overall water composition, in terms of elements of potential concern, is in good agreement with the reports by the Willow Creek Reclamation Committee [39, 40].

Zinc concentrations ranged from 98.1 and 142.3 g/L in Eastfork and Westfork waters, respectively, to 11088.3 g/L in Cottonwood, S05CO-079-002B, suggesting significant metal loading downstream to Willow Creek. Zinc had the highest concentrations of all metals in study waters. Kimball et al. [38] similarly found East Willow Creek an insignificant source of metal loads and Zn concentrations the highest of all metals in the Willow creek watershed, with its cumulative in stream load one of the highest of all the watersheds in the rocky mountains [6671]. Zinc under oxidizing conditions and in the absence of other anions with which to react tends to remains as a dissolved cation, but in the presence of certain reactants, Zn can form sulfates, oxides/hydroxides, carbonates, and phosphates whose solubilities depend more on pH than on redox potential [62]. August sampling of Willow Creek showed SO4 ranging from below detection limits to 1.7 mmol/L and the absence of Fe and Al. Study water pH ranged from 6.5 to 7.5, with those waters with pH (Gravelbar, S05CO-079-001A-C and Cottonwood, S05CO-079-002A-C) rendering the highest Zn values of 4027.9 and 11088.3 g/kg, respectively. Total concentrations of colloidal Fe and Al in the Animas river near Silverton, Colorado, were considered limiting factors to Zn adsorption during the summer months, with Zn adsorption greater downstream in association with pH increases and during spring when Fe and Al concentrations were higher [37]. Studies that show significant Zn adsorption in the water column over a narrow range of about one pH unit near pH 7 have a large excess of absorbent, usually Fe or Al hydroxides, relative to the total Zn concentration [37, 72, 73].

Copper levels ranged from 3.6 to 10.3 g/L. Cadmium which most often exists in small quantities in Zn ores ranged from 1.0 to 77.2 g/L. Lead was detected in some samples, with the greatest amount (149.6 g/L) in one replicate of the second-day sampling of the Creek water (S05CO-079-003B-2). This higher Pb value may be due in part to the capturing of more sediment in the sample compared to the other samplings and/or reflective of time or length scales of temporal or lateral variability in the stream. In a study of streams near Creede, Kimball et al. [38] found colloidal and dissolved Pb relatively high compared with many streams affected by mine drainage, with colloidal Pb > dissolved Pb and in higher concentrations downstream from Nelson Tunnel versus upstream despite no great variation in pH or major-ion concentrations. With observed pH range of study waters exhibiting mixed solid phases (organics and inorganics) [39, 40], preferential trace metal adsorption sequence in water system may be Pb > Cu > Cd > Zn [74], explaining in part the lack of Pb in some samples of stream waters compared to Zn.

Strontium and Ba levels showed relatively high levels, ranging from 112.7 to 283.1 and 107.6 to 414.0 g/L, respectively. Strontium sources are linked to nonore minerals with its loading profiles indicating sources upstream [38]. Study waters showed relatively low Ag levels, ranging from trace to 0.35 g/L. Arsenic, Be, Cr, P, and Sb were not detected in water samples.

Geochemistry of water samples associated with sampled soils showed higher levels of trace elements compared to stream samplings, for example, Cottonwood water (S05CO-079-002A) with Zn, Cd, and Mn levels of 11088.3, 77.2, and 173.2 g/L, respectively. Cottonwood had some of the lowest pH (6.5 to 6.6) and highest EC (0.20 to 0.2 dS/m) values compared to other sites. Under some circumstances, such as increases in salinity, decreases in redox potential (oxygen deficient), or decreases in pH, metals may be desorbed from sediments and solubilize in the water [75]. These circumstances may include when the water flow is physically obstructed causing water to pool and stagnate. The stream system may have a reasonably good buffering system, with low daily variability in pH and EC. The ability to resist abrupt changes in pH tends to prevent or slow metal conversion to soluble forms, which would increase the concentration of potentially toxic trace elements.

The “reference” or “baseline” sites of Westfork and Eastfork waters had relatively low levels of trace elements compared to other sites with Zn (98.1 and 142.3 g/L, resp.) and Cd (1.0 and 2.6 g/L, resp.). Overall, Eastfork showed lower levels of trace elements compared to Westfork. These data indicate trace element concentrations are significantly different at these higher elevations compared to the sampling sites at the lower elevations.

3.3. Soils: Physical, Mineralogical, and Chemical Properties

Study soils were classified as Aquic Cryofluvents, showing redox depletions with low chroma within the upper 50 cm of surface, and having aquic conditions for some time in normal years [54]. Soil textures were primarily coarse-textured (loam to coarse sand), with the >2-mm fraction ranging from 1 to 73% (Table 6). Soils had mixed mineralogy (kaolinite, smectite, and hydroxyl-interlayer smectite) (Table 6). Smectitic clays typically show the finest particle-size distribution of naturally occurring clay minerals (particles <0.05 μm) and have electrically charged surfaces that permit the adsorption and absorption of water and other substances, whereas the kaolin group of clays includes particles with relatively large size and neutral electrical charge resulting in less potential for persistent suspension and turbidity. Smectitic clays have been linked as a primary cause of persistent turbidity (fine inorganic and organic particles suspended within water column) in the Cascade Range, North Santiam River, Willamette National Forest [76]. For August and July surface and mine water sampling, total suspended solids (TSS) ranged upwards of 15.4, 18.1, and 12.8 mg/L in main channel sites in East Willow Creek, West Willow Creek, and Mainstream Willow Creek, respectively [39].

For sampled soils, the pH (1 : 1 water) for surface materials ranged from 5.4 to 6.8 (strongly acid to neutral), with the exception of the Denuded Mound soil with very strongly acidic pH of 4.6 (Table 6). The Sedge Mound soil which was in close proximity to the Denuded Mound soil had a pH of 5.6.

Acidic conditions and low fertility may explain in part the disparity of plant vegetation in the Sedge Mound (100% base saturation by CEC-7, pH 5.6) compared to the Denuded Mound (34% base saturation and 4.6 pH). Total C and N were also lower in the Denuded Mound (1.33 and 0.185%, resp.,) compared to the Sedge Mound (2.43 and 0.246%, resp.,) (Table 6). Additionally, the Sedge Mound compared to other soils showed the highest clay content (19.2%) and highest amount of available P (26.3 mg/kg) as measured by Bray P-1 (Table 6).

The Salty soil showed pH 5.7 with a relatively high saturated paste EC (8.00 dS/m) (Table 6), with the water soluble anion sulfate (168.1 mmol(−)/L) in greatest amount (data not shown). Typically, these water soluble cations and anions would show greater ion balance, suggesting some cations present in significant amounts were not measured by this method, with potential candidates of Fe, Cu, and Zn sulfates. Melanterite (hydrated Fe sulfate) is one of only a few wate-soluble sulfate minerals, forming white to green encrustations and crystal aggregates. The primary source of Fe for melanterite is Fe sulfates such as pyrite, marcasite, and chalcopyrite. The chemical processes involved in mining that lead to widespread acidification of water bodies and deposition of heavy metals have been studied extensively, with some researchers [77] finding that microbes coat the pyrite surface with exudates that dramatically speed up oxidation and production of sulfuric acid instead of protecting the surfaces from oxidation as had been previously expected.

3.4. Soils: Total Extractable Major and Trace Elements

Total concentrations of major elements showed good agreement among soils. The relative trend among major elements showed , with good agreement among soils (Table 7). This trend was also in relatively good agreement with a dataset of 486 pedons from across the US, inclusive of both anthropogenic and nonanthropogenic soils such as those affected by mining activities, with data showing [6].

Relatively high amounts of K in these coarse-textured soils were related to naturally high levels of potassium feldspars. Optical analysis showed the same trend for predominant minerals among all soils as follows: Chert > Plagioclase Feldspar > Quartz > Potassium Feldspar.

Total Fe and Al ranged from 25 to 38 and 52 to 74 g/kg, respectively. These values exceed the geometric means for Fe and Al ( and  g/kg, resp.) as determined by Burt [6]. Dithionite-citrate and ammonium oxalate extractions [7882] showed “active” Fe, Al, and Mn components, with Fed 1.0% to 1.6%, Alo 0.03% to 0.09%, and Mnd trace to 0.2% (Table 6). The Fed and Alo are generally considered good measures of the total pedogenic Fe and Al in soils, respectively [79]. The Fed/Fet ratios in soils showed free Fe oxides comprise 38% to 46% of the total Fe pool. Data suggest study soils serve as seasonal sources of Fe- and Al-rich colloids to stream waters. The interface between soil surface horizons and streams becomes tightly linked during high-flow events, resulting in increases in stream concentrations and depletions of soil constituents [83].

Total concentrations of trace elements were significantly greater in soils than in water samples. Relatively low concentrations of metals in the Hei River, China despite the contamination to soils and sediments were related to water flow rates high enough to dilute releases from the stream bank soils and sediments [84]. Total trace element (TE) concentrations in study soils showed significant short-range spatial variability, ranging from 2819 in Sedge Mound to 19274 mg/kg in Salty, and seeming to follow no directional pattern along the stream (Table 8). Most of this variability was explained by differences in Zn and Pb levels among soils, varying by factors of nearly 20 and 30, respectively.

Trace elements from anthropogenic sources have been found to be more mobile than those of native origin [8587]. A common method to quantify the release and transport of metals in surface–contaminated soils is the distribution of these contaminants within the soil profile [88]. Total concentrations were higher in surface layers of study soils compared to underlying materials, with Spruce showing 10991 versus 9925 mg/kg, respectively, and Railroad 11117 versus 7893 mg/kg, respectively. These data are in agreement with Burt [6] finding the median value of most trace elements (Cd, Cu, Ni, Pb, Zn, and Hg) significantly higher in surface horizons of anthropogenic soils () compared with those in surface horizons of nonanthropogenic soils.

Similar to the water data, there was a wide suite of trace elements present in the soils across the sampling area. The relative trend among trace elements showed , with good agreement among soils. Geochemistry of soils does not parallel the geochemistry of the stream waters. Relatively high Pb amounts in the soils compared to the water samples may be reflective of its relatively low solubility in water and low mobility in soils, accumulating primarily on the surface where its increasing presence can begin to affect soil microflora [89].

The Denuded Mound and Sedge Mound soils were demonstrative of the vegetative disparities along the stream bank. The Denuded Mound and Sedge Mound showed significant differences in total trace element concentrations (15317 versus 2819). Soil pH for the Denuded Mound and Sedge Mound was 4.6 versus 5.6, respectively. These results would tentatively suggest that a correlation may exist between metal levels and soil pH and vegetation or lack thereof along the stream. In a study of smelter–contaminated soils (Cu, Zn, Pb, and Cd) in Anaconda and Deer Lodge Valley, Montana, Burt et al. [49] found that low pH, alteration in species composition in plant communities, and reduction in plant growth and vigor often occur in conjunction with soil contamination, but that these easily measurable and observable properties alone could not be used as field indicators of total trace element concentrations. Observable differences in these Montana plant communities were determined to serve more appropriately as criteria to determine impact classes (e.g., severe, moderate, and slight) for soils affected by contamination but not correlated to total trace element content [49]. Soil pH is useful information in terms of bioavailability of trace elements [9092].

3.5. Soils: Chemical Fractionations

Percent recovery ranged from 104 to 140%, indicating relatively good agreement between SUMSEQ and SUMTE. There was also good reproducibility between the two replications of the sequential fractionation (data not shown). Precision and accuracy of the sequential extraction procedure are generally considered inherently good, with the limiting factor the inherent heterogeneity of the specimen [10]. Study soils were located on the floodplain with colluvial and alluvial deposits which create a more heterogeneous medium for soil development (e.g., horizonation, morphology, and texture) [93]. For more details, see Table 9.

In general, study soils showed the OX (15 to 68%), SS/CAR (23 to 47%), and RES (7 to 25%) as the predominant fractions over the EX (3 to 22%), OM/S (5 to 14%), and WS (1 to 16%) fractions. Element–specific bioavailable fractions (WS and EX) showed Pb (0.29% to 4.60% and 0.94% to 24.27%, resp., see Table 10), Cd (0.99% to 23.37% and 5.28% to 40.18%, resp., see Table 11), Zn (0.77% to 24.39% and 4.14% to 32.67%, respectively, see Table 12), and As (0.09% to 2.15% and 0.39% to 2.73%, resp., see Table 14). In contrast, the Montana smelter soils [11] had lower total metal concentrations (1836 to 3605 mg/kg), smaller WS (0.02% to 1.0%) and EX (0.2% to 9.3%) fractions, larger RES (14.3% to 70.1%) fraction, and smaller potentially bioavailable fractions (Pb, Cd, and Zn), exception being the greater EX Cu fraction in the Montana soils (4.67% to 30%) compared to Willow Creek (1.01% to 5.32%, Table 13). Similarly, the Montana smelter soils showed the OX fraction (29.7% to 39.8%) as the predominant fraction, with this fraction exceeding the RES fraction in surface layers in both studies. These data suggest that some anthropogenic inputs were added, as well as later transformed to, OX forms. Anthropogenic metals, initially present as poorly crystalline phases, may oxidize in surface soil horizons. Oxides have a scavenging action for trace elements disproportionate to their own concentrations, and as such constitute important sources of potentially bioavailable trace elements [10]. Overall, results suggest the Willow Creek soils are more contaminated than the Montana smelter soils, showing not only greater total metal levels but also an overall greater pool of potentially bioavailable trace elements, specifically those of concern (e.g., Cd, Pb, and Zn).

Soils in both studies have been affected by over 100 years of anthropogenic inputs, with the Willow Creek soils from mining and the Montana study soils from smelting of ores, primarily copper. While both sites have been inactive for a number of years, the geochemical processes are more active and impacts on–going in Willow Creek compared to the Montana study soils, related to site–specific differences, for example, topography, geology, climate, soil properties, and contaminant type. The Willow Creek soils are located on a floodplain whereas the Montana soils are on uplands with steep slopes (e.g., 30%), suggesting more active processes of metal transport to and potential accumulation in the Willow Creek soils. Smelter–derived particles have been shown to rapidly alter soon after deposition through weathering processes, leading to a release of metals in soils [94]. The aquic conditions of the soils along Willow Creek provide for a more dynamic system of metal release and further transport into the soil–water environment compared to the aridic soils in Montana. The higher pH Montana soils (derived from calcareous alluvium) provide for a greater buffering against acidification, resulting in a less soluble and therefore less mobile trace element fraction compared to the strongly acid to neutral pH Willow Creek soils (derived from ashflow tuffs). Soil pH and redox conditions are considered important factors affecting the chemistry of metals in soils and their uptake by organisms, with alternating aerobic and anaerobic conditions leading to changes in both factors thereby affecting most of those processes regulating metal speciation in soils [9599].

All study soils showed high levels of contamination, with the Salty and Denuded Mound soils being the most contaminated showing the greatest SUMTE (19274 and 15317 mg/kg, resp.) and some of the greatest WS (16 and 8%, resp.) and EX (16 and 22%, resp.) fractions. The Salty soil with pH 5.7 and EC 8.00 dS/m had the highest Cd levels of total, WS, and EX fractions, with 84, 19.4, and 21.1 mg/kg, respectively. In general, soluble Cd levels are greater in moderately acid, oxidizing conditions, increasing 10 fold under these conditions compared to reducing conditions at all pH levels [75]. Cadmium has also been shown to be affected by salinity, with increased salinity supporting increased Cd concentration in solution [75]. The Salty soil also had the highest levels of Zn with 12134, 2959.0, and 2850.1 mg/kg in the total, WS, and EX fractions, respectively.

Metal mobility was evaluated using the mobility factor (MF) [59, 100], representing the proportion of metals in the water soluble and exchangeable fractions and calculated as follows: where x = element.

The MF was assessed for the highly contaminated salty and denuded soils as follows: . Values ranged as follows: 48.7 to 55.3% Cd; 44.7 to 47.9% Zn; 1.4 to 27.6% Pb; 1.9 to 2.8% Cu; and 0.5 to 0.6% As. These MFs are in relatively good agreement with other studies [101, 102], with Cd and Zn showing the highest mobility. Low CEC, sandy nature, and high permeability of study soils indicate inherently low sorption capacity for Cd and Zn, and as such high mobility. Generally, in all study soils, Cd appears mainly in the SS/CAR form, whereas Zn in the SS/CAR, OX, and OM/S forms. The relatively high MFPb in the Denuded soil may be due in part to having the lowest OX fraction (15%) compared to other study soils, and similarly a low OM fraction (9%). Both of these fractions (Fe–Mn oxide and organic) have been shown to adsorb or complex Pb in less available forms in soils [103105]. Generally, in all study soil, Pb appears mainly in the SS/CAR and OX forms. The low MFAs for the salty and denuded soils may be due in part to its greater association with the OX and RES fractions over the more soluble fractions. These results show good agreement with other studies investigating As forms and distribution in soils [106109].

4. Conclusions

Historic mining activities and impacts from these activities have resulted in significant fluvial contamination of sampled soils and water in the Willow Creek floodplain. Thorough understanding of the factors affecting this contamination requires a good knowledge of the local geology, geomorphology, hydrology, climate, soil properties, and contaminant amount, type, and form. There is a wide suite of trace elements present in both soils and waters. Trace element forms in water are strongly affected by site-specific water quality factors such as pH, EC, and other dissolved constituents. Total concentrations of trace elements in water showed significant variability among sites, with much of this disparity explained by Zn levels. Soil properties such as pH and redox conditions affect chemical forms and potential reactivity. Total trace element concentrations in study soils showed significant short-range spatial variability, and seemed to follow no directional pattern along the stream. In general, study soils showed the OX, SS/CAR, and RES as the predominant fractions over the EX, OM/S, and WS. The relative distribution of trace elements within the more soluble fraction (e.g., WS and EX) of soils is similar to their respective total concentrations, suggesting a relationship between total trace elements and the concentration of soluble, mobile fraction. Largest element–specific bioavailable fractions (WS and EX) were found in Cd, Zn, and Pb.

The expected outcome of the Willow Creek study and similar studies of anthropogenic soils conducted by USDA Soil Survey (e.g., [6, 7, 11, 4951, 110]) is to better characterize important soils in specified–project areas, and more generally establish important relationships with other soil properties as well as with other soils both uncontaminated and contaminated. Our knowledge of the various soil processes affecting the elemental amounts and distribution in US soils and their relationships with other soil properties will enhance the understanding of the fate and transport of anthropogenic elements, thereby expanding the utility and application of soil survey knowledge in areas of environmental concern (e.g., urban, mine spoil reclamation, munitions disposal, smelter emissions, and agricultural waste applications).