Trench Subduction Research Articles

Sedimentology and provenance of newly identified Upper Cretaceous trench basin strata, Dênggar, southern Tibet: Implications for development of the Eurasian margin prior to India–Asia collision

AbstractTrench basins preserved along the southern margin of the Lhasa Terrane, Tibet, are sedimentologic records of convergent margin processes preceding Cenozoic India–Asia collision. We present new sedimentologic, petrographic and geochronologic data from the Rongmawa Formation and surrounding strata near Dênggar, Tibet, to determine depositional environment, provenance and age. Depositional ages range from ca. 92 to 87 Ma and lithofacies are consistent with deposition by low‐ and high‐density turbidity currents and suspension settling of pelagic detritus in a deep‐marine, trench basin setting. Sandstone modal analyses and U–Pb geochronology indicate that trench basin detritus in this region was derived from the Lhasa Terrane. We interpret that the Cretaceous subduction trench received detritus from an axial sediment dispersal system that transported sediment from headwaters in the central‐southern Lhasa terrane near Lhasa City directly to the trench and then flowed westwards parallel to the trench. The preservation of trench basin strata deposited during Late Cretaceous time compared with the lack of trench deposits prior to ca. 92 Ma and after ca. 80 Ma suggests the margin experienced a period of significant accretion during this interval. In addition, deposition of trench basin strata occurred during Late Cretaceous adakitic magmatism and high‐temperature metamorphism, which are hypothesized to be explained by subduction of an oceanic ridge or subduction zone retreat and related upper plate extension along the southern margin of the Lhasa terrane. Subduction of an oceanic ridge may provide a mechanism to potentially erode forearc basin strata and promote increased sediment delivery directly to the trench.

Forearc structures and deformation along the Manila Trench

The Manila Trench subduction zone is characterized by varying tectonic structures from north to south. Analyses of new seismic reflection data and bathymetric data indicate distinct morphological and deformational patterns in the forearc region. Differences in the nature of the subducting oceanic lithosphere (i.e. seafloor relief related to seamounts and ridges, sediment supply, reactivated features and faults associated with the South China Sea opening) cause along-strike heterogeneity in the Manila Trench and the Luzon forearc region. The northern segment is classified as an accretionary margin while the southern segment is mainly an erosive margin (with a narrow, steep, and often eroded frontal wedge). This sharp contrast is attributed to abundant sediment supply to the trench in the north and to the highly eroded frontal wedge in the south due to scarce sediment supply. The southern trench segment is prone to submarine slope failures and mass wasting processes. The 17°N latitude boundary also separates the forearc basin into the North Luzon Trough and the West Luzon Trough. Associated with this is the initiation of the Scarborough Seamount Chain (South China Sea extinct spreading ridge) subduction at 16°N latitude. A combination of forearc uplift and submarine mass movements attributed to subduction of bathymetric highs near and south of 17°N latitude produced the Stewart Bank. Seamount and other seafloor spreading features induced complex responses arising to diverse forearc architectures in the southern segment. Seamounts and other seafloor spreading-related features in the subducting slab induce slope steepening and significant vertical movement in the frontal wedge and the forearc region, respectively. Deformation associated with the subduction is overprinted by shearing related to the Philippine Fault Zone splay faults with the frontal wedge shortening associated with the ongoing subduction, further complicating the forearc development in the trench and marine forearc region.

Applying a Deep Learning Algorithm to Tsunami Inundation Database of Megathrust Earthquakes

AbstractWe explore recent developments in computer science on deep learning to estimate high‐resolution tsunami inundation from a quick low‐resolution computation result. Deep network architecture is capable of storing large information acquired via a training/learning process by pairing low‐ and high‐resolution deterministic simulation results from precalculated hypothetical scenarios. In a real case, with a real‐time source estimate and linear simulation computed on a relatively low‐grid resolution, optimized network parameters can be used to rapidly and accurately transform low‐resolution simulation outputs to higher‐resolution grids. We generate 532 source scenarios of interplate earthquakes (Mw 8–9) along the Japan Trench subduction zone to simulate tsunamis and utilize a precalculated library for deep learning. To test the proposed method, we consider a realistic source model inferred from static ground displacements and offshore tsunami data presumably available in real‐time and compare forecasted inundation heights against observations in Rikuzentakata and Otsuchi cities associated with the 2011 Tohoku‐oki tsunami. Our results show that the proposed method exhibits comparable accuracy to the conventional physics‐based simulation but achieves approximately 90% reduction of real‐time computational efforts. Thus, it has good potential as a future tsunami forecasting algorithm.

Evolution of the South-Iberian paleomargin: From hyperextension to continental subduction

Three crustal-scale sequentially restored cross-sections along the Central External Betic Cordillera enable us to pinpoint the role of extreme crustal thinning, tectonic inheritance and salt tectonics during the evolution of the South-Iberian paleomargin from Mesozoic rifting to Cenozoic continental subduction. This interpretation puts forward the existence of a former highly extended Mesozoic rift system, namely the Subbetic Basin. The basin comprises three main domains that coincide with the classical zoning of the External Betics: (i) The Intermediate units in the northern basin margin represent primary thick minibasins separated by salt walls; (ii) the External Subbetic units attest an inflated salt body in a central trough; and (iii) the Median and Internal Subbetic units in the southern basin margin are analogous to a group of primary minibasins with intercalations of submarine volcanic rocks overlying hyperextended crust. Shortening led to the extrusion of a salt canopy and the development of secondary minibasins. Furthermore, the tectonic inversion of the Subbetic Basin created two major thrusts sheets detached at the Upper Triassic evaporites, accommodating a shortening of ≈100–145 km at early Eocene-middle Miocene time. Compression reshaped the hyperextended domain into a subduction trench, and the former transfer faults into tear faults. Ongoing shortening led to the subduction of the South-Iberian paleomargin beneath the Alborán Domain lithosphere. Our results provide a new perspective on the inversion of segmented rift basins partly floored by mantle rocks and call for a major reconsideration on the tectonic architecture and evolution of the Betic-Rif Cordillera.

Geochronological, Geochemical and Sr–Nd–Pb Isotope Characteristics of the Meydan Ophiolite, SE Turkey: Petrogenesis and Implications for Mesozoic Tectonic Evolution

Gabbros, sheeted dikes and volcanic rocks of the Meydan ophiolite in the South-east Anatolian orogenic belt were examined to determine their U–Pb zircon ages, whole-rock, mineral chemistry and Nd‒Sr–Pb isotope compositions. The new U–Pb dating of zircon separates yielded a 88.3 ± 1.3 Ma crystallization age for gabbroic rocks of Meydan ophiolite. The ophiolite includes tectonites, ultramafic-mafic cumulates, isotropic gabbros, isolated mafic dikes, sheeted dikes, plagiogranites and rare basalts. The sheeted dikes and volcanics of the Meydan ophiolite exhibit near-flat rare earth element (REE) patterns, flat high field strength element trends and negative Nb anomalies that are similar to other SSZ–type eastern Mediterranean ophiolites. The whole-rock and mineral chemistry of the Meydan ophiolite/volcanic rocks suggests that they formed a supra-subduction zone and the dismembered ophiolite metamorphic sole formed beneath the ophiolite in a subduction trench setting. Isotopic ratios, i.e., 87Sr/86Sr versus 206Pb/204Pb and 143Nd/144Nd versus 87Sr/86Sr, show only limited variations, whereas 143Nd/144Nd versus 206Pb/204Pb, 87Sr/86Sr ratios are much more variable. The isotopic ratios are indicative of mixed magmatic contributions from DMM/MORB and EMII reservoirs that probably reflect the overall recycling of mantle wedge, pelagic sediment, oceanic crust and continentally derived sediment. The Meydan ophiolite, inferred arc volcanics and accretionary melange are all intruded by Late Cretaceous (87.7 ± 0.47 to 80.1 ± 0.46 Ma; zircon U–Pb) granitic rocks which are evaluated as the product of continental margin arc magmatism. The available regional comparisons suggest that the rocks studied were derived from the Berit Ocean to the north of the Bitlis Puturge continent, in which the Komurhan, Ispendere and Goksun ophiolites also formed during the Late Cretaceous.

Hypocenter and Magnitude Analysis of Aftershocks of the 2018 Lombok, Indonesia, Earthquakes Using Local Seismographic Networks

Abstract The island of Lombok in Indonesia is located between the Indo-Australian and Eurasian subduction trenches and the Flores back-arc thrust, making it vulnerable to earthquakes. On 29 July 2018, a significant earthquake Mw 6.4 shook this region and was followed by series of major earthquakes (Mw>5.8) on 5, 9, and 19 August, which led to severe damage in the northern Lombok area. In this study, we attempt to reveal the possible cause of the sequences of the 2018 Lombok earthquakes based on aftershock monitoring data. Twenty stations were deployed to record earthquake waveform data from 4 August to 9 September 2018. In total, 3259 events were identified using 28,728 P- and 20,713 S-wave arrival times during the monitoring. The aftershock hypocenters were determined using a nonlinear approach and relocated using double-difference method. The moment magnitude (Mw) of each event was determined by fitting the displacement spectrum amplitude using a Brune-type model. The magnitudes of the aftershocks range from Mw 1.7 to 6.7. The seismicity pattern reveals three clusters located in the Flores oceanic crust, which fit well with the occurrences of the four events with Mw>6. We interpret these events as the main rupture area of the 2018 Lombok earthquake sequence. Furthermore, an aseismic zone in the vicinity of Rinjani extending toward the northwestern part of Lombok was observed. We propose that the crust in this area has elevated temperatures and is highly fractured thus inhibiting the generation of large earthquakes. The aseismic nature is therefore an artifact of the detection threshold of our network (Mw 4.6).

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.

Earthquake slip surfaces identified by biomarker thermal maturity within the 2011 Tohoku-Oki earthquake fault zone

Extreme slip at shallow depths on subduction zone faults is a primary contributor to tsunami generation by earthquakes. Improving earthquake and tsunami risk assessment requires understanding the material and structural conditions that favor earthquake propagation to the trench. We use new biomarker thermal maturity indicators to identify seismic faults in drill core recovered from the Japan Trench subduction zone, which hosted 50 m of shallow slip during the Mw9.1 2011 Tohoku-Oki earthquake. Our results show that multiple faults have hosted earthquakes with displacement ≥ 10 m, and each could have hosted many great earthquakes, illustrating an extensive history of great earthquake seismicity that caused large shallow slip. We find that lithologic contrasts in frictional properties do not necessarily determine the likelihood of large shallow slip or seismic hazard.

Structural-morphological and sedimentary features of forearc slope off Miyagi, NE Japan: implications for development of forearc basins and plumbing systems

Multibeam (MB) and subbottom profile (SBP) data along the forearc slope off Miyagi, Japan, are combined to investigate structural-morphological and sedimentary structures along the forearc slope of the Japan Trench subduction zone. In addition to dip-oriented slope gullies, the MB data image a nearly dip-perpendicular slope trough bounded by a fault scarp landward of the trench-slope break. Seaward of the trench-slope break, the subbottom mostly contains normal faults that dip in opposite directions. The SBP data show not only unconformity and sliding surfaces, but also underfilled (forearc trough) and filled structures that may reflect the most recent forearc subsidence and basin filling. We propose these forearc trough and associated filling structures may indicate the earliest developments of modern isolated basins. Thus, a model of forearc basin development from slope gully to slope trough, isolated basins and forearc basin, transferred by active structures, could be proposed. Observed seafloor liquefaction in the SBP data may represent a near-surface seep structure of a plumbing system that enables hydrosphere-mantle fluid migration plausibly activated by earthquake cycles.

Subduction of an extinct rift and its role in the formation of submarine landslides in NW South America

Abstract On the eastern margin of the Panama Basin, the Nazca oceanic plate converges towards the continental plate of South America at approximately 53 mm a −1 . Subduction processes are accompanied by the presence of anomalous bathymetric elements including the Sandra Ridge. This east–west-orientated ridge is catalogued as an aborted rift derived from a magmatic spreading axis that was active between 12 and 9 Ma. Seismic activity within this structure is considered evidence of fault reactivation and tectonism. Once the structure reached the subduction trench several submarine landslides were triggered. Run-out lengths of these submarine landslides are perpendicular to the convergence of the structure with some units spreading and forming a wide fan that reaches tens of kilometres to the north and south of the trench. The area affected by the three main landslides varies between 130 and 300 km 2 approximately, with relatively superficial earthquakes (&lt;33 km) and with magnitudes that reach up to M w 7.2. The morphology of the landslides suggests a retrogressive nature with younger events proximal to shore. This paper presents estimates of the age of these landslides and discusses sources of uncertainty regarding these times of occurrence.

Subduction-driven heterogeneity of the lithospheric mantle beneath the Cathaysia block, South China

To evaluate the modification of the lithospheric mantle by the Paleo-Pacific Ocean subduction beneath the Cathaysia block, we report major and trace elements, and Li isotopes of olivine from Niutoushan, Mingxi, Xilong and Jiande peridotite xenoliths in South China. Spoon-shaped, convex-upward and highly LREE-enriched chondrite-normalized REE patterns of clinopyroxene in the Niutoushan, Mingxi and Xilong xenoliths reflect the interplay between partial melting and subsequent enrichment by metasomatism. Lithium isotope systematics indicate that at least two different metasomatic events affected these peridotites. High Li concentrations in clinopyroxene (up to 59.1 ppm) are ascribed to diffusive Li influx from silicate melts. Olivine grains with lower diffusivity of Li have much higher Li concentrations (0.56–9.76 ppm) and lighter δ7Li values (−52.75 to +10.49‰) than those in the Jiande and other peridotite xenoliths worldwide, showing a prominent carbonatite-metasomatic signature occurring before silicate metasomatism. In addition, increasing Li concentrations of olivine and clinopyroxene with increasing distance from the present subduction trench is consistent with those of low δ26Mg basalts in South China, indicating that the metasomatic agent was related to the westward subduction of the Paleo-Pacific slab. The Jiande lherzolite xenoliths further away from the trench show ‘mantle-like’ Li features and may represent upwelling asthenospheric materials accreted to the continent. Therefore, both types of metasomatism as well as upwelling asthenosphere could have played an important role in the modification of the chemical composition of the mantle beneath the Cathaysia block.

The influence of faulting style on the size-distribution of global earthquakes

We derive a unifying formulation, reliable at all scales, linking Anderson's faulting theory with the earthquake size-distribution, whose exponent is known as the b-value. Anderson's theory, introduced in 1905, related fault orientation to stress conditions. Independently, laboratory measurements on acoustic emissions have established that the applied differential stress controls their b-value. Our global survey revealed that observed spatial variations of b are controlled by different stress regimes, generally being lower in compressional (subduction trenches and continental collisional systems) and higher in extensional regimes (oceanic ridges). This confirmed previous observations that the b-value depends on the rake angle of focal mechanisms. Using a new plunge/dip-angles-based b-value analysis, we also identified further systematic influences of faulting geometry: steep normal faults (also typical of the oldest subduction zones) experience the highest proportion of smaller events, while low-angle thrust faults (typical of youngest subduction zones) undergo proportionally larger, more hazardous, events, differently from what would be expected by only allowing for rake-angle dependency. To date, however, no physical model has ever been proposed to explain how earthquakes size-distribution, differential stress and faulting styles relate to each other. Here, we propose and analytically derive a unifying formulation for describing how fault orientation and differential stresses determine b-value. Our formulation confirms that b-values decay linearly with increasing differential stress, but it also predicts a different dip-dependent modulation according to the tectonic environment, opening up new ways of assessing a region's seismic hazard.

Maximum magnitude of subduction earthquakes along the Japan-Kuril-Kamchatka trench estimated from seismic moment conservation

SUMMARYWe estimated the maximum magnitude of earthquakes in the Japan-Kuril-Kamchatka trench subduction zone with a method based on the conservation of seismic moment and the record of interplate seismicity from 1977 to 2017. The key point of this method is to base calculations on the tectonic moment rate instead of the total seismic moment rate. We modeled a seismic-moment-frequency distribution for the Japan-Kuril-Kamchatka trench on the basis of the truncated Gutenberg–Richter (G–R) law, the formula published by Utsu in 1974, the gamma distribution, and the tapered G–R law. We estimated the maximum magnitude along the Japan-Kuril-Kamchatka trench as ∼10 under the truncated G–R law and ∼11 under Utsu's formula, although the latter may be an overestimate. Therefore, the 2011 Tohoku earthquake, of moment magnitude 9.2, may not be the largest possible event in this area. The recurrence interval for magnitude 10 events based on the truncated G–R law is 4000 yr. Although these two models perform equally well in terms of Akaike Information Criterion, the range of the 95 per cent confidence level is consistently narrower for the truncated G–R law than for Utsu's formula. The estimated maximum magnitude depends not only on the model used, but also on the parameters that constitute the tectonic moment. It is essential to accumulate more seismic data and achieve more precise estimates of tectonic moment to improve estimates of maximum magnitude.

The sediment‐fill of Pago Pago Bay (Tutuila Island, American Samoa): New insights on the sediment record of past tsunamis

AbstractExtensive bathymetric and two‐dimensional seismic surveys have been carried out and cores collected in Pago Pago Bay (Tutuila, American Samoa) in order to describe and gain a better understanding of the sediment fill of the bay, which was affected by the 2009 South Pacific Tsunami. Eight sedimentary units were identified over the volcanic bedrock. The basal transgressive unit displays retrograding onlaps towards the shore, whereas the overlying seven aggradational layers alternate between four draping units and three pinching out seaward units. ‘Core to seismic’ correlation reveals that draping units are composed of homogeneous silts, while pinching out units are dominated by very coarse coral fragments showing fresh cuts, mixed with Halimeda plates. The basal unit is attributed to transgressive sedimentation in response to flooding of the bay after the last glacial maximum, followed by the upper aggradational units corresponding to highstand sedimentation. The changeovers in these upper units indicate an alternation between low‐energy silt units and high‐energy coral debris units interpreted as tsunami‐induced deposits. The 14C dating reveals that high‐energy sedimentation units can last up to approximately 2000 years while low‐energy sedimentation units can last up to approximately 1000 years. This alternation, deposited during the last highstand, may be explained by cycles of tectonic activity and quiescence of the Tonga Trench subduction, which is the main source of tsunamigenic earthquakes impacting the Samoan archipelago. In the uppermost silt unit, only the geochemical signature of the terrestrial input of the 2009 SPT backwash deposits was detected between 7 cm and 9 cm depth. Hence, Pago Pago Bay offers a unique sediment record of Holocene bay‐fill under the impact of past tsunamis intermittently during the last 7000 years.

Vertical motion of Phuket Island (1994–2018) due to the Sumatra-Andaman mega-thrust earthquake cycle: impact on sea-level and consequences for coral reefs

Macro-tidal coral reefs are particularly sensitive to medium to long-term changes in sea-level. Vertical motions of the seabed contribute to both lower or higher relative sea-level changes, particularly in tectonic plate boundary deformation zones along active subduction trenches. Phuket Island in Southern Thailand is subject to both horizontal and vertical land deformations during the seismic cycle of mega thrust earthquakes along the Sumatra and Andaman trenches. The relative sea-level changes in this region were historically monitored using the space geodetic techniques GPS and satellite altimetry alongside the traditional tide-gauge measurements over a period of almost 25 years. The GPS results show that the south of Thailand is still undergoing post-seismic deformations from the 2004 Mw 9.2 Sumatra-Andaman earthquake, after a significant change in the vertical motion of Phuket: from stable quasi-linear uplift at 2.5 ± 0.2 mm/yr, to (temporary) non-linear subsidence rates of 5–10 mm/yr in the past 14 quake aftermath years. The satellite altimetry data estimates the absolute sea-level rise in the Andaman Sea around the island at 3.9 ± 0.5 mm/yr. Therefore relative sea-level changes in Phuket appear to have been slightly positive until the end of 2004, followed by a significant increase averaging ~9 mm/yr that accumulated in 12 ± 1 cm by the end of 2018. Tide-gauge data gives a very similar result. The increased sea-level rise from 2005 onwards, due to tectonic land subsidence, correlates with the enhanced shallow coral reef growth and swift recoveries from bleaching events reported by long-term coral studies in this area.

Shallow Low‐Frequency Tremor in the Northern Japan Trench Subduction Zone

AbstractWe found offshore shallow low‐frequency tremor in the northern Japan Trench subduction zone, using records from the Seafloor Observation Network for Earthquakes and Tsunamis along the Japan Trench. The tremor aligned with a narrow zone, which follows the 10‐ to 25‐km depth contours of the plate interface and occurred in two distinct areas, off Tokachi and Sanriku, separated by a gap. In the 1.75 years examined, two episodes lasting for a few weeks were identified off Tokachi. During these episodes, tremor occurred concurrently with swarms of very low frequency earthquakes and migrated laterally about 120 km along the trench strike. In contrast, tremor off Sanriku lasted shorter (half a day to 5 days) and repeated at shorter intervals (1–4 months). The variation in tremor behavior such as duration, recurrence, migration, and the association with very low frequency earthquakes suggests along‐strike heterogeneities in frictional properties at the shallow plate interface.

Cretaceous termination of subduction at the Zealandia margin of Gondwana: The view from the paleo-trench

During the Permian to Cretaceous, Zealandia occupied a position on the proto-Pacific-facing, convergent margin of Gondwana. Subduction on this margin ceased somewhere between ~105 Ma and perhaps 70 Ma, but the timing of this tectonic transition remains controversial. Resolution of this uncertainty is important for tectonic reconstructions of the southwest Pacific and for global plate-tectonic models. Here, we revisit the problem by reference to new stratigraphic and geochemical data from the East Coast Basin of New Zealand, which occupied a position adjacent to and possibly superimposed on the relict Cretaceous subduction trench at the time subduction ceased; this basin is expected to preserve direct structural and stratigraphic evidence for or against Late Cretaceous subduction.In western parts of the East Coast Basin – the “Western Sub-belt” – a conspicuous unconformity separates undoubted accretionary prism rocks of the Torlesse Composite Terrane from younger Cretaceous “cover” units. Strata beneath and overlying this unconformity vary in age from place to place, but abundant paleontological data and published ages of detrital zircons (some reinterpreted herein), indicate that exposed Torlesse rocks are nowhere younger than ~100 Ma. Overlying Zealandia Megasequence “cover” strata are mostly younger than ~110 Ma. Between ~110 and ~85 Ma, the entire length of the Western Sub-belt reveals complex stratigraphic architecture of relatively small-scale depocentres subjected to alternating episodes of subsidence and local uplift and erosion. There is widespread evidence for compression on the Western Sub-belt over this period, although the overall amount of shortening is likely to be relatively modest. In contrast, the Eastern Sub-belt preserves a record of near-continuous and apparently simple deposition over the same interval of time. We see little and somewhat equivocal evidence for extension in either sub-belt through the Late Cretaceous. Superimposed on this general pattern, we observe discrete, widespread or basin-wide, tectonic events at 95 Ma, and within the intervals ~86–83 Ma and ~83–81 Ma, indicated by the presence of unconformities and sedimentary facies changes. Importantly, all these events can be correlated with coeval events recorded elsewhere in Zealandia, suggesting that the East Coast Basin was subject to the same tectonic regime as the rest of Zealandia and shared a common history during the mid- to Late Cretaceous.Integrating these findings with data from elsewhere in Zealandia, we argue that there is diverse evidence to indicate that subduction finished along the New Zealand segment of the Gondwana margin by 100 Ma. That said, the causes of on-going, Late Cretaceous compression in the East Coast Basin and in other parts of Zealandia are not well constrained, although several possible explanations are plausible. There are few analogue, abandoned subduction margins described in the literature; however, the situation in the South Shetland Trench shows strong similarities with patterns and complexities observed in the East Coast Basin.

The Thermal Effects of Plate‐Bending‐Related Thickening of the Oceanic Crustal Aquifer in the Nankai Trough and Japan Trench Subduction Zones

AbstractFor the Nankai Trough and Japan Trench subduction zones, we examine the combined thermal effects of lateral heat exchange by fluid circulation in an oceanic crustal aquifer and the thickening of that aquifer due to plate‐bending‐related faulting. Faults induced by the bending of a plate entering a subduction zone are hypothesized to increase the depth over which vigorous hydrothermal circulation can redistribute heat in the oceanic crust. Previous 1‐D (vertical) thermal models have examined how aquifer thickening can mine heat from deep in the crust seaward of the trench. Here we construct 2‐D thermal models that include aquifer thickening and lateral heat exchange in the aquifer. We vary the maximum aquifer thickness and the landward extent to which hydrothermal circulation persists within the subducted crust. For the Nankai margin, models most consistent with heat flux data require vigorous fluid circulation extending up to 150 km landward of the trench; aquifer thickening is permitted but not required for models to be consistent with the heat flux observations. Conversely, for the Japan Trench, preferred models include substantial aquifer thickening (maximum aquifer thickness of 1.8–5 km); vigorous circulation extending landward of the trench is permitted but is not required. For hot subduction zones, the thermal effects of aquifer thickening are modest relative to the large lateral advective heat redistribution. For cold subduction zones, where small lateral temperature gradients limit the amount of lateral heat redistribution, aquifer thickening can be the dominant process generating thermal anomalies.

Chapter 10 Sedimentary geochemistry used to infer the provenance of Permian–Triassic marine sandstones related to the SE Gondwana active continental margin, South Island, New Zealand

Abstract Major, trace and rare earth element data for sandstones and conglomerates from the Mid-Permian–Mid-Triassic Maitai Group are compared with other tectonostratigraphic units, using discrimination diagrams and comparisons with potential source terranes. Maitai Group sandstones reveal a mainly ophiolitic–oceanic-arc source during the Mid-Permian, followed by a mixed continental margin-arc–terrigenous source during the Late Permian. Latest Permian–Early Triassic sandstones mainly came from little-evolved continental margin-arc extrusives, tending to more evolved (but variable) during the Triassic. Source volcanism of the Murihiku Terrane sandstones was magmatically evolved relative to the Maitai Group generally (except during the Late Triassic). The Maitai Group and Murihiku Terrane are restored as proximal and more distal parts, respectively, of the SE Gondwana forearc basin. The localized Willsher Group shows some Maitai Group affinities. Sandstones in two melanges that formed in an outer forearc–subduction trench setting mainly indicate a mixed terrigenous–continental margin-arc source, similar to the Late Permian Maitai Group. The Caples Terrane, a Triassic accretionary prism, received detritus from little-evolved, to evolved continental margin-arc volcanics and terrigenous sources. Much of the arc-related material in all units is compatible with derivation from the latest Permian–Triassic Median Batholith, or a lateral equivalent along the SE Gondwana active margin.

Chapter 7 Patuki and Croisilles melanges in South Island, New Zealand: genesis related to Permian subduction–accretion processes

Abstract The late Early Permian ( c. 278–270 Ma) supra-subduction zone (SSZ) Dun Mountain ophiolite is bordered to the east by the Pataki Melange and to the NE by the Croisilles Melange. In the south, the ophiolite passes into a dismembered incipient oceanic arc (Otama Complex). The above units represent an oceanic forearc generated above a west-dipping subduction zone. Terrigenous sediment reached the subduction trench after the Mid-Permian(?) docking of the oceanic forearc with the long-lived SE Gondwana active continental margin. Mixed terrigenous–volcaniclastic turbidites accumulated in the trench prior to and during melange accretion. Fragments of the overriding oceanic forearc (and incipient arc, locally) detached and mixed to form melange and broken formation. Despite some individual features (e.g. of the basalt chemistry), the Patuki and Croisilles melanges are interpreted as originally representing a single Permian trench–accretionary complex. The more distal (easterly) part was sliced into the adjacent accretionary complex of the Caples Terrane to form the Croisilles Melange (and equivalent Greenstone Melange) probably after the Triassic. The South Island melanges exemplify accretionary processes in which igneous and sedimentary rocks were detached from the overriding plate by subduction–erosion, together with accretion, including seamount material from the subducting oceanic plate, with implications for melanges elsewhere.