USGS Science⌕ Search

SEARCH · USGS Science

Results for “Frontiers in Microbiology”

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

Fully contained laboratory earthquakes: The effect of asperity aspect ratio and free surfaces

Corner frequency ( f c ) and seismic moment ( M 0 ) are key parameters derived from seismic signals that are used to characterize earthquake stress drop, rupture area, and slip. These parameters are also affected by fault geometry and boundary conditions. However, the systematic study of these effects in laboratory settings has been challenging. This study presents laboratory earthquake experiments that examine how rupture dynamics are influenced by (a) the aspect ratio of rectangular PMMA velocity-weakening (VW) asperities surrounded by the Teflon velocity-strengthening (VS) patches, and (b) whether the sides of a VW asperity are confined with VS patches or are free surfaces. We found that increasing confinement by reducing free surfaces or increasing the VW asperity aspect ratio stabilizes fault slip, so that higher normal stress is required to transition from aseismic to seismic slip. Increased confinement and high aspect ratios also reduced M 0 and increased f c , both of which were determined from the radiated seismic waves. M 0 and f c were primarily controlled by the shorter dimension of the VW asperity. Analysis of high-frequency acoustic emission signals revealed that ruptures on high-aspect-ratio VW asperities propagated more unidirectionally, whereas ruptures on square VW asperities were more complex. Further, the high-aspect-ratio asperities were more likely to be eroded by surrounding VS regions while low-aspect-ratio asperities were more likely to rupture into the VS surroundings. These results demonstrate that both the confinement from surrounding stable areas and the geometry of the seismogenic patch can affect rupture nucleation, propagation, and seismic source characteristics.

JGR Solid Earth↗
USGS2024

From exploration to production: Understanding the development dynamics of lithium mining projects

Recently, there has been considerable recent controversy whether current and new lithium mines will be able to supply the rapidly growing needs of the electromobility transition. Mineral exploration projects are typically active for many years, and only some become operational mines. From exploration to production, the projects go through several stages of characterisation and evaluation. At each stage, decisions are made by companies and stakeholders to advance, continue or stop the project. This is a complex process, and even projects with very similar geological and technical characteristics may take very different trajectories, depending on external factors such as global market conditions and local regulatory environments. The present study investigates the dynamics of this process for lithium exploration projects. A global database of 397 lithium projects was compiled, covering their progression through major development stages between 2004 and 2022. Ordinal logistic regression was used for the statistical analysis of this data. Different explanatory variables were tested, including economic, geological, technical, and geographic factors, to identify the best predictors for project progress at each development stage. The results suggest an essential role for lithium carbonate prices, and a variable role for other factors at each stage. Critically, the already elapsed lead time and project economics, which are traditionally considered important for the prediction of the start-up of individual mines, do not appear to be relevant in all cases. The results provide important insights into the dynamics of lithium supply and may eventually allow more realistic forecasts to be made for future lithium market dynamics.

Resources Policy↗
USGS2025

Multi-lithofacies alluvial characterization via airborne electromagnetic-borehole fusion using ordinary interval kriging and geologic constraints

Airborne electromagnetic (AEM) survey provides extensive spatial coverage and detailed resolution at the near surface and can be used to develop hydrogeological models. However, utilization of AEM data is not straightforward because AEM resistivity is an indirect measurement for inferring sediment types. This study develops an ordinary interval kriging (OIK) algorithm and a resistivity-to-multi-lithofacies (R2ML) data fusion workflow for multi-lithofacies alluvial characterization. OIK utilizes irregular interval data to construct three-dimensional (3D) resistivity fields from one-dimensional inverted AEM resistivity models. The R2ML workflow maps the resistivity field generated from OIK into a multi-facies lithological model, incorporating geologic constraints derived from well logs and geological observations. The numerical and real-world cases demonstrate that OIK is computationally efficient, accounts for 3D anisotropy, and minimizes the smoothing effect, thereby preserving resistivity contrasts and reducing interpolation uncertainty. The methodology is applied to lithologic characterization of the Mississippi River Valley alluvial aquifer (MRVA) in the Shellmound area, Mississippi, U.S. A frequency-domain AEM survey was conducted to support groundwater studies for the managed aquifer recharge (MAR) to the MRVA. The resulting lithological model, including four types of lithofacies—clay, very fine sand, fine-medium sands, and graveliferous sands, illustrates the geomorphological processes of the MRVA and implies potential MAR. The alignment between the lithological model and existing geological and hydrogeological investigations demonstrates that OIK and R2ML workflow effectively capture the subsurface architecture of the MRVA. The methods have broad applicability for characterizing alluvial aquifers through AEM-borehole data fusion, supporting sustainable groundwater management.

Mississippi↗
USGS2025

Dendroseismological investigation of redwood trees along the North Coast section of the San Andreas Fault

Sequoia sempervirens (coast redwood) tree rings have the potential to annually resolve late-Holocene earthquakes on the northern San Andreas Fault based on direct (e.g., physical damage) and indirect (e.g., co-seismic environmental change) impacts, but scarcity of suitable samples and challenges crossdating this long-lived species have limited progress. More precise dating of the pre-1906 (penultimate) earthquake can improve hazard assessment and understanding of rupture segmentation. We target old trees (maximum >815 yr) along the North Coast section of the fault (increment cores via rope-climbing, 11 living trees; plunge cuts, 23 stumps) and employ complementary disturbance detection methods including radial-growth averaging (tree- and series-level), cataloging anatomical indicators (e.g., traumatic resin ducts, TRD), and dating structural components (e.g., reiterated trunks, leans). Multi-centennial ring-width chronologies at Fort Ross (1569−2023) and Gualala (1397−2023) promote continued study with incomplete crossdating limiting utilization of some series. Growth pulses (reductions, releases) and TRD dispersed across the record reflect dynamic environments that obfuscate detection of earthquake signals. The 1906 earthquake did not leave strong signatures on most trees, and when it did, within-tree response varied from normal presentation to discoloration, TRD, and missing rings. Synchrony of indicators at both locations identified 1678−1680 (6 of 15 trees) and 1698−1700 (8 of 16 trees) as the strongest disturbances among dated rings in the time range of the penultimate earthquake, peaking at 1698 (15.7 % of possible growth and anatomical indicators), but the triggering mechanisms for these events are unknown.

California↗
USGS2025

Statewide cumulative human health risk assessment of inorganics-contaminated groundwater wells, Montana, USA

Across the United States, rural residents rely on unregulated and generally unmonitored private wells for drinking water, which may pose serious health risks due to unrecognized contaminants. We assessed the nature, degree, and spatial distribution of cumulative health risks from inorganic contaminants in groundwater. Our analysis included nearly 84,000 data points from 6500+ wells, across 51 of Montana's 98 watersheds, using a public groundwater database. We compared a drinking water screening level cumulative risk assessment (CRA) for inorganics based on the U.S. Environmental Protection Agency (EPA) protective health thresholds (Maximum Contaminant Level Goals, Health Advisories [MCLG-HAs]) to a CRA based on EPA public supply enforceable Maximum Contaminant Levels (MCLs). Based on median concentrations of 19 inorganics (antimony, arsenic, barium, beryllium, boron, cadmium, chromium, copper, fluoride, manganese, molybdenum, nickel, nitrate, lead, selenium, strontium, thallium, uranium, zinc), 75% of watersheds had MCLG-HA-based cumulative risk values > 1.0; arsenic and uranium contributed the most risk, followed by strontium, fluoride, manganese and boron. Hence, this screening level (Tier I) CRA indicated widespread potential for unrecognized human health risk to private well users from inorganic contaminants considering both carcinogenic and non-carcinogenic risks. Sensitivity analysis showed that benchmarks applied (MCLG-HAs versus MCLs) exerted the largest control on results. Our findings identify priority regions for Tier 2 risk assessments to elucidate local sources and distributions of geogenic versus anthropomorphic contaminants. Our study is the first statewide assessment of cumulative health risk from groundwater that we are aware of, and results support increased statewide drinking water education and testing to reduce human health risks from contaminated private well water.

Montana↗
USGS2025

Seismic intensities, ground motions and basin amplification from the 2022 Mw 6.4 Ferndale, California intraslab earthquake

On December 20, 2022, a Mw 6.4 earthquake occurred at a depth of 18 km within the subducting Gorda plate in the Mendocino Triple Junction (MTJ), one of the most seismically active regions of the contiguous United States, causing widespread damage to local communities. Here we document the seismic intensities, ground motions, and basin amplification effects recorded by this earthquake across Humboldt County as part of an ongoing scientific effort to understand subduction zone earthquake hazards. Modified Mercalli Intensity (MMI) values from our post-earthquake field survey report shaking intensities as high as VIII (Severe). Strong ground motion data from 54 seismic stations were processed to calculate amplitude and frequency content parameters across Northern California. The maximum calculated geometric mean PGA and PGV are 1053 cm/s 2 and 52 cm/s (respectively), both recorded within the 3 km deep Eel River sedimentary basin. Comparisons with four published Ground Motion Prediction Equations indicate that PGA and PGV measurements align with expected attenuation-distance patterns for intraslab earthquakes of this nature. Within the Eel River Basin, ground motions for frequencies above 1 Hz are amplified, with respect to reference stations, by factors greater than 2. Our findings suggest that peak ground motions were mainly caused by sedimentary basin site-effects within the Eel River basin, although rupture directivity may have also increased ground motion amplitudes locally. While the Ferndale area is primarily impacted by shallow intraslab earthquakes, our results also raise questions about site-response and basin amplification hazards from a potential megathrust earthquake. The high seismicity rates of the southernmost Cascades call for stronger regional preparedness and improved strategies to mitigate the effects of such a large-scale disaster.

California↗
USGS2025

A scientific vision and roadmap for earthquake rupture forecast developments, a USGS perspective

We articulate a scientific vision and roadmap for the development of improved Earthquake Rupture Forecast models, which are one of the two main modeling components used in modern seismic hazard and risk analysis. One primary future objective is to provide fully time-dependent models that include both elastic rebound and spatiotemporal clustering nationwide, which is particularly important for shorter-term hazard and risk considerations (e.g., earthquake insurance products). We also discuss the importance and perennial challenges associated with quantifying epistemic uncertainties, including those associated with deformation-model slip rates, un-quantified sampling errors with respect to off-fault seismicity, and any spatial covariances. The need for more physics-based approaches is also emphasized, as is the benefit of adding model valuation (quantifying usefulness) to our verification and validation protocols. Given the multidisciplinary and system-level nature of this activity, modular design is critical. Future updates will also draw from best-available science by both the United States Geological Survey and the external community. The primary goal of this paper is to highlight plans that guide research and facilitate community engagement with model development, especially with respect to lowering the entry barrier for early career scientists and engineers. The paper is written so readers can focus on the sections that interest them most (see table of contents), with the Introduction and Discussion providing a stand-alone overview and summary.

Bulletin of the Seismological Society of America↗
USGS2025

Geologic framework and Holocene sand thickness offshore of Seven Mile Island, New Jersey

The U.S. Geological Survey assessed the Quaternary evolution of Seven Mile Island, New Jersey, to quantify coastal sediment availability, which is crucial for establishing sediment budgets, understanding sediment dispersal, and managing coastlines. This report presents preliminary interpretations of seismic profiles, maps of Holocene sand thickness from the shoreline to 2 kilometers offshore, and tables quantifying the volume of available sediment along the coastal margin based on data collected during 2021 and 2022. The results reveal spatial variability in the thickness and cross-shore extent of Holocene sand. The study area was separated into northern, central, and southern zones by using underlying stratigraphy and geomorphic features. The characteristics and spatial extent of the Holocene sand deposit indicate that hydrodynamic processes contribute to its spatial variability. Northern Seven Mile Island contains the thickest deposits of Holocene sand that were formed by sediment bypass around the Townsends Inlet ebb-tidal delta. Specifically, swash bars have welded to the updrift end of Seven Mile Island and have formed thick deposits of Holocene sand that thicken landward and taper seaward. Despite their thickness, these deposits have the smallest cross-shore extent; therefore, northern Seven Mile Island has the smallest volume of Holocene sand of the three geomorphic zones. Central Seven Mile Island has the thinnest Holocene sand deposits because this section of the barrier island is outside the influence of ebb-tidal deltas. Southern Seven Mile Island has the greatest volumes of Holocene sand because of increased accommodation and deposition adjacent to the Hereford Inlet ebb-tidal delta. Even though tidal inlets exert variable influence on the three geomorphic zones, sediment is distributed fairly uniformly within each geomorphic zone; each of the three zones contains 31.05–36.48 percent of the volume of available Holocene sand.

New Jersey↗
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↗
USGS2024

Fairweather transform boundary Oligocene to present orogenesis: Fairweather Range vertical extrusion and rotation of the Yakutat microplate at ca. 3 Ma

Oblique-slip along transform fault boundaries is often partitioned between a strike-slip system and thrust faults that accommodate contraction. However, topography along the Yakutat-North American transform (Fairweather fault), is asymmetric with low-terrain above active thrusts on the western, Yakutat side of the transform and high topography on the continental side with peaks >4500 m (Mount Fairweather: 4671 m) to the west of the Border Ranges fault, limited recorded earthquakes >M4, and no apparent reverse faults to generate the highest terrain. In this study we compile, for the first time, published U-Pb zircon, 40 Ar/ 39 Ar and K-Ar (hornblende, muscovite, and biotite) and U-Th/He and fission-track (zircon and apatite) bedrock ages (109) from 75 samples to investigate the exhumation history of the Fairweather Range region, complemented by a published detrital sample (ZFT and AFT) and 13 new 40 Ar/ 39 Ar (hornblende, biotite, and K-feldspar) ages on 9 bedrock samples from both sides of the Fairweather fault. Additionally, we examined published seismicity and geodetic data of the Fairweather region and assessed if plate paleo-vectors correlate with the cooling history of the Fairweather Range. Cooling age, seismic, and block-motion patterns indicate the Fairweather Range has been vertically extruded between the Fairweather and the Border Ranges faults as a coherent block since ca. 25 Ma. The pre-6 Ma Pacific plate motion (N30°W) aligns with the N33°W strike of the Fairweather Fault whereas a hypothetical pre-6 Ma Yakutat microplate paleo-vector of (N39°W) does not: indicating a post-6 Ma timing for Yakutat microplate counter-clockwise rotation (9°). We infer that rotation and impingement of the Yakutat microplate along the Fairweather fault at ca. 3 Ma led to the development of the Fairweather restraining bend and increased cooling rates. The resultant thickened Fairweather welt and the ∼30 km thick southeast end of the Yakutat microplate compounded double-indenter tectonics into Alaska's southeast convergent corner

Tectonophysics↗
USGS2025

Component identification of solid biomass fuels using reflected light microscopy: Interlaboratory study 2

As nations transition toward sustainable energy systems , biomass has become a vital component of global energy portfolios. Derived from organic materials such as wood, agricultural residues, forestry byproducts, and organic waste, biomass is a renewable energy source with significant environmental and economic benefits. Responsible biomass energy production can improve waste management, reduce emissions of greenhouse gases, and mitigate environmental pollution. However, as the diversity of biomass-derived fuels increases, robust quality assessment methods are essential to ensure their efficiency, safety, and minimal environmental impact. Reflected light microscopy (RLM) is one such technique with the potential to complement conventional physico-chemical analyses by enabling a rapid identification of material constituents and impurities. To refine this methodology and evaluate the reproducibility of solid biomass component identification using RLM, an interlaboratory study (ILS) was conducted. The study involved the recognition of 58 components across 45 photomicrographs, with the participation of 65 scientists and students from 25 countries. The participants faced high difficulty identifying some of the marked components, and as a result, the percentage of correct answers ranged from 19.0 % to 98.3 %, with an average correct identification rate of 62.7 %. The most challenging aspects of the identification process included distinguishing between woody and non-woody (agro) biomass, accurately identifying petroleum-derived materials, and differentiating agro biomass from inorganic matter. The results suggest that while RLM is an important tool for characterizing solid biomass, further development of methodology guidelines and training are necessary to enhance its effectiveness. Future research should prioritize preparing detailed, image-rich, microscopic morphological descriptions of biomass fuel components, which could improve the accuracy and reliability of using RLM in biomass fuel characterization.

International Journal of Coal Geology↗
USGS2026

Precipitation-based flood-inundation maps for the East Fork Little Blue River and tributaries at Lee’s Summit, Missouri, 2024

The U.S. Geological Survey, in cooperation with the City of Lee’s Summit, Missouri, assessed flooding of the East Fork Little Blue River and tributaries for varying precipitation magnitudes and durations, varying antecedent runoff conditions, and projected climate-change conditions. The precipitation scenarios were used to develop a library of flood-inundation maps for a 2.95-mile reach of the East Fork Little Blue River and tributaries within the city. A two-dimensional U.S. Army Corps of Engineers Hydrologic Engineering Center–River Analysis System (HEC–RAS; ver. 6.5) rain-on-grid model was calibrated to selected runoff events representing a range of antecedent runoff conditions and hydrologic responses. Lowest adjacent grades for structures within the nearby study area were incorporated into the terrain, and depth grids and water-surface elevation grids were developed for the study area. Simulated velocities at selected bridge locations were also developed from the model. The model was calibrated using water-surface elevation data collected from water-level loggers (pressure transducers) and streamflow measurements and water-surface elevation measurements made at a reference point during runoff events. The calibrated HEC–RAS model was used to simulate streamflows from design rainfall events of 15-minute to 24-hour durations and ranging from a 100- to 0.1-percent annual exceedance probability (1-year to 1,000-year recurrence intervals). Flood-inundation maps were produced for depths at a reference location of 3 to 16 feet, or a depth exceeding the 0.1-percent annual exceedance probability interval precipitation. The results of each precipitation duration-frequency value were represented by a 1-foot-increment inundation map based on the generated peak streamflow from that rainfall event and the corresponding water-surface elevation at the East Fork Little Blue River reference location. Within the HEC–RAS model, 240 scenarios were developed from the design rainfall events with each of 3 antecedent conditions. Additional scenarios were created to simulate the effects of projected precipitation scenarios on the 100-year recurrence interval, 24-hour storm and the 100-year recurrence interval, 6-hour storm. All simulation results were assigned to a flood-inundation map condition based on the generated peak flow and corresponding water-surface elevation at the East Fork Little Blue River reference location. The flood-inundation maps are shown on a web mapping application made available to the public through the City of Lee’s Summit (hyperlink will be added when available). The flood-inundation maps are tied to real-time precipitation data obtained from the Automated Surface Observing System weather station at the Lee’s Summit Municipal Airport, accessible at https://mesonet.agron.iastate.edu/request/download.phtml?network=MO_ASOS . The availability of these maps, along with information regarding observed rainfall, could help provide emergency management personnel and residents with information that is critical for flood-response activities, such as evacuations and road closures, and for postflood recovery efforts.

Missouri↗