Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Frontiers in Ecology and Evolution”

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,673 records · Page 46Linked to original sources
USGS1976

Effect of deicing chemicals on the hydrologic environment in Massachusetts; evaluation of surface-water data collection

The objective of this study is to develop predictive relationships that can be used to describe the impact of highway deicing salts on the hydrologic environment. A method, presented in an earlier report, for estimating yearly mean chloride concentrations from estimated or actual runoff and salt-application data was applied to the period 1971 through1974 and the results compared with yearly mean chloride concentrations computed from records of specific conductance. Estimates were found to be from minus 63 percent to plus 56 percent in error when the error was expressed as the difference between estimated and computed values as a percentage of the computed values. Chloride concentrations used in this study are computed from specific conductance/chloride concentration relationships and records of specific conductance. Preliminary results from graphic analyses and least-squares regression analysis show that relationships between measured values of specific conductance and chloride concentration have correlation coefficients ranging from 0.30 to 0.97. Analysis of streamflow for all major dissolved constituents is recommended with the purpose of attempting to describe the variations in the specific conductance/chloride concentration relationships.

Massachusetts↗
USGS2024

Declining groundwater storage expected to amplify mountain streamflow reductions in a warmer world

Groundwater interactions with mountain streams are often simplified in model projections, potentially leading to inaccurate estimates of streamflow response to climate change. Here, using a high-resolution, integrated hydrological model extending 400 m into the subsurface, we find groundwater an important and stable source of historical streamflow in a mountainous watershed of the Colorado River. In a warmer climate, increased forest water use is predicted to reduce groundwater recharge resulting in groundwater storage loss. Losses are expected to be most severe during dry years and cannot recover to historical levels even during simulated wet periods. Groundwater depletion substantially reduces annual streamflow with intermittent conditions predicted when precipitation is low. Expanding results across the region suggests groundwater declines will be highest in the Colorado Headwater and Gunnison basins. Our research highlights the tight coupling of vegetation and groundwater dynamics and that excluding explicit groundwater response to warming may underestimate future reductions in mountain streamflow.

Nature Water↗
USGS2025

Comparison of hydrologic data and water budgets between 2003–08 and 2018–23 for the eastern part of the Arbuckle-Simpson aquifer, south-central Oklahoma

The Arbuckle-Simpson aquifer is divided spatially into three parts (eastern, central, and western). The largest groundwater withdrawals are from the eastern part of the Arbuckle-Simpson aquifer, which provides water to approximately 39,000 people in Ada and Sulphur, Oklahoma, and surrounding areas. The Arbuckle-Simpson aquifer, including the eastern part, is designated a sole source aquifer for its service area. Based primarily on data collected between 2003 and 2008, a series of comprehensive hydrologic studies of the Arbuckle-Simpson aquifer was published to provide the information necessary to perform groundwater-flow model simulations so that the Oklahoma Water Resources Board could determine how much water could be withdrawn from the aquifer while maintaining flow to springs and streams. As part of the Phase 1 studies, an aquifer water budget was developed from a numerical model for the period 2003–08. For this report, Phase 1 refers to the 2003–08 data collection period, although for some of the analyses, data collected prior to 2003 were used to inform model development work. Allocation of water from this aquifer was then established by the Oklahoma Water Resources Board in 2013. Additional well-spacing rules were also established by the Oklahoma Water Resources Board for sensitive sole source groundwater basins. To determine how the water budget for the eastern part of the Arbuckle-Simpson aquifer has changed over time, recently collected hydrologic data (2018–23) were compared to data collected during 2003–08. The analysis of changes in the aquifer water budget from 2003–08 to 2018–23 could help resource managers better understand changes in the overall balance of water in storage and the potential effects on streamflow, changes in groundwater levels, and the effects of different water uses in the aquifer area on available water in the eastern part of the Arbuckle-Simpson aquifer and streams overlying the eastern part of the Arbuckle-Simpson aquifer.

Oklahoma↗
USGS2025

Rapid fault healing from cementation controls the dynamics of deep slow slip and tremor

Despite its status as one of the most important discoveries in geophysics, the physical mechanism(s) responsible for slow slip events (SSEs) are not well understood. Here, we synthesize observations of deep SSEs in the Cascadia Subduction Zone and argue that rapid, cohesive fault strengthening may control the dynamics of deep SSEs. Cohesive strength is frequently ignored in constitutive laws used to describe fault rheology in numerical simulations of earthquakes and SSEs alike. To demonstrate its importance, we perform and analyze a suite of petrological experiments that simulate fault healing under representative pressure and temperature conditions. We show that significant cohesive strength recovery caused by dissolution-precipitation processes occurs on timescales of just a few hours. Together, our experimental and observational results support the idea that cohesion is a key component of fault strength under SSE conditions and highlight the need for its inclusion in both future experiments and numerical models of fault slip.

British Columbia, Washington↗
USGS2026

Quantitative mineral resource assessment of lithium pegmatite deposits in the southern Appalachian orogen

The first quantitative mineral resource assessment for undiscovered lithium pegmatite deposits in the southern Appalachian region of the United States was conducted. Permissive tracts for lithium pegmatite deposits were delineated by integrating lithological, tectonic, geochemical, geophysical and mineral occurrence data. Lithium pegmatite prospectivity of the tracts was ranked with simplified mappable criteria, including proximity to Paleozoic felsic intrusions and major lithotectonic structures, stream sediment geochemical anomalies, and pegmatite occurrence data. The geospatial data and permissive tracts were used to estimate the number of undiscovered lithium pegmatite deposits. These estimates were integrated into probabilistic simulations along with a new global lithium pegmatite grade and tonnage dataset to quantify potential contained undiscovered lithium resources. An economic filter was applied to convert the probabilistic estimates of contained lithium into recoverable material. The identified lithium pegmatite resources for the Carolina Lithium and Kings Mountain deposits, North Carolina, contain 1589 thousand tons (kt) of Li 2 O. The median contained undiscovered resource for the southern Appalachian orogen was estimated to be 2240 kt Li 2 O. At 90% confidence, the region contains at least 130 kt Li 2 O, and 10,700 kt at 10% confidence. After applying economic filters, the median recoverable contained resource was 1430 kt Li 2 O, corresponding to approximately 201 years of current lithium imports for consumption in the United States. North and South Carolina are likely to contain most of these resources. Coarse data resolution and intra-state variations in the geological data contribute to uncertainty of undiscovered lithium pegmatite resources. Continued efforts to harmonize disparate geospatial datasets with updated or new information can improve the accuracy and precision of estimated undiscovered lithium pegmatite resources in the study area and at broader scales.

Alabama, Georgia, Maryland, North Carolina, South ↗
USGS2026

Hydrogeologic conceptualization of the Arroyos subarea and Northeast Air Force Research Laboratory subarea of the Air Force Research Laboratory, Edwards Air Force Base, Antelope Valley, California, 2017–21

Research and testing within the Air Force Research Laboratory at Edwards Air Force Base in the southwestern Mojave Desert has been conducted since the 1950s. This testing requires large quantities of groundwater for cooling and operations involving hazardous constituents that have entered the groundwater. To better understand the potential movement of contaminants through the groundwater system, the U.S. Geological Survey entered into a cooperative agreement with the U.S. Air Force Civil Engineer Center with the purpose of (1) gaining a more detailed understanding of the hydraulic connection between the unconsolidated basin-fill sediment and the underlying weathered and crystalline bedrock and (2) investigating faults and their potential effect on groundwater flow. This study was specific to the Arroyos subarea and Northeast Air Force Research Laboratory subarea of the Air Force Research Laboratory. The refined hydrogeologic characterization was derived by evaluating surface geologic maps, a regional gravity model, new subsurface geophysical surveys, lithology data from boreholes, and groundwater-level data. Refinement of the hydrogeologic characterization of the study area involved developing a map of basin-fill sediment thickness, estimating the geometry and depth to bedrock, determining the presence of weathered bedrock material, and identifying the location of faults and subsurface structures. The thickness of basin-fill sediment was derived from a regional gravity model. Basin-fill sediment ranges from less than 20 feet (ft) thick to about 150 ft thick and increases substantially toward a deep structural basin in the northwestern part of the study area. The top of bedrock estimated from the new surface geophysical surveys is generally deeper than the top of bedrock from the regional gravity model. These differences in interpreted bedrock tops indicate that there may be a zone of weathered bedrock in the subsurface, overlying crystalline bedrock. Additional subsurface data evaluating the geometry and thickness of weathered bedrock are necessary for a comprehensive evaluation of weathered bedrock in the study area. Previously unmapped extensions of the Leuhman Fault, Spring Fault Zone, and Boron Gate Fault were identified in the new surface geophysical surveys. Evaluation of groundwater flow between bedrock and basin-fill sediment, and the effects of faults on groundwater flow, was done using groundwater-level data from 2020. The groundwater table was in crystalline bedrock throughout most of the Arroyos and Northeast Air Force Research Laboratory subareas, with a few exceptions where it was in basin-fill sediment and weathered bedrock. Basin-fill sediment, and likely weathered bedrock, was partially saturated in parts of the study area where basin-fill sediment was as thick as 150 ft. Weathered bedrock material may be present near the Spring Fault Zone, and if so, was likely saturated in 2020. The groundwater-level data near and within Spring Fault Zone indicate that the Spring Fault Zone acts as a vertical-along-strike conduit for, and as an across-strike impediment to, groundwater flow. In addition, steep groundwater-level gradients in the southeastern part of the study area indicate that a subsurface structure affecting groundwater flow could be related to the Leuhman Fault. Sufficient groundwater-level data were not available to conclusively determine if the Arroyos, Rich, and Boron Gate Faults, and their inferred extensions from the geophysical surveys, affect the movement or direction of groundwater flow.

California↗
USGS2026

Stepovers and beyond: Structural control of The Geysers geothermal system and the broader Clear Lake region

Fault geometry exerts a first-order control on geothermal systems by governing stress localization, fracture development, and permeability, yet in complex fault networks or broader shear zones, the relative influence of individual geometric features is often difficult to resolve. In the northern California Coast Ranges, The Geysers geothermal field is commonly interpreted to occur within a releasing stepover, although no single, clearly defined stepover is identified in published studies. To investigate the structural controls on The Geysers and the broader Clear Lake region, a two-dimensional elastic boundary element model is developed to evaluate spatial patterns of dilational strain associated with progressively more complete fault geometries. Model results show that dilation in the region is not controlled by a single structure but instead reflects the combined effects of multiple interacting fault elements. Three primary controls are identified: (1) opposing bends in the regional strike-slip fault system, including a releasing bend along the Maacama fault; (2) the southern fault tip of the Collayomi fault, which generates a prominent dilational lobe beneath the southern Geysers; and (3) a releasing stepover between the Collayomi fault and the Geyser Peak–Mercuryville–Big Sulphur Creek fault system, inferred to collectively behave as a right-lateral shear zone bounding the western margin of The Geysers. Predicted dilational strain magnitudes are sufficient to localize permeability between faults. These results highlight that incorporating complete fault networks and bedrock geological mapping can enhance geothermal assessments and provide a transferable framework for evaluating structurally controlled permeability in tectonically active regions.

California↗
USGS2024

Global variability of the composition and temperature at the 410-km discontinuity from receiver function analysis of dense arrays

Seismic boundaries caused by phase transitions between olivine polymorphs in Earth's mantle provide thermal and compositional markers that inform mantle dynamics. Seismic studies of the mantle transition zone often use either global averaging with sparse arrays or regional sampling from a single dense array. The intermediate approach of this study utilizes many densely spaced seismic arrays distributed around the globe. We systematically compute teleseismic P-to-S receiver functions for each seismic array and invert for the 1-D seismic velocity structure of the mantle transition zone beneath each array to facilitate a comparison between densely sampled regions. We stack 3,600 receiver functions on average at 67 arrays in total. The stack is used in a probabilistic inversion to estimate the mantle transition zone interface depths and velocities beneath each array. We focus on the 410-km discontinuity (410) because it is a prominent seismic interface that is clearly linked to a single mineral phase transition between olivine and wadsleyite. The depths and velocity contrasts of the 410 are mapped to temperatures and compositions using mineral physics constraints. The depth of the 410 ranges from ∼405–440 km, which is consistent with a ∼360 K temperature range in a dry mantle and a ∼260 K temperature range in a wet mantle (2 wt. % water). The Vs contrast across the 410 ranges from ∼2.5–8 %, which is consistent with ∼20–70 vol. % olivine composition in a dry mantle and ∼25–80 vol. % in a wet mantle. The bulk composition of the upper mantle near the 410-km discontinuity is typically considered to be well-mixed because there is no thermodynamic impediment to convection at the olivine to wadsleyite phase transition. However, the wide range of inferred olivine content from our study suggests that there are large lateral variations in the bulk composition of the upper mantle near the 410-km discontinuity.

Earth and Planetary Science Letters↗
USGS2026

Evaluation of submersible pressure transducers for streamflow monitoring in small streams

Compact streamgages requiring minimal infrastructure and equipped with submersible pressure transducers (PTs) are increasingly used to monitor small streams, yet disparate implementations obscure their accuracy under real-world conditions. This study isolated instrumentation-derived uncertainty in stage monitoring by co-locating various combinations of commercial vented and unvented PTs with seven U.S. Geological Survey (USGS) reference streamgages on small streams. Multi-year PT stage records from compact streamgages, collected and corrected following USGS protocols, were compared to concurrent reference observations. Vented PTs demonstrated an average measurement uncertainty of ±0.005 m and mean absolute percent error (MAPE) of ±0.2%. Unvented PTs exhibited higher uncertainty, averaging ±0.009 m and ±0.4% MAPE. Although both sensor types had stage errors up to ±0.3 m, 95% of vented PT and unvented PT observations were within 0.01 and 0.02 m of reference stage, respectively. Analysis of additional unvented PTs revealed stage errors of up to ±1.2% when using barometric sensors within 15 km of the in-water sensor. Propagation of stage error to discharge using reference rating models resulted in cumulative discharge MAPEs of ±4.5% for vented and ±5.5% for unvented PTs. These findings highlight PTs as practical alternatives to reference instrumentation when deployed with standardized procedures, potentially expanding access to reliable streamflow data.

conterminous United States↗
USGS2026

Top Elevation of Glacial Till and Thickness of the Big Sioux Aquifer Delineated From Electrical Resistivity Tomography Surveys Near Sioux Falls, South Dakota, 2022 and 2025

The City of Sioux Falls, South Dakota, requested the U.S. Geological Survey perform electrical resistivity surveys on three parcels of land north of the city. Electrical resistivity data were collected along a total of 22 transects during March 14–18, 2022, and November 17–21, 2025. Results from electrical resistivity surveys were used to delineate the top of glacial till deposits for the purpose of characterizing the Big Sioux aquifer near the city. Delineating geologic contacts provides important information on groundwater storage, flow dynamics, well design and placement, contaminant transport, groundwater–surface-water interactions, and regional water modeling. The top elevation of glacial till and the thickness of the Big Sioux aquifer varied among the three survey areas. The interpreted top elevation of glacial till in the North survey area decreases from east to west toward a slough, with elevations ranging from 1,403 to 1,418 feet (ft). The estimated thickness of the Big Sioux aquifer in the North survey area increased from east to west, with thicknesses ranging from 23 to 38 ft. The top elevation of glacial till in the Well 72 survey area generally decreases from northwest to southeast. Top elevations of the glacial till in the Well 72 survey area ranged from 1,400 to 1,409 ft along the southern end of transect W72_2. The estimated thickness of the Big Sioux aquifer in the Well 72 survey area was greatest along a southeast to northwest trending channel, with thicknesses ranging from 28 to 40 ft. The top elevation of glacial till in the Nose survey area generally decreases west toward the Big Sioux River. Top elevations of the glacial till in the Nose survey area ranged from 1,362 to 1,395 ft. The estimated thickness of the Big Sioux aquifer in the Nose survey area ranged from 33 to 70 ft.

South Dakota↗
USGS2025

A review of post-wildfire adaptations of surface-water-quality models: Synthesis, gaps, and opportunities

As wildfires increasingly affect water-supply watersheds, the demand for models to predict water-quality responses is increasing. This work reviews and synthesizes existing post-wildfire applications of water-quality models in the context of geographic and ecohydrological distribution, hydrologic and water-quality response process representation, model parameterization, model and input data scales, model calibration data availability, as well as calibration and performance evaluation approaches. Emphasis is placed on models that simulate water-quality output, rather than sediment and erosional response as the primary focus. Here, identified gaps and opportunities to advance the post-wildfire application of water-quality models include: 1. applying models in under-represented geographic and ecohydrologic regions, 2. simulating multiple streamflow generation mechanisms, including groundwater, with an emphasis on shifting dominant flow pathways as the landscape recovers following wildfire, 3. adding studies that include the simulation of metals, 4. incorporating more biogeochemical and in-stream processes to model applications, 5. applying finer spatial and temporal resolution of precipitation data input as well as finer spatial resolution hydrologic response units, 6. implementing fully distributed grid or element models or finer resolution response units to capture burn severity heterogeneity, 7. collecting enhanced water-quality data for model calibration and validation, 8. conducting model-intercomparison studies, and 9. developing model parameter value guidance in post-wildfire applications. These identified gaps and opportunities may assist users in deciding on key processes and approaches to consider in modeling post-wildfire water-quality conditions.

Science of the Total Environment↗
USGS2026

Color-Shaded Relief Bathymetry of Lake Chelan, Washington

In 2023, the U.S. Geological Survey (USGS) completed a 3-week field operation in Lake Chelan, Washington, the third deepest lake in the conterminous United States, to acquire sediment cores, high-resolution bathymetry data, and chirp sub-bottom profiles. The high-resolution bathymetry displayed in this map publication was collected using a Norbit Winghead i67-curved array multibeam system; raw data were cleaned and gridded (3 by 3 meters [m]) in the Teledyne CARIS software. In some locations where the slopes of Lake Chelan are near-vertical, multibeam returns were sparse, resulting in artifacts resembling holes that appear as grey dots within the color-shaded bathymetric relief. Although the 2023 bathymetric mapping of Lake Chelan resulted in similar morphometric characteristics to previous hydrographic campaigns (table 1; Kendra and Singleton, 1987; Pelletier and others, 1989), these new data provide a much more detailed understanding of the two physiographic basins that comprise Lake Chelan: the deeper (as deep as 445.8 m [fig. 1]; previously reported as 453 m) and narrower (1–2 kilometers [km] wide) Lucerne Basin to the north, and the shallower (as deep as 127.8 m [fig. 1]), and wider (2–3 km wide) Wapato Basin to the south. With this detailed imaging of Lake Chelan, we find pervasive morphological evidence of slope failure and mass-transport deposits (MTDs) in at least 89 locations within the lake and 63 subaqueous deltas (fig. 2). The deeper Lucerne Basin is dominated by steep slopes and numerous subaqueous deltas fed by steep fluvial catchments, whereas the Wapato Basin contains numerous landslide scarps, scars, and debris aprons along the more gently sloping lake margins, and abundant relict kettle structures from the Pleistocene deglaciation. Lake Chelan’s geomorphic, geophysical, and sediment data contain important records of glacial retreat and advance, regional volcanic eruptions, and earthquakes that generate ground motions severe enough to induce landslides recorded in the sediment as MTDs. The epicenter of the 1872 Washington State Earthquake (M 6.5–7.2; Bakun and others, 2002; Brocher and others, 2017; Sherrod and others, 2021) is along the Spencer Canyon Fault, near the southern end of Lake Chelan. Shaking from this earthquake was strong enough to affect local aquifers, induce liquefaction of unconsolidated sediments, produce large seiches (Brocher and others, 2018), and likely cause landslides within Lake Chelan.

Washington↗
USGS2025

Pre-eruptive characteristics of “suspect” silicic magmas in Carlin-type Au-forming systems

World-class Carlin-type Au deposits hosted in sedimentary rock were formed when profuse Eocene silicic magmatism swept across northern Nevada in response to arc migration. Carlin-type Au deposits formed along with porphyry/skarn Cu-Mo-W-Au deposits, epithermal Ag-Au deposits, and distal disseminated Ag-Au deposits. But unlike these other Au-bearing deposits that have clear associations with igneous intrusions, Carlin-type ore deposits appear to have formed distant from concealed plutons, and their origin remains controversial. Despite decades of abundant geophysical, geochronological, and geochemical studies suggesting the involvement of magmas, concrete evidence for magmatic involvement is still lacking. Consequently, the involvement of contemporaneous igneous systems remains inferred based on age, proximity, and variable isotopic, geochemical, and geophysical clues. A recent synthesis of deposit models postulates that Carlin-type Au deposits are intrusion-related, but that the causative magmas reside deeper (∼6–12 km) than in typical porphyry and peripheral systems (∼3–5 km), meaning that Carlin-type deposits are perhaps more distal expressions of igneous intrusions. We investigate a collection of “suspect” magmatic systems over a ∼7 m.y. timespan (∼41–34 Ma) that are contemporaneous with and near known Carlin-type ore deposits. We report results of a multifaceted array of in situ geochemical analyses (FTIR, EMP, SHRIMP-RG, LA-ICP-MS) of quartz-hosted melt inclusions, biotite, and quartz to better characterize the pre-eruptive characteristics of these magmas. We also report results of thermobarometry and thermodynamic phase equilibria modeling to help place constraints on magmatic reservoir depths and processes. Rather than a single “flavor” of silicic magma, we observe a surprisingly broad compositional spectrum of rhyolites, with one end of the spectrum exhibiting more arc-like (I-type) characteristics and the other end displaying more post-subduction, thick-crust extensional (A-type) characteristics. This broad compositional spectrum suggests a more complex picture of silicic crustal magmatism operating over a narrow span of time during slab rollback. Despite this spectrum, magmatic systems in this study are consistently ferroan and generally peraluminous, which we interpret as an expression of the relatively elevated geotherm at the time and incorporation of variable amounts of highly peraluminous metasedimentary crustal components. The silicic magma spectrum encompasses a range of mineralization associations, including subduction-related Cu-Mo-W-Au-Ag and post-subduction, thick-crust extensional rare-metal Mo-Sn-W-F-Be-Ag-Au, consistent with the prolific and diverse array of ore deposits that formed during this time. Carlin-type Au deposition appears to be associated with nearly the entire magmatic spectrum. This apparent indifference to silicic magma “flavor” would seem to imply that if magmas are involved in Carlin-type Au deposit genesis, they perhaps do not need to be compositionally specialized and/or possibly are only relevant as heat sources driving circulation to remobilize and redistribute metals.

Nevada↗
USGS2024

Development of USGS NSHMs: Do small changes in hazard imply small changes in risk?

One of the flagship products from the U.S. Geological Survey (USGS) is the National Seismic Hazard Model (NSHM). Since 1976, the NSHM has been periodically updated to reflect newly published earthquake science and provide probabilistic estimates of seismic hazard for the United States. During each update cycle, alternative models are deliberated, analyzed, and documented through logic trees and their corresponding logic tree branch weights. For example, the decision to modify a logic tree branch weight may be influenced by sensitivity analyses of the logic tree branches in their effects on the mean hazard. However, do small changes in traditional measures of hazard imply small changes in risk? In this study, we make use of two update cycles of the USGS NSHMs and a National Bridge Inventory (NBI) from the Federal Highway Administration (FHWA) to explore the preceding question. Specifically, we first identify geographic locations in the conterminous United States in which the change in hazard from one cycle to another is relatively small. Next, we model the seismic risk to highway bridges for these locations and for each update cycle, while simultaneously distinguishing low hazard environments from high hazard environments. These data enable quantitative analysis of how much changes in risk can be expected given small changes in hazard, investigating the importance of risk models in decision-making during development of the USGS NSHMs.

Conference Paper↗
USGS2024

Debris avalanches in the northern California Coast Range triggered by plate boundary earthquakes

Determining the timing and cause for ancient hillslope failures proves difficult in the western United States, yet critical as it ties directly into groundmotion estimates for hazardous events. This knowledge gap is important to confront as hillslope failures are candidates to be triggered by earthquakes along active plate boundaries. We identify two prehistoric, i.e., preinstrumental history, debris avalanches (3 10^6 and 6 10^6 m^3) in the Coast Ranges of northern California. These debris avalanches are well suited for studying the timing (to the exact year) and trigger as the densely forested landscape enables effective use of dendrochronology and high-resolution radiocarbon to pinpoint time of failure, and the legacy of landslide studies in the region provides context for evaluating a seismic trigger as the most likely triggering mechanism. Neither debris avalanche was triggered by the largest precipitation event of the instrumental record in 1964 CE, and the two debris avalanche sites are physiographically suited to accommodate topographic amplification of seismic shaking. Through a suite of geologic and dendrochronologic evidence, we establish the likely time of failure of the two debris avalanches as 1906 CE and 933 CE. The failure of the younger landslide coincides with the 1906 CE San Francisco earthquake on the Pacific/North American plate boundary, while the older landslide dates to 933 CE falling within the broad age range (850-966 CE) of a likely megathrust earthquake on the Cascadia subduction zone as recorded in coastal marshes in Humboldt Bay, California. The precise age on the 933 CE debris avalanche could trim the uncertainty on a Cascadia subduction zone event from 116 years to one year. Utilizing debris avalanche records from sites suitable to record seismic shaking improves understanding of plate-boundary earthquake timing and extent of shaking.

Bulletin of the Seismological Society of America↗
USGS2024

Limited preservation of strike-slip surface displacement in the geomorphic record

Offset geomorphic markers are commonly used to interpret slip history of strike-slip faults and have played an important role in forming earthquake recurrence models. These data sets are typically analyzed using cumulative probability methods to interpret average amounts of slip in past earthquakes. However, interpretation of the geomorphic record to infer surface slip history is complicated by slip variability, measurement uncertainty, and modification of offset features in the landscape. To investigate how well geomorphic data record surface slip, we use offset measurements from recent strike-slip surface ruptures ( n = 39), faults with geomorphic evidence of multiple strike-slip earthquakes ( n = 29), and synthetic slip distributions with added noise ( n > 10,000) to examine the constraints of the geomorphic record and the underlying assumptions of the cumulative offset probability distribution analysis method. We find that the geomorphic record is unlikely to resolve more than two paleo-slip distributions, except in specific cases with low slip variability, high slip-per-event, and semiarid climate. In cases where site-specific conditions allow for interpretation of more than two earthquakes, lateral extrapolation along a fault is not straightforward because on-fault displacement and distributed deformation may be spatially variable in each earthquake. We also find that average slip in modern earthquakes is adequately recovered by probability methods, but the reported prevalence of strike-slip faults with characteristic slip history is not supported by geomorphic data. We also propose updated methods to interpret slip history and construct uncertainty bounds for paleo-slip distributions.

Journal of Geophysical Research: Solid Earth↗
USGS2024

Are the horizontal-to-vertical spectral ratios of earthquakes and microtremors the same?

We consider the similarities and differences between earthquake and microtremor horizontal‐to‐vertical spectral ratios (eHVSR and mHVSR, respectively) using a dataset of 161 sites in southern California. Quantitative comparisons are made in terms of the eHVSR and mHVSR lognormal median curves, as well as the frequencies and amplitudes associated with the fundamental‐ and higher‐mode resonances where present. The results show only 58% of the eHVSR–mHVSR pairs agree in terms of their median curve and only 25% of the eHVSR–mHVSR pairs agree in terms of shared resonances, which increases to 68% if flat HVSRs are considered equivalent. Furthermore, while the shared resonances match very well in terms of frequency (root mean square error, RMSE, <0.11 Hz), the amplitudes of those resonances do not agree (RMSE >1.6). These findings demonstrate that while eHVSR and mHVSR agree at some sites, they are not equivalent at all sites. To investigate if the agreement between eHVSR and mHVSR could be related to features of the microtremor data, earthquake recordings, and/or the site conditions, three machine learning (ML) models at varying levels of interpretability are presented. The ML models—which include multivariate logistic regression, gradient‐boosted trees, and support vector machines—show only partial success at using site‐specific data to predict whether eHVSR and mHVSR will likely agree in terms of their median curve (accuracy of 78%) and number of resonances (accuracy of 84%). Therefore, we conclude that while eHVSR and mHVSR can be quite similar in terms of resonant frequencies at some sites, they are not identical at all sites. Furthermore, preliminary evidence shows that the agreement of eHVSR and mHVSR can be predicted a priori given features of the microtremor measurements, earthquake recordings, and site conditions, although a larger dataset will be necessary for developing a robust predictive model.

California↗
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↗