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.

Paleotsunami research is critically important in helping to understand low frequency but large tsunamis. It is also useful for tsunami risk assessments. However, there are many remaining research questions and issues. To help address these questions and issues, post-tsunami geological surveys, which are carried out soon after an event has occurred, are a key to understanding the sedimentary characteristics of tsunamis as modern analogues. Post-tsunami geological surveys are also useful for helping in the recovery process of the local government and the people. In order to perform an efficient post-tsunami geological survey, issues and points that all researchers should be aware of are summarized here. This paper is also useful as a survey manual for paleotsunami and post-tsunami geological surveys both for experts and for non-experts.

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 MW 8.3-8.6 tsunami earthquake component of the MW 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.

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. We rely on the geologic record of great earthquakes and tsunamis, often preserved in the stratigraphy of coastal estuaries, to determine the size, frequency, and location of prehistoric events. At the Cascadia subduction zone, the 1700 CE earthquake caused a meter or more of sudden coastal subsidence and a tsunami that not only inundated coastal areas in Cascadia but was also recorded in Japan. Computer modeling implies that the earthquake occurred around 9 p.m. on 26 January 1700. Although sandy tsunami deposits from the 1700 earthquake are widespread throughout Cascadia, few sites have been cored extensively enough to confidently identify the inland extent of sandy tsunami deposits. We used over 200 sediment cores from the Salmon River estuary in Oregon to map the extent of sandy tsunami deposits from 1700. We then ran numerical models of sediment transport driven by a variety of earthquake generation sources to determine that earthquakes causing at least 0.8 m of subsidence are needed to generate a modeled tsunami capable of recreating the 1700 tsunami deposits observed in cores. This study demonstrates that tsunami deposit mapping and modeling methods can be used to improve our understanding of the impacts of past great earthquakes and tsunamis. Models of sediment transport test earthquake and tsunami sources by comparing the modeled and observed distribution of sandy deposits Sandy deposits can define the minimum inundation of the last major tsunami, around 1700 CE, generated at the Cascadia subduction zone To match tsunami deposit data, models require an earthquake that caused >= 0.8 m of coastal coseismic subsidence at the Salmon River

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.

A San Francisco Estuary core was analysed at high resolution to assess its component stratigraphic signatures of the Anthropocene in the form of non-native species, Hg, spheroidal carbonaceous particles, δ13Corg, δ15N, radiogenic materials, and heavy metals. Time series analysis of the core using Ti data provides a chronology to depth 167 cm into the 1960s. Below this, to depth 230 cm, the lowermost part of the core may extend to the 1950s or potentially a little earlier. The earliest anthropogenic marker recorded in the core is the excursion in Hg (beginning at 190 cm) which may denote the early 1960s and is the closest stratigraphic marker in the core to the proposed mid-20th century timing for the onset of the Anthropocene. Biostratigraphical signatures of non-native species arriving in the 1970s–1980s are widespread key markers and are significant tools for the correlation of Anthropocene deposits across the estuary. The absence of signals that indicate pre-1950s deposits precludes the use of the core to mark the Holocene–Anthropocene boundary. However, the core provides an important reference section to demonstrate the palaeontological distinctiveness of Anthropocene series deposits.

First posted March 31, 2023 For additional information, contact: Pacific Coastal and Marine Science CenterU.S. Geological Survey2885 Mission St.Santa Cruz, CA 95060 Bathymetric change analyses document historical patterns of sediment deposition and erosion, providing valuable insight into the sediment dynamics of coastal systems, including pathways of sediment and sediment-bound contaminants. In 2014 and 2015, the Office for Coastal Management, in partnership with the National Oceanic and Atmospheric Administration (NOAA) Office of Coastal Management, provided funding for new bathymetric surveys of large portions of San Francisco Bay. A total of 93 bathymetric surveys were conducted during this 2-year period, using a combination of interferometric sidescan and multibeam sonar systems. These data, along with recent NOAA, U.S. Geological Survey (USGS), U.S. Army Corps of Engineers, and private contractor surveys collected from 1999 to 2020 (hereinafter referred to as 2010s), were used to create the most comprehensive bathymetric digital elevation models (DEMs) of San Francisco Bay since the 1980s. Comparing DEMs created from these 2010s surveys with USGS DEMs created from NOAA's 1971–1990 (hereinafter referred to as 1980s) surveys provides information on the quantities and patterns of erosion and deposition in San Francisco Bay during the 9 to 47 years between surveys. This analysis reveals that in the areas surveyed in both the 1980s and 2010s, the bay floor lost about 34 million cubic meters of sediment since the 1980s. Results from this study can be used to assess how San Francisco Bay has responded to changes in the system, such as sea-level rise and variation in sediment supply from the Sacramento-San Joaquin Delta and local tributaries, and supports the creation of a new, system-wide sediment budget. This report provides data on the quantities and patterns of sediment volume change in San Francisco Bay for ecosystem managers that are pertinent to various sediment-related issues, including restoration of tidal marshes, exposure of legacy contaminated sediment, and strategies for the beneficial use of dredged sediment.

First posted March 7, 2023 For additional information, contact: Director,Water Resources, Earth System Processes DivisionU.S. Geological Survey345 Middlefield RoadMenlo Park, California, 94025 The South Bay Salt Pond Restoration Project (SBSPRP) encompasses over 6,000 hectares of former salt production ponds along the south edge of the San Francisco Bay and represents the largest wetland restoration effort on the west coast of North America. A series of studies associated with Phase 1 (2010–2018) restoration activities that are focused on a historically mercury contaminated slough and series of ponds within the restoration area have recently been completed. This report brings together the key findings of these loosely coordinated studies and integrates the results into a more comprehensive and holistic product that informs future restoration activities associated with the SBSPRP and elsewhere. This report documents key findings associated with the breach of pond A6: (1) a short-term spike in slough fish (Mississippi silverside) total mercury concentration in lower Alviso Slough; (2) a short-term spike in surface-water particulate total mercury in lower Alviso Slough; (3) significant sediment scour in Alviso Slough adjacent to and downstream of the breach points; (4) a decrease in surface-sediment methylmercury (as a percentage of total mercury) in lower Alviso Slough; (5) the transport of 70 kilograms per year of sediment-associated total mercury into pond A6 during the first 2 years following the breach but with much of this coming from outside of Alviso Slough, presumably from the nearby shallows, Guadalupe Slough, and the larger southern San Francisco Bay area; and (6) a slowing of bed sediment erosion in lower Alviso Slough 3–5 years after the breaching of pond A6. Although this report is not intended to be prescriptive in terms of the next steps the SBSPRP should or should not take, the totality of the findings presented provide critical process-level information regarding the extent and the duration of spikes in mercury levels in water, sediment, fish, and birds, which appeared to result from the two management actions under study. Thus, these results can be used to anticipate similar ecosystem responses associated with similar management actions that may be considered in the future. We also conclude this report by highlighting unanswered questions associated with mercury dynamics as it relates to the restoration project, and possible future directions for research.

Modern microfossil distributions reflect site‐specific habitats and provide an opportunity to assess sediment transport pathways in the nearshore environment. When applied to overwash deposits in the geological record, they provide insight into sediment provenance and transport, factors important for understanding patterns of frequency and intensity of past storms and tsunamis. Modern distribution studies are rare and often the first established ones occur immediately after an overwash event as part of a post‐event field survey. This is problematic because it is unclear what effect overwash events have on nearshore microfossil assemblages and what time interval is necessary for them to return to pre‐event conditions. This study documents the impacts of Hurricane Irma on nearshore sediments off the coast of Anegada, British Virgin Islands, using distributions of Homotrema rubrum , an encrusting foraminifer with a defined provenance in coral reefs. At four sampling intervals spanning two years, from six months pre‐Hurricane Irma to eighteen months after, surface sediment was collected from three transects on the northern and southern shores of the island. Partitioning Around Medoids cluster analysis revealed that Hurricane Irma introduced an influx of well‐preserved fragments into the reef flat and made the sediments more uniform, limiting the foraminifer’s utility as a known sediment transport indicator. The mixing of sediments along the two northern transects (reef proximal) persisted for seven to eighteen months before returning to near pre‐hurricane conditions. However, the southern transect (absence of reef), where Homotrema rubrum concentrations are significantly less, failed to recover within the time period assessed by this study, indicating a variable recovery period between Atlantic Ocean and Caribbean Sea facing shorelines. Results from this study suggest that a waiting period of at least eighteen months after a major storm is recommended before collecting surface sediment from the nearshore environments of reef‐dominated coastlines.

On the evening of 15 January 2022, the Hunga Tonga-Hunga Ha’apai volcano1 unleashed a violent underwater eruption, blanketing the surrounding land masses in ash and debris2,3. The eruption generated tsunamis observed around the world. An event of this type last occurred in 1883 during the eruption of Krakatau4, 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. January 2022 saw the first observations of a tsunami resulting from a large emergent volcanic eruption (Hunga Tonga) captured using modern instrumentation, with broad implications for hazard management in similar geophysical settings.

Periods of very shallow water (water depth in the order of 10 cm) occur daily on tidal flats because of the propagation of tides over very gently sloping beds, leading to distinct morphodynamical phenomena. To improve the understanding of the characteristics of velocity and suspended sediment concentration (SSC) surges and their contribution to sediment transport and local bed changes during periods of very shallow water, measurements of near-bed flow, and SSC were carried out at two cross-shore locations on an intertidal flat along the Jiangsu coast, China. Furthermore, the role of surges in local resuspension and morphological change was explored. Results indicate that flow and SSC surges occurred at both stations during very shallow water periods. On the lower intertidal flat, flood surges were erosive, while weaker surges on the middle intertidal flat were not. Surges on lower intertidal flats resulted in local resuspension and strong turbidity, contributing up to 25% of the onshore-suspended sediment flux during flood tides, even though they last only 10% of the flood duration. When surges travel across the flats, conditions change from erosional to depositional. Velocity surges on the middle intertidal flat were too weak to resuspend bed sediment, and the associated SSC surges were produced by advection.

The lack of historic seismicity along the Cascadia subduction zone (CSZ) fuels the uncertainty of the megathrust rupture process and resulting hazards in the US and Canadian Pacific Northwest. Diverse types of paleoseismic proxies uniquely record different aspects of the earthquake rupture process. Thus, margin-wide integration of paleoseismic proxy records leverages these unique transcripts to provide a richer understanding of past earthquake variability in time and space. In a series of virtual workshops, the USGS Powell Center research group on CSZ earthquake science and hazards have compiled paleoseismic records of coseismic coastal land-level change, tsunami inundation, and ground shaking proxies such as landslides, liquefaction, and turbidites. To consistently integrate different datasets, we constructed ranking schemes to quasi-quantitatively evaluate the quality of geochronologic age control and certainty of earthquake evidence at 49 coastal land level change sites, 81 tsunami sites, 76 terrestrial paleoseismic shaking sites, and 13 marine paleoseismic shaking sites. At each site, we evaluated how well geochronologic samples and results (14C dates, stratigraphy, geological context etc.) constrain the timing of candidate earthquake events. Similarly, we applied an evidence quality ranking scheme to assess how certain the observations of coastal land-level change, tsunami inundation, and strong ground motion were the result of a CSZ megathrust earthquake source. We also focused on historic records of non-local events, such as the 1964 Alaska earthquake and the 1992 Cape Mendocino earthquake, to identify the characteristics of distal subduction zone earthquakes and small local events, respectively, in the paleoseismic record. In the future we plan to identify the likely rupture characteristics of major CSZ earthquakes during the Holocene and to synthesize paleoseismic, geophysical, and instrumental datasets to identify the existence and persistence of rupture boundaries along the subduction zone.

ABSTRACT 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 laminae with a high concentration of heavy minerals (“black sand laminae”) were observed in three excavated trenches located on a washover terrace and fan on the east side of Fire Island. The mineral composition of these laminae reveals the presence of high quantities of magnetite, ilmenite (as opaque minerals), and garnet (as the main translucent mineral). These heavy-mineral-enriched laminae 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 less dense opaque mineral (ilmenite, 4.7) and almandine (4.3), the most common transparent heavy mineral. The number of laminae 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 the greatest number of laminae and higher heavy-mineral concentrations is located farthest 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 deposition of heavy-mineral layers is presented based on geomorphological and sedimentological evidence allowing the definition of a density threshold.

Organic geochemistry is commonly used in environmental studies. In tsunami research, however, its applications are in their infancy and it is still rarely used. We present results for two types of organic geochemical markers, biomarkers and anthropogenic markers, present in deposits left by 2011 Tohoku-oki tsunami on the Sendai Plain, Japan. As the tsunami inundated the coastal lowland up to 4.85 km inland, sediments from various sources were eroded, transported and deposited. This led to the distribution of biomarkers from different sources across the Sendai Plain creating a unique geochemical signature in the tsunami deposits. The tsunami also caused destruction along the Sendai coastline, leading to the release of large quantities of environmental pollutants (e.g., fossil fuels, tarmac, pesticides, plastics, etc.) that were distributed across the inundated area. These anthropogenic markers, represented by three main compound groups (polycyclic aromatic hydrocarbons, pesticides, and halogenated compounds), were preserved in tsunami deposits (at least until 2013, prior to land clearing). Their concentrations differed significantly from the pre- and post-tsunami background contamination levels. Organic proxy concentrations can differ for sand and mud deposits due to various factors (e.g., preservation, dilution, microbial alteration). However, it can be concluded that anthropogenic markers and biomarkers have the potential to be a valuable proxy for future studies of recent tsunami deposits because of their high source specificity and relatively good preservation potential providing information about sediment sources and transport pathways (e.g., marine source, evidence of backwash).

Far‐field tsunami deposits observed in the Kahana Valley, O‘ahu, Hawai‘i (USA), were investigated for their organic‐geochemical content. During short high‐energy events, (tsunamis and storms) organic and chemical components are transported with sediment from marine to terrestrial areas. This study investigates the use of anthropogenic based organic geochemical compounds (such as polycyclic aromatic hydrocarbons, pesticides and organochlorides) as a means to identify tsunami deposits. Samples were processed by solid–liquid extraction and analyzed using gas chromatography–mass spectrometry. A total of 21 anthropogenic marker compounds were identified, of which 11 compounds were selected for detailed analysis. Although the tsunami deposits pre‐date industrial activity in Hawai‘i by several hundred years, distinct changes were found in the concentrations of anthropogenic marker compounds between sandy tsunami deposits and the surrounding mud/peat layers, which may help in identifying tsunami deposits within cores. As expected, low overall concentrations of anthropogenic markers and pollutants were observed due to the lack of industrial input‐sources and little anthropogenic environmental impact at the study site. This geochemical characterization of tsunami deposits shows that anthropogenic markers have significant potential as another high‐resolution, multi‐proxy method for identifying tsunamis in the sedimentary record.