USGS Science⌕ Search

SEARCH · USGS Science

Results for “Canadian Journal of Earth Sciences”

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 83Linked to original sources
USGS2025

Flow heterogeneity controls dissolution dynamics in topologically complex rocks

Rock dissolution is a common subsurface geochemical reaction affecting pore space properties, crucial for reservoir stimulation, carbon storage, and geothermal energy. Predictive models for dissolution remain limited due to incomplete understanding of the mechanisms involved. We examine the influence of flow, transport, and reaction regimes on mineral dissolution using 29 time-resolved data from 3D rocks. We find that initial pore structure significantly influences the dissolution pattern, with reaction rates up to two orders of magnitude lower than batch conditions, given solute and fluid-solid boundary constraints. Flow unevenness determines the location and rate of dissolution. We propose two models describing expected dissolution patterns and effective reaction rates based on dimensionless metrics for flow, transport, and reaction. Finally, we analyze feedback between evolving flow and pore structure to understand conditions that regulate/reinforce dissolution hotspots. Our findings underscore the major impact of flow arrangement on reaction-front propagation and provide a foundation for controlling dissolution hotspots.

Geophysical Research Letters↗
USGS2026

Magnitude conversion relations create substantial differences in seismic hazard models

Earthquake catalogs are essential data inputs for seismic hazard modeling. Because earthquake magnitudes are reported in a variety of types (e.g., local magnitudes and moment magnitudes), magnitude conversion relationships must be used to convert the different magnitude types present in a catalog to a uniform magnitude type to avoid biases in the hazard computation. However, these conversion relationships are often uncertain and have been shown to sometimes perform poorly. Here, we investigate the sensitivity of the gridded seismicity component of the National Seismic Hazard Model (NSHM) to the catalog conversion equations in the Eastern United States. In the 2023 NSHM, magnitudes of various types were converted to moment magnitudes using equations developed by the Central and Eastern United States Seismic Source Characterization for Nuclear Facilities (CEUS‐SSCn), based on least‐squares (LS) regressions made using data from a catalog containing events up through 2008. We recompute these equations using events in the Advanced National Seismic System Comprehensive Earthquake Catalog with multiple magnitudes from 2000 to 2023. Although we prefer the use of orthogonal regressions for our datasets, LS regressions produce broadly similar results, with both approaches exhibiting large deviations from the CEUS‐SSCn conversions, especially at smaller magnitudes. We compare the spatial distribution of annual rates using three different models: (1) the 2023 NSHM conversions, (2) our updated conversions, and (3) no conversions. We find that the choice of conversions leads to substantial differences in the rate forecasts, which can greatly impact the seismic hazard model, particularly in regions with low‐seismicity rates such as the Eastern United States, where the hazard is dominated by gridded seismicity rather than a fault model.

Seismological Research Letters↗
USGS2025

Estimating agricultural irrigation water consumption for the High Plains aquifer region with integrated energy- and water-balance evapotranspiration modeling approaches

Estimation of irrigation water use provides essential information for the management and conservation of agricultural water resources. Conventionally, water use data are created based on reports and surveys from water users, whereas manual records may not be complete due to lacking flow meters, measurement gaps, inconsistent methods across regions, and time- and cost-consuming data processing. Alternatively, spatially explicit estimation of irrigation water use can be conducted efficiently using remote sensing evapotranspiration (ET) modeling approaches. In this study, we created a gridded blue water evapotranspiration (BWET) dataset to estimate historical irrigation water consumption (1986 – 2020) in the croplands across the United States High Plains aquifer region. The BWET data were generated by integrating an energy-balance ET model [Operational Simplified Surface Energy Balance model (SSEBop)] and a water-balance ET model [Vegetation ET model (VegET)]. BWET in croplands indicates crop consumptive use of irrigation water extracted from surface water and groundwater resources. The BWET estimates were compared with reported irrigation water use data for all counties within the aquifer region. The results revealed high agreement between growing season (May – September) BWET and annual water withdrawal at county level. Specifically, correlation coefficients of volumetric BWET and water withdrawal were 0.90 and 0.96, respectively, for the entire aquifer region and western Kansas. The timeseries of BWET and water withdrawal showed similar temporal trends and high covariations. The BWET estimates were systematically lower than the water withdrawal measurements, which was primarily attributed to blue water losses in the irrigation system. The irrigation efficiency, calculated as the ratio of BWET to water withdrawal depth, was 0.57 and 0.74 for the entire aquifer region and western Kansas, respectively. This study demonstrates the capability of using satellite-based ET models (e.g., SSEBop and VegET) to efficiently estimate crop water consumption and evaluate irrigation efficiency at landscape, county, and regional scales.

Colorado, Kansas, Nebraska, New Mexico, Oklahoma, ↗
USGS2025

Hydraulic conductivity and transmissivity estimates from slug tests in wells within the Mississippi Alluvial Plain, Arkansas and Mississippi, 2020

During the spring and summer of 2020, the U.S. Geological Survey conducted single-well slug tests on selected observation wells within the Mississippi Alluvial Plain in Arkansas and Mississippi to estimate hydraulic conductivity and transmissivity values for the Mississippi River Valley alluvial and middle Claiborne aquifers. Well and aquifer data were collected from field measurements, well-construction reports, and published aquifer-thickness information. A total of 324 slug-in and slug-out tests were conducted on 48 wells by using mechanical slugs to displace the water column and submersible pressure transducers to record changes in water levels in the wells. Hydraulic conductivity of the aquifers in which the wells are screened was estimated by curve fitting the water-level-change data using aquifer test analysis software. Estimates of aquifer transmissivity were made by multiplying the estimated hydraulic conductivity value by the aquifer thickness at well locations. Mean hydraulic conductivity estimates for 44 observation wells screened in the Mississippi River Valley alluvial aquifer range from 3 to 401 feet per day, and mean transmissivity estimates range from 285 to 80,559 feet squared per day. Mean hydraulic conductivity estimates for four observation wells screened in units of the middle Claiborne aquifer range from 0.14 to 183 feet per day, and mean transmissivity estimates range from 55 to 67,913 feet squared per day. The results from these tests can be used to improve the understanding of water availability and groundwater migration, to refine groundwater models, and to ultimately provide stakeholders and decisionmakers better information for management of the groundwater resources within the Mississippi Alluvial Plain.

Arkansas, Mississippi↗
USGS2025

Persistent deep long-period seismicity near the Lassen Volcanic Center

Deep long-period (DLP) earthquakes have been observed at many volcanic settings around the world and linked to the magmatic processes that drive volcanic unrest. At the Lassen Volcanic Center (LVC) of the Cascade arc, limited detection of DLP activity hinders classification of anomalous seismic behavior and its relationship to the LVC magmatic system. This study uses a template matching approach with seismic data from temporary nodal and permanent stations to detect and locate DLP earthquakes near the LVC between 2017 and 2024. Within the DLP catalog of 611 events, a transition occurs from scattered (2017–2020) to oscillatory (2021–2024) occurrence rates. During the oscillatory period, regional earthquakes with large amplitude velocity waveforms observed near the LVC are associated with abrupt changes in DLP occurrence rates. Continued monitoring of DLP activity has the potential to better define the processes that drive volcanic unrest at the LVC in the future.

California↗
USGS2024

New methodology for assessing underground natural gas storage resources – Example from Michigan Basin, United States

Energy consumption in the United States (U.S.) and across the world is shifting away from traditional fossil fuels like coal and oil, and towards natural gas and renewable sources, including hydrogen. Because gas demand is typically greatest during cold seasons and renewable sources sometimes produce variable supplies, it is important to store energy for use when demand exceeds supply. Whereas batteries and tanks typically store energy above ground, geologic (underground) storage may be able to retain much greater quantities of energy over much longer durations (e.g., Matos et al. 2019; Buursink et al. 2023). Consequently, the U.S. Geological Survey (USGS) is developing a methodology to locate new geologic energy storage sites with an initial focus on depleted hydrocarbon reservoirs followed by estimates of associated pore space or gas storage capacity. Depleted hydrocarbon fields in Michigan Basin in the U.S. already host large quantities of injected natural gas (about 681 billion cubic feet of seasonal gas capacity in 42 gas storage facilities; U.S. Energy Information Administration 2024). Therefore, we used this basin as an example to assess potential new storage resources alongside existing storage capacity and carbon dioxide sequestration resources (e.g., Katz and Coats 1968; Haagsma et al. 2020). In Michigan Basin gas storage is taking place in reservoirs within several key hydrocarbon producing formations with different lithologies, including the informal Stray sandstone of the Michigan Formation and the pinnacle reefs of the Niagara Group (Gautier et al. 1995; Swezey et al. 2015). The assessment methodology introduced here consists of two steps relying on both well information and reservoir data.

Michigan↗
USGS2025

Structural controls on splay fault rupture dynamics during Cascadia megathrust earthquakes

Great subduction earthquakes ( M w ≥ 8.0) can generate devastating tsunamis by rapidly displacing the seafloor and overlying water column. These potentially tsunamigenic seafloor offsets result from coseismic fault slip and deformation beneath or within the accretionary wedge. The mechanics of these shallow rupture phenomena and their dependence on subduction zone properties remain unresolved, partly due to the sparsity of offshore observations of shallow megathrust earthquake deformation. Here, we analyze how offshore structure influences shallow rupture mechanics and slip partitioning using 3D dynamic earthquake simulations of the Cascadia subduction zone (CSZ) megathrust with and without variably dipping seaward- or landward-vergent splay faults in the wedge that sole into the megathrust. Resulting tradeoffs between splay and megathrust slip reveal structural controls on rupture partitioning, with greater splay slip leading to less shallow megathrust slip updip. Gently dipping and seaward-vergent splays host more slip than those with steeper, landward-vergent splays. To isolate the underlying mechanisms, we compare models with Andersonian and plunging principal stresses. Results suggest distinct static and dynamic processes control the dip- and vergence-dependence of splay rupture: static (mis)alignment relative to far-field tectonic loading favors slip on more optimally oriented, shallowly dipping splay faults. In contrast, dynamic stress interactions of an updip-propagating megathrust rupture front with the free surface and potential branch faults favor forward branching onto seaward-vergent splays and inhibit backward branching onto landward-vergent splays. Resulting seafloor displacements suggest splay fault structure may influence coseismic tsunami source processes, highlighting the importance of dynamically viable rupture scenarios in subduction hazard assessments.

Cascadia subduction zone↗
USGS2025

Water-resources inventory and assessment at Katahdin Woods and Waters National Monument

The U.S. Geological Survey, in cooperation with the National Park Service, prepared a water-resources inventory and assessment for Katahdin Woods and Waters National Monument (KAWW). This compilation includes published and publicly accessible hydrologic data and resource assessments of streams, rivers, ponds, lakes, wetlands, vernal pools, and groundwater in and near KAWW. It also includes reports and datasets summarizing attributes of KAWW’s hydrologic infrastructure, such as stream crossings, dams, wastewater discharge plants, groundwater monitoring wells, and U.S. Geological Survey streamflow-gaging stations. Descriptions of data and details of current limitations in available datasets are included. Wetland, groundwater, streamflow, and water-quality information are all limited. Hydrography data are available; however, there are limited ground-truth data. Accurate streamlines within KAWW were developed from light detection and ranging (lidar) as a part of this work. Hydrologic infrastructure information is available from multiple sources; however, differences exist among the datasets. Datasets are summarized in appendix 1.

Maine↗
USGS2024

The 3D National Topography Model Call for Action—Part 1. The 3D Hydrography Program

The U.S. Geological Survey is initiating the 3D Hydrography Program (3DHP), the first systematic remapping of the Nation’s surface waters since the original 1:24,000-scale topographic mapping program was active from 1947 to 1992. Building on decades of experience maintaining the National Hydrography Dataset (NHD), the Watershed Boundary Dataset (WBD), and the NHDPlus High Resolution (NHDPlus HR), the 3DHP will completely refresh the Nation’s hydrography data and improve discovery and sharing of water-related data. The design of the 3DHP is based on the results of a study that estimated that the fully implemented program would have the potential to provide more than $1 billion in benefits to Federal, State, Tribal, Territorial, and local governments and to private and nonprofit organizations every year, in addition to myriad societal benefits. The 3DHP would directly support better decision making regarding water resources by providing more accurate, complete, and integrated information than is currently available. The 3DHP datasets will include a three-dimensional (3D) hydrography network generated from and integrated with elevation data from the 3D Elevation Program (3DEP) to better represent stream gradients and channel conditions, along with waterbodies, hydrologic units, hydrologically enhanced elevation and other surfaces, and more consistent and accurate attributes. The 3DHP datasets will inherit key attributes of the NHD, WBD, and NHDPlus HR, and they also will include new attributes and links to other data such as the U.S. Fish and Wildlife Service National Wetlands Inventory, groundwater data, and engineered hydrologic systems such as stormwater networks. The 3DHP will be designed to provide a set of open and interoperable web-based tools, maps, and data catalogs, creating a robust system for users to reference their information about water; the system elements are collectively referred to as the “infostructure.” The 3DHP and the infostructure can provide a foundational geospatial underpinning for the Internet of Water, a community-based effort to modernize tools and technologies to share water data. As proposed, the 3DHP would begin providing products and services to the public in 2024.

Circular↗
USGS2026

Quantifying methane emissions from a rich fen with uncrewed aircraft systems in boreal Alaska

Thawing of permafrost in northern latitudes is accelerating, potentially releasing substantial amounts of methane (CH 4 ) as forested permafrost plateaus transition into wetlands. This ecosystem shift alters the carbon exchange between the soil and atmosphere, influencing the permafrost-carbon feedback. Monitoring these changes may require measurement platforms operating across varied spatial and temporal scales. Recent advancements in small uncrewed aircraft systems (sUAS) enable high resolution CH 4 flux quantification in remote, complex terrains; however, comparisons with established methods such as eddy covariance flux towers remain limited. We used a hexacopter sUAS to quantify CH 4 emissions from the Alaska Peatland Experiment, a wetland within the Bonanza Creek Experimental Forest. Using an ensemble of methods to define the background CH 4 concentration, along with near surface emissions from soil chambers, helped constrain our flux estimates. The sUAS method yielded an average flux of 0.0077 ± 0.0019 mol s −1 CH 4 , within a factor of two concurrent tower-derived total source flux estimates (0.0036 ± 0.00042 mol s −1 CH 4 ). To assess spatial drivers of observed fluxes, we conducted a 2D footprint analysis and overlaid the results with high-resolution hyperspectral land cover classification, quantifying vegetative contributions within each footprint. This revealed higher fen representation in sUAS measurements (73.8%) than in tower footprints (58.8%), and lower tussock meadow representation (15.6% and 30.3%, respectively). These differences were consistent with known variation in vegetation-specific CH 4 emissions. Our results highlight that combining footprint modeling with land cover characterization can enhance interpretations of CH 4 fluxes and guide cross-platform comparisons.

JGR Atmospheres↗
USGS2026

Uranium—Deposits, production and resources, market dynamics, and supply chain risks

Introduction Interest in nuclear power for the generation of electricity has risen with the increase in the need for more diverse baseload power, enhanced energy security, and the development of new technologies, such as small modular reactors (SMRs), which could provide power for remote areas, industrial applications, and artificial intelligence (AI) data centers. In 2024, the U.S. Department of Energy received $2.7 billion in congressional funding to bolster the domestic uranium production and nuclear fuel supply chain and address reliance on imports from foreign suppliers. In 2025, the U.S. Government issued several Executive and Secretary’s orders aimed at revitalizing the U.S. nuclear sector. If SMRs are to be as widely deployed in the United States and worldwide as envisioned, demand for uranium (nuclear reactor fuel) will likely increase. After the Fukushima nuclear accident in 2011, the market spot price of uranium began a decline, followed by a decrease in U.S. and global uranium exploration and mine development expenditures that led to a uranium supply deficit until 2020, when prices started to recover, prompting a resurgence in uranium exploration and development. In January of 2024, the uranium spot price rose to a 17-year high $106 (U.S. dollars) per pound of U3O8 (triuranium oxide, commonly known as “yellowcake”), which is expected to increase uranium exploration, mine development, and uranium production domestically and worldwide.

Fact Sheet↗
USGS2024

Aboveground carbon stocks across a hydrological gradient: Ghost forests to non-tidal freshwater forested wetlands

Upper estuarine forested wetlands (UEFWs) play an important role in the sequestration of atmospheric carbon (C), which is facilitated by their position at the boundary of terrestrial and maritime environments but threatened by sea level rise. This study assessed the change in aboveground C stocks along the estuarine–riverine hydrogeomorphic gradient spanning salt-impacted freshwater tidal forested wetlands to freshwater forested wetlands in seasonally tidal and nontidal landscape positions. Standing stocks of C in forested wetlands were measured along two major coastal river systems, the Winyah Bay in South Carolina and the Savannah River in Georgia (USA), replicating and expanding a previous study to allow the assessment of change over time. Aboveground C stocks on these systems averaged 172.9 Mg C ha −1 , comparable to those found in UEFWs across the globe and distinct from the terrestrial forested ecosystems they are often considered to be a part of during large-scale C inventory efforts. Groundwater salinity conditions as low as 1.3 ppt were observed in conjunction with losses of aboveground C. When viewed in context alongside expected sea level rise and corresponding saltwater intrusion estimates, these data suggest a marked decrease in aboveground C stocks in forested wetlands situated in and around tidal estuaries.

Georgia, South Carolina↗
USGS2024

Earthquake-triggered ground-failure inventory associated with the M7.1 2018 Southcentral Alaska earthquake

The 30 November 2018, magnitude (Mw) 7.1 earthquake in Southcentral Alaska triggered substantial landslides, liquefaction, and ground cracking throughout the region, resulting in widespread geotechnical damage to buildings and infrastructure. Despite a challenging reconnaissance and remote-sensing environment, we constructed a detailed digital inventory of ground failure associated with the event from several sources. Sources included information derived from remotely sensed data, and data compiled from literature, social media postings, and earthquake damage information compiled by local, state, and federal agencies. Each instance of ground failure within the inventory contains information on the location and type of observed ground failure, and the methods and data used to document the occurrence. Where high-quality data, such as LIDAR or satellite imagery, were available and showed the ground-failure instance clearly, the extent is mapped as a polygon or polyline. All other locations are mapped as points. There are a total of 886 ground-failure instances documented within the inventory (400 landslides, 286 liquefaction features, and 200 features unattributed to specific processes). A semi-quantitative confidence scheme is used to describe mapping certainty associated with each ground-failure feature. This inventory represents a relatively moderate ground-failure-triggering event that occurred in a subarctic environment. This data paper describes the content within the inventory, the inventory data collection procedures, and limitations of the data. Events of this type are not often documented in detail; thus, adding the inventory data to the US Geological Survey Open Repository of Earthquake-Triggered Ground-Failure Inventories further diversifies the datasets available to the scientific community to be used to better understand and model earthquake-triggered ground failure.

Alaska↗
USGS2026

Evidence for organic and igneous carbon sources in the Mesoproterozoic Adirondack Highlands graphite district, New York

This study examined the origin of graphite mineralization in the Mesoproterozoic Adirondack Highlands. Flake graphite was mined from the historic mining district around Ticonderoga, New York, USA, from the 1850s to the 1920s. Most of the production and reserves from the district are disseminated graphite from feldspathic quartz schists (e.g., Dixon-American Graphite Co. mine, Faxon property, Flake Graphite Company). Protolith sediments of the schists were deposited in the Trans-Adirondack backarc basin at ca. 1.25 Ga and were metamorphosed to high grade at 1.18–1.15 Ga during the accretionary Shawinigan orogeny and at 1.09–1.02 Ga during the Ottawan phase of the Grenvillian orogeny. Major elements and carbon isotopes of graphite from schists (δ 13 C = −28.1‰ to −5.4‰ relative to Vienna Peedee belemnite) are consistent with derivation from organic carbon in protolith sandstones that devolatilized during metamorphism. Some small mines and prospects in the district mined graphitic calcite marbles and clinopyroxene ± scapolite skarns associated with a pegmatite suite dated by U-Pb in zircon to 1.04–1.01 Ga (e.g., Lead Hill and Crown Point deposits), and rare graphite veins are observed at several deposits. Hydrothermal flake graphite from pegmatites, skarns, and veins are distinct from graphite in schist-hosted deposits and have δ 13 C values that cluster around −7‰ to −5‰, consistent with deposition from igneous CO 2 . Covariation of carbon and oxygen stable isotopes and U-Pb zircon geochronology are not compatible with precipitation from fluids derived from country rock metapelites or marbles. Some graphite districts elsewhere show compelling evidence for remobilization and redeposition of sedimentary-derived carbon as hydrothermal graphite, sometimes as thick veins producing world-class, high-grade deposits. In the Adirondacks, there is evidence for two generations of graphite mineralization with distinct carbon sources: syngenetic graphite with an organic origin in Shawinigan quartz schists and epigenetic graphite with an igneous origin in Ottawan skarns and pegmatites.

New York↗
USGS2025

Reply to, “Comment on ‘The 1886 Charleston, South Carolina, earthquake: Relic railroad offset reveals rupture,’ by Roger Bilham and Susan E. Hough”

We welcome this opportunity to respond to Pratt et al. (2024) (hereinafter P24). Bilham and Hough (2023) proposed a “first-cut” elastic deformation model for the 1886 earthquake, a quantitative source model constrained by identified coseismic constraints. A key observation was the measurement of a lateral offset of a railroad line south of Summerville, leading to a model with predominately dextral slip and minor convergence, from which we concluded that active faulting had raised the Penholoway Marine Terrace >6 m since ∼770 ka. P24 questioned these constraints and proposed an alternative rupture model with predominantly reverse slip. This alternative model is neither consistent with coseismic constraints nor with other geophysical data. In a revised model presented here, we recognize that uplift of the Penholoway Terrace is confined to the eastern edge of the terrace, which we conclude results from active folding and tectonic transpression centered on the dextral fault that offset the railroad in 1886.

South Carolina↗
USGS2024

Automated, near real-time ground-motion processing at the U.S. Geological Survey

We describe automated ground‐motion processing software named gmprocess that has been developed at the U.S. Geological Survey (USGS) in support of near‐real‐time earthquake hazard products. Because of the open‐source development process, this software has benefitted from the involvement and contributions of a broad community and has been used for a wider range of applications than was initially envisioned. Here, we give an overview and introduction to the software, including how it has leveraged other open‐source libraries. We highlight some key features that gmprocess provides, compare response spectra calculated with the automated processing approach of gmprocess to the response spectra provided by the Next Generation Attenuation projects, and summarize projects that have utilized gmprocess. These use‐cases demonstrate that this software development effort has been successfully leveraged in earthquake research activities both within and outside the USGS.

Seismological Research Letters↗
USGS2025

Bølling-Allerød productivity in the subarctic Pacific driven by seasonal upwelling

The Bølling-Allerød deglacial event is marked by high diatom productivity and opal deposition throughout the subarctic Pacific. This opal could either constitute a strengthened biological pump and thus carbon sequestration, or a weakened biological pump and release of marine-sequestered CO 2 to the atmosphere. We quantify silicic acid supply at IODP Site U1340 in the Bering Sea using biogenic opal and δ 30 Si of Coscinodiscus , a diatom genus. These records, along with diatom environmental indicators, suggest the Bølling-Allerød had high silicic acid availability related to a shift from stratification to seasonal upwelling dynamics. We thus propose the primary cause of the high productivity event was increased macronutrient supply from vertical exchange that injected old, nutrient-rich, CO 2 -rich waters into the surface. Enhanced CO 2 release from the subarctic Pacific may help explain critical intervals of CO 2 rise that occur at the onsets of the Bølling-Allerød and PreBoreal.

Geophysical Research Letters↗
USGS2025

Seismic intensities, ground motions and basin amplification from the 2022 Mw 6.4 Ferndale, California intraslab earthquake

On December 20, 2022, a Mw 6.4 earthquake occurred at a depth of 18 km within the subducting Gorda plate in the Mendocino Triple Junction (MTJ), one of the most seismically active regions of the contiguous United States, causing widespread damage to local communities. Here we document the seismic intensities, ground motions, and basin amplification effects recorded by this earthquake across Humboldt County as part of an ongoing scientific effort to understand subduction zone earthquake hazards. Modified Mercalli Intensity (MMI) values from our post-earthquake field survey report shaking intensities as high as VIII (Severe). Strong ground motion data from 54 seismic stations were processed to calculate amplitude and frequency content parameters across Northern California. The maximum calculated geometric mean PGA and PGV are 1053 cm/s 2 and 52 cm/s (respectively), both recorded within the 3 km deep Eel River sedimentary basin. Comparisons with four published Ground Motion Prediction Equations indicate that PGA and PGV measurements align with expected attenuation-distance patterns for intraslab earthquakes of this nature. Within the Eel River Basin, ground motions for frequencies above 1 Hz are amplified, with respect to reference stations, by factors greater than 2. Our findings suggest that peak ground motions were mainly caused by sedimentary basin site-effects within the Eel River basin, although rupture directivity may have also increased ground motion amplitudes locally. While the Ferndale area is primarily impacted by shallow intraslab earthquakes, our results also raise questions about site-response and basin amplification hazards from a potential megathrust earthquake. The high seismicity rates of the southernmost Cascades call for stronger regional preparedness and improved strategies to mitigate the effects of such a large-scale disaster.

California↗