Description: Historical records indicate frequent seismic activity along the north-east Caribbean plate boundary over the past 500 years, particularly on the island of Hispaniola. We use accounts of historical earthquakes to assign intensities and the intensity assignments for the 2010 Haiti earthquakes to derive an intensity attenuation relation for Hispaniola. The intensity assignments and the attenuation relation are used in a grid search to find source locations and magnitudes that best fit the intensity assignments.Here we describe a sequence of devastating earthquakes on the Enriquillo fault system in the eighteenth century. An intensity magnitude MI 6.6 earthquake in 1701 occurred near the location of the 2010 Haiti earthquake, and the accounts of the shaking in the 1701 earthquake are similar to those of the 2010 earthquake. A series of large earthquakes migrating from east to west started with the 18 October 1751 MI 7.4–7.5 earthquake, probably located near the eastern end of the fault in the Dominican Republic, followed by the 21 November 1751 MI 6.6 earthquake near Port-au-Prince, Haiti, and the 3 June 1770 MI 7.5 earthquake west of the 2010 earthquake rupture. The 2010 Haiti earthquake may mark the beginning of a new cycle of large earthquakes on the Enriquillo fault system after 240 years of seismic quiescence. The entire Enriquillo fault system appears to be seismically active; Haiti and the Dominican Republic should prepare for future devastating earthquakes.

Description: Author Posting. © American Geophysical Union, 2006. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geophysical Research Letters 33 (2006): L11307, doi:10.1029/2006GL026125.

Description: We have established for the first time a size frequency distribution for carbonate submarine slope failures. Using detailed bathymetry along the northern edge of the carbonate platform north of Puerto Rico, we show that the cumulative distribution of slope failure volumes follows a power-law distribution. The power-law exponent of this distribution is similar to those for rock falls on land, commensurate with their interpreted failure mode. The carbonate volume distribution and its associated volume-area relationship are significantly different from those for clay-rich debris lobes in the Storegga slide, Norway. Coupling this relationship with tsunami simulations allows an estimate of the maximum tsunami runup and the maximum number of potentially damaging tsunamis from landslides to the north shore of Puerto Rico.

Description: Author Posting. © American Geophysical Union, 2004. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research 109 (2004): B12310, doi:10.1029/2004JB003031.

Description: Strike-slip faults in the forearc region of a subduction zone often present significant seismic hazard because of their proximity to population centers. We explore the interaction between thrust events on the subduction interface and strike-slip faults within the forearc region using three-dimensional models of static Coulomb stress change. Model results reveal that subduction earthquakes with slip vectors subparallel to the trench axis enhance the Coulomb stress on strike-slip faults adjacent to the trench but reduce the stress on faults farther back in the forearc region. In contrast, subduction events with slip vectors perpendicular to the trench axis enhance the Coulomb stress on strike-slip faults farther back in the forearc, while reducing the stress adjacent to the trench. A significant contribution to Coulomb stress increase on strike-slip faults in the back region of the forearc comes from “unclamping” of the fault, i.e., reduction in normal stress due to thrust motion on the subduction interface. We argue that although Coulomb stress changes from individual subduction earthquakes are ephemeral, their cumulative effects on the pattern of lithosphere deformation in the forearc region are significant. We use the Coulomb stress models to explain the contrasting deformation pattern between two adjacent segments of the Caribbean subduction zone. Subduction earthquakes with slip vectors nearly perpendicular to the Caribbean trench axis is dominant in the Hispaniola segment, where the strike-slip faults are more than 60 km inland from the trench. In contrast, subduction slip motion is nearly parallel to the Caribbean trench axis along the Puerto Rico segment, where the strike-slip fault is less than 15 km from the trench. This observed jump from a strike-slip fault close to the trench axis in the Puerto Rico segment to the inland faults in Hispaniola is explained by different distributions of Coulomb stress in the forearc region of the two segments, as a result of the change from the nearly trench parallel slip on the Puerto Rico subduction interface to the more perpendicular subduction slip beneath Hispaniola. The observations and modeling suggest that subduction-induced strike-slip seismic hazard to Puerto Rico may be smaller than previously assumed but the hazard to Hispaniola remains high.

Description: J.L. was supported by the National Science Foundation through grant NSF-EAR0003888.

Description: This paper is not subject to U.S. copyright. The definitive version was published in Journal of Geophysical Research 116 (2011): B12318, doi:10.1029/2011JB008497.

Description: We evaluate the long-term seismic activity of the North-American/Caribbean plate boundary from 500 years of historical earthquake damage reports. The 2010 Haiti earthquakes and other earthquakes were used to derive regional attenuation relationships between earthquake intensity, magnitude, and distance from the reported damage to the epicenter, for Hispaniola and for Puerto Rico and the Virgin Islands. The attenuation relationship for Hispaniola earthquakes and northern Lesser Antilles earthquakes is similar to that for California earthquakes, indicating a relatively rapid attenuation of damage intensity with distance. Intensities in Puerto Rico and the Virgin Islands decrease less rapidly with distance. We use the intensity-magnitude relationships to systematically search for the location and intensity magnitude MI which best fit all the reported damage for historical earthquakes. Many events occurred in the 20th-century along the plate-boundary segment from central Hispaniola to the NW tip of Puerto Rico, but earlier events from this segment were not identified. The remaining plate boundary to the east to Guadeloupe is probably not associated with M 〉 8 historical subduction-zone earthquakes. The May 2, 1787 earthquake, previously assigned an M 8–8.25, is probably only MI 6.9 and could be located north, west or SW of Puerto Rico. An MI 6.9 earthquake on July 11, 1785 was probably located north or east of the Virgin Islands. We located MI 〈 8 historical earthquakes on April 5, 1690, February 8, 1843, and October 8, 1974 in the northern Lesser Antilles within the arc. We speculate that the December 2, 1562 (MI 7.7) and May 7, 1842 (MI 7.6) earthquakes ruptured the Septentrional Fault in northern Hispaniola. If so, the recurrence interval on the central Septentrional Fault is ∼300 years, and only 170 years has elapsed since the last event. The recurrence interval of large earthquakes along the Hispaniola subduction segment is likely longer than the historical record. Intra-arc M ≥ 7.0 earthquakes may occur every 75–100 years in the 410-km-long segment between the Virgin Islands and Guadeloupe.

Description: This paper is not subject to U.S. copyright. The definitive version was published in Geophysical Research Letters 39 (2012): L10304, doi:10.1029/2012GL051485.

Description: Kinematic similarities between the Sumatra and Puerto Rico Trenches highlight the potential for a mega-earthquake along the Puerto Rico Trench and the generation of local and trans-Atlantic tsunamis. We used the horizontal components of continuous GPS (cGPS) measurements from 10 sites on NE Caribbean islands to evaluate strain accumulation along the North American (NA) – Caribbean (CA) plate boundary. These sites move westward and slightly northward relative to CA interior at rates ≤2.5 mm/y. Provided this motion originates in the subduction interface, the northward motion suggests little or no trench-perpendicular thrust accumulation and may in fact indicate divergence north of Puerto Rico, where abnormal subsidence, bathymetry, and gravity are observed. The Puerto Rico Trench, thus, appears unable to generate mega-earthquakes, but damaging smaller earthquakes cannot be discounted. The westward motion, characterized by decreasing rate with distance from the trench, is probably due to eastward motion of CA plate impeded at the plate boundary by the Bahamas platform. Two additional cGPS sites in Mona Passage and SW Puerto Rico move to the SW similar to Hispaniola and unlike the other 10 sites. That motion relative to the rest of Puerto Rico may have given rise to seismicity and normal faults in Mona Rift, Mona Passage, and SW Puerto Rico.

Description: © The Author(s), 2010. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Natural Hazards 63 (2012): 51-84, doi:10.1007/s11069-010-9622-6.

Description: 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–3.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–1800 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.

Description: The work was supported in part by the Nuclear Regulatory Commission under its project N6480, a tsunami-hazard assessment for the eastern United States.

Description: © The Author(s), 2017. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Geosphere 13 (2017): 301-368, doi:10.1130/GES01356.1.

Description: Extraordinary marine inundation scattered clasts southward on the island of Anegada, 120 km south of the Puerto Rico Trench, sometime between 1200 and 1480 calibrated years (cal yr) CE. Many of these clasts were likely derived from a fringing reef and from the sandy flat that separates the reef from the island’s north shore. The scattered clasts include no fewer than 200 coral boulders, mapped herein for the first time and mainly found hundreds of meters inland. Many of these are complete colonies of the brain coral Diploria strigosa. Other coral species represented include Orbicella (formerly Montastraea) annu­laris, Porites astreoides, and Acropora palmata. Associated bioclastic carbonate sand locally contains articulated cobble-size valves of the lucine Codakia orbicularis and entire conch shells of Strombus gigas, mollusks that still inhabit the sandy shallows between the island’s north shore and a fringing reef beyond. italicmbricated limestone slabs are clustered near some of the coral boulders. In addition, fields of scattered limestone boulders and cobbles near sea level extend mainly southward from limestone sources as much as 1 km inland. Radiocarbon ages have been obtained from 27 coral clasts, 8 lucine valves, and 3 conch shells. All these additional ages predate 1500 cal yr CE, all but 2 are in the range 1000–1500 cal yr CE, and 16 of 22 brain coral ages cluster in the range 1200–1480 cal yr CE. The event marked by these coral and mollusk clasts likely occurred in the last centuries before Columbus (before 1492 CE). The pre-Columbian deposits surpass Anegada’s previously reported evidence for extreme waves in post-Columbian time. The coarsest of the modern storm deposits consist of coral rubble that lines the north shore and sandy fans on the south shore; neither of these storm deposits extends more than 50 m inland. More extensive overwash, perhaps by the 1755 Lisbon tsunami, is marked primarily by a sheet of sand and shells found mainly below sea level beneath the floors of modern salt ponds. This sheet extends more than 1 km southward from the north shore and dates to the interval 1650–1800 cal yr CE. Unlike the pre-Columbian deposits, it lacks coarse clasts from the reef or reef flat; its shell assemblage is instead dominated by cerithid gastropods that were merely stirred up from a marine pond in the island’s interior. In their inland extent and clustered pre-Columbian ages, the coral clasts and associated deposits suggest extreme waves unrivaled in recent millennia at Anegada. Bioclastic sand coats limestone 4 m above sea level in areas 0.7 and 1.3 km from the north shore. A coral boulder of nearly 1 m3 is 3 km from the north shore by way of an unvegetated path near sea level. As currently understood, the extreme flooding evidenced by these and other clasts represents either an extraordinary storm or a tsunami of nearby origin. The storm would need to have produced tsunami-like bores similar to those of 2013 Typhoon Haiyan in the Philippines. Normal faults and a thrust fault provide nearby tsunami sources along the eastern Puerto Rico Trench.

Description: Published 2018. This article is a U.S. Government work and is in the public domain in the USA. The definitive version was published in Journal of Geophysical Research: Solid Earth 123 (2018): 4223-4242, doi:10.1002/2017JB014770.

Description: The Pacific/North America (PA/NA) plate boundary between Vancouver Island and Alaska is similar to the PA/NA boundary in California in its kinematic history and the rate and azimuth of current relative motion, yet their deformation styles are distinct. The California plate boundary shows a broad zone of parallel strike slip and thrust faults and folds, whereas the 49‐mm/yr PA/NA relative plate motion in Canada and Alaska is centered on a single, narrow, continuous ~900‐km‐long fault, the Queen Charlotte Fault (QCF). Using gravity analysis, we propose that this plate boundary is centered on the continent/ocean boundary (COB), an unusual location for continental transform faults because plate boundaries typically localize within the continental lithosphere, which is weaker. Because the COB is a boundary between materials of contrasting elastic properties, once a fault is established there, it will probably remain stable. We propose that deformation progressively shifted to the COB in the wake of Yakutat terrane's northward motion along the margin. Minor convergence across the plate boundary is probably accommodated by fault reactivation on Pacific crust and by an eastward dipping QCF. Underthrusting of Pacific slab under Haida Gwaii occurs at convergence angles 〉14°–15° and may have been responsible for the emergence of the archipelago. The calculated slab entry dip (5°–8°) suggests that the slab probably does not extend into the asthenosphere. The PA/NA plate boundary at the QCF can serve as a structurally simple site to investigate the impact of rheology and composition on crustal deformation and the initiation of slab underthrusting.

Description: Author Posting. © American Geophysical Union, 2005. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Eos 86 (2005): 90, doi:10.1029/2005EO090004.

Description: RayGUI 2.0 is a new version of RayGUI, a graphical user interface (GUI) to the seismic travel time modeling program of Zelt and Smith [1992]. It represents a significant improvement over the previous version of RayGUI (RayGUI 1.04; Loss et al.[1998a,1998b]).