USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Hydrology, Regional Studies”

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,662 records · Page 45Linked to original sources
USGS2024

Cross section N–N' through the Valley and Ridge province of the southern Appalachian basin, from Greene County, west-central Alabama, to Bibb County, central Alabama

Introduction Geologic cross section N–N′ is the sixth in a series of geologic cross sections constructed by the U.S. Geological Survey to document and improve understanding of the geologic framework and petroleum systems of the Appalachian basin. Cross section N–N′ provides a regional view of the structural and stratigraphic framework of the Appalachian basin in the Valley and Ridge province in western and central Alabama; it spans approximately 69 miles (mi) (111 kilometers [km]). This geologic cross section is a companion to geologic cross sections E–E′ , D–D′ , C–C′ , I–I′ , and A–A′ that are located approximately 350 to 550 mi (563 to 885 km) to the northeast. Cross section N–N' complements earlier geologic cross sections through the Alabama part of the Appalachian basin. Although some of the other cross sections show more structural and stratigraphic detail, they are of more limited extent geographically and stratigraphically. Cross section N–N′ contains information that is useful for evaluating energy resources in the Appalachian basin. Although the Appalachian basin petroleum systems are not shown on the cross section, many of their key elements (such as source rocks, reservoir rocks, seals, and traps) can be inferred from lithologic units, unconformities, and geologic structures shown on the cross section. Other aspects of petroleum systems (such as the timing of petroleum generation and petroleum migration pathways) may be evaluated by burial history, thermal history, and fluid flow models based on what is shown on the cross section. In addition, cross section N–N′ may be used as a reconnaissance tool to identify plausible geologic structures and strata for the subsurface storage of liquid waste or for the sequestration of carbon dioxide.

Alabama↗
USGS2025

A comprehensive geologic framework of the National Crustal Model for seismic hazard studies in the conterminous United States

A three-dimensional (3D) geologic framework has been developed for the conterminous United States (U.S.) as part of the U.S. Geological Survey National Crustal Model to enhance seismic hazard modeling. The geologic framework is created from geologic maps and multiple subsurface geologic unit boundaries including the base of the Miocene, Cenozoic, Phanerozoic, and the Mohorovičić discontinuity. Modifications are made to surficial geologic maps to remove discontinuities across state and country borders. The subsurface distribution of rock type and age is extrapolated from the surface, seeded with subsurface geologic information, and constrained by a map of basement geology. The framework provides the basis for estimates of subsurface seismic velocity and density that is needed to improve estimates of earthquake ground shaking and seismic hazard. The present framework greatly expands and updates a previously published 3D geologic framework of the western part of the U.S. that was itself a first-of-its-kind digital 3D portrayal of the nation.

conterminous United States↗
USGS2026

Investigating the influence of climate and volcanic surface aging on fluvial erosion: A case study of Réunion Island, Indian Ocean

Precipitation is one of the dominant drivers of landscape erosion and evolution; however, the effects of typical rainfall compared with less frequent, high-magnitude precipitation events on erosion remain unclear. Volcanic islands are ideal locations to study such phenomena due to their simple geometries, nontectonic construction, and strong spatiotemporal rainfall gradients. However, spatial variation in surface age, created during their construction, often complicates their degradation histories by introducing temporal changes in erosion rates as drainage networks develop. Réunion Island (western Indian Ocean) presents a clear example of this, with an east–west gradient in both surface age and mean annual precipitation, as well as infrequent cyclones that alter the background rainfall pattern. In this study, we analyze the effects of surface age, average rainfall, and rainfall variability on basin development and fluvial erosion across the island. We calculate basin-averaged values of basin morphology, age, precipitation, river discharge, eroded volumes, and erosion rates, and use these to analyze the dominant drivers of landscape evolution through a series of correlation analyses. Our results indicate a temporal dependence on the influence of precipitation, with young surfaces being dominantly eroded by high-rainfall events and older surfaces eroded by mean annual rainfall patterns. Furthermore, we show that drainage development of shield volcanoes follows similar trends to other volcano types, and suggest that surface permeability and groundwater structure are important controls on runoff-driven erosion on shield volcanoes. These results add new components to the question of how precipitation impacts erosion.

Réunion Island↗
USGS2025

Using public participatory geographic information systems (PPGIS) to explore uses and values for Mojave Trails National Monument, California

Many people ascribe a variety of values to public lands and waters, but some values are more difficult to assess and quantify than others. Public participatory geographic information systems (PPGIS) are tools that have been used to help quantify and map the public’s diverse values for a landscape. This work describes the first known Office of Management and Budget–approved use of PPGIS by a Department of the Interior bureau. The U.S. Geological Survey developed an internet-based application to aid in gathering PPGIS data, called Values Mapping for Planning in Regional Ecosystems (VaMPIRE). Further, this work describes the first pilot of the VaMPIRE application in coordination with the Bureau of Land Management to collect spatial data and other survey data regarding the public’s uses of and values for locations within Mojave Trails National Monument. We emailed the link to the VaMPIRE application to an interested party email list in 2024 with 207 valid emails and received 74 responses; we also received 47 responses from members of an off-roading social media group. Of the list of 16 value options, recreation was the most popular value for the monument, followed by wilderness and inspirational. Over 1,000 points were placed throughout the monument, indicating locations people use or value, with the locations spread throughout the entire monument. Additionally, most survey respondents stated their ability to receive benefits in locations they mapped would not change in response to a hypothetical scenario related to recreational facility development. This report describes exploratory results from the first use of the VaMPIRE tool in Mojave Trails National Monument and includes reflections on how the process went and considerations for future use of VaMPIRE.

California↗
USGS2026

Long-term intermittent connection between the western Snake River Plain and Columbia basin: A two-phased incision history of Hells Canyon

For more than a century, researchers have debated the route of the Snake River across the northern Cordillera and U.S. Pacific Northwest, including the associated incision history of Hells Canyon. Here, we use detrital zircon U-Pb provenance analysis of Miocene strata upstream and downstream of Hells Canyon to constrain the evolution in fluvial pathways over time. Downstream of Hells Canyon, we show that the ca. 10−6 Ma Clarkston Heights gravel was dominantly sourced from the nearby Salmon and Clearwater rivers but also with substantial input from drainages that now flow into the western Snake River Plain. These same western Snake River Plain drainages were contributing to time-equivalent Lake Idaho strata, upstream of Hells Canyon, indicating a Hells Canyon fluvial connection between the western Snake River Plain and Columbia basin before 6 Ma. Results from Lake Idaho strata suggest a two-phased lacustrine history in the western Snake River Plain: Before ca. 6.0 Ma, lake strata were derived from local tributaries, suggesting that the western Snake River Plain was isolated from the modern upper Snake River system. Between ca. 4.3 Ma and 2.2 Ma, the source area expanded to include eastern tributaries in conjunction with migration of high-standing topography of the Yellowstone hotspot. Our results challenge the long-held hypothesis of “capture” of the Snake River through Hells Canyon. Instead, we infer that Hells Canyon was a long-established route for outflow of lakes occupying the western Snake River Plain, possibly intermittently, during a ca. 10−2 Ma phase of slow regional incision, followed by rapid incision in Hells Canyon starting ca. 2 Ma in conjunction with erosion of the topographic barrier impounding Lake Idaho.

Idaho, Oregon↗
USGS2024

Post-glacial stratigraphy and late Holocene record of great Cascadia earthquakes in Ozette Lake, Washington, USA

Ozette Lake is an ~100-m-deep coastal lake located along the outer coast of the Olympic Peninsula (Washington, USA); it is situated above the locked portion of the northern Cascadia megathrust but also relatively isolated from active crustal faults and intraslab earthquakes. Here we present a suite of geophysical and geological evidence for earthquake-triggered mass transport deposits (MTDs) and related turbidite deposition in Ozette Lake since ca. 14 ka. Comprehensive high-resolution bathymetry data, seismic reflection profiles, and sediment cores are used to characterize the post-glacial stratigraphic framework and examine paleoseismic evidence in the lacustrine sediments. Stacked sequences of MTDs along the steep eastern flanks of the lake appear to grade basin-ward from thick, chaotic, blocky masses to thin, parallel-bedded turbidite beds. The discrete turbidite event layers are separated by fine-grained (silt and clay) lake sedimentation. The event layers are observed throughout the lake, but the physical characteristics of the deposits vary considerably depending on proximity to primary depocenters, steep slopes, and subaqueous deltas. A total of 30–34 event deposits are observed in the post-glacial record. Radiometric dating was used to reconstruct a detailed sedimentation history over the last ~5.5 k.y., develop an age model, and estimate the recurrence (365–405 yr) for the most recent 12 event layers. Based on sedimentological characteristics, temporal overlap with other regional paleoseismic chronologies, and recurrence estimates, at least 10 of the dated event layers appear to be sourced from slope failures triggered by intense shaking during megathrust ruptures; the recurrence interval for these 10 events is 440–560 yr. Thus, Ozette Lake contains one of the longest and most robust geological records of repeated shaking along the northern Cascadia subduction zone.

Washington↗
USGS2024

Recent applications of the USGS National Crustal Model for Seismic Hazard Studies

The U.S. Geological Survey is developing the National Crustal Model (NCM) for seismic hazard studies to facilitate modeling site, path, and source components of seismic hazard across the conterminous United States. The NCM is composed of a 1km grid of geophysical profiles, extending from the Earth’s surface into the upper mantle. It is constructed from a threedimensional (3D) geologic framework and geophysical rules that use (1) a petrologic and mineral physics database; (2) a 3D temperature model; and (3) a calibrated rock type- and age-dependent porosity model. Parameters needed to estimate site response for existing ground motion models (GMMs), including the time-averaged velocity in the upper 30 meters (VS30), the depths to 1.0 and 2.5 km/s shear-wave velocity (Z1.0 and Z2.5), and sediment thickness, can be computed from the NCM. As GMMs continue to improve in the future, other metrics could also be extracted or derived from the NCM, such as fundamental period, site attenuation (ko), a fully frequency-dependent site response function, or 3D geophysical volumes for wavefield simulations. Application of the NCM may also benefit other aspects of seismic hazard analysis, including better accounting for path-dependent attenuation and geometric spreading, more accurate estimation of earthquake source properties such as hypocentral location and stress drop, and calculation of crustal strength profiles that inform estimates of the base of seismicity.

conterminous United States↗
USGS2024

A global assessment of SAOCOM-1 L-band stripmap data for InSAR characterization of volcanic, tectonic, cryospheric, and anthropogenic deformation

SAOCOM-1 is an L-band (23.5 cm) synthetic aperture radar (SAR) constellation made up of two satellites launched in 2018 and 2020 by Comisión Nacional de Actividades Espaciales (CONAE, Argentina). In this contribution, we present a global summary of interferometric SAR (InSAR) observations of ground deformation with SAOCOM-1 stripmap data for tracking volcanic, tectonic, glacier, and anthropogenic deformation. These examples include: 1) episodes of unrest at volcanoes in the Aleutian Islands, Southern Andes, and Italy, with line-of-sight (LOS) deformation from 4 cm/yr in InSAR time series to ~70 cm in interferograms; 2) dike intrusions in Hawai’i; 3) earthquakes in the Andean fold and thrust belt and the East Anatolian fault; 4) ice flow of the Southern Patagonia icefield; and 5) subsidence due to lithium brine extraction in the Salar de Atacama basin (northern Chile). Comparisons between SAOCOM-1, ALOS-2 SM3, Sentinel-1, and TerraSAR-X/ TanDEM-X/PAZ (TSX/TDX/PAZ) mean velocities from InSAR time series show a 1:1 ± 3% correlation in the LOS velocity, which highlights the high accuracy of SAOCOM-1 data. The minimum deformation that we measured in individual interferograms is 4 ± 0.6 cm. One limitation of SAOCOM-1 is the lack of a global acquisition program, which reduces its global and broader applications. Considering the repeat periods, background observation program, and lack of a controlled orbital tube, the best suited targets for SAOCOM-1 InSAR are two. First, volcanoes that deform with secular rates located in vegetated regions in mid- and high-latitudes, and/or that undergo transient episodes of fast deformation in which C-band coherence is lost quickly. Second, glaciers where coherence can be sustained during the repeat period of eight days.

IEEE Transactions on Geoscience and Remote Sensing↗
USGS1999

Selected aquatic biological investigations in the Great Salt Lake basins, 1875-1998, National Water-Quality Assessment Program

This report summarizes previous investigations of aquatic biological communities, habitat, and contaminants in streams and selected large lakes within the Great Salt Lake Basins study unit as part of the U.S. Geological Survey’s National Water-Quality Assessment Program (NAWQA). The Great Salt Lake Basins study unit is one of 59 such units designed to characterize water quality through the examination of chemical, physical, and biological factors in surface and ground waters across the country. The data will be used to aid in the planning, collection, and analysis of biological information for the NAWQA study unit and to aid other researchers concerned with water quality of the study unit. A total of 234 investigations conducted during 1875-1998 are summarized in this report. The studies are grouped into three major subjects: (1) aquatic communities and habitat, (2) contamination of streambed sediments and biological tissues, and (3) lakes. The location and a general description of each study is listed. The majority of the studies focus on fish and macroinvertebrate communities. Studies of algal communities, aquatic habitat, riparian wetlands, and contamination of streambed sediment or biological tissues are less common. Areas close to the major population centers of Salt Lake City, Provo, and Logan, Utah, are generally well studied, but more rural areas and much of the Bear River Basin are lacking in detailed information, except for fish populations.

Idaho, Utah, Wyoming↗
USGS2025

Overview of The SCEC/USGS Community Stress Drop Validation Study using the 2019 Ridgecrest earthquake sequence

We present initial findings from the ongoing Community Stress Drop Validation Study to compare spectral stress‐drop estimates for earthquakes in the 2019 Ridgecrest, California, sequence. This study uses a unified dataset to independently estimate earthquake source parameters through various methods. Stress drop, which denotes the change in average shear stress along a fault during earthquake rupture, is a critical parameter in earthquake science, impacting ground motion, rupture simulation, and source physics. Spectral stress drop is commonly derived by fitting the amplitude‐spectrum shape, but estimates can vary substantially across studies for individual earthquakes. Sponsored jointly by the U.S. Geological Survey and the Statewide (previously, Southern) California Earthquake Center our community study aims to elucidate sources of variability and uncertainty in earthquake spectral stress‐drop estimates through quantitative comparison of submitted results from independent analyses. The dataset includes nearly 13,000 earthquakes ranging from M 1 to 7 during a two‐week period of the 2019 Ridgecrest sequence, recorded within a 1° radius. In this article, we report on 56 unique submissions received from 20 different groups, detailing spectral corner frequencies (or source durations), moment magnitudes, and estimated spectral stress drops. Methods employed encompass spectral ratio analysis, spectral decomposition and inversion, finite‐fault modeling, ground‐motion‐based approaches, and combined methods. Initial analysis reveals significant scatter across submitted spectral stress drops spanning over six orders of magnitude. However, we can identify between‐method trends and offsets within the data to mitigate this variability. Averaging submissions for a prioritized subset of 56 events shows reduced variability of spectral stress drop, indicating overall consistency in recovered spectral stress‐drop values.

California↗
USGS2026

Effect of mineral deposit data on predictions from the three-part approach to quantitative mineral resource assessment—A study of 16 previous U.S. Geological Survey assessments

The three-part approach to quantitative mineral resource assessment requires information about the properties of undiscovered mineral deposits in an assessment area. These properties are unknown, so the properties of discovered mineral deposits of the same mineral deposit type are used instead. In the three-part approach, these discovered mineral deposits come from around the world, and their properties constitute the pooled data for that mineral deposit type. Alternatively, these discovered mineral deposits could come from the assessment area, and their properties constitute the tract data for that mineral deposit type. Tract data may be more representative of the undiscovered mineral deposits in the assessment area than the pooled data. The goal of this study was to determine whether resource predictions using pooled data are equivalent to resource predictions using tract data. To this end, 16 previous U.S Geological Survey assessments were studied. For each assessment, resources were predicted for one undiscovered mineral deposit in the assessment area. One set of predictions used pooled data, and another used tract data. The two sets of predictions were compared with an equivalence test, using the six assessment statistics that are commonly reported for mineral resource assessments. Practical equivalence is the condition that two corresponding assessment statistics are within a factor of 1.5 of one another. For each of 2 assessments, all 6 assessment statistics were practically equivalent. For both assessments, the assessment statistics from the pooled data, relative to the corresponding assessment statistics from the tract data, ranged from 1.30 times smaller to 1.03 times larger. For each of 14 assessments, 1 or more of the 6 assessment statistics were not practically equivalent. The assessment statistics from the pooled data, relative to the corresponding assessment statistics from the tract data, ranged from 26.6 times smaller to 5.53 times larger. The use of pooled data has been a standard procedure in the three-part approach since at least 1986. The 16 assessments in this study are not a representative sample of those prior assessments that used pooled data. So, it is inappropriate to use the study results to infer whether pooled data affected the resource predictions for those prior assessments.

Scientific Investigations Report↗
USGS2024

SCEC/USGS Community Stress Drop Validation Study: How spectral fitting approaches influence measured source parameters

Spectral source parameters used to estimate an earthquake’s stress drop (Δσ) can vary significantly across measurement approaches. The Statewide California Earthquake Center/U.S. Geological Survey Community Stress-Drop Validation Study was initiated to compare source parameter estimates, focusing initially on a dataset from the 2019 Ridgecrest earthquake sequence. As part of that validation effort, here we focus on one potential source of uncertainty: whether spectral fitting approaches alone, applied to a common set of spectra from the 2019 Ridgecrest sequence result in different source parameter estimates. By using a common set of benchmark spectra analyzed across a consistent frequency band of 1–40 Hz, we eliminate many sources of variability. A subgroup of validation study participants volunteered to estimate the low-frequency displacement (Ω0) and corner frequency ( f c ) by fitting a smooth function to benchmark displacement spectra. Participants used linear- or log-sampled spectra, assumed a Brune or Boatwright spectral model, and applied different misfit criteria. We compare 17 approaches used to estimate Ω0, f c , and Δσ for 54 earthquake spectra. Our results reveal that 35% of events have Δσ estimates within a factor of two, whereas others exhibit variations exceeding an order of magnitude. The variability in Ω0 and f c can largely be attributed to whether a spectrum is consistent with the smooth function of an idealized simple crack model. The trade-off between Ω0 and f c may be more pronounced when using linearly sampled spectra, as higher frequency spectral bumps control the fits. As expected, methods that assumed a Boatwright model tended to have lower Ω0 and somewhat higher f c compared to those assuming a Brune model, although resulting Δσ estimates are similar. When compared to the overall validation study results, the fitting approach alone may account for between 5% and 90% (25% on average) of the total variability in spectral Δσ.

California↗
USGS1976

Water-quality study of a reach of the Merced River in Yosemite National Park and vicinity, California, April 1973 through September 1974

A 25-mile (40-kilometre) reach of the Merced River was studied to provide the National Park Service with water-quality data prior to the operation of a new sewage-treatment plant. Data were collected periodically at four Merced River and two inflow stations. Discharge varied seasonally with highest flow occurring at the time of spring snowmelt runoff. Lowest flow occurred during late summer, coinciding with high park visitor use and elevated water temperature. Water-quality variables investigated were discharge, temperature, major chemical constituents, selected trace metals, specific conductance, suspended sediment, turbidity, dissolved oxygen, pH, alkalinity, carbon dioxide, total organic carbon, nitrogen, phosphorus, algal growth potential, coliform bacteria, periphyton, and benthic invertebrates. Diel measurements were made at some stations to delineate daily fluctuations in selected water-quality variables. Results indicate that water in the reach of the Merced River studied was of good quality. The river had low concentrations of dissolved solids (8 to 41 milligrams per litre) and low alkalinity values (2 to 24 milligrams per litre of alkalinity as CaCO 3 ). Mean concentrations of nitrogen and phosphorus were highest in the Merced River at Rancheria Flat. Algal growth potential tests showed that phosphorus was limiting algal growth at Happy Isles Bridge, El Capitan Bridge, and South Fork Merced River, whereas nitrogen was limiting growth at Rancheria Flat. Input of excessive nitrogen and phosphorus concentrations in the Merced River during the summer low-flow period, and mixing with warm tributary inflow could cause excessive algal growth and production. The algal types found were diatoms, green algae, and blue-green algae. Biomass/chlorophyll a ratios in the Merced River were low (0.3 to 7.7). Benthic invertebrates consisted chiefly of Diptera, Ephemeroptera, and Plecoptera.

California↗
USGS2024

A partially nonergodic ground-motion model for Fourier amplitude spectra for the San Francisco Bay area, California, USA

We develop a partially nonergodic ground-motion model (GMM) for Fourier amplitude spectra for the San Francisco Bay Area, California, USA, using the Bayless and Abrahamson (2019) GMM as a reference ergodic GMM and developing location-dependent adjustments to the predicted median and variance. We compile regional ground-motion data from moment magnitude (𝑀 w ) >3 earthquakes occurring during 2000–2022 for which magnitude information is available in the U.S. Geological Survey Comprehensive Catalog (Guy et al., 2015). The data set predominantly consists of records from 𝑀 w 3.5–4.5 earthquakes but includes three well-recorded 𝑀 w > 5 events. Ground-motion residuals are evaluated using the time-averaged shear-wave velocity in the top 30 m (𝑉 S30 ) from the California-specific map of Thompson et al. (2018) and basin-depth site parameters from the seismic velocity model of Aagaard and Hirakawa (2021). The 𝑉 S30 dependence and basin-depth scaling of the reference ergodic GMM of Bayless and Abrahamson (2019) are evaluated and modified with the updated data set. We compute maps of site adjustments using a varying-coefficient model that considers the spatial correlation structure and uncertainties at each observation location. The spatial covariance model is developed using ground-motion residuals that are standardized by the uncertainty model, which allows for consideration of the aleatory variability in developing the site adjustments. The covariance model is fit considering the means and standard deviations of the site terms at all locations. The use of partially nonergodic median adjustments results in modified variance components of the within-event variability. Due to the low number of large-magnitude earthquakes that control seismic hazard in the data set, we do not modify between-event variance; however, we present adjustments to site-to-site variability for use in partially nonergodic hazard assessments.

California↗
USGS2025

Phase 1 technical implementation plan for the expansion of the ShakeAlert earthquake early warning system to Alaska

Executive Summary The conference report accompanying the fiscal year (FY) 2022 Consolidated Appropriations Act (Public Law 117–103) for the U.S. Department of the Interior and related agencies directed the U.S. Geological Survey (USGS) to “work with the State of Alaska to develop an implementation plan to be completed within two years in order to put ShakeAlert/Earthquake Early Warning in Alaska” (p. 29). Congress included $1 million in the FY 2022 appropriation to conduct this effort. The USGS Earthquake Hazards Program, along with partner organizations, has developed the ShakeAlert earthquake early warning (EEW) system for the West Coast, which currently operates in California, Oregon, and Washington. The purpose of the system and its alert delivery partners is to reduce the impact of earthquakes and save lives and property by delivering ShakeAlert-powered alerts that are transmitted to the public via mass notification technologies, and by providing more detailed data streams to institutional users and commercial service providers to trigger automated, user-specific, protective actions. ShakeAlert was designed in such a way that it could be expanded to other U.S. regions with high earthquake risk, after the build-out of seismic and geodetic networks to support ShakeAlert in a specified region is completed and the necessary funding is secured for long-term operation and maintenance. When an earthquake occurs, seismic waves radiate from the rupturing fault like waves on a pond. It is these waves that people feel as earthquake shaking and that can cause damage to structures. Using networks of ground-motion sensors and sophisticated computer algorithms, ShakeAlert can detect an earthquake seconds after it begins, calculate its location and magnitude, and estimate the resulting intensity of shaking. Early warnings of impending shaking are then sent to people and systems that may experience damaging shaking, allowing them to take appropriate protective actions. Depending on the user’s distance from the earthquake, alerts may be delivered before, during, or after the arrival of strong shaking. There will almost always be a region near the earthquake epicenter where alerts arrive after damaging shaking has begun. The ShakeAlert system updates its ground-motion estimates as an earthquake grows larger. In response to the FY 2022 congressional direction, the USGS worked with the State of Alaska to devise this implementation plan for ShakeAlert expansion to Alaska. The USGS engaged with the Alaska Division of Homeland Security and Emergency Management (DHS&EM) and the Alaska Division of Geological and Geophysical Surveys (DGGS). A cooperative agreement was awarded to the Alaska Earthquake Center (AEC) at the University of Alaska Fairbanks (UAF) for their contributions to the plan and their work coordinating with other networks in Alaska. The USGS engaged with the Alaska Seismic Hazards Safety Commission (ASHSC) throughout the process. The USGS also held a series of Alaska stakeholder engagements. The process of developing the implementation plan was facilitated by contracted staff from Corner Alliance, which is a government consulting firm. This implementation plan describes the details and estimates the costs for a Phase 1 expansion of the ShakeAlert system to Alaska. A geographically limited Phase 1 goal was chosen that covers the highest risk and most populated areas of Alaska. The areas proposed encompass the State’s main population centers and 90 percent of the State’s population. This Phase 1 design is considered very challenging and ambitious from the viewpoint of network operators. The lessons learned if this plan is implemented could be used to consider subsequent phases to expand EEW beyond Phase 1 in Alaska in the future. ShakeAlert is built on the foundation of the sensor networks and data processing infrastructure of the USGS-led Advanced National Seismic System (ANSS). This implementation plan calls for a total of 450 high-quality, real-time EEW-capable ANSS seismic stations in Alaska: 270 new stations, 160 upgraded stations, and 20 existing stations. These seismic station numbers are based on a station spacing of 10 kilometers (km) in urban areas, 20 km in seismic source areas that endanger population centers, and 40 km in other areas. The associated costs also include support for some EEW-capable global navigation satellite system (GNSS) stations, with a focus on improving warnings for large subduction zone earthquakes. For effective EEW, ShakeAlert requires low-latency, high-availability, robust telemetry links to deliver continuous, real-time data from field stations to the data centers. The Alaska data processing hardware infrastructure would follow the general design for fail-safe operation that is used for the ShakeAlert system on the West Coast. The ShakeAlert architecture uses two independent layers: the production layer for earthquake processing and the alert layer to make alerting decisions and serve alerts to users. This implementation plan includes two geographically separated data centers in Alaska, each with two fully independent production and alert layers using the same system design developed for the West Coast. As of March 2024, the ShakeAlert system is at version 3.0.1, with more advanced versions in the development and testing pipeline. ShakeAlert originally used two algorithms to determine the location and magnitude of earthquakes using seismic data. A third algorithm that can calculate very large magnitudes of very large earthquakes with geodetic data was added in March 2024. ShakeAlert publishes several data and alert products to meet the needs of different users. All messages include the location of the earthquake, either as a point or a line, and its magnitude. Ground-shaking estimates are published in two forms, as ground-motion contours and a map grid. Providing adequate warning time for strong shaking (the “target threshold”) requires sending alerts at a threshold lower than that strong shaking level (the “alert threshold”). The thresholds for public alerting in Alaska would be a joint USGS and State decision. To have the greatest benefit, ShakeAlert-powered alerts would be delivered to institutional users and individuals by all practical pathways. The USGS alert layer can support thousands of institutional users and alert redistributors, but the USGS does not have the mission nor the infrastructure and expertise to perform mass notifications to the public or implement automatic actions for end users of the alerts. To meet this need, ShakeAlert recruits private sector “technology enablers” that have the necessary expertise to develop end-user implementations using EEW alerts with the goal of stimulating an EEW industry. Earthquake early warning alerts are useless if people do not know how to respond to them. Although the alert messages include instructions about what to do (drop, cover, and hold on), alerts are more effective if people have been trained in advance. Messages about ShakeAlert’s capabilities, limitations, and benefits could be integrated with existing earthquake education programs, including State-run programs. Therefore, ShakeAlert would coordinate with both public and private partners and stakeholders through various partnerships and agreements to accomplish consistent and ongoing public earthquake hazard education. The estimated capital cost of completing the computing infrastructure and sensor networks for the Phase 1 ShakeAlert expansion to Alaska is approximately $66 million in 2024 dollars. The annual operation and maintenance cost of the completed system is estimated to be $12 million per year in 2024 dollars when fully built out.

Alaska↗
USGS2025

Combining multisite tsunami and deformation modeling to constrain slip distributions for the 1700 C.E. Cascadia earthquake

A major earthquake ruptured the Cascadia subduction zone (CSZ) on 26 January 1700. Key paleoseismic evidence associated with this event include tsunami deposits, stratigraphic evidence of coastal coseismic subsidence, written Japanese records of a tsunami unaccompanied by earthquake shaking, and margin‐wide turbidites found offshore and in lacustrine environments. Despite this wealth of independent clues, important details about this event remain unresolved. Dating uncertainties do not conclusively establish whether the proxies are from one earthquake or a sequence of them, and we have limited knowledge of the likely slip distributions of the event or events. Here, we use a catalog of 37,500 candidate synthetic ruptures between M w 7.8 and 9.2 and simulate their resulting coseismic deformation and tsunami inundation. Each model is then compared against estimated Japan tsunami arrivals, regional coastal subsidence records, and local paleotsunami deposits mapped at six different coastal marshes and one coastal lake along the CSZ. We find that seven full‐margin ruptures with a median magnitude of M w 9.1 satisfy all three constraints. We favor one M w 9.11 model that best matches all site paleoseismic observations and suggests that the Cascadia megathrust slipped up to ∼30 m and must have shallow geodetic coupling. We also find that some sequences composed of three or four ruptures can still satisfy the observations, yet no sequences of two ruptures can. Sequences are differentiated into three groups based on whether they contain a mainshock rupture located in the south (>44° N) or further north. All sequences contain unruptured portions of the megathrust and most contain mainshocks with peak slip above 40 m. The fit of the geologic evidence from sequences is poor in comparison to single‐event models. Therefore, sequences are generally less favored compared to full‐margin events.

British Columbia, California, Oregon, Washington↗
USGS2026

The value of electronic tagging and tracking studies for understanding fish–hypoxia interactions

Objective Hypoxia (i.e., low dissolved oxygen availability) is a natural phenomenon but can also be induced by human activities (e.g., nutrient enrichment from runoff). Given that dissolved oxygen is essential for aquatic life, periods of hypoxia tend to have negative consequences (e.g., sublethal disturbances, mortality) for most fishes. Extensive laboratory research has documented hypoxia thresholds and physiological and behavioral outcomes for a variety of freshwater and marine fishes and there is also an extensive body of fieldwork assessing population-level responses (e.g., survival, distribution). Yet, studying how individual fish respond to hypoxia in the wild has proved challenging; electronic tagging and tracking tools (e.g., biotelemetry, biologging) have made it easier to study individual responses to hypoxia and complement other tools like hydroacoustics that tend to focus on population-level responses. Methods We review what has been learned from contemporary studies that employ electronic tagging and tracking tools to understand how fish respond to hypoxia in the wild. Topics explored include identifying and validating hypoxia thresholds, habitat compression, connectivity, mortality, physiological and bioenergetic consequences, and extreme weather conditions. We also consider what we have learned about evaluating various management plans for hypoxia and reflect on how electronic tags have also been used in aquaculture systems and in hybrid studies that combine laboratory and field research. We highlight fishes in Lake Erie as a research narrative to demonstrate how electronic tagging and tracking have markedly improved our understanding of a longstanding hypoxia issue. Results Our synthesis revealed that electronic tagging and tracking have provided critical information on how fish respond to hypoxia in the field, revealing complex trade-offs and compensatory mechanisms as well as cryptic hypoxia-induced mortality. Beyond just illuminating space use and mortality, tags equipped with various sensors are revealing how fish deal with hypoxia in real time in terms of physiology, bioenergetics, and behavior. Conclusions As electronic tagging and tracking methods experience further innovation and are increasingly applied to understand the effects of hypoxia on fish, we expect more unanticipated findings about the effects of hypoxia on fish in all aquatic ecosystems, which will strengthen our ability to manage and mitigate hypoxia. Combining tools and approaches (e.g., lab and field) is perhaps the best way to generate comprehensive understanding.

Transactions of the American Fisheries Society↗
USGS2025

Deterministic physics-based earthquake sequence simulators match empirical ground-motion models and enable extrapolation to data poor regimes: Application to multifault multimechanism ruptures

We use the deterministic earthquake simulator RSQSim to generate complex sequences of ruptures on fault systems used for hazard assessment. We show that the source motions combined with a wave propagation code create surface ground motions that fall within the range of epistemic uncertainties for the Next Generation Attenuation‐West2 set of empirical models. We show the model is well calibrated where there are good data constraints, and has good correspondence in regions with fewer data constraints. We show magnitude, distance, and mechanism dependence all arising naturally from the same underlying friction. The deterministic physics‐based approach provides an opportunity for better understanding the physical origins of ground motions. For example, we find that reduced stress drops in shallow layers relative to constant stress drop with depth lead to peak ground velocities in the near field that better match empirical models. The simulators may also provide better extrapolations into regimes that are poorly empirically constrained by data because physics, rather than surface shaking data parameterizations, is underlying the extrapolations. Having shown the model is credible, we apply it to a problem where observations are lacking. We examine the case of crustal faults above a shallow subduction interface seen to break coseismically in simulations of the New Zealand fault system. These types of events were left out of consideration in the most recent New Zealand national seismic hazard model due to the modeling complexity and lack of observational data to constrain ground‐motion models (GMMs). Here, we show that in the model, by breaking up the coseismic crustal and interface rupturing fault motions into two separate subevents, and then recombining the resulting ground‐motion measures in a square‐root‐of‐sum‐of‐squares incoherent manner, we reproduce well the ground‐motion measures from the full event rupture. This provides a new method for extrapolating GMMs to more complex multifault ruptures.

Seismological Research Letters↗