Geology Reports⌕ Search

Geology topics

SeanPaul La Selle

Publications and source records attributed to SeanPaul La Selle.

At least 19 recordsLinked to original sources
USGS2026

Sedimentological reconstruction of hurricane storm surges from a carbonate reef coastline at Anegada, British Virgin Islands

Anegada is a small, low-lying carbonate island in the British Virgin Islands that is highly vulnerable to marine inundation from frequent hurricanes. A major challenge in forecasting the timing, intensity, and impacts of future hurricanes to islands such as Anegada is the limited short-term instrumental record that does not capture the full spectrum of inundation events from the frequent but low-magnitude to rare, but high-magnitude events. Washover deposits, which are transported onshore by hurricane-related storm surges and preserved within the onshore stratigraphy, can provide insight into the recurrence of these events. Here, we use grain size distributions, X-ray fluorescence (XRF) data, and an inverse sediment transport model to characterize a series of anomalous sand deposits from shallow hypersaline ponds to determine the variability in storm surges impacting Anegada over at least the last 70 years. The analyses were successfully applied to modern analogues of storms impacting Anegada's southern shoreline (Hurricanes Donna in 1960 at Category 3; Earl in 2010 at Category 4; and Irma in 2017 at Category 5). These deposits are identified by coarse-grained material and elevated marine geochemical signals (Ca, Sr, ln(Br/Ti)). Inverse modelled flow parameters (average flow depth = 0.89–1.46 m relative to the reef flat, average flow speed = 2.95–3.78 m/s) agree with observational data and eyewitness accounts, and the deposit associated with Hurricane Donna indicates the greatest impacts occurred along the southern shoreline (modelled flow depth = 1.24–1.46 m relative to the reef flat, flow speed = 3.49–3.78 m/s). Using the modern analogues as a basis for comparison, we infer the origin of an older sand deposit to be from Hurricane Dog in 1950, when it tracked by Anegada as a Category 3 storm. This methodology provides a transferable framework for interpreting storm recurrence in carbonate-dominated coastlines worldwide that could be used to inform hurricane hazard assessments and risk management efforts.

Anegada↗
USGS2026

The USGS Rotating X-ray Computed Tomography (RXCT) Coral-Core Archive: Scope, access, and standardization

We announce the U.S. Geological Survey (USGS) Rotating X-ray Computed Tomography (RXCT) Coral-Core Archive, a digital resource derived from ~400 coral reef cores curated at the USGS Pacific and St. Petersburg Coastal and Marine Science Centers. The archive de-livers calibrated 3-D image volumes that enable reproducible values of skeletal density, linear extension, and calcification from decadal-to-centennials scale records of coral growth and bioerosion. Cross-study comparability within the archive is supported by a unified RXCT workflow that minimizes imaging artifacts. This includes rejecting im-age-intensity–density regressions with r² < 0.95, back-calculating standard densities to verify a ±10% target precision, and confirming that band-averaged density values fall within published species- and site-specific ranges. Given global coral-reef decline and the rarity of physical archives, we release data under FAIR (Findable, Accessible, Interopera-ble, Reusable) principles. Calibrated imagery and scan metadata are distributed through CoralCache/CoralCT for analysis (DeCarlo et al., 2025), while core locations and collection metadata are published through the USGS Geologic Core and Sample Database (DOI: 10.5066/F7319TR3) with links to CT slice imagery in a USGS ScienceBase repository (DOI: 10.5066/P139Y9H4). This archive provides a powerful dataset for evaluating environ-mental controls on coral growth, establishing restoration baselines, and improving coastal-hazard assessments in the face of global coral-reef declines.

American Somoa, Guam, Florida, Hawaii, U.S. Virgin↗
USGS2025

Increased flood exposure in the Pacific Northwest following earthquake-driven subsidence and sea-level rise

Climate-driven sea-level rise is increasing the frequency of coastal flooding worldwide, exacerbated locally by factors like land subsidence from groundwater and resource extraction. However, a process rarely considered in future sea-level rise scenarios is sudden (over minutes) land subsidence associated with great (>M8) earthquakes, which can exceed 1 m. Along the Washington, Oregon, and northern California coasts, the next great Cascadia subduction zone earthquake could cause up to 2 m of sudden coastal subsidence, dramatically raising sea level, expanding floodplains, and increasing the flood risk to local communities. Here, we quantify the potential expansion of the 1% floodplain (i.e., the area with an annual flood risk of 1%) under low (~0.5 m), medium (~1 m), and high (~2 m) earthquake-driven subsidence scenarios at 24 Cascadia estuaries. If a great earthquake occurred today, floodplains could expand by 90 km 2 (low), 160 km 2 (medium), or 300 km 2 (high subsidence), more than doubling the flooding exposure of residents, structures, and roads under the high subsidence scenario. By 2100, when climate-driven sea-level rise will compound the hazard, a great earthquake could expand floodplains by 170 km 2 (low), 240 km 2 (medium), or 370 km 2 (high subsidence), more than tripling the flooding exposure of residents, structures, and roads under the high subsidence scenario compared to the 2023 floodplain. Our findings can support decision-makers and coastal communities along the Cascadia subduction zone as they prepare for compound hazards from the earthquake cycle and climate-driven sea-level rise and provide critical insights for tectonically active coastlines globally.

California, Oregon, Washington↗
USGS2025

Linking tidal-creek sediment fluxes to vertical sediment accretion in a restored salt marsh

Despite growing interest and investment in salt-marsh restoration, relatively few marshes subjected to restoration efforts have been systematically monitored to assess physical restoration trajectory or success. In south San Francisco Bay, CA, USA, where 83% of wetlands were lost via human manipulation, the largest wetland restoration effort on the U.S. west coast is currently underway, restoring approximately 6,000 ha of former salt-production ponds to mixed habitats. The Whale’s Tail–Cargill Mitigation salt-marsh complex in south San Francisco Bay has a century-long history of drainage, industrial use as salt-production ponds, and subsequent restoration and recovery. Restoration of the 20-ha Cargill Mitigation Marsh was initiated in the late 1990s when the levee surrounding the subsided, former salt-production pond was breached in two locations, enabling conversion back to salt-marsh habitat in the subsequent decades. Here we present time-series measurements of sediment fluxes in the primary tidal creek entering the salt-marsh complex, which are compared to decadal-scale sedimentation patterns determined from repeat elevation surveys and cores collected at the study site. All three methods show net sediment import to the restored marsh. The greatest equivalent sedimentation rates occurred early in the restoration, with generally decreasing rates through time. The long-term average, as determined from cores and expressed as a vertical sedimentation rate, is approximately 1.8 cm yr -1 . Rates from the elevation data are between 1.4 and 2.6 cm yr -1 , with higher rates earlier in the restoration. The most recent estimates, computed from time-series instrument deployments, indicate seasonal variability in sediment import. Annualized rates are lower in winter, approximately 0.1 cm yr -1 , and higher in summer, approximately 1.7 cm yr -1 . Although our measured long-term equivalent sedimentation rates are considerably greater than the current local relative sea-level rise (SLR) of 0.3 cm yr -1 , an increase in SLR or decrease in available suspended sediment would threaten the ability of the marsh to keep pace with SLR and avoid drowning in the future.

California↗
USGS2025

A 700-year rupture sequence of great eastern Aleutian earthquakes from tsunami modeling of stratigraphic records

Great Aleutian underthrusting earthquakes produced destructive tsunamis impacting Hawaiʻi in 1946 and 1957. Prior modeling of the 1957 tsunami deposit and runup records on eastern Aleutian and Hawaiian Islands jointly with tide-gauge observations across the Pacific Ocean constrained a rupture model with shallow slip up to 26 m along 600 km of the plate boundary. Here we implement this modeling approach to older deposits and show alternating deep and shallow megathrust slip up to 26, 32, and 22 m for great earthquakes along the same segment in the 18 th , 15 th , and 14 th centuries. All three modeled prehistoric Aleutian earthquakes produce tsunami inundation in Hawaiʻi with the most severe, 14 th century event having impacts exceeding the 1957 event. The along-dip variability of these four ruptures spanning seven centuries provides insights on earthquake cycles for engineering design and hazard assessment. The 15 th century and 1957 rupture models provide evidence for recurrence of tsunami earthquakes, which can produce disproportionately large tsunamis for a given moment magnitude due to reduced rigidity in the shallow megathrust. The 14 th and 18 th century events likely ruptured deeper regions that did not slip in 1957, suggesting potential for corresponding deeper failure in the next great eastern Aleutian earthquake.

Alaska, Hawaii↗
USGS2025

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

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

British Columbia, California, Oregon, Washington↗
USGS2024

Testing megathrust rupture models using tsunami deposits

The 26 January 1700 CE Cascadia subduction zone earthquake ruptured much of the plate boundary and generated a tsunami that deposited sand in coastal marshes from northern California to Vancouver Island. Although the depositional record of tsunami inundation is extensive in some of these marshes, few sites have been investigated in enough detail to map the inland extent of sand deposition and depict variability in tsunami deposit thickness and grain size. We collected 129 cores in marshes of the Salmon River estuary in Oregon and reanalyzed 114 core logs from a 1987–88 study that mapped the inland extent of circa 1700 CE sandy tsunami deposits. The ca. 1700 CE tsunami deposit in the Salmon River estuary is easily recognized in cores ≤1 m deep in which a buried marsh peat is overlain by a well sorted sand bed with a sharp lower contact that thins and fines inland. We use tsunami deposit data and models of sandy tsunami sediment transport (using Delft3D-FLOW) to test 15 rupture models that could represent a ca. 1700 CE earthquake. At least 12–16 m of slip offshore of the Salmon River, which results in 0.8–1.0 m of coastal coseismic subsidence, is required to match the ca. 1700 CE sand deposit's inland extent, which is consistent with models of heterogeneous megathrust slip in ca. 1700 CE. Our methods of detailed tsunami deposit mapping, combined with sediment transport modeling, can be used to test models of megathrust ruptures and their tsunamis to potentially improve earthquake and tsunami hazard assessments.

JGR Earth Surface↗
USGS2024

A great tsunami earthquake component of the 1957 Aleutian Islands earthquake

The great 1957 Aleutian Islands earthquake ruptured ∼1200 km of the plate boundary along the Aleutian subduction zone and produced a destructive tsunami across Hawaiʻi. Early seismic and tsunami analyses indicated that large megathrust fault slip was concentrated in the western Aleutian Islands, but tsunami waves generated by slip in the west cannot explain the large observed runup in Hawaiʻi far to the southeast. Recently mapped 1957 geologic deposits on eastern Aleutian Islands suggest occurrence of very large nearby slip. Jointly modeling tsunami runup along the eastern Aleutian and Hawaiian Islands together with tide gauge recordings across the Pacific resolves 12-26 m shallow slip along 600 km of the eastern Aleutian Islands in addition to modest, deeper western slip inferred from seismic records. The eastern near-trench slip results in an M W 8.3-8.6 tsunami earthquake component of the M W 8.6-8.8 rupture, comparable in size to the adjacent 1946 Aleutian tsunami earthquake to the east. The reexamination of the 1957 rupture confirms the tsunami hazards posed by the eastern Aleutian subduction zone to Hawaiʻi and lays the groundwork for investigation of large prehistoric earthquakes through modeling tsunami runup inferred from stratigraphic observations to constrain their rupture processes.

Alaska↗
USGS2023

Storm and tsunami overwash sediment transport inferred from recent deposits

Overwash deposits from storms and tsunamis record information about sediment transport and flow that can be used to inform hazard assessments. Here we explore deposits from two extreme wave events: (1) the 2012 Hurricane Sandy, a Category 5 hurricane that is the largest storm in the Atlantic basin on historical record, and (2) the 2011 Tohoku-oki tsunami, created by a 9.0 Mw earthquake, that was up to 20 m high at the coast.

Conference Paper↗
USGS2023

Witnessing history: Comparison of a century of sedimentary and written records in a California protected area

We use a combination of proxy records from a high-resolution analysis of sediments from Searsville Lake and adjacent Upper Lake Marsh and historical records to document over one and a half centuries of vegetation and socio-ecological change—relating to logging, agricultural land use change, dam construction, chemical applications, recreation, and other drivers—on the San Francisco Peninsula. A relatively open vegetation with minimal oak ( Quercus ) and coast redwood ( Sequoia sempervirens ) in the late 1850s reflects widespread logging and grazing during the nineteenth century. Forest and woodland expansion occurred in the early twentieth century, with forests composed of coast redwood and oak, among other taxa, as both logging and grazing declined. Invasive species include those associated with pasturage ( Rume x , Plantago ), landscape disturbance ( Urtica , Amaranthaceae), planting for wood production and wind barriers ( Eucalyptus ), and agriculture. Agricultural species, including wheat, rye, and corn, were more common in the early twentieth century than subsequently. Wetland and aquatic pollen and fungal spores document a complex hydrological history, often associated with fluctuating water levels, application of algaecides, raising of Searsville Dam, and construction of a levee. By pairing the paleoecological and historical records of both lakes, we have been able to reconstruct the previously undocumented impacts of socio-ecological influences on this drainage, all of which overprinted known climate changes. Recognizing the ecological manifestations of these impacts puts into perspective the extent to which people have interacted with and transformed the environment in the transition into the Anthropocene.

California↗
USGS2023

The Searsville Lake Site (California, USA) as a candidate Global boundary Stratotype Section and Point for the Anthropocene series

Cores from Searsville Lake within Stanford University’s Jasper Ridge Biological Preserve, California, USA, are examined to identify a potential GSSP for the Anthropocene: core JRBP2018-VC01B (944.5 cm-long) and tightly correlated JRBP2018-VC01A (852.5 cm-long). Spanning from 1900 CE ± 3 years to 2018 CE, a secure chronology resolved to the sub-annual level allows detailed exploration of the Holocene-Anthropocene transition. We identify the primary GSSP marker as first appearance of 239,240 Pu (372–374 cm) in JRBP2018-VC01B and designate the GSSP depth as the distinct boundary between wet and dry season at 366 cm (6 cm above the first sample containing 239,240 Pu) and corresponding to October-December 1948 CE. This is consistent with a lag of 1–2 years between ejection of 239,240 Pu into the atmosphere and deposition. Auxiliary markers include: first appearance of 137 Cs in 1958; late 20th-century decreases in δ 15 N; late 20th-century elevation in SCPs, Hg, Pb, and other heavy metals; and changes in abundance and presence of ostracod, algae, rotifer and protozoan microfossils. Fossil pollen document anthropogenic landscape changes related to logging and agriculture. As part of a major university, the Searsville site has long been used for research and education, serves users locally to internationally, and is protected yet accessible for future studies and communication about the Anthropocene.

California↗
USGS2022

Diverse tsunamigenesis triggered by the Hunga Tonga-Hunga Ha’apai eruption

On the evening of 15 January 2022, the Hunga Tonga-Hunga Ha’apai volcano 1 unleashed a violent underwater eruption, blanketing the surrounding land masses in ash and debris . The eruption generated tsunamis observed around the world. An event of this type last occurred in 1883 during the eruption of Krakatau , and thus we have the first observations of a tsunami from a large emergent volcanic eruption captured with modern instrumentation. Here we show that the explosive eruption generated waves through multiple mechanisms, including: (1) air–sea coupling with the initial and powerful shock wave radiating out from the explosion in the immediate vicinity of the eruption; (2) collapse of the water cavity created by the underwater explosion; and (3) air–sea coupling with the air-pressure pulse that circled the Earth several times, leading to a global tsunami. In the near field, tsunami impacts are strongly controlled by the water-cavity source whereas the far-field tsunami, which was unusually persistent, can be largely described by the air-pressure pulse mechanism. Catastrophic damage in some harbours in the far field was averted by just tens of centimetres, implying that a modest sea level rise combined with a future, similar event would lead to a step-function increase in impacts on infrastructure. Piecing together the complexity of this event has broad implications for coastal hazards in similar geophysical settings, suggesting a currently neglected source of global tsunamis.

Nature↗
USGS2022

River floodplain abandonment and channel deepening coincide with the onset of clear-cut logging in a coastal California redwood forest

Changes in both land use and climate can alter the balance of transport capacity and sediment supply in rivers. Hence, the primary driver of recent incision or aggradation in alluvial channels is often unclear. The San Lorenzo River on the central coast of California is one location where both climate and land use—specifically, clear-cut forestry of coastal redwoods—could explain recent vertical incision and floodplain abandonment. At our field site on the San Lorenzo, we estimate the magnitude of recent incision using both the ratio of bankfull to critical Shields numbers ( &#x3C4; bf * / &#x3C4; c * "> 𝜏 ∗ 𝑏𝑓 / 𝜏 ∗ 𝑐 ) and the geomorphically effective discharge, calculated from historical gauge data. The Shields number ratio suggests that the normalized bankfull stress of the San Lorenzo River is in the upper 1–2% of West Coast rivers, and the effective discharge corresponds to flow depths ~2–4 m below current bankfull conditions. Radiocarbon ages from detrital charcoal in floodplain sediment reveal active floodplain deposition during the 1600s and possibly into the 1800s, constraining the timing of incision to the last few centuries. Multiple hanging tributaries above the mainstem San Lorenzo River, along with patterns in vegetation on terrace surfaces, corroborate our estimates of the magnitude and timing of incision. Taken together, our findings suggest that floodplain abandonment in this reach was mainly due to methods employed during logging that increased shear stress on the channel bed and reduced sediment storage capacity. We suggest that direct channel modifications in rivers can counterbalance increases in sediment delivery due to clear-cutting, resulting in channel incision rather than aggradation. Today, a young, lower surface appears to be forming adjacent to the San Lorenzo River, which we interpret as an incipient floodplain that is in equilibrium with modern sediment supply and transport capacity.

California↗
USGS2021

Can modeling the geologic record contribute to constraining the tectonic source of the AD 1755 Great Lisbon earthquake?

The precise location of the seismic source of 1755 CE Great Lisbon earthquake is still uncertain. The aim of this work is to use an onland sedimentary record in southern Portugal to test and validate seismic sources for the earthquake. To achieve this, tsunami deposit thicknesses from over 150 cores collected at Salgados in southern Portugal were compared to the results of a tsunami sediment transport model (Delft3D-FLOW) that simulates tsunami propagation, inundation, erosion, and deposition. Five different hypothetical seismic sources were modeled with varying bed roughness coefficients to assess how well they reproduced observed patterns of tsunami deposit thicknesses and dune. Modeled and observed historical tsunami arrival times were also used to test different earthquake sources. Based on these comparisons, three modeled earthquake sources were able to reproduce the observed data, suggesting they should be regarded as somewhat more likely sources for the 1755 earthquake in contrast to four other modeled sources. The fault closest to shore (Marquês de Pombal) yielded the best correlations between model and observations.

Lisbon↗
USGS2021

Onshore flow characteristics of the 1755 CE Lisbon tsunami: Linking forward and inverse numerical modeling

The 1755 CE Lisbon earthquake triggered the largest historical tsunami ever impacting the Atlantic coasts of Europe. Despite recent efforts to better understand this event, there are still unanswered questions about the location of its epicenter and whether physical and historical evidences are in agreement. Inverse modeling using tsunami sediments can be applied to quantify onshore flow characteristics. Forward numerical modeling is also a powerful tool capable of simulating tsunami hydrodynamics and the induced sediment transport. This work presents novel results from a combination of inverse and forward modeling to assess tsunami characteristics onshore. The study site is located on the Portuguese southern coast, at the Salgados lowland where inverse modeling was performed using TsuSedMod ( Jaffe and Gelfenbaum, 2007 ) based on data retrieved from sediment samples. Forward modeling, including tsunami generation and propagation, was performed using the FLOW module of Delft3D suite model. Onshore topography was corrected for the 1755 CE scenario based on extensive tsunami sedimentary deposit thickness data. The tsunami source was chosen based on recent results from the authors that pointed to a good correlation between modeled and field tsunami data for the Marques de Pombal Fault (MPF), Horseshoe Fault (HSF) and a hypothetical scenario represented by a simple combination between Gorringe Bank and Horseshoe Fault (Scenario 1 - SC1). Results from inverse model show tsunami onshore average speed varying from 7.3 up to 9.3 m/s and shear velocities from 0.52 up to 0.66 m/s. Forward modeling results show a wide variation according to the seismic source and tsunami onshore velocities can range from around 7 m/s when considering MPF to even an absence of inundation (SC1). The good agreement between both modeling approaches estimating tsunami velocity confirms the potential of numerical modeling coupled with geological records to improve the understanding of historical events.

Marine Geology↗
USGS2020

Selective sediment transport during Hurricane Sandy on Fire Island (New York, USA): Inferences from heavy-mineral assemblages

In October 2012, Hurricane Sandy caused severe erosion on beaches and dunes of Fire Island (New York, USA). Major shoreline changes occurred with erosional dominance in the upper shoreline and aggradation in the lowermost section of the beach due to the deposition of eroded upper beach and dune sediment. Sand laminations with a high concentration of heavy minerals (“black sand laminations”) were observed in three excavated trenches located on a washover terrace and fan on the east side of Fire Island. The mineralogical composition of these laminations reveals the presence of high quantities of magnetite, ilmenite (as opaque minerals), and garnet (as main translucent mineral). These heavy mineral enriched laminations were formed as waves eroded and transported sand from the primary dune and smaller relict dunes under specific hydrodynamic conditions that promoted grain sorting according to differences in size and specific gravity. Based on the concentrations of certain heavy minerals, the threshold for primarily density-driven sorting probably lies between the specific gravity of a less dense mineral (ilmenite, 4.7) and Garnet almandine (4.3), the most common transparent heavy mineral. The number of laminations and concentrations of heavy minerals vary between trenches on the overwash terrace and appear to be controlled by their distance from sediment sources. The trench with a greatest number of laminations and higher heavy mineral concentrations is located furthest from the main dune but is just 10 m inland of a relict dune that acts as the primary source of sediment of the washover deposit in this trench. A conceptual model for heavy mineral layer deposition is presented based on geomorphological and sedimentological evidences allowing the definition of a density threshold.

New York↗
USGS2020

Sedimentary evidence of prehistoric distant-source tsunamis in the Hawaiian Islands

Over the past 200 years of written records, the Hawaiian Islands have experienced tens of tsunamis generated by earthquakes in the subduction zones of the Pacific "Ring of Fire" (e.g., Alaska-Aleutian, Kuril-Kamchatka, Chile, and Japan). Mapping and dating anomalous beds of sand and silt deposited by tsunamis in low-lying areas along Pacific coasts, even those distant from subduction zones, is critical for assessing tsunami hazard throughout the Pacific basin. We searched for evidence of tsunami inundation using stratigraphic and sedimentologic analyses of potential tsunami deposits beneath present and former Hawaiian wetlands, coastal lagoons, and river floodplains. Coastal wetland sites on the islands of Hawai΄i, Maui, O΄ahu, and Kaua΄i were selected based on historical tsunami runup, numerical inundation modeling, proximity to sandy source sediments, degree of historical wetland disturbance, and breadth of prior geologic and archaeologic investigations. We interpret sand beds containing marine calcareous sediment within peaty and/or muddy wetland deposits on the north and northeastern shores of Kaua΄i, O΄ahu, and Hawai΄i as tsunami deposits. At some sites, deposits of the 1946 and 1957 Aleutian tsunamis are analogs for deeper, older probable tsunami deposits. Radiocarbon-based age models date sand beds from three sites to ~700-500 cal yr B.P., which overlaps ages for tsunami deposits in the eastern Aleutian Islands that record a local subduction zone earthquake (Witter et al., 2016; Witter et al., 2018). The overlapping modeled ages for tsunami deposits at our sites support a plausible correlation with an eastern Aleutian earthquake source for a large prehistoric tsunami in the Hawaiian Islands.

Hawaii↗