The Coles Hill Uranium Deposit, Virginia, USA: Geology, Geochemistry, Geochronology, and Genetic Model

The Coles Hill uranium deposit, with an indicated resource of about 130 Mlb of U3O8, is the largest unmined uranium deposit in the United States. The deposit is hosted in the Taconian (approx. 480–450 Ma) Martinsville igneous complex, which consists of the Ordovician Leatherwood Granite (granodiorite) and the Silurian Rich Acres Formation (diorite). The host rock was metamorphosed to orthogneiss during the Alleghanian orogeny (approx. 325–260 Ma), when it also underwent dextral strike-slip movement along the Brookneal shear zone. During the Triassic, extensional tectonics led to the development of the Dan River Basin that lies east of Coles Hill. The mineralized zone is hosted in brittle structures in the footwall of the Triassic Chatham fault that forms the western edge of the basin. Within brittle fracture zones, uranium silicate and uranium-bearing fluorapatite with traces of brannerite form veins and breccia-fill with chlorite, quartz, titanium oxide, pyrite, and calcite. Uranium silicates also coat and replace primary titanite, zircon, ilmenite, and sulfides. Sodium metasomatism preceded and accompanied uranium mineralization, pervasively altering host rock and forming albite from primary feldspar, depositing limpid albite rims on igneous feldspar, altering titanite to titanium oxide and calcite, and forming riebeckite. Various geothermometers indicate temperatures of less than ~200°C during mineralization. In situ U-Pb analyses of titanite, Ti-oxide, and apatite, along with Rb/Sr and U/Pb isotope systematics of whole-rock samples, resolve the timing of geologic processes affecting Coles Hill. The host Leatherwood Granite containing primary euhedral titanite is dated at 450 to 445 Ma, in agreement with previously obtained ages from zircon in the Martinsville igneous complex. A regional metamorphic event at 330 to 310 Ma formed anhedral titanite and some apatite, reequilibrated whole-rock Rb/Sr and U-Pb isotopes, and is interpreted to have coincided with movement along the Brookneal shear zone. During shearing and metamorphism, primary refractory uranium-bearing minerals including titanite, zircon, and uranothorite were recrystallized, and uranium was liberated and mixed locally with hematite, clay, and other fine-grained minerals. Uranium mineralization was accompanied by a metasomatic episode between 250 and 200 Ma that reset the Rb-Sr and U-Pb isotope systems and formed titanium oxide and apatite that are associated and, in places, intimately intergrown with uranium silicate dating mineralization. This event coincides with rifting that formed the Dan River Basin and was a precursor to the breakup of Pangea. The orientation of late-stage tectonic stylolites is compatible with their formation during Late Triassic to Early Jurassic basin inversion, postdating the main stage of uranium mineralization and effectively dating mineralization as Mesozoic. Based on the close spatial and temporal association of uranium with apatite, we propose that uranium was carried as a uranyl-phosphate complex. Uranium was locally reduced by coupled redox reactions with ferrous iron and sulfide minerals in the host rock, forming uranium silicates. The release of calcium during sodium metasomatic alteration of primary calcic feldspar and titanite in the host rock initiated successive reactions in which uranium and phosphate in mineralizing fluids combined with calcium to form U-enriched fluorapatite. Based on the deposit mineralogy, oxygen isotope geochemistry, and trace element characteristics of uranium silicate and gangue minerals, the primary mineralizing fluids likely included connate and/or meteoric water sourced from the adjacent Dan River Basin. High heat flow related to Mesozoic rifting may have driven these (P-Na-F-rich) fluids through local aquifers and into basin margin faults, transporting uranium from the basin or mobilizing uranium from previously formed U minerals in the Brookneal shear zone, or from U-enriched older basement rock.

Supplemental Material: A robust age model for the Cryogenian Pocatello Formation of southeastern Idaho (northwestern USA) from tandem in situ and isotope dilution U-Pb dating of volcanic tuffs and epiclastic detrital zircons

<div>Figure S1. Thin-section photomicrographs for lithologies of the Bannock Volcanic Member and Scout Mountain Member of the Pocatello Formation exposed at Scout Mountain, Idaho. Sample numbers are illustrated on the stratigraphic section of Figure 4. Field of view is 24 mm × 40 mm. Tables S1–S3: LA-ICPMS U-Pb isotope and trace element concentration data. Table S4: CA-IDTIMS U-Pb isotope data.<br></div>

Supplemental Material: A robust age model for the Cryogenian Pocatello Formation of southeastern Idaho (northwestern USA) from tandem in situ and isotope dilution U-Pb dating of volcanic tuffs and epiclastic detrital zircons

<div>Figure S1. Thin-section photomicrographs for lithologies of the Bannock Volcanic Member and Scout Mountain Member of the Pocatello Formation exposed at Scout Mountain, Idaho. Sample numbers are illustrated on the stratigraphic section of Figure 4. Field of view is 24 mm × 40 mm. Tables S1–S3: LA-ICPMS U-Pb isotope and trace element concentration data. Table S4: CA-IDTIMS U-Pb isotope data.<br></div>

“Non-invasive” portable laser ablation sampling for lead isotope analysis of archaeological silver: a comparison with bulk and in situ laser ablation techniques

Portable laser ablation is tested for Pb isotope analysis of silver, yielding results comparable to bulk analyses. This low impact method removes the need for invasive sampling, opening the way for analysis of archaeological silver worldwide.

Investigating influences on the Pb pseudo-isochron using three-dimensional mantle convection models with a continental reservoir

For mid-ocean ridge basalts (MORBs) and ocean island basalts (OIBs), measurements of Pb isotope ratios show broad linear correlations with a certain degree of scatter. In 207Pb/204Pb - 206Pb/204Pb space, the best fit line defines a pseudo-isochron age (τPb) of ~1.9 Gyr.Previous modelling suggests a relative change in the behaviours of U and Pb between 2.25-2.5 Ga, resulting in net recycling of HIMU (high U/Pb) material in the latter part of Earth's history, to explain the observed τPb. However, simulations in which fractionation is controlled by a single set of partition coefficients throughout the model runs fail to reproduce τPb and the observed scatter in Pb isotope ratios. We build on these models with 3D mantle convection simulations including parameterisations for melting, U recycling from the continents and preferential removal of Pb from subducted oceanic crust.We find that both U recycling after the great oxygenation event (GOE) and Pb extraction after the onset of plate tectonics, are required in order to fit the observed gradient and scatter of both the 207Pb/204Pb - 206Pb/204Pb and 208Pb/204Pb - 206Pb/204Pb arrays. Unlike much previous work, our model does not require accumulations of subducted oceanic crust to persist at the CMB for long periods of time in order to match geochemical observations.

The volcanic and magmatic evolution of Tongariro volcano, New Zealand

<p>Detailed mapping studies of Quaternary stratovolcanoes provide critical frameworks for examining the long-term evolution of magmatic systems and volcanic behaviour. For stratovolcanoes that have experienced glaciation, edifice-forming products also act as climate-proxies from which ice thicknesses can be inferred at specific points in time. One such volcano is Tongariro, which is located in the southern Taupō Volcanic Zone of New Zealand’s North Island. This study presents the results of new detailed mapping, geochronological and geochemical investigations on edifice-forming materials to reconstruct Tongariro’s volcanic and magmatic history which address the following questions: (1) Does ice coverage on stratovolcanoes influence eruptive rates and behaviour (or record completeness)? (2) What is the relationship between magmatism, its expression (i.e. volcanism) and external but related processes such as tectonics? (3) How are intermediate-composition magmas assembled and what controls their diversity? (4) What are the relative proportions of mantle-derived and crust-derived materials in intermediate composition arc magmas? (5) Do genetic relationships exist between andesite and rhyolite magmas in arc settings? Samples from 250 new field localities in under-examined areas of Tongariro were analysed for major oxide, trace element and Sr-Nd-Pb isotope compositions. Analyses were performed on whole-rock, groundmass and xenolith samples. The stratigraphic framework for these geochemical data was established from field observations and 29 new 40Ar/39Ar age determinations, which were synthesised with volume estimates and petrographic observations for all Tongariro map units. Mapping results divide Tongariro into 36 distinct map units (at their greatest level of subdivision) which were organised into formations and constituent members. New 40Ar/39Ar age determinations reveal continuous eruptive activity at Tongariro from at least 230 ka to present, including during glacial periods. This adds to the discovery of an inlier of old basaltic-andesite (512 ± 59 ka) on Tongariro’s NW sector that has an unclear source vent. Hornblende-phyric andesite boulders, mapped into the Tupuna Formation (new), yield the oldest 40Ar/39Ar age determination (304 ± 11 ka) for materials confidently attributed to Tongariro. Tupuna Formation andesites are correlated with Turakina Formation debris flows that were deposited between 349 to 309 ka in the Wanganui Basin, ~100 km south of Tongariro, which indicates that Ruapehu did not exist at this time, at least not in its current form. Tongariro has a total edifice volume of ~90 km3, 19 km3 of which is represented by exposed mapped units. The total ringplain volume immediately adjacent to Tongariro contains ~60 km3 of material. The volume of exposed glacial deposits are no more than 1 km3. During periods of major ice coverage, edifice-building rates on Tongariro were only 17-21 % of edifice-building rates during warmer climatic periods. Because shifts in edifice-building rates do not coincide with changes in erupted compositions, differences in edifice-building rates reflect a preservation bias. Materials erupted during glacial periods were emplaced onto ice masses and conveyed to the ringplain as debris, which explains reduced preservation rates at these times. MgO concentrations in Tongariro stratigraphic units with ages between 230 and 0 ka display successive and irregular cyclicity that occurs over ~10-70 kyr intervals, which reflect episodes of enhanced mafic magma replenishment. During these cycles, more rapid (≤10 kyr) increases in MgO concentrations to ≥5-9 wt% are followed by gradual declines to ~2-5 wt%, with maxima at ~230, ~160, ~117, ~88, ~56, ~35, ~17.5 ka and during the Holocene. Contemporaneous variations in Tongariro and Ruapehu magma compositions (e.g. MgO, Rb/Sr, Sr-Nd-Pb isotope ratios) for the 200-0 ka period coincide with reported zircon growth model-ages in Taupō magmas. This contemporaneity reflects regional tectonic processes that have externally regulated and synchronised the timings of elevated mafic replenishment episodes versus periods of prolonged crustal residence at each of these volcanoes. Isotopic Sr-Nd-Pb data from metasedimentary xenoliths, groundmass separates and whole-rock samples indicate that two or three separate metasedimentary terranes (in the upper 15 km of the crust) were assimilated into Tongariro magmas. These are the Kaweka terrane and the Waipapa or Pahau terranes (or both). Subhorizontal juxtapositioning of these terranes is indicated by the coexistence of multiple terranes in the same eruptive units. Paired whole-rock and groundmass (interstitial melt) samples have effectively equal Sr-Nd-Pb isotope ratios for the complete range of Tongariro compositions. Despite intra-crystal isotopic heterogeneities that are likely widespread, the new data show that crystal fractionation and assimilation occur in approximately equal balance for essentially all Tongariro eruptives. Assimilated country rock accounts for 22-31 wt% of the average Tongariro magma. Initial evolution from a Kakuki basalt-type to a Tongariro Te Rongo Member basaltic-andesite reflects the addition of 17 % assimilated metasedimentary basement with a mass assimilation rate/mass crystal fractionation rate ratio—a.k.a. ‘r value’ of 1.8-3.5. Subsequent evolution from a Te Rongo Member basaltic-andesite to other Tongariro eruptive compositions represents 5-14 % more assimilated crust (r values of ~0.1-1.0). Magma evolution from high (>1) to lower (0.1-1.0) r values can explain the dearth of andesitic melt inclusions in (bulk) andesite magmas observed globally. High relative assimilation rates characterise rapid evolution from basalt to basaltic-andesite bulk compositions which contain andesitic interstitial melts. Thus, andesitic melt inclusions have a reduced chance of being preserved in crystals which can explain their low representation in global datasets.</p>