SEARCH · Geology Reports
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.
A U.S. Geological Survey program to monitor and assess channel reconfiguration activity is described. A data base available on the world wide web will enable land-management agencies and other interested parties to evaluate the long-term success of various channel reconfiguration projects. A demonstration project on the Lake Fork of the Gunnison River, Colorado, illustrates the program objectives and approach.
Many subsurface processes involve transitions in granular material states, from arrested to creeping to flowing. Experiments and frameworks for idealized systems reveal that granular fabrics develop during shearing, co-evolve with applied stress, and govern such transitions. We use microtomography to test whether fabrics at two San Andreas fault sites reflect slip history and whether idealized frameworks extend to nature. Near-surface sediments within the fault zone transition between deformation patterns over the seismic cycle, including bulk/localized grain re-arrangements, individual grain fracturing, and localized zones of fracturing. Aseismic and co-seismic shearing produce distinct preferred grain orientations. Co-seismic fabrics can be preserved after centuries of aseismic strain, aseismic fabrics may be overprinted, and grain size and coordination number influence the fabrics. Idealized frameworks, namely anisotropic critical state theory, frictional jamming, and material memory, can explain our observations, and fault zone sediments likely undergo cycles of memory creation and erasure that influence rigidity spatiotemporally.
We present a method to update the geospatial liquefaction model used by the U.S. Geological Survey’s near‐real‐time ground failure product with subsurface geotechnical data. The geospatial model estimates liquefaction probability from peak ground velocity (via ShakeMap) and geospatial susceptibility proxies. In many regions, additional information relevant to constraining liquefaction likelihood is also available, including surface geology maps and subsurface geotechnical measurements. There is currently no mechanism to use these data in the ground failure product liquefaction model, even though these data could provide more precise constraints on spatial variations in the lithologic character of the soil (surface geology) and direct measurements of the subsurface mechanical properties that affect liquefaction occurrence and severity (geotechnical measurements). In this study, we develop a method to integrate these data with the geospatial model and assess how these data can improve regional‐scale predictions. We develop a Bayesian updating framework and apply it to the 1989 magnitude 6.9 Loma Prieta, California, earthquake, for which mapped observations are available to evaluate performance. We constrain the Bayesian framework with 373 Northern California cone penetration tests and liquefaction susceptibility classes based on the mapped surface geology. This Bayesian model incorporates geotechnical information into the geospatial model and more accurately predicts liquefaction occurrences than the geospatial model, while sacrificing less accuracy in terms of predicting the absence of liquefaction than the geotechnical model. In future applications, this approach could be adapted to update other geospatial models using locally available subsurface data.
The 2011 Mw9.1 Tohoku-oki earthquake is exceptional among great earthquakes for having peak slip of ~50-70 m on the shallowest portion of the plate boundary megathrust. International Ocean Discovery Program Expedition 405 drilled multiple holes through the megathrust in the large slip region and at a Pacific Plate input site. The megathrust preferentially develops at the top or base of the pelagic clay in the input section where pronounced contrasts in physical properties are present. This results in a narrow, weak fault located at a major mechanical contact between frontal prism mud and subducted clay. Localization imposed by the input section enhances the tendency for shallow seismic slip, showing the Japan Trench may be more susceptible to ruptures with large shallow slip than margins without weak clays.
The electrical grid relies on minerals with unique properties to transform, transmit, and convert electricity. Transmission needs are growing—to support the expansion of data centers and link diverse energy resources into the grid. This expansion will require greater quantities of minerals. Although some minerals are produced in the United States, the United States relies on imports for others. This infographic shows the key minerals that support energy transmission, and how reliant the United States is on imports for these minerals.
Foreward This document contains design specifications of a navigational standard for the Moon, including a Lunar Transverse Mercator system, a Lunar Polar Stereographic system, a Lunar Grid Reference System, and a unique coordinate structure, Artemis Condensed Coordinates, for Artemis mission navigation and lunar surface science. The National Aeronautics and Space Administration (NASA) Artemis campaign seeks to place humans on the Moon for the first time since the Apollo missions. Early Artemis missions are heavily focused on the lunar south pole, which promises to return valuable data on the Moon’s geologic record, amongst other mission objectives. Coordinate systems in use today for the lunar south pole provides crew members on the surface neither an efficient nor intuitive means to communicate their position and orientation. A novel grid coordinate system, the Lunar Grid Reference System, is proposed to address these concerns for use in real-time extravehicular activity operations on the lunar surface. The many stakeholders involved in the Artemis missions will need a common system to communicate position and orientation while astronauts are operating on the lunar surface. To that end, Artemis crew members will need that system to be efficient and intuitive to promote efficient extravehicular activity timelines and reduce confusion. In the context of this document, these characteristics are addressed on the design of lunar coordinate systems: Efficient.—The number of characters required to communicate a location within a desired precision level in both local and global contexts, and how many steps are required for a recipient or sender to interpret a location. Intuitive.—How well the system aligns with human perceptual abilities, and whether the system yields distances that have the same relationship to actual lunar surface distance in all directions from the point where a person is located. Technological systems are currently being investigated to supplement the crew members’ ability to locate and orient themselves and other assets on the lunar surface; however, it is unlikely that those systems will be fully operational for the first few landed missions. Even with future positional aids, crew members will still need an efficient and intuitive means to communicate position and orientation. In addition, if technological systems fail, the crews will require land navigation skills and have maps available, thus providing further motivation for a crew-centric coordinate system. The contents of this U.S. Geological Survey (USGS) document detail a comprehensive framework for standardizing lunar crewed surface navigation within NASA and outlines the protocols, methods, and designs necessary for achieving consistency and interoperability across relevant space mission teams and lunar surface navigators. Key components of this document include designs of map projections, projected coordinate reference systems (Lunar Transverse Mercator and Lunar Polar Stereographic systems), and a grid system (Lunar Grid Reference System and Artemis Condensed Coordinates) for the Moon. The work proposed in this document seeks to accomplish something similar to the National Geospatial-Intelligence Agency (NGA) document SIG 0012 (NGA, 2014a), but for using grid systems for the Moon. This report incorporates initial feedback and input from NASA’s Artemis Geospatial Data Team, NASA’s Flight Operations Directorate, National Geodetic Survey, USGS Astrogeology Science Center, and NGA and is intended to serve as a resource for all involved with the Artemis missions, as well as for engineers designing and operating lunar infrastructure.
We investigate the microtextural, microchemical, and isotopic effects of late-stage ductile deformation in quartzite mylonites and kyanite–muscovite–quartz veins from the Raft River shear zone (Utah). Quartz microstructures record pervasive disequilibrium, expressed by unannealed features including undulatory extinction, deformation lamellae, and poorly defined fabrics, typical of waning deformation in shear zones. Microchemical disequilibrium is best preserved in kyanite quartzites, where CL-zoned kyanite records repeated fractures, overgrowth, and mineral precipitation, and in muscovite from muscovite-poor quartzite mylonites that shows minor-element zoning consistent with syn-deformational overgrowth on detrital cores. In contrast, muscovite from muscovite-rich kyanite quartzites exhibits minimal chemical zoning. These microchemical variations correlate with 40 Ar/ 39 Ar age systematics. Chemically zoned muscovite preserves variable single-step ages, including ∼150 Ma ages, reflecting retention of detrital cores. In contrast, syndeformational muscovite consistently yields Miocene ages, indicating recrystallization and new growth below argon closure temperatures that reset inherited isotopic signatures. Similar trends are observed in quartzite mylonites, where increasing quartz recrystallization and stronger crystallographic preferred orientations occur toward deeper structural levels. Together, these observations indicate increasing retrograde deformation and recrystallization with depth in the Raft River shear zone and demonstrate that strain-driven recrystallization exerts a first-order control on muscovite ages. We suggest that apparent thermochronologic gradients in retrograde shear zones may reflect recrystallization gradients rather than temperature gradients. Where cooling limits the thermal driving force for recrystallization, isotopic relics are preserved, and local deformation and fluid availability control re-equilibration. Consequently, isotopic disequilibrium—particularly in the 40 Ar/ 39 Ar system—may be the rule rather than the exception of retrograde tectonic environments.
The Bureau of Economic Geology at the University of Texas at Austin (UT-BEG) is leading the Gulf Coast Carbon Ore, Rare Earth, and Critical Minerals (CORE-CM) Initiative to assess the potential to produce critical minerals (CMs), including rare earth elements (REEs) from coal, coal ash, and produced water related to oil and gas production, and related materials (alumina processing waste [red mud], heavy mineral sands, graphite, and zeolite) within the Gulf Coast Basin. This project represents the first phase in a long-term program and provides reconnaissance data that will be foundational for future work by assessing resources and suggesting plans to be conducted in future work and expanding stakeholder engagement. The project includes several tasks designed to identify, characterize, and assess several necessary aspects for development of CMs and REEs in the Gulf Coast Basin.
The Southern California Borderland (SCB) is a topographically complex region on the active continental margin that hosts varied hardground habitats, including ferromanganese (FeMn) crusts and phosphorites, marine minerals being considered for resource extraction. The SCB is influenced by seasonal upwelling and terrestrial inputs, and has a well-defined oxygen-minimum zone (OMZ). We analyzed megafaunal community composition, density, and diversity across 41 ROV video transects at 10 SCB hard substrate sites spanning depth (378–2765 m), temperature (1.79–7.96 °C), and oxygen (3.91–105.67 μM) gradients, with varying rock types (FeMn crusts, phosphorites, other rocks). We counted 32,426 individuals representing 146 unique taxa over a total area of 21,935 m 2 . Echinodermata contributed 48 % of the total abundance; Cnidaria 24 %; Porifera 13 %; Annelida 6 %; Arthropoda 6 %; Chordata 2 %; Mollusca, Foraminifera, and Hemichordata <1 %. Megafauna communities showed clear heterogeneity, with density, diversity, and community composition varying among transects within sites and across sites. Rock type recorded no influence on megafaunal density, but habitats that included FeMn crusts had higher diversity and distinct taxonomic composition. Location, depth, and oxygen were the most important variables explaining variance among megafaunal communities, with distinct communities formed at deeper depths and within the OMZ. This study provides baseline information about the megafauna inhabiting SCB hardgrounds, particularly those associated with minerals considered for their resource potential. It extends existing documentation of such communities to greater depths, reveals increased representation of sponges and corals below 1000 m, and highlights the variability within and among deep-sea megafaunal communities, offering context for informed policy decisions.
We are leading an effort to develop international standards and guidelines for curating, disseminating, and validating simulated ground-motion data. This effort is organized as a working group within the Consortium of Organizations for Strong Motion Observation Systems (COSMOS). In 2022 we held online workshops on 7 and 8 June and 20 October. The first workshop focused on curating and disseminating simulated ground-motion data. The second workshop focused on validating simulated ground-motions for engineering applications. About 100 people participated in each of the workshops with strong representation from North America, Europe, and western Asia. In the coming year, we intend to form a technical committee to draft the international guidelines and standards while continuing to engage stakeholders from across the globe through online and in-person meetings. The key points from the workshops include: (1) Numerous groups are generating simulated earthquake ground motions and making them openly available; however there is very little coordination among groups to provide consistent interfaces for searching and retrieving data; (2) Standardizing interfaces for metadata and data access should consider agile approaches that can adapt to changing capabilities and user needs while building upon existing efforts; (3) A distributed architecture is preferred by most participants to allow institutions to host and manage their own data while broadcasting their holdings to a combined catalog; (4) Validation of ground-motion simulations applies to the entire workflow for simulating earthquake ground-motions, including the rupture model, seismic velocity model, and seismic wave propagation software; (5) Metrics to evaluate the validation are application dependent; more research is needed to tie structural response characteristics to ground-motion characteristics; and (6) Validation results should provide a clear, transparent, and quantitative assessment of the simulated ground motions.
Reductions in snow have left many streams in northern New Mexico dry or with very low flows during the summer. Base flow, defined as the contribution of groundwater to streamflow, can sustain streamflow during the summer and during drought conditions. Beaver dam analogs can be used to “slow the flow,” or increase the infiltration of rain-based runoff to replenish shallow groundwater reservoirs and increase base flow and potentially increase summertime streamflow.
Climate change is increasing sulfate export and changing wetland extent in mountain regions. These changes may increase microbially mediated production of the neurotoxic substance methylmercury due to enhanced sulfate metabolism in mountain environments. Here, we assess methylmercury concentrations and formation rates across high-elevation wetlands in the Colorado Rocky Mountains. We also investigate sulfate controls on methylmercury production within subalpine peatlands by amending soils with sulfate to mimic increased stream export of sulfate from the alpine zone and measuring methylmercury formation rates for different sulfate treatments. We found that subalpine peatlands have statistically significant higher methylmercury concentrations and formation rates compared to alpine, mineral-soil wetlands. Methylmercury production in subalpine peatlands also increased significantly ( p < 0.05) following sulfate additions; the highest rates occurred in sediments with intermediate extractable sulfate concentrations (∼0.60–1.4 mg sulfate g −1 dry soil). Our study is the first to identify soil sulfate-related thresholds for methylmercury production and sulfate-limitation of methylmercury production in subalpine peatlands. These findings highlight important linkages between climate-driven mineral weathering and mercury cycling in mountain regions globally.
On July 9, 1926, monitoring well Oe-151 at Woodgate, New York (USGS ID 433112075091501) recorded its first groundwater data. Since then, the well has provided water data nearly continuously and has now reached a 100-year milestone for data collection. The well is part of the U.S. Geological Survey (USGS) Climate Response Network (CRN), which is a national network of wells selected to monitor natural groundwater conditions. Well Oe-151 is the first well in the network to reach a 100-year Centennial milestone.
Near-surface seismic velocity structure may significantly impact the intensity, duration, and frequency content of ground shaking during an earthquake. In this study, we compile 649 shear wave velocity (Vs) profiles throughout the U.S. Pacific Northwest and southern British Columbia (PNW) and use these measured profiles to develop a representative soil velocity model for four major Holocene soil provinces: Puget Lowlands, Willamette Valley, fill and alluvium, and `other' soils. The resulting soil velocity model shows good agreement to measured data for a wide range of site conditions, with variability between different geologic domains reflecting fundamental differences in depositional environments. We then show that using this regional soil velocity model in simulations of the 2001 M6.8 Nisqually, Washington earthquake improves the fit to observed high-frequency (≥ 0.5 Hz) ground motions in the Puget Sound region compared to simulations that do not incorporate shallow (≤ 200 m) seismic velocity structure. Overall, this work shows that incorporating localized soil velocity profiles into seismic velocity models is important for accurately estimating high-frequency ground motion and regional seismic hazard in earthquake simulations. Future earthquake simulations and hazard studies in the PNW could incorporate these soil velocity profiles to capture the region's distinct site response characteristics.
We evaluate, select, and describe the ground-motion models (GMMs) used in the 2025 update of the U.S. National Seismic Hazard Model (NSHM) for Puerto Rico and the U.S. Virgin Islands (PRVI). We identify the most appropriate models that align with GMM selection criteria for use in the PRVI region to improve the accuracy of seismic hazard assessments. The update incorporates globally applicable GMMs suited for the active crustal and subduction earthquakes in the region. We include region-specific adjustments to these GMMs derived from local site response analyses derived from ground motion records. The unadjusted and regionally-corrected GMMs are combined to create a robust model for predicting median ground motion. The model integrates epistemic uncertainty through a median ground motion logic tree that accounts for variations in magnitude and distance. This study compares the GMMs selected for the 2025 PRVI NSHM, including both as-provided and regionally adjusted NGA-West2 and NGA-Subduction models, with those used in the 2003 PRVI NSHM. We evaluate how changes in model selection, weighting, aleatory variability, and epistemic uncertainty influence seismic hazard estimates. Trends with distance, magnitude, and spectral period are analyzed to evaluate how the scaling behavior of the newer GMMs differs from that of earlier models. Relative to the GMMs used in the 2003 NSHM for this region, the 2025 models generally predict lower ground motions. Comparisons with additional GMMs indicate that the adjustments applied for PRVI are consistent with regional-specific modifications developed elsewhere globally. The increase in aleatory variability and epistemic uncertainty in the 2025 update results in a notable increase in hazard levels from these wider uncertainty bounds. These changes can result in as much as a 10%–20% variation in probabilistic ground motion at the 2% in 50 years exceedance level for hazard maps computed across the region for representative site classes and periods.
Landslides can dam rivers and require rapid response to mitigate catastrophic outburst floods. Here, we present a workflow to map landslide dam formation susceptibility at a regional scale. We define a probabilistic function that combines river valley width and landslide volume to efficiently determine the likelihood of a landslide dam or “damability”. We combine damability values with landslide susceptibility to estimate landslide dam susceptibility. The valley width measurements are automated using a new elevation threshold-based algorithm. Landslide volume is represented as a statistical distribution from mapped landslides. We validate and apply our approach to the Oregon Coast Range, USA and find that 36 % of river stretches exceed a dam potential threshold; these are in river headwaters and steeper terrain, which in this case correlate with more resistant lithologies. We also estimate volumes of the potential dammed lakes and find that most rivers with high dam susceptibility are less likely to impound large lakes because they have low drainage areas. However, widespread susceptibility, and the potential impacts from exceptionally large landslides, suggest that this hazard should be considered in the Pacific Northwest. The damability function workflow can ingest new data and be applied more broadly to assess future landslide dam hazards.
Emergency assessments of postfire debris-flow hazards that are performed by the U.S. Geological Survey (USGS) provide estimates of debris-flow likelihood and rainfall triggering conditions that are used for evaluating and managing runoff-generated debris-flow hazards in recently burned areas throughout the western United States. Although the immediate postfire period, within roughly one year after fire, is typically the most susceptible to runoff-generated debris flows, the hazard evolves in time and space as the burned area recovers. The recovery trajectory a given burned area will take depends on local climate and weather and can be difficult to predict. Some burned areas recover quickly, whereas others experience debris flows for multiple years after fire. As a result, extending our ability to update debris-flow likelihood estimates and rainfall thresholds based on observed recovery of the burned area would be beneficial. We present a method for multi-year runoff-generated debris-flow hazard assessment that leverages the USGS “M1” debris-flow likelihood model and integrates updated, satellite-derived, normalized burn ratio data to estimate vegetation recovery. We predict recovery-aware rainfall thresholds and validate them against a multi-year debris-flow hazard prediction and could be adapted for use with other debris-flow models that incorporate burn severity data.
It has been previously proposed that a megasplay fault within the Cascadia accretionary wedge, spanning from offshore Vancouver Island to Oregon, has the potential to slip during a future Cascadia subduction zone earthquake. This hypothetical fault has major implications for tsunami size and arrival times and is included in disaster-planning scenarios currently in use in the region. This hypothesis is evaluated in this study using CASIE21 deep-penetrating and U.S. Geological Survey high-resolution seismic reflection profiles. We map changes in wedge structural style and seismic character to identify the inner-outer wedge transition zone where a megasplay fault has been previously hypothesized to exist and evaluate evidence for active faulting within this zone. Our results indicate that there is not an active, through-going megasplay fault in Cascadia, but instead, the structure and activity of faulting at the inner-outer wedge transition zone is highly variable and segmented along strike, consistent with the segmentation of other physical and mechanical properties in Cascadia. Wedge sedimentation, plate dip, and subducting topography are proposed to play a major role in controlling megasplay fault development and evolution. Incorporating updated megasplay fault location, geometry, and activity into modeling of Cascadia earthquakes and tsunamis could help better constrain associated hazards.