USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of the American Water Resources Association”

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,691 records · Page 41Linked to original sources
USGS2026

Long-term intermittent connection between the western Snake River Plain and Columbia basin: A two-phased incision history of Hells Canyon

For more than a century, researchers have debated the route of the Snake River across the northern Cordillera and U.S. Pacific Northwest, including the associated incision history of Hells Canyon. Here, we use detrital zircon U-Pb provenance analysis of Miocene strata upstream and downstream of Hells Canyon to constrain the evolution in fluvial pathways over time. Downstream of Hells Canyon, we show that the ca. 10−6 Ma Clarkston Heights gravel was dominantly sourced from the nearby Salmon and Clearwater rivers but also with substantial input from drainages that now flow into the western Snake River Plain. These same western Snake River Plain drainages were contributing to time-equivalent Lake Idaho strata, upstream of Hells Canyon, indicating a Hells Canyon fluvial connection between the western Snake River Plain and Columbia basin before 6 Ma. Results from Lake Idaho strata suggest a two-phased lacustrine history in the western Snake River Plain: Before ca. 6.0 Ma, lake strata were derived from local tributaries, suggesting that the western Snake River Plain was isolated from the modern upper Snake River system. Between ca. 4.3 Ma and 2.2 Ma, the source area expanded to include eastern tributaries in conjunction with migration of high-standing topography of the Yellowstone hotspot. Our results challenge the long-held hypothesis of “capture” of the Snake River through Hells Canyon. Instead, we infer that Hells Canyon was a long-established route for outflow of lakes occupying the western Snake River Plain, possibly intermittently, during a ca. 10−2 Ma phase of slow regional incision, followed by rapid incision in Hells Canyon starting ca. 2 Ma in conjunction with erosion of the topographic barrier impounding Lake Idaho.

Idaho, Oregon↗
USGS2026

USGS Geochron Database

Introduction Geochronology helps us understand Earth’s history by determining when important events, like volcanic eruptions, the rise of mountains, the formation of mineral resources, and changes in the landscape, happened. Geochronological data directly support geologic mapping and can inform decisions about geologic hazard mitigation, natural resource management, and infrastructure resilience. The U.S. Geological Survey (USGS) Geochron database provides access to more than 300,000 published, publicly available age measurements from more than 40,000 geological samples. This database is the result of a collaborative effort with State geological surveys and geoscientists from across the globe. The USGS Geochron database is the most comprehensive collection of geochronological data available for the United States. Users can view data through an interactive map explorer, download datasets, and integrate data into geospatial software or other analysis tools.

Fact Sheet↗
USGS2024

Preliminary ground and airborne-based geophysical mapping and modelling of an active hydrothermal system at Mammoth Lakes, California

Mammoth Lakes, California hosts a productive hydrothermal system within the seismically active south moat of Long Valley Caldera. Surficial evidence of the shallow hydrothermal system includes discrete zones of tree-kill dispersed between Shady Rest Park and the Casa Diablo Geothermal Power Plant (40 MW), as well as east of the power plant. The tree-kill areas are associated with elevated diffuse CO₂ emissions, heated ground, hydrothermal alteration, diffuse soil H₂S emissions, and gas vents. Previous mapping delineates prominent north and northwest trending structures within the south moat along the southwestern edge of the resurgent dome that may accommodate gas and fluid flow at the Shady Rest Park and Basalt Canyon Tree Kill Areas (SRTKA and BCTKA, respectively). Both tree-kill areas are also located along contacts between resurgent rhyolite, mafic lavas, and surficial deposits which may provide additional pathways for gas and fluid migration in the shallow subsurface. Characterizing structure and lithology using geophysical anomalies is critical to determining primary structural controls on the hydrothermal system and the extent of subsurface alteration at these sites. We conducted ground and airborne-based potential field geophysical surveys to map gravity and magnetic anomalies. These anomalies are then used to model subsurface geology, structure, and hydrothermal alteration. Here we present our preliminary geophysical mapping and modelling results at both tree-kill locations. Gravity and magnetic data suggest complex structural intersections are coincident with heated ground and gas emissions at the SRTKA and BCTKA. Hydrothermal systems are often observed or interpreted to exploit fault intersections which can serve as highly permeable pathways for hydrothermal fluid and gas discharge, enabling economic geothermal energy production. Geophysical mapping and modelling are an effective means of investigating such structural complexity at Mammoth Lakes due to the presence of unidentified and concealed structures.

California↗
USGS2026

Site response models based on geometric parameters for southern California sedimentary basins

Site response in sedimentary basins is influenced by complex three-dimensional (3D) features, including trapping of seismic waves, focusing of seismic energy and basin resonance. Current ground motion models (GMMs) incorporate basin effects using one-dimensional parameters like V S30 and shear wave velocity isosurface depths, which are limited in capturing lateral and 3D effects. To address these limitations, we develop seismic site response models based on novel parameters that represent multi-dimensional properties of the Los Angeles Basin (LAB) geometry and shear wave velocity. We define a basin shape for the LAB using depth to subsurface geologic interfaces associated with the oldest sedimentary deposits (depth to a particular shear wave velocity horizon, i.e., 1.5 km/s - z 1.5 ) and the depth to the crystalline basement ( z cb ) which are determined using geologic cross sections and community seismic velocity model profiles. We explore a suite of geometric descriptors computed for the LAB and southern California, from which three parameters with the greatest predictive potential are selected and evaluated using empirical ground motion residual analyses in combination with the Boore et al. GMM. The results demonstrate that the zonal heterogeneity index ( ), standard deviation of the absolute difference between z 1.5 and z cb ( ) and standard deviation of z cb ( ) each provide a reduction in site-to-site variability ( ϕ S2S ) of empirical GMMs. The reduction in ϕ S2S is period-dependent, with average decreases of 3%, 26% and 6% for , , and , respectively. Although these reductions are modest from an engineering application perspective, they are statistically significant, underscoring the inherent difficulty in fully characterising complex basin effects. Collectively, these findings indicate that the inclusion of basin-specific geometric parameters yields measurable, albeit incremental, improvements in site response prediction and establishes a framework for the progressive refinement of seismic hazard characterisation within sedimentary basins.

California↗
USGS2026

Software to support remote sensing of river discharge based on critical flow theory

Water resource management requires accurate observations of streamflow but standard field methods for measuring river discharge ( Q ) are costly and can be hazardous for equipment and personnel. Remote sensing has become a viable alternative, but many image-based techniques require field data for calibration and depth and velocity can seldom be mapped with a single sensor. A new approach based on critical flow theory, in contrast, allows both of these attributes to be inferred from readily available image data. This technique only pertains to sites with standing waves, called undular hydraulic jumps (UHJs), but a recent investigation demonstrated its potential to provide accurate discharge estimates. This paper introduces software designed to facilitate Inferring Q from UHJs Identified in River Images (InQUIRI). The package includes modules for retrieving data from image servers, making the measurements of wavelength and width required to calculate discharge, inferring a representative wavelength from a profile digitized along a wave train, combining multiple estimates to obtain an ensemble median discharge, and assessing accuracy via comparison to gage records from the U.S. Geological Survey. By making these steps easier to implement, InQUIRI enables users to apply the workflow to a variety of UHJ-containing images. Accumulating more case studies, some successful and others less so, would help constrain the range of applicability of the critical flow approach and foster development of refined guidelines for selecting and measuring waves. The software described herein could play an important role in promoting informed use of this new technique for non-contact streamflow measurement.

Arizona, Colorado, New Mexico, Utah↗
USGS2026

The mineral chemistry networks of tin and tungsten reflect metallogenic eras of the Mesozoic

Continental remobilization is a crucial driver of metallogenesis and the formation of ore deposits. Some of the world’s largest mineral deposits of the economically valuable elements tin (Sn), tungsten (W), gold (Au), copper (Cu), lead (Pb), and zinc (Zn) formed during the Mesozoic Era. Additionally, the chemistry and distribution of the elements Sn and W have been investigated in previous studies to understand planetary formation and differentiation processes. These two elements are largely co-located during certain South China Mesozoic metallogenic events but are not co-located during other time periods in the same regions. Here, we investigated the mineral chemistry network similarities and dissimilarities of Sn and W to understand their mineral formation and distribution during the Mesozoic Era and throughout Earth history. Mineral chemistry network community detection analysis and electronegativity associations among mineral constituent elements of Sn minerals and W minerals indicate that the elements have similar chemistry among their oxide minerals. However, Sn forms a much wider range of minerals that also contain S compared to W, which occurs in a limited number of S-containing minerals. The divergent constituent element interactions among S-containing Sn minerals and W minerals reflect the redox sensitivity and importance of oxygen (O) fugacity in Sn mineral formation. Conversely, extensive W mineral deposits are known to form at both high and low O fugacities. The similarities and differences between the mineral chemistry networks of Sn and W reflect the mineral distribution of the two elements in the Sn-W mineralization event from 160 to 139 Ma vs. the Sn–uranium (U) mineralization event from 125 to 98 million years ago (Ma). The mineral chemistry and distribution of Mesozoic Sn and W deposits illustrate the contrasting importance of redox and O fugacity on the mineral formation of different elements, and the dynamic crustal evolution that took place during this period of Earth history.

Geosciences↗
USGS2012

Development of vegetation based soil quality indices for mineralized terrane in arid and semi-arid regions

Soil quality indices (SQIs) are often management driven and attempt to describe key relationships between above- and below-ground parameters. In terrestrial systems, indices that were initially developed and modified for agroecosystems have been applied to non-agricultural systems in increasing number. We develop an SQI in arid and semi-arid ecosystems of the Western US impacted by different types of geologic mineralization using the relationship between vegetation community parameters and soil abiotic and biotic properties. We analyze these relations in soils associated with three different mineralization types: podiform chromite, Cu/Mo porphyry, and acid-sulfate gold vein systems at four different sites in California and Nevada. Soil samples were collected from undisturbed soils in both mineralized and nearby unmineralized substrates as well as from waste rock and tailings. Aboveground net primary productivity (ANPP), canopy cover and shrub density were measured for the vegetative communities. Minimum data sets were developed based on correlations between the soil and vegetation parameters, refined using principal components analysis, scored using non-linear functions, and combined into an overall SQI. The indices are comprised of one or two microbial parameters and three to six abiotic parameters, the latter consisting of nutrients and metals. Given the preliminary development of this approach, the parameters and combinations to arrive at an SQI for a given site cannot at this time be correlated or compared with that of another site. This SQI approach provides a means of quantifying disturbed ecosystem recovery resulting from mining, and could be applied to other disturbances in a way that readily distills the information for potential use by land managers. However, severely disturbed areas with little to no aboveground biomass, such as unreclaimed tailings, have likely crossed an ecological threshold that precludes the use of this type of monitoring tool.

California↗
USGS2025

Temporal and spatial comparison of coal mine ventilation methane emissions and mitigation quantified using PRISMA satellite data and on-site measurements

Emission monitoring at the facility level (bottom-up, BU) is key for accurate reporting of coal mine methane (CMM) emissions. Recent advances in aerial and satellite observations (top-down, TD) indicate that these methods have the potential to support CMM emissions monitoring and reporting of, as well as track the effectiveness of, mitigation actions. That said, studies have shown discrepancies between BU and TD estimations. Performing TD monitoring with concurrent BU measurements at the same mines may help address the observed discrepancies and improve quantification methods to narrow the gap between BU and TD data. This paper presents a comparison of concurrent BU-TD methane emission fluxes monitored from two ventilation shafts at a southwestern Virginia, USA, longwall mine to complement the existing body of studies on satellite-based monitoring of coal mines by incorporating continuous ground-based monitoring with concurrent TD monitoring of methane emission fluxes to address the gap and provide valuable insights into temporal emission patterns. The shafts were monitored on multiple dates between 2020 and 2023. BU monitoring was performed at the exhaust fans, while TD quantifications used PRISMA hyperspectral satellite data and two different wind reanalysis datasets (i.e., GEOS-FP and ERA5). This mine also offered a unique opportunity for BU-TD comparisons before and after ventilation air methane (VAM) oxidizer operation, which was installed at one of the shafts in 2022. The results showed that TD-estimated mean fluxes were generally lower than BU data, which were attributed to quantification difficulties associated with the low albedo caused by heavy vegetation and the terrain of the area. However, despite the discrepancies in mean emissions, both the interquartile range and the data range of the distributions generally overlapped, and the estimates correctly showed the emission trends.

Virginia↗
USGS2024

Are the horizontal-to-vertical spectral ratios of earthquakes and microtremors the same?

We consider the similarities and differences between earthquake and microtremor horizontal‐to‐vertical spectral ratios (eHVSR and mHVSR, respectively) using a dataset of 161 sites in southern California. Quantitative comparisons are made in terms of the eHVSR and mHVSR lognormal median curves, as well as the frequencies and amplitudes associated with the fundamental‐ and higher‐mode resonances where present. The results show only 58% of the eHVSR–mHVSR pairs agree in terms of their median curve and only 25% of the eHVSR–mHVSR pairs agree in terms of shared resonances, which increases to 68% if flat HVSRs are considered equivalent. Furthermore, while the shared resonances match very well in terms of frequency (root mean square error, RMSE, <0.11 Hz), the amplitudes of those resonances do not agree (RMSE >1.6). These findings demonstrate that while eHVSR and mHVSR agree at some sites, they are not equivalent at all sites. To investigate if the agreement between eHVSR and mHVSR could be related to features of the microtremor data, earthquake recordings, and/or the site conditions, three machine learning (ML) models at varying levels of interpretability are presented. The ML models—which include multivariate logistic regression, gradient‐boosted trees, and support vector machines—show only partial success at using site‐specific data to predict whether eHVSR and mHVSR will likely agree in terms of their median curve (accuracy of 78%) and number of resonances (accuracy of 84%). Therefore, we conclude that while eHVSR and mHVSR can be quite similar in terms of resonant frequencies at some sites, they are not identical at all sites. Furthermore, preliminary evidence shows that the agreement of eHVSR and mHVSR can be predicted a priori given features of the microtremor measurements, earthquake recordings, and site conditions, although a larger dataset will be necessary for developing a robust predictive model.

California↗
USGS2026

Computation of regional groundwater budgets for the Virginia Coastal Plain aquifer system

Computation of detailed groundwater flow budgets for subdivisions of Virginia’s Coastal Plain aquifer system has enabled quantification and more thorough understanding of groundwater flow within this important water resource. A zone budget analysis conducted on previously published groundwater models of the Virginia Coastal Plain and Virginia Eastern Shore shows that groundwater conditions vary substantially throughout the Coastal Plain aquifer system due to local variations in hydrogeology and historical and ongoing variations in groundwater use and management. Decades of substantial groundwater withdrawal from the Coastal Plain aquifer system have fundamentally altered groundwater flow from pre-development conditions. Rates of sustainable withdrawal are limited because the downward groundwater flow rate into confined aquifers supplying groundwater is a relatively small portion of the total groundwater water budget for the aquifer system. Analyses of groundwater budgets from the Virginia Coastal Plain model show that groundwater flow is generally outward from the surficial aquifer to rivers and coastal water bodies and downward through a series of underlying aquifers and confining units to the Potomac aquifer, which is the deepest aquifer and the source of most groundwater withdrawals. Downward flow into the Potomac aquifer currently is estimated to be only 7 percent of total net precipitation-derived net recharge at the land surface but makes up about 66 percent of inflow to the aquifer in Virginia, with much of the remaining inflow occurring laterally from areas outside of defined groundwater budget regions in Virginia. For several decades prior to 2010, high rates of withdrawal from the Potomac aquifer resulted in substantial decline in groundwater storage in the aquifer and in most overlying aquifers and confining units. From 2010 to 2025, rates of withdrawal substantially lower than the historical maximum have resulted in small net increases in groundwater storage in the confined aquifer system for most regions of the Virginia Coastal Plain. Nevertheless, for the same period, groundwater storage for the entire model domain continues to incrementally decline, indicating that storage recovery in Virginia is offset by a continued decrease in storage in areas beneath the Chesapeake Bay or in adjacent areas of Maryland and North Carolina. Withdrawals from the Potomac aquifer have induced substantial downward flow which is a large part of groundwater budgets for confined aquifers such as the Potomac. Downward groundwater flow continues under current conditions, but because vertical flow rates are a function of the difference between water pressure in the upper surficial systems and lower confined units, those rates are lower than those in earlier decades as the confined water levels partially recover from larger groundwater withdrawals in the past. Geographically, groundwater flow is generally inward from perimeter regions of the Virginia Coastal Plain toward central regions with the largest withdrawal rates. Estimated groundwater inflow from coastal regions could be contributing to saltwater intrusion, though that was not measured directly in this study. Analyses of groundwater budgets from the Virginia Eastern Shore peninsula, a geographic region of the Virginia Coastal Plain, show that groundwater flow for that isolated aquifer system is generally outward from the surficial aquifer to coastal water bodies and downward into the confined Yorktown-Eastover aquifer system, which is the source of most withdrawals. Downward groundwater flow into the confined Yorktown-Eastover aquifer system is estimated to be less than 2 percent of total recharge and less than 9 percent of net recharge at the water table but makes up over 93 percent of all inflow to the confined aquifer system. Decades of substantial but relatively consistent groundwater withdrawals have induced greater downward flow rates into the confined aquifer system but also have resulted in loss of groundwater from storage. Currently, estimated storage loss accounts for slightly under 7 percent of withdrawals from the confined aquifer system. The current withdrawal rate from the confined Yorktown-Eastover system is near the highest reported rate for the Eastern Shore, which means that the storage depletion is expected to continue, even though groundwater levels appear to be relatively stable. Estimated groundwater flow rates upward from the confining unit underlying the Yorktown-Eastover system and small rates of inflow from coastal water bodies underscore ongoing concerns about up-coning and lateral intrusion of salty groundwater.

EarthArXiv↗
USGS2025

Combining multisite tsunami and deformation modeling to constrain slip distributions for the 1700 C.E. Cascadia earthquake

A major earthquake ruptured the Cascadia subduction zone (CSZ) on 26 January 1700. Key paleoseismic evidence associated with this event include tsunami deposits, stratigraphic evidence of coastal coseismic subsidence, written Japanese records of a tsunami unaccompanied by earthquake shaking, and margin‐wide turbidites found offshore and in lacustrine environments. Despite this wealth of independent clues, important details about this event remain unresolved. Dating uncertainties do not conclusively establish whether the proxies are from one earthquake or a sequence of them, and we have limited knowledge of the likely slip distributions of the event or events. Here, we use a catalog of 37,500 candidate synthetic ruptures between M w 7.8 and 9.2 and simulate their resulting coseismic deformation and tsunami inundation. Each model is then compared against estimated Japan tsunami arrivals, regional coastal subsidence records, and local paleotsunami deposits mapped at six different coastal marshes and one coastal lake along the CSZ. We find that seven full‐margin ruptures with a median magnitude of M w 9.1 satisfy all three constraints. We favor one M w 9.11 model that best matches all site paleoseismic observations and suggests that the Cascadia megathrust slipped up to ∼30 m and must have shallow geodetic coupling. We also find that some sequences composed of three or four ruptures can still satisfy the observations, yet no sequences of two ruptures can. Sequences are differentiated into three groups based on whether they contain a mainshock rupture located in the south (>44° N) or further north. All sequences contain unruptured portions of the megathrust and most contain mainshocks with peak slip above 40 m. The fit of the geologic evidence from sequences is poor in comparison to single‐event models. Therefore, sequences are generally less favored compared to full‐margin events.

British Columbia, California, Oregon, Washington↗
USGS2024

Precariously balanced rocks in northern New York and Vermont, U.S.A.: Ground-motion constraints and implications for fault sources

Precariously balanced rocks (PBRs) and other fragile geologic features have the potential to constrain the maximum intensity of earthquake ground shaking over millennia. Such constraints may be particularly useful in the eastern United States (U.S.), where few earthquake‐source faults are reliably identified, and moderate earthquakes can be felt at great distances due to low seismic attenuation. We describe five PBRs in northern New York and Vermont—a region of elevated seismic hazard associated with historical seismicity. These boulders appear to be among the most fragile PBRs in the region, based on reports from hobbyists. The PBRs are glacial erratics, best evidenced by glacial striations on bedrock pedestals. The pedestals themselves are locally high knobs, often situated on regionally high topography; this setting limits soil development and indicates that any outwash deposits were likely ephemeral. As a result, PBR ages can be reliably established by the retreat of the last continental ice sheet, ∼15–13 ka. To quantify the fragility of the PBRs, we surveyed them with ground‐based light detection and ranging and calculated geometric parameters from the point clouds, field observations, and seismic responses. Preliminary validation of the 2023 time‐independent U.S. National Seismic Hazard Model (NSHM) shows that the existence of PBRs is generally consistent with the median site‐specific hazard curves. Only the Blue Ridge Road site suggests a modest reduction in hazard. To visualize the ensemble of data, we mapped the minimum permissible distance to potential source faults around each PBR site as a function of source magnitude by using the ground‐motion models from the 2023 NSHM. Viewed in this manner, our data are consistent with potential M ∼6.5 earthquake‐source faults in many parts of the Lake Champlain Valley and northern Adirondack Mountains. Our work illustrates a potential pathway for better constraining earthquake‐source faults in regions of cryptic faults.

Bulletin of the Seismological Society of America↗
USGS2021

New York and Landsat

From the iconic skyline of New York City to the forested landscapes of the Adirondack Mountains and the countryside of the Allegheny Plateau, the State of New York is overflowing with diversity and life. Bordered by the Atlantic Ocean on the east and two of the Great Lakes to the north and west, New York has more than 7,600 lakes, ponds, and reservoirs and more than 70,000 miles of rivers and streams. New York’s stewardship of its freshwater resources is fundamental to the health and well-being of all who work at, reside in, and visit the State’s landmarks and places. Harmful algal blooms in the State’s waterbodies are a growing concern and threaten the health of the region and its inhabitants. Images and data from Landsat satellites continue to provide critical information to scientists, public health officials, and resource managers who are studying the effects and risks of the problem.

New York↗
USGS2026

A practical decision tool for marine bird mortality assessments

Given the rise in anthropogenic, environmental, and disease events contributing to marine bird mortality, there is a critical need to improve the rigor of mortality assessments. Deficits in data collection and mortality estimation can hinder a manager’s ability to document the scale of events and assess population level impacts. Therefore, to inform decisions required during activities, such as conservation status assessments or harvest management, organizations may choose to incorporate mortality assessments into response plans. Resources, capacity, and assets to assess mortality vary across jurisdictions (federal, state, Indigenous, local, etc.), and clear guidance to support mortality assessments is often unavailable or not clearly addressed. Here, we present a decision support tool to help managers identify and evaluate survey options to assess bird mortality in a diverse array of scenarios. The objective of the decision tool is to improve data collection and availability, which will increase the ability to estimate mortality robustly, given situation-specific attributes and constraints. This decision tool is designed to guide the response when a mortality event is initially encountered and offers suggestions for assessment and reporting procedures in the absence of other guidance or to complement existing protocols. The decision tool is also meant to inform decision making for response determination and resource allocation. The tool facilitates examination of options for further assessment and monitoring, which users determine by examining questions pertaining to species prioritization, determination of mortality minimum spatial extent, and the potential magnitude of impacts on affected species. Finally, identification of appropriate survey methods that address imperfect detection when a complete census is not possible are determined by exploring location, spatial and temporal extent, and the type of species affected. Ultimately, this decision tool aims to facilitate and improve the standardization of mortality assessments, equipping managers with a practical resource to navigate the decision-making process for marine bird mortality estimation.

Ornithological Applications↗
USGS2026

Preliminary bedrock geologic map of the Port Henry quadrangle, Essex County, New York, and Addison County, Vermont

Introduction The bedrock geology of the 7.5-minute Port Henry quadrangle 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. 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 Ga (giga-annum). In the Port Henry quadrangle, the AMCG metaigneous rocks (Yhg, Ygb, Yanw) intruded older, mostly metasedimentary rocks of the Grenville Complex during the middle to late Shawinigan orogeny (~1,160–1,150 Ma [mega-annum]). All rocks were subsequently metamorphosed to upper amphibolite to granulite facies conditions during the 1,080–1,050 Ma Ottawan orogeny. New mapping reveals four periods of deformation: (1) D1 produced rarely preserved isoclinal folds in the paragneiss and marble and predates AMCG magmatism. (2) Subsequent D2 deformation produced the dominant gneissic fabric preserved in the rock, recumbent folding, and deformed all the Proterozoic units in the map area. Syn- to late-D2 felsic magmatism resulted in the regionally extensive Lyon Mountain Granite Gneiss, which hosts numerous magnetite ore bodies. (3) Mylonitic extensional shear zones and core complex formation marked the beginning of D3 deformation. Protracted D3 deformation resulted in F3 upright folding, dome and basin formation, pegmatite intrusion, 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. (4) D4 created northeast- and northwest-trending local high-grade ductile shear zones and boudinage, northwest-trending regional kilometer (km)-wide ductile shear zones, and crosscutting granitic pegmatite dikes. The development of the late-stage regional shear zones (D4) was likely due to the continuation of extensional doming and uplift from upper amphibolite facies conditions at the end of the Ottawan orogeny. The majority of iron-ore deposits in the Port Henry and adjacent Witherbee quadrangles are in the hanging wall of these extensional shear zones. In the Port Henry quadrangle, the km-wide Cheney Mountain shear zone is the result of D4 deformation. Kilometer-scale lineaments readily observed in lidar data are Ediacaran mafic dikes and Phanerozoic brittle faults. The Paleozoic rocks are part of the Early Cambrian to Late Ordovician carbonate bank on the ancient margin of Laurentia. The approximately 1-km-thick Cambrian to Ordovician stratigraphy records a transition from synrift clastics to passive-margin peritidal carbonate buildups to gradually deeper-water subtidal- to shelf-carbonates during foreland basin development associated with the Taconic orogeny. The Paleozoic rocks are weakly folded and block faulted. Large areas of the Champlain Valley are covered by undifferentiated glacial deposits, some of which contain mapped landslides. The map also shows waste rock piles and tailings from 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 Adirondacks, provide a context for historical iron mines in the eastern Adirondacks, and update the stratigraphy of the Champlain Valley in New York and Vermont. This Open-File Report includes a bedrock geologic map; a description of map units; a correlation of map units; and a geographic information system database that includes bedrock geologic units, faults, outcrops, and structural geologic information.

New York, Vermont↗
USGS2026

Distributed faulting of the northern West Napa Fault Zone in Napa Valley, California

Mapped surface ruptures from the 24 August 2014 M w 6.0 South Napa earthquake in the Napa Valley, California, show a 2‐km‐wide zone of distributed faulting in the southern and central West Napa fault zone (WNFZ). In the northern WNFZ at Hendry Winery (HW), however, the mapped 2014 surface ruptures encompass an ∼100‐m‐wide zone, implying significant narrowing of the near‐surface fault zone to the north. We present a tomographic shear‐wave velocity ( ⁠⁠ V S ) model and guided‐wave data that indicate the northern WNFZ is at least 400‐m wide, with multiple near‐surface fault traces. Our V S model shows that the 2014 surface ruptures are underlain by discrete low‐velocity zones (LVZs), and coincident guided‐wave data show that the LVZs carry fault‐zone guided waves. If nearby (<500 m) mapped faults to the east of HW are part of the WNFZ, the entire WNFZ is more than 1 km wide in the northern Napa Valley. WNFZ guided waves travel up to 38% slower than S body waves, and low‐strain guided‐wave shaking is up to five times stronger than the associated body‐wave shaking. Our data suggest that guided waves, traveling along distributed faults, may result in an increased shaking hazard over a 1‐km‐wide area of the northern Napa Valley during future significant earthquakes. In places, the 2014 surface ruptures were difficult to find one year after the earthquake, and paleoseismic trenching showed only weak evidence for faulting, which may not have been identified in trenches if the locations of the 2014 surface ruptures had not been previously mapped ( Prentice et al. , 2015 ). Guided‐wave and V S tomography data, however, show strong evidence for faulting beneath the 2014 surface ruptures and at locations to the east. Although paleoseismic trenching is the gold standard for identifying near‐surface faulting, methods such as peak ground velocities of guided waves may better identify immature near‐surface fault traces.

California↗
USGS2025

Structural analysis of brittle-plastic shear zones in the Sangre de Cristo Range, southern Colorado USA: Superposition of Rio Grande rift extension on Laramide contraction

The Sangre de Cristo Range in southern Colorado exposes some of the deepest Cenozoic structural levels in the Rocky Mountain region, including mylonitic shear zones associated with both the Laramide orogeny and Rio Grande rift. We investigated the relation between Laramide contraction and Rio Grande rift extension with detailed geologic mapping, kinematic analysis, and geochronometry in a 50 km 2 area centered on the Independence Mine shear zone (IMSZ). The 15−100-m-thick IMSZ is one of several shallowly to moderately (∼45° ± 20°) W-SW−dipping brittle-plastic shear zones along the western flank of the range. These shear zones display microstructural evidence of initiation as top-NE contractional mylonite zones, consistent with regional Laramide kinematics, which have been pervasively overprinted by shear fabrics indicating top-SW extensional reactivation. Both top-NE and top-SW shear fabrics involve cataclasis and quartz dislocation creep, although top-SW shear is more commonly localized along phyllosilicate-lined shear bands. Shear zones are hosted predominately within Proterozoic gneiss, and contain abundant chlorite and white mica derived from alteration of hornblende and feldspar, which indicates that weakening driven by fluid reactions played an important role in localizing strain. Extensional overprinting appears to be most pervasive along more steeply dipping portions of shear zones and where secondary phyllosilicates form an interconnected weak phase, which suggests that reactivation was primarily controlled by geometry and rheological contrasts inherited from contraction. One top-SW shear zone adjacent to the IMSZ cuts a late Oligocene gabbro stock, and monazite grains synkinematic with top-SW shear in the IMSZ yielded late Oligocene to Early Miocene U-Th-Pb dates that correspond with initiation of the Rio Grande rift. Reactivation of weak reverse faults may represent an important structural control during initial extension in the middle crust, prior to slip along the high-angle Sangre de Cristo normal fault system.

Colorado↗
USGS2025

Critical mineral inventory of select IOA-IOCG deposits, southwestern USA

Critical minerals are necessary for modern technology and strategic purposes. Their increasing importance requires finding new and nontraditional resources. Samples of ore, altered, and unaltered host rock were collected from 26 iron mines and prospects in California, Nevada, and Utah to assess the potential of these deposits to host economic quantities of different critical minerals. Geochemical analyses were conducted by 61 element ICP-OES-MS sodium peroxide fusion and major elements determined by WDXRF. These deposits concentrated many critical minerals beyond what is found in average upper crustal abundances, such as Sb, As, Bi, Co, Ga, Mg, Mn, Ni, Nb, Pd, REE, Sc, Te, Sn, Ti, W, V, and Zn. However, most of these are not concentrated enough in the ore to be considered as economic resources. Those critical minerals that are enriched enough in some of these deposits to possibly be considered as by-product commodities are Ni, REE, V, and potentially Co and Ga. These enrichments were not uniform, with REE more likely to be enriched in IOA deposits, whereas Co, Ga, Ni, and V could be found enriched in either IOA or IOCG deposits.

California, Nevada, Utah↗