USGS Science⌕ Search

SEARCH · USGS Science

Results for “earthquake hazards”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,675 records · Page 93Linked to original sources

Displaying seismic deaggregation: The importance of the various sources

Seismic hazard deaggregation has become a standard part of probabilistic seismic hazard assessment (PSHA). The first product of PSHA is calculation of the likely severity of ground motion at a given range of annual probability levels, and this is extremely important for seismic design of structures to be built at the site under examination. However, for full analysis of proposed structural designs, engineers also need to examine scenario events to produce detailed time histories. To select such scenarios, a deaggregation of the hazard is performed, whereby the details of sources that contribute to the annual frequency of exceeding specified levels of ground motion, or P exc , are identified. A common format for such a deaggregation is shown in Figure 1 . This relates to the 475-year peak ground acceleration (pga) at Wellington, New Zealand (41.28°S 174.77°E), and shows the distribution in magnitude and distance of sources that contribute to P exc . Return period is approximately the reciprocal of P exc . Stiff soil site conditions ( Standards New Zealand 2004 ) were assumed. The analysis in Figure 1 used the interim version of the updated seismic hazard model for New Zealand ( Stirling et al. 2007 ), with the attenuation function developed by McVerry et al. (2007). Based on a Poisson time dependence model, a return period of 475 years corresponds to a 10% probability of exceedance in 50 years. From Figure 1 , it is apparent that for this site the main contribution to ground motion of this severity is from earthquakes of magnitude about 7.6 less than 10 km from the site (blue), and there is another strong contribution from larger events in the distance range 10 to 20 km (red). These correspond to the Wellington and Wairarapa faults, respectively (see Table 1 ). There are other events less than 10 km from the site and small contributions from other sources. At this site the major contributions are from specific faults nearby, which are readily identified. At sites where there is significant background seismicity, however, the plot will be much more complicated and not so easy to interpret. Figure 1 deaggregates probabilistic pga at the site; other parameters are also commonly deaggregated in the same way, in particular response spectral acceleration at a variety of natural periods. But the figure has a major shortcoming in that it represents only one return period; to obtain a full appreciation of the various contributing sources it is necessary to perform a succession of analyses to cover the full range of return periods.

Seismological Research Letters↗

Sediment thickness and ground motion site amplification along the United States Atlantic and Gulf Coastal Plains

Past and present research on earthquake ground motions along the Atlantic and Gulf Coastal Plains and Mississippi Embayment show significant period-dependent site response that is not presently accounted for in ground motion models. These deviations are strongly correlated with the thickness of Mesozoic and younger syn- and post-rift sediments. With the recent incorporation of deep basin depth measurements in the U.S. Geological Survey National Seismic Hazard Model for select regions in the western United States, we move toward a similar analysis in the greater Coastal Plains region by constructing a sediment thickness model and considering three new site response models conditioned on sediment thickness. As of the preparation of this conference paper, we have performed a preliminary evaluation of the Chapman and Guo model and find that the predicted ratio between pseudo-spectral accelerations for the Coastal Plains relative to the continental interior are broadly consistent with our independent dataset.

Atlantic and Gulf Coastal Plains↗

Preliminary analysis of the earthquake (MW 8.1) and tsunami of April 1, 2007, in the Solomon Islands, southwestern Pacific Ocean

On April 1, 2007, a destructive earthquake (Mw 8.1) and tsunami struck the central Solomon Islands arc in the southwestern Pacific Ocean. The earthquake had a thrust-fault focal mechanism and occurred at shallow depth (between 15 km and 25 km) beneath the island arc. The combined effects of the earthquake and tsunami caused dozens of fatalities and thousands remain without shelter. We present a preliminary analysis of the Mw-8.1 earthquake and resulting tsunami. Multichannel seismic-reflection data collected during 1984 show the geologic structure of the arc's frontal prism within the earthquake's rupture zone. Modeling tsunami-wave propagation indicates that some of the islands are so close to the earthquake epicenter that they were hard hit by tsunami waves as soon as 5 min. after shaking began, allowing people scant time to react.

Science of Tsunami Hazards↗

GEM Building Taxonomy (Version 2.0)

This report documents the development and applications of the Building Taxonomy for the Global Earthquake Model (GEM). The purpose of the GEM Building Taxonomy is to describe and classify buildings in a uniform manner as a key step towards assessing their seismic risk, Criteria for development of the GEM Building Taxonomy were that the Taxonomy be relevant to seismic performance of different construction types; be comprehensive yet simple; be collapsible; adhere to principles that are familiar to the range of users; and ultimately be extensible to non-buildings and other hazards. The taxonomy was developed in conjunction with other GEM researchers and builds on the knowledge base from other taxonomies, including the EERI and IAEE World Housing Encyclopedia, PAGER-STR, and HAZUS. The taxonomy is organized as a series of expandable tables, which contain information pertaining to various building attributes. Each attribute describes a specific characteristic of an individual building or a class of buildings that could potentially affect their seismic performance. The following 13 attributes have been included in the GEM Building Taxonomy Version 2.0 (v2.0): 1.) direction, 2.)material of the lateral load-resisting system, 3.) lateral load-resisting system, 4.) height, 5.) date of construction of retrofit, 6.) occupancy, 7.) building position within a block, 8.) shape of the building plan, 9.) structural irregularity, 10.) exterior walls, 11.) roof, 12.) floor, 13.) foundation system. The report illustrates the pratical use of the GEM Building Taxonomy by discussing example case studies, in which the building-specific characteristics are mapped directly using GEM taxonomic attributes and the corresponding taxonomic string is constructed for that building, with "/" slash marks separating attributes. For example, for the building shown to the right, the GEM Taxonomy string is: DX 1 /MUR+CLBRS+MOCL 2 /LWAL 3 / DY/MUR+CLBRS+MOCL/LWAL/YPRE:1939 4 /HEX:2 5 /RES 6 / 7 / 8 /IRRE 9 /10/RSH3+RWO2 11 /FW 12 / 13 / which can be read as (1) Direction = [DX or DY] (the building has the same lateral load-resisting system in both directions); (2) Material = [Unreinforced Masonry + solid fired clay bricks + cement: lime mortar]; (3) Lateral Load-Resisting System = [Wall]; (4) Date of construction = [pre-1939]; (5) Heaight = [exactly 2 storeys]; (6) Occupancy = [residential, unknown type]; (7) Building Position = [unknown = no entry]; (8) Shape of building plan = [unknown = no entry]; (9) Structural irregularity = [regular]; (10) Exterior walls = [unknown = no entry]; (11) Roof = [Shape: pitched and hipped, Roof covering: clay tiles, Roof system material: wood, Roof system type: wood trusses]; (12) Floor = [Floor system: Wood, unknown]; (13) Foundation = [unknown = no entry]. Mapping of GEM Building Taxonomy to selected taxonomies is included in the report -- for example, the above building would be referenced by previous structural taxonomies as: PAGER-STR as UFB or UFB4, by the World Housing Encyclopedia as 7 or 8 and by the European Macroseismic Scale (98) as M5. The Building Taxonomy data model is highly flexible and has been incorporated within a relational database architecture. Due to its ability to represent building typologies using a shorthand form, it is also possible to use the taxonomy for non-database applications, and we discuss possible application of adaptation for Building Information Modelling (BIM) systems, and for the insurance industry. The GEM Building Taxonomy was independently evaluated and tested by the Earthquake Engineering Research Institute (EERI), which received 217 TaxT reports from 49 countries, representing a wide range of building typologies, including single and multi-storey buildings, reinforced and unreinforced masonry, confined masonry, concrete, steel, wood, and earthern buildings used for residential, commercial, industrial, and educational occupancy. Based on these submissions and other feedback, the EERI team validated that the GEM Building Taxonomy is highly functional, robust and able to describe different buildings around the world. The GEM Building Taxonomy is accompanied by supplementary resources. All terms have been explained in a companion online Glossary, which provides both text and graphic descriptions. The Taxonomy is accompanied by TaxT, a computer application that enables a user record information about a building or a building typology using the attributes of the GEM Building Taxonomy v2.0. TaxT can generate a taxonomy string and enable a user to generate a report in PDF format which summarizes the attribute values (s)he has chosen as representative of the building typology under consideration. The report concludes with recommendations for future development of the GEM Building Taxonomy. Appendices provide the detailed GEM Building Taxonomy tables and additional resource, as well as mappings to other taxonomies.

GEM Technical Report↗

Material properties and triggering mechanisms of an andesitic lava dome collapse at Shiveluch Volcano, Kamchatka, Russia, revealed using the finite element method

Shiveluch volcano (Kamchatka, Russia) is an active andesitic volcano with a history of explosive activity, dome extrusion, and structural collapse during the Holocene. The most recent major (> 1 km 3 ) dome collapse occurred in November 1964, producing a ~ 1.5 km 3 debris avalanche that traveled over 15 km from the vent and triggered a phreatic explosion followed by a voluminous (~ 0.8 km 3 ) eruption of juvenile pyroclastic material. Seismic records suggest that the collapse was likely triggered by a magnitude 5.1 earthquake associated with the ascent of magma into the edifice. The geomechanical properties of the pre-1964 dome are unknown; accordingly, the mechanics of the collapse are poorly understood. This project employs numerical slope stability modeling using the finite element method to constrain probable ranges of geomechanical properties for the materials involved in the collapse, considering earthquake loading as the most likely triggering mechanism. Model results show good agreement with the 1964 collapse geometry considering Geological Strength Index and horizontal pseudo-static seismic coefficient ranges of 30–60 and 0.05–0.15 g, respectively, representing variably fractured and altered dome rocks under moderate earthquake loading, confirming that ground acceleration alone could have triggered the dome collapse. Deep-seated rotational sliding is the dominant failure mode, but local extension within the dome during failure appears to play an important role in the development of the collapse. The findings of this work allow for better forward modeling of potential future collapses, the results of which can be incorporated into regional hazard and risk assessments.

Shiveluch Volcano↗

Modeling non-structural strategies to reduce pedestrian evacuation times for mitigating local tsunami threats in Guam

Reducing the potential for loss of life from local tsunamis is challenging for emergency managers given the need for self-protective behavior of at-risk individuals within brief windows of time to evacuate. There has been considerable attention paid to discussing the use of tsunami vertical-evacuation structures for areas where there may be insufficient time to evacuate. This strategy may not be feasible or needed for at-risk populations in island communities for multiple reasons. We examine the influence of three non-structural interventions (reducing departure delays, increasing travel speeds, and managing vegetation to create new paths) that may improve the evacuation potential for at-risk individuals in island communities and use the United States territory of Guam as our case study. We model pedestrian travel times out of a modeled inundation zone for a local tsunami generated by a M w 8.3 earthquake within the Mariana subduction zone. Evacuation-modeling results indicate that reducing departure delays has a larger impact than increasing travel speeds or creating evacuation corridors through heavy brush on reducing the number of at-risk individuals with insufficient time to evacuate. Travel times to safety are less than wave-arrival times for almost all at-risk individuals in the tsunami-hazard zone if one assumes all three interventions are implemented.

Guam↗

Helping coastal communities at risk from tsunamis: the role of U.S. Geological Survey research

In 1946, 1960, and 1964, major tsunamis (giant sea waves usually caused by earthquakes or submarine landslides) struck coastal areas of the Pacific Ocean. In the U.S. alone, these tsunamis killed hundreds of people and caused many tens of millions of dollars in damage. Recent events in Papua New Guinea (1998) and elsewhere are reminders that a catastrophic tsunami could strike U.S. coasts at any time. The USGS, working closely with NOAA and other partners in the National Tsunami Hazard Mitigation Program, is helping to reduce losses from tsunamis through increased hazard assessment and improved real-time warning systems.

Fact Sheet↗

Multiscale seismic imaging of active fault zones for hazard assessment: A case study of the Santa Monica fault zone, Los Angeles, California

High‐resolution seismic reflection profiles at two different scales were acquired across the transpressional Santa Monica Fault of north Los Angeles as part of an integrated hazard assessment of the fault. The seismic data confirm the location of the fault and related shallow faulting seen in a trench to deeper structures known from regional studies. The trench shows a series of near‐vertical strike‐slip faults beneath a topographic scarp inferred to be caused by thrusting on the Santa Monica fault. Analysis of the disruption of soil horizons in the trench indicates multiple earthquakes have occurred on these strike‐slip faults within the past 50 000 years, with the latest being 1000 to 3000 years ago. A 3.8-km-long, high‐resolution seismic reflection profile shows reflector truncations that constrain the shallow portion of the Santa Monica Fault (upper 300 m) to dip northward between 30° and 55°, most likely 30° to 35°, in contrast to the 60° to 70° dip interpreted for the deeper portion of the fault. Prominent, nearly continuous reflectors on the profile are interpreted to be the erosional unconformity between the 1.2 Ma and older Pico Formation and the base of alluvial fan deposits. The unconformity lies at depths of 30–60 m north of the fault and 110–130 m south of the fault, with about 100 m of vertical displacement (180 m of dip‐slip motion on a 30°–35° dipping fault) across the fault since deposition of the upper Pico Formation. The continuity of the uncomformity on the seismic profile constrains the fault to lie in a relatively narrow (50 m) zone, and to project to the surface beneath Ohio Avenue immediately south of the trench. A very high‐resolution seismic profile adjacent to the trench images reflectors in the 15 to 60 m depth range that are arched slightly by folding just north of the fault. A disrupted zone on the profile beneath the south end of the trench is interpreted as being caused by the deeper portions of the trenched strike‐slip faults where they merge with the thrust fault.

Geophysics↗

Paleoseismology of the Sawtooth fault and implications for fault behavior in the epicentral region of the 2020 Mw 6.5 Stanley, Idaho, earthquake

The 2020 moment magnitude ( M w ⁠⁠ ) 6.5 Stanley, Idaho, earthquake raised questions about the history and extent of complex faulting in the northwestern Centennial Tectonic Belt (CTB) and its relation to the Sawtooth normal fault and Eocene Trans‐Challis fault system (TCFS). To explore faulting in this area, we excavated a paleoseismic trench across the Sawtooth fault along the western margin of the CTB, and compared an early Holocene (9.1 ± 2.1 ka, 1δ ⁠ ) rupture at the site with lacustrine paleoseismic data and fault mapping in the 2020 epicentral region. We find: (1) a history of partial to full rupture of the Sawtooth fault ( M w ⁠ 6.8–7.4), (2) that shorter ruptures ( M w ≤ ⁠⁠6.9 ) are likely along distributed and discontinuous faults in the epicentral region, (3) that this complex system that hosted the 2020 earthquake is not directly linked to the Sawtooth fault, (4) that the northeast‐trending TCFS likely plays a role in controlling fault length and rupture continuity for adjacent faults, and (5) that parts of the TCFS may facilitate displacement transfer between normal faults that accommodate crustal extension and rotation. Our results help unravel complex faulting in the CTB and imply that relict structures can help inform regional seismic hazard assessments.

Idaho↗

The third hans cloos lecture. Urban landslides: Socioeconomic impacts and overview of mitigative strategies

As a result of population pressures, hillsides in the world's urban areas are being developed at an accelerating rate. This development increases the risk for urban landslides triggered by rainfall or earthquake activity. To counter this risk, four approaches have been employed by landslide managers and urban planners: (1) restricting development in landslide-prone areas; (2) implementing and enforcing excavation, grading, and construction codes; (3) protecting existing developments by physical mitigation measures and (4) developing and installing monitoring and warning systems. Where they have been utilized, these approaches generally have been effective in reducing the risk due to landslide hazards. In addition to these practices, landslide insurance holds promise as a mitigative measure by reducing the financial impact of landslides on individual property owners. Until recently, however, such insurance has not been widely available and, where it is available, it is so expensive that it has been little used. ?? Springer-Verlag 2006.

Bulletin of Engineering Geology and the Environmen↗

Timescales of surface faulting preservation in low-strain intraplate regions from landscape evolution modeling and the geomorphic and historical record

Large surface-rupturing intraplate earthquakes in stable continental regions (SCRs) are uncommon globally and have recurrence intervals of thousands to hundreds of thousands of years based on the paleoseismic and geomorphic record, challenging accurate active fault identification in these regions. To constrain the timescales of preservation for scarps created by surface ruptures from dip-slip earthquakes, we use a two-dimensional scarp diffusion model for typical intraplate settings and explore which parameters influence fault scarp preservation. These parameters include the coseismic vertical surface offset, the recurrence interval of similar magnitude earthquakes, diffusivity (as a proxy for mean annual precipitation rate), and the erodibility of the surficial material. We constrain parameter ranges from a compilation of historical surface ruptures in intraplate settings in a variety of climates, including the Central and Eastern United States, Australia, Europe, Central Asia (Mongolia, China), India, and West Africa. The timescales of scarp preservation from landscape evolution modeling agree well with observations of scarp preservation in low-strain SCR and intraplate tectonic settings, with some notable exceptions for Australian scarps. We find that the erodibility of the surficial material and earthquake recurrence interval have a stronger effect on the timescales of scarp preservation than diffusivity or coseismic vertical surface offset. Our model results may aid in identifying and characterizing subtle, slow-moving active faults in low-strain SCR and intraplate tectonic settings for different tectonic, geomorphic, and climatic characteristics. Accurate fault locations and characterization from the landscape record has implications for both probabilistic seismic and fault displacement hazard analyses.

Journal of Geophysical Research Solid Earth↗

Spatiotemporal clustering of great earthquakes on a transform fault controlled by geometry

Minor changes in geometry along the length of mature strike-slip faults may act as conditional barriers to earthquake rupture, terminating some and allowing others to pass. This hypothesis remains largely untested because palaeoearthquake data that constrain spatial and temporal patterns of fault rupture are generally imprecise. Here we develop palaeoearthquake event data that encompass the last 20 major-to-great earthquakes along approximately 320 km of the Alpine Fault in New Zealand with sufficient temporal resolution and spatial coverage to reveal along-strike patterns of rupture extent. The palaeoearthquake record shows that earthquake terminations tend to cluster in time near minor along-strike changes in geometry. These terminations limit the length to which rupture can grow and produce two modes of earthquake behaviour characterized by phases of major ( M w 7–8) and great ( M w > 8) earthquakes. Physics-based simulations of seismic cycles closely resemble our observations when parameterized with realistic fault geometry. Switching between the rupture modes emerges due to heterogeneous stress states that evolve over multiple seismic cycles in response to along-strike differences in geometry. These geometric complexities exert a first-order control on rupture behaviour that is not currently accounted for in fault-source models for seismic hazard.

South Island↗

Seismic hazard mapping of California considering site effects

In this paper, we have combined the U.S. Geological Survey's National Seismic Hazard Maps model with the California geologic map showing 17 generalized geologic units that can be defined by their VS30. We regrouped these units into seven VS30 values and calculated a probabilistic seismic hazard map for the entire state for each VS30 value. By merging seismic hazard maps based on the seven different V S30 values, a suite of seismic hazard maps was computed for 0.2 and 1.0 s spectral ordinates at 2% probability of exceedance (PE) in 50 years. The improved hazards maps explicitly incorporate the site effects and their spatial variability on ground motion estimates. The spectral acceleration (SA) at 1.0 s map of seismic shaking potential for California has now been published as California Geological Survey Map Sheet 48.

California↗

Will the 1990’s be a decade of increasingly destructive natural disasters?

Today, there is a considerable body of knowledge about natural hazards that laso enables us to devise effective means to limit the damage they cause. Yet, the reality of the situation is that disasters are increasing in number worldwide, and their social and economic impacts are becoming more nad more difficult to hear.

Earthquakes & Volcanoes (USGS)↗

Community variations in population exposure to near-field tsunami hazards as a function of pedestrian travel time to safety

Efforts to characterize population exposure to near-field tsunami threats typically focus on quantifying the number and type of people in tsunami-hazard zones. To develop and prioritize effective risk-reduction strategies, emergency managers also need information on the potential for successful evacuations and how this evacuation potential varies among communities. To improve efforts to properly characterize and differentiate near-field tsunami threats among multiple communities, we assess community variations in population exposure to tsunamis as a function of pedestrian travel time to safety. We focus our efforts on the multiple coastal communities in Grays Harbor and Pacific Counties (State of Washington, USA), where a substantial resident and visitor population is threatened by near-field tsunamis related to a potential Cascadia subduction zone earthquake. Anisotropic, path-distance modeling is conducted to estimate travel times to safety and results are merged with various population data, including residents, employees, public venues, and dependent-care facilities. Results suggest that there is substantial variability among communities in the number of people that may have insufficient time to evacuate. Successful evacuations may be possible in some communities assuming slow-walking speeds, are plausible in others if travel speeds are increased, and are unlikely in another set of communities given the large distances and short time horizon. Emergency managers can use these results to prioritize the location and determine the most appropriate type of tsunami risk-reduction strategies, such as education and training in areas where evacuations are plausible and vertical-evacuation structures in areas where they are not.

Washington↗

Tearing the terroir: Details and implications of surface rupture and deformation from the 24 August 2014 M6.0 South Napa earthquake, California

The M w 6.0 South Napa earthquake of 24 August 2014 caused slip on several active fault strands within the West Napa Fault Zone (WNFZ). Field mapping identified 12.5 km of surface rupture. These field observations, near-field geodesy and space geodesy, together provide evidence for more than ~30 km of surface deformation with a relatively complex distribution across a number of subparallel lineaments. Along a ~7 km section north of the epicenter, the surface rupture is confined to a single trace that cuts alluvial deposits, reoccupying a low-slope scarp. The rupture continued northward onto at least four other traces through subparallel ridges and valleys. Postseismic slip exceeded coseismic slip along much of the southern part of the main rupture trace with total slip 1 year postevent approaching 0.5 m at locations where only a few centimeters were measured the day of the earthquake. Analysis of airborne interferometric synthetic aperture radar data provides slip distributions along fault traces, indicates connectivity and extent of secondary traces, and confirms that postseismic slip only occurred on the main trace of the fault, perhaps indicating secondary structures ruptured as coseismic triggered slip. Previous mapping identified the WNFZ as a zone of distributed faulting, and this was generally borne out by the complex 2014 rupture pattern. Implications for hazard analysis in similar settings include the need to consider the possibility of complex surface rupture in areas of complex topography, especially where multiple potentially Quaternary-active fault strands can be mapped.

California↗

Pedestrian flow-path modeling to support tsunami evacuation and disaster relief planning in the U.S. Pacific Northwest

Successful evacuations are critical to saving lives from future tsunamis. Pedestrian-evacuation modeling related to tsunami hazards primarily has focused on identifying areas and the number of people in these areas where successful evacuations are unlikely. Less attention has been paid to identifying evacuation pathways and population demand at assembly areas for at-risk individuals that may have sufficient time to evacuate. We use the neighboring coastal communities of Hoquiam, Aberdeen, and Cosmopolis (Washington, USA) and the local tsunami threat posed by Cascadia subduction zone earthquakes as a case study to explore the use of geospatial, least-cost-distance evacuation modeling for supporting evacuation outreach, response, and relief planning. We demonstrate an approach that uses geospatial evacuation modeling to (a) map the minimum pedestrian travel speeds to safety, the most efficient paths, and collective evacuation basins, (b) estimate the total number and demographic description of evacuees at predetermined assembly areas, and (c) determine which paths may be compromised due to earthquake-induced ground failure. Results suggest a wide range in the magnitude and type of evacuees at predetermined assembly areas and highlight parts of the communities with no readily accessible assembly area. Earthquake-induced ground failures could obstruct access to some assembly areas, cause evacuees to reroute to get to other assembly areas, and isolate some evacuees from relief personnel. Evacuation-modeling methods and results discussed here have implications and application to tsunami-evacuation outreach, training, response procedures, mitigation, and long-term land use planning to increase community resilience.

Washington↗

Repeated coseismic uplift of coastal lagoons above the Patton Bay Splay Fault System, Montague Island, Alaska, USA

Coseismic slip on the Patton Bay splay fault system during the 1964 M w 9.2 Great Alaska Earthquake contributed to local tsunami generation and vertically uplifted shorelines as much as 11 m on Montague Island in Prince William Sound (PWS). Sudden uplift of 3.7–4.3 m caused coastal lagoons along the island's northwestern coast to gradually drain. The resulting change in depositional environment from marine lagoon to freshwater muskeg created a sharp, laterally continuous stratigraphic contact between silt and overlying peat. Here, we characterize the geomorphology, sedimentology, and diatom ecology across the 1964 earthquake contact and three similar prehistoric contacts within the stratigraphy of the Hidden Lagoons locality. We find that the contacts signal instances of abrupt coastal uplift that, within error, overlap the timing of independently constrained megathrust earthquakes in PWS—1964 Common Era, 760–870 yr BP, 2500–2700 yr BP, and 4120–4500 yr BP. Changes in fossil diatom assemblages across the inferred prehistoric earthquake contacts reflect ecological shifts consistent with repeated draining of a lagoon system caused by >3 m of coseismic uplift. Our observations provide evidence for four instances of combined megathrust-splay fault ruptures that have occurred in the past ∼4,200 years in PWS. The possibility that 1964-style combined megathrust-splay fault ruptures may have repeated in the past warrants their consideration in future seismic and tsunami hazards assessments.

Alaska↗