The Bigrlyi Tabular Sandstone-Hosted Uranium–Vanadium Deposit, Ngalia Basin, Central Australia

: The Bigrlyi deposit is a tabular, sandstone-hosted, uranium–vanadium deposit of Carboniferous age located in the Ngalia Basin of central Australia. The deposit is similar to the continental, ﬂuvial Saltwash-type of sandstone-hosted U-V deposits which are well known from the Colorado Plateau, USA. Most mineralization at Bigrlyi occurs as thin, multiple-stacked, stratiform lenses at the base of ﬂuvial channels near the contact between a grey sandstone succession and a hematitic, purple–red sandstone succession. grading into vanadian illite, the interlayer mineral corrensite, and altered detrital biotite. The V-Fe–oxyhydroxide minerals montroseite, haggite and doloresite, and altered detrital Fe-Ti oxides, are minor V-hosts. Mineralized zones correlate with enrichments in Se, Li, Ba, Be, Mo, Mg and Fe, and elevated Se / S ratios are characteristic of U-mineralized zones. Petrographic studies show that a heterogeneous mixture of variably mineralized lithic clasts is present; in the same rock, some clasts are Fe-rich and only weakly U-V-mineralized, while other clasts are strongly V- and / or U-mineralized. These observations point to mineralization processes that did not take place in-situ in the host sandstone at the site of deposition as required by conventional groundwater models. Lead isotope results provide evidence of the open-system mobility of radiogenic elements in parts of the deposit. In V-bearing zones, radiogenic Pb contents were found to be unsupported by current U levels, suggesting that over time U has been mobilized from these zones and redistributed, resulting in U-enrichment in other parts of the deposit. Mobility pathways were likely open over time from early in the history of the Bigrlyi deposit. A hybrid mineralization model, involving an interplay between solution-precipitation processes, detrital transport and post-depositional U remobilization, is proposed for Bigrlyi. Ferrous-ion-bearing clay minerals and pyrite are considered to be the most likely primary reductants / adsorbents, while the deposit is lacking carbonaceous matter. using an extensive dataset of drillhole-derived data that includes downhole geochemical assays and gamma-ray logs. Calcium has been used as a proxy for calcite to better understand the distribution of cemented sandstone. Drillhole assay data for V 2 O 5 , U 3 O 8 and Ca were compiled and composited to 1 meter intervals downhole, producing some 25,300 records. A plot of U 3 O 8 versus V 2 O 5 concentrations on this basis shows that the Bigrlyi deposit is dominated by vanadium-rich but uranium-poor intervals, with 90% of the assayed meter intervals being of this type, but overall U and V are well correlated.


Introduction
The Bigrlyi deposit is a tabular, sandstone-hosted, uranium-vanadium deposit of Carboniferous age located in the Ngalia Basin of central Australia [1]. Deposits of this style are distinct from

Bigrlyi Deposit
The Bigrlyi U-V deposit occurs over 11 km of strike length within steeply upturned beds of the MES; it was discovered by a ground scintillometer survey in 1973 as part of a targeted uranium exploration program in the Ngalia Basin [25]. Most U-V mineralization at Bigrlyi occurs as thin, multistacked, stratiform lenses in close proximity to the contact between a pyrite-bearing, grey sandstone unit and an underlying purple-to-red, variably oxidized, hematite-bearing sandstone. Additional mineralization occurs higher in the stratigraphy, near the upper contact between grey sandstone and overlying red-beds. Mineralized lenses are up to a few meters wide, and are typically ~50-100 m in The Ngalia Basin was affected by the long-lived intraplate Alice Springs Orogeny (ASO), which involved several tectonic movements over a period of 130 Ma from the Late Ordovician to Late Carboniferous [16]. Tectonic events in the late Neoproterozoic and early ASO led to the isolation of the Ngalia Basin from other elements of the original Centralian Superbasin [19][20][21]. Late in the ASO, the basin was affected by a phase of N-S shortening that led to basin inversion, associated thrusting and folding, and uplift of the adjacent Proterozoic Arunta basement [16,19,22]. The Carboniferous Mount Eclipse Sandstone (MES) was deposited with angular unconformity on older basin units as a consequence of this event [16]. A detrital zircon U-Pb age spectrum from the MES [23] indicates that the major source of detritus was likely~1770 Ma age granites, including the proximal Carrington Suite and Yalgoolgarrie Granite, which outcrop to the north and east of Bigrlyi. Such sources are consistent with measured paleocurrent directions [18]. The MES is typically folded, forming broad synclinal and anticlinal structures, and is often tilted vertically or overturned at the basin margins and along internal basin structures; it is the only known host for U-V mineralization in the basin. The present-day maximum thickness of the MES is 2400 m [14,24].

Bigrlyi Deposit
The Bigrlyi U-V deposit occurs over 11 km of strike length within steeply upturned beds of the MES; it was discovered by a ground scintillometer survey in 1973 as part of a targeted uranium exploration program in the Ngalia Basin [25]. Most U-V mineralization at Bigrlyi occurs as thin, multi-stacked, stratiform lenses in close proximity to the contact between a pyrite-bearing, grey sandstone unit and an underlying purple-to-red, variably oxidized, hematite-bearing sandstone. Additional mineralization occurs higher in the stratigraphy, near the upper contact between grey sandstone and overlying red-beds. Mineralized lenses are up to a few meters wide, and are typically~50-100 m in length with a distinct SW plunge. The Bigrlyi deposit has estimated mineral resources of approximately 10 kt U 3 O 8 at a grade of 0.13%, and 20 kt V 2 O 5 at a grade of 0.11%, for a 0.05% cut-off level for both metal oxides [26].
The Bigrlyi deposit is divided into four major sub-deposits (from west to east): Anomaly-2 (A2), Anomaly-4 (A4), Anomaly-7/9 (A7/9) and Anomaly-13/15 (A15) (Figure 1). The A4 and A15 sub-deposits host most of the Bigrlyi uranium and both are predominantly reduced deposits in which the bounding impermeable units are carbonate-cemented beds and shales. A7/9 is a smaller vanadium-rich sub-deposit, and A2 is oxidized to a >100 m depth and is carnotite-bearing. The sub-deposits are cut by a number of oblique faults with displacements of up to about 40 m, which have disrupted, and in places, structurally thickened the mineralized zones. However, there is no evidence for structural control of mineralization.
The grey sandstone unit that hosts the majority of mineralization is a fluvial channel sequence dominated by medium-grained, immature arkosic sandstone with common interbeds of siltstone and shale and variable contents of calcite cement. Lithic clasts of mudstone and siltstone are locally common. The underlying variably oxidized purple-to-red sandstone, which is characterized by the presence of patchy, hematite-bearing mottles, is sedimentologically indistinguishable from grey sandstone and only hosts minor mineralization. The combined thickness of both units decreases from~280 m at A2 to the west to <100 m at A15 in the east, and grey sandstone pinches out entirely east of A15. Fining upward cycles are well developed, typically defined by a basal erosional surface, which is bound by cross-stratified channel sandstone containing abundant mudstone 'rip-up' clasts and scattered pebbles of quartzite; trough cross-beds are common. Both, the grey and purple-red units are characterized by~1-10 m thick fining upward cycles with the tops of the cycles comprising fine micaceous sandstone overlain by either green siltstone (in grey sandstone) or red-brown siltstone (in purple-red sandstone), which commonly contain carbonate (calcrete) nodules and lenses at the very top of the cycle. The fine-grained units are interpreted to be abandoned channel tops [27].
The depositional environment is interpreted as a mature alluvial fan system comprising stacked fluvial channel sequences with occasional abandoned channel tops. The variation in color and mineralogy is related to the water saturation level in the channel fill due to the ephemeral nature of deposition. As such, the grey sandstone shows the development of early diagenetic kaolinite relative to illite in a fully water saturated channel, while hematite-bearing channels reflect variable water saturation with ongoing oxidation. Prolonged periods of aridity led to calcite precipitation due to evaporation [27].
Impressions and casts of fossil plant stems, trunks and other debris are locally present, especially in the finer grained channel sandstones [14]. However, unlike the Colorado Plateau U-V deposits, preserved carbonaceous matter is uncommon at Bigrlyi and is rarely found despite some earlier reports to the contrary [25]. Further east in the Ngalia Basin, carbonaceous and coaly shales have been intersected by shallow drillholes into the MES, but these units are not strongly mineralized [14,28].

Materials and Methods
Sampling and portable X-Ray Fluorescence (pXRF) spot analyses of drill core using an Olympus-InnovX Delta instrument were carried out during on-site visits to Energy Metals' Bigrlyi exploration camp and core yard. However, although standards were monitored, no matrix-specific calibration was carried out and values should only be used as indicators of relative abundance.
Polished thin sections were prepared for a range of samples of various mineralization styles and were characterized at CSIRO's micro-analytical facilities in Perth, Australia. A Zeiss AxioImager semi-automated petrographic microscope was used for the routine acquisition of plane-polarized, cross-polarized and reflected light imagery with a resolution of 2.5 µm. Selected sections were scanned with a Bruker Tornado M4 X-Ray Fluorescence (XRF) element mapper using a spot resolution of 25 µm.
To provide detailed petrographic, mineralogical and element distribution information, polished sections of drill core were examined using a Phillips XL40 Scanning Electron Microscope (SEM) with an attached Bruker energy dispersive X-ray spectroscopic system (EDS) and cathodoluminescence (CL) detector, and a Zeiss Field Emission Gun Scanning Electron Microscope (FEG-SEM) equipped with a Bruker EDS system. Backscattered electron (BSE) images, EDS element maps, CL images and composites of these were used to characterize the samples.
Multi-element geochemical data for Bigrlyi samples were extracted from Energy Metals' extensive in-house validated exploration database. Whole-rock major and trace element geochemical analyses (101 samples) were undertaken by commercial laboratories ALS Global and Labwest Minerals Analysis using low detection limit four-acid digest and lithium-borate-fusion ICP-MS methods monitored by certified standards. Organic carbon (240 samples) was determined by acid leach/combustion furnace analysis.
Lead isotopic compositions (63 samples) were determined by commercial laboratories ALS Global and Labwest Minerals Analysis using four-acid digest HR-ICP-MS methods, in which special analytical procedures were used to eliminate 204 Hg mass overlap. Accuracy and precision were monitored by certified Pb-isotope standards and by Energy Metals' internal standards. Pb-isotope ratios measured for Energy Metals' EMST22 standard over several sample batches gave values (± 1σ precision) of 206 Pb/ 204 Pb 187 ± 3, 207 Pb/ 204 Pb 25.3 ± 0.2 and 208 Pb/ 204 Pb 40.4 ± 0.5. X-ray Diffraction (XRD) patterns of 22 powdered samples were collected on whole-rock geochemistry sample pulps with a Bruker D4 Endeavor instrument fitted with a Co tube, Fe filter, and a Lynxeye position-sensitive detector. The measured 2-theta range was 5-90 • (or 50 • ), with a step size of 0.02 • and a divergence slit of 1 • .
Supporting data and further information is provided in the Supplementary Materials.

Features of the Host Sandstone
The host sandstone is dominantly a feldspathic sandstone (subarkose) with~10-25% K-feldspar content, grading into arkose and lithic-rich variants. Textures range from grain-supported to matrix-supported with a variable content of clay matrix, lithic clasts and carbonate cement. The dominant detrital minerals are quartz and K-feldspar. The quartz grains are angular to sub-angular and partially interlocked. The K-feldspar grains, which consist mainly of microcline, are also angular and often fractured. Other detrital minerals include variably degraded or chloritized biotite, muscovite, roscoelite grading into vanadian illite, chlorite, dolomite and rare plagioclase. Heavy mineral detrital grains comprise zircon, monazite, apatite, rutile/leucoxene, Fe-Ti±V-oxides, tourmaline and garnet.
Lithic clasts are present in a variety of forms. Scattered rounded pebbles and cobbles are common throughout the MES, and consist dominantly of quartzite with less abundant sandstone, granite, quartz-muscovite schist and phyllite, limestone and dolostone. They are likely derived from fully lithified, underlying sedimentary units of the Ngalia Basin, as well as the adjacent Proterozoic basement. Lithic 'rip-up' clasts of fine-grained sedimentary rocks, comprising mudstone (shale), siltstone, sandstone, occasional micritic carbonate and composites of these, are abundant at the base of fining-upward sedimentary cycles. They are typically of irregular or elongate form due to apparent soft-sediment deformation, and may be 'smeared out' and grade into clay-rich matrix. The clasts Minerals 2020, 10, 896 6 of 29 range from a few millimeters up to 40 cm in size and are likely derived from partly lithified sediments within the fluvial system. Rip-up clasts of fine-grained sedimentary rocks are frequently the focus of, and have acted as sites for, U-V mineralization.
The clay-rich matrix, grading into mudstone rip-up lithic clasts, may be a significant component of the host sandstones. A detailed study of the clay mineralogy has shown that a complex mixture of discrete clay species is present. The clay minerals include kaolinite, chlorite, illite and smectite, together with the interstratified species corrensite (chlorite-smectite) and illite-smectite [18,27,29]. The abundances of kaolinite and corrensite usually dominate over illite and chlorite in grey sandstone.
Pyrite is the main Fe-sulfide and is common in the clay-rich matrix, where it typically occurs as fine grains of globular form, as framboidal aggregates, discrete subhedral to euhedral grains, and as zoned or skeletal forms; pyrite is also associated with degraded biotite, where it may decorate cleavage planes, and as a partial replacement product of detrital Fe-Ti oxides. Marcasite is occasionally present as interstitial aggregates intergrown with pyrite, but at Bigrlyi there is no apparent spatial relationship between marcasite and mineralization, as is characteristic of the ore zones of some Colorado Plateau U-V deposits, for example [30].
Calcite is an authigenic mineral occurring as interstitial calcite cement. The cement is variable and patchy in distribution, but in some cases may flood and occlude all available pore space. Ore minerals are sometimes found enclosed by carbonate cement, indicating they formed prior to early diagenetic cementation. A detailed description and discussion of the diagenetic history can be found in Schmid and Quigley [27].

Distribution of U-V Across the Deposit
The distributions of uranium and vanadium, as well as calcium, have been investigated across the deposit using an extensive dataset of drillhole-derived data that includes downhole geochemical assays and gamma-ray logs. Calcium has been used as a proxy for calcite to better understand the distribution of cemented sandstone. Drillhole assay data for V 2 O 5 , U 3 O 8 and Ca were compiled and composited to 1 meter intervals downhole, producing some 25,300 records. A plot of U 3 O 8 versus V 2 O 5 concentrations on this basis ( Figure 2) shows that the Bigrlyi deposit is dominated by vanadium-rich but uranium-poor intervals, with 90% of the assayed meter intervals being of this type, but overall U and V are well correlated. the deposit using an extensive dataset of drillhole-derived data that includes downhole geochemical assays and gamma-ray logs. Calcium has been used as a proxy for calcite to better understand the distribution of cemented sandstone. Drillhole assay data for V2O5, U3O8 and Ca were compiled and composited to 1 meter intervals downhole, producing some 25,300 records. A plot of U3O8 versus V2O5 concentrations on this basis ( Figure 2) shows that the Bigrlyi deposit is dominated by vanadiumrich but uranium-poor intervals, with 90% of the assayed meter intervals being of this type, but overall U and V are well correlated. Figure 2. Density data plot showing the variation in V2O5 versus U3O8 assay values for 1 meter composited drillhole intervals covering the entire Bigrlyi deposit (25,300 meter-intervals in total). The vanadium-rich halo is associated with vanadium levels in the hundreds of ppm (up to about 2000 ppm V2O5) but <200 ppm U3O8.
Visualization of this data in three dimensions shows that the deposit is defined by a large halo of V mineralization (Figure 3), which surrounds U-rich intervals, and occupies a volume more than three times larger than the U-mineralized volume. The halo is characterized by V2O5 values in the hundreds of ppm range (median 200 ppm V2O5), but with low U contents. For U-mineralized intervals above the 100 ppm U3O8 level, the median V2O5:U3O8 weight ratio for all Bigrlyi mineralization is 3.8. Visualization of this data in three dimensions shows that the deposit is defined by a large halo of V mineralization (Figure 3), which surrounds U-rich intervals, and occupies a volume more than three times larger than the U-mineralized volume. The halo is characterized by V 2 O 5 values in the hundreds of ppm range (median 200 ppm V 2 O 5 ), but with low U contents. For U-mineralized intervals above the 100 ppm U 3 O 8 level, the median V 2 O 5 :U 3 O 8 weight ratio for all Bigrlyi mineralization is 3.8. The Bigrlyi deposit contains a significant number of meter intervals of elevated V concentrations, with over 5% of the dataset falling in the 0.2% V 2 O 5 to 12% V 2 O 5 range. The maximum meter-interval U concentration is 4.5% U 3 O 8 . The Bigrlyi deposit contains a significant number of meter intervals of elevated V concentrations, with over 5% of the dataset falling in the 0.2% V2O5 to 12% V2O5 range. The maximum meter-interval U concentration is 4.5% U3O8. In terms of spatial variation, the distribution of U relative to V varies with each sub-deposit along the strike; in particular, the A4 sub-deposit is characterized by low V:U and is host to nearly 60% of the total uranium in the Bigrlyi deposit, whereas the adjacent A7/9 sub-deposit is V-rich ( Figure 3) with high V:U, and hosts over 30% of vanadium in Bigrlyi. Typical features of the U-V-Ca distribution modeled in 3D (Figure 4) support the drillhole observations showing that the mineralization is stratiform in nature, and largely confined within steeply-dipping grey sandstone, with U mineralization concentrated near the upper and lower contacts; the V-rich halo within grey sandstone and the surrounding envelope of calcite-cemented sandstone is evident. In terms of spatial variation, the distribution of U relative to V varies with each sub-deposit along the strike; in particular, the A4 sub-deposit is characterized by low V:U and is host to nearly 60% of the total uranium in the Bigrlyi deposit, whereas the adjacent A7/9 sub-deposit is V-rich ( Figure 3) with high V:U, and hosts over 30% of vanadium in Bigrlyi. Typical features of the U-V-Ca distribution modeled in 3D ( Figure 4) support the drillhole observations showing that the mineralization is stratiform in nature, and largely confined within steeply-dipping grey sandstone, with U mineralization concentrated near the upper and lower contacts; the V-rich halo within grey sandstone and the surrounding envelope of calcite-cemented sandstone is evident. In terms of spatial variation, the distribution of U relative to V varies with each sub-deposit along the strike; in particular, the A4 sub-deposit is characterized by low V:U and is host to nearly 60% of the total uranium in the Bigrlyi deposit, whereas the adjacent A7/9 sub-deposit is V-rich ( Figure 3) with high V:U, and hosts over 30% of vanadium in Bigrlyi. Typical features of the U-V-Ca distribution modeled in 3D (Figure 4) support the drillhole observations showing that the mineralization is stratiform in nature, and largely confined within steeply-dipping grey sandstone, with U mineralization concentrated near the upper and lower contacts; the V-rich halo within grey sandstone and the surrounding envelope of calcite-cemented sandstone is evident.  Because stratigraphic correlations across fluvial channels in an alluvial fan are difficult to make due to the lateral migration of the channels, a distinctive playa marker horizon (i.e., calcareous mudstone with nodular calcrete) was used as the datum for correlation. The horizon lies stratigraphically above the mineralization, and was identified based on gamma ray data and core logging. The distribution of mineralization for a transect from the A4 to A7 sub-deposits, corrected for stratigraphic level, shows that the bulk of mineralization is hosted in the basal interval of the stacked channel system about 120 to 150 m below the deposit-wide marker horizon ( Figure 5). due to the lateral migration of the channels, a distinctive playa marker horizon (i.e., calcareous mudstone with nodular calcrete) was used as the datum for correlation. The horizon lies stratigraphically above the mineralization, and was identified based on gamma ray data and core logging. The distribution of mineralization for a transect from the A4 to A7 sub-deposits, corrected for stratigraphic level, shows that the bulk of mineralization is hosted in the basal interval of the stacked channel system about 120 to 150 m below the deposit-wide marker horizon ( Figure 5).

Distribution of U-V in Drill Core and Its Relation to Sedimentology
The distribution of U-V mineralization was investigated in the context of sedimentary facies variations and cyclicity. The fining upward cycles are capped by calcite-cemented, fine-grained sandstone and siltstone at the cycle tops, which creates a relative permeability and porosity barrier. In order to show the association between facies and U-V mineralization, a series of spot geochemical measurements was undertaken on drill core using a portable XRF (pXRF).

Drill Core-B07120
A set of pXRF measurements was taken across a mineralized zone over a core interval of 14 m ( Figure 6). The interval is associated with intense calcite cementation enveloping the mineralization. Visual observations were confirmed by the geochemical measurements, which show elevated Ca (20- Figure 5. Distribution of mineralization for a transect from the A4 to A7 sub-deposits corrected for stratigraphic level using the upper marker bed (grey dashed line). Uranium abundance in yellow, vanadium abundance in black.

Distribution of U-V in Drill Core and Its Relation to Sedimentology
The distribution of U-V mineralization was investigated in the context of sedimentary facies variations and cyclicity. The fining upward cycles are capped by calcite-cemented, fine-grained sandstone and siltstone at the cycle tops, which creates a relative permeability and porosity barrier. In order to show the association between facies and U-V mineralization, a series of spot geochemical measurements was undertaken on drill core using a portable XRF (pXRF).

Drill Core-B07120
A set of pXRF measurements was taken across a mineralized zone over a core interval of 14 m ( Figure 6). The interval is associated with intense calcite cementation enveloping the mineralization. Visual observations were confirmed by the geochemical measurements, which show elevated Ca (20-23% as calcite cement) on either side of the mineralization ( Figure 6). The mineralized zone is defined by elevated V (4510 ppm) and U (3973 ppm), as well as Fe (1.6%) and Al (6.8%) concentrations, due to the abundance of phyllosilicates.

Drill Core-B07129
A set of pXRF measurements was taken over a ~40 cm interval within a mineralized zone across an erosional contact ( Figure 7). The interval below the contact is a medium-to fine-grained grey sandstone with a very thin micaceous top (contact). The measurement 10 cm below the contact is weakly V-mineralized (367 ppm), with relatively low U (44 ppm). The base of the overlying cycle comprising pink-grey arkosic sandstone shows a sharp increase in U-V mineralization (U 2184 ppm, V 809 ppm). Approximately 10 cm above the contact, the sandstone is black and sooty in appearance due to strong enrichment in V (5851 ppm), while U concentrations are unchanged from the interval below. It is noticeable that the Fe concentrations also increase with V, which correlates with a relative enrichment of phyllosilicates, as seen visually ( Figure 7).

Drill Core-B07129
A set of pXRF measurements was taken over a~40 cm interval within a mineralized zone across an erosional contact ( Figure 7). The interval below the contact is a medium-to fine-grained grey sandstone with a very thin micaceous top (contact). The measurement 10 cm below the contact is weakly V-mineralized (367 ppm), with relatively low U (44 ppm). The base of the overlying cycle comprising pink-grey arkosic sandstone shows a sharp increase in U-V mineralization (U 2184 ppm, V 809 ppm). Approximately 10 cm above the contact, the sandstone is black and sooty in appearance due to strong enrichment in V (5851 ppm), while U concentrations are unchanged from the interval below. It is noticeable that the Fe concentrations also increase with V, which correlates with a relative enrichment of phyllosilicates, as seen visually ( Figure 7).

Mineralogy and Petrology
To investigate U-V mineral associations, mineralized and non-mineralized intervals were studied in more detail. Powder XRD results show that samples of host sandstone without mineralization are composed of quartz, K-feldspar (microcline), muscovite, biotite, illite, kaolinite, calcite, and chlorite/corrensite, with traces of dolomite, pyrite and hematite (in partly oxidized sandstone). Samples dominated by U-V mineralization contain in addition to the above coffinite, uraninite, montroseite and roscoelite. Samples dominated by V>U mineralization were also found to contain, in addition to montroseite, the vanadium oxide-hydroxide minerals haggite and doloresite (see Data Appendix Figure S1). The mineralogy was further refined by petrographic investigations.

Uranium-Dominated U-V Mineralization
Uraninite and coffinite are the predominant uranium minerals, and are mainly concentrated in clay-rich lithic clasts and the clay-rich matrix ( Figure 8a). QEMSCAN (automated mineralogical analysis) studies of mineralized samples indicate that coffinite is generally the dominant species, with the uraninite to coffinite ratio varying from about 1:1 to 1:4 [31]. Uraninite typically occurs as finegrained, irregular or globular aggregates, as disseminations intergrown with clays ( Figure 8a,b), and with chloritized/degraded biotite in association with montroseite ( Figure 8c). Uraninite may also rim and partly replace detrital rutile and Fe-Ti-V oxide grains (Figure 8a,f), and may be found as a replacement product of early-formed pyrite. Coffinite characteristically encrusts the margins of quartz grains (Figure 8a,b) and degraded biotite; in more porous samples it coats grain contacts and internal fractures within K-feldspar (Figure 8d,e). Coffinite may encase or replace fine globulartextured uraninite leaving disc-shaped remnants ( Figure 8b). Galena is a common fine-grained inclusion in coffinite, which likely reflects the exsolved former radiogenic Pb solid-solution content.
Detrital Fe-Ti-oxide grains have undergone a process of replacement by V and U ( Figure 8f). Element mapping shows the pervasive replacement of Fe by V on the grain margins and along internal fractures; there is growth of uraninite crystallites on the grain rims, and infiltration of fine uraninite into V-enriched areas along micro-scale fractures.
Under cathodoluminescence (CL) imaging, radiation-induced damage is noticeable on Kfeldspar grains, and to a lesser extent on quartz grains, where zones of quenched CL emissions on grain rims and internal fractures are related to the proximity of uranium minerals (Figure 8e). Of particular interest in Figure 8e is the presence of obvious radiation-induced damage affecting Kfeldspar and quartz, where no uranium mineralization is currently visible (left-hand part of Figure  8e).

Mineralogy and Petrology
To investigate U-V mineral associations, mineralized and non-mineralized intervals were studied in more detail. Powder XRD results show that samples of host sandstone without mineralization are composed of quartz, K-feldspar (microcline), muscovite, biotite, illite, kaolinite, calcite, and chlorite/corrensite, with traces of dolomite, pyrite and hematite (in partly oxidized sandstone). Samples dominated by U-V mineralization contain in addition to the above coffinite, uraninite, montroseite and roscoelite. Samples dominated by V>U mineralization were also found to contain, in addition to montroseite, the vanadium oxide-hydroxide minerals haggite and doloresite (see Supplementary Materials Figure S1). The mineralogy was further refined by petrographic investigations.

Uranium-Dominated U-V Mineralization
Uraninite and coffinite are the predominant uranium minerals, and are mainly concentrated in clay-rich lithic clasts and the clay-rich matrix (Figure 8a). QEMSCAN (automated mineralogical analysis) studies of mineralized samples indicate that coffinite is generally the dominant species, with the uraninite to coffinite ratio varying from about 1:1 to 1:4 [31]. Uraninite typically occurs as fine-grained, irregular or globular aggregates, as disseminations intergrown with clays ( Figure 8a,b), and with chloritized/degraded biotite in association with montroseite ( Figure 8c). Uraninite may also rim and partly replace detrital rutile and Fe-Ti-V oxide grains (Figure 8a,f), and may be found as a replacement product of early-formed pyrite. Coffinite characteristically encrusts the margins of quartz grains (Figure 8a,b) and degraded biotite; in more porous samples it coats grain contacts and internal fractures within K-feldspar (Figure 8d,e). Coffinite may encase or replace fine globular-textured uraninite leaving disc-shaped remnants (Figure 8b). Galena is a common fine-grained inclusion in coffinite, which likely reflects the exsolved former radiogenic Pb solid-solution content.
Detrital Fe-Ti-oxide grains have undergone a process of replacement by V and U (Figure 8f). Element mapping shows the pervasive replacement of Fe by V on the grain margins and along internal fractures; there is growth of uraninite crystallites on the grain rims, and infiltration of fine uraninite into V-enriched areas along micro-scale fractures.
Under cathodoluminescence (CL) imaging, radiation-induced damage is noticeable on K-feldspar grains, and to a lesser extent on quartz grains, where zones of quenched CL emissions on grain rims and internal fractures are related to the proximity of uranium minerals (Figure 8e). Of particular interest in Figure 8e is the presence of obvious radiation-induced damage affecting K-feldspar and quartz, where no uranium mineralization is currently visible (left-hand part of Figure 8e).

Vanadium-Dominated U-V Mineralization
Vanadium-uranium mineralization at Bigrlyi is primarily associated with the clay-rich matrix, with clay-rich lithic clasts, and with detrital micas and chlorite imbedded in the matrix. The main V- BSE-EDS micrographs illustrating U-mineral textures.
(a) Fine, disseminated uraninite/coffinite (pink) resides within Fe-clay-rich matrix and partly replaces Fe-Ti-V oxides (green); coffinite also replaces the margins of quartz grains; (b) on a quartz-mica grain contact, coffinite (light grey) replaces quartz and invades mica cleavage planes; globular uraninite has been partly replaced by coffinite leaving disc-shaped remnants of about 5 microns diameter; (c) uraninite crystallites (red) decorate the cleavage planes of altered biotite (dark blue) in association with bladed montroseite crystals (light blue-green); (d) coffinite (pink) fills fractures in K-feldspar and replaces the margins of some quartz grains; fine pyrite (light blue) occurs separately from coffinite mainly within the clay-rich matrix and as thin coatings on some grain boundaries.

Vanadium-Dominated U-V Mineralization
Vanadium-uranium mineralization at Bigrlyi is primarily associated with the clay-rich matrix, with clay-rich lithic clasts, and with detrital micas and chlorite imbedded in the matrix. The main V-host minerals are Fe-V-bearing clays and chlorite, including roscoelite grading into vanadian illite, corrensite and smectite, plus altered biotite, where vanadium occurs as a solid-solution component. Vanadium is also hosted in the V-Fe-oxyhydroxide minerals montroseite, doloresite and haggite, and may be present as a component of altered Fe-Ti oxides and rutile/leucoxene. Texturally, montroseite occurs within the clay-rich matrix as disseminated bladed crystals grading into massive aggregates of fine-grained, bladed montroseite up to a few millimeters across. Montroseite also occurs as crystal aggregates replacing altered/chloritized detrital biotite (Figure 8c). Montroseite is locally associated with pyrite, uraninite/coffinite and rutile/leucoxene. Occasional zoned pyrite grains intergrown with montroseite are found imbedded within the matrix.
QEMSCAN studies of V-U-mineralized samples confirm that vanadiferous phyllosilicates are the dominant V-hosting species. In a representative V>U ore material containing 0.5% V 2 O 5 and 0.1% U 3 O 8 from the A4 sub-deposit, the ratio of vanadiferous aluminosilicates to montroseite was found to be~3:1, with about half the vanadium being associated with V-chlorite/corrensite [31]. In agreement with petrographic observations, V-bearing detrital biotite is more abundant in the coarse fractions, whereas the opposite is true for chlorite/corrensite and montroseite. QEMSCAN inter-mineral correlations indicate that montroseite has a strong association with both vanadiferous biotite and chlorite.
On the thin-section scale, optical and XRF elemental mapping studies of V>U-mineralized sandstones show that lithic rip-up clasts may host a significant proportion of U-V mineralization. Mineralized samples from drill core B07270 (A15 sub-deposit), for example, contain a texturally and chemically heterogeneous mixture of different kinds of deformed lithic clasts, including mudstone, siltstone and fine-to-coarse grained sandstone ( Figure 9). The clasts are heterogeneous with respect to the distribution of mineralizing elements, and U-V-Fe concentrations are highly variable; some clasts are Fe-rich and only weakly U-V-mineralized, while other clasts are strongly V-and/or U-mineralized. A proportion of clasts are chemically zoned and/or are composites containing variably enriched Fe, V and U domains (Figure 10a). Some clasts are also distinct in that they are rich in heavy mineral content or may contain patches of calcite cement, features that reflect their diverse origins. Compared to the lithic clasts, the sandstone matrix in these specimens is minimally mineralized; U-V-Fe mineralization is scattered amongst interstitial clay pockets that likely represent disaggregated particles abraded from the larger clasts.
In the example shown in Figure 9, the various lithic clasts may be assigned to different types based on their contents of mineralizing elements and textures. These observations are significant because they demonstrate that, at least in parts of the deposit, U-V mineralization did not to any large extent take place in-situ in the host sandstone, but was introduced in a variety of mostly fine-grained, clay-rich pre-existing sediments. The diversity of clasts present indicates that the source sediment was variable in nature and in mineralization type. Based on the observations made, they were likely poorly consolidated sediments that were 'ripped-up' and re-deposited down-channel in the stacked fluvial channel system. As discussed above, lithic clasts containing chlorite and biotite may be variably mineralized, which is also observed petrographically, as shown in Figure 10a, where a V-rich lithic clast is in contact with a composite clast containing a U-rich domain. Detail from another V-rich clast shows that microcrystalline pyrite is associated with altered biotite, while uraninite is present as disseminated grains imbedded in V-chlorite (Figure 10b). In some lithic clasts, detrital angular Kfeldspar, quartz and dolomite grains, and minor Ti-V-Fe oxides, have been engulfed by pore-filling, V-bearing chlorite, which is overgrown by calcite cement in places. Zones that are strongly calcitecemented typically comprise angular grains of quartz and K-feldspar, and remnant mineralized domains floating in the calcite cement, implying calcite cementation occurred before compaction (Figure 10c). Some V-chlorite-rich clasts, and the intergranular matrix of other clasts, are variably mineralized and contain aggregates of montroseite (Figure 10d) or other oxyhydroxide V-minerals (doloresite, haggite). Radiation-induced damage of K-feldspar is apparent adjacent to, but not necessarily in contact with, U mineralization (Figure 10d), but is absent in calcite-cemented detrital grains lacking clay mineral content. As discussed above, lithic clasts containing chlorite and biotite may be variably mineralized, which is also observed petrographically, as shown in Figure 10a, where a V-rich lithic clast is in contact with a composite clast containing a U-rich domain. Detail from another V-rich clast shows that microcrystalline pyrite is associated with altered biotite, while uraninite is present as disseminated grains imbedded in V-chlorite (Figure 10b). In some lithic clasts, detrital angular K-feldspar, quartz and dolomite grains, and minor Ti-V-Fe oxides, have been engulfed by pore-filling, V-bearing chlorite, which is overgrown by calcite cement in places. Zones that are strongly calcite-cemented typically comprise angular grains of quartz and K-feldspar, and remnant mineralized domains floating in the calcite cement, implying calcite cementation occurred before compaction (Figure 10c). Some V-chlorite-rich clasts, and the intergranular matrix of other clasts, are variably mineralized and contain aggregates of montroseite (Figure 10d) or other oxyhydroxide V-minerals (doloresite, haggite). Radiation-induced damage of K-feldspar is apparent adjacent to, but not necessarily in contact with, U mineralization (Figure 10d), but is absent in calcite-cemented detrital grains lacking clay mineral content. Pyrite may be present in different generations in mineralized samples. Early-formed microcrystalline skeletal pyrite is commonly imbedded in the V-clay-rich matrix, and overgrown and partly replaced by uraninite (Figure 10f). Late-formed pyrite may occur as large aggregates occupying pore space and coating grain boundaries. In some instances, pyrite aggregates are associated with strongly U-V-mineralized clay clasts, and are enclosed by calcite cement (Figure 10e). In some V-rich samples that lack uranium or evidence of previously residing uranium, a different detrital provenance is recognized among the lithic clasts. The clasts are composed of rounded quartz grains and minor K-feldspar without radiation-induced damage bonded by silica Pyrite may be present in different generations in mineralized samples. Early-formed microcrystalline skeletal pyrite is commonly imbedded in the V-clay-rich matrix, and overgrown and partly replaced by uraninite (Figure 10f). Late-formed pyrite may occur as large aggregates occupying pore space and coating grain boundaries. In some instances, pyrite aggregates are associated with strongly U-V-mineralized clay clasts, and are enclosed by calcite cement (Figure 10e).
In some V-rich samples that lack uranium or evidence of previously residing uranium, a different detrital provenance is recognized among the lithic clasts. The clasts are composed of rounded quartz grains and minor K-feldspar without radiation-induced damage bonded by silica cement (Figure 11a,b). Pore spaces and fractures are occupied by V-chlorite and calcite cement. The clasts likely represent recycled silicified sandstone or quartzite sourced from older stratigraphic units. On the other hand, some V-rich samples containing detrital, immature, interlocked K-feldspar and quartz grains show intense radiation-induced damage, but lack any significant uranium content (Figure 11c,d). This is further illustrated in Figure 11e,f which show a K-feldspar grain isolated by calcite cement with cracks and fractures that are filled by pyrite and V-bearing clay minerals (Figure 11f), but uranium minerals are absent. Radiation-induced damage, however, is visible in fractures and along the edges of the grain (Figure 11f), suggesting it is inherited and hints at a former significant uranium content in a grain of granitic provenance.
Minerals 2020, 10, 896 15 of 28 cement (Figure 11a,b). Pore spaces and fractures are occupied by V-chlorite and calcite cement. The clasts likely represent recycled silicified sandstone or quartzite sourced from older stratigraphic units. On the other hand, some V-rich samples containing detrital, immature, interlocked K-feldspar and quartz grains show intense radiation-induced damage, but lack any significant uranium content (Figure 11c,d). This is further illustrated in Figure 11e,f which show a K-feldspar grain isolated by calcite cement with cracks and fractures that are filled by pyrite and V-bearing clay minerals ( Figure  11f), but uranium minerals are absent. Radiation-induced damage, however, is visible in fractures and along the edges of the grain (Figure 11f), suggesting it is inherited and hints at a former significant uranium content in a grain of granitic provenance.

Geochemical Signature of U-V Mineralization
The whole-rock major and trace-element geochemistry of the Bigrlyi deposit has been investigated in order to better understand the geochemical signatures related to U-V mineralization and to further test the mineralogical correlations observed in drill core.
Major and trace element variations across U-V-mineralized zones have been explored using the isocon plot method after Grant [32]. Originally, isocon plots were used to characterize metasomatic alteration, showing mass transfer between mineralized and non-mineralized zones, if the host-rock was geochemically identical. The Bigrlyi sandstones are not metasomatically altered, however similar principles can be applied. A non-mineralized sample, of similar grain size and bulk major element contents including Si, Al, K and Ca, has been used as the host-rock reference sample. Ratios of barren versus mineralized samples were calculated, and plotted for a transect across the mineralized zone from 191.3 to 200.3 m in drill core B07120 from the A4 sub-deposit (Figure 12a-d). The concentrations of Se, Li, Ba, Be, Mo, Mg and Fe, in relation to these elements in the non-mineralized sample, show marked enrichments across the mineralized interval, whereas S is strongly depleted.

Geochemical Signature of U-V Mineralization
The whole-rock major and trace-element geochemistry of the Bigrlyi deposit has been investigated in order to better understand the geochemical signatures related to U-V mineralization and to further test the mineralogical correlations observed in drill core.
Major and trace element variations across U-V-mineralized zones have been explored using the isocon plot method after Grant [32]. Originally, isocon plots were used to characterize metasomatic alteration, showing mass transfer between mineralized and non-mineralized zones, if the host-rock was geochemically identical. The Bigrlyi sandstones are not metasomatically altered, however similar principles can be applied. A non-mineralized sample, of similar grain size and bulk major element contents including Si, Al, K and Ca, has been used as the host-rock reference sample. Ratios of barren versus mineralized samples were calculated, and plotted for a transect across the mineralized zone from 191.3 to 200.3 m in drill core B07120 from the A4 sub-deposit (Figure 12a-d). The concentrations of Se, Li, Ba, Be, Mo, Mg and Fe, in relation to these elements in the non-mineralized sample, show marked enrichments across the mineralized interval, whereas S is strongly depleted. A set of multi-element geochemical data from barren to highly mineralized sandstone, categorized by mineralization type, was used to further explore chemical variance via the Principal Component Analysis (PCA) method. PCA allows complex multivariate data to be displayed in a lower dimensional space to simplify and aid interpretation. From inspection of an initial correlation A set of multi-element geochemical data from barren to highly mineralized sandstone, categorized by mineralization type, was used to further explore chemical variance via the Principal Component Analysis (PCA) method. PCA allows complex multivariate data to be displayed in a lower dimensional space to simplify and aid interpretation. From inspection of an initial correlation matrix, it was found that certain elements were strongly inter-correlated, reflecting similar geochemical behavior. For the purposes of the analysis, these elements were grouped together as follows: chalcophile-related, Se + As + Mo + Sb; monazite/apatite-related, La + Ce + Y (REE) + P; heavy mineral-related, Zr + Nb + Ti + Th + Sn; mafic/ultramafic-related, Cr + Ni + Co + Cu; and biotite/chlorite/corrensite-related, Fe + Mg. The data were log-transformed, and based on the correlation matrix, the principle components were extracted using the Fast MCD algorithm. Table 1 provides a summary of key element abundances and ratios for each mineralization category. Further data are provided in Supplementary Data Table S1. The PCA analysis shows that the elements Ca, Sr and Mn behave similarly and have high variance vectors that clearly relate to the variable abundance of calcite cement in the sandstones ( Figure 13). Aluminum shows the opposite behavior due to the lower contents of clays and K-feldspar in calcite-cemented sandstone. The U vector is accompanied by a similar vector for Pb, reflecting radiogenic lead production from U decay. Other U-associated vectors include the chalcophile elements Se + As + Mo + Sb. The Ba and Be vectors also show a strong relationship with U (see also Table 1). Of interest is the decoupling seen between the Se + As + Mo + Sb and S vectors, with Se/S ratios, for example, showing an increase of two orders of magnitude in strongly U-mineralized samples compared to barren sandstone ( Table 1). The V vector is accompanied by a similar vector for Li, which is presumably a result of the preferential concentration of Li in the phyllosilicate minerals that host most of the V. Strongly V-mineralized samples are influenced by the Fe + Mg vector, which reflects the importance of chlorite and related minerals as hosts for V. The vector for K is anticorrelated with V, which likely relates to the enhanced alteration and chloritization of detrital biotite and lower K-feldspar-to-clay-mineral ratios in V-mineralized samples; this is also seen in the lower K/Al ratios of mineralized samples (Table 1).

Organic Carbon
In view of the importance of organic carbon as a reductant in tabular U-V sandstone-hosted deposits of the Colorado Plateau, and its apparent rarity in petrographic studies of the Bigrlyi deposit, the Bigrlyi geochemical dataset was interrogated for organic carbon analyses from samples of barren and mineralized sandstone and siltstone/shale ( Table 2). The results confirm observations, showing that only low levels of organic carbon are present, with a maximum value of 0.46% C recorded from a weakly mineralized shale. Although the shales have slightly higher organic carbon contents than the sandstones, no association between the content of carbonaceous matter and mineralization is evident ( Table 2).

Organic Carbon
In view of the importance of organic carbon as a reductant in tabular U-V sandstone-hosted deposits of the Colorado Plateau, and its apparent rarity in petrographic studies of the Bigrlyi deposit, the Bigrlyi geochemical dataset was interrogated for organic carbon analyses from samples of barren and mineralized sandstone and siltstone/shale ( Table 2). The results confirm observations, showing that only low levels of organic carbon are present, with a maximum value of 0.46% C recorded from a weakly mineralized shale. Although the shales have slightly higher organic carbon contents than the sandstones, no association between the content of carbonaceous matter and mineralization is evident ( Table 2).

Pb Isotopes
The Pb isotopes 206 Pb and 207 Pb are the ultimate radiogenic daughter products of radioactive decay in the 238 U and 235 U decay chains, respectively. In uranium deposits such as Bigrlyi that are old enough to allow sufficient ingrowth of radiogenic Pb, there is a strong correlation between the abundances of Pb and U ( Figure 13 and Table 1). Because U and Pb have different chemical behaviors, they may be affected to different extents by post-depositional processes such as oxidative leaching. The abundance of radiogenic Pb relative to U may therefore be used to provide an insight into the differential mobility of these elements over time since deposit formation. If the mineralized system remained closed and undisturbed since formation, i.e., U and Pb remained immobile, then the amount of radiogenic Pb present will be proportional to the uranium content of the rock and the time elapsed since the mineralization event; in such cases, the U-Pb isotopic system will exhibit concordant behavior. If open system processes disturbed the ore body and mobilized U or Pb post-mineralization, then the U-Pb isotopic system will exhibit discordant behavior. Measurements of the Pb-isotopic composition of a sample can be used to quantify the mobility (gain or loss) of U relative to Pb. As discussed, in the Bigrlyi deposit there is petrographic evidence for post-mineralization U mobility, which is expected to have disturbed the U-Pb isotopic system. To explore the nature of these post-mineralization processes, a whole-rock Pb-isotope dataset, that covers the range of U-V mineralization types, has been compiled (Supplementary Data Table S2). Standard equations were used to calculate the abundance of radiogenic Pb and a common Pb isotope correction was applied to calculate the concentrations of the individual radiogenic isotopes ( 206 Pb*, 207 Pb* and 208 Pb*). Other key isotope ratios were then determined. For the Bigrlyi deposit, Th is of negligible abundance relative to U, and therefore 208 Pb*, the daughter of the Th-decay chain, is not a significant component of this system.
The whole-rock U-Pb results show that the majority (70%) of the data points are significantly discordant, and plot above and below the concordia curve in 207 Pb*/ 206 Pb* versus 238 U/ 206 Pb* space; the scatter of data points confirms that parts of the Bigrlyi deposit have experienced time-integrated gains or losses of U and/or Pb from the system (Table S2). Significantly, however, 11 of the data points plot within a few percent of the concordia curve, do not show evidence for significant U gain or loss, and using the 238 U-206 Pb chronometer give a geologically reasonable, broad timing of mineralization in the range 282 to 335 Ma (i.e., mid-to-late Carboniferous) at the 2σ level ( 238 U/ 206 Pb* 20.5 ± 1.8). This provides confidence in the quality and internal consistency of the dataset.
The amount of radiogenic Pb present in a sample was used to calculate the expected uranium content of that sample, assuming 300 Ma of evolution of the U-Pb decay chains and closed-system behavior (Supplementary Materials Table S2). Any gain or loss of U relative to Pb in a sample can be expressed as a percent deviation of the present-day uranium content from the expected uranium concentration. As determined from reproducibility of the internal standard, the precision of the percentage deviation is ±8% absolute at the 2σ level, so that any deviation above this level is expected to be of significance. As a simple rule of thumb, the ratio of U to total radiogenic Pb (U/Pb*), as listed in Table S2, is indicative of closed system behavior when that value is near 20, and deviates to lower or higher values where there has been mobility of radiogenic elements (i.e., open system behavior).

Pb-Isotope Constraints on U-Mobility and Open System Behavior
Petrographic observations showing radiation-induced damage of detrital minerals in the absence of associated uranium provide evidence for the mobility of uranium, and hence for open system behavior for radiogenic elements over time. Since V is likely to be largely immobile as a component of aluminosilicate minerals, the mobilization of U relative to V and its re-precipitation could be a significant process in upgrading U-rich parts of the deposit following sedimentary deposition. This view is supported and reinforced by the Pb-isotope results, which provide evidence for time-integrated U mobility across the deposit. The deviation between measured U, and U calculated on the basis of radiogenic Pb content, has been explored as a function of V/U ratio in Figure 14, with the data points categorized by mineralization type as per Table 1. The following interpretations can be made: • About 30 % of the Pb-isotope dataset, including many of the strongly U-mineralized samples, show closed system behaviors within error; presumably these mineralization sites were protected from U-and/or Pb-mobility by impermeable barriers soon after formation, preventing the infiltration of oxidized fluids post-mineralization; • About 45% of the Pb-isotope dataset shows evidence for time-integrated U-loss relative to Pb (i.e., these samples have radiogenic Pb concentrations that are unsupported by current U levels). This category mainly comprises samples from the vanadium halo and from strongly V-mineralized zones; • The remaining 25% of the Pb-isotope dataset show evidence for time-integrated U-gain relative to Pb (i.e., these samples have U levels that are unsupported by their current radiogenic Pb contents).
Minerals 2020, 10, 896 20 of 28 14, with the data points categorized by mineralization type as per Table 1. The following interpretations can be made: • About 30 % of the Pb-isotope dataset, including many of the strongly U-mineralized samples, show closed system behaviors within error; presumably these mineralization sites were protected from U-and/or Pb-mobility by impermeable barriers soon after formation, preventing the infiltration of oxidized fluids post-mineralization; • About 45% of the Pb-isotope dataset shows evidence for time-integrated U-loss relative to Pb (i.e., these samples have radiogenic Pb concentrations that are unsupported by current U levels). This category mainly comprises samples from the vanadium halo and from strongly Vmineralized zones; • The remaining 25% of the Pb-isotope dataset show evidence for time-integrated U-gain relative to Pb (i.e., these samples have U levels that are unsupported by their current radiogenic Pb contents). Figure 14. Percent deviation of uranium concentration from that expected for 300 Ma of closed system uranium decay, based on the concentration of present-day radiogenic Pb plotted versus sample V/U ratio. Data points are categorized by mineralization type as in Table 1. Samples showing timeintegrated U-loss relative to Pb (i.e., samples with radiogenic Pb unsupported by current U levels) plot below the 0% line, and mainly comprise those from the vanadium halo and from strongly Vmineralized zones.
In Table 3, some examples of downhole variation are used to illustrate the scale of radiogenic element mobility. The data show that variations between closed system and open system behavior may fluctuate downhole over scales, from ~0.4 m up to several meters' width, with U-loss relative to Pb being most common in intervals of high V/U. Note also that the lithic-clast-rich, V>U-mineralized sample B07270-227 m, previously described (Figure 9), shows strictly closed-system behavior, confirming this sample as an example of unmodified mineralized sandstone.  Figure 14. Percent deviation of uranium concentration from that expected for 300 Ma of closed system uranium decay, based on the concentration of present-day radiogenic Pb plotted versus sample V/U ratio. Data points are categorized by mineralization type as in Table 1. Samples showing time-integrated U-loss relative to Pb (i.e., samples with radiogenic Pb unsupported by current U levels) plot below the 0% line, and mainly comprise those from the vanadium halo and from strongly V-mineralized zones.
In Table 3, some examples of downhole variation are used to illustrate the scale of radiogenic element mobility. The data show that variations between closed system and open system behavior may fluctuate downhole over scales, from~0.4 m up to several meters' width, with U-loss relative to Pb being most common in intervals of high V/U. Note also that the lithic-clast-rich, V>U-mineralized sample B07270-227 m, previously described (Figure 9), shows strictly closed-system behavior, confirming this sample as an example of unmodified mineralized sandstone. The Pb-isotope results suggest that over time U has been mobilized from the vanadium halo and from strongly V-mineralized zones and redistributed to other parts of the deposit, resulting in zones of later U-enrichment. These mobility pathways were presumably open from early in the history of the Bigrlyi deposit, although the earliest events, occurring immediately post-deposition, may not necessarily be apparent due to the time required for the ingrowth of sufficient radiogenic Pb to be detectable. The timing and history of U-Pb mobility is not known, but a number of events in the post-depositional period are considered likely in response to tectonic movements or climatic/hydrological changes. Open system processes do not appear to have significantly affected the Bigrlyi deposit within the last 1 Ma or so, since the U-decay chains in unoxidized, sub-surface parts of the Bigrlyi deposit are in secular equilibrium within error [33], and the ratios of 238 U to short-lived daughter radionuclides, including 210 Pb, have not been disturbed in recent times.

Nature and Timing of Bigrlyi U-V Mineralization
In this study, an integrated sedimentological, mineralogical and geochemical investigation of the distribution of U-V mineralization in the Bigrlyi deposit has been presented. A summary of key observations and interpretations is outlined below.
Sedimentologically, U-V mineralization is concentrated in the basal parts of stacked fluvial channels within the locally reduced facies of an immature, arkosic sandstone (grey sandstone). Uranium-rich and/or V-rich mineralization is mainly concentrated near redox boundaries at the top and base of the host grey sandstone, and a broad halo of lower grade V mineralization extends through the grey sandstone unit. Minor remnant U mineralization and broader zones of V-enrichment occur in the underlying and overlying oxidized, hematitic sandstone units. Along the strike there are notable differences in the bulk V:U ratios of the various Bigrlyi sub-deposits, with the A−7/9 sub-deposit being V-rich and U-poor. Within sub-deposit variation of V:U ratios is partly the result of the remobilization of U, as Pb-isotope data and petrographic evidence show. Pb-isotope studies and unsupported radiation-damage-related features in detrital K-feldspar show that fluid pathways did exist within ore zones, allowing access for oxidized fluids and consequent uranium mobilization. Over time, open system element mobility affected some areas of the deposit more than others, including the V-halo and V-rich zones, which show U-loss relative to Pb, with the mobilized U being re-concentrated in other zones.
Mineralogically, U-V mineralization is strongly associated with phyllosilicate minerals in lithic rip-up clasts and in the clay-rich matrix of the host reduced sandstone. The clay-rich lithic clasts are comprised of mostly fine-grained, partly lithified sediments (mudstone, siltstone and fine sandstone) that are texturally and chemically heterogeneous. Some lithic clasts are rich in V, some are rich in U, and some may be zoned with U-and V-rich domains; other clasts are unmineralized and contain Fe-rich phyllosilicates and pyrite. The lithic rip-up clasts were likely derived from a diverse range of heterogeneously mineralized, reduced, mainly fine-grained sediments, which were originally deposited up-channel of their present depositional site. Fluvial channel stacking caused by high-energy events was a result of erosion and scouring of older channels and overbank deposits, and the downstream redeposition of the eroded materials.
The U-V mineralization was in place prior to the precipitation of early diagenetic, groundwater-related calcite cements. The lack of porosity caused by the pore-occluding calcite cement prevented the later infiltration of fluids and oxidative mobilization in large parts of the deposit.
Vanadium occurs dominantly as a V (III) solid-solution component of Fe-bearing phyllosilicates (roscoelite/illite, chlorite/corrensite and altered biotite), which are present as constituents of lithic rip-up clasts, the intergranular matrix, or gradational zones between the two. The Fe and V-bearing phyllosilicates are compositionally variable and heterogeneously distributed, while the weakly mineralized zones may contain an abundance of Fe-bearing-only phyllosilicates. Vanadium-oxyhydroxide minerals such as montroseite and haggite are less important hosts for V. Montroseite typically occurs as a replacement phase of degraded detrital biotite, or forms aggregates of bladed crystals within the clay-rich matrix or lithic clasts. Detrital Fe-Ti minerals contain variable amounts of V due to alteration and partial replacement of Fe by V, but the timing of the alteration process is not known.
Uranium in the unoxidized Bigrlyi sub-deposits is present as coffinite and uraninite, with the former dominant. Uraninite commonly occurs as fine-grained disseminations and globular aggregates in association with phyllosilicates, including degraded detrital biotite, and is texturally an early-formed phase. Coffinite is generally a later phase present as massive aggregates replacing uraninite, as encrustations and overgrowths on grain boundaries of detrital quartz, and as fracture-fillings in K-feldspar; the presence of galena inclusions suggests it was recrystallized in later events, releasing and re-precipitating contained radiogenic Pb. Where uranium mineralization is present, pyrite is typically less abundant than in unmineralized zones.
Geochemically, U-V mineralization is characterized by the enrichment of trace Li accompanied by elevated Fe + Mg, relating to the abundance of phyllosilicates; it is possible that the association of Li with mineralizing elements is a source-related feature, reflecting a common provenance of detrital phyllosilicates. Potassium concentrations and K/Al ratios are lower in the mineralized zone, likely due to the higher abundance of phyllosilicates relative to K-feldspar. In strongly U-mineralized zones chalcophile elements, such as Se, are significantly enriched relative to sulfur, resulting in characteristically elevated Se/S ratios; the effect is present but less pronounced in V>U-mineralized zones. Uranium mineralization is also correlated with elevated Be and Ba, but it is unclear whether this association is related to a common source-region provenance or to preferential concentration during the mineralization process. Lastly, V mineralization is poorly correlated and decoupled from the mafic/ultramafic elements, suggesting a mafic igneous source for V is less likely than a black shale or hydrothermal source.

Implications for Mineralization Models
As petrographic observations suggest, vanadium was largely introduced, as detrital clasts or disseminated V (III) -bearing phyllosilicates (roscoelite, chlorite, biotite), into the Bigrlyi host sandstone as detritus eroded from sites located up-channel and in the exposed basement of the hinterland. It is probable that localized diagenetic decomposition of the V-bearing minerals resulted in the formation of vanadium oxyhydroxides, such as montroseite (V (III) -bearing), haggite (V (III) ,V (IV) -bearing) and doloresite (V (IV) -bearing), prior to calcite cementation and under mostly reducing conditions. Therefore, V mineralization at Bigrlyi did not take place in-situ in the host sandstone by processes involving the transport of soluble V in groundwaters and its precipitation at the site of deposition.
Further, as early-formed U-minerals are intimately associated with Fe-V-phyllosilicates, such as those occurring in heterogenous U-enriched lithic clasts, it is proposed that at least a portion of Bigrlyi U was also originally introduced into the host sandstone as detrital components. Original up-channel sites associated with mineralization most likely comprised reduced, clay-rich sediment intervals; at the site of deposition, prolonged reducing conditions preserved the mineralized zones.
While parts of the deposit remained under reducing conditions due to the lack of porosity and permeability (chemical compaction by calcite-cemented beds, and mechanical compaction), at least some zones were open to fluid exchange. Evidence provided from both petrography and Pb-isotope studies indicate that fluid access pathways existed over long time periods within the V-halo and within other V-rich mineralized domains of the Bigrlyi deposit. Such pathways permitted, over time, the variable access of oxidized fluids into the grey sandstone host, resulting in mobilization and the loss of U from the V-bearing domains, and its re-concentration elsewhere in the deposit, while V remained essentially fixed as an immobile component of aluminosilicate minerals. Patches of V mineralization, associated with weak remnant U mineralization located in oxidized sandstone stratigraphically above and below the main host, hint at the more widespread mobilization of U outside ore zones.
The presence of earlier sites of U-V mineralization located in up-channel sediments echo the multiple migration-accretion hypothesis of Gruner [34] and the pre-concentration model of Finch and Davis [35]. What kind of processes could be responsible for mineralization in the earlier deposited sediments? Gruner proposed that the "Concentration of large deposits could proceed by several stages of oxidation-solution-migration-accretion, a kind of recycling action" involving multiple cycles of oxidative dissolution of U and V and the re-precipitation of these elements by interaction with carbonaceous matter or hydrogen sulfide. Finch and Davis [35] suggested that "uranium could have been pre-concentrated during deposition of overbank boggy mudstone", again emphasizing a role for carbonaceous matter as a reductant. Neither of these models include a role for detrital materials, and in the case of Bigrlyi the low organic carbon content would rule out a significant role for carbonaceous matter as a reductant. Yet at Bigrlyi, vanadium and some uranium were pre-concentrated in detrital components. A hybrid model involving an interplay between solution-precipitation processes, detrital transport, pre-concentration and post-depositional U remobilization may be the most appropriate model to accommodate all the features of the Bigrlyi deposit.
The ultimate source of Bigrlyi vanadium is not known, but a black shale source or a hydrothermal source (such as the roscoelite-related alteration identified in some alkaline magmatic systems, e.g., [36]) in the hinterland are possible. Original detrital components, as well as oxidatively leached V (V) from such a source, could have both been contributors to the V source inventory. Uranium was most likely leached from uraniferous granitoids and pegmatites in the basement Arunta terrain, and transported as stable U (VI) carbonate complexes. In fact, polycrystalline quartz-K-feldspar clasts showing radiation-induced damage features, but lacking uranium content are found in the host sandstone, and likely represent uranium-leached detritus from the granitic source region.
The primary reductants (and adsorbents) were likely to be Fe(II)-bearing clay minerals and pyrite; Fe(II) is known to be a thermodynamically favored abiotic reducer of U (VI) [37], and pyrite is an efficient scavenger of soluble U (VI) [38]. There may also be a role for V (III) -bearing minerals as reductants.
The key element selenium is mobile as soluble selenite anions in oxidized waters, and is precipitated under reducing conditions as metal selenides. In roll-front deposits, Se is closely associated with the U orebody due to these elements having similar redox behaviors [39]. In Bigrlyi, elevated Se/S ratios in U-mineralized zones may be the result of pyrite replacement, S mobilization and the liberation of pyrite-derived Se, followed by selenide and U mineral co-precipitation. The proposed mineralization model for Bigrlyi is envisaged as a solution-precipitation process combined with the transport of U-V-mineralized detritus, as follows. The erosion of a hinterland vanadium source results in the introduction of Fe-V-bearing detrital phyllosilicate-rich materials into the fluvial system. Vanadium-bearing minerals located in sites above a fluctuating water table may be oxidatively mobilized and locally re-precipitated by reaction with Fe-rich detrital minerals, in a process of 'vanadization'. Oxidative leaching of uraniferous granites and pegmatites in the hinterland provides the main source of soluble U, which is reduced and precipitated at the same depositional sites as host Fe-V-bearing phyllosilicates. These initial sites (immature fluvial channels) were eroded and scoured out in flood events, and a mix of semi-lithified mineralized sediments were re-deposited downstream. Mineralized detrital materials are eventually concentrated in a reduced aquifer downstream in the fluvial channel system. Uranium is mobilized and reprecipitated as coffinite in various later events that cause U to be redistributed from V-rich domains and from partly oxidized parts of the stratigraphy. The aquifer is isolated, and mineralization is protected by chemical barriers and mechanical compaction. The deposit is affected by tectonism associated with basin shortening and inversion, and the host stratigraphy is exhumed and tilted vertically on the northern margin of the basin, a possible driver for U redistribution early in the deposit history.
The preservation and mineralization potential in the Bigrlyi deposit is controlled by the deposition of detrital phyllosilicates as clasts or matrix clays, calcite cementation driven by evaporation, and fluvial channel architecture [27]. An exploration implication of the mineralization model is the recognition that additional favorable sites in which U-V minerals were concentrated likely existed upstream and downstream of the present depositional site. Unfortunately, at Bigrlyi, uplift on the basin margin has likely eroded upstream portions of the fluvial channel system, but other mineralized sites within the fluvial channel system might be traceable using a combination of palaeocurrent direction indicators and geophysical techniques that detect pyritic sandstone.

Comparison with Colorado Plateau U-V Deposits
The Colorado Plateau U-V deposits have been known for over 120 years, and controls on mineralization were first discussed by Wright [40] and subsequently studied in the late 1950s and early 1990s. The most important tabular style deposits occur within the Jurassic Saltwash Member of the Morrison Formation, and other deposits are known in Triassic and Cretaceous sandstones. As tat Bigrlyi, tabular ore bodies are associated with fluvial channels in arkosic to micaceous sandstones. Unlike at Bigrlyi, mineralization here is often associated with an abundance organic plant matter [41,42]; organic matter-poor deposits are known but are less common. Furthermore, U-V mineralization does not always occur together; for example, the Jurassic Entrada Sandstone is dominated by V mineralization with only minor U [40].
Vanadium minerals in the Colorado Plateau deposits commonly comprise hydrated oxides (e.g., montroseite, haggite, doloresite), vanadates (e.g., carnotite, tyuyamunite, pascoite) and phyllosilicates (roscoelite, V-hydrous mica, V-chlorite) and other rare minerals [10,41,43]. Haggite and montroseite are common in unoxidized ore, while doloresite is an indicator for the onset of oxidation. In the Saltwash-type deposits, the highest concentrations of V occur together with U mineralization, but V is distributed over a larger volume, and as with the Bigrlyi deposit, forms a halo around the deposit [13]. Other element distributions related to mineralization include the co-enrichment of S, Se, Mo and Cr in the mineralized zone [13], and the formation of authigenic calcite and dolomite stratigraphically above mineralization [7].
The genesis of Saltwash-type tabular U-V deposits is highly debated, particularly with regard to the relative timing of mineralization, the nature of diagenetic processes, the inferred fluid source(s), and the role of reductants including organic matter [4,12,13,44]. Northrop and Goldhaber [13] argue that mineralization formed at the interface between a relatively stagnant, evaporite-derived basinal brine and flowing meteoric water. On the other hand, Hansley and Spirakis [12] argue that the abundance and distribution of amorphous organic matter, in combination with diagenetic alteration processes of variable intensity, are crucial in the genesis of Saltwash-type deposits. Other works that invoke the importance of organic matter include Northrop and Goldhaber [13] and Wanty et al. [7].
The published work on the Colorado Plateau U-V deposits in general shows that these deposits have similar characteristics to those described for Bigrlyi, including ore mineralogy, vanadium distribution and abundance of carbonate cements engulfing ore zones, and the distribution of kaolinite [13].
However, there are differences, such as the lack of organic matter and sulfur depletion rather than enrichment; the timing and source of V mineralization is also interpreted differently. Northrop and Goldhaber [13] suggested that hypersaline brine played an important role in dolomite precipitation and ore formation. By contrast, the Bigrlyi carbonate mineral is predominantly evaporation-related calcite.
Uranium mineralization at Bigrlyi is strongly associated with detrital Fe-V-clays, which likely acted as reductants, but this has not been demonstrated in the Colorado Plateau deposits. Although V-bearing chlorite and chlorite-smectite (corrensite) are present in the Colorado Plateau deposits, they only appear to occur in the mineralized intervals, where they were interpreted to be authigenic phases formed as a result of U mineralization [13].
Based mainly on the Colorado Plateau deposits, Finch and Davis [35] outlined nine diagnostic characteristics for Sandstone-hosted uranium deposits, including the Saltwash type. As considered below, it is instructive to compare how these characteristics might or might not apply to the Bigrlyi deposit. In the future, these characteristics may need to be modified to encompass a wider geographical range of deposits, and additional investigations of the Colorado Plateau deposits may be necessary to better clarify particular criteria. The lack of modern published petrographic/core imagery studies and geochemical data hinders any further detailed comparison with the Colorado Plateau deposits (Table 4). Table 4. Comparison Bigrlyi Deposit-Finch and Davis Sandstone-hosted U-V Deposit Criteria.

Criteria Bigrlyi
Host-rock unit is Silurian or younger in age.

Yes, Carboniferous
Host-rock unit was deposited most commonly in fluvial and lacustrine environments in continental settings. Yes Provenance for sedimentation was commonly granitic or felsic volcanic terranes, either of which provided a plausible source for uranium.

Yes, granitic
Fossil carbonized plant matter or humic material is commonly present.
Only plant fossil imprints are present at Bigrlyi, organic matter preservation was poor, and concentrations of organic carbon are low.
Uranium concentrations were controlled by sedimentary features rather than tectonic fracture structures. Yes, fluvial system architecture controls were important at Bigrlyi The host-rock units are those with good regional transmissivity, and the deposits are localized where sandstone/mudstone ratios are near 1:1.
Good transmissivity aquifer characteristics were presumably true for the pre-diagenetic sediments deposited at Bigrlyi, and bounding shale units would have been important in confining the original aquifer in grey sandstone. Although mudstone lithic clasts are important sites for mineralization, overall, mudstone is not as abundant as sandstone in mineralized domains. Loss of permeability and porosity occurred due to early diagenetic carbonate cementation.
Low-temperature groundwaters were the mineralizing solutions, rather than high-temperature hypogene fluids.
Yes, but earlier-formed detrital components are also present, some could have been of original hydrothermal origin.
The ore minerals are epigenetic even though mineralization was commonly part of diagenesis.
No, mineralization was largely pre-diagenetic and there were detrital components.
Mineralization took place in rocks having original low-angle basin-ward dips, and in most cases the deposits were preserved because of only slight increases in regional dip. More severe tectonic events have redistributed ore in some places.
Yes and no, vertical tilting of the host sediments and deformation at the basin margins did not appear to affect the preservation of ore domains at Bigrlyi, but may have initiated and been a driver for U redistribution.

Conclusions
In this contribution, the key features of the petrography, mineralogy and geochemistry of the Bigrlyi U-V sandstone-hosted deposit have been presented and discussed. Although the Bigrlyi deposit shares many of the features that characterize Saltwash-type tabular sandstone U-V deposits, as exemplified by those from the Colorado Plateau, there are a number of differences. In particular, the recognition of a role for detrital components in hosting U-V mineralization has led to the development of a new mineralization model. Following the multiple migration-accretion hypothesis of Gruner [34], the proposed mineralization model is envisaged as involving solution-precipitation processes, detrital transport, and U-V mineral preconcentration. It is hoped this work will encourage the further study and comparison of tabular, sandstone U-V deposits from around the world.
Author Contributions: All authors contributed to conceptualization, writing, visualization. All authors have read and agreed to the published version of the manuscript.
Funding: This research received no external funding.