Trench Axis Research Articles

Seismic imaging of Northwest Pacific and East Asia: New insight into volcanism, seismogenesis and geodynamics

Recent studies on seismic imaging of the Northwest Pacific and East Asian region are reviewed. High-resolution tomographic images reveal significant lateral heterogeneities in the crust and upper mantle, which shed new light on interplate and intraplate volcanism, earthquake mechanism, and mantle dynamics. Significant recent advances in seismic imaging are tomographic inversions for 3-D distribution of seismic anisotropy and attenuation in the crust and upper mantle, which provide important new information on lithospheric deformation, mantle convection, melts and fluids associated with plate subductions and mantle dynamics. The results show that intraplate volcanism in NE Asia is caused by hot and wet upwelling flows in the big mantle wedge above the stagnant Pacific slab in the mantle transition zone. The subducting Pacific and Philippine Sea slabs exhibit mainly trench-parallel fast-velocity directions (FVDs) of azimuthal anisotropy, which may reflect frozen-in lattice-preferred orientation of aligned anisotropic minerals and/or shape-preferred orientation, such as transform faults formed at the mid-ocean ridge and normal faults produced at the outer-rise area near the trench axis. Trench-normal FVDs are generally revealed in the mantle wedge under the volcanic front and back-arc region, which may reflect corner flows in the mantle wedge due to the plate subduction and dehydration. Trench-normal FVDs also appear in the subslab mantle, which may reflect asthenospheric shear deformation associated with the overlying slab subduction. The nucleation of large earthquakes, including both crustal and megathrust events, is controlled by structural heterogeneities in and around seismogenic faults. The cause of deep-focus earthquakes is a subject of ongoing debate. A recent example is the 2015 Bonin deep earthquake (M 7.9, ~670 km depth), which was probably caused by joint effects of several factors, including the Pacific slab's fast deep subduction, slab tearing and thermal variation, stress changes and phase transformations in the slab, and complex interactions between the slab and ambient mantle. The key to future advances in seismic imaging is instrumentation. Gradual deployment of seismometers on seafloors and in those less instrumented land areas will be the most important task for seismologists from now. For this purpose, international cooperation of seismologists in different countries will be necessary

Tsunami damage estimation in Esmeraldas, Ecuador using fragility functions

<abstract> <p>The current study investigated the probable impact from a tsunami to a populated area located along the northwest ecuadorian coast, specifically in the key oil-industrial city of esmeraldas. a numerical tsunami simulation was performed considering the seismological and tectonic aspects of the area. The damage probability was calculated using fragility functions (ffs). Briefly, 16 cases of source models with slightly different fault parameters were tested, where one was selected as the worst scenario of tsunami inundation. This scenario was a hypothetic earthquake case (mw 8.7) located in front of esmeraldas city, approximately 100 km offshore along the ecuador—colombia trench, with three shallow fault segments (top depth of 10 km), a strike aligned with the trench axis, a middle dip angle of 28°, assuming large slips of 5 to 15 m, and a rake angle of 90°. The results from the numerical simulation were comparable to a similar study previously conducted and with those of historically documented data. The tsunami damage estimation using FFs resulted in estimated damages of 50% and 44% in exposed buildings and population, respectively. Results also showed that the most impacted areas were located next to the coastal shoreline and river. tourism, oil exports, and port activities, in general, would be affected in this scenario; thus, compromising important industries that support the national budget. Results from this study would assist in designing or improving tsunami risk reduction strategies, disaster management, use of coastal zones, and planning better policies.</p> </abstract>

The First 10 Million Years of Rear‐Arc Magmas Following Backarc Basin Formation Behind the Izu Arc

AbstractIODP Site U1437 is located in the Izu rear‐arc region, approximately 330 km west of the Izu‐Bonin trench axis. The oldest four units (Units IV through Unit VII) include volcaniclastic sediment and in situ hyaloclastites. They have ages of about 6–15 Ma, shortly after cessation of Shikoku backarc basin opening. Three magma types are identified by their distinct geochemistry; they are similar types to those found in the modern Izu arc (Rear Arc Seamount Chain [RASC]‐type, Rift‐type, and volcanic front [VF]‐type). RASC‐type has the most enriched Nd and Hf isotope and fluid‐immobile trace element ratios and dominates from 9 to 6 Ma. Rift‐type, dominant from 15 to 9 Ma, is similar to VF‐type in Nd‐Hf isotopes but has the least radiogenic Sr and Pb, and intermediate La/Yb and Nb/Yb, indicating a more fertile mantle source and less hydrous slab component than VF‐type. Less common and randomly distributed VF‐type sediments have the most radiogenic Sr and Pb, and the highest Ba/(Th, LREE [light rare earth element]) ratios, and are interpreted to be distally derived. The genesis of mafic Unit VII samples (~15 Ma) was modeled using the Arc Basalt Simulator. Results are most similar to those for basalts in the modern rift environment indicating the addition of ~1% of a melt‐rich slab component generated at ~125 km, to a Philippine Sea Plate ambient mantle that was more depleted than DMM (depleted MORB mantle). The initial post‐Shikoku basin magmatism in the Izu rear‐arc generated Rift‐type magmas for about 6 million years before the distinctive RASC‐type magmatism began, which then became increasingly enriched.

Calcite twinning fabrics along the Middle America trench, Costa Rica and the Motagua sinistral fault, Honduras and Jamaica: Tectonic implications for the Caribbean plate

The Cretaceous oceanic crust of the Caribbean plate has a long and complex tectonic history as an indentor plate bounded by transcurrent faults (N and S) and subduction margins (W [east slab dip] and E [west slab dip]) with no strain data. Miocene forearc sediments and accreted ophiolites are exposed along the central Pacific coast of Costa Rica and preserve a variety of folding and faulting styles. These clastic sediments, deposited in the forearc as the Middle America trench evolved, contain calcite cements and are cross cut by calcite veins. Analysis of the mechanical twins in these two calcite groups (8 samples, n = 267 twins) at Punta Dominical and the Nicoya Peninsula ophiolite preserve a consistent sub-horizontal shortening strain that is oriented parallel to the current trench axis, as are most of the accretion-related fold axes. Differential stress magnitudes for the calcite twinning strain average −450 bars, and there is little evidence of a twinning strain overprint (low NEVs). Obduction along the Middle America trench involves, in part, a component of trench-parallel horizontal shortening.The Motagua fault originates in the Middle America trench and strikes east as the northern boundary of the Caribbean plate as a sinistral margin. The Roatan Islands are a series of islands composed of south-dipping rocks along the southern side of the Motagua fault. Recent, deformed clastic sediments are drapped over gneisses, schists (40Ar-39Ar biotite age of 14 Ma), conodont-free marbles and ophiolitic fragments of presumed Cretaceous age. Twinning studies (n = 6 samples, n = 425 twins) of these marbles preserve a N–S, horizontal shortening strain, normal to the sinistral plate boundary. Differential stress magnitudes for the calcite twinning strain average −36 MPa. Only in a combined sample (n = 189 grains) is a strain overprint (high NEVs; n = 47) present, which records horizontal shortening strain at 132°, 5°. Jamaica is also along the eastern Motagua fault (Cayman trough) and twinning strain studies in Eocene limestones and younger calcite veins (n = 4 samples; n = 136 twins) preserve a horizontal shortening strain normal to the sinistral plate boundary.Twinning strain results are compared to regional GPS, focal mechanism and SKS anisotropy results around the Caribbean plate, especially along the Middle America trench on the western margin of Costa Rica with some new insights into subduction dynamics.

Characteristic activities of slow earthquakes in Japan.

Slow earthquakes are a recently discovered phenomenon that mainly occur updip and downdip of the seismogenic zones of great earthquakes along the subducting plate interface. The spatiotemporal activity of various slow earthquakes occurring in the Nankai subduction zone is characterized by along-strike heterogeneity and along-dip systematic changes. Various slow earthquakes are horizontally distributed at their own depths and along-strike segments can be observed with respect to this distribution downdip of the locked zone; however, slow and great earthquakes occur in the same depth range near the Nankai Trough and Japan Trench axes. The frequently observed spatiotemporal interactions between different slow earthquakes can be attributed to their sensitivity and the stress transfer of the surrounding areas. This stress transfer is expected to extend to the adjacent sections in the seismogenic zone. Therefore, precise monitoring of slow earthquakes is important for future evaluations of great earthquakes, which requires the long-term maintenance and continuous improvement of the high-quality observation networks.

Slow Earthquakes Illuminating Interplate Coupling Heterogeneities in Subduction Zones

AbstractSlow earthquakes are mainly distributed in the vicinity of seismogenic zones of megathrust earthquakes and relationships between both types of earthquakes are expected. We examined the activity of very low frequency earthquakes (VLFEs), classified as one type of slow earthquake, around Japan because they have the potential to clarify detailed spatiotemporal slip behaviors at the plate boundaries. The distribution of the shallow VLFE activity rate is heterogeneous along trench axes and exhibits an anticorrelation relationship with the spatial distribution of the interplate coupling ratio, whereas deep VLFEs are distributed only in weakly coupled areas, and the spatial variation of the activity rate is small. Furthermore, VLFEs are mainly hosted by low seismic velocity anomalies. Thus, slow earthquakes can be triggered by decreased effective stress due to the high pore fluid pressure within regions with weak interplate coupling, and their activity can be an indicator of interplate slip behavior.

Correlation of frontal prism structures and slope failures near the trench axis with shallow megathrust slip at the Japan Trench

Since the giant 2011 Tohoku earthquake and tsunami, much research has focused on the distribution of coseismic slip at shallow depths during this subduction megathrust event. Here we present seismic images obtained in the immediate vicinity of the trench axis, that show thrust faults and fold-and-thrust type deformation structures near the epicenter of the 2011 Tohoku earthquake where the large coseismic slip has been inferred, and chaotic structure and the absence of thrust faults in northern and southern source areas. Seismic profiles show evidence of slope failures of the trench inner wall in a proposed tsunami source region around 39°–40° N, where the slips estimated from previous studies are in disagreement. Our results show that structural characteristics in the trench axis may be related to the occurrence of shallow megathrust slip and tsunamigenesis in the Japan Trench.

Glycerol dialkyl glycerol tetraethers in surface sediments from three Pacific trenches: Distribution, source and environmental implications

Glycerol dialkyl glycerol tetraethers (GDGTs) have been widely used to elucidate sources of sediment total organic carbon (TOC), past temperature and presence of methanogenesis in diverse environments. However, their applicability to hadal trenches with their unique deposition dynamics remains unknown. Here, we analyzed GDGTs and their stable isotope values and content of TOC in surface sediments from the Kermadec Trench region (KT; 6080–10010 m), New Britain Trench region (NBT; 1553–8931 m), and Atacama Trench region (AT; 2560–8085 m). These regions are at very different distances from terrestrial sources and have varying net primary productivity (NPP) in the waters above them. The GDGT concentration was highly variable (54.5–2416 μg g−1 TOC) within and between trench regions and was not directly related to local NPP or apparent terrestrial inputs. This finding is presumably due to complex deposition dynamics within the trench interior. Isoprenoidal GDGTs (isoGDGTs; 75.4–99.1%) were dominant over branched GDGTs (brGDGTs, 0.91–24.6%) in all samples, leading to low levels of Branched versus Isoprenoidal Tetraether (BIT) index (0.01–0.27). Thus, sediment TOC is mainly derived from marine sources. However, compared to adjacent non-hadal sites, trench axis sites have a higher BIT index, lower acyclic hexa-/pentamethylated brGDGT and lower δ13C, supporting relative enrichment of terrestrial organic carbon at the trench axis. The application of TetraEther indeX of tetraethers consisting of 86 carbon atoms (TEX86) resulted in sea surface temperature (SST) estimates of 18.9–23.7 °C in the KT, 28.6–30.2 °C in the NBT, and 17.9–20.4 °C in the AT. The close agreement between TEX86-SST with observed in situ SST suggests that isoGDGTs are not selectively degraded during the transport towards the hadal realm, and that TEX86 from hadal settings robustly records an integrated regional SST signal.

Seismicity of the Incoming Plate and Forearc Near the Mariana Trench Recorded by Ocean Bottom Seismographs

AbstractEarthquakes near oceanic trenches are important for studying incoming plate bending and updip thrust zone seismogenesis, yet are poorly constrained using seismographs on land. We use an ocean bottom seismograph (OBS) deployment spanning both the incoming Pacific Plate and the forearc to study seismicity near the Mariana Trench. The yearlong deployment in 2012–2013 consisted of 20 broadband OBSs and 5 suspended hydrophones, with an additional 59 short period OBSs and hydrophones recording for 1 month. We locate 1,692 earthquakes using a nonlinear method with a 3D velocity model constructed from active source profiles and surface wave tomography results. Events occurring seaward of the trench occur to depths of ~35 km below the seafloor, and focal mechanisms of the larger events indicate normal faulting corresponding to plate bending. Significant seismicity emerges about 70 km seaward from the trench, and the seismicity rate increases continuously towards the trench, indicating that the largest bending deformation occurs near the trench axis. These plate‐bending earthquakes occur along faults that facilitate the hydration of the subducting plate, and the lateral and depth distribution of earthquakes is consistent with low‐velocity regions imaged in previous studies. The forearc is marked by a heterogeneous distribution of low magnitude (<5 Mw) thrust zone seismicity, possibly due to the rough incoming plate topography and/or serpentinization of the forearc. A sequence of thrust earthquakes occurs at depths ~10 km below seafloor and within 20 km of the trench axis, demonstrating that the megathrust is seismically active nearly to the trench.

Evidence of giant earthquakes and tsunamis of the seventeenth-century type along the southern Kuril subduction zone, eastern Hokkaido, northern Japan: a review

Abstract Large earthquakes of magnitude 7–8 occur at intervals of about 70–100 years in the southern Kuril Trench off eastern Hokkaido, northern Japan. However, especially large giant earthquakes (M w c. 8.8) have recurred at 300–500 year intervals, and the most recent of these was in the early seventeenth century. Surveys on the Pacific coast of eastern Hokkaido have documented field evidence that includes sharp contacts of peat and mud layers, indicative of coseismic coastal uplift, and large tsunami sand layers. Several proposed fault models of the seventeenth-century giant earthquake have been tested with numerical tsunami simulations. The best of these models incorporates a multi-segment interplate earthquake (M w 8.8) with large slip near the trench axis. The seventeenth-century giant earthquake and the 1611 Keicho Sanriku tsunami appear to have had different sources because varve counts from Lake Harutori suggest that the seventeenth-century giant earthquake occurred at c. 1637. Moreover, a consideration of regional tectonics suggests that convergence of the Kuril forearc sliver controls the current topography of eastern Hokkaido. Further progress in understanding the seismicity of this region will benefit from the combined insights of both geologists and seismologists.

Plate geometry model and seismicity in the northern Ryukyu subduction zone, Japan, deduced from amphibious seismic observations

In this study, the geometry of the subducted Philippine Sea Plate interface beneath the northern Ryukyu Trench subduction zone was estimated based on seismic tomography and focal mechanism estimations using passive data obtained from a temporary amphibious seismic network and permanent land stations. Using relocated low-angle thrust-type earthquakes and information obtained from previous studies, we constrained the plate geometry from the trench axis to a depth of 60 km with uncertainties of <2 km between the depths of 15 and 30 km. The estimated plate geometry model appears to be smooth; there is no clear evidence of the existence of large subducted seamounts or ridges in the forearc of the study area. The active earthquake region near the plate boundary was primarily limited to the depth range of 15–30 km. These upper and lower limits may be controlled by the temperature along the plate boundary and the existence of serpentinized mantle in the hanging-wall plate, respectively. Moreover, we determined that the active region of interplate earthquakes is clearly separated from the area dominated by short-term slow slip events and low-frequency tremors, suggesting that suitable frictional conditions along the plate boundary for regular and slow earthquakes are different. There are two areas that were nearly seismically quiet, indicating either full interplate coupling or full creeping conditions. On comparing the historical record in the seismic catalog, tsunami records, and the very-low frequency earthquake distribution, one of the areas has a plausible tsunamigenic potential for an M∼8 event. The spatial pattern of seismicity shows obvious variations between the northern and southern halves of the study area; intraslab earthquakes were dominant in the northern half, whereas interplate earthquakes were dominant in the southern half. Furthermore, the lower limit of the seismogenic zone in the northern half was ∼5–10 km deeper than that in the southern half. These regionalities indicate that differences in age and/or temperature or degree of hydration in the north–south direction exist in the slab.

Geophysical Modelling of the Middle America Trench using GMT

The study is focused on the geomorphological analysis of the Guatemala Trench, East Pacific Ocean. Research goal is to find geometric variations in western and eastern flanks of the trench and correlation of the submarine geomorphology with geologic settings and seismicity through numerical and graphical modelling. Methods include GMT based analysis of the bathymetry, geomorphic shape and surface trends in topography and gravity grids. Dataset contains raster grids on bathymetry, gravity, geoid and geological layers. Technical workflow is following: 1) Bathymetric mapping by modules ('grdcut', 'grdimage'); 2) Datasets visualizing and analysis, 3) Topographic and gravimetric surface modelling by ASCII data; 4) Cartographic mapping ('psbasemap', 'psxy', 'grdcontour') 4) 3D-mesh modelling; 5) Automatically digitized orthogonal cross-stacked profiles (‘grdtrack’); 6) Visualizing curvature trends (‘trend1d’); 7) statistical histograms. Results reveal unevenness in the structure of the submarine landforms. Modelling cross-section profiles highlighted depth variation at different parts of the transects and seafloor segments. Geomorphic structure has straight shape form of the slopes with steep oceanward forearc. Its geometry has steep and strait shape which correlates with seismicity. Depth samples vary: -3000 to -6200 m, seafloor is 3-5 km wide. Trench has symmetric accurate 'V- shaped' geometric form for the segment of -30 to -30 km. Oceanward side slope increases towards continental shelf, left flank stabilizes in depths at -50 km from the trench axis, deepens at -50 to 0 m, decreases at -3,200 to -5,800. Oceanward flank gradient varies: slope steepness of 35,63° at 0-20 km, 42,17° at 20-40 km, 44° at 40-60km. Left flank has 43,24° at 0-24 km, 28,33° at 24-40 km, 14° at 40-60 km, 1,13° at 60-200 km. Marine vertical gravity correlates with tectonic slab contours (20-30 mGal on the trench slopes). The study demonstrated effective GMT-based framework with a multi-disciplinary scope that combined cartographic methods of modelling deep-sea trench profiles with geological, bathymetric and geophysical analysis. Technical application of the advanced cartographic solutions and high quality mapping by GMT demonstrated its functionality for data analysis and effectiveness of the geoinformation processing. Proposed techniques and workflow methodology provided framework for geological mapping and can be applied in further similar research.

Event Stratigraphy in a Hadal Oceanic Trench: The Japan Trench as Sedimentary Archive Recording Recurrent Giant Subduction Zone Earthquakes and Their Role in Organic Carbon Export to the Deep Sea

Hadal trenches are the deepest places on Earth and are important foci for natural carbon sequestration. Much of the sedimentary sequences that accumulate within hadal trenches have been linked to widespread slope sediment remobilisation events, triggered by subduction zone earthquakes. Therefore, hadal trench deposits may provide valuable insights into the hazards posed by large earthquakes and their implications for the carbon cycle. Despite this strong societal relevance, no studies to date have provided the necessary coverage to understand the spatial and temporal variations of earthquake-triggered deposition along a hadal trench axis. We address these issues by integrating high-resolution bathymetry and subbottom profiler data, and sediment cores acquired over the entire hadal trench axis of the Japan Trench. We identify around 40 isolated trench-fill basins along the trench axis of the Japan Trench that document 115 sediment remobilization event deposits. We map the spatio-temporal distribution of the acoustically-transparent event deposit bodies imaged in subbottom profiler data from the trench-fill basins. Using radiocarbon dating, slope failure deposits identified from subbotom profiles and sediment coring were shown to be co-eval with major historic earthquake (e.g., AD2011 Mw9.0–9.1 Tohoku-oki, AD1454 Mw≥8.4 Kyotoku, and AD869 Mw≥8.6 Jogan events). Furthermore, the lower part of the acoustically-imaged stratigraphic succession in isolated basins along the Japan Trench also documents several thick acoustically-transparent bodies that relate to older events. These identifications of event deposits allow quantitative constraints of along-strike variation of sediment volumes redistributed by episodic events along the entire trench axis, revealing that the total volumes of event deposits triggered by different historic large earthquakes are highly variable. We conclude that at least 7 Tg (1012 g) of organic carbon remobilized from surficial slope sediments is exported to the hadal axis of Japan Trench in the last 2,000 years by giant earthquakes. These findings highlight the significance of seismo-tectonic events for the long-term carbon cycle in hadal trenches and societal implications.

Intra-trench variations in flexural bending of the subducting Pacific Plate along the Tonga-Kermadec Trench

We conducted a detailed analysis of along-trench variations in the flexural bending of the subducting Pacific Plate at the Tonga-Kermadec Trench. Inversions were conducted to obtain best-fitting solutions of trench-axis loadings and variations in the effective elastic plate thickness for the analyzed flexural bending profiles. Results of the analyses revealed significant along-trench variations in plate flexural bending: the trench relief (W0) of 1.9 to 5.1 km; trench-axis vertical loading (V0) of −0.5×1012 to 2.2×1012 N/m; axial bending moment (M0) of 0.1×1017 to 2.2×1017 N; effective elastic plate thickness seaward of the outer-rise region (\(T_{\rm{e}}^{\rm{M}}\)) of 20 to 65 km, trench-ward of the outer-rise (\(T_{\rm{e}}^{\rm{m}}\)) of 11 to 33 km, and the transition distance of 20 to 95 km. The Horizon Deep, the second greatest trench depth in the world, has the greatest trench relief (w0 of 5.1 km) and trench-axis loading (V0 of 2.2×1012 N/m); these values are only slightly smaller than that of the Challenger Deep (W0 of 5.7 km and V0 of 2.9×1012 N/m) and similar to that of the Sirena Deep (W0 of 5.2 km and V0 of 2.0×1012 N/m) of the Mariana Trench, suggesting that these deeps are linked to great flexural bending of the subducting plates. Analyses using three independent methods, i.e., the \(T_{\rm{e}}^{\rm{M}}\)/\(T_{\rm{e}}^{\rm{m}}\) inversion, the flexural curvature/yield strength envelope analysis, and the elasto-plastic bending model with normal faults, all yielded similar average Te reduction of 28%–36% and average Te reduction area \({S_{\Delta {T_{\rm{e}}}}}\) of 1 195–1 402 km2 near the trench axis. The calculated brittle yield zone depth from the flexural curvature/yield strength envelope analysis is also consistent with the distribution of the observed normal faulting earthquakes. Comparisons of the Manila, Philippine, Tonga-Kermadec, Japan, and Mariana Trenches revealed that the average values of \(T_{\rm{e}}^{\rm{M}}\) and \(T_{\rm{e}}^{\rm{m}}\) both in general increase with the subducting plate age.

Petit-spot volcanoes on the oldest portion of the Pacific plate

The formation of a petit-spot volcano involves the ascent of an asthenospheric melt to the outboard part of a plate subduction zone where flexure of the plate has taken place in the outer trench swell. On the Pacific Plate, the only previously known examples of such volcanoes were in the vicinity of the Japan Trench. We describe here a newly discovered petit-spot submarine volcano that formed in one of the oldest parts of the Pacific Plate among a cluster of small conical knolls to the southeast of Minamitorishima (Marcus) Island. Geochronological data indicate that this petit-spot volcanic eruption occurred less than 3 million years ago. The volcano erupted on the eastern slope of the outer rise, and it is at an unusually large distance from the Mariana trench axis (~800 km) because this rise, near the older part of the Pacific Plate, is wider than those of other subduction systems.

Flexural bending curvature and yield zone of subducting plates

ABSTRACT We quantify flexural deformation of subducting oceanic plates at a global array of 15 ocean trenches, using a new approach of modelling spatial variations in flexural bending shape and curvature. The investigated trenches are chosen to represent a diverse range of subducting plate age of 24–150 Ma, including the Middle America, Peru, Chile, West and East Aleutian, Sumatra, North and South Philippine, Tonga, Kermadec, Kuril, Japan, Izu-Bonin, and South and North Mariana. The studied trenches show systematic intra- and inter-trench variability in the calculated flexural bending curvature, stress distribution, extensional brittle yield zone, and effective elastic plate thickness Te of the subducting plates. We find that subducting plate age is a critical factor controlling the bending curvature and the corresponding extensional yield zone. The width, depth, and area of the extensional yield zone are all calculated to increase systematically with the subducting plate age. The newly-developed curvature analysis can yield continuously varying apparent Te(x) from the trench axis to outer rise. The calculated extensional yield zones from the curvature analysis are in general consistent with the observed normal faulting earthquakes of magnitude≥6.0 at the 15 global trenches. Our analyses also reveal that the five deepest regions of the global trenches, i.e. the Challenger, Horizon, Serena, and Scholl Deeps and Galathea Depth, are associated with relatively large flexural bending and calculated yield zones comparing to their respective adjacent trench segments. The Serena and Challenger Deeps of the S. Mariana trench are calculated to have the largest flexural bending among the five deeps.

1D-velocity structure and seismotectonics of the Ecuadorian margin inferred from the 2016 Mw7.8 Pedernales aftershock sequence

On April 16th 2016 a Mw 7.8 earthquake ruptured the central coastal segment of the Ecuadorian subduction zone. Shortly after the earthquake, the Instituto Geofisico de la Escuela Politecnica Nacional of Ecuador, together with several international institutions deployed a dense, temporary seismic network to accurately categorize the post-seismic aftershock sequence. Instrumentation included short-period and broadband sensors, along with Ocean Bottom Seismometers. This deployment complemented the permanent Ecuadorian seismic network and recorded the developing aftershock sequence for a period of one year following the main-shock. A subset of 345 events with ML > 3.5, were manually picked in the period of May to August 2016, providing highly accurate P- and S-onset times. From this catalogue, a high-quality dataset of 227 events, with an azimuthal gap <200°, are simultaneously inverted for, obtaining the minimum 1D velocity model for the rupture region, along with hypocentral locations and station corrections. We observe an average Vp/Vs of 1.82 throughout the study region, with relatively higher Vp/Vs values of 1.95 and 2.18 observed for the shallowest layers down to 7.5 km. The high relative Vp/Vs ratio (1.93) of the deeper section, between 30 km and 40 km, is attributed to dehydration and serpentinization processes. For the relocated seismicity distribution, clusters of events align perpendicular to the trench, and crustal seismicity is also evidenced, along with earthquakes located close to the trench axis. We also compute Regional Moment Tensors to analyze the different sources of seismicity after the mainshock. Aside from thrust events related to the subduction process, normal and strike-slip mechanisms are detected. We suggest that the presence of subducting seamounts coming from the Carnegie Ridge act as erosional agents, helping to create a scenario which promotes locking and allows seismicity to extend up to the trench, along zones of weakness activated after large earthquakes.

Sediment delivery and depositional processes along the eastern Aleutian Trench

The Aleutian Trench contains vast volumes of clastic sediment sourced mainly from the glacial erosion of Alaska and controlled by allocyclic influences. GLORIA (Geologic Long Range Inclined Asdic) and multi-channel seismic surveys conducted by the US Geological Survey revealed details of both the geometry and distribution of the trench sediments and their associated supply systems. Milankovitch-induced climatic cycles cause the episodic release of sediment transported as both unconfined and confined turbid flows, the latter routed along submarine channels that cross the abyssal plain of the Gulf of Alaska and discharge into the trench. Small volumes of sediment enter the trench directly from glacier snouts at the transform trough-end of the trench and from the Chugach Mountains in the north. A major channel system extends along the trench axis, which acts as a conduit for turbid sediment flows. The channel also reworks and redistributes material supplied from the accretionary prism and the oceanic plate. Subduction of seamount chains coupled with the westward-increasing obliquity of plate collision and accretionary prism-sourced slump, which may be seismically triggered, segment the trench environment, dividing the axial channel system into discrete flow patterns. This effect of compartmentalization creates petrographically distinct segments of fill partly owing to the primary source landscape of accreted terranes formed of ancient island arcs. This serves to emphasize the cyclic nature of Earth processes.

Sources, Degradation, and Transport of Organic Matter in the New Britain Shelf‐Trench Continuum, Papua New Guinea

AbstractHadal trenches are considered as depocenters for organic matter and hotspots for microbial diagenetic activity. Here, we explore the sources, degradation, and transport of organic matter in the shelf‐trench continuum using seven short sediment cores collected along two transects with water depths ranging between 1,553 and 8,901 m in the New Britain Trench area, Papua New Guinea. Carbon isotopic compositions (δ13C) and radiocarbon contents (Δ14C) of sedimentary organic matter accompanied by total organic carbon/total nitrogen ratios suggest an important contribution from the preaged soil organic matter mixed by the marine algae and terrigenous C3 vascular plants. In addition, the trench axis sites are characterized by elevated accumulation of terrigenous organic materials. Rates of organic matter mineralization approximated by dissolved inorganic carbon fluxes at the sediment‐water interface reveal an approximately threefold higher rate at the trench axis sites than the abyssal sites. 210Pbxs profiles and burial of carbonate (up to 50%) at both trench axis sites reflect recent occurrence of mass‐wasting events possibly induced by earthquakes, which is responsible for the transport of preaged, terrigenous organic matter to the trench bottom. This tectonically triggered deposition event, which was also shown to occur in the Japan Trench, the Tonga Trench, and probably in many other trenches, is likely to efficiently transport terrigenous organic matter to hadal trenches, thereby underpinning the importance of terrigenous organic matter burial in hadal trenches for the ocean organic carbon budget. Furthermore, we hypothesize that hadal trenches may host a distinct microbial community that is capable of feeding on old, refractory terrigenous organic matter.

Crustal dense blocks in the fore-arc and arc region of Chilean ranges and their role in the magma ascent and composition: Breaking paradigms in the Andean metallogeny

We use Intermediate-to-short wavelength (20–150 km) gravity anomalies to investigate the role played by crustal scale dense blocks in the genesis of Cu ore deposits along the Andes. These dense blocks are widespread along the margin and located at the fore-arc, arc, and back-arc domains, including depths in the range of 12–18 km, and thickness of 5–8 km. The interplay of these dense blocks with long-lived Trans-Lithospheric-Faults (TLF) controls the emplacement, and most likely, provides the required fertile fluids for the Cu ore genesis. Dense blocks represent impervious crustal domains, surrounded by damage zones at the intersection with TLF systems. Concentration of seismicity, large Vp/Vs ratios, low gravity domains, active volcanism, and structural fabric, demonstrate that these damage zones are in fact the most likely places for magma as well as mineralized flow ascent and emplacement. Basic flow models in porous media show a flow path that migrates horizontally at the base of dense/impervious blocks until the high permeability edge flanks, where pressure gradients provide the ideal conditions for the upward and focusing hydrothermal fluid ascent. Although not fully addressed in this paper, evidence suggests that the interaction of subducting-related fluids and the mafic-dense block might produce a fertile fluid with a concentration of Cu.In contrast with the long-standing view of metallogenic belts of the same age and parallel to the trench axis, we propose a multi-age metallogenic cluster, along the periphery of dense blocks, and organized oblique to the margin, following the trend of a given TLF.