Reconstructing hydrodynamic flow parameters of the 1700 tsunami at Cannon Beach, Oregon, USA

Coastal communities in the western United States face risks of inundation by distant tsunamis that propagate across the Pacific Ocean as well as local tsunamis produced by great (M_w?>?8) earthquakes on the Cascadia subduction zone. In 1964, the M_w 9.2 Alaska earthquake launched a Pacific-wide tsunami that flooded Cannon Beach, a small community (population 1640) in northwestern Oregon, causing over $230,000 in damages. However, since the giant 2004 Indian Ocean tsunami, the 2010 Chile tsunami and the recent 2011 Tohoku-Oki tsunami, renewed concern over potential impacts of a Cascadia tsunami on the western US has motivated closer examination of the local hazard. This study applies a simple sediment transport model to reconstruct the flow speed of the most recent Cascadia tsunami that flooded the region in 1700 using the thickness and grain size of sand layers deposited by the waves. Sedimentary properties of sand from the 1700 tsunami deposit provide model inputs. The sediment transport model calculates tsunami flow speed from the shear velocity required to suspend the quantity and grain size distribution of the observed sand layers. The model assumes a steady, spatially uniform tsunami flow and that sand settles out of suspension forming a deposit when the flow velocity decreases to zero. Using flow depths constrained by numerical tsunami simulations for Cannon Beach, the sediment transport model calculated flow speeds of 6.5?C7.6?m/s for sites within 0.6?km of the beach and higher flow speeds (~8.8?m/s) for sites 0.8?C1.2?km inland. Flow speed calculated for sites within 0.6?km of the beach compare well with maximum velocities estimated for the largest tsunami simulation. The higher flow speeds calculated for the two sites furthest landward contrast with much lower maximum velocities (<3.8?m/s) predicted by numerical simulations. Grain size distributions of sand layers from the most distal sites are inconsistent with deposition from sediment falling out of suspension. We infer that rapid deceleration in tsunami flow and convergences in sediment transport formed unusually thick deposits. Consequently, higher flow speeds calculated by the sediment model probably overestimate the actual wave speed at sites furthest inland.     

Coastal subsidence in Oregon,USA, during the giant Cascadia earthquake of AD 1700

Quantitative estimates of land-level change during the giant AD 1700 Cascadia earthquake along the Oregon coast are inferred from relative sea-level changes reconstructed from fossil foraminiferal assemblages preserved within the stratigraphic record. A transfer function, based upon a regional training set of modern sediment samples from Oregon estuaries, is calibrated to fossil assemblages in sequences of samples across buried peat-mud and peat-sand contacts marking the AD 1700 earthquake. Reconstructions of sample elevations with sample-specific errors estimate the amount of coastal subsidence during the earthquake at six sites along 400 km of coast. The elevation estimates are supported by lithological, carbon isotope, and faunal tidal zonation data. Coseismic subsidence at Nehalem River, Nestucca River, Salmon River, Alsea Bay, Siuslaw River and South Slough varies between 0.18 m and 0.85 m with errors between 0.18 m and 0.32 m. These subsidence estimates are more precise, consistent, and generally lower than previous semi-quantitative estimates. Following earlier comparisons of semi-quantitative subsidence estimates with elastic dislocation models of megathrust rupture during great earthquakes, our lower estimates for central and northern Oregon are consistent with modeled rates of strain accumulation and amounts of slip on the subduction megathrust, and thus, with a magnitude of 9 for the AD 1700 earthquake.     

Recurrence intervals of major paleotsunamis as calibrated by historic tsunami deposits in three localities: Port Alberni, Cannon Beach, and Crescent City, along the Cascadia margin, Canada and USA

We review geologic records of both historic and prehistoric tsunami inundations at three widely separated localities that experienced significant damage from the 1964 Alaskan tsunami along the Cascadia margin. The three localities are Port Alberni, Cannon Beach, and Crescent City, representing, respectively, the north, central, and south portions of the study area (1,000 km in length). The geologic records include anomalous sand sheets from marine surges that are hosted in supratidal peaty mud deposits. Paleotsunami sand sheets that exceed the thickness, continuity and/or extent of the 1964 historic tsunami are counted as major paleotsunami inundations. Major paleotsunamis (6–7 in number) at each locality during the last 3,000 years demonstrate mean recurrence intervals of 450–540 years, and within-cluster intervals (three events each) of 270–460 years. It has been 313 years since the last major paleotsunami from a great Cascadia earthquake in AD 1700. We compare the dated sequences of major paleotsunami inundations to the nearest regional records of coastal coseismic subsidence in Willapa Bay in the central margin, Waatch/Neah Bay in the northern margin, and Coquille in the southern margin. Similar numbers of events from both types of records suggest that the major paleotsunamis are locally derived (near-field) from ruptures of the Cascadia margin megathrust fault zone, rather than from transoceanic tsunamis (far-field) originating at other subduction zones around the Pacific Rim. Given the catastrophic hazard of the near-field Cascadia margin tsunamis, we propose a basic rule for reminding the general public of the need for self-initiated evacuation following a great earthquake at the Cascadia margin.     

Confidence levels for tsunami-inundation limits in northern Oregon inferred from a 10,000-year history of great earthquakes at the Cascadia subduction zone

To explore the local tsunami hazard from the Cascadia subduction zone we (1) evaluate geologically reasonable variability of the earthquake rupture process, (2) specify 25 deterministic earthquake sources, and (3) use resulting vertical coseismic deformations for simulation of tsunami inundation at Cannon Beach, Oregon. Maximum runup was 9–30 m (NAVD88) from earthquakes with slip of ~8–38 m and M _w ~8.3–9.4. Minimum subduction zone slip consistent with three tsunami deposits was 14–15 m. By assigning variable weights to the source scenarios using a logic tree, we derived percentile inundation lines that express the confidence level (percentage) that a Cascadia tsunami will not exceed the line. Ninety-nine percent of Cascadia tsunami variation is covered by runup ≤30 m and 90% ≤16 m with a “preferred” (highest weight) value of ~10 m. A hypothetical maximum-considered distant tsunami had runup of ~11 m, while the historical maximum was ~6.5 m.     

Great earthquakes of variable magnitude at the Cascadia subduction zone

Comparison of histories of great earthquakes and accompanying tsunamis at eight coastal sites suggests plate-boundary ruptures of varying length, implying great earthquakes of variable magnitude at the Cascadia subduction zone. Inference of rupture length relies on degree of overlap on radiocarbon age ranges for earthquakes and tsunamis, and relative amounts of coseismic subsidence and heights of tsunamis. Written records of a tsunami in Japan provide the most conclusive evidence for rupture of much of the plate boundary during the earthquake of 26 January 1700. Cascadia stratigraphic evidence dating from about 1600 cal yr B.P., similar to that for the 1700 earthquake, implies a similarly long rupture with substantial subsidence and a high tsunami. Correlations are consistent with other long ruptures about 1350 cal yr B.P., 2500 cal yr B.P., 3400 cal yr B.P., 3800 cal yr B.P., 4400 cal yr B.P., and 4900 cal yr B.P. A rupture about 700-1100 cal yr B.P. was limited to the northern and central parts of the subduction zone, and a northern rupture about 2900 cal yr B.P. may have been similarly limited. Times of probable short ruptures in southern Cascadia include about 1100 cal yr B.P., 1700 cal yr B.P., 3200 cal yr B.P., 4200 cal yr B.P., 4600 cal yr B.P., and 4700 cal yr B.P. Rupture patterns suggest that the plate boundary in northern Cascadia usually breaks in long ruptures during the greatest earthquakes. Ruptures in southernmost Cascadia vary in length and recurrence intervals more than ruptures in northern Cascadia.     

Relative sea-level change inferred from intertidal sediments beneath marshes on Vancouver Island, British Columbia

Relative sea-level change at the time of, and since, the most recent great earthquake at the Cascadia subduction zone is estimated from intertidal sediments at three marshes on western Vancouver Island, British Columbia. We compare the elevation of the pre-earthquake surface, which is marked by a tsunami sand sheet, with the modern depositional elevation range of the sediment type upon which the sand was deposited. At a site south of the Nootka fault zone, which is the northern boundary of the subducting Juan de Fuca plate, tidal mud overlies the pre-earthquake marsh surface. The stratigraphy at this site indicates 0.2–1.6 m of coseismic submergence and 1.1 m of subsequent emergence. In contrast, two sites to the north lack obvious stratigraphic evidence for coseismic land-level change and record between 0.1 and 1.7 m of post-earthquake submergence. These results indicate a difference in tectonic environment across the Nootka fault zone and suggest that plate-boundary rupture during the last great Cascadia earthquake probably did not extend north of central Vancouver Island.     

GPR studies over the tsunami affected Karaikal beach,Tamil Nadu,south India

In this study, results of GPR profiling related to mapping of subsurface sedimentary layers at tsunami affected Karaikal beach are presented . A 400 MHz antenna was used for profiling along 262 m stretch of transect from beach to backshore areas with penetration of about 2.0 m depth (50 ns two-way travel time). The velocity analysis was carried out to estimate the depth information along the GPR profile. Based on the significant changes in the reflection amplitude, three different zones are marked and the upper zone is noticed with less moisture compared to other two (saturated) zones. The water table is noticed to vary from 0.5 to 0.75 m depth (12–15 ns) as moving away from the coastline. Buried erosional surface is observed at 1.5 m depth (40–42 ns), which represents the limit up to which the extreme event acted upon. In other words, it is the depth to which the tsunami sediments have been piled up to about 1.5 m thickness. Three field test pits were made along the transect and sedimentary sequences were recorded. The sand layers, especially, heavy mineral layers, recorded in the test pits indicate a positive correlation with the amplitude and velocity changes in the GPR profile. Such interpretation seems to be difficult in the middle zone due to its water saturation condition. But it is fairly clear in the lower zone located just below the erosional surface where the strata is comparatively more compact. The inferences from the GPR profile thus provide a lucid insight to the subsurface sediment sequences of the tsunami sediments in the Karaikal beach.     

Sedimentology of the December 26, 2004, Sumatra tsunami deposits in eastern India (Tamil Nadu) and Kenya

The December 26, 2004 Sumatra tsunami caused severe damage at the coasts of the Indian ocean. We report results of a sedimentological study of tsunami run-up parameters and the sediments laid down by the tsunami at the coast of Tamil Nadu, India, and between Malindi and Lamu, Kenya. In India, evidence of three tsunami waves is preserved on the beaches in the form of characteristic debris accumulations. We measured the maximum run-up distance at 580 m and the maximum run-up height at 4.85 m. Flow depth over land was at least 3.5 m. The tsunami deposited an up to 30 cm thick blanket of moderately well to well-sorted coarse and medium sand that overlies older beach deposits or soil with an erosional unconformity. The sand sheet thins inland without a decrease of grain-size. The deposits consist frequently of three layers. The lower one may be cross-bedded with foresets dipping landward and indicating deposition during run-up. The overlying two sand layers are graded or parallel-laminated without indicators of current directions. Thus, it remains undecided whether they formed during run-up or return flow. Thin dark laminae rich in heavy minerals frequently mark the contacts between successive layers. Benthic foraminifera indicate an entrainment of sediment by the tsunami from water depths less than ca. 30 m water depth. On the Indian shelf these depths are present at distances of up to 5 km from the coast. In Kenya only one wave is recorded, which attained a run-up height of 3 m at a run-up distance of ca. 35 m from the tidal water line at the time of the tsunami impact. Only one layer of fine sand was deposited by the tsunami. It consists predominantly of heavy minerals supplied to the sea by a nearby river. The sand layer thins landward with a minor decrease in grain-size. Benthic foraminifera indicate an entrainment of sediment by the tsunami from water depths less than ca. 30 m water depth, reaching down potentially to ca. 80 m. The presence of only one tsunami-related sediment layer in Kenya, but three in India, reflects the impact of only one wave at the coast of Kenya, as opposed to several in India. Grain-size distributions in the Indian and Kenyan deposits are mostly normal to slightly positively skewed and indicate that the detritus was entrained by the tsunami from well sorted pre-tsunami deposits in nearshore, swash zone and beach environments.     

Tsunami impact to Washington and northern Oregon from segment ruptures on the southern Cascadia subduction zone

This paper explores the size and arrival of tsunamis in Oregon and Washington from the most likely partial ruptures of the Cascadia subduction zone (CSZ) in order to determine (1) how quickly tsunami height declines away from sources, (2) evacuation time before significant inundation, and (3) extent of felt shaking that would trigger evacuation. According to interpretations of offshore turbidite deposits, the most frequent partial ruptures are of the southern CSZ. Combined recurrence of ruptures extending ~490 km from Cape Mendocino, California, to Waldport, Oregon (segment C) and ~320 km from Cape Mendocino to Cape Blanco, Oregon (segment D), is ~530 years. This recurrence is similar to frequency of full-margin ruptures on the CSZ inferred from paleoseismic data and to frequency of the largest distant tsunami sources threatening Washington and Oregon, ~M _w 9.2 earthquakes from the Gulf of Alaska. Simulated segment C and D ruptures produce relatively low-amplitude tsunamis north of source areas, even for extreme (20 m) peak slip on segment C. More than ~70 km north of segments C and D, the first tsunami arrival at the 10-m water depth has an amplitude of <1.9 m. The largest waves are trapped edge waves with amplitude ≤4.2 m that arrive ≥2 h after the earthquake. MM V–VI shaking could trigger evacuation of educated populaces as far north as Newport, Oregon for segment D events and Grays Harbor, Washington for segment C events. The NOAA and local warning systems will be the only warning at greater distances from sources.     

Luminescence dating of buried cobble surfaces from sandy beach ridges: a case study from Denmark

Here we investigate the use of optically stimulated luminescence (OSL) for dating cobbles from the body of successive beach ridges and compare cobble surface‐derived ages to standard quartz OSL ages from sand. Between four and eight cobbles and sand samples (age control) were dated with the luminescence method, taken from the modern beach and from beach ridges on the south and north extremes of a prograding spit on the westernmost coast of Lolland, Denmark. Luminescence‐depth profiles perpendicular to the surfaces of the cobbles show that the feldspar infrared signals stimulated at 50 °C were fully reset to various depths into the cobbles prior to final deposition; as a result, the equivalent doses determined from close to the surface of such cobbles can be used to calculate burial ages. Beach‐ridge burial ages given by the average of ages of individual cobbles taken from the same site are consistent, within errors, with the ages derived from the sand samples. Cobble‐ and sand‐derived ages show that the southernmost beach ridge at Albuen was formed around 2 ka ago, indicating that this sandy spit is younger than other coastal systems in Denmark. The agreement between ages derived from clasts and from standard quartz OSL in this study confirms that, even in the absence of sandy sediments, we can reliably date sites using OSL by targeting larger clasts. In addition, the record of prior light exposure contained in the shape of the cobbles’ luminescence‐depth profile removes one of the major uncertainties (i.e. the degree of signal reset prior to burial) in the luminescence dating of high latitude sites.     

Occurrence of soft sediment deformation at Dive Agar beach, west coast of India: possible record of the Indian Ocean tsunami (2004)

We describe here a sequence of soft sediment deformation (SSD) structures at Dive Agar beach near Srivardhan in the west coast of India. The ~120-cm-thick sediment package is represented by a basal undeformed sand (layer A) sharply cut by ~30-cm-thick intermixed beach sand and terrigenous sand (layer B1) followed by complex load structures and convolutions (8?C15?cm) within a coarse sandy layer (B2). The layer B2 is scoured by terrigenous sand (layer C1) which is capped with a silty mud layer (C2). The entire sequence (B2?CC1?CC2) is intruded by sand dykes originating from the lower layer B1. This sediment package is further overlain by a heavy mineral reach marine sand (layer D) with liquefactions long axes inclined southward as a result of forceful long-shore drift. The profile ends up with coarse-grained, poorly sorted sand including angular clasts of terrigenous outwash deposits indicating return of distal inundations. Intense deformation (liquefaction) is restricted to the heavy mineral-rich marine and the intermixed sands (layers B2 and D), whereas the terrigenous sand layers show scoured bases with oscillatory and pebbly tops. The presence of complex load structures injecting into the underlying layers, the top-truncated sand dykes, macro-thrust faults, scouring, and inclusion of coral fragments can explain it as a record of tsunami in the west coast. Occurrence of un-decayed consumer plastic material within the deformed layers suggests it as one of the most recent tsunami events (i.e., 2004 IOT), the only reported event after 1945 in the west coast. Alternative marine and terrigenous sands are characteristic of tsunami run-up and backwash deposits, while the dimensions of SSDs may be related to the <2?m magnitude (height) of the 2004 IOT at Dive Agar.     

Indirect economic losses of drought under future projections of climate change: a case study for Spain

This study presents the results of numerical simulations of the 2004 Indian Ocean earthquake and tsunami in the Bay of Lhok Nga (northwestern coast of Sumatra, Indonesia) integrating sediment erosion and deposition. We investigate the transport of sediment both by suspension and by bedload under different scenarii of long breaking dispersive waves through a series of numerical experiments. The earthquake source model used by Koshimura et al. (Coast Eng J 51:243–273, 2008) with a 25-m dislocation better reproduces the wave travel time, flow depth and inundation area than the other models tested. The model reproduces realistically the pronounced coastal retreat in the northern part of Lhok Nga Bay (retreat ranging between 50 and 150 m), where Paris et al. (Geomorphology 104:59–72, 2009) estimated a mean retreat of 80 m. There is also a good agreement between the simulated area of coastal retreat (195,400 m²) and the field observations (203,200 m²). The simulation may underestimate the volume of tsunami deposits (611,700 m³ vs. 500,000–1,000,000 m³ estimated by Paris et al. (2009). The model fully reproduces the observed thickness of tsunami deposits when considering both bedload and suspension, even if bedload transport dominates. Limitations are due to micro-scale topographic, anthropic features (which are not always represented by the DEM) and the amount of debris which may influence flow dynamics and sediment transport.     

Tsunami hazard assessment of Canada

We present a preliminary probabilistic tsunami hazard assessment of Canadian coastlines from local and far-field, earthquake, and large submarine landslide sources. Analyses involve published historical, palaeotsunami and palaeoseismic data, modelling, and empirical relations between fault area, earthquake magnitude, and tsunami run-up. The cumulative estimated tsunami hazard for potentially damaging run-up (≥1.5 m) of the outer Pacific coastline is ~40–80 % in 50 years, respectively one and two orders of magnitude greater than the outer Atlantic (~1–15 %) and the Arctic (<1 %). For larger run-up with significant damage potential (≥3 m), Pacific hazard is ~10–30 % in 50 years, again much larger than both the Atlantic (~1–5 %) and Arctic (<1 %). For outer Pacific coastlines, the ≥1.5 m run-up hazard is dominated by far-field subduction zones, but the probability of run-up ≥3 m is highest for local megathrust sources, particularly the Cascadia subduction zone; thrust sources further north are also significant, as illustrated by the 2012 Haida Gwaii event. For Juan de Fuca and Georgia Straits, the Cascadia megathrust dominates the hazard at both levels. Tsunami hazard on the Atlantic coastline is dominated by poorly constrained far-field subduction sources; a lesser hazard is posed by near-field continental slope failures similar to the 1929 Grand Banks event. Tsunami hazard on the Arctic coastline is poorly constrained, but is likely dominated by continental slope failures; a hypothetical earthquake source beneath the Mackenzie delta requires further study. We highlight areas susceptible to locally damaging landslide-generated tsunamis, but do not quantify the hazard.     

Distribution and sedimentary characteristics of tsunami deposits along the Cascadia margin of western North America

Tsunami deposits have been found at more than 60 sites along the Cascadia margin of Western North America, and here we review and synthesize their distribution and sedimentary characteristics based on the published record. Cascadia tsunami deposits are best preserved, and most easily identified, in low-energy coastal environments such as tidal marshes, back-barrier marshes and coastal lakes where they occur as anomalous layers of sand within peat and mud. They extend up to a kilometer inland in open coastal settings and several kilometers up river valleys. They are distinguished from other sediments by a combination of sedimentary character and stratigraphic context. Recurrence intervals range from 300–1000 years with an average of 500–600 years. The tsunami deposits have been used to help evaluate and mitigate tsunami hazards in Cascadia. They show that the Cascadia subduction zone is prone to great earthquakes that generate large tsunamis. The inclusion of tsunami deposits on inundation maps, used in conjunction with results from inundation models, allows a more accurate assessment of areas subject to tsunami inundation. The application of sediment transport models can help estimate tsunami flow velocity and wave height, parameters which are necessary to help establish evacuation routes and plan development in tsunami prone areas.     

The Most Recent Cascadia Earthquake and Native American Narratives

The Cascadia subduction zone fault lies just off the Pacific coast of the USA and Canada. Although this fault has been seismically inactive over the written history of the Cascadia region, it has the potential to produce catastrophic earthquakes and tsunamis. A variety of dating methods have been used to show that the most recent Cascadia earthquake occurred in 1700. Among these methods is an informal analysis of oral traditions handed down by Native American peoples that appear to refer to a major earthquake in this region. A central difficulty in analyzing these narratives quantitatively is their use of a generation and other qualitative measures of time that have no fixed lengths. Here, these narratives are analyzed under an explicit statistical model of the lengths of these measures. The results raise a question about the previous conclusion that these narratives all refer to the most recent Cascadia earthquake.

     

Storegga tsunami sand in peat below the Tapes beach ridge at Harøy, western Norway, and its possible relation to an early Stone Age settlement

One of the early problems with the Storegga tsunami deposit was how to distinguish it from deposits of the midHolocene (Tapes) transgression. An excavation on Harøy, an island on the outermost western coast of Norway, shows a distinct, clean sand bed embedded in peat and clearly separated from the overlying Tapes beach deposits. This sand bed continues in the peat landwards of the beach ridge for at least 60 m. Radiocarbon dates of the peat show that the sand was deposited some time between 6900 and 7700 yr BP. The sedimentary structures of the bed, the ¹⁴C dates, and the fact that this is the only sand bed in the peat, suggest that the sand bed was deposited by a short-lived event, the Storegga tsunami. On the neighbouring island, Fjørtoft, a Stone Age settlement, dated to 7500 yr BP, was discovered in the early 1970s. The settlement was found underneath a sand bed that later had been covered by the Tapes beach ridge deposits. When discovered, the sand covering the settlement was inferred as eolian sand. However, this investigation shows that the Storegga tsunami deposited a widespread sand bed on the land surface around this time with a similar grain size distribution to eolian sand. It is therefore suggested that the sand bed covering this settlement was deposited from the Storegga tsunami. Both the stratigraphy and ¹⁴C dates demonstrate that the Tapes transgression maximum was reached well after the Storegga tsunami on Harøy, between 6500 and 6100 yr BP.     

A Late Pleistocene coarse-grained spit-platform sequence in northern Jylland, Denmark

Several well-preserved Late Pleistocene spit systems occur uplifted in northern Jylland, Denmark. Their present-day morphological expression allows detailed study of spit growth patterns while the internal sedimentological organisation can be examined in a series of pits distributed along the length of the spits. Two characteristic vertical sequences are recognized in the systems. The first (Sequence I) consists of a giant-scale cross-bedded foreset unit, overlain by topset and beach units, while the second (Sequence II) consists of the foreset unit overlain by bar-trough and beach units. The two sequence types pass laterally into each other with a short overlap zone. They can be interpreted in terms of Meistrell's (1966, 1972) model for spit-platform growth based on scaled wave tank experiments. The giant-scale cross-bedded unit corresponds to prograding of a coarse-grained subaqueous spit-platform while the topset, bar-trough and beach units reflect the growth of the subaerial spit. The alternation between sequence I and II reflects the inversely related growth of the spit and platform structures: when the rate of subaqueous platform progradation declines, the subaerial spit grows uniformly, and when the platform progrades uniformly spit growth declines. The model is probably only valid for relatively coarse-grained systems because only these deposits would have a relatively steep front. The water depth in which the spit system progrades and thus bottom topography, determines the thickness of the giant-scale cross-bedded foreset unit because the water depth over the top of the platform is relatively constant. If the water is less than a few metres deep the spit-platform is not developed as seen where the Late Pleistocene spit systems prograded over elevations of the sea bottom. Conversely, the correct recognition of spit-platform sequences allows precise determination of sea-level and water depth at the time of formation. Finally, the model adds one further mode of formation of giant-scale cross-bedding to those already known from fluvial transverse, lateral and point bars, subtidal sand waves and Gilbert deltas.     

Brisbane Airport: an alluvial landscape veiled by marine sediments

At Brisbane Airport, the construction of a diversion channel for Kedron Brook exposed a former beach, low cliff and sand spit, which, with their associated sediments and acid sulfate soils, demonstrate a postglacial high sea-level 1.3 – 1.4 m above present mean sea-level. The beach appears to date from 4000 to 5000 y BP. It varies in level where it lies above soft ground; these variations, and sag depressions that follow buried streamlines, indicate sediment consolidation since withdrawal of the sea from the former shore. Most of the area consists of former estuarine deposits, mangrove and saline marshes, and stranded tidal flats on which acid sulfate soils are widely developed. The modern landforms mostly reproduce subsurface features, to the extent that the surface relief replicates the landscape transgressed by the sea 7000 years ago. A small rise of sea-level possibly to +0.65 m occurred about 2000 – 3000 years ago. Foredunes near the present shore that are related to a slightly lower level 1000 – 500 years ago (?0.25 m) are currently subject to wave erosion.     

Geomorphic and stratigraphic evidence for an unusual tsunami or storm a few centuries ago at Anegada, British Virgin Islands

Waters from the Atlantic Ocean washed southward across parts of Anegada, east-northeast of Puerto Rico, during a singular event a few centuries ago. The overwash, after crossing a fringing coral reef and 1.5?km of shallow subtidal flats, cut dozens of breaches through sandy beach ridges, deposited a sheet of sand and shell capped with lime mud, and created inland fields of cobbles and boulders. Most of the breaches extend tens to hundreds of meters perpendicular to a 2-km stretch of Anegada??s windward shore. Remnants of the breached ridges stand 3?m above modern sea level, and ridges seaward of the breaches rise 2.2?C3.0?m high. The overwash probably exceeded those heights when cutting the breaches by overtopping and incision of the beach ridges. Much of the sand-and-shell sheet contains pink bioclastic sand that resembles, in grain size and composition, the sand of the breached ridges. This sand extends as much as 1.5?km to the south of the breached ridges. It tapers southward from a maximum thickness of 40?cm, decreases in estimated mean grain size from medium sand to very fine sand, and contains mud laminae in the south. The sand-and-shell sheet also contains mollusks??cerithid gastropods and the bivalve Anomalocardia??and angular limestone granules and pebbles. The mollusk shells and the lime-mud cap were probably derived from a marine pond that occupied much of Anegada??s interior at the time of overwash. The boulders and cobbles, nearly all composed of limestone, form fields that extend many tens of meters generally southward from limestone outcrops as much as 0.8?km from the nearest shore. Soon after the inferred overwash, the marine pond was replaced by hypersaline ponds that produce microbial mats and evaporite crusts. This environmental change, which has yet to be reversed, required restriction of a former inlet or inlets, the location of which was probably on the island??s south (lee) side. The inferred overwash may have caused restriction directly by washing sand into former inlets, or indirectly by reducing the tidal prism or supplying sand to post-overwash currents and waves. The overwash happened after A.D. 1650 if coeval with radiocarbon-dated leaves in the mud cap, and it probably happened before human settlement in the last decades of the 1700s. A prior overwash event is implied by an inland set of breaches. Hypothetically, the overwash in 1650?C1800 resulted from the Antilles tsunami of 1690, the transatlantic Lisbon tsunami of 1755, a local tsunami not previously documented, or a storm whose effects exceeded those of Hurricane Donna, which was probably at category 3 as its eye passed 15?km to Anegada??s south in 1960.     

Luminescence dating of Holocene spit deposits: An example from Skagen Odde, Denmark

Clemmensen, L. B. & Murray, A. S. 2009: Luminescence dating of Holocene spit deposits: An example from Skagen Odde, Denmark. Boreas , 10.1111/j.1502-3885.2009.00110.x. ISSN 0300-9483.
Skagen Odde is a large, active spit system in northern Denmark that started to form about 7200 years ago. Models for spit growth have usually relied on radiocarbon-dating of swale peat (Martørv). In this study, we date the spit deposits at three sites directly using Optically Stimulated Luminescence (OSL) to obtain supplementary age control on spit development. The spit deposits consist of a lowermost succession of shoreface, beach and backshore aeolian deposits topped by a swale peat and followed by an uppermost succession of aeolian sand sheet and dune deposits. The ages of the shallow marine, beach and backshore aeolian deposits at the main study site are indistinguishable, implying good resetting of the shallow marine deposits; the average age of 4640±250 years compares well with earlier model predictions based on radiocarbon-dating of swale peat. Aeolian sand extracted from the uppermost part of the swale peat at this site provides OSL ages of between 1600 and 2500 years, in good agreement with a calibrated AMS age from the same level of 2330–2200 years. The uppermost aeolian succession consists of two units separated by a thin palaeosol, and the aeolian units have OSL ages of about 1500 years and younger than 130 years. Lowermost spit deposits at the two additional sites have average ages of 5010±240 and 3730±190, respectively, supporting the existing chronology for spit growth based on radiocarbon-dating.