USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Petrology”

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,659 records · Page 79Linked to original sources
USGS2025

An unexplained tsunami: Was there megathrust slip during the 2020 Mw7.6 Sand Point, Alaska, earthquake?

On October 19, 2020, the M w 7.6 Sand Point earthquake struck south of the Shumagin Islands in Alaska. Moment tensors indicate the earthquake was primarily strike-slip, yet the event produced an enigmatic tsunami that was larger and more widespread than expected for an earthquake of that magnitude and mechanism. Using a suite of hydrodynamic, seismic, and geodetic modeling techniques, we explore plausible causes of the tsunami. We find that strike-slip models consistent with the moment tensor orientation cannot produce the observed tsunami. Hydrodynamic inversion of sea surface deformation from deep ocean and tide gauge data suggest seafloor deformation more closely matches a megathrust, rather than a strike-slip, source. Static slip inversions, using sea level and Global Navigation Satellite System data, allow for a portion of co-seismic megathrust slip that can explain tsunamigenesis. Combining all available geophysical datasets to model the kinematic rupture, we show that considerable, relatively slow, megathrust slip is allowable in the Shumagin segment, concurrent with strike-slip faulting. We hypothesize that the slow megathrust rupture does not contribute much seismic radiation allowing it to previously go unnoticed with traditional seismic monitoring.

Alaska↗
USGS2025

Hydrogeology and groundwater quality in the Snake River alluvial aquifer at Jackson Hole Airport, Wyoming, 2011–20

The Snake River alluvial aquifer underlying the Jackson Hole Airport (JHA) in northwest Wyoming is an important source of water used for domestic, commercial, and irrigation purposes by the airport and nearby residents. The U.S. Geological Survey, in response to previously identified water-quality concerns in the area, monitored and evaluated changes in hydrogeologic characteristics and groundwater-quality conditions of the alluvial aquifer during 2011–20. During that period, the Jackson Hole Airport made several changes that potentially improved water quality at and downgradient from the airport. Well, water level, and hydrogeologic data were collected from the alluvial aquifer to identify hydrogeologic characteristic and groundwater quality changes. Additionally, results of statistical tests were applied to water-quality results to evaluate trends in selected physical properties and constituent concentrations with time. The trends of those data show that water quality did improve overall during the study period compared to previously collected data. Presumably, these trends are in response to the changes in the aircraft deicing/anti-icing fluid (ADAF) formulation used by the JHA, the many JHA infrastructure improvements made during 2011–20, the degradation of existing ADAFs in subsurface soils and groundwater, or some combination of these possibilities.

Wyoming↗
USGS2024

Using surficial geologic maps, vegetation, and monitoring to address erosion impacts from grazing in Channel Islands National Park, California

Employing a map-unit classification scheme based on geomorphic process and age, the U.S. Geological Survey (USGS) mapped Quaternary surficial deposits of the five islands comprising Channel Islands National Park (CHIS), California, as no such maps previously existed. Mapping was motivated through an agreement with the National Park Service (NPS) to aid natural resource assessments, including post-grazing disturbance recovery and identification of mass wasting and tectonic hazards. The resulting detailed (1:12,000 scale) maps portray areas of upland erosional transport processes and alluvial, fluvial, eolian, beach, marine terrace, mass wasting, and mixed depositional processes at a scale commensurate with the process signal on the landscape. Detailed, GPS-located observations of sedimentology, pedogenic development, landscape position, and geomorphology constrain map unit identification. In addition, map-unit boundary delineation was determined in a GIS using a variety of high-resolution (sub-meter) aerial imagery and lidar-based DEMs complemented by derivative raster products. The GIS geodatabase schema used is GeMS compliant. Resolvable map-unit widths are as small as 10-15 meter. Locational accuracy of unit boundaries is estimated to be as low as 10 m or better. Detailed mapping is at a scale allowing determination of geomorphic transport processes operating on the landscape and conveys geologic information to land managers at a resolution commensurate with park resources. Map unit ages are constrained by relative soil chronostratigraphy, vertebrate and invertebrate fossils, and materials collected for dating control using radiogenic, luminescence, and cosmogenic radionuclide techniques. Associated topical studies examined connections between geology, surface water, and vegetation distribution and restoration.

California↗
USGS2025

U.S. Geological Survey National Groundwater Climate Response Network

What is the U.S. Geological Survey National Groundwater Climate Response Network? The U.S. Geological Survey’s (USGS) National Water Monitoring Network is a network of networks that includes the Climate Response Network (CRN). The CRN is a network of wells selected to monitor the effects of climate variability, including droughts, on groundwater systems within the United States, Puerto Rico, and the U.S. Virgin Islands. The primary purpose of the CRN is to monitor these effects on groundwater levels in unconfined or near-surface confined aquifers that are minimally affected by pumping or other human-induced stresses.

Fact Sheet↗
USGS2024

Stream nitrate dynamics driven primarily by discharge and watershed physical and soil characteristics at intensively monitored sites: Insights from deep learning

We developed a suite of models using deep learning to make hindcast predictions of the 7‐day average backward‐looking nitrate concentration at 46 predominantly agricultural sites across the midwestern and eastern United States. The models used daily observations of discharge and meteorological variables and watershed attributes describing anthropogenic modification to hydrology, nitrogen application, climate, groundwater, land use, watershed physiographic attributes, and soils. Across all sites, discharge and watershed soil and physiographic attributes showed a strong influence on model performance. Analysis of drivers across sites revealed considerable regional differences related to controlling processes such as groundwater contributions. We tested several ways to pool data across sites to develop accurate models and make the most effective use of available data. Single‐site models, in which models are trained and tested at a single location, showed generally strong predictive performance (median Kling‐Gupta Efficiency = 0.66), and accuracy at poorly performing sites could be improved by grouping sites with similar characteristics. Developing a single model for all sites reduced performance at several locations with distinct characteristics, suggesting that there is a threshold of dissimilarity beyond which more data does not improve the model. While many deep learning studies have shown that national or even global models can outperform local models, it is not clear that this is true for water quality constituents. This study demonstrates how data can be combined effectively, using deep learning to develop accurate and interpretable models of instream nitrate at sites where varying processes are responsible for changes in nitrate concentration.

Water Resources Research↗
USGS2025

Induced earthquakes are generally not tidally triggered in Oklahoma and Kansas

Human-induced earthquakes occur along critically stressed faults as injected wastewater simultaneously heightens fluid pressure and pushes faults to failure. We investigate the possibility that small stresses imposed by Earth tides could trigger earthquakes in the induced seismicity region of Oklahoma and Kansas from 2011 to 2018. We decluster a catalog consisting of ∼110,000 earthquakes using three methods (Reasenberg, nearest-neighbor distance, and phase-bin). We find no significant tidal earthquake triggering using Schuster's p -value test for the declustered catalogs as a whole. We search for localized triggering using discretized space-time cells and find ∼0–6% of cells have significant tidal triggering which is close to what is randomly expected (5%) and indicates there is an insignificant amount of tidal triggering for the full study region. One area that has significant p -values across multiple time windows, ∼2014–2016 is ∼15 km from a region of large wastewater injection volume. It is possible that localized tidal triggering occurs for this time and area because faults remain critically stressed and are particularly susceptible to slip under the small stress load from the semidiurnal tide. Possible explanations for the lack of tidal triggering in our broader study are that the pre-seismic stressing rate in the earthquake nucleation area is faster than the tidal stressing rate (∼3 kPa/day), faults are not close enough to critically stressed to be affected by tidal forcing, and that nucleation occurs over longer periods than the tides considered in this study (∼1, ∼14 days). Fluid injection could be the source of a higher pre-seismic stress rate.

Kansas, Oklahoma↗
USGS2024

Seismicity and anisotropic imaging reveal an active detachment beneath the northern Alaska Range foothills

North of the Denali Fault, the collision between the Yakutat block with North America is accommodated by a fold-thrust belt giving rise to the northern Alaska Range foothills. At the western end, the Kantishna Hills anticline hosts prominent microseismicity and surface deformation, interpreted as active folding of the Kantishna Hills anticline above a midcrustal detachment. We test for this detachment by using anisotropy-aware receiver functions to image fabric contrasts within the crust in context with seismicity. Seismic stations near the crest of the Kantishna Hills anticline and near its southern flank show a single strong contrast in dipping fabric at depths of 12–13 km near microseismicity clustering depths, consistent with a detachment plane beneath the fold. A minimum b -value at 10–13 km depth is consistent with seismicity on the detachment, compatible with the imaged anisotropic contrast, while off-fault seismicity is shallower and deeper with smaller magnitudes. South-dipping imbricate thrusts in schist characterize the northern Alaska Range foothills structure. This supports our interpretation of the observed anisotropy as reflecting SSW–SSE-dipping foliation above a detachment at ∼10–13 km depth that exploits existing crustal weaknesses along subtle fabric contrasts observed in the seismically quiescent region north of the actively deforming belt.

Alaska↗
USGS2024

Subducting plate structure and megathrust morphology from deep seismic imaging linked to earthquake rupture segmentation at Cascadia

The origin of rupture segmentation along subduction zone megathrusts and linkages to the structural evolution of the subduction zone are poorly understood. Here, regional-scale seismic imaging of the Cascadia margin is used to characterize the megathrust spanning ~900 km from Vancouver Island to the California border, across the seismogenic zone to a few tens of kilometers from the coast. Discrete domains in lower plate geometry and sediment underthrusting are identified, not evident in prior regional plate models, which align with changes in lithology and structure of the upper plate and interpreted paleo-rupture patches. Strike-slip faults in the lower plate associated with oblique subduction mark boundaries between regions of distinct lower plate geometry. Their formation may be linked to changes in upper plate structure across long-lived upper plate faults. The Juan de Fuca plate is fragmenting within the seismogenic zone at Cascadia as the young plate bends beneath the heterogeneous upper plate resulting in structural domains that coincide with paleo-rupture segmentation.

Oregon, Washington↗
USGS2024

The U.S. Geological Survey Ohio Water Microbiology Laboratory

The U.S. Geological Survey Ohio Water Microbiology Laboratory is a part of the Ohio-Kentucky-Indiana Water Science Center. The mission of the laboratory is to provide microbiological data of public health significance from surface waters, groundwaters, and sediments for a variety of study objectives. The laboratory conducts internal projects, works with external cooperators, and assists U.S. Geological Survey offices and National programs. The laboratory offers guidance, study design, and data interpretation expertise to collaborators, all following rigorous quality control and quality assurance procedures.

Fact Sheet↗
USGS2025

Light absorbing particles deposited to snow cover across the Upper Colorado River Basin, Colorado Rocky Mountains, 2013-16: Interannual variations from multiple natural and anthropogenic sources

Atmospheric particulate matter (PM) as light-absorbing particles (LAPs) deposited to snow cover can result in early onset and rapid snow melting, challenging management of downstream water resources. We identified LAPs in 38 snow samples (water years 2013–2016) from the mountainous Upper Colorado River basin by comparing among laboratory-measured spectral reflectance, chemical, physical, and magnetic properties. Dust sample reflectance, averaged over the wavelength range of 0.35–2.50 μm, varied by a factor of 1.9 (range, 0.2300–0.4444) and was suppressed mainly by three components: (a) carbonaceous matter measured as total organic carbon (1.6–22.5 wt. %) including inferred black carbon, natural organic matter, and carbon-based synthetic, black road-tire-wear particles, (b) dark rock and mineral particles, indicated by amounts of magnetite (0.11–0.37 wt. %) as their proxy, and (c) ferric oxide minerals identified by reflectance spectroscopy and magnetic properties. Fundamental compositional differences were associated with different iron oxide groups defined by dominant hematite, goethite, or magnetite. These differences in iron oxide mineralogy are attributed to temporally varying source-area contributions implying strong interannual changes in regional source behavior, dust-storm frequency, and (or) transport tracks. Observations of dust-storm activity in the western U.S. and particle-size averages for all samples (median, 25 μm) indicated that regional dust from deserts dominated mineral-dust masses. Fugitive contaminants, nevertheless, contributed important amounts of LAPs from many types of anthropogenic sources.

Colorado↗
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↗
USGS2025

A method to obtain remotely sensed grain size distributions from nonplanar granular deposits

Constraining the grain size distribution of granular deposits with complex surfaces is difficult with existing approaches. Field and laboratory techniques are time consuming and limited by the maximum grain size that laboratories can accommodate. In this study, we present a new method to identify the coarse fraction of the grain size distribution at a debris-flow fan deposit surveyed with terrestrial laser scanning (TLS) in Glenwood Canyon, Colorado, USA. This method is a novel grain segmentation algorithm developed for application to point cloud data of deposits with complex surfaces and angular grains ranging in size from centimeters to a meter. This approach combines an existing random forest machine learning method with a novel iterative clustering algorithm. We compared the grain size distribution from our algorithm with a Wolman pebble count conducted in the field, and found a root mean squared error of less than 2 cm from the 5th to 95th percentile of the grain size distribution of grains ranging from cobble to boulder sized (6.3–78 cm in our application). Finally, we compared our new algorithm with an existing open-source grain segregation algorithm, and our method outperformed the selected alternative when applied to the debris-flow deposit point cloud.

Colorado↗
USGS2025

Dynamic feedbacks between river meandering and landsliding in northwestern Washington glacial terraces

Landsliding in river valleys poses unique risks for cascading hazards and can damage infrastructure and cause fatalities. In postglacial valleys, many landslides are posited to occur in relation to lateral river erosion, but the dynamics of fluvial-hillslope interactions are not well understood. Here, we investigate a section of the Nooksack River in western Washington State where the channel is flanked by landslide-prone glacial terraces similar to those that failed in the 2014 State Route 530 “Oso” landslide. We map 216 landslides through time across 17 aerial imagery data sets (1933–2022) and analyze them in relation to river meandering and curvature. We observe dynamic feedbacks between lateral river meandering and valley-adjacent landsliding. Terrace lateral retreat rates of up to 25 m/year owing to combined fluvial erosion and slope failure occur on pinned, outer meander bends immediately downstream from peaks in river curvature (>0.0075 1/m); these locations are predisposed to both shallow and deep-seated landslides. Deep-seated landslides extending 17%–32% of the active valley width into the floodplain can displace the river away from the floodplain margin and change the channel planform. River-displacing landslides relocate meanders up- or downstream, thereby conditioning the location of subsequent landslides. This conceptual model of coupled landslide-driven meander displacement and valley-adjacent landsliding is exemplified across western Washington river systems. The distance between up- and downstream valley-adjacent landsliding scales with valley width, meander wavelength, and terrace height. Our results can advance our understanding of the river-hillslope interface in landscape evolution and can be used to inform hazard management in river corridors.

Washington↗
USGS2024

Coastal Science Navigator companion guide—Discover the U.S. Geological Survey coastal science products you need

The Coastal Science Navigator is an online gateway to a wide variety of U.S. Geological Survey (USGS) coastal change hazards-related information, data, and tools relevant to stakeholders’ scientific and decision-making needs. The products within the Coastal Science Navigator provide data related to past, present, and future threats to our coastlines. The filter search allows users to see all available products and identify relevant options by time scale, geographic scope, coastal hazard theme, and other filters. The guided search suggests products based on users’ answers to a short series of questions. A comprehensive summary is available for each product. The idea for the Coastal Science Navigator arose in 2020 in response to stakeholder feedback identifying the need for a central source for USGS coastal science information. It was published in July 2023 and initially included 55 products. Regular updates are planned to integrate other existing and new products. This guide introduces some of the many coastal change hazards-related products available through the USGS. In it, we showcase the products included in the Coastal Science Navigator’s initial publication in July 2023. While it is not representative of all the information, tools, and data available, we hope it serves as a compelling snapshot of what the USGS has to offer and encourages you to explore the Coastal Science Navigator to discover more of the products you need. To navigate this guide, the products have been organized by the time scale they are best suited for—past, present, or future—although many products cover multiple time scales. An additional section features software, one of the many product types available as filters within the Coastal Science Navigator. Other products include downloadable data, websites, and geonarratives (web pages that combine text, images, and interactive maps into narratives you can scroll through). Featured geographic scopes are also highlighted within this guide, detailing some of the many regions in which the USGS conducts research and illustrating another way to filter products within the Coastal Science Navigator.

Circular↗
USGS2026

Magnitude conversion relations create substantial differences in seismic hazard models

Earthquake catalogs are essential data inputs for seismic hazard modeling. Because earthquake magnitudes are reported in a variety of types (e.g., local magnitudes and moment magnitudes), magnitude conversion relationships must be used to convert the different magnitude types present in a catalog to a uniform magnitude type to avoid biases in the hazard computation. However, these conversion relationships are often uncertain and have been shown to sometimes perform poorly. Here, we investigate the sensitivity of the gridded seismicity component of the National Seismic Hazard Model (NSHM) to the catalog conversion equations in the Eastern United States. In the 2023 NSHM, magnitudes of various types were converted to moment magnitudes using equations developed by the Central and Eastern United States Seismic Source Characterization for Nuclear Facilities (CEUS‐SSCn), based on least‐squares (LS) regressions made using data from a catalog containing events up through 2008. We recompute these equations using events in the Advanced National Seismic System Comprehensive Earthquake Catalog with multiple magnitudes from 2000 to 2023. Although we prefer the use of orthogonal regressions for our datasets, LS regressions produce broadly similar results, with both approaches exhibiting large deviations from the CEUS‐SSCn conversions, especially at smaller magnitudes. We compare the spatial distribution of annual rates using three different models: (1) the 2023 NSHM conversions, (2) our updated conversions, and (3) no conversions. We find that the choice of conversions leads to substantial differences in the rate forecasts, which can greatly impact the seismic hazard model, particularly in regions with low‐seismicity rates such as the Eastern United States, where the hazard is dominated by gridded seismicity rather than a fault model.

Seismological Research Letters↗
USGS2026

Channel morphology and large wood control postfire debris-flow erosion and deposition

Runoff-generated debris flows are a known response to wildfire, and accurately predicting the volume of these debris flows is important for estimating the magnitude of downstream hazards. Prior data collection efforts have focused on debris-flow volume measurements at catchment outlets, but few studies have considered how erosion and deposition modulate the volume of sediment arriving at catchment outlets. This study takes advantage of a high-resolution dataset to document the factors that control the total debris-flow volume reaching the catchment outlet during a fatal postfire debris flow. Using pre- and post-event airborne lidar, satellite imagery and field mapping, we found that a postfire debris flow in the Black Hollow catchment in northern Colorado eroded 136,000 ± 30,000 m 3 and redeposited 27,000 ± 7,500 m 3 in the main channel. Most of the in-channel deposition (52% by volume) occurred where a confined channel reach transitioned to an unconfined channel reach downstream, allowing the flow to widen and deposit material. Wood jams played multiple roles in the debris-flow dynamics, both nucleating deposition (25% of the deposit volume was stored behind wood jams) and exacerbating erosion (50% of the total erosion occurred downstream from a wood dam break). The remaining deposition occurred due to spatial changes in channel slope as well as deposition observed at newly formed channel bars. Using these data in this study, we identified topographic and vegetation metrics that can be used (pre-event) to anticipate where deposition may occur in channels prior to a debris flow.

Colorado↗
USGS2025

Estimating agricultural irrigation water consumption for the High Plains aquifer region with integrated energy- and water-balance evapotranspiration modeling approaches

Estimation of irrigation water use provides essential information for the management and conservation of agricultural water resources. Conventionally, water use data are created based on reports and surveys from water users, whereas manual records may not be complete due to lacking flow meters, measurement gaps, inconsistent methods across regions, and time- and cost-consuming data processing. Alternatively, spatially explicit estimation of irrigation water use can be conducted efficiently using remote sensing evapotranspiration (ET) modeling approaches. In this study, we created a gridded blue water evapotranspiration (BWET) dataset to estimate historical irrigation water consumption (1986 – 2020) in the croplands across the United States High Plains aquifer region. The BWET data were generated by integrating an energy-balance ET model [Operational Simplified Surface Energy Balance model (SSEBop)] and a water-balance ET model [Vegetation ET model (VegET)]. BWET in croplands indicates crop consumptive use of irrigation water extracted from surface water and groundwater resources. The BWET estimates were compared with reported irrigation water use data for all counties within the aquifer region. The results revealed high agreement between growing season (May – September) BWET and annual water withdrawal at county level. Specifically, correlation coefficients of volumetric BWET and water withdrawal were 0.90 and 0.96, respectively, for the entire aquifer region and western Kansas. The timeseries of BWET and water withdrawal showed similar temporal trends and high covariations. The BWET estimates were systematically lower than the water withdrawal measurements, which was primarily attributed to blue water losses in the irrigation system. The irrigation efficiency, calculated as the ratio of BWET to water withdrawal depth, was 0.57 and 0.74 for the entire aquifer region and western Kansas, respectively. This study demonstrates the capability of using satellite-based ET models (e.g., SSEBop and VegET) to efficiently estimate crop water consumption and evaluate irrigation efficiency at landscape, county, and regional scales.

Colorado, Kansas, Nebraska, New Mexico, Oklahoma, ↗
USGS2025

Hydraulic conductivity and transmissivity estimates from slug tests in wells within the Mississippi Alluvial Plain, Arkansas and Mississippi, 2020

During the spring and summer of 2020, the U.S. Geological Survey conducted single-well slug tests on selected observation wells within the Mississippi Alluvial Plain in Arkansas and Mississippi to estimate hydraulic conductivity and transmissivity values for the Mississippi River Valley alluvial and middle Claiborne aquifers. Well and aquifer data were collected from field measurements, well-construction reports, and published aquifer-thickness information. A total of 324 slug-in and slug-out tests were conducted on 48 wells by using mechanical slugs to displace the water column and submersible pressure transducers to record changes in water levels in the wells. Hydraulic conductivity of the aquifers in which the wells are screened was estimated by curve fitting the water-level-change data using aquifer test analysis software. Estimates of aquifer transmissivity were made by multiplying the estimated hydraulic conductivity value by the aquifer thickness at well locations. Mean hydraulic conductivity estimates for 44 observation wells screened in the Mississippi River Valley alluvial aquifer range from 3 to 401 feet per day, and mean transmissivity estimates range from 285 to 80,559 feet squared per day. Mean hydraulic conductivity estimates for four observation wells screened in units of the middle Claiborne aquifer range from 0.14 to 183 feet per day, and mean transmissivity estimates range from 55 to 67,913 feet squared per day. The results from these tests can be used to improve the understanding of water availability and groundwater migration, to refine groundwater models, and to ultimately provide stakeholders and decisionmakers better information for management of the groundwater resources within the Mississippi Alluvial Plain.

Arkansas, Mississippi↗