Recurrence Intervals for Great Earthquakes of the Past 3,500 Years at Northeastern Willapa Bay, Washington

Plate-boundary earthquakes have occurred repeatedly in the past several thousand years at the Cascadia subduction zone, where they are widely recorded by buried marsh and forest soils beneath estuarine wetlands. This report adds to previous accounts of such soils along the Pacific coast of southern Washington State. Our new evidence comes from outcrop surveys, diatom analyses, and radiocarbon dating of soils exposed about 10 km apart in banks of the Niawiakum and Willapa Rivers, tidal arms of Willapa Bay. This new evidence clarifies the timing of great (magnitude 8 or larger) earthquakes during the past 3,500 years at this part of the Cascadia subduction zone. All the surveyed outcrops display buried soils that probably record tectonic subsidence during earthquakes. Although in many cased we cannot rule out every alternative to such coseismic subsidence, we found no buried soil that is better explained by stream migration, storm, river flood, sea-level rise, barrier breaching, or sediment compaction. Each of the soils probably represents a marsh or forest that suddenly became a tidal flat and consequently was buried by tidal mud. In nearly every case the soils have too much lateral continuity and too little relief to record cutting and filling by tidal streams. Vascular plant fossils preserved within and above many of the soils show that storms or floods, if unaccompanied by lasting submergence, cannot account for burial of the soils. Where their remains are preserved, plants that had lived on the soils belong to species indicative of high parts of tidal marshes, or of tidal swamps or uplands. By contrast, the main vascular -plant species preserved in mud above the soils is Triglochin maritimum, a colonist of saltwater mudflats in southern coastal Washington. Assemblages of diatoms verify vascular-plant evidence for lasting submergence and show that such submergence is recorded also by buried soils with which vascular-plant fossils are not preserved. In an outcrop along the Niawiakum River, for example, diatoms show that each of six successive soils represents a high marsh or an upland, and that mud above each of these soils represents an intertidal or subtidal mudflat. Fossils further show that the change from high marsh or upland to mudflat probably happened too fast to have resulted from a gradual rise in sea level. For every soil studied for fossil diatoms, diatom assemblages imply that this change happened suddenly, without transition through low marsh. Gradual sea-level rise is further precluded for some soils by remains of vascular plants that had lived on buried soils. These remains include stems and leaves that were surrounded by mudflat deposits before they had time to decompose. Other alternatives to coseismic subsidence can be discounted as well. Submergence from breaching of a baymouth barrier is ruled out by soils of brackish marshes that require tidal connection with the sea. Localized settlement from earthquake-induced compaction of unconsolidated Holocene deposits does not explain the presence of buried soils directly above well-consolidated Pleistocene deposits. Evidence against alternative explanation at one outcrop can be extrapolated to other outcrops by radiocarbon and stratigraphic correlation of buried soils. Plate-boundary earthquakes probably account for all the subsidence events. The plate boundary is the only recognized fault common to all areas having evidence of coseismic subsidence in southern coastal Washington. Although some of these areas coincide with mapped late Cenozoic synclines, where coseismic subsidence might accompany earthquakes on faults in the North American plate, others are outside such synclines. The coseismically subsided areas include part of welt of Eocene basement rock 40 km long and 15 km wide, herein termed the South Bend antiform. This structural high, the largest in southern coastal Washington, has probably grown in late Cenozoic time. The South Bend Antiform should subside during plate-boundary earthquakes that flex the North American plate throughout southern coastal Washington. By contrast, the antiform might grow upward during an upper-plate earthquake that produces subsidence only in late Cenozoic synclines that flank the antiform. With these expectations in mind, we compared outcrops off the South Bend antiform (Niawiakum River) with outcrops on the antiform (Willapa River). We found no difference in sense of timing of earthquake-induced changes in land level. Like the surveyed outcrops along the Niawiakum River, the surveyed outcrops along the Willapa River contain buried soils indicative of earthquake-induced subsidence. the longest stratigraphic sequence exposed along both streams contain six or seven soils less than 3,500 years old, and these sequences correlate with one another on the basis of soil horizons, spruce roots, and radiocarbon ages. Though every buried soil identified in the surveyed outcrops probably records a plate-boundary earthquake is widely recorded by a buried soil. One buried soil was obliterated at many sites through centuries of decomposition, probably because shallow burial left the soil high in the weathering profile of a succeeding soil. Such decomposition also destroyed most of the organic matter associated with two other soils. In addition, a widespread soil might not be available for burial if an earthquake occurs too soon after its predecessor for much rebuilding of tidal marshes. Despite these limitations, the composite record from buried soils in large outcrops of northeastern Willapa Bay probably includes every earthquake in the area during the past 3,500 years that caused at least .5 meters of widespread coseismic subsidence and followed the preceding great earthquake by more than a century. Willapa Bay's earthquake history in the past 3,500 years probably includes seven events, each comprising a single rupture or multiple contiguous ruptures on the Washington part of the Cascadia plate boundary. Each event probably included at least one great earthquake, as judged from likely rupture widths inferred from modern geophysical evidence, likely rupture lengths inferred from coastwise correlation of buried soils, extensive sea-floor displacement inferred from a tsunami in Japan, and seismic-moment release deduced from plate motions and average recurrence intervals. The history begins with three events between 3,300-3,500 years ago and 2,400-2,800 years ago (ranges include estimated 95-percent confidence interval). The next recorded events occurred 1,500-1,700 and 1.130-1,350 years ago. They were followed by a poorly dated event that probably occurred before 900 years ago and may have been associated with rupture in inland faults in the North American plate 1,000-1,100 years ago. The most recent of the events happened close to 300 years ago, probably in January 1700. The six intervals between events in this inferred history average 500-540 years but range from about one to three centuries to about a millennium. The earliest two intervals sum to 540-1,100 years. The first of them may have been the longer, as judged from spruce roots that may indicate prolonged interseismic emergence. Next came an interval of 700-1,300 years, when spruce forests spread onto emerging tidal marshes and decomposition largely destroyed an underlying buried soil. This exceptionally long interval was followed by two short ones that together spanned no more than 800 years. The most recent complete interval, marked by another spreading of spruce forests and decomposition of earlier buried soils, lasted 600-1,000 years. This pattern of long and short recurrence intervals at Willapa Bay may match the pattern of intervals between turbidity currents in Cascadia Channel , on the abyssal sea floor 200 km off the central Oregon coast. Previous work showed that these currents largely originated at submarine canyon heads about 50 km west of Willapa Bay, and that great earthquakes may have generated 13 currents in the past 7,500 years. Although pelagic layers between the turbidites have been interpreted as evidence for recurrence intervals of fairly uniform duration, borrows in the turbidites suggest variability in recurrence intervals inferred from buried soils at Willapa Bay.

[1]  A. Moore,et al.  A Tsunami About 1000 Years Ago in Puget Sound, Washington , 1992, Science.

[2]  J. Pizzuto,et al.  Sea-level rise acceleration and the drowning of the Delaware Bay coast at 1.8 ka , 1993 .

[3]  W. Thatcher The earthquake deformation cycle at the Nankai Trough, southwest Japan , 1984 .

[4]  S. R. Holdahl,et al.  Postseismic crustal uplift near Anchorage, Alaska , 1977 .

[5]  B. F. Atwater,et al.  Sudden, probably coseismic submergence of Holocene trees and grass in coastal Washington State , 1991 .

[6]  Steven N. Ward,et al.  The 1960 Chile earthquake: inversion for slip distribution from surface deformation , 1990 .

[7]  H. Clifton Discrimination between Subtidal and Intertidal Facies in Pleistocene Deposits, Willapa Bay, Washington , 1983 .

[8]  R. Minor,et al.  Earthquake-Induced Subsidence and Burial of Late Holocene Archaeological Sites, Northern Oregon Coast , 1996, American Antiquity.

[9]  A. Nelson,et al.  Diatom zonation in southern Oregon tidal marshes relative to vascular plants, foraminifera, and sea level , 1993 .

[10]  Joe W. Johnson Characteristics and Behavior of Pacific Coast Tidal Inlets , 1973 .

[11]  Preliminary estimates of the strength of prehistoric shaking in the Columbia River valley and the southern half of coastal Washington, with emphasis for a Cascadia subduction zone earthquake about 300 years ago , 1995 .

[12]  C. B. Davis,et al.  Salt tolerance of five California salt marsh plants. , 1970 .

[13]  John Adams,et al.  Paleoseismicity of the Cascadia Subduction Zone: Evidence from turbidites off the Oregon‐Washington Margin , 1990 .

[14]  H. C. Cowles Submerged Forests , 1915, Botanical Gazette.

[15]  A. Jennings,et al.  Foraminiferal assemblage zones in Oregon tidal marshes; relation to marsh floral zones and sea level , 1992 .

[16]  J. Clague,et al.  Detection of Large Prehistoric Earthquakes in the Pacific Northwest by Microfossil Analysis , 1994, Science.

[17]  Thomas H. Heaton,et al.  Source characteristics of hypothetical subduction earthquakes in the northwestern United States , 1986 .

[18]  I. Shennan,et al.  Analysis and interpretation of Holocene sea-level data , 1983, Nature.

[19]  D. Chen,et al.  AIR ENTRY IN SALT MARSH SEDIMENTS H. F. HEMOND , 1990 .

[20]  M. Stein,et al.  Long-term earthquake clustering: A 50,000-year paleoseismic record in the Dead Sea Graben , 1996 .

[21]  H. Kelsey,et al.  Topographic form of the Coast Ranges of the Cascadia Margin in relation to coastal uplift rates and plate subduction , 1994 .

[22]  Late-Holocene sea-level record in a core from the Puget Lowland , 1987 .

[23]  H. Kanamori The energy release in great earthquakes , 1977 .

[24]  Kenji Satake,et al.  Time and size of a giant earthquake in Cascadia inferred from Japanese tsunami records of January 1700 , 1996, Nature.

[25]  S. Kennedy,et al.  Subrecent Intertidal Diatoms from St. Catherines Island, Georgia: Taphonomic Implications , 1989 .

[26]  J. C. Savage,et al.  Strain accumulation in western Washington , 1991 .

[27]  Herb Dragert,et al.  Continuous GPS monitoring of elastic strain in the Northern Cascadia Subduction Zone , 1995 .

[28]  R. Wells,et al.  Cenozoic evolution of the continental margin of Oregon and Washington , 1991 .

[29]  M. Umitsu,et al.  Stratigraphic Radiocarbon and Diatom Indicato rs of Sudden Submergence along the Nesselle River Mouth, Washington, USA , 1995 .

[30]  G. Griggs,et al.  Sedimentation in Cascadia Deep-Sea Channel , 1970 .

[31]  A. Nelson Holocene Tidal-Marsh Stratigraphy in South-Central Oregon-Evidence for Localized Sudden Submergence in the Cascadia Subduction Zone , 1992 .

[32]  K. Fujimoto,et al.  Late Holocene Sea-Level Fluctuations and Mangrove Forest Formation on Ponape Island, Micronesia , 1990 .

[33]  R. Bucknam,et al.  Abrupt Uplift Within the Past 1700 Years at Southern Puget Sound, Washington , 1992, Science.

[34]  Garry C. Rogers,et al.  An assessment of the megathrust earthquake potential of the Cascadia subduction zone , 1988 .

[35]  C. Scholz,et al.  On the mechanism of seismic decoupling and back arc spreading at subduction zones , 1995 .

[36]  R. Yeats,et al.  Oblique strike-slip faulting of the central Cascadia , 1997 .

[37]  D. Johnson The New England-Acadian shoreline , 1926 .

[38]  R. Hyndman,et al.  Seismicity and rates of relative motion on the plate boundaries of Western North America , 1983 .

[39]  E. Oele,et al.  The Quaternary history of the North Sea , 1979 .

[40]  B. F. Atwater,et al.  Evidence for Great Holocene Earthquakes Along the Outer Coast of Washington State , 1987, Science.

[41]  H. Kanamori,et al.  VARIABLE RUPTURE MODE OF THE SUBDUCTION ZONE ALONG THE ECUADOR-COLOMBIA COAST , 1982 .

[42]  K. Berryman,et al.  Holocene paleoseismicity in the fold and thrust belt of the Hikurangi subduction zone, eastern North Island, New Zealand , 1989 .

[43]  M. Darienzo,et al.  Stratigraphic Evidence for Great Subduction-zone Earthquakes at Four Estuaries in Northern Oregon, U.S.A. , 1994 .

[44]  H. Streif The Results of Stratigraphical and Facial Investigations in the Coastal Holocene of Woltzeten/Ostfriesland, Germany , 1972 .

[45]  H. Kelsey,et al.  Late Quaternary tectonic deformation in the Cape Arago-Bandon Region of coastal Oregon as deduced from wave-cut platforms , 1990 .

[46]  J. Steers,et al.  Holocene Sea Level Changes in the Netherlands , 1962 .

[47]  J. Clague,et al.  Tsunami deposits beneath tidal marshes on Vancouver Island, British Columbia , 1994 .

[48]  Emery I. Balazs,et al.  Geodetic evidence for aseismic subduction of the Juan de Fuca Plate , 1979 .

[49]  J. Day,et al.  Estimating shallow subsidence in microtidal salt marshes of the southeastern United States: Kaye and Barghoorn revisited , 1995 .

[50]  Hiromichi Tsuji,et al.  Silent fault slip following an interplate thrust earthquake at the Japan Trench , 1997, Nature.

[51]  Eric L. Geist,et al.  Source parameters controlling the generation and propagation of potential local tsunamis along the Cascadia margin , 1996 .

[52]  E M Scott,et al.  An Overview of All Three Stages of the International Radiocarbon Intercomparison , 1990, Radiocarbon.

[53]  M. Ando Source mechanisms and tectonic significance of historical earthquakes along the nankai trough, Japan , 1975 .

[54]  J. Clague,et al.  Amount of subsidence during a late Holocene earthquake—evidence from fossil tidal marsh foraminifera at Vancouver Island, west coast of Canada , 1995 .

[55]  John J. Clague,et al.  Summary of Coastal Geologic Evidence for past Great Earthquakes at the Cascadia Subduction Zone , 1995 .

[56]  Daniel M. Davis,et al.  GREAT THRUST EARTHQUAKES AND ASEISMIC SLIP ALONG THE PLATE BOUNDARY OF THE MAKRAN SUBDUCTION ZONE , 1992 .

[57]  James C. Savage,et al.  Mechanism of the Chilean Earthquakes of May 21 and 22, 1960 , 1970 .

[58]  L. Denys Diatom analysis of coastal deposits: methodological aspects , 1984 .

[59]  J. Clague,et al.  Radiocarbon evidence for extensive plate-boundary rupture about 300 years ago at the Cascadia subduction zone , 1995, Nature.

[60]  H. Kanamori,et al.  A moment magnitude scale , 1979 .

[61]  G. Carver,et al.  Late Holocene Tectonics and Paleoseismicity, Southern Cascadia Subduction Zone , 1992, Science.

[62]  Alfred C. Redfield,et al.  Development of a New England Salt Marsh , 1972 .

[63]  H. Kelsey Late Quaternary deformation of marine terraces on the Cascadia Subduction Zone near Cape Blanco, Oregon , 1990 .

[64]  James T. Morris,et al.  Emission of gaseous carbon dioxide from salt-marsh sediments and its relation to other carbon losses , 1986 .

[65]  A. Hine,et al.  Coastal storm deposition: Salt-marsh response to a severe extratropical storm, March 1993, west-central Florida , 1995 .

[66]  Kelin Wang,et al.  Thermal constraints on the zone of major thrust earthquake failure: The Cascadia Subduction Zone , 1993 .

[67]  G. Griggs,et al.  Deep-sea sedimentation and sediment-fauna interaction in Cascadia Channel and on Cascadia Abyssal Plain☆ , 1969 .

[68]  Derald G. Smith,et al.  Evidence for eight great earthquake-subsidence events detected with ground-penetrating radar, Willapa barrier, Washington , 1996 .

[69]  Daniel E. Bunker,et al.  Tree-ring dating the 1700 Cascadia earthquake , 1997, Nature.

[70]  J. C. Savage,et al.  Tide gage measurements of uplift along the south coast of Alaska , 1991 .

[71]  M. E. Mackay,et al.  Transverse structural trends along the Oregon convergent margin: Implications for Cascadia earthquake potential and crustal rotations , 1992 .

[72]  K. Walker,et al.  Tidal marsh stratigraphy, sea-level change and large earthquakes, i: a 5000 year record in washington, U.S.A.☆ , 1996 .

[73]  John J. Clague,et al.  Tsunami Deposits Beneath Tidal Marshes on Northwestern Vancouver Island, British Columbia , 1997, Quaternary Research.

[74]  Manuel G. Bonilla,et al.  Effects of the earthquake of March 27, 1964, on the Alaska Railroad: Chapter D in The Alaska earthquake, March 27, 1964: effects on transportation, communications, and utilities , 1970 .

[75]  M. Torresan A review and comparison of the hydrophotometer and pipette methods in the analysis of fine-grained sediment , 1987 .

[76]  P. Vincent,et al.  Present‐day vertical deformation of the Cascadia Margin, Pacific Northwest, United States , 1994 .

[77]  J. Clague,et al.  Optical Dating of Tsunami-Laid Sands , 1996, Quaternary Research.

[78]  High-precision decadal calibration of the radiocarbon time scale, AD 1950-2500 BC. , 1986 .

[79]  M. Darienzo,et al.  Episodic tectonic subsidence of Late Holocene salt marshes, northern Oregon Central Cascadia Margin , 1990 .

[80]  M. Stuiver,et al.  High-Precision Calibration of the Radiocarbon Time Scale, AD 1950–500 BC , 1986, Radiocarbon.

[81]  K. Berryman Age, height, and deformation of Holocene marine terraces at Mahia Peninsula, Hikurangi Subduction Margin, New Zealand , 1993 .

[82]  I. Shennan Flandrian sea‐level changes in the Fenland. I: The geographical setting and evidence of relative sea‐level changes , 1986 .

[83]  J. Stein,et al.  Earthquake-induced burial of archaeological sites along the southern Washington coast about A.D. 1700 , 1996 .

[84]  C. Lomnitz Major earthquakes and tsunamis in Chile during the period 1535 to 1955 , 1970 .

[85]  C. Goldfinger,et al.  Forearc Deformation and Great Subduction Earthquakes: Implications for Cascadia Offshore Earthquake Potential , 1995, Science.

[86]  A. Meyerson Pollen and paleosalinity analyses from a Holocene tidal marsh sequence, Cape May County, New Jersey , 1972 .

[87]  M. Stuiver,et al.  Radiocarbon test of earthquake magnitude at the Cascadia subduction zone , 1991, Nature.

[88]  David R. Brillinger,et al.  A more precise chronology of earthquakes produced by the San Andreas fault in southern California , 1989 .

[89]  C. Hitchcock,et al.  Flora of the Pacific Northwest: An Illustrated Manual , 1973 .

[90]  Kelin Wang,et al.  Current deformation and the width of the seismogenic zone of the northern Cascadia subduction thrust , 1994 .

[91]  E. Hemphill-Haley Occurrences of Recent and Holocene intertidal diatoms (Bacillariophyta) in northern Willapa Bay, Washington , 1993 .