SEARCH · USGS Science

Results for “Science of the Total Envionrment”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

1,687 records · Page 55Linked to original sources

USGS2026

Methodology for quantitative spatial sensitivity analysis of volcanic geodetic networks

Introduction This report introduces a methodology for assessing the state of the U.S. Geological Survey Volcano Observatories’ geodetic monitoring networks that measure how volcanoes deform or change shape. This new method uses a model-based approach that considers the uniqueness of the instrument environments at each volcano. This report focuses on simplified volcanic sources, is independent of the shape or size of the volcano, or the network geometry, and thus highlights the strengths and potential vulnerabilities of each volcano’s geodetic network in an actionable visual format. This analysis can help observatories to make informed decisions about whether volcanoes have an adequate level of geodetic monitoring and indicate where improvements are needed.

Lassen Peak, Mount Shasta↗

USGS2025

The δ13C signature of dissolved organic and inorganic carbon reveals complex carbon transformations within a salt marsh

Coastal wetlands have high rates of atmospheric CO 2 uptake, which is subsequently respired back to the atmosphere, stored as organic matter within flooded, anoxic soils, or exported to the coastal ocean. Transformation of fixed carbon occurs through a variety of subsurface aerobic and anaerobic microbial processes, and results in a large inventory of dissolved carbon. Carbon source and the roles of aerobic respiration, sulfate reduction, and methane cycling were evaluated within salt marsh peat and the underlying sandy subterranean estuary. There is a large increase in dissolved inorganic carbon (DIC, 7,350 ± 3,900 μmol L −1 ), dissolved organic carbon (DOC, 1,040 ± 1,480 μmol L −1 ) and CH 4 (14.5 ± 33.3 μmol L −1 ) within the marsh porewaters compared to creek waters. Alkalinity production (5,730 ± 2,170 μeq L −1 ) and sulfate removal (1,810 ± 1,970 μmol L −1 ) indicate anaerobic respiration, however, relative contributions from the various decomposition pathways cannot be identified due to overlapping geochemical signatures. The δ 13 C of the DOC (−29.0 ± 3.7‰) and DIC (−11.2 ± 1.1‰) produced within the marsh differed from the bulk soil organic matter δ 13 C (−14.5 ± 0.2‰). We explore a variety of mechanisms that could result in co-occurring depleted δ 13 C-DOC and enriched δ 13 C-DIC compared to the bulk soil organic carbon pool and salt marsh vegetation, including selective mineralization, production of δ 13 C-depleted bacterial biomass, and methane-derived DOC. While important questions remain about carbon cycling pathways, we found evidence of a cryptic methane cycle. Alteration of the δ 13 C of carbon species complicates source attribution in solid and dissolved phases and careful consideration should be used when carbon is partitioned between in situ salt marsh production and external marine and terrestrial sources.

Massachusetts↗

USGS1996

Technique for estimating magnitude and frequency of peak flows in Maryland

A convenient and reliable technique for estimating flood magnitudes is required for effective flood-plain management and for the efficient design of bridges, culverts, embankments, and flood-protection structures. Methods are presented for estimating peak-flow magnitudes of selected frequencies, ranging from 2 to 500 years, for all nontidal drainage basins in Maryland. The methods were developed by generalized least-squares regression techniques using data from 219 gaged basins in and near Maryland. The State is divided into five hydrologic regions: the Appalachian Plateaus and Allegheny Ridges region, the Blue Ridge and Great Valley region, the Piedmont region, the Western Coastal Plain region, and the Eastern Coastal Plain region. These regions correspond to the physiographic provinces of the State, with the exceptions that (1) the Coastal Plain Province is divided into two hydrologic regions, and (2) there is no distinct hydrologic region corresponding to the Valley and Ridge Province as it is divided into its constituent Allegheny Ridges and Great Valley subdivisions. Sets of equations for calculating peak discharges based on physical basin characteristics are provided for each of the regions. Based on the peak-flow equations, methods for estimating peak flows are presented for ungaged and gaged streams in Maryland. The methods and equations are supported by generalized least-squares analysis of basin and flood-frequency characteristics data from 219 drainage basins in and near Maryland. estimates for each of the five regions are calculated using combinations of the fol-lowing basin characteristics: drainage area, forest cover, basin relief, carbonate rock coverage, storage, and runoff-curve number. Drainage area contributes to the estimate in all five study regions. Carbonate rock coverage is used only in the Blue Ridge and Great Valley region. Storage and runoff-curve number are used solely in the Eastern Coastal Plain region. All other basin characteristics are used in two or more regions. Standard errors of estimate for the regression equations range from 19 to 31 percent in the Appalachian Plateaus and Allegheny Ridges region, 34 to 47 percent in the Blue Ridge and Great Valley region, 33 to 48 percent in the Piedmont region, 45 to 64 percent in the Western Coastal Plain region, and 36 to 42 percent in the Eastern Coastal Plain region.

Maryland↗

USGS2025

Anatectic origin of Mississippian spodumene-bearing pegmatites in western Maine during orogenic plateau collapse

Spodumene pegmatites are an important lithium source, but the processes and tectonic settings in which they form are poorly understood. The Rumford pegmatite district surrounding Plumbago Mountain, western Maine, is host to numerous spodumene pegmatites, including the Plumbago North pegmatite (a world-class spodumene resource). Competing petrogenetic models for these spodumene pegmatites include (1) highly fractionated melts of the Mooselookmeguntic igneous complex and (2) anatexis. We tested these hypotheses by constraining the geologic, magmatic, metamorphic, and tectonic history of the Plumbago Mountain area with detailed geologic mapping and U-(Th)-Pb geochronology. The Silurian Rangeley Formation records initial isoclinal folding prior to, and contact-related metamorphism synchronous with, the intrusion of the 417 ± 4 Ma Plumbago Mountain pluton. Peak amphibolite facies metamorphism and crustal melting occurred during the ca. 410 to 400 Ma Acadian orogeny. Pulsed emplacement of the Mooselookmeguntic igneous complex occurred between ca. 389 and 356 Ma. Cassiterite U-Pb dates of spodumene pegmatites (333–327 Ma) are ≥23 m.y. younger than nearby granitic plutons, strongly arguing against the fractional crystallization model. Metamorphic monazite and xenotime (346–328 Ma) and 330 to 308 Ma 40 Ar/ 39 Ar hornblende dates indicate metamorphism coeval with spodumene pegmatite emplacement, supporting anatectic models. Reheating, anatexis, and spodumene pegmatite emplacement occurred during collapse of the 380 to 330 Ma Acadian orogenic plateau. Lithium enrichment may be linked to one or more stages of partial melting of metasedimentary and plutonic rocks during the formation, tenure, and collapse of the Acadian altiplano and emphasizes the role of anatexis in producing spodumene pegmatites of economic significance.

Maine↗

USGS2025

A new groundwater energy transport model for the MODFLOW hydrologic simulator

Heat transport in the subsurface is an important aspect of research related to the effects of a warming climate on ecological services (i.e., cold-water refugia); the development of geothermal resources for energy banking schemes (i.e., aquifer thermal energy storage [ATES]); and the effects of temperature on other aspects of groundwater quality, such as nutrient cycling. Historically, simulation of heat transport using the MODFLOW groundwater simulator and related codes was performed by scaling the input parameters of a solute-transport model to emulate heat transport. However, that approach required additional pre- and post-processing of input and output and could not account for the variation in effective thermal storage and transport properties during transient, unsaturated flow, for example. True heat-transport capabilities in the context of MODFLOW were first introduced in a variant called USG-Transport. More recently, a new groundwater energy-transport (GWE) model type has been added to MODFLOW 6, the core version of the MODFLOW hydrologic simulator. GWE supports the simulation of heat transport on structured or unstructured grids as well as within and between features of advanced packages that represent streams, lakes, multi-aquifer wells, and the unsaturated zone. GWE is integrated within MODFLOW 6 and is accessible through the FloPy Python package and the MODFLOW 6 application programming interface (API). An example simulation demonstrates conduction between grid cells through both the water and the solid aquifer material, including thermal bleeding from saturated overburden cells into a groundwater flow field.

Groundwater↗

USGS2026

Microtextural characteristics of adularia in banded quartz veins from the Midas low-sulfidation epithermal deposit, Nevada

High-grade ores at the Miocene Midas low-sulfidation epithermal deposit in northern Nevada are confined to crustiform quartz veins containing abundant adularia. Micro-X-ray fluorescence elemental mapping reveals that adularia is a common gangue mineral occurring in colloform bands, bands showing bladed textures, and bands with dendritic terminations. The adularia aggregates have delicate shapes and are comprised of stacked, submillimeter crystals hosted by fine-grained quartz. The textural evidence suggests that the adularia aggregates originally formed within a gel-like, noncrystalline silica matrix, which subsequently transformed into quartz. This indicates that the adularia did not precipitate in open space along the vein walls. Correlative microscopy, involving scanning electron microscopy-based automated mineralogy and optical petrography, demonstrates that bands containing abundant adularia are not the primary host to ore minerals. The ore minerals occur in different bands within the crustiform veins, implying that adularia and ore mineral precipitation did not always occur simultaneously. It is hypothesized here that fluid flow at Midas involved intermittent short-lived events of fluid flashing, causing rapid solute supersaturation in the liquid. During each flashing event, different amounts of vapor were produced along a given vein. Compositional differences between adjacent bands in the crustiform quartz veins may, therefore, be linked to variations in the amount of vapor formed during each flash event.

Nevada↗

USGS1998

Technique for simulating peak-flow hydrographs in Maryland

The efficient design and management of many bridges, culverts, embankments, and flood-protection structures may require the estimation of time-of-inundation and (or) storage of floodwater relating to such structures. These estimates can be made on the basis of information derived from the peak-flow hydrograph. Average peak-flow hydrographs corresponding to a peak discharge of specific recurrence interval can be simulated for drainage basins having drainage areas less than 500 square miles in Maryland, using a direct technique of known accuracy. The technique uses dimensionless hydrographs in conjunction with estimates of basin lagtime and instantaneous peak flow. Ordinary least-squares regression analysis was used to develop an equation for estimating basin lagtime in Maryland. Drainage area, main channel slope, forest cover, and impervious area were determined to be the significant explanatory variables necessary to estimate average basin lagtime at the 95-percent confidence interval. Qualitative variables included in the equation adequately correct for geographic bias across the State. The average standard error of prediction associated with the equation is approximated as plus or minus (+/-) 37.6 percent. Volume correction factors may be applied to the basin lagtime on the basis of a comparison between actual and estimated hydrograph volumes prior to hydrograph simulation. Three dimensionless hydrographs were developed and tested using data collected during 278 significant rainfall-runoff events at 81 stream-gaging stations distributed throughout Maryland and Delaware. The data represent a range of drainage area sizes and basin conditions. The technique was verified by applying it to the simulation of 20 peak-flow events and comparing actual and simulated hydrograph widths at 50 and 75 percent of the observed peak-flow levels. The events chosen are considered extreme in that the average recurrence interval of the selected peak flows is 130 years. The average standard errors of prediction were +/- 61 and +/- 56 percent at the 50 and 75 percent of peak-flow hydrograph widths, respectively.

Maryland↗

USGS2026

Global performance of remote sensing-based and reanalysis-driven models to estimate open water evaporation

Evaporation plays an essential role in the water cycle, influencing local and regional climates while directly impacting water availability in lakes. However, directly measuring evaporation over water bodies remains challenging due to the high costs of installing and maintaining the required in situ instrumentation. Although several remote sensing algorithms have been providing evaporation estimates, the lack of a global validation hinders our understanding of their relative uncertainties and performances across different regions. Here, we analyze the performance of a suite of models that leverage satellite data and meteorological reanalysis to estimate evaporation over lakes worldwide. We compare 3 remote sensing-based models, 1 reanalysis-driven model and 1 ensemble approach, using in situ observations from 27 lakes representing a diverse range of geographic and climatic regions. Our results demonstrate that, overall, the ensemble outperformed any individual model in terms of accuracy, with a RMSE and a bias of 1.3 and 0.3 mm day −1 , respectively. These findings highlight the benefits of using an ensemble approach to estimate open water evaporation with satellite-based models at the global scale, leveraging the unique strengths of each model. For the individual models, differences in the representation of heat storage changes and advection effects led to lower values of RMSE and bias, depending on the location and depth of the lakes. This study sets the path for future improvement of open water evaporation algorithms globally, while remote sensing techniques are proven satisfactory to monitoring of water loss in lakes globally, an essential step toward effective large-scale water resources management.

Water Resources Research↗

USGS2026

The effect of temperature, flow, density and disease on smallmouth bass (Micropterus dolomieu) early life growth rate

Individual growth rates of fish are influenced by gradually changing environmental conditions and rapid-onset events. However, there is a paucity of information on drivers of spatiotemporal variability in growth rates of fishes across large spatial extents. Using a 36-year dataset of smallmouth bass ( Micropterus dolomieu ) length-at-age data across 10 river reaches in Pennsylvania, USA ( n = 54,068 individuals), we estimated annual early life growth rates (mm year −1 ) and assessed the effects of relative abundance, summer temperature, summer streamflow rate and a disease outbreak on growth rate. We found that temperature had a positive effect on growth rate, varied spatially, and had a stronger effect in diseased reaches. There was also evidence for an impact of disease on early life growth through density-dependent mechanisms, with growth rates increasing from ∼85 to 101 mm·year −1 following disease-related decreases in abundance. This study adds to our understanding of the factors that shape early life growth, including rapid-onset events like disease which can have immediate and lasting effects on growth rate trajectory.

Pennsylvania↗

USGS2025

Preliminary field report of landslide hazards following Hurricane Helene

Executive Summary This report reflects our knowledge regarding the widespread landslide activity associated with Hurricane Helene observed during the U.S. Geological Survey’s (USGS) mission assignment to North Carolina in October 2024. The material in this report was originally prepared for the Federal Emergency Management Agency under mission assignment DR-4827-NC. The data and commentary in this report are reflective of a report provided to the Federal Emergency Management Agency (FEMA) on October 18, 2024, as well as information provided in briefings at the Buncombe County Emergency Operations Center. The report has been modified for public dissemination. This assessment was based on systematic visual examination and mapping of landslide locations from aerial and satellite imagery, visual and photographic observations from low-level helicopter overflights and conversations with local landslide experts from the North Carolina Geological Survey and Appalachian Landslide Consultants PLLC, and more than 50 years of combined landslide hazard professional experience of the mission-assigned field team. No systematic field investigations were done by the USGS. While responding to the event, the USGS did not identify any landslides that posed an immediate major threat to recovery personnel in parts of nine counties in North Carolina (Avery, Buncombe, Henderson, McDowell, Mitchell, Polk, Rutherford, Watauga, and Yancey); however, threats from renewed landslide activity may remain heightened in localized areas for months or even years. Known areas of the most abundant landslide occurrence include Bat Cave, Lake Lure, Chimney Rock, Swannanoa, Black Mountain, Fairview, steep areas in Asheville, and the Blue Ridge Parkway. The USGS shared detailed locations of known landslides with the Emergency Operations Centers. The thousands of landslide scars on hillsides and landslide deposits on flatter ground may present some threat to recovery activities. Soil and rocks will continue to erode from newly exposed landslide scars and may pose a threat to people and infrastructure who are immediately nearby. In general, the steeper and taller the landslide scar, the greater the potential threat. This threat is heightened during periods of rainfall and increases with the duration and intensity of rainstorms. Very heavy rainfall, or repeated rainfall events during short periods, could also initiate new landslides on steep slopes. Excavation of landslide deposits, particularly excavation of those deposits directly adjacent to steep slopes, may also pose a threat to nearby people and equipment. An interagency collaborative mapping effort led by the USGS that informed this assessment identified 1,155 landslide locations by the October 2024 briefings, but that number increased to 2,217 in a final reviewed version of the locations published in January 2025. Locations were mapped from satellite imagery, fixed-wing and helicopter surveys, media and social media, and field reports in the 3 weeks following the passage of the remnants of Hurricane Helene. USGS products outlined in this report are publicly available and include geotagged photographs from aerial reconnaissance, hazard models, an interactive view of mapped landslide locations, and landslide safety and education resources.

North Carolina, South Carolina, Tennessee, Virgini↗

USGS2025

Refinements to the attenuated propagation of local earthquake shaking (APPLES) ground-motion-based earthquake early warning algorithm

We refined the Attenuated ProPagation of Local Earthquake Shaking (APPLES) ground-motion-based earthquake early warning (EEW) approach, and directly compare APPLES performance with that of the source-characterization-based U.S. ShakeAlert EEW system for a suite of historical earthquakes in the U.S. West Coast and Japan. APPLES is an extension of the Propagation of Local Undamped Motion (PLUM) algorithm in which observed shaking intensity at seismic stations is used to forward-predict intensity distributions to surrounding areas using an attenuation model derived from an intensity prediction equation. We test new configuration options within APPLES, such as using the second highest estimated ground motion rather than the maximum, to better match median ground-motion observations and reduce alerts for small magnitude earthquakes, both of which are key alerting priorities within ShakeAlert. We evaluate these configurations alongside ShakeAlert by comparing the ground-motion estimation accuracy and available warning times relative to station observations and ShakeMap distributions. Our preferred APPLES configuration produces accurate ground-motion estimates and corresponds better with median observations compared to ShakeAlert’s estimates. This preferred configuration substantially reduces alert issuance for M < 5.0 earthquakes compared to the previous APPLES configuration, and alert-release criteria can further restrict alerts to primarily M ≥ 5.5 earthquakes without requiring magnitude estimation. Prioritizing matching median-observed ground motions may reduce APPLES warning times compared to configurations that were tuned to avoid missed alerts (such as those that use the maximum estimated ground motions), which can lead to shorter warning times compared to ShakeAlert for the same alert threshold. However, station-based warning time assessments demonstrate that APPLES can outperform ShakeAlert for high target thresholds. APPLES is a simple, independent EEW approach that may improve the robustness of EEW for the West Coast of the U.S.

California, Oregon, Washington↗

USGS2024

Using the horizontal-to-vertical spectral ratio method to estimate thickness of the Barry Arm landslide, Prince William Sound, Alaska

Conducting detailed investigations of large landslides is difficult, especially in the subsurface, largely due to environmental factors such as steep slopes, difficult access, and numerous objective hazards. These factors have made it challenging to accurately estimate the depth to the failure surface of the Barry Arm landslide, a large (roughly 10 8 cubic meters), deep-seated bedrock landslide in Prince William Sound, Alaska, recognized in 2019. The landslide has exhibited accelerated movement in recent years and poses a potential tsunamigenic hazard if rapid failure occurs. Failure surface depth, equivalent to landslide thickness, is a necessary metric for landslide-volume calculations and associated tsunami wave models. In this report, we used seismic noise recorded by a seismometer located on the Barry Arm landslide in Alaska to calculate the horizontal-to-vertical spectral ratio (HVSR) to investigate the site fundamental frequency ( f 0 ) and depth of the failure surface. To ensure that observed peak frequencies in the spectral ratio were related to the underlying stratigraphy (and not caused by other noise sources like nearby glaciers, topographic resonance, weather, or human activities), we also calculated HVSRs using earthquake signals, HVSRs at other seismic stations within a 2.5-kilometer radius, and a standard spectral ratio between the landslide station and other sites. We observed multiple peaks in the landslide HVSR curves at 1.5 hertz (Hz), 4–5 Hz, and 7–11 Hz. The frequencies of these peaks were consistent at the landslide site through time and across methods and were dissimilar to those identified at other seismic stations in the area, making it unlikely the peaks were caused by local noise. Directional HVSRs calculated at 15-degree intervals showed amplification of the higher frequency peaks in the direction parallel to slip, indicating two-dimensional site effects. We used the distinct frequency peaks in the seismic record to develop a 4-layer conceptual model of the landslide wherein the top of the deepest layer represents the primary failure surface, or the boundary between damaged (mobile) and undamaged material. We inverted Rayleigh wave ellipticity curves within this 4-layer configuration with constraints on S-wave velocity and layer thickness based on analogous material properties identified in the literature. This was necessary absent any site-specific subsurface S-wave velocity data. The best-fitting models indicate a mean slope-normal depth to the failure surface of 188 (±9) meters (m), with additional stratigraphic boundaries at 4 and 20 m below ground surface, potentially representing layered motion. These results agree with and improve upon ranges estimated by previous studies and can support future modeling and assessment efforts at Barry Arm.

Alaska↗

USGS2024

The advantages of electric-grade heat storage in long closed-loop wells

Though long (e.g., tens of kilometers) closed-loop geothermal wells are under consideration for the extraction of subsurface heat, these wells might also serve as an efficient energy storage mechanism for electricity generation. Using a semi-analytic model solution, the potential for electric-grade heat storage as a function of ambient temperature (e.g., corresponds to depth of loop), well length, well diameter, and flow rate is evaluated. The following simplified cases for comparison of standard closed loop operation with energy storage operation were considered: [1] constant flow rate and constant temperature (90 °C) injection to represent the standard closed-loop base-case; and [2] constant flow rate and annual cycle sinusoidal temperature (90-150 °C) to represent seasonal (summer) charging while solar resources are peak. Initial temperatures for all scenarios considered herein are a uniform 175 °C. Electric-grade heat is assumed to be delivered whenever temperatures at the extraction point exceed 90 °C. For calculation purposes only, if temperature falls below 100 °C, it is assumed that heat delivered is sub-economic, so no electricity would be produced. For all scenarios, temperatures at the extraction well asymptotically approach the flow-weighted average injection temperature, but energy storage scenarios exhibited a damped time-varying signal that diminishes in magnitude with length of the loop. The asymptotic approach depends on initial temperatures in the rock and the heat extraction rate (a function of well diameter and flowrate). This analysis demonstrates that shorter closed loops can produce more electricity over time than longer closed-loops previously proposed for electricity production over typical engineering design lifetimes (e.g., 30 years). Although only a high-temperature scenario is considered herein, rock that is initially below boiling temperature would not host a standard closed-loop resource, but injection of hot water seasonally would asymptotically heat this low-temperature system to temperatures capable of electricity production. In other words, regardless of initial temperatures, closed loops could be used to store electricity with no critical minerals in the geothermal battery.

Conference Paper↗

USGS2024

Biogeochemical and community ecology responses to the wetting of non-perennial streams

Transitions between dry and wet hydrologic states are the defining characteristic of non-perennial rivers and streams, which constitute the majority of the global river network. Although past work has focused on stream drying characteristics, there has been less focus on how hydrology, ecology and biogeochemistry respond and interact during stream wetting. Wetting mechanisms are highly variable and can range from dramatic floods and debris flows to gradual saturation by upwelling groundwater. This variation in wetting affects ecological and biogeochemical functions, including nutrient processing, sediment transport and the assembly of biotic communities. Here we synthesize evidence describing the hydrological mechanisms underpinning different types of wetting regimes, the associated biogeochemical and organismal responses, and the potential scientific and management implications for downstream ecosystems. This combined multidisciplinary understanding of wetting dynamics in non-perennial streams will be key to predicting and managing for the effects of climate change on non-perennial ecosystems.

Nature Water↗

USGS2000

Biodegradation of chlorinated ethenes at a karst site in middle Tennessee

This report presents results of field and laboratory investigations examining the biodegradation of chlorinated ethenes in a karst aquifer contaminated with trichloroethylene (TCE). The study site, located in Middle Tennessee, was selected because of the presence of TCE degradation byproducts in the karst aquifer and available site hydrologic and chlorinated-ethene information. Additional chemical, biological, and hydrologic data were gathered to evaluate whether the occurrence of TCE degradation byproducts in the karst aquifer was the result of biodegradation within the aquifer or simply transport into the aquifer. Geochemical analysis established that sulfate-reducing conditions, essential for reductive dechlorination of chlorinated solvents, existed in parts of the contaminated karst aquifer. Other areas of the aquifer fluctuated between anaerobic and aerobic conditions and contained compounds associated with cometabolism, such as ethane, methane, ammonia, and dissolved oxygen. A large, diverse bacteria population inhabits the contaminated aquifer. Bacteria known to biodegrade TCE and other chlorinated solvents, such as sulfate-reducers, methanotrophs, and ammonia-oxidizers, were identified from karst-aquifer water using the RNA-hybridization technique. Results from microcosms using raw karst-aquifer water found that aerobic cometabolism and anaerobic reductive-dechlorination degradation processes were possible when appropriate conditions were established in the microcosms. These chemical and biological results provide circumstantial evidence that several biodegradation processes are active in the aquifer. Additional site hydrologic information was developed to determine if appropriate conditions persist long enough in the karst aquifer for these biodegradation processes to be significant. Continuous monitoring devices placed in four wells during the spring of 1998 indicated that pH, specific conductance, dissolved oxygen, and oxidation-reduction potentials changed very little in areas isolated from active ground-water flow paths. These stable areas in the karst aquifer had geochemical conditions and bacteria conducive to reductive dechlorination of chlorinated ethenes. Other areas of the karst aquifer were associated with active ground-water flow paths and fluctuated between anaerobic and aerobic conditions in response to rain events. Associated with this dynamic environment were bacteria and geochemical conditions conducive to cometabolism. In summary, multiple lines of evidence developed from chemical, biological, and hydrologic data demonstrate that a variety of biodegradation processes are active in this karst aquifer.

Tennessee↗

USGS2026

Archean to Mesozoic–Cenozoic seismic crustal structure: Implications for geological and biological evolution

We use >4500 measurements of crustal structure to investigate the seismic structure of continental crust, Archean to Mesozoic–Cenozoic. The mean crustal thickness of continents, including their margins, is 36.5 km. We find that Archean, Paleoproterozoic, and Mesoproterozoic crust have similar mean crustal thickness (41 km), seismic velocities, Vp/Vs ratio, and density. Crusts of these ages span 3 Ga, from 4.0 to 1.0 Ga, and their similarity of physical properties (thickness, density, and seismic velocities) suggests that the process of crustal formation may have been similar during this time period, which covers 66% of Earth's history. Neoproterozoic crustal properties differ significantly from Mesoproterozoic and older crusts. The mean thickness of Neoproterozoic crust is 32 km, some 9 km thinner than the mean thickness of Archean, Paleoproterozoic, and Mesoproterozoic crust. A lithospheric root with a thickness of 150–200 km underlies Archean, Paleoproterozoic, and Mesoproterozoic crusts, and this root resists lithospheric rifting and crustal extension and thinning. Based on these observations, particularly the lithospheric thickness, we conclude that Archean, Paleoproterozoic, and Mesoproterozoic lithosphere are unique and together form the stable cratonic nuclei, defined as the thick (150–250 km), long-lived cores of continents. Higher mantle temperatures during the Archean, Paleoproterozoic, and Mesoproterozoic may have played a key role in the formation of the thick lithospheric roots. The second key finding is that the similarities of mean seismic properties indicate that the process of crustal formation operating in the Archean continued during the Paleoproterozoic and Mesoproterozoic. The thin (mean value < 135 km) lithospheric root beneath Neoproterozoic and younger crust may be related to the steady decrease in mantle temperature through time. Neoproterozoic and Paleozoic crust have similar physical properties, and these eras are characterized by pronounced biodiversification, including the renowned Garden of Ediacara, the Cambrian Explosion, and the Great Ordovician Biodiversification. Mesozoic–Cenozoic crust is the most diverse and reflects the current tectonic and magmatic processes of crustal formation.

Book chapter↗

USGS2024

Deep lithospheric controls on surface deformation and seismicity around the East Anatolian Fault Zone and A3 Triple Junction

The East Anatolian Fault Zone (EAFZ) is a plate-bounding strike-slip fault capable of hosting large earthquakes, as demonstrated by the extremely damaging February 2023 M w 7.8 and M w 7.7 mainshocks of the Kahramanmaraş earthquake sequence. Deformation related to this boundary, part of the Anatolia-Arabia-Africa (A 3 ) Triple Junction, is diffuse, as was shown by part of this earthquake sequence occurring on a northern splay of the EAFZ (the Sürgü-Çardak Fault Zone; SCFZ). Controls on surface deformation are commonly linked to stress in the brittle upper crust, but the complex deformation and seismicity patterns in this region may also reflect deeper processes, such as variations in the location and extent of the strong Arabian Plate lithospheric mantle. Seismic tomography indicates that the Arabian Plate underthrusts Anatolia as far north as the SCFZ and extends as far west as the central Adana Basin, coincident with a zone of relatively deep (>30 km) strike-slip seismogenesis that has produced M w >6 earthquakes. By investigating the relationship between deformation since the inception of the EAFZ (ca. 5 Ma), seismic structure, and seismicity, we infer that the SCFZ will become the future SE boundary of the Anatolian Plate as part of the evolving A 3 Triple Junction.

GSA Today↗

USGS2024

Integration of rupture directivity models for the US National Seismic Hazard Model

Several rupture directivity models (DMs) have been developed in recent years to describe the near-source spatial variations in ground motion amplitudes related to propagation of rupture along the fault. We recently organized an effort towards incorporating these directivity effects into the USGS National Seismic Hazard Model (NSHM), by first evaluating the community's work and potential methods to implement directivity adjustments into probabilistic seismic hazard analysis (PSHA). Guided by this evaluation and comparison among the considered DMs, we selected an approach that can be readily implemented into the USGS hazard software, that provides an azimuthally varying adjustment to the median ground motion and its aleatory variability. This method allows assessment of the impact on hazard levels and provides a platform to test the DM amplification predictions using a generalized coordinate system, necessary for consistent calculation of source-to-site distance terms for complex ruptures. We give examples of the directivity-related impact on hazard, progressing from a simple, hypothetical rupture, to more complex fault systems, composed of multiple rupture segments and sources. The directivity adjustments were constrained to strike-slip faulting, where DMs have good agreement. We find that rupture directivity adjustments using a simple median and aleatory adjustment approach can impact hazard both from a site perspective and on a regional scale, increasing shaking off the end of the fault trace up to 30--40\% and potentially reducing it for sites along strike. Statewide hazard maps of California show that the change in shaking along major faults can be a factor to consider for assessing long-period (>ls) near-source effects within the USGS NSHM going forward, reaching up to 10--20\%. Finally, we suggest consideration of minimum parameter ranges and baseline requirements as future DMs are developed to minimize single approach adaptations, to enable more consistent application within both ground motion and hazard studies.

Earthquake Spectra↗