Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Great Lakes Research”

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,663 records · Page 3Linked to original sources
USGS2025

Refining the earthquake history of south-central Alaska through lake records

The Alaska–Aleutian subduction zone (AASZ) is one of the world's most seismically active plate boundaries and the source of the 1964 M w 9.2 Great Alaska earthquake–the second largest instrumentally recorded earthquake in the world. Understanding the nature and frequency of such earthquakes is necessary for seismic and tsunami hazard assessment, but instrumental and historical records that span less than 150 years are too short to allow a statistically reliable analysis of earthquake recurrence times. This calls for studies of evidence of past earthquakes, extending the earthquake catalog further back in time. Subduction-zone paleoseismology in south-central Alaska is predominantly based on coastal evidence of land-level changes and tsunamis generated by megathrust earthquakes and preserved in the geological record. A complementary approach is lacustrine paleoseismology, which is still a relatively young discipline in Alaska. However, globally, lake basins are well-established high-resolution and continuous recorders of paleoseismic activity along subduction zones, relying on the identification of underwater landslide deposits and turbidites generated by seismic shaking. As a result, lake basins not only register ground shaking from megathrust earthquakes, but also from intraslab and crustal earthquakes, which are typically not accompanied by significant land-level changes. In this review paper, we combine coastal and lacustrine paleoseismology approaches to refine the south-central Alaskan earthquake history by comparing the paleoseismic records from two lakes (i.e., Eklutna Lake, located in the Chugach Mountain Range, and Skilak Lake, situated on the Kenai Peninsula) with the coastal and crustal earthquake catalog in Alaska. The resulting age ranges of all known megathrust earthquakes involving the Alaskan megathrust between the Kodiak and Prince William Sound (PWS) sections are more precise and accurate for the last 1.3 kyrs BP than the previously published age ranges from coastal records. As a result, this study supports the following key conclusions: (1) The 1964 CE earthquake was an exceptionally strong and unique event in the last 2000 years, rupturing the PWS, Kenai, Barren Islands, and Kodiak sections simultaneously. (2) The high-resolution and seasonal markings of the varved lake records now disentangle for the first time closely timed earthquakes, which was not possible based on the coastal evidence alone. (3) No persistent megathrust rupture boundaries exist. So, the possibility of a full rupture of the entire eastern AASZ, from PWS to Semidi cannot be excluded. (4) The rupture pattern in the eastern AASZ reveals superimposed cycles of multi-asperity ruptures (1964 earthquake) and clustered complementary partial ruptures, or rupture cascades. (5) The PWS section hosts the largest asperity in the eastern AASZ. (6) The shaking record of megathrust earthquakes indicates a time-dependent (quasiperiodic) behavior for the study area, but the observation of complementary clusters means that the hazards will not drop to zero but instead may even increase for a neighboring section. (7) The time-independent behavior of intraplate earthquakes implies that the intraslab hazard did not decrease following the 2016 and 2018 earthquakes. This study utilizes an integrated approach for subduction zone paleoseismology as a solution for unraveling recurrence and rupture patterns in Alaska, which can be applied worldwide.

Alaska↗
USGS2025

Bathymetric and acoustic-backscatter mapping of Lake Sammamish, Washington, during USGS field activity 2021-656-FA

In 2021, scientists from the U.S. Geological Survey (USGS) acquired bathymetric and acoustic-backscatter data within Lake Sammamish, Washington. Mapping was completed to find evidence of past earthquakes such as underwater landslides and is part of a larger USGS project to understand the overall geologic hazards history of the Cascadia Margin region. The survey was conducted using the USGS research vessel Parke Snavely outfitted with an interferometric sidescan-sonar system for swath mapping. Three estimates of bathymetric uncertainty were calculated from the post-processed bathymetric data. This report provides a summary of the mapping mission.

Washington↗
USGS2025

Multi-scale predictors of Northern Long-eared Bat (Myotis septentrionalis) occupancy in the United States

Historically, Myotis septentrionalis (Northern Long eared Bat) was among the most common forest-interior species in North America. Largely due to high mortality from white-nose syndrome, this species has experienced severe population declines across its range. To create an updated species distribution map representing summer occupancy probabilities from 2017 to 2022, we integrated stationary acoustic data with live-capture data from the database of the North American Bat Monitoring Program into a multi-scale, multi-method occupancy modeling framework. Our results provide data-driven predictions with quantified uncertainty for summer occupancy probabilities for Northern Long-eared Bats at 2 spatial scales across the range of the species, while also accounting for inherent observation biases (e.g., imperfect detection).

Journal of North American Bat Research↗
USGS1976

Solute transport and modeling of water quality in a small stream

An injection of chloride, sodium, and stable strontium was made at a constant rate for 3 hours into Uvas Creek, Santa Clara County, Calif., to determine the mass transport processes in a small stream. Five observation points were selected within a 610-metre reach of the stream below the injection site. Water samples were collected at the observation points during and immediately after the injection. A mathematical model of the stream was obtained by solving analytically and optimally the one-dimensional mass transport equation of the solutes in the stream. Comparison of field results with a simplified mathematical model indicates the dominance of convection in the behavior of sodium and chloride. The concentration of chloride and sodium can be closely simulated by the model. However, strontium cannot be well represented by the simplified model, which contains a first-order decay-type sink.

California↗
USGS2026

Synthesizing a twelve-year sediment trap time series of planktic foraminiferal flux in the Gulf of America (Mexico)

Sediment trap time series provide powerful frameworks for testing hypotheses about planktic foraminiferal assemblage composition, seasonality, and geochemical responses to environmental variability, all of which are central to improving paleoceanographic reconstructions. We present high resolution foraminiferal assemblage data from a long-running sediment trap (2008–2020) in the northern Gulf of America (Mexico). This study summarizes the species composition, seasonality, and size distribution of the fifteen most abundant species of planktic foraminifera, which account for 98% of total flux. Foraminiferal flux peaks in winter and reaches a minimum in summer, following the seasonal pattern of primary production in the northern Gulf. Winter assemblages are dominated by non-spinose taxa, whereas spinose taxa prevail during summer. Across nearly all species, average monthly test size covaries with temperature, independent of seasonal flux trends. Notably, Trilobatus sacculifer and Neogloboquadrina dutertrei show a significant decline in relative abundance and flux during 2017–2020, nearly disappearing from the assemblage.

Journal of Foraminiferal Research↗
USGS2026

Global performance of remote sensing-based and reanalysis-driven models to estimate open water evaporation

Evaporation plays an essential role in the water cycle, influencing local and regional climates while directly impacting water availability in lakes. However, directly measuring evaporation over water bodies remains challenging due to the high costs of installing and maintaining the required in situ instrumentation. Although several remote sensing algorithms have been providing evaporation estimates, the lack of a global validation hinders our understanding of their relative uncertainties and performances across different regions. Here, we analyze the performance of a suite of models that leverage satellite data and meteorological reanalysis to estimate evaporation over lakes worldwide. We compare 3 remote sensing-based models, 1 reanalysis-driven model and 1 ensemble approach, using in situ observations from 27 lakes representing a diverse range of geographic and climatic regions. Our results demonstrate that, overall, the ensemble outperformed any individual model in terms of accuracy, with a RMSE and a bias of 1.3 and 0.3 mm day −1 , respectively. These findings highlight the benefits of using an ensemble approach to estimate open water evaporation with satellite-based models at the global scale, leveraging the unique strengths of each model. For the individual models, differences in the representation of heat storage changes and advection effects led to lower values of RMSE and bias, depending on the location and depth of the lakes. This study sets the path for future improvement of open water evaporation algorithms globally, while remote sensing techniques are proven satisfactory to monitoring of water loss in lakes globally, an essential step toward effective large-scale water resources management.

Water Resources Research↗
USGS2025

Colored shaded-relief bathymetry and acoustic backscatter of Lake Sammamish, Washington

Evidence of strong earthquakes (such as underwater landslides and associated deposits) may be recorded within the lacustrine sediments of Pacific Northwest lakes. The floor of Lake Sammamish, Wash., an approximately 11 kilometer (6.8 mile) long, 2 kilometer (1.2 mile) wide, and 35 meter (114.8 feet) deep lake located in a populated region just east of Seattle, was mapped by the U.S. Geological Survey in November of 2021 to search for evidence of past earthquakes. Mapping was conducted using a SWATHplus-M 234-kHz interferometric side-scan sonar system was pole-mounted on the U.S. Geological Survey research vessel Parke Snavely, and the system collected full-coverage bathymetric and acoustic backscatter data, which were processed to 2-meter spatial resolution. Two maps were created, a colored shaded-relief bathymetric map showing lake floor morphology (sheet 1), and an acoustic-backscatter map showing backscatter intensities (sheet 2). The results may then be utilized together to investigate past earthquake activity.

Washington↗
USGS2026

Decadal-scale characteristics of natural and anthropogenic dune morphology along North Carolina Barrier Islands (SE Atlantic Coast)

Dunes offer protection to coastal communities from storms, yet quantitative research characterizing human influences on dune systems across decades and 10–100 km spatial scales constitutes a known data gap. This study analyzed a multidecadal archive of dune crest elevations and positions to map mesoscale spatiotemporal trends in foredune morphology across the barrier island coastline of North Carolina, SE Atlantic coast. Results are summarized within four categories of coastline representing a range of coastal management intensities, integrating beach width and dune lateral accommodation space to assess the comparative influence of differing intensities of anthropogenic intervention on dune morphology. These categories account for the presence of structures, roads, and alongshore connectivity to anthropogenic sediment sources (beach nourishment). These categories show distinct dune and beach system characteristics. Dunes are largest in the areas with the longest legacy of coastal management. As categories decrease in development intensity, dune and beach system characteristics generally depart from the large, static dunes and narrow beaches consistent with hold-the-line management. Where small or transgressive dunes exist on heavily developed coasts, they coincide with relatively strong hydrodynamic forcing or sediment-limiting geology. There is evidence of indirect coastal management influence on undeveloped barrier islands with alongshore connectivity to developed barriers. While this analysis was conducted at the mesoscale, the results provide researchers with the framework to make insights at smaller scales than those demonstrated in this study.

North Carolina↗
USGS2026

Evaluation of best management practices at an edge-of-field site in the Eagle Creek watershed, Ohio, 2012–20

In 2010, the U.S. Geological Survey worked in partnership with the Great Lakes Restoration Initiative and the Natural Resources Conservation Service to identify farm fields in priority watersheds—watersheds critically important to Great Lake health. Intensive best management practices were implemented in five States to test the efficacy of best management practices, also referred to as agricultural conservation practices, to reduce sediment and nutrient runoff on agricultural fields. These fields were chosen as representatives of the priority watersheds because their farming practices and geographical conditions were common among farms in those watersheds. A farm site located near Findlay, Ohio, in the Eagle Creek watershed used two conservation practices, cover crop (Natural Resources Conservation Service practice standard 340) and variable rate technology (VRT; Natural Resources Conservation Service practice standard 590), which represent common practices in the Maumee priority watershed. This study monitored surface runoff and subsurface tile runoff at the Eagle Creek watershed site. The effects of cover crop and VRT nutrient application on sediment and nutrient runoff were assessed from October 2012 to September 2020. Cover crops were applied in fall of 2016, 2017, and 2018, coinciding with VRT applications. The following parameters were analyzed as part of the study: runoff discharge; peak discharge; and concentrations and loads of total phosphorus, particulate phosphorus, dissolved reactive phosphorus, and suspended sediment. High peak runoff events disproportionally affected nutrient concentrations and loads in runoff. A threshold of the 85th percentile of the peak discharge was chosen to assess the effect of cover crop and VRT nutrient application on nutrient mitigation at each gage. Runoff events below the 85th percentile were considered small, and runoff events over the 85th percentile were considered large. Concentrations from collected water-quality samples were used for the analysis; estimated concentrations were excluded from the statistical evaluation. All surface runoff parameters during small runoff events, except surface runoff volume, were significantly lower in the period with cover crops and VRT implementation. Contrarily, subsurface runoff parameters during small events were not statistically different between the periods with and without cover crops and VRT implementation. Surface and subsurface parameters during large events were also not statistically different between the periods with and without cover crops and VRT implementation. These results indicate cover crops and VRT may improve water quality during small runoff events, but additional best management practices that mitigate large runoff events may lead to greater water-quality improvement given the contribution of large events to overall losses.

Ohio↗
USGS2024

Exploring the dynamic interactions between the Southern San Andreas Fault and a normal fault under the Salton Sea

We investigate the dynamic interactions between the Southern San Andreas Fault (SSAF) and a proximal normal fault (NF) beneath the Salton Sea in southern California. The NF, positioned near the SSAF terminus at Bombay Beach, exhibits 11–15 displacement events across 14 stratigraphic sequences, with a range of 0.2–1.4 m of vertical offset since ∼2–3 ka. Notably, four of these events may align temporally with SSAF earthquakes, raising questions about the possible interplay between the two faults. Utilizing dynamic rupture models, we analyze the coseismic interactions between the SSAF and NF, addressing under what conditions the SSAF induces slip on the NF. Our findings reveal that a suite of SSAF ruptures, particularly those propagating from north to south, can trigger slip on the normal fault and replicate observed vertical offsets. If the SSAF extends beneath the Salton Sea, earthquakes originating south of the NF intersection are less likely to trigger normal fault slip, although we cannot exclude this possibility. Some SSAF ruptures do not trigger discernible slip on the NF, rendering such events undetectable in the stratigraphic record. Our research contributes toward discussions regarding the seismic hazard in southern California, shedding light on the interplay between the SSAF and NF.

California↗
USGS2026

Identifying overlap between native fish movements and a possible seasonal Grass Carp deterrent in the Sandusky River, Ohio, USA

Riverine fishes are vulnerable to habitat fragmentation caused by interrupted access to vital habitats. Fragmentation using artificial structures can be used to limit the spread of invasive fishes by inhibiting movement between critical habitats. Thus, tradeoffs exist between maintaining connectivity for native species and restricting movement of invasive fishes. A nonphysical deterrent has been proposed to prevent invasive Grass Carp ( Ctenopharyngodon idella ) from reaching spawning habitats in the Sandusky River, a tributary to Lake Erie, during late spring–summer. However, this timing overlaps with spawning periods of a suite of native fishes. We used acoustic telemetry and ichthyoplankton surveys to examine movement and spawning activity of multiple native species and to assess potential impacts of deterrent operation on reproduction. All species used the proposed deterrent location during the Grass Carp spawning period, and their movements corresponded with larval fish collections, indicating upstream spawning. These findings suggest the deterrent could unintentionally restrict native species’ access to essential habitats. Balancing habitat connectivity and invasive species control will likely depend on species-specific responses and careful timing of deterrent operation.

Ohio↗
USGS2024

Earthquake cycle mechanics during caldera collapse: Simulating the 2018 Kīlauea eruption

In multiple observed caldera-forming eruptions, the rock overlying a draining magma reservoir dropped downward along ring faults in sequences of discrete collapse earthquakes. These sequences are analogous to tectonic earthquake cycles and provide opportunities to examine fault mechanics and collapse eruption dynamics over multiple events. Collapse earthquake cycles have been studied with zero-dimensional slider-block models, but these do not account for the complicated interplay between fluid and elastic dynamics or for factors such as the heterogeneous fault properties and non-vertical ring fault geometries often inferred at volcanoes. We present two-dimensional axisymmetric mafic piston-like collapse earthquake cycle models that include rate-and-state friction, fully-dynamic elasticity, and compressible viscous fluid magma flow. We demonstrate that collapse earthquake intervals and magnitudes are highly sensitive to inertial effects, evolving stress fields, fault geometry, and depth-varying fault friction. Given the consistent earthquake cycles observed in most eruptions, this suggests that ring faults can quickly stabilize and often become nearly vertical at depth. We use the well-monitored 2018 collapse sequence at Kı̄lauea as a case study. Our model can produce many features of Kı̄lauea seismic and geodetic observations, except for a significant amount of interseismic slip, which cannot be readily explained with simple rate-and-state friction parameterizations.

Hawaii↗
USGS2025

3D viscoelastic models of slip-deficit rate along the Cascadia subduction zone

Interseismic deformation in the Pacific Northwest is constrained by the horizontal crustal velocity field derived from the Global Positioning System (GPS) in addition to vertical rates derived from GPS, leveling, and tide gauge measurements. Such measurements were folded in to deformation models of fault slip rates as part of the 2023 National Seismic Hazard Model (NSHM) update. Here I build upon one of the contributing models, the viscoelastic earthquake-cycle model of Pollitz [2022]. This model permits inclusion of effects of time-dependent viscoelastic relaxation within earthquake cycles (i.e., ‘ghost transients’) and laterally variable elastic and/or ductile material properties. I lever-age these capabilities to incorporate the Cascadia megathrust into Western U.S.-wide deformation models in which crustal fault slip rates are estimated simultaneously with slip deficit rates along the interplate boundary between the descending Juan de Fuca plate and North American plate. This effort includes construction of a margin-wide model of viscoelastic structure founded on the Slab 2.0 model and probes different models of the ductile properties of the surrounding oceanic asthenosphere, continental lower crust, and mantle asthenosphere. This results in new estimates of the distribution of slip deficit rate along the ∼ 1000 km long margin, highlights the importance of correcting for glacial-isostatic adjustment effects, and permits assessment of sensitivity of results to assumed ductile properties.

California, Oregon, Washington↗
USGS2000

Storm impact scale for barrier islands

A new scale is proposed that categorizes impacts to natural barrier islands resulting from tropical and extra-tropical storms. The proposed scale is fundamentally different than existing storm-related scales in that the coupling between forcing processes and the geometry of the coast is explicitly included. Four regimes, representing different levels of impact, are defined. Within each regime, patterns and relative magnitudes of net erosion and accretion are argued to be unique. The borders between regimes represent thresholds defining where processes and magnitudes of impacts change dramatically. Impact level 1 is the 'swash' regime describing a storm where runup is confined to the foreshore. The foreshore typically erodes during the storm and recovers following the storm; hence, there is no net change. Impact level 2 is the 'collision' regime describing a storm where the wave runup exceeds the threshold of the base of the foredune ridge. Swash impacts the dune forcing net erosion. Impact level 3 is the 'overwash' regime describing a storm where wave runup overtops the berm or, if present, the foredune ridge. The associated net landward sand transport contributes to net migration of the barrier landward. Impact level 4 is the 'inundation' regime describing a storm where the storm surge is sufficient to completely and continuously submerge the barrier island. Sand undergoes net landward transport over the barrier island; limited evidence suggests the quantities and distance of transport are much greater than what occurs during the 'overwash' regime.

Journal of Coastal Research↗
USGS2024

Afterslip and creep in the rate-dependent framework: Joint inversion of borehole strain and GNSS displacements for the Mw 7.1 Ridgecrest earthquake

The elusive transition toward afterslip following an earthquake is challenging to capture with typical data resolution limits. A dense geodetic network recorded the Mw 7.1 Ridgecrest earthquake, including 16 Global Navigation Satellite System (GNSS) stations and 3 borehole strainmeters (BSM). The sub-nanostrain precision and sub-second sampling rate of BSMs bridges a gap between conventional seismologic and geodetic methods, exemplified by atypical postseismic shear strain reversals observed at nearfield (<2 km) station B921 that remain unexplained. We jointly invert GNSS displacements and BSM strains for coseismic and postseismic slip spanning hours to months over 7 independent periods. Cosiesmically, our model resolves the largest slip magnitudes of up to 6.6 m on the mainshock rupture plane, with similar patterns to other inferred slip distributions. The foreshock fault appears to slip coincidently with mainshock, revealing potential asperities activated during the preceding Mw 6.4 event. Postseismically, the best-fitting models adhere to mechanical rate-and-state expectations of logarithmically decaying slip adjacent to the coseismic rupture terminus, and where deep rheologic conditions favor creep. Most spatial variation occurs in the early postseismic timeframe (<1–2 weeks), with evidence for regional rheologic control and static stress dependence. Triggered creep on the neighboring Garlock Fault unexpectedly persists for >178 days—further highlighting the importance of fault networks in postseismic stress redistribution, critical to assessing future hazard.

Journal of Geophysics Research↗
USGS2024

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

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

Journal of Geophysical Research: Solid Earth↗
USGS2024

RegionGrow3D: A deterministic analysis for characterizing discrete three-dimensional landslide source areas on a regional scale

Regional-scale characterization of shallow landslide hazards is important for reducing their destructive impact on society. These hazards are commonly characterized by (a) their location and likelihood using susceptibility maps, (b) landslide size and frequency using geomorphic scaling laws, and (c) the magnitude of disturbance required to cause landslides using initiation thresholds. Typically, this is accomplished through the use of inventories documenting the locations and triggering conditions of previous landslides. In the absence of comprehensive landslide inventories, physics-based slope stability models can be used to estimate landslide initiation potential and provide plausible distributions of landslide characteristics for a range of environmental and forcing conditions. However, these models are sometimes limited in their ability to capture key mechanisms tied to discrete three-dimensional (3D) landslide mechanics while possessing the computational efficiency required for broad-scale application. In this study, the RegionGrow3D (RG3D) model is developed to broadly simulate the area, volume, and location of landslides on a regional scale (≥1,000 km 2 ) using 3D, limit-equilibrium (LE)-based slope stability modeling. Furthermore, RG3D is incorporated into a susceptibility framework that quantifies landsliding uncertainty using a distribution of soil shear strengths and their associated probabilities, back-calculated from inventoried landslides using 3D LE-based landslide forensics. This framework is used to evaluate the influence of uncertainty tied to shear strength, rainfall scenarios, and antecedent soil moisture on potential landsliding and rainfall thresholds over a large region of the Oregon Coast Range, USA.

Journal of Geophysical Research: Earth Surface↗
USGS2024

Evidence for low effective stress within the crust of the subducted Gorda plate from the 2022 December Mw 6.4 Ferndale earthquake sequence

Stress levels on and adjacent to megathrust faults at seismogenic depths remain a key but difficult to constrain parameter for assessing seismic hazard in subduction zones. Although strong ground motions have been observed to be generated from distinct, high-stress regions on the downdip end of the megathrust rupture areas in many great earthquakes, we lack direct constraints on the stress level in the lower seismogenic portion of the Cascadia megathrust. On 2022 December 20, a M w 6.4 strike-slip earthquake occurred near Ferndale, California in southern Cascadia and likely ruptured the Gorda slab crust in the lower seismogenic portion, providing an opportunity to assess the stress level in this region. Here, we relocate the Ferndale mainshock and the first two weeks of aftershocks using a high-resolution 3-D velocity model and estimate rupture dimensions, directivity, and stress drop for several M w 4-5 aftershocks and recent earthquakes. The aftershocks define a strike-slip fault in the slab crust striking ENE, consistent with the mainshock focal mechanism. The orientation of this fault is about 45 ° off the ideally oriented fault plane given the stress state in the slab. The aftershock zone is extensive and broad in the forward direction of the mainshock rupture but still constrained within the volume of high Vp/Vs within the slab crust. Our stress drop estimates are generally lower for M w 4-5 earthquakes located within the slab crust compared to those a few km deeper in the slab mantle. Combined, our results support a relatively low effective stress level in the vicinity of the megathrust in the lower portion of the seismogenic zone in southern Cascadia, likely due to elevated fluid pressures. Consequently, the ground motion in the onshore region above this low-stress seismogenic portion in southern Cascadia may not be as intense as that observed during great earthquakes in other subduction zones.

California↗