Tungsten (W) became a new emerging contaminant of concern for the United States (U.S.) Environmental Protection Agency (EPA) in 2002 after multiple studies correlated elevated tungsten concentrations with adverse health effects and environmental risks(EPA 2014). Three childhood leukemia clusters appeared between 1997 and 2003 in Fallon, NV, Sierra Vista, AZ, and Grove, CA, prompting extensive testing and exposure assessments of environmental contaminants in drinking and surface water, surface soils, residential indoor dust, and air (Rubin et al. 2007; P. R. Sheppard and Witten 2003; Paul R. Sheppard et al. 2012; Dermatas et al. 2006; Larson and Dermatas 2009; Strigul, Koutsospyros, and Christodoulatos 2010; Strigul et al. 2005; Koutsospyros et al. 2006; Tuna et al. 2012a). Biological testing done by the Center for Disease Control (CDC) reported high levels of tungsten in urine (as high as 15 times the national average) from residents near Fallon, NV (CDC 2003), as well as high tungsten concentrations in trees from all three communities (P. R. Sheppard and Witten 2003). All three towns are located near anthropogenic sources of tungsten including active or inactive tungsten mines, processing operations, and/or military bases. These events prompted the National Center for Environmental Health (NCEH) and CDC to designate tungsten for further study under the National Toxicity Program (ATSDR 2005). While the events in Fallon, NV and Sierra Vista, AZ brought awareness to the threat of tungsten toxicity, to date there are no drinking water standards or regulations for the discharge of tungsten into the air, surface water, groundwater, or soils in the United States (EPA 2021).
The market for tungsten has developed quickly since its official discovery in 1783 and tungsten is now recognized as one of 50 minerals deemed “critical to U.S national security and the economy” by the Energy Act of 2020 (Nassar and Fortier 2021). Tungsten ore is used to produce tungsten alloys and tungsten carbide which have a wide range of applications including welding, high-temperature cutting/drilling technologies, aerospace industries, filaments in light bulbs, and in manufacturing of electrodes and components used in touchscreen products (EPA 2014; U.S. Geological Survey 2021; Kabata-Pendias and Mukherjee 2007). While tungsten has been used in munitions since the Second World War, production ramped up in the mid-1990s following a military initiative to go “green” by substituting toxic lead and radioactive uranium with an “environmentally friendly” alternative thought to be tungsten (EPA 2014), (Jay L. Clausen et al. 2011), (Jay L. Clausen et al. 2007), (Dermatas et al. 2006), (Bostick et al. 2018), (Koutsospyros et al. 2006). China is currently the world’s largest producer of tungsten controlling more than 80% of the world’s supply but has effectively discontinued their exports in recent years, prompting tungsten to be moved to the top of the “endangered list” due to the limited supply (BGS 2016). There have been no active U.S commercial tungsten mines in operation since 2015 (USGS n.d.), however in 1955, the U.S. boasted more than 740 operating tungsten mines (Shedd 2002) and according to USGS scientist Jeff Mauk, “the United States may have future production potential from U.S. projects that are currently in the advanced exploration stage,”(Burton and Demas 2021) . The economic importance of tungsten and the growing concern about limited supply gives an urgency to the research on the environmental fate and transport of tungsten. High concentrations of W have the potential to flow and seep into natural waters and soils from these legacy tungsten mines as well as future tungsten mines. Identifying the controls on tungsten mobility at these abandoned mines is a vital first step towards mitigation.
The purpose of this study is to identify if adsorption is the main control on W contamination mobility in surficial soils and surface water from exposed tailings and waste piles at the Boriana and Black Pearl abandoned tungsten mines, in western, Arizona. To determine whether adsorption is the main control on W mobility, we must identify controls on tungsten adsorption in this field setting. The ultimate goal of this research is to provide a framework we can build upon to predict and prevent the mobilization of tungsten contamination from anthropogenic sources in the environment.
The previous and long held notion that tungsten is immobile and environmentally “green” (Kerwien 1996; Koutsospyros et al. 2006) is not supported by current research. Studies conducted at small arms firing ranges that utilized “green” ammunition found tungsten contamination in groundwater, subsurface and surficial soil, and plants Jay L. Clausen et al. (2007); Jay L. Clausen et al. (2011); A. J. Bednar et al. (2006); A. J. Bednar et al. (2009); Bostick et al. (2018); Seiler, Stollenwerk, and Garbarino (2005); Tuna et al. (2012a); Datta et al. (2017); Johannesson et al. (2013). Elevated tungsten levels attributed to nearby active mining operations have been found in soils, the crops grown in the contaminated soils in rice fields in China (Silvetti et al. 2014; Lin et al. 2014; Bob Wilson and Pyatt 2006; B. Wilson and Pyatt 2009), as well as in soils and plants near inactive mining operations in southwestern France and northwest England.
Elevated tungsten levels observed in anthropogenically contaminated surficial soils occur as two major aqueous species: metallic tungsten W(0) and oxide species W(VI) (Bostick et al. 2018; Jay L. Clausen and Korte 2009; Seiler, Stollenwerk, and Garbarino 2005; Strigul et al. 2005; Tuna et al. 2012b; Hobson et al. 2020). Bostick et al. (2018) reports that 98% of tungsten in contaminated soils from a munitions firing range is present as W(VI) after only a few years of exposure to the elements. This rapid environmental oxidation of metallic W(0) to highly soluble W(VI) species suggests W(VI) is the redox state of interest when considering transport and fate (Hobson et al. 2020; Datta et al. 2017; Gustafsson 2003). While aqueous speciation of tungstate W(VI) is complicated and geochemical properties of polymer species seem to vary from those of monomeric tungstate (Davantes et al. 2017; A. J. Bednar et al. 2006; Jay L. Clausen et al. 2011; Petruzzell and Pedron 2021; Hur and Reeder 2016), recent studies have concluded the fundamental controls on tungsten speciation include pH, concentration of tungsten, and presence of other solutes such as silica and phosphorus (Wasylenki et al. 2020; Bostick et al. 2018; Hobson et al. 2020; Datta et al. 2017; Johannesson et al. 2013; A. J. Bednar et al. 2006). The monomeric tungstate species (WO42-) observed at low concentrations and high pH (>8 pH) appear stable and dominate in solution while polymeric tungstate species observed at high concentrations and low pH (<5 pH) appear to dominate solution (Bostick et al. 2018; Johannesson et al. 2013).
Previous reviewed literature from the past two decades has progressed our knowledge on the controls of tungsten mobility significantly however the question of applicability of these results to natural field settings remains a challenge . Tungsten concentrations reported for non-polluted soils in the U.S range from 0.5 to 5.0 mg/kg, from 0.01 to 0.05 µg/L in non-polluted stream waters, and from 100 to 200 mg/kg in soils surrounding tungsten ore-processing plants (Kabata-Pendias and Mukherjee 2007). Concentrations in the majority of these experiments are 1 to 3 orders or magnitude higher than what would be typically observed in a field setting for non-polluted soils and 8 orders of magnitude higher for surficial water (1420 mg/kg (Jay L. Clausen and Korte 2009), 515 mg/L & 700 mg/kg (Strigul et al. 2005), 7080 mg/kg (A. J. Bednar et al. 2009), 50 mg/kg (Johannesson et al. 2013)). Simple experiments to determine the aqueous speciation and sorption affinities of tungsten in a laboratory at field relative tungsten concentrations so far do not address the effects of continual water flow or other variability introduced in natural environmental settings (Wasylenki et al. 2020). Field testing completed at munitions testing facilities are not only representative of extreme examples of concentrated contamination in soils (up to 5100 mg/kg [W]) but also only represent one anthropogenic source of tungsten contamination (Jay L. Clausen et al. 2011; A. J. Bednar et al. 2009; Bostick et al. 2018). Methods and procedures for the collection, preparation, and analysis of tungsten in the environment are still under development by U.S federal regulatory agencies (EPA United States Environmental Protection Agency 2017). As such, the effectiveness and reproducibility of results from various laboratory testing procedures is another challenge. Tungsten can be measured using ICP-MS, and multiple studies utilize modified EPA acid digestion methods to extract bulk concentrations of tungsten from soil and water (B. Wilson and Pyatt 2009; Hobson et al. 2020; Jay L. Clausen et al. 2007; Zimmerman and Weindorf 2010; A. J. Bednar et al. 2009; Cutler and Stillings 2011; Dold 2003; Assam and Stillings 2007; Joseph Taggart and Editor, n.d.; Cao et al. 2020; Li et al. 2019; Wenzel et al. 2001; Gómez-Álvarez et al. 2015; Violante et al. 2010; Tessier, Campbell, and Bisson 1979). In addition to the variations of acid digestions for bulk extraction, many studies include variations of sequential extraction procedures (SEP) which aim to test the controls of W mobility using the relationship between mobility and W fractionation (Bostick et al. 2018; Hobson et al. 2020; Zimmerman and Weindorf 2010; Cutler and Stillings 2011; Joseph Taggart and Editor, n.d.; Li et al. 2019; Wenzel et al. 2001; Tessier, Campbell, and Bisson 1979; Gaudino et al. 2007; Kalyvas, Gasparatos, and Massas 2018; Petruzzelli and Pedron 2017; Gómez Ariza et al. 2000; A. J. Bednar et al. 2010). Most SEPs follow a similar fractionation order generally cited from Tessier et al. (1979) beginning with exchangeable, then carbonate bound, followed by manganese and iron oxide bound, organic matter, and finally residual. However, many studies publish drastically different and often conflicting results even when utilizing the same SEP (Li et al. 2019; Cao et al. 2020; Gómez Ariza et al. 2000; Alan and Kara 2019; Hall, Jefferson, and Michel 1988; Claff et al. 2010; Violante et al. 2010; Miller, Martens, and Zelazny 1986; Mesko et al. 2013).
The commonalities between the Boriana and Black Pearl mines make them comparable study sites and will potentially yield corroborating results. Both the Boriana and Black Pearl mines operated c. 1914-1956 with a majority of mining operations occurring below ground and milling and processing operations above ground. These mines produced ore with an average recovery percentage of wolframite (WO3) at 70.0% (Dale 1961; Mines 1991; Stein, Hannah, and America. 1990). At present, mine waste and tailings piles are barren and exposed with no coverings or repository closure methods in place. Both mine localities are characterized as semi-arid regions averaging less than 17 inches of precipitation annually and include ephemeral streams as drainage channels that discharge into the intermittent streams of Mackenzie Creek (Boriana) and Boulder Creek (Black Pearl). I hypothesize these drainage channels act as preferential pathways for tungsten and other heavy metal contaminants to mobilize in the dissolved phase via surface runoff and precipitation events to downgradient streams and surrounding watersheds. Surficial soil and water samples were collected in October 2021 and March 2022 (Maps 1 & 2) from the Black Pearl and Boriana Mines and pseudo bulk soil analysis was completed on a subset of soil samples from the Black Pearl Mine using an inductively coupled mass spectrometer (ICP-MS) in April and June 2022. The dataset used for visualization in this report comprises of a subset of these samples.
Map 1. Sample location map of the Black Pearl Mine.
Map 2. Sample location map of the Boriana Mine.
Fieldwork and Sample Collection Fieldwork consists of the collection of surficial soil samples at the Boriana and Black Pearl Mines from the source points (tailings and waste piles), upgradient from source points, downgradient from the source points within the drainage channels at 20-25m intervals to the intermittent streams of Boulder and Mackenzie Creek, in the intermittent streams of Boulder and Mackenzie Creek above, between, and below the entry points of drainage channels, and downgradient from source points outside the drainage channels. Water samples and water quality measurements using a YSI multi-parameter probe will be collected from Boulder and Mackenzie Creek above, between, and below the entry points of drainage channels. Soil samples are collected in air-tight plastic bags while water samples are collected in 15mL plastic centrifuge tubes and stored in a cool dark place. To prevent cross contamination in the field, sample gloves are switched out between each sample. Following field collection, grain size analysis is done on the soil samples using plastic #40 and #200 sieves to determine if there is a relationship between grain size and W adsorption during transport.
Laboratory Analysis Laboratory work includes pseudo bulk concentration analysis via acid digestion according to modified US EPA Method 3050b as described in Bednar et al. (A. J. Bednar et al. 2010, 2009) and the SEP described in Hobson et al. (Hobson et al. 2020). To complete pseudo bulk analysis, soil samples will be digested in concentrated nitric and phosphoric acid with heat and then in hydrogen peroxide to leach secondary precipitates deposited on the sediment grains during runoff, filtered through 0.2µm polypropylene syringe filters, dried, then dissolved in 2% nitric and trace hydrofluoric acid (HF), diluted, and analyzed using an Inductively Coupled Plasma Mass Spectrometer (ICP-MS) to determine W, Fe, and Mn concentrations. In conjunction with W, Fe, Mn, Ni, Cu, Zn, Mo, Pb, and U may also be analyzed in select samples to determine any correlations between other heavy metals and W adsorption from these mine tailings. Water samples are filtered (same filters as soil samples), diluted with 2% nitric acid and trace HF, and analyzed using on the ICP-MS for W, Fe, and Mn concentrations. Bulk soil concentrations are extracted using the following calculation. \[\frac{MetalAmount(ng)}{BulkSoil(g)} = ICP-MS(ppb) * \frac{QuadSample}{D1SampleAliquot} * \frac{Digested Sample}{Bulk Sample} \]
Figure 1. Concentration plots of W vs all other heavy metals at Black Pearl Mine.
Figure 2. Concentration plots of Fe vs all other heavy metals at Black Pearl Mine.
Figure 3. Concentration plots of Mn vs all other heavy metals at Black Pearl Mine.
Map 3. Heavy metal concentrations at the Black Pearl Mine.
Map 4. Ratio of heavy metals to Fe concentrations at the Black Pearl Mine.
Map 5. Ratio of heavy metals to Mn concentrations at the Black Pearl Mine.
Finally, we can statistically compute these relationships using principal component analysis (PCA). In figure 4, we can see there is a clear and positively correlated relationship between Fe, W, and Mo and they have the strongest effect on PC1 while Mn is inversely correlated and strongly effects PC1. I also feel confident interpreting from the PCA that Pb, Cu, and U are correlated with a stronger effect on PC2. There also seems to be some minor correlation with Zn and Mn which we also observe in Figure 3.
Figure 4. PCA of soil concentrations at the Black Pearl Mine.
Thus far, I think the results I have trend with what we expect to see geochemically at the Black Pearl Mine. While there are huge amounts of error in this dataset, the largest source being my own initial laboratory work, I have made tremendous progress developing my laboratory procedures the past few months that should significantly cut down on that error. After I have more completed laboratory testing on more samples however, I believe PCA as well as the spatial visualization will be extremely useful tools for my research.