Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Rock Products”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,423 records · Page 79Linked to original sources

Paleozoic magmatism and porphyry Cu-mineralization in an evolving tectonic setting in the North Qilian Orogenic Belt, NW China

The NWW-striking North Qilian Orogenic Belt records the Paleozoic accretion–collision processes in NW China, and hosts Paleozoic Cu–Pb–Zn mineralization that was temporally and spatially related to the closure of the Paleo Qilian-Qinling Ocean. The Wangdian Cu deposit is located in the eastern part of the North Qilian Orogenic Belt, NW China. Copper mineralization is spatially associated with an altered early Paleozoic porphyritic granodiorite, which intruded tonalites and volcaniclastic rocks. Alteration zones surrounding the mineralization progress outward from a potassic to a feldspar-destructive phyllic assemblage. Mineralization consists mainly of quartz-sulfide stockworks and disseminated sulfides, with ore minerals chalcopyrite, pyrite, molybdenite, and minor galena and sphalerite. Gangue minerals include quartz, orthoclase, biotite, sericite, and K-feldspar. Zircon LA-ICPMS U–Pb dating of the ore-bearing porphyritic granodiorite yielded a mean 206 Pb/ 238 U age of 444.6 ± 7.8 Ma, with a group of inherited zircons yielding a mean U–Pb age of 485 ± 12 Ma, consistent with the emplacement age (485.3 ± 6.2 Ma) of the barren precursor tonalite. Rhenium and osmium analyses of molybdenite grains returned model ages of 442.9 ± 6.8 Ma and 443.3 ± 6.2 Ma, indicating mineralization was coeval with the emplacement of the host porphyritic granodiorite. Rhenium concentrations in molybdenite (208.9–213.2 ppm) suggest a mantle Re source. The tonalities are medium-K calc-alkaline. They are characterized by enrichment of light rare-earth elements (LREEs) and large-ion lithophile elements (LILEs), depletion of heavy rare-earth elements (HREEs) and high-field-strength elements (HFSEs), and minor negative Eu anomalies. They have ε Hf ( t ) values in the range of +3.6 to +11.1, with two-stage Hf model ages of 0.67–1.13 Ga, suggesting that the ca. 485 Ma barren tonalites were products of arc magmatism incorporating melts from the mantle wedge and the lithosphere. In contrast, the 40-m.y.-younger ore-bearing porphyritic granodiorite is sub-alkaline and peraluminous. They are enriched in LREEs and LILEs, depleted in HFSEs, and show weak negative Eu anomalies. They display ε Hf ( t ) values of captured or inherited zircons in the range of +8.5 to +10.0, and younger two-stage Hf model ages of 0.78 Ga and 0.86 Ga, similar to those of ca. 485 Ma tonalite. The ca. 445 Ma zircons have ε Hf ( t ) values of −2.1 to +9.9, with two-stage Hf model ages of 0.75–1.27 Ga. Moreover, they have relatively high oxygen fugacity than that of the precursor barren tonalite. The ca. 445 Ma magmas at Wangdian thus formed in a subduction setting, and incorporated melts from the subduction-modified lithosphere that had previously been enriched by additions of chalcophile and siderophile element-rich materials by the earlier magmatism and metasomatism during the Paleo Qilian-Qinling Ocean subduction event.

North Qilian Orogenic Belt↗

Coupled arsenotrophy in a hot spring photosynthetic biofilm at Mono Lake, California

Red-pigmented biofilms grow on rock and cobble surfaces present in anoxic hot springs located on Paoha Island in Mono Lake. The bacterial community was dominated (∼ 85% of 16S rRNA gene clones) by sequences from the photosynthetic Ectothiorhodospira genus. Scraped biofilm materials incubated under anoxic conditions rapidly oxidized As(III) to As(V) in the light via anoxygenic photosynthesis but could also readily reduce As(V) to As(III) in the dark at comparable rates. Back-labeling experiments with 73 As(V) demonstrated that reduction to 73 As(III) also occurred in the light, thereby illustrating the cooccurrence of these two anaerobic processes as an example of closely coupled arsenotrophy. Oxic biofilms also oxidized As(III) to As(V). Biofilms incubated with [ 14 C]acetate oxidized the radiolabel to 14 CO 2 in the light but not the dark, indicating a capacity for photoheterotrophy but not chemoheterotrophy. Anoxic, dark-incubated samples demonstrated As(V) reduction linked to additions of hydrogen or sulfide but not acetate. Chemoautotrophy linked to As(V) as measured by dark fixation of [ 14 C]bicarbonate into cell material was stimulated by either H 2 or HS − . Functional genes for the arsenate respiratory reductase ( arrA ) and arsenic resistance ( arsB ) were detected in sequenced amplicons of extracted DNA, with about half of the arrA sequences closely related (∼98% translated amino acid identity) to those from the family Ectothiorhodospiraceae . Surprisingly, no authentic PCR products for arsenite oxidase ( aoxB ) were obtained, despite observing aerobic arsenite oxidation activity. Collectively, these results demonstrate close linkages of these arsenic redox processes occurring within these biofilms. Oxyanions of the group 15 element arsenic, arsenate [As(V)] and arsenite [As(III)], have been known for millennia to be potent poisons. Despite its well-established toxicity to life, the phenomenon of arsenic resistance was discovered whereby some microorganisms maintain an otherwise “normal” existence in the presence of high concentrations of As(V) or As(III) (17, 29, 31). More recently it has become recognized that certain representatives from the bacterial and archaeal domains can actually exploit the electrochemical potential of the As(V)/As(III) redox couple (+130 mV) to gain energy for growth. This can be achieved either by employing As(III) as an autotrophic electron donor or by using As(V) as a respiratory electron acceptor (18, 21, 34). The latter phenomenon, although most commonly associated with chemoheterotrophy, can also employ inorganic substances like sulfide or H 2 . Indeed, As(V)-respiring anaerobes displaying a capacity for chemoautotrophy with these electron donors have been isolated and described (5, 7, 16). We recently reported that photoautotrophy is supported by As(III) in anoxic biofilms located in hot springs on Paoha Island in Mono Lake, CA (15). This process represented a novel means of As(III) oxidation achieved via anoxygenic photosynthesis occurring in certain photosynthetic bacteria (i.e., Ectothiorhodospira ) and possibly within some cyanobacteria as well (e.g., “ Oscillatoria ”). Whether or not a microbial habitat is overtly oxic or anoxic, or temporally shifts between these two states over a diel cycle, critical energy linkages between aerobes and anaerobes have long been known for the biogeochemical cycles of key elements, such as sulfur, iron, and nitrogen. Most prominently studied is the case of nitrogen, whereby an ecological coupling exists between the processes of nitrification and denitrification (9, 10, 28). The former process provides energy to aerobic nitrifiers, while the latter process consumes the nitrate produced by this reaction, thereby meeting the energy needs of the denitrifiers. For arsenic, the detection of both As(III) oxidation and As(V) reduction in oxic and anoxic incubations of freshly collected periphyton suggested that an analogous coupled process may also occur for this element (12). Similarly, several uncontaminated soils in Japan displayed a capacity for either As(V) reduction or As(III) oxidation upon arsenic oxyanion amendment and whether they were incubated under oxic or anoxic conditions (39). A defined coculture consisting of an aerobic As(III) oxidizer (strain OL1) and an anaerobic As(V) respirer (strain Y5) was shown to function in this fashion under manipulated laboratory conditions of oxygen tension (26). We pursued the phenomenon of coupled arsenic metabolism further by using materials collected from the hot spring biofilms in Mono Lake, but we focused on examination of the cycling of arsenic under anoxic conditions. In this paper we report results obtained by manipulated incubations of red-pigmented biofilms found in the hot springs of Paoha Island. Preliminary community characterizations of these biofilms show that they are dominated by B acteria from the genus Ectothiorhodospira but also harbor an assemblage of Archaea related to the Halobacteriacaea . Incubation results have demonstrated the presence of the following arsenic metabolic activities: respiratory As(V) reduction, photosynthetic anaerobic As(III) oxidation, and aerobic As(III) oxidation, along with the ecophysiological conditions under which they occur. Surprisingly, we were unable to obtain authentic PCR products for arsenite oxidase genes ( aoxB ), despite observing aerobic As(III) oxidation activity. These biofilms serve as a model system for how anaerobic cycling of arsenic can be sustained with oxidation of As(III) by anoxygenic photosynthesis coupled to regeneration of this electron donor via dissimilatory As(V) reduction. The significance that such a light-driven anaerobic ecosystem may have played in the Archean Earth is discussed.

California↗

Cierco Pb-Zn-Ag vein deposits: Isotopic and fluid inclusion evidence for formation during the mesozoic extension in the pyrenees of Spain

The Cierco Pb-Zn-Ag vein deposits, located in the central Pyrenees of Spain, crosscut Paleozoic metasedimentary rocks and are in close proximity to Hercynian granodiorite dikes and plutons. Galena and sphalerite in the deposits have average delta 34 S values of - 4.3 and - 0.8 per mil (CDT), respectively. Coexisting mineral pairs give an isotopic equilibration temperature range of 89 degrees to 163 degrees C which overlaps with the 112 degrees to 198 degrees C range obtained from primary fluid inclusions. Coexisting quartz has a delta 18 O value of 19 + or - 1 per mil (VSMOW). The fluid which deposited these minerals is inferred to have had delta 18 O (sub H 2 ) O and delta 34 S (sub H 2 ) S values of 5 + or - 1 and -1 + or - 1 per mil, respectively. Chemical and microthermometric analyses of fluid inclusions in quartz and sphalerite indicate salinities of 3 to 29 wt percent NaCl equiv with Na (super +) and Ca (super 2+) as the dominant cations in solution. The Br/Cl and I/Cl ratios differ from those characteristic of magmatic waters and pristine seawater, but show some similarity to those observed in deep ground waters in crystalline terranes, basinal brines, and evaporated seawater. Barite, which postdates the sulfides, spans isotopic ranges of 13 to 21 per mil, 10 to 15 per mil, and 0.7109 to 0.7123 for delta 34 S, delta 18 O, and 87 Sr/ 86 Sr, respectively. The three parameters are correlated providing strong evidence that the barites are products of fluid mixing.We propose that the Cierco deposits formed along an extensional fault system at the margin of a marine basin during the breakup of Pangea at some time between the Early Triassic and Early Cretaceous. Sulfide deposition corresponded to an upwelling of hydrothermal fluid from the Paleozoic basement and was limited by the amount of metals carried by the fluid. Barite deposition corresponded to the waning of upward flow and the collapse of sulfate-rich surface waters onto the retreating hydrothermal plume. Calcite precipitated late in the paragenesis as meteoric or marine waters descended into the fault system, possibly during a regression in the overlying basin.There are other deposits resembling Cierco elsewhere in the Iberian peninsula. Taken as a group, they are evidence that hydrothermal circulation systems were widespread during Mesozoic extension. Differences among the deposits can be related to the fact that H 2 S and other solutes had local and variable sources.

Economic Geology↗

Aeolian dust in Colorado Plateau soils: Nutrient inputs and recent change in source

Aeolian dust (windblown silt and clay) is an important component in arid-land ecosystems because it may contribute to soil formation and furnish essential nutrients. Few geologic surfaces, however, have been characterized with respect to dust-accumulation history and resultant nutrient enrichment. We have developed a combination of methods to identify the presence of aeolian dust in arid regions and to evaluate the roles of this dust in ecosystem processes. Unconsolidated sandy sediment on isolated surfaces in the Canyonlands region of the Colorado Plateau differs greatly in mineralogical and chemical composition from associated bedrock, mainly aeolian sandstone. Detrital magnetite in the surficial deposits produces moderately high values of magnetic susceptibility, but magnetite is absent in nearby bedrock. A component of the surficial deposits must be aeolian to account for the abundance of magnetite, which formed originally in far-distant igneous rocks. Particle-size analysis suggests that the aeolian dust component is typically as much as 20–30%. Dust inputs have enriched the sediments in many elements, including P, Mg, Na, K, and Mo, as well as Ca, at sites where bedrock lacks calcite cement. Soil-surface biologic crusts are effective dust traps that apparently record a change in dust sources over the past several decades. Some of the recently fallen dust may result from human disturbance of land surfaces that are far from the Canyonlands, such as the Mojave Desert. Some land-use practices in the study area have the potential to deplete soil fertility by means of wind-erosion removal of aeolian silt. Many studies have addressed the presence of aeolian dust in soils and surficial deposits in deserts to provide important geologic and ecologic information bearing on landscape dynamics (1–16). From this body of work, we have improved understanding about: ( i ) current and past sources and flux of dust, hence changing conditions of dust emission; ( ii ) the genesis of desert soils; ( iii ) the influences of aeolian silt and clay on water-infiltration rates in soil; ( iv ) the evolution of desert surfaces (such as desert pavement) relevant to surface stability, as well as the distribution of surface and subsurface water; and ( v ) interrelations among aeolian dust, distribution of plants and soil crust, rain-water runoff, and productivity. Nevertheless, we lack fundamental knowledge about the accumulation history of aeolian dust on most landscapes and about how to discriminate between contributions from parent material and aeolian dust to the biotic system. Although many different methods provide clear evidence for aeolian input into soils (4, 6, 8, 9, 11, 13, 17, 18), ecosystem studies would benefit from rapid assessment of aeolian components in soils. Here we demonstrate that magnetic methods, which characterize the type and distribution of iron oxide minerals, can be applied to detect the presence of aeolian dust in young surficial sediments and soils over large arid-land areas. The magnetic results, combined with geochemical and textural analyses, form a basis for understanding the influence of fine-grained aeolian inputs on soil fertility of the central Colorado Plateau, Utah.

Utah↗

Petroleum system modeling capabilities for use in oil and gas resource assessments

Summary: Petroleum resource assessments are among the most highly visible and frequently cited scientific products of the U.S. Geological Survey. The assessments integrate diverse and extensive information on the geologic, geochemical, and petroleum production histories of provinces and regions of the United States and the World. Petroleum systems modeling incorporates these geoscience data in ways that strengthen the assessment process and results are presented visually and numerically. The purpose of this report is to outline the requirements, advantages, and limitations of one-dimensional (1-D), two-dimensional (2-D), and three-dimensional (3-D) petroleum systems modeling that can be applied to the assessment of oil and gas resources. Primary focus is on the application of the Integrated Exploration Systems (IES) PetroMod? software because of familiarity with that program as well as the emphasis by the USGS Energy Program on standardizing to one modeling application. The Western Canada Sedimentary Basin (WCSB) is used to demonstrate the use of the PetroMod? software. Petroleum systems modeling quantitatively extends the 'total petroleum systems' (TPS) concept (Magoon and Dow, 1994; Magoon and Schmoker, 2000) that is employed in USGS resource assessments. Modeling allows integration of state-of-the-art analysis techniques, and provides the means to test and refine understanding of oil and gas generation, migration, and accumulation. Results of modeling are presented visually, numerically, and statistically, which enhances interpretation of the processes that affect TPSs through time. Modeling also provides a framework for the input and processing of many kinds of data essential in resource assessment, including (1) petroleum system elements such as reservoir, seal, and source rock intervals; (2) timing of depositional, hiatus, and erosional events and their influences on petroleum systems; (3) incorporation of vertical and lateral distribution and lithologies of strata that compose the petroleum systems; and (4) calculations of pressure-volume-temperature (PVT) histories. As digital data on petroleum systems continue to expand, the models can integrate these data into USGS resource assessments by building and displaying, through time, areas of petroleum generation, migration pathways, accumulations, and relative contributions of source rocks to the hydrocarbon components. IES PetroMod? 1-D, 2-D, and 3-D models are integrated such that each uses the same variables for petroleum systems modeling. 1-D burial history models are point locations, mainly wells. Maps and cross-sections model geologic information in two dimensions and can incorporate direct input of 2-D seismic data and interpretations using various formats. Both 1-D and 2-D models use data essential for assessments and, following data compilation, they can be completed in hours and retested in minutes. Such models should be built early in the geologic assessment process, inasmuch as they incorporate the petroleum system elements of reservoir, source, and seal rock intervals with associated lithologies and depositional and erosional ages. The models can be used to delineate the petroleum systems. A number of 1-D and 2-D models can be constructed across a geologic province and used by the assessment geologists as a 3-D framework of processes that control petroleum generation, migration, and accumulation. The primary limitation of these models is that they only represent generation, migration, and accumulation in two dimensions. 3-D models are generally built at reservoir to basin scales. They provide a much more detailed and realistic representation of petroleum systems than 1-D or 2-D models because they portray more fully the temporal and physical relations among (1) burial history; (2) lithologies and associated changes through burial in porosity, permeability, and compaction; (3) hydrodynamic effects; and (4) other parameters that influence petroleum gen

Open-File Report↗

Geophysical Studies Based on Gravity and Seismic Data of Tule Desert, Meadow Valley Wash, and California Wash Basins, Southern Nevada

Gravity and seismic data from Tule Desert, Meadow Valley Wash, and California Wash, Nevada, provide insight into the subsurface geometry of these three basins that lie adjacent to rapidly developing areas of Clark County, Nevada. Each of the basins is the product of Tertiary extension accommodated with the general form of north-south oriented, asymmetrically-faulted half-grabens. Geophysical inversion of gravity observations indicates that Tule Desert and Meadow Valley Wash basins are segmented into subbasins by shallow, buried basement highs. In this study, basement refers to pre-Cenozoic bedrock units that underlie basins filled with Cenozoic sedimentary and volcanic units. In Tule Desert, a small, buried basement high inferred from gravity data appears to be a horst whose placement is consistent with seismic reflection and magnetotelluric observations. Meadow Valley Wash consists of three subbasins separated by basement highs at structural zones that accommodated different styles of extension of the adjacent subbasins, an interpretation consistent with geologic mapping of fault traces oblique to the predominant north-south fault orientation of Tertiary extension in this area. California Wash is a single structural basin. The three seismic reflection lines analyzed in this study image the sedimentary basin fill, and they allow identification of faults that offset basin deposits and underlying basement. The degree of faulting and folding of the basin-fill deposits increases with depth. Pre-Cenozoic units are observed in some of the seismic reflection lines, but their reflections are generally of poor quality or are absent. Factors that degrade seismic reflector quality in this area are rough land topography due to erosion, deformed sedimentary units at the land surface, rock layers that dip out of the plane of the seismic profile, and the presence of volcanic units that obscure underlying reflectors. Geophysical methods illustrate that basin geometry is more complicated than would be inferred from extrapolation of surface topography and geology, and these methods aid in defining a three-dimensional framework to understand groundwater storage and flow in southern Nevada.

Open-File Report↗

Eruptive stratigraphy of the Tatara-San Pedro complex, 36°S, sourthern volcanic zone, Chilean Andes: Reconstruction method and implications for magma evolution at long-lived arc volcanic centers

The Quaternary Tatara–San Pedro volcanic complex (36°S, Chilean Andes) comprises eight or more unconformity-bound volcanic sequences, representing variably preserved erosional remnants of volcanic centers generated during ∼930 ky of activity. The internal eruptive histories of several dominantly mafic to intermediate sequences have been reconstructed, on the basis of correlations of whole-rock major and trace element chemistry of flows between multiple sampled sections, but with critical contributions from photogrammetric, geochronologic, and paleomagnetic data. Many groups of flows representing discrete eruptive events define internal variation trends that reflect extrusion of heterogeneous or rapidly evolving magma batches from conduit–reservoir systems in which open-system processes typically played a large role. Long-term progressive evolution trends are extremely rare and the magma compositions of successive eruptive events rarely lie on precisely the same differentiation trend, even where they have evolved from similar parent magmas by similar processes. These observations are not consistent with magma differentiation in large long-lived reservoirs, but they may be accommodated by diverse interactions between newly arrived magma inputs and multiple resident pockets of evolved magma and/or crystal mush residing in conduit-dominated subvolcanic reservoirs. Without constraints provided by the reconstructed stratigraphic relations, the framework for petrologic modeling would be far different. A well-established eruptive stratigraphy may provide independent constraints on the petrologic processes involved in magma evolution—simply on the basis of the specific order in which diverse, broadly cogenetic magmas have been erupted. The Tatara–San Pedro complex includes lavas ranging from primitive basalt to high-SiO 2 rhyolite, and although the dominant erupted magma type was basaltic andesite (∼52–55 wt % SiO 2 ) each sequence is characterized by unique proportions of mafic, intermediate, and silicic eruptive products. Intermediate lava compositions also record different evolution paths, both within and between sequences. No systematic long-term pattern is evident from comparisons at the level of sequences. The considerable diversity of mafic and evolved magmas of the Tatara–San Pedro complex bears on interpretations of regional geochemical trends. The variable role of open-system processes in shaping the compositions of evolved Tatara–San Pedro complex magmas, and even some basaltic magmas, leads to the conclusion that addressing problems such as arc magma genesis and elemental fluxes through subduction zones on the basis of averaged or regressed reconnaissance geochemical datasets is a tenuous exercise. Such compositional indices are highly instructive for identifying broad regional trends and first-order problems, but they should be used with extreme caution in attempts to quantify processes and magma sources, including crustal components, implicated in these trends.

Andes Mountains, Tatara-San Pedro Volcanic Complex↗

Hydrogeology, simulated ground-water flow, and ground-water quality, Wright-Patterson Air Force Base, Ohio

Ground water is the primary source of water in the Wright-Patterson Air Force Base area. The aquifer consists of glacial sands and gravels that fill a buried bedrock-valley system. Consolidated rocks in the area consist of poorly permeable Ordovician shale of the Richmondian stage, in the upland areas, the Brassfield Limestone of Silurian age. The valleys are filled with glacial sediments of Wisconsinan age consisting of clay-rich tills and coarse-grained outwash deposits. Estimates of hydraulic conductivity of the shales based on results of displacement/recovery tests range from 0.0016 to 12 feet per day; estimates for the glacial sediments range from less than 1 foot per day to more than 1,000 feet per day. Ground water flow from the uplands towards the valleys and the major rivers in the region, the Great Miami and the Mad Rivers. Hydraulic-head data indicate that ground water flows between the bedrock and unconsolidated deposits. Data from a gain/loss study of the Mad River System and hydrographs from nearby wells reveal that the reach of the river next to Wright-Patterson Air Force Base is a ground-water discharge area. A steady-state, three-dimensional ground-water-flow model was developed to simulate ground-water flow in the region. The model contains three layers and encompasses about 100 square miles centered on Wright-Patterson Air Force Base. Ground water enters the modeled area primarily by river leakage and underflow at the model boundary. Ground water exits the modeled area primarily by flow through the valleys at the model boundaries and through production wells. A model sensitivity analysis involving systematic changes in values of hydrologic parameters in the model indicates that the model is most sensitive to decreases in riverbed conductance and vertical conductance between the upper two layers. The analysis also indicates that the contribution of water to the buried-valley aquifer from the bedrock that forms the valley walls is about 2 to 4 percent of the total ground-water flow in the study area. Ground waters in the vicinity of Wright-Patterson Air Force Base can be classified into two compositional groups on the basis of their chemical composition: calcium magnesium bicarbonate-type and sodium chloride-type waters. Calcium magnesium bicarbonate-type waters are found in the glacial deposits and the Brassfield Limestone, whereas the sodium chloride waters are exclusively associated with the shales. Equilibrium speciation calculations indicate that ground water of the glacial drift aquifer is in equilibrium with calcite, dolomite, and chalcedony, but is undersaturated with respect to gypsum and fluorite. Waters from the shales are slightly supersaturated with respect to calcite, dolomite, and siderite but are undersaturated with respect to chalcedony. Simple-mass balance calculations treating boron as a conservative species indicate that little (< 5 percent) or no recharge from the shales to the glacial drift aquifer takes place. Data on the stable isotopes of oxygen and hydrogen indicate a meteoric origin for all ground water beneath Wright-Patterson Air Force Base, but the data were inconclusive with respect to identification of distinct isotopic differences between water collected from the glacial drift and bedrock aquifers. Tritium concentrations used to distinguish waters having a pre-and post-1953 recharge component indicate that most water entered the glacial drift aquifer after 1953. This finding indicates that recharge from shallow to deep parts (greater than 150 feet) of the aquifer takes place over time intervals of a few years or decades. However, the fact that some deep parts of the glacial aquifer did not contain measurable tritium indicates that ground-water flow from recharge zones to these parts of the aquifer takes decades or longer.

Water-Resources Investigations Report↗

Incremental assembly and prolonged consolidation of Cordilleran magma chambers--Evidence from the Southern Rocky Mountain volcanic field

Recent inference that Mesozoic Cordilleran plutons grew incrementally during >10 6 yr intervals, without the presence of voluminous eruptible magma at any stage, minimizes close associations with large ignimbrite calderas. Alternatively, Tertiary ignimbrites in the Rocky Mountains and elsewhere, with volumes of 1&ndash;5 &times; 10 3 km 3 , record multistage histories of magma accumulation, fractionation, and solidification in upper parts of large subvolcanic plutons that were sufficiently liquid to erupt. Individual calderas, up to 75 km across with 2&ndash;5 km subsidence, are direct evidence for shallow magma bodies comparable to the largest granitic plutons. As exemplified by the composite Southern Rocky Mountain volcanic field (here summarized comprehensively for the first time), which is comparable in areal extent, magma composition, eruptive volume, and duration to continental-margin volcanism of the central Andes, nested calderas that erupted compositionally diverse tuffs document deep composite subsidence and rapid evolution in subvolcanic magma bodies. Spacing of Tertiary calderas at distances of tens to hundreds of kilometers is comparable to Mesozoic Cordilleran pluton spacing. Downwind ash in eastern Cordilleran sediments records large-scale explosive volcanism concurrent with Mesozoic batholith growth. Mineral fabrics and gradients indicate unified flow-age of many pluton interiors before complete solidification, and some plutons contain ring dikes or other textural evidence for roof subsidence. Geophysical data show that low-density upper-crustal rocks, inferred to be plutons, are 10 km or more thick beneath many calderas. Most ignimbrites are more evolved than associated plutons; evidence that the subcaldera chambers retained voluminous residua from fractionation. Initial incremental pluton growth in the upper crust was likely recorded by modest eruptions from central volcanoes; preparation for caldera-scale ignimbrite eruption involved recurrent magma input and homogenization high in the chamber. Some eroded calderas expose shallow granites of similar age and composition to tuffs, recording sustained postcaldera magmatism. Plutons thus provide an integrated record of prolonged magmatic evolution, while volcanism offers snapshots of conditions at early stages. Growth of subvolcanic batholiths involved sustained multistage open-system processes. These commonly involved ignimbrite eruptions at times of peak power input, but assembly and consolidation processes continued at diminishing rates long after peak volcanism. Some evidence cited for early incremental pluton assembly more likely records late events during or after volcanism. Contrasts between relatively primitive arc systems dominated by andesitic compositions and small upper-crustal plutons versus more silicic volcanic fields and associated batholiths probably reflect intertwined contrasts in crustal thickness and magmatic power input. Lower power input would lead to a Cascade- or Aleutian-type arc system, where intermediate-composition magma erupts directly from middle- and lower-crustal storage without development of large shallow plutons. Andean and southern Rocky Mountain&ndash;type systems begin similarly with intermediate-composition volcanism, but increasing magma production, perhaps triggered by abrupt changes in plate boundaries, leads to development of larger upper-crustal reservoirs, more silicic compositions, large ignimbrites, and batholiths. Lack of geophysical evidence for voluminous eruptible magma beneath young calderas suggests that near-solidus plutons can be rejuvenated rapidly by high-temperature mafic recharge, potentially causing large explosive eruptions with only brief precursors.

Southern Rocky Mountains↗

The Schwarzwalder uranium deposit, III: Alteration, vein mineralization, light stable isotopes, and genesis of the deposit

The Schwartzwalder uranium deposit formed at 69.3 + or - 1.1 m.y. in a complex fracture system during the inception of Laramide uplift of the Front Range in Colorado. Geologic and isotopic evidence demonstrates that the ore-forming fluids were in chemical equilibrium with the Proterozoic metavolcanic and metasedimentary host-rock terrane at depth and that the metals, sulfur, and carbonate deposited in the fractures were derived from the metamorphic rocks. The data are not consistent with chemical contributions from an unrecognized magma or from the overlying Phanerozoic sedimentary rocks. The protoliths for the metamorphic rocks were submarine volcanic rocks and related volcanogenic exhalative iron-formations and chert. Water trapped along the basement faults and in the regolith between the basement and the overlying Paleozoic sedimentary rocks interacted with the metavolcanic rocks to produce isotopically heavy fluids containing high concentrations of dissolved metals and carbonate. Calculated delta 18 O values for this fluid range from 4.3 to 8.2 per mil, indicating a low water/rock ratio in the source terrane.Two stages of alteration and three stages of vein mineralization are recorded in the Schwartzwalder deposit. At the onset of Laramide faulting, fluids migrated along the fracture systems to zones of low hydraulic potential. These fluids contained CO 2 and had a metastably large K/Na ratio; they altered the gneissic wall rocks to a carbonate-sericite assemblage, adding K (super +) and CO 2 and removing SiO 2 with little or no change in volume. As the fractures continued to open, CO 2 was evolved from the fluids, increasing the pH and superimposing a hematite-adularia alteration assemblage on the earlier alteration.The veins record three stages of mineralization, the second of which generated the high-grade uranium veins. Evidence for the stage I sulfide-carbonate mineralization is poorly preserved, but isotopic and temperature data from this stage are consistent with a trend in fluid composition culminating in stage II pitchblende deposition. Sudden, large movements along the faults caused episodic evolution of CO 2 from the fluid. This loss of CO 2 decreased the solubilities of carbonates and adularia and the stabilities of sulfur species in solution. Uranyl carbonate complexes dissociated and sulfur species in solution likely reduced the uranyl ions to produce stage II pitchblende. Carbonate, adularia, and sulfides dominated the vein mineralogy after deposition of pitchblende. Progressively lower delta 18 O values in vein carbonates suggest the mixing of cooler, less evolved, perhaps meteoric, waters during the later stages of mineralization. However, the fluid pressure remained high, as indicated by explosion breccias and inward collapse features which formed as fault movements produced sudden decreases in the confining pressure. Only the stage III carbonate-iron disulfide assemblage in the major postore segment of the Illinois fault may be the product of meteoric water alone.

Economic Geology↗

Katmai volcanic cluster and the great eruption of 1912

In June 1912, the world's largest twentieth century eruption broke out through flat-lying sedimentary rocks of Jurassic age near the base of Trident volcano on the Alaska Peninsula. The 60 h ash-flow and Plinian eruptive sequence excavated and subsequently backfilled with ejecta a flaring funnel-shaped vent since called Novarupta. The vent is adjacent to a cluster of late Quaternary stratocones and domes that have released about 140 km 3 of magma in the past 150 k.y. Although the 1912 vent is closest to the Trident group and is also close to Mageik and Griggs volcanoes, it was the summit of Mount Katmai, 10 km east of Novarupta, that collapsed during the eruption to form a 5.5 km 3 caldera. Many earthquakes, including 14 in the range M 6−7, took place during and after the eruption, releasing 250 times more seismic energy than the 1991 caldera-forming eruption of the Philippine volcano, Pinatubo. The contrast in seismic behavior may reflect the absence of older caldera faults at Mount Katmai, lack of upward (subsidence opposing) magma flow owing to lateral magma withdrawal in 1912, and the horizontally stratified structure of the thick shale-rich Mesozoic basement. The Katmai caldera compensates for only 40% of the 13 km 3 of 1912 magma erupted, which included 7–8 km 3 of slightly zoned high-silica rhyolite and 4.5 km 3 of crystal-rich dacite that grades continuously into 1 km 3 of crystal-rich andesite. We have now mapped, sampled, and studied the products of all 20 components of the Katmai volcanic cluster. Pyroxene dacite and silicic andesite predominate at all of them, and olivine andesite is also common at Griggs, Katmai, and Trident volcanoes, but basalt and rhyodacite have erupted only at Mount Katmai. Rhyolite erupted only in 1912 and is otherwise absent among Quaternary products of the cluster. Pleistocene products of Mageik and Trident and all products of Griggs are compositionally distinguishable from those of 1912 at Novarupta. Holocene products of Mount Martin and Trident are closer in composition to the andesite-dacite array of 1912, but they reveal consistent differences. The affinity of the 1912 suite is closest with the array of products erupted by the Southwest Katmai cone, the edifice that had produced the only pre-1912 rhyodacite as well as the largest prehistoric Plinian eruption in the cluster. It is doubtful that any 1912 magma had been stored beneath Novarupta or Trident, and there is no evidence that more than one magma chamber erupted in 1912. Despite a compositional gap separating the aphyric rhyolite from the very crystal-rich andesite-dacite continuum, isotopic and chemical affinities linking all the 1912 ejecta and the continuity of all those ejecta in magmatic temperature and oxygen fugacity suggest that the rhyolite originated principally by incremental upward expulsion of interstitial melt from subjacent andesite-dacite mush. A large reservoir of such hot crystal mush is required both as the residue of rhyolitic melt separation and as a proximate heat source to thermally sustain the nearly aphyric condition of the overlying rhyolite. A model is presented for a unitary zoned chamber beneath Mount Katmai.

Alaska↗

Simulation of groundwater flow and brine discharge to the Dolores River in the Paradox Valley, Montrose County, Colorado

Salinity, or total dissolved solids (TDS), of the Colorado River affects agricultural, municipal, and industrial water users and is an important concern in the Western United States. In the Paradox Valley of southwestern Colorado, natural discharge of sodium-chloride brine to the Dolores River from the underlying core of a salt-valley anticline accounts for about 6 percent of the salinity load to the Colorado River. Formation of the Paradox Valley began during the Miocene, and subsequent erosion exposed the Pennsylvania Paradox Formation in the core of the anticline where a cap rock, collapse features, breccia, and sodium-chloride saturated brine developed at the top of the exposed salt diapir. The discharge of brine to the Dolores River is affected by these dissolution features, along with seasonal hydrologic conditions and density-dependent flow between older dense brine and the younger fresh groundwater in the overlying alluvial aquifer. To reduce TDS concentrations in the Dolores River through the Paradox Valley, the Bureau of Reclamation has pumped brine from a series of shallow wells adjacent to the river since July 1996. The pumped brine is collected and piped to a deep disposal well where it is injected into the Mississippian Leadville Limestone at a depth of about 4,570-meters below land surface. The pumping and injection operation is collectively known as the Paradox Valley Unit (PVU), and by 2015, the PVU had substantially reduced TDS concentrations in the Dolores River by about 70 percent. Since 2019, injection-pressure limits and related seismic activity have constrained deep-well injection and thus brine pumping at the PVU. In cooperation with the Bureau of Reclamation, the U.S. Geological Survey developed a MODFLOW-6 three-dimensional, variable-density groundwater flow and TDS transport model of the Paradox Valley to evaluate the effects of PVU pumping operations on brine discharge to the Dolores River and to guide additional research. The finite-difference model grid consists of 76 rows and 48 columns oriented from northwest to southeast in alignment with valley topography and groundwater-flow directions in the near-surface freshwater alluvial aquifer. A 7-layer hydrogeologic framework was developed from existing datasets to represent the alluvial aquifer, cap rock, collapse breccia, and groundwater flow and TDS transport from the underlying Paradox Formation salt to the Dolores River. The model represents a 33-year transient calibration period from 1987 through 2020 that includes pre-PVU conditions from 1987 through June 1996 and post-PVU conditions from July 1996 through 2020. A 1,000-year simulation of groundwater flow and coupled TDS transport computed the initial conditions for the subsequent 33-year transient simulation. Observations of precipitation, streamflow, evaporation, agricultural land use, and PVU brine pumping rates were used to specify appropriate boundary conditions to the model representing time-varying recharge, tributary streamflow, groundwater underflow, evapotranspiration (ET), and PVU pumping. Values for average monthly streamflow and TDS concentration at the upstream streamgage, the Dolores River at Bedrock (USGS streamgage 09169500), were specified as model input where the Dolores River enters Paradox Valley. Observed pumping from the PVU, water levels and TDS concentrations in groundwater, and streamflow and estimated TDS concentrations at the downstream streamgage, the Dolores River near Bedrock (USGS streamgage 09171100), were calibration targets that constrained the manual calibration of model parameters representing aquifer hydraulic conductivity, storage, streambed conductance, recharge, and (ET). Two primary model-calibration targets were the match between observed and simulated TDS mass flux from PVU pumping wells and the match between estimated and simulated TDS mass flux to the Dolores River. The simulated TDS mass withdrawn by pumping wells is calculated by the model as the product of the assigned pumping rate and simulated groundwater TDS concentrations. Because actual pumping rates were assigned as simulated values, the total simulated PVU pumping for the 33-year calibration is within 0.5 percent of the observed values. However, simulated concentrations and thus mass flux of TDS withdrawn by the PVU pumping wells were consistently about 26 percent less than observed values for all the simulated time periods (33-year simulation, pre-PVU, and post-PVU). The representation of brine inflow was explored through additional modeling to evaluate the effect of the simulated brine source on groundwater TDS concentrations. Results indicated that a saturated-salt constant-flux brine source best replicated the magnitude and transient pattern observed for TDS mass flux from PVU pumping wells. The simulated TDS mass flux to the Dolores River is compared to estimates based on observed streamflow and specific conductance (SC) data for the downstream streamgage. The calibrated model provided a close fit of simulated to measured streamflow at the downstream streamgage, and the calibrated model fit to estimated TDS concentrations at the downstream streamgage was reasonable. The greatest differences between simulated and estimated values occurred during drought periods from June 2000 to March 2003, May 2012 to June 2013, and October 2013 to October 2014, when simulated TDS concentrations in the river were greater than estimated concentrations. In general, simulated TDS mass flux to the river for the pre-PVU period is in good agreement with estimated values (2-percent difference), but the model overestimated TDS mass flux to the river by about 41 percent during the post-PVU period. The model uncertainty with respect to TDS mass flux to the river indicates other processes or model parameters not well represented by the model are affecting the system, especially during drought. During model calibration, the most sensitive parameters were identified as vertical hydraulic conductivity of the alluvial aquifer, conductance of the Dolores River streambed, ET extinction depth and rate, and recharge rate. Five 5-year scenarios of conditions for 2021–25 were simulated to assist evaluation of alternative strategies to manage the discharge of brine into the Dolores River. The first scenario simulates no PVU pumping and serves as a base case for comparison to the other scenarios. Two scenarios simulate the effects of varying withdrawal timing at an annual rate about one-third less than during 2010 through 2018. During high-flow spring snowmelt runoff periods when brine discharge is naturally minimized, PVU pumping does not substantially affect salinity in the Dolores River, and comparison of these two scenarios indicates that scheduling brine withdrawals during times of low river stage is nearly as effective at reducing TDS mass flux to the river as pumping brine year-round. Cessation of pumping during periods of high river stage may be advantageous for system maintenance, brine injection, and seismic-risk reduction. The fourth scenario tested the effect of reducing irrigation-return flow on brine discharge and predicted a slight reduction of TDS mass flux to the Dolores River, but not as great a reduction as that of using the PVU to remove brine. The fifth scenario simulated 5 years of drought conditions without PVU pumping and indicates brine discharge during drought about 15 percent greater than during average hydrologic conditions. Results from scenario 5 are consistent with the calibrated model results and indicate that aquifer properties and ET processes and parameters may be affecting simulation results during drought. The Paradox Valley groundwater model provides a reasonable overall match to observed conditions in the Dolores River. The model is useful for evaluating relative differences between brine management scenarios to inform PVU operational decisions and to identify gaps in data and process understanding. Representation of the brine source, hydraulic-conductivity parameters, and recharge and ET processes were identified as potential areas for additional field and modeling research. Additional research in the Paradox Valley might include field-data collection that provides additional information on the hydrogeologic framework, groundwater levels, groundwater TDS concentrations, stream characteristics, and aquifer properties. Additional modeling efforts could benefit from applying advanced tools for model development, calibration, and visualization including parameter-estimation and sensitivity analysis. Statistical evaluation of known model uncertainties such as hydraulic conductivity, streambed conductance, representations of the brine source, recharge, and ET could improve the match between simulated and estimated TDS mass flux from PVU pumping wells and to the Dolores River further informing model predictions and system understanding for the Paradox Valley.

Colorado↗

Update of the 2008 provisional Enhanced Geothermal Systems (EGS) assessment for the Great Basin, USA

In response to the Energy Act of 2020, the U.S. Geological Survey (USGS) is updating the Enhanced Geothermal Systems (EGS) resource assessment for the Great Basin, USA. The previous 2008 provisional assessment estimated how much electricity could be generated from EGS resources of the western United States using models of electric-grade heat, models of heat extraction over time, and estimates of how much rock might be stimulated to produce viable amounts of heat. Herein, a similar conceptual strategy is applied, using updated models of heat extraction as a function of fracture spacing and well distance. Previously used reservoir heat delivery models are updated to have a dependence on fracture and well spacing, potentially improving future estimates of EGS resources as ongoing research provides a better understanding about the success of reservoir stimulation as a function of geology and location. For a range of well distances (250-1000 m) and fracture spacings (1-50 m), heat extraction efficiency ranges from 25-62%, demonstrating the importance of accounting for the most likely results of proven viable fracturing technologies. Although fracturing is important, the biggest uncertainty by far in estimating the EGS resource for the Great Basin is estimating which geologic units at what depths can be stimulated sufficiently to produce geothermal energy economically and efficiently. Uncertainties in these factors yield estimates that range over two orders of magnitude with an upper limit of ~174 terawatts-thermal (TWth) produced for 30 years from the upper 7 km of the crust. This upper limit would require significant technological advances to access most of the electric-grade resource across the Great Basin. Assuming that 1% of this estimate will be accessible in the next few decades gives a resource estimate similar to that made in the 2008 provisional assessment. These estimated EGS heat extraction rates far exceed ( greater than 100x) the natural geothermal heat production rate, thus geothermal electricity production at these rates might not be sustainable unless heat is also recharged from other sources (e.g., excess solar energy when supply exceeds demand). In addition to assessment maps and cumulative estimates, the new models of fractured reservoirs developed herein can be used to estimate steady power production given a set of fractures and well spacing, and estimates can be made for setback distances to ensure no thermal interference with nearby powerplants.

Arizona, California, Idaho, Nevada, Oregon, Utah↗

To mix or not to mix: Details of magma storage, recharge, and remobilization during the Pacheco stage at Misti Volcano, Peru (≤21–2 ka)

We investigate ten of the most recent tephra-fall deposits emplaced between ≤21–2 ka from the Pacheco stage of Misti volcano, Peru, to elucidate magma dynamics and explosive eruption triggers related to magma storage, recharge, and remobilization. Whole-rock, glass, and mineral textures and compositions indicate the presence of broadly felsic, intermediate, and mafic magmas in a chemically and thermally stratified magma storage system (Zones 1–3) that interact to differing extents prior to eruption. Intermediate magmas are defined by plagioclase + amphibole + two-pyroxenes + Fe-Ti oxides and phase equilibria indicate they formed at ~300–600 MPa and ~950–1000 °C. Intermediate magmas dominate the Pacheco stage and either erupted alone as hybridized magmas or mingled with minor volumes of cool felsic magmas (~800 °C) in which only plagioclase + Fe-Ti oxides are stable. Felsic magmas do not exclusively comprise any tephra-fall deposit emplaced during the Pacheco stage but were remobilized by recharge and mixing with intermediate magmas in order to erupt. Furthermore, felsic-hosted amphibole cognate to the intermediate magmas are reacted despite the felsic magmas being water saturated, which suggests they are staged above the amphibole stability limit (≤200 MPa). The cryptic presence of mafic magmas is indicated by high-An plagioclase cores (An 74–88 ), rare anhedral olivine (Fo 77–80 ), and possibly high Mg# augite and amphibole (up to Mg# 84 and 77, respectively). The dearth of basalt to basaltic andesite melts recorded in erupted glasses and exclusivity of high-An plagioclase to crystal cores signals mafic magmas are staged deeper in the crust than the intermediate magmas. Periodic interactions between these magmas tracked via glass compositions and crystal exchange reveal an alternation between the production of mingled magmas and their eruption shortly after a recharge event, followed by a period of homogenization and eruption of hybridized magmas. As such, we identify magma recharge as a key mechanism by which half of the explosive eruptions were triggered in the Pacheco stage. A >100 °C increase in Misti’s fumarole temperatures from 1967 to 2018 coincident with changes in fumarolic gas composition is consistent with degassing of a mafic recharge magma, signaling that Misti could produce similar explosive eruptions in the future.

Misti Volcano↗

Succession of Laramide magmatic and magmatic-hydrothermal events in the Patagonia Mountains, Santa Cruz County, Arizona

This investigation of the space-time progression of magmatism and hydrothermal activity in the Patagonia Mountains of southern Arizona is based on field and paragenetic relationships, and on U-Pb and 40 Ar/ 39 Ar geochronology of igneous and hydrothermal minerals. The Patagonia Mountains consist of Precambrian, Paleozoic, and Mesozoic sedimentary, granitic, and volcanic rocks, Laramide volcanic rocks, and a core of Laramide intrusions that comprise the Patagonia Mountains batholith. Laramide igneous rocks and adjacent Paleozoic and Mesozoic rocks contain significant porphyry Cu-Mo deposits, Mo-Cu breccia pipes, Ag replacement deposits, and numerous other Cu-Pb-Zn-Ag replacement and vein deposits. Ages of igneous and hydrothermal minerals from 20 U-Pb and 52 40 Ar/ 39 Ar determinations define four magmatic and magmatic-hydrothermal events that formed the batholith and altered parts of it and adjacent rocks; cumulatively the events span at least 16 m.y., from ~74 to 58 Ma. The oldest event of this succession includes the 74 Ma Washington Camp stock and spatially associated Cu-Pb-Zn-Ag replacement deposits in Paleozoic carbonate rocks of the Washington Camp-Duquesne district. Eruption of 73 to 68 Ma volcanic rocks in the northern part of the range was the next youngest event, which coincides temporally with replacement and vein deposits in Paleozoic carbonate rocks at the Flux mine (~71 Ma). An event at 65 to 62 Ma is marked by emplacement of small-volume quartz monzonite, granodiorite, and diorite intrusions, formation of the Ventura breccia deposit in Jurassic granite at 65 to 64 Ma, and formation of other Pb-Zn-Ag-Cu replacement and vein deposits (~62 Ma; Blue Nose and Morning Glory). The Red Mountain porphyry Cu-Mo system is hosted by ~62 Ma granodiorite and Laramide volcanic rocks (73–68 Ma) at the northern end of the batholith. It includes a deep, chalcopyrite-bornite resource (~60 Ma) that is associated with potassic and sericitic alteration and a near-surface chalcocite-enargite resource (60 Ma) that is associated with advanced, supergene-enriched argillic alteration. The youngest event includes the Sunnyside porphyry Cu-Mo system and a Cu-Mo breccia deposit at Red Hill (Four Metals mine), both of which formed in large-volume quartz monzonite, granodiorite, quartz monzonite porphyry, and quartz feldspar porphyry (~61–59 Ma). Similar to the Red Mountain system, the Sunnyside system consists of a deep chalcopyrite resource that occurs in ~60 to 59 Ma quartz feldspar porphyry, and a near-surface, slightly younger (~59–58 Ma) enargite-chalcocite-tennantite resource that occurs in quartz feldspar porphyry, quartz monzonite porphyry, and Mesozoic rocks. The Red Hill Cu-Mo breccia deposit is hosted by large-volume quartz monzonite, granodiorite, and quartz monzonite porphyry (~61–59 Ma). Discrepancies between field and paragenetic relationships and some analytic ages at Sunnyside and Red Hill preclude precise dating of mineralization stages, and may reflect disturbance of isotope systems by multiple, co-spatial to juxtaposed intrusive and hydrothermal events, and/or by unrecognized intrusions. Numerous vein and replacement deposits at the northern end of the batholith, including the Hardshell Ag resource and the Three R supergene chalcocite resource, are distal deposits of the Sunnyside and Red Mountain systems. Small, ~61 to 59 Ma Cu-Mo deposits in large-volume intrusions in the southern part of the batholith consist of hydrothermal quartz, biotite, K-feldspar, muscovite, chalcopyrite, and molybdenite. The age span of magmatic and magmatic-hydrothermal events in the Patagonia Mountains, minimally 16 m.y., is comparable to that of certain other magmatic-hydrothermal successions that contain porphyry Cu-Mo systems. Magmatic-hydrothermal events of the Wasatch-Oquirrh igneous trend, Utah, and the Boulder batholith, Montana, both span ~17 m.y. and include the Bingham and Butte porphyry Cu-Mo, vein and replacement deposits, respectively. Plutons and mineral deposits in the Pima district, Arizona, which includes the porphyry Cu-Mo deposits at Sierrita-Esperanza, Mission-Pima-San Xavier North, and Twin Buttes, formed over an interval of ~14 m.y. The diversity of igneous and hydrothermal products likely reflects evolutionary processes occurring at multiple sites in the lithosphere and at different time scales from >10 m.y. to less than the geochronologic precision currently achievable.

Arizona↗

Core Research Center

The Core Research Center (CRC) of the U.S. Geological Survey (USGS), located at the Denver Federal Center in Lakewood, Colo., currently houses rock core from more than 8,500 boreholes representing about 1.7 million feet of rock core from 35 States and cuttings from 54,000 boreholes representing 238 million feet of drilling in 28 States. Although most of the boreholes are located in the Rocky Mountain region, the geologic and geographic diversity of samples have helped the CRC become one of the largest and most heavily used public core repositories in the United States. Many of the boreholes represented in the collection were drilled for energy and mineral exploration, and many of the cores and cuttings were donated to the CRC by private companies in these industries. Some cores and cuttings were collected by the USGS along with other government agencies. Approximately one-half of the cores are slabbed and photographed. More than 18,000 thin sections and a large volume of analytical data from the cores and cuttings are also accessible. A growing collection of digital images of the cores are also becoming available on the CRC Web site Internet http://geology.cr.usgs.gov/crc/.

General Information Product↗

Reconnaissance of the chemical quality of surface waters of the Neches River basin, Texas

The kinds and quantities of minerals dissolved in the surface water of the Neches River basin result from such environmental factors as geology, streamflow patterns and characteristics, and industrial influences. As a result of high rainfall in the basin, much of the readily soluble material has been leached from the surface rocks and soils. Consequently, the water in the streams is usually low in concentrations of dissolved minerals and meets the U.S. Public Health Service drinking-water standards. In most streams the concentration of dissolved solids is less than 250 ppm (parts per million). The Neches River drains an area of about 10,000 square miles in eastern Texas. From its source in southeast Van Zandt County the river flows in a general southeasterly direction and empties into Sabine Lake, an arm of the Gulf of Mexico. In the basin the climate ranges from moist subhumid to humid, and the average annual rainfall ranges from 46 inches is the northwest to more than 52 inches in the southeast. Annual runoff from the basin has averaged 11 inches; however, runoff rates vary widely from year to year. The yearly mean discharge of the Neches River at Evadale has ranged from 994 to 12,720 cubic feet per second. The rocks exposed in the Neches River basin are of the Quaternary and Tertiary Systems and range in age from Eocene to Recent. Throughout most of the basin the geologic formations dip generally south and southeast toward the gulf coast. The rate of dip is greater than that of the land surface; and as a result, the older formations crop out to the north of the younger formations. Water from the outcrop areas of the Wilcox Group and from the older formations of the Claiborne Group generally has dissolved-solids concentrations ranging from 100 to 250 ppm; water from the younger formations has concentrations less than 100 ppm. The northern half of the basin has soft water, with less than 60 ppm hardness. The southern half of .the basin has very soft water, usually with less than 30 ppm hardness. The chloride concentrations are less than 20 ppm in surface water in the southern half of the basin and usually range from 20 to 100 ppm in the northern half of the basin. Concentrations greater than 100 ppm are found only where pollution is occurring. The Neches River basin has an abundance of surface water, but uneven distribution of runoff makes storage projects necessary to provide dependable water supplies. The principal existing reservoirs, with the exception of Striker Creek Reservoir, contain water of excellent quality. Chemical-quality data for the Striker Creek drainage area indicate that its streams are affected by .the disposal of brines associated with oil production. Sam Rayburn Reservoir began impounding water in 1965. The water impounded should prove of acceptable quality for most uses, but municipal and industrial wastes released into the Angelina River near Lufkin may have a degrading effect on the quality of the water, especially during extended periods of low flows. Water available for storage at the many potential reservoir sites will be of good quality; but, if the proposed salt-water barrier is to impound acceptable water, the disposal of oilfield brine into Pine Island Bayou should be discontinued.

Water Supply Paper↗

Porphyry copper assessment of British Columbia and Yukon Territory, Canada: Chapter C in Global mineral resource assessment

The U.S. Geological Survey does regional, national, and global assessments of resources (mineral, energy, water, biologic) to provide science in support of land management and decision making. Mineral resource assessments provide a synthesis of available information about where mineral deposits are known and suspected to be in the Earth&rsquo;s crust, which commodities may be present, and estimates of amounts of resources that may be present in undiscovered deposits. Canada is an important source of copper, consistently ranking as one of the top 10 world producers during the past decade (2000&ndash;2010). The preponderance of this production has been from porphyry-copper-type deposits in the western Canadian Cordillera. A probabilistic mineral resource assessment of undiscovered resources associated with porphyry copper deposits in western Canada was completed as part of a global mineral resource assessment. The purpose of the assessment was to (1) compile a database of known deposits and significant prospects, (2) delineate permissive areas (tracts) for undiscovered porphyry copper deposits that may be present in the upper kilometer (minimally) of the Earth&rsquo;s crust, and (3) provide probabilistic estimates of amounts of copper (Cu), molybdenum (Mo), gold (Au), and silver (Ag) that could be contained in undiscovered porphyry copper deposits in the tracts. The study was done by the U.S. Geological Survey (USGS) in collaboration with geologists from the British Columbia Geological Survey, Yukon Geological Survey, and industry consultants. The database of known deposits and significant prospects includes an inventory of mineral resources in 89 known porphyry copper (and 2 related copper-bearing polymetallic vein) ore zones, representing 50 porphyry copper deposits, and lists key characteristics of 280 additional porphyry copper and related copper-bearing prospects, as indicated by currently available exploration results, which also are summarized. Resource and exploration and development activity are updated with information current through April 2010. The delineation of permissive tracts and probabilistic estimation of resources in undiscovered deposits were done using the USGS three-part mineral resource assessment approach. Permissive tracts are defined in accordance with descriptive models for porphyry copper deposits to include igneous rocks and known deposits and prospects within magmatic arcs related to convergent plate-margin boundary zones. Frequency distributions of total tonnages and average grades of thoroughly explored deposits were used as models for undiscovered deposits and include a new grade and tonnage model for calc-alkaline porphyry Cu&plusmn;Mo&plusmn;Au deposits in western Canada. Five permissive tracts for the occurrence of porphyry copper deposits were delineated: 2 island-arc tracts, 1 tract of transitional, mixed island-arc and continental arc affinities, and 2 continental arc tracts. In permissive tract 003pCu2001, calc-alkaline igneous rocks of Middle Triassic to Late Jurassic age in accreted island-arc terranes of the Intermontane belt are assessed for calc-alkaline porphyry Cu&plusmn;Mo&plusmn;Au deposits. The area of this tract is 175,250 km 2 . In 12 known deposits, the total reported tonnage of ore is 8,100 million metric tons (Mt) containing 24.6 Mt copper. An estimated 6.9 undiscovered deposits contain a calculated mean of 8.9 Mt copper and a median of 6.9 Mt copper. The spatial density for the 18.9 known plus estimated undiscovered deposits in this tract is approximately 11 deposits per 100,000 km 2 . In permissive tract 003pCu2002, alkaline igneous rocks of Middle Triassic to Late Jurassic age within the Intermontane accreted island-arc terranes are assessed for alkaline porphyry Cu-Au deposits. The area of this tract is 109,290 km 2 . In 12 known deposits the total reported tonnage of ore is 6,440 Mt, containing 20.9 Mt copper. An estimated 7 undiscovered deposits contain a calculated mean of 22 Mt copper and a median of 13 Mt copper. The spatial density for the 19 known plus estimated undiscovered deposits in this tract is approximately 17 deposits per 100,000 km 2 . In permissive tract 003pCu2003, calc-alkaline igneous rocks of Late Triassic to Early Cretaceous age within the accreted Insular terranes of mixed island-arc and continental arc affinities are assessed for calc-alkaline porphyry Cu&plusmn;Mo&plusmn;Au deposits. The area of this tract is 58,360 km 2 . The total tonnage of ore reported in the 2 known deposits is 1,160 Mt containing 3.17 Mt copper. An estimated 2.3 undiscovered deposits contain a calculated mean of 3 Mt copper and a median of 1.9 Mt copper. The spatial density for the 4.3 known plus estimated undiscovered deposits in this tract is approximately 7 deposits per 100,000 km 2 . In permissive tract 003pCu2004, calc-alkaline igneous rocks in continental magmatic arcs of Jurassic to Eocene age are assessed for porphyry Cu&plusmn;Mo&plusmn;Au deposits. The area of this tract is 639,500 km 2 . The total tonnage of ore reported for the 23 known deposits is 6,520 Mt containing 17.9 Mt copper. An estimated 9.6 undiscovered deposits contain a calculated mean of 13 Mt copper and a median of 11 Mt copper. The spatial density for the 32.6 known deposits plus the estimated undiscovered deposits in this tract is approximately 5 deposits per 100,000 km 2 . In permissive tract 003pCu2005, calc-alkaline igneous rocks in continental magmatic arcs of Oligocene to Pliocene age are assessed for porphyry Cu&plusmn;Mo&plusmn;Au deposits. The area of this tract is 32,840 km 2 . The total tonnage of ore reported for the 1 known deposit is 44.8 Mt containing 0.224 Mt copper. An estimated 1.4 undiscovered deposits contain a calculated mean of 1.8 Mt copper and a median of 0.72 Mt copper. The spatial density for the 2.4 known plus estimated undiscovered deposits in this permissive tract is approximately 7 deposits per 100,000 km 2 . Western Canada has been thoroughly explored for porphyry copper deposits. The total estimated copper contained in known deposits is about 66.8 Mt (based on 2010 data), as compared to a 49 Mt mean of estimated copper in undiscovered deposits and a 34 Mt median of estimated copper in undiscovered deposits. The copper contained in known porphyry copper deposits represents about 58 percent of the total of known and undiscovered porphyry copper deposits (based on mean values). About 86 percent of the increase in estimated copper resources between 1993 and 2009 resulted from the discovery of extensions to known deposits. Nevertheless, exploration for undiscovered deposits continues, especially in and around significant prospects and in parts of permissive tracts that are mostly hidden beneath younger volcanic, sedimentary, or vegetated surficial cover.

British Columbia, Yukon Territory↗