Geology topics

Niamh Cahill

Publications and source records attributed to Niamh Cahill.

10 recordsLinked to original sources

USGS2022

Reproducibility and variability of earthquake subsidence estimates from saltmarshes of a Cascadia estuary

We examine fossil foraminiferal assemblages from 20 sediment cores to assess sudden relative sea-level (RSL) changes across three mud-over-peat contacts at three salt marshes in northern Humboldt Bay, California (~44.8°N, -124.2°W). We use a validated foraminiferal-based Bayesian transfer function to evaluate the variability of subsidence stratigraphy at a range of 30-6000 m across an estuary. We use the consistency in RSL reconstructions to support estimates of coseismic subsidence from megathrust earthquakes. To assess the variability of subsidence estimates, we analyzed: nine examples of the 1700 CE earthquake (average of 0.64 ±0.14 m subsidence; range of 0.24 ±0.27 to 1.00 ±0.44 m), five examples of the ca. 875 cal a BP earthquake (average of 0.43 ±0.16 m; range of 0.41 ±0.36 to 0.48 ±0.39 m), and six examples of the ca. 1120 cal a BP earthquake (average of 0.70±0.18 m; range of 0.47 ±0.36 to 0.80 ±0.49 m). Our subsidence estimate results suggest ~±0.3 m of within-site (intra-site) variability, which is consistent with previous research. We also identify inconsistencies between sites (inter-site) at northern Humboldt Bay greater than one-sigma uncertainties, driven by variable foraminiferal assemblages in the mud overlying the 1700 CE subsidence contact. Therefore, we recommend at least two quantitative microfossil reconstructions across the same stratigraphic sequence from different marsh sites within an estuary to account for estimate variability and provide increased confidence in vertical coseismic deformation estimates. Our results have broad implications for quantitative, microfossil-based reconstructions of coseismic subsidence at temperate coastlines globally.

California↗

USGS2021

Timing and amount of southern Cascadia earthquake subsidence over the past 1700 years at northern Humboldt Bay, California, USA

Stratigraphic, lithologic, foraminiferal, and radiocarbon analyses indicate that at least four abrupt mud-over-peat contacts are recorded across three sites (Jacoby Creek, McDaniel Creek, and Mad River Slough) in northern Humboldt Bay, California, USA (∼44.8°N, −124.2°W). The stratigraphy records subsidence during past megathrust earthquakes at the southern Cascadia subduction zone ∼40 km north of the Mendocino Triple Junction. Maximum and minimum radiocarbon ages on plant macrofossils from above and below laterally extensive (>6 km) contacts suggest regional synchroneity of subsidence. The shallowest contact has radiocarbon ages that are consistent with the most recent great earthquake at Cascadia, which occurred at 250 cal yr B.P. (1700 CE). Using Bchron and OxCal software, we model ages for the three older contacts of ca. 875 cal yr B.P., ca. 1120 cal yr B.P., and ca. 1620 cal yr B.P. For each of the four earthquakes, we analyze foraminifera across representative mud-over-peat contacts selected from McDaniel Creek. Changes in fossil foraminiferal assemblages across all four contacts reveal sudden relative sea-level (RSL) rise (land subsidence) with submergence lasting from decades to centuries. To estimate subsidence during each earthquake, we reconstructed RSL rise across the contacts using the fossil foraminiferal assemblages in a Bayesian transfer function. The coseismic subsidence estimates are 0.85 ± 0.46 m for the 1700 CE earthquake, 0.42 ± 0.37 m for the ca. 875 cal yr B.P. earthquake, 0.79 ± 0.47 m for the ca. 1120 cal yr B.P. earthquake, and ≥0.93 m for the ca. 1620 cal yr B.P. earthquake. The subsidence estimate for the ca. 1620 cal yr B.P. earthquake is a minimum because the pre-subsidence paleoenvironment likely was above the upper limit of foraminiferal habitation. The subsidence estimate for the ca. 875 cal yr B.P. earthquake is less than (<50%) the subsidence estimates for other contacts and suggests that subsidence magnitude varied over the past four earthquake cycles in southern Cascadia.

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

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↗

USGS2017

Reconstructing Common Era relative sea-level change on the Gulf Coast of Florida

To address a paucity of Common Era data in the Gulf of Mexico, we reconstructed ~ 1.1 m of relative sea-level (RSL) rise over the past ~ 2000 years at Little Manatee River (Gulf Coast of Florida, USA). We applied a regional-scale foraminiferal transfer function to fossil assemblages preserved in a core of salt-marsh peat and organic silt that was dated using radiocarbon and recognition of pollution, 137 Cs and pollen chronohorizons. Our proxy reconstruction was combined with tide-gauge data from four nearby sites spanning 1913–2014 CE. Application of an Errors-in-Variables Integrated Gaussian Process (EIV-IGP) model to the combined proxy and instrumental dataset demonstrates that RSL fell from ~ 350 to 100 BCE, before rising continuously to present. This initial RSL fall was likely the result of local-scale processes (e.g., silting up of a tidal flat or shallow sub-tidal shoal) as salt-marsh development at the site began. Since ~ 0 CE, we consider the reconstruction to be representative of regional-scale RSL trends. We removed a linear rate of 0.3 mm/yr from the RSL record using the EIV-IGP model to estimate climate-driven sea-level trends and to facilitate comparison among sites. This analysis demonstrates that since ~ 0 CE sea level did not deviate significantly from zero until accelerating continuously from ~ 1500 CE to present. Sea level was rising at 1.33 mm/yr in 1900 CE and accelerated until 2014 CE when a rate of 2.02 mm/yr was attained, which is the fastest, century-scale trend in the ~ 2000-year record. Comparison to existing reconstructions from the Gulf coast of Louisiana and the Atlantic coast of northern Florida reveal similar sea-level histories at all three sites. We explored the influence of compaction and fluvial processes on our reconstruction and concluded that compaction was likely insignificant. Fluvial processes were also likely insignificant, but further proxy evidence is needed to fully test this hypothesis. Our results indicate that no significant Common Era sea-level changes took place on the Gulf and southeastern Atlantic U.S. coasts until the onset of modern sea-level rise in the late 19th century.

Florida↗

USGS2017

Extended late Holocene relative sea-level histories for North Carolina, USA

We produced ∼3000-year long relative sea-level (RSL) histories for two sites in North Carolina (USA) using foraminifera preserved in new and existing cores of dated salt-marsh sediment. At Cedar Island, RSL rose by ∼2.4 m during the past ∼3000 years compared to ∼3.3 m at Roanoke Island. This spatial difference arises primarily from differential GIA that caused late Holocene RSL rise to be 0.1–0.2 mm/yr faster at Roanoke Island than at Cedar Island. However, a non-linear difference in RSL between the two study regions (particularly from ∼0 CE to ∼1250 CE) indicates that additional local- to regional-scale processes drove centennial-scale RSL change in North Carolina. Therefore, the Cedar Island and Roanoke Island records should be considered as independent of one another. Between-site differences on sub-millennial timescales cannot be adequately explained by non-stationary tides, sediment compaction, or local sediment dynamics. We propose that a period of accelerating RSL rise from ∼600 CE to 1100 CE that is present at Roanoke Island (and other sites north of Cape Hatteras at least as far as Connecticut), but absent at Cedar Island (and other sites south of Cape Hatteras at least as far as northeastern Florida) is a local-to regional-scale effect of dynamic ocean and/or atmospheric circulation.

North Carolina↗

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

Sea-level change during the last 2500 years in New Jersey, USA

Relative sea-level changes during the last ∼2500 years in New Jersey, USA were reconstructed to test if late Holocene sea level was stable or included persistent and distinctive phases of variability. Foraminifera and bulk-sediment δ 13 C values were combined to reconstruct paleomarsh elevation with decimeter precision from sequences of salt-marsh sediment at two sites using a multi-proxy approach. The additional paleoenvironmental information provided by bulk-sediment δ 13 C values reduced vertical uncertainty in the sea-level reconstruction by about one third of that estimated from foraminifera alone using a transfer function. The history of sediment deposition was constrained by a composite chronology. An age–depth model developed for each core enabled reconstruction of sea level with multi-decadal resolution. Following correction for land-level change (1.4 mm/yr), four successive and sustained (multi-centennial) sea-level trends were objectively identified and quantified (95% confidence interval) using error-in-variables change point analysis to account for age and sea-level uncertainties. From at least 500 BC to 250 AD, sea-level fell at 0.11 mm/yr. The second period saw sea-level rise at 0.62 mm/yr from 250 AD to 733 AD. Between 733 AD and 1850 AD, sea level fell at 0.12 mm/yr. The reconstructed rate of sea-level rise since ∼1850 AD was 3.1 mm/yr and represents the most rapid period of change for at least 2500 years. This trend began between 1830 AD and 1873 AD. Since this change point, reconstructed sea-level rise is in agreement with regional tide-gauge records and exceeds the global average estimate for the 20th century. These positive and negative departures from background rates demonstrate that the late Holocene sea level was not stable in New Jersey.

New Jersey↗