Geology topics

Uri S. ten Brink

Publications and source records attributed to Uri S. ten Brink.

At least 19 recordsLinked to original sources

USGS2026

Fault sources of the 2020 southern Puerto Rico seismic sequence and the regional tectonic structure

The 2020 southwestern Puerto Rico seismic sequence occurred in an area not known for seismic activity and included multiple normal and strike‐slip focal mechanisms. The largest earthquake, M w 6.4, caused significant damage and permanent shoreline subsidence. The locations of the ruptured faults have been explored extensively using seismological and geodetic evidence, but no consensus has been reached. Here, we report the results of a high‐resolution seismic reflection survey that identified a north‐verging 20 km long curved normal fault spanning the upper reach of the Guayanilla canyon, which we call the Guayanilla Canyon fault. Part of the fault likely ruptured during the M w 6.4. The concave‐to‐the‐north fault bounds a south‐dipping ∼1 km deep asymmetric basin. This long‐lived structure, identified also on sediment isopach and magnetic anomaly maps, appears to control the formation of the Guayanilla canyon system, the only submarine drainage system in southern Puerto Rico, suggesting a relationship among seismic activity, geomorphology, and reef development. Our seismic reflection data do not support the existence of previously published major seismogenic faults in the area, such as the north‐northeast–south‐southwest Guayanilla fault zone, the Northern Boqueron Bay‐Punta Montalva fault, or a northwest–southeast strike‐slip fault, suggested to have ruptured as part of the M w 6.4 earthquake. We identified smaller faults in the survey area, some of which can be correlated with seafloor expressions, and/or with changes in sediment thickness and magnetic anomalies, and/or with moderate‐sized earthquakes. Multiple seismic images of Investigator fault, a prominent mid‐slope feature parallel to the south coast of Puerto Rico, show buried low‐angle thrusts and an abrupt change in the attitude of the overlying sediments, leading us to interpret this feature as the expression of a boundary between the now‐inactive upper accretionary prism of Muertos trough and a buried back‐arc basin that overlies the arc basement of Puerto Rico.

Puerto Rico↗

USGS2026

Rapid cooling of magma-rich margins exemplified by the Eastern North American Margin

Thermal subsidence allows sediment accumulation in passive continental margins, modulating basin development, lithospheric dynamics and Earth-system processes. Classic thermal cooling models underestimate the lithospheric subsidence needed to explain thick early post-rift sedimentary sequences in magma-rich margins. Here we isolate the processes driving this anomalous subsidence at the Eastern North American Margin, a type-example for magma-rich margins. We construct an analytical model simultaneously accounting for non-uniform lithospheric extension, magmatic additions, and variable thermal diffusivity, which controls the rate at which heat spreads through rock. We then calibrate the model parameters against observed sediment and crustal thicknesses across the Baltimore Canyon Trough. Anomalous subsidence can be reconciled if early post-rift lithospheric heat transport is up to 60% faster than typically assumed, particularly at the distal margin edge. We propose that this increase represents hydrothermal fluid circulation enhancing lithospheric cooling during early post-rift periods, offering a new perspective on continental margin evolution.

Eastern North American Margin↗

USGS2026

Erratum to earthquake recurrence estimates for northern Caribbean faults from combinatorial optimization

In the original article (ten Brink and Geist, 2025), we derived on‐fault magnitude–frequency relationships with uncertainties from a combinatorial optimization method that assigns individual earthquakes from a regional catalog to the different faults. Our model used a magnitude versus area relationship from Shaw (2023). His relationship was published with an error, M=f(1/√A)⁠, instead of M=f(1/A) (Shaw, 2025). Here, we recalculate the magnitude–frequency relationship with uncertainties using the Shaw (2025) corrected formula. The remaining parameters of the model, the model itself, and the number of repeated optimization runs have not changed. Figures 1 and 2 show the corrected results and replace figures 3 and 4 in the original article.

The Seismic Record↗

USGS2025

Crustal structure along and surrounding the Hawaiian Islands: Volcanic island construction across scarred oceanic crust

The Hawaiian Ridge is a classic example of an intraplate volcanic island chain emplaced on oceanic lithosphere. We seek to constrain both the deformation from island loading around the Hawaiian Ridge and the influence of the oceanic lithosphere, including the Molokaʻi fracture zone (MFZ), on off-axis volcanic emplacement. To examine these processes, we conducted a marine geophysical experiment in 2018 that included the acquisition of eight multichannel seismic reflection lines and used a 6,600 in 3 tuned air gun array and an ultra-long hydrophone streamer cable towed behind R/V Marcus G. Langseth across and around the southern Hawaiian Ridge. We image both the top of igneous oceanic crust and the Moho, and we observe significant variations in regional crustal structure and thickness variations, Moho characteristics, the locations of the buried MFZ, and the flexure of the Pacific oceanic lithosphere under the Hawaiian Ridge. We observe up to ∼4.5 km of deflection of the igneous oceanic crust in response to the volcanic load with sediment thickness increasing to ∼3–3.4 km near the ridge but not correlating with the deflection. A systematic difference in igneous oceanic crustal thickness is found north (average ∼5.2 km) and south (average ∼5.9 km), of the MFZ that also correspond to changes in Moho characteristics. The MFZ itself is associated with the largest crustal thickness variations (∼3.7–7.6 km). Magma ascent through these different crustal regions may account for some variations of magmatic flux to the surface along the Hawaiian Ridge.

Hawaii↗

USGS2025

Kilometers-scale subsidence of the inner Puerto Rico Trench wall since the Pleistocene

We have observed a thin veneer of Pleistocene-age carbonate rocks on the southern side of Mona Block, a seamount northwest of Puerto Rico, above a water depth of 1800 m. Mona Block is located within the inner wall of the Puerto Rico Trench. A similar veneer of carbonate rocks (corals and pavement) was encountered during a deep submergence vehicle (DSV) Alvin dive in 1976 on the north side of Mona Block at depths of 3650–2900 m. We also observed a remnant Oligocene–Pliocene carbonate platform at a depth of 2080 m on the east side of the block. Restoring the carbonate outcrops to near sea level implies that Mona Block, with a summit currently at 670 m bsl, was an island with an elevation of ≥1000 m until the mid-Pleistocene and has subsided and tilted to the north thereafter. We suggest that the observed subsidence was caused by the foundering or retreating North American slab, which gave rise to the anomalously deep Puerto Rico Trench. Our estimated trench wall subsidence varies nonlinearly from ~800 m in the distal part to ≥3700 m near the trench. Carbonate platform subsidence north of Puerto Rico is older (mid-Pliocene), suggesting that subsidence and, by inference, trench collapse had propagated westward with time. Mona Block is thus a yardstick for studying mantle processes associated with slab retreat and foundering.

Geology↗

USGS2025

Evidence for petit-spot volcanism in the Puerto Rico Trench

Petit-spot volcanism occurs in intraplate settings along the outer rise of subducting plates. Here we present evidence for petit-spot type of volcanism from multibeam bathymetry and backscatter data from the Puerto Rico Trench (PRT). It is the first report of such volcanism in the Atlantic basin. Up to 34 possible petit-spot volcanic cones are mapped in the eastern section of the PRT, with an average diameter of 950 m and a mean height of 92 m. The seamounts are in a region of the PRT that has the longest and highest bending faults, up to 50.6 km long and 1.49 km high. Some of the volcanic cones are associated with fault scarps. The larger bending fault scarps in the region of the volcanic cones suggest the stress state in the eastern PRT is favorable for allowing asthenospheric melt to escape to the surface.

Puerto Rico Trench↗

USGS2025

Earthquake recurrence estimates for northern Caribbean faults from combinatorial optimization

We use combinatorial optimization to find the optimal spatial distribution of random samples of earthquakes (≥6.5) that minimize the misfit in target slip rates for all faults in the northeast Caribbean, and we derive magnitude-frequency relationships with uncertainties for these faults. Slip rates for many faults are derived from GPS block models, not direct measurements, because of their underwater locations. Predicted recurrence rates for eastern Hispaniola and Puerto Rico Trench faults are 220-450 yr for M7 and 3-5 kyr for M8, with maximum feasible magnitude of M8.2. The most frequent earthquakes with magnitudes ≥7.0 are predicted on the large upper plate strike-slip faults, Enriquillo (EF) and Septentrional Fault, commensurate with the historical record. Calais et al. (2023) suggested that shortening in western Hispaniola is accommodated on the offshore Jérémie and onshore Malpasse faults north and south of EF, instead of on terrestrial faults in western Hispaniola and EF. Because of our system-modeling approach, such a configuration predicts less frequent earthquakes on EF and on western Hispaniola and Muertos convergent zones. Recurrence times of a few 100s yr for M6.7 earthquakes is predicted on the submerged faults in Mona Passage, and infrequent M>7 earthquakes are predicted on the Virgin Islands faults.

Lesser Antilles↗

USGS2024

Field geology under the sea with a remotely operated vehicle: Mona Rift, Puerto Rico

We implemented concepts of field geology at great ocean depths by constructing virtual outcrops from a string of overlapping video frames collected by remotely operated vehicles (ROVs). This lower-cost alternative to drilling boreholes allows stratigraphic extension into the offshore and regional interpretation of marine seismic profiles. The imagery was collected along a dive transect on the western wall of Mona Rift, a deep and narrow rift northwest of Puerto Rico, between water depths of 1560 m and 3927 m. The northern coast of Puerto Rico and its large offshore area are underlain by a mid-Eocene and younger forearc basin topped by a thick carbonate platform. There are no drill holes offshore, and tying seismic lines across the shoreline there is problematic. We describe our virtual outcrop and constrain its age and stratigraphy using seven rock samples collected by ROV and compare the outcrop's stratigraphy to deep boreholes and outcrops on land. Our formation descriptions and ages agree, for the most part, with those on land, but we identified a 100-m-thick section that is represented on land by an unconformity. Our stratigraphic interpretation indicates lateral variations in formation thicknesses and establishes a cross-section for additional sampling of the Eocene–Pliocene geology. It also suggests that Mona Rift has formed since the mid-Pliocene. The presence or absence of ferromanganese (Fe-Mn) crust on rocks along the transect may be correlated with the smoothness of the rock surface.

Mona Rift, Puerto Rico↗

USGS2024

Controls on the stratigraphic architecture of the US Atlantic margin: Processes forming the accommodation space

Accommodation space governs the spatial and temporal distributions of sediments in continental margins. Mapping the sedimentation patterns, therefore, offers insights into the solid-Earth processes that shape accommodation space. We assembled an unprecedented amount of seismic and borehole data along the Eastern North American Margin and used it to divide the margin's sedimentary package into eight chronostratigraphic intervals, identifying temporal shifts in depocenters under the continental shelf, slope, and rise. The Jurassic depocenters follow the syn-rift structure and its thermal subsidence loci. The Long Island Platform is the only margin segment where the early post-rift sediment thickness matches subsidence predictions from uniform-stretching models, whereas in Georges Bank Basin (GBB) and Baltimore Canyon Trough (BCT), sediment thickness is 1.5–3 times higher than predicted, pointing to other factors at play. A margin-wide Jurassic transient shoulder uplift is inferred from the occurrence of stratigraphic onlaps above thinned crust. Unlike the Jurassic, the Cretaceous and Cenozoic depocenters disregard the inherited subsidence pattern. The accommodation space over the shelf and coastal plain during the Cretaceous was affected by regional isostatic compensation of the sedimentary loads accumulated on the shelf and rise. Accommodation space development in the GBB was interrupted during the Cretaceous after the margin crossed the Great Meteor Hotspot track, resulting in a widespread permanent uplift, erosion, and sediment redistribution. The distribution of anomalous Neogene subsidence in the BCT challenges previous suggestions of mantle dynamic control on the accommodation space and favors flexural downwarping of the shelf by sediment accumulation on the rise.

Journal of Geophysical Research: Solid Earth↗

USGS2024

Geologic input databases for the 2025 Puerto Rico – U.S. Virgin Islands National Seismic Hazard Model update: Crustal faults component

The last National Seismic Hazard Model (NSHM) for Puerto Rico and the U.S. Virgin Islands (PRVI) was published in 2003. In advance of the 2025 PRVI NSHM update, we created three geologic input databases to summarize new onshore and offshore fault source information in the northern Caribbean region between 62°–70° W and 16°–21° N. These databases, of fault sections, fault‐zone polygons, and geologic estimates of fault activity (fault‐slip rate and earthquake recurrence intervals) at specific sites, document updates to fault parameters used in prior seismic hazard models in PRVI. Fault sources were reviewed from published studies since 2003, which document substantial changes to the understanding of fault location, geometry, or activity. New fault section sources were added for features that meet the criteria of (1) length ≥7 km, (2) unequivocal evidence of recurrent tectonic Quaternary activity, and (3) documentation that is publicly available in a peer‐reviewed source. In addition, we revised several broad areal sources, such as the Mona and Anegada extensional zones. The 2003 model included three fault sections and two fault‐zone polygons (areal sources). These databases include 35 fault sections, 6 fault‐zone polygons, and 51 earthquake geology sites. To characterize fault activity rates, slip‐rate bins were assigned based on landscape expression and paleoseismic trench observations for faults without published slip‐rate sites. Additional fault sources were evaluated but not included in these databases due to a lack of published information about fault location, geometry, or recurrent Quaternary activity. The PRVI NSHM 2025 geologic input databases describe crustal faulting; the geometries and coupling of Puerto Rico subduction zone and Muertos Trough models are considered in a separate database. Updates to the fault sections, fault‐zone polygons, and earthquake geology databases can help inform the location and recurrence rate of damaging earthquakes in the PRVI NSHM implementation.

Puerto Rico, U.S. Virgin Islands↗

USGS2024

Photogrammetry of the deep seafloor from archived unmanned submersible exploration dives

Large amounts of video images have been collected for decades by scientific and governmental organizations in deep (>1000 m) water using manned and unmanned submersibles and towed cameras. The collected images were analyzed individually or were mosaiced in small areas with great effort. Here, we provide a workflow for utilizing modern photogrammetry to construct virtual geological outcrops hundreds or thousands of meters in length from these archived video images. The photogrammetry further allows quantitative measurements of these outcrops, which were previously unavailable. Although photogrammetry had been carried out in recent years in the deep sea, it had been limited to small areas with pre-defined overlapping dive paths. Here, we propose a workflow for constructing virtual outcrops from archived exploration dives, which addresses the complicating factors posed by single non-linear and variable-speed vehicle paths. These factors include poor navigation, variable lighting, differential color attenuation due to variable distance from the seafloor, and variable camera orientation with respect to the vehicle. In particular, the lack of accurate navigation necessitates reliance on image quality and the establishment of pseudo-ground-control points to build the photogrammetry model. Our workflow offers an inexpensive method for analyzing deep-sea geological environments from existing video images, particularly when coupled with rock samples.

Journal of Marine Science and Engineering↗

USGS2024

Modeled flooding by tsunamis and a storm versus observed extent of coral erratics on Anegada, British Virgin Islands— Further evidence for a great Caribbean earthquake six centuries ago

Models of near-field tsunamis and an extreme hurricane provide further evidence for a great precolonial earthquake along the Puerto Rico Trench. The models are benchmarked to brain-coral boulders and cobbles on Anegada, 125 km south of the trench. The models are screened by their success in flooding the mapped sites of these erratics, which were emplaced some six centuries ago. Among 25 tsunami scenarios, 19 have megathrust sources and the rest posit normal faulting on the outer rise. The modeled storm, the most extreme of 15 hurricanes of category 5, produces tsunami-like bores from surf beat. In the tsunami scenarios, simulated flow depth is 1 m or more at all the clast sites, and 2 m or more at nearly all, given either a megathrust rupture 255 km long with 7.5 m of dip slip and M8.45, or an outer-rise rupture 130 km long with 11.4 m of dip slip and M8.17. By contrast, many coral clasts lie beyond the reach of simulated flooding from the extreme hurricane. The tsunami screening may underestimate earthquake size by neglecting trees and shrubs that likely impeded both the simulated flows and the observed clasts; and it may overestimate earthquake size by leaving coastal sand barriers intact. The screening results broadly agree with those from previously published tsunami simulations. In either successful scenario, the average recurrence interval spans thousands of years, and flooding on the nearest Caribbean shores begins within a half-hour.

Anegada↗

USGS2023

Crustal structure across the central Dead Sea Transform and surrounding areas: Insights into tectonic processes in continental transforms

New geophysical profiles across the central Dead Sea Transform (DST) near the Sea of Galilee, Israel, and surrounding highlands, augmented by static stress modeling, allow us to study continental transform plate deformation. The DST separates a ∼10 km thick sedimentary column above a thinned (16–23 km) crust to the west from a ∼7 km column above a ∼30-km thick crust to the east. Crustal thinning starts under the DST, as observed also farther south, indicating that the DST is indeed located along the boundary between the Arabian plate and its continental margin. Moho step here is gradual. The DST's eastern shoulder dips westward toward the DST unlike the upward flexed shoulder observed farther south, perhaps delineating the northern limit of a thinner and hotter lithosphere. The shape of the Sea of Galilee is modeled as an asymmetric pull-apart basin formed by a left-lateral stepover of 2.6 km between slightly divergent and underlapping strike-slip fault strands dipping 70° to the west. Reflection data indicate that these strands are not connected. Several fault traces within the Sea of Galilee have previously been suggested to carry part of the relative plate motion. However, given slip along the main DST faults, Coulomb stress will increase only on fault portions in the northern part of the lake, in accord with the geographical distribution of seismicity, suggesting that these faults are likely secondary. Mismatch between the DST strand locations in the geophysical profiles and the subsidence model, may reflect temporal changes in fault geometry.

Tectonics↗

USGS2023

Seafloor observations eliminate a landslide as the source of the 1918 Puerto Rico Tsunami

The 11 October 1918 devastating tsunami in northwest Puerto Rico had been used as an example for earthquake‐induced landslide tsunami hazard. Three pieces of evidence pointed to a landslide as the origin of the tsunami: the discovery of a large submarine landslide scar from bathymetry data collected by shipboard high‐resolution multibeam sonar, reported breaks of submarine cable within the scar, and the fit of tsunami models to flooding observations. Newly processed seafloor imagery collected by remotely operated vehicle (ROV) show, however, pervasive Fe–Mn crust (patina) on the landslide walls and floor, indicating that the landslide scar is at least several hundred years old. C 14 "> 14 C dates of sediment covering the landslide floor verify this interpretation. Although we have not searched the region systematically for an alternative tsunami source, we propose a possible source—a two‐segment normal‐fault rupture along the eastern wall of Mona rift. The proposed fault location matches the published normal faults with steep bathymetry and is close to the International Seismological Center–Global Earthquake Model catalog locations of the 1918 mainshock and aftershocks. The ROV observations further show fresh vertical slickensides and rock exposure along the proposed fault trace. Hydrodynamic models from an M w "> M w 7.2 earthquake rupture along the eastern wall of the rift faithfully reproduce the reported tsunami amplitudes, polarities, and arrival times. Our analysis emphasizes the value of close‐up observations and physical samples to augment remote sensing data in natural hazard studies.

Puerto Rico↗

USGS2023

Late Pleistocene-Holocene age and stratigraphy of the Currituck Slide Complex, U.S. mid-Atlantic continental slope: Implications for landslide triggering

Considerable effort has been made to link submarine slope failures to changes in local and global-scale environmental conditions, in order to assess landslide hazard probability. Here we provide the first radiocarbon dates of hemipelagic sediment overlying mass transport deposits and inferred failure surfaces of the Currituck Slide Complex (CSC), a prominent landslide scar on the U.S. mid-Atlantic continental slope. The dates, taken from both the upper and lower scars of the complex, constrain the age of the last major failure event to 13,835 and 16,020 years BP. Time correlation of the hemipelagic sediments across the landslide scar and proximal deposit suggests a single failure of both the upper and lower parts of this 160 km 3 volume. A higher rate of sediment supply from the periglacial Appalachian Mountains and from glacial melt-water pulses, with an exposed continental shelf at that time, may have enlarged a shelf-edge delta at the site of the CSC, which may have facilitated or triggered failure. A smaller landslide at the southern edge of the complex with a less well-defined geomorphologic footprint is dated at 5500 BP and possibly represents the reshaping of the seafloor around the CSC triggered by low-frequency earthquakes on the nearby continental margin.

Marine Geology↗

USGS2022

Quantifying permanent uplift due to lithosphere-hotspot interaction

Vertical motions that accompany the passage of the lithosphere over a mantle hotspot can shed light on the nature of the hotspot and its effect on the lithosphere. However, quantifying the temporal vertical and spatial extent, is challenging due to the paucity of evidence in the geological record. Here, we utilize dense seismic and well data covering the intersection of the Great Meteor Hotspot (GMH) track with the U.S. Atlantic continental margin to constrain the surface expression of the hotspot passage under the lithosphere. The continuous sedimentary record of the eastern North American margin during its passage over the hotspot allows determination of the timing, magnitude, width and rate of denudation. We find that a ∼300 km wide region was denuded by up to 850 m between ∼97 and 86 Ma, ∼10 m.y. after the passage of the GMH. Stratigraphic relationships suggest a decaying rock uplift rate with time and no subsequent sagging. The broad, long-lasting, and delayed uplift was modeled as a surface manifestation of either sub-lithospheric mantle depletion, permanently eroded base of the continental lithosphere, or intrusions of depleted magma. We consider sub-lithospheric depletion to be the most likely cause, based on seismic imaging results.

Tectonics↗

USGS2022

Systematic mapping of the ocean-continent transform plate boundary of the Queen Charlotte fault system, southeastern Alaska and western British Columbia—A preliminary bathymetric terrain model

In 2015, U.S. Geological Survey scientists in collaboration with scientists from other institutions began a study of the Queen Charlotte fault—the first systematic study of the fault in more than three decades. The primary goal of the study was to gain a better understanding of the earthquake, tsunami, and underwater-landslide hazards throughout southeastern Alaska, as well as gather data to develop geologic models that can be applied to similar plate boundaries around the globe, such as the San Andreas fault system in southern California, the Alpine fault in New Zealand, and the North Anatolian fault in Turkey. A bathymetric terrain model was compiled from six different multibeam surveys of the previously unmapped Queen Charlotte fault offshore of southeastern Alaska and Haida Gwaii archipelago.

Alaska, British Columbia↗

USGS2022

Mature diffuse tectonic block boundary revealed by the 2020 southwestern Puerto Rico seismic sequence

Distributed faulting typically tends to coalesce into one or a few faults with repeated deformation. The progression of clustered medium-sized (≥Mw4.5) earthquakes during the 2020 seismic sequence in southwestern Puerto Rico (SWPR), modeling shoreline subsidence from InSAR, and sub-seafloor mapping by high-resolution seismic reflection profiles, suggest that the 2020 SWPR seismic sequence was distributed across several short intersecting strike-slip and normal faults beneath the insular shelf and upper slope of Guayanilla submarine canyon. Multibeam bathymetry map of the seafloor shows significant erosion and retreat of the shelf edge in the area of seismic activity as well as slope-parallel lineaments and submarine canyon meanders that typically develop over geological time. The T -axis of the moderate earthquakes further matches the extension direction previously measured on post early Pliocene (∼>3 Ma) faults. We conclude that although similar deformation has likely taken place in this area during recent geologic time, it does not appear to have coalesced during this time. The deformation may represent the southernmost part of a diffuse boundary, the Western Puerto Rico Deformation Boundary, which accommodates differential movement between the Puerto Rico and Hispaniola arc blocks. This differential movement is possibly driven by the differential seismic coupling along the Puerto Rico—Hispaniola subduction zone. We propose that the compositional heterogeneity across the island arc retards the process of focusing the deformation into a single fault. Given the evidence presented here, we should not expect a single large event in this area but similar diffuse sequences in the future.

Puerto Rico↗