Possible Involvement of Permian Phosphoria Formation Oil as a Source of REE and Other Metals Associated with Complex U-V Mineralization in the Northern Bighorn Basin?

The origin of V, U, REE and other metals in the Permian Phosphoria Formation have been speculated and studied by numerous scientists. The exceptionally high concentrations of metals have been interpreted to reflect fundamental transitions from anoxic to oxic marine conditions. Much of the oil in the Bighorn Basin, is sourced by the Phosphoria Formation. Two of the top 10 producing oil fields in Wyoming are located approximately 50 km west of two abandoned U-V mining districts in the northern portion of the basin. These fields produce from basin margin anticlinal structures from Mississippian age reservoir rock. Samples collected from abandoned U-V mines and prospects hosted in Mississippian aged paleokarst in Montana and Wyoming have hydrocarbon residue present and contain anomalous high concentrations of many metals that are found in similar concentrations in the Phosphoria Formation. As, Hg, Mo, Pb, Tl, U, V and Zn, often metals of environmental concern occur in high concentrations in Phosphoria Formation samples and had values ranging from 30–1295 ppm As, 0.179–12.8 ppm Hg, 2–791 ppm Mo, <2–146 ppm Pb, 10–490 ppm Tl, 907–86,800 ppm U, 1240–18,900 ppm V, and 7–2230 ppm Zn, in mineralized samples from this study. The REE plus Y composition of Madison Limestone- and limestone breccia hosted-bitumen reflect similar patterns to both mineralized samples from this study and to U.S. Geological Survey rock samples from studies of the Phosphoria Formation. Geochemical, mineralogical and field data were used to investigate past theories for mineralization of these deposits to determine if U present in home wells and Hg content of fish from rivers on the proximal Crow Indian Reservation may have been derived from these deposits or related to their mode of mineralization.


Introduction
The original impetus for this study was to try and determine the source of Hg and Pb in the Bighorn River. The Bighorn River is listed as a 303d impaired water body downstream of the Crow Indian Reservation due to elevated Pb and Hg and the source of contamination is unknown [1]. In a 2015 United States Geological Survey (USGS) study, fish tissue from the Bighorn Reservoir had the highest Hg levels observed in the entire Yellowstone River Basin (0.676 µg/g wet weight) [2]. The Bighorn River and its tributaries drain both the Pryor Mountain U-V Mining District (PMD), Montana (MT) and the Little Mountain U-V Mining District (LMD), Wyoming (WY) (Figure 1).

Stratigraphy
Rocks in the region range from Archean to Quaternary in age ( Figure 2). Approximately 760-915 m of Paleozoic rock overlie Archean gneiss and schist in the study areas [27]. The main formations of exposure relevant to the U-V deposits in the mining districts include limestones and dolomites of the Mississippian Madison Group (Madison Formation in WY). Extensive regional paleokarst features at the top of the Madison Group locally host collapse breccia [28,29], which are host structures of the U-V deposits [5,6,10,11]. The horizon is characterized by solution cavities along bedding planes, fractures and joints, and collapse breccia filled with angular limestone fragments of variable size (from several cm to several meters across) and reddish clay silt from the overlying Amsden Formation [6]. The in filling with Amsden Formation imparts a red to purple matrix to many of the breccia units.

Structural Setting
During the Laramide Orogeny , the central Rocky Mountain region experienced a general NE-SW (40 • -55 • ) transpressional strain regime as a result of the low-angle subduction of the Farallon plate at its western continental margin [31,32]. Beginning in the late Cretaceous, Laramide shortening was accommodated by uplift of broad basement arches that eventually isolated the Bighorn Basin ( Figure 3) [31]. The Bighorn Basin has a generally NW-SE sinuous axis [26]. Smaller scale folds in the form of anticlines and synclines formed along the basin margins and served to accommodate regional shortening [31][32][33]. These folds host most of the oil fields of the Bighorn Basin [34]. The structures hosting the U-V mining districts and the Little Sheep Mountain anticline are all folds created during this time. Structurally the Pryor Mountains are divided into N and S segments by the Sage Creek fault zone, the eastern extension of the Nye-Bowler lineament. The lineament is interpreted to be the expression of a left-lateral wrench fault at depth and defines the northern end of the Bighorn Basin [26,34,35]. The U-V deposits in the PMD occur in the southern portion of the Pryor Mountains, MT. The location of the mining districts is outlined by the red rectangle. Map modified from Blackstone [26] and Stahl [36].

Mineralization and Alteration
The main difference in the ore of the two districts is the silica content. The silica content of the PMD was high and the LMD was low [5,6,11]. Abundant silica observed in the PMD is concentrated in and around siliceous tectonic hydrothermal breccia (THB). These were emplaced in fractures along the crest of the Gypsum Creek anticline which forms the southern portion of the Big Pryor Mountain block ( Figure 1). The U-V ore minerals of MT/WY are easy to recognize in outcrop due to their bright yellow color ( Figure 4A-G). Boxwork is common in both districts ( Figure 4H,I). At many of the deposits, the host rock appears fresh and unaltered within a few feet of the mineralization. The wall rock in the PMD has the most alteration of the three areas studied including bleaching, liesegang banding, silicification, hematite-, and limonite-staining. Tectonic Hydrothermal Breccia, PMD Siliceous THB only observed in the PMD, are relatively small, 0.5-5 m wide, hematite-and limonite-stained dike-like features ( Figure 5A). Tectonic breccia were sampled to determine the isotope markers (stable C, O and radiogenic Sr) [37]. Although bleaching does not extend far from the THB, the wall rock is generally bleached where the breccia cut through. Wall rock around some of the features is liesegang banded ( Figure 5B). The matrices of the THB are composed of cryptocrystalline silica and local microscopic Fe-oxide inclusions. Unmineralized THB have barite, local fluorite, minor oxides, and sulfides along fracture surfaces and in vugs. The mineralized THB have U-V minerals filling in vugs and fractures in addition to the other minerals.

Field Work
Samples were chosen to test theories proposed for the origin of these deposits. A Thermo Scientific™ RadEye™ Personal Radiation Detector (Thermo Fisher Scientific Inc., Erlangen, Germany) was used to identify elevated radiation sources in the field. Rock samples included breccia clasts and matrix samples, altered and unaltered host rock, generations of cements from veins, and mineralized samples from dumps and outcrops. A bentonite sample was collected from an operating bentonite mine approximately 20 km southeast of the LMD to use to test the theory that U leached from volcanic ash that once covered the region. Dried bitumen was collected from fractures in Mm limestone and limestone breccia in two caves as well as from fractures in a limestone solution breccia in the Permian Goose Egg Formation exposed in the river cut of Little Sheep Mountain ( Figure 1). The bitumen samples were collected to compare the REE and Y composition of our U-V mineralized samples to test the hydrothermal theory, which is an association with oil and brines from the BHB. Two speleothem calcite samples from Big Pryor Mountain were collected to compare REE and assay analyses with calcite cement of breccia samples and calcite fracture fill to test the meteoric theory.

Mineralogical and Geochemical Characterization
Work performed to complete analyses included: standard petrography, cathodoluminescence (CL), scanning electron microscopy (SEM) utilizing energy dispersive X-ray spectrometry (EDS) and backscatter electron imaging (BSE), powder X-ray diffraction (XRD) spectroscopy and minor fluid inclusion work. Most of the analytical work was performed using instrumentation in the Imaging and Chemical Analysis Laboratory (ICAL) facility at Montana State University (MSU), Bozeman, MT, USA. XRD work was performed on a SCINTAG X1 Diffraction Spectrometer (Scintag Inc., West Palm Beach, FL, USA) and computer-aided mineral identification system at ICAL. SEM work was conducted at ICAL using a JEOL JSM-6100 (JEOL USA Inc., Peabody, MA, USA). CL was performed at the Department of Earth Sciences (Montana State University, Bozeman, MT, USA) using a RELIOTRON CL luminoscope instrument (RELION Industries LLC, Bedford, MA, USA) attached to a Nikon microscope (Nikon Corporation, Tokyo, Japan). Assays were performed by two different commercial laboratories: American Analytical Services (AAS) Inc., in Osborn, ID, USA and an Australian Laboratory Services (ALS) Global branch out of Reno, NV, USA.

Sample Preparation
Samples were prepared for chemical assay and REE analyses at MSU. Most samples were broken into chips and then pulverized in a Plattner's diamond mortar and pestle set while a few required micro drilling with a diamond tipped dremel bit. Some of the samples for REE analysis were pulverized in a shatter box. The mineralized U-V breccia were selected to try to get mainly the yellow U-V minerals, samples included part of the breccia host as the U-V minerals were disseminated throughout the samples. Fluorite was microdrilled. Barite was collected as individual crystals in the field and pulverized at MSU. Thin sections and some hand samples were stained for carbonate identification using a mixture of alizarin red and potassium cyanide solution [38].

Results, Mineralogy and Geochemistry
Samples examined by multiple methods revealed the episodic nature of mineralization in both districts. Bitumen was observed at all scales in this study; in outcrop along fractures and bedding planes, petrographically in interparticle porosity, and as inclusions in cement and minerals. Oxides, hydroxides and sulfides were closely associated with U-V minerals in both districts. The white card petrographic method, placing a piece of white paper beneath a petrographic slide, differentiated bitumen from sulfides and oxides [39]. It was also useful for viewing U-V minerals. CL revealed zoning in carbonates, fluorites, and highlighted U-V minerals and Fe-V oxides. The EDS and BSE modes of the SEM revealed compositional variations in minerals. The Secondary Electron Image (SEI) mode of the SEM aided in viewing morphology of minerals with the same composition. BSE mode discriminates between different phases as a function of the mean atomic number [40]. Minerals with larger Z appear brighter than lighter ones. This mode was useful for seeing zonation in some of the sulfides and for identifying elements present in minerals. Portable XRF and EDS detected Au (Appendix A, Figure A1) and platinum group metals (PGM) in some of the samples. A few samples were sent in for assays for these elements (Table A1) The primary ore minerals in the two districts are the Ca-uranyl vanadate minerals Tyuyamunite (Ca(UO 2 ) 2 (VO 4 ) 2 ·5-8H 2 O) and metatyuyamunite (Ca(UO 2 ) 2 (VO 4 ) 2 ·3-5H 2 O). These were the only uranyl vanadates reported in the literature and were reconfirmed in both districts in this study by XRD and EDS. Two additional metal-uranyl vanadates were identified by XRD and EDS in only the PMD in this study [41]. These minerals are francevillites in the Curienite-Francevillite Series (Pb(UO 2 ) 2 (VO 4 ) 2 ·5H 2 O to Ba(UO 2 ) 2 (VO 4 ) 2 ·5H 2 O) [42]. Additionally, Carnotite, a K-uranyl vanadate (K 2 (UO 2 ) 2 (VO 4 ) 2 ·3H 2 O), [43], was identified only by EDS in this study in the PMD ( Figure A2). U-K grains with occur locally as 2-10 µm grains in dark purple fluorite ( Figure A3). EDS detected thallium (Tl) in an unidentified uranyl vanadate from a prospect site, in the PMD ( Figure A4). A mineral with the same morphology was found with tyuyamunite from the Dandy mine, which had the highest Tl (490 ppm) concentration of the study in analyses performed by ALS (Appendix A, Table A2). The only Tl-uranyl vanadate found in Mindat.org is thallian carnotite (K,Tl) 2 (UO 2 ) 2 [VO 4 ] 2 ·3H 2 O [44]. This mineral did not contain K and may be a Tl-bearing phase of curientite, with Tl substituting for some of the Pb because Tl and Pb have similar atomic and covalent radii: Tl (1.71.6 Å and 1.45 ± 7 Å) and Pb (1.75 Å and 1.46 ± 5 Å) [45]. EDS detected Tl in another uranyl vanadate from a sample collected from a stockpile area near the Dandy mine ( Figure A5). The weight percent of the elements are similar to francevillite; the mineral may be a Tl-bearing francevillite with Tl substituting for Pb.
Minerals were cross-referenced with information in the volume Reviews in Mineralogy, Uranium: Mineralogy, Geochemistry and the Environment [45] and internet sources such as Mindat.org, handbookofmineralogy.org, rruff.info and webmineral.com. The U-V minerals locally occur together. In the mines that have abundant barite, the Ba phase followed by the (Ba,Pb) phase, then the (Ba,Tl) phase of francevillite occur together with metatyuyamunite and tyuyamunite in the deposits and probably reflect the changes in available metals in the fluids sourcing the deposits. Carnotite is found only in trace amounts. Metatyuyamunite and tyuyamunite, the only uranyl vanadate minerals originally identified in the early literature make up approximately 80%-90% of the U-V minerals in the PMD and are the only U-V minerals identified in the LMD. All of these minerals are indistinguishable in hand sample. The Tl-bearing uranyl vanadates have similar morphologies to the francevillites observed in this study. The morphologies of some of the metal-uranyl vanadates observed in this study are shown in Figure 6.

Gangue Minerals Associated with Hydrothermal Dolomitization
Tectonic hydrothermal brecciation with associated hydrothermal dolomite [46] and calcite has been documented in the Madison Group in the Bighorn Basin and nearby Wind River Basin in oil reservoir characterization studies [47][48][49][50][51]. The silicic THB, in this district are the first reported in the Bighorn Basin. Gangue minerals in both deposits associated with hydrothermal dolomitization include sulfides, oxides, hydroxides, hydrothermal "saddle" dolomite [52], and herkimer style, double terminated, quartz [53]. Pyrobitumen and bitumen, though not minerals, are also associated with hydrothermal dolomitization have been observed in this study. Fluorite and barite, only observed in the PMD, are also associated with hydrothermal dolomitization [46]. Other gangue minerals in the PMD that may have hydrothermal origins are halloysite and witherite, identified by XRD and EDS in this study [54,55] Round pyrobitumen preserved in late-stage calcite ( Figure 7E) may serve as an analogue for spherical shaped vugs found in cements of both districts. Vugs likely formed after hydrocarbon or volatiles were released or possibly dissolved in other episodes of fluid migration. Fluids introduced either hydrothermally or later from surface waters may have biodegraded the hydrocarbons by processes such as "water washing" when the breached structure was invaded with water. Spherical vugs with sulfides ( Figure 7F), barite, clay coatings ( Figure 7I,J), witherite, fluorite, Fe-V oxides ( Figure 7I-N), or a combination of these minerals are found in samples from the PMD.
In mineralized silicified THB of the PMD, vugs lined with Fe-V oxides often have U-V minerals in the centers surrounded by bleached quartz rims ( Figure 7M). Spherical, cubic, and hexagonal vugs are present in quartz cements of breccia in the PMD ( Figure 7F,I,K,L). Cubic and hexagonal vugs may reflect minerals that were dissolved out, possibly magnetite, apatite, or barite.
In the LMD, tyuyamunite or metatyuyamunite crystals that radiate from round vugs in calcite cements are closely associated with Fe-oxides and or Fe-hydroxides ( Figure 7O). Irregular shaped inclusions (possibly the remaining residuum after bacteria biodegradation of oil) were observed in calcite cement from LMD samples. Some of the inclusions had U-V minerals growing into them ( Figure 7P); these were closely associated with Fe-oxides or Fe-hydroxides.
Fe-V oxides, identified by EDS ( Figures A6-A9) in samples are found in association with bitumen. EDS spectra of the Fe-V oxides show a high percentage of Fe which varies from around 32% to 70% weight percent. The spectra from the PMD samples also show moderately strong V, Si, Al, O, and no S peaks, suggesting that the metals are iron oxides (magnetite, goethite, or hematite) coated with silicate minerals (probably clays).
Zn was detected in Fe-V oxide samples from the LMD ( Figure A9). The Fe-V oxides may represent an isomorphous substitution of trace elements with Fe-oxides and hydroxides rather than a specific mineral phase. Using EDS, Zn, Cu, and Ni were detected in a mineral from the East Pryor mine area ( Figure A10) likely reflecting metals derived from hydrocarbons based on chemical analyses of the bitumen from this study (Table A2).

Cathodoluminescence
CL revealed the episodic nature of calcite and fluorite. It highlighted U-V minerals and Fe-V oxides showing the close relationship of bitumen or hydrocarbons in the samples. The U-V minerals and Fe-V oxides fluoresce bright green due to V; these minerals occur along fractures and interparticle pore space ( Figure 8A). REE are common activators in fluorite [56,57]. A strong blue CL signal is generally attributed to the LREE fraction, with Eu 2+ , Sm 3+ , and Dy 3+ , serving as the main activators [56] while Dy 3+ and Tb 3+ are usually the activators in HREE environments, which impart a lilac color [56]. Hexaoctahedral fluorite from the Old Glory mine, PMD, alternates between blue and lilac colored zones ( Figure 8B). Spheres of Fe-V oxides fluoresce green in the interparticle porosity of calcite and fluorite ( Figure 8C,D). In transmitted, plane polarized light, the spheres are clearly visible in the same view of Figure 8D and in a standard petrographic image in quartz ( Figure 8E). Luminescence in carbonates is generally due to Mn 2+ in the crystal and quenching by Fe 2+ [58]. The calcite cements from this study fluoresce bright red orange alternating with dark zones, which is indicative of episodic Fe-rich hydrothermal sourced fluids [59]. REE-activated luminescence in sedimentary calcite is much lower than hydrothermal influenced carbonates, which often crystalize from REE-enriched solutions [58,59]. A corroded and fractured quartz veinlet from the Marie mine, PMD, reflects episodic fluid interactions. Calcite and bitumen filled in the pore spaces of the quartz veinlet. Inclusions in the bitumen along the edges of the quartz luminesce green, likely indicating V content ( Figure 8G). In transmitted light view of the same image the quartz appears white ( Figure 8H). Sulfides are also common in the interparticle porosity of calcite and dolomite breccia in both mining districts and in fluorite in the PMD. The sulfides appear black in CL images ( Figure 8I).

Fluid Inclusions
Fluid inclusions in several samples from the PMD were analyzed by FLUID INC., Denver, CO, USA. A barite sample from the Swamp Frog mine contained yellow fluorescing primary inclusions; these were too small to study. Hydrocarbon (HC) inclusions are usually yellow and usually imply some sort of low API gravity [60] e.g., relatively heavy liquids-consistent with the HCs generated from Phosphoria source rocks for oil of the Bighorn Basin [61]. Jim Reynolds (pers. comm.) noted that the quartz and fluorite he examined from our study area appeared hydrothermal in origin, typical of what he observes in low-temperature (<200 • C) epithermal systems. The quartz and fluorite samples with H 2 O fluid inclusion assemblages had highly variable liquid to vapor ratios due to necking and were too small to obtain reliable homogenization temperatures. This is typical for low-temperature quartz in hydrothermal environments, because of the short-lived pulses of warm fluids, there is not enough geological time for maturation of the inclusions at low (<200 • C) temperatures [61].

Rare Earth Elements REE Plus Y
Spidergram plots of REE plus Y normalized to NASC (North American Shale Composite) were produced from the data ( Figure 9) using the values of Gromet et al. [62]. Plots were split up by similar REE plus Y patterns. The data from the bitumen, bentonite, and the average upper continental crust [63] are plotted together ( Figure 9A). This plot shows the bitumen have a low concentration of REE compared to the average continental crust and the bentonite.
Most of the patterns plotted from the data of mineralized and gangue minerals resemble portions of the bitumen patterns. The bentonite sample data plotted with the bitumen and average crustal abundance data has the largest negative Eu anomaly of the data. The Permian Goose Egg hosted bitumen differs from the Mississippian Madison (Mm) Group hosted bitumen. It has lower REE concentration, no Eu, no Lu, and a flatter, irregular HREE and Y pattern. Many of the patterns from the samples in this study reflect portions of the Mm hosted bitumen patterns ( Figure 9B-E). These patterns were split up based on the steepness of the LREE patterns and similar patterns to each other. Some samples from both districts have patterns similar to both the Permian Goose Egg hosted and Mm hosted bitumen ( Figure 9F,G). Other samples have patterns more similar to the Permian Goose Egg hosted bitumen than to the Mm hosted bitumen ( Figure 9H,I). The Tensleep Sandstone (TS53) hosted quartz ( Figure 9H) differs the most of the quartz samples. It lacks Eu and has a LREE composition very similar to the Permian Goose Egg carbonate hosted bitumen, but differs from the bitumen and other quartz samples ( Figure 9F,H,I) with a positive sloping HREE pattern.
All the samples have greater REE plus Y concentrations than the bitumen except for two speleothem calcite samples from Big Pryor Mountain ( Figure 9H). The speleothem calcite have the least concentration of REE of all the samples reflecting the typical low REE concentrations in natural waters [63]. The fluorites are the most enriched in REE plus Y of the gangue minerals and are more enriched than many of the U-V mineralized samples ( Figure 9C,D). The Sandra Mine U-V samples are the most enriched in REE plus Y ( Figure 9C).
The REE and Y pattern for barite from the Swamp Frog mine differs most from other data patterns in this study though is similar to the HREE portion of the Goose Egg bitumen pattern between Dy and Tm ( Figure 9I). The pattern is similar to the quartz with pyrite from the Lisbon mine sample pattern ( Figure 9I). The Old Glory mine sample (solid black line and yellow filled symbols) is similar to the barite sample from Ho-Yb in the HREE portion of the pattern. Both the Old Glory sample and the barite sample are somewhat enriched in Sm relative to the other REE in each sample. The Old Glory sample has Eu where the rest of the samples in this group do not.
Three U-V breccia data patterns from the Dandy, Old Glory and Marie mines do not match the bitumen patterns. These patterns have convex light and middle REE with flat to irregular HREE patterns ( Figure 9J). The patterns in this set as well as those from Figure 9I and the Tensleep Sandstone hosted quartz ( Figure 9H) may reflect a contribution from a different source, possibly the inferred intrusion to the west of the PMD ( Figure 1) [5,6,16], or some other fractionation process during mineral formation.  [20][21][22][64][65][66][67][68][69][70][71][72][73][74][75][76][77]. The deposit has been a major phosphate resource for agricultural fertilizer [70]. Early studies were to evaluate the economic resources of this deposit which is enriched in many metals, define its origin, and aid in exploration for similar deposits [70]. REE and trace element data were used in many of the studies to define the origin [69][70][71]73,74]. Later studies focused on reaction pathways, transport, and fate of potentially toxic trace elements (Se, Cd, Tl, Hg, As, etc.) associated with the occurrence, development, and societal use of phosphate [71,72]. A recent study has shown that phosphorite deposits including the Phosphoria Formation are enriched in REE that would be easily recoverable as a byproduct of phosphate production [75]. The REE occur in the francolite (Ca 5 (PO 4 ) 3 F) portion of the Phosphoria phosphorite rock, a carbonate-rich fluorapatite, which is beneficiated or physically concentrated and dissolved in sulfuric acid to produce fertilizer [76]. In the study, dilute H 2 SO 4 and HCl, extracted nearly 100% of the total REE content of the phosphorite samples in leaching experiments [75].  [66,67,69,70,74,75,77]. Red filled squares are sites that had samples with greater than 100 ppm concentration of Tl. The map is adapted from Sheldon, 1963 [65] and locations from the different studies.
For this study, data from several studies carried out by the USGS and the Wyoming State Geological Survey of the Phosphoria Formation in ID, WY and MT [64][65][66][67][73][74][75][76][77] was examined. Figure 10 shows the locations of samples sites and the approximate location of this study; map shows approximate location of the Permian age rock. Also noted are the locations of samples that had greater than 100 ppm Tl. Not all of the studies included Tl as one of the analytes.
Data from analyses of rock samples studied by the USGS from MT, ID, and WY was normalized to NASC (North American Shale Composite) using the values of Gromet et al. [62]. Two main patterns are present in the USGS data: (1) a strong negative Ce anomaly with HREE patterns showing slightly elevated Gd and Y ( Figure 11A) and (2) patterns with little or no Ce anomaly and flatter HREE and Y patterns compared to the average crustal abundance pattern ( Figure 11B). Most of the data reviewed (approximately 90%) in the USGS studies had the strong negative Ce anomaly. Six sample analyses from this study plotted with data from the USGS studies show similar REE and Y concentrations and patterns ( Figure 11C). Samples from the Sandra mine are comparable in concentration of HREE and Y to the most enriched of the USGS samples ( Figure 11C).
Fractionation of Ce in the marine environment has been attributed to the oxidation of Ce 3+ to Ce 4+ and subsequent removal by suspended Fe and Mn particles in the water column or solid phases residing on the sea floor [75,[78][79][80]. This results in REE patterns with the distinctive negative Ce anomaly characteristic of seawater [75,[78][79][80]. The negative Ce anomaly is often associated with the precipitation of Mn nodules [80]. For an example, Mn nodules occur in a fine grained argillaceous sandstone bed of the Franson Tongue Member of the Park City Formation, Dillon, MT [66]. The depletion of Ce and the REE content of the U.S. Geological Survey samples, including the Permian nodules, were attributed mainly to the apatite fraction present in the samples as the mineral francolite [69][70][71]75]. The bitumen samples in our study lack the Ce anomaly, although they display similar HREE and Y patterns to the Phosphoria rock sample REE patterns. The bitumen REE patterns likely represent the fractionation of the REE as a group due to the maturation and migration of the hydrocarbons and brines from the Phosphoria Formation as the depth of the overlying sediments increased in the Bighorn Basin. REE are mobilized in hot reducing conditions typical of deep basin brines [63].  [20,21,[64][65][66][67][68][69][70][71][72][73][74][75]. Analyses of the Mn nodule used in Figure 11 above from a U.S. Geological Survey study near Dillon, MT had low total REE but high compositions of other metals: 400 ppm As, 1200 ppm Cd, 2200 ppm Co, 70 ppm Mo, 13,000 ppm Ni, 1 ppm Te, 1000 ppm Tl, 700 ppm V, and 25,000 ppm Zn [66].

Concentration of Metals in the Phosphoria Formation and in the U-V Deposits
U-V samples from this study contain anomalously high concentrations of many of the same metals supporting derivation from this source ( Table 2). The bitumen do not have high concentrations but do have many of the elements of concern in this study. All these elements would not stoichiometrically fit into the structure of the minerals in this study. As hydrocarbons and brines migrated into the structures hosting the deposits, it is likely many of these metals were incorporated as inclusions in minerals or breccia cement. The most enriched PMD sample has similar, slightly higher concentration of REE and Y from Sm-Ho to the most enriched USGS sample. The only PMD sample with a negative Eu anomaly is similar to a phosphatic carbonate from the USGS Lakeridge core, WY that also has a negative Eu anomaly.

Phosphoria Oil Migration and Oilfields in the Bighorn Basin
Originally, Phosphoria Formation-sourced hydrocarbons accumulated in the Bighorn Basin in stratigraphic traps created primarily by up dip facies change, pinch out and truncation of reservoir carbonates, and by uneven Phosphoria Goose Egg truncation of the underlying Pennsylvanian Tensleep Sandstone [23,24]. The timing of the first migration of Phosphoria sourced oil ranges from Jurassic to late cretaceous in time [81]. The Cottonwood Creek oil field in the southeastern Bighorn Basin produces from a stratigraphic trap (Figure 12). The oil produced in the field is from a Phosphoria Formation reservoir with a seal formed by impermeable Permian Goose Egg shales and anhydrites [82], similar to the original pre-Laramide stratigraphic traps. Seven of the top ten producing oilfields in WY are in the Bighorn Basin. These are highlighted in red in Figure 12. The primary oil production in the northern Bighorn Basin was from Mississippian age reservoir rock in two of the top seven fields show in Figure 12. These fields, Elk Basin and Garland, may serve as analogues for the structures hosting the U-V deposits in the two districts. Both Elk Basin and Garland produce heavy oil and significant hydrogen sulfide is associated with the Elk Basin oilfield [40].

Origin of Bitumen and Heavy Oil
Most heavy oil and natural bitumen is thought to be the residue of formerly light oil that has lost its light-molecular-weight components through degradation by bacteria, water-washing, and evaporation [83,84]. Meteoric water is the commonly accepted vehicle for bringing dissolved oxygen and microbes into contact with oil reservoirs, either along faults or fractures or by aquifers in hydrodynamic connection with the surface [84]. Generally, this is under aerobic conditions at depths of about 1524 m or less and temperatures below 80 • C [83]. Hydrocarbons are oxidized completely to carbon dioxide and water without intermediate products by many bacteria [84]. This results in the remaining residue of solid bitumen and asphalts which are enriched in heavy N, S, O compounds [84]. Some average chemical variations in conventional oil and bitumen that have been observed in studies from basins around the world show an increase in the metal concentration in bitumen compared to conventional oil. Based on studies from approximately 8500 deposits there is an increase in the metal concentration in bitumen compared to conventional oil [83]. For example conventional oils have average concentrations of approximately 1 ppm Pb, 16 ppm V, 6 ppm Fe, 1 ppm Al, 8 ppm Ni, and 0.5 ppm Cu. and bitumen has concentrations of approximately 5 ppm Pb, 335 ppm V, 4290 ppm Fe, 21,000 ppm Al, 89 ppm Ni and 44 ppm Cu [83]. Hg however, is less concentrated in bitumen 0.019 ppm compared to an average of about 19 ppm in conventional oil [83]. This may suggest that Hg is being mobilized in the biodegradation process. The combination of aerobic and anaerobic heterotrophic bacteria metabolic processes and subsequent water washing likely results in mobilization of some metals during biodegradation and enrichment or a sink for others.

Style of Breaching of Oil Traps in the Study Areas
We suggest that the structures hosting the U-V deposits likely had Phosphoria sourced oil in structural traps until they were breached. Breaching of the structures in each district was not synchronous or of the same style. The folds hosting the deposits are of different scales, maturity, and elevation. We suggest that breaching and subsequent U mineralization in the PMD was earlier than in the LMD based on elevation and structural maturity of the folds.
The PMD is approximately 1000 m higher elevation than the LMD. The Gypsum Creek anticline, the main fold hosting the U-V deposits in the PMD, was breached as the Crooked Creek fault propagated up through its core extending beyond the Madison Group.
The folds hosting the deposits in the Little Mountain area are smaller in scale and less structurally mature than those in the PMD in the sense that the master faults did not fully propagate through the structures. The Porcupine Creek anticline, one of the main structures hosting the U-V deposits in the LMD, was breached by the superposition of the ancestral Bighorn River across its axis, forming part of the Bighorn Canyon.

Paragenesis of U Minerals in the Two Districts
The paragenesis of U minerals in the two districts is slightly different due to structure, elevation and position in the Bighorn Basin. The initial evolution of the U-V mineralization in both districts is related to the migration of the oil and brines from the Bighorn Basin into structural traps as folding proceeded in the two areas. Simplified conditions for mobilizing metals in a basin environment are increased temperatures, lower pHs, and increased salinities. U would be stable in a reduced U 4+ state in the basin environment and would likely form F and Cl complexes, which stabilize U 4+ in solution [85]. Abundant F would have been available from the apatite of the Permian Phosphoria Formation. These complexes could have migrated with oil and been in solution in the traps until the structure was breached. Once oxidized meteoric water was available U +4 would oxidize to U +6 to stabilize the uranyl ion, UO 2 2+ , and its complexes [86,87]. The uranyl ions would bind with available metals and precipitate out as various uranyl minerals. All the identified U minerals in the districts are uranyl minerals with the Ca-uranyl vanadates being the most common. The uranyl silicates, phosphates, and vanadates are relatively insoluble and require Si, P, and V 5+ , respectively to precipitate [86]. These elements can be derived from a variety of sources. In many U deposits, metal-uranyl vanadates are derived from the oxidation of U +4 minerals such as uraninite [85][86][87]. In this study, however, no U +4 minerals have been observed. This suggests that the metal-uranyl vanadates are primary (did not oxidize from a reduced U mineral such as uraninite). We suggest that these metals were derived from the biodegradation of the oil once the basin fluids mixed with the meteoric groundwater after the structures were breached.

Mineral Paragenesis PMD
The two mining districts are hosted in the same paleokarst horizon with nearly identical ore minerals but differ in gangue minerals; notably fluorite, barite and abundant quartz present only in the PMD. Two episodes of siliceous fluid migration are apparent in the PMD. The first significant siliceous fluid migration was the emplacement of silica cemented THB. We suggest these were emplaced as a result of earthquake induced seismic pumping during the Laramide Orogeny that hydrofractured the host rock [47,88,89]. The resulting fractures and brecciation were healed by the rapid precipitation of quartz in the earlier episodes during subsequent and sudden decrease in the partial pressure of CO 2 [88]. As the master faults propagated along the same faults and the process would be repeated forming new or following previous fractures [88]. The earlier siliceous THB had abundant Fe-oxides with trace amounts of V, Tl, Pb, and Ti present in the matrix ( Figure 13A). Early barite was associated with the first pulses of siliceous THB evident by the inclusions of barite in the siliceous cement of the THB ( Figure 13B). Later barite precipitated along fractures and in vugs ( Figure 13C). Quartz, barite and fluorite also precipitated after the emplacement of the siliceous THB in fractures and vugs, often accompanied by U-V minerals. These minerals likely represent the breaching of the structure by the full propagation of the fault through the Madison Group. The oxidized meteoric ground waters likely mixed with the oils and brines that were present in the structure, migrated down into fractures, and precipitated the ore and gangue minerals in the vugs and fractures of the THB and along other fractures in collapse paleokarst in the area ( Figure 13E-G). Oxidized fluids migrating down into fractures occupied by the THB could explain the hematite-and limonite-staining that resulted from the oxidation of the Fe-V oxides present in the THB. The meteoric sources may have been lakes, ponds and or an ancestral river that occupied the area during the time the structure was breached or later as the Bighorn Basin was exhumed.
The quartz cement in the THB may be related to the inferred intrusion. It may have formed in a similar timeframe to the Sliderock volcano and the Lodgepole intrusive, part of an extensive NW-trending volcanic field that lies astride the folds and faults of the Nye-Bowler lineament [90,91] (Figure 3, laccoliths). The alignment of these volcanic intrusives along the Nye-Bowler lineament strongly suggests the fault zone controlled the emplacement of the volcanic centers of the field [90]. Thus, it would not be unreasonable to have another intrusion along the same NW-trending lineament that formed because of low-angle subduction of the Farallon plate. The igneous products of the Sliderock Mountain stratovolcano were interpreted to be the products of arc-related magmatism based on geochemistry [90]. The inferred intrusion may account for the Si, Ba, Au, and other metals present in the minerals of the PMD and not present in the LMD. Evidence supporting the inferred intrusion near the PMD include: High aeromagnetic anomaly identified in 1982 southwest of the Big Pryor Mountain block by the Department of Energy [16] and supplemental geophysical work by the USGS in 1985 [5,6]. Gold detected in samples from this study by commercial assays (Appendix A, Table A1) and using the SEM. An abandoned gold lode deposit west of the PMD may have been associated with an intrusion (Figure 1, green symbol). High concentrations of Tl (up to 490 ppm) in the PMD. Tl is not common but is commonly associated with Carlin-type gold deposits that have a close proximity to Eocene volcanic centers or plutons [92][93][94][95][96]. Thus, the Tl and Au detected in the PMD could support an intrusion. Proximity of the Late Cretaceous 78-75 Ma. diorite to granodiorite Sliderock stratovolcano and the Lodgepole intrusive [90,91] that is suggested to have been controlled by the E-W trending Nye Bowler lineament. Abundant fluorite and barite in the PMD. High concentrations of Ba in the PMD samples with the highest values excluding the barite samples ranging from (5460-6550 ppm). High silica content of the tectonic breccias, highly fractured zones with white "bull" quartz and a few outcrops or lag deposits of the same, and high silica content (95%) in reported ore during the years that the mines were operating [5].

Mineral Paragenesis LMD
The LMD deposits are hosted in smaller folds than the Gypsum Creek and East Pryor Mountain anticlines. The close association of bitumen and Fe-oxides and hydroxides with the U-V minerals in the LMD suggests that oilfield basin brines were present in the structures until the down cutting of the Bighorn River breached the Porcupine Creek anticline. Gypsum, dogtooth calcite, bitumen, sulfides and speleothem type occurrences from waste piles and outcrops suggest that hypogene sulfuric acid speleogenesis [18,19] may have contributed to concentrating ore forming minerals in the LMD. Though underground mines are no longer safely accessible descriptions of the caves encountered and actually mined out in the LMD [6,7] support this interpretation. REE data show good correlation with the Mm hosted bitumen ( Figure 9E). Hypogenic caves likely developed after oxygenated water and microbes were introduced to the system. Hydrogen sulfide that was produced by microbes consuming the light fraction of the hydrocarbons likely contributed sulfuric acid speleogenesis initiation in the LMD and to the mobilization of the metals needed to form the U-V minerals. The bitumen that was the residuum of the microbial process as well as the caves support this interpretation. The presence of Fe-and S-reducing bacteria in the Lower Kane Cave system [18] as well as U and V present in the mud and water of the geothermal artesian spring may serve as an analogue for the LMD. The Lower Kane Cave is actively forming by sulfuric acid speleogenesis. It has an artesian geothermal spring emanating from fractures at the core of the anticline that are mixing basin derived fluids with the groundwater at the level of the Bighorn River.
Timing of U mineralization and active sulfuric acid speleogenesis in the LMD likely coincided with the timing of formation of Spence Cave in Sheep Mountain, WY, upper Kane Cave in Little Sheep Mountain, WY and the Horsethief cave formed near the mining district. Sheep Mountain and Little Sheep Mountain are anticlines that have been superposed by the Bighorn River. Each of these caves is proposed to have developed by hypogene sulfuric acid speleogenesis [92]. Windblown sand and volcanic ash that were deposited in the Spence and Horsethief caves by wind were dated using cosmogenic nuclide concentrations [92]. A cosmogenic 26 Al/ 10 Be burial date for eolian sand in the entrance to Spence Cave yielded a burial date of 0.31 ± 0.19 million years. Tephrochronology of ash in the Horsethief Cave is the Lava Creek B fallout ash, erupted from the Yellowstone Plateau volcanic field ca. 0.64 million years [96]. These dates may provide minimum estimates for the timing of cave development in the northern Bighorn Basin and for the U-V mineralization in the LMD.
The LMD lacks the strong silica component, barite, and fluorite that the PMD has. If the oilfield brines that originally migrated into structures had F from the abundant apatite it would seem likely that there should have been some present also in the LMD. This may reflect more F and Ba available from the inferred intrusion in the PMD and not as much concentration in either Ba or F to form barite and fluorite in the LMD. The lack of siliceous THB in the LMD may suggest the district was structurally separated from the fluids that formed the early THB in the PMD. It may also be due to less sampling in the LMD than in the PMD due to accessibility, though no barite or fluorite was observed in this study or documented in previous studies. Ba was detected in high (4300 ppm) concentration (Appendix A, No. 19) in an unmineralized breccia from the LMD though barite was not observed in the field, petrographically or using the SEM.
Although the paragenesis of each deposit examined varies, the geologic events and processes affecting Madison Limestone in the two districts have some consistent similarities. The origin of the deposits reflects a close but complex relationship with the migration of hydrocarbons, hydrology, and the exhumation of the Bighorn Basin.

Future Directions
After elevated U was identified in home water wells on the Crow Reservation collaborative research related to U contaminated drinking water on the Navajo and Crow Indian Reservations was established with the University of New Mexico, the University of New Mexico and Montana State University. All three Universities are working on projects related to U and the health of people on reservations. Part of the research being conducted at the University of Montana, is the development of silica polyamine composites to aid in removal of U and other toxic metals from contaminated drinking water.

Selective Solid Phase Extraction of U from Contaminated Water
In an effort to address the long-standing health hazards imposed by abandoned uranium mines and mills, naturally occurring U, and trace metal contamination research is currently being conducted at the University of Montana using Silica Polyamine Composites (SPCs). SPCs filter, isolate, and remove unwanted metals. A variety of other methods such as precipitation, solvent extraction, electrolysis and ion exchange have been employed to remove dissolved metals from aqueous samples. Most of these have disadvantages such as poor removal efficiency, high cost, and generation of secondary pollutants and ineffectiveness for low concentration removal [93]. We are testing contaminated water from wells on the Navajo Reservation and using these SPCs as a potential method of remediation through various bench scale studies. Our aim is to eventually produce point of source filters that are economical, efficient and available to people on American Indian reservations and other areas in need of ways to clean contaminated drinking water.
Since the mid-20th century amorphous silica gels have been a favored solid-state matrix for applications in chromatography, catalysis, colloid chemistry and as drying agents [94]. SPCs are solid phase hybrid materials that consist of an inorganic matrix combined with an organic polymer developed by the Rosenberg Research Group at the University of Montana [95,96]. These hybrid materials bring together the best properties of each component in order to enhance functionality. These silica based organic-inorganic materials offer a rigid matrix with high porosity and good thermal stability [97]. Furthermore, the polar nature of the silica polyamine surface also makes for better mass-transfer kinetics in the case of aqueous solutions, and the polyamine can be easily modified with metal selective functional groups [98,99]. SPCs act as chelating agents for a range of metal ions and can also be useful for toxic metal immobilization and disposal [99]. They have led to a wide range of metal-selective materials that in most cases provide a high loading, rigid, surface that does not shrink or swell [97]. SPCs are currently being used in medium to large-scale recovery of metals and the remediation of toxic anionic and cationic pollutants [96,97].
Acid stripping is the primary method of extracting immobilized metals from SPCs. Without the acid strip, metals can remain on the SPC and can be used as a medium for long-term disposal [99]. These composites have been previously found to maintain activity over more than 7000 test cycles that consist of treating a packed column of SPCs with a metal ion solution and then acid stripping. These composites showed no visible signs of degradation and negligible loss in metal ion adsorption capacity under a variety of conditions [95].
SPCs are materials that are designed to bond with a specified metal and extract the metal from aqueous solutions for later recovery [95][96][97]100,101]. SPCs have the ability to be modified with a variety of metal selective ligands which is vital to an effective remediation processes. Previously the Rosenberg research group has synthesized a number of different SPCs that are selective for Cu 2+ , divalent transition metals, Fe 3+ , Ni 2+ , Ga 3+ , divalent and trivalent metals ( Figure 14). In addition, earlier studies have evaluated the ability of SPCs to separate REE from mine leachates from Mt. Weld in Laverton, Australia. These studies separated Ce(III), La(III), Nd(III), Sm(III), and Pr(III) from the solution matrix that contained Fe(III), Mn(II), Ca(II), Mg(II), Al(III), Zn(II), Ti(IV), and other lower concentration trace metals [102]. This separation required first the use of WP-4 followed by BPAP. The WP-4 removed the large excess of iron and BPAP selectively removed the entire family of lanthanides from mixtures of transition metals, alkali and alkaline earth metals at low pH [102]. This study demonstrates the feasibility to separate REE metals from a metal matrix like that of mine leachates. Thus, if the concentration of REE or other desirable metals was economical they might be recovered using SPCs possibly offsetting the cost of remediation efforts.

Elements of Environmental Concern
Alluvial low lands of the Crow Reservation, MT, located along the Bighorn and Little Bighorn Rivers are where U was found to exceed EPA's MCL in some tested home wells [3]. The source of the U has not been determined but it is likely that material eroded from the Wolf, Big Horn and Pryor Mountain ranges may have had U minerals present either from the U-V deposits, or thin Permian Phosphoria sediments. If the U was derived from erosion of materials associated with the U-V deposits or from Phosphoria Formation sediments then there is a possibility other metals of environmental concern e.g., As, Cd, Hg, Pb, Se, and Tl may be present. Tl is not found in many rock types and Tl minerals are rare. The average crustal abundance of Tl is 0.75 ppm [103]. Tl is one of the metals found in anomalous high concentrations in both the Phosphoria Formation and in the U-V deposits of this study. All the mineralized samples from this study contained anomalous concentrations (10-490 ppm) Tl. Tl is more toxic to mammals than Cd, Pb, Cu, or Zn [104,105]. Tl is associated with hydrothermal sourced fluids in sediment hosted gold deposits [106][107][108][109]. Tl in these deposits was likely mobilized and transported with hydrocarbons and brines though could be related to the inferred intrusion that might also account for the Au and Ag present in the U-V deposits. Tl disperses readily and accumulates in soil, water transport and food crops such as cabbage, affecting the health of people [110]. Tl forms highly volatile complexes with halides and can enter the crystal lattices of sulfides and sulfosalts [111]. Drinking water wells on the Crow Reservation, MT, were analyzed for Tl because U was found some of the wells tested and might be associated with the eroded material sourced from deposits in WY and MT. Initial well water samples (n = 29) from the Crow Reservation tested for Tl at the University of New Mexico did not show elevated concentrations of Tl. The results ranged from 0.01 to 0.36 µg/L well below the maximum contaminant level (MCL) of 0.002 mg/L and the EPA's MCLG (goal) of 0.0005 mg/L [112].

Conclusions and Future Work
U-V mineralization is complex in the two districts. Our proposed model for paragenesis is a simplified sequence of mineralization. Both districts are hosted in the same paleokarst horizon but are in slightly different structural positions in the Bighorn Basin. The structures in both mining districts likely held hydrocarbons in structural traps until they were breached. This conclusion is based on observations of bitumen in the samples both in the field and in the lab using standard petrography, SEM, and CL techniques.
We conclude that the source of the many of the metals in the deposits of this study was likely Permian Phosphoria Formation oil introduced through fluid migration in fractures as folding proceeding during the later stages of the Laramide Orogeny. The process of episodic tectonic hydrothermal brecciation accompanied folding and provided both a mechanism and fluid source for metals in both districts. This conclusion is based on similar REE and other trace element geochemistry of: (1) the bitumen from Little Sheep Mountain; (2) the U-V breccia and gangue minerals from this study and (3) the data from Phosphoria Formation studies conducted by the USGS.
The exhumation of the Bighorn Basin affected each area differently and mineralization style and gangue minerals reflect this difference. U was a later element in both districts where it followed the Fe-V oxides. Fe-V oxides were emplaced with the earliest siliceous THB in the PMD. The episodic nature of hydrothermal tectonic brecciation and the interaction of meteoric groundwater has created a complex mineralization sequence. The model has components of both hydrothermal and meteoric sources. It is likely that most of the V was hydrothermally introduced while the U was probably derived from both hydrothermal and meteoric sources. Especially as U 6+ is very mobile in water and is still being mobilized evident by U-V minerals found on the surfaces of burned wooden supports of worked prospects in the PMD, on surfaces of reclaimed materials, along sprung stylolites and along open fractures. The triuranium carbonate minerals are generally the most soluble [85] and some are still present in the deposits. As surface water percolates through the deposits these minerals are likely dissolved and reprecipitated as the relatively insoluble metal-uranyl vanadates when encountering a V source such as material from sprung stylolites. Thus, surface waters flowing through these deposits probably do not contribute many of these elements to the hydrologic system. This study does not support a direct source of U, Hg or Pb being derived from the U-V deposits as a source of contamination in the Bighorn River or home wells on the reservation. Though material derived from the erosion of the deposits may be a contributing factor.
Hg may be entering the hydrologic system through the geothermal Kane Cave system that flows directly into the Bighorn River above the Bighorn Reservoir area. Drops of oil were described in S.J. Egemeier's 1973 dissertation which helped to identify sulfuric acid speleogenesis as a cave forming process and suggested that the ore deposits might have formed through this process. Sulfate reducing bacteria present in the Kane cave system [19] may be processing the sulfate (SO 4 2− ) and also taking up mercury in its inorganic form and converting it to methylmercury through metabolic processes. The conversion of inorganic mercury to methylmercury increases the toxicity. These methylmercury-containing bacteria may be consumed by the next higher level in the food chain, or the bacteria may excrete the methylmercury to the water where it can quickly adsorb to plankton, which are also consumed by the next level in the food chain. Thus, future work is warranted to investigate the microbial processes related to biodegradation of hydrocarbons which may show that Hg is indeed being introduced to the hydrologic system through these processes. The geochemical data from Table A2 shows that the mineralized U-V samples (No. 1-13) that were assayed had the highest Hg (0.179-12.8 ppm) and the select gangue minerals that were assayed for Hg all had detectable Hg. These included quartz sample No. 35 (0.095 ppm); dolomite breccia sample, No. 34 (0.02 ppm); two fluorites No. 28 (0.618 ppm) and No. 29 (1.05 ppm); and calcite sample No. 26 (0.034 ppm). The low concentration of Hg remaining in the bitumen in this study (<0.01-0.01 ppm samples [39][40][41] suggests that Hg was likely mobilized somewhere between the oil and bitumen phases, assuming there was a higher initial concentration in the oil. This is supported by the study [83] described earlier in the text in the section on heavy oil and bitumen. Hg and the other metals including the REE's do not fit in the structures of most of the minerals that contain them. This suggests that the Hg, REE and other metals were likely incorporated as Permian Phosphoria Formation sourced hydrocarbon or brine fluid inclusions in the minerals. It is likely that Hg and Pb may be entering the hydrologic system from atmospheric anthropogenic sources as well. Future work might include: (1) tracing anthropogenic Hg and Pb input using stable Hg and Pb isotope ratios in sediments; (2) studies on how bacteria metabolize the metals of interest (Hg and Pb) using oils from the Garland and Elk Basin oil fields that are producing from the Mississippian and Pennsylvanian reservoir rocks in the northern Bighorn Basin; and (3) using portable DNA analyzers in some of the hydrothermal springs entering the Bighorn River area to identify different bacteria that may be metabolizing the metals.
Water quality on many reservations will undoubtedly continue to be an issue. Research will be continued to develop a way to remove U and other metals from the water to ensure the most important need on the Navajo, Crow, and other reservations, clean water for good health. Author Contributions: Anita L. Moore-Nall: Researched the area's geology, did the field work, submitted samples, interpreted the geological and geochemical data, produced most of the figures and is primary author on paper; Ranalda L. Tsosie: Researched, did the chemistry and wrote the Selective Solid Phase Extraction of Uranium portion of the paper; reviewed drafts of the paper; contributed to the manuscript.

Conflicts of Interest:
The authors declare no conflict of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

Appendix A
Five samples were analyzed for platinum group metals (PGM) and Au by Fire Assay with ICP finish by American Analytical Services (AAS), Osburn, ID, USA. The samples included two U-V mineralized breccia, two barites, and a fluorite. Table A1 shows these results. Values are in ppm. The Nevada branch of Australian Laboratory Services (ALS) Global analyzed samples for REE and trace elements. A total of 41 samples were analyzed: 18 U-V mineralized breccia, 6 limestone and dolomite host rock, 5 calcite, 1 bentonite, 1 barite, 5 fluorite, and 4 quartz were analyzed for REE by by Lithium Borate Fusion and ICP-MS (Fused bead, acid digestion and ICP-MS) in a 30 element package (ME-MS81). Three samples containing bitumen were analyzed for REE by Aqua Regia and ICP-AES/ICP-MS (ME-MS41r) a 51-element package. Twenty two samples were analyzed for base metals by a four acid digestion and ICP-MS (ME-4ACD81) to test for Pb and 17 samples were analyzed for trace element Hg by Aqua Regia and ICP-MS (Hg-MS42). Additionally 17 mineralized samples were analyzed by by XRF (V,U-XRF10 ) for ore grade % U and 3 for ore grade %V as these samples were above the maximum detection level using ME-MS81. Trace element data is presented in Table A2 and REE data is presented in Table A3. Figure A1. EDS spectrum of Au mineral from Lisbon mine sample. The Y-axis shows the counts in 10 3 (number of X-rays received and processed by the detector) and the X-axis shows the energy level of those counts in kilo-electron volts (keV).   Appendix B Figure A2. Image A in figure shows Old Glory mine silicified breccia sample with U-V minerals in a vug. Image B is a higher magnification of Image A. The EDS spectrum was collected from + in Image B. Carnotite is indicated based on the empirical formula K 2 (UO 2 ) 2 (VO 4 ) 2 ·3(H 2 O) which has the weight percentage for elements K (8.67%) V (11.29%) U (59.86%) [113] outlined in red in the PUzaf results. Sample was coated with Ir. The Ir peak is denoted by the *Ir and arrow. Ir coating was subtracted out before the normalization was performed by the program utilized by the instrument.    [114]. The sample may be a Tl-francevillite with Tl substituting for Pb. Sample was coated with Ir. The Ir peak is denoted by the *Ir and arrow. Ir coating was subtracted out before the normalization was performed by the program utilized by the instrument. Figure A6. FeV oxides with Tl, Pb and Ti from the Old Glory mine, PMD. The large Si peak is likely picked up by the Si matrix. Davidite, a U REE mineral in the Crichtonite Group, a group of about 13 complex metal titanates, identified by early workers using optical methods, was described as occurring in disseminations in and interstitial to a fluorite matrix [5,14]. Fe has the largest weight percent in the oxides of this sample. A much larger percentage of Ti would be needed to have a Crichtonite Group mineral, though these may occur locally in other samples. The FeV oxides show up well in image B, BSE, 1000× magnification. EDS collection site + is shown in image A, SEI 2000× magnification. Petrographic image C shows inclusions of dark oxides in a quartz matrix and vein urefrom the Old Glory mine, PMD, xpl. Sample was coated with Ir. The Ir peak is denoted by the *Ir and arrow. Ir coating was subtracted out before the quantification was performed by the program utilized by the instrument. Figure A7. EDS spectrum collected from + from Fe-V oxides in sample from the Old Glory mine, PMD. The elements V, Fe, Tl, Pb, Al and O were detected. These oxides have a lighter colored rim indicating element(s) with larger Z. The rims of the oxides may have the Tl and or Pb. Sample was coated with Ir. The Ir peak is denoted by the *Ir and arrow. Ir coating was subtracted out before the quantification was performed by the program utilized by the instrument. Figure A8. EDS spectrum collected from + in image above. As was determined by the peak at its characteristic K a energy of 10.53 keV and Dy from the peak at 1.293 keV where it partially overlaps with the As L a line. EDS analysis suggests that native iron is the phase located at the red dot; however, seems unlikely due to the associated minerals, and may be an artifact of thin section preparation. This work was performed at Montana Tech-CAMP by Gary Wyss. Figure A9. EDS spectrum collected from + in the BSE image above. Elements detected were C, Ca, Fe, V, Ti, Zn, As, Tl, Si, Cl and O. The sample was coated with Ir. Cl is likely from the epoxy for mounting the sample. Sample is a polished thin-section from a limestone breccia, Leo Incline area, LMD. The U-V minerals occur in close association with Fe-oxides and or hydroxides and bitumen. The black area in the image is an inclusion of bitumen. The red orange Fe-oxides, some with darker Fe-V rims show up well in plain polarized light as do the yellow to yellow green U-V minerals. Figure A10. EDS spectrum collected from the red plus symbol + in the BSE image for unidentified mineral from the East Pryor mine, PMD, that did not have U. The elements Zn, Ni, Cu, and O were detected. This mineral likely formed from the metals derived from hydrothermal fluids associated with Permian Phosphoria Formation sourced oil.