Imaging megathrust zone and Yakutat/Pacific plate interface in Alaska subduction zone

NASA Astrophysics Data System (ADS)

Kim, Y.; Abers, G. A.; Li, J.; Christensen, D. H.; Calkins, J. A.

2012-12-01

We image the subducted slab underneath a 450 km long transect of the Alaska subduction zone. Dense stations in southern Alaska are set up to investigate (1) the geometry and velocity structure of the downgoing plate and their relation to slab seismicity, and (2) the interplate coupled zone where the great 1964 (magnitude 9.3) had greatest rupture. The joint teleseismic migration of two array datasets (MOOS, Multidisciplinary Observations of Onshore Subduction, and BEAAR, Broadband Experiment Across the Alaska Range) based on teleseismic receiver functions (RFs) using the MOOS data reveal a shallow-dipping prominent low-velocity layer at ~25-30 km depth in southern Alaska. Modeling of these RF amplitudes shows a thin (3-6.5 km) low-velocity layer (shear wave velocity less than 3 km/s), which is ~20-30% slower than normal oceanic crustal velocities, between the subducted slab and the overriding North America plate. The observed low-velocity megathrust layer (with Vp/Vs ratio exceeding 2.0) may be due to a thick sediment input from the trench in combination of elevated pore fluid pressure in the channel. The subducted crust below the low-velocity channel has gabbroic velocities with a thickness of 11-15 km. Both velocities and thickness of the low-velocity channel abruptly increase as the slab bends in central Alaska, which agrees with previously published RF results. Our image also includes an unusually thick low-velocity crust subducting with a ~20 degree dip down to 130 km depth at approximately 200 km inland beneath central Alaska. The unusual nature of this subducted segment has been suggested to be due to the subduction of the Yakutat terrane. Subduction of this buoyant crust could explain the shallow dip of the thrust zone beneath southern Alaska. We also show a clear image of the Yakutat and Pacific plate subduction beneath the Kenai Peninsula, and the along-strike boundary between them at megathrust depths. Our imaged western edge of the Yakutat terrane, at ~30-42 km depth in the central Kenai along the megathrust, aligns with the western end of the geodetically locked patch with high slip deficit, and coincides with the boundary of aftershock events from the 1964 earthquake. It seems plausible that this sharp change in the nature of the downgoing plate controls the slip distribution of great earthquakes on this plate interface.

A thermo-mechanical model of horizontal subduction below an overriding plate

NASA Astrophysics Data System (ADS)

van Hunen, Jeroen; van den Berg, Arie P.; Vlaar, Nico J.

2000-10-01

Subduction of young oceanic lithosphere cannot be explained by the gravitational driving mechanisms of slab pull and ridge push. This deficiency of driving forces can be overcome by obduction of an actively overriding plate, which forces the young plate either to subduct or to collide. This mechanism leads to shallow flattening of the slab as observed today under parts of the west coast of North and South America. Here this process is examined by means of numerical modeling. The convergence velocity between oceanic and continental lithospheric plates is computed from the modeling results, and the ratio of the subduction velocity over the overriding velocity is used as a diagnostic of the efficiency of the ongoing subduction process. We have investigated several factors influencing the mechanical resistance working against the subduction process. In particular, we have studied the effect of a preexisting lithospheric fault with a depth dependent shear resistance, partly decoupling the oceanic lithosphere from the overriding continent. We also investigated the lubricating effect of a 7 km thick basaltic crustal layer on the efficiency of the subduction process and found a log-linear relation between convergence rate and viscosity prefactor characterizing the strength of the oceanic crust, for a range of parameter values including values for basaltic rocks, derived from empirical data. A strong mantle fixes the subducting slab while being overridden and prevents the slab from further subduction in a Benioff style. Viscous heating lowers the coupling strength of the crustal interface between the converging plates with about half an order of magnitude and therefore contributes significantly to the subduction process. Finally, when varying the overriding velocity from 2.5 to 10 cm yr -1, we found a non-linear increase of the subduction velocity due to the presence of non-linear mantle rheology. These results indicate that active obduction of oceanic lithosphere by an overriding continental lithosphere is a viable mechanism for shallow flat subduction over a wide range of model parameters.

Small-scale Forearc Structure from Residual Bathymetry and Vertical Gravity Gradients at the Cocos-North America Subduction Zone offshore Mexico

NASA Astrophysics Data System (ADS)

Garcia, E. S. M.; Ito, Y.

2017-12-01

The subduction of topographic relief on the incoming plate at subduction zones causes deformation of the plate interface as well as the overriding plate. Whether the resulting geometric irregularities play any role in inhibiting or inducing seismic rupture is a topic of relevance for megathrust earthquake source studies. A method to discern the small-scale structure at subduction zone forearcs was recently developed by Bassett and Watts (2015). Their technique constructs an ensemble average of the trench-perpendicular topography, and the removal of this regional tectonic signal reveals the short-wavelength residual bathymetric anomalies. Using examples from selected areas at the Tonga, Mariana, and Japan subduction zones, they were able to link residual bathymetric anomalies to the subduction of seamount chains, given the similarities in wavelength and amplitude to the morphology of seamounts that have yet to subduct. We focus here on an analysis of forearc structures found in the Mexico segment of the Middle America subduction zone, and their potential mechanical interaction with areas on the plate interface that have been previously identified as source regions for earthquake ruptures and aseismic events. We identified several prominent residual bathymetric anomalies off the Guerrero and Oaxaca coastlines, mainly in the shallow portion of the plate interface and between 15 and 50 kilometers away from the trench axis. The residual amplitude of these bathymetric anomalies is typically in the hundreds of meters. Some of the residual bathymetric anomalies offshore Oaxaca are found landward of seamount chains on the incoming Cocos Plate, suggesting that these anomalies are associated with the prior subduction of seamounts at the margin. We also separated the residual and regional components of satellite-based vertical gravity gradient data using a directional median filter to isolate the possible gravity signals from the seamount edifices.

Seismic evidence for deep fluid circulation in the overriding plate of subduction zones

NASA Astrophysics Data System (ADS)

Tauzin, B.; Reynard, B.; Bodin, T.; Perrillat, J. P.; Debayle, E.

2015-12-01

In subduction zones, non-volcanic tremors are associated with fluid circulations (Obara, 2002). Their sources are often located on the interplate boundary (Rogers and Dragert, 2003; Shelly et al, 2006; La Rocca, 2009), consistent with fluids released by the dehydration of subducted plates (Hacker et al., 2003). Reports of tremors in the overriding continental crust of several subduction zones in the world (Kao et al., 2005; Payero et al., 2008; Ide, 2012) suggest fluid circulation at shallower depths but potential fluid paths are poorly documented. Here we obtained seismic observations from receiver functions that evidence the close association between the shallow tremor zone, electrical conductivity, and tectonic features of the Cascadia overriding plate. A seismic discontinuity near 15 km depth in the crust of the overriding North American plate is attributed to the Conrad discontinuity. This interface is segmented, and its interruption is spatially correlated with conductive regions and shallow swarms of seismicity and non-volcanic tremors. These observations suggest that shallow fluid circulation, tremors and seismicity are controlled by fault zones limiting blocks of accreted terranes in the overriding plate (Brudzinski and Allen, 2007). These zones constitute fluid "escape" routes that may contribute unloading fluid pressure on the megathrust. Obara, K. (2002). Science, 296, 1679-1681. Rogers, G., & Dragert, H. (2003). Science, 300, 1942-1943. Shelly, D. R., et al. (2006). Nature, 442, 188-191. La Rocca, M., et al. (2009). Science, 323, 620-623. Kao, H., et al. (2005). Nature, 436, 841-844. Payero, J. S., et al. (2008). Geophysical Research Letters, 35. Ide, S. (2012). Journal of Geophysical Research: Solid Earth, 117. Brudzinski, M. R., & Allen, R. M. (2007). Geology, 35, 907-910.

Serpentinites and Boron Isotope Evidence for Shallow Fluid Transfer Across Subduction Zones

NASA Astrophysics Data System (ADS)

Scambelluri, M.; Tonarini, S.

2012-04-01

In subduction zones, fluid-mediated chemical exchanges between subducting plates and overlying mantle dictate volatile and incompatible element cycles in earth and influence arc magmatism. One of the outstanding issues is concerned with the sources of water for arc magmas and mechanisms for its slab-to-mantle wedge transport. Does it occur by slab dehydration at depths directly beneath arc front, or by hydration of fore-arc mantle and subsequent subduction of the hydrated mantle? Historically, the deep slab dehydration hypothesis had strong support, but it appears that the hydrated mantle wedge hypothesis is gaining ground. At the center of this hypothesis are studies of fluid-mobile element tracers in volatile-rich mantle wedge peridotites (serpentinites) and their subducted high-pressure equivalents. Serpentinites are key players in volatile and fluid-mobile element cycles in subduction zones. Their dehydration represents the main event for fluid and element flux from slabs to mantle, though direct evidence for this process and identification of dehydration environments have been elusive. Boron isotopes are known markers of fluid-assisted element transfer during subduction and can be the tracers of these processes. Until recently, the altered oceanic crust has been considered the main 11B reservoir for arc magmas, which largely display positive delta11B. However, slab dehydration below fore-arcs transfers 11B to the overlying hydrated mantle and leaves the residual mafic crust very depleted in 11B below sub-arcs. The 11B-rich composition of serpentinites candidate them as the heavy B carriers for subduction. Here we present high positive delta11B of Alpine high-pressure (HP) serpentinites recording subduction metamorphism from hydration at low gades to eclogite-facies dehydration: we show a connection among serpentinite dehydration, release of 11B-rich fluids and arc magmatism. In general, the delta11B of these rocks is heavy (16‰ to + 24‰ delta11B). No B loss and no 11B fractionation occurs in these rocks with progressive burial: their high B and 11B compositions demonstrate that initially high budgets acquired during shallow hydration are transferred and released to fluids at arc magma depths, providing the high-boron component requested for arcs. Interaction of depleted mantle-wedge with de-serpentinization fluids and/or serpentinite diapirs uprising from the slab-mantle interface thus provide an efficient self-consistent mechanism for water and B transfer to many arcs. The boron compositions documented here for Erro-Tobbio serpentinites are unexpected for slabs, deputed to loose much B and 11B during subduction dehydration. Their isotopic compositions can be achieved diluting through the mantle the subduction-fluids released during shallow dehydration (30 km) of a model slab. Moreover their delta11B is close to values measured in Syros eclogite blocks, hosted in mélanges atop of the slab and metasomatized by uprising subduction-fluids. The nature of serpentinizing fluids and the fluid-transfer mechanism in Erro-Tobbio is further clarified integrating B isotopes with O-H and Sr isotopic systems. Low deltaD (-102‰), high delta18O (8‰) of early serpentinites suggest low-temperature hydration by metamorphic fluids. 87Sr/86Sr ranges from 0.7044 to 0.7065 and is lower than oceanic serpentinites formed from seawater. Our data indicate that alteration occurred distant from mid-ocean ridges: we propose metamorphic environments like the slab-mantle interface or the fore-arc mantle fed by B- and 11B-rich slab fluids. We therefore provide field-based evidence for delivery of water and 11B at sub-arcs by serpentinites formed by subduction-fluid infiltration in mantle rocks atop of the slab since the early stages of burial, witnessing shallow fluid transfer across the subduction zone.

An image of P- to S-wave velocity ratios in the forearc of the Central Andean subduction zone

NASA Astrophysics Data System (ADS)

Wasja Bloch, Nikolai; Kummerow, Jörn; Wigger, Peter; Shapiro, Serge

2014-05-01

The ratio of seismic P- to S-wave velocities (the Vp/Vs ratio) of a given rock volume is a sensitive proxy for the detection of fluids and melts. In subduction regimes it has often been inferred from seismic tomography and been used, e.g., to detect pathways of ascending melt above the seismogenic zone, where tomographic methods have their highest resolution. We present Vp/Vs ratios that were computed using only seismic arrival time observations following the approach of Lin and Shearer (2007). This approach has its highest sensitivity in the source volume of a set of nearby seismic events and is hence particularly well suited to directly probe the plate interface. We present data from a temporary local network of short period seismometers that was in operation in the forearc of the Central Andean subduction zone at 21° S between 2005 and 2012. From this database we were able to localize 3253 seismic events (Ml ~0.5--4) with high precision, yielding a detailed image of the seismicity distribution in this region. Seismicity is pervasive within the entire crust of the South American continental plate and exhibits three distinct bands in the subducting slab, the lowermost one being located in the lithospheric mantle of the subducting plate. The highest concentration of seismic events is found in the contact zone between the continental and the oceanic lithosphere at depths between 30 and 50 km. We group seismic events into approximately 100 subsets of nearby events that origin from the same geological structure. For about half of these subsets we are able to extract a reliable local Vp/Vs ratio. In the middle continental crust, Vp/Vs ratios show slightly enhanced values (~1.75). In the lower continental crust towards the plate interface they tend to increase from this value updip and decrease downdip. At the plate interface itself, we observe higher Vp/Vs ratios (>1.8) at shallower depths (between 20 and 40 km). Downdip (40--60 km depth) Vp/Vs ratios decrease to rather typical values (~1.75). The same trend is observed in the lowermost band of mantle seismicity in the subducting slab. Below 80 km depth, where mineral transitions toward the eclogite facies are expected to occur, Vp/Vs ratios tend to be low (<1.75). The consistently high Vp/Vs ratios in the shallow part of the subducting slab hint at the presence of fluids in the porespace of the subducting lithosphere there. In the deeper part, downdip variations of Vp/Vs may be attributed to mineral phase transitions due to the changing P-T-conditions along the subduction pathway.

Shallow seismicity patterns in the northwestern section of the Mexico Subduction Zone

NASA Astrophysics Data System (ADS)

Abbott, Elizabeth R.; Brudzinski, Michael R.

2015-11-01

This study characterizes subduction related seismicity with local deployments along the northwestern section of the Mexico Subduction Zone where 4 portions of the plate interface have ruptured in 1973, 1985, 1995, and 2003. It has been proposed that the subducted boundary between the Cocos and Rivera plates occurs beneath this region, as indicated by inland volcanic activity, a gap in tectonic tremor, and the Manzanillo Trough and Colima Graben, which are depressions thought to be associated with the splitting of the two plates after subduction. Data from 50 broadband stations that comprised the MARS seismic array, deployed from January 2006 to June 2007, were processed with the software program Antelope and its generalized source location algorithm, genloc, to detect and locate earthquakes within the network. Slab surface depth contours from the resulting catalog indicate a change in subduction trajectory between the Rivera and Cocos plates. The earthquake locations are spatially anti-correlated with tectonic tremor, supporting the idea that they represent different types of fault slip. Hypocentral patterns also reveal areas of more intense seismic activity (clusters) that appear to be associated with the 2003 and 1973 megathrust rupture regions. Seismicity concentrated inland of the 2003 rupture is consistent with slip on a shallowly dipping trajectory for the Rivera plate interface as opposed to crustal faulting in the overriding North American plate. A prominent cluster of seismicity within the suspected 1973 rupture zone appears to be a commonly active portion of the megathrust as it has been active during three previous deployments. We support these interpretations by determining focal mechanisms and detailed relocations of the largest events within the 1973 and inland 2003 clusters, which indicate primarily thrust mechanisms near the plate interface.

Preliminary deformation model for National Seismic Hazard map of Indonesia

DOE Office of Scientific and Technical Information (OSTI.GOV)

Meilano, Irwan; Gunawan, Endra; Sarsito, Dina

Preliminary deformation model for the Indonesia’s National Seismic Hazard (NSH) map is constructed as the block rotation and strain accumulation function at the elastic half-space. Deformation due to rigid body motion is estimated by rotating six tectonic blocks in Indonesia. The interseismic deformation due to subduction is estimated by assuming coupling on subduction interface while deformation at active fault is calculated by assuming each of the fault‘s segment slips beneath a locking depth or in combination with creeping in a shallower part. This research shows that rigid body motion dominates the deformation pattern with magnitude more than 15 mm/year, except inmore » the narrow area near subduction zones and active faults where significant deformation reach to 25 mm/year.« less

Variations in fluid transport and seismogenic properties in the Lesser Antilles subduction zone: constraints from joint active-source and local earthquake tomography

NASA Astrophysics Data System (ADS)

Paulatto, M.; Laigle, M.; Charvis, P.; Galve, A.

2015-12-01

The degree of coupling and the seismogenic properties of the plate interface at subduction zones are affected by the abundance of slab fluids and subducted sediments. High fluid input can cause high pore-fluid pressures in the subduction channel and decrease coupling leading to aseismic behaviour. Constraining fluid input and transfer is therefore important for understanding plate coupling and large earthquake hazard, particularly in places where geodetic and seismological constraints are scarce. We use P-wave traveltimes from several active source seismic experiments and P- and S-wave traveltimes from shallow and intermediate depth (< 150 km) local earthquakes recorded on a vast amphibious array of OBSs and land stations to recover the Vp and Vp/Vs structure of the central Lesser Antilles subduction zone. Our model extends between Martinique and Antigua from the prism to the arc and from the surface to a depth of 160 km. We find low Vp and high Vp/Vs ratio (> 1.80) on the top of the slab, at depths of up to 100 km. We interpret this high Vp/Vs ratio anomaly as evidence of elevated fluid content either as free fluids or as bound fluids in hydrated minerals (e.g. serpentinite). The strength and depth extent of the anomaly varies strongly from south to north along the subduction zone and correlates with variations in forearc morphology and with sediment input constrained by multi-channel seismic reflection profiles. The anomaly is stronger and extends to greater depth in the south, offshore Martinique, where sediment input is elevated due to the vicinity of the Orinoco delta. The gently dipping forearc slope observed in this region may be the result of weak coupling of the plate interface. A high Vp/Vs ratio is also observed in the forearc likely indicating a fractured and water-saturated overriding plate. On the other hand the anomaly is weaker and shallower offshore Guadeloupe, where sediment input is low due to subduction of the Barracuda ridge. Here a strong plate coupling is likely responsible for uplifting the inner forearc and formation of the Karukera spur. We infer that variations in plate coupling modulated by slab fluid transport and release are a major factor in determining the distribution of seismic slip in the Lesser Antilles subduction zone.

Metasomatism, Fluid Overpressure and Brecciation at the Slab-Mantle Interface: Insights from the Livingstone Fault, New Zealand

NASA Astrophysics Data System (ADS)

Tarling, M.; Smith, S. A. F.; Scott, J.

2017-12-01

Juxtaposition of mantle peridotite and serpentinite against quartzofeldspathic and mafic schists occurs along the shallow slab-mantle interface in some subduction zones. This part of the subduction interface has been invoked as a possible source region of episodic tremor and slow slip, yet geological observations of fault zone structures and chemical reactions pertinent to this region are quite rare. The >1000 km long Livingstone Fault in New Zealand is a superbly exposed fault zone that provides a suitable analogue (both in terms of scale and the rock types involved) for the shallow slab-mantle interface. The fault is characterized by a foliated and highly sheared serpentinite mélange tens to several hundreds of meters wide that separates (partially serpentinised) peridotites from quartzofeldspathic schists. Talc- and tremolite-forming metasomatic reactions occurred along the margins of the mélange and around entrained pods due to mixing of serpentinite with silica- and calcium-rich fluids derived from the adjacent quartzofeldspathic schist. The metasomatic reactions generated significant volumes of water at the melange-schist contact that became trapped between the two relatively impermeable fault zone lithologies. On the schist side of the contact, brittle faulting was promoted by the formation of a laterally-continuous silicified zone up to tens of metres wide. On the melange side, a zone up to tens of metres wide of `crackle-breccias' containing veined stockworks of tremolite indicates periodic increases of pore pressure sufficient to cause hydraulic fracture of serpentinite. The crackle-breccias are multi-generational indicating that this process was episodic. Sr and Nd isotope data and permeability calculations suggest that the episodic brecciation process was critical to the transfer of fluids across the melange. Our observations suggest that fluid-producing metasomatic reactions along the shallow slab-mantle interface may contribute to the tremor signal by triggering brecciation events and promoting brittle failure in serpentinite and schist.

Subduction Zone Dewatering at the Southern End of New Zealand's Hikurangi Margin - Insights from 2D Seismic Tomography

NASA Astrophysics Data System (ADS)

Crutchley, G. J.; Klaeschen, D.

2016-12-01

The southern end of New Zealand's Hikurangi subduction margin is characterised by highly-oblique convergence as it makes a southward transition into a right-lateral transform plate boundary. Long-offset seismic data that cross part of the offshore portion of this transition zone give new insight into the nature of the margin. We have carried out two-dimensional pre-stack depth migrations with an iterative reflection tomography to update the velocity field on two seismic lines in this area. The depth-migrated sections show much-improved imaging of faulting within the wedge, and the seismic velocities themselves give clues about the distribution of gas and/or overpressured regions at the plate boundary and within the overlying wedge. A fascinating observation is a major splay fault that has been (or continues to be) a preferred dewatering pathway through the wedge, evidenced by a thermal anomaly that has left its mark on the overlying gas hydrate layer. Another interesting observation is a thick and laterally extensive low velocity zone beneath the subduction interface, which might have important implications for the long-term mechanical stability of the interface. Our on-going work on these data is focused on amplitude versus offset analysis in an attempt to better understand the nature of the subduction interface and also the shallower gas hydrate system. This study is an example of how distinct disturbances of the gas hydrate system can provide insight into subduction zone fluid flow processes that are important for understanding wedge stability and ultimately earthquake hazard.

Plateau subduction, intraslab seismicity and the Denali Volcanic Gap

NASA Astrophysics Data System (ADS)

Bostock, M. G.; Chuang, L. Y.; Wech, A.; Plourde, A. P.

2017-12-01

Tectonic tremors in Alaska (USA) are associated with subduction of the Yakutat plateau, but their origins are unclear due to lack of depth constraints. We have processed tremor recordings to extract low-frequency earthquakes (LFEs), and generated a set of six LFE waveform templates via iterative network matched filtering and stacking. The timing of impulsive P (compressional) wave and S (shear) wave arrivals on template waveforms places LFEs at 40-58 km depth, near the upper envelope of intraslab seismicity and immediately updip of increased levels of intraslab seismicity. S waves at near-epicentral distances display polarities consistent with shear slip on the plate boundary. We compare characteristics of LFEs, seismicity, and tectonic structures in central Alaska with those in warm subduction zones, and propose a new model for the region's unusual intraslab seismicity and the enigmatic Denali volcanic gap (i.e., an area of no volcanism where expected). We argue that fluids in the Yakutat plate are confined to its upper crust, and that shallow subduction leads to hydromechanical conditions at the slab interface in central Alaska akin to those in warm subduction zones where similar LFEs and tremor occur. These conditions lead to fluid expulsion at shallow depths, explaining strike-parallel alignment of tremor occurrence with the Denali volcanic gap. Moreover, the lack of double seismic zone and restriction of deep intraslab seismicity to a persistent low-velocity zone are simple consequences of anhydrous conditions prevailing in the lower crust and upper mantle of the Yakutat plate.

Plateau subduction, intraslab seismicity, and the Denali (Alaska) volcanic gap

USGS Publications Warehouse

Chuang, Lindsay Yuling; Bostock, Michael; Wech, Aaron; Plourde, Alexandre

2018-01-01

Tectonic tremors in Alaska (USA) are associated with subduction of the Yakutat plateau, but their origins are unclear due to lack of depth constraints. We have processed tremor recordings to extract low-frequency earthquakes (LFEs), and generated a set of six LFE waveform templates via iterative network matched filtering and stacking. The timing of impulsive P (compressional) wave and S (shear) wave arrivals on template waveforms places LFEs at 40–58 km depth, near the upper envelope of intraslab seismicity and immediately updip of increased levels of intraslab seismicity. S waves at near-epicentral distances display polarities consistent with shear slip on the plate boundary. We compare characteristics of LFEs, seismicity, and tectonic structures in central Alaska with those in warm subduction zones, and propose a new model for the region’s unusual intraslab seismicity and the enigmatic Denali volcanic gap (i.e., an area of no volcanism where expected). We argue that fluids in the Yakutat plate are confined to its upper crust, and that shallow subduction leads to hydromechanical conditions at the slab interface in central Alaska akin to those in warm subduction zones where similar LFEs and tremor occur. These conditions lead to fluid expulsion at shallow depths, explaining strike-parallel alignment of tremor occurrence with the Denali volcanic gap. Moreover, the lack of double seismic zone and restriction of deep intraslab seismicity to a persistent low-velocity zone are simple consequences of anhydrous conditions prevailing in the lower crust and upper mantle of the Yakutat plate.

Seismicity and structure of Nazca Plate subduction zone in southern Peru

NASA Astrophysics Data System (ADS)

Lim, H.; Kim, Y.; Clayton, R. W.; Thurber, C. H.

2016-12-01

We define subducting plate geometries in the Nazca subduction zone by (re)locating intra-slab earthquakes in southern Peru (2-18°S) and taking previously published converted phase analysis results, to clarify the slab geometry and inferred relationships to the seismicity. We also provide both P- and S-wave velocities of the subducting Nazca Plate and mantle wedge portions close to the slab using double-difference tomography (Zhang and Thurber, 2003) to understand upper plate volcanism and subduction process. A total of 492 regional earthquakes from August 2008 to February 2013 recorded from the dense seismic array (PeruSE, 2013) are selected for the relocation and tomography. The relocated seismicity shows a smooth contortion in the slab-dip transition zone for 400 km between the shallow (25°)-to-flat dipping interface in the north and 40°-dipping interface in the south. We find a significant slab-dip difference (up to 10°) between our results and previously published slab models along the profile region sampling the normal-dip slab at depth (>100 km). Robust features in both P- and S-wave tomography inversions are dipping low-velocity slabs down to 100 km transitioning to higher-velocities at 100-140 km in both flat slab and dipping slab regions. Differences in the velocities of the mantle wedge between the two regions may indicate different hydration states in the wedge.

Strong interseismic coupling, fault afterslip, and viscoelastic flow before and after the Oct. 9, 1995 Colima-Jalisco earthquake: continuous GPS measurements from Colima, Mexico

USGS Publications Warehouse

Azua, B.M.; DeMets, C.; Masterlark, Timothy

2002-01-01

Continuous GPS measurements from Colima, Mexico during 4/93-6/01, bracketing the Oct. 9, 1995 M = 8.0 Colima-Jalisco earthquake, provide new constraints on Rivera plate subduction mechanics. Modeling of margin-normal strain accumulation before the earthquake suggests the Rivera-North America subduction interface was fully locked. Transient postseismic motion from 10/ 95-6/97 is well fit by a model that includes logarithmically-decaying fault afterslip, elastic strain from shallow fault relocking, and possibly a minor viscoelastic response, but is fit poorly by models that assume a dominant Maxwell viscoelastic response of the lower crust and upper mantle, independent of the assumed viscosities. Landward, margin-normal motion since mid-1997 is parallel to but ??? 75% slower than the pre-seismic velocity. Afterslip alone fails to account for this slowdown. The viscoelastic response predicted by a FEM correctly resolves the remaining velocity difference within the uncertainties. Both processes thus offset elastic strain accumulating from the relocked subduction interface.

Disparate Tectonic Settings of Devastating Earthquakes in Mexico, September 2017

NASA Astrophysics Data System (ADS)

Li, J.; Chen, W. P.; Ning, J.

2017-12-01

Large earthquakes associated with thrust faulting along the plate interface typically pose the highest seismic risk along subduction zones. However, both damaging earthquakes in Mexico of September 2017 are notable exceptions. The Tehuantepec event on the 8th (Mw 8.1) occurred just landward of the trench but is associated with normal faulting, akin to the large (Ms 8) historical event of 1931 that occurred about 200 km to the northwest along this subduction zone. The Puebla earthquake (on the 19th, Mw 7.1) occurred almost 300 km away from the trench where seismic imaging had indicated that the flat-lying slab steepens abruptly and plunges aseismically into the deep mantle. Here we show that both types of tectonic settings are in fact common along a large portion of the Mexican subduction zone, thus identifying source zones of potentially damaging earthquakes away from the plate interface. Additionally, modeling of broadband waveforms made clear that another significant event (Mw 6.1) on the 23rd, is associated with shallow normal faulting in the upper crust, not directly related to the two damaging earthquakes.

Exhumed Blueschists and Eclogites: Hotter Than the Average Model

NASA Astrophysics Data System (ADS)

Penniston-Dorland, S.; Kohn, M. J.; Manning, C. E.

2014-12-01

The maximum-pressure P-T conditions (Pmax-T) of exhumed subduction-related metamorphic rocks are compared to estimates of P-T conditions predicted by computational thermal models of subduction systems. While the range of proposed models encompasses most Pmax-T, most models are 200-400°C too cold. In general, discrepancies are greatest for Pmax < 2 GPa where only a few of the highest-T modeled paths overlap petrologic observations. Comparison among published models suggests several possible explanations for these differences. Variables that affect temperatures within the subduction zone include the timing of subduction initiation relative to metamorphism, age of the subducting oceanic crust, rate of convergence, and the dip angle of the subducting plate. An additional factor is whether subducted material is constrained to move coherently with the incoming plate or whether it convects within the plate interface. Higher temperatures are predicted for relatively young subducting crust, slow convergence rates, and shallow subduction dips. Simulations in which material from the subducted slab decouples from the slab and rises buoyantly into an overlying weak layer (e.g. hydrated mantle) also result in higher temperatures for exhumed oceanic crust. Our compilation and comparison suggest either that most models are missing one or more important controls on heat sources and heat transfer or that exhumed blueschists and eclogites are more buoyant than typical subducted rocks.

Drainage Asperities on Subduction Megathrusts

NASA Astrophysics Data System (ADS)

Sibson, R. H.

2012-12-01

Geophysical observations coupled with force-balance analyses suggest that the seismogenic shear zone interface of subduction megathrusts is generally fluid-overpressured to near-lithostatic values (λv = Pf/σv > 0.9) below the forearc hanging-wall, strongly modulating the profile of frictional shear resistance. Fluid sources include the accretionary prism at shallow levels and, with increasing depth, metamorphic dehydration of material entrained within the subduction shear zone together with progressive metamorphism of oceanic crust in the downgoing slab. Solution transfer in fine-grained material contained within the deeper subduction shear zone (150 < T < 350°C) likely contributes to hydrothermal sealing of fractures. A dramatic difference may therefore exist between low prefailure permeability surrounding the megathrust and high postfailure fracture permeability along the rupture zone and adjacent areas of aftershock activity. Observed postseismic changes in the velocity structure of the fore-arc hanging-wall led Husen and Kissling (2001) to propose massive fluid loss across the subduction interface following the 1995 Antofagasta, Chile, Mw8.0 megathrust rupture. Such trans-megathrust discharges represent a variant of 'fault-valve' action in which the subduction interface itself acts as a seal trapping overpressured fluids derived from metamorphic dehydration beneath. In low-permeability assemblages the maximum sustainable overpressure is limited by the activation or reactivation of brittle faults and fractures under the prevailing stress state. Highest overpressures tend to occur at low differential stress in compressional stress regimes. Loci for fluid discharge are likely determined by stress heterogeneities along the megathrust (e.g. the hangingwall of the rupture at its downdip termination). Discharge sites may be defined by swarm aftershocks defining activated fault-fracture meshes. However, fluid loss across a subduction interface will be enhanced when the stress-state in the forearc hanging-wall switches from compressional reverse-slip faulting before failure to extensional normal-slip faulting postfailure, as occurred during the 2011 Mw9.0 Tohoku megathrust rupture. Mean stress and fault-normal stress then change from being greater than vertical stress prefailure, to less than vertical stress postfailure. Postfailure reductions in overpressure are expected from a combination of poroelastic effects and fluid loss through fault-fracture networks, enhancing vertical permeability. Mineralised fault-fracture meshes in exhumed fore-arc assemblages (e.g. the Alaska-Juneau Au-quartz vein swarm) testify to the episodic discharge of substantial volumes of hydrothermal fluid (< tens of km3). Localized drainage from the subduction interface shear zone increases frictional strength significantly, giving rise to a postfailure strength asperities. Anticipated strength increases from such fluid discharge depends on the magnitude of the drop in overpressure but are potentially large (< hundreds of MPa). Time to the subsequent failure is then governed by reaccumulation of fluid overpressure as well as shear stress along the subduction interface.

Implications of loading/unloading a subduction zone with a heterogeneously coupled interface

NASA Astrophysics Data System (ADS)

Herman, M. W.; Furlong, K. P.; Govers, R. M. A.

2017-12-01

Numerical models of subduction zones with appropriate physical properties may help understand deformation throughout great earthquake cycles, as well as associated observations such as the distribution of smaller magnitude megathrust earthquakes and surface displacements. Of particular interest are displacements near the trench, where tsunamis are generated. The patterns of co-seismic strain release in great megathrust earthquakes depend on the frictional coupling of the plate interface prior to the event. Geodetic observations during the inter-seismic stage suggest that the plates are fully locked at asperities surrounded by zones of apparent partial coupling. We simulate the accumulation (and release) of elastic strain in the subduction system using a finite element model with a relatively simple geometry and material properties. We demonstrate that inter-seismic apparent partial coupling can be dominantly explained by a distribution of completely locked asperities and zero friction elsewhere. In these models, the interface up-dip of the locked zone (< 15 km depth) accumulates large slip deficit even if its coefficient of friction is zero, as might be inferred from the scarcity of megathrust earthquakes shallower than 15 km in global earthquake catalogs. In addition, the upper plate above a low-friction shallow megathrust accumulates large displacements with little internal strain, potentially leading to large co-seismic block displacements (low displacement gradients) of the near-trench seafloor like those observed following the 2011 Mw 9.0 Tohoku earthquake. This is also consistent with anomalously low co-seismic frictional heating of the shallow megathrust indicated by borehole heat flow measurements after the Tohoku event. Our models also yield insights into slip partitioning throughout multiple earthquake cycles. In smaller ruptures, fault slip is inhibited by nearby locked zones; in subsequent multi-segment ruptures, the rest of this slip deficit may be released, producing significantly larger slip than might be expected based on historical earthquake magnitudes. Finally, because low-friction areas around asperities accumulate some slip deficit but may not rupture co-seismically, these regions may be the primary locations of afterslip following the rupture of the locked patch.

Chilean megathrust earthquake recurrence linked to frictional contrast at depth

NASA Astrophysics Data System (ADS)

Moreno, M.; Li, S.; Melnick, D.; Bedford, J. R.; Baez, J. C.; Motagh, M.; Metzger, S.; Vajedian, S.; Sippl, C.; Gutknecht, B. D.; Contreras-Reyes, E.; Deng, Z.; Tassara, A.; Oncken, O.

2018-04-01

Fundamental processes of the seismic cycle in subduction zones, including those controlling the recurrence and size of great earthquakes, are still poorly understood. Here, by studying the 2016 earthquake in southern Chile—the first large event within the rupture zone of the 1960 earthquake (moment magnitude (Mw) = 9.5)—we show that the frictional zonation of the plate interface fault at depth mechanically controls the timing of more frequent, moderate-size deep events (Mw < 8) and less frequent, tsunamigenic great shallow earthquakes (Mw > 8.5). We model the evolution of stress build-up for a seismogenic zone with heterogeneous friction to examine the link between the 2016 and 1960 earthquakes. Our results suggest that the deeper segments of the seismogenic megathrust are weaker and interseismically loaded by a more strongly coupled, shallower asperity. Deeper segments fail earlier ( 60 yr recurrence), producing moderate-size events that precede the failure of the shallower region, which fails in a great earthquake (recurrence >110 yr). We interpret the contrasting frictional strength and lag time between deeper and shallower earthquakes to be controlled by variations in pore fluid pressure. Our integrated analysis strengthens understanding of the mechanics and timing of great megathrust earthquakes, and therefore could aid in the seismic hazard assessment of other subduction zones.

Shoreline-crossing shear-velocity structure of the Juan de Fuca plate and Cascadia subduction zone from surface waves and receiver functions

NASA Astrophysics Data System (ADS)

Janiszewski, Helen; Gaherty, James; Abers, Geoffrey; Gao, Haiying

2017-04-01

The Cascadia subduction zone (CSZ) is the site of the onshore-offshore Cascadia Initiative, which deployed seismometers extending from the Juan de Fuca ridge to the subduction zone and onshore beyond the volcanic arc. This array allows the unique opportunity to seismically image the evolution and along-strike variation of the crust and mantle of the entire CSZ. We compare teleseismic receiver functions, ambient-noise Rayleigh-wave phase velocities in the 10-20 s period band, and earthquake-source Rayleigh-wave phase velocities from 20-100 s, to determine shear-velocity structure in the upper 200 km. Receiver functions from both onshore and shallow-water offshore sites provide constraints on crustal and plate interface structure. Spectral-domain fitting of ambient-noise empirical Green's functions constrains shear velocity of the crust and shallow mantle. An automated multi-channel cross-correlation analysis of teleseismic Rayleigh waves provides deeper lithosphere and asthenosphere constraints. The amphibious nature of the array means it is essential to examine the effect of noise variability on data quality. Ocean bottom seismometers (OBS) are affected by tilt and compliance noise. Removal of this noise from the vertical components of the OBS is essential for the teleseismic Rayleigh waves; this stabilizes the output phase velocity maps particularly along the coastline where observations are predominately from shallow water OBS. Our noise-corrected phase velocity maps reflect major structures and tectonic transitions including the transition from high-velocity oceanic lithosphere to low-velocity continental lithosphere, high velocities associated with the subducting slab, and low velocities beneath the ridge and arc. We interpret the resulting shear-velocity model in the context of temperature and compositional variation in the incoming plate and along the strike of the CSZ.

Shoreline-Crossing Shear-Velocity Structure of the Juan de Fuca Plate and Cascadia Subduction Zone from Surface Waves and Receiver Functions

NASA Astrophysics Data System (ADS)

Janiszewski, H. A.; Gaherty, J. B.; Abers, G. A.; Gao, H.

2016-12-01

The Cascadia subduction zone (CSZ) is the site of the onshore-offshore Cascadia Initiative, which deployed seismometers extending from the Juan de Fuca ridge to the subduction zone and onshore beyond the volcanic arc. This array allows the unique opportunity to seismically image the evolution and along-strike variation of the crust and mantle of the entire CSZ. We compare teleseismic receiver functions, ambient-noise Rayleigh-wave phase velocities in the 10-20 s period band, and earthquake-source Rayleigh-wave phase velocities from 20-100 s, to determine shear-velocity structure in the upper 200 km. Receiver functions from both onshore and shallow-water offshore sites provide constraints on crustal and plate interface structure. Spectral-domain fitting of ambient-noise empirical Green's functions constrains shear velocity of the crust and shallow mantle. An automated multi-channel cross-correlation analysis of teleseismic Rayleigh waves provides deeper lithosphere and asthenosphere constraints. The amphibious nature of the array means it is essential to examine the effect of noise variability on data quality. Ocean bottom seismometers (OBS) are affected by tilt and compliance noise. Removal of this noise from the vertical components of the OBS is essential for the teleseismic Rayleigh waves; this stabilizes the output phase velocity maps particularly along the coastline where observations are predominately from shallow water OBS. Our noise-corrected phase velocity maps reflect major structures and tectonic transitions including the transition from high-velocity oceanic lithosphere to low-velocity continental lithosphere, high velocities associated with the subducting slab, and low velocities beneath the ridge and arc. We interpret the resulting shear-velocity model in the context of temperature and compositional variation in the incoming plate and along the strike of the CSZ.

Subduction-Related Structure in the Mw 9.2, 1964 Megathrust Rupture Area Offshore Kodiak Island, Alaska

NASA Astrophysics Data System (ADS)

Krabbenhoeft, A.; von Huene, R.; Klaeschen, D.; Miller, J. J.

2016-12-01

Some of the largest earthquakes worldwide, including the 1964 9.2 Mw megathrust earthquake, occurred in Alaskan subduction zones. To better understand rupture processes and their mechanisms, we relate seafloor morphology from multibeam and regional bathymetric compilations with sub-seafloor images and seismic P-wave velocity structures. We re-processed legacy multichannel seismic (MCS) data including shot- and intra-shotgather interpolation, multiple removal and Kirchhoff depth migration. These images even reveal the shallow structure of the subducting oceanic crust. Traveltime tomography of a coincident vintage (1994) wide angle dataset reveals the P-wave velocity distribution as well as the deep structure of the subducting plate to the ocean crust Moho. The subducting oceanic crust morphology is rough and partly hidden by a thick sediment cover that reaches 3 km depth at the trench axis. Bathymetry shows two major contrasting upper plate morphologies: the shallow dipping lower slope consists of trench-parallel ridges that form the accreted prism whereas the steep rough middle and upper slopes are composed of competent older rock.Thrust faults are distributed across the entire slope, some of which connect with the subducted plate interface. A subtle change in seafloor gradient from the lower to the middle slope coincides with a thrust fault zone marking the boundary between the margin framework and the frontal prism. It corresponds to the most prominent lateral increase in seismic P-wave velocities, 25 km landward of the trench axis.Major thrusts in several MCS-lines are correlated with bathymetric data, showing their > 100 km lateral extent, which might also be tsunamigenic paths of earthquake rupture from the seismogenic zone to the seafloor.

Investigation of complex slow slip behavior along the Hikurangi subduction zone with earthquake simulator RSQSim

NASA Astrophysics Data System (ADS)

Colella, H.; Ellis, S. M.; Williams, C. A.

2015-12-01

The Hikurangi subduction zone (New Zealand) is one of many subudction zones that exhibit slow slip behavior. Geodetic observations along the Hikurangi subduction zone are unusual in that not only does the subduction zone exhibit periodic slow slip events at "typical" subduction-zone depths of 25-50 km along the southern part of the margin, but also much shallower depths of 8-15 km along the northern part of the margin. Furthermore, there is evidence for interplay between slow slip events at these different depth ranges (between the deep and shallow events) along the central part of the margin, and some of the slow slip behavior is observed along regions of the interface that were previously considered locked, which raises questions about the slip behavior of this region. This study employs the earthquake simulator, RSQSim, to explore variations in the effective normal stress (i.e., stress after the addition of pore fluid pressures) and the frictional instability necessary to generate the complex slow slip events observed along the Hikurangi margin. Preliminary results suggest that to generate slow slip events with similar recurrence intervals to those observed the effective normal stress (MPa) is 3x higher in the south than the north, 6-9MPa versus 2-3MPa, respectively. Results also suggest that, at a minimum, that some overlap along the central margin must exist between the slow slip sections in the north and south to reproduce the types of slip events observed along the Hikurangi subduction zone. To further validate the results from the simulations, Okada solutions for surface displacements will be compared to geodetic solution to more accurately constrain the areas in which slip behavior varies and the cause(s) for the variation(s).

Triggered Slow Slip and Afterslip on the Southern Hikurangi Subduction Zone Following the Kaikōura Earthquake

NASA Astrophysics Data System (ADS)

Wallace, Laura M.; Hreinsdóttir, Sigrún; Ellis, Susan; Hamling, Ian; D'Anastasio, Elisabetta; Denys, Paul

2018-05-01

The 2016 MW7.8 Kaikōura earthquake ruptured a complex sequence of strike-slip and reverse faults in New Zealand's northeastern South Island. In the months following the earthquake, time-dependent inversions of Global Positioning System and interferometric synthetic aperture radar data reveal up to 0.5 m of afterslip on the subduction interface beneath the northern South Island underlying the crustal faults that ruptured in the earthquake. This is clear evidence that the far southern end of the Hikurangi subduction zone accommodates plate motion. The MW7.8 earthquake also triggered widespread slow slip over much of the subduction zone beneath the North Island. The triggered slow slip included immediate triggering of shallow (<15 km), short (2-3 weeks) slow slip events along much of the east coast, and deep (>30 km), long-term (>1 year) slow slip beneath the southern North Island. The southern Hikurangi slow slip was likely triggered by large (0.5-1.0 MPa) static Coulomb stress increases.

Interactions between strike-slip earthquakes and the subduction interface near the Mendocino Triple Junction

NASA Astrophysics Data System (ADS)

Gong, Jianhua; McGuire, Jeffrey J.

2018-01-01

The interactions between the North American, Pacific, and Gorda plates at the Mendocino Triple Junction (MTJ) create one of the most seismically active regions in North America. The earthquakes rupture all three plate boundaries but also include considerable intraplate seismicity reflecting the strong internal deformation of the Gorda plate. Understanding the stress levels that drive these ruptures and estimating the locking state of the subduction interface are especially important topics for regional earthquake hazard assessment. However owing to the lack of offshore seismic and geodetic instruments, the rupture process of only a few large earthquakes near the MTJ have been studied in detail and the locking state of the subduction interface is not well constrained. In this paper, first, we use the second moments inversion method to study the rupture process of the January 28, 2015 Mw 5.7 earthquake on the Mendocino transform fault that was unusually well recorded by both onshore and offshore strong motion instruments. We estimate the rupture dimension to be approximately 6 km by 3 km corresponding to a stress drop of ∼4 MPa for a crack model. Next we investigate the frictional state of the subduction interface by simulating the afterslip that would be expected there as a result of the stress changes from the 2015 earthquake and a 2010 Mw 6.5 intraplate earthquake within the subducted Gorda plate. We simulate afterslip scenarios for a range of depths of the downdip end of the locked zone defined as the transition to velocity strengthening friction and calculate the corresponding surface deformation expected at onshore GPS monuments. We can rule out a very shallow downdip limit owing to the lack of a detectable signal at onshore GPS stations following the 2010 earthquake. Our simulations indicate that the locking depth on the slab surface is at least 14 km, which suggests that the next M8 earthquake rupture will likely reach the coastline and strong shaking should be expected there.

How geometry and structure control the seismic radiation : spectral element simulation of the dynamic rupture of the Mw 9.0 Tohoku earthquake

NASA Astrophysics Data System (ADS)

Festa, G.; Vilotte, J.; Scala, A.

2012-12-01

The M 9.0, 2011 Tohoku earthquake, along the North American-Pacific plate boundary, East of the Honshu Island, yielded a complex broadband rupture extending southwards over 600 km along strike and triggering a large tsunami that ravaged the East coast of North Japan. Strong motion and high-rate continuous GPS data, recorded all along the Japanese archipelago by the national seismic networks K-Net and Kik-net and geodetic network Geonet, together with teleseismic data, indicated a complex frequency dependent rupture. Low frequency signals (f< 0.1 Hz) inverted from seismic, geodetic and tsunami data, evidenced an extremely compact region of large slip (between 30 to 50 meters), extending along-dip over about 100 km, between the hypocenter and the trench, and 150 to 200 km along strike. This slip asperity was likely the cause of the localized tsunami source and of the large amplitude tsunami waves. High-frequency signals (f>0.5 Hz) were instead generated close to the coast in the deeper part of the subduction zone, by at least four smaller size asperities, with possible repeated slip, and were mostly the cause for the ground shaking felt in the Eastern part of Japan. The deep origin of the high-frequency radiation was also confirmed by teleseismic high frequency back projection analysis. Intermediate frequency analysis showed a transition between the shallow and deeper part of the fault, with the rupture almost confined in a small stripe containing the hypocenter before propagating southward along the strike, indicating a predominant in-plane rupture mechanism in the initial stage of the rupture itself. We numerically investigate the role of the geometry of the subduction interface and of the structural properties of the subduction zone on the broadband dynamic rupture and radiation of the Tohoku earthquake. Based upon the almost in-plane behavior of the rupture in its initial stage, 2D non-smooth spectral element dynamic simulations of the earthquake rupture propagation are performed including the non planar and kink geometry of the subduction interface, together with bi-material interfaces taking into account rapid and large variations of the impedance properties along the subduction interfaces and dynamic normal stress coupling. Based on a number of tomographic studies of the NE Japan subduction zone at different space, evidencing a high-velocity "toe" mantle wedge, and wide-angle reflection and refraction studies, supporting a non planar geometry of the subduction interface with at least two strong bending or kink features, we constrain the subduction geometry and the structural properties of the subduction zone model along an off-Miyagi profile. Through several simulations, we investigate possible structural control on the broadband rupture process of the Tohoku earthquake, in terms of the rupture velocity, seismic radiation and slip/stress distribution along the subduction interface. We Explored the influence of initial stress and interface behavior to capture the main features of the rupture and its radiation pattern. Implications for the broad band strong motion observation are discussed, together with implications for the seismic cycle and future earthquake nucleation.

Source time functions of large Mexican subduction earthquakes, morphology of the Benioff Zone, age of the plate, and their tectonic implications

NASA Astrophysics Data System (ADS)

Singh, S. K.; Mortera, F.

1991-12-01

We study source parameters of large, shallow Mexican subduction zone earthquakes (95°W to 106°W) which occurred between 1928 and 1986 by modeling the P waves recorded on Galitzin-Wilip seismograph in DeBilt (DBN), Holland. For post-1962 events the source parameters retrieved from DBN seismograms alone agree well with those obtained from long-period World-Wide Standardized Seismograph Network records, giving us confidence in our results for pre-1962 events. All earthquakes are shallow (H˜10 to 20 km). With few exceptions the sources in Oaxaca (95°W to 99°W) are very simple. To the northwest of 99°W they are simple as well as complex. The ratio of surface wave to body wave seismic moment (Mos/MoP), which is a measure of long- to short-period radiation, is smaller in Oaxaca (˜ 1.5±0.5) than in the regions northwest of 99°W (˜3.1±1.3). These results suggest a change in the plate interface characteristics near 99°W. The sharp change in the rupture mode and the intersection of the O'Gorman Fracture Zone (OFZ) with the trench occur near 99°W. Two strike-slip events offshore, close to OFZ, suggest a segmentation of the subducting plate near 99°W. The age of the plate near the trench in Oaxaca is not well known; it is possible that it does not increase continuously from northwest to southeast in the region but jumps across 99°W. If so, then the older age of the subducted plate southeast of 99°W may be the cause of the distinct rupture mode of the Oaxaca earthquakes. The length of the Benioff zone, which is greatest below Oaxaca ( ≈ 400 km) and decreases toward the northwest, can be explained by the correlation between the length of the subducted slab and the product of the lithosphere age and convergence rate. The relative complexity of sources, the weaker background seismicity, and the lesser number of aftershocks northwest of Oaxaca may be explained by a stronger interface coupling resulting from subduction of younger oceanic slabs (˜5 to 13 m.y. old) in this region. This, however, explains neither larger Mos/MoP values northwest of Oaxaca nor the low stress drop estimates obtained from the analysis of near-field strong-motion data for the Michoacan earthquake of 1985, both of which indicate weaker coupling of the interface. Thus the issue of whether subduction of very young plates (≤ 10 m.y. old) results in strong or weak coupling remains unsolved from the presently available Mexican data.

Segmentation of plate coupling, fate of subduction fluids, and modes of arc magmatism in Cascadia, inferred from magnetotelluric resistivity

USGS Publications Warehouse

Wannamaker, Philip E.; Evans, Rob L.; Bedrosian, Paul A.; Unsworth, Martyn J.; Maris, Virginie; McGary, R. Shane

2014-01-01

Five magnetotelluric (MT) profiles have been acquired across the Cascadia subduction system and transformed using 2-D and 3-D nonlinear inversion to yield electrical resistivity cross sections to depths of ∼200 km. Distinct changes in plate coupling, subduction fluid evolution, and modes of arc magmatism along the length of Cascadia are clearly expressed in the resistivity structure. Relatively high resistivities under the coasts of northern and southern Cascadia correlate with elevated degrees of inferred plate locking, and suggest fluid- and sediment-deficient conditions. In contrast, the north-central Oregon coastal structure is quite conductive from the plate interface to shallow depths offshore, correlating with poor plate locking and the possible presence of subducted sediments. Low-resistivity fluidized zones develop at slab depths of 35–40 km starting ∼100 km west of the arc on all profiles, and are interpreted to represent prograde metamorphic fluid release from the subducting slab. The fluids rise to forearc Moho levels, and sometimes shallower, as the arc is approached. The zones begin close to clusters of low-frequency earthquakes, suggesting fluid controls on the transition to steady sliding. Under the northern and southern Cascadia arc segments, low upper mantle resistivities are consistent with flux melting above the slab plus possible deep convective backarc upwelling toward the arc. In central Cascadia, extensional deformation is interpreted to segregate upper mantle melts leading to underplating and low resistivities at Moho to lower crustal levels below the arc and nearby backarc. The low- to high-temperature mantle wedge transition lies slightly trenchward of the arc.

Tectonic slicing of subducting oceanic crust along plate interfaces: Numerical modeling

NASA Astrophysics Data System (ADS)

Ruh, J. B.; Le Pourhiet, L.; Agard, Ph.; Burov, E.; Gerya, T.

2015-10-01

Multikilometer-sized slivers of high-pressure low-temperature metamorphic oceanic crust and mantle are observed in many mountain belts. These blueschist and eclogite units were detached from the descending plate during subduction. Large-scale thermo-mechanical numerical models based on finite difference marker-in-cell staggered grid technique are implemented to investigate slicing processes that lead to the detachment of oceanic slivers and their exhumation before the onset of the continental collision phase. In particular, we investigate the role of the serpentinized subcrustal slab mantle in the mechanisms of shallow and deep crustal slicing. Results show that spatially homogeneous serpentinization of the sub-Moho slab mantle leads to complete accretion of oceanic crust within the accretionary wedge. Spatially discontinuous serpentinization of the slab mantle in form of unconnected patches can lead to shallow slicing of the oceanic crust below the accretionary wedge and to its deep slicing at mantle depths depending on the patch length, slab angle, convergence velocity and continental geothermal gradient. P-T paths obtained in this study are compared to natural examples of shallow slicing of the Crescent Terrane below Vancouver Island and deeply sliced crust of the Lago Superiore and Saas-Zermatt units in the Western Alps.

Release of mineral-bound water prior to subduction tied to shallow seismogenic slip off Sumatra

NASA Astrophysics Data System (ADS)

Hüpers, Andre; Torres, Marta E.; Owari, Satoko; McNeill, Lisa C.; Dugan, Brandon; Henstock, Timothy J.; Milliken, Kitty L.; Petronotis, Katerina E.; Backman, Jan; Bourlange, Sylvain; Chemale, Farid; Chen, Wenhuang; Colson, Tobias A.; Frederik, Marina C. G.; Guèrin, Gilles; Hamahashi, Mari; House, Brian M.; Jeppson, Tamara N.; Kachovich, Sarah; Kenigsberg, Abby R.; Kuranaga, Mebae; Kutterolf, Steffen; Mitchison, Freya L.; Mukoyoshi, Hideki; Nair, Nisha; Pickering, Kevin T.; Pouderoux, Hugo F. A.; Shan, Yehua; Song, Insun; Vannucchi, Paola; Vrolijk, Peter J.; Yang, Tao; Zhao, Xixi

2017-05-01

Plate-boundary fault rupture during the 2004 Sumatra-Andaman subduction earthquake extended closer to the trench than expected, increasing earthquake and tsunami size. International Ocean Discovery Program Expedition 362 sampled incoming sediments offshore northern Sumatra, revealing recent release of fresh water within the deep sediments. Thermal modeling links this freshening to amorphous silica dehydration driven by rapid burial-induced temperature increases in the past 9 million years. Complete dehydration of silicates is expected before plate subduction, contrasting with prevailing models for subduction seismogenesis calling for fluid production during subduction. Shallow slip offshore Sumatra appears driven by diagenetic strengthening of deeply buried fault-forming sediments, contrasting with weakening proposed for the shallow Tohoku-Oki 2011 rupture, but our results are applicable to other thickly sedimented subduction zones including those with limited earthquake records.

How much does geometry of seismic sources matter in tsunami modeling? A sensitivity analysis for the Calabrian subduction interface

NASA Astrophysics Data System (ADS)

Tonini, R.; Maesano, F. E.; Tiberti, M. M.; Romano, F.; Scala, A.; Lorito, S.; Volpe, M.; Basili, R.

2017-12-01

The geometry of seismogenic sources could be one of the most important factors concurring to control the generation and the propagation of earthquake-generated tsunamis and their effects on the coasts. Since the majority of potentially tsunamigenic earthquakes occur offshore, the corresponding faults are generally poorly constrained and, consequently, their geometry is often oversimplified as a planar fault. The rupture area of mega-thrust earthquakes in subduction zones, where most of the greatest tsunamis have occurred, extends for tens to hundreds of kilometers both down dip and along strike, and generally deviates from the planar geometry. Therefore, the larger the earthquake size is, the weaker the planar fault assumption become. In this work, we present a sensitivity analysis aimed to explore the effects on modeled tsunamis generated by seismic sources with different degrees of geometric complexities. We focused on the Calabrian subduction zone, located in the Mediterranean Sea, which is characterized by the convergence between the African and European plates, with rates of up to 5 mm/yr. This subduction zone has been considered to have generated some past large earthquakes and tsunamis, despite it shows only in-slab significant seismic activity below 40 km depth and no relevant seismicity in the shallower portion of the interface. Our analysis is performed by defining and modeling an exhaustive set of tsunami scenarios located in the Calabrian subduction and using different models of the subduction interface with increasing geometrical complexity, from a planar surface to a highly detailed 3D surface. The latter was obtained from the interpretation of a dense network of seismic reflection profiles coupled with the analysis of the seismicity distribution. The more relevant effects due to the inclusion of 3D complexities in the seismic source geometry are finally highlighted in terms of the resulting tsunami impact.

Volatile transfer and recycling at convergent margins: Mass-balance and insights from high-P/T metamorphic rocks

NASA Astrophysics Data System (ADS)

Bebout, Gray E.

The efficiency with which volatiles are deeply subducted is governed by devolatilization histories and the geometries and mechanisms of fluid transport deep in subduction zones. Metamorphism along the forearc slab-mantle interface may prevent the deep subduction of many volatile components (e.g., H2O, Cs, B, N, perhaps As, Sb, and U) and result in their transport in fluids toward shallower reservoirs. The release, by devolatilization, and transport of such components toward the seafloor or into the forearc mantle wedge, could in part explain the imbalances between the estimated amounts of subducted volatiles and the amounts returned to Earth's surface. The proportion of the initially subducted volatile component that is retained in rocks subducted to depths greater than those beneath magmatic arcs (>100 km) is largely unknown, complicating assessments of deep mantle volatile budgets. Isotopic and trace element data and volatile contents for the Catalina Schist, the Franciscan Complex, and eclogite-facies complexes in the Alps (and elsewhere) provide insight into the nature and magnitude of fluid production and transport deep in subduction zones and into the possible effects of metamorphism on the compositions of subducting rocks. Compatibilities of the compositions of the subduction-related rocks and fluids with the isotopic and trace element compositions of various mantle-derived materials (igneous rocks, xenoliths, serpentinite seamounts) indicate the potential to trace the recycling of rock and fluid reservoirs chemically and isotopically fractionated during subduction-zone metamorphism.

Fault zone structure and seismic reflection characteristics in zones of slow slip and tsunami earthquakes

NASA Astrophysics Data System (ADS)

Bell, Rebecca; Henrys, Stuart; Sutherland, Rupert; Barker, Daniel; Wallace, Laura; Holden, Caroline; Power, William; Wang, Xiaoming; Morgan, Joanna; Warner, Michael; Downes, Gaye

2015-04-01

Over the last couple of decades we have learned that a whole spectrum of different fault slip behaviour takes place on subduction megathrust faults from stick-slip earthquakes to slow slip and stable sliding. Geophysical data, including seismic reflection data, can be used to characterise margins and fault zones that undergo different modes of slip. In this presentation we will focus on the Hikurangi margin, New Zealand, which exhibits marked along-strike changes in seismic behaviour and margin characteristics. Campaign and continuous GPS measurements reveal deep interseismic coupling and deep slow slip events (~30-60 km) at the southern Hikurangi margin. The northern margin, in contrast, experiences aseismic slip and shallow (<10-15 km) slow slip events (SSE) every 18-24 months with equivalent moment magnitudes of Mw 6.5-6.8. Updip of the SSE region two unusual megathrust earthquakes occurred in March and May 1947 with characteristics typical of tsunami earthquakes. The Hikurangi margin is therefore an excellent natural laboratory to study differential fault slip behaviour. Using 2D seismic reflection, magnetic anomaly and geodetic data we observe in the source areas of the 1947 tsunami earthquakes i) low amplitude interface reflectivity, ii) shallower interface relief, iii) bathymetric ridges, iv) magnetic anomaly highs and in the case of the March 1947 earthquake v) stronger geodetic coupling. We suggest that this is due to the subduction of seamounts, similar in dimensions to seamounts observed on the incoming Pacific plate, to depths of <10 km. We propose a source model for the 1947 tsunami earthquakes based on geophysical data and find that extremely low rupture velocities (c. 300 m/s) are required to model the observed large tsunami run-up heights (Bell et al. 2014, EPSL). Our study suggests that subducted topography can cause the nucleation of moderate earthquakes with complex, low velocity rupture scenarios that enhance tsunami waves, and the role of subducted rough topography in seismic hazard should not be under-estimated. 2D seismic reflection data along the northern Hikurangi margin also image thick (c. 2 km) high-amplitude reflectivity zones (HRZ) coinciding broadly with the source areas of shallow SSEs. The HRZ may be the result of high-fluid content within subduction sediments, suggesting fluids may exert an important control on the generation of SSEs by reducing effective stress (Bell et al. 2010, GJI). However, this hypothesis remains untested. In this presentation, using synthetic models, we will discuss planned future applications of an advanced seismic imaging technique called Full-waveform inversion, integrated with drilling, at subduction margins like Hikurangi to recover fault physical properties at high-resolution in 3D to examine the properties of heterogeneous fault zones.

How long-term dynamics of sediment subduction controls short-term dynamics of seismicity

NASA Astrophysics Data System (ADS)

Brizzi, S.; van Zelst, I.; van Dinther, Y.; Funiciello, F.; Corbi, F.

2017-12-01

Most of the world's greatest earthquakes occur along the subduction megathrust. Weak and porous sediments have been suggested to homogenize the plate interface and thereby promote lateral rupture propagation and great earthquakes. However, the importance of sediment thickness, let alone their physical role, is not yet unequivocally established. Based on a multivariate statistical analysis of a global database of 62 subduction segments, we confirm that sediment thickness is one of the key parameters controlling the maximum magnitude a megathrust can generate. Moreover, Monte Carlo simulations highlighted that the occurrence of great earthquakes on sediment-rich subduction segments is very unlikely (p-value≪0.05) related to pure chance. To understand how sediments in the subduction channel regulate earthquake size, this study extends and demystifies multivariate, spatiotemporally limited data through numerical modeling. We use the 2D Seismo-Thermo-Mechanical modeling approach to simulate both the long- and short-term dynamics of subduction and related seismogenesis (van Dinther et al., JGR, 2013). These models solve for the conservation of mass, momentum and energy using a visco-elasto-plastic rheology with rate-dependent friction. Results show that subducted sediments have a strong influence on the long-term evolution of the convergent margin. Increasing the sediment thickness on the incoming plate from 0 to 6 km causes a decrease of slab dip from 23° to 10°. This, in addition to increased radiogenic heating, extends isotherms, thereby widening the seismogenic portion of the megathrust from 80 to 150 km. Consequently, over tens of thousands of years, we observe that the maximum moment magnitude of megathrust earthquakes increases from 8.2 to 9.2 for these shallower and warmer interfaces. In addition, we observe more and larger splay faults, which could enhance vertical seafloor displacements. These results highlight the primary role of subducted sediments in controlling the seismogenic behavior of subduction megathrusts. Furthermore, they illustrate the distinct impact long-term subduction dynamics can have on short-term seismogenesis. This impact likely surpasses the effect a lower static friction would have on seismogenesis, especially for the maximum magnitude of subduction segments.

Seismicity and geodynamics in the central part of the Vanuatu Arc

NASA Astrophysics Data System (ADS)

Baillard, C.; Crawford, W. C.; Ballu, V.; Regnier, M. M.; Pelletier, B.; Garaebiti, E.

2013-12-01

The Vanuatu Arc (VA) in the southwest Pacific ocean (167°E, 13-20°S), is highly seismically active, with more than 35 events of magnitude Mw ≥ 7 since 1973 (USGS catalog). The geodynamics are dominated by the east-dipping subduction of the Australian Plate under the North Fiji Basin microplate. Convergence rates are estimated to be between 130 and 170 mm/yr, except in the central part of the VA where convergence slows to 30-40 mm/yr. This slowing appears to be the result of blockage by the subducting d'Entrecastaux ridge. To quantify the tectonics of this blocked section, we deployed 30 seismometers in 2008-2009 and 8 GPS stations since 2008, in the forearc region of the central VA. The seismometers recorded more than 100 events/day. Detailed analysis of the earthquake catalog reveals: 1) a seismic gap between 40 and 60 km deep under the two largest islands of the VA (Santo and Malekula); 2) subduction plane and intraplate faulting within the down-going plate; and 3) reduced activity beneath Malekula island , perhaps indicating a locked patch on the subduction plane. We infer the geometry of the subduction interface by combining our catalog with unpublished data from the 2000 Santo Mw 6.9 earthquake and aftershocks and the USGS and Global CMT catalogs. The subduction interface appears to be composed of two different panels: a shallow one with a small dip angle and a deeper one with higher dip starting at a depth of ~50 km. We compare finite-element modeling of these panels to the geodetic data to test the connectedness of the two panels and their degree of locking.

The IODP NanTroSEIZE Transect: Accomplishments and Future Plans

NASA Astrophysics Data System (ADS)

Tobin, H. J.; Kinoshita, M.; Araki, E.; Byrne, T. B.; Kimura, G.; McNeill, L. C.; Moore, G. F.; Saffer, D. M.; Underwood, M.; Saito, S.

2009-12-01

The Nankai Trough Seismogenic Zone Experiment (NanTroSEIZE) is a decade-long project to investigate the processes and properties that determine the nature of frictional locking, creep and other fault behavior governing seismogenic rupture and tsunamigenesis on a major plate boundary where great subduction earthquakes occur. The main goal of the science plan is to sample and instrument the key faults in several locations across the transition from those dominated by frictionally stable, aseismic processes vs. those hypothesized to be frictionally locked (seismogenic) faults of the megathrust system. The transect includes primary drill sites from the incoming plate, across the outer accretionary complex of the lower slope, to the Kumano forearc basin and underlying up-dip end of the likely locked plate interface. The scale of this project required a division into multiple stages of operations, spanning a number of years and IODP expeditions. From September 2007 through October 2009, the NanTroSEIZE science team has achieved many of its primary goals during 5 expeditions. Completed drill sites to date include penetrations ranging from ~200 m to ~1600 m below the sea floor that have documented the faults and wall rocks of both the frontal thrust and out-of-sequence splay faults in the accretionary system, the sedimentary section of the subducting plate, and the thick forearc basin sedimentary record and underlying older subduction complex in the hanging wall of the main plate interface. Major results include characterization of: the fault zone geology, strain localization, and physical properties shallower than ~ 1 km, the distribution of ambient (and paleo-) stress orientations across the transect, the absence of evidence for focused fluid channeling along the principal shallow fault systems, and the tectonic history of the subduction system. Extensive downhole measurements and a 2-ship VSP have further documented stress, pressure, rock strength, and elastic properties around the boreholes. The first temporary long-term monitoring instruments are now in place in one sealed borehole, recording pore pressure and temperature. The most ambitious aspect of the NanTroSEIZE project remains for the now-scheduled next stage: drilling to ~ 7000 m below the sea bed across the faults of the main plate boundary, then placing long-term monitoring instruments into both deep and shallow sealed borehole observatories - all to test hypotheses of locking, strain accumulation, and interseismic fault processes.

Subduction and vertical coastal motions in the eastern Mediterranean

NASA Astrophysics Data System (ADS)

Howell, Andy; Jackson, James; Copley, Alex; McKenzie, Dan; Nissen, Ed

2017-10-01

Convergence in the eastern Mediterranean of oceanic Nubia with Anatolia and the Aegean is complex and poorly understood. Large volumes of sediment obscure the shallow structure of the subduction zone, and since much of the convergence is accommodated aseismically, there are limited earthquake data to constrain its kinematics. We present new source models for recent earthquakes, combining these with field observations, published GPS velocities and reflection-seismic data to investigate faulting in three areas: the Florence Rise, SW Turkey and the Pliny and Strabo Trenches. The depths and locations of earthquakes reveal the geometry of the subducting Nubian plate NE of the Florence Rise, a bathymetric high that is probably formed by deformation of sediment at the surface projection of the Anatolia-Nubia subduction interface. In SW Turkey, the presence of a strike-slip shear zone has often been inferred despite an absence of strike-slip earthquakes. We show that the GPS-derived strain-rate field is consistent with extension on the orthogonal systems of normal faults observed in the region and that strike-slip faulting is not required to explain observed GPS velocities. Further SW, the Pliny and Strabo Trenches are also often interpreted as strike-slip shear zones, but almost all nearby earthquakes have either reverse-faulting or normal-faulting focal mechanisms. Oblique convergence across the trenches may be accommodated either by a partitioned system of strike-slip and reverse faults or by oblique slip on the Aegean-Nubia subduction interface. The observed late-Quaternary vertical motions of coastlines close to the subduction zone are influenced by the interplay between: (1) thickening of the material overriding the subduction interface associated with convergence, which promotes coastal uplift; and (2) subsidence due to extension and associated crustal thinning. Long-wavelength gravity data suggest that some of the observed topographic contrasts in the eastern Mediterranean are supported by mantle convection. However, whether the convection is time dependent and whether its pattern moves relative to Nubia are uncertain, and its contribution to present-day rates of vertical coastal motions is therefore hard to constrain. The observed extension of the overriding material in the subduction system is probably partly related to buoyancy forces arising from topographic contrasts between the Aegean, Anatolia and the Mediterranean seafloor, but the reasons for regional variations are less clear.

Paleogeodetic records of seismic and aseismic subduction from central Sumatran microatolls, Indonesia

USGS Publications Warehouse

Natawidjaja, D.H.; Sieh, K.; Ward, S.N.; Cheng, H.; Edwards, R. Lawrence; Galetzka, J.; Suwargadi, B.W.

2004-01-01

We utilize coral microatolls in western Sumatra to document vertical deformation associated with subduction. Microatolls are very sensitive to fluctuations in sea level and thus act as natural tide gauges. They record not only the magnitude of vertical deformation associated with earthquakes (paleoseismic data), but also continuously track the long-term aseismic deformation that occurs during the intervals between earthquakes (paleogeodetic data). This paper focuses on the twentieth century paleogeodetic history of the equatorial region. Our coral paleogeodetic record of the 1935 event reveals a classical example of deformations produced by seismic rupture of a shallow subduction interface. The site closest to the trench rose 90 cm, whereas sites further east sank by as much as 35 cm. Our model reproduces these paleogeodetic data with a 2.3 m slip event on the interface 88 to 125 km from the trench axis. Our coral paleogeodetic data reveal slow submergence during the decades before and after the event in the areas of coseismic emergence. Likewise, interseismic emergence occurred before and after the 1935 event in areas of coseismic submergence. Among the interesting phenomenon we have discovered in the coral record is evidence of a large aseismic slip or "silent even" in 1962, 27 years after the 1935 event. Paleogeodetic deformation rates in the decades before, after, and between the 1935 and 1962 events have varied both temporally and spatially. During the 25 years following the 1935 event, submergence rates were dramatically greater than in prior decades. During the past four decades, however, rates have been lower than in the preceding decades, but are still higher than they were prior to 1935. These paleogeodetic records enable us to model the kinematics of the subduction interface throughout the twentieth century. Copyright 2004 by the American Geophysical Union.

Seismic variability of subduction thrust faults: Insights from laboratory models

NASA Astrophysics Data System (ADS)

Corbi, F.; Funiciello, F.; Faccenna, C.; Ranalli, G.; Heuret, A.

2011-06-01

Laboratory models are realized to investigate the role of interface roughness, driving rate, and pressure on friction dynamics. The setup consists of a gelatin block driven at constant velocity over sand paper. The interface roughness is quantified in terms of amplitude and wavelength of protrusions, jointly expressed by a reference roughness parameter obtained by their product. Frictional behavior shows a systematic dependence on system parameters. Both stick slip and stable sliding occur, depending on driving rate and interface roughness. Stress drop and frequency of slip episodes vary directly and inversely, respectively, with the reference roughness parameter, reflecting the fundamental role for the amplitude of protrusions. An increase in pressure tends to favor stick slip. Static friction is a steeply decreasing function of the reference roughness parameter. The velocity strengthening/weakening parameter in the state- and rate-dependent dynamic friction law becomes negative for specific values of the reference roughness parameter which are intermediate with respect to the explored range. Despite the simplifications of the adopted setup, which does not address the problem of off-fault fracturing, a comparison of the experimental results with the depth distribution of seismic energy release along subduction thrust faults leads to the hypothesis that their behavior is primarily controlled by the depth- and time-dependent distribution of protrusions. A rough subduction fault at shallow depths, unable to produce significant seismicity because of low lithostatic pressure, evolves into a moderately rough, velocity-weakening fault at intermediate depths. The magnitude of events in this range is calibrated by the interplay between surface roughness and subduction rate. At larger depths, the roughness further decreases and stable sliding becomes gradually more predominant. Thus, although interplate seismicity is ultimately controlled by tectonic parameters (velocity of the plates/trench and the thermal regime), the direct control is exercised by the resulting frictional properties of the plate interface.

3D dynamics of crustal deformation driven by oblique subduction: Northern and Central Andes

NASA Astrophysics Data System (ADS)

Schütt, Jorina M.; Whipp, David M., Jr.

2017-04-01

The geometry and relative motion of colliding plates will affect how and where they deform. In oblique subduction systems, factors such as the dip angle of the subducting plate and the convergence obliquity, as well as the presence of weak zones in the overriding plate, all influence how oblique convergence is partitioned onto various fault systems in the overriding plate. The partitioning of strain into margin-normal slip on the plate-bounding fault and horizontal shearing on a strike-slip system parallel to the margin is mainly controlled by the margin-parallel shear forces acting on the plate interface and the strength of the continental crust. While these plate interface forces are influenced by the dip angle of the subducting plate (i.e., the length of plate interface in the frictional domain) and the obliquity angle between the normal to the plate margin and the plate convergence vector, the strength of the continental crust in the upper plate is strongly affected by the presence or absence of weak zones such as regions of arc volcanism, pre-existing fault systems, or boundaries of stronger crustal blocks. In order to investigate which of these factors are most important in controlling how the overriding continental plate deforms, we compare results of lithospheric-scale 3D numerical geodynamic experiments from two regions in the north-central Andes: the Northern Volcanic Zone (NVZ; 5°N - 3°S) and adjacent Peruvian Flat Slab Segment (PFSS; 3°S -14°S). The NVZ is characterized by a 35° subduction dip angle with an obliquity angle of about 40°, extensive volcanism and significant strain partitioning in the continental crust. In contrast, the PFSS is characterized by flat subduction (the slab flattens beneath the continent at around 100 km depth for several hundred kilometers), an obliquity angle of about 20°, no volcanism and minimal strain partitioning. The plate geometry and convergence obliquity for these regions are incorporated in 3D (1600 x 1600 x 160 km) numerical experiments of oceanic subduction beneath a continent, focusing on the conditions under which strain partitioning occurs in the continental plate. In addition to different slab geometries and obliquity angles, we consider the effect of a continental crustal of uniform strength (friction angle Φ=15^°) versus one including a weak zone in the continental crust (Φ=4^°) that runs parallel to the margin. Results of our experiments show that the obliquity angle has the largest effect on initiating strain partitioning, as expected based on strain partitioning theory, but strain partitioning is clearly enhanced by the presence of a continental weakness. Margin-parallel mass transport velocities in the continental sliver are similar to the values observed in the NVZ (about 1 cm/year) in models with a continental weakness and twice as high as those without. In addition, a shallower subduction angle results in formation of a wider continental sliver. Based upon our results, the lack of strain partitioning observed in the PFSS results from both a low convergence obliquity and lack of a weak zone in the continent, even though the shallow subduction should make strain partitioning more favorable.

Seismic imaging of lithospheric discontinuities and continental evolution

NASA Astrophysics Data System (ADS)

Bostock, M. G.

1999-09-01

Discontinuities in physical properties within the continental lithosphere reflect a range of processes that have contributed to craton stabilization and evolution. A survey of recent seismological studies concerning lithospheric discontinuities is made in an attempt to document their essential characteristics. Results from long-period seismology are inconsistent with the presence of continuous, laterally invariant, isotropic boundaries within the upper mantle at the global scale. At regional scales, two well-defined interfaces termed H (˜60 km depth) and L (˜200 km depth) of continental affinity are identified, with the latter boundary generally exhibiting an anisotropic character. Long-range refraction profiles are frequently characterized by subcontinental mantle that exhibits a complex stratification within the top 200 km. The shallow layering of this package can behave as an imperfect waveguide giving rise to the so-called teleseismic Pn phase, while the L-discontinuity may define its lower base as the culmination of a low velocity zone. High-resolution, seismic reflection profiling provides sufficient detail in a number of cases to document the merging of mantle interfaces into lower continental crust below former collisional sutures and magmatic arcs, thus unambiguously identifying some lithospheric discontinuities with thrust faults and subducted oceanic lithosphere. Collectively, these and other seismic observations point to a continental lithosphere whose internal structure is dominated by a laterally variable, subhorizontal layering. This stratigraphy appears to be more pronounced at shallower lithospheric levels, includes dense, anisotropic layers of order 10 km in thickness, and exhibits horizontal correlation lengths comparable to the lateral dimensions of overlying crustal blocks. A model of craton evolution which relies on shallow subduction as a principal agent of craton stabilization is shown to be broadly compatible with these characteristics.

Strong Matrix & Weak Blocks: Evolutionary Inversion of Mélange Rheological Relationships During Subduction and Its Implications for Seismogenesis

NASA Astrophysics Data System (ADS)

Clarke, A. P.; Vannucchi, P.; Ougier-Simonin, A.; Morgan, J. P.

2017-12-01

Subduction zone interface layers are often conceived to be heterogeneous, polyrheological zones analogous to exhumed mélanges. Mélanges typically contain mechanically strong blocks within a weaker matrix. However, our geomechanical study of the Osa Mélange, SW Costa Rica shows that this mélange contains blocks of altered basalt which are now weaker in friction than their surrounding indurated volcanoclastic matrix. Triaxial deformation experiments were conducted on samples of both the altered basalt blocks and the indurated volcanoclastic matrix at confining pressures of 60 and 120 MPa. These revealed that the volcanoclastic matrix has a strength 7.5 times that of the altered basalt at 60 MPa and 4 times at 120 MPa, with the altered basalt experiencing multi-stage failure. The inverted strength relationship between weaker blocks and stronger matrix evolved during subduction and diagenesis of the melange unit by dewatering, compaction and diagenesis of the matrix and cataclastic brecciation and hydrothermal alteration of the basalt blocks. During the evolution of this material, the matrix progressively indurated until its plastic yield stress became greater than the brittle yield stress of the blocks. At this point, the typical rheological relationship found within melanges inverts and melange blocks can fail seismically as the weakest links along the subduction plate interface. The Osa Melange is currently in the forearc of the erosive Middle America Trench and is being incorporated into the subduction zone interface at the updip limit of seismogenesis. The presence of altered basalt blocks acting as weak inclusions within this rock unit weakens the mélange as a whole rock mass. Seismic fractures can nucleate at or within these weak inclusions and the size of the block may limit the size of initial microseismic rock failure. However, when fractures are able to bridge across the matrix between blocks, significantly larger rupture areas may be possible. While this mechanism is a promising candidate for the updip limit of the unusually shallow seismogenic zone beneath Osa, it remains to be seen whether analogous evolutionary strength-inversions control the updip limit of other subduction seismogenic zones.

An Evaluation of Proposed Mechanisms of Slab Flattening in Central Mexico

NASA Astrophysics Data System (ADS)

Skinner, Steven M.; Clayton, Robert W.

2011-08-01

Central Mexico is the site of an enigmatic zone of flat subduction. The general geometry of the subducting slab has been known for some time and is characterized by a horizontal zone bounded on either side by two moderately dipping sections. We systematically evaluate proposed hypotheses for shallow subduction in Mexico based on the spatial and temporal evidence, and we find no simple or obvious explanation for the shallow subduction in Mexico. We are unable to locate an oceanic lithosphere impactor, or the conjugate of an impactor, that is most often called upon to explain shallow subduction zones as in South America, Japan, and Laramide deformation in the US. The only bathymetric feature that is of the right age and in the correct position on the conjugate plate is a set of unnamed seamounts that are too small to have a significant effect on the buoyancy of the slab. The only candidate that we cannot dismiss is a change in the dynamics of subduction through a change in wedge viscosity, possibly caused by water brought in by the slab.

Detailed Image of the Subducting Plate and Upper mantle Seismic Discontinuities in the Mariana Subduction Zone

NASA Astrophysics Data System (ADS)

Tibi, R.; Wiens, D. A.; Shiobara, H.; Sugioka, H.; Yuan, X.

2006-12-01

We use P-to-S converted teleseismic phases recorded at island and ocean bottom stations in Mariana to image the subducting plate and the upper mantle seismic discontinuities in the Mariana subduction zone. The land and seafloor stations which operated from June 2003 to May 2004, were deployed within the framework of the MARGINS Subduction Factory experiment of the Mariana system. The crust in the sudducting plate is observed at about 80--90 km depth beneath the islands of Saipan, Tinian and Rota. For most of the island stations, a low velocity layer is imaged in the forearc at depth between about 20 and 60 km, with decreasing depths toward the arc. The nature of this feature is not yet clear. We found evidence for double seismic discontinuities at the base of the transition zone near the Mariana slab. A shallower discontinuity is imaged at depths of ~650--715 km, and a deeper interface lies at ~740-- 770 km depth. The amplitudes of the seismic signals suggest that the shear velocity contrasts across the two features are comparable. These characteristics support the interpretation that the discontinuities are the results of the phase transformations in olivine (ringwoodite to post-spinel) and garnet (ilminite to perovskite), respectively, for the pyrolite model of mantle composition.

Imaging megathrust zone and Yakutat/Pacific plate interface in Alaska subduction zone

NASA Astrophysics Data System (ADS)

Kim, Y.; Abers, G. A.; Li, J.; Christensen, D. H.; Calkins, J. A.

2013-05-01

We image the subducted slab underneath a 450 km long transect of the Alaska subduction zone. Dense stations in southern Alaska are set up to investigate (1) the geometry and velocity structure of the downgoing plate and their relation to slab seismicity, and (2) the interplate coupled zone where the great 1964 (magnitude 9.3) had greatest rupture. The joint teleseismic migration of two array datasets (MOOS, Multidisciplinary Observations of Onshore Subduction, and BEAAR, Broadband Experiment Across the Alaska Range) based on teleseismic receiver functions (RFs) using the MOOS data reveal a shallow-dipping prominent low-velocity layer at ~25-30 km depth in southern Alaska. Modeling of these RF amplitudes shows a thin (<6.5 km) low-velocity layer (shear wave velocity of ~3 km/s), which is ~20-30% slower than normal oceanic crustal velocities, between the subducted slab and the overriding North American plate. The observed low-velocity megathrust layer (with P-to-S velocity ratio (Vp/Vs) exceeding 2.0) may be due to a thick sediment input from the trench in combination of elevated pore fluid pressure in the channel. The subducted crust below the low-velocity channel has gabbroic velocities with a thickness of 11-12 km. Both velocities and thickness of the low-velocity channel abruptly increase as the slab bends in central Alaska, which agrees with previously published RF results. Our image also includes an unusually thick low-velocity crust subducting with a ~20 degree dip down to 130 km depth at approximately 200 km inland beneath central Alaska. The unusual nature of this subducted segment has been suggested to be due to the subduction of the Yakutat terrane. We also show a clear image of the Yakutat and Pacific plate subduction beneath the Kenai Peninsula, and the along-strike boundary between them at megathrust depths. Our imaged western edge of the Yakutat terrane, at 25-30 km depth in the central Kenai along the megathrust, aligns with the western end of the geodetically locked patch with high slip deficit, and coincides with the boundary of aftershock events from the 1964 earthquake. It seems plausible that this sharp change in the nature of the downgoing plate controls the slip distribution of great earthquakes on this plate interface.

Spatiotemporal Relationship between Shallow Slow Slip and Repeating Earthquakes in the Northern Hikurangi Subduction Margin, New Zealand

NASA Astrophysics Data System (ADS)

Shaddox, H. R.; Schwartz, S. Y.; Todd, E. K.; Sheehan, A.; Yarce, J.; Nakai, J. S.

2017-12-01

Repeating earthquakes (REs), events with similar waveforms and locations, have been used to estimate aseismic slip-rate variations and their presence has been used as a proxy for slow slip. Here, we look for REs in the area of a well-recorded slow slip event (SSE) using ocean-bottom data from the Hikurangi Ocean Bottom Investigation of Tremor and Slow Slip (HOBITSS) experiment. The HOBITSS array of absolute pressure gauges (APG) and ocean-bottom seismometers (OBS) was deployed offshore of Gisborne, New Zealand from May 2014 - June 2015 above shallow (<15 km), episodic (every 18-24 months) SSEs in the northern Hikurangi subduction margin. The HOBITSS network recorded a Mw 6.8 SSE in September/October 2014. To find repeating earthquakes associated with the 2014 SSE we apply matched filtering using the HOBITSS catalog and OBS stations. Events are initially located in Antelope (dbgenloc or dblocsat2), relocated with either NonLinLoc or BayesLoc, and then relocated again (where possible) using GrowClust, a relative relocation algorithm that uses both differential travel times and waveform cross-correlations. Events with final locations within 5 km of the plate interface in the area of the SSE ( 70x100 km area, <2-12 km depth) are used as templates (85 total); few of these plate boundary events occur within areas of large slip. We define RE pairs using a minimum average cross-correlation coefficient of 0.95 at 3+ stations. Only one family of REs (2 events) is found within peak slip ( 20 cm at 8 km depth), days before slip conclusion. One family (2 events) is found in the southern portion of slip ( 15 cm at 7 km depth), two weeks after the SSE. All other REs found (10 families, 30 events) are collocated with a subducting seamount at the northeast edge of slip (<10 cm). One of these families contains two events that occur days before slip termination; the remainder occur days to weeks after the SSE. First-motion focal mechanisms are determined for each family and variability is evident between families; these events may be related to a fracture network surrounding the subducting seamount rather than occurring on the plate interface. Overall, we only find REs toward the end and after the well-recorded 2014 SSE. Thus, in some environments REs may not be reliable indicators of shallow SSEs.

An investigation of deformation and fluid flow at subduction zones using newly developed instrumentation and finite element modeling

NASA Astrophysics Data System (ADS)

Labonte, Alison Louise

Detecting seafloor deformation events in the offshore convergent margin environment is of particular importance considering the significant seismic hazard at subduction zones. Efforts to gain insight into the earthquake cycle have been made at the Cascadia and Costa Rica subduction margins through recent expansions of onshore GPS and seismic networks. While these studies have given scientists the ability to quantify and locate slip events in the seismogenic zone, there is little technology available for adequately measuring offshore aseismic slip. This dissertation introduces an improved flow meter for detecting seismic and aseismic deformation in submarine environments. The value of such hydrologic measurements for quantifying the geodetics at offshore margins is verified through a finite element modeling (FEM) study in which the character of deformation in the shallow subduction zone is determined from previously recorded hydrologic events at the Costa Rica Pacific margin. Accurately sensing aseismic events is one key to determining the stress state in subduction zones as these slow-slip events act to load or unload the seismogenic zone during the interseismic period. One method for detecting seismic and aseismic strain events is to monitor the hydrogeologic response to strain events using fluid flow meters. Previous instrumentation, the Chemical Aqueous Transport (CAT) meter which measures flow rates through the sediment-water interface, can detect transient events at very low flowrates, down to 0.0001 m/yr. The CAT meter performs well in low flow rate environments and can capture gradual changes in flow rate, as might be expected during ultra slow slip events. However, it cannot accurately quantify high flow rates through fractures and conduits, nor does it have the temporal resolution and accuracy required for detecting transient flow events associated with rapid deformation. The Optical Tracer Injection System (OTIS) developed for this purpose is an electronic flow meter that can measure flow rates of 0.1 to >500 m/yr at a temporal resolution of 30 minutes to 0.5 minutes, respectively. Test deployments of the OTIS at cold seeps in the transpressional Monterey Bay demonstrated the OTIS functionality over this range of flow environments. Although no deformation events were detected during these test deployments, the OTIS's temporally accurate measurements at the vigorously flowing Monterey Bay cold seep rendered valuable insight into the plumbing of the seep system. In addition to the capability to detect transient flow events, a primary functional requirement of the OTIS was the ability to communicate and transfer data for long-term real-time monitoring deployments. Real-time data transfer from the OTIS to the desktop was successful during a test deployment of the Nootka Observatory, an acoustically-linked moored-buoy system. A small array of CAT meters was also deployed at the Nootka transform-Cascadia subduction zone triple junction. Four anomalous flow rate events were observed across all four meters during the yearlong deployment. Although the records have low temporal accuracy, a preliminary explanation for the regional changes in flow rate is made through comparison between flow rate records and seismic records. The flow events are thought to be a result of a tectonic deformation event, possibly with an aseismic component. Further constraints are not feasible given the unknown structure of faulting near the triple junction. In a final proof of concept study, I find that use these hydrologic instruments, which capture unique aseismic flow rate patterns, is a valuable method for extracting information about deformation events on the decollement in the offshore subduction zone margin. Transient flow events observed in the frontal prism during a 1999--2000 deployment of CAT meters on the Costa Rica Pacific margin suggest episodic slow-slip deformation events may be occurring in the shallow subduction zone. The FEM study to infer the character of the hypothetical deformation event driving flow transients verify that indeed, a shallow slow-slip event can reproduce the unique flow rate patterns observed. Along (trench) strike variability in the rupture initiation location, and bidirectional propagation, is one way to explain the opposite sign of flow rate transients observed at different along-strike distances. The larger question stimulated by this dissertation project, is: What are the controls on fault mechanics in offshore subduction zone environments? It appears the shallow subduction zone plate interface doesn't behave solely in response to frictional properties of the sediment lining the decollement. Shallow episodic slip at the Costa Rica Pacific margin and further north off Nicaragua, where a slow earthquake broke through the shallow 'stable-sliding' zone and resulted in a tsunami, are potentially conceived through the normally faulted incoming basement topography. Scientists should seek to map out the controls of faulting mechanics, whatever they may be, at all temporal and spatial scales in order to understand these dynamic subduction zone systems. The quest to understanding these controls, in part, requires the characterization of aseismic and seismic strain occurring over time and space. The techniques presented in this dissertation advance scientists' capability for quantifying such strains. With the new instrumentation presented here, long-term real-time observatory networks on the seafloor, and modeling for characterization of deformation events, the pieces of the subduction zone earthquake cycle puzzle may start to come together.

The 2017/09/08 Mw 8.2 Tehuantepec, Mexico Earthquake: A Large but Compact Dip-Slip Faulting Event Severing the Slab

NASA Astrophysics Data System (ADS)

Hjorleifsdottir, V.; Iglesias, A.; Suarez, G.; Santoyo, M. A.; Villafuerte, C. D.; Ji, C.; Franco-Sánchez, S. I.; Singh, S. K.; Cruz-Atienza, V. M.; Ando, R.

2017-12-01

The Mw 8.2 September 8 earthquake occurred in the middle of the "Tehuantepec Gap", a segment of the Mexican subduction zone that has no historical mentions of a large earthquake. It was, however, not the expected subduction megathrust earthquake, but rather an intraplate, normal faulting event, in the subducting oceanic Cocos plate. The earthquake rupture initiated at a depth of 50 km and propagated NW on a near-vertical plane, breaking towards the surface. Most of the slip was concentrated in the distance range 30-100 km from the hypocenter and at depth between 15 and 50 km, with maximum slip of 15m. The earthquake seems to have broken the entire lithosphere, estimated to be 35 km thick. The strike of the fault is about 20 degrees oblique to the trench but aligned with the existing fabric on the incoming oceanic plate, suggesting a structural control by preexisting intraslab fractures and activation by the extensional stress due to the slab bending and pulling. Aftershocks occurred along the fault plane during the first day after the event, with activation of other parallel structures within the subducting plate, towards the east, as well as in upper plate, in the following days. Coulomb stress modeling suggests that the stress on the plate interface above the rupture was significantly increased where shallow thrust aftershoks took place, and reduced updip of the earthquake. There are several other examples of large intraslab normal faulting earthquakes, near the downdip edge (1931 Mw 7.8 and 1999 Mw 7.5, Oaxaca) or directly below (1997 Mw 7.1, Michoacan) the coupled plate interface, along the Mexican subduction zone. The possibility of events of similar magnitude to the 2017 earthquake occurring close to the coastline, all along this part of the subduction zone, cannot be ruled out.

Deformation in the mantle wedge associated with Laramide flat-slab subduction

NASA Astrophysics Data System (ADS)

Behr, W. M.; Smith, D.

2013-12-01

Early Tertiary crustal deformation preserved ~1500 km from the plate boundary in the western U.S. is considered by most to be related to a narrow segment of shallow Farallon-slab subduction, similar to the modern Pampean flat-slab of the central Andes. Evidence that the slab shallowed enough to penetrate several hundred kilometers inboard of the plate boundary includes a) shearing off of lithosphere and underplating of schists derived from the accretionary wedge beneath the volcanic arc; b) a cessation of arc magmatism and eastward sweeping of the magmatic front; and c) mid-Tertiary eruptions as far east as the Four Corners region of serpentinized ultramafic microbreccia (SUM) sourced from very cold, hydrated mantle lithosphere. Included within the SUM diatremes are eclogites interpreted to represent fragments of the slab itself and/or remnants of older rock from the mantle wedge metasomatized and recrystallized to eclogite along the top of the slab. Also included within the SUM diatremes are deformed peridotites that represent pieces of the variably hydrated mantle wedge as well as tectonically eroded and entrained fragments of the plate interface. These include weakly deformed to strongly foliated tectonites, spectacularly sheared mylonites and ultramylonites, and cataclasites, formed at temperatures ranging from 500-650°C. Some of the deformed samples contain hydrous minerals, including antigorite, chlorite, and/or tremolite/pargasite that were formed in-situ prior to or during deformation. We investigate the rheological and seismic properties of the peridotite samples using detailed microstructural and petrological analyses. Initial EBSD data indicate that an antigorite-bearing mylonite exhibits a B-type olivine LPO, whereas an ultramylonite that lacks hydrous minerals exhibits an A-type olivine LPO. This is consistent with experimental data that indicate B-type LPOs form under hydrous conditions; and it suggests that these rocks record a transition from trench-parallel to trench-perpendicular seismic anisotropy, as commonly observed in the mantle wedge above active subduction zones. We also show that the deformation within these sheared peridotites can be used to estimate the magnitude of shear stress along the contact between the Farallon slab and the overlying North American lithosphere. Shear stresses along the plate interface were moderate to high (~40 MPa), allowing a strong degree of interplate coupling, consistent with the stress transfer required to deform the upper plate and produce the basement-cored uplifts characteristic of the Laramide orogeny (e.g. the Rocky Mountains). These results place important natural constraints on flat-slab subduction mechanics. Schematic representation of Laramide flat-slab subduction (modified from Humphreys et al., 2003, Int. Geo. Rev.). The mantle inclusions examined here are sourced from the mantle wedge above the slab and from a serpentinite melange along the slab interface.

Role of H2O in Generating Subduction Zone Earthquakes

NASA Astrophysics Data System (ADS)

Hasegawa, A.

2017-03-01

A dense nationwide seismic network and high seismic activity in Japan have provided a large volume of high-quality data, enabling high-resolution imaging of the seismic structures defining the Japanese subduction zones. Here, the role of H2O in generating earthquakes in subduction zones is discussed based mainly on recent seismic studies in Japan using these high-quality data. Locations of intermediate-depth intraslab earthquakes and seismic velocity and attenuation structures within the subducted slab provide evidence that strongly supports intermediate-depth intraslab earthquakes, although the details leading to the earthquake rupture are still poorly understood. Coseismic rotations of the principal stress axes observed after great megathrust earthquakes demonstrate that the plate interface is very weak, which is probably caused by overpressured fluids. Detailed tomographic imaging of the seismic velocity structure in and around plate boundary zones suggests that interplate coupling is affected by local fluid overpressure. Seismic tomography studies also show the presence of inclined sheet-like seismic low-velocity, high-attenuation zones in the mantle wedge. These may correspond to the upwelling flow portion of subduction-induced secondary convection in the mantle wedge. The upwelling flows reach the arc Moho directly beneath the volcanic areas, suggesting a direct relationship. H2O originally liberated from the subducted slab is transported by this upwelling flow to the arc crust. The H2O that reaches the crust is overpressured above hydrostatic values, weakening the surrounding crustal rocks and decreasing the shear strength of faults, thereby inducing shallow inland earthquakes. These observations suggest that H2O expelled from the subducting slab plays an important role in generating subduction zone earthquakes both within the subduction zone itself and within the magmatic arc occupying its hanging wall.

Seismic‐wave attenuation determined from tectonic tremor in multiple subduction zones

USGS Publications Warehouse

Yabe, Suguru; Baltay, Annemarie S.; Ide, Satoshi; Beroza, Gregory C.

2014-01-01

Tectonic tremor provides a new source of observations that can be used to constrain the seismic attenuation parameter for ground‐motion prediction and hazard mapping. Traditionally, recorded earthquakes of magnitude ∼3–8 are used to develop ground‐motion prediction equations; however, typical earthquake records may be sparse in areas of high hazard. In this study, we constrain the distance decay of seismic waves using measurements of the amplitude decay of tectonic tremor, which is plentiful in some regions. Tectonic tremor occurs in the frequency band of interest for ground‐motion prediction (i.e., ∼2–8  Hz) and is located on the subducting plate interface, at the lower boundary of where future large earthquakes are expected. We empirically fit the distance decay of peak ground velocity from tremor to determine the attenuation parameter in four subduction zones: Nankai, Japan; Cascadia, United States–Canada; Jalisco, Mexico; and southern Chile. With the large amount of data available from tremor, we show that in the upper plate, the lower crust is less attenuating than the upper crust. We apply the same analysis to intraslab events in Nankai and show the possibility that waves traveling from deeper intraslab events experience more attenuation than those from the shallower tremor due to ray paths that pass through the subducting and highly attenuating oceanic crust. This suggests that high pore‐fluid pressure is present in the tremor source region. These differences imply that the attenuation parameter determined from intraslab earthquakes may underestimate ground motion for future large earthquakes on the plate interface.

Effect of Sediments on Rupture Dynamics of Shallow Subduction Zone Earthquakes and Tsunami Generation

NASA Astrophysics Data System (ADS)

Ma, S.

2011-12-01

Low-velocity fault zones have long been recognized for crustal earthquakes by using fault-zone trapped waves and geodetic observations on land. However, the most pronounced low-velocity fault zones are probably in the subduction zones where sediments on the seafloor are being continuously subducted. In this study I focus on shallow subduction zone earthquakes; these earthquakes pose a serious threat to human society in their ability in generating large tsunamis. Numerous observations indicate that these earthquakes have unusually long rupture durations, low rupture velocities, and/or small stress drops near the trench. However, the underlying physics is unclear. I will use dynamic rupture simulations with a finite-element method to investigate the dynamic stress evolution on faults induced by both sediments and free surface, and its relations with rupture velocity and slip. I will also explore the effect of off-fault yielding of sediments on the rupture characteristics and seafloor deformation. As shown in Ma and Beroza (2008), the more compliant hanging wall combined with free surface greatly increases the strength drop and slip near the trench. Sediments in the subduction zone likely have a significant role in the rupture dynamics of shallow subduction zone earthquakes and tsunami generation.

On the role of subducting oceanic plateaus in the development of shallow flat subduction

NASA Astrophysics Data System (ADS)

van Hunen, Jeroen; van den Berg, Arie P.; Vlaar, Nico J.

2002-08-01

Oceanic plateaus, aseismic ridges or seamount chains all have a thickened crust and their subduction has been proposed as a possible mechanism to explain the occurrence of flat subduction and related absence of arc magmatism below Peru, Central Chile and at the Nankai Trough (Japan). Their extra compositional buoyancy could prohibit the slab from sinking into the mantle. With a numerical thermochemical convection model, we simulated the subduction of an oceanic lithosphere that contains an oceanic crustal plateau of 18-km thickness. With a systematic variation, we examined the required physical parameters to obtain shallow flat subduction. Metastability of the basaltic crust in the eclogite stability field is of crucial importance for the slab to remain buoyant throughout the subduction process. In a 44-Ma-old subducting plate, basalt must be able to survive a temperature of 600-700 °C to keep the plate buoyant sufficiently long to cause a flat-slab segment. We found that the maximum yield stress in the slab must be limited to about 600 MPa to allow for the necessary bending to the horizontal. Young slabs show flat subduction for larger parameter ranges than old slabs, since they are less gravitationally unstable and show less resistance against bending. Hydrous weakening of the mantle wedge area and lowermost continent are required to allow for the necessary deformation of a change in subduction style from steep to flat. The maximum flat slab extent is about 300 km, which is sufficient to explain the observed shallow flat subduction near the Nankai Trough (Japan). However, additional mechanisms, such as active overthrusting by an overriding continental plate, need to be invoked to explain the flat-slab segments up to 500 km long below Peru and Central Chile.

Widespread afterslip and triggered slow slip events following the M7.8 Kaikoura earthquake, New Zealand

NASA Astrophysics Data System (ADS)

Wallace, L. M.; Hreinsdottir, S.; Hamling, I. J.; D'Anastasio, E.; Bartlow, N. M.

2017-12-01

Just after midnight on 14 Nov 2016 (NZ Local time), the M7.8 Kaikoura earthquake ruptured a complex sequence of strike-slip and reverse faults over an approximately 150 km length in the northeastern South Island of New Zealand (Hamling et al., 2017, Science). In the months following the earthquake, time-dependent inversions of InSAR observations and continuous and semi-continuous GPS measurements reveal up to 0.5 m of afterslip on the subduction interface beneath the northern South Island underlying the region of large coseismic slip on crustal faults in the M7.8 earthquake. The geodetic data also require significant afterslip on a subset of the crustal faults that ruptured in the earthquake, including the Needles, Jordan Thrust, and Kekerengu faults. Our best-fitting models also suggest significant afterslip on an offshore reverse fault, in a similar position to one inferred by Clark et al. (2017, EPSL) from coseismic coastal uplift data. The M7.8 earthquake also triggered widespread slow slip occurring over much of the Hikurangi subduction zone beneath the North Island. Immediately following the earthquake, continuous GPS sites operated by GeoNet (www.geonet.org.nz) along the North Island's east coast (above the Hikurangi subduction zone) detected several to 30 mm of eastward motion over the two-week period immediately following the M7.8 event. These sites are located 350-650 km from the M7.8 earthquake. Such large eastward motion along the North Island's east coast following the earthquake is consistent with the initiation of a large slow slip event along the shallow, offshore portion of the Hikurangi subduction zone. In addition to shallow slow slip (<15 km depth) triggered offshore the east coast, we also observe deeper slow slip (>30 km depth) triggered in the Kapiti region at the southern Hikurangi margin. The Kapiti SSE was still ongoing as of August 2017, although we expect it to finish before the end of 2017. Given the large distance of the shallow east coast SSE from the M7.8 earthquake, we suggest that the shallow SSE was more likely to be triggered by dynamic stress changes, while the deeper SSEs closer to the Mw 7.8 were more likely triggered by static stress changes.

Anatomy of an ancient subduction interface at 40 km depth: Insights from P-T-t-d data, and geodynamic implications (Dent Blanche, Western Alps)

NASA Astrophysics Data System (ADS)

Angiboust, Samuel; Glodny, Johannes; Oncken, Onno; Chopin, Christian

2014-05-01

An exhumed metamorphic suture zone over 40 km long is exposed in the Dent Blanche Region of the Western Alps belt, along the Swiss-Italian border. In this region, the metasediment-bearing ophiolitic remnants of the Liguro-Piemontese ocean (Tsaté complex) are overthrusted by a continental, km-sized complex (Dent Blanche Tectonic System: DBTS) of Austro-Alpine affinity. The DBTS represents a strongly deformed composite terrane with independent tectonic slices of continental and oceanic origin. In order to better understand the nature and the geodynamic meaning of the shear zone at the base of the DBTS (Dent Blanche Thrust, DBT) we re-evaluated the pressure-temperature-time-deformation (P-T-t-d) history of these two units using modern thermobarometric tools, Rb/Sr deformation ages and field relationships. Our results show that the Tsaté complex is formed by a stack of km-thick calcschists-bearing tectonic slices, having experienced variable maximum burial temperatures of between 360°C and 490°C at depths of ca. 25-40 km, between 41 Ma and 37 Ma. The Arolla gneissic mylonites constituting the base of the DBTS experienced a continuous record of protracted high-pressure (12-14 kbar), top-to-NW D1 deformation at 450-500°C between 43 and 55 Ma. Some of these primary, peak metamorphic fabrics have been sheared (top-to-SE D2) and backfolded during exhumation and collisional overprint (20 km depth, 35-40 Ma) leading to the regional greenschist facies retrogression particularly prominent within Tsaté metasediments. The final juxtaposition of the DBTS with the Tsaté complex occurred between 350 and 500°C during this later, exhumation-related D2 event. Although some exhumation-related deformation partially reworked D1 primary features, we emphasize that the DBT can be viewed as a remnant of the Alpine early Eocene blueschist-facies subduction interface region. The DBT therefore constitutes the deeper equivalent of some shallower portions of the Alpine subduction interface exposed 200 km eastwards in eastern Switzerland (e.g. Bachmann et al., 2009). Our results shed light on deep (25-45 km) subduction zone structures and dynamics and are therefore of major interest for geophysical studies imaging the plate interface region in active subduction zones.

Tectonic controls on earthquake size distribution and seismicity rate: slab buoyancy and slab bending

NASA Astrophysics Data System (ADS)

Nishikawa, T.; Ide, S.

2014-12-01

There are clear variations in maximum earthquake magnitude among Earth's subduction zones. These variations have been studied extensively and attributed to differences in tectonic properties in subduction zones, such as relative plate velocity and subducting plate age [Ruff and Kanamori, 1980]. In addition to maximum earthquake magnitude, the seismicity of medium to large earthquakes also differs among subduction zones, such as the b-value (i.e., the slope of the earthquake size distribution) and the frequency of seismic events. However, the casual relationship between the seismicity of medium to large earthquakes and subduction zone tectonics has been unclear. Here we divide Earth's subduction zones into over 100 study regions following Ide [2013] and estimate b-values and the background seismicity rate—the frequency of seismic events excluding aftershocks—for subduction zones worldwide using the maximum likelihood method [Utsu, 1965; Aki, 1965] and the epidemic type aftershock sequence (ETAS) model [Ogata, 1988]. We demonstrate that the b-value varies as a function of subducting plate age and trench depth, and that the background seismicity rate is related to the degree of slab bending at the trench. Large earthquakes tend to occur relatively frequently (lower b-values) in shallower subduction zones with younger slabs, and more earthquakes occur in subduction zones with deeper trench and steeper dip angle. These results suggest that slab buoyancy, which depends on subducting plate age, controls the earthquake size distribution, and that intra-slab faults due to slab bending, which increase with the steepness of the slab dip angle, have influence on the frequency of seismic events, because they produce heterogeneity in plate coupling and efficiently inject fluid to elevate pore fluid pressure on the plate interface. This study reveals tectonic factors that control earthquake size distribution and seismicity rate, and these relationships between seismicity and tectonic properties may be useful for seismic risk assessment.

Kinematic variables and water transport control the formation and location of arc volcanoes.

PubMed

Grove, T L; Till, C B; Lev, E; Chatterjee, N; Médard, E

2009-06-04

The processes that give rise to arc magmas at convergent plate margins have long been a subject of scientific research and debate. A consensus has developed that the mantle wedge overlying the subducting slab and fluids and/or melts from the subducting slab itself are involved in the melting process. However, the role of kinematic variables such as slab dip and convergence rate in the formation of arc magmas is still unclear. The depth to the top of the subducting slab beneath volcanic arcs, usually approximately 110 +/- 20 km, was previously thought to be constant among arcs. Recent studies revealed that the depth of intermediate-depth earthquakes underneath volcanic arcs, presumably marking the slab-wedge interface, varies systematically between approximately 60 and 173 km and correlates with slab dip and convergence rate. Water-rich magmas (over 4-6 wt% H(2)O) are found in subduction zones with very different subduction parameters, including those with a shallow-dipping slab (north Japan), or steeply dipping slab (Marianas). Here we propose a simple model to address how kinematic parameters of plate subduction relate to the location of mantle melting at subduction zones. We demonstrate that the location of arc volcanoes is controlled by a combination of conditions: melting in the wedge is induced at the overlap of regions in the wedge that are hotter than the melting curve (solidus) of vapour-saturated peridotite and regions where hydrous minerals both in the wedge and in the subducting slab break down. These two limits for melt generation, when combined with the kinematic parameters of slab dip and convergence rate, provide independent constraints on the thermal structure of the wedge and accurately predict the location of mantle wedge melting and the position of arc volcanoes.

Imaging the transition from Aleutian subduction to Yakutat collision in central Alaska, with local earthquakes and active source data

USGS Publications Warehouse

Eberhart-Phillips, D.; Christensen, D.H.; Brocher, T.M.; Hansen, R.; Ruppert, N.A.; Haeussler, Peter J.; Abers, G.A.

2006-01-01

In southern and central Alaska the subduction and active volcanism of the Aleutian subduction zone give way to a broad plate boundary zone with mountain building and strike-slip faulting, where the Yakutat terrane joins the subducting Pacific plate. The interplay of these tectonic elements can be best understood by considering the entire region in three dimensions. We image three-dimensional seismic velocity using abundant local earthquakes, supplemented by active source data. Crustal low-velocity correlates with basins. The Denali fault zone is a dominant feature with a change in crustal thickness across the fault. A relatively high-velocity subducted slab and a low-velocity mantle wedge are observed, and high Vp/Vs beneath the active volcanic systems, which indicates focusing of partial melt. North of Cook Inlet, the subducted Yakutat slab is characterized by a thick low-velocity, high-Vp/Vs, crust. High-velocity material above the Yakutat slab may represent a residual older slab, which inhibits vertical flow of Yakutat subduction fluids. Alternate lateral flow allows Yakutat subduction fluids to contribute to Cook Inlet volcanism and the Wrangell volcanic field. The apparent northeast edge of the subducted Yakutat slab is southwest of the Wrangell volcanics, which have adakitic composition consistent with melting of this Yakutat slab edge. In the mantle, the Yakutat slab is subducting with the Pacific plate, while at shallower depths the Yakutat slab overthrusts the shallow Pacific plate along the Transition fault. This region of crustal doubling within the shallow slab is associated with extremely strong plate coupling and the primary asperity of the Mw 9.2 great 1964 earthquake. Copyright 2006 by the American Geophysical Union.

Latent heat effects of the major mantle phase transitions on low-angle subduction

NASA Astrophysics Data System (ADS)

van Hunen, Jeroen; van den Berg, Arie P.; Vlaar, Nico J.

2001-08-01

Very low to zero shallow dip angles are observed at several moderately young subduction zones with an active trenchward moving overriding plate. We have investigated the effects of latent heat for this situation, where mantle material is pushed through the major mantle phase transitions during shallow low-angle subduction below the overriding plate. The significance of the buoyancy forces, arising from the latent heat effects, on the dynamics of the shallowly subducting slab is examined by numerical modeling. When a 32-Ma-old slab is overridden with 2.5 cm/yr by a continent, flat subduction occurs with a 4-5 cm/yr convergence rate. When latent heat is included in the model, forced downwellings cause a thermal anomaly and consequently thermal and phase buoyancy forces. Under these circumstances, the flat slab segment subducts horizontally about 350 km further and for about 11 Ma longer than in the case without latent heat, before it breaks through the 400-km phase transition. The style of subduction strongly depends on the mantle rheology: increasing the mantle viscosity by one order of magnitude can change the style of subduction from steep to shallow. Similarly, an overriding velocity of less than 1 cm/yr leads to steep subduction, which gradually changes to flat subduction when increasing the overriding velocity. However, these model parameters do not change the aforementioned effect of the latent heat, provided that low-angle subduction occurs. In all models latent heat resulted in a substantial increase of the flat slab length by 300-400 km. Varying the olivine-spinel transition Clapeyron slope γ from 1 to 6 MPa/K reveals a roughly linear relation between γ and the horizontal length of the slab. Based on these results, we conclude that buoyancy forces due to latent heat of phase transitions play an important role in low-angle subduction below an overriding plate.

Seismic subduction of the Nazca Ridge as shown by the 1996-97 Peru earthquakes

USGS Publications Warehouse

Spence, W.; Mendoza, C.; Engdahl, E.R.; Choy, G.L.; Norabuena, E.

1999-01-01

By rupturing more than half of the shallow subduction interface of the Nazca Ridge, the great November 12, 1996 Peruvian earthquake contradicts the hypothesis that oceanic ridges subduct aseismically. The mainshock's rupture has a length of about 200 km and has an average slip of about 1.4 m. Its moment is 1.5 x 1028 dyne-cm and the corresponding M(w) is 8.0. The mainshock registered three major episodes of moment release as shown by a finite fault inversion of teleseismically recorded broadband body waves. About 55% of the mainshock's total moment release occurred south of the Nazca Ridge, and the remaining moment release occurred at the southern half of the subduction interface of the Nazca Ridge. The rupture south of the Nazca Ridge was elongated parallel to the ridge axis and extended from a shallow depth to about 65 km depth. Because the axis of the Nazca Ridge is at a high angle to the plate convergence direction, the subducting Nazca Ridge has a large southwards component of motion, 5 cm/yr parallel to the coast. The 900-1200 m relief of the southwards sweeping Nazca Ridge is interpreted to act as a 'rigid indenter,' causing the greatest coupling south of the ridge's leading edge and leading to the large observed slip. The mainshock and aftershock hypocenters were relocated using a new procedure that simultaneously inverts local and teleseismic data. Most aftershocks were within the outline of the Nazca Ridge. A three-month delayed aftershock cluster' occurred at the northern part of the subducting Nazca Ridge. Aftershocks were notably lacking at the zone of greatest moment release, to the south of the Nazca Ridge. However, a lone foreshock at the southern end of this zone, some 140 km downstrike of the mainshock's epicenter, implies that conditions existed for rupture into that zone. The 1996 earthquake ruptured much of the inferred source zone of the M(w) 7.9-8.2 earthquake of 1942, although the latter was a slightly larger earthquake. The rupture zone of the 1996 earthquake is immediately north of the seismic gap left by the great earthquakes (M(w) ~8.8-9.1) of 1868 and 1877. The M(w) 8.0 Antofagasta earthquake of 1995 occurred at the southern end of this great seismic gap. The M(w) 8.2 deep-focus Bolivian earthquake of 1994 occurred directly downdip of the 1868 portion of that gap. The recent occurrence of three significant earthquakes on the periphery of the great seismic gap of the 1868 and 1877 events, among other factors, may signal an increased seismic potential for that zone.

Revisiting the physical characterisitics of the subduction interplate seismogenic zones

NASA Astrophysics Data System (ADS)

Heuret, Arnauld; Lallemand, Serge; Funiciello, Francesca; Piromallo, Claudia

2010-05-01

Based on the Centennial earthquake catalog, the revised 1964-2007 EHB hypocenters catalog and the 1976-2007 CMT Harvard catalog, we have extracted the hypocenters, nodal planes and seismic moments of worldwide subduction earthquakes for the 1900-2007 period. For the 1976-2007 period, we combine the focal solutions provided by Harvard and the revised hypocenters from Engdahl et al. (1998). Older events are extracted from the Centennial catalogue (Engdahl and Villasenor, 2002) and they are used to estimate the cumulated seismic moment only. The selection criteria for the subduction earthquakes are similar to those used by Mc Caffrey (1994), i.e., we test if the focal mechanisms are consistent with 1/ shallow thrust events (depth > 70 km, positive slips, and at least one nodal plane gets dip < 45°), and, 2/ the plate interface local geometry and orientation (one nodal plane is oriented toward the volcanic arc, the azimuth of this nodal plane ranges between ± 45° with respect to the trench one, its dip ranges between ± 20° with respect to the slab one and the epicentre is located seaward of the volcanic arc). Our study concerns segments of subduction zones that fit with estimated paleoruptures associated with major events (M > 8). We assume that the seismogenic zone coincides with the distribution of 5.5 < M < 7 subduction earthquakes. We provide a map of the interplate seismogenic zones for 80% of the trench systems including dip, length, downdip and updip limits, we revisit the statistical study done by Pacheco et al. (1993) and test some empirical laws obtained for example by Ruff and Kanamori (1980) in light of a more complete, detailed, accurate and uniform description of the subduction interplate seismogenic zone. Since subduction earthquakes result from stress accumulation along the interplate and stress depends on plates kinematics, subduction zone geometry, thermal state and seismic coupling, we aim to isolate some correlations between parameters. The statistical analysis reveals that: 1- vs, the subduction velocity is the first order controlling parameter of seismogenic zone variability, both in term of geometry and seismic behaviour; 2- steep dip, large vertical extent and narrow horizontal extent of the seismogenic zone are associated to fast subductions, and cold slabs, the opposite holding for slow subductions and warm slabs; the seismogenic zone usually ends in the fore-arc mantle rather than at the upper plate Moho depth; 3- seismic rate () variability is coherent with the geometry of the seismogenic zone:  increases with the dip and with the vertical extent of the seismogenic zone, and it fits with vs and with the subducting plate thermal state; 4- mega-events occurrence determines the level of seismic energy released along the subduction interface, whatever  is; 5- to some extent, the potential size of earthquakes fits with vs and with the seismogenic zone geometry, but second order controlling parameters are more difficult to detect; 6- the plate coupling, measured through Upper Plate Strain, is one possible second order parameter: mega-events are preferentially associated to neutral subductions, i.e. moderate compressive stresses along the plate interface; high plate coupling (compressive UPS) is thought to inhibit mega-events genesis by enhancing the locking of the plate interface and preventing the rupture to extend laterally. This research was supported as part of the Eurohorcs/ESF — European Young Investigators Awards Scheme (resp. F.F.), by funds from the National Research Council of Italy and other National Funding Agencies participating in the 3rd Memorandum of Understanding, as well as from the EC Sixth Framework Programme.

Linking Serpentinite Geochemistry with Possible Alteration and Evolution of Supra-Subduction Wedge Mantle

NASA Astrophysics Data System (ADS)

Scambelluri, M.; Cannaò, E.; Agostini, S.; Gilio, M.

2016-12-01

Serpentinites are able to transport and release volatiles and fluid-mobile elements (FME) found in arc magmas. Constraining the trace element compositions of these rocks and of fluids released by de-serpentinization improves our knowledge of mass transfer from subduction zones to volcanic arcs, and of the role of slab and wedge mantle in this global process. Studies of high-pressure ultramafic rocks exhumed from plate interface settings reveal the fluid/rock interactions atop the slab and the processes that can affect the mantle wedge. Alpine eclogite-facies antigorite serpentinite (Voltri Massif) and fully de-serpentinized meta-peridotite (Cima di Gagnone) are enriched in sediment-derived As, Sb, U, Pb before peak dehydration. Their Sr, Pb and B isotopic compositions are reset during prograde (forearc) interaction with slab fluids. The eclogitic garnet and olivine from the Cima di Gagnone metaperidotite trap primary inclusions of the fluid released during breakdown of antigorite and chlorite. The inclusions display FME enrichments (high Cl, S; variable Cs, Rb, Ba, B, Pb, As, Sb) indicating element release from rocks to fluids during dehydration under subarc conditions. Our studies show that serpentinized mantle rocks from subduction zones sequester FME from slab fluids and convey these components and radiogenic isotopes into the mantle wedge upon dehydration. The geochemical processes revealed by such plate-interface rocks can apply to the supra-subduction mantle. Shallow element release from slabs to mantle wedge, downdrag of this altered mantle and its subsequent (subarc) dehydration transfers crust-derived FMEs to the arc magma sources without the need of concomitant subarc dehydration/melting of metasedimentary slab components. The slab signature detected in arc lavas can thus result from geochemical mixing of sediment, oceanic crust and ultramafic reservoirs into altered wedge-mantle rocks, rather than being attributed to multiple fluids.

Large-scale displacement following the 2016 Kaikōura earthquake

NASA Astrophysics Data System (ADS)

Wang, T.; Peng, D.; Barbot, S.; Wei, S.; Shi, X.

2017-12-01

The 2016 Mw 7.9 Kaikōura earthquake occurred near the southern termination of the Hikurangi subduction system, where a transition from subduction to strike-slip motion dominates the pre-seismic strain accumulation. Dense spatial coverage of the GPS measurements and large amount of Interferometric Synthetic Aperture Radar (InSAR) images provide valuable constraints, from the near field to the far field, to study how the slip is distributed among the subduction interface and the overlying fault system before, during and after the earthquake. We extract time-series deformation from the New Zealand continuous GPS network, and SAR images acquired from Japanese ALOS-2 and European Sentinel-1A/B satellites to image the surface deformation related to the 2016 Kaikōura earthquake. Both GPS and InSAR data, which cover the entire New Zealand region, show that the co-seismic and post-seismic deformations are distributed in an extraordinary large area, as far as to the north tip of the North Island. Based on a coseismic slip model derived from seismic and geodetic observations, we calculate the stress perturbation incurred by the earthquake. We explore a range of possibilities of friction laws and rheology via a linear combination of strain rate in finite volumes and slip velocity on ruptured faults. We obtain the slip distribution that can best explain our geodetic measurements using outlier-insensitive hierarchical Bayesian model, to better understand different mechanisms behind the localized shallow after slip and distributed deformation. Our results indicate that complex interactions between the subduction interface and the overlying fault system play an important role in causing such large-scale deformation during and after the earthquake event.

A quantitative analysis of global intermediate and deep seismicity

NASA Astrophysics Data System (ADS)

Ruscic, Marija; Becker, Dirk; Le Pourhiet, Laetitita; Agard, Philippe; Meier, Thomas

2017-04-01

The seismic activity in subduction zones around the world shows a large spatial variabilty with some regions exhibiting strong seismic activity down to depths of almost 700km while in other places seismicity terminates at depths of about 200 or 300 km. Also the decay of the number of seismic events or of the seismic moment with depth is more pronounced in some regions than in others. The same is true for the variability of the ratio of large to small events (the b-value of the Gutenberg-Richter relation) that is varying with depth. These observations are often linked to parameters of the downgoing plate like age or subduction velocity. In this study we investigate a subset of subduction zones utilizing the revised ISC catalogue of intermediate and deep seismicity to determine statistical parameters well suited to describe properties of intermediate deep and deep events. The seismicity is separated into three depth intervals from 50-175km, 175-400km and >400km based on the depth at which the plate contact decouples, the observed nearly exponential decay of the event rate with depth and the supposed depth of phase transition at 410 km depth where also an increase of the event number with depth is observed. For estimation of the b-value and the exponential decay with depth, a restriction of the investigated time interval to the period after 1997 produced significantly better results indicating a globally homogeneous magnitude scale with the magnitude of completeness of about Mw 5. On a global scale the b-value decreases with depth from values of about 1 at 50-175km to values of slightly below 0.8 for events below 400km. Also, there is a slight increase of the b-value with the age of the subducting plate. These changes in the b-value with depth and with age may indicate a varying fragmentation of the slab. With respect to the ratio of the seismic moment between deeper and shallower parts of the subduction zones a dependence on the age is apparent with older slabs exhibiting higher ratios indicating stronger hydration of older slabs and consequently stronger seismic activity at depth in older and thicker slabs. Furthermore, older slabs show the tendency to larger b-values. This indicates stronger fragmentation of older slabs favoring smaller events. Between 50 km and 300 km depth, seismicity in subduction zones decays nearly exponentially with depth. However, the majority of subduction zones show between about 60 km and 100 km lower seismic activity than expected by an exponential decay. This observation correlates well with findings from petrological studies that rocks are rarely scraped off from the downgoing plate at these depths indicating low seismic coupling and low stresses at the plate interface in a depth range below the seismogenic zone and above 100 km depth were dehydration reactions become virulent. Interestingly, the percentage of this deficit becomes larger with plate age for event frequency (reduced number of events), but decreases for moment release (events have larger magnitudes). It is observed that the forearc high is located above the plate interface with reduced seismic coupling. The forearc high is thus an indication of upward directed return flow along the seismically decoupled plate interface. In addition, it is found that the topography of the forearc high is larger above shallow dipping slabs. A correlation of the depth dependent seismic behavior with the subduction or trench velocity is not observed for the investigated subduction zones. Plate age seems to be the dominating factor for properties of intermediate deep and deep seismicity.

How the Slip Distribution Complexities Control the Tsunami Scenarios: a Sensitivity Analysis for the Hellenic and Calabrian Subduction Interfaces.

NASA Astrophysics Data System (ADS)

Scala, A.; Murphy, S.; Herrero, A.; Maesano, F. E.; Lorito, S.; Romano, F.; Tiberti, M. M.; Tonini, R.; Volpe, M.; Basili, R.

2017-12-01

Recent giant tsunamigenic earthquakes (Sumatra 2004, Chile 2010, Tohoku 2011) have confirmed that the complexity of seismic slip distributions may play a fundamental role in the generation and the amplitude of the tsunami waves. In particular, big patches of large slip on the shallower part of the subduction zones, as well as slow rupture propagation within low rigidity areas, can contribute to increase the tsunamigenic potential thus generating devastating coastal inundation. In the Mediterranean Sea, some subduction structures can be identified, such as the Hellenic Arc at the boundary between the African and Aegean plates, and the Calabrian Arc between the European and African plates. We have modelled these areas using discretized high-resolution 3D fault geometries with realistic variability of the strike and dip angles. In particular, the latter geometries have been constrained from the analysis of a dense network of seismic reflection profiles and the seismicity of the areas. To study the influence of different rigidity conditions, we compare the tsunami scenarios deriving from homogeneous slip to those obtained from depth-dependent slip distributions at different magnitudes. These depth-dependent slip distributions are obtained by imposing a variability with depth of both shear modulus and seismic rate, and the conservation of the dislocation over the whole subduction zone. Furthermore, we generate along the Hellenic and Calabrian arc subduction interfaces an ensemble of stochastic slip distributions using a composite source model technique. To mimic either single or multiple asperity source models, the distribution of sub-events whose sum produces the stochastic slip, are distributed based on a PDF, defined as the combination of either one or more Gaussian functions. Tsunami scenarios are then generated from this ensemble in order to address how the position of the main patch of slip can affect the tsunami amplitude along the coast.

Cascadia Seismicity Related to Seamount Subduction as detected by the Cascadia Initiative Amphibious Data

NASA Astrophysics Data System (ADS)

Morton, E.; Bilek, S. L.; Rowe, C. A.

2016-12-01

Unlike other subduction zones, the Cascadia subduction zone (CSZ) is notable for the absence of detected and located small and moderate magnitude interplate earthquakes, despite the presence of recurring episodic tremor and slip (ETS) downdip and evidence of pre-historic great earthquakes. Thermal and geodetic models indicate that the seismogenic zone exists primarily, if not entirely, offshore; therefore the perceived unusual seismic quiescence may be a consequence of seismic source location in relation to land based seismometers. The Cascadia Initiative (CI) amphibious community seismic experiment includes ocean bottom seismometers (OBS) deployed directly above the presumed locked seismogenic zone. We use the CI dataset to search for small magnitude interplate earthquakes previously undetected using the on-land sensors alone. We implement subspace detection to search for small earthquakes. We build our subspace with template events from existing earthquake catalogs that appear to have occurred on the plate interface, windowing waveforms on CI OBS and land seismometers. Although our efforts will target the entire CSZ margin and full 4-year CI deployment, here we focus on a previously identified cluster off the coast of Oregon, related to a subducting seamount. During the first year of CI deployment, this target area yields 293 unique detections with 86 well-located events. Thirty-two of these events occurred within the seamount cluster, and 13 events were located in another cluster to the northwest of the seamount. Events within the seamount cluster are separated into those whose depths place them on the plate interface, and a shallower set ( 5 km depth). These separate event groups track together temporally, and seem to agree with a model of seamount subduction that creates extensive fracturing around the seamount, rather than stress concentrated at the seamount-plate boundary. During CI year 2, this target area yields >1000 additional event detections.

The 2012 August 27 Mw7.3 El Salvador earthquake: expression of weak coupling on the Middle America subduction zone

NASA Astrophysics Data System (ADS)

Geirsson, Halldor; LaFemina, Peter C.; DeMets, Charles; Hernandez, Douglas Antonio; Mattioli, Glen S.; Rogers, Robert; Rodriguez, Manuel; Marroquin, Griselda; Tenorio, Virginia

2015-09-01

Subduction zones exhibit variable degrees of interseismic coupling as resolved by inversions of geodetic data and analyses of seismic energy release. The degree to which a plate boundary fault is coupled can have profound effects on its seismogenic behaviour. Here we use GPS measurements to estimate co- and post-seismic deformation from the 2012 August 27, Mw7.3 megathrust earthquake offshore El Salvador, which was a tsunami earthquake. Inversions of estimated coseismic displacements are in agreement with published seismically derived source models, which indicate shallow (<20 km depth) rupture of the plate interface. Measured post-seismic deformation in the first year following the earthquake exceeds the coseismic deformation. Our analysis indicates that the post-seismic deformation is dominated by afterslip, as opposed to viscous relaxation, and we estimate a post-seismic moment release one to eight times greater than the coseismic moment during the first 500 d, depending on the relative location of coseismic versus post-seismic slip on the plate interface. We suggest that the excessive post-seismic motion is characteristic for the El Salvador-Nicaragua segment of the Central American margin and may be a characteristic of margins hosting tsunami earthquakes.

Zoned mantle convection.

PubMed

Albarède, Francis; Van Der Hilst, Rob D

2002-11-15

We review the present state of our understanding of mantle convection with respect to geochemical and geophysical evidence and we suggest a model for mantle convection and its evolution over the Earth's history that can reconcile this evidence. Whole-mantle convection, even with material segregated within the D" region just above the core-mantle boundary, is incompatible with the budget of argon and helium and with the inventory of heat sources required by the thermal evolution of the Earth. We show that the deep-mantle composition in lithophilic incompatible elements is inconsistent with the storage of old plates of ordinary oceanic lithosphere, i.e. with the concept of a plate graveyard. Isotopic inventories indicate that the deep-mantle composition is not correctly accounted for by continental debris, primitive material or subducted slabs containing normal oceanic crust. Seismological observations have begun to hint at compositional heterogeneity in the bottom 1000 km or so of the mantle, but there is no compelling evidence in support of an interface between deep and shallow mantle at mid-depth. We suggest that in a system of thermochemical convection, lithospheric plates subduct to a depth that depends - in a complicated fashion - on their composition and thermal structure. The thermal structure of the sinking plates is primarily determined by the direction and rate of convergence, the age of the lithosphere at the trench, the sinking rate and the variation of these parameters over time (i.e. plate-tectonic history) and is not the same for all subduction systems. The sinking rate in the mantle is determined by a combination of thermal (negative) and compositional buoyancy and as regards the latter we consider in particular the effect of the loading of plates with basaltic plateaux produced by plume heads. Barren oceanic plates are relatively buoyant and may be recycled preferentially in the shallow mantle. Oceanic plateau-laden plates have a more pronounced negative buoyancy and can more easily founder to the very base of the mantle. Plateau segregation remains statistical and no sharp compositional interface is expected from the multiple fate of the plates. We show that the variable depth subduction of heavily laden plates can prevent full vertical mixing and preserve a vertical concentration gradient in the mantle. In addition, it can account for the preservation of scattered remnants of primitive material in the deep mantle and therefore for the Ar and (3)He observations in ocean-island basalts.

Resolution Study of Marine CSEM Imaging of Subduction Zones

NASA Astrophysics Data System (ADS)

Gustafson, C.; Key, K.

2016-12-01

Marine controlled source electromagnetic (CSEM) data allow us to image seafloor electrical resistivity from which we can constrain the porosity and fluid content of the subsurface. In subduction zones, CSEM data can be used to constrain geologic structure, hydrogeology and fluid-tectonic processes. The scales of features we are interested in recovering with CSEM data range from large-scale features such as the incoming tectonic plate and subducting slab, to the narrow dipping plate boundary interface where slip occurs, to thin faults that cut the overriding forearc crust and shallow fluid seeps and mounds on the seafloor. Thus electrical structure is expected to vary on scales ranging from scales of meters to tens of kilometers. CSEM data collected by Scripps at the Middle America Trench in 2010 is the first and to-date the only application of the method for studying a subduction zone. The results from this pioneering data set highlight the types of new discoveries that are possible with CSEM data, such as imaging conductive bending faults and a water-rich channel of subducting sediments. In this work we explore the magnitude and scale of 2D resistivity structures that can be resolved with CSEM data through a suite of synthetic inversion studies. We build resistivity models that are representative of various known and hypothesized subduction zone plate boundary structures. We generate synthetic noisy data for these models and invert them using the freely available MARE2DEM inversion code. We compare the recovered models to the original models in order to determine which resistivity structures may be successfully identified using CSEM. We explore the potential effects of receiver spacing, frequency bandwidth and system noise levels on the ability of CSEM to recover these different subduction zone structures.

Slab1.0: A three-dimensional model of global subduction zone geometries

NASA Astrophysics Data System (ADS)

Hayes, Gavin P.; Wald, David J.; Johnson, Rebecca L.

2012-01-01

We describe and present a new model of global subduction zone geometries, called Slab1.0. An extension of previous efforts to constrain the two-dimensional non-planar geometry of subduction zones around the focus of large earthquakes, Slab1.0 describes the detailed, non-planar, three-dimensional geometry of approximately 85% of subduction zones worldwide. While the model focuses on the detailed form of each slab from their trenches through the seismogenic zone, where it combines data sets from active source and passive seismology, it also continues to the limits of their seismic extent in the upper-mid mantle, providing a uniform approach to the definition of the entire seismically active slab geometry. Examples are shown for two well-constrained global locations; models for many other regions are available and can be freely downloaded in several formats from our new Slab1.0 website, http://on.doi.gov/d9ARbS. We describe improvements in our two-dimensional geometry constraint inversion, including the use of `average' active source seismic data profiles in the shallow trench regions where data are otherwise lacking, derived from the interpolation between other active source seismic data along-strike in the same subduction zone. We include several analyses of the uncertainty and robustness of our three-dimensional interpolation methods. In addition, we use the filtered, subduction-related earthquake data sets compiled to build Slab1.0 in a reassessment of previous analyses of the deep limit of the thrust interface seismogenic zone for all subduction zones included in our global model thus far, concluding that the width of these seismogenic zones is on average 30% larger than previous studies have suggested.

How to Simulate the Interplate Domain in Thermo-mechanical Experiments of Subduction ? Critical Effects of Resolution and Rheology, and Consequences on Wet Mantle Melting

NASA Astrophysics Data System (ADS)

Arcay, D.

2017-12-01

Oceanic plate subduction implies tight interactions between converging lithospheres and surrounding mantle. Plate-mantle couplings can be modeled using thermo-chemical codes of mantle convection. But how to model correctly with a continuous fluid approach the subduction interface, characterised by strong and localised discontinuities? The present study aims at better deciphering the different mechanisms involved in the functioning of the subduction interplate, simply modeled by a weak crust layer, free to evolve. Pseudo-brittle and non-Newtonian behaviours are modelled. This study shows first that the numerical resolution is critical. If the subducting plate is 100 Myr old, subduction occurs for any crust strength. The stiffer the crust is, the shallower the interplate down-dip extent is and the hotter the fore-arc base is. Conversely, imposing a very weak subduction channel leads to an extreme mantle wedge cooling and inhibits mantle melting in wet conditions. If the incoming plate is 20 Myr old, subduction occurs only if the crust is either stiff and denser than the mantle, or weak and buoyant. These conditions lead notably to (1) fore-arc lithosphere cooling, and (2) partial or complete hindrance of wet mantle melting. Finally, subduction plane dynamics is intimately linked to the regime of subduction-induced corner flow: either focussed towards the mantle wedge tip and strongly warming the subduction plate, or, diffuse and favoring global cooling by the lengthening of the subduction plane. The thermal states simulated within the mantle wedge are compared with observations to decipher the best rheological ranges modelling the subduction channel. Two intervals of crustal activation energy are underlined: 345-385 kJ/mol to reproduce the slab surface temperature range inferred from geothermometry, and 415-455 kJ/mol to reproduce the hot mantle wedge core suggested by seismic tomographies. As these two intervals do not overlap, an extra process involved in subduction dynamics is needed. A moderate mantle viscosity reduction, caused by metasomatism in the mantle wedge, is proposed. From these results, it can be inferred that the subduction channel down-dip extent should vary with the subduction setting, consistently with the worldwide variability of sub-arc depths of the subducting plate surface.

Subduction Initiation under Unfavorable Conditions and New Fault Formation

NASA Astrophysics Data System (ADS)

Mao, X.; Gurnis, M.; May, D.

2017-12-01

How subduction initiates with unfavorable dipping lithospheric heterogeneities is an important and rarely studied topic. We build a geodynamic model starting with a vertical weak zone for the Puysegur incipient subduction zone (PISZ). A true free surface is tracked in pTatin3D, based on the Arbitrary Lagrangian Eulerian (ALE) finite element method, and is used to follow the dynamic mantle-surface interaction and topographic evolution. A simplified surface process, based on linear topography diffusion, is implemented. Density and free water content for different phase assemblages are gained by referring to precalculated 4D (temperature, pressure, rock type and total water content) phase maps using Perplex. Darcy's law is used to migrate free water, and a linear water weakening is applied to the mantle material. A new visco-elastic formulation called Elastic Viscous Stress Splitting (EVSS) method is also included. Our predictions fit the morphology of the Puysegur Trench and Ridge and the deformation history on the overriding plate. We show a new thrust fault forms and evolves into a smooth subduction interface, and the preexisting weak zone becomes a vertical fault inboard of the thrust fault during subduction initiation, which explains the two-fault system at PISZ. Our model suggests that the PISZ may not yet be self-sustaining. We propose that the Snares Trough is caused by plate coupling differences between shallower and deeper parts, the tectonic sliver between two faults experiences strong rotation, and low density materials accumulate beneath the Snares trough. Extended models show that with favorable dipping heterogeneities, no new fault forms, and subduction initiates with smaller resisting forces.

Rupture Process During the Mw 8.1 2017 Chiapas Mexico Earthquake: Shallow Intraplate Normal Faulting by Slab Bending

NASA Astrophysics Data System (ADS)

Okuwaki, R.; Yagi, Y.

2017-12-01

A seismic source model for the Mw 8.1 2017 Chiapas, Mexico, earthquake was constructed by kinematic waveform inversion using globally observed teleseismic waveforms, suggesting that the earthquake was a normal-faulting event on a steeply dipping plane, with the major slip concentrated around a relatively shallow depth of 28 km. The modeled rupture evolution showed unilateral, downdip propagation northwestward from the hypocenter, and the downdip width of the main rupture was restricted to less than 30 km below the slab interface, suggesting that the downdip extensional stresses due to the slab bending were the primary cause of the earthquake. The rupture front abruptly decelerated at the northwestern end of the main rupture where it intersected the subducting Tehuantepec Fracture Zone, suggesting that the fracture zone may have inhibited further rupture propagation.

Anomalous Seismic Radiation in the Shallow Subduction Zone Explained by Extensive Poroplastic Deformation in the Overriding Wedge

NASA Astrophysics Data System (ADS)

Hirakawa, E. T.; Ma, S.

2012-12-01

The deficiency of high-frequency seismic radiation from shallow subduction zone earthquakes was first recognized in tsunami earthquakes (Kanamori, 1972), which produce larger tsunamis than expected from short-period (20 s) surface wave excitation. Shallow subduction zone earthquakes were also observed to have unusually low energy-to-moment ratios compared to regular subduction zone earthquakes (e.g., Newman and Okal, 1998; Venkataraman and Kanamori, 2004; Lay et al., 2012). What causes this anomalous radiation and how it relates to large tsunami generation has remained unclear. Here we show that these anomalous observations can be due to extensive poroplastic deformation in the overriding wedge, which provides a unifying interpretation. Ma (2012) showed that the pore pressure increase in the wedge due to up-dip rupture propagation significantly weakens the wedge, leading to widespread Coulomb failure in the wedge. Widespread failure gives rise to slow rupture velocity and large seafloor uplift (landward from the trench) in the case of a shallow fault dip. Here we extend this work and demonstrate that the large seafloor uplift due to the poroplastic deformation significantly dilates the fault behind the rupture front, which reduces the normal stress on the fault and increases the stress drop, slip, and rupture duration. The spectral amplitudes of the moment-rate time function is significantly less at high frequencies than those from elastic simulations. Large tsunami generation and deficiency of high-frequency radiation are thus two consistent manifestations of the same mechanism (poroplastic deformation). Although extensive poroplastic deformation in the wedge represents a significant portion of total seismic moment release, the plastic deformation is shown to act as a large energy sink, leaving less energy to be radiated and leading to low energy-to-moment ratios as observed for shallow subduction zone earthquakes.

Are diamond-bearing Cretaceous kimberlites related to shallow-angle subduction beneath western North America?

NASA Astrophysics Data System (ADS)

Currie, C. A.; Beaumont, C.

2009-05-01

The origin of deep-seated magmatism (in particular, kimberlites and lamproites) within continental plate interiors remains enigmatic in the context of plate tectonic theory. One hypothesis proposes a relationship between kimberlite occurrence and lithospheric subduction, such that a subducting plate releases fluids below a continental craton, triggering melting of the deep lithosphere and magmatism (Sharp, 1974; McCandless, 1999). This study provides a quantitative evaluation of this hypothesis, focusing on the Late Cretaceous- Eocene (105-50 Ma) kimberlites and lamproites of western North America. These magmas were emplaced along a corridor of Archean and Proterozoic lithosphere, 1000-1500 km inboard of the plate margin separating the subducting Farallon Plate and continental North America Plate. Kimberlite-lamproite magmatism coincides with tectonic events, including the Laramide orogeny, shut-down of the Sierra Nevada arc, and eastward migration of volcanism, that are commonly attributed to a change in Farallon Plate geometry to a shallow-angle trajectory (<25° dip). Thermal-mechanical numerical models demonstrate that rapid Cretaceous plate convergence rates and enhanced westward velocity of North America result in shallow-angle subduction that places the Farallon Plate beneath the western edge of the cratonic interior of North America. This geometry is consistent with the observed continental dynamic subsidence that lead to the development of the Western Interior Seaway. The models also show that the subducting plate has a cool thermal structure, and subducted hydrous minerals (serpentine, phengite and phlogopite) remain stable to more than 1200 km from the trench, where they may break down and release fluids that infiltrate the overlying craton lithosphere. This is supported by geochemical studies that indicate metasomatism of the Colorado Plateau and Wyoming craton mantle lithosphere by an aqueous fluid and/or silicate melt with a subduction signature. Through Cretaceous shallow-angle subduction, the Farallon Plate was in a position to mechanically and chemically interact with North American craton lithosphere at the time of kimberlite-lamproite magmatism, making the subduction hypothesis a viable mechanism for the genesis of these magmas. REFERENCES: McCandless, T.E., Proceedings of the 7th International Kimberlite Conference, v.2, pp.545-549, 1999; Sharp, W.E., Earth Planet. Sci. Lett., v.21, pp.351-354, 1974.

New Insights on the Structure of the Cascadia Subduction Zone from Amphibious Seismic Data

NASA Astrophysics Data System (ADS)

Janiszewski, Helen Anne

A new onshore-offshore seismic dataset from the Cascadia subduction zone was used to characterize mantle lithosphere structure from the ridge to the volcanic arc, and plate interface structure offshore within the seismogenic zone. The Cascadia Initiative (CI) covered the Juan de Fuca plate offshore the northwest coast of the United States with an ocean bottom seismometer (OBS) array for four years; this was complemented by a simultaneous onshore seismic array. Teleseismic data recorded by this array allows the unprecedented imaging of an entire tectonic plate from its creation at the ridge through subduction initiation and back beyond the volcanic arc along the entire strike of the Cascadia subduction zone. Higher frequency active source seismic data also provides constraints on the crustal structure along the plate interface offshore. Two seismic datasets were used to image the plate interface structure along a line extending 100 km offshore central Washington. These are wide-angle reflections from ship-to-shore seismic data from the Ridge-To-Trench seismic cruise and receiver functions calculated from a densely spaced CI OBS focus array in a similar region. Active source seismic observations are consistent with reflections from the plate interface offshore indicating the presence of a P-wave velocity discontinuity. Until recently, there has been limited success in using the receiver function technique on OBS data. I avoid these traditional challenges by using OBS constructed with shielding deployed in shallow water on the continental shelf. These data have quieter horizontals and avoid water- and sediment-multiple contamination at the examined frequencies. The receiver functions are consistently modeled with a velocity structure that has a low velocity zone (LVZ) with elevated P to S-wave velocity ratios at the plate interface. A similar LVZ structure has been observed onshore and interpreted as a combination of elevated pore-fluid pressures or metasediments. This new offshore result indicates that the structure may persist updip indicating the plate interface may be weak. To focus more broadly on the entire subduction system, I calculate phase velocities from teleseismic Rayleigh waves from 20-100 s period across the entire onshore-offshore array. The shear-wave velocity model calculated from these data can provide constrains on the thermal structure of the lithosphere both prior to and during subduction of the Juan de Fuca plate. Using OBS data in this period band requires removal of tilt and compliance noise, two types of water-induced noise that affect long period data. To facilitate these corrections on large seismic arrays such as the CI, an automated quality control routine was developed for selecting noise windows for the calculation of the required transfer functions. These corrections typically involve either averaging out transient signals, which requires the assumption of stationarity of the noise over the long periods of time, or laborious hand selection of noise segments. This new method calculates transfer functions based on daily time series that exclude transient signals, but allows for the investigation of long-term variation over the course of an instrument's deployment. I interpret these new shoreline-crossing phase velocity maps in terms of the tectonics associated with the Cascadia subduction system. Major findings include that oceanic plate cooling models do not explain the velocities observed beneath the Juan de Fuca plate, that slow velocities in the forearc appear to be more prevalent in areas modeled to have experienced high slip in past Cascadia megathrust earthquakes, and along strike variations in phase velocity reflect variations in arc structure and backarc tectonics.

Frictional behavior of carbonate-rich incoming sediment in the Hikurangi subduction zone

NASA Astrophysics Data System (ADS)

Rabinowitz, H. S.; Savage, H. M.; Carpenter, B.; Ikari, M.; Collettini, C.

2017-12-01

In recent years, the traditional view of the seismogenic zone has been challenged by observations of a range of seismic behaviors both above and below the depths previously considered capable of nucleating earthquakes. The Hikurangi trench is one of the few subduction zones where this transitional seismic behavior has been observed at the shallowest portions of the subduction zone, providing an opportunity to investigate the mechanical controls on seismic behavior through measurements of directly sampled sediment. To this end, an IODP cruise (March-May, 2018; Exp. 375) will recover sample from the faults that participate in this shallow seismic behavior. In order to obtain preliminary frictional characterization of the sedimentary inputs to the Hikurangi Trench, we conducted deformation experiments on samples from an ocean drill core through the incoming sediments (ODP Site 1124). The sedimentary package subducting at Hikurangi contains carbonate-rich lithologies, which have been shown to be more frictionally unstable (velocity-weakening, high healing rates) than the clays that comprise the majority of the sedimentary inputs to global subduction zones. Such frictional properties could promote seismic behavior in the shallower reaches of the subduction zone. We focus on a section of ODP Site 1124 which has a carbonate content of 40 wt% to investigate the effect of this lithology. Samples were saturated with distilled water mixed with 35 g/l sea salt. Velocity-stepping and slide-hold-slide tests were performed in multiple biaxial and triaxial deformation apparatus to investigate a range of pressures, temperatures and velocities relevant to the shallow subduction zone (σeff = 1-150 MPa, sliding velocities of 1.7 nm/s-300 μm/s, hold times of 1-1000 s, and T = 20-100 ºC). We observe transitions from velocity-strengthening to velocity-weakening behavior over these conditions which could contribute to shallow seismic behavior in the Hikurangi trench.

Various mechanisms to induce present-day shallow flat subduction and implications for the younger Earth: a numerical parameter study

NASA Astrophysics Data System (ADS)

van Hunen, Jeroen; van den Berg, Arie P.; Vlaar, Nico J.

2004-08-01

Shallow flat subduction is a relatively common feature at present-day subduction zones. Several mechanisms to explain this feature have been proposed, and can be subdivided into three groups: overthrusting of the subducting plate, subduction of a plume-generated oceanic plateau, and slab suction forces. We developed a numerical model to investigate these mechanisms and tested it through a comparison of the model results with the observations of the Peru flat slab where all three mechanisms seem to be contributing. The ratio of contributions of overthrusting continent to plateau subduction is in the range of 1:1 to 1:2, and the role of plate suction forces is likely to be significant. By applying the overthrusting continent and plateau subduction mechanisms separately, we were able to determine the most important model parameters for each of the mechanisms. Overthrusting easily results in flat subduction, and the flat slab length is primarily a function of slab age, overriding plate motion and mantle viscosity. An oceanic plateau is much less likely to cause flat subduction, and favorable conditions for flat subduction include a young slab age, long-lived plateau buoyancy after subduction, a strong mantle, and addition of slab suction forces that are large enough to further reduce the subduction dip angle, once the plateau initiates this flattening. Furthermore, we found that even though today flat subduction can be explained with the dominant model parameters within a reasonable range, for a slightly hotter, younger Earth, these flat subduction conditions are much less favorable, and so this style of subduction was probably not present in the past. This contradicts earlier predictions that flat subduction was a more wide-spread phenomenon in the early stages of plate tectonics in a younger earth.

Characterizing Shallow Seismicity at the Western End of the Middle America Trench

NASA Astrophysics Data System (ADS)

Abbott, E. R.; Brudzinski, M. R.

2011-12-01

The Middle America Trench along southwestern Mexico marks the subduction of both Cocos and Rivera plates. A wide range of seismic activity is seen all along this trench including great earthquakes with short (50-100 y) cycles, abundant microseismicity, prominent earthquake afterslip, recurring interseismic slow slip, and bands of non-volcanic tremor. Despite the fact that each of these different fault behaviors should be controlled by stress on the plate interface, no reliable relationship has been found between these phenomena as of yet. This study focuses on characterizing seismicity at the western end of the subduction zone where 4 portions of the plate interface have ruptured in 1973, 1985, 1995, and 2003. The subducted boundary between the Cocos and Rivera plates occurs beneath this region, indicated by the Manzanillo Trough, Colima Graben, inland volcanic activity, and a curious gap in tremor activity. Our data was collected by the MARS seismic array, which consists of about 50 three-component broadband seismometers deployed across Jalisco, Colima and Michoacán from January 2006 to June 2007, covering an along-strike distance of ~400 km. 18 months of data from this array was processed with Antelope for hypocentral locations of shallow (<30 km) earthquakes. To confirm the reliability of the automated locations, analyst refinement was performed on the first ~700 events, revealing little change in location and a similar clustering of events. Compilation of the resulting hypocenters reveals clusters that appear to be associated with the 2003 and 1973 megathrust earthquakes. While there are some events within the 2003 Tecomán earthquake rupture zone, more events are found inland and directly northward. Modeling of geodetic data following the neighboring 1995 Colima-Jalisco earthquake showed significant afterslip immediately downdip from that event, and there are also geodetic signatures consistent with afterslip following the 2003 event such that seismicity patterns in 2006-2007 may be influenced by ongoing afterslip. Seismicity may be concentrated north of the 2003 event as opposed to more broadly covering the region immediately downdip as it appears to follow the western edge of the Colima Graben. The prominent cluster of seismicity within the suspected 1973 rupture zone is curious both in that there is ongoing megathrust related seismicity at this point in the earthquake cycle and that it seems to match a cluster of aftershocks recorded in the days immediately after the 1973 mainshock. Finally, in comparison with observed non-volcanic tremor in the region, shallow seismicity appears to be most prominent where there are notable gaps in tremor distribution indicating that shallow earthquakes are anticorrelated with tremor locations.

P and S wave attenuation tomography of the Japan subduction zone

NASA Astrophysics Data System (ADS)

Wang, Zewei; Zhao, Dapeng; Liu, Xin; Chen, Chuanxu; Li, Xibing

2017-04-01

We determine the first high-resolution P and S wave attenuation (Q) tomography beneath the entire Japan Islands using a large number of high-quality t∗ data collected from P and S wave velocity spectra of 4222 local shallow and intermediate-depth earthquakes. The suboceanic earthquakes used in this study are relocated precisely using sP depth phases. Significant landward dipping high-Q zones are revealed clearly, which reflect the subducting Pacific slab beneath Hokkaido and Tohoku, and the subducting Philippine Sea (PHS) slab beneath SW Japan. Prominent low-Q zones are visible in the crust and mantle wedge beneath the active arc volcanoes in Hokkaido, Tohoku, and Kyushu, which reflect source zones of arc magmatism caused by fluids from the slab dehydration and corner flow in the mantle wedge. Our results also show that nonvolcanic low-frequency earthquakes (LFEs) in SW Japan mainly occur in the transition zone between a narrow low-Q belt and its adjacent high-Q zones right above the flat segment of the PHS slab. This feature suggests that the nonvolcanic LFEs are caused by not only fluid-affected slab interface but also specific conditions such as high pore pressure which is influenced by the overriding plate.

Peru Subduction Zone Seismic Experiment (PeruSZE): Preliminary Results From a Seismic Network Between Mollendo and Lake Titicaca, Peru.

NASA Astrophysics Data System (ADS)

Guy, R.; Stubailo, I.; Skinner, S.; Phillips, K.; Foote, E.; Lukac, M.; Aguilar, V.; Tavera, H.; Audin, L.; Husker, A.; Clayton, R.; Davis, P. M.

2008-12-01

This work describes preliminary results from a 50 station broadband seismic network recently installed from the coast to the high Andes in Peru. UCLA's Center for Embedded Network Sensing (CENS) and Caltech's Tectonic Observatory are collaborating with the IRD (French L'Institut de Recherche pour le Developpement) and the Institute of Geophysics, in Lima Peru in a broadband seismic experiment that will study the transition from steep to shallow slab subduction. The currently installed line has stations located above the steep subduction zone at a spacing of about 6 km. In 2009 we plan to install a line of 50 stations north from this line along the crest of the Andes, crossing the transition from steep to shallow subduction. A further line from the end of that line back to the coast, completing a U shaped array, is in the planning phase. The network is wirelessly linked using multi-hop network software designed by computer scientists in CENS in which data is transmitted from station to station, and collected at Internet drops, from where it is transmitted over the Internet to CENS each night. The instrument installation in Peru is almost finished and we have been receiving data daily from 10 stations (out of total 50) since June 2008. The rest are recording on-site while the RF network is being completed. The software system provides dynamic link quality based routing, reliable data delivery, and a disruption tolerant shell interface for managing the system from UCLA without the need to travel to Peru. The near real-time data delivery also allows immediate detection of any problems at the sites. We are building a seismic data and GPS quality control toolset that would greatly minimize the station's downtime by alerting the users of any possible problems.

Fluid and mass transfer at subduction interfaces-The field metamorphic record

NASA Astrophysics Data System (ADS)

Bebout, Gray E.; Penniston-Dorland, Sarah C.

2016-01-01

The interface between subducting oceanic slabs and the hanging wall is a structurally and lithologically complex region. Chemically disparate lithologies (sedimentary, mafic and ultramafic rocks) and mechanical mixtures thereof show heterogeneous deformation. These lithologies are tectonically juxtaposed at mm to km scales, particularly in more intensely sheared regions (mélange zones, which act as fluid channelways). This juxtaposition, commonly in the presence of a mobile fluid phase, offers up huge potential for mass transfer and related metasomatic alteration. Fluids in this setting appear capable of transporting mass over scales of kms, along flow paths with widely varying geometries and P-T trajectories. Current models of arc magmatism require km-scale migration of fluids from the interface into mantle wedge magma source regions and implicit in these models is the transport of any fluids generated in the subducting slab along and ultimately through the subduction interface. Field and geochemical studies of high- and ultrahigh-pressure metamorphic rocks elucidate the sources and compositions of fluids in subduction interfaces and the interplay between deformation and fluid and mass transfer in this region. Recent geophysical studies of the subduction interface - its thickness, mineralogy, density, and H2O content - indicate that its rheology greatly influences the ways in which the subducting plate is coupled with the hanging wall. Field investigation of the magnitude and styles of fluid-rock interaction in metamorphic rocks representing "seismogenic zone" depths (and greater) yields insight regarding the roles of fluids and elevated fluid pore pressure in the weakening of plate interface rocks and the deformation leading to seismic events. From a geochemical perspective, the plate interface contributes to shaping the "slab signature" observed in studies of the composition of arc volcanic rocks. Understanding the production of fluids with hybridized chemical/isotopic compositions could improve models aimed at identifying the relative contributions of end-member rock reservoirs through analyses of arc volcanic rocks. Production of rocks rich in hydrous minerals, along the subduction interface, could stabilize H2O to great depths in subduction zones and influence deep-Earth H2O cycling. Enhancement of decarbonation reactions and dissolution by fluid infiltration facilitated by deformation at the interface could influence the C flux from subducting slabs entering the sub-arc mantle wedge and various forearc reservoirs. In this paper, we consider records of fluid and mass transfer at localities representing various depths and structural expressions of evolving paleo-interfaces, ranging widely in structural character, the rock types involved (ultramafic, mafic, sedimentary), and the rheology of these rocks. We stress commonalities in styles of fluid and mass transfer as related to deformation style and the associated geometries of fluid mobility at subduction interfaces. Variations in thermal structure among individual margins will lead to significant differences in not only the rheology of subducting rocks, and thus seismicity, but also the profiles of devolatilization and melting, through the forearc and subarc, and the element/mineral solubilities in any aqueous fluids or silicate melts that are produced. One key factor in considering fluid and mass transfer in the subduction interface, influencing C cycling and other chemical additions to arcs, is the uncertain degree to which sub-crustal ultramafic rocks in downgoing slabs are hydrated and release H2O-rich fluids.

Stress Transfer Processes during Great Plate Boundary Thrusting Events: A Study from the Andaman and Nicobar Segments

NASA Astrophysics Data System (ADS)

Andrade, V.; Rajendran, K.

2010-12-01

The response of subduction zones to large earthquakes varies along their strike, both during the interseismic and post-seismic periods. The December 26, 2004 earthquake nucleated at 3° N latitude and its rupture propagated northward, along the Andaman-Sumatra subduction zone, terminating at 15°N. Rupture speed was estimated at about 2.0 km per second in the northern part under the Andaman region and 2.5 - 2.7 km per second under southern Nicobar and North Sumatra. We have examined the pre and post-2004 seismicity to understand the stress transfer processes within the subducting plate, in the Andaman (10° - 15° N ) and Nicobar (5° - 10° N) segments. The seismicity pattern in these segments shows distinctive characteristics associated with the outer rise, accretionary prism and the spreading ridge, all of which are relatively better developed in the Andaman segment. The Ninety East ridge and the Sumatra Fault System are significant tectonic features in the Nicobar segment. The pre-2004 seismicity in both these segments conform to the steady-state conditions wherein large earthquakes are fewer and compressive stresses dominate along the plate interface. Among the pre-2004 great earthquakes are the 1881 Nicobar and 1941 Andaman events. The former is considered to be a shallow thrust event that generated a small tsunami. Studies in other subduction zones suggest that large outer-rise tensional events follow great plate boundary breaking earthquakes due to the the up-dip transfer of stresses within the subducting plate. The seismicity of the Andaman segment (1977-2004) concurs with the steady-state stress conditions where earthquakes occur dominantly by thrust faulting. The post-2004 seismicity shows up-dip migration along the plate interface, with dominance of shallow normal faulting, including a few outer rise events and some deeper (> 100 km) strike-slip faulting events within the subducting plate. The September 13, 2002, Mw 6.5 thrust faulting earthquake at Diglipur (depth: 21 km) and the August 10, 2009, Mw 7.5 normal faulting earthquake near Coco Island (depth: 22 km), within the northern terminus of the 2004 rupture are cited as examples of the alternating pre and post earthquake stress conditions. The major pre and post 2004 clusters were associated with the Andaman Spreading Ridge (ASR). In the Nicobar segment, the most recent earthquake on June 12, 2010, Mw 7.5 (focal depth: 35 km) occurred very close to the plate boundary, through left lateral strike-slip faulting. A segment that does not feature any active volcanoes unlike its northern and southern counterparts, this part of the plate boundary has generated several right lateral strike-slip earthquakes, mostly on the Sumatra Fault System. The left-lateral strike-slip faulting associated with the June 12 event on a nearly N-S oriented fault plane consistent with the trend of the Ninety East ridge and the occasional left-lateral earthquakes prior to the 2004 mega-thrust event suggest the involvement of the Ninety East ridge in the subduction process.

Slab2 - Updated Subduction Zone Geometries and Modeling Tools

NASA Astrophysics Data System (ADS)

Moore, G.; Hayes, G. P.; Portner, D. E.; Furtney, M.; Flamme, H. E.; Hearne, M. G.

2017-12-01

The U.S. Geological Survey database of global subduction zone geometries (Slab1.0), is a highly utilized dataset that has been applied to a wide range of geophysical problems. In 2017, these models have been improved and expanded upon as part of the Slab2 modeling effort. With a new data driven approach that can be applied to a broader range of tectonic settings and geophysical data sets, we have generated a model set that will serve as a more comprehensive, reliable, and reproducible resource for three-dimensional slab geometries at all of the world's convergent margins. The newly developed framework of Slab2 is guided by: (1) a large integrated dataset, consisting of a variety of geophysical sources (e.g., earthquake hypocenters, moment tensors, active-source seismic survey images of the shallow slab, tomography models, receiver functions, bathymetry, trench ages, and sediment thickness information); (2) a dynamic filtering scheme aimed at constraining incorporated seismicity to only slab related events; (3) a 3-D data interpolation approach which captures both high resolution shallow geometries and instances of slab rollback and overlap at depth; and (4) an algorithm which incorporates uncertainties of contributing datasets to identify the most probable surface depth over the extent of each subduction zone. Further layers will also be added to the base geometry dataset, such as historic moment release, earthquake tectonic providence, and interface coupling. Along with access to several queryable data formats, all components have been wrapped into an open source library in Python, such that suites of updated models can be released as further data becomes available. This presentation will discuss the extent of Slab2 development, as well as the current availability of the model and modeling tools.

Magnitude and location of historical earthquakes in Japan and implications for the 1855 Ansei Edo earthquake

USGS Publications Warehouse

Bakun, W.H.

2005-01-01

Japan Meteorological Agency (JMA) intensity assignments IJMA are used to derive intensity attenuation models suitable for estimating the location and an intensity magnitude Mjma for historical earthquakes in Japan. The intensity for shallow crustal earthquakes on Honshu is equal to -1.89 + 1.42MJMA - 0.00887?? h - 1.66log??h, where MJMA is the JMA magnitude, ??h = (??2 + h2)1/2, and ?? and h are epicentral distance and focal depth (km), respectively. Four earthquakes located near the Japan Trench were used to develop a subducting plate intensity attenuation model where intensity is equal to -8.33 + 2.19MJMA -0.00550??h - 1.14 log ?? h. The IJMA assignments for the MJMA7.9 great 1923 Kanto earthquake on the Philippine Sea-Eurasian plate interface are consistent with the subducting plate model; Using the subducting plate model and 226 IJMA IV-VI assignments, the location of the intensity center is 25 km north of the epicenter, Mjma is 7.7, and MJMA is 7.3-8.0 at the 1?? confidence level. Intensity assignments and reported aftershock activity for the enigmatic 11 November 1855 Ansei Edo earthquake are consistent with an MJMA 7.2 Philippine Sea-Eurasian interplate source or Philippine Sea intraslab source at about 30 km depth. If the 1855 earthquake was a Philippine Sea-Eurasian interplate event, the intensity center was adjacent to and downdip of the rupture area of the great 1923 Kanto earthquake, suggesting that the 1855 and 1923 events ruptured adjoining sections of the Philippine Sea-Eurasian plate interface.

Anatomy of the western Java plate interface from depth-migrated seismic images

NASA Astrophysics Data System (ADS)

Kopp, H.; Hindle, D.; Klaeschen, D.; Oncken, O.; Reichert, C.; Scholl, D.

2009-11-01

Newly pre-stack depth-migrated seismic images resolve the structural details of the western Java forearc and plate interface. The structural segmentation of the forearc into discrete mechanical domains correlates with distinct deformation styles. Approximately 2/3 of the trench sediment fill is detached and incorporated into frontal prism imbricates, while the floor sequence is underthrust beneath the décollement. Western Java, however, differs markedly from margins such as Nankai or Barbados, where a uniform, continuous décollement reflector has been imaged. In our study area, the plate interface reveals a spatially irregular, nonlinear pattern characterized by the morphological relief of subducted seamounts and thicker than average patches of underthrust sediment. The underthrust sediment is associated with a low velocity zone as determined from wide-angle data. Active underplating is not resolved, but likely contributes to the uplift of the large bivergent wedge that constitutes the forearc high. Our profile is located 100 km west of the 2006 Java tsunami earthquake. The heterogeneous décollement zone regulates the friction behavior of the shallow subduction environment where the earthquake occurred. The alternating pattern of enhanced frictional contact zones associated with oceanic basement relief and weak material patches of underthrust sediment influences seismic coupling and possibly contributed to the heterogeneous slip distribution. Our seismic images resolve a steeply dipping splay fault, which originates at the décollement and terminates at the sea floor and which potentially contributes to tsunami generation during co-seismic activity.

Anatomy of the western Java plate interface from depth-migrated seismic images

USGS Publications Warehouse

Kopp, H.; Hindle, D.; Klaeschen, D.; Oncken, O.; Reichert, C.; Scholl, D.

2009-01-01

Newly pre-stack depth-migrated seismic images resolve the structural details of the western Java forearc and plate interface. The structural segmentation of the forearc into discrete mechanical domains correlates with distinct deformation styles. Approximately 2/3 of the trench sediment fill is detached and incorporated into frontal prism imbricates, while the floor sequence is underthrust beneath the d??collement. Western Java, however, differs markedly from margins such as Nankai or Barbados, where a uniform, continuous d??collement reflector has been imaged. In our study area, the plate interface reveals a spatially irregular, nonlinear pattern characterized by the morphological relief of subducted seamounts and thicker than average patches of underthrust sediment. The underthrust sediment is associated with a low velocity zone as determined from wide-angle data. Active underplating is not resolved, but likely contributes to the uplift of the large bivergent wedge that constitutes the forearc high. Our profile is located 100 km west of the 2006 Java tsunami earthquake. The heterogeneous d??collement zone regulates the friction behavior of the shallow subduction environment where the earthquake occurred. The alternating pattern of enhanced frictional contact zones associated with oceanic basement relief and weak material patches of underthrust sediment influences seismic coupling and possibly contributed to the heterogeneous slip distribution. Our seismic images resolve a steeply dipping splay fault, which originates at the d??collement and terminates at the sea floor and which potentially contributes to tsunami generation during co-seismic activity. ?? 2009 Elsevier B.V.

Asymmetric Subductions in an Asymmetric Earth: Geodynamics and Numerical Modeling

NASA Astrophysics Data System (ADS)

Dal Zilio, L.; Ficini, E.; Doglioni, C.; Gerya, T.

2016-12-01

The driving mechanism of plate tectonics is still controversial. Moreover, mantle kinematics is still poorly constrained due to the limited information available on its composition, thermal state, and physical parameters. The net rotation of the lithosphere, or so-called W-ward drift, however, indicates a decoupling of the plates relative to the underlying asthenosphere at about 100-200 km depth in the Low-Velocity Zone and a relative "E-ward" mantle counterflow. This mantle flow can account for a number of tectonic asymmetries on subduction dynamics such as steep versus shallow slab dip, diverging versus converging subduction hinge, low versus high topography of mountain belts, etc. This asymmetry is generally interpreted to reflect the age-dependent negative buoyancy of the subducting lithosphere. However, slab dip is insensitive to the age of the lithosphere. Here we investigate the role of mantle flow in controlling subduction dynamics using a high-resolution rheologically consistent two-dimensional numerical modeling. Results show the evolution of a subducting oceanic plate beneath a continent: when the subducting plate is dipping in opposite direction with respect to the mantle flow, the slab is sub-vertically deflected by the mantle flow, thus leading the coeval development of a back-arc basin. In contrast, agreement between mantle flow and dipping of the subducting slab relieves shallow dipping subduction zone, which in turn controls the development of a pronounced topography. Moreover, this study confirms that the age of the subducting oceanic lithosphere (i.e. its negative buoyancy) has a second order effect on the dip angle of the slab and, more generally, on subduction dynamics. Our numerical experiments show strong similarities to the observed evolution of subduction zone worldwide and demonstrate that the possibility of a horizontal mantle flow is universally valid.

Fluid Pressure in the Shallow Plate Interface at the Nankai Trough Subduction Zone

NASA Astrophysics Data System (ADS)

Tobin, H. J.; Saffer, D.

2003-12-01

The factors controlling the occurrence, magnitude, and other characteristics of great earthquakes is a fundamental outstanding question in fault physics. Pore fluid pressure is perhaps the most critical yet poorly known parameter governing the strength and seismogenic character of plate boundary faults, but unfortunately cannot be directly inferred through available geophysical sensing methods. Moreover, true in situ fluid pressure has proven difficult to measure even in boreholes. At the Nankai Trough, several hundred meters of sediment are subducted beneath the frontal portion of the accretionary prism. The up-dip portion of the plate interface is therefore hosted in these fine-grained marine sedimentary rocks. ODP Leg 190 and 196 showed that these rapidly-loaded underthrust sediments are significantly overpressured near the deformation front. Here, we attempt to quantitatively infer porosity, pore pressure, and effective normal stress at the plate interface at depths currently inaccessible to drilling. Using seismic reflection interval velocity calibrated at the boreholes to porosity, we quantitatively infer pore pressure to ˜ 20 km down-dip of the deformation front, to a plate interface depth of ˜ 6 km. We have developed a Nankai-specific velocity-porosity transform using ODP cores and logs. We use this function to derive a porosity profile for each of two down-dip seismic sections extracted from a 3-D dataset from the Cape Muroto region. We then calculate pore fluid pressure and effective vertical (fault-normal) stress for the underthrust sediment section using a compaction disequilibrium approach and core-based consolidation test data. Because the pore fluid pressure at the fault interface is likely controlled by that of the top of the underthrust section, this calculation represents a quantitative profile of effective stress and pore pressure at the plate interface. Results show that seismic velocity and porosity increase systematically downdip in the underthrust section, but the increase is suppressed relative to that expected from normally consolidating sediments. The computed pore pressure increases landward from an overpressure ratio (λ * = hydrostatic pressure divided by the lithostatic overburden) of ˜ 0.6 at the deformation front to ˜ 0.77 where sediments have been subducted 15 km. The results of this preliminary analysis suggest that a 3-dimensional mapping of predicted effective normal stress in the seismic data volume is possible.

Empirical ground-motion relations for subduction-zone earthquakes and their application to Cascadia and other regions

USGS Publications Warehouse

Atkinson, G.M.; Boore, D.M.

2003-01-01

Ground-motion relations for earthquakes that occur in subduction zones are an important input to seismic-hazard analyses in many parts of the world. In the Cascadia region (Washington, Oregon, northern California, and British Columbia), for example, there is a significant hazard from megathrust earthquakes along the subduction interface and from large events within the subducting slab. These hazards are in addition to the hazard from shallow earthquakes in the overlying crust. We have compiled a response spectra database from thousands of strong-motion recordings from events of moment magnitude (M) 5-8.3 occurring in subduction zones around the world, including both interface and in-slab events. The 2001 M 6.8 Nisqually and 1999 M 5.9 Satsop earthquakes are included in the database, as are many records from subduction zones in Japan (Kyoshin-Net data), Mexico (Guerrero data), and Central America. The size of the database is four times larger than that available for previous empirical regressions to determine ground-motion relations for subduction-zone earthquakes. The large dataset enables improved determination of attenuation parameters and magnitude scaling, for both interface and in-slab events. Soil response parameters are also better determined by the data. We use the database to develop global ground-motion relations for interface and in-slab earthquakes, using a maximum likelihood regression method. We analyze regional variability of ground-motion amplitudes across the global database and find that there are significant regional differences. In particular, amplitudes in Cascadia differ by more than a factor of 2 from those in Japan for the same magnitude, distance, event type, and National Earthquake Hazards Reduction Program (NEHRP) soil class. This is believed to be due to regional differences in the depth of the soil profile, which are not captured by the NEHRP site classification scheme. Regional correction factors to account for these differences are proposed for Cascadia and Japan. The results of this study differ significantly from previous analyses based on more limited data and have important implications for seismic-hazard analysis. The ground-motion relations predict that a great megathrust earthquake (M ???8) at a fault distance of about 100 km would produce pseudoacceleration (PSA), 5% damped, horizontal component on soil sites of about 110 cm/sec2 at 0.5 Hz, 660 cm/sec2 at 2.5 Hz, and 410 cm/sec2 at 5 Hz, with a peak ground acceleration of about 180 cm/ sec2 . These damaging levels of motion would be experienced over a very large area, corresponding to a rectangular area about 300 km wide by 500 km long. Large in-slab events (M 7.5) would produce even higher PSA values within 100 km of the fault, but the in-slab motions attenuate much more rapidly with distance. Thus the hazard posed by moderate to large in-slab events such as the 2001 Nisqually earthquake is modest compared to that of a Cascadia megathrust earthquake of M ???8, in terms of the area that would experience damaging levels of ground motion.

Structure of the tsunamigenic plate boundary and low-frequency earthquakes in the southern Ryukyu Trench

PubMed Central

Arai, Ryuta; Takahashi, Tsutomu; Kodaira, Shuichi; Kaiho, Yuka; Nakanishi, Ayako; Fujie, Gou; Nakamura, Yasuyuki; Yamamoto, Yojiro; Ishihara, Yasushi; Miura, Seiichi; Kaneda, Yoshiyuki

2016-01-01

It has been recognized that even weakly coupled subduction zones may cause large interplate earthquakes leading to destructive tsunamis. The Ryukyu Trench is one of the best fields to study this phenomenon, since various slow earthquakes and tsunamis have occurred; yet the fault structure and seismic activity there are poorly constrained. Here we present seismological evidence from marine observation for megathrust faults and low-frequency earthquakes (LFEs). On the basis of passive observation we find LFEs occur at 15–18 km depths along the plate interface and their distribution seems to bridge the gap between the shallow tsunamigenic zone and the deep slow slip region. This suggests that the southern Ryukyu Trench is dominated by slow earthquakes at any depths and lacks a typical locked zone. The plate interface is overlaid by a low-velocity wedge and is accompanied by polarity reversals of seismic reflections, indicating fluids exist at various depths along the plate interface. PMID:27447546

Linkages Between the Megathrust and Upper-plate Deformation: Lessons From the Deformational Dichotomy of the 2016 Kaikoura New Zealand Earthquake

NASA Astrophysics Data System (ADS)

Furlong, K. P.; Herman, M. W.

2017-12-01

Following the 2016 Mw 7.8 Kaikoura earthquake, the nature of the coseismic rupture was unclear. Seismological and tsunami evidence pointed to significant involvement of the subduction megathrust, while geodetic and field observations pointed to a shallow set of intra-crustal faults as the main participants during the earthquake. It now appears that the Kaikoura earthquake produced synchronous faulting on the plate boundary subduction interface - the megathrust - and on a suite of crustal faults above the rupture zone in the overlying plate. This Kaikoura-style earthquake, involving synchronous ruptures on multiple components of the plate boundary, may be an important mode of plate boundary deformation affecting seismic hazard along subduction zones. Here we propose a model to explain how these upper-plate faults are loaded during the periods between megathrust earthquakes and subsequently can rupture synchronously with the megathrust. Between megathrust earthquakes, horizontal compression, driven by plate convergence, locks the upper-plate faults, particularly those at higher angles to the convergence direction and the oblique plate motion of the subducting Pacific plate deforms the upper-plate in bulk shear. During the time interval of megathrust rupture, two things happen which directly affect the stress conditions acting on these upper-plate faults: (1) slip on the megathrust and the associated `rebound' of the upper plate reduces the compressive or normal stress acting on the upper plate faults, and (2) the base of the upper plate faults (and the upper plate itself) is decoupled from the slab in the region above rupture area. The reduction in normal stress acting on these faults increases their Coulomb Stress state to strongly favor strike-slip fault slip, and the basal decoupling of the upper plate allows it to undergo nearly complete stress recovery in that region; enabling the occurrence of very large offsets on these faults - offsets that exceed the slip on the plate interface. With these results it is clear that the 2016 Kaikoura NZ earthquake represents a mode of subduction zone rupture that must be considered in other regions.

Tectono-seismic characteristics of faults in the shallow portion of an accretionary prism

NASA Astrophysics Data System (ADS)

Hirono, Tetsuro; Ishikawa, Tsuyoshi

2018-01-01

To understand the tectono-seismic evolution of faults in the shallow part of a subduction-accretion system, we examined major faults in a fossil accretionary prism, the Emi Group (Hota Group), Boso Peninsula, Japan, by performing multiple structural, geochemical, and mineralogical analyses. Because the strata are relatively shallow (burial depth, 1-4 km), early stage deformation related to subduction, accretion, and uplifting processes is well preserved in three dominant fault zones. On the basis of both previous findings and our geochemical and mineralogical results, we inferred that early stage faulting in a near-trench setting under high pore fluid pressure and second stage faulting at relatively deep along subduction corresponded to aseismic deformations, as shown by velocity strengthening characteristics; and during late stage faulting, probably in association with accretion and uplift processes, a high-temperature fluid, revealed by a geochemical temperature proxy, triggered fault weakening by a thermal pressurization mechanism, and potentially led to the generation of a tsunami.

15 Years Of Ecuadorian-French Research Along The Ecuadorian Subduction Zone

NASA Astrophysics Data System (ADS)

Charvis, P.

2015-12-01

The Ecuadorian segment of the Nazca/South America subduction zone is an outstanding laboratory to study the seismic cycle. Central Ecuador where the Carnegie ridge enters the subduction marks a transition between a highly coupled segment that hosted one of the largest seismic sequence during the 20thcentury and a ~1200-km long weakly coupled segment encompassing southern Ecuador and northern Peru. A shallow dipping subduction interface and a short trench-coast line distance ranging from 45 to 80 km, together with La Plata Island located only 33 km from the trench axis, allow to document subduction processes in the near field with an exceptional resolution. Since 2000, a close cooperation between the Institute of Geophysics (Quito), INOCAR (Oceanographic Institute of the Ecuadorian Navy) with French groups allowed us to conduct up to 6 marine geophysics cruises to survey the convergent margin and jointly develop dense GPS and seismological networks. This fruitful collaboration now takes place in the framework of an International Joint Laboratory "Earthquakes and Volcanoes in the Northern Andes" (LMI SVAN), which eases coordinating research projects and exchanges of Ecuadorian and French scientists and students. This long-term investigation has already provided a unique view on the structure of the margin, which exhibits a highly variable subduction channel along strike. It allowed us to evidence the contrast between creeping and coupled segments of subduction at various scale, and the existence of large continental slivers whose motion accommodates the obliquity of the Nazca/South America convergence. Finally, we could evidence the first Slow Slip Events (SSE) that oppositely to most SSE documented so far, are accompanied with intense micro-seismicity. The recent support of the French National Research Agency and the Ecuadorian Agency for Sciences and Technology (Senescyt) will enable us to integrate the already obtained results, in an attempt to develop an earthquake forecast model for Ecuador.

What role did the Hikurangi subduction zone play in the M7.8 Kaikoura earthquake?

NASA Astrophysics Data System (ADS)

Wallace, L. M.; Hamling, I. J.; Kaneko, Y.; Fry, B.; Clark, K.; Bannister, S. C.; Ellis, S. M.; Francois-Holden, C.; Hreinsdottir, S.; Mueller, C.

2017-12-01

The 2016 M7.8 Kaikoura earthquake ruptured at least a dozen faults in the northern South Island of New Zealand, within the transition from the Hikurangi subduction zone (in the North Island) to the transpressive Alpine Fault (in the central South Island). The role that the southern end of the Hikurangi subduction zone played (or did not play) in the Kaikoura earthquake remains one of the most controversial aspects of this spectacularly complex earthquake. Investigations using near-field seismological and geodetic data suggest a dominantly crustal faulting source for the event, while studies relying on teleseismic data propose that a large portion of the moment release is due to rupture of the Hikurangi subduction interface beneath the northern South Island. InSAR and GPS data also show that a large amount of afterslip (up to 0.5 m) occurred on the subduction interface beneath the crustal faults that ruptured in the M7.8 earthquake, during the months following the earthquake. Modeling of GPS velocities for the 20 year period prior to the earthquake indicate that interseismic coupling was occurring on the Hikurangi subduction interface beneath the northern South Island, in a similar location to the suggested coseismic and postseismic slip on the subduction interface. We will integrate geodetic, seismological, tsunami, and geological observations in an attempt to balance the seemingly conflicting views from local and teleseismic data regarding the role that the southern Hikurangi subduction zone played in the earthquake. We will also discuss the broader implications of the observed coseismic and postseismic deformation for understanding the kinematics of the southern termination of the Hikurangi subduction zone, and its role in the transition from subduction to strike-slip in the central New Zealand region.

Seismic and aseismic slip on the ``uncoupled'' Tonga subduction megathrust

NASA Astrophysics Data System (ADS)

Beavan, R. J.; Wang, X.; Bevis, M. G.; Kautoke, R'

2010-12-01

The Tonga subduction zone has been a type example of a weakly coupled subduction interface since soon after the birth of plate tectonics. Yet in the September 2009 double earthquake, the northern Tonga subduction interface failed in a great Mw 8 earthquake that was probably dynamically triggered by a Mw 8 extensional intraplate earthquake in the outer trench slope region of the incoming Pacific Plate. There are some discrepancies between models of the September 2009 doublet derived from seismic data and those derived from geodetic and DART tsunami data, in particular about which fault plane failed in the intraplate earthquake. In this presentation we explore how well the geodetic and tsunami data can be fit using the alternative fault plane. We also present new GPS data that show the subduction interface is continuing to slip faster than its 1996-2005 “long-term” rate, and we speculate on what this means for the mechanisms by which interplate slip is accommodated at the Tonga subduction zone.

Modelling the interplate domain in thermo-mechanical simulations of subduction: Critical effects of resolution and rheology, and consequences on wet mantle melting

NASA Astrophysics Data System (ADS)

Arcay, Diane

2017-08-01

The present study aims at better deciphering the different mechanisms involved in the functioning of the subduction interplate. A 2D thermo-mechanical model is used to simulate a subduction channel, made of oceanic crust, free to evolve. Convergence at constant rate is imposed under a 100 km thick upper plate. Pseudo-brittle and non-Newtonian behaviours are modelled. The influence of the subduction channel strength, parameterized by the difference in activation energy between crust and mantle (ΔEa) is investigated to examine in detail the variations in depth of the subduction plane down-dip extent, zcoup . First, simulations show that numerical resolution may be responsible for an artificial and significant shallowing of zcoup if the weak crustal layer is not correctly resolved. Second, if the age of the subducting plate is 100 Myr, subduction occurs for any ΔEa . The stiffer the crust is, that is, the lower ΔEa is, the shallower zcoup is (60 km depth if ΔEa = 20 kJ/mol) and the hotter the fore-arc base is. Conversely, imposing a very weak subduction channel (ΔEa > 135 J/mol) leads there to an extreme mantle wedge cooling and inhibits mantle melting in wet conditions. Partial kinematic coupling at the fore-arc base occurs if ΔEa = 145 kJ/mol. If the incoming plate is 20 Myr old, subduction can occur under the conditions that the crust is either stiff and denser than the mantle, or weak and buoyant. In the latter condition, cold crust plumes rise from the subduction channel and ascend through the upper lithosphere, triggering (1) partial kinematic coupling under the fore-arc, (2) fore-arc lithosphere cooling, and (3) partial or complete hindrance of wet mantle melting. zcoup then ranges from 50 to more than 250 km depth and is time-dependent if crust plumes form. Finally, subduction plane dynamics is intimately linked to the regime of subduction-induced corner flow. Two different intervals of ΔEa are underlined: 80-120 kJ/mol to reproduce the range of slab surface temperature inferred from geothermometry, and 10-40 kJ/mol to reproduce the shallow hot mantle wedge core inferred from conditions of last equilibration of near-primary arc magmas and seismic tomographies. Therefore, an extra process controlling mantle wedge dynamics is needed to satisfy simultaneously the aforementioned observations. A mantle viscosity reduction, by a factor 4-20, caused by metasomatism in the mantle wedge is proposed. From these results, I conclude that the subduction channel down-dip extent, zcoup , should depend on the subduction setting, to be consistent with the observed variability of sub-arc depths of the subducting plate surface.

Seismicity and Structure of the Incoming Pacific Plate Subducting into the Japan Trench off Miyagi

NASA Astrophysics Data System (ADS)

Obana, K.; Fujie, G.; Kodaira, S.; Takahashi, T.; Yamamoto, Y.; Sato, T.; Yamashita, M.; Nakamura, Y.; Miura, S.

2015-12-01

Stresses within the oceanic plate in trench axis and outer-rise region have been characterized by shallow extension and deep compression due to the bending of the plate subducting into the trench. The stress state within the incoming/subducting oceanic plate is an important factor not only for the occurrence of shallow intraplate normal-faulting earthquakes in the trench-outer rise region but also the hydration of the oceanic plate through the shallow normal faults cutting the oceanic lithosphere. We investigate seismic velocity structure and stress state within the incoming/subducting Pacific Plate in the Japan Trench based on the OBS aftershock observations for the December 2012 intraplate doublet, which consists of a deep reverse faulting (Mw 7.2) and a shallow normal faulting (Mw 7.2) earthquake, in the Japan Trench off Miyagi. Hypocenter locations and seismic velocity structures were estimated from the arrival time data of about 3000 earthquakes by using double-difference tomography method (Zhang and Thurber, 2003). Also, focal mechanisms were estimated from first motion polarities by using the program HASH by Hardebeck and Shearer (2002). The results show that the earthquakes occurred mainly within the oceanic crust and the uppermost mantle. The deepest event was located at a depth of about 60 km. Focal mechanisms of the earthquakes shallower than a depth of 40 km indicate normal-faulting with T-axis normal to the trench. On the other hand, first motion polarities of the events at depths between 50 and 60 km can be explained a reverse faulting. The results suggest that the neutral plane of the stress between shallow extension and deep compression locates at 40 to 50 km deep. Seismic velocity structures indicate velocity decrease in the oceanic mantle toward the trench. Although the velocity decrease varies with locations, the results suggest the bending-related structure change could extend to at least about 15 km below the oceanic Moho in some locations.

Geological perspectives of shallow slow earthquakes deduced from deformation in subduction mélanges

NASA Astrophysics Data System (ADS)

Ujiie, K.; Saishu, H.; Kinoshita, T.; Nishiyama, N.; Otsubo, M.; Ohta, K.; Yamashita, Y.; Ito, Y.

2017-12-01

Shallow (< 15 km depth) slow earthquakes are important to understand, as they occur along the subduction thrust where devastating tsunamis are generated. Geophysical studies have revealed that shallow slow earthquakes are not restricted to specific temperature conditions and depths but occur in regions of high fluid pressure. In the Nankai subduction zone, the shallow slow slip appears to trigger tremor and very-low-frequency-earthquake. However, the geologic perspectives for shallow slow earthquakes remain enigmatic. The Makimine mélange in the Late Cretaceous Shimanto accretionary complex of southwest Japan was formed during the subduction of young oceanic plate. Within the mélange, the quartz-filled veins and viscous shear zones are concentrated in the zones of 10 to 60 m-thick. The veins consist of shear veins showing low-angle thrust or normal faulting mechanisms and extension veins parallel or at high angle to mélange foliation. The geometrical relationship between shear and extension veins indicates that shear slip and tensile fracturing occur by small differential stress under elevated fluid pressure. The shear and extension veins typically show crack-seal textures defined by the solid inclusions bands. The time scale of each crack-seal event, which is determined from the quartz kinetics considering inclusion band spacing and vein length, is a few years. The shear slip increments estimated from the spacing of inclusions bands at dilational jogs are 0.1 mm. The viscous shear is accommodated by pressure solution creep and consistently shows low-angle thrust shear sense. These geologic features are suggested to explain seismogenic environment for shallow slow earthquakes. The shear veins and viscous shear zones showing low-angle thrust faulting mechanism could represent episodic tremor and slip, while the shear veins showing low-angle normal faulting mechanism may represent the tremor that occurred after the passage of slow slip front.

Time-Varying Upper-Plate Deformation during the Megathrust Subduction Earthquake Cycle

NASA Astrophysics Data System (ADS)

Furlong, Kevin P.; Govers, Rob; Herman, Matthew

2015-04-01

Over the past several decades of the WEGENER era, our abilities to observe and image the deformational behavior of the upper plate in megathrust subduction zones has dramatically improved. Several intriguing inferences can be made from these observations including apparent lateral variations in locking along subduction zones, which differs from interseismic to coseismic periods; the significant magnitude of post-earthquake deformation (e.g. following the 20U14 Mw Iquique, Chile earthquake, observed on-land GPS post-EQ displacements are comparable to the co-seismic displacements); and incompatibilities between rates of slip deficit accumulation and resulting earthquake co-seismic slip (e.g. pre-Tohoku, inferred rates of slip deficit accumulation on the megathrust significantly exceed slip amounts for the ~ 1000 year recurrence.) Modeling capabilities have grown from fitting simple elastic accumulation/rebound curves to sparse data to having spatially dense continuous time series that allow us to infer details of plate boundary coupling, rheology-driven transient deformation, and partitioning among inter-earthquake and co-seismic displacements. In this research we utilize a 2D numerical modeling to explore the time-varying deformational behavior of subduction zones during the earthquake cycle with an emphasis on upper-plate and plate interface behavior. We have used a simplified model configuration to isolate fundamental processes associated with the earthquake cycle, rather than attempting to fit details of specific megathrust zones. Using a simple subduction geometry, but realistic rheologic layering we are evaluating the time-varying displacement and stress response through a multi-earthquake cycle history. We use a simple model configuration - an elastic subducting slab, an elastic upper plate (shallower than 40 km), and a visco-elastic upper plate (deeper than 40 km). This configuration leads to an upper plate that acts as a deforming elastic beam at inter-earthquake loading times and rates with a viscously relaxed regime at depths greater than 40 km. Analyses of our preliminary model results lead to the following: 1. Co-seismic stress transfer from the unloading elastic layer (shallow) into an elastically loading visco-elastic layer (deeper) - extends ~ 100 km inboard of locked zone. This stress transfer affects both coseismic and post-seismic surface displacements. 2. Post-seismic response of upper plate involves seaward motion for initial 10-20 years (~ 2 Maxwell times) after EQ. This occurs in spite of there being no slip on locked plate boundary - i.e. this is not plate boundary after-slip but rather is a consequence of stress relaxation in co-seismically loaded visco-elastic layer. However standard inversions of the surface displacement field would indicate significant after-slip along the locked plate interface. 3. By approximately 80 years (8 Maxwell times) system has returned to simple linear displacement pattern - the expected behavior for a shortening elastic beam. Prior to that time, the surface (observable) displacement pattern changes substantially over time and would result in an apparent temporal variation in coupling - from near-zero coupling to fully locked over ~ 80 years post-earthquake. These preliminary results indicate that care is needed in interpreting observed surface displacement fields from GPS, InSAR, etc. during the interseismic period. temporal variations in crustal deformation observed in regions such as the recent Tohoku, Maule, and Iquique megathrust events which are ascribed to fault plane after-slip may in fact reflect processes associated with re-equilibration of the visco-elastic subduction system.

A silent slip event on the deeper Cascadia subduction interface.

PubMed

Dragert, G; Wang, K; James, T S

2001-05-25

Continuous Global Positioning System sites in southwestern British Columbia, Canada, and northwestern Washington state, USA, have been moving landward as a result of the locked state of the Cascadia subduction fault offshore. In the summer of 1999, a cluster of seven sites briefly reversed their direction of motion. No seismicity was associated with this event. The sudden displacements are best explained by approximately 2 centimeters of aseismic slip over a 50-kilometer-by-300-kilometer area on the subduction interface downdip from the seismogenic zone, a rupture equivalent to an earthquake of moment magnitude 6.7. This provides evidence that slip of the hotter, plastic part of the subduction interface, and hence stress loading of the megathrust earthquake zone, can occur in discrete pulses.

Rheological Properties of Natural Subduction Zone Interface: Insights from "Digital" Griggs Experiments

NASA Astrophysics Data System (ADS)

Ioannidi, P. I.; Le Pourhiet, L.; Moreno, M.; Agard, P.; Oncken, O.; Angiboust, S.

2017-12-01

The physical nature of plate locking and its relation to surface deformation patterns at different time scales (e.g. GPS displacements during the seismic cycle) can be better understood by determining the rheological parameters of the subduction interface. However, since direct rheological measurements are not possible, finite element modelling helps to determine the effective rheological parameters of the subduction interface. We used the open source finite element code pTatin to create 2D models, starting with a homogeneous medium representing shearing at the subduction interface. We tested several boundary conditions that mimic simple shear and opted for the one that best describes the Grigg's type simple shear experiments. After examining different parameters, such as shearing velocity, temperature and viscosity, we added complexity to the geometry by including a second phase. This arises from field observations, where shear zone outcrops are often composites of multiple phases: stronger crustal blocks embedded within a sedimentary and/or serpentinized matrix have been reported for several exhumed subduction zones. We implemented a simplified model to simulate simple shearing of a two-phase medium in order to quantify the effect of heterogeneous rheology on stress and strain localization. Preliminary results show different strength in the models depending on the block-to-matrix ratio. We applied our method to outcrop scale block-in-matrix geometries and by sampling at different depths along exhumed former subduction interfaces, we expect to be able to provide effective friction and viscosity of a natural interface. In a next step, these effective parameters will be used as input into seismic cycle deformation models in an attempt to assess the possible signature of field geometries on the slip behaviour of the plate interface.

A Benchmarking setup for Coupled Earthquake Cycle - Dynamic Rupture - Tsunami Simulations

NASA Astrophysics Data System (ADS)

Behrens, Joern; Bader, Michael; van Dinther, Ylona; Gabriel, Alice-Agnes; Madden, Elizabeth H.; Ulrich, Thomas; Uphoff, Carsten; Vater, Stefan; Wollherr, Stephanie; van Zelst, Iris

2017-04-01

We developed a simulation framework for coupled physics-based earthquake rupture generation with tsunami propagation and inundation on a simplified subduction zone system for the project "Advanced Simulation of Coupled Earthquake and Tsunami Events" (ASCETE, funded by the Volkswagen Foundation). Here, we present a benchmarking setup that can be used for complex rupture models. The workflow begins with a 2D seismo-thermo-mechanical earthquake cycle model representing long term deformation along a planar, shallowly dipping subduction zone interface. Slip instabilities that approximate earthquakes arise spontaneously along the subduction zone interface in this model. The absolute stress field and material properties for a single slip event are used as initial conditions for a dynamic earthquake rupture model.The rupture simulation is performed with SeisSol, which uses an ADER discontinuous Galerkin discretization scheme with an unstructured tetrahedral mesh. The seafloor displacements resulting from this rupture are transferred to the tsunami model with a simple coastal run-up profile. An adaptive mesh discretizing the shallow water equations with a Runge-Kutta discontinuous Galerkin (RKDG) scheme subsequently allows for an accurate and efficient representation of the tsunami evolution and inundation at the coast. This workflow allows for evaluation of how the rupture behavior affects the hydrodynamic wave propagation and coastal inundation. We present coupled results for differing earthquake scenarios. Examples include megathrust only ruptures versus ruptures with splay fault branching off the megathrust near the surface. Coupling to the tsunami simulation component is performed either dynamically (time dependent) or statically, resulting in differing tsunami wave and inundation behavior. The simplified topographical setup allows for systematic parameter studies and reproducible physical studies.

Transitory signal in the Copiapó region of the Chilean subduction interface

NASA Astrophysics Data System (ADS)

Klein, E.; Vigny, C.; Zigone, D.; Duputel, Z.; Doubre, C.; Meneses, G.

2017-12-01

After the occurrence of the last 3 major Chilean earthquakes of Maule (2010, Mw8.8), Iquique (2014, Mw8.1) and Illapel (2015, Mw8.3), the North-Central Chile region of Atacama has now great odds of being the scene of the next great megathrust event along the Chilean subduction zone. In this region, a clear gap of deep seismicity is observed and the historical seismicity reveals a seismic cycle of about a century, after the occurrence of two Mw 8.5, in 1819 and 1922. Consistently, the coupling pattern derived from interseismic GPS velocities, estimated between 2010 and 2014, highlights two highly coupled segments of Atacama and Chañaral, separated by the Low Coupled Zone of Baranquilla. Yet, the last two years of survey GPS observations show a clear change of trend compared to interseismic velocities, at all sites located between Baranquilla (27.5°S) and Chañaral (26.5°S). We estimate the offsets between the 2016 positions extrapolated from interseismic velocities and measured positions, which reveal a horizontal westward diverging motion spatially coherent over the whole area and of several centimeters in amplitude. Additionally, a non-negligible vertical, uplifting offset is observed. The first analysis of data coming from the rare broadband stations of the region reveals some signals that could be associated with non-volcanic tremors but that are not occurring exactly at the same time as the GPS anomaly. South of this area, the Copiapó Ridge enters the subduction and repeated swarms have occurred offshore in 1973, 1976, 2006. Nevertheless, the pattern observed by GPS suggests an inland and deep source, which cannot directly be associated with shallow seismic swarms nor aseismic slip on the interface.

Deep Structure of Northern Apennines Subduction Orogen (Italy) as Revealed by a Joint Interpretation of Passive and Active Seismic Data

NASA Astrophysics Data System (ADS)

Piana Agostinetti, Nicola; Faccenna, Claudio

2018-05-01

The Apennines is a well-studied orogeny formed by the accretion of continental slivers during the subduction of the Adriatic plate, but its deep structure is still a topic of controversy. Here we illuminated the deep structure of the Northern Apennines belt by combining results from the analysis of active seismic (CROP03) and receiver function data. The result from combining these two approaches provides a new robust view of the structure of the deep crust/upper mantle, from the back-arc region to the Adriatic subduction zone. Our analysis confirms the shallow Moho depth beneath the back-arc region and defines the top of the downgoing plate, showing that the two plates separate at depth about 40 km closer to the trench than reported in previous reconstructions. This spatial relationship has profound implications for the geometry of the shallow subduction zone and of the mantle wedge, by the amount of crustal material consumed at trench.

Elastic Properties of Subduction Zone Materials in the Large Shallow Slip Environment for the Tohoku 2011 Earthquake: Laboratory data from JFAST Core Samples

NASA Astrophysics Data System (ADS)

Jeppson, T.; Tobin, H. J.

2014-12-01

The 11 March 2011 Tohoku-Oki earthquake (Mw=9.0) produced large displacements of ~50 meters near the Japan Trench. In order to understand earthquake propagation and slip stabilization in this environment, quantitative values of the real elastic properties of fault zones and their surrounding wall rock material is crucial. Because elastic and mechanical properties of faults and wallrocks are controlling factors in fault strength, earthquake generation and propagation, and slip stabilization, an understanding of these properties and their depth dependence is essential to understanding and accurately modeling earthquake rupture. In particular, quantitatively measured S-wave speeds, needed for estimation of elastic properties, are scarce in the literature. We report laboratory ultrasonic velocity measurements performed at elevated pressures, as well as the calculated dynamic elastic moduli, for samples of the rock surrounding the Tohoku earthquake principal fault zone recovered by drilling during IODP Expedition 343, Japan Trench Fast Drilling Project (JFAST). We performed measurements on five samples of gray mudstone from the hanging wall and one sample of underthrust brown mudstone from the footwall. We find P- and S-wave velocities of 2.0 to 2.4 km/s and 0.7 to 1.0 km/s, respectively, at 5 MPa effective pressure. At the same effective pressure, the hanging wall samples have shear moduli ranging from 1.4 to 2.2 GPa and the footwall sample has a shear modulus of 1.0 GPa. While these values are perhaps not surprising for shallow, clay-rich subduction zone sediments, they are substantially lower than the 30 GPa commonly assumed for rigidity in earthquake rupture and propagation models [e.g., Ide et al., 1993; Liu and Rice, 2005; Loveless and Meade, 2011]. In order to better understand the elastic properties of shallow subduction zone sediments, our measurements from the Japan Trench are compared to similar shallow drill core samples from the Nankai Trough, Costa Rica, Cascadia, and Barbados ridge subduction zones. We find that shallow subduction zone sediments in general have similarly low rigidity. These data provide important ground-truth values that can be used to parameterize fault slip models addressing the problem of shallow, tsunamigenic propagation of megathrust earthquakes.

Locking of the Chile subduction zone controlled by fluid pressure before the 2010 earthquake

NASA Astrophysics Data System (ADS)

Moreno, Marcos; Haberland, Christian; Oncken, Onno; Rietbrock, Andreas; Angiboust, Samuel; Heidbach, Oliver

2014-04-01

Constraints on the potential size and recurrence time of strong subduction-zone earthquakes come from the degree of locking between the down-going and overriding plates, in the period between large earthquakes. In many cases, this interseismic locking degree correlates with slip during large earthquakes or is attributed to variations in fluid content at the plate interface. Here we use geodetic and seismological data to explore the links between pore-fluid pressure and locking patterns at the subduction interface ruptured during the magnitude 8.8 Chile earthquake in 2010. High-resolution three-dimensional seismic tomography reveals variations in the ratio of seismic P- to S-wave velocities (Vp/Vs) along the length of the subduction-zone interface. High Vp/Vs domains, interpreted as zones of elevated pore-fluid pressure, correlate spatially with parts of the plate interface that are poorly locked and slip aseismically. In contrast, low Vp/Vs domains, interpreted as zones of lower pore-fluid pressure, correlate with locked parts of the plate interface, where unstable slip and earthquakes occur. Variations in pore-fluid pressure are caused by the subduction and dehydration of a hydrothermally altered oceanic fracture zone. We conclude that variations in pore-fluid pressure at the plate interface control the degree of interseismic locking and therefore the slip distribution of large earthquake ruptures.

Geodesy cannot presently detect the up-dip limit of frictional locking on megathrusts

NASA Astrophysics Data System (ADS)

Almeida, R. V.; Lindsey, E. O.; Bradley, K.; Hubbard, J.; Sathiakumar, S.; Malick, R.; Hill, E.

2017-12-01

Most discussions of interseismic behavior on megathrust faults focus on whether they are frictionally locked or creeping. Unfortunately, many geodetic studies of subduction zone megathrusts equate fault coupling with frictional locking. This comparison is not appropriate, as one reflects the physical properties of the fault, and the other reflects the kinematics of the fault. Much of the uncertainty about slip behavior is because in subduction zones, the shallow part of the fault is far from land, and therefore creep is not detectable by land-based GPS. Published coupling maps of subduction zone megathrusts often assume a low coupling ratio near the trench, updip from fully coupled regions. Yet, if the megathrust attains a coupling ratio of 1 anywhere on the fault (i.e., the hanging wall is moving with the same velocity as the footwall), a lower value of coupling updip of this location requires interseismic extension at a rate proportional to the decrease (Wang and Dixon, 2004). We argue that the shallow region of megathrusts lie in updip stress shadows, and do not (except under rare circumstances) experience appropriate driving forces to cause significant creep during the interseismic period. Therefore it may not be possible to determine whether these regions are frictionally locked by examining interseismic geodetic records. We demonstrate this effect using a boundary element model with rate-strengthening friction and a simplified subduction zone geometry. We show that a coupling value of zero at the trench is physically unrealistic even if only a small portion of the downdip fault zone is locked. The maximum creep at the trench depends on the width of the transition of the frictionally locked zone, but should be small (<30% of plate rate) under most circumstances. During the interseismic period, even if the shallow parts of megathrusts are frictionally unlocked, creep is likely smaller than the resolution of current seafloor geodetic techniques (which is currently in the range of cms/yr). These results have important implications for various aspects of subduction studies, including physical limits on geodetic coupling inversions, the hazard posed by slip on shallow decollements (tsunamigenic or otherwise), the seismotectonic interpretation of shallow seismicity, and the utility of seafloor geodetic measurements.

Ground motion models used in the 2014 U.S. National Seismic Hazard Maps

USGS Publications Warehouse

Rezaeian, Sanaz; Petersen, Mark D.; Moschetti, Morgan P.

2015-01-01

The National Seismic Hazard Maps (NSHMs) are an important component of seismic design regulations in the United States. This paper compares hazard using the new suite of ground motion models (GMMs) relative to hazard using the suite of GMMs applied in the previous version of the maps. The new source characterization models are used for both cases. A previous paper (Rezaeian et al. 2014) discussed the five NGA-West2 GMMs used for shallow crustal earthquakes in the Western United States (WUS), which are also summarized here. Our focus in this paper is on GMMs for earthquakes in stable continental regions in the Central and Eastern United States (CEUS), as well as subduction interface and deep intraslab earthquakes. We consider building code hazard levels for peak ground acceleration (PGA), 0.2-s, and 1.0-s spectral accelerations (SAs) on uniform firm-rock site conditions. The GMM modifications in the updated version of the maps created changes in hazard within 5% to 20% in WUS; decreases within 5% to 20% in CEUS; changes within 5% to 15% for subduction interface earthquakes; and changes involving decreases of up to 50% and increases of up to 30% for deep intraslab earthquakes for most U.S. sites. These modifications were combined with changes resulting from modifications in the source characterization models to obtain the new hazard maps.

Seismic structure of subducted Philippine Sea plate beneath the southern Ryukyu arc by receiver function and local earthquakes tomography

NASA Astrophysics Data System (ADS)

Nakamura, M.

2012-12-01

Seismic coupling of the Ryukyu subduction zone is assumed to be weak from the lack of historical interplate large earthquakes. However, recent investigation of repeating slow slip events (Heki & Kataoka, 2008), shallow low frequency earthquakes (Ando et al., 2012), and source of 1771 Yaeyama mega-tsunami (Nakamura, 2009), showed that the interplate coupling is not weak in the south of Ryukyu Trench. The biannually repeating SSEs (Mw=6.5) occur at the depth of 20-40 km on the upper interface of the subducted Philippine Sea plate beneath Yaeyama region, where earthquake swarm occurred on 1991 and 1992. To reveal the relation among the crustal structure, earthquake swarms, and occurrence of slow slip events (SSE), local earthquake tomography and receiver function (RF) analysis was computed in the southwestern Ryukyu arc. A tomographic inversion was used to determine P and S wave structures beneath Iriomote Island in the southwestern Ryukyu region for comparison with the locations of the SSE. The seismic tomography (Thurber & Eberhart-Phillips, 1999) was employed. The P- and S- wave arrival time data picked manually by Japan Meteorological Agency (JMA) are used. The 6750 earthquakes from January 2000 to July 2012 were used. For the calculation of the receiver function, the 212 earthquakes whose magnitudes are over 6.0 and epicentral distances are between 30 and 90 degrees were selected. The teleseicmic waveforms observed at two short-period seismometers of the JMA, and one broadband seismometer of F-net of National Research Institute for Earth Science and Disaster Prevention were used. The water level method (the water level is 0.01) is applied to original waveforms. Assuming that each later phase in a RF is the wave converted from P to S at a depth, I transformed the time domain RF into the depth domain one along each ray path in a reference velocity model. The JMA2001 velocity model is used in this study. The results of tomography show that the low Vp and high Vp/Vs anomalies are distributed along the hypocenters in the subducted slab. The plate interface is about 10 km above the slab earthquakes from the trace of negative RF amplitude. The slab earthquakes are distributed along the trace of positive RF amplitude. Therefore the slab earthquakes occur near the oceanic Moho of the PHS. The fault depth of the SSEs corresponds to the plate interface within 5 km. The fault-planes of the SSE are located above the low Vp and high Vp/Vs zone. Assuming that the difference between high Vp/Vs and low Vp/Vs originates to the fluid contents, this would be interpreted that the fluids from the subducted oceanic crust cannot be transported upward and is trapped at the plate interface. The observed strong S-wave reflector (Nakamura, 2001) in the upper interface of the subducted plate also supports the idea. The top of the faults of the SSEs connects to the cluster of earthquake swarms in the lower crust. This suggests that the trapped fluids are transported upward along the faults, accumulates in the lower crust, and induce the swarm of micro-earthquakes in the lower crust.

Stress analysis of the Mw 7.4 Armería, Colima, Mexico earthquake of 22 January 2003

NASA Astrophysics Data System (ADS)

Vargas-Bracamontes, D.; Nunez-Cornu, F. J.

2012-12-01

On 22 January 2003 a shallow Mw 7.4 earthquake occurred off the Pacific coast of the state of Colima. This event struck near the towns of Tecomán and Armería in western Mexico where a diffuse triple junction between the North American, Cocos and Rivera plates makes the local tectonic setting highly complex. This earthquake is the largest during the twenty-first century in the area. Some seismic studies of this earthquake indicate that this event occurred on a continental intraplate reverse fault, suggesting that the shock and its aftershocks represent partial accommodation of deformation in the continental crust caused by oblique subduction. In contrast, other works propose that the 2003 Armería earthquake was due to faulting along the subduction interface between the Rivera and North American Plates. We assess the suggested sources of this earthquake in terms of stress models that consider the controversial geometrical features that characterize this tectonic area. Also, we explore the implications for seismic hazard that this event could have caused in the Colima region.

How fault geometry controls earthquake magnitude

NASA Astrophysics Data System (ADS)

Bletery, Q.; Thomas, A.; Karlstrom, L.; Rempel, A. W.; Sladen, A.; De Barros, L.

2016-12-01

Recent large megathrust earthquakes, such as the Mw9.3 Sumatra-Andaman earthquake in 2004 and the Mw9.0 Tohoku-Oki earthquake in 2011, astonished the scientific community. The first event occurred in a relatively low-convergence-rate subduction zone where events of its size were unexpected. The second event involved 60 m of shallow slip in a region thought to be aseismicaly creeping and hence incapable of hosting very large magnitude earthquakes. These earthquakes highlight gaps in our understanding of mega-earthquake rupture processes and the factors controlling their global distribution. Here we show that gradients in dip angle exert a primary control on mega-earthquake occurrence. We calculate the curvature along the major subduction zones of the world and show that past mega-earthquakes occurred on flat (low-curvature) interfaces. A simplified analytic model demonstrates that shear strength heterogeneity increases with curvature. Stress loading on flat megathrusts is more homogeneous and hence more likely to be released simultaneously over large areas than on highly-curved faults. Therefore, the absence of asperities on large faults might counter-intuitively be a source of higher hazard.

Along-Strike Electrical Conductivity Variations in the Incoming Plate and Shallow Forearc of the Cascadia Subduction Zone

NASA Astrophysics Data System (ADS)

Key, K.; Bedrosian, P.; Egbert, G. D.; Livelybrooks, D.; Parris, B. A.; Schultz, A.

2015-12-01

The Magnetotelluric Observations of Cascadia using a Huge Array (MOCHA) experiment was carried out to study the nature of the seismogenic locked zone and the down-dip transition zone where episodic tremor and slip (ETS) originates. This amphibious magnetotelluric (MT) data set consists of 8 offshore and 15 onshore profiles crossing from just seaward of the trench to the western front of the Cascades, with a north-south extent spanning from central Oregon to central Washington. The 71 offshore stations and the 75 onshore stations (red triangles in the image below) fit into the broader context of the more sparsely sampled EarthScope MT transportable array (black triangles) and other previous and pending MT surveys (other symbols). These data allows us to image variations in electrical conductivity along distinct segments of the Cascadia subduction zone defined by ETS recurrence intervals. Since bulk conductivity in this setting depends primarily on porosity, fluid content and temperature, the conductivity images created from the MOCHA data offer unique insights on fluid processes in the crust and mantle, and how the distribution of fluid along the plate interface relates to observed variations in ETS behavior. This abstract explores the across- and along-strike variations in the incoming plate and the shallow offshore forearc. In particular we examine how conductivity variations, and the inferred fluid content and porosity variations, are related to tectonic segmentation, seismicity and deformation patterns, and arc magma variations along-strike. Porosity inferred in the forearc crust can be interpreted in conjunction with active and passive seismic imaging results and may provide new insights on the origin of recently observed extremely high heat flow values. A companion abstract (Parris et al.) examines the deeper conductivity structure of the locked and ETS zones along the plate interface in order to identify correlations between ETS occurrence rates and inferred fluid concentrations.

Slab-pull and slab-push earthquakes in the Mexican, Chilean and Peruvian subduction zones

NASA Astrophysics Data System (ADS)

Lemoine, A.; Madariaga, R.; Campos, J.

2002-09-01

We studied intermediate depth earthquakes in the Chile, Peru and Mexican subduction zones, paying special attention to slab-push (down-dip compression) and slab-pull (down-dip extension) mechanisms. Although, slab-push events are relatively rare in comparison with slab-pull earthquakes, quite a few have occurred recently. In Peru, a couple slab-push events occurred in 1991 and one slab-pull together with several slab-push events occurred in 1970 near Chimbote. In Mexico, several slab-push and slab-pull events occurred near Zihuatanejo below the fault zone of the 1985 Michoacan event. In central Chile, a large M=7.1 slab-push event occurred in October 1997 that followed a series of four shallow Mw>6 thrust earthquakes on the plate interface. We used teleseismic body waveform inversion of a number of Mw>5.9 slab-push and slab-pull earthquakes in order to obtain accurate mechanisms, depths and source time functions. We used a master event method in order to get relative locations. We discussed the occurrence of the relatively rare slab-push events in the three subduction zones. Were they due to the geometry of the subduction that produces flexure inside the downgoing slab, or were they produced by stress transfer during the earthquake cycle? Stress transfer can not explain the occurence of several compressional and extensional intraplate intermediate depth earthquakes in central Chile, central Mexico and central Peru. It seemed that the heterogeneity of the stress field produced by complex slab geometry has an important influence on intraplate intermediate depth earthquakes.

Slow slip events in the early part of the earthquake cycle

NASA Astrophysics Data System (ADS)

Voss, Nicholas K.; Malservisi, Rocco; Dixon, Timothy H.; Protti, Marino

2017-08-01

In February 2014 a Mw = 7.0 slow slip event (SSE) took place beneath the Nicoya Peninsula, Costa Rica. This event occurred 17 months after the 5 September 2012, Mw = 7.6, earthquake and along the same subduction zone segment, during a period when significant postseismic deformation was ongoing. A second SSE occurred in the middle of 2015, 21 months after the 2014 SSE and 38 months after the earthquake. The recurrence interval for Nicoya SSEs was unchanged by the earthquake. However, the spatial distribution of slip for the 2014 event differed significantly from previous events, having only deep ( 40 km) slip, compared to previous events, which had both deep and shallow slip. The 2015 SSE marked a return to the combination of deep plus shallow slip of preearthquake SSEs. However, slip magnitude in 2015 was nearly twice as large (Mw = 7.2) as preearthquake SSEs. We employ Coulomb Failure Stress change modeling in order to explain these changes. Stress changes associated with the earthquake and afterslip were highest near the shallow portion of the megathrust, where preearthquake SSEs had significant slip. Lower stress change occurred on the deeper parts of the plate interface, perhaps explaining why the deep ( 40 km) region for SSEs remained unchanged. The large amount of shallow slip in the 2015 SSE may reflect lack of shallow slip in the prior SSE. These observations highlight the variability of aseismic strain release rates throughout the earthquake cycle.