USGS Science⌕ Search

SEARCH · USGS Science

Results for “Integrated Environmental Assessment and Management”

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.

1,661 records · Page 76Linked to original sources
USGS2000

Natural attenuation of chlorinated volatile organic compounds in ground water at Area 6, Naval Air Station Whidbey Island, Washington

Natural attenuation is a viable alternative to pump and treat for meeting remediation objectives in the vicinity of the southern contaminant plume. The combination of historically low contaminant concentrations in ground water, a landfill cap that limits source area contributions, favorable conditions for degradation of VC, and a relatively long downgradient distance to potential receptors are all favorable for natural attenuation as a remediation alternative. Natural attenuation could effectively meet all but one remediation goal that extraction wells PW- 2, PW-4, PW-6, PW-7, PW-8 and PW-9 are currently being employed to meet. The goal of preventing migration of all VC across the site boundary could not be met by natural attenuation. Some VC would migrate south of the Navy boundary, but the potential for subsequent VC mineralization downgradient of the base and the existing institutional controls would result in minimal additional risk from using natural attenuation. In the western contaminant plume natural attenuation is not currently a viable alternative to pump and treat for meeting remediation objectives. There is a possibility that rates for reductive dechlorination of TCE and TCA could increase substantially if the plume was allowed to migrate beneath the Oak Harbor landfill, but there are not enough data to be certain of such an increase. TCA and TCE concentrations at the leading edge of the western plume need to be reduced to at least 25 and 40 µg/L, respectively, to be protective of potential downgradient receptors. Source area TCA and TCE concentrations have decreased substantially over the past ten years, and the extraction wells PW-3 and PW-5 in particular are removing a significant mass of contaminants from ground water, so natural attenuation may be a viable alternative for the western plume in the future. There would be some possible side benefits of using natural attenuation as an alternative to pump and treat in the southern contaminant plume. The first would be that the resulting decrease in the amount of treated water that would need to be recharged in the swale north of the landfill would result in less off-base migration of contamination across the western site boundary. The second benefit would be a substantial reduction in the amount of dissolved iron and manganese being extracted from the shallow aquifer and run through the treatment system. Removing that source of operation and maintenance problems would result in more effective containment and removal of contamination in the western contaminant plume. The most critical data gap identified in this evaluation is the paucity of contaminant chemistry information downgradient of the Navy boundary in the vicinity of the southern contaminant plume. Without such data, the behavior of the plume and the protectiveness of natural attenuation to downgradient receptors cannot be verified, and field attenuation rates for VC cannot be determined directly. The existing long-term monitoring plan would need to be reviewed and revised if natural attenuation is selected as a remedy for the contamination in the southern plume. In particular, additional performance monitoring wells may be required downgradient of the property.

Washington↗
USGS2025

Smectite-rich horizons in Inceptisols trigger shallow landslides in tropical granitic terranes

Puerto Rico was affected by >70,000 landslides in the wake of 2017 Hurricane Maria, and landslide prevalence was especially high in the Utuado region in the Cordillera Central. Landslide density was highest where soil parent material is granodiorite; landslide slip surfaces tended to be shallow (<60 cm), and often were mobilized rapidly and with long runout distances. This study combines field observations with soil mineralogy (bulk and clay fractions), soil geochemistry (bulk fraction), and soil strength as determined by field cone penetrometer testing (CPT) to test the hypothesis that clay-rich subsoil horizons function as slip planes when water-logged. Soil pits were excavated to depths of ∼200 cm in Ultisols on an undulating plateau and to ∼100 cm in Inceptisols on steep slopes (36-43 o ) that flank the plateau and cone penetrometer tests (CPT) were done within 2 m of the soil pit. Six pits were located adjacent to scarps from previous landslides, enabling analysis of soil profiles downward through extrapolated slip surfaces. Results from X-ray diffraction (XRD), X-ray fluorescence (XRF) and thermogravimetric analysis (TGA) indicate that soils are heterogeneous, often with subsoil horizons enriched in clay minerals and immobile elements (Al, Fe, Ti). Inceptisols on steep slopes often contain smectite-rich horizons at 30–60 cm depth that appear to function as slip surfaces; in other Inceptisols, such horizons are not present and landslide susceptibility is potentially lower. In Ultisols, soil mineralogy is dominated by kaolinite with minor halloysite, and increased kaolinite content at soil depths ≥80 cm at some sites suggests potential slip surfaces enhancing probability of landslides. The origin of clay-rich horizons appears to be (1) fractures in granodiorite that facilitate water flow and leaching, accelerating mineral dissolution during early weathering stages, and (2) smectite-rich buried soils under permeable colluvium likely deposited by a prior mass wasting event. Where clay-rich layers occur beneath more-permeable horizons, rapid infiltration then absorption of water in clay-rich subsoil horizons causes decreased shear strength and increased landslide susceptibility.

Puerto Rico↗
USGS2025

Evidence for offset of Cretaceous plutons by the Tintina fault in eastern Alaska: Implications for regional metallogeny

Cretaceous magmatism in eastern interior Alaska is voluminous, but temporally and spatially diverse – suggestive of varying sources and drivers. More than 150 new U-Pb zircon and more than 500 geochemical analyses of Cretaceous plutonic units allow for the grouping of distinct plutonic suites. Magmatism was continuous from 120-66 Ma but can be grouped into temporally distinct pulses from ca. 115-100 Ma, 100-90 Ma, and 75-66 Ma. Geochemical diversity occurs during each pulse, further distinguishing multiple suites. Diverse metallogenic epochs are strongly correlated to pluton chemistry. Mineralization is largely absent prior to 108 Ma. From 108-100 Ma, plutonism is coeval with sparse, but notable Au-quartz veins with variable Bi, As, W, and Mo. From 100-90 Ma, intrusion-related mineralization zones from Au-Cu(-Bi) and U-Th in the northwest to central Au- Bi-As-Te(-W), and Mo-W to the southeast. Porphyry style Cu-Mo(-Au) occurrences occur with the latest Cretaceous plutons emplaced from 75-66 Ma . Restoration of ~450 km of dextral movement on the Tintina fault and comparison of metallogenic and geochemical characteristics of Alaska plutons suggest 100-90 Ma plutons may be the continuation of the metallogenically significant Tombstone, Mayo, and Tungsten suites from the Yukon.

Alaska↗
USGS2026

A case study of maximum depositional ages from terrestrial sandstones near the Cretaceous—Paleogene transition, western Williston Basin, USA

We present a new Bayesian method for deriving maximum depositional ages from detrital K-feldspar using total-fusion 40 Ar/ 39 Ar data. Individual analyses weighted by K/Ca ratio, age uncertainty, and percent radiogenic argon prioritize a result to come from accurately measured volcanic sources. Results from this method applied to sandstones from the Hell Creek region of northeastern Montana show that detrital K-feldspar maximum depositional ages align closely with detrital U-Pb zircon ages as well as tephra based chronostratigraphic constraints from both 40 Ar/ 39 Ar and U-Pb data. An age model informed by maximum depositional ages and available tephra data provides an estimate of 1.74 ± 1.04 Ma for the duration of deposition for the Hell Creek Formation. Combined age model and maximum depositional age data suggest ≤ 1.84 Ma of missing time is represented by the Hell Creek and Fox Hills formational disconformity, coinciding with the final regression of the Western Interior seaway.

Montana↗
USGS2025

Grammar to graph—An approach for semantic transformation of annotations to triples

Data annotation is the process of labeling data to show the outcome that a related data model should predict. In this study, annotation data were transformed into semantic graph triples, mainly for use with the Resource Description Framework (RDF), a type of entity-relationship-attribute data model for graph databases. The transformation of annotation data to semantic graph triples provides complex linguistic meaning with data handling advantages such as reduced data storage needs, improved logical specification of relations between objects, and reusable classes and properties that support logic and inference. A grammar-based framework in graph form supports user questions and queries. The words defining approximately 334 topographic feature types compiled by the U.S. Geological Survey were tokenized as units of analysis and grouped by part of speech. Their dependency relations were identified for this study using natural language processing libraries. Dependency concepts are used as structured semantic relations among part-of-speech classes. Tokens, units equivalent to words, form instances of classes and were quantified within a tabular output format using PostgreSQL data storage software. Table data were logically aligned as triples following a mapping file and stored with an ontology file using Ontop virtual triplestore software. A grammar ontology schema for the data was synchronized to match queries whose results validated the graph’s structure. The text analysis produced 8 part-of-speech classes of content words for object representations and 4 classes of function words for operational applications. Dependency relations formed 27 ontology properties for topographic subgraph structures. Token occurrences shaped overall ontology salience and formed a lexicon of syntactic terms for subgraph objects and properties. The schema ontology of class and property population shapes formed the lexicon of English terms. SPARQL Protocol and RDF Query Language (SPARQL) was used with the lexicon to conform data to RDF guidelines. This study confirms the hypothesis that although linguistic logic varies from description logic, its approximation applies to ontology design. Property and query use case patterns extracted from the analysis support queries concerning complex topographic relations and patterns normally embedded within text definitions. The method used in this study could be applied to text forms in other domains, such as survey notes.

Scientific Investigations Report↗
USGS2024

Earthquake effects surveyed during the nineteenth century as ecological features of Chinookan tidelands

Lasting effects of a Cascadia earthquake in 1700 were documented during surveys of Chinookan tidelands near the mouth of the Columbia River between 1805 and 1868. The effects resemble estuarine consequences, near Anchorage, of the 1964 Alaska earthquake: fatal drowning of subsided meadows and forests by post-earthquake tides, rebirth of marshes and forests through post-earthquake sedimentation, and uplift. Chinookan remains of killed forests were recorded by James Graham Cooper, John J. Lowell, and Cleveland Rockwell. Cooper, attached to a railroad survey and the Smithsonian Institution, wrote of redcedar stumps and trunks standing dead in tidal marshes of Shoalwater (now Willapa) Bay. Two such snags served as bearing trees for Lowell as he platted a Shoalwater Bay township under contract with the General Land Office. Rockwell, of the US Coast Survey, flecked landward edges of tidal flats west of Astoria with symbols that evoke remains of a bygone spruce forest. The Lewis and Clark Expedition, while in that area in 1805–1806, mapped and puzzled over tideland vegetation that post-1700 succession helps explain.

Alaska, British Columbia, California, Oregon, Wash↗
USGS2025

Earthquake magnitude and source parameter estimation with a distributed acoustic sensing dataset in the Gorda subduction zone

Distributed acoustic sensing (DAS) systems offer a cost‐effective way to create large‐scale strainmeter arrays for seismological applications using fiber‐optic cables. DAS‐based strain measurements are known to be influenced by various factors, bringing into question their general reliability for accurate earthquake characterization. A 15‐km‐long DAS deployment in northern California was operational within 3 days of the 2022 M w 6.4 Ferndale earthquake and ran continuously throughout the aftershock sequence. We utilize these aftershock data to validate DAS‐based strain measurements in two ways. We first test the accuracy of DAS‐based magnitude estimates from peak dynamic strains by comparing them with magnitude and attenuation scaling relations derived independently from traditional borehole strainmeter (BSM) data. We demonstrate that DAS‐based magnitudes are comparable to BSM‐based magnitudes when corrections for variations in site response along the fiber‐optic cable are properly made. Magnitude errors are spatially correlated, potentially because of factors such as finite‐fault effects (e.g., stress drop) or more complex, unmodeled path attenuation or because of wave propagation effects in heterogeneous media. We then apply more advanced source characterization methodology to the DAS data using a time‐domain empirical Green’s function (EGF) deconvolution approach to measure details of the moment rate history. The EGF approach using DAS data depends on careful treatment of distorting factors such as anthropogenic sources of noise and optical phase wrapping but successfully isolates source spectra for moderate‐magnitude earthquakes: source spectral ratios obtained from DAS data, broadband seismometer data, and BSM data in the same region show consistent results, revealing differences in directivity and spectral shape among earthquakes. Although further research is needed to refine source‐time‐function estimation techniques for DAS data, particularly for larger magnitude events, these case studies demonstrate the clear potential of DAS for earthquake source characterization.

California↗
USGS2025

A high-resolution 3-D P-wave velocity structure of the south-central Cascadia subduction zone from wide-angle shore-crossing seismic refraction data

This study addresses a significant gap in understanding the features of the south-central Cascadia subduction zone, a region characterized by complex geologic, tectonic, and seismic transitions both offshore and onshore. Unlike other segments along this margin, this area lacks a 3-D velocity model to delineate its structural and geological features on a fine scale. To address this void, we developed a high-resolution 3-D P-wave velocity model using active source seismic data from ship-borne seismic shots recorded on temporary and permanent onshore seismic stations and ocean-bottom seismometers. Our model shows velocity variations across the region with distinct velocity-depth profiles for the Siletz, Franciscan, and Klamath terranes in the overlying plate. We identified seaward dipping high-velocity static backstops associated with the Siletz and Klamath terranes, situated near the shoreline and further inland, respectively. Regions of reduced crustal velocity are associated with crustal faults. Moreover, there is significant along-strike depth variation in the subducting slab, which is about 4 km deeper near the thick, dense Siletz terrane and becomes shallower near the predominantly less-dense Franciscan terrane. This highlights a sudden tectonic and geologic transition at the southern boundary of the Siletz terrane. Our velocity model also indicates slightly increased hydration, though still minimal, in both the oceanic crust and the upper mantle of the subducting plate compared to other parts of the margin.

British Columbia, California, Oregon, Washington↗
USGS2024

(Re)discovering the seismicity of Antarctica: A new seismic catalog for the southernmost continent

We apply a machine learning (ML) earthquake detection technique on over 21 yr of seismic data from on‐continent temporary and long‐term networks to obtain the most complete catalog of seismicity in Antarctica to date. The new catalog contains 60,006 seismic events within the Antarctic continent for 1 January 2000–1 January 2021, with estimated moment magnitudes (⁠Mw ⁠) between −1.0 and 4.5. Most detected seismicity occurs near Ross Island, large ice shelves, ice streams, ice‐covered volcanoes, or in distinct and isolated areas within the continental interior. The event locations and waveform characteristics indicate volcanic, tectonic, and cryospheric sources. The catalog shows that Antarctica is more seismically active than prior catalogs would indicate, examples include new tectonic events in East Antarctica, seismic events near and around the vicinity of David Glacier, and many thousands of events in the Mount Erebus region. This catalog provides a resource for more specific studies using other detection and analysis methods such as template matching or transfer learning to further discriminate source types and investigate diverse seismogenic processes across the continent.

Seismological Research Letters↗
USGS2025

Sand provenance boundary in the Mu Us Sandy Land of northern China

Desert dunes are often assumed to have uniform mineral compositions due to extensive mixing during lateral transport, which complicates provenance studies. The Mu Us Sandy Land in north-central China, near the East Asian summer monsoon precipitation boundary, experiences a wetter climate than most deserts. Climate wetting as a result of a warming climate, and the ‘Sand Control Project’ implemented by the Chinese government to decrease the lateral movement of sand dunes in this area provide an opportunity to study surface processes of sand production and transport. Previous studies using zircon U Pb geochronology and heavy mineral composition indicate distinct sand sources for the Mu Us Sandy Land: local basement-derived middle Yellow River sediments and recycled dried-up lacustrine sediments for the eastern part, and northeastern Tibetan Plateau-derived upper Yellow River sediments for the western part. However, zircons and heavy minerals only represent trace amounts of the bulk mineralogy within a sand dune, so broader provenance analysis targeting common minerals is essential. We focus on the area near the proposed provenance difference boundary between the western and eastern Mu Us Sandy Land, using comprehensive sampling and provenance techniques to confirm distinct provenance characteristics and delineate the provenance boundary. Our findings reveal that sand from most of the Mu Us Sandy Land originated from the erosion of local basement by the middle Yellow River and recycled local dried-up lacustrine sediments, whereas the southwestern corner and the neighboring western-central Chinese Loess Plateau received sediments from the distal northeastern Tibetan Plateau.

Mu Us Sandy Land↗
USGS2025

Airborne geophysical analysis to decipher salinization for coastal Louisiana

Coastal Louisiana is known for saltwater intrusion that threatens wetlands, aquifers, and rivers. However, the extent of saltwater intrusion is not well understood. This study develops an innovative framework with airborne electromagnetic (AEM) data to map chloride concentration distributions for wetlands in the Mississippi River deltaic plain and Chenier plain as well as for the Mississippi River Valley alluvial aquifer (MRVA) and Chicot aquifer. Moreover, the framework maps chloride concentrations along the Mississippi River and Atchafalaya River. Key components in the framework include the establishment of resistivity-to-chloride concentration transformation, 3D resistivity architecture building through geostatistics, and the employment of a lithologic model. The transformation functions correlate AEM resistivity data with porewater salinity measurements and groundwater and river chloride samples. The results show that AEM data reliably infers soil water chloride concentrations and correlates well with the distribution of various marsh types. AEM data reveals extensive saltwater presence at depth and near the coast, originating from salt domes and the Gulf of Mexico, respectively. The saltwater upconing pattern in the Chicot aquifer is likely due to excessive groundwater withdrawals. The AEM data also confirms a distinct tongue of saltwater intruding into the Atchafalaya Basin from the Gulf. The AEM data helps to identify faults that are obscured or eroded at the surface, which appear as leaky barriers in the subsurface where dramatic changes in chloride concentration are apparent. Finally, this study uses the AEM data to infer the presence of an extensive seawater wedge in the Mississippi River and Atchafalaya River.

Louisiana↗
USGS2026

Deep groundwater total dissolved solids mapping in the Dakota Group, Williston Basin, USA

Growing concern about the quantity of available freshwater around the world has led to interest in surveying groundwater total dissolved solids (TDS) below water well depths. Deep TDS has not been systematically mapped, and there is much to learn about the distribution and controls on deeper groundwater. In sedimentary basins across the United States, groundwater resources often overlie hydrocarbon resources, providing an opportunity to use borehole geophysical data collected for hydrocarbons to characterize groundwater and pore space resources. This study adapts a recently developed subsurface geostatistical and geophysical modeling approach to continuously map groundwater TDS, porosity, and temperature in the Dakota Group of the Williston Basin—an undercharacterized regional aquifer system overlying deeper hydrocarbon reservoirs. Groundwater TDS in the Dakota Group ranges from approximately 4800 to 26,900 mg/L. TDS patterns are stratified with higher TDS in the lower and upper Dakota Group, and relatively lower TDS in the middle Dakota Group. The lower TDS in the middle zone may represent a preferential regional flow path for lower-TDS meteoric recharge from the west. The alternating pattern of TDS may also be evidence of higher-TDS inflows into the Dakota Group from underlying and potentially from overlying aquifers. Porosity is lower near the center of the Williston Basin and tends to be higher to the east, which may be related to grain size distributions. The new regional TDS and porosity modeling serves as a quantitative reference for water users and provides supporting evidence for hypotheses on Dakota Group recharge.

Montana, North Dakota, South Dakota↗
USGS2026

Water table rise sustains carbon release from soils in wetland-dominated landscapes: An intact soil core study

Wetland-dominated landscapes influence carbon cycling through their potential to act as both carbon sinks and sources. Wetlands in low-relief landscapes have dynamic terrestrial-aquatic interfaces that change seasonally with variable surface water and groundwater levels. However, few studies have directly quantified dissolved organic matter (DOM) release and greenhouse gas (CO 2 , CH 4 ) fluxes from wetland soils along terrestrial-aquatic interfaces as they are seasonally re-saturated by groundwater. To estimate groundwater-mediated soil DOM and gas fluxes, we performed laboratory simulations of vertical groundwater rise on intact soil cores collected from four Delmarva bay wetlands located in the Mid-Atlantic United States. At each wetland, one core was collected from within the wetland basin and the other from the transitional zone at the basin edge. Cores were re-saturated with groundwater over 15 days and then kept fully saturated for an additional 25 days. Source groundwater, soil porewater, and exfiltrated surface water samples were collected and analyzed for pH, ORP, DOM concentration, and DOM optical indices. In both the wetland and transition zone cores, porewater DOM concentrations increased over the wet up and were sustained during prolonged saturation. Optical indices shifted from recently produced, microbial-like signatures towards aromatic, terrestrial-like signatures. Fluxes of CO 2 decreased as the duration of soil saturation increased and soil cores switched from CH 4 sinks to sources upon full soil core saturation. Results indicate that groundwater rise sustains carbon mobilization from soils in wetland-dominated landscapes, emphasizing the need to understand how climate-driven changes to groundwater dynamics may affect carbon fluxes along terrestrial-aquatic interfaces.

Delaware, Maryland↗
USGS2026

Bedrock geologic map of the Eagle Lake quadrangle, Essex County, New York

The bedrock geology of the 7.5-minute Eagle Lake quadrangle, Essex County, New York, consists of deformed and metamorphosed Mesoproterozoic gneisses of the Adirondack Highlands unconformably overlain by weakly deformed lower Paleozoic sedimentary rocks of the Champlain Valley. The Mesoproterozoic rocks occur on the eastern edge of the Adirondack Highlands and represent an extension of the Grenville Province of Laurentia. Granulite facies Mesoproterozoic paragneiss, marble, and amphibolite hosted the emplacement of an anorthosite-mangerite-charnockite-granite (AMCG) suite, now exposed mostly as orthogneiss, at approximately 1.18–1.15 giga-annum (Ga, billion years before present). The earliest of four phases of deformation (D1) predated AMCG magmatism and is characterized by gneissosity, rarely preserved F1 isoclinal folds, and migmatite in the paragneiss host rocks. A sample of hornblende quartz syenite from the AMCG suite, collected from an abandoned railroad cut on Old Furnace Road, yielded a U-Pb zircon age of 1,149±10 million years before present. D2 deformation produced a composite penetrative gneissosity, migmatite, and isoclinal F2 folds. Towards the end of D2, felsic magmatism (including the regionally extensive Lyon Mountain Granite Gneiss, abbreviated “LMG”) spread by penetrative migration as semiconcordant alkali feldspar granite sheets subparallel to S2 into the previously deformed lithologies. The LMG crystallized at approximately 1.15 to 1.14 Ga and displays synkinematic F2 folds thus constraining the time of D2 deformation. Exhumation of the Marcy anorthosite began during D3 along a mylonitic extensional detachment, as a type of core complex. Protracted D3 produced F3 folds exhibited in regional domes and basins, such as the Hammondville antiform, reactivation of the S2 foliation, partial melting, metamorphism, metasomatism, iron ore remobilization, and intrusion of magnetite-bearing pegmatite both as layer-parallel sills and crosscutting dikes. D4 created NE- and NW-trending boudinage, local high-grade ductile shear zones, and crosscutting granitic pegmatite dikes. Kilometer (km)-scale lineaments readily observed in lidar data are Ediacaran mafic dikes and Phanerozoic brittle faults. Lower Paleozoic rocks are part of the Early Cambrian to Late Ordovician great American carbonate bank on the ancient margin of Laurentia. The Potsdam Sandstone preserves the Cambrian stratigraphy in outliers above the Great Unconformity. The Paleozoic rocks are weakly folded and block faulted. Parts of the quadrangle are covered by undifferentiated glacial deposits, but much of the quadrangle contains only a variably thick, veneer of unmapped glacial till over significant areas of exposed bedrock. The map also shows waste rock piles and locations of historical mining operations. This study was undertaken to improve our understanding of the bedrock geology in the Adirondack Highlands, establish a modern framework for 1:24,000-scale bedrock geologic mapping in the Adirondack Mountains, and provide a modern context for historical mines. This Scientific Investigations Map of the Eagle Lake 7.5-minute quadrangle consists of a map sheet, an explanatory pamphlet, and a geographic information system database that includes bedrock geologic units, faults, outcrops, and structural geologic information. The map sheet includes a bedrock geologic map, a correlation of map units, a description of map units, an explanation of map symbols, and two cross sections. The explanatory pamphlet includes a discussion of the geology.

New York↗
USGS1977

A multi-channel digital telemetry system for low frequency geophysical data

An inexpensive general purpose digital telemetry system for collection of low frequency geophysical data from U.S. Geological Survey instruments (eg. tilt, strain, gravity, creep, water level, radon, magnetic field, resistivity, telluric current, temperature, etc.) has been designed and built. This system provides data for a more general interactive data acquisition, retrieval and analysis system. The field stations are self-contained, battery operated and housed in weather proof containers. Each accepts up to 15 analog data inputs in the range of -5 to +5 volts. The dynamic range is 70db. The units transmit information as FSK (Frequency Shift Keyed) tones onto either a phone line or radio link with up to 150 transmitters sharing one line. The average power consumption is 0.06 nR watts where n is the 1 number of input channels transmitted and R is the sample rate in minutes -1 . The central receiver-recorder unit accepts and decodes the FSK tones and converts, formats and records the digital data together with time information and station identification on IBM combatible magnetic tape. The digital data are also converted and recorded in analog form for visual monitoring.

Open-File Report↗
USGS2025

Pre-eruptive characteristics of “suspect” silicic magmas in Carlin-type Au-forming systems

World-class Carlin-type Au deposits hosted in sedimentary rock were formed when profuse Eocene silicic magmatism swept across northern Nevada in response to arc migration. Carlin-type Au deposits formed along with porphyry/skarn Cu-Mo-W-Au deposits, epithermal Ag-Au deposits, and distal disseminated Ag-Au deposits. But unlike these other Au-bearing deposits that have clear associations with igneous intrusions, Carlin-type ore deposits appear to have formed distant from concealed plutons, and their origin remains controversial. Despite decades of abundant geophysical, geochronological, and geochemical studies suggesting the involvement of magmas, concrete evidence for magmatic involvement is still lacking. Consequently, the involvement of contemporaneous igneous systems remains inferred based on age, proximity, and variable isotopic, geochemical, and geophysical clues. A recent synthesis of deposit models postulates that Carlin-type Au deposits are intrusion-related, but that the causative magmas reside deeper (∼6–12 km) than in typical porphyry and peripheral systems (∼3–5 km), meaning that Carlin-type deposits are perhaps more distal expressions of igneous intrusions. We investigate a collection of “suspect” magmatic systems over a ∼7 m.y. timespan (∼41–34 Ma) that are contemporaneous with and near known Carlin-type ore deposits. We report results of a multifaceted array of in situ geochemical analyses (FTIR, EMP, SHRIMP-RG, LA-ICP-MS) of quartz-hosted melt inclusions, biotite, and quartz to better characterize the pre-eruptive characteristics of these magmas. We also report results of thermobarometry and thermodynamic phase equilibria modeling to help place constraints on magmatic reservoir depths and processes. Rather than a single “flavor” of silicic magma, we observe a surprisingly broad compositional spectrum of rhyolites, with one end of the spectrum exhibiting more arc-like (I-type) characteristics and the other end displaying more post-subduction, thick-crust extensional (A-type) characteristics. This broad compositional spectrum suggests a more complex picture of silicic crustal magmatism operating over a narrow span of time during slab rollback. Despite this spectrum, magmatic systems in this study are consistently ferroan and generally peraluminous, which we interpret as an expression of the relatively elevated geotherm at the time and incorporation of variable amounts of highly peraluminous metasedimentary crustal components. The silicic magma spectrum encompasses a range of mineralization associations, including subduction-related Cu-Mo-W-Au-Ag and post-subduction, thick-crust extensional rare-metal Mo-Sn-W-F-Be-Ag-Au, consistent with the prolific and diverse array of ore deposits that formed during this time. Carlin-type Au deposition appears to be associated with nearly the entire magmatic spectrum. This apparent indifference to silicic magma “flavor” would seem to imply that if magmas are involved in Carlin-type Au deposit genesis, they perhaps do not need to be compositionally specialized and/or possibly are only relevant as heat sources driving circulation to remobilize and redistribute metals.

Nevada↗
USGS2026

The geologic history of the Chehalis Forearc Basin, Washington State, USA

The Chehalis basin is located between the Cascade arc and the Coast Range in southwest Washington State. It consists of sedimentary and volcanic rocks deposited throughout the Cenozoic and is underlain by the Siletzia terrane, a thick, large igneous province accreted in the Eocene. Here, we constrain evolution of the Chehalis basin depocenter using isochore maps derived from isostatic gravity anomalies, borehole data, and stratigraphy for several time periods: the Eocene (47.6–36.8 Ma), the Oligocene (36.8–20 Ma), and the Neogene (20–0 Ma). Our results suggest that local subsidence is driven by deformation on faults that bound and intersect the basin. We see northward depocenter migration, interpreted as a shift in deformation from the northwest striking Cedar Creek fault to the west striking Doty fault. We interpret these data in terms of the long-term north-south shortening and clockwise rotation of the Cascadia forearc. During the Eocene, the Cedar Creek fault was preferentially aligned with north-south shortening, but became less active when clockwise rotation brought the Doty fault into east-west alignment with the stress field. Sediment accumulation rates decreased from 196 m/Myr in the Paleogene (∼40–20 Ma) to 27 m/Myr in the Neogene (∼20–0 Ma) as Miocene uplift of the Coast Range limited accommodation space and shifted the depositional environment from marine to fluvial. Our results are consistent with the geologic evolution of the Portland and Tualatin basins to the south and reveal a depositional and structural history uniquely shaped by clockwise rotation of the Cascadia forearc.

Washington↗
USGS2025

Flow heterogeneity controls dissolution dynamics in topologically complex rocks

Rock dissolution is a common subsurface geochemical reaction affecting pore space properties, crucial for reservoir stimulation, carbon storage, and geothermal energy. Predictive models for dissolution remain limited due to incomplete understanding of the mechanisms involved. We examine the influence of flow, transport, and reaction regimes on mineral dissolution using 29 time-resolved data from 3D rocks. We find that initial pore structure significantly influences the dissolution pattern, with reaction rates up to two orders of magnitude lower than batch conditions, given solute and fluid-solid boundary constraints. Flow unevenness determines the location and rate of dissolution. We propose two models describing expected dissolution patterns and effective reaction rates based on dimensionless metrics for flow, transport, and reaction. Finally, we analyze feedback between evolving flow and pore structure to understand conditions that regulate/reinforce dissolution hotspots. Our findings underscore the major impact of flow arrangement on reaction-front propagation and provide a foundation for controlling dissolution hotspots.

Geophysical Research Letters↗