USGS Science⌕ Search

SEARCH · USGS Science

Results for “Environmental Modelling & Software”

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.

834 records · Page 33Linked to original sources
USGS2025

Effect of land cover type on 3D deformation recovery from synthetically deformed high-resolution satellite optical imagery

The limits of detection for earthquake surface deformation in the spatial domain have improved with advances in remote sensing imagery data availability, resolution, and analysis. Sub-pixel correlation and digital elevation model (DEM) differencing from sub-meter, earthquake-spanning satellite optical imagery has enhanced surface rupture mapping and deformation measurements. However, knowledge of measurement accuracy and uncertainty is limited. To address this, we construct orthophotos and digital elevation models (DEMs) from repeat high resolution (∼0.5 m) satellite optical imagery along two sections of the Garlock fault, California with clear fault geomorphology and differing land cover. We deform later sets of DEMs and images with synthetic earthquakes containing both diffuse and discrete horizontal and vertical displacements. Sub-pixel image correlation and DEM differencing demonstrate how vegetation degrades recovered displacement accuracy. In barren land cover, horizontal displacements are detectable to an expected ∼1/10th-pixel size. With shrubs, trees, and grass, detectable displacements increase to >1/2-pixel size, and filtering results by correlation score and using elevation values as input rather than image values improves accuracy. Vertical displacement detection thresholds remain lower in vegetation, at >1-pixel size. Higher slope angles degrade displacement recovery, worsened by vegetation. Diminishing seasonal separation improves accuracy over vegetated regions, though not to the level achieved in barren environments. These results will inform research and operational efforts on the utility of high resolution satellite optical imagery for detecting deformation in varied land cover. Furthermore, they reveal where alternative measurements, such as from LiDAR or radar interferometry, are required to mitigate the effects of vegetation and capture fine-scale crustal deformation.

California↗
USGS2024

Effects of stochastically-simulated near-fault ground motions on soil liquefaction

The scarcity of historically recorded near-fault ground motions poses a challenge to systematically understanding the influence of near-fault effects on various types of seismic demands for engineering purposes. In particular, the current state of knowledge of the influence of ground-shaking intensity on soil liquefaction and its consequences does not specifically account for the effects of near-fault ground motion characteristics. In this study, the influence of near-fault ground motions on liquefaction triggering and lateral spreading are investigated using non-linear modeling of a hypothetical liquefiable soil column in the finite-element computational platform OpenSees subjected to simulated ground motion time series that represent strong earthquake shaking in the near field. The simulated ground motion time series and resulting datasets are based on a parametric stochastic model and are developed for a range of source and path parameters to represent a realistic variability of ground motion characteristics. Dependencies between ground motion intensity measures (IMs) and liquefaction demand parameters are investigated for near-fault pulse and nonpulse-like ground motion sets. Evolutionary IMs, such as cumulative absolute velocity (CAV) and the time-varying magnitude-adjusted peak ground acceleration (PGAM), are considered in developing liquefaction triggering probability density functions. Post-liquefaction triggering responses such as lateral spreading displacements are examined in relation to PGA M and CAV. The ground motion simulations are validated by comparing their liquefaction-capacity PGA M fragilities and post-triggering CAV vulnerability relationships to historical records from the 1994 Northridge earthquake in California, USA. Finally, a path forward for future studies that includes finding systematic differences in the IM-liquefaction demand relationships between near-fault and far-field stochastic ground motion sets is outlined.

Conference Paper↗
USGS2026

Ductile and brittle Rio Grande Rift deformation in Oligocene granite records a two-stage rift history in southern Colorado

The timing and nature of early deformation in the Rio Grande Rift remains poorly constrained. We present evidence for the earliest structural signature of rift extension in the Sangre de Cristo Range, southern Colorado, based on new geologic mapping, structural analysis, rock magnetic data, and thermochronology. These analyses focus on the ~30.0 Ma granite of Chokecherry Canyon, which hosts discrete low-angle mylonitic shear zones and a distributed, gently SW-dipping protomylonitic fabric. Incremental stretching axes, stretching lineations, and Kmax magnetic lineations plunge gently WSW. Quartz microstructures and crystallographic orientations indicate dominantly coaxial strain in the protomylonite and general shear in the discrete shear zones. Quartz c-axis opening-angle thermometry suggests deformation at ~420–540°C. Thermal modeling of ⁴⁰Ar/³⁹Ar K-feldspar data indicates rapid post magmatic cooling below the brittle–plastic transition, supporting shear-zone formation immediately after emplacement. Slow cooling from ~20–13 Ma was followed by renewed rapid cooling at ~13 Ma, interpreted as the onset of extensional exhumation along the Sangre de Cristo Fault System. These results show that extension in the northern Rio Grande Rift was active by ~30 Ma, earlier than previously recognized. We propose a two-stage model for northern Rio Grande Rift evolution: Stage I (30–23 Ma) records ENE–WSW extension localized in low-angle mylonitic shear zones associated with mid-crustal intrusions; Stage II (≤18 Ma) reflects brittle high-angle normal faulting, focused exhumation, and rift narrowing. Stage I magmatism and deformation along the western range front likely established crustal weaknesses that guided later fault development.

Colorado↗
USGS2025

Projecting stream water quality using Weighted Regression on Time, Discharge, and Season (WRTDS): An example with drought conditions in the Delaware River Basin

Future water availability depends on understanding the responses of constituent concentrations to hydrologic change. Projecting future water quality remains a methodological challenge, particularly when using discrete observations with limited temporal resolution. This study introduces Weighted Regression on Time, Discharge, and Season for Projection (WRTDS-P), a novel, computationally efficient method that enables the projection of daily stream water quality under varying hydrologic conditions using commonly available discrete monitoring data. WRTDS-P model performance was validated using 39 sites in the Delaware River Basin (DRB) and four key constituents: specific conductance (SC), nitrate (NO 3 − ), magnesium (Mg 2+ ) and calcium (Ca 2+ ). Projections were tested against holdout data from the final 1 to 5 years of each time series, demonstrating robust predictive capability, with median Nash-Sutcliffe efficiencies of 0.67 for SC, 0.56 for NO 3 − , 0.65 for Ca 2+ , and 0.79 for Mg 2+ . Model uncertainty was correlated with indicators of hydrologic or geochemical mass-sinks, such as groundwater storage and adsorption in wetland soils. Drought scenario analyses for SC used ranges of reduced discharge including flows from the 1965 drought of record. Scenarios predicted widespread increases of SC, especially in southern DRB streams where baseline SC levels are already elevated. Fractional increases of SC were more uniformly distributed, indicating potential risk to sensitive ecosystems. Notably, drought-induced SC increases were positively correlated with interannual SC trends, indicating that hydrologic extremes could exacerbate ongoing salinization. This work provides a transferable and interpretable framework for projecting future water quality and assessing hydrologic risk to water resources and aquatic ecosystems.

Delaware, New Jersey, New York, Pennsylvania↗
USGS2025

Dynamic baseflow storage estimates and the role of topography, geology and evapotranspiration on streamflow recession characteristics in the Neversink Reservoir Watershed, New York

Estimates of dynamic groundwater volumes supplying baseflow to streams are important for water availability projections during extended periods of drought. The primary goals of this study were to provide dynamic storage volume estimates, inferred from streamflow recession analysis, for baseflow regimes within seven gaged catchments within the Neversink Reservoir Watershed (NRW), a critical municipal water source for New York City. Additionally, geomorphological properties, surficial geology and hydro-meteorological processes were quantified and described in relation to time and spatially variable recession behaviour and storage estimates across the NRW. To explore these relationships, we (1) evaluated seasonal trends in streamflow recession behaviour in relation to modelled potential evapotranspiration (PET) and catchment runoff rates, (2) derived empirical streamflow models for cool-season runoff using both linear and nonlinear reservoir assumptions for baseflow and (3) calculated metrics related to the geology and geomorphology of each catchment and compared these metrics to area normalised baseflow dynamic storage estimates. Results show that baseflow recession behaves as a nonlinear reservoir, and applying linear groundwater reservoir assumptions may underestimate the total dynamic storage volumes compared to what would be predicted for a nonlinear reservoir. Increases in PET caused decreases in storage conditions that resulted in increased recession rates and nonlinearity in streamflow recession during the growing season. Additionally, we found that while no single physical catchment characteristic solely predicted catchment storage dynamics, sediment volume and stream gradients were stronger predictors of normalised storage volumes than catchment surface area or surface topography alone. Within the NRW, catchments with the highest sediment volume exhibited the lowest recession rates and higher dynamic storage volumes, while the smallest catchment, mostly devoid of sediment, had the fastest recession rate and lowest dynamic storage volume.

New York↗
USGS2025

Ultramafic float rocks at Jezero crater (Mars): Excavation of lower crustal rocks or mantle peridotites by impact cratering?

Based on observation and data from meteorites and in situ scientific missions, experiments as well as models, the Martian mantle is assumed to share some compositional and mineralogical affinity with the terrestrial mantle. However, there might be subtle differences like the Martian mantle being more ferroan. Yet, we do not have any direct analysis of a Martian mantle rock to confirm this assumption. NASA’s Perseverance rover found olivine-rich boulder-sized float rocks on the upper Jezero fan (Mars). These boulders have an ultramafic composition and their mineralogy is dominantly composed of Fo 73±3 olivine with high-Mg orthopyroxene, Cr-rich Ti-Fe oxides and minor plagioclase and high-Ca pyroxene. Microtextural and petrological analysis reveals that these minerals crystallized at equilibrium. In addition, these boulders are different from all the bedrocks analyzed by Perseverance along its traverse which are crustal igneous rocks and sediments. Comparing our data to Martian meteorites and available Mars bulk silicate models (BSM), we discuss that these boulders could represent primitive melts and/or lower crustal material, and we specifically hypothesize that they could be mantle peridotites. We propose that these putative mantle rocks could have been excavated by the succession of impacts from the shallow mantle or lower crust in the Isidis region where Jezero crater is located. These olivine-rich boulders could thereby constitute the first direct analysis of a Martian mantle rock.

Earth and Planetary Science Letters↗
USGS2026

Monazite and xenotime U-Pb geochronology and thermometry of the Blue Ridge and Inner Piedmont of North Carolina: Implications for the thermal-metamorphic evolution of the southern Appalachian metamorphic “core”

The southern Appalachian orogen preserves a complex distribution of metamorphism and deformation varying in timing, magnitude, and spatial extent. These complexities give rise to disparate interpretations for southern Appalachian tectonic evolution, which complicates the testing and interpretation of tectonic models in this system. New monazite (Mnz) and xenotime (Xtm) laser ablation split stream (LASS) analyses alongside Mnz-Xtm thermometry in the orogenic core in the eastern Blue Ridge (EBR), western Inner Piedmont (WIP), and Cat Square terranes (CST) of North Carolina yield new constraints that define distinct pro- and retrograde metamorphic events. The EBR preserves two prograde thermal events: the Taconic (∼470-440 Ma, >660°C) and Neoacadian (∼380-340 Ma, 600–700°C), separated by a period of cooling (exhumation?) and followed by garnet breakdown from 339 to 329 Ma. Evidence of pervasive Neoacadian ductile deformation in the EBR is largely limited to the Brevard fault zone (BFZ), indicating that a major rheological gradient existed across the BFZ during the Neoacadian and early Alleghanian. Southeast of the BFZ, in the WIP and CST, monazite data define a protracted Neoacadian evolution from early mineral growth at ∼405 Ma at ∼450–600°C to >700°C at ∼360 Ma, followed by early Alleghanian retrograde metamorphism and deformation (<345 Ma, 350–500°C). These constraints, together with previously reported thermobarometric data, define a P-T-t evolution for the WIP and CST consistent with Neoacadian crustal flow, while the coeval presence of a thermal-rheological boundary along the BFZ further supports a model of Neoacadian crustal “escape” flow within the orogen.

North Carolina↗
USGS2024

Post-wildfire debris flows

Post-wildfire debris flows pose severe hazards to communities and infrastructure near and within recently burned mountainous terrain. Intense heat of wildfires changes the runoff characteristics of a watershed by combusting the vegetative canopy, litter, and duff, introducing ash into the soil and creating water repellant soils. Following wildfire, rainfall on bare ground is less able to infiltrate into the fire-altered soils and overland flow is less impeded by vegetation. Rainfall runoff in recently burned areas can erode hillslopes owing to the removal of soil binding organic matter near the soil surface by fire. In channels, loose, dry-ravel deposits composed of sand and gravel are readily entrained by concentrated runoff in channels. Entrainment of soil on hillslopes and in channels bulks up the sediment concentration of the rainfall runoff to generate debris flows capable of transporting boulders and large woody debris. Post-wildfire debris flows can be triggered by rainfall conditions that would typically produce little runoff during unburned conditions. The primary rainfall trigger for post-wildfire debris flows is high intensity rainfall during short duration convective rainstorms or periods of high rainfall intensity embedded within a long-duration frontal storm. Numerous observations of debris flows triggered by storms lasting less than an hour following periods of little to no rainfall indicate that antecedent rainfall is not a requirement for initiation of post-wildfire debris flows. Post-wildfire debris-flow hazard assessment entails estimating probability and magnitude of debris flows in the burned area, estimating debris-flow runout and intensity, and defining rainfall intensity-duration thresholds for debris-flow initiation. In the United States, probability and magnitude is estimated using empirically derived models largely based on data collected in southern California. The models provide maps to identify watersheds and drainage paths where post-wildfire hazards are most pronounced. Rainfall intensity-duration thresholds can be incorporated into flood hazard forecasting tools. Currently, work is underway to identify how to best implement debris-flow runout models in burned areas with efficiency and accuracy. Post-wildfire debris flows have been a long-recognized process in the Transverse Ranges of southern California; however, climate change is driving more frequent wildfires to burn more mountainous terrain throughout the western United States and worldwide. As a result, post-wildfire debris flows are becoming a more common threat in areas where they were once infrequent. As the threat of post-wildfire debris flow expands into new areas, evaluating the hazard becomes challenging because the degree to which wildfire increases debris-flow susceptibility varies from region to region. This chapter summarizes the knowledge to date for evaluating post-wildfire debris-flow susceptibility and hazard assessment. We summarize the characteristics of wildfire burn severity, topography, underlying soil and geology, and rainfall conditions that contribute to making a watershed most likely to produce post-wildfire debris flows. Methods for hazard assessment in the United States and other countries are summarized. We highlight knowledge gaps for how post-wildfire debris-flow susceptibility varies throughout the western United States and worldwide and identify research needs to improve hazard assessment methods in different geographies.

Book chapter↗
USGS2026

Rapid seismic and infrasound assessment of large landslides: A case study from Denali National Park and Preserve (Alaska)

Large, rapid landslides are a global hazard that can occur in remote, mountainous areas. Eyewitness reports of landslides and satellite imagery can often be limited or delayed, particularly during inclement weather. However, landslide-generated seismic and infrasound (low-frequency atmospheric sound) waves can be remotely detected in near real-time. This information can significantly expedite characterization and possible landslide response activities. Here, we highlight these capabilities using a > 4 million m³ ice–rock avalanche in Denali National Park and Preserve (Alaska). This event was detected via a landslide-specific seismic location and volume estimation algorithm deployed in Alaska, and — notably — by standard earthquake monitoring systems. Following rapid detection of this event, we combined its seismic and infrasound dataset with optical, synthetic aperture radar, and oblique aerial imagery, multitemporal digital elevation models, and a numerical flow model to reconstruct its failure timeline and dynamics. We apply array processing to infrasound signals traveling > 250 km and find that two precursory events occurred minutes prior to the main failure. We use long-period seismic signals to infer the force exerted by the landslide on the Earth and constrain the rheological parameters of our numerical flow simulation with this result and deposit morphology. The main failure produced a steeply-dipping impulsive initial downward force and reached speeds exceeding 60 m/s. This impulsive force generated relatively strong seismic body waves, which contributed to the earthquake system detection. This large, remote Alaska landslide underscores the key value of seismic and infrasound analysis for rapid landslide assessment and motivates efforts to further operationalize these approaches.

Alaska↗
USGS2026

High-resolution magnetic survey using an unoccupied aerial vehicle to constrain buried lava flow geometry, volume, and eruptive history of Little Cones, Crater Flat, Nevada

Magnetic surveys are an important tool used to augment geologic mapping in distributed volcanic fields. Using magnetic anomalies, it is possible to model the geometry of shallowly buried volcanic features, such as conduits, sills, and lava flows. This subsurface mapping is important for understanding eruption dynamics and emplacement of lava flows, and it sometimes reveals buried volcanoes no longer visible at the surface. These data are critical to better interpret the numbers, styles, and magnitudes of eruptions in distributed volcanic fields and their associated volcanic hazards. New advances in unoccupied aerial vehicles (UAVs) offer an attractive middle range of resolution and aerial coverage between ground-based magnetic surveys and aeromagnetic surveys. Here, we present the results of a UAV fluxgate magnetic survey of the Little Cones, Nevada, scoria cones, which have been the target of previous ground and aeromagnetic surveys. The magnetic anomalies at Little Cones are of interest because the surrounding alluvium conceals lava flows that erupted from Little Cones, making it very difficult to understand the volume and morphology of lava flows from geologic mapping alone. Nonlinear inversion of UAV-collected magnetic data were used to model the thickness and morphology of buried Little Cones’ lava flows with higher precision than achieved previously. The sequence of events and calculated flow characteristics are then interpreted. The total volume of Little Cones, including concealed lava flows, is approximately 0.016 cubic kilometer, and the initial sheet flow erupted in less than 24 hours. The findings presented herein demonstrate that UAV-based magnetic surveys are a reliable method of data collection and an efficient alternative to other survey methods, facilitating development of a three-dimensional perspective of distributed volcanic fields.

Nevada↗
USGS2024

Characterization of change in tree cover state and condition over the conterminous United States

Variability in the effects of disturbances and extreme climate events can lead to changes in tree cover over time, including partial or complete loss, with diverse ecological consequences. It is therefore critical to identify in space and time the change processes that lead to tree cover change. Studies of change are often hampered by the lack of data capable of consistently detecting different types of change. Using the Landsat satellite record to create a long time-series of land cover and land cover change, the U.S. Geological Survey Land Change Monitoring Assessment and Projection (LCMAP) project has made an annual time series of land cover across the conterminous United States for the period 1985 to 2018. Multiple LCMAP products analyzed together with map validation reference plots provide a robust basis for understanding tree cover change. In LCMAP (Collection 1.2), annual change detection is based on harmonic model breaks calculated at each Landsat pixel from the Continuous Change Detection and Classification (CCDC) algorithm. The results showed that the majority of CCDC harmonic model breaks (signifying change) indicated partial tree cover loss (associated with management practices such as tree cover thinning) as compared to complete tree cover loss (associated with practices like clearcut harvest or fire disturbance). Substantially fewer occurrences of complete tree cover loss were associated with change in land cover state. The area of annual tree cover change increased after the late 1990s and stayed high for the rest of the study period. The reference data showed that tree harvest dominated across the conterminous United States. The majority of tree cover change occurred in evergreen forests. Large estimates of disturbance-related tree cover change indicated that tree cover loss may have previously been underreported due to omission of partial tree cover loss in prior studies. This has considerable implications for forest carbon accounting along with tracking ecosystem goods and services.

Forests↗
USGS2025

Implications of physics-based M9 ground motions on liquefaction-induced damage in the Cascadia Subduction Zone: Looking forward and backward

Given the likelihood of future M9 Cascadia Subduction Zone (CSZ) earthquakes, various estimates of the resulting, regional ground motions have been made, including a suite of 30 physics-based simulations that reflect key modeling uncertainties. However, because the last CSZ interface rupture occurred in 1700 CE, the shaking expected in such an event is especially uncertain, as are the impacts to the built and living environments. Like other coseismic impacts, soil liquefaction poses a significant threat and must be considered by any scenario study used to inform planning and response, or to focus mitigation resources. Liquefaction is also notable for its potential to “ground truth” ground-motion estimates, given that its presence or absence in the geologic record can provide constraint on the intensities of shaking in past events. It is thus an important phenomenon looking both forward and backward. Accordingly, using recent physics-based simulations, this study (1) predicts liquefaction in M9 CSZ ruptures at 400 locations in Oregon, Washington, and British Columbia (BC) using an array of cone-penetration-test based models and (2) uses paleoliquefaction evidence at ten sites spanning from Southern Oregon to Vancouver, BC to constrain possible ground-motion intensities experienced in the 1700 CE earthquake. The forward predictions indicate that liquefaction in M9 events could be pervasive in the region and affect numerous population hubs, with the potential for damage across hundreds of square kilometers. The backward analyses suggest that 1700 CE ground-motion intensities may have been less than expected from M9 simulations in some northern portions of the CSZ (e.g. Seattle), given the paucity of 1700 CE liquefaction evidence in these areas. Ultimately, further discovery and analysis of CSZ paleoliquefaction, or lack thereof, will confirm or modify this possibility and the conclusions drawn herein.

Oregon, Washington↗
USGS2025

The 3D Elevation Program—Supporting the economy of Massachusetts

Introduction Massachusetts extends from the mountains of the Appalachian system in the west of the State to the sandy beaches and rocky shorelines of the Atlantic coast in the east. Inland topographic data support a wide range of important activities, including geologic mapping, transportation planning, forest and wildlife management, quantifying ecological services, water supply protection, commonwealth-wide infrastructure planning, local site planning, and flood-plain management. Nearshore bathymetry can be used to support coastal portions of the Commonwealth by addressing the combined threats of ocean warming, strong storm surge, and rising sea levels. The maintenance and (or) expansion of Massachusetts ports (for instance, Boston, New Bedford) and Cape Cod sediment management depends upon the accurate mapping of bathymetry and the frequent influx of sediment and redeposition. Critical applications that address the broad range of requirements depend on light detection and ranging (lidar) data that provide a highly detailed three-dimensional (3D) model of the Earth’s surface and aboveground features. The 3D Elevation Program (3DEP) is managed by the U.S. Geological Survey (USGS) in partnership with Federal, State, Tribal, U.S. territorial, and local agencies to acquire consistent lidar coverage at quality level 2 or better to meet the many needs of the Nation and Massachusetts. The status of available and in-progress 3DEP baseline lidar data in Massachusetts is shown in figure 1. 3DEP baseline lidar data include quality level 2 or better, 1-meter or better digital elevation models, and lidar point clouds, and must meet the Lidar Base Specification version 1.2 ( https://www.usgs.gov/3dep/lidarspec ) or newer requirements. The National Enhanced Elevation Assessment identified user requirements and conservatively estimated that availability of lidar data would result in at least $1.23 million in new benefits annually to Massachusetts. The top 10 Massachusetts business uses for 3D elevation data, which are based on the estimated annual conservative benefits of 3DEP, are shown in table 2.

Massachusetts↗
USGS2026

Simulation-based scenario ShakeMaps for large magnitude (MW6.5+) crustal earthquakes on the Seattle, Tacoma, and southern Whidbey Island faults, Washington, USA

Scenario ground‐motion maps based on empirical ground‐motion models (GMMs) provide a rapid and generally reliable means of estimating the amplitude and distribution of earthquake shaking. However, because GMMs are designed for broad applicability, they often rely on simplified representations of Earth structure, which can limit their accuracy in regions with complex source, path, and site effects. This can substantially impact the accuracy of predicted shaking in areas like western Washington State, where deep, interconnected basin structure exerts a strong influence on seismic‐wave propagation. In this study, we present a new suite of simulation‐based scenario ShakeMaps that characterize ground shaking from large‐magnitude ( ⁠ M W 6.5–7.5) earthquakes on the Seattle, Tacoma, and southern Whidbey Island faults. These maps are developed using results from recent 3D wave propagation simulations ( Stone et al. , 2022 , 2023 , 2025 ) that incorporate realistic rupture geometries, variable slip distributions, and a regional 3D seismic velocity model with shallow soils. Broadband ground motions are estimated by combining the low‐frequency (<1 Hz) deterministic seismograms from these studies with high‐frequency (1–10 Hz) stochastic seismograms. Simulated ground motions are corrected to account for the enforced minimum shear‐wave velocity and nonlinear site response. The resulting ShakeMaps represent median ground‐shaking estimates derived from multiple rupture scenarios with varying slip distributions and hypocenter locations for each fault. To extend ShakeMap coverage beyond the simulation domain (i.e., into eastern Washington, northern Oregon, and southwestern British Columbia), we scale GMM‐based ground‐motion estimates using amplification patterns observed in the simulations. These new ShakeMaps reveal the substantial influence of deep basin structure on shaking intensity, underscoring the importance of considering crustal structure complexity in regional hazard assessments for the Pacific Northwest.

Washington↗
USGS2026

Characterizing precariously balanced rocks (PBRs) in the eastern United States for estimating maximum past earthquake ground motions

Fragile geologic features (FGFs) can be used to estimate an upper bound to the strength of earthquake ground shaking in the past thousands to tens of thousands of years. FGFs could be especially useful in the eastern United States, where few active faults are available for paleoseismic studies and where low seismic attenuation means that an earthquake could damage FGFs at large distances. This article describes a pilot study of one form of FGF, precariously balanced rocks (PBRs), in the Blue Ridge of Virginia, ∼160 km from the central Virginia seismic zone (CVSZ) and ∼100 km from Virginia’s Giles County seismic zone (GCSZ). The main goal is to develop and refine methods for studying PBRs in the eastern United States. Photogrammetric 3D models of the rocks were used to locate their centers of mass and delineate their basal contacts, with the latter corroborated by wrapping a wire around the bases to verify the circumferences. Relations between hill geometry and amplification were used to estimate topographic amplification. Seismometers on the PBRs provided the frequencies, amplitudes, and directions of rocking after giving the rocks gentle pushes. The levels of ground shaking likely to topple the rocks were compared with the median hazard curves from the U.S. Geological Survey’s National Seismic Hazard Model and were used to map magnitudes of earthquakes in the nearby region likely to topple the rocks. Results indicate that an M w ≥ 7.5 earthquake is needed in the CVSZ or near the 2020 Sparta, North Carolina, earthquake to topple the two most fragile PBRs and that an earthquake of ~ M w 7–7.5 in the GCSZ is needed. The fragilities and estimated ages of the PBRs (20–40 ka) are thus consistent with the absence of such events in the historical record in the nearby region.

eastern United States↗
USGS2024

Precariously balanced rocks in northern New York and Vermont, U.S.A.: Ground-motion constraints and implications for fault sources

Precariously balanced rocks (PBRs) and other fragile geologic features have the potential to constrain the maximum intensity of earthquake ground shaking over millennia. Such constraints may be particularly useful in the eastern United States (U.S.), where few earthquake‐source faults are reliably identified, and moderate earthquakes can be felt at great distances due to low seismic attenuation. We describe five PBRs in northern New York and Vermont—a region of elevated seismic hazard associated with historical seismicity. These boulders appear to be among the most fragile PBRs in the region, based on reports from hobbyists. The PBRs are glacial erratics, best evidenced by glacial striations on bedrock pedestals. The pedestals themselves are locally high knobs, often situated on regionally high topography; this setting limits soil development and indicates that any outwash deposits were likely ephemeral. As a result, PBR ages can be reliably established by the retreat of the last continental ice sheet, ∼15–13 ka. To quantify the fragility of the PBRs, we surveyed them with ground‐based light detection and ranging and calculated geometric parameters from the point clouds, field observations, and seismic responses. Preliminary validation of the 2023 time‐independent U.S. National Seismic Hazard Model (NSHM) shows that the existence of PBRs is generally consistent with the median site‐specific hazard curves. Only the Blue Ridge Road site suggests a modest reduction in hazard. To visualize the ensemble of data, we mapped the minimum permissible distance to potential source faults around each PBR site as a function of source magnitude by using the ground‐motion models from the 2023 NSHM. Viewed in this manner, our data are consistent with potential M ∼6.5 earthquake‐source faults in many parts of the Lake Champlain Valley and northern Adirondack Mountains. Our work illustrates a potential pathway for better constraining earthquake‐source faults in regions of cryptic faults.

Bulletin of the Seismological Society of America↗
USGS2026

The petrogenesis of Þingmúli volcano, East Fjords, Iceland

In this work we revisit Þingmúli volcano (Þ = Th), a classic locality known as an example of a complete tholeiitic differentiation. Þingmúli is a ~ 9.5 Ma extinct central volcano located in the East Fjords of Iceland, in which the whole compositional spectrum from basalt to rhyolites have erupted. These volcanic products have been previously considered as petrogenetically related by an ideal fractionation trend, regardless any temporal relationship or volumetric considerations. Here we report new whole-rock geochemistry, mineral chemistry, isotope analyses, estimation of residence times of the different eruptive deposits, and an update of the original petrogenetic model. Our results highlight that an enriched source, likely spinel lherzolites, generated transitional-alkaline basaltic melts after 15–20% of partial melting at depths of 40–45 km. Many of these basaltic melts erupted at various stages of the volcano's history, while others remained longer in the volcanic plumbing system. These evolved by fractional crystallisation into basaltic andesite magmas with a residence time of ~5 years based on the crystal size distribution of the plagioclase population. Isotopic differences between the basalts/basaltic andesites ( 87 Sr/ 86 Sr ~ 0.7034; 143 Nd/ 144 Nd ~ 0.51315) and the erupted rhyolites ( 87 Sr/ 86 Sr ~ 0.7037; 143 Nd/ 144 Nd ~ 0.51304) indicate that the latter are not petrogenetically related to the former. Therefore, instead of a fractional crystallisation mechanism to generate the rhyolites, we propose the partial melting of ignimbrite layers located beneath the volcano. The broad range of trace element concentrations in andesites and dacites and their different isotopic values compared to the basalts strongly suggest that these magmas have been generated by magma mixing between basaltic and rhyolitic melts, similar to modern day Icelandic volcanoes such as Hekla. These results highlight the need to revisit previously studied Icelandic classic localities and reassess their traditionally proposed petrogenetic models.

Þingmúli volcano, East Fjords↗
USGS2025

Summary of selenium in the lower Gunnison River Basin, Colorado—Information and data gaps

The Cretaceous Mancos Shale is a geologic source of selenium in the lower Gunnison River Basin. Natural weathering processes and human activity mobilize selenium from the Mancos Shale and derived materials, and surface water, groundwater, and sediment all affect the transport of selenium from source areas to receiving streams and biota. Selenium accumulates through the aquatic food chain, and its toxic effects can result in invertebrate mortality and mortality, decreased reproduction, and deformities to fish and birds. The Bureau of Reclamation, in cooperation with the State of Colorado and Gunnison River Basin water users, is implementing a Selenium Management Program to reduce selenium concentrations in the lower Gunnison River Basin of Colorado. Goals of the Selenium Management Program are to (1) achieve compliance with the State of Colorado chronic aquatic-life standard for dissolved selenium (4.6 micrograms per liter) in the Gunnison River near Grand Junction, Colorado; (2) sufficiently improve water-quality conditions to assist in the recovery of endangered species in the Gunnison and Colorado Rivers by reducing selenium concentrations; and (3) support continued water uses in the basin. Many previous studies have contributed to the understanding of selenium in the environment; however, monitoring and research data gaps exist in the lower Gunnison River Basin. The purpose of this report is to summarize information regarding selenium in the lower Gunnison River Basin and describe strategies for scientific research and monitoring to potentially improve understanding of selenium sources; processes affecting the mobilization, transport, and fate of selenium; and the effects of selenium-mitigation projects in the lower Gunnison River Basin. Monitoring and research data gaps discussed in this report include geologic mapping and geochemical source characterization, long-term and ongoing monitoring of the surface-water and groundwater networks, developing and refining statistical models, characterizing selenium on suspended sediment, modeling selenium in the food web, evaluating best management practices, and more.

Colorado↗