Selenium attenuation via reductive precipitation in unsaturated waste rock as a control on groundwater impacts in the Idaho phosphate patch

Phosphate mining in southeastern Idaho has historically resulted in the release of dissolved metals and inorganics to groundwater and surface water, primarily due to leachate from waste rock in backfilled pits and overburden storage piles. Selenium (Se) is of particular concern due to its high concentration in leachate and its limited attenuation downgradient of source zones under oxic conditions. Assessments of potential groundwater/surface water impacts from waste rock typically involve laboratory characterization using saturated and unsaturated flow columns packed with waste rock. In this study, we compare the results of saturated and unsaturated column tests with groundwater quality data from the Mountain Fuel, Champ, South and Central Rasmussen Ridge Area (SCRRA), Smoky Canyon, Ballard, Henry, and Enoch Valley Mines, to understand the release and attenuation of Se in different geochemical environments. Column studies and field results demonstrate that the ratio of aqueous Se to aqueous sulfate (Se:SO₄ ratio) is a useful metric for understanding Se release and attenuation, where the extent of sulfate reduction is much less than Se reduction. Comparison of dissolved Se and sulfate results suggests that the net leachability of Se from unsaturated waste rock is variable. Overall, Se concentrations in groundwater directly beneath waste rock dumps is not as high as would be predicted from unsaturated columns. Lower Se:SO₄ ratios are observed immediately beneath waste rock dumps and backfilled pits relative to areas receiving shallow waste rock runoff. It is hypothesized that Se released in the oxic upper portions of the waste rock is subsequently attenuated via reductive precipitation at depth in unsaturated, low-oxygen portions of the waste rock. This highlights an important mechanism by which Se may be naturally attenuated within waste rock piles prior to discharge to groundwater and surface water. These results have important implications for mining practices in the region. A better understanding of Se dynamics can help drive waste rock management during active mining and capping/water management options during post-mining reclamation.     

In situ formation of welded tuff-like textures in the carapace of a voluminous silicic lava flow,Owyhee County,SW Idaho

 The Badlands rhyolite, on the Owyhee Plateau of southwestern Idaho, can be demonstrated to be a large lava flow on the basis of its geometry of large and small flow lobes, its well-exposed near-vent features, and its response to pre-existing topography. However, samples of the dense upper vitrophyre of the unit reveal a range of annealed fragmental textures, including material which closely resembles the compressed, welded glass shards which are characteristic of ignimbrites. Formation of these tuff-like textures involved processes probably common to emplacement of most silicic lava flow units. Decompression upon extrusion causes inflation of pumice at the surface of the lava flow; some of this pumice is subsequently comminuted, producing loose bubble-wall shards, bits of pumice, chips of dense glass, and fragments of phenocrysts. This debris sifts down around loose blocks and into open fractures deeper in the flow, where it can be reheated, compressed, and annealed to varying degrees. The end result is a dense vitrophyre layer (beneath the true upper, non-welded carapace breccia) which can be extremely texturally heterogeneous, with areas of flow-foliated lava occurring very near lava which in many aspects looks like welded ignimbrite, complete with flattened pumices. Identical textures in other silicic units have been cited by previous workers as evidence that those units erupted as pyroclastic flows which then underwent sufficient rheomorphism to create a flow-foliated rock which otherwise appears to be lava. The textures described herein indicate that lava flows can come to mimic rheomorphic ignimbrites, at least at scales ranging from thin sections to outcrops. Voluminous silicic units with scattered fragmental textures, but with otherwise lava-like features, are probably true effusive lava flows. Received: January 30, 1995 / Accepted: January 22, 1996     

Chemical and physical properties affecting strontium distribution coefficients of surficial-sediment samples at the Idaho National Engineering and Environmental Laboratory, Idaho

 The U.S. Geological Survey and Idaho State University, in cooperation with the U.S. Department of Energy, conducted a study to determine strontium distribution coefficients (K_ds) of surficial sediments at the Idaho National Engineering and Environmental Laboratory (INEEL). Batch experiments using synthesized aqueous solutions were used to determine K_ds, which describe the distribution of a solute between the solution and solid phase, of 20 surficial-sediment samples from the INEEL. The K_ds for the 20 surficial-sediment samples ranged from 36 to 275 ml/g. Many properties of both the synthesized aqueous solutions and sediments used in the experiments also were determined. Solution properties determined were initial and equilibrium concentrations of calcium, magnesium, and strontium, pH and specific conductance, and initial concentrations of potassium and sodium. Sediment properties determined were grain-size distribution, bulk mineralogy, whole-rock major-oxide and strontium and barium concentrations, and Brunauer-Emmett-Teller (BET) surface area. Solution and sediment properties were correlated with strontium K_ds of the 20 surficial sediments using Pearson correlation coefficients. Solution properties with the strongest correlations with strontium K_ds were equilibrium pH and equilibrium calcium concentration correlation coefficients, 0.6598 and –0.6518, respectively. Sediment properties with the strongest correlations with strontium K_ds were manganese oxide (MnO), BET surface area, and the >4.75-mm-grain-size fraction correlation coefficients, 0.7054, 0.7022, and –0.6660, respectively. Effects of solution properties on strontium K_ds were interpreted as being due to competition among similarly charged and sized cations in solution for strontium-sorption sites; effects of sediment properties on strontium K_ds were interpreted as being surface-area related. Multivariate analyses of these solution and sediment properties resulted in r² values of 0.8071 when all five properties were used and 0.8043 when three properties, equilibrium pH, MnO, and BET surface area, were used. Received: 30 November 1998 · Accepted: 16 February 1999     

Modeling changes in rill erodibility and critical shear stress on native surface roads

This study investigated the effect of cumulative overland flow on rill erodibility and critical shear stress on native surface roads in central Idaho. Rill erodibility decreased exponentially with increasing cumulative overland flow depth; however, critical shear stress did not change. The study demonstrated that road erodibility on the studied road changes over the course of one or more consecutive overland flow events. Therefore, model simulations that fail to take into consideration this change will probably over-estimate sediment yields. An exponential function describing the relationship between rill erodibility and cumulative overland flow depth is presented as a basis for future model development for simulating erosion on native surface roads. Copyright © 2008 John Wiley & Sons, Ltd.     

Detrital zircon U–Pb reconnaissance of the Franciscan subduction complex in northwestern California

In northwestern California, the Franciscan subduction complex has been subdivided into seven major tectonostratigraphic units. We report U-Pb ages of ≈2400 detrital zircon grains from 26 sandstone samples from 5 of these units. Here, we tabulate each unit’s interpreted predominant sediment source areas and depositional age range, ordered from the oldest to the youngest unit. (1) Yolla Bolly terrane: nearby Sierra Nevada batholith (SNB); ca. 118 to 98 Ma. Rare fossils had indicated that this unit was mostly 151–137 Ma, but it is mostly much younger. (2) Central Belt: SNB; ca. 103 to 53 Ma (but poorly constrained), again mostly younger than previously thought. (3) Yager terrane: distant Idaho batholith (IB); ca. 52 to 50 Ma. Much of the Yager’s detritus was shed during major core complex extension and erosion in Idaho that started 53 Ma. An Eocene Princeton River–Princeton submarine canyon system transported this detritus to the Great Valley forearc basin and thence to the Franciscan trench. (4) Coastal terrane: mostly IB, ±SNB, ±nearby Cascade arc, ±Nevada Cenozoic ignimbrite belt; 52 to <32 Ma. (5) King Range terrane: dominated by IB and SNB zircons; parts 16–14 Ma based on microfossils. Overall, some Franciscan units are younger than previously thought, making them more compatible with models for the growth of subduction complexes by progressive accretion. From ca. 118 to 70 Ma, Franciscan sediments were sourced mainly from the nearby Sierra Nevada region and were isolated from southwestern US and Mexican sources. From 53 to 49 Ma, the Franciscan was sourced from both Idaho and the Sierra Nevada. By 37–32 Ma, input from Idaho had ceased. The influx from Idaho probably reflects major tectonism in Idaho, Oregon, and Washington, plus development of a through-going Princeton River to California, rather than radical changes in the subduction system at the Franciscan trench itself.     

Geochemistry and stratigraphic correlation of basalt lavas beneath the Idaho Chemical Processing Plant, Idaho National Engineering Laboratory

 Thirty-nine samples of basaltic core were collected from wells 121 and 123, located approximately 1.8 km apart north and south of the Idaho Chemical Processing Plant at the Idaho National Engineering Laboratory. Samples were collected from depths ranging from 15 to 221 m below land surface for the purpose of establishing stratigraphic correlations between these two wells. Elemental analyses indicate that the basalts consist of three principal chemical types. Two of these types are each represented by a single basalt flow in each well. The third chemical type is represented by many basalt flows and includes a broad range of chemical compositions that is distinguished from the other two types. Basalt flows within the third type were identified by hierarchical K-cluster analysis of 14 representative elements: Fe, Ca, K, Na, Sc, Co, La, Ce, Sm, Eu, Yb, Hf, Ta, and Th. Cluster analyses indicate correlations of basalt flows between wells 121 and 123 at depths of approximately 38–40 m, 125–128 m, 131–137 m, 149–158 m, and 183–198 m. Probable correlations also are indicated for at least seven other depth intervals. Basalt flows in several depth intervals do not correlate on the basis of chemical compositions, thus reflecting possible flow margins in the sequence between the wells. Multi-element chemical data provide a useful method for determining stratigraphic correlations of basalt in the upper 1–2 km of the eastern Snake River Plain. Received: 16 February 1996 · Accepted: 1 April 1996     

Re-Os systematics and geochemistry of cobaltite (CoAsS) in the Idaho cobalt belt,Belt-Purcell Basin,USA: Evidence for middle Mesoproterozoic sediment-hosted Co-Cu sulfide mineralization with Grenvillian and Cretaceous remobilization

We report the first study of the Re-Os systematics of cobaltite (CoAsS) using disseminated grains and massive sulfides from samples of two breccia-type and two stratabound deposits in the Co-Cu-Au Idaho cobalt belt (ICB), Lemhi subbasin to the Belt-Purcell Basin, Idaho, USA. Using a ¹⁸⁵Re + ¹⁹⁰Os spike solution, magnetic and non-magnetic fractions of cobaltite mineral separates give reproducible Re-Os analytical data for aliquot sizes of 150 to 200 mg. Cobaltite from the ICB has highly radiogenic ¹⁸⁷Os/¹⁸⁸Os ratios (17–45) and high ¹⁸⁷Re/¹⁸⁸Os ratios (600–1800) but low Re and total Os contents (ca. 0.4–4 ppb and 14–64 ppt, respectively). Containing 30 to 74% radiogenic ¹⁸⁷Os, cobaltite from the ICB is amenable to Re-Os age determination using the isochron regression approach.Re-Os data for disseminated cobaltite mineralization in a quartz-tourmaline breccia from the Haynes-Stellite deposit yield a Model 1 isochron age of 1349 ± 76 Ma (2σ, n = 4, mean squared weighted deviation MSWD = 2.1, initial ¹⁸⁷Os/¹⁸⁸Os ratio = 4.7 ± 2.2). This middle Mesoproterozoic age is preserved despite a possible metamorphic overprint or a pulse of metamorphic-hydrothermal remobilization of pre-existing cobaltite that formed along fold cleavages during the ca. 1190–1006 Ma Grenvillian orogeny. This phase of remobilization is tentatively identified by a Model 3 isochron age of 1132 ± 240 Ma (2σ, n = 7, MSWD = 9.3, initial ¹⁸⁷Os/¹⁸⁸Os ratio of 9.0 ± 2.9) for cobaltite in the quartz-tourmaline breccia from the Idaho zone in the Blackbird mine.All Mesoproterozoic cobaltite mineralization in the district was affected by greenschist- to lower amphibolite-facies (garnet zone) metamorphism during the Late Jurassic to Late Cretaceous Cordilleran orogeny. However, the fine- to coarse-grained massive cobaltite mineralization from the shear zone-hosted Chicago zone, Blackbird mine, is the only studied deposit that has severely disturbed Re-Os systematics with evidence for a linear trend of mixing with (metamorphic?) fluids.The new Re-Os ages and extremely high initial ¹⁸⁷Os/¹⁸⁸Os ratios of cobaltite reported here favor a magmatic-hydrothermal genetic model for a multi-stage REE-Y-Co-Cu-Au mineralization occurring at ca. 1370 to 1349 Ma, and related to the emplacement of the Big Deer Creek granite pluton at ca. 1377 Ma. In our model, deposition of paragenetically early xenotime and gadolinite was followed by an influx of Mesoproterozoic evaporitic brines and magmatic-hydrothermal fluids containing metals and reduced sulfur derived from mafic and oceanic island-arc Archean to Paleoproterozoic rocks in the Laurentian basement. Cobaltite mineralization occurred upon cooling of these fluids at an inferred temperature of 300 °C or below.     

Scale Dependence of Species Richness Maps*

Mapping patterns of species richness is a longstanding tradition in biogeography and more recently in conservation planning. This paper describes the effects of sampling unit size on patterns of vertebrate richness across landscapes in the Intermountain Sagebrush and Rocky Mountain Forest ecoregions in Idaho. Variability of richness decreased with increasing size of the spatial sampling units at both sites, with greater overall variability in the forested site. Richness in the sagebrush region was best explained by alpha diversity, whereas beta diversity accounted for more of the variability in richness in the forested site.     

Modeling of terracette‐hillslope soil moisture as a function of aspect,slope and vegetation in a semi‐arid environment

In the semi‐arid western United States, water availability plays a defining role in land use. Soil moisture, vegetation, and microtopography are key variables in the hydrologic function of these ecosystems. Previous research has not addressed the influence of site‐specific aspect, vegetation, or slope gradient on terracette soil moisture patterns in semi‐arid rangelands. Therefore, the objectives of this study were to: (1) assess the influence of terracette site aspect, vegetation cover, and slope on soil moisture; (2) conceptualize conditions at the hillslope scale given terracette morphology; and (3) estimate the extent of terracettes at a regional scale. The Simultaneous Heat and Water (SHAW) model was used to simulate soil water dynamics of terracettes given variations in site conditions. These results were coupled with time‐of‐flight laser scans to quantify terracette bench and riser percent‐area, and statewide assessments of terracette extent using digital orthoimagery and a geographical information system (GIS). Modeling results indicated site aspect had minimal influence (±0.005 m³ m^?3) on terracette soil moisture. Vegetation, represented as leaf area index (LAI), had the single‐most influential effect on terracette volumetric water content (θ _v) demonstrated by an inverse relationship of LAI to mean terracette hillslope θ _v; and slope increases of ≥15% on northern azimuths increased mean θ _v which contrasted with southern azimuths for similar slope increases. Laser scanning results indicated bench width and riser length could be estimated from mean site slope (R ² = 0.82 risers and R ² = 0.93 benches). Aerial orthoimagery/GIS assessments estimated >159 000 ha of terracettes throughout the State of Idaho, with >41 000 ha (~26%) occurring on lands managed as grazing allotments. These findings provide an increased understanding of rangeland hydrologic processes as influenced by cattle density, vegetation, and terracettes which can aide land managers in their selection and application of best management practices on these lands. Copyright © 2017 John Wiley & Sons, Ltd.     

Holocene vegetation, fire and climate history of the Sawtooth Range, central Idaho, USA

The paucity of low- and middle-elevation paleoecologic records in the Northern Rocky Mountains limits our ability to assess current environmental change in light of past conditions. A 10,500-yr-long vegetation, fire and climate history from Lower Decker Lake in the Sawtooth Range provides information from a new region. Initial forests dominated by pine and Douglas-fir were replaced by open Douglas-fir forest at 8420 cal yr BP, marking the onset of warmer conditions than present. Presence of closed Douglas-fir forest between 6000 and 2650 cal yr BP suggests heightened summer drought in the middle Holocene. Closed lodgepole pine forest developed at 2650 cal yr BP and fires became more frequent after 1450 cal yr BP. This shift from Douglas-fir to lodgepole pine forest was probably facilitated by a combination of cooler summers, cold winters, and more severe fires than before. Five drought episodes, including those at 8200 cal yr BP and during the Medieval Climate Anomaly, were registered by brief intervals of lodgepole pine decline, an increase in fire activity, and mistletoe infestation. The importance of a Holocene perspective when assessing the historical range of variability is illustrated by the striking difference between the modern forest and that which existed 3000 yr ago.