- Open Access
Micro-spectroscopic investigation of selenium-bearing minerals from the Western US Phosphate Resource Area
© American Institute of Physics 2005
- Received: 14 September 2004
- Accepted: 08 December 2004
- Published: 28 January 2005
Mining activities in the US Western Phosphate Resource Area (WPRA) have released Se into the environment. Selenium has several different oxidation states and species, each having varying degrees of solubility, reactivity, and bioavailability. In this study we are investigating the speciation of Se in mine-waste rocks. Selenium speciation was determined using bulk and micro-x-ray absorption spectroscopy (XAS), as well as micro-x-ray fluorescence mapping. Rocks used for bulk-XAS were ground into fine powders. Shale used for micro-XAS was broken along depositional planes to expose unweathered surfaces. The near edge region of the XAS spectra (XANES) for the bulk rock samples revealed multiple oxidation states, with peaks indicative of Se(-II), Se(IV), and Se(+VI) species. Micro-XANES analysis of the shale indicated that three unique Se-bearing species were present. Using the XANES data together with ab initio fitting of the extended x-ray absorption fine structure region of the micro-XAS data (micro-EXAFS) the three Se-bearing species were identified as dzharkenite, a di-selenide carbon compound, and Se-substituted pyrite. Results from this research will allow for a better understanding of the biogeochemical cycling of Se in the WPRA.
- Shale Sample
- Atomic Shell
Approximately 40% of the United States phosphate reserves are located in southeast Idaho, northern Utah, and western Wyoming in an area known as the Western Phosphate Resource Area (WPRA). The phosphate ore in the WPRA is located in the Meade Peak Member of the Phos-phoria Formation dating to the Permian period. Approximately 6 × 106 Mg of ore is extracted annually from this area resulting in 36 × 106 Mg of waste material. Rex Chert and Cherty-Shale Members overlay the Phosphoria Formation and Pennsylvanian period Wells Formation (limestone) lies below the Phosphoria Formation. The Phosphoria Formation contains two ore zones. Typical thickness of the upper and lower ore zone is 6–14 m, with 23–55 m of phosphatic shales (middle waste shale) separating the two ore zones.[1, 2]
Selenium is an essential micronutrient with a narrow range between sufficiency and toxicity. In the environment, birds, hoofed animals, and aquatic species are particularly sensitive to Se, and may show chronic or acute signs of toxicity. Bioavailability is dependent on solubility, which also controls Se mobility through the ecosystem. Thus, in order to manage risks and remediation of the soils derived from the Se containing rocks in the WPRA better knowledge of Se speciation in the parent materials is necessary.
Selenium exists in multiple oxidation states, with each state having different fates within the environment. Selenides [Se(-II)] and elemental Se(0), which are reduced forms of Se, are formed under anaerobic conditions, and are not considered to be mobile or bioavailable. Selenite (Se O32-) and selenate (Se O42-) are more soluble than the reduced Se species. Selenite forms adsorption complexes on ferric oxide complexes, which are not readily bioavailable. Selenate is the most soluble Se species and is readily taken up by plants and leached through the soil profile. Distribution of selenite and selenate between the solid and solution phases within soil has been shown to be a function of pH.
Past work conducted in the WPRA identified a number of Se-bearing species in shale material, and found substantial heterogeneity. Desborough analyzed samples from a "carbon seam" at Maybe Canyon mine in Idaho by electron microprobe, x-ray diffraction, and bulk chemical analysis and concluded that the shale samples contained no pyrite, suggesting that the selenium in the rock was associated with organic material. In contrast, samples from the Dry Valley Mine (within the WPRA) contained both framboidal and angular pyrite, with Se concentrations ranging from <100 (detection limit) to 2000 mg kg-1. By using a combination of x-ray absorption near edge spectroscopy (XANES) and sequential-extraction techniques, Perkins and Foster determined that in unweathered, or minimally weathered shale samples from the WPRA, the largest fraction of Se was associated with iron sulfides. Most of the non-sulfide-associated Se was associated with organic matter and oxyhydroxides, and a small amount of Se was present in elemental form. In weathered shale samples the dominant form of Se observed was Se(IV) associated with oxyhydroxides.
Nondestructive x-ray spectroscopic methods such as x-ray absorption fine structure (XAFS) spectroscopy are powerful tools in studying Se speciation. However, the heterogeneity that exists in natural samples often makes it difficult to interpret speciation using bulk-XAFS spectroscopy. In recent years the development of microscopically resolved XAFS spectroscopy has created new opportunities for interpreting speciation and spatial heterogeneity within natural materials. Micro-XAFS techniques, for example, can be used to obtain chemical and structural information on the micrometer-scale typical of mineral grains. Manceau et al. suggest that even though the nanoscale detail cannot be seen with micrometer size probes, such as micro-XAFS spectroscopy, it can be inferred with confidence. Strawn et al. used micro-XAFS spectroscopy to identify two iron aggregates that were located within a few microns of each other in undisturbed soil samples. They were also able to determine that As and Se were more concentrated in iron oxide aggregates rather than the clay/shale matrix. Tokunaga et al. also used micro-XANES to obtain spatially resolved distribution maps of Se oxidation states in natural samples. In the research presented in this paper we use micro-XAFS spectroscopy to investigate the species of Se present in the middle waste shale from the Meade Peak Member of the Phosphoria Formation.
A. Sample collection
We collected shale samples from the Smoky Canyon mine in The Caribou National Forest in southeastern Idaho. For micro-XAFS analyses a shale rock was used that had no visible fractures, implying that the interior was unweathered. The shale was broken along depositional planes to expose fresh surfaces. Two separate surfaces were prepared for analysis (referred to as samples 1 and 2). The bulk Se concentration in the shale was 140 mg kg-1. A second shale sample with a concentration of 34 mg kg-1 was ground to fine powder for bulk XAFS spectroscopy (referred to as ground shale-Smoky Canyon).
Description of bulk rock samples used in this study collected by USGS (see Sec. II A)
Unit within mead peak member
As (mg kg-1)
Cu (mg kg-1)
Se (mg kg-1)
Middle waste shale
Upper waste shale
Middle waste shale
Upper waste shale
Selenium standards, purchased from Alfa Aesar or Sigma Aldrich, included selenium sulfide, copper selenide, zinc selenide, selenomethionine, sodium selenate, sodium selenite, orthogonal (gray) elemental selenium, selenocystine, Se-methyl-selenocysteine, and selenoethinoine. The minerals berzelianite, chalcomenite, krutaite, mandarinoite, and penroseite were purchased from the Excalibur Mineral Company (New York). Ferroselite was synthesized following the procedures described by Warren, and checked by powder x-ray diffraction to confirm its crystal structure. Monoclinic (red) elemental selenium was made by microbial reduction of selenite. Standards used for XAFS spectroscopy were diluted in boron nitride.
C. X-ray spectroscopy
Micro-XAFS spectroscopy and fluorescence microprobe mapping were done on beamline 10.3.2 at The Advanced Light Source (ALS), Lawrence Berkeley National Laboratory (Berkeley, CA). Beamline 10.3.2 is designed to map elemental distributions and collect XAFS data on small areas. A seven-element germanium detector was used to measure sample fluorescence. A monochromatic x-ray beam tunable to energies from 2.0 to 17 keV using a Si(111) crystal set is focused on the sample using a Kirkpatrick-Baez mirror pair. The shale samples were attached to the sample holder using double-sided tape and the sample holder is attached to an XY stage that is oriented 45° to the beam.
For elemental mapping micro-synchrotron x-ray fluorescence (micro-SXRF) was used. For sample 1, the elemental distribution map was collected using a beam focused to 5 × 5 μm (horizontal and vertical dimensions, respectively) rastered across the sample at 5 μm intervals. The elemental distribution map for sample 2 was collected in the same manner as sample 1, except the beam size was 10 × 10 μm and step size was 7 μm. Dwell time on each spot was 100 ms. The beam energy was tuned to 12 757 eV. The windowed (regions of interest) x-ray fluorescence counts for Ca, Cr, Cu, Fe, Se, Ti, and Zn were collected as the sample rastered through the beam. To determine the relationship between S, Fe, and Ca in sample 1 the beam was tuned to ~3988 eV with a beam spot size of 6 × 6 μm (relying on harmonics to excite the Fe and Ca K edges).
All spectra were collected at room temperature. Select Se standards were measured on beamline 10.3.2. Additional selenium standards and the powdered shale samples were measured on beamline X-11b at The National Synchrotron Light Source (NSLS) at Brookhaven National Laboratory, Upton, NY. Energy is focused at this beamline using a Si(111) monochromator. Transmission spectra from Se standards diluted in boron nitride were collected using N2-filled ionization detectors. Fluorescence spectra from the powdered reference materials were collected on beamline X-11b using a krypton purged Stern–Heald type detector. A Se(0) standard (orthorhombic Se, Sigma Aldrich) was run either between samples (beamline 10.3.2) or in situ (X-11b) for energy calibration. The XAS data were averaged (three to seven scans), normalized, and background subtracted using standard procedures in the program ATHENA. Step size through the XANES regions was 0.70 eV. Measurement of the standard deviation of the edge energies from several points of interests with the same XANES features between two separate beamtime runs indicated that the edge-energy precision was less than the step size. Thus we estimate that the calibrated data are accurate to within 0.70 eV. Extended x-ray absorption fine structure (EXAFS) data were transformed into k space (χ spectra) using the sample absorption edge. Finally the data were transformed into frequency distributions by Fourier transforming (FT) the data using square windows.
The EXAFS data were fit using theoretical backscattering amplitudes and phase shifts calculated in the program FEFF 8.0 using the unit cell and space groups for dzharkenite (FeSe2, pyrite structure) or Se-substituted pyrite. Fitting was done in WINXAS. The amplitude reduction factor was fixed to 1, which was derived from optimizing the fit of monoclinic elemental Se. Coordination numbers were fixed at values appropriate for the atomic shells from the mineral that most closely simulated the unknown EXAFS spectra. This is consistent with the approach used to fit overlapping first shells in arsenic sulfides. We also used the results from the XANES analyses to help guide the atomic-shell selections. To determine which models best simulated the unknown spectra we fit the first peak in the FT with two shells of Fe and either Se or S. Then, based on first shell fit results, we deduced which mineral the atomic distances most closely matched, added the additional atomic shells, fixed coordination numbers, and allowed the Debye–Waller factors and bond distances to vary. Fixing the coordination numbers was necessary so that the degrees of freedom did not exceed the number allowed during the fitting procedure. Additionally, to reduce the number of atomic shells necessary to fit the data we assumed that atomic shells with the same atoms that were closer than ~0.2 Å could be fit as a single shell. Bond distances derived from fitting these shells represent weighted averages, and Debye–Waller factors reflect the deviation in the inter-atomic separations. The final best fit was derived based on the reduced chi square and qualitative agreement between the dominant features in the fit and the raw x and FT spectra. As a result of the fitting we were able to compare atomic coordination environment and derived bond distances with those in theoretical minerals, and evaluate the relative local atomic disorder (the Debye–Waller factor) for each shell. Errors for bond distances were less than 3% as estimated by comparing the fit-derived parameters for ferroselite EXAFS data to the theoretical values.
Principal component analysis was done on the micro-XANES spectra to determine the number of end members that comprise the 9 micro-XANES. The minimum of the indicator (IND), which is at a minimum when the correct number of significant components is reached, was used to assess the minimum number of components necessary to reproduce the spectral set. Target transforms using the number of components necessary to define the sample space was done to determine if the components could be matched with known reference compounds. To determine whether the target transforms were acceptable, moderately acceptable, or unacceptable, the SPOIL and the chi squared values were evaluated, as well as visual confirmation that the target was well defined by the components.
A. XANES standards
Selenium K-edge energies for reference standards, based on maximum of first derivative of XANES spectra
K-edge energy (eV)
Organic Se compounds
FeSe2 (Synthetic Ferroselite)
Se (0) compounds
Elemental Se (monoclinic)
Elemental Se (orthorhomblic)
Se (IV) compounds
Se (VI) compounds
Selenium standards had absorption K-edge energies between 12 656.0 and 12 664.0 eV (Fig. 2 and Table II). In general, the reduced Se compounds have edge energies between 12 656 and 12 659.5 eV, which are significantly less than the K edges of oxidized Se species [Se(IV) = 12 662.0 eV and Se(VI) = 12 664.0 eV]. For both allotropes of zero-valent Se the edge energy is 12 658.0 eV. For the selenides the edge energies range between 12 656.0 and 12 658.0 eV. For the organic Se compounds the K-edge energies ranged from 12 657.5 to 12 659.5 eV. Although the edge energies for the standard compound classes overlap, the overall XANES for the compounds are distinct, indicating that Se-XANES is highly sensitive to oxidation state and local atomic coordination, and thus can be used to elucidate speciation.
B. Elemental distributions
Electron microprobe analyses of siltstone/shale thin sections from the Smoky Canyon Mine in the WPRA were performed by Bond. This work identified pyrite and sphalerite grains in the rock, with pyrite comprising about 0.2% (by area) of the sample studied. Selenium was primarily associated with pyrite at concentrations approaching 0.7% (by mass) in the grains, fully accounting for the bulk Se concentration in the sample. Desborough et al. and Munkers used scanning electron microscope-energy dispersive spectroscopy (SEM-EDS) analysis and documented the presence of pyrite grains with both framboidal and cubic morphology in the middle waste shale, with total Se levels in the grains ranging from 0.1% to 3% (by mass).
Perkins and Foster used SEM-EDS and electron microprobe to analyze shale samples from the Meade peak formation. They observed that Se is associated with sulfide minerals, particularly iron sulfides such as pyrite. Perkins and Foster further reported that selenium concentrations detected in pyrite were 10–100 times higher than concentrations observed in bulk-rock, with average rock concentrations exceeding 2000 mg kg-1. Positive correlations between Se and Cu, Ni, Cd, and Zn were also detected.
Although a positive correlation for enrichment of Se in iron sulfide minerals and shale has been observed, speciation and mineralogy of Se in these materials is still unconfirmed. At best it is hypothesized that Se mirrors sulfur and forms similar minerals. The most thermodynamically stable iron selenide is ferroselite, which is isostructural to marcasite. Yakolvleva et al. used SEM-EDS, powder XRD, and bulk chemical analysis to study a supergene iron sulpho-selenide deposit associated with a massive copper-zinc sulfide ore deposit. They discovered a solid solution series of pyrite-dzharkenite (FeS2–FeSe2) occurring in the weathering zone, suggesting that dzharkenite is a stable end member. Pyrite is the dominant iron sulfide mineral in the middle waste shale samples. Based on the micro-SXRF analyses of the freshly fractured rock surfaces, we hypothesize that most of the Se is associated with iron and S in the Permian sedimentary formations. Additional information on Se speciation in the shale requires spectroscopy to probe Se oxidation state and molecular coordination.
C. Bulk XANES
Despite the fact that the bulk samples have three degrees of assigned weathering, least, less, and more, there was no clear trend in Se oxidation state among the samples. It was expected that the least weathered sample would contain reduced species, the more weathered samples would contain oxidized species, and the less weathered samples would contain a combination of reduced and oxidized species. The presence of oxidized species in the least weathered sample may be an artifact from sample collection (drilling), preparation (grinding), and/or storage in aerobic environments. Alternatively the least weathered samples may actually have some unexpected weathering. Knudsen and Gunter and Herring and Grauch proposed that the enrichment of Se in the middle waste shale is a result of weathering of Se from the overlying materials and precipitation in the less permeable middle waste shale, making it possible for Se in the least weathered material to exist in the mobile, oxyanionic, Se(IV) and Se(VI) species. Due to the presence of organic matter, microbes, and lack of oxygen, diagenic elemental and metal selenides would be expected, however, the presence of some nonreduced Se cannot be ruled out. Perkins and Foster also observed Se(IV) species in weathered samples.
XANES from spot g is similar to seleniferous pyrite/marcasite as identified by Perkins and Foster. This standard was not available to make a direct comparison, however qualitatively the features in the XANES from spot g are the same as those for the Se-pyrite/marcasite standard reported in Perkins and Foster. Perkins and Foster used bulk XANES analysis to identify seleniferous pyrite/macrcasite in the middle waste shale samples from WPRA. Howard reported that when the concentration ratio of H2Se is high relative to H2S, and excess iron is present, ferroselite is ther-modynamically stable relative to elemental Se and pyrite. However, Howard[31, 36] suggested that when the concentration of H2Se is below a critical point, pyrite is stable, and Se will form either elemental Se or a Se–pyrite solid solution.
The micro-XANES spectrum from spot b, which is a spectral end member, has XANES features that most closely match the spectrum of seleno-cystine [Fig. 9(b)]. Black shales contain high concentrations of organic carbon compounds, including sulfur compounds. The similarity in the XANES from spot b with seleno-cystine is likely a result of a di-sulfide-type bond involving either two Se atoms or a Se and S atom, as in seleno-cystine, but does not necessarily indicate that seleno-cystine is the molecular form of Se in the shale. In fact, the EXAFS structure directly following the XANES region of spot b (data not shown) is clearly distinct from seleno-cystine and elemental Se. The concomitant presence of iron at this spot suggests that the Se may be present bonded into an organic molecule mixed with colloidal iron sulfide minerals, or possibly some type of iron-selenide carbon compound, such as a hydrogenase.
There are three possible mechanisms through which organic Se can form in the shales: (1) detrital Se introduced from plant materials, particularly since plants tend to shuttle Se into amino acids such as cystine/cysteine for detoxification, possibly altered by catagenic processes; (2) abiotic, nucleophillic replacement of organic functional groups by reduced Se(-II); or (3) microbial by-products that formed in the anaerobic sedimentary environment. Herring and Grauch used factor analysis to show that the Se in the Meade Peak Formation materials was highly correlated with organic C and total S, which supports our observation that an organically-bonded Se compound is present in the shale. Bond performed total Se analyses on methylene chloride (MeCl) extracts of middle waste shale from the WPRA with varying amounts of weathering. Methylene chloride is a selective extractant for semi-volatile organic compounds. The percentage of MeCl extractable Se in the bulk rock ranged from 0.4% to 2.6%, with the unweathered sample containing the most MeCl extractable Se. Additional analyses of the extracts by gas chromatography mass spec-troscopy indicated that C8 to C25 hydrocarbons were present.
Spot i has a XANES spectra similar to that of ferroselite [Fig. 9(b)]. As discussed earlier, ferroselite is considered to be a stable mineral formed when reduced Se and Fe concentrations are high relative to sulfide concentrations. Of the 12 spots measured on the two samples, only one had XANES that resembled the XANES from ferroselite. The other spots had XANES spectra indicative of mixtures of Se-substituted pyrite/marcasite and a diselenide organic compound [see Fig. 9(a)]. Although marcasite and ferroselite are isostructural, their unit cell dimensions differ, and thus should result in distinct XANES from Se-substituted pyrite/marcasite (see Table IV for theoretical interatomic distances). Dzharkenite, an iron selenide isostructural to pyrite, or achavalite, an iron selenide isostructural to mackinawite, were not available as standard references. Iron selenides are commonly found together with iron sulfides, and the geochemistry of Se is largely controlled by Fe. The kinetic and thermodynamic processes associated with diagenesis of these reduced iron selenide minerals are poorly understood, but related to temperature, redox potential, relative concentrations, and speciation. A summary reaction for the formation of ferroselite is
FeS2 + 2H2Se↔FeSe2 + 2H2S.
Howard proposed that for this reaction to occur the H2S/H2Se ratio must be less than 10-2.6, and suggested that elemental Se is a reactant. Elemental Se was not detected in our samples or the XANES spectra of the Perkins and Foster study. Grauch et al. identified elemental Se in the middle waste shale using SEM-EDS. Although Howard proposed possible reaction mechanisms for the formation of ferroselite in a pyrite environment, the actual reactions occurring in the environment are not known; highlighting the need for accurate speciation studies.
Analysis of the Se XANES indicated that there are three species present in the middle waste shale: selenium substituted pyrite/marcasite, organic diselenide compounds, and an iron selenide such as ferroselite. We found no evidence of elemental Se or oxidized Se in the freshly cleaved shale samples. Previous micro-analyses have suggested that Zn and Cu selenides are also present in the middle waste shale. We did not detect these species; however, there were positive correlations with these elements in the micro-SXRF mapping. The lack of detection of these metal selenides in the XANES data may be due to the low concentrations present relative to the iron selenides, i.e., the Zn or Cu selenides may be present as small isolates in the iron-sulfide/selenides.
E. Micro-EXAFS spectroscopy
To further identify the mineralogy of the selenium present in the middle waste shale we collected micro-EXAFS spectra on end-member spots i and g. A limited range of micro-EXAFS data were collected for spot b because it was not recognized as an end member at the time of data collection. While the near edge spectra (XANES) are sensitive to oxidation state and first shell coordination environment, EXAFS spectra are sensitive to backscattering off of multiple atomic shells, making this technique a sensitive probe for molecular composition and structure.
To further elucidate the local atomic environment surrounding Se the χ spectra for spots i and g were Fourier transformed (FT) to isolate the phase and amplitude contributions. Peaks in the FT correspond to atomic backscatterers surrounding the central Se atom (offset by a phase shift). In cases where multiple atoms exist at similar bond distances, separation of the backscatterers by FT is difficult. For example, in the mineral ferroselite the first atomic shell surrounding Se has four atoms between 2.34 and 2.5 Å (three Fe atoms and one Se atom, Table IV). In the FT of the ferroselite EXAFS spectra the peak at 2.05 Å corresponds to backscattering off of these atoms. The lack of resolution of the two atomic backscatterers in the FT illustrates the limitations of resolving molecular structures that have multiple shells with short distances from the central atom.
Fit results for sample 2 POI g and i and ferroselite and elemental Se (monoclinic) standards
ΔR (Å)e (% difference)
POI g sample 2 fit with Se-substituted pyrite model
POI i sample 2 fit with dzharkenite model
Elemental Se (monoclinic)
Theoretical bond distances and coordination numbers for Se and S-bearing minerals
2.34–2.38 ( = 2.37)
3.29–3.58 ( = 3.39)
3.74–3.91 ( = 3.85)
3.29–3.55 ( = 3.42)
3.67–3.87 ( = 3.72)
3.08–3.32 ( = 3.20)
3.43–3.61 ( = 3.48)
2.23–2.25 ( = 2.24)
2.97–3.38 ( = 3.19)
3.48–3.60 ( = 3.56)
The EXAFS from spot g was fit using theoretical backscatterers from pyrite. The model fit resulted in four atomic shells (Table III). The first atomic shell surrounding Se contains three Fe atoms at 2.25 Å. Attempts to fit S in this shell with and without iron were unsuccessful. One Se atom at 2.25 Å was fit, but did not improve the quality of the fit, suggesting it has a minimal contribution. Ideal pyrite has three Fe atoms at 2.26 Å and one S atom at 2.14 Å. Marcasite has a similar first atomic shell structure with S atoms at 2.21 and Fe between 2.23 and 2.25 Å. The third shell was fit with 12 Se atoms at 3.24 Å, which is close to the weighted average for these twelve atoms in pyrite or marcasite (Table IV), but distinctly shorter than the expected distance for dzharkenite or ferroselite. The fourth atomic shell was fit with four Fe atoms at 3.49 Å, which is most similar to the bond distance observed in pyrite. Thus we hypothesize that the mineral at spot g is a Se-substituted pyrite based on three observations: (1) the first shell Fe-bond distance is distinctly shorter than the FeSe2 polymorphs; (2) the second and higher shell backscattering atoms occur at similar bond distances as marcasite/pyrite (a slightly better match for pyrite); and (3) the similarity of the XANES to the Se-substituted marcasite/ pyrite presented by Perkins and Foster.
Fitting the EXAFS data from spot i resulted in four atomic shells. The fitted bond distances for the first four atomic shells agreed with the mineral dzharkenite as opposed to ferroselite (Tables III and IV). Specifically, in ferroselite the first shell Se atom is at a much longer bond distance than the three iron atoms, while in dzharkenite this is reversed. The third and fourth shells were fit with Se and Fe atoms, respectively, resulting in bond distances consistent with the weighted average for these atomic shells in dzharkenite. Based on the fit of the first four atomic shells having similar bond distances as the theoretical structure of dzharkenite, we hypothesize that spot i is dzharkenite and not ferroselite.
Dzharkenite is a rare mineral with only a few known occurrences. However, it has been reported as a mineral in solid solution with pyrite. Thus, the mineral is likely present in the shale as a Fe–Se–S mineral, but with Se dominating the composition such that it has a dzharkenite structure rather than a Se-substituted pyrite or marcasite structure.
The analyses of the Se micro-EXAFS spectra provided insights into the identity of two of the three end members we observed. Kulp and Pratt analyzed selenium speciation in chalk and shale collected from South Dakota and Wyoming using selective sequential extractions. Based on their results they hypothesized that the Se was present in the rocks predominantly as organic Se compounds and di-selenide minerals (e.g., ferroselite). Therefore the macroscopic observations of Kulp and Pratt support our findings. By using a molecular level technique we have extended the general observations of Kulp and Pratt and identified specific mineral and molecular speciation of Se in the shale. Given the difficulty in interpreting mineral structures due to sample heterogeneity, microscopic synchrotron-based XRF, XANES, and EXAFS spectroscopy allows for unique insights into Se speciation within heterogeneous materials such as rocks. We know of no other technique better suited for investigating the speciation of elements with bulk concentrations too low for traditional techniques.
The most commonly reported forms of Se in shale are Se-substituted pyrite,[6, 1, 31, 42] or precipitation of metal selenides. However, results of past studies are based on standard micro-analytical techniques, bulk EXAFS spectroscopy, and chemical extractions. In this study we used advanced microscopic spectroscopy to resolve the specific speciation of Se that occurs in shale. Our results suggest that there are three distinct phases in the shale that contain Se; dzharkenite, an organic di-selenide, and Se-substituted pyrite. In all of these phases Se existed in negative oxidation states, indicating that the oxidized Se present in weathered shales is a result of oxidative weathering of the reduced mineral Se to elemental Se, selenite, or selenate.
Improved knowledge of the forms of Se present in shale will allow for a better understanding of the biogeochemical cycling of Se as the Se-bearing rocks begin to weather. Oxidative breakdown of the various selenides occurring in the shale will lead to increased mobility and bioavailability of Se. These reactions are likely rate controlled. We are not aware of differences in the solubility or dissolution rates of dzharkenite, Se-substituted pyrite, and organic selenide compounds. However, by knowing the Se minerals and compounds that are relevant, further studies can be undertaken to determine how Se is released from Se-bearing mine waste. Such studies are important in assessing environmental risks associated with mining areas in the WPRA, where Se-bearing rocks have been displaced from the stable, slowly weathering strata, to the surficial environment in which weathering processes are much more rapid.
Funding for this project was provided by US-EPA grant ID#X970339010 and by a grant from the J.R. Simplot Company. The Advanced Light Source is supported by the Director, Office of Science, Office of Basic Energy Sciences, Materials Sciences Division, of the U.S. Department of Energy under Contract No. DE-AC03-76SF00098 at Lawrence Berkeley National Laboratory. Research carried out (in whole or in part) at the National Synchrotron Light Source, Brookhaven National Laboratory is supported by the U.S. Department of Energy, Division of Materials Sciences and Division of Chemical Sciences, under Contract No. DE-AC02-98CH10886. Assistance from Brian Hart and Sirine Fakara is appreciated.
- Mislevy P, Blue WG, Strieker JA, Cook BC, Vice MJ: Reclamation of Drastically Disturbed Lands. Edited by: Barnhisel RI, Darmody RG, Daniels WL. 2000, American Society of Agronomy, Madison, WI, 41: 961.Google Scholar
- Mansfield GR: U.S. Geol Surv Prof Pap. 152. 1927, U.S. Government Print OfficeGoogle Scholar
- Knudsen AC: 2002, Ph.D. dissertation, University of IdahoGoogle Scholar
- Desborough G, DeWitt E, Jones J, Meier A, Meeker G: U.S. Geol Surv Open-File Rep No. 99-129. 1999, 20.Google Scholar
- Schrauzer GN: J Nutr. 2000, 130: 1653.Google Scholar
- Presser TS, Swain WC: Appl Geochem. 1990, 5: 703-10.1016/0883-2927(90)90066-E.View ArticleGoogle Scholar
- Zawislanski PT, Benson SM, Terberg R, Borglin SE: Environ Sci Technol. 2003, 37: 2415-10.1021/es020977z.View ArticleGoogle Scholar
- Workman SM: Reclaiming Mine Soils and Overburden in the Western United States: Analytic Parameters and Procedures. Edited by: Williams RD, Schuman GE. 1987, Soil Conservation Society of America, Ankeny, IA, 109-Report No. AG 680, 1983; S. E. Fisher, F. F. Munshower, and F. ParadyGoogle Scholar
- Hayes KF, Roe AL, Brown GE, Hodgson KO, Leckie JO, Parks GA: Science. 1987, 238: 783.View ArticleGoogle Scholar
- Desborough GA: U.S. Geol Surv Open-File Rep No. 77-341. 1977, 20.Google Scholar
- Perkins RB, Foster AL: Life Cycle of the Phosphoria Formation: From Deposition to Post-Mining Environment. Edited by: Hein JR. 2004, Elsevier, Amsterdam, 8: 251.View ArticleGoogle Scholar
- Pickering IJ, Brown GE, Tokunaga TK: Environ Sci Technol. 1995, 29: 2456.View ArticleGoogle Scholar
- Manceau A, Marcus MA, Tamura N, et al: Applications of Synchrotron Radiation in Low-Temperature Geochemistry and Environmental Sciences. Edited by: Fenter PA, Rivers ML, Sturchio NC. 2002, Geochemical Society and Mineralogical Society of America, Washington, DC, 49: 341.Google Scholar
- Strawn DG, Doner HE, Zavarin M, McHugo S: Geoderma. 2002, 108: 237-10.1016/S0016-7061(02)00133-7.View ArticleGoogle Scholar
- Tokunaga TK, Sutton SR, Bujt S, Nuessle P, Shea-McCarthy G: Environ Sci Technol. 1998, 32: 1092-10.1021/es9706559.View ArticleGoogle Scholar
- Herring JR, Desborough GA, Wilson SA, Tysdal RG, Grauch RI, Gunter ME: U.S. Geol Surv Open-File Rep No. 99-147-A. 1999, 24.Google Scholar
- Herring JR, Grauch RI, Desborough GA, Wilson SA, Tysdal RG: U.S. Geol Surv Open-File Rep No. 99-147-C. 2000, 35.Google Scholar
- Herring JR, Grauch RI, Desborough GA, Wilson SA, Tysdal RG: U.S. Geol Surv Open-File Rep No. 99-147-D. 2000, 38.Google Scholar
- Herring JR, Wilson SA, Stillings LA, Knudsen AC, Gunter ME, Tysdal RG, Grauch RI, Desborough GA, Zielinski RA: U.S. Geol Surv Open-File Rep No. 99-147-B. 2000, 34.Google Scholar
- Warren CG: Econ Geol. 1968, 63: 418.View ArticleGoogle Scholar
- Marcus MA, MacDowell AA, Celestre R, Manceau A, Miller T, Padmore HA, Sublett RE: J Synchrotron Radiat. 2004, 11: 239-10.1107/S0909049504005837.View ArticleGoogle Scholar
- Lytle FW, Greegor RB, Sandstom DR, Marques EC, Wond J, Spiro CL, Huffman GP, Huggins FE: Nucl Instrum Methods Phys Res A. 1984, 226: 542-10.1016/0168-9002(84)90077-9.View ArticleGoogle Scholar
- Newville J: J Synchrotron Radiat. 2001, 8: 322-10.1107/S0909049500016964.View ArticleGoogle Scholar
- Ankudinov AL, Ravel B, Rehr JJ, Conradson SD: Phys Rev B. 1998, 58: 7565-10.1103/PhysRevB.58.7565.View ArticleGoogle Scholar
- Ressler T: J Synchrotron Radiat. 1998, 5: 118-10.1107/S0909049597019298.View ArticleGoogle Scholar
- Bostick BC, Fendorf S: Geochim Cosmochim Acta. 2003, 67: 909-10.1016/S0016-7037(02)01170-5.View ArticleGoogle Scholar
- Beauchemin S, Hesterberg D, Beauchemin M: Soil Sci Soc Am J. 2002, 66: 83.View ArticleGoogle Scholar
- Frank P, Heman B, Carlson RMK, Hodgson KO: Inorg Chem. 1994, 33: 3794-10.1021/ic00095a024. K. Xia, F. Weesner, W. F. Bleam, P. R. Bloom, U. L. Skyllberg, and P. A. Helmke, Soil Sci. Soc. Am. J. 62, 1240 (1998).View ArticleGoogle Scholar
- Bond MM: Master of Science thesis. 1999, University of IdahoGoogle Scholar
- Munkers JP: Master of Science thesis. 2000, University of IdahoGoogle Scholar
- Howard JH: Geochim Cosmochim Acta. 1977, 41: 1665-10.1016/0016-7037(77)90176-4.View ArticleGoogle Scholar
- Yakovleva VA, Belogub EV, Novoselov KA: Miner Mag. 2003, 67: 355-10.1180/0026461036720106.View ArticleGoogle Scholar
- Knudsen AC, Gunter ME: Phosphates-Geochemical, Geobiological, and Materials Importance. Edited by: Kohn MJ, Rakovan J, Hughes JM. 2003, Mineralogical Society of America, Washington, DCGoogle Scholar
- Herring JR, Grauch RI: Life Cycle of the Phosphoria Formation. Edited by: Hein JR. 2004, Elsevier, Amsterdam, 8: 321.View ArticleGoogle Scholar
- Grauch RI, Desborough GA, Meeker GP, Foster AL, Tysdal RG, Herring JR, Lowers HA, Ball BA, Zielinski RA, Johnson EA: Life Cycle of the Phosphoria Formation: From Deposition to the Post-Mining Environment. Edited by: Hein JR. 2004, Elsevier, Amsterdam, 8: 189.View ArticleGoogle Scholar
- Coleman RG, Delevaux M: Econ Geol. 1957, 52: 499.View ArticleGoogle Scholar
- Passier HF, Lang GD: Proceedings of the Ocean Drilling Program, Scientific Results, 1998. Edited by: Orr WL, White CM. 1990, American Chemical Society, Washington, DC, 160: 249-M. L. Tuttle, C. A. Rice, and M. B. Goldhaber, in Geochemistry of Sulfur in Fossil Fuels, Vol. 429, 114 pp.Google Scholar
- Thompson JF, Smith IK, Madison JT: Sulfur in Agriculture. Edited by: Tabatabai MA. 1986, ASA, CSSA, and SSSA, Madison, WI, 27: 57.Google Scholar
- Vairavamurthy A, Mopper K: Biogenic Sulfur in the Environment. Edited by: Saltzman ES, Cooper WJ. 1989, American Chemical Society, Washington, DC, 393: 231.View ArticleGoogle Scholar
- Tracy ML, Möller G: J Assoc Off Anal Chem. 1990, 73: 404.Google Scholar
- Kulp TR, Pratt LM: Geochim Cosmochim Acta. 2004, 68: 3687-10.1016/j.gca.2004.03.008.View ArticleGoogle Scholar
- Velinsky DJ, Cutter GA: Anal Chim Acta. 1990, 235: 419-10.1016/S0003-2670(00)82102-9.View ArticleGoogle Scholar
- Wyckoff RWG: Crystal Structures. 1963, Wiley & Sons, New York, 1:Google Scholar
- Mandarino JA: Miner Record. 1996, 27: 115.Google Scholar