Chemical compositions of sediment pore waters are presented from 13 piston and gravity cores that were collected on similar to 24 Ma crust of the Cocos Plate seaward of the Middle America Trench and near the onset of crustal faulting from subduction. Cores were collected mainly within a 1.75 km2 area overlying a buried basement topographic high that supports an elevated heat flux, consistent with seawater transport within the upper volcanic crust. Systematic variations in pore water chemical profiles indicate upward seepage speeds (up to 1.7 cm yr-1 providing a net flux of 0.1 L s-1), constrain the chemical composition of the formation water within the underlying upper basaltic basement, and elucidate diagenetic reactions in the sediment. Relative to seawater, formation water has an elevated temperature (70-80 degrees C) and concentrations or values of Ca, chlorinity, Sr, Ba, Li, Fe, Mn, Si, Cs, D/H, and Mo, and lower concentrations or values of Mg, Na, sulfate, alkalinity, TCO2, K, B, F, phosphate, 87Sr/86Sr, delta 13C, delta 18O, U, and Rb. Although this site is located only 30 km from the trench axis, there is no chemical evidence for subduction-related hydrologic influences. Instead, the data are explained by a combination of seawater-basalt reactions within the upper basement and diffusive exchange with overlying sediment, as part of a shallow, ridge-flank hydrothermal system. It is unclear why this site has an elevated heat flux relative to neighboring crust, but this may result from variations in crustal properties or modification related to flexural faulting outboard of the trench. An area of elevated volcanic rock buried below sediment west of the Middle American Trench is unusually warm, with temperatures of 70-80 degrees C compared to more regional-distributed crust that is vigorously cooled (10-20 degrees C) from circulating seawater (formation water) within the upper volcanic crust. Because of the proximity to subduction-related crustal faults, sediment cores were collected from this area of warm crust and interstitial sediment pore waters were extracted and analyzed to assess if and how this area with a high heat flux is linked to subduction processes, potentially providing a source of water that serpentinizes the underlying mantle. Chemical data from sediment pore waters confirm upward seepage, the composition of crustal formation water in underlying volcanic rocks, and the extent of reactions as the formation water ascended the sediment column. Crustal formation water appears to be unaffected by nearby plate subduction. Instead, this warm and chemically altered formation water defines a ridge-flank hydrothermal system, which is distinct from hydrothermal systems to the west and south. How this hydrothermal system became isolated is unknown, but it could be a result of aging and evolving crustal properties or of plate faulting related to subduction processes. Pore water chemical profiles reveal sediment diagenesis, seepage speeds up to 1.7 cm yr-1, and discharge from 1.75 km2 of 0.1 L s-1Crustal formation water is warm (similar to 75 degrees C) and chemically altered, stemming from water-basalt reactions and diffusive exchange with pore waterThis ridge-flank hydrothermal system is hydrologically isolated from the ventilated crust to the west and the trench to the east

Abstract High‐resolution bathymetry and three‐dimensional seismic data along the Cocos Ridge reveal a 245 km2 field of ∼1–4 km in diameter seafloor depressions. The seafloor depressions are part of a two‐tiered honeycomb pattern. The lower‐tier depressions have steep faults that truncate strata with chaotic internal reflections consistent with sediment collapse into the depression. These extend into a lens shaped interval just above igneous basement. Overlying these depressions is a second broader set with rough seafloor morphology with gently dipping boundaries defined by pinch‐out stratigraphic patterns. Drilling results indicate that the lens‐shaped zones that host the deeper depressions represent anomalous regions of high porosity, low velocity, and low density within calcareous rich sediment. Analysis of nannofossils from IODP Site U1414 suggests the collapse structures formed during the late Miocene, whereas the younger shallower depressions likely formed between the early Pliocene and the Pliocene‐Pleistocene boundary. Geochemical and petrological analysis at Site U1414 suggests that hydrothermal circulation during the late Miocene led to carbonate dissolution and collapse. Following collapse, focused fluid‐flow and bottom current scouring resulted in formation of the overlying set of depressions and a honeycomb seafloor morphology. Similar sets of depressions along the Carnegie Ridge to the south support the hypothesis that two‐tiered depressions formed in response to processes that occurred broadly across the Panama Basin between the late Miocene and the Pliocene‐Pleistocene transition. Geochemical results at Site U1414, combined with geophysical data, suggest this two‐tiered system of depressions currently guides ongoing fluid outflow.

Logging data are measurements of physical properties of the formation surrounding a borehole, acquired in situ after completion of coring (wireline logging) or during drilling (Logging-While-Drilling, LWD). The range of data (resistivity, gamma radiation, velocity, density, borehole images,…) in any hole depends on the scientific objectives and operational constraints.

Abstract The Costa Rican convergent margin has been considered a type erosive margin, with erosional models suggesting average losses up to −153 km3/km/m.y. However, three‐dimensional (3D) seismic reflection and Integrated Ocean Drilling Program data collected offshore the Osa Peninsula images accretionary structures and vertical motions that conflict with the forearc basal erosion model. Here we integrate such data to do an in‐situ accounting of material transfer at the plate boundary across the outermost 10 km of the forearc, characterized by active and inactive megathrusts. Our in‐situ budget finds an approximate balance between sediment recycling via accretion and underplating, 0.7–2.3 km3/km/m.y., and basal erosion, 0.7 km3/km/m.y., while subducting sediment volumes, 7.8 km3/km/m.y., greatly outpace either material transfer volumes. These budget results differ significantly from published estimates based on simple proxies of trench axis deflection and slope subsidence. These budget results are the summation of thin incoming hemipelagic sediments that variably accrete along the deformation front, underplating of hemipelagic sediments on the upthrown‐side and basal erosion on the downthrown‐side of active plate bending faulting landward of the trench axis, and sediment subduction primarily composed of pelagic sediments.

The bottom simulating reflection (BSR) is widely observed along continental margins and is believed to mark the base of gas hydrate stability zone (BGHSZ). In some regions, double or multiple overlapping BSRs are observed, yet their formation mechanisms and geologic implications are not well understood. Here we present 3D seismic images from the 2018 NZ3D experiment that covers a 14 x 60 km(2) survey area on New Zealand's northern Hikurangi subduction margin. We observe double BSRs in five locations. Beneath the Tuaheni Basin in the mid-slope, a secondary BSR (BSR2) lies similar to 100-360 m deeper than the primary BSR (BSR1) and its 3D geometry mimics the unconformity at the base of the basin. At three thrust ridges located 18-38 km from the deformation front, BSR2 lies similar to 55-130 m below and is subparallel to BSR1. At another thrust ridge similar to 14 km from the deformation front, BSR2 forms above the BSR1, and the two BSRs converge towards the peak of the ridge. Through 3D modeling of BGHSZ and analysis of the geometry and reflection characteristics of the double BSRs, we identify three potential mechanisms for their formation (1) rapid sedimentation, (2) tectonic uplift, (3) overpressure/ heat advection caused by fluid migration. Our study demonstrates that formation of double BSRs is closely linked to subduction processes along the northern Hikurangi margin, and double BSRs may be used as indicators for areas with recent sedimentation, tectonic and/or fluid activities. (C) 2021 Elsevier B.V. All rights reserved.

Megathrust roughness and structural complexity are thought to be controls on earthquake slip at subduction zones because they result in heterogeneity in shear strength and resolved stress. However, because active megathrust faults are difficult to observe, the causes and scales of complexity are largely unknown. Here we measured the in situ properties of the megathrust of the Middle America subduction zone in a three-dimensional seismic reflection volume to determine how fault properties vary. We quantify spatial variability in the megathrust roughness, overburden and rock physical properties. Heterogeneity in the megathrust roughness exists at length scales of a few kilometres because the megathrust is dissected by active lower-plate normal faults, which offset the megathrust and renewed fault roughness. Spatial variations in the rock physical properties at the plate interface are characterized by correlation length scales of hundreds of metres. Frontal prism taper, historical seismicity and the variation in earthquake stress drop values local to the megathrust are all affected by the heterogeneity at these length scales. Both geometric and rheological complexities may therefore control the mechanical behaviour of the subduction plate interface, which includes earthquake rupture characteristics. Geometric and rheological complexities may control the mechanical behaviour of megathrusts, according to an analysis of the heterogeneity in roughness and rock properties of the Middle America megathrust from 3D seismic reflection data.

In the version of this Article originally published, in Fig. 1a the Nicoya Peninsula, Osa Peninsula and Burica Peninsula were incorrectly labelled as the Nicoya Plate, Osa Plate and Burica Plate, respectively. This has been corrected in the online versions.

Understanding the links between subducting slabs and upper-plate deformation is a longstanding goal in the field of tectonics. New 3D seismic sequence stratigraphy, mapped within the Costa Rica Seismogenesis Project (CRISP) seismic-reflection volume offshore southern Costa Rica, spatiotemporally constrains several Pleistocene outer forearc processes and provides clearer connections to subducting plate dynamics. Three significant shelf and/or slope erosional events at ca. 2.5-2.3 Ma, 1.95-1.78 Ma, and 1.78-1.19 Ma, each with notable differences in spatial extent, volume removed, and subsequent margin response, caused abrupt shifts in sedimentation patterns and rates. These shifts, coupled with observed deformation, suggest three primary mechanisms for Pleistocene shelf and slope vertical motions: (1) regional subaerial erosion and rapid subsidence linked to the southeastward Panama Fracture Zone triple-junction migration, with associated abrupt bathymetric variations and plate kinematic changes; (2) transient, kilometer-scale uplift and subsidence due to inferred subducting plate topography; and (3) progressive outer wedge shortening accommodated by landward-and seaward-dipping thrust faults and fold development due to the impinging Cocos Ridge. Furthermore, we find that the present-day wedge geometry (to within similar to 3 km along strike) has been maintained through the Pleistocene, in contrast to modeled landward margin retreat. We also observe that deformation, i.e., extension and shortening, is decoupled from net margin subsidence. Our findings do not require basal erosion, and they suggest that the vertical motions of the Costa Rican outer forearc are not the result of a particular continuous process, but rather are a summation of plate to plate changes (e.g., passage of a fracture zone triple junction) and episodic events (e.g., subducting plate topography).

Exhumed faults are rough, often exhibiting topographic corrugations oriented in the direction of slip; such features are fundamental to mechanical processes that drive earthquakes and fault evolution. However, our understanding of corrugation genesis remains limited due to a lack of in situ observations at depth, especially at subducting plate boundaries. Here we present three-dimensional seismic reflection data of the Costa Rica subduction zone that image a shallow megathrust fault characterized by corrugated, and chaotic and weakly corrugated topographies. The corrugated surfaces extend from near the trench to several kilometres down-dip, exhibit high reflection amplitudes (consistent with high fluid content/pressure) and trend 11–18° oblique to subduction, suggesting 15 to 25 mm yr−1 of trench-parallel slip partitioning across the plate boundary. The corrugations form along portions of the megathrust with greater cumulative slip and may act as fluid conduits. In contrast, weakly corrugated areas occur adjacent to active plate bending faults where the megathrust has migrated up-section, forming a nascent fault surface. The variations in megathrust roughness imaged here suggest that abandonment and then reestablishment of the megathrust up-section transiently increases fault roughness. Analogous corrugations may exist along significant portions of subduction megathrusts globally. Mature parts of the shallow megathrust beneath Costa Rica are characterized by striking corrugations that may channel fluids, according to seismic images. Nascent sections of the subduction zone plate boundary appear only weakly corrugated.

In 2011 we acquired a 3D seismic reflection volume across the Costa Rica margin NW of the Osa Peninsula to investigate the complex structure and the development of the seismogenic zone within the Costa Rican subduction zone in the vicinity of recent International Ocean Drilling Program (IODP) drilling. In contrast to previous interpretations, these newly acquired seismic images show that the margin wedge is composed of a layered fabric that is consistent with clastic sediments, similar to materials recovered from IODP drilling, that have been thrust and thickened into thrust-bounded folded sequences. These structures are consistent with a balanced sequence that has been frontally accreted in the context of an accretionary model. We interpret these sequences as sediment originally deposited on the subducting crust in a trench basin created by the southward migration of the Cocos–Nazca–Caribbean triple junction, and accreted during recent margin subduction that also accelerated with passage of the triple junction. The margin is composed of relatively rapidly accreted sediment that was added to the margin during a phase of accretion within the last ∼5 Ma that was probably preceded throughout the Neogene by periods of non-accretion or erosion.

Here we propose a new framework for forearc evolution that focuses on the potential feedbacks between subduction tectonics, sedimentation, and geomorphology that take place during an extreme event of subduction erosion. These feedbacks can lead to the creation of a “depositionary forearc,” a forearc structure that extends the traditional division of forearcs into accretionary or erosive subduction margins by demonstrating a mode of rapid basin accretion during an erosive event at a subduction margin. A depositionary mode of forearc evolution occurs when terrigenous sediments are deposited directly on the forearc while it is being removed from below by subduction erosion. In the most extreme case, an entire forearc can be removed by a single subduction erosion event followed by depositionary replacement without involving transfer of sediments from the incoming plate. We need to further recognize that subduction forearcs are often shaped by interactions between slow, long‐term processes, and sudden extreme events reflecting the sudden influences of large‐scale morphological variations in the incoming plate. Both types of processes contribute to the large‐scale architecture of the forearc, with extreme events associated with a replacive depositionary mode that rapidly creates sections of a typical forearc margin. The persistent upward diversion of the megathrust is likely to affect its geometry, frictional nature, and hydrogeology. Therefore, the stresses along the fault and individual earthquake rupture characteristics are also expected to be more variable in these erosive systems than in systems with long‐lived megathrust surfaces.

In 2011 we acquired an 11x55 km, 3-D seismic reflection volume across the Costa Rica margin, NW of the Osa Peninsula, to accurately image the subduction thrust in 3-D, to examine fault zone properties, and to infer the hydrogeology that controls fluid accumulation along the thrust. Following processing to remove water column multiples, noise, and acquisition artifacts, we constructed a 3-D seismic velocity model for Kirchhoff prestack depth migration imaging. Images of the plate boundary thrust show high-reflection amplitudes underneath the middle to lower slope that we attribute to fluid-rich, poorly drained portions of the subduction thrust. At similar to 5km subseafloor, beneath the upper slope, the plate interface abruptly becomes weakly reflective, which we interpret as a transition to a well-drained subduction thrust. Mineral dehydration during diagenesis may also diminish at 5km subseafloor to reduce fluid production and contribute to the downdip change from high to low amplitude. There is also a layered fabric and systems of both thrust and normal faults within the overriding plate that form a plumbing system. Faults commonly have fault plane reflections and are presumably fluid charged. The faults and layered fabric form three compartmentalized hydrogeologic zones: (1) a shallow NE dipping zone beneath the slope, (2) a steeply SW dipping zone beneath the shelf slope break, and (3) a NE dipping zone beneath the shelf. The more direct pathway in the middle zone drains the subduction thrust more efficiently and contributes to reduced fluid pressure, elevates effective stress, and creates greater potential for unstable coseismic slip.

The March 13th 1888 collapse of Ritter Island in Papua New Guinea is the largest known sector collapse of an island volcano in historical times. One single event removed most of the island and its western submarine flank, and produced a landslide deposit that extends at least 70 km from the headwall of the collapse scar. We have mapped and described the deposits of the debris avalanche left by the collapse using full-coverage multibeam bathymetry, side-scan sonar backscatter intensity mapping, chirp seismic-reflection profiles, TowCam photographs of the seafloor and samples from a single dredge. Applying concepts originally developed on the 1980 Mount St. Helens collapse landslide deposits, we find that the Ritter landslide deposits show three distinct morphological facies: large block debris avalanche, matrix-rich debris avalanche and distal debris flow facies. Restoring the island's land and submarine topography we obtained a volume of 4.2 km3 for the initial collapse, about 75% of which is now forming the large block facies at distances less than 12 km from the collapse scar. The matrix-rich facies volume is unknown, but large scale erosion of the marine sediment substrate yielded a minimum total volume of 6.4 km3 in the distal debris flow and/or turbidite deposits, highlighting the efficiency of substrate erosion during the later history of the landslide movement. Although studying submarine landslide deposits we can never have the same confidence that subaerial observations provide, our analysis shows that well-exposed submarine landslide deposits can be interpreted in a similar way to subaerial volcano collapse deposits, and that they can in turn be used to interpret older, incompletely exposed submarine landslide deposits. Studying the deposits from a facies perspective provides the basis for reconstructing the kinematics of a collapse event landslide; understanding the mechanisms involved in its movement and deposition; and so providing key inputs to tsunami models.