USGS Science⌕ Search

SEARCH · USGS Science

Results for “Applied and Environmental 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,657 records · Page 81Linked to original sources
USGS2025

Groundwater response to managed aquifer recharge at the Southeast Houghton Artificial Recharge Project in Tucson, Arizona

Managed aquifer recharge is a widespread practice for storing water in the subsurface as groundwater. At a managed aquifer recharge facility in southern Arizona, groundwater-level and repeat microgravity data were collected to monitor aquifer response. These data were used to inform parameter identification for an unsaturated-zone flow model used to simulate the recharge process. The facility, the Southeast Houghton Artificial Recharge Project (SHARP), consists of 3 surface basins (about 27,600 square meters [6.8 acres] total surface area) where recycled water is distributed in recharge cycles lasting several months, with dry periods in between. During the study period, December 2020–December 2022, Tucson Water (the City of Tucson’s water utility) reported 6.56×10 6 cubic meters of water (5,320 acre-feet) recharged. Monitoring included groundwater-level observations at 3 monitoring wells and repeat microgravity measurements at as many as 22 locations (some stations were destroyed between surveys). Six gravity surveys were carried out using absolute- and relative-gravity meters. Large gravity increases, more than 250 microgals, were observed during the first repeat survey, 3.5 months after the start of recharge, but only in the immediate vicinity of the recharge basins. Data show that water moved downward to the water table, and storage changes in the unsaturated zone away from the facility were likely minimal. Gravity decreased at stations more than 1 kilometer from the facility, consistent with regional groundwater-level changes. Groundwater-level increases in wells adjacent to the recharge basins began 2 months after the second repeat gravity survey, and 5.5 months after recharge began. Unsaturated-zone flow modeling was carried out using software that simulates water movement and parameter estimation. Model calibration was carried out by minimizing an objective function calculated from the differences between simulated and observed groundwater levels, and between simulated and observed repeat microgravity data. Including repeat microgravity data in the objective function reduced the uncertainty in estimated parameter values for saturated hydraulic conductivity and saturated water content. Modeling indicated that the unsaturated zone between the recharge basins and the water table does not become saturated even after 685 days of simulated infiltration. This gradual wetting may account for increasing infiltration rates over time, as hydraulic conductivity increases with increasing water content. Unsaturated-zone water content decreased rapidly between recharge cycles. Model-simulated groundwater mounding extended about 1 kilometer from the center of SHARP after the 685-day period following the onset of recharge.

Arizona↗
USGS1976

Results of phytoplankton sampling at National Stream Quality Accounting Network stations in Montana, 1975 water year

Twelve National Stream Quality Accounting Network stations were operated in Montana during the 1975 water year. The network was established for the purpose of acquiring a base of hydrologic data for use by agencies engaged in water-resources planning on a national or regional scale. Among the characteristics analyzed were phytoplankton identification and cell counts. Samples consisted of composites of equal aliquots, collected at the center of each quartile of flow, using modified suspended-sediment samplers and sediment collection techniques. Identification and counting were done using the Sedgwick-Rafter cell method. Cell counts ranged from 21 cells per millilitre at Flathead River at Flathead, British Columbia to 27,000 cells per millilitre at Yellowstone River near Sidney. The class Bacillariophyceae was most abundant in both number and variety at all sampling sites. Anabaena and Aphanizomenon of the phyllum Cyanophyta were found at six stations and two stations respectively. These two genera of blue-green algae often become abundant in enriched waters resulting in nuisance conditions.

Montana↗
USGS2024

Correction to A regime shift in sediment export from a coastal watershed during a record wet winter, California: Implications for landscape response to hydroclimatic extremes

In the referenced article, the authors would like to correct text in the first paragraph on page 2571, Figure 9 and its caption. The changes reflect an error made in the processing of the rainfall intensity-duration data used to compare storms to published debris flow triggering thresholds. The correctly processed data does not change the interpretations made in the paper but does correctly indicate that the investigated storms did not exceed the rainfall intensity – duration threshold of Cannon (1988) but did significantly exceed the debris flow triggering threshold of Wieczorek (1987).

Earth Surface Processes and Landforms↗
USGS2026

The potential impact of three-dimensional distributed slip models derived from real-time GNSS data on the performance of the ShakeAlert earthquake early warning system for slab interface earthquakes

The ShakeAlert® earthquake early warning (EEW) system is designed to warn users of imminent strong ground motion with sufficient time to take protective actions. ShakeAlert currently uses three algorithms to characterize the earthquake source. One estimates the location and magnitude using the first few seconds of the P wave and, while fast, tends to underestimate magnitude for M w 7.0+ earthquakes. A second estimates the location, orientation, length, and corresponding magnitude of a line source using observed peak ground acceleration and contributes primarily to M w 5.5+ earthquakes. The third infers earthquake magnitude from peak ground displacement measured using Global Navigation Satellite System (GNSS) data and offers nonsaturating magnitudes for M w 7.0+ earthquakes. Other EEW algorithms exist that infer temporally evolving spatially variable slip on a 3D fault surface from real‐time GNSS data, information that might enable more accurate and timely alerts in the event of large‐magnitude subduction interface earthquakes. Here, we evaluate the potential contribution of one such algorithm, BEFORES ( Minson et al. , 2014 ), to improve ShakeAlert performance through a simulated real‐time implementation of Bayesian evidence‐based fault orientation and real‐time earthquake slip (BEFORES) and other ShakeAlert algorithms using data for eight M w 7.6+ earthquakes. The test results demonstrate that BEFORES can produce well‐constrained and accurate magnitude estimates as soon as or sooner than other EEW algorithms, in turn enabling it to increase the amount of warning time users receive in many cases. However, with a modified Mercalli intensity (MMI) threshold of 3.5, which is commonly used for issuing alerts, BEFORES would tend to alert large geographic regions that did not feel strong shaking (MMI 6+). This effect can be mitigated using a higher alert threshold of MMI 4.5 without negative impact on the amount of warning time obtainable with BEFORES.

Bulletin of the Seismological Society of America↗
USGS2024

Unscrambling the Proterozoic supercontinent record of northeastern Washington State, USA

The time interval from Supercontinent Nuna assembly in the late Paleoproterozoic to Supercontinent Rodinia breakup in the Neoproterozoic is considered by some geologists to comprise the “Boring Billion,” an interval possibly marked by a slowdown in plate tectonic processes. In northeastern Washington State, USA, similar to much of western Laurentia, early workers generally thought the tectonostratigraphic framework of this interval of geologic time consisted of two major sequences, the (ca. 1480–1380 Ma) Mesoproterozoic Belt Supergroup and unconformably overlying (<720 Ma) Neoproterozoic Windermere Supergroup. However, recent research indicates that strata considered by early workers as Belt Supergroup equivalents are actually younger, and a post-Belt, pre-Windermere record is present within the <1360 Ma Deer Trail Group and <760 Ma Buffalo Hump Formation. Thus, the northeastern Washington region perhaps comprises the most complete stratigraphic record of the “Boring Billion” time interval in the northwestern United States and holds important insights into global Proterozoic supercontinent tectonic processes. In light of these exciting developments, this field guide will address the early historic economic geology and original mapping of these Proterozoic sequences in the northeastern Washington region, and from that foundation explore more recent isotopic provenance data and their regional to global context. Finally, the guide will end with a discussion of remaining questions with a goal of stimulating interest in these relatively understudied, yet important, rocks.

Washington↗
USGS2024

Updated three-dimensional temperature maps for the Great Basin, USA

As part of the periodic update of the geothermal energy assessments for the USA (e.g., last update by Williams and others, 2008), a new three-dimensional temperature map has been constructed for the Great Basin, USA. Williams and DeAngelo (2011) identified uncertainty in estimates of conductive heat flow near land surface as the largest contributor to uncertainty in previously published temperature maps. The new temperature maps incorporate new conductive heat flow estimates developed by DeAngelo and others (2023). Predicted temperatures at depth are compared with representative measurements (for conductively dominated conditions), showing good agreement under relatively simple uniform conditions. Inputs included radiogenic heat production for all layers of 1.89 μW/m 3 , effective bulk thermal conductivity of 2.7 W/m/°C for all rocks underlying sedimentary basins, and a previously published (Williams and DeAngelo, 2011) empirically driven estimate of increasing thermal conductivity with depth in sedimentary sequences. The resulting three-dimensional temperature model is published in a USGS data release associated with this manuscript (Burns and others, 2023).

Arizona, California, Idaho, Nevada, Oregon, Utah↗
USGS2024

Aboveground carbon stocks across a hydrological gradient: Ghost forests to non-tidal freshwater forested wetlands

Upper estuarine forested wetlands (UEFWs) play an important role in the sequestration of atmospheric carbon (C), which is facilitated by their position at the boundary of terrestrial and maritime environments but threatened by sea level rise. This study assessed the change in aboveground C stocks along the estuarine–riverine hydrogeomorphic gradient spanning salt-impacted freshwater tidal forested wetlands to freshwater forested wetlands in seasonally tidal and nontidal landscape positions. Standing stocks of C in forested wetlands were measured along two major coastal river systems, the Winyah Bay in South Carolina and the Savannah River in Georgia (USA), replicating and expanding a previous study to allow the assessment of change over time. Aboveground C stocks on these systems averaged 172.9 Mg C ha −1 , comparable to those found in UEFWs across the globe and distinct from the terrestrial forested ecosystems they are often considered to be a part of during large-scale C inventory efforts. Groundwater salinity conditions as low as 1.3 ppt were observed in conjunction with losses of aboveground C. When viewed in context alongside expected sea level rise and corresponding saltwater intrusion estimates, these data suggest a marked decrease in aboveground C stocks in forested wetlands situated in and around tidal estuaries.

Georgia, South Carolina↗
USGS2025

3D viscoelastic models of slip-deficit rate along the Cascadia subduction zone

Interseismic deformation in the Pacific Northwest is constrained by the horizontal crustal velocity field derived from the Global Positioning System (GPS) in addition to vertical rates derived from GPS, leveling, and tide gauge measurements. Such measurements were folded in to deformation models of fault slip rates as part of the 2023 National Seismic Hazard Model (NSHM) update. Here I build upon one of the contributing models, the viscoelastic earthquake-cycle model of Pollitz [2022]. This model permits inclusion of effects of time-dependent viscoelastic relaxation within earthquake cycles (i.e., ‘ghost transients’) and laterally variable elastic and/or ductile material properties. I lever-age these capabilities to incorporate the Cascadia megathrust into Western U.S.-wide deformation models in which crustal fault slip rates are estimated simultaneously with slip deficit rates along the interplate boundary between the descending Juan de Fuca plate and North American plate. This effort includes construction of a margin-wide model of viscoelastic structure founded on the Slab 2.0 model and probes different models of the ductile properties of the surrounding oceanic asthenosphere, continental lower crust, and mantle asthenosphere. This results in new estimates of the distribution of slip deficit rate along the ∼ 1000 km long margin, highlights the importance of correcting for glacial-isostatic adjustment effects, and permits assessment of sensitivity of results to assumed ductile properties.

California, Oregon, Washington↗
USGS2025

The influence of postfire root strength decay on shallow landslide susceptibility in western Oregon

Severe wildfires can increase the potential for postfire landslides, partly due to the loss of vegetation and root reinforcement. This study investigates the multi-year dynamics of root strength following wildfire in coniferous forests of western Oregon and evaluates implications for postfire slope stability. Root tensile strength was measured through laboratory testing of 969 roots collected from 26 test pits over four postfire timeframes: unburned conditions, 4–12 months postfire (Cedar Creek, 2022), 24–35 months postfire (Holiday Farm, 2020), and 59–67 months postfire (Eagle Creek, 2017). Variables analyzed include root diameter, tensile thread strength, stiffness, progressive tensile strength with displacement, time-dependent strength and ductility, root area ratio, and cohesion. Results indicate a 50% reduction in root strength after wildfire, with the lowest values occurring around 4 years postfire. Projections indicate it could take 10–22 years following fire for root strength to recover to 70% and 90% of the prefire original root strength, respectively. This timeline indicates a multi-year window of vulnerability for shallow landslides that reaches a maximum approximately 4 years following fire. Brittle failure behavior was observed in burned root systems compared to unburned root systems, signaling a more abrupt loss of strength at yield. Simple slope stability analyses show that reduced root strength can lead to instability on more gentle slopes and under less saturation in comparison to unburned conditions. These findings highlight the critical role of root reinforcement in postfire slope stability and the long-term implications of wildfire disturbance on landslide susceptibility.

Oregon↗
USGS2024

Flood-inundation maps for the Cuyahoga River at Jaite, Ohio, 2024

Digital flood-inundation maps for a nearly 6-mile reach of the Cuyahoga River at Jaite, Ohio, were created by the U.S. Geological Survey (USGS) in cooperation with the Northeast Ohio Regional Sewer District Board of Trustees. The maps depict estimates of the extent and depth of flooding corresponding to selected water levels (stages) at USGS streamgage 04206425 on the Cuyahoga River at Jaite, Ohio. Water-surface profiles were computed for the stream reach by using a one-dimensional steady-state step-backwater model. The hydraulic model was calibrated to the current USGS streamgage data and then used to compute 15 water-surface profiles for flood stages at 1-foot intervals referenced to the streamgage datum and ranging from 6 to 20 feet, which correspond to below “action stage” to “major flood stage” as reported by the National Weather Service. The simulated water-surface profiles were then used with a geographic information system digital elevation model derived from light detection and ranging data to delineate the areas flooded at each stage. These maps, along with current stage data from the USGS streamgage and forecasted high-flow stages from the National Weather Service, can provide emergency management personnel and residents with information that is critical for flood response activities such as evacuations and road closures, as well as for postflood recovery efforts.

Ohio↗
USGS2026

Rapid hydrothermal triggering of induced seismicity at the Coso Geothermal Field

The long-term producing Coso Geothermal Field (CGF) in California operates over 100 wells tapping into a reservoir characterized by an extensive fracture network, complex fluid pathways, and regular seismic activity. Understanding the interaction between seismicity and injection can shed important light on the hydrothermal characteristics of the field. Here, we analyze 15 years of local seismic and daily operational data from the CGF, identifying a strong correlation between short-term increase in seismicity rate and seasonal volumetric and temperature variations in the reinjected fluid. Furthermore, the seismic footprint during peak injection of colder fluids reveals a near-instantaneous response up to 2 km away from the injection well, too rapid for pore pressure diffusion alone. This short-term and distant response is observed to have directional preference, indicating structural or permeability anisotropy within the reservoir. Additionally, the seismic response correlates with the initial volumetric increase of colder fluids, but also with temperature decrease during stable injection periods, suggesting thermal effects alone can play an important role in triggering distant seismicity.

California↗
USGS2026

Towards a refined definition of “tidal influence” for the coastal wetland sciences

The astronomical tide relates to lunar (and to a lesser extent solar) influences on rising and falling water bodies, and is a key control of many coastal processes. However, tides manifest along coastal margins in different ways to muddle application beyond Sir Isaac Newton’s first descriptions. Here, we argue that tide, or specifically “tidal,” has not been adequately defined for the coastal wetland sciences to facilitate contemporary interdisciplinary distinction for rapidly changing coastlines and future climate adaptation planning. We provide a brief history of the tidal concept, and we describe tides in a diversity of uncommonly recognized marginal coastal environments. Common astronomical tide manifestations are described through hydrographic examples, then expanded to include those tidal environments exposed additionally to fluvial, meteorological (wind), or groundwater influences. Our expanded definition of “tidal” is intended to be a forward-looking construct as environments will inevitably change with sea-level rise, water extraction, and anthropogenic flow alteration. Our definition of wetland tidal influence, which we broaden, includes those “wetlands affected by astronomic water level fluctuations within coastally restricted environments causing frequent or infrequent surface inundation or groundwater tidal variation that affects the biogeochemistry of the soil (e.g., ion exchange capacity, oxygen state, mineralization processes)”. Astronomic water level fluctuations are often amplified or dampened by non-astronomical events, such as wind. Expanding the definition of what is considered tidal by the scientific community is critical for improved understanding of coastal wetland function and the suite of ecosystem services they provide people and society, including climate mitigation opportunity, management, and adaptation.

Wetlands↗
USGS2024

Projections of multiple climate-related coastal hazards for the US Southeast Atlantic

Faced with accelerating sea level rise and changing ocean storm conditions, coastal communities require comprehensive assessments of climate-driven hazard impacts to inform adaptation measures. Previous studies have focused on flooding but rarely on other climate-related coastal hazards, such as subsidence, beach erosion and groundwater. Here, we project societal exposure to multiple hazards along the Southeast Atlantic coast of the United States. Assuming 1 m of sea level rise, more than 70% of the coastal residents and US \ $ 1 trillion in property are in areas projected to experience shallow and emerging groundwater, 15 times higher than daily flooding. Storms increase flooding exposure by an order of magnitude over daily flooding, which could impact up to ~50% of all coastal residents and US \ $ 770 billion in property value. The loss of up to ~80% of present-day beaches and high subsidence rates that currently affect over 1 million residents will exacerbate flooding and groundwater hazard risks.

Florida, Georgia, North Carolina, South Carolina, ↗
USGS2025

3D Converted wave reverse time migration imaging

We describe a newly developed method for recovering high-resolution images of seismic discontinuities, such as subducting slabs, in 3D. Our method makes use of converted S → P or P → S waves observed by dense arrays of seismometers to infer the locations and relative strengths of seismic discontinuities at depth in a target region. Observed direct and converted waves are backpropagated to their times of origin. The time-reversed wavefield is then separated into its constituent P and S components via the Helmholtz decomposition, and those separated wavefields are used to compute imaging functions that characterize the locations and relative strengths of seismic discontinuities. Imaging functions may be designed to use either S → P or P → S waves, so that users can target those arrivals expected to be most dominant in a given dataset. We have previously demonstrated the efficacy of our method in two dimensions, and we now present a 3D implementation of our technique which addresses the significant computational challenges posed by the size of volumetric wavefield data in three dimensions. Through a series of synthetic examples, we demonstrate that our method is capable of recovering the fine scale structure of a subducting slab given realistic station coverage and earthquake sources. We investigate optimal seismic station geometries for our technique and explore image interpretability in regions with poor data coverage. We find that linear station geometries yield more optimal, interpretable imaging functions than collections of small arrays can. We also show that our method can successfully recover bothS → P or P → S images when realistic shear earthquake sources are used, and we explore the additional computational challenges presented by the high frequency content of S waves. Our results demonstrate the potential for our technique to recover high-resolution information about subducting slabs in real-world regions, given that relatively sparse seismic arrays with only approximately 100 stations are capable of recovering interpretable imaging functions from just a few realistic earthquake sources for multiple discontinuities at significant depth in an area of approximately 400~sq~km.

Seismica↗
USGS2024

On the provenance of field reports of the 1886 Charleston, South Carolina, earthquake: A seismo-historical whodunnit

Much of what is known about the effects of the 1886 Charleston, South Carolina, earthquake throughout the epicentral region can be attributed to meticulous field investigations by an individual with training in geology and engineering, Earle Sloan ( Clendenin, 1926 ). In a recent study, Bilham and Hough (2024) undertook a detailed analysis of the effects of the earthquake on railroads in the Charleston region, drawing heavily from Sloan’s reports. This exercise identified several inconsistencies in Sloan’s field reports, including understandable measurement imprecision, inferred data entry mistakes, and transcription errors. The study also begged the question, where was Sloan at the time of the mainshock and over the following week? And to what extent did he draw from secondhand information in compiling his reports? On this question Sloan’s reports were sometimes enigmatic, lending themselves to misinterpretation in contemporaneous as well as modern interpretations. Beyond the details that were germane for, and briefly summarized by, the studies of Bilham and Hough (2023 , 2024) , in this report we don our historical seismologist caps to chronicle Sloan’s activities following the earthquake. We summarize our inferences here for the benefit of future scholars who might attempt to retrace either Sloan’s footsteps or our own. This study also serves to highlight Sloan’s singular contributions to earthquake science, which were never published separately.

South Carolina↗
USGS2026

The role of groundwater in contributing to surface water salinization in the Upper Colorado River Basin

Freshwater salinization impacts the availability of water for human use and ecosystem needs worldwide. It has been estimated that total dissolved solids (TDS) in the Colorado River Basin cause $350 million/year in damages and substantial resources are devoted to reducing TDS loading to streams. This study describes the development and application of coupled watershed models that enable TDS source tracking through the subsurface and across the landscape at a seasonal timestep for 35 years in the Upper Colorado River Basin. Results indicate that, on average, 75% of TDS loading to streams originates as baseflow, and 50% of loading is lagged in delivery by longer than one season. Snowmelt was identified as a dominant process controlling the transport of lagged TDS to streams. This approach informs when and where TDS mitigation efforts may be effective in a watershed that serves as a critical water supply for the southwestern United States.

Arizona, Colorado, New Mexico, Utah, Wyoming↗
USGS2026

Earthquake catalog for the Fairbanks region of central Alaska, 2014–2024, based on waveform cross-correlation

The Fairbanks region of central Alaska is part of a broad zone of intraplate crustal deformation, situated north of the Denali fault and north of the ongoing collision and flat‐slab subduction of the Yakutat oceanic plateau. Seismicity in the Fairbanks region occurs both in diffuse areas as well as in well‐defined lineaments, such as the left‐lateral Salcha fault, which hosted the 1937 M 8 7.3 earthquake. Starting with the regional seismicity catalog, we perform waveform cross‐correlation, network‐matched filtering, and relative relocation to obtain an enhanced seismicity catalog over the time period 2014–2024. Based on the relocated catalog, we interpret a set of 15 fault segments, including two conjugate faults and two new faults east of the previously documented fault system. Considering the combined seismicity in the Minto and Fairbanks regions, the median depth of seismicity decreases from east (6 km) to west (20 km). Our interpreted faults provide guidance for future tectonic modeling and assessment of seismic hazards in this region.

Alaska↗
USGS1999

Effects of snowmobile use on snowpack chemistry in Yellowstone National Park, 1998

Snowmobile use in Yellowstone National Park has increased substantially in the past three decades. In areas of greatest snowmobile use, elevated levels of by-products of gasoline combustion such as ammonium and benzene have been detected in snowpack samples. Annual snowpacks and snow-covered roadways trap deposition from local and regional atmospheric emissions. Snowpack samples representing most of the winter precipitation were collected at about the time of maximum annual snow accumulation at a variety of locations in the park to observe effects of a range of snowmobile traffic levels. Concentrations of organic and inorganic compounds in snow samples from pairs of sites located directly in and off snow-packed roadways used by snowmobiles were compared to concentrations in samples collected at nearby off-road sites. Concentrations of ammonium were 2 to 5 times higher for the in-road snow compared to off-road snow for each pair of sites. Thus, concentrations decreased rapidly with distance from roadways. In addition, concentrations of ammonium, nitrate, sulfate, benzene, and toluene in snow were positively correlated with snowmobile use.

Wyoming↗