Search USGS⌕ Search

SEARCH · Search USGS

Results for “Interfaces”

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 433 records · Page 24Linked to original sources
USGS2016

A decision framework for identifying models to estimate forest ecosystem services gains from restoration

Restoring degraded forests and agricultural lands has become a global conservation priority. A growing number of tools can quantify ecosystem service tradeoffs associated with forest restoration. This evolving “tools landscape” presents a dilemma: more tools are available, but selecting appropriate tools has become more challenging. We present a Restoration Ecosystem Service Tool Selector (RESTS) framework that describes key characteristics of 13 ecosystem service assessment tools. Analysts enter information about their decision context, services to be analyzed, and desired outputs. Tools are filtered and presented based on five evaluative criteria: scalability, cost, time requirements, handling of uncertainty, and applicability to benefit-cost analysis. RESTS uses a spreadsheet interface but a web-based interface is planned. Given the rapid evolution of ecosystem services science, RESTS provides an adaptable framework to guide forest restoration decision makers toward tools that can help quantify ecosystem services in support of restoration.

Forest Ecosystems↗
USGS1994

Large-explosive source, wide-recording aperture, seismic profiling on the Columbia Plateau, Washington

Clear subsurface seismic images have been obtained at low cost on the Columbia Plateau, Washington. The Columbia Plateau is perhaps the most notorious of all 'bad-data' areas because large impedance contrasts in surface flood basalts severely degrade the seismic wavefield. This degradation was mitigated in this study via a large-explosive source, wide-recording aperture shooting method.The shooting method emphasizes the wide-angle portion of the wavefield, where Fermat's principle guarantees reverberation will not interfere with the seismic manifestations of crucial geologic interfaces. The basalt diving wave, normally discarded in standard common midpoint (CMP) seismic profiling, can be used to image basalt velocity structure via traveltime inversion. Maximum depth-penetration of the diving wave tightly constrains basalt-sediment interface depth. An arrival observed only at shot-receiver offsets greater than 15 km can be used to determine the velocity and geometry of basement via simultaneous inversion.The results from this study suggest that previous geologic hypotheses and hydrocarbon play concepts for the Columbia Plateau may have been in error.

Geophysics↗
USGS1994

An image of the Columbia Plateau from inversion of high‐resolution seismic data

We use a method of traveltime inversion of high‐resolution seismic data to provide the first reliable images of internal details of the Columbia River Basalt Group (CRBG), the subsurface basalt/sediment interface, and the deeper sediment/basement interface. Velocity structure within the basalts, delineated on the order of 1 km horizontally and 0.2 km vertically, is constrained to within ±0.1 km/s for most of the seismic profile. Over 5000 observed traveltimes fit our model with an rms error of 0.018 s. The maximum depth of penetration of the basalt diving waves (truncated by underlying low‐velocity sediments) provides a reliable estimate of the depth to the base of the basalt, which agrees with well‐log measurements to within 0.05 km (165 ft). We use image blurring, calculated from the resolution matrix, to estimate the aspect ratio of imaged velocity anomaly widths to true widths for velocity features within the basalt. From our calculations of image blurring, we interpret low velocity zones (LVZ) within the basalts at Boylston Mountain and the Whiskey Dick anticline to have widths of 4.5 and 3 km, respectively, within the upper 1.5 km of the model. At greater depth, the widths of these imaged LVZs thin to approximately 2 km or less. We interpret these linear, subparallel, low‐velocity zones imaged adjacent to anticlines of the Yakima Fold Belt to be brecciated fault zones. These fault zones dip to the south at angles between 15 to 45 degrees.

Geophysics↗
USGS2007

Seismoelectric imaging of shallow targets

We have undertaken a series of controlled field experiments to develop seismoelectric experimental methods for near-surface applications and to improve our understanding of seismoelectric phenomena. In a set of off-line geometry surveys (source separated from the receiver line), we place seismic sources and electrode array receivers on opposite sides of a man-made target (two sand-filled trenches) to record separately two previously documented seismoelectric modes: (1) the electromagnetic interface response signal created at the target and (2) the coseismic electric fields located within a compressional seismic wave. With the seismic source point in the center of a linear electrode array, we identify the previously undocumented seismoelectric direct field, and the Lorentz field of the metal hammer plate moving in the earth's magnetic field. We place the seismic source in the center of a circular array of electrodes (radial and circumferential orientations) to analyze the source-related direct and Lorentz fields and to establish that these fields can be understood in terms of simple analytical models. Using an off-line geometry, we create a multifold, 2D image of our trenches as dipping layers, and we also produce a complementary synthetic image through numerical modeling. These images demonstrate that off-line geometry (e.g., crosswell) surveys offer a particularly promising application of the seismoelectric method because they effectively separate the interface response signal from the (generally much stronger) coseismic and source-related fields. ?? 2007 Society of Exploration Geophysicists.

Geophysics↗
USGS2019

Shear-wave seismic reflection studies of unconsolidated sediments in the near surface

We have successfully applied of SH-wave seismic reflection methods to two different near-surface problems targeting unconsolidated sediments. At the former Fort Ord, where the water table is approximately 30m deep, we imaged aeolian and marine aquifer and aquitard stratigraphy to a depth of approximately 80m. We identified reflections from sand/clay and sand/silt interfaces and we mapped these interfaces along our transects. At an aggregate study site in Indiana, where the water table is at a depth of 1to2m, we imaged stratigraphy in alluvial sand and gravel, and observe a strong reflection from the 32-m-deep bedrock surface. In both cases, we exploited the high resolution potential of SH waves, their insensitivity to water content, and the possibility of reducing Love wave contamination by working along a roadway. We accomplished our results using only sledgehammer sources and simple data processing flows.

California, Indiana↗
USGS2015

Ground motion models used in the 2014 U.S. National Seismic Hazard Maps

The National Seismic Hazard Maps (NSHMs) are an important component of seismic design regulations in the United States. This paper compares hazard using the new suite of ground motion models (GMMs) relative to hazard using the suite of GMMs applied in the previous version of the maps. The new source characterization models are used for both cases. A previous paper ( Rezaeian et al. 2014 ) discussed the five NGA-West2 GMMs used for shallow crustal earthquakes in the Western United States (WUS), which are also summarized here. Our focus in this paper is on GMMs for earthquakes in stable continental regions in the Central and Eastern United States (CEUS), as well as subduction interface and deep intraslab earthquakes. We consider building code hazard levels for peak ground acceleration (PGA), 0.2-s, and 1.0-s spectral accelerations (SAs) on uniform firm-rock site conditions. The GMM modifications in the updated version of the maps created changes in hazard within 5% to 20% in WUS; decreases within 5% to 20% in CEUS; changes within 5% to 15% for subduction interface earthquakes; and changes involving decreases of up to 50% and increases of up to 30% for deep intraslab earthquakes for most U.S. sites. These modifications were combined with changes resulting from modifications in the source characterization models to obtain the new hazard maps.

Earthquake Spectra↗
USGS2024

The role of clay minerals in the preservation of organic matter in sediments of Qinghai Lake, NW China

The role of saline lake sediments in preserving organic matter has long been recognized. In order to further understand the preservation mechanisms, the role of clay minerals was studied. Three sediment cores, 25, 57, and 500 cm long, were collected from Qinghai Lake, NW China, and dissected into multiple subsamples. Multiple techniques were employed, including density fractionation, X-ray diffraction, scanning and transmission electron microscopy (SEM and TEM), total organic carbon (TOC) and carbon compound analyses, and surface area determination. The sediments were oxic near the water-sediment interface, but became anoxic at depth. The clay mineral content was as much as 36.8%, consisting mostly of illite, chlorite, and halloysite. The TEM observations revealed that organic matter occurred primarily as organic matter-clay mineral aggregates. The TOC and clay mineral abundances are greatest in the mid-density fraction, with a positive correlation between the TOC and mineral surface area. The TOC of the bulk sediments ranges from 1 to 3% with the non-hydrocarbon fraction being predominant, followed by bitumen, saturated hydrocarbon, aromatic hydrocarbons, and chloroform-soluble bitumen. The bimodal distribution of carbon compounds of the saturated hydrocarbon fraction suggests that organic matter in the sediments was derived from two sources: terrestrial plants and microorganisms/algae. Depthrelated systematic changes in the distribution patterns of the carbon compounds suggest that the oxidizing conditions and microbial abundance near the water-sediment interface promote degradation of labile organic matter, probably in adsorbed form. The reducing conditions and small microbial biomass deeper in the sediments favor preservation of organic matter, because of the less labile nature of organic matter, probably occurring within clay mineral-organic matter aggregates that are inaccessible to microorganisms. These results have important implications for our understanding of mechanisms of organic matter preservation in saline lake sediments.

Qinghai Lake↗
USGS2012

A global earthquake discrimination scheme to optimize ground-motion prediction equation selection

We present a new automatic earthquake discrimination procedure to determine in near-real time the tectonic regime and seismotectonic domain of an earthquake, its most likely source type, and the corresponding ground-motion prediction equation (GMPE) class to be used in the U.S. Geological Survey (USGS) Global ShakeMap system. This method makes use of the Flinn–Engdahl regionalization scheme, seismotectonic information (plate boundaries, global geology, seismicity catalogs, and regional and local studies), and the source parameters available from the USGS National Earthquake Information Center in the minutes following an earthquake to give the best estimation of the setting and mechanism of the event. Depending on the tectonic setting, additional criteria based on hypocentral depth, style of faulting, and regional seismicity may be applied. For subduction zones, these criteria include the use of focal mechanism information and detailed interface models to discriminate among outer-rise, upper-plate, interface, and intraslab seismicity. The scheme is validated against a large database of recent historical earthquakes. Though developed to assess GMPE selection in Global ShakeMap operations, we anticipate a variety of uses for this strategy, from real-time processing systems to any analysis involving tectonic classification of sources from seismic catalogs.

Bulletin of the Seismological Society of America↗
USGS2014

Seismic‐wave attenuation determined from tectonic tremor in multiple subduction zones

Tectonic tremor provides a new source of observations that can be used to constrain the seismic attenuation parameter for ground‐motion prediction and hazard mapping. Traditionally, recorded earthquakes of magnitude ∼3–8 are used to develop ground‐motion prediction equations; however, typical earthquake records may be sparse in areas of high hazard. In this study, we constrain the distance decay of seismic waves using measurements of the amplitude decay of tectonic tremor, which is plentiful in some regions. Tectonic tremor occurs in the frequency band of interest for ground‐motion prediction (i.e., ∼2–8 Hz) and is located on the subducting plate interface, at the lower boundary of where future large earthquakes are expected. We empirically fit the distance decay of peak ground velocity from tremor to determine the attenuation parameter in four subduction zones: Nankai, Japan; Cascadia, United States–Canada; Jalisco, Mexico; and southern Chile. With the large amount of data available from tremor, we show that in the upper plate, the lower crust is less attenuating than the upper crust. We apply the same analysis to intraslab events in Nankai and show the possibility that waves traveling from deeper intraslab events experience more attenuation than those from the shallower tremor due to ray paths that pass through the subducting and highly attenuating oceanic crust. This suggests that high pore‐fluid pressure is present in the tremor source region. These differences imply that the attenuation parameter determined from intraslab earthquakes may underestimate ground motion for future large earthquakes on the plate interface.

Cascadia↗
USGS2016

Toppling analysis of the Echo Cliffs precariously balanced rock

Toppling analysis of a precariously balanced rock (PBR) can provide insight into the nature of ground motion that has not occurred at that location in the past and, by extension, can constrain peak ground motions for use in engineering design. Earlier approaches have targeted 2D models of the rock or modeled the rock–pedestal contact using spring‐damper assemblies that require recalibration for each rock. Here, a method to model PBRs in 3D is presented through a case study of the Echo Cliffs PBR. The 3D model is created from a point cloud of the rock, the pedestal, and their interface, obtained using terrestrial laser scanning. The dynamic response of the model under earthquake excitation is simulated using a rigid‐body dynamics algorithm. The veracity of this approach is demonstrated through comparisons against data from shake‐table experiments. Fragility maps for toppling probability of the Echo Cliffs PBR as a function of various ground‐motion parameters, rock–pedestal interface friction coefficient, and excitation direction are presented. These fragility maps indicate that the toppling probability of this rock is low (less than 0.2) for peak ground acceleration (PGA) and peak ground velocity (PGV) lower than 3 m/s 2 and 0.75 m/s, respectively, suggesting that the ground‐motion intensities at this location from earthquakes on nearby faults have most probably not exceeded the above‐mentioned PGA and PGV during the age of the PBR. Additionally, the fragility maps generated from this methodology can also be directly coupled with existing probabilistic frameworks to obtain direct constraints on unexceeded ground motion at a PBR’s location.

California↗
USGS2021

The impact of 3D finite‐fault information on ground‐motion forecasting for earthquake early warning

We identify aspects of finite‐source parameterization that strongly affect the accuracy of estimated ground motion for earthquake early warning (EEW). EEW systems aim to alert users to impending shaking before it reaches them. The U.S. West Coast EEW system, ShakeAlert, currently uses two algorithms based on seismic data to characterize the earthquake’s location, magnitude, and origin time, treating it as a point or line source. From this information, ShakeAlert calculates shaking intensity and alerts locations where shaking estimates exceed a threshold. Several geodetic EEW algorithms under development would provide 3D finite‐fault information. We investigate conditions under which this information produces sufficiently better intensity estimates to potentially improve alerting. Using scenario crustal and subduction interface sources, we (1) identify the most influential source geometry parameters for an EEW algorithm’s shaking forecast, and (2) assess the intensity alert thresholds and magnitude ranges for which more detailed source characterization affects alert accuracy. We find that alert regions determined using 3D‐source representations of correct magnitude and faulting mechanism are generally more accurate than those obtained using line sources. If a line‐source representation is used and magnitude is calculated from the estimated length, then incorrect length estimates significantly degrade alert region accuracy. In detail, the value of 3D‐source characterization depends on the user’s chosen alert threshold, tectonic regime, and faulting style. For the suite of source models we tested, the error in shaking intensity introduced by incorrect geometry could reach levels comparable to the intrinsic uncertainty in ground‐motion calculations (e.g., 0.5–1.3 modified Mercalli intensity [MMI] units for MMI 4.5) but, especially for crustal sources, was often less. For subduction interface sources, 3D representations substantially improved alert area accuracy compared to line sources, and incorrect geometry parameters were more likely to cause error in calculated shaking intensity that exceeded uncertainties.

Bulletin of the Seismological Society of America↗
USGS2025

Deterministic physics-based earthquake sequence simulators match empirical ground-motion models and enable extrapolation to data poor regimes: Application to multifault multimechanism ruptures

We use the deterministic earthquake simulator RSQSim to generate complex sequences of ruptures on fault systems used for hazard assessment. We show that the source motions combined with a wave propagation code create surface ground motions that fall within the range of epistemic uncertainties for the Next Generation Attenuation‐West2 set of empirical models. We show the model is well calibrated where there are good data constraints, and has good correspondence in regions with fewer data constraints. We show magnitude, distance, and mechanism dependence all arising naturally from the same underlying friction. The deterministic physics‐based approach provides an opportunity for better understanding the physical origins of ground motions. For example, we find that reduced stress drops in shallow layers relative to constant stress drop with depth lead to peak ground velocities in the near field that better match empirical models. The simulators may also provide better extrapolations into regimes that are poorly empirically constrained by data because physics, rather than surface shaking data parameterizations, is underlying the extrapolations. Having shown the model is credible, we apply it to a problem where observations are lacking. We examine the case of crustal faults above a shallow subduction interface seen to break coseismically in simulations of the New Zealand fault system. These types of events were left out of consideration in the most recent New Zealand national seismic hazard model due to the modeling complexity and lack of observational data to constrain ground‐motion models (GMMs). Here, we show that in the model, by breaking up the coseismic crustal and interface rupturing fault motions into two separate subevents, and then recombining the resulting ground‐motion measures in a square‐root‐of‐sum‐of‐squares incoherent manner, we reproduce well the ground‐motion measures from the full event rupture. This provides a new method for extrapolating GMMs to more complex multifault ruptures.

Seismological Research Letters↗
USGS1995

Near real-time monitoring of seismic events and status of portable digital recorders using satellite telemetry

Near real-time monitoring of seismic events and status of portable 16-bit digital recorders has been established for arrays near Parkfield, Mammoth Lakes, and San Francisco, California. This monitoring system provides near real-time seismic event identification (rough location and magnitude) and a cost-effective means to maintain arrays at near 100% operational level. Principal objectives in the design of this system have been portability and low-cost telemetry. The system has been developed to use portable digital seismic recorders ( GEOS —General Earthquake Observation System) and portable data collection platforms ( DCP 's) for the Geostationary Operational Environmental Satellite (GEOS) telemetry system. Data are transferred asynchronously from the GEOS seismic system through a microprocessor-controlled interface every 10 min. The interface stores, determines priority, converts, and synchronously transfers these data to a Sutron Corp. model 8004 DCP for transmission through the GEOS satellite telemetry system. Event parameters include trigger time, peak amplitude, time of peak amplitude, and event duration. Instrument configuration parameters, transmitted at system start-up time and every 24 hr, include recording parameters, trigger parameters, GEOS software version, clock reference, and location parameter. Instrument status includes battery voltage, number of events, and percentage of tape usage. These data are transmitted as appropriate to the U.S. Geological Survey satellite downlink and computers located in Menlo Park, California, where they are processed and displayed.

Bulletin of the Seismological Society of America↗
USGS2007

Human influence on California fire regimes

Periodic wildfire maintains the integrity and species composition of many ecosystems, including the mediterranean-climate shrublands of California. However, human activities alter natural fire regimes, which can lead to cascading ecological effects. Increased human ignitions at the wildland-urban interface (WUI) have recently gained attention, but fire activity and risk are typically estimated using only biophysical variables. Our goal was to determine how humans influence fire in California and to examine whether this influence was linear, by relating contemporary (2000) and historic (1960-2000) fire data to both human and biophysical variables. Data for the human variables included fine-resolution maps of the WUI produced using housing density and land cover data. Interface WUI, where development abuts wildland vegetation, was differentiated from intermix WUI, where development intermingles with wildland vegetation. Additional explanatory variables included distance to WUI, population density, road density, vegetation type, and ecoregion. All data were summarized at the county level and analyzed using bivariate and multiple regression methods. We found highly significant relationships between humans and fire on the contemporary landscape, and our models explained fire frequency (R2 = 0.72) better than area burned (R2 = 0.50). Population density, intermix WUI, and distance to WUI explained the most variability in fire frequency, suggesting that the spatial pattern of development may be an important variable to consider when estimating fire risk. We found nonlinear effects such that fire frequency and area burned were highest at intermediate levels of human activity, but declined beyond certain thresholds. Human activities also explained change in fire frequency and area burned (1960-2000), but our models had greater explanatory power during the years 1960-1980, when there was more dramatic change in fire frequency. Understanding wildfire as a function of the spatial arrangement of ignitions and fuels on the landscape, in addition to nonlinear relationships, will be important to fire managers and conservation planners because fire risk may be related to specific levels of housing density that can be accounted for in land use planning. With more fires occurring in close proximity to human infrastructure, there may also be devastating ecological impacts if development continues to grow farther into wildland vegetation. ?? 2007 by the Ecological Society of America.

Ecological Applications↗
USGS2006

Source apportionment modeling of volatile organic compounds in streams

It often is of interest to understand the relative importance of the different sources contributing to the concentration cw of a contaminant in a stream; the portions related to sources 1, 2, 3, etc. are denoted cw,1, cw,2, cw,3, etc. Like c w, 'he fractions ??1, = cw,1/c w, ??2 = cw,2/cw, ??3 = cw,3/cw, etc. depend on location and time. Volatile organic compounds (VOCs) can undergo absorption from the atmosphere into stream water or loss from stream water to the atmosphere, causing complexities affecting the source apportionment (SA) of VOCs in streams. Two SA rules are elaborated. Rule 1: VOC entering a stream across the air/water interface exclusively is assigned to the atmospheric portion of cw. Rule 2: VOC loss by volatilization, flow loss to groundwater, in-stream degradation, etc. is distributed over cw,1 cw,2, c w,3, etc. in proportion to their corresponding ?? values. How the two SA rules are applied, as well as the nature of the SA output for a given case, will depend on whether transport across the air/water interface is handled using the net flux F convention or using the individual fluxes J convention. Four hypothetical stream cases involving acetone, methyl-tert-butyl ether (MTBE), benzene, chloroform, and perchloroethylene (PCE) are considered. Acetone and MTBE are sufficiently water soluble from air for a domestic atmospheric source to be capable of yielding cw values approaching the common water quality guideline range of 1 to 10 ??g/L. For most other VOCs, such levels cause net outgassing (F > 0). When F > 0 in a given section of stream, in the net flux convention, all of the ??j, for the compound remain unchanged over that section while cw decreases. A characteristic time ??d can be calculated to predict when there will be differences between SA results obtained by the net flux convention versus the individual fluxes convention. Source apportionment modeling provides the framework necessary for comparing different strategies for mitigating contamination at points of interest along a stream. ?? 2006 SETAC.

Environmental Toxicology and Chemistry↗
USGS2021

Improving ESRI ArcGIS performance of coastal and seafloor analysis with the Python multiprocessing module

Coastal research frequently involves the use of a GIS to analyze large areas for changes in response to major weather events, human action, and other factors. The GIS workflows used to conduct these analyses can be complex and sometimes require multiple days to complete. Long runtimes often exist even on modern high-powered workstations if the GIS software does not use parallel computing techniques, which prevents it from fully utilizing the capabilities of multicore processors. If a GIS application supports a programming interface that allows geoprocessing tools to be called from an external program, then GIS workflows can use parallel functionality embedded in that programming language to divide the load of a large workflow among multiple child processes. In ArcMap and ArcGIS Pro, this technique can be implemented by using the Python programming interface and the multiprocessing module in Python to run geoprocessing tools in child processes. This method was used in the Seafloor Elevation Change Analysis Tool (SECAT), a Python script written for ArcMap and ArcGIS Pro that calculates changes in seafloor elevation over time using two different digital elevation models. Running SECAT with between one and eight child processes on two different datasets improved execution times by at least a factor of 2.4. These results demonstrate that using the Python multiprocessing module can significantly accelerate a variety of time-consuming workflows.

Journal of Coastal Research↗
USGS1959

Geology and uranium-vanadium deposits of the slick rock district, San Miguel and Dolores counties, Colorado

Sedimentary rocks known in the Slick Rock district in southwestern Colorado range in age from Devonian (?) to Cretaceous, and aggregate about 13,000 feet in maximum thickness. Important uranium-vanadium production has come from deposits in the Salt Wash member of the Morrison formation of Late Jurassic age. The sedimentary rocks are gently folded in the Dolores and Glade anticlines and the Disappointment syncline, and are cut by the Dolores fault zone in the north part of the district and by the Glade fault zone in the south part of the district . Principal fracture sets are oriented approximately parallel to the major faults. Detrital hematite, magnetite, and ilmenite in rocks of the Morrison formation not affected by epigenetic alteration contain appreciable amounts of several of the elements found in the ore deposits . Epigenetic alteration processes have bleached large volumes of rock and largely destroyed these minerals. Such alteration is spatially associated with the Dolores fault zone. Most of the known ore deposits are in the north part of the Slick Rock district in a belt called the Dolores ore zone. The zone lies along the Dolores fault zone but is wider than the fault zone. All known deposits are associated with abundant carbonaceous plant material. Uranium-vanadium deposits in the district are chiefly tabular to lenticular and are roughly parallel to the sedimentary bedding. Some ore bodies, however, are narrow, elongate, and curve sharply across bedding; these bodies have been called "rolls" by the miners. Mineral zoning is evident in some roll bodies; carbonates, goethite (altered from pyrite), selenides, and sulfides are commonly found in concentric layers at the concave edge of rolls. This zoning, and the relationship of roll ore bodies to sedimentary structures and lithology, suggest that ore was deposited at an interface between two solutions, possibly cool connate water and a warmer ore solution. On a district scale, copper and lead are distinctly most abundant in the ore deposits within and immediately adjacent to the Dolores fault zone, and less abundant in deposits toward the edge of the zone. Uranium-vanadium deposits in the district occur only in sandstone that is considered to be epigenetically altered, and the most extensive epigenetic changes have occurred close to ore bodies. It is concluded that ground water, heated and set into circulation near the end of Cretaceous time by igneous intrusions in the La Sal and other centers on the Colorado Plateau, picked up elements from sedimentary rocks where they had been faulted and fractured, and deposited the elements at solution interfaces where accumulations of carbonaceous material provided favorable chemical conditions for precipitation.

Colorado, Utah↗
USGS1990

Geochemistry of vanadium in an epigenetic, sandstone-hosted vanadium-uranium deposit, Henry Basin, Utah

The epigenetic Tony M vanadium-uranium orebody in south-central Utah is hosted in fluvial sandstones of the Morrison Formation (Upper Jurassic). Although the deposit is mined for uranium, vanadium has a higher average abundance in the ore. Thus, the geochemistry of vanadium in the orebody was studied to characterize ore-forming processes within the inferred ground-water flow regime. Measurements of the relative amounts of V (super +3) and V (super +4) in ore minerals show that V (super +3) is more abundant. Thermodynamic calculations show that vanadium was more likely transported to the site of mineralization as V (super +4) . The ore formed as V (super +4) was reduced by hydrogen sulfide, followed by hydrolysis and precipitation of V (super +3) in oxide minerals (e.g., montroseitc or paramontroseite) or chlorite. Uranium was transported as uranyl ion (U (super +6) ), or some complex thereof, and reduced by hydrogen sulfide, forming coffinite. Detrital organic matter in the rocks served as the carbon source for sulfate-reducing bacteria. It was this bacteriogenic H 2 S that reduced the metals in the mineralization process.Possible sources for the V and U in this deposit have been identified previously. Vanadium most likely was derived from the dissolution of iron-titanium oxides, which liberated Fe as well as V. A zone of titanium-rich remnants is observed updip and up the hydrologic gradient from the deposit (M. Goldhaber and R. L. Reynolds, unpub. data). Uranium probably was derived from the overlying Brushy Basin Member of the Morrison Formation (Northrop, 1982). A preliminary age date for the deposit of 115 Ma (K. Ludwig, 1986, pers. commun.) indicates that the ore formed after deposition of the Brushy Basin Member.Previous studies have shown that the ore formed at the density-stratified interface between a basinal brine and dilute meteoric water. The mineralization processes described above occurred within the mixing zone between these two fluids. Stable isotope analyses of ore-stage dolomite show a progressively heavier carbon and oxygen isotope signature with increasing depth through an ore horizon, consistent with the two-solution interface model.

Economic Geology↗