Forearc Region Research Articles

Plate boundary reorganization in the active Banda Arc–continent collision: Insights from new GPS measurements

New GPS measurements reveal that large sections of the SE Asian Plate are progressively accreting to the edge of the Australian continent by distribution of strain away from the deformation front to forearc and backarc plate boundary segments. The study was designed to investigate relative motions across suspected plate boundary segments in the transition from subduction to collision. The oblique nature of the collision provides a way to quantify the spatial and temporal distribution of strain from the deformation front to the back arc. The 12 sites we measured from Bali to Timor included some from an earlier study and 7 additional stations, which extended the epoch of observation to ten years at many sites. The resulting GPS velocity field delineates at least three Sunda Arc–forearc regions around 500 km in strike-length that shows different amounts of coupling to the Australian Plate. Movement of these regions relative to SE Asia increases from 21% to 41% to 63% eastward toward the most advanced stages of collision. The regions are bounded by the deformation front to the south, the Flores–Wetar backarc thrust system to the north, and poorly defined structures on the sides. The suture zone between the NW Australian continental margin and the Sunda–Banda Arcs is still evolving with more than 20 mm/yr of movement measured across the Timor Trough deformation front between Timor and Australia.

Shocks in a box: An analogue model of subduction earthquake cycles with application to seismotectonic forearc evolution

We introduce and test an experimental approach to simulate elastoplastic megathrust earthquake cycles using an analogue model and apply it to study the seismotectonic evolution of subduction zones. The quasi‐two‐dimensional analogue model features rate‐ and state‐dependent elastic‐frictional plastic and viscoelastic material properties and is scaled for gravity, inertia, elasticity, friction, and viscosity. The experiments are monitored with a high‐resolution strain analysis tool based on digital image correlation (particle imaging velocimetry, PIV), providing deformation time series comparable to seismologic, geodetic, and geologic observations. In order to separate elastic and nonelastic effects inherent the experimental deformation patterns, we integrate elastic dislocation modeling (EDM) into a hybrid approach: we use the analogue earthquake slip and interseismic locking distribution as EDM dislocation input and forward model the coseismic and interseismic elastic response. The residual, which remains when the EDM prediction is subtracted from the experimental deformation pattern, highlights the accumulation of permanent deformation in the model. The setup generates analogue earthquake sequences with realistic source mechanisms and elastic forearc response and recurrence patterns and reproduces principal earthquake scaling relations. By applying the model to an accretionary‐type plate margin, we demonstrate how strain localization at the rupture peripheries may lead to a seismotectonically segmented forearc, including a tectonically stable shelf and coastal high (<20% plate convergence accommodated by internal shortening) overlying the area of large megathrust earthquake slip. Fifty to 75% of plate convergence is accommodated by internal shortening in the slope region where earthquake slip tapers out toward the trench. The inner forearc region remains undeformed and represents a basin.

Mapping the mantle wedge and interplate thrust zone of the northeast Japan arc

We detected sP depth phases from three-component, short-period seismograms from the forearc earthquakes under the Pacific Ocean recorded by the dense seismic network on the Japan Islands. Then we use the sP depth-phase data together with P and S wave arrival times to locate the suboceanic earthquakes accurately. A large number of P and S wave arrival times from the relocated suboceanic events and earthquakes under the land area are used to determine the 3-D P and S wave velocity and Poisson's ratio structures under the entire NE Japan arc from the Japan Trench to the Japan Sea. Our results exhibit strong lateral heterogeneity under the forearc region and a clear correlation between the structural heterogeneity and distribution of large earthquakes (M 7.0–8.2) occurred in the interplate thrust zone. Prominent low-velocity and high Poisson's ratio zones are visible in the mantle wedge under the volcanic front and back-arc region, which show nicely the images of arc magma caused by the slab dehydration and corner flow in the mantle wedge.

Mid-crustal horizontal shear zone in the forearc region of the mid-Cretaceous SW Japan arc, inferred from strain analysis of rocks within the Ryoke metamorphic belt

The strain geometry of a Cretaceous low- P/high- T metamorphosed accretionary complex in SW Japan is one of the key geological constraints in understanding the Cretaceous tectonics of the forearc region at the eastern margin of Eurasia. We studied the strain geometry of upper greenschist to lower amphibolite facies metamorphic rocks from the mid-Cretaceous low- P/high- T Ryoke metamorphic belt, which during the mid-Cretaceous were located at mid-crustal depths close to the volcanic front. Strain analysis focused on deformed radiolarian fossils in metachert and sandstone/chert clasts within metapelite. The maximum stretching direction ( X-direction) of the strain ellipsoid of the schistosity-forming deformation at the peak metamorphism is oriented E-W, and the XY-plane of the strain ellipsoid developed parallel to schistosity. The metachert and metapelite exhibit plane strain to general flattening strain that probably largely resulted from the schistosity-forming deformation. Based on the results of the previous strain analyses and this study, it reveals that the schistosity-forming deformation occurred under bedding-normal shortening, resulting in general flattening strain throughout the Ryoke metamorphic belt. The X-direction of the strain ellipsoid of the schistosity-forming deformation is oriented parallel to the length of the Ryoke metamorphic belt throughout the SW Japan arc. The schistosity formed parallel to bedding planes that developed horizontally in non-metamorphosed accretionary complexes of SW Japan, indicating horizontal XY-plane of the strain ellipsoid. These observations indicate the occurrence of horizontal shear deformation within the middle crust (∼15 km depth) at the time of the low- P/high- T metamorphism. The deep-seated horizontal shear zone and the shallower incipient Median Tectonic Line (MTL) fault system could result in a crustal-scale detachment fault in the forearc region of the mid-Cretaceous SW Japan arc. After the formation of the detachment, the left-lateral pull-apart forearc basins, Onogawa-Izumi Group distributed on the north side of the MTL and the southern marginal mylonite zone might be formed by the sinistral movement along the detachment fault.

Spatial heterogeneity of the mantle wedge structure and interplate coupling in the NE Japan forearc region

We performed 3‐D seismic tomography in the forearc region of the northeastern Japan subduction zone using both onshore and offshore seismic station data. We obtained the Vp, Vs, and Vp/Vs structures around the plate boundary with high spatial resolution. The position of the plate boundary as defined by relocated hypocenters coincides with the sharp velocity boundary between the oceanic crust and the mantle wedge. The mantle wedge above the coseismic slip area of the 1978 and 2005 off‐Miyagi interplate earthquakes (M > 7) is characterized by high Vp and Vs, but low Vp/Vs. Off Fukushima, however, where large earthquakes rarely occur, we found a high Vp/Vs anomaly at the tip of the mantle wedge. The spatial distribution of serpentinized mantle wedge limits the spatial extent of the strongly coupled area on the plate boundary, and thus can explain the difference in seismic activity between the off‐Miyagi and off‐Fukushima regions.

Pb isotope variations in hydrogenetic Fe–Mn crusts from the Izu–Bonin fore-arc

Fe–Mn crusts were recovered from the western escarpment of the Bonin Ridge in the Izu–Bonin fore-arc region (dive site #824: 28.612°N, 141.803°E) at water depths of c. 2900 m using the Shinkai 6500 submersible during cruise YK 04–05. Major and trace element data and XRD mineralogy indicate that the crusts are hydrogenetic in origin. We present profiles of variations in Pb isotope composition measured in-situ by laser ablation MC-ICP-MS across two of the crusts. The isotopic variations are systematic and can be matched up between the two crusts, indicating similar growth rates. The crust Pb isotope composition rules out any local source for Pb from within the Izu–Bonin–Mariana arc system, either from hydrothermal activity or through leaching of volcanic detritus. Input of a globally well-mixed volcanic Pb component, either from aerosols or as an absorbed component on aeolian dust, has been proposed as a mechanism to explain the Pb isotope composition of Central Pacific deep water. However, the Izu–Bonin crusts are displaced to lower 206Pb/ 204Pb and higher 208Pb/ 204Pb, which requires an additional Pb source. One possibility is that as water is advected from the south, outboard of the Luzon–Ryukyu–Honshu arc system, it is progressively polluted by Pb derived by weathering and erosion of these young island arc volcanic systems. Using a constant Co-flux model, growth rates are estimated at ∼ 7–13 mm/Ma, which would suggest that these crusts provide a record of changes in the composition of deep water in the Izu–Bonin fore-arc region of the western Pacific Ocean over the last 4–8 Ma. Over this interval, the main feature has been a progressive decrease in 207Pb/ 206Pb (0.843 to 0.839) and 208Pb/ 206Pb (2.088 to 2.080) with time. The interior parts have compositions similar to those of crusts from the Izu–Bonin fore-arc, while the rims have compositions similar to crusts from the central Western Pacific.

Perkembangan Tektonik Daerah Busur Muka Selat Sunda dan Hubungannya dengan Zona Sesar Sumatera

The subduction of the Indian-Australian Plate beneath the Eurasian Plate is oblique (~45o) along Sumatra. A major zone of dextral strike slip displacement along Sumatra Island, called the Sumatra Fault, partially accommodates the oblique convergence off Sumatra. On the contrary, the subduction is normal along Java; therefore, no major fault zone can be found along the Java Island. The transition zone of the Sunda Strait fore arc, which is located between Java and Sumatra fore arc, is subject to northwest-southeast extension related to the motion of the Sumatra Fault and north-south compression because of subduction. Geophysical studies show that continuous extension dominated the Sunda Strait fore arc region. The results are interpreted as showing ongoing separation of the Sunda Strait fore arc region as the Sumatra fore arc plate has moved northwest, bounded by the Sumatra Fault. Therefore, the Sumatra Fault can also be interpreted to extend across the fore arc to the trench in the form of several graben systems. Key words: Sunda Strait, Sumatra Fault Zone, fore arc tectonic

Fault-induced seismic anisotropy by hydration in subducting oceanic plates

In subduction zones, the anisotropic fast shear-wave component is generally observed to be orientated parallel to the strike of the trench. The interpretation of this shear-wave splitting above subduction zones has, however, been controversial and none of the inferred models appears to be explain the whole range of anisotropic patterns observed worldwide. Manuele Faccenda and colleagues now show that the amount and geometry of seismic anisotropy measured in the forearc regions of subduction zones strongly depend on the preferred orientation of hydrated faults in the subducting oceanic plate. The anisotropy originates from the crystallographic preferred orientation of highly anisotropic hydrous minerals formed along steeply dipping faults and from the larger scale vertical layering consisting of dry and hydrated crust-mantle sections whose spacing is several times smaller than teleseismic wavelengths. Faults orientations and estimated delay times are consistent with the observed shear-wave splitting patterns in most subduction zones. It is shown that the amount and geometry of seismic anisotropy measured in the forearc regions of subduction zones strongly depend on the preferred orientation of hydrated faults in the subducting oceanic plate. The anisotropy originates from the crystallographic preferred orientation of highly anisotropic hydrous minerals formed along steeply dipping faults and from the larger-scale vertical layering consisting of dry and hydrated crust–mantle sections, the spacing of which is several times smaller than teleseismic wavelengths. The variation of elastic-wave velocities as a function of the direction of propagation through the Earth’s interior is a widely documented phenomenon called seismic anisotropy. The geometry and amount of seismic anisotropy is generally estimated by measuring shear-wave splitting, which consists of determining the polarization direction of the fast shear-wave component and the time delay between the fast and slow, orthogonally polarized, waves. In subduction zones, the teleseismic fast shear-wave component is oriented generally parallel to the strike of the trench1, although a few exceptions have been reported (Cascadia2 and restricted areas of South America3,4). The interpretation of shear-wave splitting above subduction zones has been controversial and none of the inferred models seems to be sufficiently complete to explain the entire range of anisotropic patterns registered worldwide1. Here we show that the amount and the geometry of seismic anisotropies measured in the forearc regions of subduction zones strongly depend on the preferred orientation of hydrated faults in the subducting oceanic plate. The anisotropy originates from the crystallographic preferred orientation of highly anisotropic hydrous minerals (serpentine and talc) formed along steeply dipping faults and from the larger-scale vertical layering consisting of dry and hydrated crust–mantle sections whose spacing is several times smaller than teleseismic wavelengths. Fault orientations and estimated delay times are consistent with the observed shear-wave splitting patterns in most subduction zones.

Trench‐forearc interactions reflected in the sedimentary fill of Talara basin, northwest Peru

ABSTRACTExceptional exposure of the forearc region of NW Peru offers insight into evolving convergent margins. The sedimentary fill of the Talara basin spans the Cretaceous to the Eocene for an overall thickness of 9000 m and records within its stratigraphy the complicated history of plate interactions, subduction tectonics, terrane accretion, and Andean orogeny. By the early Tertiary, extensional tectonism was forming a complex horst and graben system that partitioned the basin into a series of localized depocentres. Eocene strata record temporal transitions from deltaic and fluvial to deep‐water depositional environments as a response to abrupt, tectonically controlled relative sea‐level changes across those depocentres. Stratigraphic and provenance data suggest a direct relationship between sedimentary packaging and regional tectonics, marked by changes in source terranes at major unconformities. A sharp shift is recognized at the onset of deepwater (bathyal) sedimentation of the Talara Formation, whose sediments reflect an increased influx of mafic material to the basin, likely related to the arc region. Although the modern topography of the Amotape Mountains partially isolates the Talara basin from the Lancones basin and the Andean Cordillera to the east, provenance data suggest that the Amotape Mountains were not always an obstacle for Cordilleran sediment dispersal. The mountain belt intermittently isolated the Talara basin from Andean‐related sediment throughout the early Tertiary, allowing arc‐related sediment to reach the basin only during periods of subsidence in the forearc region, probably related to plate rearrangement and/or seamounts colliding with the trench. Intraplate coupling and/or partial locking of subduction plates could be among the major causes behind shifts from contraction to extension (and enhanced subduction erosion) in the forearc region. Eventually, collisional tectonic and terrane accretion along the Ecuadorian margin forced a major late‐Eocene change in sediment dispersal.

Strong variations in seismic anisotropy across the Hikurangi subduction zone, North Island, New Zealand

Shear wave splitting measurements from SKS, SKKS, S and ScS phases reveal the anisotropic structure of North Island, New Zealand. Waveforms were recorded on temporary and permanent broadband seismographs. Strong anisotropy varies rapidly across the Hikurangi subduction zone. These changes constrain the width of the deformation associated with the plate boundary and define three different regions of mantle flow. Stations in the fore-arc region, eastern central North Island record fast polarisations subparallel to the trench and moderate delay times of around 2.5 s, consistent with earlier studies of North Island. This is interpreted as caused by trench-parallel mantle flow beneath the subducted slab with a possible contribution of trench-parallel fossil anisotropy within the subducted slab. Delay times increase towards the extending Central Volcanic Region (CVR) with values up to 4.5 s, one of the largest delay times measured worldwide. Strong S-wave anisotropy (7.5%) through 80 km in the mantle wedge in addition to the sub-slab anisotropy can explain the observed high delay times. One mechanism to create such high anisotropy is melt segregation of partially molten rock. A frequency dependence of this mechanism could also explain persisting discrepancies between local and teleseismic shear wave splitting measurements. To the west, delay times decrease systematically. Measurements in the far back-arc region of western central North Island show no apparent splitting. This suggests that mantle wedge dynamics terminate beneath the western border of the CVR. The apparent isotropy beneath the far back-arc may be related to local small-scale mantle convection or vertical mantle flow, possibly caused by the detachment of lithosphere.

Collision of the Ganges–Brahmaputra Delta with the Burma Arc: Implications for earthquake hazard

We take a fresh look at the topography, structure and seismicity of the Ganges–Brahmaputra Delta (GBD)–Burma Arc collision zone in order to reevaluate the nature of the accretionary prism and its seismic potential. The GBD, the world's largest delta, has been built from sediments eroded from the Himalayan collision. These sediments prograded the continental margin of the Indian subcontinent by ∼ 400 km, forming a huge sediment pile that is now entering the Burma Arc subduction zone. Subduction of oceanic lithosphere with > 20 km sediment thickness is fueling the growth of an active accretionary prism exposed on land. The prism starts at an apex south of the GBD shelf edge at ∼ 18°N and widens northwards to form a broad triangle that may be up to 300 km wide at its northern limit. The front of the prism is blind, buried by the GBD sediments. Thus, the deformation front extends 100 km west of the surface fold belt beneath the Comilla Tract, which is uplifted by 3–4 m relative to the delta. This accretionary prism has the lowest surface slope of any active subduction zone. The gradient of the prism is only ∼ 0.1°, rising to ∼ 0.5° in the forearc region to the east. This low slope is consistent with the high level of overpressure found in the subsurface, and indicates a very weak detachment. Since its onset, the collision of the GBD and Burma Arc has expanded westward at ∼ 2 cm/yr, and propagated southwards at ∼ 5 cm/yr. Seismic hazard in the GBD is largely unknown. Intermediate-size earthquakes are associated with surface ruptures and fold growth in the external part of the prism. However, the possibility of large subduction ruptures has not been accounted for, and may be higher than generally believed. Although sediment-clogged systems are thought to not be able to sustain the stresses and strain-weakening behavior required for great earthquakes, some of the largest known earthquakes have occurred in heavily-sedimented subduction zones. A large earthquake in 1762 ruptured ∼ 250 km of the southern part of the GBD, suggesting large earthquakes are possible there. A large, but poorly documented earthquake in 1548 damaged population centers at the northern and southern ends of the onshore prism, and is the only known candidate for a rupture of the plate boundary along the subaerial part of the GBD–Burma Arc collision zone.

The three-dimensional thermo-mechanical signature of ridge subduction and slab window migration

We present a series of three-dimensional numerical models investigating possible thermal and mechanical effects of ridge subduction and slab window migration. Ridge subduction is a fundamental process of plate tectonics, and the geologic manifestations in the overriding lithosphere have been investigated in many convergent margin settings. Many of the geologic effects of slab window migration (e.g., anomalous high- T metamorphism in the forearc, non-arc-like magmatism in the volcanic arc) are explained primarily by the introduction of a region of hot, upwelling asthenospheric mantle in the subduction zone environment. Using end-member boundary conditions and idealized geometries, our models address the thermal and corresponding mechanical manifestations of slab window migration, and our results provide a general agreement with observations from interpreted regions of slab window migration. Our models show that protracted heating in the forearc region should result in protracted high-temperature metamorphism, associated rheologic weakening and strain partitioning and changes in the topographic uplift pattern in regions affected by slab window migration. Although these models are idealized, they do provide valuable insight into the geodynamics of slab window environments and provide a valuable starting point for further exploration.

Seismic evidence for widespread serpentinized forearc mantle along the Mariana convergence margin

Seismic imaging of subduction zones can provide constraints on mineral reactions in the slab and surrounding regions. We use P‐to‐S converted phases from teleseisms recorded at broadband stations in the Mariana Islands to image the forearc and arc regions of the Mariana convergent margin. The subducting oceanic crust is observed between 75 and 110 km depth as a thin low velocity zone overlying the subducting Moho, demonstrating that the basalt to eclogite phase transition must occur at a greater depth. A low velocity zone (LVZ), approximately 10−25 km thick, whose upper boundary is imaged at about 40−55 km depth, is detected in the forearc region of the mantle wedge along the entire margin. The anomaly is located too shallow to represent subducted oceanic crust. We interpret the LVZ as a serpentinized region in the forearc mantle, resulting from hydration by slab‐expelled water. The inferred S‐wave velocity in the LVZ of as low as ∼3.6 km/s represents a level of serpentinization of 30−50%, corresponding to a chemically bound water content of about 4−6 wt%.

Differential Late Paleozoic active margin evolution in South-Central Chile (37°S–40°S) – the Lanalhue Fault Zone

The N–S oriented Coastal Cordillera of South Central Chile shows marked lithological contrasts along strike at ∼38°S. Here, the sinistral NW–SE-striking Lanalhue Fault Zone ( nomen novum) juxtaposes Permo-Carboniferous magmatic arc granitoids and associated, frontally accreted metasediments (Eastern Series) in the northeast with a Late Carboniferous to Triassic basal-accretionary forearc wedge complex (Western Series) in the southwest. The fault is interpreted as an initially ductile deformation zone with divergent character, located in the eastern flank of the basally growing, upwarping, and exhuming Western Series. It was later transformed and reactivated as a semiductile to brittle sinistral transform fault. Rb–Sr data and fluid inclusion studies of late-stage fault-related mineralizations revealed Early Permian ages between 280 and 270 Ma for fault activity, with subsequent minor erosion. Regionally, crystallization of arc intrusives and related metamorphism occurred between ∼306 and ∼286 Ma, preceded by early increments of convergence-related deformation. Basal Western Series accretion started at >290 Ma and lasted to ∼250 Ma. North of the Lanalhue fault, Late Paleozoic magmatic arc granitoids are nearly 100 km closer to the present day Andean trench than further south. We hypothesize that this marked difference in paleo-forearc width is due to an Early Permian period of subduction erosion north of 38°S, contrasting with ongoing accretion further south, which kinematically triggered the evolution of the Lanalhue Fault Zone. Permo-Triassic margin segmentation was due to differential forearc accretion and denudation characteristics, and is now expressed in contrasting lithologies and metamorphic signatures in todays Andean forearc region north and south of the Lanalhue Fault Zone.

Evidence of strong lateral inhomogeneous structure beneath SE Carpathians and specific mantle flow patterns

The Vrancea region, in the South-Eastern Carpathians (Romania), represents a unique case among the seismic areas in the world taking into account the extreme concentration and persistence of seismicity and tectonic stress field. Our goal is to show that the particular shape of attenuation and shear wave splitting properties can be interpreted in the light of the decoupling and slab retreat processes, which hypothetically induce a specific configuration of the upper mantle flow. Delamination and break-off processes combined with retrograde motion of the slab imply lateral asymmetry in flow geometry and geotectonic properties. Particularly relevant is the strike-parallel flow localized in front of the downgoing slab, in contrast with the steeply dip flow along the slab in the back side. The local upper mantle flows around the descending seismic active body explains the contrast of the seismic wave attenuation in the back-arc region against fore-arc region and the strong anisotropy anomaly observed in the South-East Carpathians (time delays of 1.5–2 s). Outside the epicentral area, the coherent pattern of the shear wave splitting follows the trends of the deformation field as outlined by the GPS measurements, in favour of a strong coupling between the surface and mantle processes.

Pb–Nd isotopic constraints on sedimentary input into the Lesser Antilles arc system

The Lesser Antilles arc is a particularly interesting island arc because it is presently very active, it is located perpendicular to the South American continent and its chemical and isotopic compositions display a strong north–south gradient. While the presence in the south of a thick pile of sedimentary material coming from the old South American continent has long been suspected to explain the geochemical gradient, previous studies failed to demonstrate unambiguously a direct link between the arc lava compositions and the subducted sediment compositions. Here, we present new Nd, Sm, Th, U and Pb concentrations and Nd–Pb isotopic data for over 60 sediments from three sites located in the fore arc region of the Lesser Antilles arc. New data for DSDP Site 543 drill core located east of Dominica Island complement the data published by White et al. [White, W.M., Dupré, B. and Vidal, P., 1985. Isotope and trace element geochemistry of sediments from the Barbados Ridge–Demerara Plain region, Atlantic Ocean. Geochimica et Cosmochimica Acta, 49: 1875–1886.] and confirm their relatively uniform isotopic compositions (i.e., 206Pb/ 204Pb between 19.13 and 19.53). In contrast, data obtained on DSDP Site 144 located further south, on the edge of the South American Rise and on sediments from Barbados Island are much more variable ( 206Pb/ 204Pb ranges from 18.81 to 27.69). The very radiogenic Pb isotopic compositions are found in a 60 m thick black shale unit, which has no age equivalent in the Site 543 drill core. We interpret the peculiar composition of the southern sediments as being due to two factors, (a) the proximity of the South American craton, which contributes coarse grain old detrital material that does not travel far from the continental shelf, and (b) the presence of older sediments including the thick black shale unit formed during Oceanic Anoxic events 2 and 3. The north–south isotopic change known along the Lesser Antilles arc can be explained by the observed geographical changes in the composition of the subducted sediments. About 1% contamination of the mantle wedge by Site 543 sediments explains the composition of the northern islands while up to 10% sediments like those of Site 144 is required in the source of the southern island lavas. The presence of black shales in the subducted pile provides a satisfactory explanation for the very low Δ8/4 values that characterize the Lesser Antilles arc.

Neogene subsidence and stratigraphy of the Simeulue forearc basin, Northwest Sumatra

Interpretation of new multi-channel seismic data reveals the Neogene subsidence history of the Simeulue forearc basin, located off northwestern Sumatra. The trench-parallel elongated basin extends over 260 km in NW–SE direction and 100 km in SW–NE direction. Maximum water depth is about 1300 m and the Neogene sedimentary fill shows a thickness of up to 5 s two-way traveltime. We identified three stages of subsidence evolution after the formation of a regional basal Neogene unconformity which are characterized by changing accumulation space, sedimentation rates, source areas and tectonic phases. An Early/Mid Miocene stage is marked by occurrence of subsidence in several half grabens forming along the western border of the basin. During the Late Miocene/Pliocene this small scale subsidence changed to a steadily subsiding trench-parallel trough. The present setup of the forearc region under the influence of strike-slip faults due to oblique subduction is active at least since this time as evidenced by wrench faulting. At the end of this stage subsidence expanded significantly eastwards drowning a large carbonate platform that evolved in the then shallow southern and eastern parts of the basin. Presently, the central part of the Simeulue basin is subject to inversion, probably related to reactivation of Early/Mid Miocene half grabens.

Seismic velocity structure around the Hyuganada region, Southwest Japan, derived from seismic tomography using land and OBS data and its implications for interplate coupling and vertical crustal uplift

The Hyuganada region, a forearc region of Southwest Japan, is characterized by several interesting geological and geophysical features, i.e., significant aseismic crustal uplift of ∼120 m during the past ∼120 thousand years at the Miyazaki Plain, negative free-air gravity anomalies with the maximum magnitude of −130 mgal, and relatively less cohesive interplate coupling compared with that for off the Shikoku and Kii Peninsula. In order to examine the causes of these observations, we determined a detailed three-dimensional seismic velocity structure based on the seismic data observed by ocean bottom seismometers (OBS) and land stations. P- and S-wave tomographic velocity structures clearly indicate the subducting slab and also the zones of high Poisson's ratio at 25–35 km depth along the coastline of the northeastern part of the Hyuganada. The region with high Poisson's ratio may correspond to the serpentinized mantle wedge as suggested for other mantle wedges, and appears to be coincident with the zone for observed aseismic slips such as the slow-slip and after-slip events. Also, the detection may be related to a relatively weak interplate coupling in the Hyuganada region. The tomographic structures also indicate low velocity zones with a horizontal scale comparable to the Kyushu-Palau Ridge in and around the subducting slab. If we assume that the low velocity zones correspond to the subducted Kyushu-Palau Ridge, then the predicted gravity anomaly due to the density contrast between the low velocity zones and the surrounding region can explain about 60% of the gravity anomaly in the Hyuganada region. The buoyancy is probably an important factor for the crustal uplift observed in the Miyazaki Plain, the steep bending of the subducting slab and the normal fault-type earthquakes around the Hyuganada region.

Synthesis and tectonic interpretation of the westernmost Paleozoic Variscan orogen in southern Mexico: From rifted Rheic margin to active Pacific margin

Paleozoic rocks in southern Mexico occur in two terranes, Oaxaquia (Oaxacan Complex) and Mixteca (Acatlán Complex) that appear to record: (1) Ordovician rifting on the southern margin of the Rheic Ocean, (2) passive drifting with Amazonia during the Silurian, (3) Devonian-Permian subduction beneath southern Mexico producing an arc complex that was partially removed by subduction erosion, subjected to HP metamorphism and Mississippian extrusion into the upper plate, followed by reestablishment of a Permian arc. In the Oaxaquia terrane, the 920–1300 Ma basement is unconformably overlain by a ∼ 200 m uppermost Cambrian–lowest Ordovician shelf sequence containing Gondwanan fauna (Tiñu Formation), unconformably overlain by 650 m of shallow marine-continental Carboniferous sedimentary rocks containing a Midcontinent (USA) fauna. In the Mixteca terrane, the low-grade Paleozoic sequence is composed of: (a) a ?Cambrian–Ordovician clastic sequence intruded by ca. 480–440 Ma bimodal, rift-related igneous rocks; and (b) a latest Devonian-Permian shallow marine sequence (> 906 m) consisting of metapsammites, metapelites and tholeiitic mafic volcanic rocks. High pressure (HP) metamorphic rocks in the Mixteca terrane consists of: (i) a Cambro-Ordovician rift-shelf intruded by bimodal rift-related intrusions that are similar to the low-grade rocks; (ii) periarc ultramafic rocks, and (iii) arc and MORB rocks. The Ordovician granitoids contain concordant inherited zircons that range in age from ca. 900 to 1300 Ma, indicating a source in the Oaxacan Complex. Concordant ages of detrital zircons in both the low- and high-grade Cambro-Ordovician metasedimentary rocks indicate a provenance in local Ordovician plutons and/or ca. 1 Ga Oaxacan basement, and distal northwestern Gondwana sources with a unique source in the 900–750 Ma Goiás magmatic arc within the Brasiliano orogen. These data combined with the rift-related nature of the Cambro-Ordovician rocks are most consistent with an origin along the southern margin of the Rheic Ocean. Latest Devonian-Permian deposition was synchronous with Mississippian extrusion of the HP rocks into the upper plate during extensional deformation. The HP Cambro-Orodivician rift-shelf rocks are inferred to have originated in the forearc region of the upper plate that was removed by subduction erosion, carried down the subduction channel, and then extruded into the upper plate in the middle of the Mixteca terrane. The presence of arc-related rocks in the HP assemblage suggests that the arc complex was also removed, whereas the MORB rocks may have been derived from the subducting slab. Detrital zircons in the Carboniferous rocks of both the Mixteca and Oaxaquia terranes contain Devonian detrital zircons and volcanic clasts that are inferred to have come from the removed Devonian arc on the western margin of the Mixteca terrane and/or from exhumed HP rocks. During the Permian, arc-related intrusions in both the Mixteca and Oaxaquia terranes were accompanied by dextral transtensional deformation and deposition of clastic rocks containing Permian detrital zircons and carbonates in periarc, pull-apart basins. Empirical relationships between the dip of the Benioff zone and the widths of arc and forearc indicate that the Permian trench lay beneath the eastern edge of the Mesozoic Guerrero terrane.

Olivine fabric transitions and shear wave anisotropy in the Ryukyu subduction system

Recent shear wave splitting measurements from the fore-arc region of the Ryukyu subduction system show large magnitude (0.3–1.6 s) trench-parallel splitting in both local and teleseismic phases. The similarity of splitting parameters associated with shallow local-S and teleseismic phases suggests that the source of anisotropy is located in the fore-arc mantle. One explanation for this pattern of shear wave splitting involves a transition from commonly observed high-temperature olivine fabrics with flow-parallel seismically fast directions to a flow-normal B-type olivine fabric in the cold fore-arc mantle of the Ryukyu wedge. We test the B-type fabric hypothesis by comparing observed splitting parameters to those predicted from geodynamic models that incorporate olivine fabric development. The distribution of olivine fabric is calculated with high-resolution thermomechanical models of the Ryukyu subduction zone that include realistic slab geometry and an experimentally based wet olivine rheology. We conclude that B-type fabric can explain the magnitude and trench-parallel orientation of deep local-S phases that sample the core of the fore-arc mantle. However, our calculations show that B-type fabric alone cannot account for large magnitude trench-parallel splitting associated with teleseismic phases that sample the shallow tip of the fore-arc mantle. Alternative models for trench-parallel teleseismic splitting in the shallow tip of the fore-arc mantle involve the addition of crustal or slab anisotropy and highly anisotropic foliated antigorite serpentinite.