Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “National Hazards Observer”

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.

At least 91 records · Page 5Linked to original sources
USGS2024

Learning from a high-severity fire event—Conditions following the 2018 Carr Fire at Whiskeytown National Recreation Area

The 2018 Carr Fire burned more than 90 percent of Whiskeytown National Recreation Area, with much of the park burning at high severity. California yellow pine and mixed conifer forests are not well adapted to large, high-severity fires, and forest recovery after these events may be problematic. Large, high-severity fire patches pose difficulties for recruitment with interiors that are long distances from potential seed trees and may develop fuel structures that can promote further high-severity fire. This report details patterns of forest structure derived from field plots measured 2–3 years after the Carr Fire, providing a characterization of immediate fire effects. We coupled these observations with remotely sensed information, including data collected from unoccupied aircraft system surveys. The remotely sensed data were used to depict erosion after the Carr Fire as well as to create a high-resolution land cover classification map, a debris flow risk map and hazard assessment, and a post-fire canopy vegetation loss map. Results indicated high levels of tree mortality after the Carr Fire, including high-value old growth forest stands, supporting remotely sensed assessments of fire severity. The high-resolution tree mortality model also aligned well with other remotely sensed estimates of immediate burn severity. Results of the land cover classification illustrated the high percentage of dead vegetation remaining in the understory and canopy 8 months post-fire. Changes in vegetation height identified areas with canopy vegetation loss from 1- to 8-months post-fire. Pairing the post-fire debris accumulation with debris flow probabilities may identify high-risk debris flow areas. The results of this study will help inform future decisions concerning wildland fire and vegetation management strategies at Whiskeytown National Recreation Area and are broadly relevant for management in the aftermath of large, high-severity fires in mixed, dry coniferous forests in the western United States.

California↗
USGS2013

Final report and archive of the swath bathymetry and ancillary data collected in the Puerto Rico Trench region in 2002 and 2003

In 2002 and 2003, the U.S. Geological Survey (USGS), in cooperation with the National Oceanic and Atmospheric Administration (NOAA), conducted three exploration cruises that mapped for the first time the morphology of the entire tectonic plate boundary stretching from the Dominican Republic in the west to the Lesser Antilles in the east, a distance of approximately 700 kilometers (430 miles). Observations from these three exploration cruises, coupled with computer modeling and published Global Positioning System (GPS) results and earthquake focal mechanisms, have provided new information that is changing the evaluation of the seismic and tsunami hazard from this plate boundary. The observations collected during these cruises also contributed to the basic understanding of the mechanisms that govern plate tectonics, in this case, the creation of the island of Puerto Rico and the deep trench north of it. Results of the sea floor mapping have been an important component of the study of tsunami and earthquake hazards to the northeastern Caribbean and the U.S. Atlantic coast off the United States.

Puerto Rico↗
USGS2017

A fault‐based model for crustal deformation in the western United States based on a combined inversion of GPS and geologic inputs

We develop a crustal deformation model to determine fault‐slip rates for the western United States (WUS) using the Zeng and Shen (2014) method that is based on a combined inversion of Global Positioning System (GPS) velocities and geological slip‐rate constraints. The model consists of six blocks with boundaries aligned along major faults in California and the Cascadia subduction zone, which are represented as buried dislocations in the Earth. Faults distributed within blocks have their geometrical structure and locking depths specified by the Uniform California Earthquake Rupture Forecast, version 3 (UCERF3) and the 2008 U.S. Geological Survey National Seismic Hazard Map Project model. Faults slip beneath a predefined locking depth, except for a few segments where shallow creep is allowed. The slip rates are estimated using a least‐squares inversion. The model resolution analysis shows that the resulting model is influenced heavily by geologic input, which fits the UCERF3 geologic bounds on California B faults and ±one‐half of the geologic slip rates for most other WUS faults. The modeled slip rates for the WUS faults are consistent with the observed GPS velocity field. Our fit to these velocities is measured in terms of a normalized chi‐square, which is 6.5. This updated model fits the data better than most other geodetic‐based inversion models. Major discrepancies between well‐resolved GPS inversion rates and geologic‐consensus rates occur along some of the northern California A faults, the Mojave to San Bernardino segments of the San Andreas fault, the western Garlock fault, the southern segment of the Wasatch fault, and other faults. Off‐fault strain‐rate distributions are consistent with regional tectonics, with a total off‐fault moment rate of 7.2 × 10 18 "> 7.2 × 10 18 and 8.5 × 10 18    N · m / year "> 8.5 × 10 18 N ⋅ m / year for California and the WUS outside California, respectively.

Bulletin of the Seismological Society of America↗
USGS2019

Toward a national coastal hazard forecast of total water levels

Storm surge and large waves combine to erode beaches, cause marsh and coral decay, and inundate low-elevation areas, resulting in hazards to coastal communities and loss of natural resources. The USGS, in collaboration with NOAA, is developing a real-time system to provide ∼ 6-day forecasts of total water levels (TWLs) combining tides, storm surge, and wave runup. TWL is compared with dune elevations along sandy coastlines to obtain regional flood and coastal hazard forecasts. The predictions are available online (https://coastal.er.usgs.gov/hurricanes/research/twlviewer/) and provide intuitive representations of regional coastal risk. Current efforts include 1) assessing skill by comparing forecasts with observations along sandy coastlines and 2) developing methods more appropriate for non-sandy coastlines (e.g., rocky, reef, marsh, cliff, and built shorelines). The vision of this work is to provide a tool for short- to medium-term coastal hazard forecasts to inform planners and emergency responders along nearly the entire US coastline.

Conference Paper↗
USGS2010

Topographic Science

The mission of U.S. Geological Survey (USGS) Earth Resources Observation and Science (EROS) Center Topographic Science is to establish partnerships and conduct research and applications that facilitate the development and use of integrated national and global topographic datasets. Topographic Science includes a wide range of research and applications that result in improved seamless topographic datasets, advanced elevation technology, data integration and terrain visualization, new and improved elevation derivatives, and development of Web-based tools. In cooperation with our partners, Topographic Science is developing integrated-science applications for mapping, national natural resource initiatives, hazards, and global change science. http://topotools.cr.usgs.gov/.

Fact Sheet↗
USGS2020

A 36-year record of rock avalanches in the Saint Elias Mountains of Alaska, with implications for future hazards

Glacial retreat and mountain-permafrost degradation resulting from rising global temperatures have the potential to impact the frequency and magnitude of landslides in glaciated environments. Several recent events, including the 2015 Taan Fiord rock avalanche, which triggered a tsunami with one of the highest wave runups ever recorded, have called attention to the hazards posed by landslides in regions like southern Alaska. In the Saint Elias Mountains, the presence of weak sedimentary and metamorphic rocks and active uplift resulting from the collision of the Yakutat and North American tectonic plates create landslide-prone conditions. To differentiate between the typical frequency of landsliding resulting from the geologic and tectonic setting of this region, and landslide processes that may be accelerated due to changes in climate, we used Landsat imagery to create an inventory of rock avalanches in a 3700 km2 area of the Saint Elias Mountains. During the period from 1984-2019, we identified 220 rock avalanches with a mean recurrence interval of 60 days. We compared our landslide inventory with a catalog of M ≥ 4 earthquakes to identify potential coseismic events, but only found three possible earthquake-triggered rock avalanches. We observed a distinct temporal cluster of 41 rock avalanches from 2013 through 2016 that correlated with above average air temperatures (including the three warmest years on record in Alaska, 2014-2016); this cluster was similar to a temporal cluster of recent rock avalanches in nearby Glacier Bay National Park and Preserve. The majority of rock avalanches initiated from bedrock ridges in probable permafrost zones, suggesting that ice loss due to permafrost degradation, as opposed to glacial thinning, could be a dominant factor contributing to rock-slope failures in the high elevation areas of the Saint Elias Mountains. Although earthquake-triggered landslides have episodically occurred in southern Alaska, evidence from our study suggests that area-normalized rates of non-coseismic rock avalanches were greater during the period from 1964 to 2019, and that the frequency of these events will continue to increase as the climate continues to warm. These findings highlight the need for hazard assessments in Alaska that address changes in landslide patterns related to climate change.

Alaska↗
USGS2026

The U.S. Geological Survey 2025 Puerto Rico and U.S. Virgin Islands time-independent earthquake rupture forecast

We present the 2025 U.S. Geological Survey Puerto Rico and U.S. Virgin Islands (PRVI) time‐independent earthquake rupture forecast (ERF), developed for the 2025 update to the National Seismic Hazard Model (NSHM) for PRVI. The updated ERF improves upon a prior model from 2003, including an expanded fault inventory with slip‐rate estimates, updated seismicity catalogs, and refined subduction zone geometries and deformation models. It applies the fault‐system inversion methodology to solve for rates of ruptures on modeled faults, adapted from the 2023 NSHM (NSHM23) for the western United States, including the first application of the inversion to model rates on a U.S. subduction interface. Off‐fault and intraslab seismicity are constrained by observed seismicity and use updated methods developed for NSHM23. Uncertainties in model components are substantial, and the ERF represents epistemic uncertainties through a comprehensive logic tree consisting of 1.7 billion logic‐tree branches combined across all sources.

Puerto Rico, U.S. Virgin Islands↗
USGS2020

Applications of knowledge and predictions of atmospheric rivers

This chapter reviews how AR research is being applied in real-world situations to address issues of flood planning and emergency intervention. It includes water supply management case studies. Examples comprise five distinct sections that show how AR research is being directly applied to the challenges that water managers, dam operators, crisis-management engineers such as USACE, National Weather Service (NWS) personnel, the media, and others face. These topics include how decision-makers on the ground must iteratively alternate between forecasts and their own field observations, especially in unfolding emergency-response conditions, and the trade-offs necessitated between acting on competing priorities such as flood-risk management and water supply management. Ultimately, almost all AR studies have the potential to directly benefit the public’s need for ongoing water supply as well as for accurate weather forecasts and deployable emergency protocols for natural hazards that necessitate municipal, state, and federal government personnel to collaborate.

Book chapter↗
USGS2023

Sediment thickness map of United States Atlantic and Gulf Coastal Plain Strata, and their influence on earthquake ground motions

With the recent successful accounting of basin depth ground-motion adjustments in seismic hazard analyses for select areas of the western United States, we move toward implementing similar adjustments in the Atlantic and Gulf Coastal Plains by constructing a sediment thickness model and evaluating multiple relevant site amplification models for central and eastern United States seismic hazard analyses. We digitize and combine existing sediment thickness data sets into a composite surface that delineates the base of Cretaceous sediments under the Atlantic Coastal Plain and the base of Mesozoic sediments under the Gulf Coastal Plain. Amplification models dependent on sediment thickness, site natural period, and source-to-site path length are compared with data sets of observed ground motions to evaluate the ability of the new models to improve ground motion estimates. We find that the amplification models can account for observed trends in sediment-thickness and period-dependent residuals, but some tuning is required. For example, the model of Chapman and Guo requires a reference V S 30 , the time-averaged shear-wave velocity within 30 m of the Earth’s surface, for non-Coastal Plain sites, which we estimate to be between about 1 and 2 km/s. Along with our sediment thickness model, we estimate a velocity profile for application to the Harmon et al. site-natural-period-based model in order to best match the Chapman and Guo period dependence for a broad range of sediment thicknesses. The Next Generation of Attenuation models for the eastern United States Gulf Coast path-based adjustment models can also account for seismic attenuation in the Coastal Plain sediments and reduce the standard deviation of total residuals. If enacted in the U.S. Geological Survey National Seismic Hazard Model, these amplification models will reduce predicted short-period (<1 s) and increase predicted long-period (>1 s) ground motions in the Coastal Plains appreciably.

Earthquake Spectra↗
USGS2018

2018 update to the U.S. Geological Survey national volcanic threat assessment

When erupting, all volcanoes pose a degree of risk to people and infrastructure, however, the risks are not equivalent from one volcano to another because of differences in eruptive style and geographic location. Assessing the relative threats posed by U.S. volcanoes identifies which volcanoes warrant the greatest risk-mitigation efforts by the U.S. Geological Survey and its partners. This update of the volcano threat assessment of Ewert and others (2005) considers new research in order to determine which volcanic systems should be added or removed from the list of potentially active volcanoes, updates the scoring of active volcanoes, and updates the 24-factor hazard and exposure matrix used to create the threat ranking. The threat assessment places volcanoes into five threat categories: very low, low, moderate, high, and very high. Within all five threat categories there are changes in relative rankings of volcanoes, and in a few cases, volcanoes moved between categories owing to changes in our understanding of their hazard, unrest, and exposure factors. Scorings of hazard factors were updated for some volcanoes where new research has identified Holocene eruptive activity or clarified our understanding of Holocene eruptive history and the occurrence of particular hazards such as tephra fall or pyroclastic density currents. The most numerous scoring changes made in the threat matrix since 2005 have been made among the hazard factors, particularly those accounting for observed eruptive activity or unrest. The very low threat category underwent the greatest amount of change, dropping from 32 to 21 volcanoes, owing to better knowledge of the eruptive histories of those volcanoes. The list of 18 very high threat volcanoes determined by Ewert and others (2005) remains the same; 11 of the 18 volcanoes are located in Washington, Oregon, or California, where explosive and often snow- and ice-covered edifices can project hazards long distances to densely populated and highly developed areas. Five of the 18 very high threat volcanoes are in Alaska near important population centers, economic infrastructure, or below busy air traffic corridors. The remaining two very high threat volcanoes are on the Island of Hawaiʻi, where densely populated and highly developed areas now exist on the flanks of highly active volcanoes. The high- and moderate-threat categories are dominated by Alaskan volcanoes. In these categories the generally more active and more explosive volcanoes in Alaska can have a substantial effect on national and international aviation, and large eruptions from any of the moderate- to very-high-threat volcanoes could cause regional or national-scale disasters. This revised threat assessment includes 18 very high threat, 39 high threat, 49 moderate threat, 34 low threat, and 21 very low threat volcanoes. The total of 161 volcanoes is a decrease of 8 from the total reported by Ewert and others (2005).

Scientific Investigations Report↗
USGS2009

Geologic Map of the Carlton Quadrangle, Yamhill County, Oregon

The Carlton, Oregon, 7.5-minute quadrangle is located in northwestern Oregon, about 35 miles (57 km) southwest of Portland. It encompasses the towns of Yamhill and Carlton in the northwestern Willamette Valley and extends into the eastern flank of the Oregon Coast Range. The Carlton quadrangle is one of several dozen quadrangles being mapped by the U.S. Geological Survey (USGS) and the Oregon Department of Geology and Mineral Industries (DOGAMI) to provide a framework for earthquake- hazard assessments in the greater Portland, Oregon, metropolitan area. The focus of USGS mapping is on the structural setting of the northern Willamette Valley and its relation to the Coast Range uplift. Mapping was done in collaboration with soil scientists from the National Resource Conservation Service, and the distribution of geologic units is refined over earlier regional mapping (Schlicker and Deacon, 1967). Geologic mapping was done on 7.5-minute topographic base maps and digitized in ArcGIS to produce ArcGIS geodatabases and PDFs of the map and text. The geologic contacts are based on numerous observations and samples collected in 2002 and 2003, National Resource Conservation Service soils maps, and interpretations of 7.5-minute topography. The map was completed before new, high-resolution laser terrain mapping was flown for parts of the northern Willamette Valley in 2008.

Open-File Report↗
USGS2016

Municipal solid waste landfills harbor distinct microbiomes

Landfills are the final repository for most of the discarded material from human society and its “built environments.” Microorganisms subsequently degrade this discarded material in the landfill, releasing gases (largely CH 4 and CO 2 ) and a complex mixture of soluble chemical compounds in leachate. Characterization of “landfill microbiomes” and their comparison across several landfills should allow the identification of environmental or operational properties that influence the composition of these microbiomes and potentially their biodegradation capabilities. To this end, the composition of landfill microbiomes was characterized as part of an ongoing USGS national survey studying the chemical composition of leachates from 19 non-hazardous landfills across 16 states in the continental U.S. The landfills varied in parameters such as size, waste composition, management strategy, geography, and climate zone. The diversity and composition of bacterial and archaeal populations in leachate samples were characterized by 16S rRNA gene sequence analysis, and compared against a variety of physical and chemical parameters in an attempt to identify their impact on selection. Members of the Epsilonproteobacteria, Gammaproteobacteria, Clostridia, and candidate division OP3 were the most abundant. The distribution of the observed phylogenetic diversity could best be explained by a combination of variables and was correlated most strongly with the concentrations of chloride and barium, rate of evapotranspiration, age of waste, and the number of detected household chemicals. This study illustrates how leachate microbiomes are distinct from those of other natural or built environments, and sheds light on the major selective forces responsible for this microbial diversity.

Frontiers in Microbiology↗
USGS1967

Volcanic hazards at Mount Rainier, Washington

Mount Rainier is a large stratovolcano of andesitic rock in the Cascade Range of western Washington. Although the volcano as it now stands was almost completely formed before the last major glaciation, geologic formations record a variety of events that have occurred at the volcano in postglacial time. Repetition of some of these events today without warning would result in property damage and loss of life on a catastrophic scale. It is appropriate, therefore, to examine the extent, frequency, and apparent origin of these phenomena and to attempt to predict the effects on man of similar events in the future. The present report was prompted by a contrast that we noted during a study of surficial geologic deposits in Mount Rainier National Park, between the present tranquil landscape adjacent to the volcano and the violent events that shaped parts of that same landscape in the recent past. Natural catastrophes that have geologic causes - such as eruptions, landslides, earthquakes, and floods - all too often are disastrous primarily because man has not understood and made allowance for the geologic environment he occupies. Assessment of the potential hazards of a volcanic environment is especially difficult, for prediction of the time and kind of volcanic activity is still an imperfect art, even at active volcanoes whose behavior has been closely observed for many years. Qualified predictions, however, can be used to plan ways in which hazards to life and property can be minimized. The prediction of eruptions is handicapped because volcanism results from conditions far beneath the surface of the earth, where the causative factors cannot be seen and, for the most part, cannot be measured. Consequently, long-range predictions at Mount Rainier can be based only on the past behavior of the volcano, as revealed by study of the deposits that resulted from previous eruptions. Predictions of this sort, of course, cannot be specific as to time and locale of future events, and clearly are valid only if the past behavior is, as we believe, a reliable guide. The purpose of this report is to infer the events recorded by certain postglacial deposits at Mount Rainier and to suggest what bearing similar events in the future might have on land use within and near the park. In addition, table 2 (page 22) gives possible warning signs of an impending eruption. We want to increase man's understanding of a possibly hazardous geologic environment around Mount Rainier volcano, yet we do not wish to imply for certain that the hazards described are either immediate or inevitable. However, we do believe that hazards exist, that some caution is warranted, and that some major hazards can be avoided by judicious planning. Most of the events with which we are concerned are sporadic phenomena that have resulted directly or indirectly from volcanic eruptions. Although no eruptions (other than steam emission) of the volcano in historic time are unequivocally known (Hopson and others, 1962), pyroclastic (air-laid) deposits of pumice and rock debris attest to repeated, widely spaced eruptions during the 10,000 years or so of postglacial time. In addition, the constituents of some debris flows indicate an origin during eruptions of molten rock; other debris flows, because of their large size and constituents, are believed to have been caused by steam explosions. Some debris flows, however, are not related to volcanism at all.

Washington↗
USGS2007

Simulations of potential runout and deposition of the Ferguson rockslide, Merced River Canyon, California

INTRODUCTION An active rockslide in Merced River Canyon was first noticed on April 29, 2006 when a few rocks rolled onto Highway 140 between mileposts 103 and 104, compromising traffic on this highway and signaling the onset of renewed activity of the Ferguson rockslide. State highway 140 is one of the main entrances to Yosemite National Park and is the primary road for large commercial trucks access into the park from the west. Continued rockslide activity during 2006 built a large talus cone that covered the highway and encroached into the Merced River below it. Observations by the US Forest Service (USFS), the California Department of Transportation (CALTRANS), and the U.S. Geological Survey (USGS) confirm that the rockslide remained active through 2006 and represents a potential threat to traffic along the rerouted highway as well as to recreational users of the Merced River in the runout path below the rockslide. Delineation of the hazards posed by the Ferguson rockslide is a necessary prerequisite to mitigating them. Field observations of the rockslide, shown in the photo of Figure 1, have constrained the geometry and structure of the slide mass (Beck, 2006; Gallegos and DeGraff, 2006). Based on initial estimates by geologists from USFS, CALTRANS and the USGS, the rockslide, active in 2006 and 2007, has an area of approximately 40,000 square meters and a volume of approximately 800,000 cubic meters. Structural mapping suggests that the motion of the slide is translational along a planar bed, and that differential motion of the slide from the toe to the headwall has resulted in formation of large tension cracks that transect the slide across the slope (Beck, 2006). These indications of persistent movement were confirmed during 2006 and 2007 by GPS measurements made by the USGS at three points on the rockslide (Rick LaHusen, USGS, written communication). The larger of these cracks divide the slide into regions that moved at different rates in 2006, with the toe of the rockslide moving five to ten times faster than the middle portion or headwall part of the slide. Downslope of the main rockslide mass, a talus slope consist in of angular blocks ranging in size from 0.1 to greater than 10 meters (Gallegos and DeGraff, 2006), buries Highway 140. Both the main rockslide and the talus consist of angular blocks ranging in size from 0.1 to more than 10 meters and are composed of highly fractured phyllite, slate, and chert from the Phyllite and Chert of Hite Cove (Bateman and Krauskopf, 1987). The purpose of this report is to assess the hazard posed by the Ferguson rockslide by simulating the runout and deposition of a portion of the slide if rapid failure occurs. As discussed by Gallegos and DeGraf (2006), a runout analysis is needed to delineate slide hazards. The report is restricted to calculations of potential runout and does not address the likelihood of rapid failure. Based on discussions with Allan Gallegos (USFS), two end-member initial slide volumes were chosen: (1) the toe of the slide along boundaries defined by Tim Beck (CALTRANS) in (Beck, 2006), and (2) the entire sliding rock mass, again along boundaries defined by Tim Beck. The simulated runout of these volumes during rapid failure uses granular flow mechanics developed by Iverson (1997) and the model developed by Denliner and Iverson, (2004). This model has been thoroughly tested against experimental data and provides plausible, defensible results.

Open-File Report↗
USGS2021

Forest resistance to extended drought enhanced by prescribed fire in low elevation forests of the Sierra Nevada

Prescribed fire reduces fire hazards by removing dead and live fuels (small trees and shrubs). Reductions in forest density following prescribed fire treatments (often in concert with mechanical treatments) may also lessen competition so that residual trees might be more likely to survive when confronted with additional stressors, such as drought. The current evidence for these effects is mixed and additional study is needed. Previous work found increased tree survivorship in low elevation forests with a recent history of fire during the early years of an intense drought (2012 to 2014) in national parks in the southern Sierra Nevada. We extend these observations through additional years of intense drought and continuing elevated tree mortality through 2017 at Sequoia and Kings Canyon National Parks. Relative to unburned sites, we found that burned sites had lower stem density and had lower proportions of recently dead trees (for stems ≤47.5 cm dbh) that presumably died during the drought. Differences in recent tree mortality among burned and unburned sites held for both fir (white fir and red fir) and pine (sugar pine and ponderosa pine) species. Unlike earlier results, models of individual tree mortality probability supported an interaction between plot burn status and tree size, suggesting the effect of prescribed fire was limited to small trees. We consider differences with other recent results and discuss potential management implications including trade-offs between large tree mortality following prescribed fire and increased drought resistance.

California↗
USGS2026

Suitability of ground-motion models for seismic hazard assessment in Puerto Rico and the U.S. Virgin Islands

We perform linear mixed-effects analyses with a ground-motion dataset to evaluate how well ground-motion models (GMMs) fit active crustal, subduction interface, and subduction intraslab earthquakes in Puerto Rico and the US Virgin Islands. Most of the GMMs reproduce the spatial variation in peak ground motions with earthquake magnitude and rupture distance but predict ground motions 0.3–1.0 natural log units (35%–270%) greater than observations. Two GMMs developed for Puerto Rico that are based on ground-motion records from mostly magnitude 4–5 earthquakes do not perform as well. We attribute the period-dependent overprediction in observed ground motions to differences between the observed site response and the linear site response in the GMMs. Consequently, we developed region-specific GMMs by adjusting the period-dependent linear site response coefficients and period-dependent constant coefficients to remove most of the bias between predicted ground motions and observations. For the analysis, we compile ground-motion records and process waveforms to build a dataset with 10,127 records at 72 stations from 849 magnitude 4.0–6.4 earthquakes between 1 January 2006 and 31 March 2024. The earthquakes include active crustal, subduction interface, and subduction intraslab events. We evaluate the GMMs using the time-averaged shear wave speed in the top 30 m ( ), which we compile from site surveys and proxy values computed from horizontal to vertical spectral ratios. Site terms exhibit strong consistency across GMMs and crustal and subduction earthquakes, indicating that the linear mixed-effects analysis successfully isolates the effects of local site response. The event terms show little spatial correlation and more substantial variability than in other regions, which we attribute to uncertainties in the earthquake magnitudes. This analysis guides the selection of the GMMs for the 2025 update of the National Seismic Hazard Model for Puerto Rico and the US Virgin Islands.

Puerto Rico, U.S. Virgin Islands↗
USGS2019

Exposure and potential effects of pesticides and pharmaceuticals in protected streams of the US National Park Service southeast Region

Globally protected areas offer refugia for a broad range of taxa including threatened and endangered species. The United States National Park Service (NPS) manages public lands to preserve biodiversity, but increasing park visitation and development of surrounding landscapes increase exposure to and effects from bioactive contaminants. The risk (exposure and hazard) to NPS protected-stream ecosystems within the highly urbanized southeast region (SER) from bioactive contaminants was assessed in five systems based on 334 pesticide and pharmaceutical analytes in water and 119 pesticides in sediment. Contaminant mixtures were common across all sampled systems, with approximately 24% of the unique analytes (80/334) detected at least once and 15% (49/334) detected in half of the surface-water samples. Pharmaceuticals were observed more frequently than pesticides, consistent with riparian buffers and concomitant spatial separation from non-point pesticide sources in four of the systems. To extrapolate exposure data to biological effects space, site specific cumulative exposure-activity ratios (ΣEAR) were calculated for detected surface-water contaminants with available ToxCast data; common exceedances of a 0.001 ΣEAR effects-screening threshold raise concerns for molecular toxicity and possible, sub-lethal effects to non-target, aquatic vertebrates. The results illustrate the need for continued management of protected resources to reduce contaminant exposure and preserve habitat quality, including prioritization of conservation practices (riparian buffers) near stream corridors and increased engagement with upstream/up-gradient property owners and municipal wastewater facilities.

Alabama, Florida, Georgia, Kentucky, Mississippi, ↗
USGS2026

Wavelet Inversion for SliP (WISP): Open-source earthquake slip modeling software

Models of the spatiotemporal evolution of earthquake slip, termed finite-fault models, are a critical component of rapid earthquake and tsunami response, earthquake forecasting, seismic ground-motion estimates, and studies of earthquake kinematics. Here, we detail a newly released finite-fault modeling software, Wavelet Inversion for SliP (WISP), in use at the U.S. Geological Survey’s National Earthquake Information Center (NEIC) and available to the public. WISP version 1.1.0 allows inversion of teleseismic body and surface waves, as well as local strong-motion, static and dynamic Global Navigation Satellite System, and satellite imagery (e.g., Interferometric Synthetic Aperture Radar) observations on single or multiple planar fault segments. The software is used in NEIC rapid response of earthquakes M w ≥ 7, generally resulting in a published model within the first few hours after the event origin time. The rupture location and dimensions are then used as inputs to downstream products to estimate earthquake shaking, predict loss, and model the likelihood of secondary hazards, namely landslides and liquefaction. WISP is also used in research studies to evaluate the characteristics of complex ruptures including multifault ruptures and earthquake doublets, among others. The WISP version 1.1.0 software release is composed of Python-wrapped FORTRAN code to accomplish the inversion procedure. A simple command line interface facilitates ease of use even for those with only a cursory knowledge of Python scripting. WISP version 1.1.0 includes a Jupyter Notebook tutorial demonstrating use of the software for modeling the 2015 M w 8.3 Illapel, Chile, earthquake. In parallel with the tutorial, we demonstrate the typical usage of the WISP software using the M w 8.3 Illapel earthquake example here.

Seismological Research Letters↗