SEARCH · Geology Reports

Results for “The Earth Observer”

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

SourceResearch connectionFrom yearThrough year

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

At least 1,171 records · Page 65Linked to original sources

USGS1993

Differentiation of volcanic ash-fall and water-borne detrital layers in the Eocene Senakin coal bed, Tanjung Formation, Indonesia

The Sangsang deposit of the Eocene Senakin coal bed, Tanjung Formation, southeastern Kalimantan, Indonesia, contains 11 layers, which are thin (<5 cm) and high in ash (> 70%). These layers are characterized by their pelitic macroscopic texture. Examination of eight of the layers by scanning-electron microscopy, energy-dispersive X-ray, and X-ray diffraction analyses show that they are composed primarily of fairly well-crystallized kaolinite, much of which is vermicular. Accessory minerals include abundant Ti oxide, rare-earth element-rich Ca and A1 phosphates, quartz that luminescences in the blue color range, and euhedral to subhedral pyroxene, hornblende, zircon, and sanidine. Although this mineral suite is suggestive of volcanic ash-fall material, only the four pelitic layers in the middle of the bed are thought to be solely derived from volcanic ash-falls on the basis of diagnostic minerals, replaced glass shards, and lithostratigraphic relationships observed in core and outcrop. The three uppermost pelitic layers contain octahedral chromites, some quartz grains that luminesce in teh orange color range, and some quartz grains that contain two-phase fluid inclusions. These layers are interpreted to be derived from a combination of volcanic ash-fall material and hydrologic transport of volcaniclastic sediment. In contrast, the lowermost pelitic layer, which contains large, rounded FeMg-rich chromites, is thought to have been dominantly deposited by water. The source of the volcanic ash-fall material may have been middle Tertiary volcanism related to plate tectonic activity between Kalimantan and Sulawesi. The volcanic ash was deposited in sufficient amounts to be preserved as layers within the coal only in the northern portions of the Senakin region: the southern coal beds in the region do not contain pelitic layers.

southeastern Kalimantan↗

USGS1989

Voyager 2 at Neptune: Imaging science results

Voyager 2 images of Neptune reveal a windy planet characterized by bright clouds of methane ice suspended in an exceptionally clear atmosphere above a lower deck of hydrogen sulfide or ammonia ices. Neptune's atmosphere is dominated by a large anticyclonic storm system that has been named the Great Dark Spot (GDS). About the same size as Earth in extent, the GDS bears both many similarities and some differences to the Great Red Spot of Jupiter. Neptune's zonal wind profile is remarkably similar to that of Uranus. Neptune has three major rings at radii of 42,000, 53,000, and 63,000 kilometers. The outer ring contains three higher density arc-like segments that were apparently responsible for most of the ground-based occultation events observed during the current decade. Like the rings of Uranus, the Neptune rings are composed of very dark material; unlike that of Uranus, the Neptune system is very dusty. Six new regular satellites were found, with dark surfaces and radii ranging from 200 to 25 kilometers. All lie inside the orbit of Triton and the inner four are located within the ring system. Triton is seen to be a differentiated body, with a radius of 1350 kilometers and a density of 2.1 grams per cubic centimeter; it exhibits clear evidence of early episodes of surface melting. A now rigid crust of what is probably water ice is overlain with a brilliant coating of nitrogen frost, slightly darkened and reddened with organic polymer material. Streaks of organic polymer suggest seasonal winds strong enough to move particles of micrometer size or larger, once they become airborne. At least two active plumes were seen, carrying dark material 8 kilometers above the surface before being transported downstream by high level winds. The plumes may be driven by solar heating and the subsequent violent vaporization of subsurface nitrogen.

Science↗

USGS2022

Flood lavas on Earth, Io and Mars

Flood lavas are major geological features on all the major rocky planetary bodies. They provide important insight into the dynamics and chemistry of the interior of these bodies. On the Earth, they appear to be associated with major and mass extinction events. It is therefore not surprising that there has been significant research on flood lavas in recent years. Initial models suggested eruption durations of days and volumetric fluxes of order 10 7 m 3 s −1 with flows moving as turbulent floods. However, our understanding of how lava flows can be emplaced under an insulating crust was revolutionized by the observations of actively inflating pahoehoe flows in Hawaii. These new ideas led to the hypothesis that flood lavas were emplaced over many years with eruption rates of the order of 10 4 m 3 s −1 . The field evidence indicates that flood lava flows in the Columbia River Basalts, Deccan Traps, Etendeka lavas, and the Kerguelen Plateau were emplaced as inflated pahoehoe sheet flows. This was reinforced by the observation of active lava flows of ≥100 km length on Io being formed as tube-fed flows fed by moderate eruption rates (10 2 –10 3 m 3 s −1 ). More recently it has been found that some flood lavas are also emplaced in a more rapid manner. New high-resolution images from Mars revealed ‘platy–ridged’ flood lava flows, named after the large rafted plates and ridges formed by compression of the flow top. A search for appropriate terrestrial analogues found an excellent example in Iceland: the 1783–1784 Laki Flow Field. The brecciated Laki flow top consists of pieces of pahoehoe, not aa clinker, leading us to call this ‘rubbly pahoehoe’. Similar flows have been found in the Columbia River Basalts and the Kerguelen Plateau. We hypothesize that these flows form with a thick, insulating, but mobile crust, which is disrupted when surges in the erupted flux are too large to maintain the normal pahoehoe mode of emplacement. Flood lavas emplaced in this manner could have (intermittently) reached effusion rates of the order of 10 6 m 3 s −1 .

Journal of the Geological Society↗

USGS2001

Mapping crustal heterogeneity using Lg propagation efficiency throughout the Middle East, Mediterranean, Southern Europe and Northern Africa

In this paper we describe a technique for mapping the lateral variation of Lg characteristics such as Lg blockage, efficient Lg propagation, and regions of very high attenuation in the Middle East, North Africa, Europe and the Mediterranean regions. Lg is used in a variety of seismological applications from magnitude estimation to identification of nuclear explosions for monitoring compliance with the Comprehensive Nuclear-Test-Ban Treaty (CTBT). These applications can give significantly biased results if the Lg phase is reduced or blocked by discontinuous structure or thin crust. Mapping these structures using quantitative techniques for determining Lg amplitude attenuation can break down when the phase is below background noise. In such cases Lg blockage and inefficient propagation zones are often mapped out by hand. With our approach, we attempt to visually simplify this information by imaging crustal structure anomalies that significantly diminish the amplitude of Lg. The visualization of such anomalies is achieved by defining a grid of cells that covers the entire region of interest. We trace Lg rays for each event/ station pair, which is simply the great circle path, and attribute to each cell a value equal to the maximum value of the Lg/P-coda amplitude ratio for all paths traversing that particular cell. The resulting map, from this empirical approach, is easily interpreted in terms of crustal structure and can successfully image small blockage features often missed by analysis of raypaths alone. This map can then be used to screen out events with blocked Lg prior to performing Q tomography, and to avoid using Lg-based methods of event identification for the CTBT in regions where they cannot work. For this study we applied our technique to one of the most tectonically complex regions on the earth. Nearly 9000 earthquake/station raypaths, traversing the vast region comprised of the Middle East, Mediterranean, Southern Europe and Northern Africa, have been analyzed. We measured the amplitude of Lg relative to the P-coda and mapped the lateral variation of Lg propagation efficiency. With the relatively dense coverage provided by the numerous crossing paths we are able to map out the pattern of crustal heterogeneity that gives rise to the observed character of Lg propagation. We observe that the propagation characteristics of Lg within the region of interest are very complicated but are readily correlated with the different tectonic environments within the region. For example, clear strong Lg arrivals are observed for paths crossing the stable continental interiors of Northern Africa and the Arabian Shield. In contrast, weakened to absent Lg is observed for paths crossing much of the Middle East, and Lg is absent for paths traversing the Mediterranean. Regions that block Lg transmission within the Middle East are very localized and include the Caspian Sea, the Iranian Plateau and the Red Sea. Resolution is variable throughout the region and strongly depends on the distribution of seismicity and recording stations. Lg propagation is best resolved within the Middle East where regions of crustal heterogeneity on the order of 100 km are imaged (e.g., South Caspian Sea and Red Sea). Crustal heterogeneity is resolvable but is poorest in seismically quiescent Northern Africa.

Pure and Applied Geophysics↗

USGS2014

Maximum magnitude earthquakes induced by fluid injection

Analysis of numerous case histories of earthquake sequences induced by fluid injection at depth reveals that the maximum magnitude appears to be limited according to the total volume of fluid injected. Similarly, the maximum seismic moment seems to have an upper bound proportional to the total volume of injected fluid. Activities involving fluid injection include (1) hydraulic fracturing of shale formations or coal seams to extract gas and oil, (2) disposal of wastewater from these gas and oil activities by injection into deep aquifers, and (3) the development of enhanced geothermal systems by injecting water into hot, low-permeability rock. Of these three operations, wastewater disposal is observed to be associated with the largest earthquakes, with maximum magnitudes sometimes exceeding 5. To estimate the maximum earthquake that could be induced by a given fluid injection project, the rock mass is assumed to be fully saturated, brittle, to respond to injection with a sequence of earthquakes localized to the region weakened by the pore pressure increase of the injection operation and to have a Gutenberg-Richter magnitude distribution with a b value of 1. If these assumptions correctly describe the circumstances of the largest earthquake, then the maximum seismic moment is limited to the volume of injected liquid times the modulus of rigidity. Observations from the available case histories of earthquakes induced by fluid injection are consistent with this bound on seismic moment. In view of the uncertainties in this analysis, however, this should not be regarded as an absolute physical limit.

Journal of Geophysical Research B: Solid Earth↗

USGS2018

Factors controlling landslide frequency-area distributions

A power‐law relation for the frequency–area distribution (FAD) of medium and large landslides (e.g. tens to millions of square meters) has been observed by numerous authors. But the FAD of small landslides diverges from the power‐law distribution, with a rollover point below which frequencies decrease for smaller landslides. Some studies conclude that this divergence is an artifact of unmapped small landslides due to lack of spatial or temporal resolution; others posit that it is caused by the change in the underlying failure process. An explanation for this dilemma is essential both to evaluate the factors controlling FADs of landslides and power‐law scaling, which is a crucial factor regarding both landscape evolution and landslide hazard assessment. This study examines the FADs of 45 earthquake‐induced landslide inventories from around the world in the context of the proposed explanations. We show that each inventory probably involves some combination of the proposed explanations, though not all explanations contribute to each case. We propose an alternative explanation to understand the reason for the divergence from a power‐law. We suggest that the geometry of a landslide at the time of mapping reflects not just one single movement but many, including the propagation of numerous smaller landslides before and after the main failure. Because only the resulting combination of these landslides can be observed due to a lack of temporal resolution, many smaller landslides are not taken into account in the inventory. This reveals that the divergence from the power‐law is not necessarily attributed to the incompleteness of an inventory. This conceptual model will need to be validated by ongoing observation and analysis. Also, we show that because of the subjectivity of mapping procedures, the total number of landslides and total landslide areas in inventories differ significantly, and therefore the shapes of FADs also differ considerably.

Earth Surface Processes and Landforms↗

USGS2002

Aftershock triggering by complete Coulomb stress changes

We examine the correlation between seismicity rate change following the 1992, M 7.3, Landers, California, earthquake and characteristics of the complete Coulomb failure stress (CFS) changes (ΔCFS( t )) that this earthquake generated. At close distances the time-varying “dynamic” portion of the stress change depends on how the rupture develops temporally and spatially and arises from radiated seismic waves and from permanent coseismic fault displacement. The permanent “static” portion (ΔCFS) depends only on the final coseismic displacement. ΔCFS diminishes much more rapidly with distance than the transient, dynamic stress changes. A common interpretation of the strong correlation between ΔCFS and aftershocks is that load changes can advance or delay failure. Stress changes may also promote failure by physically altering properties of the fault or its environs. Because it is transient, ΔCFS( t ) can alter the failure rate only by the latter means. We calculate both ΔCFS and the maximum positive value of ΔCFS( t ) (peak ΔCFS( t )) using a reflectivity program. Input parameters are constrained by modeling Landers displacement seismograms. We quantify the correlation between maps of seismicity rate changes and maps of modeled ΔCFS and peak ΔCFS( t ) and find agreement for both models. However, rupture directivity, which does not affect ΔCFS, creates larger peak ΔCFS( t ) values northwest of the main shock. This asymmetry is also observed in seismicity rate changes but not in ΔCFS. This result implies that dynamic stress changes are as effective as static stress changes in triggering aftershocks and may trigger earthquakes long after the waves have passed.

Journal of Geophysical Research B: Solid Earth↗

USGS2004

Stress interaction between subduction earthquakes and forearc strike-slip faults: Modeling and application to the northern Caribbean plate boundary

Strike-slip faults in the forearc region of a subduction zone often present significant seismic hazard because of their proximity to population centers. We explore the interaction between thrust events on the subduction interface and strike-slip faults within the forearc region using three-dimensional models of static Coulomb stress change. Model results reveal that subduction earthquakes with slip vectors subparallel to the trench axis enhance the Coulomb stress on strike-slip faults adjacent to the trench but reduce the stress on faults farther back in the forearc region. In contrast, subduction events with slip vectors perpendicular to the trench axis enhance the Coulomb stress on strike-slip faults farther back in the forearc, while reducing the stress adjacent to the trench. A significant contribution to Coulomb stress increase on strike-slip faults in the back region of the forearc comes from "unclamping" of the fault, i.e., reduction in normal stress due to thrust motion on the subduction interface. We argue that although Coulomb stress changes from individual subduction earthquakes are ephemeral, their cumulative effects on the pattern of lithosphere deformation in the forearc region are significant. We use the Coulomb stress models to explain the contrasting deformation pattern between two adjacent segments of the Caribbean subduction zone. Subduction earthquakes with slip vectors nearly perpendicular to the Caribbean trench axis is dominant in the Hispaniola segment, where the strike-slip faults are more than 60 km inland from the trench. In contrast, subduction slip motion is nearly parallel to the Caribbean trench axis along the Puerto Rico segment, where the strike-slip fault is less than 15 km from the trench. This observed jump from a strike-slip fault close to the trench axis in the Puerto Rico segment to the inland faults in Hispaniola is explained by different distributions of Coulomb stress in the forearc region of the two segments, as a result of the change from the nearly trench parallel slip on the Puerto Rico subduction interface to the more perpendicular subduction slip beneath Hispaniola. The observations and modeling suggest that subduction-induced strike-slip seismic hazard to Puerto Rico may be smaller than previously assumed but the hazard to Hispaniola remains high. Copyright 2004 by the American Geophysical Union.

Puerto Rico↗

USGS2007

Sorted bed forms as self-organized patterns: 2. complex forcing scenarios

We employ a numerical model to study the development of sorted bed forms under a variety of hydrodynamic and sedimentary conditions. Results indicate that increased variability in wave height decreases the growth rate of the features and can potentially give rise to complicated, a priori unpredictable, behavior. This happens because the system responds to a change in wave characteristics by attempting to self-organize into a patterned seabed of different geometry and spacing. The new wavelength might not have enough time to emerge before a new change in wave characteristics occurs, leading to less regular seabed configurations. The new seabed configuration is also highly dependent on the preexisting morphology, which further limits the possibility of predicting future behavior. For the same reasons, variability in the mean current magnitude and direction slows down the growth of features and causes patterns to develop that differ from classical sorted bed forms. Spatial variability in grain size distribution and different types of net sediment aggradation/degradation can also result in the development of sorted bed forms characterized by a less regular shape. Numerical simulations qualitatively agree with observed geometry (spacing and height) of sorted bed forms. Also in agreement with observations is that at shallower depths, sorted bed forms are more likely to be affected by changes in the forcing conditions, which might also explain why, in shallow waters, sorted bed forms are described as ephemeral features. Finally, simulations indicate that the different sorted bed form shapes and patterns observed in the field might not necessarily be related to diverse physical mechanisms. Instead, variations in sorted bed form characteristics may result from variations in local hydrodynamic and/or sedimentary conditions.

Journal of Geophysical Research F: Earth Surface↗

USGS1997

Permeability and porosity of the Illinois UPH 3 drillhole granite and a comparison with other deep drillhole rocks

Permeability, porosity, and volumetric strain measurements were conducted on granite cores obtained at depths of 0.7 to 1.6 km from the Illinois UPH 3 drillhole at effective confining pressures from 5 to 100 MPa. Initial permeabilities were in the range of 10-17 to 10-19 m2 and dropped rapidly with applied pressure to values between 10-20 and 10-24 m2 at 100 MPa, typical of other deep granite core samples. These values are several decades lower than equivalent weathered surface granites at comparable effective confining pressures, where weathering products in cracks and pores inhibit crack closure with applied pressure. Permeabilities of the Illinois cores were inversely related to sample depth, suggesting that stress relief and thermal microfractures induced during core retrieval dominated the fluid flow. Thus these samples provide an upper bound on in situ matrix permeability values. A comparison of core permeability from UPH 3 and other deep drillholes shows that stress relief damage can often dominate laboratory permeability measurements. We conclude that it may be difficult to make meaningful estimates of in situ permeability based on either borehole samples (possible damage during retrieval) or surface-derived analogs (altered by weathering). Volumetric strain determined from porosity measurements was compared with differential strain analysis (DSA) data reported by other investigators on samples from the same depths in the drillhole. Our strain measurements (0.002 to 0.005 at 100 MPa) were nearly twice as large as the DSA values, probably because of the crack-enhancing effects of fluids present in our samples that are absent in the dry DSA cores, as well as other time-dependent deformation effects. This difference in observed strain magnitudes between the two measurement methods may be an important consideration if strain and/or porosity data from deep core samples are used in models of stress, fluid circulation, and excess fluid pressure generation in the midcrust.

Journal of Geophysical Research B: Solid Earth↗

USGS2012

Acceleration spectra for subduction zone earthquakes

We estimate the source spectra of shallow earthquakes from digital recordings of teleseismic P wave groups, that is, P + p P + s P , by making frequency dependent corrections for the attenuation and for the interference of the free surface. The correction for the interference of the free surface assumes that the earthquake radiates energy from a range of depths. We apply this spectral analysis to a set of 12 subduction zone earthquakes which range in size from M S = 6.2 to 8.1, obtaining corrected P wave acceleration spectra on the frequency band from 0.01 to 2.0 Hz. Seismic moment estimates from surface waves and normal modes are used to extend these P wave spectra to the frequency band from 0.001 to 0.01 Hz. The acceleration spectra of moderate subduction zone earthquakes, that is, earthquakes whose seismic moments are less than 10 27 dyn cm, exhibit ω-square or Brune-type spectra, while the acceleration spectra of large subduction zone earthquakes, that is, earthquakes whose seismic moments are greater than 10 27 dyn cm, exhibit intermediate slopes where ü(ω) ∝ ω 5/4 for frequencies from 0.005 to 0.05 Hz. For this set of earthquakes, spectral shape appears to be a discontinuous function of seismic moment. Using reasonable assumptions for the phase characteristics, we transform the spectral shape observed for large earthquakes into the time domain to fit Ekström's (1987) moment rate functions for the M S = 8.1 Michoacan earthquake of September 19, 1985, and the M S = 7.6 Michoacan aftershock of September 21, 1985.

Journal of Geophysical Research Solid Earth↗

USGS1996

The Border Ranges fault system in Glacier Bay National Park, Alaska: Evidence for major early Cenozoic dextral strike-slip motion

The Border Ranges fault system of southern Alaska, the fundamental break between the arc basement and the forearc accretionary complex, is the boundary between the Peninsular–Alexander–Wrangellia terrane and the Chugach terrane. The fault system separates crystalline rocks of the Alexander terrane from metamorphic rocks of the Chugach terrane in Glacier Bay National Park. Mylonitic rocks in the zone record abundant evidence for dextral strike-slip motion along north-northwest-striking subvertical surfaces. Geochronologic data together with regional correlations of Chugach terrane rocks involved in the deformation constrain this movement between latest Cretaceous and Early Eocene (~50 Ma). These findings are in agreement with studies to the northwest and southeast along the Border Ranges fault system which show dextral strike-slip motion occurring between 58 and 50 Ma. Correlations between Glacier Bay plutons and rocks of similar ages elsewhere along the Border Ranges fault system suggest that as much as 700 km of dextral motion may have been accommodated by this structure. These observations are consistent with oblique convergence of the Kula plate during early Cenozoic and forearc slivering above an ancient subduction zone following late Mesozoic accretion of the Peninsular–Alexander–Wrangellia terrane to North America.

Alaska↗

USGS2019

Hydroacoustic, meteorologic and seismic observations of the 2016 Nansen Ice Shelf calving event and iceberg formation

On 7 April 2016, the Nansen Ice Shelf (NIS) front calved into two icebergs, the first large-scale calving event in >30 years. Three hydrophone moorings were deployed seaward of the NIS in December 2015 and over the following months recorded hundreds of short duration, broadband (10-400 Hz) cryogenic signals, likely caused by fracturing of the ice-shelf. The majority of these icequakes occur between January and early March 2016, several weeks prior to the calving observed by satellite on 7 April. Barometric pressure and wind speed records show the day the icebergs drifted from the NIS coincided with the largest low-pressure storm system recorded in the previous 7 months. A nearby seismic station also shows an increase in low-frequency energy, harmonic tremor, and microseisms on 7 April. Our interpretation is the northern segment of the NIS leading edge broke free during mid-January to February, producing high acoustic energy, but the icebergs remained stationary until the combination of a strong low-pressure system, with high winds freed the icebergs. As the unpinning of Antarctic ice shelves is not a well-documented process, our observations show that storm systems may play an under-appreciated role in Antarctic ice shelf break-up.

Frontiers in Earth Science↗

USGS2007

Statistical analysis of water-quality data containing multiple detection limits II: S-language software for nonparametric distribution modeling and hypothesis testing

Analysis of low concentrations of trace contaminants in environmental media often results in left-censored data that are below some limit of analytical precision. Interpretation of values becomes complicated when there are multiple detection limits in the data-perhaps as a result of changing analytical precision over time. Parametric and semi-parametric methods, such as maximum likelihood estimation and robust regression on order statistics, can be employed to model distributions of multiply censored data and provide estimates of summary statistics. However, these methods are based on assumptions about the underlying distribution of data. Nonparametric methods provide an alternative that does not require such assumptions. A standard nonparametric method for estimating summary statistics of multiply-censored data is the Kaplan-Meier (K-M) method. This method has seen widespread usage in the medical sciences within a general framework termed "survival analysis" where it is employed with right-censored time-to-failure data. However, K-M methods are equally valid for the left-censored data common in the geosciences. Our S-language software provides an analytical framework based on K-M methods that is tailored to the needs of the earth and environmental sciences community. This includes routines for the generation of empirical cumulative distribution functions, prediction or exceedance probabilities, and related confidence limits computation. Additionally, our software contains K-M-based routines for nonparametric hypothesis testing among an unlimited number of grouping variables. A primary characteristic of K-M methods is that they do not perform extrapolation and interpolation. Thus, these routines cannot be used to model statistics beyond the observed data range or when linear interpolation is desired. For such applications, the aforementioned parametric and semi-parametric methods must be used.

Computers & Geosciences↗

USGS2022

The consequences of climate change for dryland biogeochemistry

Drylands, which cover more than 40% of Earth’s terrestrial surface, are dominant drivers of global biogeochemical cycling and home to more than one third of the human population. Climate projections predict warming, drought frequency and severity, and evaporative demand will increase in drylands at faster rates than global means. Due to extreme temperatures and high biological dependency on limited water availability, drylands are predicted to be exceptionally sensitive to climate change and, indeed, significant climate impacts are already being observed. Yet our understanding and ability to forecast climate change effects on dryland biogeochemistry and ecosystem functions lag behind many mesic systems. To improve our capacity to forecast ecosystem change, we propose focusing on the controls and consequences of two key characteristics affecting dryland biogeochemistry: i) high spatial and temporal heterogeneity in environmental conditions and ii) generalized resource scarcity. In addition to climate change, drylands are experiencing accelerating land use change. Building our understanding of dryland biogeochemistry in both intact and disturbed systems will better equip us to address the interacting effects of climate change and landscape degradation. Responding to these challenges will require a diverse, globally distributed, and interdisciplinary community of dryland experts united towards better understanding these vast and important ecosystems.

New Phytologist Foundation↗

USGS2010

Evolution of Ore Deposits and Technology Transfer Project: Isotope and Chemical Methods in Support of the U.S. Geological Survey Science Strategy, 2003-2008

Principal functions of the U.S. Geological Survey (USGS) Mineral Resources Program are providing assessments of the location, quantity, and quality of undiscovered mineral deposits, and predicting the environmental impacts of exploration and mine development. The mineral and environmental assessments of domestic deposits are used by planners and decisionmakers to improve the stewardship of public lands and public resources. Assessments of undiscovered mineral deposits on a global scale reveal the potential availability of minerals to the United States and other countries that manufacture goods imported to the United States. These resources are of fundamental relevance to national and international economic and security policy in our globalized world economy. Performing mineral and environmental assessments requires that predictions be made of the likelihood of undiscovered deposits. The predictions are based on geologic and geoenvironmental models that are constructed for the diverse types of mineral deposits from detailed descriptions of actual deposits and detailed understanding of the processes that formed them. Over the past three decades the understanding of ore-forming processes has benefited greatly from the integration of laboratory-based geochemical tools with field observations and other data sources. Under the aegis of the Evolution of Ore Deposits and Technology Transfer Project (referred to hereinafter as the Project), a 5-year effort that terminated in 2008, the Mineral Resources Program provided state-of-the-art analytical capabilities to support applications of several related geochemical tools to ore-deposit-related studies. The analytical capabilities and scientific approaches developed within the Project have wide applicability within Earth-system science. For this reason the Project Laboratories represent a valuable catalyst for interdisciplinary collaborations of the type that should be formed in the coming years for the United States to meet its natural-resources and natural-science needs. This circular presents an overview of the Project. Descriptions of the Project laboratories are given first including descriptions of the types of chemical or isotopic analyses that are made and the utility of the measurements. This is followed by summaries of select measurements that were carried out by the Project scientists. The studies are grouped by science direction. Virtually all of them were collaborations with USGS colleagues or with scientists from other governmental agencies, academia, or the private sector.

Circular↗

USGS1999

P wave crustal velocity structure in the greater Mount Rainier area from local earthquake tomography

We present results from a local earthquake tomographic imaging experiment in the greater Mount Rainier area. We inverted P wave arrival times from local earthquakes recorded at permanent and temporary Pacific Northwest Seismograph Network seismographs between 1980 and 1996. We used a method similar to that described by Lees and Crosson [1989], modified to incorporate the parameter separation method for decoupling the hypocenter and velocity problems. In the upper 7 km of the resulting model there is good correlation between velocity anomalies and surface geology. Many focal mechanisms within the St. Helens seismic zone have nodal planes parallel to the epicentral trend as well as to a north-south trending low-velocity trough, leading us to speculate that the trough represents a zone of structural weakness in which a moderate ( M 6.5–7.0) earthquake could occur. In contrast, the western Rainier seismic zone does not correlate in any simple way with anomaly patterns or focal mechanism fault planes, leading us to infer that it is less likely to experience a moderate earthquake. A ∼10 km-wide low-velocity anomaly occurs 5 to 18 km beneath the summit of Mount Rainier, which we interpret to be a signal of a region composed of hot, fractured rock with possible small amounts of melt or fluid. No systematic velocity pattern is observed in association with the southern Washington Cascades conductor. A midcrustal anomaly parallels the Olympic-Wallowa lineament as well as several other geophysical trends, indicating that it may play an important role in regional tectonics.

Journal of Geophysical Research B: Solid Earth↗

USGS2021

An eddy-resolving numerical model to study turbulent flow, sediment and bed evolution using detached eddy simulation in a lateral separation zone at the field-scale

Turbulence-resolving simulations elucidate key elements of fluid dynamics and sediment transport in fluvial environments. This research presents a feasible strategy for applying state-of-the-art computational fluid mechanics to the study of sediment transport and morphodynamic processes in lateral separation zones, which are common features in canyon rivers where massive lateral flow separation causes large-scale turbulence that controls sediment erosion and deposition. An eddy-resolving model was developed and tested at the field-scale, coupling a viscous flow and sediment transport solver using Detached Eddy Simulation techniques. A morphodynamic model was applied to the viscous flow/sediment solver to calculate erosion and deposition. A simulation of turbulence was performed at the grid resolution for a straight channel to determine the relative contributions of modeled and resolved diffusivity. The time-dependent, energetically important, correlative, non-stationary signals of the simulated quantities were captured at the lateral separation zone. Strong periodic signals featured by high amplitude were found at the separation zone, while low frequency pulsations were observed at the reattachment zone of the lateral separation zone. Interactions between the eddies and the loose bed boundaries resulted in erosion of sediment at the main channel followed by deposition at the primary eddy and eddy bars. tions elucidate key elements of fluid dynamics and sediment transport in fluvial environments. This research presents a feasible strategy for applying state-of-the-art computational fluid mechanics to the study of sediment transport and morphodynamic processes in lateral separation zones, which are common features in canyon rivers where massive lateral flow separation causes large-scale turbulence that controls sediment erosion and deposition. An eddy-resolving model was developed and tested at the field scale, coupling viscous flow and sediment transport solver using Detached Eddy Simulation (DES) techniques. A morphodynamic model was applied to the viscous flow/sediment solver to calculate erosion and deposition. A simulation of turbulence was performed at the grid resolution for a straight channel to determine the relative contributions of modeled and resolved diffusivity. The time-dependent, energetically important, correlative non-stationary signals of the simulated quantities were captured at the lateral separation zone. Strong periodic signals featured by high amplitude were found at the separation zone, while low frequency pulsations were observed at the reattachment zone of the lateral separation zone. Interactions between the eddies and the loose bed boundaries resulted in massive erosion of sediment at the main channel followed by deposition at the primary eddy and eddy bars.

Arizona↗