Geology topics

Andrea D. Hawkes

Publications and source records attributed to Andrea D. Hawkes.

12 recordsLinked to original sources

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↗

USGS2021

Complex sedimentary processes in large coastal embayments and their potential for coastal morphological and paleo tropical cyclone studies: A case study from Choctawhatchee Bay Western Florida, U.S.A

Storminess and sea-level can both have a significant impact on landforms in cyclone-prone coastal regions, although much of our understanding comes from short-timescale modern observations. This study aims to understand the variability of sediment transport and deposition in the Choctawhatchee Bay/Santa Rosa Island in the northern Gulf of Mexico, establishing the dominant sediment transport processes and morphological response of the barrier system to long-term variations in storminess and rising sea-levels. Here, we study the spatial and temporal changes in physicochemical properties of the sedimentary record of Choctawhatchee Bay to examine the character and fidelity of records of storm impacts spanning the Holocene. Proxies for marine and terrestrial conditions in the cores situated closer to the present barrier (proximal) show that sedimentation in coastal areas and marine influence of the bay during the last ~8000 yrs. were mainly determined by barrier response to the Holocene transgression and changes in storminess. In contrast, sedimentation close to the landward shore was governed by terrigenous input. The correlation of grain size and terrigenous proxies with regional hurricane records indicates that hinterland erosion by the rainfall during hurricane events is likely the dominant terrigenous sediment transport mechanism in areas close to the landward shore of the bay. These results suggest that sediment archives in large coastal deposition environments are equally suitable for sea level and cyclone modulated coastal morphological studies and paleo tropical cyclone studies, depending on the location, selected with an understanding of sedimentation processes in the vicinity.

Florida↗

USGS2021

A maximum rupture model for the central and southern Cascadia subduction zone—reassessing ages for coastal evidence of megathrust earthquakes and tsunamis

A new history of great earthquakes (and their tsunamis) for the central and southern Cascadia subduction zone shows more frequent (17 in the past 6700 yr) megathrust ruptures than previous coastal chronologies. The history is based on along-strike correlations of Bayesian age models derived from evaluation of 554 radiocarbon ages that date earthquake evidence at 14 coastal sites. We reconstruct a history that accounts for all dated stratigraphic evidence with the fewest possible ruptures by evaluating the sequence of age models for earthquake or tsunami contacts at each site, comparing the degree of temporal overlap of correlated site age models, considering evidence for closely spaced earthquakes at four sites, and hypothesizing only maximum-length megathrust ruptures. For the past 6700 yr, recurrence for all earthquakes is 370–420 yr. But correlations suggest that ruptures at ∼1.5 ka and ∼1.1 ka were of limited extent (<400 km). If so, post-3-ka recurrence for ruptures extending throughout central and southern Cascadia is 510–540 yr. But the range in the times between earthquakes is large: two instances may be ∼50 yr, whereas the longest are ∼550 and ∼850 yr. The closely spaced ruptures about 1.6 ka may illustrate a pattern common at subduction zones of a long gap ending with a great earthquake rupturing much of the subduction zone, shortly followed by a rupture of more limited extent. The ruptures of limited extent support the continued inclusion of magnitude-8 earthquakes, with longer ruptures near magnitude 9, in assessments of seismic hazard in the region.

British Columbia, Washington, Oregon, California↗

USGS2020

Minimal stratigraphic evidence for coseismic coastal subsidence during 2000 yr of megathrust earthquakes at the central Cascadia subduction zone

Lithology and microfossil biostratigraphy beneath the marshes of a central Oregon estuary limit geophysical models of Cascadia megathrust rupture during successive earthquakes by ruling out >0.5 m of coseismic coastal subsidence for the past 2000 yr. Although the stratigraphy in cores and outcrops includes as many as 12 peat-mud contacts, like those commonly inferred to record subsidence during megathrust earthquakes, mapping, qualitative diatom analysis, foraminiferal transfer function analysis, and 14 C dating of the contacts failed to confirm that any contacts formed through subsidence during great earthquakes. Based on the youngest peat-mud contact’s distinctness, >400 m distribution, ∼0.6 m depth, and overlying probable tsunami deposit, we attribute it to the great 1700 CE Cascadia earthquake and(or) its accompanying tsunami. Minimal changes in diatom assemblages from below the contact to above its probable tsunami deposit suggest that the lower of several foraminiferal transfer function reconstructions of coseismic subsidence across the contact (0.1–0.5 m) is most accurate. The more limited stratigraphic extent and minimal changes in lithology, foraminifera, and(or) diatom assemblages across the other 11 peat-mud contacts are insufficient to distinguish them from contacts formed through small, gradual, or localized changes in tide levels during river floods, storm surges, and gradual sea-level rise. Although no data preclude any contacts from being synchronous with a megathrust earthquake, the evidence is equally consistent with all contacts recording relative sea-level changes below the ∼0.5 m detection threshold for distinguishing coseismic from nonseismic changes.

British Columbia, California, Oregon, Washington↗

USGS2020

Historically unprecedented Northern Gulf of Mexico hurricane activity from 650 to 1250 CE

Hurricane Michael (2018) was the first Category 5 storm on record to make landfall on the Florida panhandle since at least 1851 CE (Common Era), and it resulted in the loss of 59 lives and $25 billion in damages across the southeastern U.S. This event placed a spotlight on recent intense (exceeding Category 4 or 5 on the Saffir-Simpson Hurricane Wind Scale) hurricane landfalls, prompting questions about the natural range in variability of hurricane activity that the instrumental record is too short to address. Of particular interest is determining whether the frequency of recent intense hurricane landfalls in the northern Gulf of Mexico (GOM) is within or outside the natural range of intense hurricane activity prior to 1851 CE. In this study, we identify intense hurricane landfalls in northwest Florida during the past 2000 years based on coarse anomaly event detection from two coastal lacustrine sediment archives. We identified a historically unprecedented period of heightened storm activity common to four Florida panhandle localities from 650 to 1250 CE and a shift to a relatively quiescent storm climate in the GOM spanning the past six centuries. Our study provides long-term context for events like Hurricane Michael and suggests that the observational period 1851 CE to present may underrepresent the natural range in landfalling hurricane activity.

Gulf of Mexico↗

USGS2020

Identifying the greatest earthquakes of the past 2000 years at the Nehalem River Estuary, Northern Oregon Coast, USA

We infer a history of three great megathrust earthquakes during the past 2000 years at the Nehalem River estuary based on the lateral extent of sharp (≤3 mm) peat-mud stratigraphic contacts in cores and outcrops, coseismic subsidence as interpreted from fossil diatom assemblages and reconstructed with foraminiferal assemblages using a Bayesian transfer function, and regional correlation of 14 C-modeled ages for the times of subsidence. A subsidence contact from 1700 CE (contact A), sometimes overlain by tsunami-deposited sand, can be traced over distances of 7 km. Contacts B and D, which record subsidence during two earlier megathrust earthquakes, are much less extensive but are traced across a 700-m by 270-m tidal marsh. Although some other Cascadia studies report evidence for an earthquake between contacts B and D, our lack of extensive evidence for such an earthquake may result from the complexities of preserving identifiable evidence of it in the rapidly shifting shoreline environments of the lower river and bay. Ages (95% intervals) and subsidence for contacts are: A, 1700 CE (1.1 ± 0.5 m); B, 942–764 cal a BP (0.7 ± 0.4 m and 1.0 m ± 0.4 m); and D, 1568–1361 cal a BP (1.0 m ± 0.4 m). Comparisons of contact subsidence and the degree of overlap of their modeled ages with ages for other Cascadia sites are consistent with megathrust ruptures many hundreds of kilometers long. But these data cannot conclusively distinguish among different types or lengths of ruptures recorded by the three great earthquake contacts at the Nehalem River estuary.

California, Oregon, Washington, British Columbia↗

USGS2019

Salt marsh ecosystem restructuring enhances elevation resilience and carbon storage during accelerating relative sea-level rise

Salt marshes respond to sea-level rise through a series of complex and dynamic bio-physical feedbacks. In this study, we found that sea-level rise triggered salt marsh habitat restructuring, with the associated vegetation changes enhancing salt marsh elevation resilience. A continuous record of marsh elevation relative to sea level that includes reconstruction of high-resolution, sub-decadal, marsh elevation over the past century, coupled with a lower-resolution 1500-year record, revealed that relative sea-level rose 1.5 ± 0.4 m, following local glacial isostatic adjustment (1.2 mm/yr). As sea-level rise has rapidly accelerated, the high marsh zone dropped 11 cm within the tidal frame since 1932, leading to greater inundation and a shift to flood- and salt-tolerant low marsh species. Once the marsh platform fell to the elevation favored by low-marsh Spartina alterniflora , the elevation stabilized relative to sea level. Currently low marsh accretion keeps pace with sea-level rise, while present day high marsh zones that have not transitioned to low marsh have a vertical accretion deficit. Greater biomass productivity, and an expanding subsurface accommodation space favorable for salt marsh organic matter preservation, provide a positive feed-back between sea-level rise and marsh platform elevation. Carbon storage was 46 ± 28 g C/m 2 /yr from 550 to 1800 CE, increasing to 129 ± 50 g C/m 2 /yr in the last decade. Enhanced carbon storage is controlled by vertical accretion rates, rather than soil carbon density, and is a direct response to anthropogenic eustatic sea-level rise, ultimately providing a negative feedback on climate warming.

Massachusetts↗

USGS2016

Relationships between diatoms and tidal environments in Oregon and Washington, USA

A new regional dataset comprising 425 intertidal diatom taxa from 175 samples from 11 ecologically diverse Oregon and Washington estuaries illustrates the importance of compiling a large modern dataset from a range of sites. Cluster analyses and detrended correspondence analysis of the diatom assemblages identify distinct vertical zones within supratidal, intertidal and subtidal environments at six of the 11 study sites, but the abundance of some of the most common species varies widely among and within sites. Canonical correspondence analysis of the regional dataset shows relationships between diatom species and tidal exposure, salinity and substratum (grain size and organic content). Correspondence analyses of local datasets show higher values of explained variation than the analysis of the combined regional dataset. Our results emphasize that studies of the autecology of diatom species require many samples from a range of modern environments to adequately characterize species–environment relationships.

Oregon, Washington↗

USGS2015

Accommodation space, relative sea level, and the archiving of paleo-earthquakes along subduction zones

The spatial variability of Holocene relative sea-level (RSL) change influences the capacities of coastal environments to accommodate a sedimentary record of paleoenvironmental change. In this study we couch a specific investigation in more general terms in order to demonstrate the applicability of the relative sea-level history approach to paleoseismic investigations. Using subsidence stratigraphy, we trace the different modes of coastal sedimentation over the course of time in the eastern Indian Ocean where RSL change evolved from rapidly rising to static from 8000 yr ago to present. Initially, the coastal sites from the Aceh, Sumatra, coastal plain, which are subject to repeated great earthquakes and tsunamis, built up a sedimentary sequence in response to a RSL rise of 1.4 mm/yr. The sequence found at 2 sites 8 km apart contained 3 soils of a mangrove origin ( Rhizophora , Bruguiera / Ceriops , Avicennia pollen, and/or intertidal foraminifera) buried by sudden submergence related to coseismic subsidence and 6 tsunami sands that contain pristine subtidal and planktic foraminifera. After 3800 cal yr B.P. (years before A.D. 1950), sea level stabilized and remained such to the present. The stable relative sea level reduced accommodation space in the late Holocene, suggesting that the continued aggradation of the coastal plain was a consequence of periodic coastal inundation by tsunamis.

Geology↗

USGS2015

Late Holocene sea- and land-level change on the U.S. southeastern Atlantic Coast

Late Holocene relative sea-level (RSL) reconstructions can be used to estimate rates of land-level (subsidence or uplift) change and therefore to modify global sea-level projections for regional conditions. These reconstructions also provide the long-term benchmark against which modern trends are compared and an opportunity to understand the response of sea level to past climate variability. To address a spatial absence of late Holocene data in Florida and Georgia, we reconstructed ~ 1.3 m of RSL rise in northeastern Florida (USA) during the past ~ 2600 years using plant remains and foraminifera in a dated core of high salt-marsh sediment. The reconstruction was fused with tide-gauge data from nearby Fernandina Beach, which measured 1.91 ± 0.26 mm/year of RSL rise since 1900 CE. The average rate of RSL rise prior to 1800 CE was 0.41 ± 0.08 mm/year. Assuming negligible change in global mean sea level from meltwater input/removal and thermal expansion/contraction, this sea-level history approximates net land-level (subsidence and geoid) change, principally from glacio-isostatic adjustment. Historic rates of rise commenced at 1850–1890 CE and it is virtually certain ( P = 0.99) that the average rate of 20th century RSL rise in northeastern Florida was faster than during any of the preceding 26 centuries. The linearity of RSL rise in Florida is in contrast to the variability reconstructed at sites further north on the U.S. Atlantic coast and may suggest a role for ocean dynamic effects in explaining these more variable RSL reconstructions. Comparison of the difference between reconstructed rates of late Holocene RSL rise and historic trends measured by tide gauges indicates that 20th century sea-level trends along the U.S. Atlantic coast were not dominated by the characteristic spatial fingerprint of melting of the Greenland Ice Sheet.

Florida↗

USGS2013

Heterogeneous rupture in the great Cascadia earthquake of 1700 inferred from coastal subsidence estimates

Past earthquake rupture models used to explain paleoseismic estimates of coastal subsidence during the great A.D. 1700 Cascadia earthquake have assumed a uniform slip distribution along the megathrust. Here we infer heterogeneous slip for the Cascadia margin in A.D. 1700 that is analogous to slip distributions during instrumentally recorded great subduction earthquakes worldwide. The assumption of uniform distribution in previous rupture models was due partly to the large uncertainties of then available paleoseismic data used to constrain the models. In this work, we use more precise estimates of subsidence in 1700 from detailed tidal microfossil studies. We develop a 3-D elastic dislocation model that allows the slip to vary both along strike and in the dip direction. Despite uncertainties in the updip and downdip slip extensions, the more precise subsidence estimates are best explained by a model with along-strike slip heterogeneity, with multiple patches of high-moment release separated by areas of low-moment release. For example, in A.D. 1700, there was very little slip near Alsea Bay, Oregon (~44.4°N), an area that coincides with a segment boundary previously suggested on the basis of gravity anomalies. A probable subducting seamount in this area may be responsible for impeding rupture during great earthquakes. Our results highlight the need for more precise, high-quality estimates of subsidence or uplift during prehistoric earthquakes from the coasts of southern British Columbia, northern Washington (north of 47°N), southernmost Oregon, and northern California (south of 43°N), where slip distributions of prehistoric earthquakes are poorly constrained.

British Columbia;Washington;Oregon;California↗

USGS2013

Modern foraminifera, δ13C, and bulk geochemistry of central Oregon tidal marshes and their application in paleoseismology

We assessed the utility of δ 13 C and bulk geochemistry (total organic content and C:N) to reconstruct relative sea-level changes on the Cascadia subduction zone through comparison with an established sea-level indicator (benthic foraminifera). Four modern transects collected from three tidal environments at Siletz Bay, Oregon, USA, produced three elevation-dependent groups in both the foraminiferal and δ 13 C/bulk geochemistry datasets. Foraminiferal samples from the tidal flat and low marsh are identified by Miliammina fusca abundances of > 45%, middle and high marsh by M. fusca abundances of < 45% and the highest marsh by Trochamminita irregularis abundances > 25%. The δ 13 C values from the groups defined with δ 13 C/bulk geochemistry analyses decrease with an increasing elevation; − 24.1 ± 1.7‰ in the tidal flat and low marsh; − 27.3 ± 1.4‰ in the middle and high marsh; and − 29.6 ± 0.8‰ in the highest marsh samples. We applied the modern foraminiferal and δ 13 C distributions to a core that contained a stratigraphic contact marking the great Cascadia earthquake of AD 1700. Both techniques gave similar values for coseismic subsidence across the contact (0.88 ± 0.39 m and 0.71 ± 0.56 m) suggesting that δ 13 C has potential for identifying amounts of relative sea-level change due to tectonics.

Oregon↗